Recognition: 2 theorem links
· Lean TheoremChiral interactions and superfluidity in the calcium isotopic chain
Pith reviewed 2026-05-16 09:41 UTC · model grok-4.3
The pith
Bayesian calibration of third-order chiral interactions leaves neutron pairing gaps in calcium isotopes too small.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We perform ab initio calculations of three-point mass differences in the odd- and even-mass 39-49Ca isotopes to probe nuclear superfluidity via empirical neutron pairing gaps. We also quantify the sensitivity of those gaps to the parameters of the interaction at mean-field level. Recent studies employing accurate chiral nuclear interactions have found these gaps to be too small. We show that experimental values can be reproduced at mean-field level by substantially increasing the attraction of the singlet S-wave two-nucleon contact interaction, but doing so induces an unphysical bound state of the di-neutron. The sensitivity of these predictions to the full calibration of the nuclear is the
What carries the argument
Bayesian posterior sampling in delta-full chiral effective field theory at third chiral order, applied to mean-field calculations of pairing gaps extracted from three-point mass differences.
Where Pith is reading between the lines
- Pairing in these nuclei likely depends on collective or correlation effects not captured at mean-field level with current forces.
- The result may constrain models of superfluidity in neutron-rich systems relevant to neutron-star crusts.
- Tests with explicit four-nucleon forces or coupled-cluster methods could show whether the gap can be recovered without the di-neutron artifact.
Load-bearing premise
Mean-field calculations combined with third-order chiral interactions and the chosen Bayesian sampling are sufficient to conclude that the pairing gap discrepancy originates outside the interaction calibration.
What would settle it
A calculation at higher chiral order or beyond mean-field that reproduces the experimental calcium pairing gaps without producing a bound di-neutron state.
Figures
read the original abstract
We perform ab initio calculations of three-point mass differences in the odd- and even-mass $^{39-49}$Ca isotopes to probe nuclear superfluidity via empirical neutron pairing gaps. We also quantify the sensitivity of those gaps to the parameters of the interaction at mean-field level. Recent studies employing accurate chiral nuclear interactions have found these gaps to be too small. We show that experimental values can be reproduced at mean-field level by substantially increasing the attraction of the singlet $S$-wave two-nucleon contact interaction, but doing so induces an unphysical bound state of the di-neutron. The sensitivity of these predictions to the full calibration of the nuclear interaction is then studied by performing Bayesian posterior sampling in a delta-full chiral effective field theory at third chiral order. We find that pairing gaps remain largely unaffected, leaving the explanation of nuclear superfluidity as a future task for improved many-body modeling and refined interactions at higher chiral orders.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript computes empirical neutron pairing gaps in the $^{39-49}$Ca chain via three-point mass differences using ab initio mean-field calculations with delta-full chiral EFT interactions at third order. Bayesian posterior sampling over the low-energy constants shows that the gaps remain largely insensitive to parameter variations within the posterior and are systematically too small relative to experiment. Reproducing the experimental gaps at this level requires unphysically strong singlet S-wave attraction that induces a bound di-neutron state. The authors conclude that the origin of nuclear superfluidity must be addressed by improved many-body methods or refined interactions at higher chiral orders.
Significance. If the central result holds, the work provides a clear demonstration that uncertainties in the N3LO chiral interaction calibration do not resolve the underprediction of pairing gaps at mean-field level. The combination of direct ab initio gap computation with Bayesian sampling isolates the discrepancy to the many-body truncation, offering a reproducible and falsifiable benchmark for future studies. This strengthens the case for prioritizing beyond-mean-field correlations or N4LO+ interactions in ab initio nuclear theory.
major comments (2)
- [Results on Bayesian sampling and mean-field calculations] The conclusion that interaction calibration is not the source of the discrepancy rests on the mean-field truncation (explicitly stated in the abstract and results section). While the paper defers resolution to improved many-body modeling, it does not quantify how the sensitivity to LECs would change under a beyond-mean-field treatment such as second-order perturbation theory or coupled-cluster with triples; this is load-bearing for the claim that the gaps 'remain largely unaffected'.
- [Section on three-point mass differences] In the section deriving the three-point mass differences, the formula is applied directly at the mean-field level without explicit inclusion of higher-order corrections to the binding energies. If these corrections are comparable in size to the pairing gaps themselves, they could alter the reported insensitivity; the manuscript should state whether such corrections are estimated to be negligible.
minor comments (2)
- [Figure showing pairing gaps] The caption of the figure displaying pairing gaps across the isotopic chain should explicitly note that the shaded bands represent the 68% credible interval from the Bayesian posterior samples.
- [Discussion of unphysical adjustments] A short statement clarifying the precise definition of the singlet S-wave contact strength (e.g., the value of C_{^1S0} that induces the di-neutron bound state) would improve reproducibility.
Simulated Author's Rebuttal
We thank the referee for the careful review and the recommendation of minor revision. We respond to the major comments point by point below.
read point-by-point responses
-
Referee: The conclusion that interaction calibration is not the source of the discrepancy rests on the mean-field truncation (explicitly stated in the abstract and results section). While the paper defers resolution to improved many-body modeling, it does not quantify how the sensitivity to LECs would change under a beyond-mean-field treatment such as second-order perturbation theory or coupled-cluster with triples; this is load-bearing for the claim that the gaps 'remain largely unaffected'.
Authors: We agree that our conclusions are drawn at the mean-field level. The manuscript demonstrates insensitivity of the pairing gaps to LEC variations specifically within the mean-field approximation using Bayesian sampling. We do not claim this insensitivity extends to beyond-mean-field methods; rather, the results highlight that the discrepancy persists even after accounting for interaction uncertainties at N3LO in mean-field. Quantifying the LEC sensitivity at higher many-body orders would require substantial additional computations (e.g., using CCSD(T)) and is beyond the scope of this work focused on mean-field. We have revised the text to more explicitly limit the scope of the 'largely unaffected' statement to the mean-field calculations performed. revision: partial
-
Referee: In the section deriving the three-point mass differences, the formula is applied directly at the mean-field level without explicit inclusion of higher-order corrections to the binding energies. If these corrections are comparable in size to the pairing gaps themselves, they could alter the reported insensitivity; the manuscript should state whether such corrections are estimated to be negligible.
Authors: The three-point mass differences are indeed computed from the mean-field binding energies without incorporating higher-order corrections. We have updated the manuscript to include a statement clarifying that higher-order many-body corrections to the binding energies are not included in the present mean-field study. We do not estimate their magnitude relative to the gaps in this work, as doing so would entail separate calculations. This omission is consistent with the mean-field focus of the paper, and any significant impact from such corrections would reinforce our call for improved many-body modeling. revision: yes
- Quantifying how the sensitivity of the pairing gaps to the LECs changes when going beyond the mean-field approximation, for example using second-order perturbation theory or coupled-cluster with triples.
Circularity Check
No significant circularity detected in derivation chain
full rationale
The paper computes three-point mass differences to extract neutron pairing gaps at mean-field level from third-order delta-full chiral interactions. It shows that reproducing experimental gaps requires unphysical strengthening of the singlet S-wave contact term that induces a di-neutron bound state. Bayesian sampling over the posterior (constrained by external calibration data) then demonstrates that the gaps remain largely insensitive across the sampled parameter space. This insensitivity is a direct numerical output of the sampling procedure, not a quantity defined in terms of itself or recovered by construction from the fit. No step matches self-definitional, fitted-input-called-prediction, or self-citation-load-bearing patterns; the central claim is independent of the inputs and defers resolution to future many-body methods or higher-order interactions.
Axiom & Free-Parameter Ledger
free parameters (1)
- low-energy constants of the chiral interaction
axioms (2)
- domain assumption Chiral effective field theory provides a systematic, order-by-order expansion of nuclear forces consistent with QCD symmetries.
- domain assumption Mean-field level treatment with the chosen many-body method captures the dominant contributions to three-point mass differences.
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We find that pairing gaps remain largely unaffected, leaving the explanation of nuclear superfluidity as a future task for improved many-body modeling and refined interactions at higher chiral orders.
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Bayesian posterior sampling in a delta-full chiral effective field theory at third chiral order
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
L. N. Cooper, R. L. Mills, and A. M. Sessler, Possible Su- perfluidity of a System of Strongly Interacting Fermions, Phys. Rev.114, 1377 (1959)
work page 1959
-
[2]
D. J. Dean and M. Hjorth-Jensen, Pairing in nuclear sys- tems: from neutron stars to finite nuclei, Rev. Mod. Phys. 75, 607 (2003)
work page 2003
- [3]
- [4]
- [5]
- [6]
- [7]
- [8]
-
[9]
E. Epelbaum, H.-W. Hammer, and U.-G. Meißner, Mod- ern Theory of Nuclear Forces, Rev. Mod. Phys.81, 1773 (2009)
work page 2009
-
[10]
R. Machleidt and D. R. Entem, Chiral effective field the- ory and nuclear forces, Phys. Rep.503, 1 (2011)
work page 2011
-
[11]
P. Ring and P. Schuck,The Nuclear Many-Body Problem (Springer-Verlag, New-York, 1980)
work page 1980
-
[12]
Schunck, ed.,Energy Density Functional Methods for Atomic Nuclei, 2053-2563 (IOP Publishing, 2019)
N. Schunck, ed.,Energy Density Functional Methods for Atomic Nuclei, 2053-2563 (IOP Publishing, 2019)
work page 2053
-
[13]
A. Ekstr¨ om, C. Forss´ en, G. Hagen, G. R. Jansen, W. Jiang, and T. Papenbrock, What is ab ini- tio in nuclear theory?, Frontiers in Physics11, 10.3389/fphy.2023.1129094 (2023)
-
[14]
Hergert, A Guided Tour ofab initioNuclear Many- Body Theory, Front
H. Hergert, A Guided Tour ofab initioNuclear Many- Body Theory, Front. in Phys.8, 379 (2020)
work page 2020
-
[15]
A. Signoracci, T. Duguet, G. Hagen, and G. R. Jansen, Ab initio bogoliubov coupled cluster theory for open-shell nuclei, Phys. Rev. C91, 064320 (2015)
work page 2015
- [16]
- [17]
-
[18]
W. G. Jiang, A. Ekstr¨ om, C. Forss´ en, G. Hagen, G. R. Jansen, and T. Papenbrock, Accurate bulk properties of nuclei froma= 2 to∞from potentials with ∆ isobars, Phys. Rev. C102, 054301 (2020)
work page 2020
-
[19]
M. Wang, W. Huang, F. Kondev, G. Audi, and S. Naimi, The ame 2020 atomic mass evaluation (ii). tables, graphs and references*, Chin. Phys. C45, 030003 (2021)
work page 2020
- [20]
-
[21]
K. Hebeler, S. K. Bogner, R. J. Furnstahl, A. Nogga, and A. Schwenk, Improved nuclear matter calculations from chiral low-momentum interactions, Phys. Rev. C 83, 031301 (2011)
work page 2011
- [22]
-
[23]
A. Scalesi,On the characterization of nuclear many-body correlations in the ab initio approach, Phd thesis, Uni- versit´ e Paris-Saclay (2024)
work page 2024
-
[24]
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
- [25]
-
[26]
Z. H. Sun, A. Ekstr¨ om, C. Forss´ en, G. Hagen, G. R. Jansen, and T. Papenbrock, Multiscale physics of atomic nuclei from first principles, Phys. Rev. X15, 011028 (2025)
work page 2025
- [27]
-
[28]
T. Miyagi, Nuhamil : A numerical code to gen- erate nuclear two- and three-body matrix elements from chiral effective field theory, Eur. Phys. J. A59, 10.1140/epja/s10050-023-01039-y (2023)
-
[29]
R. F. Garcia Ruiz, M. L. Bissell, K. Blaum, A. Ekstr¨ om, N. Fr¨ ommgen, G. Hagen, M. Ham- men, K. Hebeler, J. D. Holt, G. R. Jansen, M. Kowalska, K. Kreim, W. Nazarewicz, R. Neugart, G. Neyens, W. N¨ ortersh¨ auser, T. Papenbrock, J. Papuga, A. Schwenk, J. Simonis, K. A. Wendt, and D. T. Yor- danov, Unexpectedly large charge radii of neutron-rich calcium...
work page 2016
-
[30]
E. Caurier, K. Langanke, G. Martinez-Pinedo, F. Nowacki, and P. Vogel, Shell model description of isotope shifts in calcium, Phys. Lett. B522, 240 (2001)
work page 2001
- [31]
- [32]
-
[33]
A. Scalesi, T. Duguet, P. Demol, M. Frosini, V. Som` a, and A. Tichai, Impact of correlations on nuclear bind- ing energies: Ab initio calculations of singly and dou- bly open-shell nuclei, Eur. Phys. J. A60, 209 (2024), arXiv:2406.03545 [nucl-th]
- [34]
-
[35]
M. Kortelainen, Z. Sun, G. Hagen, W. Nazarewicz, T. Papenbrock, and P.-G. Reinhard, Universal trend of charge radii of even-even ca–zn nuclei, Phys. Rev. C105, L021303 (2022)
work page 2022
-
[36]
S. J. Novario, D. Lonardoni, S. Gandolfi, and G. Hagen, Trends of neutron skins and radii of mirror nuclei from first principles, Phys. Rev. Lett.130, 032501 (2023)
work page 2023
-
[37]
M. Frosini, T. Duguet, B. Bally, Y. Beaujeault-Taudi` ere, J. P. Ebran, and V. Som` a, In-mediumk-body reduction ofn-body operators: A flexible symmetry-conserving ap- proach based on the sole one-body density matrix, Eur. Phys. J. A57, 151 (2021), arXiv:2102.10120 [nucl-th]
- [38]
-
[39]
J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Theory of superconductivity, Phys. Rev.108, 1175 (1957)
work page 1957
- [40]
-
[41]
I. Svensson, A. Ekstr¨ om, and C. Forss´ en, Bayesian pa- rameter estimation in chiral effective field theory using the Hamiltonian Monte Carlo method, Phys. Rev. C105, 014004 (2022), arXiv:2110.04011 [nucl-th]
-
[42]
I. Svensson, A. Ekstr¨ om, and C. Forss´ en, Inference of the low-energy constants in ∆-full chiral effective field theory including a correlated truncation error, Phys. Rev. C109, 064003 (2024), arXiv:2304.02004 [nucl-th]
-
[43]
I. Svensson, A. Ekstr¨ om, and C. Forss´ en, Bayesian esti- mation of the low-energy constants up to fourth order in the nucleon-nucleon sector of chiral effective field theory, Phys. Rev. C107, 014001 (2023)
work page 2023
-
[44]
S. Wesolowski, I. Svensson, A. Ekstr¨ om, C. Forss´ en, R. J. Furnstahl, J. A. Melendez, and D. R. Phillips, Rigor- ous constraints on three-nucleon forces in chiral effec- tive field theory from fast and accurate calculations of few-body observables, Phys. Rev. C104, 064001 (2021), arXiv:2104.04441 [nucl-th]
-
[45]
B. Huet al., Ab initio predictions link the neutron skin of 208Pb to nuclear forces, Nature Phys.18, 1196 (2022), arXiv:2112.01125 [nucl-th]
- [46]
-
[47]
W. G. Jiang, C. Forss´ en, T. Dj¨ arv, and G. Hagen, Nuclear-matter saturation and symmetry energy within ∆-full chiral effective field theory, Phys. Rev. C109, L061302 (2024)
work page 2024
- [48]
- [49]
-
[50]
Kondoet al., First observation of 28O, Nature620, 965 (2023), [Erratum: Nature 623, E13 (2023)]
Y. Kondoet al., First observation of 28O, Nature620, 965 (2023), [Erratum: Nature 623, E13 (2023)]
work page 2023
-
[51]
A. F. M. Smith and A. E. Gelfand, Bayesian statistics without tears: A sampling-resampling perspective, Am. Stat.46, 84 (1992)
work page 1992
-
[52]
W. Jiang and C. Forss´ en, Bayesian probability updates using sampling/importance resampling: Applications in nuclear theory, Frontiers in PhysicsV olume 10 - 2022, 10.3389/fphy.2022.1058809 (2022)
-
[53]
T. Duguet, J.-P. Ebran, M. Frosini, H. Hergert, and V. Som` a, Rooting the EDF method into the ab initio framework: PGCM-PT formalism based on MR-IMSRG pre-processed Hamiltonians, The European Physical Journal A59, 13 (2023)
work page 2023
- [54]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.