pith. machine review for the scientific record. sign in

arxiv: 2601.16133 · v2 · submitted 2026-01-22 · ⚛️ nucl-th

Recognition: no theorem link

Computational schemes for the Magnus expansion of the in-medium similarity renormalization group

Matthias Heinz
Authors on Pith no claims yet

Pith reviewed 2026-05-16 11:55 UTC · model grok-4.3

classification ⚛️ nucl-th
keywords IMSRGMagnus expansionnuclear many-body methodsIMSRG(3)hunter-gatherer schemeground-state energiesexcitation energies
0
0 comments X

The pith

The hunter-gatherer scheme for IMSRG(3) approximations changes nuclear ground-state energies by up to 7 MeV compared to standard IMSRG(2) methods.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper examines the numerical uncertainty introduced by the hunter-gatherer scheme when solving the IMSRG flow equations in a factorized approximation to three-body operators. It reports concrete differences in energies obtained with this scheme versus conventional two-body truncations. A reader would care because the observed shifts can reach the same scale as the three-body corrections the approximation is meant to include, directly affecting the reliability of computed nuclear properties.

Core claim

The hunter-gatherer scheme for solving the IMSRG equations in the context of factorized IMSRG(3) approximations leads to differences of up to 7 MeV in ground-state energies and 0.5 MeV in excitation energies relative to conventional IMSRG(2) calculations. These discrepancies are in some cases comparable to the expected magnitude of full IMSRG(3) corrections.

What carries the argument

The hunter-gatherer scheme, an approximate method for integrating the IMSRG flow equations while capturing leading three-body operator effects at two-body cost.

If this is right

  • Ground-state energies can shift by as much as 7 MeV when the hunter-gatherer scheme is used.
  • Excitation energies can shift by up to 0.5 MeV.
  • These shifts reach the size of the IMSRG(3) corrections the scheme is intended to approximate.
  • The scheme's uncertainty must be tracked when using factorized IMSRG(3) for precision nuclear calculations.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Switching to more accurate integration methods for the flow equations could reduce these shifts at modest extra cost.
  • Repeating the comparisons for a wider set of nuclei would clarify when the scheme's errors are largest.
  • Hybrid approaches that use the hunter-gatherer scheme only when three-body contributions are small could improve overall accuracy.

Load-bearing premise

The reported energy differences arise primarily from the hunter-gatherer scheme itself rather than from basis choices, truncation artifacts, or other implementation details in the IMSRG(2) reference calculations.

What would settle it

Running identical IMSRG(2) calculations with an exact non-approximate integration scheme and comparing the resulting energies to those from the hunter-gatherer scheme would show whether the 7 MeV and 0.5 MeV differences persist.

Figures

Figures reproduced from arXiv: 2601.16133 by Matthias Heinz.

Figure 1
Figure 1. Figure 1: FIG. 1. Comparison of split Magnus and hunter-gatherer ap [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Comparison of ground-state energies and charge radii of [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Energies from the IMSRG(2) as a function of the flow [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
read the original abstract

The in-medium similarity renormalization group (IMSRG) is a popular many-body method used for computations of nuclei. It solves the many-body Schr\"odinger equation through a continuous unitary transformation of the many-body Hamiltonian. The IMSRG transformation is typically truncated at the normal-ordered two-body level, the IMSRG(2), but recently several approaches have been developed to capture the effects of normal-ordered three-body operators, the IMSRG(3). In particular, a factorized approximation to the IMSRG(3) proposes to capture the leading effects of three-body operators at the same computational cost as the IMSRG(2) approximation. This approach often employs an approximate scheme for solving the IMSRG equations, the so-called hunter-gatherer scheme. In this work, I study the uncertainty associated with this scheme. I find that the hunter-gatherer scheme differs by up to $7\,\mathrm{MeV}$ for ground-state energies and $0.5\,\mathrm{MeV}$ for excitation energies from standard IMSRG(2) approaches. These differences are in some cases comparable to the expected size of IMSRG(3) corrections.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 1 minor

Summary. The manuscript investigates computational schemes for the Magnus expansion in the in-medium similarity renormalization group (IMSRG), focusing on the hunter-gatherer scheme used within factorized IMSRG(3) approximations. Through direct numerical comparisons, it reports that this scheme produces differences of up to 7 MeV in ground-state energies and 0.5 MeV in excitation energies relative to standard IMSRG(2) calculations, and concludes that these differences can be comparable in magnitude to expected IMSRG(3) corrections.

Significance. If the reported differences can be isolated to the hunter-gatherer scheme, the work would provide a useful quantification of numerical uncertainty in approximate IMSRG solvers. This is relevant for the nuclear many-body community, as it helps assess the reliability of IMSRG(2) and factorized IMSRG(3) results when higher-order corrections are estimated. The emphasis on concrete numerical comparisons is a constructive element of the study.

major comments (1)
  1. [Abstract and numerical results section] The central claim attributes energy differences of up to 7 MeV (ground states) and 0.5 MeV (excitations) to the hunter-gatherer scheme. However, the manuscript does not explicitly demonstrate that all other variables—single-particle basis, normal-ordering truncation, flow-equation solver tolerances, and numerical convergence criteria—are held fixed between the hunter-gatherer and standard IMSRG(2) calculations. This isolation is required to support the conclusion that the differences are comparable to IMSRG(3) corrections; without it, the spread could include implementation artifacts.
minor comments (1)
  1. [Abstract] The abstract refers to 'standard IMSRG(2) approaches' without naming the specific implementations or the nuclei studied; adding these details would improve reproducibility and context.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their constructive feedback and careful reading of our manuscript. We address the major comment below and agree to revise the text to make the numerical isolation explicit.

read point-by-point responses

  1. Referee: [Abstract and numerical results section] The central claim attributes energy differences of up to 7 MeV (ground states) and 0.5 MeV (excitations) to the hunter-gatherer scheme. However, the manuscript does not explicitly demonstrate that all other variables—single-particle basis, normal-ordering truncation, flow-equation solver tolerances, and numerical convergence criteria—are held fixed between the hunter-gatherer and standard IMSRG(2) calculations. This isolation is required to support the conclusion that the differences are comparable to IMSRG(3) corrections; without it, the spread could include implementation artifacts.

    Authors: We thank the referee for highlighting this point. In the calculations presented, the single-particle basis, normal-ordering truncation (two-body level for the IMSRG(2) reference), flow-equation solver tolerances, and numerical convergence criteria were identical between the hunter-gatherer and standard IMSRG(2) runs; the only controlled difference was the use of the hunter-gatherer scheme itself. This setup was chosen precisely to isolate the scheme's effect. We acknowledge that the manuscript does not state this isolation explicitly. In the revised version we will add a dedicated paragraph in the numerical results section confirming that all other parameters were held fixed, and we will insert a brief clarifying sentence in the abstract and introduction. These changes will directly address the concern and strengthen the attribution of the observed energy differences to the hunter-gatherer scheme. revision: yes

Circularity Check

0 steps flagged

No circularity in numerical comparison study

full rationale

The paper reports direct numerical comparisons of ground-state and excitation energies between the hunter-gatherer scheme and standard IMSRG(2) approaches, with differences up to 7 MeV and 0.5 MeV. No derivations, predictions, or results reduce to inputs by construction, self-definition, or fitted parameters. The central claims rest on computational outputs from the IMSRG equations rather than any self-referential or load-bearing self-citation chain. The study is self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Review based on abstract only; no explicit free parameters, new axioms, or invented entities are stated. The work rests on standard assumptions of the IMSRG flow and normal-ordering truncations already present in the cited literature.

pith-pipeline@v0.9.0 · 5492 in / 1052 out tokens · 23400 ms · 2026-05-16T11:55:07.499926+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Ab initio calculation of symmetry-breaking observables

    nucl-th 2026-05 unverdicted novelty 7.0

    A new IMSRG variant computes ab initio anapole and Schiff moments in medium-mass nuclei, benchmarked against no-core shell model results in light systems.

Reference graph

Works this paper leans on

43 extracted references · 43 canonical work pages · cited by 1 Pith paper · 1 internal anchor

  1. [1]

    Hergert, A guided tour of ab initio nuclear many-body theory, Front

    H. Hergert, A guided tour of ab initio nuclear many-body theory, Front. Phys.8, 379 (2020)

    work page 2020

  2. [2]

    T. D. Morris, J. Simonis, S. R. Stroberg, C. Stumpf, G. Hagen, J. D. Holt, G. R. Jansen, T. Papenbrock, R. Roth, and A. Schwenk, Structure of the lightest tin isotopes, Phys. Rev. Lett.120, 152503 (2018)

    work page 2018

  3. [3]

    Gysbers, G

    P. Gysbers, G. Hagen, J. D. Holt, G. R. Jansen, T. D. Mor- ris, P. Navr´ atil, T. Papenbrock, S. Quaglioni, A. Schwenk, S. R. Stroberg, and K. A. Wendt, Discrepancy between experimental and theoretical β-decay rates resolved from first principles, Nat. Phys.15, 428 (2019)

    work page 2019

  4. [4]

    Som` a, P

    V. Som` a, P. Navr´ atil, F. Raimondi, C. Barbieri, and T. Duguet, Novel chiral Hamiltonian and observables in light and medium-mass nuclei, Phys. Rev. C101, 014318 (2020)

    work page 2020

  5. [5]

    G. B. King, L. Andreoli, S. Pastore, M. Piarulli, R. Schi- avilla, R. B. Wiringa, J. Carlson, and S. Gandolfi, Chiral effective field theory calculations of weak transitions in light nuclei, Phys. Rev. C102, 025501 (2020)

    work page 2020

  6. [6]

    S. R. Stroberg, J. D. Holt, A. Schwenk, and J. Simonis, Ab initio limits of atomic nuclei, Phys. Rev. Lett.126, 022501 (2021)

    work page 2021

  7. [7]

    B. S. Hu, W. G. Jiang, T. Miyagi, Z. H. Sun, A. Ekstr¨ om, C. Forss´ en, G. Hagen, J. D. Holt, T. Papenbrock, S. R. Stroberg, and I. Vernon, Ab initio predictions link the neutron skin of 208Pb to nuclear forces, Nat. Phys.18, 1196 (2022)

    work page 2022

  8. [8]

    Tichai, P

    A. Tichai, P. Demol, and T. Duguet, Towards heavy-mass ab initio nuclear structure: Open-shell Ca, Ni and Sn isotopes from Bogoliubov coupled-cluster theory, Phys. Lett. B851, 138571 (2024)

    work page 2024

  9. [9]

    Miyagi, X

    T. Miyagi, X. Cao, R. Seutin, S. Bacca, R. F. Garcia Ruiz, K. Hebeler, J. D. Holt, and A. Schwenk, Impact of two- body currents on magnetic dipole moments of nuclei, Phys. Rev. Lett.132, 232503 (2024)

    work page 2024

  10. [10]

    Porro, T

    A. Porro, T. Duguet, J.-P. Ebran, M. Frosini, R. Roth, and V. Som` a, Ab initio description of monopole resonances in light- and medium-mass nuclei: II. Ab initio PGCM calculations in 46Ti, 28Si and 24Mg, Eur. Phys. J. A60, 134 (2024)

    work page 2024

  11. [11]

    Miyagi, M

    T. Miyagi, M. Heinz, and A. Schwenk, Ab initio computa- tions of the fourth-order charge density moments of 48Ca and 208Pb, Phys. Lett. B872, 140032 (2026)

    work page 2026

  12. [12]

    Bonaiti, G

    F. Bonaiti, G. Hagen, and T. Papenbrock, Structure of the doubly magic nuclei 208Pb and 266Pb from ab initio computations, arXiv:2508.14217

    work page arXiv

  13. [13]

    Bonaiti, C

    F. Bonaiti, C. Balos, K. Godbey, G. Hagen, T. Papen- brock, and C. S. Woodward, Computing nuclear response functions with time-dependent coupled-cluster theory (2026), Phys. Rev. C, in press, arXiv:2510.19940

    work page arXiv 2026

  14. [14]

    Gazda, R

    D. Gazda, R. Catena, and C. Forss´ en, Ab initio nuclear response functions for dark matter searches, Phys. Rev. D95, 103011 (2017)

    work page 2017

  15. [15]

    Andreoli, V

    L. Andreoli, V. Cirigliano, S. Gandolfi, and F. Pederiva, Quantum Monte Carlo calculations of dark matter scat- tering off light nuclei, Phys. Rev. C99, 025501 (2019)

    work page 2019

  16. [16]

    C. G. Payne, S. Bacca, G. Hagen, W. Jiang, and T. Pa- penbrock, Coherent elastic neutrino-nucleus scattering on 40Ar from first principles, Phys. Rev. C100, 061304 (2019)

    work page 2019

  17. [17]

    Novario, P

    S. Novario, P. Gysbers, J. Engel, G. Hagen, G. R. Jansen, T. D. Morris, P. Navr´ atil, T. Papenbrock, and 6 S. Quaglioni, Coupled-cluster calculations of neutrinoless double-β decay in 48Ca, Phys. Rev. Lett.126, 182502 (2021)

    work page 2021

  18. [18]

    B. S. Hu, J. Padua-Arg¨ uelles, S. Leutheusser, T. Miyagi, S. R. Stroberg, and J. D. Holt, Ab initio structure factors for spin-dependent dark matter direct detection, Phys. Rev. Lett.128, 072502 (2022)

    work page 2022

  19. [19]

    Glick-Magid, C

    A. Glick-Magid, C. Forss´ en, D. Gazda, D. Gazit, P. Gys- bers, and P. Navr´ atil, Nuclear ab initio calculations of6He β-decay for beyond the Standard Model studies, Phys. Lett. B832, 137259 (2022)

    work page 2022

  20. [20]

    Cirigliano, W

    V. Cirigliano, W. Dekens, J. de Vries, S. Gandolfi, M. Hoferichter, and E. Mereghetti, Ab initio electroweak corrections to superallowed β decays and their impact on Vud, Phys. Rev. C110, 055502 (2024)

    work page 2024

  21. [21]

    Belley, J

    A. Belley, J. M. Yao, B. Bally, J. Pitcher, J. Engel, H. Hergert, J. D. Holt, T. Miyagi, T. R. Rodr´ ıguez, A. M. Romero, S. R. Stroberg, and X. Zhang, Ab initio uncer- tainty quantification of neutrinoless double-beta decay in 76Ge, Phys. Rev. Lett.132, 182502 (2024)

    work page 2024

  22. [22]

    Farren and S

    A. Farren and S. R. Stroberg, Spurious Isospin Break- ing in the In-medium Similarity Renormalization Group, arXiv:2412.10693

    work page arXiv

  23. [23]

    Gennari, M

    M. Gennari, M. Drissi, M. Gorchtein, P. Navratil, and C.-Y. Seng, Ab initio strategy for taming the nuclear- structure dependence of Vud extractions: The 10C→10B superallowed transition, Phys. Rev. Lett.134, 012501 (2025)

    work page 2025

  24. [24]

    Door, C.-H

    M. Door, C.-H. Yeh, M. Heinz, F. Kirk, C. Lyu,et al., Probing new bosons and nuclear structure with ytterbium isotope shifts, Phys. Rev. Lett.134, 063002 (2025)

    work page 2025

  25. [25]

    Heinz, M

    M. Heinz, M. Hoferichter, T. Miyagi, F. N¨ oel, and A. Schwenk, Ab initio calculations of overlap integrals for µ→e conversion in nuclei, Phys. Lett. B871, 139975 (2025)

    work page 2025

  26. [26]

    K. B. Ng, S. Foster, L. Cheng, P. Navr´ atil, and S. Malbrunot-Ettenauer, Nuclear Schiff moment of fluo- rine isotope 19F, arXiv:2507.19811

    work page arXiv

  27. [27]

    G. B. King, J. Carlson, A. R. Flores, S. Gandolfi, E. Mereghetti, S. Pastore, M. Piarulli, and R. B. Wiringa, Quantum Monte Carlo calculation of δNS in 10C using an effective field theory approach, arXiv:2509.07310

    work page arXiv

  28. [28]

    Z. Li, T. Miyagi, and A. Schwenk, Ab initio calculations of beta-decay half-lives for N = 50 neutron-rich nuclei, arXiv:2509.06812

    work page internal anchor Pith review Pith/arXiv arXiv

  29. [29]

    Kuske, T

    J. Kuske, T. Miyagi, A. Arcones, and A. Schwenk, r- process nucleosynthesis with ab initio nuclear masses around theN= 82 shell closure, arXiv:2509.19131

    work page arXiv

  30. [30]

    Tsukiyama, S

    K. Tsukiyama, S. K. Bogner, and A. Schwenk, In-medium similarity renormalization group for nuclei, Phys. Rev. Lett.106, 222502 (2011)

    work page 2011

  31. [31]

    Hergert, S

    H. Hergert, S. K. Bogner, T. D. Morris, A. Schwenk, and K. Tsukiyama, The in-medium similarity renormalization group: A novel ab initio method for nuclei, Phys. Rep. 621, 165 (2016)

    work page 2016

  32. [32]

    S. R. Stroberg, A. Calci, H. Hergert, J. D. Holt, S. K. Bogner, R. Roth, and A. Schwenk, Nucleus-dependent valence-space approach to nuclear structure, Phys. Rev. Lett.118, 032502 (2017)

    work page 2017

  33. [33]

    S. R. Stroberg, H. Hergert, S. K. Bogner, and J. D. Holt, Nonempirical interactions for the nuclear shell model: An update, Annu. Rev. Nucl. Part. Sci.69, 307 (2019)

    work page 2019

  34. [34]

    Heinz, A

    M. Heinz, A. Tichai, J. Hoppe, K. Hebeler, and A. Schwenk, In-medium similarity renormalization group with three-body operators, Phys. Rev. C103, 044318 (2021)

    work page 2021

  35. [35]

    S. R. Stroberg, T. D. Morris, and B. C. He, In-medium similarity renormalization group with flowing 3-body op- erators, and approximations thereof, Phys. Rev. C110, 044316 (2024)

    work page 2024

  36. [36]

    Heinz, T

    M. Heinz, T. Miyagi, S. R. Stroberg, A. Tichai, K. Hebeler, and A. Schwenk, Improved structure of calcium isotopes from ab initio calculations, Phys. Rev. C111, 034311 (2025)

    work page 2025

  37. [37]

    B. C. He and S. R. Stroberg, Factorized approximation to the in-medium similarity renormalization group IM- SRG(3), Phys. Rev. C110, 044317 (2024)

    work page 2024

  38. [38]

    T. D. Morris, N. M. Parzuchowski, and S. K. Bogner, Mag- nus expansion and in-medium similarity renormalization group, Phys. Rev. C92, 034331 (2015)

    work page 2015

  39. [39]

    Hebeler, S

    K. Hebeler, S. K. Bogner, R. J. Furnstahl, A. Nogga, and A. Schwenk, Improved nuclear matter calculations from chiral low-momentum interactions, Phys. Rev. C83, 031301(R) (2011)

    work page 2011

  40. [40]

    Miyagi,NuHamil: A numerical code to generate nu- clear two- and three-body matrix elements from chiral effective field theory, Eur

    T. Miyagi,NuHamil: A numerical code to generate nu- clear two- and three-body matrix elements from chiral effective field theory, Eur. Phys. J. A59, 150 (2023)

    work page 2023

  41. [41]

    Hagen, T

    G. Hagen, T. Papenbrock, M. Hjorth-Jensen, and D. J. Dean, Coupled-cluster computations of atomic nuclei, Rep. Prog. Phys.77, 096302 (2014)

    work page 2014

  42. [42]

    Heinz, Data: On computational schemes for the mag- nus expansion of the in-medium similarity renormalization group, 10.5281/zenodo.18341412 (2026)

    M. Heinz, Data: On computational schemes for the mag- nus expansion of the in-medium similarity renormalization group, 10.5281/zenodo.18341412 (2026)

    work page doi:10.5281/zenodo.18341412 2026

  43. [43]

    S. R. Stroberget al., imsrg++, https://github.com/ ragnarstroberg/imsrg(2026)

    work page 2026