Recognition: no theorem link
Mass radius and D-term of atomic nuclei in relativistic mean field theory
Pith reviewed 2026-05-13 03:17 UTC · model grok-4.3
The pith
The D-term of atomic nuclei shows local maxima and minima at magic neutron numbers, leading to kinks in their mechanical radii.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Based on relativistic mean field theory, computations of the D-term D(t=0) and related radii for dozens of spin-0 nuclei reveal that |D| exhibits local maxima and minima at magic and sub-magic neutron numbers rather than increasing monotonically with N. This behavior induces characteristic kinks in the mass, scalar, tensor, and shear radii for isotopes across the nuclear chart, demonstrating the strong sensitivity of these mechanical properties to nuclear shell structure.
What carries the argument
The D-term as the forward limit of the gravitational form factor D(t=0), extracted from the energy-momentum tensor distributions computed in relativistic mean field theory.
If this is right
- The D-term of nuclei varies non-monotonically with neutron number N, showing extrema at closed shells.
- Mass, scalar, tensor, and shear radii of nuclei display kinks corresponding to magic and sub-magic numbers.
- Mechanical properties of nuclei are strongly influenced by shell structure effects.
- These features appear consistently in multiple isotopic series including Ca, Ni, Zr, Sn, and Pb.
Where Pith is reading between the lines
- If the shell structure sensitivity holds, it may affect predictions for how nuclei respond in gravitational or high-energy scattering experiments.
- Similar computations could be extended to deformed or odd nuclei to check the robustness of the kinks.
- Experimental extraction of D-term via deeply virtual Compton scattering or other processes could test these predictions directly.
Load-bearing premise
The chosen parametrization of relativistic mean field theory faithfully reproduces the energy-momentum tensor inside the nuclei without introducing large uncontrolled model dependence.
What would settle it
Measurement of the D-term for a sequence of isotopes showing strictly monotonic increase with neutron number, without kinks at known magic numbers, would falsify the predicted sensitivity to shell structure.
Figures
read the original abstract
Based on relativistic mean field theory for atomic nuclei, we compute the mass radius and other radii associated with the energy momentum tensor for dozens of spin-0 nuclei across the nuclear chart. We also compute the D-term of these nuclei, the forward limit of the gravitational form factor $D(t=0)=D$. The dependence on the neutron number $N$ is systematically studied for calcium (Ca), nickel (Ni), zirconium (Zr), tin (Sn) and lead (Pb) isotopes. Remarkably, $|D|$ does not monotonically increase with $N$. Instead, it exhibits local maxima and minima when $N$ equals a magic number and even a sub-magic number. This results in characteristic kinks in the mass, scalar, tensor and shear radii of these isotopes. Our work for the first time elucidates the strong sensitivity of the various mechanical properties of nuclei to the nuclear shell structure.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper computes the mass radius and radii associated with the energy-momentum tensor (scalar, tensor, shear) as well as the D-term D(t=0) for dozens of spin-0 nuclei using relativistic mean field theory. It examines the neutron-number dependence along the Ca, Ni, Zr, Sn and Pb isotopic chains and reports that |D| is non-monotonic, exhibiting local maxima and minima at magic and sub-magic neutron numbers; these features produce characteristic kinks in the various radii. The central claim is that this constitutes the first demonstration of strong sensitivity of nuclear mechanical properties to shell structure.
Significance. If the reported shell-structure dependence survives scrutiny, the work would usefully extend the study of gravitational form factors and mechanical properties from nucleons to finite nuclei, potentially linking nuclear-structure observables to the energy-momentum tensor. The computational approach is standard RMFT, but the absence of cross-validation with alternative parametrizations and the lack of quantitative error estimates or comparisons to existing data reduce the immediate impact.
major comments (3)
- [Results for isotopic chains (Ca, Ni, Zr, Sn, Pb)] The central claim that the kinks in |D| and the radii demonstrate sensitivity to nuclear shell structure rests on results obtained with a single (or narrow set of) RMF parametrization(s). No comparison is shown with other standard RMF Lagrangians (e.g., NL3, PK1, or density-dependent variants) whose single-particle spectra and shell closures differ; if the extrema shift or disappear under a different parametrization, the observed non-monotonicity would be an artifact of the chosen effective interaction rather than a robust consequence of shell structure. This issue is load-bearing for the strongest claim in the abstract.
- [Abstract and computational details] The abstract states that computations were performed for dozens of nuclei, yet no validation against known experimental mass radii, charge radii, or other D-term-related quantities is provided, nor are error estimates or sensitivity analyses to RMF parameters reported. Without these, it is impossible to assess whether the reported kinks exceed the model uncertainty.
- [Formalism and D-term definition] The extraction of the D-term from the forward limit of the gravitational form factor in the RMF framework is not cross-checked against independent calculations (e.g., non-relativistic Skyrme or ab-initio methods) for even a single nucleus; such a benchmark would be required to establish that the shell-structure features are not an artifact of the mean-field approximation.
minor comments (2)
- [Abstract] The abstract and introduction should explicitly name the RMF parametrization(s) employed and state the range of nuclei studied (e.g., which specific isotopes).
- [Results] Notation for the various radii (mass, scalar, tensor, shear) should be defined once in the text and used consistently in all figures and tables.
Simulated Author's Rebuttal
We thank the referee for the careful and constructive review of our manuscript. We address each major comment below and have made revisions to the manuscript to incorporate additional analyses and validations as suggested.
read point-by-point responses
-
Referee: The central claim that the kinks in |D| and the radii demonstrate sensitivity to nuclear shell structure rests on results obtained with a single (or narrow set of) RMF parametrization(s). No comparison is shown with other standard RMF Lagrangians (e.g., NL3, PK1, or density-dependent variants) whose single-particle spectra and shell closures differ; if the extrema shift or disappear under a different parametrization, the observed non-monotonicity would be an artifact of the chosen effective interaction rather than a robust consequence of shell structure. This issue is load-bearing for the strongest claim in the abstract.
Authors: We agree that robustness across parametrizations strengthens the claim. In the revised manuscript, we have included calculations using the NL3 parametrization for the Ca, Ni, Zr, Sn, and Pb chains. The characteristic kinks at magic neutron numbers remain present, with the positions of extrema largely unchanged. This indicates that the sensitivity to shell structure is a general feature within the RMF framework. A new figure comparing the two parametrizations has been added, along with a discussion of the minor quantitative differences. revision: yes
-
Referee: The abstract states that computations were performed for dozens of nuclei, yet no validation against known experimental mass radii, charge radii, or other D-term-related quantities is provided, nor are error estimates or sensitivity analyses to RMF parameters reported. Without these, it is impossible to assess whether the reported kinks exceed the model uncertainty.
Authors: We have revised the manuscript to include direct comparisons of the computed mass radii and charge radii with experimental data for the studied isotopic chains. The model reproduces the overall trends and magic number effects well. We have also performed a sensitivity analysis by varying the RMF parameters within their typical ranges and included error estimates in the plots of |D| and the radii. These show that the kinks are robust against small parameter changes and exceed the estimated uncertainties. A new section on model validation has been added. revision: yes
-
Referee: The extraction of the D-term from the forward limit of the gravitational form factor in the RMF framework is not cross-checked against independent calculations (e.g., non-relativistic Skyrme or ab-initio methods) for even a single nucleus; such a benchmark would be required to establish that the shell-structure features are not an artifact of the mean-field approximation.
Authors: Cross-validation with other methods is desirable but challenging for the full range of nuclei considered. We have added a benchmark comparison for the D-term of 16O using our RMF approach against published Skyrme model results. The shell-structure sensitivity is qualitatively similar. For heavier nuclei, we discuss the limitations of ab-initio methods and argue that the RMF provides a consistent mean-field description across the chart. This discussion is included in the revised formalism section. revision: partial
Circularity Check
No significant circularity in RMFT-based derivation of nuclear D-terms and radii
full rationale
The paper performs direct numerical computations of energy-momentum tensor distributions, mass radii, and the D-term D(t=0) within standard relativistic mean-field theory for selected isotopes. The reported non-monotonic behavior and kinks at magic/sub-magic neutron numbers are outputs of these calculations rather than inputs used to define or fit any parameter. No self-definitional relations, fitted quantities renamed as predictions, or load-bearing self-citations appear in the provided abstract or derivation description; the central claims remain independent of the target shell-structure sensitivities.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
Electromagnetic charge radius In RMF, protons and neutrons are treated as elementary (pointlike) Dirac particles. In reality, they have their own electromagnetic form factors reflecting the internal distribution of quarks. The charge radiusr em of a nucleus is usually obtained by convoluting the proton and neutron densities with their respective electroma...
-
[2]
It is therefore a more reasonable measure of the nuclear size, especially for neutron-rich nuclei
Baryon number radius ⟨r2⟩B = 1 A Z d3r r2ρB(r).(24) In contrast to the charge radius, protons and neutrons equally contribute to⟨r 2⟩B. It is therefore a more reasonable measure of the nuclear size, especially for neutron-rich nuclei. In principle, again convolutions with the nucleon form factors are needed. But we neglect them in⟨r 2⟩B and in the other r...
-
[3]
We will use the Belinfante form which is more appropriate in the context of GFFs
Mass radius ⟨r2⟩m = 1 M Z d3r r2T 00(r).(25) While the total massM= R d3rT 00(r) is the same for both the canonical and Belinfante-improved energy density, with the weight factorr 2 the two energy densities give different results. We will use the Belinfante form which is more appropriate in the context of GFFs. 6
-
[4]
Using D, we can express the scalar radius as [8] ⟨r2⟩s =⟨r 2⟩m − 3D M 2 .(27)
Scalar radius ⟨r2⟩s = 1 M Z d3r r2T µ µ (r).(26) In the present model,T µ µ is given by (14) and represents the spin-0 part of the energy momentum tensor. Using D, we can express the scalar radius as [8] ⟨r2⟩s =⟨r 2⟩m − 3D M 2 .(27)
-
[5]
Our definition is different from the one in [60]
Tensor radius ⟨r2⟩t = 1 M Z d3r r2 T 00(r) + 1 2 Tii(r) = 3 2 ⟨r2⟩m − 1 2 ⟨r2⟩s.(28) The linear combination represents the irreducible spin-2 part of the energy momentum tensor [43]. Our definition is different from the one in [60]. Note that⟨r 2⟩t is not independent as it is given by [43] ⟨r2⟩t =⟨r 2⟩m + 3D 2M 2 .(29)
-
[6]
We now come to an important issue of the Coulomb field
Shear radius ⟨r2⟩shear = R d3r r2s(r)R d3r s(r) =− 15 4M DR d3r s(r) .(30) This is analogous to the ‘mechanical radius’ associated with the distribution of the ‘normal force’ 2 3 s(r) +p(r) [7]. We now come to an important issue of the Coulomb field. In fact, if we evaluate⟨r 2⟩m,⟨r 2⟩t,⟨r 2⟩shear and the D-term using the formulas derived in the previous ...
work page 1903
-
[7]
H. De Vries, C. W. De Jager, and C. De Vries, Atom. Data Nucl. Data Tabl.36, 495 (1987)
work page 1987
-
[8]
I. Angeli and K. P. Marinova, Atomic Data and Nuclear Data Tables99, 69 (2013)
work page 2013
- [9]
-
[10]
R. D. Woods and D. S. Saxon, Phys. Rev.95, 577 (1954)
work page 1954
-
[11]
I. Y. Kobzarev and L. B. Okun’, Zh. Eksp. Teor. Fiz.43, 1904 (1962)
work page 1904
- [12]
- [13]
- [14]
-
[15]
C. Cebulla, K. Goeke, J. Ossmann, and P. Schweitzer, Nucl. Phys. A794, 87 (2007), hep-ph/0703025
- [16]
-
[17]
V. D. Burkert, L. Elouadrhiri, and F. X. Girod, Nature557, 396 (2018)
work page 2018
-
[18]
Holographic $J/\psi$ production near threshold and the proton mass problem
Y. Hatta and D.-L. Yang, Phys. Rev. D98, 074003 (2018), 1808.02163
work page Pith review arXiv 2018
-
[19]
C. Lorc´ e, H. Moutarde, and A. P. Trawi´ nski, Eur. Phys. J. C79, 89 (2019), 1810.09837
- [20]
- [21]
- [22]
- [23]
- [24]
- [25]
-
[26]
Duran et al., Nature615, 813 (2023), 2207.05212
B. Duran et al., Nature615, 813 (2023), 2207.05212. 13
- [27]
- [28]
- [29]
- [30]
-
[31]
M. Goharipour, F. Irani, M. H. Amiri, H. Fatehi, B. Falahi, A. Moradi, and K. Azizi (MMGPDs), Nucl. Phys. B1017, 116962 (2025), 2503.08847
-
[32]
W. Broniowski and E. Ruiz Arriola, Phys. Rev. D112, 054028 (2025), 2503.09297
- [33]
- [34]
- [35]
- [36]
- [37]
- [38]
-
[39]
A. Hippel¨ ainen, N. Jokela, and M. J¨ arvinen (2026), 2601.14939
- [40]
-
[41]
A. Freese and W. Cosyn, Phys. Rev. D106, 114013 (2022), 2207.10787
- [42]
- [43]
- [44]
-
[45]
M. V. Polyakov, Phys. Lett. B555, 57 (2003), hep-ph/0210165
work page Pith review arXiv 2003
- [46]
-
[47]
S. Liuti and S. K. Taneja, Phys. Rev. C72, 034902 (2005), hep-ph/0504027
-
[48]
A. Garcia Martin-Caro, M. Huidobro, and Y. Hatta, Phys. Rev. D108, 034014 (2023), 2304.05994
-
[49]
A. Garc´ ıa Mart´ ın-Caro, M. Huidobro, and Y. Hatta, Phys. Rev. D110, 034002 (2024), 2312.12984
-
[50]
T. H. R. Skyrme, Proc. Roy. Soc. Lond. A260, 127 (1961)
work page 1961
- [51]
- [52]
- [53]
- [54]
- [55]
-
[56]
Halcrow, Skyrmions3D,https://github.com/chrishalcrow/Skyrmions3D.jl(2023)
C. Halcrow, Skyrmions3D,https://github.com/chrishalcrow/Skyrmions3D.jl(2023)
work page 2023
- [57]
-
[58]
S. B. Gudnason and P. Leask (2026), 2604.25876
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[59]
A. Garcıa Martın-Caro and C. Halcrow (2023), 2312.04335
-
[60]
J. D. Walecka, Annals Phys.83, 491 (1974)
work page 1974
-
[61]
B. D. Serot and J. D. Walecka, Phys. Lett. B87, 172 (1979)
work page 1979
-
[62]
M. Bender, P.-H. Heenen, and P.-G. Reinhard, Rev. Mod. Phys.75, 121 (2003)
work page 2003
- [63]
- [64]
- [65]
- [66]
- [67]
-
[68]
M. Varma and P. Schweitzer, Phys. Rev. D102, 014047 (2020), 2006.06602
-
[69]
G. A. Lalazissis, T. Niksic, D. Vretenar, and P. Ring, Phys. Rev. C71, 024312 (2005)
work page 2005
-
[70]
U. C. Perera, A. V. Afanasjev, and P. Ring, Phys. Rev. C104, 064313 (2021)
work page 2021
- [71]
-
[72]
K. K¨ onig et al., Phys. Rev. Lett.131, 102501 (2023), 2309.02839
-
[73]
Taniuchi et al., Nature569, 53 (2019), 1912.05978
R. Taniuchi et al., Nature569, 53 (2019), 1912.05978
- [74]
-
[75]
Y. Hatta, Phys. Rev. D109, L051502 (2024), 2311.14470
-
[76]
Y. Hagiwara, X.-B. Tong, and B.-W. Xiao, Phys. Rev. D111, L051503 (2025), 2401.12840
- [77]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.