A General Automatic Method for Optimal Construction of Matrix Product Operators Using Bipartite Graph Theory

Abstract

Constructing matrix product operators (MPO) is at the core of the modern density matrix renormalization group (DMRG) and its time dependent formulation. For DMRG to be conveniently used in different problems described by different Hamiltonians, in this work we propose a new generic algorithm to construct the MPO of an arbitrary operator with a sum-of-products form based on the bipartite graph theory. We show that the method has the following advantages: (i) It is automatic in that only the definition of the operator is required; (ii) It is symbolic thus free of any numerical error; (iii) The complementary operator technique can be fully employed so that the resulting MPO is globally optimal for any given order of degrees of freedom; (iv) The symmetry of the system could be fully employed to reduce the dimension of MPO. To demonstrate the effectiveness of the new algorithm, the MPOs of Hamiltonians ranging from the prototypical spin-boson model and Holstein model to the more complicated ab initio electronic Hamiltonian and the anharmonic vibrational Hamiltonian with sextic force field are constructed. It is found that for the former three cases, our automatic algorithm can reproduce exactly the same MPOs as the optimally hand-crafted ones already known in the literature.

31 pages, 6 figures

References in corpus (31)

• The density-matrix renormalization group in the age of matrix product states
• Thermodynamic limit of the density matrix renormalization for the spin-1 Heisenberg chain
• Unifying time evolution and optimization with matrix product states
• The Density Matrix Renormalization Group in Chemistry and Molecular Physics: Recent Developments and New Challenges
• An Efficient Matrix Product Operator Representation of the Quantum-Chemical Hamiltonian
• Efficient Tree Tensor Network States (TTNS) for Quantum Chemistry: Generalizations of the Density Matrix Renormalization Group Algorithm
• Density-matrix renormalization-group method in momentum space
• Generic Construction of Efficient Matrix Product Operators
• Finite automata for caching in matrix product algorithms
• Time-Dependent Density Matrix Renormalization Group Algorithms for Nearly Exact Absorption and Fluorescence Spectra of Molecular Aggregates at Both Zero and Finite Temperature
• Finding matrix product state representations of highly-excited eigenstates of many-body localized Hamiltonians
• Large-scale quantum-dynamics with matrix product states
• Time-dependent Density Matrix Renormalization Group Quantum Dynamics for Realistic Chemical Systems
• Second-Order Self-Consistent-Field Density-Matrix Renormalization Group
• Multi-Layer Potfit: An Accurate Potential Representation for Efficient High-Dimensional Quantum Dynamics
• Numerical Assessment for Accuracy and GPU Acceleration of TD-DMRG Time Evolution Schemes
• Time-step targeting time-dependent and dynamical density matrix renormalization group algorithms with ab initio Hamiltonians
• Computing vibrational eigenstates with tree tensor network states (TTNS)
• Targeted Excited State Algorithms
• Vibrational Density Matrix Renormalization Group
• T3NS: three-legged tree tensor network states
• Calculating vibrational spectra of molecules using tensor train decomposition
• Linear Response Theory for the Density Matrix Renormalization Group: Efficient Algorithms for Strongly Correlated Excited States
• Finite Temperature TD-DMRG for the Carrier Mobility of Organic Semiconductors
• The Thouless theorem for matrix product states and subsequent post-density matrix renormalization group methods
• Optimization of highly excited matrix product states with an application to vibrational spectroscopy
• Matrix product state formulation of the multiconfiguration time-dependent Hartree theory
• Externally-Contracted Multi-Reference Configuration Interaction Method Using a DMRG Reference Wave Function
• Ultrafast ab-initio Quantum Chemistry Using Matrix Product States
• Towards DMRG-tailored coupled cluster method in the 4c-relativistic domain
• On the generalization of the exponential basis for tensor network representations of long-range interactions in two and three dimensions

Cited by in corpus (19)

• Block2: a comprehensive open source framework to develop and apply state-of-the-art DMRG algorithms in electronic structure and beyond
• A General Charge Transport Picture for Organic Semiconductors with Nonlocal Electron-Phonon Couplings
• Low communication high performance ab initio density matrix renormalization group algorithms
• Learning tensor networks with tensor cross interpolation: new algorithms and libraries
• Non-Markovian Stochastic Schrödinger Equation: Matrix Product State Approach to the Hierarchy of Pure States
• Transcorrelated Density Matrix Renormalization Group
• Optimal Tree Tensor Network Operators for Tensor Network Simulations: Applications to Open Quantum Systems
• A comparison between the one- and two-step spin-orbit coupling approaches based on the ab initio Density Matrix Renormalization Group
• Unbiasing Fermionic Auxiliary-Field Quantum Monte Carlo with Matrix Product State Trial Wavefunctions
• Optimal Symbolic Construction of Matrix Product Operators and Tree Tensor Network Operators
• State Diagrams to determine Tree Tensor Network Operators
• QCMaquis 4.0: Multi-Purpose Electronic, Vibrational, and Vibronic Structure and Dynamics Calculations with the Density Matrix Renormalization Group
• Simulating Quantum Turbulence with Matrix Product States
• Extended Non-Markovian Stochastic Schrödinger Equation with Complex Frequency Modes for General Basis Functions
• Scaling up the transcorrelated density matrix renormalization group
• The Software Landscape for the Density Matrix Renormalization Group
• Enhanced Krylov Methods for Molecular Hamiltonians: Reduced Memory Cost and Complexity Scaling via Tensor Hypercontraction
• Multipartite Entanglement in Crossing the Quantum Critical Point
• paces: Parallelized Application of Co-Evolving Subspaces, a method for computing quantum dynamics on GPUs