DanceQ: High-performance library for number-conserving bases

Abstract

The complexity of quantum many-body problems scales exponentially with the size of the system, rendering any finite size scaling analysis a formidable challenge. This is particularly true for methods based on the full representation of the wave function, where one simply accepts the enormous Hilbert space dimensions and performs linear algebra operations, e.g., for finding the ground state of the Hamiltonian. If the system satisfies an underlying symmetry where an operator with degenerate spectrum commutes with the Hamiltonian, it can be block-diagonalized, thus reducing the complexity at the expense of additional bookkeeping. At the most basic level, required for Krylov space techniques (like the Lanczos algorithm) it is necessary to implement a matrix-vector product of a block of the Hamiltonian with arbitrary block-wavefunctions, potentially without holding the Hamiltonian block in memory. An efficient implementation of this operation requires the calculation of the position of an arbitrary basis vector in the canonical ordering of the basis of the block. We present here an elegant and powerful, multi-dimensional approach to this problem for the $U(1)$ symmetry appearing in problems with particle number conservation. Our divide-and-conquer algorithm uses multiple subsystems and hence generalizes previous approaches to make them scalable. In addition to the theoretical presentation of our algorithm, we provide DanceQ, a flexible and modern - header only - C++20 implementation to manipulate, enumerate, and map to its index any basis state in a given particle number sector as open source software under https://DanceQ.gitlab.io/danceq.

Source code available under https://gitlab.com/DanceQ/danceq, documentation here https://DanceQ.gitlab.io/danceq

References in corpus (21)

• Anyons in an exactly solved model and beyond
• The density-matrix renormalization group in the age of matrix product states
• Ultracold atomic gases in optical lattices: mimicking condensed matter physics and beyond
• Many-body localization edge in the random-field Heisenberg chain
• The ITensor Software Library for Tensor Network Calculations
• Pyrochlore Photons: The U(1) Spin Liquid in a S=1/2 Three-Dimensional Frustrated Magnet
• Massive Parallel Quantum Computer Simulator
• Polynomially filtered exact diagonalization approach to many-body localization
• Exact diagonalization: the Bose-Hubbard model as an example
• Massively parallel quantum computer simulator, eleven years later
• The $S=1/2$ Kagome Heisenberg Antiferromagnet Revisited
• Quantum Lattice Model Solver $\mathcal{H}\Phi$
• Improving Hamiltonian encodings with the Gray code
• The pyrochlore S=1/2 Heisenberg antiferromagnet at finite temperature
• A novel general mapping for bosonic and fermionic operators in Fock space
• Sublattice Coding Algorithm and Distributed Memory Parallelization for Large-Scale Exact Diagonalizations of Quantum Many-Body Systems
• Optimized implementation of the Lanczos method for magnetic systems
• Abundance of hard-hexagon crystals in the quantum pyrochlore antiferromagnet
• Polynomial filter diagonalization of large Floquet unitary operators
• Divide and conquer the Hilbert space of translation-symmetric spin systems
• Topological States of Matter in Frustrated Quantum Magnetism

Cited by in corpus (2)

• Probing Hilbert space fragmentation using controlled dephasing
• Thermodynamics of the Heisenberg antiferromagnet on the maple-leaf lattice