Wavefunction optimization at the complete basis set limit with Multiwavelets and DMRG

Abstract

The density matrix renormalization group (DMRG) is a powerful numerical technique to solve strongly correlated quantum systems: it deals well with systems which are not dominated by a single configuration (unlike Coupled Cluster) and it converges rapidly to the Full Configuration Interaction (FCI) limit (unlike truncated Configuration Interaction (CI) expansions). In this work, we develop an algorithm integrating DMRG within the multiwavelet-based multiresolution analysis (MRA). Unlike fixed basis sets, multiwavelets offer an adaptive and hierarchical representation of functions, approaching the complete basis set limit to a specified precision. As a result, this combined technique leverages the multireference capability of DMRG and the complete basis set limit of MRA and multiwavelets. More specifically, we adopt a pre-existing Lagrangian optimization algorithm for orbitals represented in the MRA domain and improve its computational efficiency by replacing the original CI calculations with DMRG. Additionally, we substitute the reduced density matrices computation with the direct extraction of energy gradients from the DMRG tensors. We apply our method to small systems such H2, He, HeH2, BeH2 and N2. The results demonstrate that our approach reduces the final energy while keeping the number of orbitals low compared to FCI calculations on an atomic orbital basis set.

References in corpus (14)

• A variational eigenvalue solver on a quantum processor
• The density-matrix renormalization group in the age of matrix product states
• Ab Initio Quantum Chemistry using the Density Matrix Renormalization Group
• Heat-bath Configuration Interaction: An efficient selected CI algorithm inspired by heat-bath sampling
• Tensor network states and algorithms in the presence of a global U(1) symmetry
• Multiresolution analysis of electronic structure: semicardinal and wavelet bases
• Finite automata for caching in matrix product algorithms
• The Elephant in the Room of Density Functional Theory Calculations
• MADNESS: A Multiresolution, Adaptive Numerical Environment for Scientific Simulation
• The MRChem multiresolution analysis code for molecular electronic structure calculations: performance and scaling properties
• Direct determination of optimal real-space orbitals for correlated electronic structure of molecules
• State Diagrams to determine Tree Tensor Network Operators
• VAMPyR -- A High-Level Python Library for Mathematical Operations in a Multiwavelets Representation
• Heat semigroup representation of Laplacian