{"id":"6d77b44e-77c0-48f8-8390-3ae7bb354924","arxiv_id":"2509.03008","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":5,"one_line_summary":"The parity doublet model combined with the new NICER radius measurement restricts the chiral invariant nucleon mass m0 to 800-860 MeV, implying it is at least 85% of the nucleon mass.","lead":"Using the small newly measured radius of neutron star PSR J0614-3329, the authors tighten the allowed value of the chiral invariant part of the nucleon mass to 800-860 MeV. This suggests most of the nucleon mass may come from gluon condensation and other non-chiral mechanisms rather than from spontaneous chiral symmetry breaking.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Lower bound m0 >= 800 MeV hinges on 1-sigma acceptance rule; at 2-sigma, m0 = 700 MeV is admitted, so the claimed 220-MeV raise is not robust.","rationale":"The reader's verdict is CONDITIONAL and already lists the 1σ acceptance rule among the sensitivities in the rationale, so we agree with the overall conditional assessment. However, the reader's 'weakest_assumption' focuses on the identification of m0 with the chirally invariant mass and on the interpolation form. We identify a more immediate, falsifiable flaw: the core quantitative claim (raising the lower bound to 800 MeV) is an artifact of adopting 1σ credible intervals as hard exclusion limits. This is especially concerning because the paper's own results show m0 = 700 MeV is consistent with PSR J0614-3329 at 2σ, so the entire '220 MeV raise' vanishes if a more permissive (and arguably more standard) threshold is used. This is not a dispute about model interpretation; it is a methodological choice that can be settled by re-analysis, and it directly controls the headline statement. We therefore regard this as the single most load-bearing concern. The interpolation form and the m0 identification remain secondary issues that would also need discussion, but they are less immediately decisive because the paper's sensitivity tests partially address the former, and the latter is a matter of model interpretation rather than a numerical inconsistency. Hence, the paper should remain CONDITIONAL, not be rejected outright, because the authors could either justify the 1σ threshold with a full statistical treatment or soften the conclusion accordingly.","tokens_in":14330,"tokens_out":6906,"duration_ms":71913,"concrete_test":"Re-run the parameter-acceptance analysis using 2σ (or 90%) credible intervals for PSR J0614-3329 (R = 10.29^{+1.01}_{-0.86} km at M = 1.44^{+0.06}_{-0.07} M_sun) and for the other NICER/GW constraints, and check whether any m0 = 700 MeV EOS (with some (H,gV)) satisfies all constraints simultaneously. If it does—as the paper's own Sec. IIIB indicates—the claimed lower bound of 800 MeV is an artifact of the 1σ acceptance rule, and the headline conclusion must be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim—that PSR J0614-3329 raises the lower bound of the chiral invariant mass to 800 MeV and thus m0 ≥ 85% of the nucleon mass—is not derived from a physical principle but from the decision to accept only models that satisfy every observational constraint within 1σ (68%) credible intervals. In Sec. IIIB, the authors state that m0 = 700 MeV produces M–R curves consistent with GW170817, PSR J0030+0451, and PSR J0740+6620 at 1σ, and with PSR J0614-3329 only at 2σ; m0 = 800 MeV satisfies all at 1σ. Therefore, the entire exclusion of m0 = 700 MeV, and hence the 'approximately 220 MeV' raise in the lower bound, is determined by the arbitrary choice of a 68% threshold. If a 90% or 2σ standard is used—as is common in NICER analyses to guard against systematic errors—m0 = 700 MeV would be admitted, lowering the bound to at most 700 MeV and reducing the inferred fraction to ~74% of the nucleon mass. The paper gives no statistical justification for the 1σ rule, nor does it propagate the quoted uncertainties of the NICER contours into the m0 constraint. The 'dramatic refinement' therefore rests on a statistical convention rather than on the astrophysical data alone.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs unified neutron-star equations of state by combining a parity doublet model (PDM) for hadronic matter up to 2n0 with an NJL-type quark matter model beyond 5n0, connected by a smooth fifth-order polynomial crossover in the baryon chemical potential. The EOS is matched to nuclear saturation properties and the symmetry-energy slope, then used to solve the TOV equations. The authors scan the chiral invariant mass m0 from 600 to 900 MeV and vary quark-matter couplings H and gV, imposing causality, the existence of 2-solar-mass stars, and current NICER/GW170817 constraints. Their central claim is that including the new small-radius measurement of PSR J0614-3329 tightens the allowed range from 580 MeV ≲ m0 ≲ 860 MeV to 800 MeV ≲ m0 ≲ 860 MeV, implying m0 ≥ 85% of the nucleon mass and hence that non-chiral mechanisms dominate the origin of nucleon mass.","tokens_in":14736,"tokens_out":6041,"duration_ms":74652,"significance":"If the result holds, the paper would provide an astrophysical bridge between neutron-star radii and the QCD decomposition of the nucleon mass—an unusual and valuable connection. The internal chain (mean-field EOS, crossover construction, TOV integration, comparison with observational contours) is coherent, and the direction of the effect (smaller radius requires larger m0) is physically clear. The authors also test sensitivity to the interpolation boundaries nL and nU, and they explicitly enforce causality and thermodynamic stability. However, the main quantitative claim rests on a hard 1-sigma acceptance rule and on an identification between an effective-model parameter and a QCD mass decomposition; both are load-bearing and neither is satisfactorily justified. The paper is therefore potentially significant but requires substantial strengthening before the abstract-level conclusion can be accepted.","major_comments":[{"comment":"The exclusion of m0 = 700 MeV, and hence the new lower bound m0 ≥ 800 MeV, is an artifact of requiring consistency with PSR J0614-3329 at 1σ. The text states that m0 = 700 MeV is consistent with GW170817, J0030+0451, and J0740+6620 at 1σ and with J0614-3329 only at 2σ. Since the quoted NICER intervals are posterior credible intervals, using the 68% level as a hard accept/reject threshold without a likelihood or an explicitly justified decision rule is statistically arbitrary. Under a 90% or 2σ standard, m0 = 700 MeV would be admitted, lowering the bound to ≤700 MeV and reducing the inferred fraction from ≥85% to about 74%. The '220 MeV raise' and the '85%' claim are therefore not robust to the choice of confidence level. The authors should report the compatibility level of each m0 value, propagate the full NICER posterior contours, or otherwise justify the 1σ cutoff.","section":"Sec. III.B, Fig. 4"},{"comment":"The paper identifies the PDM parameter m0 with the chirally invariant component of the physical nucleon mass and concludes that gluon condensation and other non-chiral mechanisms dominate nucleon mass generation. This identification is an interpretation of an effective-model parameter, not a derived equality. The mass formula m*_αj = ... with a constant m0 in the parity-doublet Lagrangian is one possible parametrization; it need not coincide numerically with the chiral-invariant component defined in lattice QCD or QCD sum rules after mean-field and quantum corrections are absorbed. Thus the radius fit can be internally valid while the 'origin of nucleon mass' conclusion does not follow. The authors should either derive the connection from a QCD-based operator decomposition or explicitly frame the conclusion as model-dependent.","section":"Eq. (3) and Sec. IV"},{"comment":"The crossover region 2n0–5n0 is modeled with a fifth-order polynomial in μB matched through second derivatives. The paper tests sensitivity to the boundary densities nL and nU, but not to the functional form or matching order. Since PSR J0614-3329 constrains radii at densities that likely include the interpolation region, the assumed crossover form is load-bearing. Another interpolation scheme that satisfies the same boundary conditions but differs in the interior (e.g., a speed-of-sound interpolation or a different thermodynamic variable) could shift the M–R relation by an amount comparable to the 220-MeV effect claimed here. A systematic comparison of interpolation forms is needed before the tightened m0 range can be considered robust.","section":"Sec. III.A, Eqs. (25)-(28)"},{"comment":"The symmetry-energy slope is fixed to L = 57.7 MeV, the central value of the empirical constraint 57.7 ± 19 MeV. Because L controls the isospin-dependent pressure and hence the radius, and because the a0 meson and ω-ρ mixing are introduced specifically to reproduce L, the m0 constraint should be tested over the empirical range of L. The paper currently gives no indication whether the allowed range 800 ≲ m0 ≲ 860 MeV shifts when L = 38.7 MeV or L = 76.7 MeV is used. This uncertainty should be propagated into the final m0 bounds.","section":"Sec. II.A and Table III"}],"minor_comments":[{"comment":"The sentence 'This model was originally constructed to include strange quark effects through KMT-type interactions' appears twice in consecutive paragraphs; one occurrence should be removed.","section":"Sec. II.A"},{"comment":"Ref. [85] appears to duplicate Ref. [83] (same paper, one with arXiv number). Please consolidate.","section":"References"},{"comment":"The upper bound m0 ≲ 860 MeV is inherited from the previous analyses cited as Refs. [42,83], not derived from the new PSR J0614-3329 data. The phrase 'narrows this range to 800 ≲ m0 ≲ 860 MeV' could be misread as a new upper bound; clarify that only the lower bound is sharpened here.","section":"Sec. III.B"},{"comment":"The interpolation polynomial is written in powers of μB, but the text says 'fifth order polynomial of μB' and then sums from i=0 to 5; this is consistent but could be stated more explicitly as degree 5 to avoid confusion with a fifth-order expansion in another variable.","section":"Eq. (25)"}],"recommendation":"major_revision","confidential_remarks":"The 'previous constraint' 580 ≲ m0 ≲ 860 MeV is taken from the authors' own previous works, and the new lower bound is obtained within the same framework. An independent implementation of the PDM-plus-interpolation EOS would increase confidence. The statistical-threshold issue is the most serious: if the authors adopt a 90% or 2σ standard, the central abstract claim (m0 ≥ 800 MeV and hence m0 ≥ 85%) is not supported by the current analysis. The paper should be revised to either justify the 1σ rule with posterior information or soften the central claim accordingly."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the paper does a legitimate, timely thing — it runs the authors' parity-doublet plus quark-hadron crossover machinery against the new PSR J0614-3329 radius and finds that the small radius prefers large chiral invariant mass m0. The direction is plausible and the modeling is transparent. But the headline bound m0 ≥ 800 MeV (≥85% of the nucleon mass) rests on a 1-σ acceptance rule that the paper never justifies, and the paper itself reports m0=700 MeV passes all other constraints at 1-σ and fails only PSR J0614 at 2-σ. At a 2-σ standard the lower bound would be about 700 MeV, or roughly 74% of the nucleon mass. So the claimed 220 MeV raise is not robust.\n\nWhat's genuinely good: the EOS construction is clear, the constraint set is current (GW170817, NICER pulsars, 2-solar-mass requirement), and the authors check how the result changes with the interpolation boundaries. The causality and maximum-mass scans over (H,gV) are handled properly, and the comparison with the χEFT band is a useful sanity check. The paper also notes explicitly that m0=600 falls outside that band.\n\nSoft spots, roughly in order. First and biggest: the 1-σ acceptance rule. The difference between m0=700 and m0=800 is exactly a 1-σ vs 2-σ decision on one measurement. No statistical motivation is given, and common practice for NICER contours—which include systematic uncertainties—would be a 90% or 2-σ criterion. Since the central quantitative claim hangs on that choice, this is the load-bearing weakness. Second: the m0 grid is coarse (600, 700, 800, 900 MeV). The true lower bound could be anywhere between 700 and 800; the '220 MeV raise' is partly a grid artifact. Third: the identification of the effective Lagrangian parameter m0 with the chirally invariant mass of the physical nucleon is a model assumption, not a derived QCD statement. It's a reasonable interpretation within the PDM, but the paper's abstract and summary push it as a fact about QCD mass generation. Fourth, minor: L and other saturation properties are fixed at their central values; no uncertainty propagation.\n\nOverall, this is a serious paper with a real gap between its strongest claim and its evidence. It deserves peer review, but the authors should either adopt a defensible confidence criterion or present the constraint as a function of that criterion. I'd bring it to a reading group and would cite it with a caveat.","headline":"Legitimate model update using the new small-radius pulsar, but the m0≥800 MeV bound is cut at an unjustified 1-σ line; at 2-σ it drops to 700 MeV (~74%), so the headline is not robust.","tokens_in":15216,"tokens_out":3594,"would_cite":true,"duration_ms":38689,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The smallest reliably measured neutron star radius, about 10.3 km for a 1.44-solar-mass pulsar, forces the chiral-invariant share of the nucleon mass above 800 MeV — at least 85% of it — in a parity-doublet model with quark-hadron crossover","keywords":["nucleon mass origin","chiral invariant mass","parity doublet model","neutron star equation of state","PSR J0614-3329","NICER radius measurement","quark-hadron crossover","gluon condensation"],"falsifier":"A future high-precision radius measurement of a ~1.4-solar-mass neutron star would settle the central claim: a central value clearly above the J0614-3329 1σ band (above ~11.3 km) dissolves the pressure that raised the m0 lower bound, and the allowed window reverts toward 580–860 MeV; a value below ~10 km pushes the required m0 past 860 MeV and clashes with the two-solar-mass constraint, forcing changes in the quark-matter sector. Independently, a lattice QCD or QCD sum-rule computation of the baryon mass in the chirally restored phase below 800 MeV would falsify the nucleon-mass-origin interpr","tokens_in":14194,"feed_emoji":"⭐","tokens_out":21878,"duration_ms":181398,"temperature":0.7,"pith_summary":"This paper uses the recent NICER measurement of PSR J0614-3329 — the smallest reliably measured neutron star radius, about 10.3 km for a 1.44-solar-mass star — to weigh in on a question usually reserved for particle physics: where does the nucleon's mass come from? The authors build unified equations of state that stitch a parity doublet hadronic model, valid up to twice nuclear saturation density, to a quark model above five times that density through a smooth crossover; the chiral invariant mass m0 — the part of the nucleon mass that would remain if chiral symmetry were restored — is treated as a free parameter, and every current astrophysical constraint is applied. Older constraints allowed 580–860 MeV for m0; adding the small radius of PSR J0614-3329 shrinks the allowed window to 800–860 MeV. If the calculation is right, over 85% of the nucleon mass comes from chiral-invariant mechanisms such as gluon condensation rather than from spontaneous chiral symmetry breaking — a direct challenge to the traditional picture.","feed_headline":"Small neutron star radius: 85% of nucleon mass is chiral-invariant","feed_subtitle":"It means under 15% of a proton's mass comes from spontaneous chiral symmetry breaking.","key_machinery":"The central object is the chiral invariant mass m0 in the parity doublet model: the parameter that sets the mass of the nucleon doublet when the chiral condensate vanishes. In the model's mass formula, m* = (1/2)[sqrt((g1+g2)^2 (σ−ja)^2 + 4m0^2) ± (g1−g2)(σ−ja)], the condensate-dependent term dies as σ → 0 while m0 remains, so m0 is the natural candidate for the mass that would survive chiral restoration. The argument chains m0 to astrophysics: larger m0 weakens the σ coupling, softens the equation of state, and shrinks the predicted radius of a 1.4-solar-mass neutron star. Carrying the chain is a unified equation of state — parity doublet hadronic matter up to 2n0, a fifth-order polynomial","core_discovery":"The central claim is that the chirally invariant component of the nucleon mass, m0, must lie in the window 800 MeV ≲ m0 ≲ 860 MeV — at least 85% of the physical nucleon mass — once the small radius of PSR J0614-3329 joins the constraint set. In the parity doublet model, the nucleon mass splits into a chiral variant piece tied to the quark condensate and the parameter m0, which persists in the chirally restored phase. Larger m0 softens the equation of state and shrinks neutron stars, so the small measured radius rules out m0 ≲ 800 MeV; conversely m0 = 900 MeV makes stars too small for the two-solar-mass pulsar PSR J0740+6620. The authors take this as evidence that gluon condensation and other","pith_inferences":["The interpretation implies a concrete prediction: a second, independent small-radius measurement of a 1.4-solar-mass neutron star should land near 10.3 km. A central value above roughly 11.3 km would reopen m0 ≈ 700 MeV; one below roughly 10 km would push m0 past 860 MeV and force revisions in the quark-matter sector (diquark coupling or vector repulsion).","The fit constrains m0 only through the model's mass formula, so a direct lattice QCD or QCD sum-rule determination of the baryon mass in the chirally restored phase below 800 MeV would falsify the nucleon-mass-origin conclusion even while leaving the equation-of-state fit intact.","The paper tests the crossover boundaries but not the interpolation shape; repeating the analysis with a Maxwell construction (first-order transition) or with a family of alternate interpolation profiles would show how much of the 220 MeV shift is forced by the radius measurement and how much is carried by the assumed smooth crossover."],"forward_implications":["Nucleon mass is mostly not a product of spontaneous chiral symmetry breaking: at most about 15% of the 940 MeV can come from the quark condensate.","The allowed window 800–860 MeV sits in a narrow band where the equation of state is soft enough to give a 10.3 km radius yet stiff enough to still reach two solar masses.","Future radius measurements of 1.4-solar-mass neutron stars become direct probes of m0: the current constraint sits near the edge of the PSR J0614-3329 contour, so a slightly different central radius would move the window.","The lower bound is robust to where the hadron-quark crossover is placed: moving the matching density between 1.5n0 and 2.5n0 shifts the lower bound only between about 790 and 800 MeV, while the upper bound responds more strongly."],"supporting_citations":[{"why":"supplies the PSR J0614-3329 radius measurement (R = 10.29 km at 1.44 M⊙) that drives the tightened m0 window.","marker":"[63]"},{"why":"establishes the earlier 580 ≲ m0 ≲ 860 MeV constraint that the present analysis refines.","marker":"[42]"},{"why":"provides the previous parity-doublet matter analysis whose m0 range is being updated.","marker":"[83]"},{"why":"introduces the parity doublet model whose mass formula defines the chiral invariant mass m0.","marker":"[21]"},{"why":"supplies the low-density parity doublet hadronic equation of state used up to 2n0.","marker":"[66]"},{"why":"supplies the NJL-type quark model with color-flavor locked pairing used above 5n0.","marker":"[48]"},{"why":"provides the unified equation-of-state construction connecting hadronic and quark matter through crossover.","marker":"[89]"},{"why":"fixes the empirical symmetry energy slope L = 57.7 MeV used to pin the model's isovector parameters.","marker":"[67]"},{"why":"provides the two-solar-mass pulsar (PSR J0740+6620) measurement that excludes the largest m0 values.","marker":"[7]"}],"fun_headline_variants":["Tiny neutron star radius: 85% of proton mass is chiral-invariant","Neutron star radius: proton's chiral-invariant mass is ≥85%","Small pulsar radius reveals nucleon mass is 85% chiral-invariant","85% of nucleon mass is chiral-invariant, tiny star confirms"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The load-bearing premise is that the parameter m0 in the parity doublet mass formula equals the chirally invariant part of the physical nucleon mass — the mass that would survive in the chirally restored phase — so that fitting neutron star radii can be read as measuring the QCD origin of nucleon mass; equally load-bearing is the assumption that the smooth polynomial crossover between hadronic (2n0) and quark (5n0) matter represents the real transition region.","fun_headline_variants_meta":{"raw":{"variants":["Tiny neutron star radius: 85% of proton mass is chiral-invariant","Neutron star radius: proton's chiral-invariant mass is ≥85%","Small pulsar radius reveals nucleon mass is 85% chiral-invariant","85% of nucleon mass is chiral-invariant, tiny star confirms"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001489,"raw_usage":{"total_tokens":5899,"prompt_tokens":912,"completion_tokens":4987,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":656,"completion_tokens_details":{"reasoning_tokens":4915}},"tokens_in":656,"tokens_out":4987,"duration_ms":38790,"temperature":1.0,"reasoning_tokens":4915,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T11:13:18.519537+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A future high-precision radius measurement of a ~1.4-solar-mass neutron star would settle the central claim: a central value clearly above the J0614-3329 1σ band (above ~11.3 km) dissolves the pressure that raised the m0 lower bound, and the allowed window reverts toward 580–860 MeV; a value below ~10 km pushes the required m0 past 860 MeV and clashes with the two-solar-mass constraint, forcing changes in the quark-matter sector. Independently, a lattice QCD or QCD sum-rule computation of the baryon mass in the chirally restored phase below 800 MeV would falsify the nucleon-mass-origin interpr","supporting_citations":[{"cited_title":"Impacts of anomaly on nuclear and neutron star equation of state based on a parity doublet model","cited_arxiv_id":"2207.05970","evidence_quote":"supplies the low-density parity doublet hadronic equation of state used up to 2n0."}],"review_version":1}