The self-consistent field method utilized for solving the Hartree-Fock (HF) problem and the closely related Kohn-Sham problem is typically thought of as one of the cheapest methods available to quantum chemists. This intuition has been developed from the numerous applications of the self-consistent field method to a large variety of molecular systems. However, as characterized by its worst-case behavior, the HF problem is NP-complete. In this work, we map out boundaries of the NP-completeness by investigating restricted instances of HF. We have constructed two new NP-complete variants of the problem. The first is a set of Hamiltonians whose translationally invariant Hartree-Fock solutions are trivial, but whose broken symmetry solutions are NP-complete. Second, we demonstrate how to embed instances of spin glasses into translationally invariant Hartree-Fock instances and provide a numerical example. These findings are the first steps towards understanding in which cases the self-consistent field method is computationally feasible and when it is not.
The self-consistent field method utilized for solving the Hartree-Fock (HF) problem and the closely related Kohn-Sham problem is typically thought of as one of the cheapest methods available to quantum chemists. This intuition has been developed from the numerous applications of the self-consistent field method to a large variety of molecular systems. However, as characterized by its worst-case behavior, the HF problem is NP-complete. In this work, we map out boundaries of the NP-completeness by investigating restricted instances of HF. We have constructed two new NP-complete variants of the problem. The first is a set of Hamiltonians whose translationally invariant Hartree-Fock solutions are trivial, but whose broken symmetry solutions are NP-complete. Second, we demonstrate how to embed instances of spin glasses into translationally invariant Hartree-Fock instances and provide a numerical example. These findings are the first steps towards understanding in which cases the self-consistent field method is computationally feasible and when it is not.
Many problems in physical chemistry involve systems that are coupled to an environment, such as a molecule interacting with a nearby surface. In these cases, electron density may flow between the system and its surroundings, leading to metastable states and other open-system phenomena that cannot be captured by conventional Hermitian quantum mechanics. Non-Hermitian quantum mechanics (NHQM) [1] provides a natural framework for the description of open systems. Hartree-Fock (HF) and Kohn-Sham (KS) methods can be extended to the non-Hermitian domain [2–7]. Here, we show that even for nominally isolated systems, within HF and KS, each electron can be viewed as coupled to the remaining electrons, which serve as an environment. This viewpoint reveals that an electron can leave or rejoin the system, producing self-consistent complex potentials. These unconventional solutions are absent in the conventional closed-system HF and KS approach. We show that these additional solutions can be physically meaningful. Our results highlight the importance of incorporating non-Hermitian perspectives into electronic structure calculations, offering new insights and extending the range of problems HF and KS can be applied to. [1] Moiseyev, N. Non-Hermitian Quantum Mechanics, Cambridge University Press (2011). [2] Zhou, Y., Ernzerhof, M. J. Chem. Phys. 136, 094105 (2012). [3] Ernzerhof, M. J. Chem. Phys. 125, 124104 (2006). [4] Zhou, Y., Ernzerhof, M. J. Phys. Chem. Lett. 3, 1916–1920 (2012). [5] Burton, H. G. A., Thom, A. J. W. J. Chem. Theory Comput. 12, 167–173 (2016). [6] Burton, H. G. A., Gross, M., Thom, A. J. W. J. Chem. Theory Comput. 14, 607–618 (2018). [7] Burton, H. G. A., Thom, A. J. W., Loos, P.-F. J. Chem. Theory Comput. 15, 4374–4385 (2019).
Many problems in physical chemistry involve systems that are coupled to an environment, such as a molecule interacting with an adjacent surface, possibly resulting in meta-stable molecular states where electron density is transferred to the surface. Such systems can be described by non-Hermitian quantum mechanics (NHQM), where the Hamiltonian includes dissipative terms. Within NHQM, one can also formulate the Hartree-Fock (HF) and Kohn-Sham (KS) methods and, as in the conventional theory, an effective independent-particle picture is employed. The crucial observation of the present work is that even for systems that are not coupled to an environment, in the HF or KS equation a single electron is coupled to a bath of the remaining electrons which can act as an environment, opening up the possibility for the exchange of current density between the one-electron and the remaining N-1 electron system. The corresponding self-consistent states represent a new uncharted space of solutions to the HF and KS equations. We show that the additional solutions can have a physical interpretation and thus extend the range of problems HF and KS can be applied to. If open-system HF and KS calculations are performed, the new class of solutions is always encountered but this has also not been noted previously.
A simple and commonly employed approximate technique with which one can examine spatially disordered systems when strong electronic correlations are present is based on the use of real-space unrestricted self-consistent Hartree-Fock wave functions. In such an approach the disorder is treated exactly while the correlations are treated approximately. In this report we critique the success of this approximation by making comparisons between such solutions and the exact wave functions for the Anderson-Hubbard model. Due to the sizes of the complete Hilbert spaces for these problems, the comparisons are restricted to small one-dimensional chains, up to ten sites, and a 4x4 two-dimensional cluster, and at 1/2 filling these Hilbert spaces contain about 63,500 and 166 million states, respectively. We have completed these calculations both at and away from 1/2 filling. This approximation is based on a variational approach which minimizes the Hartree-Fock energy, and we have completed comparisons of the exact and Hartree-Fock energies. However, in order to assess the success of this approximation in reproducing ground-state correlations we have completed comparisons of the local charge and spin correlations, including the calculation of the overlap of the Hartree-Fock wave functions with those of the exact solutions. We find that this approximation reproduces the local charge densities to quite a high accuracy, but that the local spin correlations, as represented by < S_i . S_j >, are not as well represented. In addition to these comparisons, we discuss the properties of the spin degrees of freedom in the HF approximation, and where in the disorder-interaction phase diagram such physics may be important.
Abstract. A simple and commonly employed approximate technique with which one can examine spatially disordered systems when strong electronic correlations are present is based on the use of real-space unrestricted self-consistent Hartree-Fock wave functions. In such an approach the disorder is treated exactly while the correlations are treated approximately. In this report we critique the success of this approximation by making comparisons between such solutions and the exact wave functions for the Anderson-Hubbard model. Due to the sizes of the complete Hilbert spaces for these problems, the comparisons are restricted to small one-dimensional chains, up to ten sites, and a 4x4 two-dimensional cluster, and at 1/2 filling these Hilbert spaces contain about 63,500 and 166 million states, respectively. We have completed these calculations both at and away from 1/2 filling. This approximation is based on a variational approach which minimizes the Hartree-Fock energy, and we have completed comparisons of the exact and Hartree-Fock energies. However, in order to assess the success of this approximation in reproducing ground-state correlations we have completed comparisons of the local charge and spin correlations, including the calculation of the overlap of the Hartree-Fock wave functions with those of the exact solutions. We find that this approximation reproduces the local charge densities to quite a high accuracy, but that the local spin correlations, as represented by 〈Si · Sj〉, are not as well represented. In addition to these comparisons, we discuss the properties of the spin degrees of freedom in the HF approximation, and where in the disorderinteraction phase diagram such physics may be important. PACS numbers: 71.27.+a,71.23.-k,71.10.-wCritiquing HF for the Anderson-Hubbard Model 2 1.
In quantum chemistry, self-consistent field (SCF) algorithms define a nonlinear optimization problem, with both continuous and discrete components. In this work, we derive Hartree–Fock-inspired SCF algorithms that can be exactly written as a sequence of Quadratic Unconstrained Spin/Binary Optimization problems (QUSO/QUBO). We reformulate the optimization problem as a series of MaxCut graph problems, which can be efficiently solved using semidefinite programming techniques. This procedure provides performance guarantees at each SCF step, irrespective of the complexity of the optimization landscape. We numerically demonstrate the QUBO-SCF and MaxCut-SCF methods by studying the hydroxide anion OH – and molecular Nitrogen N 2 . The largest problem addressed in this study involves a system comprised of 220 qubits (equivalently, spin–orbitals). Our results show that QUBO-SCF and MaxCut-SCF suffer much less from internal instabilities compared with conventional SCF calculations. Additionally, we show that the new SCF algorithms can enhance single-reference methods, such as configuration interaction. Finally, we explore how quantum algorithms for optimization can be applied to the QUSO problems arising from the Hartree–Fock method. Four distinct hybrid-quantum classical approaches are introduced: GAS-SCF, QAOA-SCF, QA-SCF and DQI-SCF.
I will present our recent results on the calculation of symmetry-projected wave functions [1,2] for electronic structure theory. For Hartree-Fock and spin projection, this is a 50-year old problem in quantum chemistry, going all the way back to Lowdin and his “extended HF” theory. For Hartree-Fock-Bogoliubov and number projection, our approach offers new perspectives on Antisymmetrized Geminal Power (AGP) wavefunctions that were the focus of much attention in the 1980s. In our work, all molecular symmetries (electron number, spin S 2 and Sz , point group, and complex conjugation) are deliberately broken and restored in a self-consistent variation-after-projection approach. The resulting method yields a comprehensive black-box treatment of static correlation with one-electron (mean-field) computational cost. The ensuing wave function is of high quality multireference character competitive with exact diagonalization. The method can also be applied to calculate excited states and spectral functions [3]. The curse of the thermodynamic limit and the quest for a low-cost treatment of residual correlations will also be addressed.