“…In order to redistribute charges and spins over the open subsystem and the reservoir, we will constrain the ground state wave function of the total system to attain a particular feature of the open system (in this study, those features will be Mulliken (spin)populations of atomic domains). The framework needed to impose such constraints was initially developed by Mukherji and Karplus in the context of observables like the dipole moment in unrestricted variational theories and has been significantly expanded upon since then, ,− culminating in modern day applications such as, among others, constrained density functional theory, electronegativity equalization, ,,, and restricted open-shell Hartree–Fock reformulated as a constrained unrestricted Hartree–Fock model . In order to make this paper self-contained, we will recall the necessary theoretical concepts needed to construct the framework that will allow us to constrain the (spin)population of such an open subsystem.…”