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

arXiv:1605.00940 · Physical review. B./Physical review. B

Efficient perturbation theory to improve the density matrix renormalization group

Emanuele Tirrito, Shi-Ju Ran, Andrew J. Ferris, Ian P. McCulloch, Maciej Lewenstein

Published Feb 21, 2017• 4 citations
arXivHTMLPDF

Abstract

The density matrix renormalization group (DMRG) is one of the most powerful numerical methods available for many-body systems. It has been applied to solve many physical problems, including the calculation of ground states and dynamical properties. In this work, we develop a perturbation theory of the DMRG (PT-DMRG) to greatly increase its accuracy in an extremely simple and efficient way. Using the canonical matrix product state (MPS) representation for the ground state of the considered system, a set of orthogonal basis functions $\left\{|{\ensuremath{\psi}}_{i}\ensuremath{\rangle}\right\}$ is introduced to describe the perturbations to the ground state obtained by the conventional DMRG. The Schmidt numbers of the MPS that are beyond the bond dimension cutoff are used to define these perturbation terms. The perturbed Hamiltonian is then defined as ${\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{H}}_{ij}=\ensuremath{\langle}{\ensuremath{\psi}}_{i}|\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}|{\ensuremath{\psi}}_{j}\ensuremath{\rangle}$; its ground state permits us to calculate physical observables with a considerably improved accuracy compared to the original DMRG results. We benchmark the second-order perturbation theory with the help of a one-dimensional Ising chain in a transverse field and the Heisenberg chain, where the precision of the DMRG is shown to be improved $O(10)$ times. Furthermore, for moderate $L$ the errors of the DMRG and PT-DMRG both scale linearly with ${L}^{\ensuremath{-}1}$ (with $L$ being the length of the chain). The linear relation between the dimension cutoff of the DMRG and that of the PT-DMRG at the same precision shows a considerable improvement in efficiency, especially for large dimension cutoffs. In the thermodynamic limit we show that the errors of the PT-DMRG scale with $\sqrt{{L}^{\ensuremath{-}1}}$. Our work suggests an effective way to define the tangent space of the ground-state MPS, which may shed light on the properties beyond the ground state. This second-order PT-DMRG can be readily generalized to higher orders, as well as applied to models in higher dimensions.

DOI
10.1103/physrevb.95.064110

Related papers

A Perturbative Density Matrix Renormalization Group Algorithm for Large Active Spaces

Sheng Guo, Zhendong Li, Garnet Kin‐Lic Chan

We describe a low cost alternative to the standard variational DMRG (density matrix renormalization group) algorithm that is analogous to the combination of the selected configuration interaction plus perturbation theory (SCI+PT). We denote the resulting method p-DMRG (perturbative DMRG) to distinguish it from the standard variational DMRG. p-DMRG is expected to be useful for systems with very large active spaces, for which variational DMRG becomes too expensive. Similar to SCI+PT, in p-DMRG, a zeroth-order wave function is first obtained by a standard DMRG calculation but with a small bond dimension. Then, the residual correlation is recovered by a second-order perturbative treatment. We discuss the choice of partitioning for perturbation theory, which is crucial for its accuracy and robustness. To circumvent the problem of a large bond dimension in the first-order wave function, we use a sum of matrix product states to expand the first-order wave function, yielding substantial savings in computational cost and memory. We also propose extrapolation schemes to reduce the errors in the zeroth- and first-order wave functions. Numerical results for Cr 2 with a (28e, 76o) active space and 1,3-butadiene with an (22e, 82o) active space reveal that p-DMRG provides ground state energies of a similar quality to variational DMRG with very large bond dimensions but at a significantly lower computational cost. This suggests that p-DMRG will be an efficient tool for benchmark studies in the future.

arXiv:1803.07150•Journal of Chemical Theory and Computation•Jun 21, 2018•39 citations

Scaling up the transcorrelated density matrix renormalization group

Benjamin C. Corbett, Akimasa Miyake

Explicitly correlated methods, such as the transcorrelated method which shifts a Jastrow or Gutzwiller correlator from the wave function to the Hamiltonian, are designed for high-accuracy calculations of electronic structures, but their application to larger systems has been hampered by the computational cost. We develop improved techniques for the transcorrelated density-matrix renormalization group (DMRG), in which the ground state of the transcorrelated Hamiltonian is represented as a matrix product state (MPS), and demonstrate large-scale calculations of the ground-state energy of the two-dimensional Fermi-Hubbard model. Our developments stem from three technical inventions: (i) constructing matrix product operators (MPOs) of transcorrelated Hamiltonians with low bond dimension and high sparsity, (ii) exploiting the entanglement structure of the ground states to increase the accuracy of the MPS representation, and (iii) optimizing the nonlinear parameter of the Gutzwiller correlator to mitigate the nonvariational nature of the transcorrelated method. We examine systems of size up to $12\ifmmode\times\else\texttimes\fi{}12$ lattice sites, four times larger than previous transcorrelated DMRG studies, and demonstrate that transcorrelated DMRG yields significant improvements over standard nontranscorrelated DMRG for equivalent computational effort. Transcorrelated DMRG reduces the error of the ground-state energy by $2.4\ifmmode\times\else\texttimes\fi{}$--$14\ifmmode\times\else\texttimes\fi{}$, with the smallest improvement seen for a small system at half filling and the largest improvement in a dilute closed-shell system.

arXiv:2506.07441•Physical review. B./Physical review. B•Sep 15, 2025•1 citation

Density Matrix Renormalization Group for Transcorrelated Hamiltonians: Ground and Excited States in Molecules

Ke Liao, Huanchen Zhai, Evelin Martine Corvid Christlmaier +4 more

We present the theory of a density matrix renormalization group (DMRG) algorithm which can solve for both the ground and excited states of non-Hermitian transcorrelated Hamiltonians and show applications in molecular systems. Transcorrelation (TC) accelerates the basis set convergence rate by including known physics (such as, but not limited to, the electron-electron cusp) in the Jastrow factor used for the similarity transformation. It also improves the accuracy of approximate methods such as coupled cluster singles and doubles (CCSD) as shown by recent studies. However, the non-Hermiticity of the TC Hamiltonians poses challenges for variational methods like DMRG. Imaginary-time evolution on the matrix product state (MPS) in the DMRG framework has been proposed to circumvent this problem, but this is currently limited to treating the ground state and has lower efficiency than the time-independent DMRG (TI-DMRG) due to the need to eliminate Trotter errors. In this work, we show that with minimal changes to the existing TI-DMRG algorithm, namely, replacing the original Davidson solver with the general Davidson solver to solve the non-Hermitian effective Hamiltonians at each site for a few low-lying right eigenstates, and following the rest of the original DMRG recipe, one can find the ground and excited states with improved efficiency compared to the original DMRG when extrapolating to the infinite bond dimension limit in the same basis set. An accelerated basis set convergence rate is also observed, as expected, within the TC framework.

Journal of Chemical Theory and Computation•Mar 13, 2023•33 citations

Inner Space Perturbation Theory in Matrix Product States: Replacing Expensive Iterative Diagonalization

Jiajun Ren, Yuanping Yi, Zhigang Shuai

We propose an inner space perturbation theory (isPT) to replace the expensive iterative diagonalization in the standard density matrix renormalization group theory (DMRG). The retained reduced density matrix eigenstates are partitioned into the active and secondary space. The first-order wave function and the second- and third-order energies are easily computed by using one step Davidson iteration. Our formulation has several advantages including (i) keeping a balance between the efficiency and accuracy, (ii) capturing more entanglement with the same amount of computational time, (iii) recovery of the standard DMRG when all the basis states belong to the active space. Numerical examples for the polyacenes and periacene show that the efficiency gain is considerable and the accuracy loss due to the perturbation treatment is very small, when half of the total basis states belong to the active space. Moreover, the perturbation calculations converge in all our numerical examples.

Journal of Chemical Theory and Computation•Sep 6, 2016•21 citations

Density Matrix Renormalization Group for Transcorrelated Hamiltonians: Ground and Excited States in \emph{ab initio} Systems

Ke Liao, Zhai, Huanchen, Evelin Martine Corvid Christlmaier +4 more

We present the theory of a density matrix renormalization group (DMRG) algorithm which can solve for both the ground and excited states of non-Hermitian transcorrelated Hamiltonians, and show applications in \emph{ab initio} molecular systems. Transcorrelation (TC) accelerates the basis set convergence rate by including known physics (such as, but not limited to, the electron-electron cusp) in the Jastrow factor used for the similarity transformation. It also improves the accuracy of approximate methods such as coupled cluster singles and doubles (CCSD) as shown by some recent studies. However, the non-Hermiticity of the TC Hamiltonians poses challenges for variational methods like DMRG. Imaginary-time evolution on the matrix product state (MPS) in the DMRG framework has been proposed to circumvent this problem; but this is currently limited to treating the ground state, and has lower efficiency than the time-independent DMRG (TI-DMRG), due to the need to eliminate Trotter errors. In this work, we show that with minimal changes to the existing TI-DMRG algorithm, namely replacing the original Davidson solver with the general Davidson solver to solve the non-Hermitian effective Hamiltonians at each site for a few low-lying right eigenstates, and following the rest of the original DMRG recipe, one can find the ground and excited states with improved efficiency compared to the original DMRG when extrapolating to the infinite bond dimension limit in the same basis set. Accelerated basis set convergence rate is also observed, as expected, within the TC framework.

arXiv:2211.00173•Oct 31, 2022•1 citation

An efficient density matrix renormalization group algorithm for chains with periodic boundary condition

Dayasindhu Dey, Debasmita Maiti, Manoranjan Kumar

The Density Matrix Renormalization Group (DMRG) is a state-of-the-art numerical technique for a one dimensional quantum many-body system; but calculating accurate results for a system with Periodic Boundary Condition (PBC) from the conventional DMRG has been a challenging job from the inception of DMRG. The recent development of the Matrix Product State (MPS) algorithm gives a new approach to find accurate results for the one dimensional PBC system. The most efficient implementation of the MPS algorithm can scale as O($p \times m^3$), where $p$ can vary from 4 to $m^2$. In this paper, we propose a new DMRG algorithm, which is very similar to the conventional DMRG and gives comparable accuracy to that of MPS. The computation effort of the new algorithm goes as O($m^3$) and the conventional DMRG code can be easily modified for the new algorithm.Received: 2 August 2016, Accepted: 12 October 2016; Edited by: K. Hallberg; DOI: http://dx.doi.org/10.4279/PIP.080006Cite as: D Dey, D Maiti, M Kumar, Papers in Physics 8, 080006 (2016)This paper, by D Dey, D Maiti, M Kumar, is licensed under the Creative Commons Attribution License 3.0.

arXiv:1605.09301•Papers in Physics•Nov 18, 2016•18 citations
The density matrix renormalization group (DMRG) is one of the most powerful numerical methods available for many-body systems. It has been applied to solve many physical problems, including the calculation of ground states and dynamical properties. In this work, we develop a perturbation theory of the DMRG (PT-DMRG) to greatly increase its accuracy in an extremely simple and efficient way. Using the canonical matrix product state (MPS) representation for the ground state of the considered system, a set of orthogonal basis functions {∣\ensuremathψi\ensuremath⟩}\left\{|{\ensuremath{\psi}}_{i}\ensuremath{\rangle}\right\}{∣\ensuremathψi​\ensuremath⟩} is introduced to describe the perturbations to the ground state obtained by the conventional DMRG. The Schmidt numbers of the MPS that are beyond the bond dimension cutoff are used to define these perturbation terms. The perturbed Hamiltonian is then defined as H\ifmmode~\else˜\fiij=\ensuremath⟨\ensuremathψi∣H\ifmmode^\elseˆ\fi∣\ensuremathψj\ensuremath⟩{\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{H}}_{ij}=\ensuremath{\langle}{\ensuremath{\psi}}_{i}|\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}|{\ensuremath{\psi}}_{j}\ensuremath{\rangle}H\ifmmode~\else˜\fi​ij​=\ensuremath⟨\ensuremathψi​∣H\ifmmode^\elseˆ\fi​∣\ensuremathψj​\ensuremath⟩; its ground state permits us to calculate physical observables with a considerably improved accuracy compared to the original DMRG results. We benchmark the second-order perturbation theory with the help of a one-dimensional Ising chain in a transverse field and the Heisenberg chain, where the precision of the DMRG is shown to be improved O(10)O(10)O(10) times. Furthermore, for moderate LLL the errors of the DMRG and PT-DMRG both scale linearly with L\ensuremath−1{L}^{\ensuremath{-}1}L\ensuremath−1 (with LLL being the length of the chain). The linear relation between the dimension cutoff of the DMRG and that of the PT-DMRG at the same precision shows a considerable improvement in efficiency, especially for large dimension cutoffs. In the thermodynamic limit we show that the errors of the PT-DMRG scale with L\ensuremath−1\sqrt{{L}^{\ensuremath{-}1}}L\ensuremath−1​. Our work suggests an effective way to define the tangent space of the ground-state MPS, which may shed light on the properties beyond the ground state. This second-order PT-DMRG can be readily generalized to higher orders, as well as applied to models in higher dimensions.
We describe a low cost alternative to the standard variational DMRG (density matrix renormalization group) algorithm that is analogous to the combination of the selected configuration interaction plus perturbation theory (SCI+PT). We denote the resulting method p-DMRG (perturbative DMRG) to distinguish it from the standard variational DMRG. p-DMRG is expected to be useful for systems with very large active spaces, for which variational DMRG becomes too expensive. Similar to SCI+PT, in p-DMRG, a zeroth-order wave function is first obtained by a standard DMRG calculation but with a small bond dimension. Then, the residual correlation is recovered by a second-order perturbative treatment. We discuss the choice of partitioning for perturbation theory, which is crucial for its accuracy and robustness. To circumvent the problem of a large bond dimension in the first-order wave function, we use a sum of matrix product states to expand the first-order wave function, yielding substantial savings in computational cost and memory. We also propose extrapolation schemes to reduce the errors in the zeroth- and first-order wave functions. Numerical results for Cr 2 with a (28e, 76o) active space and 1,3-butadiene with an (22e, 82o) active space reveal that p-DMRG provides ground state energies of a similar quality to variational DMRG with very large bond dimensions but at a significantly lower computational cost. This suggests that p-DMRG will be an efficient tool for benchmark studies in the future.
Explicitly correlated methods, such as the transcorrelated method which shifts a Jastrow or Gutzwiller correlator from the wave function to the Hamiltonian, are designed for high-accuracy calculations of electronic structures, but their application to larger systems has been hampered by the computational cost. We develop improved techniques for the transcorrelated density-matrix renormalization group (DMRG), in which the ground state of the transcorrelated Hamiltonian is represented as a matrix product state (MPS), and demonstrate large-scale calculations of the ground-state energy of the two-dimensional Fermi-Hubbard model. Our developments stem from three technical inventions: (i) constructing matrix product operators (MPOs) of transcorrelated Hamiltonians with low bond dimension and high sparsity, (ii) exploiting the entanglement structure of the ground states to increase the accuracy of the MPS representation, and (iii) optimizing the nonlinear parameter of the Gutzwiller correlator to mitigate the nonvariational nature of the transcorrelated method. We examine systems of size up to $12\ifmmode\times\else\texttimes\fi{}12$ lattice sites, four times larger than previous transcorrelated DMRG studies, and demonstrate that transcorrelated DMRG yields significant improvements over standard nontranscorrelated DMRG for equivalent computational effort. Transcorrelated DMRG reduces the error of the ground-state energy by $2.4\ifmmode\times\else\texttimes\fi{}$--$14\ifmmode\times\else\texttimes\fi{}$, with the smallest improvement seen for a small system at half filling and the largest improvement in a dilute closed-shell system.
We present the theory of a density matrix renormalization group (DMRG) algorithm which can solve for both the ground and excited states of non-Hermitian transcorrelated Hamiltonians and show applications in molecular systems. Transcorrelation (TC) accelerates the basis set convergence rate by including known physics (such as, but not limited to, the electron-electron cusp) in the Jastrow factor used for the similarity transformation. It also improves the accuracy of approximate methods such as coupled cluster singles and doubles (CCSD) as shown by recent studies. However, the non-Hermiticity of the TC Hamiltonians poses challenges for variational methods like DMRG. Imaginary-time evolution on the matrix product state (MPS) in the DMRG framework has been proposed to circumvent this problem, but this is currently limited to treating the ground state and has lower efficiency than the time-independent DMRG (TI-DMRG) due to the need to eliminate Trotter errors. In this work, we show that with minimal changes to the existing TI-DMRG algorithm, namely, replacing the original Davidson solver with the general Davidson solver to solve the non-Hermitian effective Hamiltonians at each site for a few low-lying right eigenstates, and following the rest of the original DMRG recipe, one can find the ground and excited states with improved efficiency compared to the original DMRG when extrapolating to the infinite bond dimension limit in the same basis set. An accelerated basis set convergence rate is also observed, as expected, within the TC framework.
We propose an inner space perturbation theory (isPT) to replace the expensive iterative diagonalization in the standard density matrix renormalization group theory (DMRG). The retained reduced density matrix eigenstates are partitioned into the active and secondary space. The first-order wave function and the second- and third-order energies are easily computed by using one step Davidson iteration. Our formulation has several advantages including (i) keeping a balance between the efficiency and accuracy, (ii) capturing more entanglement with the same amount of computational time, (iii) recovery of the standard DMRG when all the basis states belong to the active space. Numerical examples for the polyacenes and periacene show that the efficiency gain is considerable and the accuracy loss due to the perturbation treatment is very small, when half of the total basis states belong to the active space. Moreover, the perturbation calculations converge in all our numerical examples.
We present the theory of a density matrix renormalization group (DMRG) algorithm which can solve for both the ground and excited states of non-Hermitian transcorrelated Hamiltonians, and show applications in \emph{ab initio} molecular systems. Transcorrelation (TC) accelerates the basis set convergence rate by including known physics (such as, but not limited to, the electron-electron cusp) in the Jastrow factor used for the similarity transformation. It also improves the accuracy of approximate methods such as coupled cluster singles and doubles (CCSD) as shown by some recent studies. However, the non-Hermiticity of the TC Hamiltonians poses challenges for variational methods like DMRG. Imaginary-time evolution on the matrix product state (MPS) in the DMRG framework has been proposed to circumvent this problem; but this is currently limited to treating the ground state, and has lower efficiency than the time-independent DMRG (TI-DMRG), due to the need to eliminate Trotter errors. In this work, we show that with minimal changes to the existing TI-DMRG algorithm, namely replacing the original Davidson solver with the general Davidson solver to solve the non-Hermitian effective Hamiltonians at each site for a few low-lying right eigenstates, and following the rest of the original DMRG recipe, one can find the ground and excited states with improved efficiency compared to the original DMRG when extrapolating to the infinite bond dimension limit in the same basis set. Accelerated basis set convergence rate is also observed, as expected, within the TC framework.
The Density Matrix Renormalization Group (DMRG) is a state-of-the-art numerical technique for a one dimensional quantum many-body system; but calculating accurate results for a system with Periodic Boundary Condition (PBC) from the conventional DMRG has been a challenging job from the inception of DMRG. The recent development of the Matrix Product State (MPS) algorithm gives a new approach to find accurate results for the one dimensional PBC system. The most efficient implementation of the MPS algorithm can scale as O(p×m3p \times m^3p×m3), where ppp can vary from 4 to m2m^2m2. In this paper, we propose a new DMRG algorithm, which is very similar to the conventional DMRG and gives comparable accuracy to that of MPS. The computation effort of the new algorithm goes as O(m3m^3m3) and the conventional DMRG code can be easily modified for the new algorithm.Received: 2 August 2016, Accepted: 12 October 2016; Edited by: K. Hallberg; DOI: http://dx.doi.org/10.4279/PIP.080006Cite as: D Dey, D Maiti, M Kumar, Papers in Physics 8, 080006 (2016)This paper, by D Dey, D Maiti, M Kumar, is licensed under the Creative Commons Attribution License 3.0.