Where  represent the reduced density matrix of the entire system, X  represent the reduced density matrix of either the active or environment subsystem, h represent the core Hamiltonian of the whole system, X G    represent the two-electron interactions within the active or environment subsystem, and ,