Efficient simulations of Hartree--Fock equations by an accelerated gradient descent method
arXiv:2402.17843 · doi:10.1103/PhysRevE.110.055304
Abstract
We develop convergence acceleration procedures that enable a gradient descent-type iteration method to efficiently simulate Hartree--Fock equations for atoms interacting both with each other and with an external potential. Our development focuses on three aspects: (i) optimization of a parameter in the preconditioning operator; (ii) adoption of a technique that eliminates the slowest-decaying mode to the case of many equations (describing many atoms); and (iii) a novel extension of the above technique that allows one to eliminate multiple modes simultaneously. We illustrate performance of the numerical method for the 2D model of the first layer of helium atoms above a graphene sheet. We demonstrate that incorporation of aspects (i) and (ii) above into the ``plain" gradient descent method accelerates it by at least two orders of magnitude, and often by much more. Aspect (iii) -- a multiple-mode elimination -- may bring further improvement to the convergence rate compared to aspect (ii), the single-mode elimination. Both single- and multiple-mode elimination techniques are shown to significantly outperform the well-known Anderson Acceleration. We believe that our acceleration techniques can also be gainfully employed by other numerical methods, especially those handling hard-core-type interaction potentials.
arXiv commentmain text (39 pages); supplement appended (7 pages)
References in corpus (23)
- Supersolid hardcore bosons on the triangular lattice
- Phase diagram of 4He adsorbed on graphite
- 4He on a single graphene sheet
- Obtaining Wannier Functions of a Crystalline Insulator within a Hartree-Fock approach: Applications to LiF and LiCl
- Quantum phases of hardcore bosons with long-range interactions on a square lattice
- Monolayer Solid Helium-4 Clusters on Graphite
- Phase diagram of He on graphene
- Multiconfigurational Hartree-Fock theory for identical bosons in a double well
- Possible quantum liquid crystal phases of helium monolayers
- Spatially Modulated Superfluid State in Two-Dimensional He Films
- A mode elimination technique to improve convergence of iteration methods for finding solitary waves
- Iterative approaches to the self-consistent nuclear energy density functional problem. Heavy ball dynamics and potential preconditioning
- Superfluid and supersolid phases of 4He on the second layer of graphite
- Diffusion Monte Carlo Study of the O Adsorption on a Single Layer Graphene
- Two-Dimensional Bose-Hubbard Model for Helium on Graphene
- Conjugate gradient method for finding fundamental solitary waves
- Hard-core boson approach to the spin-1/2 triangular-lattice antiferromagnet CsCuCl at finite temperatures in magnetic fields higher than the saturation field
- Novel substrates for Helium adsorption: Graphane and Graphene-Fluoride
- Commensurate-incommensurate solid transition in the He monolayer on -graphyne
- Strain-induced superfluid transition for atoms on graphene
- Adsorption of Helium atoms on two-dimensional substrates
- Phase Stability in the Two dimensional Anisotropic Boson Hubbard Hamiltonian
- Imaginary-time time-dependent density functional theory for periodic systems