Volume 55, Number 3, May-June 2021
|Page(s)||763 - 787|
|Published online||05 May 2021|
Numerical solution of large scale Hartree–Fock–Bogoliubov equations
Department of Mathematics, University of California, Berkeley, CA 94720, USA
2 Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
* Corresponding author: email@example.com
Accepted: 2 November 2020
The Hartree–Fock–Bogoliubov (HFB) theory is the starting point for treating superconducting systems. However, the computational cost for solving large scale HFB equations can be much larger than that of the Hartree–Fock equations, particularly when the Hamiltonian matrix is sparse, and the number of electrons N is relatively small compared to the matrix size Nb. We first provide a concise and relatively self-contained review of the HFB theory for general finite sized quantum systems, with special focus on the treatment of spin symmetries from a linear algebra perspective. We then demonstrate that the pole expansion and selected inversion (PEXSI) method can be particularly well suited for solving large scale HFB equations. For a Hubbard-type Hamiltonian, the cost of PEXSI is at most 𝒪(Nb2) for both gapped and gapless systems, which can be significantly faster than the standard cubic scaling diagonalization methods. We show that PEXSI can solve a two-dimensional Hubbard-Hofstadter model with Nb up to 2.88 × 106, and the wall clock time is less than 100 s using 17 280 CPU cores. This enables the simulation of physical systems under experimentally realizable magnetic fields, which cannot be otherwise simulated with smaller systems.
Mathematics Subject Classification: 65 / 81
Key words: Hartree–Fock–Bogoliubov / pole expansion and selected inversion / superconductivity / Hubbard-Hofstadter
© EDP Sciences, SMAI 2021
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.