In-Medium Similarity Renormalization Group (IMSRG)

• IMSRG is an ab initio many-body framework that applies continuous unitary transformations to nuclear Hamiltonians to suppress off-diagonal couplings and decouple targeted sectors.
• It employs normal ordering with respect to a chosen reference state and allows systematic truncations like IMSRG(2) and IMSRG(3) for practical calculations in nuclear structure.
• The method has been extended to include multireference formulations, deformation, continuum and finite-temperature effects, thereby enhancing the predictive power for diverse nuclear phenomena.

The In-Medium Similarity Renormalization Group (IMSRG) is an ab initio many-body framework that applies a continuous unitary transformation to a nuclear Hamiltonian, normal-ordered with respect to a chosen reference state, in order to suppress off-diagonal couplings and decouple targeted sectors of the many-body problem. In practice, the method has been used both to extract correlated ground-state energies and to derive effective Hamiltonians and operators for valence-space calculations, while later extensions have incorporated correlated references, deformation, continuum degrees of freedom, explicit normal-ordered three-body operators, infinite matter, and finite temperature (Hergert et al., 2015, Hergert, 2016, Hergert et al., 2018).

1. Formal basis of the IMSRG

At its core, the IMSRG defines a flow of Hamiltonians

$$H(s)=U(s)H(0)U^\dagger(s), \qquad \frac{dH(s)}{ds}=[\eta(s),H(s)],$$

with anti-Hermitian generator

$$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$

The flow parameter $s$ labels a continuous family of unitarily equivalent Hamiltonians. The guiding idea is to partition the Hamiltonian into “diagonal” and “off-diagonal” pieces, $H=H^{d}+H^{od}$, and choose $\eta(s)$ so that $H^{od}(s)\to 0$ as $s\to\infty$ (Hergert et al., 2015, Hergert et al., 2018).

The in-medium aspect enters through normal ordering with respect to a finite-density reference state $|\Phi\rangle$, typically a Hartree–Fock Slater determinant in single-reference applications. In second quantization, the normal-ordered Hamiltonian is written as

$$H = E +\sum_{ij} f_{ij}\{a_i^\dagger a_j\} +\frac14\sum_{ijkl}\Gamma_{ijkl}\{a_i^\dagger a_j^\dagger a_l a_k\} +\frac1{36}\sum_{ijklmn}W_{ijklmn}\{a_i^\dagger a_j^\dagger a_k^\dagger a_n a_m a_l\},$$

where $E$, $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$0, $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$1, and $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$2 are the zero-, one-, two-, and three-body normal-ordered terms. Normal ordering folds substantial medium dependence, including contributions from three-nucleon forces, into low-rank operators (Hergert et al., 2015, Hergert et al., 2012).

For closed-shell ground-state decoupling, the off-diagonal Hamiltonian is usually defined as the set of matrix elements that couple the reference to $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$3 and $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$4 excitations,

$$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$5

so that the evolved zero-body term $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$6 approximates the correlated ground-state energy (Hergert et al., 2018, Hergert, 2016). In valence-space formulations, the decoupling target is instead a shell-model $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$7 space, and the effective Hamiltonian is the projected evolved operator,

$$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$8

which can then be diagonalized in configuration interaction (Tsukiyama et al., 2012, Hergert, 2016).

2. Truncation hierarchy, generators, and the Magnus formulation

The exact flow generates operators of arbitrarily high rank, so practical IMSRG calculations rely on systematic truncations. The standard IMSRG(2) retains normal-ordered zero-, one-, and two-body operators and discards explicit three-body terms during the flow. In conjunction with the normal-ordered two-body approximation (NO2B), this has become the workhorse formulation for medium-mass nuclei (Hergert et al., 2012, Hergert et al., 2015). Reviews of the perturbative content emphasize that IMSRG(2) reproduces the complete energy through third order in MBPT and a subset of fourth-order contributions, while remaining polynomial-scaling (Hergert et al., 2015).

A full IMSRG(3) keeps the normal-ordered three-body operator $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$9 explicitly. The 2021 extension to IMSRG(3) established that the method is fourth-order complete in MBPT in the NO2B setting and benchmarked a hierarchy of lower-cost approximations, including IMSRG(3)-MP4, IMSRG(3)-$s$0, IMSRG(3)-$s$1, and IMSRG(3)-$s$2 (Heinz et al., 2021). A later calcium study used the IMSRG(3)-$s$3 truncation in realistic spaces and found that it yields smaller many-body uncertainties and a significantly better description of the first $s$4 excitation energy of $s$5Ca than IMSRG(2), while leaving problematic calcium isotope shifts largely unresolved (Heinz et al., 2024).

Generator choice governs both numerical behavior and the character of the decoupling. The Wegner generator,

$s$6

guarantees monotonic suppression of off-diagonal couplings but can produce stiff differential equations (Hergert et al., 2015). White-type generators use ratios of off-diagonal matrix elements to energy denominators and typically yield faster, less stiff flows. Imaginary-time generators replace the denominator by $s$7, while Brillouin generators are particularly natural for correlated multireference states because they avoid an explicit diagonal/off-diagonal partition (Hergert et al., 2014, Hergert et al., 2018). In deformed single-reference IMSRG, the White generator is used together with the off-diagonal set $s$8 in the deformed basis (Yuan et al., 2022).

The Magnus formulation rewrites the unitary transformation as

$s$9

so that

$H=H^{d}+H^{od}$0

and similarly for any operator $H=H^{d}+H^{od}$1. This reformulation reduces memory requirements, permits low-order ODE integrators without loss of accuracy in the transformed observables, and makes consistent operator evolution essentially routine once $H=H^{d}+H^{od}$2 has been obtained (Morris et al., 2015, Yuan et al., 2022).

3. Reference states, multireference formulations, and symmetry restoration

The simplest IMSRG implementation uses a single Slater determinant reference and is therefore most natural for closed-shell nuclei. Open-shell and collective systems require more elaborate references. The multireference IMSRG (MR-IMSRG) generalizes normal ordering and Wick’s theorem to correlated reference states by introducing irreducible density matrices $H=H^{d}+H^{od}$3 and has been used with particle-number-projected Hartree–Fock–Bogoliubov references to describe even calcium and nickel isotopes with chiral NN+3N interactions (Hergert et al., 2014, Hergert, 2016).

The generalized formalism expands the Hamiltonian as

$H=H^{d}+H^{od}$4

with cumulants such as

$H=H^{d}+H^{od}$5

entering both the flow equations and the Brillouin generator (Yao et al., 2018, Hergert et al., 2014). In MR-IMSRG(2), terms nonlinear in $H=H^{d}+H^{od}$6 and those involving $H=H^{d}+H^{od}$7 may be dropped depending on the reference and implementation, whereas in IMSRG+GCM calculations inclusion of $H=H^{d}+H^{od}$8 in the Brillouin generator is necessary for convergence in $H=H^{d}+H^{od}$9Ti (Yao et al., 2018).

A distinct strategy for open-shell nuclei is to break rotational symmetry at the reference level. The deformed IMSRG (D-IMSRG) uses an axially deformed Hartree–Fock Slater determinant in an $\eta(s)$0-scheme basis, thereby constructing a single-reference IMSRG for even-even open-shell systems. The corresponding reference preserves axial, parity, and time-reversal symmetries while breaking spherical rotational symmetry. Static rotational correlations are then approximated by angular-momentum projection of the HF state,

$\eta(s)$1

with leading-order correction

$\eta(s)$2

This avoids the much more expensive full projection of the correlated IMSRG state (Yuan et al., 2022).

Another extension uses symmetry-restored generator-coordinate-method references. IMSRG+GCM starts from deformed HFB states projected on good angular momentum and particle number,

$\eta(s)$3

and then applies MR-IMSRG(2) with a Brillouin generator. In pf-shell benchmarks this improves low-lying spectra relative to either spherical-reference IMSRG or GCM with unevolved operators, and it provides a consistent route to evolved $\eta(s)$4 operators via the PI and PF prescriptions (Yao et al., 2018).

Variant	Reference state	Primary target
Single-reference IMSRG	HF Slater determinant	Closed-shell ground states
MR-IMSRG	Projected HFB or correlated reference	Open-shell ground states
VS-IMSRG	Core reference plus valence-space decoupling	Shell-model Hamiltonians
D-IMSRG	Axially deformed HF Slater determinant	Deformed even-even nuclei
IMSRG+GCM	Symmetry-restored GCM states	Collective spectra and $\eta(s)$5

4. Evolved operators, excited states, and transition observables

A central technical advantage of the IMSRG is that observables can be evolved consistently with the Hamiltonian. In Magnus form,

$\eta(s)$6

and matrix elements are evaluated in the decoupled reference or in post-IMSRG many-body states (Yuan et al., 2022, Parzuchowski et al., 2017). This is essential for charge radii, moments, transition strengths, and weak matrix elements.

For charge radii, D-IMSRG uses the intrinsic decomposition

$\eta(s)$7

with $\eta(s)$8, $\eta(s)$9, and $H^{od}(s)\to 0$0 in that implementation (Yuan et al., 2022). A more recent calcium study used

$H^{od}(s)\to 0$1

together with a corrected spin-orbit term proportional to $H^{od}(s)\to 0$2 (Heinz et al., 2024).

Excited states can be accessed either by deriving effective interactions for subsequent configuration-interaction diagonalization or by working directly with the decoupled Hamiltonian. EOM-IMSRG formulates

$H^{od}(s)\to 0$3

with a ladder operator truncated, for example, at the $H^{od}(s)\to 0$4 level. TDA-IMSRG and valence-restricted TDA-IMSRG instead engineer additional decouplings so that a Tamm–Dancoff diagonalization becomes exact in a chosen sector, though EOM-IMSRG is more robust numerically (Parzuchowski et al., 2016).

For electromagnetic observables, the consistent IMSRG evolution of tensor operators was implemented in both EOM-IMSRG(2,2) and VS-IMSRG(2). Magnetic dipole observables are generally in reasonable agreement with experiment, whereas the more collective electric quadrupole and octupole observables are significantly underpredicted, often by over an order of magnitude, indicating missing physics at the present truncation level (Parzuchowski et al., 2017). In the pf-shell benchmark for neutrinoless double-beta decay, IMSRG+GCM gives better spectra than either spherical-reference IMSRG or GCM with unevolved operators, while the improvement in the $H^{od}(s)\to 0$5 matrix element is slight in the single-shell setup; the dominant collective isoscalar-pairing correlations must be included explicitly in the GCM coordinates (Yao et al., 2018).

5. Finite nuclei: closed shells, valence spaces, deformation, and continuum

Applications to finite nuclei span closed-shell energies, open-shell spectroscopy, deformation, and open quantum systems. Early chiral NN+3N calculations for closed shells up to $H^{od}(s)\to 0$6Ni showed good agreement with experiment in $H^{od}(s)\to 0$7He and the closed-shell oxygen isotopes, while calcium and nickel isotopes were somewhat overbound; comparison with coupled-cluster and importance-truncated no-core shell model established the accuracy of IMSRG(2) and the role of the NO2B approximation (Hergert et al., 2012). MR-IMSRG calculations across even calcium and nickel isotopes subsequently demonstrated the importance of chiral 3N interactions for correct shell evolution and dripline trends, with agreement at shell closures comparable to CC results obtained with the same Hamiltonians (Hergert et al., 2014).

The valence-space IMSRG provides a nonperturbative derivation of shell-model Hamiltonians. The original open-shell application derived effective $H^{od}(s)\to 0$8- and $H^{od}(s)\to 0$9-shell Hamiltonians for $s\to\infty$0Li and $s\to\infty$1O, obtaining a spectrum for $s\to\infty$2Li in very good agreement with ab initio results and establishing a direct alternative to diagrammatic effective-interaction theory for shell-model studies (Tsukiyama et al., 2012). Later reviews emphasized successful spectroscopy in the lower $s\to\infty$3 shell and access to intrinsically deformed states through shell-model diagonalization (Hergert, 2016).

For deformed nuclei, D-IMSRG calculations with NNLO$s\to\infty$4 showed that projected HF corrections lower the energies of $s\to\infty$5Be by about $s\to\infty$6–$s\to\infty$7 MeV, with specific corrections

$s\to\infty$8

and gains of about $s\to\infty$9–$|\Phi\rangle$0 MeV in Ne and Mg near the $|\Phi\rangle$1 island of inversion. In that study, both D-IMSRG and VS-IMSRG predicted the even-even Be dripline at $|\Phi\rangle$2Be, whereas experiment places it at $|\Phi\rangle$3Be; the omission of continuum coupling was explicitly identified as a likely source of the discrepancy (Yuan et al., 2022).

Continuum degrees of freedom can be incorporated directly through the Berggren basis. The Gamow IMSRG reformulates the flow for a complex-symmetric Hamiltonian, using a complex-orthogonal similarity transformation in a basis that satisfies the Berggren completeness relation. In $|\Phi\rangle$4O it reproduces the observed resonant structure, including the triplet of $|\Phi\rangle$5, $|\Phi\rangle$6, and $|\Phi\rangle$7 states around $|\Phi\rangle$8 MeV. In $|\Phi\rangle$9C it predicts a pronounced halo density and demonstrates that the neutron $$H = E +\sum_{ij} f_{ij}\{a_i^\dagger a_j\} +\frac14\sum_{ijkl}\Gamma_{ijkl}\{a_i^\dagger a_j^\dagger a_l a_k\} +\frac1{36}\sum_{ijklmn}W_{ijklmn}\{a_i^\dagger a_j^\dagger a_k^\dagger a_n a_m a_l\},$$0 continuum is essential: the matter radius increases from $$H = E +\sum_{ij} f_{ij}\{a_i^\dagger a_j\} +\frac14\sum_{ijkl}\Gamma_{ijkl}\{a_i^\dagger a_j^\dagger a_l a_k\} +\frac1{36}\sum_{ijklmn}W_{ijklmn}\{a_i^\dagger a_j^\dagger a_k^\dagger a_n a_m a_l\},$$1 to $$H = E +\sum_{ij} f_{ij}\{a_i^\dagger a_j\} +\frac14\sum_{ijkl}\Gamma_{ijkl}\{a_i^\dagger a_j^\dagger a_l a_k\} +\frac1{36}\sum_{ijklmn}W_{ijklmn}\{a_i^\dagger a_j^\dagger a_k^\dagger a_n a_m a_l\},$$2 fm with NNLO$$H = E +\sum_{ij} f_{ij}\{a_i^\dagger a_j\} +\frac14\sum_{ijkl}\Gamma_{ijkl}\{a_i^\dagger a_j^\dagger a_l a_k\} +\frac1{36}\sum_{ijklmn}W_{ijklmn}\{a_i^\dagger a_j^\dagger a_k^\dagger a_n a_m a_l\},$$3 and from $$H = E +\sum_{ij} f_{ij}\{a_i^\dagger a_j\} +\frac14\sum_{ijkl}\Gamma_{ijkl}\{a_i^\dagger a_j^\dagger a_l a_k\} +\frac1{36}\sum_{ijklmn}W_{ijklmn}\{a_i^\dagger a_j^\dagger a_k^\dagger a_n a_m a_l\},$$4 to $$H = E +\sum_{ij} f_{ij}\{a_i^\dagger a_j\} +\frac14\sum_{ijkl}\Gamma_{ijkl}\{a_i^\dagger a_j^\dagger a_l a_k\} +\frac1{36}\sum_{ijklmn}W_{ijklmn}\{a_i^\dagger a_j^\dagger a_k^\dagger a_n a_m a_l\},$$5 fm with NNLO$$H = E +\sum_{ij} f_{ij}\{a_i^\dagger a_j\} +\frac14\sum_{ijkl}\Gamma_{ijkl}\{a_i^\dagger a_j^\dagger a_l a_k\} +\frac1{36}\sum_{ijklmn}W_{ijklmn}\{a_i^\dagger a_j^\dagger a_k^\dagger a_n a_m a_l\},$$6 when the $$H = E +\sum_{ij} f_{ij}\{a_i^\dagger a_j\} +\frac14\sum_{ijkl}\Gamma_{ijkl}\{a_i^\dagger a_j^\dagger a_l a_k\} +\frac1{36}\sum_{ijklmn}W_{ijklmn}\{a_i^\dagger a_j^\dagger a_k^\dagger a_n a_m a_l\},$$7-wave is promoted from discrete states to a Berggren continuum (Hu et al., 2019).

Recent IMSRG(3)-$$H = E +\sum_{ij} f_{ij}\{a_i^\dagger a_j\} +\frac14\sum_{ijkl}\Gamma_{ijkl}\{a_i^\dagger a_j^\dagger a_l a_k\} +\frac1{36}\sum_{ijklmn}W_{ijklmn}\{a_i^\dagger a_j^\dagger a_k^\dagger a_n a_m a_l\},$$8 calculations for $$H = E +\sum_{ij} f_{ij}\{a_i^\dagger a_j\} +\frac14\sum_{ijkl}\Gamma_{ijkl}\{a_i^\dagger a_j^\dagger a_l a_k\} +\frac1{36}\sum_{ijklmn}W_{ijklmn}\{a_i^\dagger a_j^\dagger a_k^\dagger a_n a_m a_l\},$$9Ca sharpened the method’s quantitative status. For $E$0Ca, the first $E$1 excitation in VS-IMSRG(2) is $E$2 MeV, and the IMSRG(3)-$E$3 correction is $E$4 MeV, moving the result toward the experimental $E$5 MeV and improving the description of the $E$6 shell closure. By contrast, charge-radius corrections are small: for $E$7Ca the single-reference radius changes from $E$8 fm by only $E$9 fm, and the large experimental increase in $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$00 remains unexplained (Heinz et al., 2024).

6. Infinite matter and finite-temperature extensions

The IMSRG has also been extended beyond finite nuclei. A 2025 development formulated IMSRG(2) directly in homogeneous infinite nuclear matter using a plane-wave basis around a filled Fermi sea, normal ordering at finite density, and a White generator with Epstein–Nesbet denominators. The reference occupations are those of the Fermi sphere, and the Hamiltonian is written as

$$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$01

For pure neutron matter with the soft NNLO$$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$02 interaction, IMSRG(2) and CCD agree very well up to $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$03 for NN only and up to $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$04 when 3N forces are included through NO2B. For a harder N$$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$05LO interaction with $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$06 MeV, pronounced discrepancies emerge among MBPT3, MBPT4, and IMSRG(2), especially above $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$07, while in symmetric nuclear matter even soft interactions already show large differences among MBPT, CCD, and IMSRG(2) near and above saturation density (Zhen et al., 5 Jan 2025).

A finite-temperature extension generalizes the reference from a zero-temperature Slater determinant to a thermal Hartree–Fock ensemble, with occupations

$$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$08

in the grand canonical ensemble. The IMSRG(2) flow equations retain their form, but the integer occupations are replaced by thermal factors. Because a straightforward White generator can become ill-behaved when occupations appear in both numerators and denominators, the finite-temperature implementation adopts an arctan-regularized White generator with occupation-dependent denominators and preserves the IMSRG flow equation

$$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$09

in the thermal setting (Smith et al., 2024).

Benchmarking on an exactly solvable pairing-plus-particle–hole model showed that FT-IMSRG substantially improves over FT-HF across temperatures, with relative errors typically below the percent level for weak to moderate couplings and the largest deviations at intermediate $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$10. Canonical and grand canonical implementations differ more strongly in small systems, but those differences diminish rapidly with particle number, which suggests that the grand canonical ensemble is the practical choice for realistic nuclear applications (Smith et al., 2024).

7. Computational strategies, uncertainty quantification, and outlook

Because IMSRG calculations are dominated by tensor contractions in the two- and three-body sectors, algorithmic developments have become a major subfield. Importance truncation for IMSRG(2) selects only “important” two-body matrix elements according to a measure $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$11 and leaves the rest out of the flow. The most effective strategy in the benchmark study was to truncate only the $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$12 and $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$13 blocks, using either

$$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$14

or

$$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$15

together with a perturbative third-order correction $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$16 as an error indicator. For soft Hamiltonians and large single-particle spaces, compression ratios up to $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$17 are feasible with sub-MeV discrepancies, while hard interactions produce larger and more conservative error estimates (Hoppe et al., 2021).

Machine-learning acceleration has also been explored. IMSRG-Net learns the Magnus operator $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$18 directly from a short training window of the flow, using a loss function that combines data fidelity for $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$19 and its derivative-based generator proxy. In benchmarks for $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$20O and $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$21Ca, training on ten points up to $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$22 sufficed to reproduce full IMSRG(2) ground-state energies with typical absolute differences below $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$23 keV and charge radii with typical absolute differences of order $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$24 fm. The same strategy was successfully applied to the second step of VS-IMSRG decoupling (Yoshida, 2023).

The numerical treatment of the Magnus expansion itself remains nontrivial. A 2026 analysis of the hunter-gatherer scheme showed that, although it reduces storage to two Magnus operators, it can differ from standard IMSRG(2) approaches by up to $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$25 for ground-state energies and $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$26 for excitation energies, with the discrepancies in some cases comparable to expected IMSRG(3) corrections. The paper therefore recommends split Magnus or direct-flow schemes when the target accuracy is at the level of IMSRG(3) improvements, especially in valence-space applications with large decouplings (Heinz, 22 Jan 2026).

Uncertainty quantification is now increasingly explicit. In calcium, IMSRG(3)-$$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$27 corrections imply practical IMSRG(2) uncertainty estimates of approximately $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$28–$$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$29 for correlation energies, $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$30–$$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$31 for charge radii, and $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$32–$$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$33 for neutron skins for soft interactions such as 1.8/2.0 (EM) (Heinz et al., 2024). Earlier MR-IMSRG studies estimated beyond-(2) effects in two-neutron separation energies at the few-hundred-keV level by comparison to CCSD and CR-CC(2,3), while NO2B overbinding in calcium and nickel was found to be below $$\eta(s)=\frac{dU(s)}{ds}U^\dagger(s)=-\eta^\dagger(s).$$34 in total energies (Hergert et al., 2014). These benchmarks suggest that many-body truncation, omitted residual and induced three-body terms, continuum coupling, and deficiencies of current chiral Hamiltonians contribute at comparable levels depending on the observable.

The emerging picture is that the IMSRG is no longer a single approximation but a family of controlled transformations. Single-reference, multireference, valence-space, deformed, Gamow, finite-temperature, and infinite-matter variants all preserve the same central idea—continuous in-medium decoupling—while differing in reference choice, generator, and operator-rank truncation. Current developments point toward fuller symmetry restoration, more systematic treatment of explicit and induced three-body operators, broader continuum embeddings, and tighter integration of uncertainty quantification with the evolving hierarchy from IMSRG(2) to IMSRG(3) and beyond (Hergert et al., 2018, Heinz et al., 2021, Heinz et al., 2024).