{"id":"d9908be5-97dc-43cc-8c9e-50e6202db165","arxiv_id":"2502.05513","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":0.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A review of neutron star equations of state which finds hadronic and hard quark models still plausible, and soft NJL quark models disfavored.","lead":"This preprint reviews the equations of state used for neutron stars, comparing hadronic, hybrid, and quark matter models against pulsar, X-ray, and gravitational wave observations. It concludes that stiffer equations of state remain viable, while soft NJL-type quark models struggle to explain two-solar-mass neutron stars.","discovery_kind":"review","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Unsourced comparative tables are the load-bearing pillar of Section 7; if tabulated model values are misreported or unrepresentative, the conclusion that soft NJL models are disfavored while RMF/FRG/V-QCD remain viable loses support.","rationale":"The reader's weakest-assumption analysis identifies the compilation accuracy of Tables 1-4 as the critical premise for the review's conclusion. My independent reading agrees: Section 7's ranking of models as 'soft and inconsistent' versus 'stiff and viable' is exactly the content of those tables, and nothing else in the manuscript provides a quantitative derivation of those ranges. The tables lack citations, the extraction methodology is absent, and several quoted values are already internally problematic (e.g., the text admits NJL can be stiffened, while the table labels all NJL as < 2.0 M_sun). I therefore consider this the single most load-bearing concern. The proposed concrete test would settle it by checking whether the tabulated entries are traceable and accurate; until then, the review cannot be used as a reliable reference for the comparative claim. I do not see a separate concern that would change the verdict in a different direction, so the reader's CONDITIONAL verdict remains appropriate. The broken equations and reference errors are real but secondary: they weaken the review's pedagogical and technical value, not the comparative conclusion directly, unless they signal the same lack of quality control that makes the tables suspect.","tokens_in":13829,"tokens_out":2288,"duration_ms":25780,"concrete_test":"Construct a provenance table mapping every row of Tables 1-4 to a specific published EoS parameterization (e.g., QHC21 from Kojo et al. 2022; NJL with vector coupling G_V from Otto et al. 2020; GM1/DD2 RMF from Xia et al. 2022). Then recompute the Section 7 comparison using only rows whose source can be verified and whose quoted values match the source. If any unverifiable or misquoted row changes the classification of a model family as above or below the 2.0 M_sun threshold or above/below the GW170817 Lambda_1.4 constraint, the conclusion must be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim in Section 7 rests on the comparative ranking of EoS models presented in Tables 1-4 and Section 6.1. These tables contain no inline citations, and the manuscript does not describe how the ranges were extracted from the literature. Without source-to-entry traceability, the conclusion is unverifiable. This is not a cosmetic issue: the conclusion explicitly contrasts 'traditional NJL models' (M_max < 2.0 M_sun, Lambda_1.4 > 500) with stiffer models. Yet Section 3.3.1 itself notes that NJL models with vector repulsion or diquark pairing can stiffen the EoS and support > 2 M_sun. If the table entries omit these variants, the soft-vs-stiff dichotomy is a strawman rather than a fair comparison. The problem also appears elsewhere: Equation (3), the tidal deformability formula, is broken, and the hyperonic RMF Lambda_1.4 range (> 500) in Table 3 conflicts with many published hyperonic RMF models that include repulsive hyperon couplings and yield lower Lambda. No observation in the paper independently verifies these numbers, and the reference list contains duplicated and mislabeled entries (e.g., [14]'s page/article data appears garbled). The review's conclusion may be consistent with the broad field consensus, but the manuscript as written does not provide the evidence needed to support it.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript is a review of neutron-star equations of state, classifying models into hadronic (RMF, non-relativistic), hybrid (QHC, first-order), and quark-matter (NJL, FRG, V-QCD) categories, and comparing their predictions against mass, radius, tidal deformability, and nuclear-physics constraints. It summarizes historical developments, observational constraints from NICER, GW170817, pulsar timing, and future prospects with ET, CE, SKA, and heavy-ion facilities. The central conclusion (Section 7) is that models predicting a soft EoS, such as traditional NJL models, are largely disfavored by current observations, while stiffer models such as RMF, FRG, and V-QCD remain viable, and that hybrid first-order transitions may manifest through mass twins.","tokens_in":14067,"tokens_out":3148,"duration_ms":31331,"significance":"If the comparative tables and equations were properly sourced and corrected, this review could serve as a broad, accessible overview of the current EoS landscape for nonspecialists. The paper covers the major model classes and observables, and it identifies real tensions in the field (e.g., the hyperon puzzle, phase-transition signatures). The central conclusion is consistent with the broad community consensus, and the paper does not introduce circular reasoning or ad hoc inventions. However, the manuscript's load-bearing comparative analysis currently rests on unsourced tables and incomplete equations, so the paper as written does not itself provide the evidence needed to support its own conclusion. The review has potential but requires substantial technical repair before it can be considered reliable.","major_comments":[{"comment":"The Tolman-Oppenheimer-Volkoff equation is printed as '= −' with no right-hand side, so the central structure equation of neutron-star modeling is not actually stated. The standard dP/dr = −G(ε+P)(m+4πr^3P)/(r(r−2Gm)) form should be given with all variables defined.","section":"4.1, Eq. (1)"},{"comment":"The tidal deformability formula is printed as 'Λ = k2 5', which is dimensionally meaningless without the radius and mass dependence. The correct expression Λ = (2/3) k2 (R/M)^5 must be used, and the Love number k2 should be defined.","section":"4.2, Eq. (3)"},{"comment":"Several key model equations are incomplete fragments: Eq. (9) for RMF pressure is missing the baryon pressure sum and the correct meson-field coefficients, Eq. (15) defines PNJL only as a minus sign times an undefined sum, Eq. (16) is '= Tr' with no argument, and Eq. (17) is 'ϵ = T4' without the Nc and λ factors. A reader cannot follow the model predictions from these expressions; they must be written out correctly or removed.","section":"5.1, Eq. (9) and 5.3, Eqs. (15)-(17)"},{"comment":"The comparative tables are the empirical basis for the Section 7 conclusion, but they carry no inline citations and the manuscript does not describe how the model ranges were extracted or selected from the literature. This is not a cosmetic issue: Section 3.3.1 explicitly states that NJL models with vector repulsion or diquark pairing can stiffen the EoS and support masses above 2 M⊙, yet Table 1 lists 'NJL Model < 2.0' without distinguishing variants. Similarly, Table 3 assigns 'Hyperonic RMF > 500' for Λ1.4, which conflicts with published hyperonic RMF models that include repulsive hyperon couplings and yield lower tidal deformabilities. Each table entry needs a source or at least a clearly stated selection criterion; otherwise the soft-vs-stiff ranking in the conclusion is unverifiable.","section":"6.1, Tables 1-4"},{"comment":"The crossover and Gibbs-construction pressure formulas use weight functions w(ρ) and χ that are not defined. The text should specify the range and physical meaning of these volume fractions, and state how they are computed in the Maxwell versus Gibbs constructions.","section":"5.2, Eqs. (12) and (14)"}],"minor_comments":[{"comment":"The historical account of the Oppenheimer-Volkoff derivation is repeated almost verbatim in two consecutive paragraphs; the duplication should be removed.","section":"Section 2.1"},{"comment":"Several references are garbled or duplicated: [14] appears to share the title of [13] and is assigned an incorrect journal identifier (Nature Communications, 2019, 14: 8352 is the same as [13]); [24] duplicates [7] with 'Submitted'; and [18]'s volume/page data (Phys. Rev. D 101, 103006, 2020) matches [16], not the cited Read et al. article. The reference list needs systematic checking.","section":"References"},{"comment":"There are many typographical artifacts, including 'sufficient' (with a ligature), odd spacing in section headings such as 'M o d e l' and 'R e l a t i v i s t i c M e a n - F i e l d', and inconsistent dashes in ranges. These should be cleaned up in the revised version.","section":"Throughout"},{"comment":"The quoted pressure range P(2ρ0) ≈ 50–80 MeV/fm3 is stated without a citation; since this is an empirical constraint, a source should be provided.","section":"Section 4.3, Eq. (6)"},{"comment":"The NICER radius constraint '12.2 km ≤ R1.4 ≤ 13.7 km' is presented without specifying the posterior credible interval or the exact NICER analysis used; the paper should state the confidence level and cite the corresponding analysis.","section":"Section 4.1, Eq. (2)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript has the structure of a review but currently reads as a rough draft: several equations are empty fragments, the comparisons tables lack sources, and the reference list contains duplicated entries. The Section 7 conclusion is broadly consistent with community consensus, so the paper is salvageable, but the authors need to either properly source the comparative tables and correct the equations or substantially narrow the claims to what can be verified from the text. I would recommend a thorough technical revision before considering publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick verdict: this is a review, not a research paper, and the field already has better reviews. The bar for a new one is the quality of the synthesis, and on that bar this manuscript falls short in mechanics but not in coverage.\n\nWhat's actually good: the categorization into hadronic, hybrid, and quark-matter models is sensible, and the comparison tables, if they could be trusted, would give a newcomer a fast map of where current models sit against NICER and GW170817. The section on observational constraints covers the right measurements. The conclusion — that soft NJL models are disfavored while RMF/FRG/V-QCD stay viable — matches the broad consensus, and the review correctly notes that NJL with vector repulsion or diquark pairing can get above 2 solar masses, which is the main nuance.\n\nThe soft spots are substantial. Equations are broken: Eq. (1) is just '= −', Eq. (3) is 'Λ = k2 5' with the radius and mass factors missing, Eq. (9) has missing pieces, and Eqs. (15)-(17) are fragments. For a review whose purpose is to explain these relations to new readers, this is not cosmetic. The tables in Section 6 have no inline citations, so a reader cannot check where the numbers came from. That is the load-bearing part: the whole Section 7 ranking rests on those numbers. And the reference list has at least one duplicated/garbled entry ([14] looks like a bad copy of a Nature Communications paper). I also think the 'hyperonic RMF > 500' entry in Table 3 is suspiciously high unless they are cherry-picking a specific parameterization; many hyperonic RMF models with repulsive hyperon couplings give lower Lambda.\n\nNone of this means the central claim is wrong. The consensus is likely correct. But the manuscript as written does not supply the evidence needed to check.\n\nI'd send it back for major revision: fix the equations, add a source column or per-entry citations to the tables, and clean up the reference list. A well-executed version of this review would be worth reading, especially for students. As is, I wouldn't cite it, but it deserves a serious referee.","headline":"A well-structured but sloppy review whose load-bearing tables and equations are not yet trustworthy enough to support its (otherwise plausible) conclusion.","tokens_in":14587,"tokens_out":2141,"would_cite":false,"duration_ms":21627,"reading_group":"maybe","serious_thinker":"unclear","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This review argues that soft quark-matter equations of state are now ruled out by neutron-star observations, while stiffer models remain viable.","keywords":["neutron star equation of state","first-order phase transition","quark matter","hadronic matter","hybrid EoS models","tidal deformability","mass-radius relation","multi-messenger observations"],"falsifier":"Suppose a future X-ray timing measurement found the radius of a $1.4\\,M_\\odot$ neutron star to be 11.5 km with an uncertainty of 0.3 km. That would sit outside the bands quoted for RMF, FRG, and V-QCD (12.0 to 13.2 km) and inside the band quoted for NJL (11.0 to 11.8 km), directly contradicting the review's conclusion that soft quark models are ruled out. Conversely, a radius near 13.0 km at the same mass would confirm the stiff side.","tokens_in":13594,"feed_emoji":"⭐","tokens_out":8918,"duration_ms":85801,"temperature":0.7,"pith_summary":"This review surveys the competing equations of state proposed for neutron-star interiors and asks which ones survive the current astrophysical data. Its central conclusion, reached in Section 7, is that the soft end of the model space is effectively closed: traditional NJL quark-matter models predict maximum masses below two solar masses and tidal deformabilities above the range allowed by the 2017 binary neutron star merger, while stiffer frameworks such as RMF, FRG, and V-QCD remain consistent with the measured high-mass pulsars, X-ray radius bounds, and tidal deformability limits. The paper reaches this by sorting models into hadronic, hybrid, and pure quark categories and comparing their predicted maximum mass, radius, and tidal deformability against the same set of observational constraints. If the comparison holds, the open question is not whether any quark-matter phase exists, but which stiff description is right.","feed_headline":"Soft neutron-star equations of state ruled out, stiff ones survive","feed_subtitle":"Review finds the data favor RMF, FRG, and V-QCD models over traditional NJL quark matter.","key_machinery":"The central object is the equation of state itself, the relation between pressure and energy density inside a neutron star, because everything observable is computed from it through the Tolman-Oppenheimer-Volkoff equations and the tidal deformability $\\Lambda$, which measures how much the star deforms under its companion's gravity. The review's load-bearing machinery is its comparison tables: each model is reduced to three numbers, maximum mass, radius at $1.4\\,M_\\odot$, and tidal deformability $\\Lambda_{1.4}$, and those numbers are checked against the same observational windows. In addition, the distinction between a smooth hadron-quark crossover and a first-order Maxwell or Gibbs transition is what generates the mass-twin and post-merger-oscillation signatures the review highlights.","core_discovery":"In the paper's own framing, no equation of state is definitively confirmed, but the data have already removed a whole class. Models that produce a soft equation of state, exemplified by traditional NJL quark matter, cannot support the $\\ge 2\\,M_\\odot$ pulsars and predict $\\Lambda_{1.4} > 500$, in tension with the GW170817 constraint; models that produce a stiff equation of state, such as RMF, FRG, and V-QCD, support $>2\\,M_\\odot$ stars and give $\\Lambda_{1.4}$ in the observed range. Hybrid crossover models, notably QHC21, also remain viable, and first-order phase transitions are not excluded, with mass twins and post-merger oscillations named as the signatures that could confirm them. The conclusion is that continued observations will mainly discriminate among the surviving stiff models rather than between stiff and soft physics.","pith_inferences":["The review treats soft versus stiff as the main axis, but a model that is stiff at low and high density with a brief softening in between could evade the current windows; such a non-monotonic equation of state is not explicitly considered in the comparison.","Because the tables summarize other groups' published numbers, a re-derivation directly from the primary sources would be the fastest way to test the review's ranking; no new observation is needed.","If soft NJL models are truly dead, heavy-ion experiments that show a strong softening of nuclear matter at $2$ to $4$ times nuclear saturation density would force the field to revisit the conclusion, since terrestrial and astrophysical data would then point in opposite directions.","The paper's implied ordering could be sharpened into a statistical model-selection test: fit each surviving model to the same mass, radius, and tidal deformability data and compute evidence ratios; the review does not do this, but its tables provide the needed ranges."],"forward_implications":["If the conclusion holds, any viable quark-matter model must include strong repulsive interactions, such as vector couplings or FRG-type density-dependent couplings, to reach $2\\,M_\\odot$.","First-order phase transitions remain possible only if they occur at sufficiently high density, above about $2$ to $4$ times nuclear saturation, so that the softer branch does not violate the tidal deformability bound.","Mass twins become a sharp test: V-QCD and first-order hybrid models predict them, while RMF and FRG do not; a twin-star discovery would settle the nature of the transition.","The next generation of gravitational-wave detectors would be measuring the fine structure of the surviving stiff equations of state, not deciding between stiff and soft bulk behavior.","A radius measurement of a $1.4\\,M_\\odot$ star closer to 12 km would favor RMF-like models, while one closer to 13 km would favor V-QCD."],"supporting_citations":[],"fun_headline_variants":["Soft neutron-star EoS ruled out, stiff survive","Neutron-star data favor stiff EoS, kill soft","Data eliminate soft equations of state for neutron stars","Soft neutron-star EoS dead, stiff alive"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The weakest point is the tabulated comparison: Tables 1 through 4 list each model's maximum mass, radius, and tidal deformability without inline citations or a description of how the entries were extracted, so if even one entry is misreported the ranking of soft versus stiff models and the paper's conclusion would shift.","fun_headline_variants_meta":{"raw":{"variants":["Soft neutron-star EoS ruled out, stiff survive","Neutron-star data favor stiff EoS, kill soft","Data eliminate soft equations of state for neutron stars","Soft neutron-star EoS dead, stiff alive"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000802,"raw_usage":{"total_tokens":3539,"prompt_tokens":975,"completion_tokens":2564,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":591,"completion_tokens_details":{"reasoning_tokens":2501}},"tokens_in":591,"tokens_out":2564,"duration_ms":19586,"temperature":1.0,"reasoning_tokens":2501,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T19:01:17.607188+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Suppose a future X-ray timing measurement found the radius of a $1.4\\,M_\\odot$ neutron star to be 11.5 km with an uncertainty of 0.3 km. That would sit outside the bands quoted for RMF, FRG, and V-QCD (12.0 to 13.2 km) and inside the band quoted for NJL (11.0 to 11.8 km), directly contradicting the review's conclusion that soft quark models are ruled out. Conversely, a radius near 13.0 km at the same mass would confirm the stiff side.","supporting_citations":[],"review_version":1}