{"id":"cca2432a-eccb-48b9-89af-b635f62b96ce","arxiv_id":"2509.10589","paper_version":2,"verdict":"REJECT","confidence":"LOW","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":6,"one_line_summary":"A hybrid neutron-star model with a first-order hadron-quark phase transition predicts twin-star branches, reduced maximum masses, and 200-400 Hz post-merger gravitational-wave shifts, claimed to be detectable with next-generation detectors.","lead":"This preprint models neutron stars and their collisions with equations of state that include a sudden transition to quark matter, and lists several observable signals such as double-radius stars and shifted gravitational-wave frequencies. It claims these signals should be testable with next-generation detectors, but most quantitative results are asserted rather than derived.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Lattice QCD constraints are never actually coupled to the hybrid EoS construction; if the susceptibilities of Sec. II.A are varied, the TOV results likely don't change, removing the first-principles basis for the claimed signatures.","rationale":"I read the paper in good faith. The authors clearly state in Sec. IV.D that extrapolating lattice QCD from mu_B/T < 3 to neutron star conditions is model-dependent, which is honest. However, the more serious issue is structural: the paper never shows how the lattice data enter the EoS. The susceptibilities in Sec. II.A are computed and parameterized but are not used in Sec. II.B. This is not a matter of consensus; it is an internal gap between the stated framework and its implementation. If the coupling is absent, the signatures are not anchored in lattice QCD, and the central claim is overstated. A sensitivity test would settle this. If the authors can show a term like G_V(T,chi_2) or a matching condition, the concern is resolved. As submitted, the load-bearing numbers are unverifiable, so I do not object to the reader's REJECT, though the issue is fixable.","tokens_in":12872,"tokens_out":7534,"duration_ms":81157,"concrete_test":"Implement the EoS from Eqs. (7)-(11) as described and solve the TOV equations. Then repeat with the lattice susceptibility parameters A2,B2,... in Eqs. (3)-(5) shifted by their quoted 1-sigma errors (or, more directly, set the lattice input to zero) while holding the PNJL/NJL parameters fixed. If the mass-radius relation and maximum masses are unchanged to numerical precision, the lattice QCD input is not load-bearing and the abstract's 'constrained by lattice QCD' claim is not supported. Additionally, scan the manuscript for any term involving chi_n in the thermodynamic potential; absence confirms the gap.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The manuscript's central claim is that the hybrid EoS are 'constrained by lattice QCD' and therefore produce theoretically grounded signatures of QCD phase transitions. But Sec. II.A presents lattice susceptibilities chi_2, chi_4, chi_6 and a resummed pressure (Eq. 6) that never reappear in the EoS construction of Sec. II.B. The hadronic pressure uses chiral EFT and PNJL (with no lattice coupling); the quark phase is an NJL model fitted to vacuum pion properties (Sec. II.B.2); the Maxwell/Gibbs coexistence conditions (Eqs. 12-19) involve only the two pressures and chemical potentials. No equation links Eqs. (3)-(6) to Eqs. (7)-(11). The only mention is a qualitative sentence in Sec. II.B.1 that the PNJL model is 'constrained by the available lattice QCD data,' without any fitted parameters or matching conditions. Sec. IV.D further concedes that the finite-T lattice results require 'model-dependent extrapolations to neutron star conditions.' Thus the advertised first-principles anchor is not operational: the twin-star radii, maximum mass reductions, and post-merger frequency shifts are properties of the specific NJL/PNJL parameter sets, not of lattice QCD. Without demonstrating the coupling, the central claim collapses to 'a generic hybrid EoS produces generic signatures,' which is not what the abstract promises.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a framework for constructing hybrid hadron-quark equations of state (EoS) for neutron stars and neutron-star mergers, combining chiral effective field theory, PNJL/NJL quark models, Maxwell/Gibbs constructions, and lattice QCD data at finite mu_B. It then applies these EoS to TOV solutions and merger hydrodynamics, claiming signatures: twin-star branches with radius differences 0.5-2.0 km, reductions in M_max by 0.2-0.4 M_sun, post-merger gravitational-wave frequency shifts of 200-400 Hz, enhanced neutrino emission, and kilonova changes, all while remaining 'marginally consistent' with GW170817, NICER, and pQCD constraints. The abstract emphasizes that the EoS are 'constrained by lattice QCD.' The manuscript also discusses speed-of-sound behavior and universal-relation violations.","tokens_in":13318,"tokens_out":3802,"duration_ms":45939,"significance":"If the claimed signatures were derived from a demonstrably lattice-QCD-anchored EoS framework, the paper would be valuable for guiding next-generation detectors such as Einstein Telescope and Cosmic Explorer. The paper covers a broad and important set of observables (twin stars, post-merger GW f2 shifts, neutrino cooling, kilonova properties) and lists concrete, falsifiable predictions. It also correctly identifies limitations of current lattice and nuclear-matter calculations. However, the central advertised connection to lattice QCD is not made operational, and several quantitative predictions are asserted without the corresponding computational details. As written, the paper is better described as a phenomenological survey of hybrid-EoS effects than as a lattice-QCD-constrained derivation.","major_comments":[{"comment":"The lattice QCD input is never coupled to the EoS construction. Equations (3)-(6) give the lattice susceptibilities and the resummed pressure, but no equation in Sec. II.B uses them: the hadronic PNJL pressure and the quark NJL pressure are parameterized independently, and the coexistence conditions (12)-(19) involve only the two pressures and chemical potentials. The statement in Sec. II.B.1 that the PNJL model is 'constrained by the available lattice QCD data' is not supported by any fit, matching condition, or comparison. Consequently, the advertised first-principles anchor is not operational: the twin-star radii, mass reductions, and frequency shifts reported in Sec. III are properties of the specific NJL/PNJL parameter sets, not consequences of the lattice QCD analysis in Sec. II.A. This directly undermines the abstract's claim that the hybrid EoS are 'constrained by lattice QCD.'","section":"Sec. II.A vs. II.B"},{"comment":"The post-merger gravitational-wave results are reported without the computational details needed to assess them. No numerical relativity code, grid resolution, initial data, binary mass ratio, or equation-of-state table treatment is specified. The quoted frequencies f2 = 2150 ± 50 Hz and f2 = 1850 ± 80 Hz have no stated provenance for the central values or the error bars. Equation (28) is the quadrupole formula, which is not an adequate wave extraction for strong-field post-merger remnants. The detectability statement (SNR ≥ 15 at 100 Mpc with ET/CE) is unsupported: no noise curves, injection procedure, or detector-sensitivity calculation is presented. These are load-bearing for the paper's central claim (iii).","section":"Sec. III.C and Sec. II.D.2"},{"comment":"The reported uncertainties in the neutron-star structure results have no stated origin. For example, M_max(H) = 2.15 ± 0.08 M_sun and R = 13.2 ± 0.4 km are quoted in Sec. III.A without explaining whether these reflect EoS-parameter variations, lattice-fit errors, or some other source. The same applies to the twin-star radius differences. Additionally, Sec. IV.A concludes that the viable parameter space requires rho_trans ≲ 2.5 rho_sat, but the baseline model in Sec. III.A uses rho_trans = 2.8 rho_sat, and the main quoted signatures (Delta R = 1.2-2.0 km, M_max reduction) are for that baseline. The paper does not reconcile this tension, leaving unclear which numbers are actually consistent with the multimessenger constraints.","section":"Sec. III.A and Sec. IV.A"},{"comment":"The neutrino and kilonova predictions are asserted without the corresponding modeling details. Sec. III.D quotes neutrino luminosity enhancement factors of 2-5 and specific luminosities (e.g., 2 x 10^52 to 8 x 10^52 erg/s) but gives no neutrino transport scheme, no thermal evolution calculation, and no microphysical rates beyond a list of processes. Sec. III.E quotes ejecta masses, electron fractions, and light-curve shifts (0.3-0.8 magnitudes) without stating the ejecta model, the radiative-transfer solver, or the nuclear heating model. Since the abstract lists enhanced neutrino emission as a key signature, these omissions are load-bearing.","section":"Sec. III.D and Sec. III.E"},{"comment":"The paper itself concedes in Sec. IV.D that the lattice-QCD results require 'model-dependent extrapolations to neutron star conditions.' This admission is in tension with the abstract's strong claim that the framework is 'constrained by lattice QCD at finite temperature and baryon chemical potential up to mu_B/T < 3.' The extrapolation from mu_B/T < 3 to the cold, high-density regime relevant to neutron stars is a nontrivial step, and the paper does not demonstrate how the error bars or the qualitative conclusions change if that extrapolation is unreliable. At minimum, the claim should be softened to reflect the model-dependence the authors themselves identify.","section":"Sec. IV.D"}],"minor_comments":[{"comment":"The beta-equilibrium conditions include a neutrino chemical potential mu_nu_e. For cold neutron stars in beta equilibrium, neutrinos have left and mu_nu = 0; the presence of mu_nu_e is only relevant in trapped-neutrino or hot proto-neutron-star contexts. Please clarify which regime is being used.","section":"Sec. II.B.1, Eqs. (8)-(10)"},{"comment":"The parametrization in Eqs. (3)-(5) lists Gaussian fits with central values and errors, but it is unclear how these errors are propagated (or not) to any later result. A brief statement would help.","section":"Sec. II.A"},{"comment":"The phrase 'conformal bound violations' should be used carefully: c_s^2 = 1/3 is the conformal value, not a rigorous upper bound. The text later refers to c_s^2 > 1/3 as a violation, which is standard, but the initial wording could be misread.","section":"Sec. III.B"},{"comment":"The paper refers to figures (Figs. 1-6) with informative captions, but the actual figure panels are not discussed in the text with quantitative comparisons (e.g., how the hybrid curve in Fig. 4 is obtained from specific parameter choices). Please make the figures and their parameters explicit.","section":"General"},{"comment":"Some citations appear incomplete or imprecise: for example, Ref. [4] gives a journal page that may not correspond to the cited paper, and Ref. [16] has a typo in the quotation marks. A careful reference check is needed.","section":"References"}],"recommendation":"reject","confidential_remarks":"The manuscript has a broad scope and touches on many important observables, but the central claim of lattice-QCD-constrained hybrid EoS is not supported by any operational coupling between the lattice section and the EoS construction. The post-merger, neutrino, and kilonova results are presented without the numerical or modeling details needed for verification. These are not presentation-level issues; they affect the validity of the abstract's claims. In my view the paper would require major new content (actual matching to lattice data, full merger simulations, and uncertainty propagation) before it could be seriously considered. I therefore recommend reject rather than major_revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: the paper's headline quantitative claims are not verifiable as submitted. Post-merger frequencies, damping times, signal-to-noise ratios, neutrino luminosities, and kilonova properties are quoted without the simulations, error analysis, or parameter choices that produce them. The error bars on M_max and f2 have no stated provenance. Second, the stress-test is correct: the lattice QCD section is decorative. Equations (3)–(6) present susceptibilities and a resummed pressure, but nothing in the hybrid EoS construction (Section II.B) uses them. The hadronic phase is chiral EFT plus PNJL, the quark phase is NJL fitted to vacuum properties, and the Maxwell/Gibbs conditions only involve the two pressures and chemical potentials. No equation links the lattice inputs to the EoS. The claim of a lattice-constrained framework is therefore not operational.\n\nWhat the paper does well: it correctly identifies a real question, and it is transparent about its limitations. Section IV.D explicitly concedes that the finite-T lattice results require model-dependent extrapolations to neutron star conditions. The observation that strong first-order transitions are only marginally consistent with GW170817 and NICER if the onset density is pushed below ~2.5 rho_sat is worth stating, even if the paper gets there by moving the transition density after seeing the tension. The literature overview is serviceable, though the four headline signatures already appear in the cited literature: twin stars and frequency shifts in Refs. [39–41], detectability in Ref. [52], and merger evidence in Ref. [53]. The genuinely new part is narrow: a specific hybrid EoS family, but the paper does not provide enough detail to evaluate it.\n\nThe soft spots are the missing derivations and the post hoc tuning. These are fixable in principle: supply the parameter sets, the TOV and hydrodynamics code, the error analysis, and the simulation setup. As it stands, the load-bearing numbers are unverifiable, and the advertised first-principles foundation is not demonstrated.\n\nWho should read this: anyone wanting a compact survey of the multimessenger constraints and the standard hybrid EoS framework might get some value. A researcher looking for reliable predictions should not rely on it. It does not deserve a serious referee in its current form; I would desk reject with an invitation to resubmit with code, data, and derivations.","headline":"Unverifiable numbers and a decorative lattice QCD section sink an otherwise timely roadmap paper.","tokens_in":13769,"tokens_out":3364,"would_cite":false,"duration_ms":38341,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A first-order hadron-quark phase transition inside neutron stars would create twin-star branches, lower the maximum mass by 0.2–0.4 solar masses, and shift post-merger gravitational waves by 200–400 Hz, all within reach of next-generation d","keywords":["QCD phase transition","hybrid equation of state","neutron star","twin stars","gravitational waves","lattice QCD","quark matter","multimessenger constraints"],"falsifier":"Measure the radius of a 1.4 solar-mass neutron star to better than 0.5 km; a value above about 12.5 km would exclude the predicted hybrid branch at ~11.1 km. Alternatively, a high-signal post-merger gravitational-wave observation of a binary neutron star at ~100 Mpc that shows no 200–400 Hz drop in the dominant frequency within 20 ms after merger would rule out a strong first-order transition.","tokens_in":12774,"feed_emoji":"🌟","tokens_out":6604,"duration_ms":60453,"temperature":0.7,"pith_summary":"This paper argues that if a strong first-order phase transition from hadronic matter to deconfined quark matter occurs inside neutron stars, it will leave a distinct set of multimessenger signatures: twin-star radius gaps of 0.5–2 km, a 0.2–0.4 solar mass drop in maximum mass, a delayed 200–400 Hz shift in post-merger gravitational waves, and a factor 2–5 increase in neutrino luminosity. Building hybrid equations of state that join chiral effective field theory at nuclear densities to perturbative QCD at high density, with the intermediate region anchored to lattice QCD data at small baryon chemical potential, the paper finds these signatures are marginally consistent with current GW170817 and NICER constraints. If the central claim is right, next-generation detectors like the Einstein Telescope and Cosmic Explorer will be able to detect quark cores in merging neutron stars, turning QCD phase transition physics into an observable astrophysical subject.","feed_headline":"Quark transition leaves a 300 Hz fingerprint in neutron-star mergers","feed_subtitle":"Twin star gaps, neutrino glow, and a 300 Hz gravitational-wave drop become testable with next-gen detectors.","key_machinery":"The central object is the hybrid equation of state built by interpolating between three anchors: chiral effective field theory at densities up to about twice saturation, lattice QCD susceptibilities up to baryon chemical potential over temperature of 3, and perturbative QCD at asymptotic densities. The hadron-quark transition is implemented through Maxwell and Gibbs constructions with a parameter that interpolates between them. The key mechanism is the softening of the equation of state in the coexistence region: it creates a second stable branch of compact stars (twins), lowers the maximum mass, delays the post-merger gravitational-wave frequency drop, and boosts neutrino cooling.","core_discovery":"The central claim is that a strong first-order hadron-quark transition, implemented through Maxwell or Gibbs constructions in a hybrid equation of state, produces a specific, coherent set of observable signatures rather than a single anomaly. For a transition onset near 2.8–3.2 times nuclear saturation density, the maximum neutron star mass drops from about 2.15 to 2.05 solar masses, twin branches appear with radius differences up to 2 km, the post-merger gravitational-wave frequency shifts downward by 200–400 Hz roughly 10–20 ms after merger, and neutrino emission brightens by a factor of 2–5. The paper claims these signatures survive comparison with current observations—though only in a re","pith_inferences":["Beyond the paper: A confirmed 200–400 Hz post-merger frequency drop would not only indicate quark deconfinement but also constrain the location of the QCD critical point, effectively using astrophysical observations as a complementary probe to heavy-ion collision experiments.","Beyond the paper: The 'marginal consistency' with current data suggests the transition must occur at relatively low density; this narrows the allowed parameter space and could be sharpened with existing NICER data on PSR J0740+6620, potentially ruling out the hybrid branch without waiting for next-generation detectors.","Beyond the paper: The paper's thermal treatment in merger simulations is simplified; including non-equilibrium effects and magnetic fields could alter the predicted 10–20 ms delay time, which is a concrete extension to test.","Beyond the paper: A null detection by Einstein Telescope with a large sample of mergers would push the transition to be either weak, at higher density, or absent, effectively setting an upper bound on the latent heat and transition strength."],"forward_implications":["Twin star radius differences of 0.5–2.0 km are resolvable by next-generation X-ray timing missions and marginal with current NICER data, offering a direct test of the quark-core hypothesis.","Post-merger gravitational-wave frequency shifts of 200–400 Hz are detectable with signal-to-noise above 10 in Einstein Telescope and Cosmic Explorer for sources at ~100 Mpc, providing a clean observational signature.","Hybrid models predict 10–30% less dynamical ejecta and more neutron-rich ejecta, making kilonovae fainter and faster with enhanced production of heavy r-process elements such as rare-earth species.","Universal relations (I-Love-Q, binary Love relations) are violated by 3–15% when quark cores are present, giving an equation-of-state-independent diagnostic for phase transitions.","The speed of sound in neutron star cores would show c_s^2 below 0.5c^2 with transient conformal bound violations around 2–4 times saturation density, testable through future radius and tidal measurements."],"fun_headline_variants":["Neutron-star merger waves reveal quark transition: 200–400 Hz drop","Neutron star twins and 300 Hz drop: signs of quark matter","Next-gen detectors may catch quark matter via merger wave shifts","Neutrino glow and twin stars hint at QCD transition in neutron stars"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The predictions stand or fall with the assumption that lattice QCD results at baryon chemical potential up to about three times the temperature can be extrapolated by Taylor expansion and resummation to the cold, high-density conditions inside neutron stars—an extrapolation the paper concedes is model-dependent.","fun_headline_variants_meta":{"raw":{"variants":["Neutron-star merger waves reveal quark transition: 200–400 Hz drop","Neutron star twins and 300 Hz drop: signs of quark matter","Next-gen detectors may catch quark matter via merger wave shifts","Neutrino glow and twin stars hint at QCD transition in neutron stars"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000865,"raw_usage":{"total_tokens":3630,"prompt_tokens":830,"completion_tokens":2800,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":574,"completion_tokens_details":{"reasoning_tokens":2720}},"tokens_in":574,"tokens_out":2800,"duration_ms":22969,"temperature":1.0,"reasoning_tokens":2720,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T18:06:04.051322+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the radius of a 1.4 solar-mass neutron star to better than 0.5 km; a value above about 12.5 km would exclude the predicted hybrid branch at ~11.1 km. Alternatively, a high-signal post-merger gravitational-wave observation of a binary neutron star at ~100 Mpc that shows no 200–400 Hz drop in the dominant frequency within 20 ms after merger would rule out a strong first-order transition.","supporting_citations":[],"review_version":1}