REVIEW 3 major objections 4 minor 21 references
The breakdown scale of pionless effective field theory in the three-nucleon sector
T0 review · 3 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read The paper infers the breakdown scale of pionless effective field theory from neutron-deuteron cross sections, finding a next-to-next-to-leading-order mode near 100 MeV and a combined two- and three-nucleon posterior close to the pion mass.
desk verdict The nd-only breakdown-scale estimate is solid; the combined np+nd posterior overreaches on an unexamined independence assumption. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the order-by-order expansion $\sigma_{\mathrm{nd}}^{(\nu)}(k) = \sigma_{\mathrm{ref}}(k) \sum_{i=0}^{\nu} c_i(k) [f(k)/M_{\mathrm{hi}}]^i$, where $\sigma_{\mathrm{ref}}(k)$ is the leading-order cross section, $c_i(k)$ are dimensionless EFT coefficients, and $Q(k)=f(k)/M_{\mathrm{hi}}$ is the expansion parameter built from a soft-max-interpolated soft scale $f(k)$. Under the likelihood, the coefficients $c_i$ are treated as independent, identically distributed draws from a normal prior of natural size, so the posterior for $M_{\mathrm{hi}}$ follows from how natural the observed coefficients $c_1$ and $c_2$ look for each trial value of $M_{\mathrm{hi}}$. The cross sections themselves come from solving the pionless EFT three-body integral equations for the neutron-deuteron $t$-matrix, with a three-nucleon interaction fitted to the triton binding energy at leading and next-to-leading order and to the neutron-deuteron scattering length at next-to-next-to-leading order.
What would settle it
Compute the next-order (N3LO) coefficient $c_3(k)$ at the same three relative momenta: if $|c_3|$ is not of natural size when $M_{\mathrm{hi}}$ is set near 100 MeV, or if $c_1$, $c_2$, and $c_3$ show strong correlations across orders, the iid-naturalness premise fails and the posterior mode near 100 MeV would not survive a correlation-aware analysis.
Extended reading notes
Core claim
The central claim is that Bayesian inference from the order-by-order neutron-deuteron total cross sections yields a posterior for the pionless effective field theory breakdown scale $M_{\mathrm{hi}}$ whose next-to-next-to-leading-order mode sits near 100 MeV, with a 68% highest-posterior-density interval of [63, 117] MeV, and that using the neutron-proton posterior as a prior in a sequential update produces a combined posterior with its mode very close to the pion mass, $m_\pi \approx 138$ MeV, at both next-to-leading and next-to-next-to-leading order. The authors read this as evidence that the inferred breakdown scale is not a single universal number but depends on the process probed, since the neutron-deuteron value is lower than the neutron-proton value near 190 MeV, and that the combined two- and three-nucleon evidence places the breakdown of pionless effective field theory close to the expected pion scale.
Load-bearing premise
The load-bearing premise is that the EFT expansion coefficients at the three chosen momenta are independent, identically distributed draws from a normal prior of natural size; if correlations across orders or momenta are significant, the factorized likelihood is invalid and the inferred breakdown-scale posterior is distorted.
Editorial extensions
If this is right
- If the inferred posterior is right, pionless effective field theory is quantitatively predictive for neutron-deuteron total cross sections up to momenta of order 100 MeV, and its breakdown in the three-nucleon sector is close to the expected pion scale.
- The breakdown scale is process-dependent: the neutron-deuteron value is lower than the neutron-proton value, so statements about the breakdown scale must specify the observable and kinematic regime under study.
- Combining two- and three-nucleon information through a sequential Bayesian update locates a common breakdown near $m_\pi \approx 138$ MeV, giving a concrete momentum at which explicit pion degrees of freedom become necessary.
- The next-to-next-to-leading-order posterior is robust to changes in the prior and in the selection of the three relative momenta, whereas the next-to-leading-order posterior is more sensitive to those choices.
- Including more than three momenta artificially narrows the posteriors because the independent-identically-distributed assumption for the coefficients is eventually violated by the finite correlation length of effective field theory truncation errors.
Reading between the lines
- A testable extension of the paper's own speculation: if the lower neutron-deuteron breakdown scale reflects off-shell probing of the two-body $t$-matrix near the pion cut at $m_\pi/2$, then repeating the inference with an explicitly pionful description should push the mode back above 138 MeV.
- The sequential update assumes conditional independence of the neutron-deuteron and neutron-proton likelihoods; a Gaussian-process treatment that models correlations across orders and momenta would likely widen the combined posterior and could move its mode away from exactly the pion mass.
- Applying the same point-wise Bayesian estimate to spin-dependent two- and three-nucleon observables, where the neglected next-to-next-to-leading-order $S$-$D$ operator matters, would provide a consistency check on whether the total-cross-section inference is biased by that omission.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper computes neutron-deuteron (nd) total cross sections in pionless effective field theory (EFT) order by order through next-to-next-to-leading order (N2LO), fitting the three-nucleon couplings to the triton binding energy and the nd scattering length. It then applies the pointwise Bayesian framework of Melendez et al. to infer a posterior for the EFT breakdown scale M_hi from the order-by-order convergence pattern, using the NLO and N2LO predictions at three selected momenta. The nd-only posterior at N2LO has a mode near 100 MeV with a 68% credible interval of [63,117] MeV. Combining this with a previous np-based posterior via a sequential Bayesian update that assumes conditional independence yields a combined posterior with a mode close to the pion mass (near 138 MeV). The paper concludes that the inferred breakdown scale depends on the process and that the combined analysis sharpens the inference.
Significance. If the nd-only inference is reliable, it provides a quantitative benchmark for the pionless EFT breakdown scale in the three-nucleon sector and a nontrivial cross-system consistency check. The paper includes useful robustness checks (prior choice, hyperprior variation, momentum selection) and reports explicit overlap coefficients between the nd and np posteriors. The nd-only result appears well supported by the displayed posteriors and the stated robustness tests. However, the combined np+nd claim is not supported at the same level of evidence: it depends on a conditional-independence assumption that the paper itself does not test, and the abstract presents the combined result without the caveat that the authors state in the final section. The combined claim should therefore be either substantiated with a correlation treatment or clearly downgraded in the presentation.
major comments (3)
- [Sec. III, Eq. (17)] The combined posterior p(Mhi|σnd,σnp,I) relies on the factorization p(σnd,σnp|Mhi,I) = p(σnd|Mhi,I)p(σnp|Mhi,I), i.e., conditional independence of the truncation errors for the two observables. This is not justified: both nd and np cross sections are computed from the same pionless EFT truncated at the same orders, sharing the same two-body LECs and the same naturalness prior. The paper's own statement in Sec. III that including more momenta 'will significantly violate the iid assumption' applies equally to adding a correlated observable. The reported overlap coefficient η=0.23 at N2LO shows that the nd and np posteriors barely overlap; multiplying such incompatible likelihoods produces the apparent 'sharpening' near mπ by construction, not as an independent confirmation. A concrete test is needed, for example a Gaussian-process model for the joint truncation error across the two observables, or at least a sensitivity study that drops the product and reports the nd-only and np-only results with an explicit measure of tension.
- [Abstract and Sec. IV] The abstract states that the combined analysis 'further sharpens the inference, placing the mode close to the pion mass scale' without any caveat, whereas Sec. IV concedes that this combined result 'might change when accounting for correlations.' This is a load-bearing overstatement because the combined result is the primary evidence for the 'close to the pion mass' claim. The abstract should either carry the same caveat or the combined result should be presented as a speculative consistency check rather than a sharpened inference.
- [Sec. III, paragraph after Fig. 5] The text notes that the combined posterior 'results from averaging two distributions on opposite sides of Mhi ≈ mπ.' This reveals that the combined mode near 138 MeV is an artifact of multiplying two non-overlapping posteriors, not a consequence of additional information. The sequential update using p(Mhi|σnp,I) as a prior is described as 'almost equivalent' to the product in Eq. (17), but the equivalence depends on the np prior being proportional to a likelihood, which is only approximately true and is not demonstrated in the manuscript. The robustness of the combined claim should be assessed against the possibility that the two datasets are partially correlated, which could substantially shift the mode and reduce the perceived sharpening.
minor comments (4)
- [Sec. III, Fig. 3 caption and text] The sentence 'the N2LO mode exhibits is at a greater value' contains a grammatical error; it should read 'the N2LO mode is at a greater value' or similar.
- [Sec. III] The word 'sligthly' should be 'slightly' in the paragraph on momentum selection.
- [Eq. (4) and surrounding text] The symbol '⊙' is described as elementwise multiplication, but its action on the vectors and matrices in Eq. (4) could be clarified, especially since the other product symbol '⊗' is defined by an integral in Eq. (3).
- [Sec. IV] The statement that the nd-inferred M_hi is 'slightly lower' than the np result is in tension with the reported 68% intervals not overlapping at N2LO; consider wording such as 'noticeably lower' or explicitly noting the non-overlap in the summary.
Circularity Check
No significant circularity: the M_hi posterior is inferred from the pionless EFT's own order-by-order convergence pattern, and the self-cited np posterior is independent two-body evidence.
full rationale
The derivation of p(M_hi|sigma_nd,I) is self-contained. Equations (14)-(15) define a likelihood on the order-by-order EFT cross-section predictions through the expansion coefficients c_i(k) with Q(k)=f(k)/M_hi, and the observed coefficients are extracted from the actual LO/NLO/N2LO nd calculations. The posterior is therefore a genuine Bayesian update on the model's convergence pattern rather than a rearrangement of the prior. The np posterior from Ref. [7], used as a prior in the combined analysis, is independent two-body evidence: it is conditioned on np scattering predictions, not on any nd quantity fitted in this paper, and the sequential update is explicitly performed 'assuming independence' with the caveat that the combined result 'might change when accounting for correlations.' An unverified conditional-independence assumption is a statistical robustness concern, not a case of a claim reducing to its inputs by construction. Likewise, the statement that including many momenta 'will significantly violate the iid assumption' is an openly stated limitation of the likelihood model, not a concealed circular step. No quoted equation reduces to a fitted parameter renamed as a prediction: the three-nucleon LECs are fitted to the triton binding energy and nd scattering length, but M_hi is not one of those fitted constants; it is a hyperparameter inferred from the naturalness of the remaining order-by-order coefficients.
Assumptions & free parameters
free parameters (3)
- M_hi (breakdown scale) =
mode ~100 MeV (nd), ~138 MeV (combined)
- bar c^2 (EFT coefficient variance) =
Integrated out with inverse-chi-square prior nu0=2, tau0^2=1
- H_LO, H_NLO, H_N2LO (three-nucleon couplings) =
Not reported numerically; calibrated to Bt=8.48 MeV and ad=0.65 plus or minus 0.04 fm
assumptions (6)
- domain assumption The pionless EFT expansion sigma^(nu) = sigma_ref sum_{i=0}^{nu} c_i Q^i with Q=f(k)/M_hi correctly describes the order-by-order cross sections.
- domain assumption The EFT coefficients c_i(k) are iid draws from a normal prior with natural variance.
- domain assumption The np and nd predictions are conditionally independent given M_hi, enabling the sequential Bayesian update.
- domain assumption Neglect of the N2LO S-D coupling operator does not materially change the total nd cross section.
- ad hoc to paper The soft scale f(k) is represented by a sixth-order soft-max function smoothing over the deuteron binding momentum.
- domain assumption Log-uniform prior over M_hi in (m_pi/40, 40 m_pi), with robustness to a uniform prior checked.
Cite this review
Pith. "Pith review of The breakdown scale of pionless effective field theory in the three-nucleon sector." pith.science (2026). https://pith.science/paper/L3NCFMUK
@misc{pith2026250708700,
author = {Pith},
title = {Pith review of: The breakdown scale of pionless effective field theory in the three-nucleon sector},
year = {2026},
howpublished = {\url{https://pith.science/paper/L3NCFMUK}},
note = {Machine review of arXiv:2507.08700}
}
read the original abstract
We make order-by-order predictions of neutron-deuteron total cross sections up to next-to-next-to-leading order in pionless effective field theory. Using Bayesian methods, we infer a posterior distribution for the breakdown scale. The result shows a mode near 100 MeV, and a combined analysis with neutron-proton scattering further sharpens the inference, placing the mode close to the pion mass scale, consistent with the expected range of pionless effective field theory.
Figures
Reference graph
Works this paper leans on
-
[7]
A. Ekstr¨ om and L. Platter, Phys. Lett. B 860, 139207 (2025), arXiv:2409.08197 [nucl-th]
arXiv 2025
-
[1]
E. Epelbaum, H.-W. Hammer, and U.-G. Meissner, Rev. Mod. Phys. 81, 1773 (2009), arXiv:0811.1338 [nucl-th]
arXiv 2009
-
[2]
R. Machleidt and D. R. Entem, Phys. Rept. 503, 1 (2011), arXiv:1105.2919 [nucl-th]
arXiv 2011
-
[3]
H. W. Hammer, S. K¨ onig, and U. van Kolck, Rev. Mod. Phys. 92, 025004 (2020), arXiv:1906.12122 [nucl-th]
arXiv 2020
-
[4]
van Kolck, Nuclear Physics A 645, 273 (1999)
U. van Kolck, Nuclear Physics A 645, 273 (1999)
work page 1999
-
[5]
D. B. Kaplan, M. J. Savage, and M. B. Wise, Phys. Lett. B 424, 390 (1998), arXiv:nucl-th/9801034
arXiv 1998
-
[6]
D. B. Kaplan, M. J. Savage, and M. B. Wise, Nuclear Physics B 534, 329 (1998)
work page 1998
-
[8]
J. M. Bub, M. Piarulli, R. J. Furnstahl, S. Pastore, and D. R. Phillips, (2024), arXiv:2408.02480 [nucl-th]
arXiv 2024
Show all 21 references
-
[9]
P. F. Bedaque, H. W. Hammer, and U. van Kolck, Phys. Rev. Lett. 82, 463 (1999), arXiv:nucl-th/9809025
1999 arXiv
-
[10]
P. F. Bedaque, H. W. Hammer, and U. van Kolck, Nucl. Phys. A 676, 357 (2000), arXiv:nucl-th/9906032
2000 arXiv
-
[11]
Efimov, Phys
V. Efimov, Phys. Rev. C 44, 2303 (1991)
1991
-
[12]
Margaryan, R
A. Margaryan, R. P. Springer, and J. Vanasse, Phys. Rev. C 93, 054001 (2016), arXiv:1512.03774 [nucl-th]. 6
2016 arXiv
-
[13]
P. F. Bedaque, G. Rupak, H. W. Griesshammer, and H.- W. Hammer, Nucl. Phys. A 714, 589 (2003), arXiv:nucl- th/0207034
2003
-
[14]
J. H. Hetherington and L. H. Schick, Phys. Rev. 137, B935 (1965)
1965
- [15]
-
[16]
W. Dilg, L. Koester, and W. Nistler, Phys. Lett. B 36, 208 (1971)
1971
-
[17]
J. A. Melendez, R. J. Furnstahl, D. R. Phillips, M. T. Pratola, and S. Wesolowski, Phys. Rev. C 100, 044001 (2019), arXiv:1904.10581 [nucl-th]
2019 arXiv
-
[18]
R. J. Furnstahl, N. Klco, D. R. Phillips, and S. Wesolowski, Phys. Rev. C 92, 024005 (2015), arXiv:1506.01343 [nucl-th]
2015 arXiv
-
[19]
J. A. Melendez, S. Wesolowski, and R. J. Furnstahl, Phys. Rev. C 96, 024003 (2017), arXiv:1704.03308 [nucl- th]
2017 arXiv
-
[20]
Svensson, A
I. Svensson, A. Ekstr¨ om, and C. Forss´ en, Phys. Rev. C 109, 064003 (2024), arXiv:2304.02004 [nucl-th]
2024 arXiv
-
[21]
P. J. Millican, R. J. Furnstahl, J. A. Melendez, D. R. Phillips, and M. T. Pratola, Phys. Rev. C 110, 044002 (2024), arXiv:2402.13165 [nucl-th]
2024 arXiv
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.