Effective Mass in Quantum Hadrodynamics-I and its Impact on the Equation of State of Neutron Matter

Eko Tri Sulistyani; Ghitha Nadhira Azka Rahiemy; Pekik Nurwantoro

REVIEW 2 major objections 5 minor 19 references

Effective Mass in Quantum Hadrodynamics-I and its Impact on the Equation of State of Neutron Matter

T0 review · 2 major / 5 minor · reviewed 2026-08-04 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read The paper confirms that Walecka's self-consistent effective-mass equation is numerically robust, and that variations in meson–nucleon couplings propagate through the effective mass to control the stiffness of the QHD-I neutron-matter equati

desk verdict Clean, correct reproduction of Walecka's QHD-I effective mass calculation; no new physics, but a solid pedagogical and code baseline. read the letter →

arxiv 2509.10176 v1 pith:4XOLGDPL submitted 2025-09-12 nucl-th astro-ph.HE

Ghitha Nadhira Azka Rahiemy , Eko Tri Sulistyani , Pekik Nurwantoro This is my paper

classification nucl-thastro-ph.HE

keywords QuantumHadrodynamicseffectivenucleonmassrelativisticmean-fieldtheoryneutronmatterequationofstateGaussianquadratureWaleckamodelnuclearsaturationstar

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper re-derives the scalar-density integral of the Walecka (QHD-I) model step by step and validates the self-consistent equation for the effective nucleon mass M* by comparing 20-point Gaussian quadrature with the analytic integral. The numerical and analytic curves overlap nearly perfectly, so the paper treats this as confirmation that Walecka's M* relation is correct and robust. The authors then vary three parameter sets and show that the resulting differences in M* translate directly into differences in saturation density, compressibility, and the stiffness of the neutron-matter equation of state. The larger point is that M* functions as the bridge that carries microscopic scalar/vector meson dynamics into the macroscopic pressure–energy relation used for neutron-star modeling.

What carries the argument

The load-bearing object is the self-consistent effective nucleon mass M* = M - g_s φ_0, which enters the scalar density integral (Eq. 2). The paper pairs the analytic evaluation (Eq. 3), obtained by integration by parts and a hyperbolic substitution, with a 20-point Gauss–Legendre quadrature and a Secant root-finder for the fixed-point equation f(M*)=0. M* then feeds the QHD-I energy density and pressure formulas, so the same machinery shows how coupling choices propagate into EoS stiffness.

What would settle it

Evaluate the scalar-density integral at k_F ≈ 6 fm^-1 with a 50- or 100-point Gauss–Legendre rule and compare it to the analytic formula (Eq. A27); if the relative difference exceeds the stated 10^-3 tolerance, or if the M* fixed point shifts by more than 10^-3 M, the numerical robustness claim collapses. A second check is to run the Secant solver from a different starting guess and verify that the same root is found.

Watch

Extended reading notes

Core claim

The central claim is that the effective nucleon mass M* = M - g_s φ_0 is correctly and robustly reproduced by Walecka's analytic scalar-density integral, and that M* is what channels parameter-set differences into the equation of state. Using 20-point Gaussian quadrature and a Secant iteration for the self-consistency condition, the paper reports near-perfect overlap with the analytic M* curve for neutron matter (γ=2) and symmetric nuclear matter (γ=4). It finds Walecka, Serot-Walecka, and RHA0 parameter sets give systematically different M*(k_F), binding energies, and pressures, with all QHD-I equations of state staying stiff at high density and approaching the causal limit p=ε.

Load-bearing premise

The validation claim rests on the untested assumption that the 20-point Gaussian quadrature and the Secant iteration converge to the true solution within the stated 10^-3 tolerance, so the reported near-perfect overlap with the analytic result is asserted rather than quantified.

Editorial extensions

If this is right

The re-derivation and numerical agreement confirm that the standard QHD-I self-consistency relation is dependable enough to serve as a baseline for extended models.
The computed drop of M*/M with density reproduces the characteristic saturation dip, with neutron matter remaining unbound in all three parameter sets.
The Serot-Walecka set produces the steepest pressure rise and RHA0 the most gradual, so parameter choice directly alters the predicted stiffness of neutron-star matter.
All QHD-I equations of state are stiff and approach the causal limit p=ε at high density, so matching observed neutron-star masses and radii requires physics beyond QHD-I.
The same numerical framework can be reused, with minimal modification, for extended RMF models that include nonlinear couplings, ρ mesons, or β-equilibrated matter.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

Because the paper reports near-perfect overlap without error bars or a convergence study, a direct test of quadrature-order dependence would settle whether the stated 10^-3 tolerance is actually met; the analytic formula makes this check trivial.
The sensitivity of M* to couplings implies a potentially tight link between neutron-star radius and mass observations and the scalar/vector coupling constants, a constraint the paper gestures at but does not quantify.
One could differentiate the fixed-point equation to obtain an analytic dM*/dC_s² and thereby predict the ordering of EoS stiffness without solving the full EoS, an extension the paper leaves implicit.
The framing of M* as a microscopic–macroscopic bridge suggests that future RMF studies should report M*(k_F) as a primary diagnostic, since it seems to encode most of the parameter dependence.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Short version: this is a careful re-implementation of the Walecka QHD-I effective mass equation, and it checks out. The math is right, the numerics match the analytic integral, and the code is on GitHub. But there is no new physics here—the central integral is Walecka 1974, and the EoS is the standard QHD-I result from Serot and Walecka. The only genuinely new element is the 20-point Gaussian quadrature routine, which is routine computational work.

What the paper does well: the step-by-step re-derivation in Appendix A is algebraically correct (the square on M* in the denominator is missing in Eq. A1, but that's clearly a typo, since Eq. A2 and the final result have it), and the numerical validation is honest. The authors properly note that once the couplings are fixed, everything else is a prediction of the model, and they don't overclaim the parameter-set comparison as new. The GitHub availability is a plus—someone can rerun this easily.

Soft spots: The 'impact' framing oversells the parameter sensitivity study. Figures 5 and 6 are restatements of known QHD-I behavior; Diener's thesis and Serot-Walecka already show this. The convergence claim for 20-point Gaussian quadrature is asserted, not quantified—no error bars or convergence study, so the 'near-perfect overlap' in Fig. 3 is visual, not measured. That's a minor presentational issue, not a correctness risk: the integral is smooth, and 20-point quadrature gives machine precision in practice. Also, the paper does not compare to beta-equilibrated matter or discuss the symmetry energy, which limits the astrophysical punch, but the authors acknowledge that QHD-I is a baseline.

Bottom line: this is a solid reproduction, not a contribution to the research frontier. It should not be cited as a source for new physics, but it would be a fine pedagogical reference for someone learning RMF and numerical evaluation of the scalar density integral.

For peer review: if it landed on my desk at a journal that accepts pedagogical or code-baseline papers, I'd send it out—the derivation is correct, the code is available, and the honest framing makes it a useful check on the literature. For a research journal that expects new results, it's a desk reject.

Referee Report

2 major / 5 minor

Summary. The manuscript revisits the QHD-I (Walecka) effective nucleon mass M*, re-derives the scalar density integral in Appendix A (Eq. A27), validates the integral numerically with 20-point Gaussian quadrature, and then uses three parameter sets (Walecka, Serot-Walecka, RHA0) to compute M*/M, binding energy, and the equation of state for neutron matter and symmetric nuclear matter. The central claim is that Walecka's self-consistent M* equation is correct and numerically robust, and that parameter-set variations propagate through M* to control the stiffness of the equation of state.

Significance. If the claims hold, the paper confirms a standard QHD-I result and provides a clearly written derivation plus openly available code. The physics is not new: the integral is textbook material and the equation-of-state expressions are those of standard QHD-I. The value of the paper lies in its verification/pedagogical role and in establishing a numerical baseline for possible RMF extensions. Strengths include the explicit analytic reduction in Appendix A, the availability of the code, and the use of standard EoS forms. The reported agreement between numerical and analytic results is visually good, but the paper lacks a quantitative convergence study, and one central table definition is inconsistent.

major comments (2)

[Table 1; Eq. (2)] Table 1 defines C_s^2 ≡ g_s^2 (M^2/m_s^2)^2 and C_v^2 ≡ g_v^2 (M^2/m_v^2)^2, but the listed numerical values (266.9, 195.7, etc.) are the standard QHD-I couplings C_s^2 = g_s^2 (M/m_s)^2 and C_v^2 = g_v^2 (M/m_v)^2. As printed, the definition changes the conversion between C_i^2 and g_i^2/m_i^2 used in Eq. (2) by an extra factor (M^2/m_s^2) or (M^2/m_v^2), so the quantitative results cannot be reproduced without resolving this inconsistency. Please correct the definition and state explicitly which convention was used in the code.
[Sec. 2.2; Fig. 3] The paper asserts that a 20-point Gaussian quadrature gives convergence to the analytic result and that the overlap is 'almost perfect', but no convergence study, error table, or residual check is reported. Since the central claim includes numerical robustness, the authors should report the maximum relative error between the quadrature and Eq. (3) over the k_F range, a convergence test in quadrature order, and the residual of f(M*) after the Secant iteration. Without these, 'near-perfect overlap' remains a visual assertion rather than a quantified validation.

minor comments (5)

[Eq. (A1)] The denominator under the square root in Eq. (A1) is written as sqrt(k^2 + M*), missing the square on M*. It should be sqrt(k^2 + M*^2), as used correctly in Eq. (A27).
[Eq. (6)] The displayed baryon density formula is written as ρ_B = γ/(2π^2) ∫_0^{k_F} dk = γ/(6π^2) k_F^3. If the integrand is literally dk, the equality is wrong; the integrand should be k^2 dk. The final expression is correct, but the displayed equation should be fixed.
[Sec. 3.3; Fig. 8] The text says the logarithmic pressure-energy density relation spans about 10^13.4–10^14.5 g/cm^3, but Fig. 8 shows Log10 E from approximately 12 to 16.5. Please clarify the intended density range.
[Figures captions] Several figure captions contain OCR/typing errors, e.g. 'Walecka anal tic', 'Effecti e Mass', 'Ne tron Matter', and 'N clear Matter'. These should be corrected.
[Sec. 3.2] The phrase 'systematic exploration of how variations in meson masses and couplings influence' overstates the analysis, since only three combined parameter sets are compared and the individual roles of each coupling are not isolated. Rephrasing as 'comparison of representative parameter sets' would be more accurate.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the paper re-derives and numerically validates a known integral; the central claims are consequences of an independently published model, not fits to the paper's own outputs.

full rationale

The derivation chain is self-contained and non-circular. The self-consistent effective mass relation (Eq. 2) is taken from Walecka (ref [7]), and the analytic evaluation (Eq. 3) is re-derived in Appendix A by standard integration by parts and substitution (Eqs. A1-A27). The numerical evaluation uses 20-point Gaussian quadrature to integrate the same scalar-density integral; matching the analytic expression is a genuine check of algebraic prefactors, not a fit. M* is obtained by solving Eq. (4) with the Secant method, and no parameter is fitted in this paper; the three parameter sets (Walecka, Serot-Walecka, RHA0) are inherited from prior fits to nuclear saturation or Dirac-Hartree benchmarks, as the paper explicitly states ('Once these two parameters are fixed, all other quantities in the model are predictions'). The EoS results (Eqs. 6-8) are standard QHD-I expressions evaluated with those inputs, and the paper's 'impact' statements describe direct propagation of input couplings through M*, which is a model consequence rather than a prediction extracted from the paper's own outputs. There is no self-citation chain: refs [7] and [16] are external prior work, and the present authors do not cite themselves for any load-bearing premise. The only weakness—absence of a formal convergence study for the 20-point quadrature—is a presentational issue, not circularity.

Assumptions & free parameters 3 free parameters · 5 assumptions · 0 invented entities

The paper fits no new parameters; all coupling constants are inherited from prior fits. Its central validation claim rests on standard calculus and the QHD-I mean-field approximation. The quadrature order is an ad hoc numerical choice without a published convergence test.

free parameters (3)

C_s^2, C_v^2 (Walecka set) = 266.9, 195.7
Fitted in Walecka (1974) to nuclear saturation; used as input in Table 1.
C_s^2, C_v^2 (Serot-Walecka set) = 357.4, 273.8
From Serot and Walecka (1992); used as input in Table 1.
C_s^2, C_v^2 (RHA0 set) = 228.0, 147.5
RHA0 set from Serot and Walecka (1992); used as input in Table 1.

assumptions (5)

domain assumption QHD-I mean-field approximation: scalar and vector meson fields are replaced by their expectation values
Invoked in Eqs. (1)-(2) and (6)-(8); standard for the model but a physical approximation.
standard math Standard integral identities: integration by parts and x = a sinh(theta) substitution for the J integral
Used in Appendix A (Eqs. A3-A26) to evaluate the scalar density integral.
domain assumption Pure neutron matter degeneracy gamma = 2 and symmetric matter gamma = 4
Stated in Sec. 2.1 following Walecka [7]; all neutron matter figures use gamma = 2.
domain assumption Parameter sets from prior literature (Walecka, Serot-Walecka, RHA0) are valid inputs for dense matter
Table 1; the sensitivity analysis trusts these fits without reassessment at high density.
ad hoc to paper A 20-point Gaussian quadrature is sufficient to converge the scalar density integral
Adopted in Sec. 2.2 as a compromise; no convergence test is reported to justify this order.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Effective Mass in Quantum Hadrodynamics-I and its Impact on the Equation of State of Neutron Matter." pith.science (2026). https://pith.science/paper/4XOLGDPL

@misc{pith2026250910176,
  author       = {Pith},
  title        = {Pith review of: Effective Mass in Quantum Hadrodynamics-I and its Impact on the Equation of State of Neutron Matter},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/4XOLGDPL}},
  note         = {Machine review of arXiv:2509.10176}
}

read the original abstract

The effective nucleon mass (M^*) plays a central role in Quantum Hadrodynamics-I (QHD-I), linking scalar meson interactions at the microscopic level to the macroscopic properties of dense nuclear matter. In this work, we re-derive the scalar density integral in detail and validate it numerically using Gaussian quadrature. The numerical and analytic results are found to be in excellent agreement, confirming the robustness of both approaches. We then investigate the sensitivity of M^* to different parameter sets, highlighting its strong influence on nuclear saturation, compressibility, and the resulting equation of state (EoS). The analysis shows that variations in meson-nucleon couplings propagate directly into differences in pressure and energy density, affecting the stiffness of the EoS. While QHD-I produces characteristically stiff EoS, the effective mass evaluation provides a transparent framework for connecting microscopic meson dynamics to macroscopic neutron star properties. These findings underline the relevance of M^* as a microscopic-macroscopic bridge and demonstrate the utility of numerical methods for extending relativistic mean-field models in nuclear astrophysics.

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

19 extracted references · 7 canonical work pages

[1]

Technical report, Technische Univ

Antic, S., al.: Generalized relativistic mean-ﬁeld model with non-linear deriva- tive nucleon-meson couplings for nuclear matter and ﬁnite nuclei. Technical report, Technische Univ. Darmstadt (Germany). Fachbereich Physik (2018)

2018
[2]

Relativistic mean-field theory applied to the study of neutron star properties

Diener, J.P.W.: Relativistic mean-ﬁeld the- ory applied to the study of neutron star properties. PhD thesis, Stellenbosch University (2008). https://doi.org/10.48550/ ARXIV.0806.0747

work page Pith review arXiv doi:10.48550/arxiv.0806.0747 2008
[3]

Cambridge Mono- graphs on Mathematical Physics

Frishman, Y., Sonnenschein, J.: Non- Perturbative Field Theory: From Two Dimensional Conformal Field Theory to QCD in Four Dimensions. Cambridge Mono- graphs on Mathematical Physics. Cambridge University Press, Cambridge, UK (2010)

2010
[4]

Dupuis, N., Canet, L., Eichhorn, A., Metzner, W., Pawlowski, J.M., Tissier, M., Wschebor, N.: The nonperturbative functional renor- malization group and its applications. Phys. Rep. 910, 1–114 (2021) https://doi.org/10. 1016/j.physrep.2021.01.001

2021
[5]

Abbott, R., Detmold, W., Illa, M., Par- reño, A., Perry, R.J., Romero-López, F., Shanahan, P.E., Wagman, M.L.: Qcd con- straints on isospin-dense matter and the nuclear equation of state. Phys. Rev. Lett. 134(1), 011903 (2025) https://doi.org/10. 1103/PhysRevLett.134.011903

2025
[6]

Neutron stars and Constraints for the Equation of State of Dense Matter

Kumar, R., Dexheimer, V., Jahan, J.: Neu- tron stars and Constraints for the Equation of State of Dense Matter. arXiv preprint arXiv:2503.23413 (2025). https://doi.org/10. 48550/ARXIV.2503.23413

work page Pith review arXiv doi:10.48550/arxiv.2503.23413 2025
[7]

Walecka, J.D.: A theory of highly con- densed matter. Ann. Phys. 83, 491–529 (1974) https://doi.org/10.1016/0003- 4916(74)90208-5

doi:10.1016/0003- 1974
[8]

Plenum Press, New York, NY (1986)

Negele, J.W., Serot, B.D., Vogt, E., Walecka, J.D.: The Relativistic Nuclear Many- body Problem. Plenum Press, New York, NY (1986). https://www.osti.gov/biblio/ 8 5443763

1986

Show all 19 references

[9]

Erkelenz, K.: Current status of the relativis- tic two-nucleon one boson exchange poten- tial. Phys. Rep. 13, 191–258 (1974) https:// doi.org/10.1016/0370-1573(74)90008-8

1974 doi
[10]

Holinde, K.: Two-nucleon forces and nuclear matter. Phys. Rep. 68, 121– 188 (1981) https://doi.org/10.1016/0370- 1573(81)90188-5

1981 doi
[11]

Wydawnictwo Uniwer- sytetu Śląskiego, Katowice, Poland (2007)

Bednarek, I.: Relativistic Mean Field: Mod- els of Neutron Stars. Wydawnictwo Uniwer- sytetu Śląskiego, Katowice, Poland (2007)

2007
[12]

: Constraining a relativis- tic mean ﬁeld model using neutron star mass–radius measurements i: nucleonic mod- els

Huang, C., et al. : Constraining a relativis- tic mean ﬁeld model using neutron star mass–radius measurements i: nucleonic mod- els. Mon. Not. R. Astron. Soc. 529(4), 4650– 4665 (2024) https://doi.org/10.1093/mnras/ stae844

2024 doi
[13]

Pais, H., Ivanytskyi, O., Providência, C.: Landau parameters and entrainment matrix of cold stellar matter: eﬀect of the sym- metry energy and strong magnetic ﬁelds. J. Cosmol. Astropart. Phys. 2022(04), 024 (2022) https://doi.org/10.1088/1475-7516/ 2022/04/024

2022 doi
[14]

arXiv preprint arXiv:2407.03113 (2024)

Scurto, L., Carvalho, V., Pais, H., Providên- cia, C.: Assessing the joint eﬀect of tem- perature and magnetic ﬁeld on the neu- tron star equation of state. arXiv preprint arXiv:2407.03113 (2024). https://doi.org/10. 48550/ARXIV.2407.03113

2024 doi arXiv
[15]

Wang, X., Li, J., Fang, J., Pais, H., Providên- cia, C.: Pasta phases in neutron stars under strong magnetic ﬁelds. Phys. Rev. D 105(6), 063004 (2022) https://doi.org/10. 1103/PhysRevD.105.063004

2022
[16]

In: Ainsworth, T.L., Campbell, C.E., Clements, B.E., Krotscheck, E

Serot, B.D., Walecka, J.D.: Relativistic nuclear many-body theory. In: Ainsworth, T.L., Campbell, C.E., Clements, B.E., Krotscheck, E. (eds.) Recent Progress in Many-Body Theories: Volume 3, pp. 49–92. Springer, Boston, MA (1992). https://doi. org/10.1007/978-1-4615-3466-2_5

1992 doi
[17]

Lalazissis, G.A., König, J., Ring, P.: New parametrization for the lagrangian density of relativistic mean ﬁeld theory. Phys. Rev. C 55(1), 540–543 (1997) https://doi.org/10. 1103/PhysRevC.55.540

1997
[18]

Long, W., Meng, J., Giai, N.V., Zhou, S.- G.: New eﬀective interactions in relativistic mean ﬁeld theory with nonlinear terms and density-dependent meson-nucleon coupling. Phys. Rev. C 69(3), 034319 (2004) https:// doi.org/10.1103/PhysRevC.69.034319

2004 doi
[19]

Piekarewicz, J.: Validating relativistic mod- els of nuclear structure against theoretical, experimental, and observational constraints. Phys. Rev. C 76(6), 064310 (2007) https:// doi.org/10.1103/PhysRevC.76.064310 9

2007 doi

Pith tools

Reviewed August 4, 2026 · model on record in the stance chip above.

[1] [1]

Technical report, Technische Univ

Antic, S., al.: Generalized relativistic mean-ﬁeld model with non-linear deriva- tive nucleon-meson couplings for nuclear matter and ﬁnite nuclei. Technical report, Technische Univ. Darmstadt (Germany). Fachbereich Physik (2018)

2018

[2] [2]

Relativistic mean-field theory applied to the study of neutron star properties

Diener, J.P.W.: Relativistic mean-ﬁeld the- ory applied to the study of neutron star properties. PhD thesis, Stellenbosch University (2008). https://doi.org/10.48550/ ARXIV.0806.0747

work page Pith review arXiv doi:10.48550/arxiv.0806.0747 2008

[3] [3]

Cambridge Mono- graphs on Mathematical Physics

Frishman, Y., Sonnenschein, J.: Non- Perturbative Field Theory: From Two Dimensional Conformal Field Theory to QCD in Four Dimensions. Cambridge Mono- graphs on Mathematical Physics. Cambridge University Press, Cambridge, UK (2010)

2010

[4] [4]

Dupuis, N., Canet, L., Eichhorn, A., Metzner, W., Pawlowski, J.M., Tissier, M., Wschebor, N.: The nonperturbative functional renor- malization group and its applications. Phys. Rep. 910, 1–114 (2021) https://doi.org/10. 1016/j.physrep.2021.01.001

2021

[5] [5]

Abbott, R., Detmold, W., Illa, M., Par- reño, A., Perry, R.J., Romero-López, F., Shanahan, P.E., Wagman, M.L.: Qcd con- straints on isospin-dense matter and the nuclear equation of state. Phys. Rev. Lett. 134(1), 011903 (2025) https://doi.org/10. 1103/PhysRevLett.134.011903

2025

[6] [6]

Neutron stars and Constraints for the Equation of State of Dense Matter

Kumar, R., Dexheimer, V., Jahan, J.: Neu- tron stars and Constraints for the Equation of State of Dense Matter. arXiv preprint arXiv:2503.23413 (2025). https://doi.org/10. 48550/ARXIV.2503.23413

work page Pith review arXiv doi:10.48550/arxiv.2503.23413 2025

[7] [7]

Walecka, J.D.: A theory of highly con- densed matter. Ann. Phys. 83, 491–529 (1974) https://doi.org/10.1016/0003- 4916(74)90208-5

doi:10.1016/0003- 1974

[8] [8]

Plenum Press, New York, NY (1986)

Negele, J.W., Serot, B.D., Vogt, E., Walecka, J.D.: The Relativistic Nuclear Many- body Problem. Plenum Press, New York, NY (1986). https://www.osti.gov/biblio/ 8 5443763

1986