Particle-Breaking Hartree-Fock Theory for Open Molecular Systems

In this work we present the particle-breaking Hartree-Fock (PBHF) model which is a mean-field approach to open molecular systems. The interaction of a system with the environment is parametrized through a particle-breaking term in the molecular Hamiltonian. The PBHF wave function is constructed through an exponential unitary transformation of a Slater determinant with a given number of electrons. We consider only the closed-shell formalism. The parametrization results in a linear combination of Slater determinants with different numbers of electrons, i.e., the PBHF wave function is not an eigenfunction of the number operator. As a result, the density matrix may have fractional occupations which are, because of the unitary parametrization, always between 0.0 and 2.0. The occupations are optimized simultaneously with the orbitals, using the trust-region optimization procedure. In the limit of a particle-conserving Hamiltonian, the PBHF optimization will converge to a standard Hartree-Fock wave function. We show that the average number of electrons may be decreased or increased depending on whether the particle-breaking term affects occupied or virtual orbitals.