pith. machine review for the scientific record. sign in

arxiv: 2605.12023 · v1 · submitted 2026-05-12 · ⚛️ nucl-th · hep-ph

Recognition: no theorem link

Mass radius and D-term of atomic nuclei in relativistic mean field theory

Authors on Pith no claims yet

Pith reviewed 2026-05-13 03:17 UTC · model grok-4.3

classification ⚛️ nucl-th hep-ph
keywords D-termnuclear radiirelativistic mean field theoryshell structureenergy-momentum tensormagic numbersgravitational form factorsatomic nuclei
0
0 comments X

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.

The paper applies relativistic mean field theory to calculate the mass radius and radii derived from the energy-momentum tensor, along with the D-term, for many spin-zero nuclei. It examines how these quantities change with neutron number in isotopic chains of calcium, nickel, zirconium, tin, and lead. The key finding is that the absolute value of the D-term does not grow steadily with more neutrons but instead peaks and dips at magic and sub-magic numbers. This pattern produces visible kinks in the various radii, highlighting how nuclear shell structure influences the mechanical properties of nuclei.

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

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

  • 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

Figures reproduced from arXiv: 2605.12023 by Makoto Oka, Tomohiro Oishi, Yoshitaka Hatta.

Figure 1
Figure 1. Figure 1: FIG. 1: D-terms of nuclear isotopes as a function of the neutron number [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Charge radius (top), mass radius (middle) and shear radius (bottom) of Ca, Ni, Zr, Sn and Pb isotopes. [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Differences [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Top: Pressure and shear distributions, [PITH_FULL_IMAGE:figures/full_fig_p011_4.png] view at source ↗
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.

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

3 major / 2 minor

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)
  1. [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.
  2. [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.
  3. [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)
  1. [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).
  2. [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

3 responses · 0 unresolved

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
  1. 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

  2. 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

  3. 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

0 steps flagged

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

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review; specific free parameters, axioms, and invented entities cannot be identified without the full manuscript methods and equations sections.

pith-pipeline@v0.9.0 · 5456 in / 987 out tokens · 73501 ms · 2026-05-13T03:17:24.472611+00:00 · methodology

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

77 extracted references · 77 canonical work pages · 1 internal anchor

  1. [1]

    In reality, they have their own electromagnetic form factors reflecting the internal distribution of quarks

    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. [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. [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. [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. [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. [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 ...

  7. [7]

    De Vries, C

    H. De Vries, C. W. De Jager, and C. De Vries, Atom. Data Nucl. Data Tabl.36, 495 (1987)

  8. [8]

    Angeli and K

    I. Angeli and K. P. Marinova, Atomic Data and Nuclear Data Tables99, 69 (2013)

  9. [9]

    Z. Sun, K. A. Beyer, Z. A. Mandrykina, I. A. Valuev, C. H. Keitel, and N. S. Oreshkina, Phys. Rev. Lett.135, 163002 (2025), 2504.19977

  10. [10]

    R. D. Woods and D. S. Saxon, Phys. Rev.95, 577 (1954)

  11. [11]

    I. Y. Kobzarev and L. B. Okun’, Zh. Eksp. Teor. Fiz.43, 1904 (1962)

  12. [12]

    Pagels, Phys

    H. Pagels, Phys. Rev.144, 1250 (1966)

  13. [13]

    M. V. Polyakov and P. Schweitzer, Int. J. Mod. Phys. A33, 1830025 (2018), 1805.06596

  14. [14]

    Goeke, J

    K. Goeke, J. Grabis, J. Ossmann, M. V. Polyakov, P. Schweitzer, A. Silva, and D. Urbano, Phys. Rev. D75, 094021 (2007), hep-ph/0702030

  15. [15]

    Cebulla, K

    C. Cebulla, K. Goeke, J. Ossmann, and P. Schweitzer, Nucl. Phys. A794, 87 (2007), hep-ph/0703025

  16. [16]

    J.-H. Jung, U. Yakhshiev, and H.-C. Kim, J. Phys. G41, 055107 (2014), 1310.8064

  17. [17]

    V. D. Burkert, L. Elouadrhiri, and F. X. Girod, Nature557, 396 (2018)

  18. [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

  19. [19]

    Lorc´ e, H

    C. Lorc´ e, H. Moutarde, and A. P. Trawi´ nski, Eur. Phys. J. C79, 89 (2019), 1810.09837

  20. [20]

    I. V. Anikin, Phys. Rev. D99, 094026 (2019), 1902.00094

  21. [21]

    D. E. Kharzeev, Phys. Rev. D104, 054015 (2021), 2102.00110

  22. [22]

    K. A. Mamo and I. Zahed, Phys. Rev. D106, 086004 (2022), 2204.08857

  23. [23]

    Fujita, Y

    M. Fujita, Y. Hatta, S. Sugimoto, and T. Ueda, PTEP2022, 093B06 (2022), 2206.06578

  24. [24]

    Won, J.-Y

    H.-Y. Won, J.-Y. Kim, and H.-C. Kim, Phys. Rev. D106, 114009 (2022), 2210.03320

  25. [25]

    X.-Y. Wang, C. Dong, and Q. Wang, Phys. Rev. D106, 056027 (2022), 2206.11644

  26. [26]

    Duran et al., Nature615, 813 (2023), 2207.05212

    B. Duran et al., Nature615, 813 (2023), 2207.05212. 13

  27. [27]

    Y. Guo, X. Ji, Y. Liu, and J. Yang, Phys. Rev. D108, 034003 (2023), 2305.06992

  28. [28]

    D. C. Hackett, D. A. Pefkou, and P. E. Shanahan, Phys. Rev. Lett.132, 251904 (2024), 2310.08484

  29. [29]

    B. Wang, F. He, G. Wang, T. Draper, J. Liang, K.-F. Liu, and Y.-B. Yang (χQCD), Phys. Rev. D109, 094504 (2024), 2401.05496

  30. [30]

    Cao, F.-K

    X.-H. Cao, F.-K. Guo, Q.-Z. Li, and D.-L. Yao, Nature Commun.16, 6979 (2025), 2411.13398

  31. [31]

    Goharipour, F

    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. [32]

    Broniowski and E

    W. Broniowski and E. Ruiz Arriola, Phys. Rev. D112, 054028 (2025), 2503.09297

  33. [33]

    Fujii, M

    D. Fujii, M. Kawaguchi, and M. Tanaka, Phys. Lett. B866, 139559 (2025), 2503.09686

  34. [34]

    Tanaka, D

    M. Tanaka, D. Fujii, and M. Kawaguchi, Phys. Rev. D112, 054048 (2025), 2507.21220

  35. [35]

    Stegeman and R

    R. Stegeman and R. Zwicky, JHEP03, 184 (2026), 2508.18537

  36. [36]

    Sugimoto and T

    S. Sugimoto and T. Tsukamoto (2025), 2503.19492

  37. [37]

    S. Nair, C. Mondal, S. Xu, X. Zhao, and J. P. Vary (BLFQ), Phys. Rev. D112, 114001 (2025), 2506.07554

  38. [38]

    A. Sain, P. Choudhary, B. Gurjar, C. Mondal, D. Chakrabarti, and A. Mukherjee, Phys. Rev. D111, 094011 (2025), 2503.12574

  39. [39]

    Hippel¨ ainen, N

    A. Hippel¨ ainen, N. Jokela, and M. J¨ arvinen (2026), 2601.14939

  40. [40]

    Fukushima and T

    K. Fukushima and T. Uji (2026), 2603.11704

  41. [41]

    Freese and W

    A. Freese and W. Cosyn, Phys. Rev. D106, 114013 (2022), 2207.10787

  42. [42]

    He and I

    F. He and I. Zahed, Phys. Rev. C109, 045209 (2024), 2310.12315

  43. [43]

    He and I

    F. He and I. Zahed, Phys. Rev. C110, 014312 (2024), 2401.09318

  44. [44]

    Cosyn, A

    W. Cosyn, A. Freese, and A. Sosa (2026), 2602.18298

  45. [45]

    M. V. Polyakov, Phys. Lett. B555, 57 (2003), hep-ph/0210165

  46. [46]

    Guzey and M

    V. Guzey and M. Siddikov, J. Phys. G32, 251 (2006), hep-ph/0509158

  47. [47]

    Liuti and S

    S. Liuti and S. K. Taneja, Phys. Rev. C72, 034902 (2005), hep-ph/0504027

  48. [48]

    Garcia Martin-Caro, M

    A. Garcia Martin-Caro, M. Huidobro, and Y. Hatta, Phys. Rev. D108, 034014 (2023), 2304.05994

  49. [49]

    Garc´ ıa Mart´ ın-Caro, M

    A. Garc´ ıa Mart´ ın-Caro, M. Huidobro, and Y. Hatta, Phys. Rev. D110, 034002 (2024), 2312.12984

  50. [50]

    T. H. R. Skyrme, Proc. Roy. Soc. Lond. A260, 127 (1961)

  51. [51]

    Braaten and L

    E. Braaten and L. Carson, Phys. Rev. D38, 3525 (1988)

  52. [52]

    Carson, Nucl

    L. Carson, Nucl. Phys. A535, 479 (1991)

  53. [53]

    O. V. Manko, N. S. Manton, and S. W. Wood, Phys. Rev. C76, 055203 (2007), 0707.0868

  54. [54]

    D. T. J. Feist, P. H. C. Lau, and N. S. Manton, Phys. Rev. D87, 085034 (2013), 1210.1712

  55. [55]

    S. B. Gudnason and C. Halcrow, JHEP08, 117 (2022), 2202.01792

  56. [56]

    Halcrow, Skyrmions3D,https://github.com/chrishalcrow/Skyrmions3D.jl(2023)

    C. Halcrow, Skyrmions3D,https://github.com/chrishalcrow/Skyrmions3D.jl(2023)

  57. [57]

    A. X. L. Freire, C. Adam, A. Garc´ ıa Mart´ ın-Caro, and D. G. D´ ıaz (2025), 2510.25841

  58. [58]

    S. B. Gudnason and P. Leask (2026), 2604.25876

  59. [59]

    Garcıa Martın-Caro and C

    A. Garcıa Martın-Caro and C. Halcrow (2023), 2312.04335

  60. [60]

    J. D. Walecka, Annals Phys.83, 491 (1974)

  61. [61]

    B. D. Serot and J. D. Walecka, Phys. Lett. B87, 172 (1979)

  62. [62]

    Bender, P.-H

    M. Bender, P.-H. Heenen, and P.-G. Reinhard, Rev. Mod. Phys.75, 121 (2003)

  63. [63]

    Niksic, D

    T. Niksic, D. Vretenar, and P. Ring, Prog. Part. Nucl. Phys.66, 519 (2011), 1102.4193

  64. [64]

    Niksic, N

    T. Niksic, N. Paar, D. Vretenar, and P. Ring, Comput. Phys. Commun.185, 1808 (2014), 1403.4039

  65. [65]

    Brockmann and H

    R. Brockmann and H. Toki, Phys. Rev. Lett.68, 3408 (1992)

  66. [66]

    Ji, Front

    X. Ji, Front. Phys. (Beijing)16, 64601 (2021), 2102.07830

  67. [67]

    J. F. Donoghue, B. R. Holstein, B. Garbrecht, and T. Konstandin, Phys. Lett. B529, 132 (2002), [Erratum: Phys.Lett.B 612, 311–312 (2005)], hep-th/0112237

  68. [68]

    Varma and P

    M. Varma and P. Schweitzer, Phys. Rev. D102, 014047 (2020), 2006.06602

  69. [69]

    G. A. Lalazissis, T. Niksic, D. Vretenar, and P. Ring, Phys. Rev. C71, 024312 (2005)

  70. [70]

    U. C. Perera, A. V. Afanasjev, and P. Ring, Phys. Rev. C104, 064313 (2021)

  71. [71]

    Naito, T

    T. Naito, T. Oishi, H. Sagawa, and Z. Wang, Phys. Rev. C107, 054307 (2023), 2209.02857

  72. [72]

    K¨ onig et al., Phys

    K. K¨ onig et al., Phys. Rev. Lett.131, 102501 (2023), 2309.02839

  73. [73]

    Taniuchi et al., Nature569, 53 (2019), 1912.05978

    R. Taniuchi et al., Nature569, 53 (2019), 1912.05978

  74. [74]

    Hagen and T

    G. Hagen and T. Papenbrock, Nature569, 49 (2019)

  75. [75]
  76. [76]

    Hagiwara, X.-B

    Y. Hagiwara, X.-B. Tong, and B.-W. Xiao, Phys. Rev. D111, L051503 (2025), 2401.12840

  77. [77]

    He and I

    F. He and I. Zahed, Phys. Rev. C112, 014303 (2025), 2407.09991