Quantifying Delocalization and Static Correlation Errors by Imposing (Spin)Population Redistributions through Constraints on Atomic Domains

Spin migration in density functional theory: Energy, potential, and density perspectives

Spin is a fundamental property of any many-electron system. The ability of density functional theory to accurately predict the physical properties of a system, while varying its spin, is crucial for describing magnetic materials and high-spin molecules, spin flips, and magnetization and demagnetization processes. Within density functional theory, when using various exchange-correlation approximations, the exact dependence of the energy on the spin often deviates from the exact constant or piecewise-linear behavior, which is directly related to the problem of strong (static) correlation and challenges the description of molecular dissociation. In this paper, we study the behavior of the energy, the frontier Kohn-Sham (KS) and generalized KS (GKS) orbitals, the KS potentials, and the electron density, with respect to fractional spin, in different atomic systems. We analyze seven standard exchange-correlation functionals and find two main scenarios of deviation from the expected exact results. We clearly recognize a jump in the frontier orbital energies upon spin variation in the exact exchange and in hybrid functionals, as well as the related plateau in the corresponding KS potential, when using the optimized effective potential method within the KS scheme. When calculations are performed using the GKS approach, no jumps are observed, as expected. Moreover, we demonstrate that for high-spin systems, a full three-dimensional treatment is crucial; the spherical approximation commonly used in atoms causes a qualitative deviation. Our results are instrumental for the assessment of the quality of existing approximations from a new perspective and for the development of advanced functionals with sensitivity to magnetic properties.

Tilted-Plane Structure of the Energy of Finite Quantum Systems

The piecewise linearity condition on the total energy with respect to the total magnetization of finite quantum systems is derived using the infinite-separation-limit technique. This generalizes the well-known constancy condition, related to static correlation error, in approximate density functional theory. The magnetic analog of Koopmans' theorem in density functional theory is also derived. Moving to fractional electron count, the tilted-plane condition is derived, lifting certain assumptions in previous works. This generalization of the flat-plane condition characterizes the total energy surface of a finite system for all values of electron count N and magnetization M. This result is used in combination with tabulated spectroscopic data to show the flat-plane structure of the oxygen atom, among others. We find that derivative discontinuities with respect to electron count sometimes occur at noninteger values. A diverse set of tilted-plane structures is shown to occur in d-orbital subspaces, depending on chemical coordination. General occupancy-based total-energy expressions are demonstrated thereby to be necessarily dependent on the symmetry-imposed degeneracies.

Extending Conceptual Density Functional Theory toward First-Order Reduced Density Matrices: An Open Subsystems Viewpoint on the Fukui Matrix

As a matrix extension of the Fukui function, a reactivity descriptor grounded within Conceptual Density Functional Theory, the Fukui matrix extends Frontier Molecular Orbital Theory to correlated regimes with its eigendecomposition in Fukui occupations and Fukui naturals. Despite successful applications, the questions remain as to whether replacing a quantity derived from a purely density-based framework by its matrix extension is theoretically well-founded and what chemical information is contained in the corresponding eigendecomposition. In this study, we show that the matrix extension of the Fukui function is only well-defined if one also generalizes the external potential to become nonlocal, leading to the introduction of Conceptual First-Order Reduced Density Matrix Functional Theory. By interpreting the Anderson impurity model from an interacting open subsystem perspective, we show how Fukui occupations and Fukui naturals reflect the influence of an increasing (static) correlation and which characteristic patterns we should expect within a molecular context. This study represents a step in generalizing Conceptual Density Functional Theory beyond its density-based perspective.

Uncovering phase transitions that underpin the flat-planes in the tilted Hubbard model using subsystems and entanglement measures

The failure of many approximate electronic structure methods can be traced to their erroneous description of fractional charge and spin redistributions in the asymptotic limit toward infinity, where violations of the flat-plane conditions lead to delocalization and static correlation errors. Although the energetic consequences of the flat-planes are known, the underlying quantum phase transitions that occur when (spin)charge is redistributed have not been characterized. In this study, we use open subsystems to redistribute (spin)charges in the tilted Hubbard model by imposing suitable Lagrange constraints on the Hamiltonian. We computationally recover the flat-plane conditions and quantify the underlying quantum phase transitions using quantum entanglement measures. The resulting entanglement patterns quantify the phase transition that gives rise to the flat-plane conditions and quantify the complexity required to accurately describe charge redistributions in strongly correlated systems. Our study indicates that entanglement patterns can uncover those phase transitions that have to be modeled accurately if the delocalization and static correlation errors of approximate methods are to be reduced.

Analyzing the Behavior of Spin Phases in External Magnetic Fields by Means of Spin-Constrained States

During molecular dissociation in the presence of an external uniform magnetic field, electrons flip their spin antiparallel to the magnetic field because of the stabilizing influence of the spin Zeeman operator. Although generalized Hartree-Fock descriptions furnish the optimal mean-field energetic description of such bond-breaking processes, they are allowed to break Ŝ_z symmetry, leading to intricate and unexpected spin phases and phase transitions. In this work, we show that the behavior of these molecular spin phases can be interpreted in terms of spin phase diagrams constructed by constraining states to target expectation values of projected spin. The underlying constrained states offer a complete electronic characterization of the spin phases and spin phase transitions, as they can be analyzed using standard quantum chemical tools. Because the constrained states effectively span the entire phase space, they could provide an excellent starting point for post-Hartree-Fock methods aimed at gaining more electron correlation or regaining spin symmetry.