REVIEW 3 major objections 4 minor 10 references
Neutron skins of heavy nuclei and tidal deformability of neutron star
T0 review · 3 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read This paper argues that two relativistic mean-field parameter sets fitted only to laboratory nuclei also reproduce neutron-skin measurements and the GW170817 tidal-deformability bound, letting a single energy functional connect nuclear…
desk verdict Useful proceedings summary of the author's own G3/IOPB-I results, but the GW170817 comparison overstates the agreement: two of three models exceed the quoted 90% upper limit. 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 load-bearing object is the extended relativistic mean-field (ERMF) energy-density functional: a Lagrangian in which nucleons couple to $\sigma$, $\omega$, $\rho$, and $\delta$ mesons and the photon, with nonlinear self-interactions and an $\omega$-$\rho$ cross-coupling. The paper uses two parameter sets, G3 and IOPB-I, fitted by simulated annealing to eight spherical nuclei together with nuclear-matter saturation properties. The $\omega$-$\rho$ cross-coupling controls the density dependence of the symmetry energy and therefore ties the neutron-skin thickness to the isovector part of the equation of state; that same equation of state, fed into the general-relativistic stellar-structure equations with a small metric perturbation, gives the dimensionless tidal deformability $\Lambda = 2k_2/(3C^5)$. One set of coupling constants thus determines both the skin and the tidal response.
What would settle it
A future $^{208}$Pb neutron-skin measurement with uncertainty below 0.02 fm whose central value falls outside the 0.16-0.23 fm range spanned by G3 and IOPB-I would break the claimed compatibility with the $\Delta r_{np}\le 0.25$ fm limit. A gravitational-wave bound placing $\Lambda_{1.4}$ below about 460 would likewise exclude the G3 equation of state, and a deformed-plus-pairing calculation for $^{238}$U that moves its skin by more than a few hundredths of a femtometer would test the spherical approximation.
Extended reading notes
Core claim
The central claim is that the calibrated energy-density functionals G3 and IOPB-I, built from an extended relativistic mean-field Lagrangian with $\sigma$, $\omega$, $\rho$, and $\delta$ mesons plus nonlinear couplings, reproduce the measured neutron-skin thicknesses of 26 stable nuclei from $^{40}$Ca to $^{238}$U. For $^{208}$Pb the models give $\Delta r_{np}=0.180$ fm (G3) and $0.221$ fm (IOPB-I), inside the PREX-II uncertainty band and below the upper limit $\Delta r_{np}\le 0.25$ fm obtained from the skin-tidal correlation. The paper further claims that the resulting tidal deformabilities of a 1.4-solar-mass neutron star, $\Lambda_{1.4}=461$ (G3), $622$ (FSUGarnet), and $681$ (IOPB-I), fall within the GW170817 90% credible bounds, while the stiffer NL3 equation of state is excluded.
Load-bearing premise
The comparison rests on the reliability of the G3 and IOPB-I calibrations—fits to eight spherical nuclei plus nuclear-matter saturation properties—and on the assumption that ignoring deformation and superfluidity for all 26 nuclei, including $^{238}$U, does not materially shift the predicted neutron-skin thicknesses.
Editorial extensions
If this is right
- A single energy-density functional, without astrophysical tuning, can account for the antiproton neutron-skin data across the nuclear chart and for the GW170817 tidal bound.
- The $\Delta r_{np}\le 0.25$ fm limit for $^{208}$Pb is satisfied by G3 (0.180 fm) and IOPB-I (0.221 fm), so the PREX-II result does not force an unusually stiff symmetry energy.
- The stiffer NL3 force is excluded by GW170817, while G3 gives the lowest $\Lambda_{1.4}$ (461) of the models considered, making it the most constrained by future merger observations.
- The reproduction of the nearly linear skin-versus-asymmetry band means the isovector part of the functional—chiefly the $\omega$-$\rho$ cross-coupling—is the quantity that future neutron-rich nuclei measurements will sharpen.
Reading between the lines
- An implicit consequence is that the antiproton skin data and the GW170817 tidal measurement are not fully independent confirmations: both respond to the same isovector terms of the functional, so agreement with both tests the internal consistency of one fitted symmetry-energy density.
- Repeating the calculation for deformed nuclei such as $^{238}$U with quadrupole deformation and pairing would show whether the spherical-skin approximation shifts the predicted values by more than a few hundredths of a femtometer.
- Because $\Lambda \propto R^5$, the spread from $\Lambda_{1.4}=461$ (G3) to $681$ (IOPB-I) corresponds to only a few percent in radius; a future precision radius measurement of a 1.4-solar-mass neutron star would cleanly separate the two parameter sets.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This manuscript uses the extended relativistic mean-field models G3 and IOPB-I to predict the neutron-skin thickness Δrnp for 26 nuclei from 40Ca to 238U, compares the results with antiproton data and with the PREX-II 208Pb result, and computes the dimensionless tidal deformability Λ1.4 (and the binary weighted combination) for comparison with GW170817. The main reported numbers are Δrnp(208Pb) = 0.180 fm (G3) and 0.221 fm (IOPB-I), and Λ1.4 = 461.03 (G3), 622.06 (FSUGarnet), and 680.79 (IOPB-I). The paper concludes that both new parameter sets are consistent with the experimental skin data and with GW170817.
Significance. The NST part of the paper is a genuinely useful, non-circular test: the G3 and IOPB-I couplings were calibrated in refs [1,2] to eight spherical nuclei and to saturation properties, so the skin and tidal predictions are not fitted to the observables used for comparison. The G3 208Pb value and the IOPB-I value fall inside the PREX-II error band, and G3 tracks the antiproton systematics reasonably well. The tidal part, however, is not supported as stated: under the GW170817 re-analysis upper limit Λ1.4 ≤ 580 quoted by the paper itself, IOPB-I (681) and FSUGarnet (622) are excluded. The central claim therefore needs correction, though it is fixable within the paper's scope.
major comments (3)
- [Section 3, last paragraph] The text quotes both the first-analysis bound Λ ≤ 800 [6] and the re-analysis 90% upper limit Λ1.4 = 580 [7], and then states that the values 680.79 (IOPB-I), 622.06 (FSUGarnet), and 461.03 (G3) are consistent with the 90% credible intervals. Under the stricter re-analysis bound of 580, only 461.03 satisfies the constraint; 622.06 and 680.79 do not. The Conclusion's statement that "the calculated value Λ1.4 of a neutron star is consistent with the recent observation of GW170817" is therefore valid only for G3, or only if the weaker first-analysis bound is explicitly adopted. Please revise the comparison and the conclusion accordingly.
- [Section 3, NST comparison paragraph] The paper claims that both G3 and IOPB-I are consistent with the upper limit Δrnp ≤ 0.25 fm derived in [5] from the tidal deformability bound Λ1.4 ≤ 580. This is internally inconsistent for IOPB-I, whose own Λ1.4 = 680.79 violates the very tidal limit used to derive the skin bound. The consistency statement should be restricted to G3, or the skin-vs-tidal correlation should be recomputed using the actual IOPB-I Λ value.
- [Section 3, first paragraph] The NST calculations are performed "without considering the nuclear deformation and superfluidity," yet Fig. 1(a) compares the models with data for deformed nuclei such as 238U and 232Th. No estimate is given for the error introduced by this approximation. Because the claimed agreement spans the entire isotope chain, please quantify the expected shift (or cite deformed RMF calculations) or restrict the conclusion to nuclei for which the spherical approximation is justified.
minor comments (4)
- [Abstract and Section 1] The abstract and the introduction contain the typo "GW1701817"; it should read "GW170817."
- [Section 3, last paragraph] The sentence "From the GW170817, the values of Λ≤800 ... and Λ = 190+390−120 ... are within the 90% credible intervals which are consistent with ..." mixes bounds and model predictions; please rewrite it so that each model's Λ1.4 is explicitly compared with each bound.
- [Abstract and Section 3] The phrase "dimensional tidal deformability" should be "dimensionless tidal deformability," since Eq. (3) defines the dimensionless quantity Λ.
- [Fig. 1(b) caption] The caption says the figure is "adopted from [2]"; please clarify whether the Λ values quoted in the text are newly computed here or reproduced from [2], and whether the figure is original or reprinted with permission.
Circularity Check
No circularity: NST and Λ1.4 are computed from externally calibrated G3/IOPB-I models and checked against independent data.
full rationale
The derivation chain is not circular. The G3 and IOPB-I parameter sets are taken from refs [1,2]; their coupling constants were fitted to properties of eight spherical nuclei and nuclear-matter saturation constraints, not to the neutron-skin or tidal-deformability observables reported here. NST is computed directly as ∆rnp = Rn − Rp, and Λ1.4 is obtained by solving the TOV and perturbing equations, so neither target observable is an input to the model calibration. The comparisons with antiproton NST data and with the GW170817 bounds are external checks, not fitted outputs. The only caveats are: (i) the author of this paper is also an author of the papers that introduced G3 and IOPB-I, so there is self-citation; however, the cited model results are independently falsifiable against the data quoted here, so this is not load-bearing circularity; and (ii) the paper's statement that all three quoted Λ1.4 values are consistent with the GW170817 re-analysis is numerically questionable for IOPB-I (680.79) and FSUGarnet (622.06) if the 90% upper limit is 580. That is a consistency/correctness concern, not a circularity: the models' outputs are not defined in terms of, nor fitted to, the bounds with which they are compared.
Assumptions & free parameters
free parameters (2)
- G3 coupling constants =
Calibrated in ref [2]; not listed here
- IOPB-I coupling constants =
Calibrated in ref [2]; not listed here
assumptions (4)
- domain assumption The ERMF Lagrangian in Eq.(1), including nonlinear sigma terms, omega-rho coupling, and the delta meson, is a valid effective theory for finite nuclei and neutron star matter.
- domain assumption The mean-field approximation can safely neglect Fock terms and vacuum fluctuations once coupling constants are fitted.
- ad hoc to paper Nuclear deformation and superfluidity are negligible for the NST of all 26 nuclei.
- standard math The TOV equations and the tidal Love number formula (Eq.3) describe neutron star tides.
Cite this review
Pith. "Pith review of Neutron skins of heavy nuclei and tidal deformability of neutron star." pith.science (2026). https://pith.science/paper/FSV7J5OM
@misc{pith2026190802909,
author = {Pith},
title = {Pith review of: Neutron skins of heavy nuclei and tidal deformability of neutron star},
year = {2026},
howpublished = {\url{https://pith.science/paper/FSV7J5OM}},
note = {Machine review of arXiv:1908.02909}
}
abstract
In this paper, I have discussed the numerical predictions for the neutron-skin thickness (NST) of various finite nuclei starting from $^{40}$Ca to $^{238}$U using recently developed effective relativistic mean-field models G3 and IOPB-I \cite{G3, IOPB}. The calculated results are compared with the PREX-II data, and the experiment has been done with antiprotons at CERN. Further, I have also calculated the dimensional tidal deformability of a canonical neutron star 1.4$M_\odot$ and compared it with the recent observation of GW1701817.
Figures
Reference graph
Works this paper leans on
-
[2]
Bharat Kumar, S. K. Patra, and B. K. Agrawal, Phys. Rev. C 97 (2018) 045806
work page 2018
-
[6]
P . B. Abbott et al., Phys. Rev. Lett. 119 (2017) 161101
work page 2017
-
[7]
P . B. Abbott et al., Phys. Rev. Lett. 121 (2018) 161101
work page 2018
-
[5]
F. J. Fattoyev, J. Piekarewicz, and C. J. Horowitz, Phys. Rev. Lett. 120 (2018) 172702
work page 2018
-
[1]
Bharat Kumar, S. K. Singh, B. K. Agrawal, and S. K. Patra, N ucl. Phys. A 966 (2017) 197
work page 2017
-
[3]
BA TES Laboratory at MIT, http: //bateslab.mit.edu/
- [4]
-
[8]
R. J. Furnstahl, B. D. Serot and H. B. Tang, Nucl. Phys. A 598 (1996) 539; R. J. Furnstahl, B. D. Serot and H. B. Tang, Nucl. Phys. A 615 (1997) 441
work page 1996
Show all 10 references
-
[9]
Trzci´ nska, J
A. Trzci´ nska, J. Jastrzebski, P . Lubi´ nski, F. J. Hartmann, R. Schmidt, T. von Egidy, and B. Klos, Phys. Rev. Lett. 87 (2001) 082501; J. Jastrze ¸bski, A. Trzci´ nska, P . Lubi´ nski, B. Klos, F. J. Hartmann, T. von Egidy, S. Wycech, Int. J. Mod. Phys. E 13 (2004) 343
2001
-
[10]
Lackey, Ryan N
Tanja Hinderer, Benjamin D. Lackey, Ryan N. Lang, and Jocelyn S. Read, Phys. Rev. D81 (2010) 123016. 4
2010
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.