Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Density matrix renormalization group calculations and entanglement optimization in quantum chemistry

Andrej Antalík

Published Sep 14, 2015• 0 citations
OpenAlex

Abstract

Density matrix renormalization group is a powerful numerical approach originating in solid state physics. Since its introduction to quantum chemis- try it has been successfully applied to many challenging problems and proved especially efficient for systems with strong multi-reference character. As the method is not invariant to the change of orbital ordering, we implemented the automatic orbital ordering optimization procedure built around the entan- glement analysis and employing the Fiedler method. A technique for active space selection based on one-orbital entropy is introduced. We followed the study by Berardo et al. in which EOM-CC methods were applied to excited states of TiO2 clusters. Four lowest-lying singlet excited states of monomer, dimer and trimer were calculated by means of the DMRG method. In case of the monomer there was good agreement for two states, but the calculations suffered from insufficient treatment of dynamic correlation. For the dimer our results were in very good agreement with EOM-CCSDT results and for the trimer further study would be appropriate. 1

Related papers

On the Ordering of Sites in the Density Matrix Renormalization Group using Quantum Mutual Information

Mazen Ali

The density matrix renormalization group (DMRG) of White 1992 remains to this day an integral component of many state-of-the-art methods for efficiently simulating strongly correlated quantum systems. In quantum chemistry, QC-DMRG became a powerful tool for ab initio calculations with the non-relativistic Schrödinger equation. An important issue in QC-DMRG is the so-called ordering problem -- the optimal ordering of DMRG sites corresponding to electronic orbitals that produces the most accurate results. To this end, a commonly used heuristic is the grouping of strongly correlated orbitals as measured via quantum mutual information. In this work, we show how such heuristics can be directly related to minimizing the entanglement entropy of matrix product states and, consequently, to the truncation error of a fixed bond dimension approximation. Key to establishing this link is the strong subadditivity of entropy. This provides rigorous theoretical justification for the orbital ordering methods and suggests alternate ordering criteria.

arXiv:2103.01111•Mar 1, 2021•1 citation

High-performance ab initio density matrix renormalization group method: Applicability to large-scale multireference problems for metal compounds

Yuki Kurashige, Takeshi Yanai

This article presents an efficient and parallelized implementation of the density matrix renormalization group (DMRG) algorithm for quantum chemistry calculations. The DMRG method as a large-scale multireference electronic structure model is by nature particularly efficient for one-dimensionally correlated systems, while the present development is oriented toward applications for polynuclear transition metal compounds, in which the macroscopic one-dimensional structure of electron correlation is absent. A straightforward extension of the DMRG algorithm is proposed with further improvements and aggressive optimizations to allow its application with large multireference active space, which is often demanded for metal compound calculations. Special efficiency is achieved by making better use of sparsity and symmetry in the operator and wave function representations. By accomplishing computationally intensive DMRG calculations, the authors have found that a large number of renormalized basis states are required to represent high entanglement of the electron correlation for metal compound applications, and it is crucial to adopt auxiliary perturbative correction to the projected density matrix during the DMRG sweep optimization in order to attain proper convergence to the solution. Potential energy curve calculations for the Cr(2) molecule near the known equilibrium precisely predicted the full configuration interaction energies with a correlation space of 24 electrons in 30 orbitals [denoted by (24e,30o)]. The energies are demonstrated to be accurate to 0.6mE(h) (the error from the extrapolated best value) when as many as 10,000 renormalized basis states are employed for the left and right DMRG block representations. The relative energy curves for [Cu(2)O(2)](2+) along the isomerization coordinate were obtained from DMRG and other correlated calculations, for which a fairly large orbital space (32e,62o) is modeled as a full correlation space. The DMRG prediction nearly overlaps with the energy curve from the coupled cluster with singles, doubles, and perturbative triple [CCSD(T)] calculations, while the multireference complete active space self-consistent field (CASSCF) calculations with the small reference configuration (8e,8o) are found to overestimate the biradical character of the electronic state of [Cu(2)O(2)](2+) according to the one-electron density matrix analysis.

The Journal of Chemical Physics•Jun 19, 2009•193 citations

Advances in the density matrix renormalization group method for use in quantum chemistry

Dominika Zgid

Despite the success of modern quantum chemistry in predicting properties of organic molecules, the treatment of inorganic systems, which have many close lying states, remains out of quantitative reach for current methods. To treat non-dynamic correlation, we take advantage of the density matrix renormalization group (DMRG) method that has become very successful in the field of solid state physics. We present a detailed study of the DMRG method, and we pay special attention to the evolution of the understanding behind the mathematical structure of the DMRG wave function. Our primary goal is to develop a density matrix renormalization group self--consistent--field (DMRG-SCF) approach, analogous to the complete active space self--consistent field (CASSCF) method, but dealing with large active spaces that are too demanding for the full configuration interaction (FCI) method. \nAs a first step towards such a DMRG-SCF procedure, we present a spin-adapted DMRG algorithm designed to target spin- and spatial-symmetry states that are hard to obtain while using an unrestricted algorithm. \nOur next step is a modification of the DMRG algorithm to obtain decreasing energy at every step during the sweep. This monotonically convergent DMRG scheme lets us obtain the two-body density matrix as a by--product of the existing procedure without any additional cost in storage. Additionally, the two-body density matrix produced at convergence using this scheme is free from the N-representability problem that is present when the two--body density matrix is produced with the two-site DMRG scheme without additional storage cost. Finally, taking advantage of the modifications developed herein, we present results obtained from our DMRG-SCF method. Lastly, we discuss possible ways of incorporating dynamical correlation into the DMRG scheme, in order to obtain a modern multireference approach.

•May 1, 2008•1 citation

Optimizing the density-matrix renormalization group method using quantum information entropy

Örs Legeza, Jenő Sólyom

In order to optimize the ordering of the lattice sites in the momentum space and quantum chemistry versions of the density-matrix renormalization group (DMRG) method we have studied the separability and entanglement of the target state for the one-dimensional Hubbard model and various molecules. By analyzing the behavior of von Neumann entropy we have found criteria that help to fasten convergence. An initialization procedure has been developed which maximizes the Kullback-Leibler entropy and extends the active space in a dynamical fashion. The dynamically extended active space procedure reduces significantly the effective system size during the first half-sweep and accelerates the speed of convergence of momentum space DMRG and quantum chemistry DMRG to a great extent. The effect of lattice site ordering on the number of block states to be kept during the RG procedure is also investigated.

Physical review. B, Condensed matter•Nov 21, 2003•319 citations

Multi-reference Epstein–Nesbet perturbation theory with density matrix renormalization group reference wavefunction

Yinxuan Song, Yifan Cheng, Yingjin Ma +1 more

Abstract The accurate electronic structure calculation for strongly correlated chemical systems requires an adequate description for both static and dynamic electron correlation, and is a persistent challenge for quantum chemistry. In order to account for static and dynamic electron correlations accurately and efficiently, in this work we propose a new method by integrating the density matrix renormalization group (DMRG) method and multi-reference second-order Epstein–Nesbet perturbation theory (ENPT2) with a selected configuration interaction (SCI) approximation. Compared with previous DMRG-based dynamic correlation methods, the DMRG-ENPT2 method extends the range of applicability, allowing us to efficiently calculate systems with very large active space beyond 30 orbitals. We demonstrate this by performing calculations on H 2 S with an active space of (16e, 15o), hexacene with an active space of (26e, 26o) and 2D H 64 square lattice with an active space of (42e, 42o).

arXiv:1909.11954•Electronic Structure•Feb 4, 2020•25 citations

Ab initio Density Matrix Renormalization Group Theory for Multireference Quantum Chemistry

Takeshi Yanai

A major challenge in molecular electronic structure theory arises from chemical systems where the mean-field picture based upon a single electronic configuration fails to properly describe mechanics of many electrons. Such complex electronic states should be accounted for using a correlated superposition of multiple determinants and are referred to as multireference in quantum chemistry. Application domain of traditional multireference methods is severely limited to systems containing a small number of correlated electrons because of their exponential complexity. In recent years, much attention has been drawn to the density matrix renormalization group (DMRG) theory as a promising alternative. It has been introduced in ab initio quantum chemistry calculations and shown to be a highly-scalable multireference method that can handle full quantum degrees of freedom for the strongly-correlated wavefunction with unprecedented large configuration space. In this article, we provide an explicative introduction to the underlying theory and formalism built into the DMRG method and give a brief overview of our methodological extensions that combine DMRG and active-space correlation models.

Molecular Science•Jan 1, 2014•0 citations
Density matrix renormalization group is a powerful numerical approach originating in solid state physics. Since its introduction to quantum chemis- try it has been successfully applied to many challenging problems and proved especially efficient for systems with strong multi-reference character. As the method is not invariant to the change of orbital ordering, we implemented the automatic orbital ordering optimization procedure built around the entan- glement analysis and employing the Fiedler method. A technique for active space selection based on one-orbital entropy is introduced. We followed the study by Berardo et al. in which EOM-CC methods were applied to excited states of TiO2 clusters. Four lowest-lying singlet excited states of monomer, dimer and trimer were calculated by means of the DMRG method. In case of the monomer there was good agreement for two states, but the calculations suffered from insufficient treatment of dynamic correlation. For the dimer our results were in very good agreement with EOM-CCSDT results and for the trimer further study would be appropriate. 1
The density matrix renormalization group (DMRG) of White 1992 remains to this day an integral component of many state-of-the-art methods for efficiently simulating strongly correlated quantum systems. In quantum chemistry, QC-DMRG became a powerful tool for ab initio calculations with the non-relativistic Schrödinger equation. An important issue in QC-DMRG is the so-called ordering problem -- the optimal ordering of DMRG sites corresponding to electronic orbitals that produces the most accurate results. To this end, a commonly used heuristic is the grouping of strongly correlated orbitals as measured via quantum mutual information. In this work, we show how such heuristics can be directly related to minimizing the entanglement entropy of matrix product states and, consequently, to the truncation error of a fixed bond dimension approximation. Key to establishing this link is the strong subadditivity of entropy. This provides rigorous theoretical justification for the orbital ordering methods and suggests alternate ordering criteria.
This article presents an efficient and parallelized implementation of the density matrix renormalization group (DMRG) algorithm for quantum chemistry calculations. The DMRG method as a large-scale multireference electronic structure model is by nature particularly efficient for one-dimensionally correlated systems, while the present development is oriented toward applications for polynuclear transition metal compounds, in which the macroscopic one-dimensional structure of electron correlation is absent. A straightforward extension of the DMRG algorithm is proposed with further improvements and aggressive optimizations to allow its application with large multireference active space, which is often demanded for metal compound calculations. Special efficiency is achieved by making better use of sparsity and symmetry in the operator and wave function representations. By accomplishing computationally intensive DMRG calculations, the authors have found that a large number of renormalized basis states are required to represent high entanglement of the electron correlation for metal compound applications, and it is crucial to adopt auxiliary perturbative correction to the projected density matrix during the DMRG sweep optimization in order to attain proper convergence to the solution. Potential energy curve calculations for the Cr(2) molecule near the known equilibrium precisely predicted the full configuration interaction energies with a correlation space of 24 electrons in 30 orbitals [denoted by (24e,30o)]. The energies are demonstrated to be accurate to 0.6mE(h) (the error from the extrapolated best value) when as many as 10,000 renormalized basis states are employed for the left and right DMRG block representations. The relative energy curves for [Cu(2)O(2)](2+) along the isomerization coordinate were obtained from DMRG and other correlated calculations, for which a fairly large orbital space (32e,62o) is modeled as a full correlation space. The DMRG prediction nearly overlaps with the energy curve from the coupled cluster with singles, doubles, and perturbative triple [CCSD(T)] calculations, while the multireference complete active space self-consistent field (CASSCF) calculations with the small reference configuration (8e,8o) are found to overestimate the biradical character of the electronic state of [Cu(2)O(2)](2+) according to the one-electron density matrix analysis.
Despite the success of modern quantum chemistry in predicting properties of organic molecules, the treatment of inorganic systems, which have many close lying states, remains out of quantitative reach for current methods. To treat non-dynamic correlation, we take advantage of the density matrix renormalization group (DMRG) method that has become very successful in the field of solid state physics. We present a detailed study of the DMRG method, and we pay special attention to the evolution of the understanding behind the mathematical structure of the DMRG wave function. Our primary goal is to develop a density matrix renormalization group self--consistent--field (DMRG-SCF) approach, analogous to the complete active space self--consistent field (CASSCF) method, but dealing with large active spaces that are too demanding for the full configuration interaction (FCI) method. \nAs a first step towards such a DMRG-SCF procedure, we present a spin-adapted DMRG algorithm designed to target spin- and spatial-symmetry states that are hard to obtain while using an unrestricted algorithm. \nOur next step is a modification of the DMRG algorithm to obtain decreasing energy at every step during the sweep. This monotonically convergent DMRG scheme lets us obtain the two-body density matrix as a by--product of the existing procedure without any additional cost in storage. Additionally, the two-body density matrix produced at convergence using this scheme is free from the N-representability problem that is present when the two--body density matrix is produced with the two-site DMRG scheme without additional storage cost. Finally, taking advantage of the modifications developed herein, we present results obtained from our DMRG-SCF method. Lastly, we discuss possible ways of incorporating dynamical correlation into the DMRG scheme, in order to obtain a modern multireference approach.
In order to optimize the ordering of the lattice sites in the momentum space and quantum chemistry versions of the density-matrix renormalization group (DMRG) method we have studied the separability and entanglement of the target state for the one-dimensional Hubbard model and various molecules. By analyzing the behavior of von Neumann entropy we have found criteria that help to fasten convergence. An initialization procedure has been developed which maximizes the Kullback-Leibler entropy and extends the active space in a dynamical fashion. The dynamically extended active space procedure reduces significantly the effective system size during the first half-sweep and accelerates the speed of convergence of momentum space DMRG and quantum chemistry DMRG to a great extent. The effect of lattice site ordering on the number of block states to be kept during the RG procedure is also investigated.
Abstract The accurate electronic structure calculation for strongly correlated chemical systems requires an adequate description for both static and dynamic electron correlation, and is a persistent challenge for quantum chemistry. In order to account for static and dynamic electron correlations accurately and efficiently, in this work we propose a new method by integrating the density matrix renormalization group (DMRG) method and multi-reference second-order Epstein–Nesbet perturbation theory (ENPT2) with a selected configuration interaction (SCI) approximation. Compared with previous DMRG-based dynamic correlation methods, the DMRG-ENPT2 method extends the range of applicability, allowing us to efficiently calculate systems with very large active space beyond 30 orbitals. We demonstrate this by performing calculations on H 2 S with an active space of (16e, 15o), hexacene with an active space of (26e, 26o) and 2D H 64 square lattice with an active space of (42e, 42o).
A major challenge in molecular electronic structure theory arises from chemical systems where the mean-field picture based upon a single electronic configuration fails to properly describe mechanics of many electrons. Such complex electronic states should be accounted for using a correlated superposition of multiple determinants and are referred to as multireference in quantum chemistry. Application domain of traditional multireference methods is severely limited to systems containing a small number of correlated electrons because of their exponential complexity. In recent years, much attention has been drawn to the density matrix renormalization group (DMRG) theory as a promising alternative. It has been introduced in ab initio quantum chemistry calculations and shown to be a highly-scalable multireference method that can handle full quantum degrees of freedom for the strongly-correlated wavefunction with unprecedented large configuration space. In this article, we provide an explicative introduction to the underlying theory and formalism built into the DMRG method and give a brief overview of our methodological extensions that combine DMRG and active-space correlation models.