Numerical solution of large scale Hartree-Fock-Bogoliubov equations

Abstract

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 $N_b$. 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 $\Or(N_b^2)$ 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 $N_b$ up to $2.88\times 10^6$, and the wall clock time is less than $100$ s using $17280$ CPU cores. This enables the simulation of physical systems under experimentally realizable magnetic fields, which cannot be otherwise simulated with smaller systems.

arXiv comment25 pages, 9 figures

References in corpus (18)

• Realization of the Hofstadter Hamiltonian with ultracold atoms in optical lattices
• From the Kondo Regime to the Mixed-Valence Regime in a Single-Electron Transistor
• Density matrix embedding: A simple alternative to dynamical mean-field theory
• Ground-state phase diagram of the square lattice Hubbard model from density matrix embedding theory
• Chebyshev-BdG: an efficient numerical approach to inhomogeneous superconductivity
• Accelerating Atomic Orbital-based Electronic Structure Calculation via Pole Expansion and Selected Inversion
• Pole-based approximation of Fermi-Dirac function
• Efficient Numerical Self-consistent Mean-field Approach for Fermionic Many-body Systems by Polynomial Expansion on Spectral Density
• Majorana corner modes with solitons in an attractive Hubbard-Hofstadter model of cold atom optical lattices
• Variational Principle of Bogoliubov and Generalized Mean Fields in Many-Particle Interacting Systems
• Minimizers for the Hartree-Fock-Bogoliubov Theory of Neutron Stars and White Dwarfs
• Many-body computations by stochastic sampling in Hartree-Fock-Bogoliubov space
• Topological Phase Transition in the Hofstadter-Hubbard Model
• Minimax rational approximation of the Fermi-Dirac distribution
• Phases of attractive spin-imbalanced fermions in square lattices
• Attractive Hofstadter-Hubbard model with imbalanced chemical and vector potentials
• BCS theory of time-reversal-symmetric Hofstadter-Hubbard model
• Proximity induced pseudogap in mesoscopic superconductor/normal-metal bilayers