{"id":"ea29f525-6d25-4eaf-afba-ddb6e18a63eb","arxiv_id":"2607.23253","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Using the Landau instead of the Dirac effective nucleon mass in neutron-star cooling calculations cools massive stars faster, changing predicted surface temperatures by ~0.03–0.06 dex.","lead":"Inside neutron stars, nucleons have two different 'effective masses' depending on how they are defined, and this paper shows the choice changes predicted cooling. The effect is comparable to observational errors in well-measured cooling stars, so neutron-star cooling calculations should specify which mass they use.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paper never derives which effective mass enters the imported emissivity formulas; if only the Landau mass belongs there, the 'Dirac mass' cooling branch is unphysical and the claimed 0.03–0.06 dex systematic uncertainty is overstated.","rationale":"Reader's weakest assumption is exactly the load-bearing point. The paper is otherwise a competent sensitivity study: the Walecka derivation of m_D and m_L is clean, the analytic relation m_L = sqrt(k_F^2 + m_D^2) is correct, and the NSCool simulations with TM1e/TM1m provide independent support that EOS-level effective mass differences alter cooling curves. The problem is the interpretation of those simulations as an 'uncertainty' between definitions. The imported emissivity formulas have a definite microscopic content; one cannot swap m_D and m_L in them without derivation. Since the non-relativistic formulas' phase-space/density-of-states factors are governed by the Landau mass, the Dirac branch is the likely suspect. However, this does not warrant REJECT: the qualitative claim that effective-mass treatment affects cooling is plausible and partly supported by the EOS comparison; the quantitative 0.03–0.06 dex uncertainty needs to be re-derived. CONDITIONAL remains the appropriate verdict, contingent on the authors either deriving the mass content of the emissivities or rephrasing the paper as a comparison of two ad hoc prescriptions rather than a physical systematic uncertainty. Thus no verdict adjustment.","tokens_in":18651,"tokens_out":6900,"duration_ms":60540,"concrete_test":"Re-derive Eq. (36) and Eq. (42) from the RMF/Walecka Lagrangian: compute the neutrino-emission amplitude with quasiparticle energies ε_k = sqrt(k^2 + m_D^2), μ* = ε_{k_F}, and density of states D(ε_F) = m_L k_F^2/π^2, then take the low-temperature/non-relativistic limit. If the resulting prefactor is ∝ m_L^4 (bremsstrahlung) and ∝ m_L,n m_L,p (DU) — as the density-of-states factors imply — the m_D curves in Figs. 2–3 are unphysical; if it is ∝ m_D^4 and m_D^2, the m_L curves are unphysical. A faster analytical check: take the non-relativistic limit of the relativistic DU expression Eq. (41) with ε_i = sqrt(p_i^2 + m_D^2) and show whether it reduces to Eq. (42) with m*_N m*_p = m_L,n m_L,p or m_D^2.","verdict_should_be":"UNCHANGED","load_bearing_attack":"At Sec. III A the authors import non-relativistic emissivities, Eq. (36) for nn bremsstrahlung and Eq. (42) for direct Urca, and then announce: 'In the numerical comparison below, we evaluate this expression using either m*_N = m_D or m*_N = m_L.' Nowhere do they derive from the Walecka Lagrangian which mass should appear. This is not a harmless choice. In their own RMF derivation, the two masses play different roles: ε_k = sqrt(k^2 + m_D^2) sets the single-particle spectrum, while the density of states and specific heat are controlled by m_L = sqrt(k_F^2 + m_D^2), via D(ε_F) = m_L k_F^2/π^2 (Eq. 30) and c_V ∝ m_L (Eq. 31). The non-relativistic emissivities were derived in a Fermi-liquid context where m* in phase-space factors is the quasiparticle mass, i.e., the Landau mass. Using m_D in those factors is therefore not a legitimate alternative prescription but a category error: it replaces a density-of-states quantity by a scalar self-energy mass. The 'Dirac' cooling branch is thus either internally inconsistent (if c_V remains m_L while the emissivity uses m_D) or unphysical (if c_V is also switched to m_D, contradicting Eq. 31). The central conclusion that the m_D vs m_L choice is a non-negligible systematic uncertainty therefore is not secured by the paper's calculation; it may be an artifact of applying the wrong mass to formulas that already encode m_L.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper addresses the role of the nucleon effective mass in neutron-star cooling, focusing on the difference between the Dirac mass m_D and the Landau mass m_L in relativistic mean-field (Walecka) models. The authors derive m_D and m_L at finite temperature, show that m_L = sqrt(k_F^2 + m_D^2), and obtain the low-temperature specific heat c_V proportional to m_L. They then evaluate standard neutrino emissivities for neutron-neutron bremsstrahlung and direct Urca processes by substituting either m_D or m_L for the effective nucleon mass, solve isothermal cooling equations, and perform NSCool simulations with the TM1e and TM1m equations of state. The central claim is that the choice between Dirac and Landau effective masses is a non-negligible systematic uncertainty in neutron-star cooling, changing predicted surface temperatures by roughly 0.03-0.06 dex for 1.4-2.0 solar-mass stars.","tokens_in":19060,"tokens_out":17608,"duration_ms":157023,"significance":"If the central comparison were justified, the result would be valuable: many cooling codes use a single effective mass without distinguishing Dirac and Landau definitions, and the paper offers a clean RMF derivation showing that the two masses differ substantially at high density. The low-temperature expansion and the derivation of c_V proportional to the Landau mass are sound and clearly presented. The use of the public NSCool code and two realistic EOS models is a welcome check on the simplified Walecka-model results. However, the load-bearing step — evaluating imported emissivity formulas with either m_D or m_L as interchangeable 'prescriptions' — is not derived from the underlying microphysics. If only one of the two masses is the correct input for those formulas, the claimed uncertainty is an artifact, not a physical systematic effect. The qualitative direction of the effect is also forced by the algebraic relation m_L > m_D, so the paper's quantitative claim is not yet secured.","major_comments":[{"comment":"This is the central issue. The paper never derives from the Walecka Lagrangian which effective mass enters the non-relativistic bremsstrahlung emissivity. In the RMF framework, the quasiparticle dispersion near the Fermi surface is governed by v_F = k_F/m_L (Eq. 14), and the density of states and specific heat are controlled by m_L (Eqs. 30-31). The non-relativistic emissivity Eq. (36) is a standard Fermi-liquid result; its m*^4 factor is the quasiparticle mass, i.e. the Landau mass. Substituting m_D into that factor is therefore not an alternative prescription but a category error. Moreover, the cooling curves labeled 'Dirac' in Fig. 2 combine an emissivity evaluated with m_D with a heat capacity evaluated with m_L, so they are internally inconsistent. As a result, the claimed 0.03-0.06 dex systematic uncertainty is not established by the calculation. The authors must either derive the","section":"Sec. III A, Eq. (36) and the sentence 'In the numerical comparison below, we evaluate this expression using either m*_N"},{"comment":"The same problem appears in the direct Urca treatment. The paper states that m*_N and m*_p in the DU emissivity 'are identified with' either Landau or Dirac masses, but it does not derive which identification is correct. In the relativistic formula Eq. (41), the single-particle energies ϵ_n and ϵ_p already contain m_D, and the role of the m*_N m*_p term must be established from the weak-current matrix element and the phase-space integral. If m* in the Leinson formula is the Dirac mass, then the Landau branch is unphysical; if m* is the Landau mass, then the Dirac branch is unphysical. The non-relativistic reduction Eq. (42) differs by a factor (m_L/m_D)^2 between the two prescriptions, so this is not a harmless ambiguity. A derivation, not a choice, is needed before these curves can be interpreted as a physical uncertainty. The same criticism applies to the NSCool runs in Sec. III B that","section":"Sec. III A, Eqs. (41)-(42)"},{"comment":"Even if both mass prescriptions were taken at face value, the qualitative result would be a mathematical consequence of Eq. (15): m_L = sqrt(k_F^2 + m_D^2) > m_D, and all emissivities increase with the effective mass. Thus the statement that 'uncertainties in the effective masses provide a non-negligible contribution to the systematic uncertainty' requires demonstrating that the community actually uses both definitions, or that the correct mass cannot be fixed by existing many-body theory. The paper does not provide such evidence. The TM1e/TM1m comparison is a legitimate model-parameter uncertainty, but it concerns the value of the Dirac mass at saturation, not the Dirac-Landau definitional ambiguity. The paper should separate these two distinct claims and either derive the correct mass or substantially weaken the conclusion.","section":"Sec. III A and Conclusion"}],"minor_comments":[{"comment":"The text says 'Panel (a) denotes the temperature evolution ... while panel (b) shows E_NN normalized by m_psi^3/T^8', but the panels are reversed: panel (a) is the emissivity and panel (b) is the temperature. Also, the sentence 'using the Landau mass enhances the emissivity ... and results in a slower decrease of the temperature' appears to be wrong: a larger emissivity with c_V ∝ m_L gives a cooling rate ∝ m_L^3, so the Landau mass should cool faster (steeper temperature decrease), not slower.","section":"Fig. 2 and text after it"},{"comment":"There is a typo: 'Compared to the Laudau mass' should read 'Landau mass'.","section":"Sec. II C"},{"comment":"The phrase 'somehow theoretical treatment' is awkward and should be rephrased.","section":"Sec. III A"},{"comment":"The step in which the neutron and proton masses are assumed equal (m*_N = m*_p = m_D) should be stated more prominently; this is an approximation that may not be justified for neutron-rich matter, where the proton fraction is small.","section":"Sec. III A, Eq. (41)"},{"comment":"The entry 'C Cabibbo factor' is ambiguous; please write, e.g., 'Cabibbo factor C' or define C explicitly.","section":"Table I"}],"recommendation":"major_revision","confidential_remarks":"The reader's stress-test concern is well-founded and is the decisive issue. The manuscript's RMF thermodynamics and c_V derivation are sound, but the central comparison of m_D and m_L in the emissivities is not derived. The authors need to either provide a microscopic derivation of which effective mass enters each emissivity or reframe the paper as a study of the sensitivity of cooling curves to the mass prescription, with the caveat that only one prescription is physically correct. Without this, the abstract's claim of a non-negligible systematic uncertainty is overreaching. The TM1e/TM1m comparison is a useful positive element that should be emphasized more, as it is a well-defined physical model difference."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper does something useful: it lays out the Dirac and Landau effective masses in the Walecka model, derives the low-temperature specific heat, and then asks which mass should go into standard neutrino emissivity formulas. It also goes beyond toy models by running NSCool with the TM1e and TM1m equations of state, which differ mainly in the effective mass. The result that cooling curves shift by ~0.03–0.06 dex is a concrete, quantitative statement about a real source of uncertainty. The TM1e-vs-TM1m comparison is, on its own, a legitimate demonstration that effective-mass differences matter for cooling predictions. That part of the paper deserves credit.\n\nThe soft spot is the one the reader flagged. The emissivity formulas are imported from Fermi-liquid literature, where the mass in the phase-space factors is the quasiparticle mass—the Landau mass. The paper never derives from the Walecka Lagrangian that the Dirac mass is a legitimate substitute in those formulas. Since the paper's own derivation shows the density of states and specific heat are controlled by m_L = sqrt(k_F^2 + m_D^2), using m_D in the emissivity while keeping c_V ∝ m_L is an inconsistent hybrid. The Dirac-mass cooling branch may not correspond to any physically realizable model. The qualitative claim that the m_D/m_L distinction is a systematic uncertainty is therefore overstated: the comparison may be between the correct formula and an incorrect one. Also, the Fig. 2 discussion has a caption/text swap, and the statement that the Landau mass \"results in a slower decrease of the temperature\" is backwards—larger emissivity means faster cooling. The claim that most previous cooling simulations do not distinguish the masses is also too broad without a systematic survey.\n\nStill, the paper is not incoherent. The TM1e/TM1m result stands regardless of the Dirac-vs-Landau pairing, and the authors are upfront about limitations. The Walecka-model derivation of c_V and the low-T expansion is clean and standard. The main issue is fixable: derive which mass enters the emissivity, or reframe the paper as an exploratory sensitivity study rather than a definite uncertainty estimate. With that revision, the quantitative conclusions would be much more solid.\n\nI would send this to peer review, not desk reject it, because the question is relevant and the TM1e/TM1m cooling comparison is a real contribution. A good referee could force the authors to tighten the derivation and clean up the presentation. I'd bring it to a reading group as a case study in how effective masses can be mishandled in a popular framework.","headline":"A useful sensitivity study showing effective-mass definitions can shift NS cooling curves, but the central Dirac-vs-Landau comparison is under-derived: the Dirac-mass emissivity branch is likely unphysical.","tokens_in":19553,"tokens_out":4072,"would_cite":true,"duration_ms":38772,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["97.60.Jd","26.60.-c"],"model":"deepseek-v4-flash","headline":"The choice between the Dirac and Landau nucleon effective masses shifts neutron-star cooling predictions by 0.03–0.06 dex in surface temperature, an ambiguity most simulations ignore.","keywords":["neutron star cooling","effective nucleon mass","Dirac mass","Landau mass","relativistic mean-field theory","Walecka model","neutrino emissivity","direct Urca process"],"falsifier":"Derive the non-relativistic neutron-neutron bremsstrahlung emissivity from the Walecka Lagrangian to leading order in 1/m_ψ and read off whether the prefactor is m_D or m_L; if the prefactor is m_D, the faster-cooling Landau-mass curves in the paper are unphysical and the claimed systematic uncertainty collapses to whichever mass actually appears.","tokens_in":18523,"feed_emoji":"⭐","tokens_out":13972,"duration_ms":113469,"temperature":0.7,"pith_summary":"The paper argues that neutron-star cooling calculations carry a systematic uncertainty that most simulations ignore: which of two definitions of the nucleon effective mass is inserted into the neutrino-emissivity formulas. In a relativistic mean-field (Walecka) description, the Dirac mass — the in-medium mass reduced by the scalar meson field — falls with density, while the Landau mass — Fermi momentum divided by Fermi velocity — rises at high density. Feeding these two prescriptions into the standard emissivities for neutron-neutron bremsstrahlung and the direct Urca process changes the predicted surface temperature of a 1.4–2.0 solar-mass neutron star by roughly 0.03–0.06 dex, comparable to or larger than the observational error on well-studied cooling neutron stars. The paper checks this with two realistic equations of state that differ only in effective-mass properties and finds broadly the same behavior. The authors note that other cooling uncertainties may be larger, but effective masses are among the hardest to constrain experimentally, so the ambiguity should be folded into any cooling-based constraint on dense matter.","feed_headline":"Effective-mass choice shifts neutron-star cooling by 0.03-0.06 dex","feed_subtitle":"Most cooling simulations ignore this distinction; it rivals the observational error on well-studied neutron stars.","key_machinery":"The load-bearing identity is m_L = k_F / v_F = sqrt(k_F^2 + m_D^2), which ties the quasiparticle dispersion at the Fermi surface to the density of states D(ε) = k_F^2/(π^2 v_F) = 3 m_L n / k_F^2 and thus to the specific heat c_V = π^2 T m_L n / k_F^2. Because the same Landau mass controls the heat capacity on the left side of the cooling equation and the neutrino-emissivity prefactors on the right side (scaling as m*^4 for bremsstrahlung and as m*_N m*_p for direct Urca), swapping m_D and m_L changes both sides of the cooling equation at high density. The Dirac mass, by contrast, is the scalar-mean-field-modified mass term m_D = m_ψ − g_σ σ, which decreases with density. The paper pushes bot","core_discovery":"The central claim is that the Dirac mass and the Landau mass are physically distinct effective masses in relativistic mean-field models, and that the choice between them changes neutron-star cooling predictions at a level that matters for observations. In the Walecka model the Dirac mass m_D = m_ψ − g_σ σ falls with density, while the Landau mass m_L = (k_F^2 + m_D^2)^{1/2} = k_F/v_F rises at high density because the Fermi momentum grows. The standard emissivities for neutron-neutron bremsstrahlung and the direct Urca process, evaluated with one mass or the other, produce cooling curves that differ by roughly 0.03–0.06 dex in surface temperature for 1.4–2.0 solar-mass stars, with the Landau-","pith_inferences":["A first-principles reduction of the bremsstrahlung and direct-Urca matrix elements from the Walecka Lagrangian would likely single out one of the two masses; if it is uniquely the Landau mass, the paper's 0.03–0.06 dex band is an upper envelope of prescription choices, not a symmetric statistical error.","The same two-mass distinction should enter other Fermi-surface quantities — thermal conductivity, shear viscosity, and superfluid pairing gaps — so the effective-mass ambiguity may also affect neutron-star thermal relaxation, r-mode damping, and pulsar-timing glitch modeling, not just cooling curves.","Because the Dirac/Landau gap grows with density, a single precisely observed cooling sequence of a neutron star near the maximum mass would discriminate between the prescriptions more sharply than the current ensemble of lower-mass cooling stars."],"forward_implications":["Cooling curves computed with the Landau mass cool faster than those computed with the Dirac mass, and the gap widens as neutron-star mass increases.","The effective-mass ambiguity shifts predicted surface temperatures by about 0.03 dex for a 1.4-solar-mass star and about 0.06 dex for a 2.0-solar-mass star, comparable to or larger than the observational error on well-studied cooling neutron stars.","Below the direct-Urca threshold the effect still appears through the slow neutron-neutron bremsstrahlung channel, whose emissivity scales with the fourth power of the chosen effective mass.","The effect persists with realistic relativistic mean-field equations of state that differ only in effective-mass properties, so it is not an artifact of the simple Walecka model."],"fun_headline_variants":["Dirac vs Landau mass alters neutron-star cooling curves","Effective mass choice shifts neutron-star cooling by 0.03–0.06 dex","Most cooling simulations ignore a key nucleon mass distinction","Neutron-star cooling depends on which effective mass you pick","Dirac and Landau masses give different cooling rates for neutron stars"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The paper treats m_D and m_L as interchangeable plug-ins for m*_N in the non-relativistic emissivity formulas (Sec. III A) without deriving from the Lagrangian which one the phase-space factors actually require; if only one prescription is physical, the claimed 0.03–0.06 dex uncertainty is not a band but an error in one of the two curves.","fun_headline_variants_meta":{"raw":{"variants":["Dirac vs Landau mass alters neutron-star cooling curves","Effective mass choice shifts neutron-star cooling by 0.03–0.06 dex","Most cooling simulations ignore a key nucleon mass distinction","Neutron-star cooling depends on which effective mass you pick","Dirac and Landau masses give different cooling rates for neutron stars"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000885,"raw_usage":{"total_tokens":3678,"prompt_tokens":782,"completion_tokens":2896,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":526,"completion_tokens_details":{"reasoning_tokens":2807}},"tokens_in":526,"tokens_out":2896,"duration_ms":18993,"temperature":1.0,"reasoning_tokens":2807,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-31T23:57:21.345781+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Derive the non-relativistic neutron-neutron bremsstrahlung emissivity from the Walecka Lagrangian to leading order in 1/m_ψ and read off whether the prefactor is m_D or m_L; if the prefactor is m_D, the faster-cooling Landau-mass curves in the paper are unphysical and the claimed systematic uncertainty collapses to whichever mass actually appears.","supporting_citations":[],"review_version":1}