pith. sign in

arxiv: 2607.00681 · v2 · pith:DEFMQK6Tnew · submitted 2026-07-01 · ⚛️ nucl-th · hep-lat

Elastic deuteron-deuteron scattering within Nuclear Lattice Effective Field Theory

Helen Meyer , Serdar Elhatisari , Fabian Hildenbrand , Ulf-G. Mei{\ss}ner This is my paper

Pith reviewed 2026-07-03 18:07 UTC · model grok-4.3

classification ⚛️ nucl-th hep-lat
keywords deuteron-deuteron scatteringnuclear lattice effective field theoryeffective range parametersquintet channeladiabatic projection methodchiral interactionsbig-bang nucleosynthesis
0
0 comments X

The pith

Nuclear lattice EFT yields a deuteron-deuteron scattering length of 12.96 fm with stronger quintet-channel repulsion.

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

The paper computes low-energy deuteron-deuteron scattering in the spin-quintet 5S2 channel using nuclear lattice effective field theory. It combines next-to-next-to-next-to-leading order chiral interactions implemented via wavefunction matching with the adiabatic projection method. Two stabilization procedures address small norm-matrix eigenvalues at large Euclidean projection times and produce consistent phase shifts. The extracted values show more negative phase shifts and a substantially larger scattering length than previous calculations, indicating stronger effective repulsion. These results supply the first nuclear-lattice benchmark for deuteron-deuteron scattering and a foundation for future coupled-channel studies of deuteron-induced reactions.

Core claim

The Coulomb-modified effective-range analysis of the computed phase shifts gives 5a_dd = (12.96 ± 0.26) fm and 5r_dd = (3.62 ± 0.79) fm. The phase shifts are more negative and the scattering length substantially larger than in previous calculations, corresponding to a stronger effective repulsion in the 5S2 channel.

What carries the argument

Adiabatic projection method stabilized by Tikhonov regularization or projection onto well-resolved norm eigenmodes, applied to N3LO chiral interactions.

If this is right

  • The calculation supplies the first nuclear-lattice benchmark for deuteron-deuteron scattering.
  • It establishes a basis for future coupled-channel calculations of deuteron-induced reactions relevant to big-bang nucleosynthesis.
  • The two stabilization procedures produce consistent Coulomb-subtracted phase shifts within statistical and numerical uncertainties.
  • The larger scattering length reflects stronger effective repulsion in the 5S2 channel relative to earlier work.

Where Pith is reading between the lines

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

  • The stronger repulsion may alter predicted rates for deuteron-involved reactions in big-bang nucleosynthesis.
  • The lattice approach could be extended to other partial waves or to include explicit electromagnetic effects beyond the Coulomb subtraction used here.
  • Direct comparison with low-energy dd scattering data would test whether the reported effective-range parameters hold outside the lattice framework.

Load-bearing premise

The two stabilization procedures correctly remove the effect of small norm-matrix eigenvalues at large Euclidean projection time without introducing systematic bias into the extracted phase shifts.

What would settle it

An independent calculation or direct experimental measurement of the quintet deuteron-deuteron scattering length that differs significantly from 12.96 fm would falsify the central result.

Figures

Figures reproduced from arXiv: 2607.00681 by Fabian Hildenbrand, Helen Meyer, Serdar Elhatisari, Ulf-G. Mei{\ss}ner.

Figure 1
Figure 1. Figure 1: Threshold energy of two non-interacting deuterons on the lattice as [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Coulomb-subtracted phase shift of the projected single-channel [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Fitting the modified effective range expansion (ERE) to eq. (28). project by providing computing time on the GCS Supercom￾puters JUWELS and JUPITER at Jülich Supercomputing Cen￾tre (JSC) and the support of the project EXOTIC by the JSC by dedicated HPC time provided on the JURECA DC GPU par￾tition. Furthermore, the authors gratefully acknowledge the computing time provided on the high-performance computer … view at source ↗
read the original abstract

We calculate low-energy deuteron-deuteron scattering in the spin-quintet $^{5}S_2$ channel using nuclear lattice effective field theory. The calculation combines chiral interactions at next-to-next-to-next-to-leading order, implemented through wavefunction matching, with the adiabatic projection method. Because the radial cluster basis develops small norm-matrix eigenvalues at large Euclidean projection time, we investigate two stabilization procedures: Tikhonov regularization and projection onto well-resolved norm eigenmodes. The two procedures yield consistent Coulomb-subtracted phase shifts within their statistical and numerical uncertainties. A Coulomb-modified effective-range analysis gives ${}^5a_{dd} = (12.96 \pm 0.26)\,\mathrm{fm}$ and ${}^5r_{dd} = (3.62 \pm 0.79)\,\mathrm{fm}$. The phase shifts are more negative, and the scattering length is substantially larger than in previous calculations, corresponding to a stronger effective repulsion in the $^{5}S_2$ channel. These results provide a first nuclear-lattice benchmark for deuteron-deuteron scattering and establish a basis for future coupled-channel calculations of the deuteron-induced reactions relevant to big-bang nucleosynthesis.

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 / 0 minor

Summary. The manuscript computes low-energy elastic deuteron-deuteron scattering in the spin-quintet 5S2 channel within nuclear lattice effective field theory. It employs N3LO chiral interactions via wavefunction matching together with the adiabatic projection method on the lattice. Two stabilization procedures (Tikhonov regularization and projection onto well-resolved norm eigenmodes) are applied to the radial cluster norm matrix at large Euclidean times; the resulting Coulomb-subtracted phase shifts are fitted to a Coulomb-modified effective-range expansion, producing 5a_dd = (12.96 ± 0.26) fm and 5r_dd = (3.62 ± 0.79) fm. These values are larger (more repulsive) than those from previous calculations and are presented as a first nuclear-lattice benchmark for dd scattering relevant to big-bang nucleosynthesis.

Significance. If the extracted phase shifts prove free of systematic bias, the work supplies the first lattice-EFT benchmark for deuteron-deuteron scattering and a foundation for future coupled-channel calculations of deuteron-induced reactions. The use of N3LO interactions, direct numerical solution of the lattice Schrödinger equation, and internal consistency between two stabilization methods constitute clear technical strengths.

major comments (1)
  1. [Stabilization procedures paragraph (abstract and methods)] Stabilization procedures paragraph (abstract and methods): The central claim of a substantially larger 5a_dd rests on the Coulomb-subtracted phase shifts obtained after Tikhonov regularization and projection onto well-resolved norm eigenmodes. Mutual statistical consistency of the two procedures does not exclude a shared systematic bias in the low-momentum regime; no external validation (exactly solvable test case, third independent regulator, or controlled limit with vanishing regularization) is described that would bound any such bias below the quoted 0.26 fm uncertainty.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and for the constructive comments. We address the major comment on the stabilization procedures below.

read point-by-point responses

  1. Referee: The central claim of a substantially larger 5a_dd rests on the Coulomb-subtracted phase shifts obtained after Tikhonov regularization and projection onto well-resolved norm eigenmodes. Mutual statistical consistency of the two procedures does not exclude a shared systematic bias in the low-momentum regime; no external validation (exactly solvable test case, third independent regulator, or controlled limit with vanishing regularization) is described that would bound any such bias below the quoted 0.26 fm uncertainty.

    Authors: We agree that external validation would strengthen the bound on possible shared systematic bias. The two procedures are mathematically distinct (one damps small eigenvalues via a penalty term while the other discards them by mode selection), and their agreement within uncertainties already reduces the likelihood of a common bias, but this does not fully exclude it. In the revised manuscript we will add a controlled test on an exactly solvable two-body problem to quantify the residual bias from each procedure and confirm it lies below the reported 0.26 fm uncertainty. revision: yes

Circularity Check

0 steps flagged

Direct numerical lattice computation; no reduction by construction

full rationale

The paper computes deuteron-deuteron phase shifts by solving the lattice Schrödinger equation with N3LO chiral interactions (via wavefunction matching) and the adiabatic projection method, followed by two stabilization procedures on the norm matrix whose mutual consistency is reported. The effective-range parameters are then obtained by a standard fit to the extracted Coulomb-subtracted phase shifts. No equation in the manuscript equates a derived quantity to an input parameter by definition, renames a fitted result as a prediction, or imports a uniqueness theorem from the authors' prior work that would force the central result. The calculation is self-contained against external benchmarks and receives score 0.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The calculation rests on N3LO chiral EFT interactions whose low-energy constants were fitted in prior work, plus the assumption that the chosen stabilization methods preserve the physical scattering observables.

free parameters (1)
  • low-energy constants of N3LO chiral interaction
    Determined by fits to nucleon-nucleon data in earlier publications; used here via wavefunction matching.
axioms (2)
  • domain assumption Chiral effective field theory at N3LO provides an accurate description of low-energy nuclear forces when implemented on the lattice
    Invoked to generate the two-body interactions for the deuteron-deuteron system.
  • domain assumption The adiabatic projection method extracts the correct asymptotic scattering information from the lattice wave functions
    Central to obtaining phase shifts from the Euclidean-time evolution.

pith-pipeline@v0.9.1-grok · 5760 in / 1454 out tokens · 29648 ms · 2026-07-03T18:07:44.174449+00:00 · methodology

Review history (2 revisions) →

discussion (0)

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

Reference graph

Works this paper leans on

42 extracted references · 42 canonical work pages · 13 internal anchors

  1. [1]

    R. V . Wagoner, W. A. Fowler, F. Hoyle, On the Synthesis of elements at very high temperatures, Astrophys. J. 148 (1967) 3–49. doi:10.1086/149126

    work page doi:10.1086/149126 1967

  2. [2]

    R. H. Cyburt, B. D. Fields, K. A. Olive, T.- H. Yeh, Big Bang Nucleosynthesis: 2015, Rev. Mod. Phys. 88 (2016) 015004. arXiv:1505.01076, doi:10.1103/RevModPhys.88.015004

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1103/revmodphys.88.015004 2015

  3. [3]

    Pitrou, A

    C. Pitrou, A. Coc, J.-P. Uzan, E. Vangioni, Resolving con- clusions about the early Universe requires accurate nu- clear measurements, Nature Rev. Phys. 3 (4) (2021) 231–

    work page 2021

  4. [4]

    arXiv:2104.11148, doi:10.1038/s42254-021-00294- 6

    work page doi:10.1038/s42254-021-00294-

  5. [5]

    Burns, T

    A.-K. Burns, T. M. P. Tait, M. Valli, PRyMordial: the first three minutes, within and beyond the standard model, Eur. Phys. J. C 84 (1) (2024) 86. arXiv:2307.07061, doi:10.1140/epjc/s10052-024-12442-0

    work page doi:10.1140/epjc/s10052-024-12442-0 2024

  6. [6]

    Y . Xu, S. Goriely, A. Jorissen, G. Chen, M. Arnould, Databases and tools for nuclear astrophysics applica- tions BRUSsels Nuclear LIBrary (BRUSLIB), Nuclear Astrophysics Compilation of REactions II (NACRE 6 II) and Nuclear NETwork GENerator (NETGEN), As- tron. Astrophys. 549 (2013) A106. arXiv:1212.0628, doi:10.1051/0004-6361/201220537

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1051/0004-6361/201220537 2013

  7. [7]

    Pitrou, A

    C. Pitrou, A. Coc, J.-P. Uzan, E. Vangioni, Precision big bang nucleosynthesis with improved Helium-4 predic- tions, Phys. Rept. 754 (2018) 1–66. arXiv:1801.08023, doi:10.1016/j.physrep.2018.04.005

    work page doi:10.1016/j.physrep.2018.04.005 2018

  8. [8]

    T. A. Lähde, U.-G. Meißner, Nuclear Lattice Effective Field Theory: An introduction, V ol. 957, Springer, 2019. doi:10.1007/978-3-030-14189-9

    work page doi:10.1007/978-3-030-14189-9 2019

  9. [9]

    Meier, W

    W. Meier, W. Glöckle, Elastic scattering 2 H(d, d) 2 H below 360 keV, Nucl. Phys. A 255 (1975) 21–34. doi:10.1016/0375-9474(75)90144-X

    work page doi:10.1016/0375-9474(75)90144-x 1975

  10. [10]

    Wildermuth, W

    K. Wildermuth, W. McClure, Cluster Representa- tions of Nuclei, Springer Berlin, Heidelberg, 1966. doi:10.1007/BFb0045472

    work page doi:10.1007/bfb0045472 1966

  11. [11]

    Wildermuth, Y

    K. Wildermuth, Y . C. Tang, A Unified Theory of the Nucleus, Vieweg+Teubner Verlag, Wiesbaden, 1977. doi:10.1007/978-3-322-85255-7

    work page doi:10.1007/978-3-322-85255-7 1977

  12. [12]

    R. A. Malfliet, J. A. Tjon, Solution of the Faddeev equations for the triton problem using local two par- ticle interactions, Nucl. Phys. A 127 (1969) 161–168. doi:10.1016/0375-9474(69)90775-1

    work page doi:10.1016/0375-9474(69)90775-1 1969

  13. [13]

    I. N. Filikhin, S. L. Yakovlev, Microscopic calculation of low-energy deuteron-deuteron scattering on the basis of the cluster-reduction method, Phys. Atom. Nucl. 63 (2000). doi:10.1134/1.855624

    work page doi:10.1134/1.855624 2000

  14. [14]

    S. L. Yakovlev, I. N. Filikhin, Calculations of scatter- ing lengths in four nucleon system on the basis of clus- ter reduction method for Yakubovsky equations (1 1997). arXiv:nucl-th/9701020

    work page internal anchor Pith review Pith/arXiv arXiv 1997

  15. [15]

    H. M. Hofmann, G. M. Hale, Microscopic calcula- tion of the He-4 system, Nucl. Phys. A 613 (1997) 69–106. arXiv:nucl-th/9608046, doi:10.1016/S0375- 9474(96)00418-6

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1016/s0375- 1997

  16. [16]

    Kellermann, H

    H. Kellermann, H. M. Hofmann, C. Elster, Gaussian Pa- rameterization of a Meson TheoreticalNNPotential for Microscopic Nuclear Structure Calculations, Acta Phys. Austriaca 7 (1989) 31–53

    work page 1989

  17. [17]

    H. M. Hofmann, G. M. Hale, He-4 can experiments serve as a database for determining the three-nucleon force?, Phys. Rev. C 77 (2008) 044002. arXiv:nucl-th/0512065, doi:10.1103/PhysRevC.77.044002

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1103/physrevc.77.044002 2008

  18. [18]

    R. B. Wiringa, V . G. J. Stoks, R. Schiavilla, An Accu- rate nucleon-nucleon potential with charge independence breaking, Phys. Rev. C 51 (1995) 38–51. arXiv:nucl- th/9408016, doi:10.1103/PhysRevC.51.38

    work page doi:10.1103/physrevc.51.38 1995

  19. [19]

    B. S. Pudliner, V . R. Pandharipande, J. Carlson, S. C. Pieper, R. B. Wiringa, Quantum Monte Carlo calculations of nuclei with A<=7, Phys. Rev. C 56 (1997) 1720–1750. arXiv:nucl-th/9705009, doi:10.1103/PhysRevC.56.1720

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1103/physrevc.56.1720 1997

  20. [20]

    J. F. Carew, Deuteron-deuteron elastic and three- and four- body breakup scattering using the Faddeev-Yakubovskii equations, Phys. Rev. C 103 (1) (2021) 014002. doi:10.1103/PhysRevC.103.014002

    work page doi:10.1103/physrevc.103.014002 2021

  21. [21]

    J. F. Carew, Variational bounds in n-particle scatter- ing using the faddeev–yakubovskii equations: Deuteron- deuteron s=2 scattering, Few-Body Systems 55 (2014) 171–190. doi:10.1007/s00601-014-0844-0

    work page doi:10.1007/s00601-014-0844-0 2014

  22. [22]

    Modern Theory of Nuclear Forces

    E. Epelbaum, H.-W. Hammer, U.-G. Meißner, Modern Theory of Nuclear Forces, Rev. Mod. Phys. 81 (2009) 1773–1825. arXiv:0811.1338, doi:10.1103/RevModPhys.81.1773

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1103/revmodphys.81.1773 2009

  23. [23]

    Elhatisari, et al., Wavefunction matching for solv- ing quantum many-body problems, Nature 630 (8015) (2024) 59–63

    S. Elhatisari, et al., Wavefunction matching for solv- ing quantum many-body problems, Nature 630 (8015) (2024) 59–63. arXiv:2210.17488, doi:10.1038/s41586- 024-07422-z

    work page doi:10.1038/s41586- 2024

  24. [24]

    Hildenbrand, S

    F. Hildenbrand, S. Elhatisari, U.-G. Meißner, H. Meyer, Z. Ren, A. Herten, M. Bode, Lattice Calculation of the Sn Isotopes near the Proton Dripline, Phys. Rev. Lett. 136 (6) (2026) 062501. arXiv:2509.08579, doi:10.1103/n7nt- s64t

    work page doi:10.1103/n7nt- 2026

  25. [25]

    Z. Ren, S. Elhatisari, U.-G. Meißner, Ab Initio Study of the Radii of Oxygen Isotopes, Phys. Rev. Lett. 135 (15) (2025) 152502. arXiv:2506.02597, doi:10.1103/y6s2- 43ym

    work page doi:10.1103/y6s2- 2025

  26. [26]

    Ab initio alpha-alpha scattering

    S. Elhatisari, D. Lee, G. Rupak, E. Epelbaum, H. Krebs, T. A. Lähde, T. Luu, U.-G. Meißner, Ab initio alpha-alpha scattering, Nature 528 (2015) 111. arXiv:1506.03513, doi:10.1038/nature16067

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1038/nature16067 2015

  27. [27]

    N. Li, S. Elhatisari, E. Epelbaum, D. Lee, B.-N. Lu, U.-G. Meißner, Neutron-proton scattering with lattice chiral effective field theory at next-to-next-to-next-to- leading order, Phys. Rev. C 98 (4) (2018) 044002. arXiv:1806.07994, doi:10.1103/PhysRevC.98.044002

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1103/physrevc.98.044002 2018

  28. [28]

    Elhatisari, T

    S. Elhatisari, T. A. Lähde, D. Lee, U.-G. Meißner, T. V onk, Alpha-alpha scattering in the Multi- verse, JHEP 02 (2022) 001. arXiv:2112.09409, doi:10.1007/JHEP02(2022)001

    work page doi:10.1007/jhep02(2022)001 2022

  29. [29]

    Elhatisari, F

    S. Elhatisari, F. Hildenbrand, U.-G. Meißner, Ab ini- tio lattice study of neutron–alpha scattering with chi- ral forces at N3LO, J. Phys. G 52 (12) (2025) 125102. arXiv:2507.08495, doi:10.1088/1361-6471/ae2145

    work page doi:10.1088/1361-6471/ae2145 2025

  30. [30]

    H. Tong, S. Elhatisari, U.-G. Meißner, Ab initio calcula- tion of hyper-neutron matter, Sci. Bull. 70 (2025) 825–

    work page 2025

  31. [31]

    arXiv:2405.01887, doi:10.1016/j.scib.2025.01.008. 7

    work page doi:10.1016/j.scib.2025.01.008 2025

  32. [32]

    Lee, Lattice Effective Field Theory Simulations of Nuclei, Ann

    D. Lee, Lattice Effective Field Theory Simulations of Nuclei, Ann. Rev. Nucl. Part. Sci. 75 (1) (2025) 109–128. arXiv:2501.03303, doi:10.1146/annurev-nucl- 101918-023343

    work page doi:10.1146/annurev-nucl- 2025

  33. [33]

    Ab initio $\alpha$-$\alpha$ scattering with high-fidelity chiral interactions

    A. Sarkar, S. Elhatisari, T. A. Lähde, U.-G. Meißner, Ab initioα-αscattering with high-fidelity chiral interactions (6 2026). arXiv:2606.28987

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  34. [34]

    M. Pine, D. Lee, G. Rupak, Adiabatic projection method for scattering and reactions on the lattice, Eur. Phys. J. A 49 (2013) 151. arXiv:1309.2616, doi:10.1140/epja/i2013- 13151-3

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1140/epja/i2013- 2013

  35. [35]

    Nucleon-deuteron scattering using the adiabatic projection method

    S. Elhatisari, D. Lee, U.-G. Meißner, G. Rupak, Nucleon- deuteron scattering using the adiabatic projection method, Eur. Phys. J. A 52 (6) (2016) 174. arXiv:1603.02333, doi:10.1140/epja/i2016-16174-2

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1140/epja/i2016-16174-2 2016

  36. [36]

    Elhatisari, Adiabatic projection method with Euclidean time subspace projection, Eur

    S. Elhatisari, Adiabatic projection method with Euclidean time subspace projection, Eur. Phys. J. A 55 (8) (2019)

    work page 2019

  37. [37]

    arXiv:1906.01046, doi:10.1140/epja/i2019-12844-9

    work page doi:10.1140/epja/i2019-12844-9 1906

  38. [38]

    A. N. Tikhonov, V . Y . Arsenin, Solutions of Ill-Posed Problems, V . H. Winston & Sons, Washington, DC, 1977, distributed by Halsted Press, New York

    work page 1977

  39. [39]

    D. E. Hilt, D. W. Seegrist, U. S. F. Service, P. Northeast- ern Forest Experiment Station (Radnor, Ridge, a computer program for calculating ridge regression estimates, V ol. no.236, Upper Darby, Pa, Dept. of Agriculture, Forest Ser- vice, Northeastern Forest Experiment Station, 1977

    work page 1977

  40. [40]

    P. C. Hansen, Rank-Deficient and Discrete Ill-Posed Prob- lems, Society for Industrial and Applied Mathematics,

    work page

  41. [41]

    doi:10.1137/1.9780898719697

    work page doi:10.1137/1.9780898719697

  42. [42]

    B.-N. Lu, T. A. Lähde, D. Lee, U.-G. Meißner, Precise determination of lattice phase shifts and mixing angles, Phys. Lett. B 760 (2016) 309–313. arXiv:1506.05652, doi:10.1016/j.physletb.2016.06.081. 8 Appendix A. Euclidean time extrapolation of phase shift results For each individual momentum point we calculate the phase shift as described in section 2.3 ...

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1016/j.physletb.2016.06.081 2016