Optimization of highly excited matrix product states with an application to vibrational spectroscopy

Abstract

Configuration-interaction-type calculations on electronic and vibrational structure are often the method of choice for the reliable approximation of many-particle wave functions and energies. The exponential scaling, however, limits their application range. An efficient approximation to the full configuration interaction solution can be obtained with the density matrix renormalization group (DMRG) algorithm without a restriction to a predefined excitation level. In a standard DMRG implementation, however, excited states are calculated with a ground-state optimization in the space orthogonal to all lower lying wave function solutions. A trivial parallelization is therefore not possible and the calculation of highly excited states becomes prohibitively expensive, especially in regions with a high density of states. Here, we introduce two variants of the density matrix renormalization group algorithm that allow us to target directly specific energy regions and therefore highly excited states. The first one, based on shift-and-invert techniques, is particularly efficient for low-lying states, but is not stable in regions with a high density of states. The second one, based on the folded auxiliary operator, is less efficient, but more accurate in targeting high-energy states. We apply the algorithm to the solution of the nuclear Schroedinger equation, but emphasize that it can be applied to the diagonalization of general Hamiltonians as well, such as the electronic Coulomb Hamiltonian to address X-ray spectra. In combination with several root-homing algorithms and a stochastic sampling of the determinant space, excited states of interest can be adequately tracked and analyzed during the optimization. We demonstrate that we can accurately calculate prominent spectral features of large molecules such as the sarcosyn-glycine dipeptide.

arXiv comment26 pages, 7 figures, 3 tables

References in corpus (23)

• The density-matrix renormalization group in the age of matrix product states
• The density-matrix renormalization group
• From density-matrix renormalization group to matrix product states
• A class of ansatz wave functions for 1D spin systems and their relation to DMRG
• The density matrix renormalization group for ab initio quantum chemistry
• Tensor product methods and entanglement optimization for ab initio quantum chemistry
• An Efficient Matrix Product Operator Representation of the Quantum-Chemical Hamiltonian
• A Strictly Single-Site DMRG Algorithm with Subspace Expansion
• Obtaining highly-excited eigenstates of many-body localized Hamiltonians by the density matrix renormalization group
• Perfect Sampling with Unitary Tensor Networks
• Finding matrix product state representations of highly-excited eigenstates of many-body localized Hamiltonians
• New Approaches for ab initio Calculations of Molecules with Strong Electron Correlation
• Time-step targeting time-dependent and dynamical density matrix renormalization group algorithms with ab initio Hamiltonians
• Targeted Excited State Algorithms
• $\sigma$-SCF: A Direct Energy-targeting Method To Mean-field Excited States
• Vibrational Density Matrix Renormalization Group
• Construction of CASCI-type wave functions for very large active spaces
• Error estimates for extrapolations with matrix-product states
• Calculating vibrational spectra of molecules using tensor train decomposition
• Externally-Contracted Multi-Reference Configuration Interaction Method Using a DMRG Reference Wave Function
• Efficient reconstruction of CASCI-type wave functions for a DMRG state using quantum information theory and genetic algorithm
• Obtaining highly excited eigenstates of the localized XX chain via DMRG-X
• Comparison of different eigensolvers for calculating vibrational spectra using low-rank, sum-of-product basis functions

Cited by in corpus (26)

• The Density Matrix Renormalization Group in Chemistry and Molecular Physics: Recent Developments and New Challenges
• Large-scale quantum-dynamics with matrix product states
• Hardware Efficient Quantum Algorithms for Vibrational Structure Calculations
• Computing vibrational eigenstates with tree tensor network states (TTNS)
• A General Automatic Method for Optimal Construction of Matrix Product Operators Using Bipartite Graph Theory
• Simulating X-ray absorption spectra with CASSCF linear response methods
• A tensor network view of multilayer multiconfiguration time-dependent Hartree methods
• Excited-state DMRG made simple with FEAST
• Transcorrelated Density Matrix Renormalization Group
• Nuclear-Electronic All-Particle Density Matrix Renormalization Group
• Electron Dynamics with the Time-Dependent Density Matrix Renormalization Group
• Quantum Proton Effects from Density Matrix Renormalization Group Calculations
• Flexible DMRG-based framework for anharmonic vibrational calculations
• SCINE -- Software for Chemical Interaction Networks
• Minimal matrix product states and generalizations of mean-field and geminal wavefunctions
• Multicomponent Heat-Bath Configuration Interaction with the Perturbative Correction for the Calculation of Protonic Excited States
• Computing vibrational energy levels by solving linear equations using a tensor method with an imposed rank
• Ultrafast Charge Migration Dynamics in Enol Keto Tautomerization Monitored with a Local Soft-X-Ray Probe
• Vibrational Entanglement through the Lens of Quantum Information Measures
• QCMaquis 4.0: Multi-Purpose Electronic, Vibrational, and Vibronic Structure and Dynamics Calculations with the Density Matrix Renormalization Group
• Computing excited eigenstates using inexact Lanczos methods and tree tensor network states
• Trotter simulation of vibrational Hamiltonians on a quantum computer
• The Software Landscape for the Density Matrix Renormalization Group
• A Synthesis of Hidden Subgroup Quantum Algorithms and Quantum Chemical Dynamics
• Modal Backflow Neural Quantum States for Anharmonic Vibrational Calculations
• Accurate, full-dimensional computations of thousands of complex vibrational eigenstates with tree tensor network states