---
title: 'DSRG-MRPT3: Third-Order MR Perturbation Theory'
url: https://www.emergentmind.com/topics/dsrg-mrpt3
type: topic
---

# DSRG-MRPT3: Third-Order MR Perturbation Theory

Searching arXiv for recent and foundational papers on DSRG-MRPT3 and related MR-DSRG methods.
DSRG-MRPT3 is the third-order multireference perturbation theory based on the driven similarity renormalization group (DSRG). In this formulation, dynamical correlation is added to a multiconfigurational reference through an \(s\)-dependent similarity transformation of the Hamiltonian, yielding a method that is intruder free, size consistent, non-iterative with \(\mathcal{O}(N^6)\) scaling, and inclusive of reference relaxation effects [1701.02011]. Within the broader MR-DSRG hierarchy, it is the third-order perturbative counterpart to iterative linearized schemes such as MR-LDSRG(2) [1602.05667].

## 1. Definition and formal framework

The DSRG framework applies a continuous unitary transformation to the electronic Hamiltonian,
\[
\bar{H}(s) = e^{-\hat{A}(s)} \hat{H}\, e^{\hat{A}(s)},
\qquad
\hat{A}(s)=\hat{T}(s)-\hat{T}^\dagger(s),
\]
where \(\hat{T}(s)\) is an internally contracted excitation operator defined with respect to a multiconfigurational reference [2405.06077]. The essential DSRG condition is not a differential flow equation but a driven algebraic constraint on the off-diagonal part of the transformed Hamiltonian,
\[
[\bar{H}(s)]_{\mathrm{od}}=\hat{R}(s),
\]
with \(\hat{R}(s)\) a source operator chosen to damp couplings associated with problematic denominators [1602.05667].

The flow parameter \(s\) has units of \(E_h^{-2}\) and defines an effective energy cutoff \(\Lambda=s^{-1/2}\) in the underlying SRG interpretation [1602.05667, 1406.0114]. In perturbative expressions, amplitudes acquire the regularized factor
\[
\frac{1-e^{-s\Delta^2}}{\Delta},
\]
so contributions from near-degenerate external configurations remain finite rather than diverging [1701.02011, 1602.05667]. This regularization is the central reason DSRG-based multireference perturbation theories are described as intruder resistant.

Historically, the single-reference DSRG was introduced as an algebraic alternative to SRG differential equations, replacing a flow of ordinary differential equations by polynomial equations for the transformed Hamiltonian [1406.0114]. The multireference extension, and later the first nonperturbative MR-LDSRG(2) model, recast this structure in generalized normal ordering relative to a complete active space reference, preserving connectedness and size-extensive behavior while avoiding the redundancy of projective internally contracted formulations [1602.05667].

## 2. Perturbative construction of third-order MR-DSRG

DSRG-MRPT3 is obtained by a perturbative analysis of the MR-LDSRG(2) equations truncated through third order in the perturbation [1701.02011]. The reference is a CASSCF wavefunction,
\[
\ket{\Phi_0} = \sum_{\mu=1}^{d} c_\mu \ket{\Phi^\mu},
\]
with orbitals partitioned into core, active, and virtual spaces, and semicanonicalized so that the generalized Fock matrix is diagonal within each block [1701.02011]. The zeroth-order Hamiltonian is chosen as
\[
\hat{H}^{(0)} = E_0 + \hat{F}^{(0)},
\]
which yields simple denominators and ensures that commutators with \(\hat{H}^{(0)}\) contribute only non-diagonal one- and two-body terms [1701.02011].

Within this perturbative expansion, the transformed Hamiltonian is expressed as
\[
\bar{H}(s)=\bar{H}^{(0)}(s)+\bar{H}^{(1)}(s)+\bar{H}^{(2)}(s)+\bar{H}^{(3)}(s)+\cdots,
\]
and the unrelaxed third-order energy is
\[
E_{\mathrm{u}}^{[3]}(s)=\sum_{n=0}^{3} E_{\mathrm{u}}^{(n)}(s).
\]
The first-order energy vanishes for the chosen partitioning, while the second- and third-order terms are built from commutators of \(\tilde{H}^{(1)}(s)\), \(\tilde{H}^{(2)}(s)\), and the first- and second-order cluster operators [1701.02011]. The formal structure mirrors coupled-cluster perturbation theory, but every denominator is renormalized by the DSRG damping functions.

A distinctive aspect of DSRG-MRPT3 is reference relaxation. After constructing the third-order transformed Hamiltonian truncated to one- and two-body operators,
\[
\bar{H}^{[3]}_{1,2}(s)=\sum_{n=0}^{3}\bar{H}^{(n)}_{1,2}(s),
\]
one diagonalizes this effective Hamiltonian in the CAS space to obtain a relaxed energy \(E^{[3]}(s)\) [1701.02011]. This one-step relaxation is inexpensive relative to the dominant external-space contractions and was found to reduce nonparallelity errors systematically.

The original DSRG-MRPT3 paper emphasized four formal properties: it is intruder free, size consistent, non-iterative with \(\mathcal{O}(N^6)\) scaling, and includes reference relaxation effects [1701.02011]. On bond dissociation curves of \(\mathrm{F}_2\), \(\mathrm{H}_2\mathrm{O}_2\), \(\mathrm{C}_2\mathrm{H}_6\), and \(\mathrm{N}_2\), its nonparallelism errors were consistent with CASPT3 and MRCISD, and showed significant improvements over DSRG second-order multireference perturbation theory [1701.02011].

## 3. Spin-free and relativistic formulations

A major formal extension is the spin-free MR-DSRG built from the ensemble normal ordering of Mukherjee and Kutzelnigg [2106.07097]. The ensemble averages over all microstates for a given total spin quantum number and is invariant with respect to SU(2) transformations, so the equations can be rewritten entirely in spin-free quantities that closely resemble spin-adapted closed-shell coupled-cluster equations [2106.07097]. In this framework, perturbation theory up to third order and iterative singles-doubles MR-DSRG variants were benchmarked on thirty-three first-row diatomic molecules, and focal point analysis for \([\mathrm{Fe}(\mathrm{H}_2\mathrm{O})_6]^{2+}\) and \([\mathrm{Fe}(\mathrm{NH}_3)_6]^{2+}\) showed that third-order perturbative corrections are essential to achieve reasonably converged energetics [2106.07097]. The resulting spin splittings were predicted to be \(-35.7\) and \(-17.1\) kcal mol\(^{-1}\), respectively [2106.07097].

DSRG-MRPT3 has also been extended to a fully relativistic four-component setting. The 4c-SA-DSRG-MRPT3 method starts from state-averaged 4c-CASSCF references and a Dirac–Coulomb–Breit Hamiltonian under the no-pair approximation, with molecular spinors and spin-orbit coupling built directly into the one- and two-electron operators [2405.06077]. In benchmarks on second- to fourth-row \(p\)-block atoms, 4c-SA-DSRG-MRPT3 generally yielded smaller mean absolute errors than 4c-CASSCF, 4c-CASPT2, and 4c-MR-CISD+Q, while showing reduced sensitivity to the flow parameter relative to 4c-SA-DSRG-MRPT2 [2405.06077]. The best overall MAE occurred near \(s=0.35\;E_h^{-2}\), and the method achieved sub-wavenumber accuracy for boron and carbon and errors below \(2\ \mathrm{cm}^{-1}\) for selenium and bromine splittings [2405.06077]. For the OH radical, MRPT3 also markedly improved spectroscopic constants relative to the reference and PT2 levels [2405.06077].

## 4. Large-active-space implementations

Recent work has pushed DSRG-MRPT3 into active spaces that were previously inaccessible for internally contracted multireference methods. A DMRG-based implementation avoids the explicit construction of high-order reduced density matrices by forming matrix-product-state compressed intermediates, allowing DSRG second- and third-order perturbation theories to be applied to dodecacene with an active space of 50 electrons in 50 orbitals [2503.01299]. This active space was described as the largest employed to date within the framework of internally contracted multireference formalism [2503.01299].

In the oligoacene series from naphthalene to dodecacene, DMRG-DSRG yielded a best estimate of the vertical singlet-triplet gap of dodecacene of \(0.22\ \mathrm{eV}\), in excellent agreement with the linearized adiabatic connection result of \(0.24\ \mathrm{eV}\) [2503.01299]. For zeaxanthin, all DSRG schemes predicted the vertical excited-state ordering
\[
2\, ^1 A_g^- < 1\, ^1 B_u^+ < 1\, ^1 B_u^-,
\]
a result of direct interest for carotenoid spectroscopy [2503.01299]. For the \(\mathrm{Cr}_2\) potential energy curve, both the equilibrium and shoulder regions were reasonably reproduced by linearized DSRG with one- and two-body operators [2503.01299].

The technical significance of this implementation lies in the way it reconciles internally contracted perturbation theory with matrix-product-state references. In conventional formulations, scalar terms in the BCH expansion require the 3-RDM, which becomes a memory bottleneck for large \(L\). The compressed-intermediate strategy replaces explicit 3-RDM construction by MPS objects built from amplitude-weighted annihilation strings, reducing memory requirements by at least a factor of \(L\) while preserving high numerical accuracy for DSRG-PT2 and PT3 energies [2503.01299]. This development suggests that DSRG-MRPT3 is not limited to traditional CAS references and can be integrated naturally with DMRG-SCF.

## 5. Spectroscopy, excited states, and ionization variants

DSRG-MRPT3 has become a practical parent Hamiltonian for several spectroscopic theories. In the XABOOM benchmark of K-edge \(1s\to\pi^\ast\) excitations, third-order corrections significantly improved the accuracy of GAS-DSRG absolute excitation energies, reducing the mean absolute deviation from experimental values to \(0.32\ \mathrm{eV}\) [2212.04369]. In the same study, DSRG-MRPT2 systematically underestimated absolute excitation energies, whereas DSRG-MRPT3 largely removed the systematic shift and proved more robust to active-space truncation and intruder-like behavior in challenging systems such as ozone and glyoxylic acid [2212.04369].

For valence ionization, the IP-EOM-DSRG formalism combines an EOM treatment of ionized states with three parent methods: DSRG-MRPT2, DSRG-MRPT3, and MR-LDSRG(2) [2506.13693]. The EOM step scales as \(\mathcal{O}(N^5)\) with basis size, and benchmarks on small molecules, radicals, and stretched geometries showed that all three IP-EOM-DSRG variants accurately reproduce vertical ionization potentials and spectroscopic constants, with the DSRG-MRPT3 and MR-LDSRG(2) versions outperforming several state-of-the-art multireference methods of comparable or higher cost [2506.13693]. In that work, EOM-DSRG-PT3 was also less sensitive to the flow parameter than EOM-DSRG-PT2 over \(s\in[0.5,1.0]\;E_h^{-2}\) [2506.13693].

An analogous core-ionization extension, CVS-IP-EOM-DSRG, was formulated for X-ray photoelectron spectra [2509.21646]. There the EOM step scales as \(\mathcal{O}(N^4)\) relative to basis size, and although all three parent Hamiltonians accurately predicted vertical core-ionization energies, only the DSRG-MRPT3 and MR-LDSRG(2) levels reliably captured the full dissociation behavior and reproduced the experimental vibrational structure of XPS spectra [2509.21646]. These developments indicate that DSRG-MRPT3 functions not only as a stand-alone energy model but also as a compact Hermitian effective Hamiltonian for excited-state and ionization theories.

## 6. Relation to adjacent methods, strengths, and limitations

DSRG-MRPT3 is best understood as occupying a middle ground between second-order multireference perturbation theories and iterative MR-DSRG models. Relative to DSRG-MRPT2, third-order corrections improve accuracy, narrow error distributions, and reduce sensitivity to the flow parameter in both relativistic and nonrelativistic applications [2405.06077, 1701.02011]. Relative to MR-LDSRG(2), DSRG-MRPT3 sacrifices iterative resummation in exchange for a non-iterative algorithm of lower practical cost, yet in many benchmarks it approaches the accuracy of the iterative method [2506.13693, 2503.01299].

Compared with CASPT2, NEVPT2, MR-CISD(+Q), and related approaches, the distinguishing feature of DSRG-MRPT3 is the renormalization of denominators by the flow parameter \(s\), rather than by a level shift or by a specially constructed zeroth-order Hamiltonian [1701.02011, 2212.04369]. In relativistic benchmarks, 4c-CASPT2 and 4c-MR-CISD+Q consistently underestimated spin-orbit splittings, whereas 4c-DSRG-MRPT3 yielded better agreement with experiment over a wide range of \(s\) [2405.06077]. In core-excitation benchmarks, DSRG-MRPT3 improved absolute energies substantially relative to PT2 and achieved benchmark-quality performance for organic K-edge spectra [2212.04369].

Its limitations are equally well established. Accuracy depends critically on the quality of the CASSCF or DMRG reference and on the chosen active space [2506.13693, 2503.01299]. Results still depend on the flow parameter, even if the dependence is weaker than in PT2 [2405.06077, 2506.13693]. The standard perturbative and linearized formulations truncate the transformed Hamiltonian to one- and two-body operators, so omitted three-body terms can become significant in strongly correlated regimes or at stretched geometries [1701.02011]. In EOM applications, limitations of the operator manifold also matter: the 2h1p ionization ansatz poorly describes states with strong 3h2p character regardless of the parent MR-DSRG level [2506.13693]. These caveats motivate ongoing work on triples, higher-body corrections, and renormalized internally contracted MRCC variants that reinterpret DSRG as a unitary multireference coupled-cluster theory and adapt its flow-equation regularization to nonunitary settings [2405.16139, 2005.10132].

In contemporary electronic-structure theory, DSRG-MRPT3 is therefore a mature third-order multireference perturbation model with a clear formal identity: it combines a multiconfigurational reference, internally contracted excitations, SRG-style denominator renormalization, and a relaxed effective Hamiltonian. Its demonstrated range now extends from bond dissociation and diradicals to core spectroscopy, spin-orbit splittings, XPS, large-active-space DMRG references, and ionized states, while preserving the methodological attributes that motivated its original introduction [1701.02011, 2503.01299, 2506.13693].

Source: https://www.emergentmind.com/topics/dsrg-mrpt3