CheMPS2: a free open-source spin-adapted implementation of the density matrix renormalization group for ab initio quantum chemistry

Abstract

The density matrix renormalization group (DMRG) has become an indispensable numerical tool to find exact eigenstates of finite-size quantum systems with strong correlation. In the fields of condensed matter, nuclear structure and molecular electronic structure, it has significantly extended the system sizes that can be handled compared to full configuration interaction, without losing numerical accuracy. For quantum chemistry (QC), the most efficient implementations of DMRG require the incorporation of particle number, spin and point group symmetries in the underlying matrix product state (MPS) ansatz, as well as the use of so-called complementary operators. The symmetries introduce a sparse block structure in the MPS ansatz and in the intermediary contracted tensors. If a symmetry is non-abelian, the Wigner-Eckart theorem allows to factorize a tensor into a Clebsch-Gordan coefficient and a reduced tensor. In addition, the fermion signs have to be carefully tracked. Because of these challenges, implementing DMRG efficiently for QC is not straightforward. Efficient and freely available implementations are therefore highly desired. In this work we present CheMPS2, our free open-source spin-adapted implementation of DMRG for ab initio QC. Around CheMPS2, we have implemented the augmented Hessian Newton-Raphson complete active space self-consistent field method, with exact Hessian. The bond dissociation curves of the 12 lowest states of the carbon dimer were obtained at the DMRG(28 orbitals, 12 electrons, D${}_{\mathsf{SU}(2)}$=2500)/cc-pVDZ level of theory. The contribution of $1s$ core correlation to the $X^1{\Sigma }_g^{+}$ bond dissociation curve of the carbon dimer was estimated by comparing energies at the DMRG(36o, 12e, D${}_{\mathsf{SU}(2)}$=2500)/cc-pCVDZ and DMRG-SCF(34o, 8e, D${}_{\mathsf{SU}(2)}$=2500)/cc-pCVDZ levels of theory.

16 pages, 13 figures

References in corpus (19)

• The density-matrix renormalization group in the age of matrix product states
• Renormalization algorithms for Quantum-Many Body Systems in two and higher dimensions
• From density-matrix renormalization group to matrix product states
• The radical character of the acenes: A density matrix renormalization group study
• Orbital Optimization in the Density Matrix Renormalization Group, with applications to polyenes and β-carotene
• A spin-adapted Density Matrix Renormalization Group algorithm for quantum chemistry
• Continuous Matrix Product States for Quantum Fields
• Simulating Strongly Correlated Quantum Systems with Tree Tensor Networks
• Multireference Correlation in Long Molecules with the Quadratic Scaling Density Matrix Renormalization Group
• Entanglement Measures for Single- and Multi-Reference Correlation Effects
• Efficient Tree Tensor Network States (TTNS) for Quantum Chemistry: Generalizations of the Density Matrix Renormalization Group Algorithm
• Accurate ab initio spin densities
• Four-Component Density Matrix Renormalization Group
• Complete-Graph Tensor Network States: A New Fermionic Wave Function Ansatz for Molecules
• Targeted Excited State Algorithms
• Construction of CASCI-type wave functions for very large active spaces
• Linear Response Theory for the Density Matrix Renormalization Group: Efficient Algorithms for Strongly Correlated Excited States
• Longitudinal static optical properties of hydrogen chains: finite field extrapolations of matrix product state calculations
• Optimizing Hartree-Fock orbitals by the density-matrix renormalization group

Cited by in corpus (58)

• Recent developments in the PySCF program package
• The density matrix renormalization group for ab initio quantum chemistry
• Tensor product methods and entanglement optimization for ab initio quantum chemistry
• A practical guide to density matrix embedding theory in quantum chemistry
• Reference Energies for Double Excitations
• An Efficient Matrix Product Operator Representation of the Quantum-Chemical Hamiltonian
• Efficient description of strongly correlated electrons with mean-field cost
• Block2: a comprehensive open source framework to develop and apply state-of-the-art DMRG algorithms in electronic structure and beyond
• Density Matrix Renormalization Group with Efficient Dynamical Electron Correlation Through Range Separation
• Projected seniority-two orbital optimization of the Antisymmetric Product of one-reference orbital Geminal
• Spin-adapted Matrix Product States and Operators
• DMRG-SCF study of the singlet, triplet, and quintet states of oxo-Mn(Salen)
• Matrix Product State applications for the ALPS project
• Low communication high performance ab initio density matrix renormalization group algorithms
• Linearized Coupled Cluster Correction on the Antisymmetric Product of 1 reference orbital Geminals
• DMRG-CASPT2 study of the longitudinal static second hyperpolarizability of all-trans polyenes
• Computing vibrational eigenstates with tree tensor network states (TTNS)
• A time-dependent formulation of multi-reference perturbation theory
• T3NS: three-legged tree tensor network states
• Open source Matrix Product States: Opening ways to simulate entangled many-body quantum systems in one dimension
• Multi-reference algebraic diagrammatic construction theory for excited states: General formulation and first-order implementation
• Investigation of metal-insulator like transition through the ab initio density matrix renormalization group approach
• Multi-State Pair-Density Functional Theory
• Quantum Entanglement in Carbon-Carbon, Carbon-Phosphorus and Silicon-Silicon Bonds
• Stochastic multi-reference perturbation theory with application to linearized coupled cluster method
• Self-Consistent Embedding of Density-Matrix Renormalization Group Wavefunctions in a Density Functional Environment
• A tensor network view of multilayer multiconfiguration time-dependent Hartree methods
• Analytical Gradient Theory for Strongly Contracted (SC-) and Partially Contracted (PC-) N-Electron Valence State Perturbation Theory (NEVPT2)
• Second-Order Multi-Reference Algebraic Diagrammatic Construction Theory for Photoelectron Spectra of Strongly Correlated Systems
• Ultrafast ab-initio Quantum Chemistry Using Matrix Product States
• A call to arms: making the case for more reusable libraries
• Combined unitary and symmetric group approach applied to low-dimensional spin systems
• Exact NMR simulation of protein-size spin systems using tensor train formalism
• Multireference Algebraic Diagrammatic Construction Theory for Excited States: Extended Second-Order Implementation and Benchmark
• Projector quantum Monte Carlo with matrix product states
• A beginner's guide to non-abelian iPEPS for correlated fermions
• The Python-based Simulations of Chemistry Framework (PySCF)
• An efficient method for strongly correlated electrons in two-dimensions
• Analysis of Electron Correlation Effects in Strongly Correlated Systems (${\mathrm{N}}_2$ and ${\mathrm{N}}_2^{+}$) by applying DMRG and Quantum Information Theory
• Density-Matrix Renormalization Group Algorithm with Multi-Level Active Space
• The Effect of Geometry, Spin and Orbital Optimization in Achieving Accurate, Fully-Correlated Results for Iron-Sulfur Cubanes
• Challenges for variational reduced-density-matrix theory with three-particle N-representability conditions
• A comparison between the one- and two-step spin-orbit coupling approaches based on the ab initio Density Matrix Renormalization Group
• The Three-Legged Tree Tensor Networks with SU(2)- and molecular point group symmetry
• Towards large-scale restricted active space calculations inspired by the Schmidt decomposition
• Unveiling Intrinsic Many-Body Complexity by Compressing Single-Body Triviality
• Data-driven Refinement of Electronic Energies from Two-Electron Reduced-Density-Matrix Theory
• The seniority quantum number in Tensor Network States
• Computing vibrational energy levels by solving linear equations using a tensor method with an imposed rank
• 2500 vibronic eigenstates of the NO${}_3$ radical
• Laplace-transformed multi-reference second-order perturbation theories in the atomic and active molecular orbital basis
• Inversion symmetry of singular values and a new orbital ordering method in tensor train approximations for quantum chemistry
• Quantification of electron correlation effects - Quantum Information Theory versus Method of Increments
• 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
• Rapidly convergent quantum Monte Carlo using a Chebyshev projector
• Effective Hamiltonians from Spin Adapted Configuration Interaction
• The Software Landscape for the Density Matrix Renormalization Group