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

Notes on density matrix perturbation theory

Truflandier, LA, Dianzinga, RM, David R. Bowler

Published Oct 28, 2020• 4 citations
OpenAlex

Abstract

Density matrix perturbation theory (DMPT) is known as a promising alternative to the Rayleigh-Schrödinger perturbation theory, in which the sum-over-states (SOS) is replaced by algorithms with perturbed density matrices as the input variables. In this article, we formulate and discuss three types of DMPT, with two of them based only on density matrices: the approach of Kussmann and Ochsenfeld [J. Chem. Phys. 127, 054103 (2007)] is reformulated via the Sylvester equation and the recursive DMPT of Niklasson and Challacombe [Phys. Rev. Lett. 92, 193001 (2004)] is extended to the hole-particle canonical purification (HPCP) from Truflandier et al. [J. Chem. Phys. 144, 091102 (2016)]. A comparison of the computational performances shows that the aforementioned methods outperform the standard SOS. The HPCP-DMPT demonstrates stable convergence profiles but at a higher computational cost when compared to the original recursive polynomial method.

Related papers

Notes on density matrix perturbation theory

Lionel A. Truflandier, Rivo M. Dianzinga, David R. Bowler

Density matrix perturbation theory (DMPT) is known as a promising alternative to the Rayleigh-Schrödinger perturbation theory, in which the sum-over-states (SOS) is replaced by algorithms with perturbed density matrices as the input variables. In this article, we formulate and discuss three types of DMPT, with two of them based only on density matrices: the approach of Kussmann and Ochsenfeld [J. Chem. Phys. 127, 054103 (2007)] is reformulated via the Sylvester equation and the recursive DMPT of Niklasson and Challacombe [Phys. Rev. Lett. 92, 193001 (2004)] is extended to the hole-particle canonical purification (HPCP) from Truflandier et al. [J. Chem. Phys. 144, 091102 (2016)]. A comparison of the computational performances shows that the aforementioned methods outperform the standard SOS. The HPCP-DMPT demonstrates stable convergence profiles but at a higher computational cost when compared to the original recursive polynomial method.

arXiv:2007.04739•The Journal of Chemical Physics•Oct 23, 2020•5 citations

Density Matrix Perturbation Theory

Anders M. N. Niklasson, Matt Challacombe

An orbital-free quantum perturbation theory is proposed. It gives the response of the density matrix upon variation of the Hamiltonian by quadratically convergent recursions based on perturbed projections. The technique allows treatment of embedded quantum subsystems with a computational cost scaling linearly with the size of the perturbed region, O(N(pert.)), and as O(1) with the total system size. The method allows efficient high order perturbation expansions, as demonstrated with an example involving a 10th order expansion. Density matrix analogs of Wigner's 2n+1 rule are also presented.

arXiv:cond-mat/0311591•Physical Review Letters•May 14, 2004•78 citations

Susceptibility formulation of density matrix perturbation theory

Anders M. N. Niklasson, Adela Habib, Joshua Finkelstein +1 more

Density matrix perturbation theory based on recursive Fermi-operator expansions provides a computationally efficient framework for time-independent response calculations in quantum chemistry and materials science. From a perturbation in the Hamiltonian, we can calculate the first-order perturbation in the density matrix, which then gives us the linear response in the expectation values for some chosen set of observables. We present an alternative, dual formulation, where we instead calculate the static susceptibility of an observable, which then gives us the linear response in the expectation values for any number of different Hamiltonian perturbations. We show how the calculation of the susceptibility can be performed with the same expansion schemes used in recursive density matrix perturbation theory, including generalizations to fractional occupation numbers and self-consistent linear response calculations, i.e., similar to density functional perturbation theory. As with recursive density matrix perturbation theory, the dual susceptibility formulation is well suited for numerically thresholded sparse matrix algebra, which has linear scaling complexity for sufficiently large sparse systems. Similarly, the recursive computation of the susceptibility also seamlessly integrates with the computational framework of deep neural networks used in artificial intelligence (AI) applications. This integration enables the calculation of quantum response properties that can leverage cutting-edge AI-hardware, such as NVIDIA Tensor Cores or Google Tensor Processing Units. We demonstrate performance for recursive susceptibility calculations using NVIDIA Graphics Processing Units and Tensor Cores.

arXiv:2409.17033•The Journal of Chemical Physics•Dec 16, 2024•0 citations

Perturbation Theory and the Sum of Squares

Matthew B. Hastings

The sum-of-squares (SoS) hierarchy is a powerful technique based on semi-definite programming that can be used for both classical and quantum optimization problems. This hierarchy goes under several names; in particular, in quantum chemistry it is called the reduced density matrix (RDM) method. We consider the ability of this hierarchy to reproduce weak coupling perturbation theory for three different kinds of systems: spin (or qubit) systems, bosonic systems (the anharmonic oscillator), and fermionic systems with quartic interactions. For such fermionic systems, we show that degree-$4$ SoS (called $2$-RDM in quantum chemsitry) does not reproduce second order perturbation theory but degree-$6$ SoS ($3$-RDM) does (and we conjecture that it reproduces third order perturbation theory). Indeed, we identify a fragment of degree-$6$ SoS which can do this, which may be useful for practical quantum chemical calculations as it may be possible to implement this fragment with less cost than the full degree-$6$ SoS. Remarkably, this fragment is very similar to one studied by Hastings and O'Donnell for the Sachdev-Ye-Kitaev (SYK) model.

arXiv:2205.12325•May 24, 2022•1 citation

Regular article Equation for the direct determination of the density matrix: Time-dependent density equation and perturbation theory*

Hiroshi Nakatsuji

The density equation proposed previously for the direct determination of the density matrix, i.e. for the wave mechanics without wave, is extended to a time- dependent case. The time-dependent density equation has been shown to be equivalent to the time-dependent Schrodinger equation so long as the density matrix, included as a self-contained variable, is N-representable. Formally, it is obtainable from the previous time- independent equation by replacing the energy E with ih @ @t . The perturbation theory formulas for the density equation have also been given for both the time- dependent and time-independent cases.

•Jan 1, 1999•0 citations

Quantum-mechanical perturbation theory

John P. Killingbeck

Modern computers have made possible the evaluation of higher order terms of perturbation series. Apart from the empirical work involving numerical calculations, much formal mathematical work has been done on the convergence properties of perturbation series and of Pade approximants which represent them. Some of the modern work in perturbation theory is reviewed in the context of traditional results from the theory of real and complex variables. The two major versions of time-independent perturbation theory, the Rayleigh-Schrodinger (RS) and Brillouin-Wigner (BW) theories, are compared and contrasted, and the alternative techniques for evaluating the terms in the energy series are discussed. The sum-over-states method, the differential equation method and the variational principle method are treated, with emphasis on their inter-relations.

Reports on Progress in Physics•Sep 1, 1977•163 citations
Density matrix perturbation theory (DMPT) is known as a promising alternative to the Rayleigh-Schrödinger perturbation theory, in which the sum-over-states (SOS) is replaced by algorithms with perturbed density matrices as the input variables. In this article, we formulate and discuss three types of DMPT, with two of them based only on density matrices: the approach of Kussmann and Ochsenfeld [J. Chem. Phys. 127, 054103 (2007)] is reformulated via the Sylvester equation and the recursive DMPT of Niklasson and Challacombe [Phys. Rev. Lett. 92, 193001 (2004)] is extended to the hole-particle canonical purification (HPCP) from Truflandier et al. [J. Chem. Phys. 144, 091102 (2016)]. A comparison of the computational performances shows that the aforementioned methods outperform the standard SOS. The HPCP-DMPT demonstrates stable convergence profiles but at a higher computational cost when compared to the original recursive polynomial method.
Density matrix perturbation theory (DMPT) is known as a promising alternative to the Rayleigh-Schrödinger perturbation theory, in which the sum-over-states (SOS) is replaced by algorithms with perturbed density matrices as the input variables. In this article, we formulate and discuss three types of DMPT, with two of them based only on density matrices: the approach of Kussmann and Ochsenfeld [J. Chem. Phys. 127, 054103 (2007)] is reformulated via the Sylvester equation and the recursive DMPT of Niklasson and Challacombe [Phys. Rev. Lett. 92, 193001 (2004)] is extended to the hole-particle canonical purification (HPCP) from Truflandier et al. [J. Chem. Phys. 144, 091102 (2016)]. A comparison of the computational performances shows that the aforementioned methods outperform the standard SOS. The HPCP-DMPT demonstrates stable convergence profiles but at a higher computational cost when compared to the original recursive polynomial method.
An orbital-free quantum perturbation theory is proposed. It gives the response of the density matrix upon variation of the Hamiltonian by quadratically convergent recursions based on perturbed projections. The technique allows treatment of embedded quantum subsystems with a computational cost scaling linearly with the size of the perturbed region, O(N(pert.)), and as O(1) with the total system size. The method allows efficient high order perturbation expansions, as demonstrated with an example involving a 10th order expansion. Density matrix analogs of Wigner's 2n+1 rule are also presented.
Density matrix perturbation theory based on recursive Fermi-operator expansions provides a computationally efficient framework for time-independent response calculations in quantum chemistry and materials science. From a perturbation in the Hamiltonian, we can calculate the first-order perturbation in the density matrix, which then gives us the linear response in the expectation values for some chosen set of observables. We present an alternative, dual formulation, where we instead calculate the static susceptibility of an observable, which then gives us the linear response in the expectation values for any number of different Hamiltonian perturbations. We show how the calculation of the susceptibility can be performed with the same expansion schemes used in recursive density matrix perturbation theory, including generalizations to fractional occupation numbers and self-consistent linear response calculations, i.e., similar to density functional perturbation theory. As with recursive density matrix perturbation theory, the dual susceptibility formulation is well suited for numerically thresholded sparse matrix algebra, which has linear scaling complexity for sufficiently large sparse systems. Similarly, the recursive computation of the susceptibility also seamlessly integrates with the computational framework of deep neural networks used in artificial intelligence (AI) applications. This integration enables the calculation of quantum response properties that can leverage cutting-edge AI-hardware, such as NVIDIA Tensor Cores or Google Tensor Processing Units. We demonstrate performance for recursive susceptibility calculations using NVIDIA Graphics Processing Units and Tensor Cores.
The sum-of-squares (SoS) hierarchy is a powerful technique based on semi-definite programming that can be used for both classical and quantum optimization problems. This hierarchy goes under several names; in particular, in quantum chemistry it is called the reduced density matrix (RDM) method. We consider the ability of this hierarchy to reproduce weak coupling perturbation theory for three different kinds of systems: spin (or qubit) systems, bosonic systems (the anharmonic oscillator), and fermionic systems with quartic interactions. For such fermionic systems, we show that degree-444 SoS (called 222-RDM in quantum chemsitry) does not reproduce second order perturbation theory but degree-666 SoS (333-RDM) does (and we conjecture that it reproduces third order perturbation theory). Indeed, we identify a fragment of degree-666 SoS which can do this, which may be useful for practical quantum chemical calculations as it may be possible to implement this fragment with less cost than the full degree-666 SoS. Remarkably, this fragment is very similar to one studied by Hastings and O'Donnell for the Sachdev-Ye-Kitaev (SYK) model.
The density equation proposed previously for the direct determination of the density matrix, i.e. for the wave mechanics without wave, is extended to a time- dependent case. The time-dependent density equation has been shown to be equivalent to the time-dependent Schrodinger equation so long as the density matrix, included as a self-contained variable, is N-representable. Formally, it is obtainable from the previous time- independent equation by replacing the energy E with ih @ @t . The perturbation theory formulas for the density equation have also been given for both the time- dependent and time-independent cases.
Modern computers have made possible the evaluation of higher order terms of perturbation series. Apart from the empirical work involving numerical calculations, much formal mathematical work has been done on the convergence properties of perturbation series and of Pade approximants which represent them. Some of the modern work in perturbation theory is reviewed in the context of traditional results from the theory of real and complex variables. The two major versions of time-independent perturbation theory, the Rayleigh-Schrodinger (RS) and Brillouin-Wigner (BW) theories, are compared and contrasted, and the alternative techniques for evaluating the terms in the energy series are discussed. The sum-over-states method, the differential equation method and the variational principle method are treated, with emphasis on their inter-relations.