{"id":"e8038f27-0dd0-4609-b9b4-bcb842a982a8","arxiv_id":"2412.07758","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Hybrid stars with color-superconducting quark matter can reproduce the mass and spin of pulsar J0952-0607, whereas the hadronic-only model cannot, and the revised Kepler-frequency relation tightens radius bounds.","lead":"The paper models rapidly spinning neutron stars with a quark matter core and shows that such hybrid stars can explain the heaviest fast pulsar J0952-0607, while a purely hadronic model cannot. It also revises the empirical relation between spin frequency, mass and radius for hybrid stars, yielding new radius constraints.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"J0952-0607's stability in the favored hybrid model relies on an assumed constant T/W>0.08 f-mode threshold, not a calculation for a 2.35 solar-mass star.","rationale":"The reader's weakest_assumption focuses on the model family and parameter ranges of the quark-matter EoS. That concern is legitimate but somewhat generic: any study using a specific effective model is conditional on that model. The more directly load-bearing point is the paper's own treatment of non-axisymmetric stability, because the central claim requires J0952-0607 to be a stable equilibrium in the favored hybrid sequence. The paper explicitly assumes a mass-independent T/W=0.08 threshold for the f-mode instability. This assumption is used both to mark the stable region in Fig. 12 and to carve out the allowed parameter region in Fig. 13. If the threshold for a 2.35 solar-mass hybrid star is materially different, the conclusion that J0952-0607 is describable by the hybrid EoS would need revision, and the comparison with the hadronic DD2npY-T EoS would lose its force. The proposed check is concrete and directly tests the assumption. I do not see evidence of arithmetic or code-level errors, and the RNS framework is standard. The paper also gives credit for using a widely used code and for presenting a detailed parameter scan. However, the stability boundary is a physical calculation that should be verified for the specific high-mass configuration, not merely scaled from lower-mass results.","tokens_in":27237,"tokens_out":7300,"duration_ms":76258,"concrete_test":"Run a general-relativistic linear perturbation analysis on the RNS background for the rotating 2.35 solar-mass hybrid star with ηV=0.452, ηD=0.780 at 709.21 Hz, computing the l=m=2 f-mode complex frequency and the resulting instability threshold T/W_crit, including shear and bulk viscosity estimates. If T/W_crit is below the star's T/W from the equilibrium model, the central claim fails; if T/W_crit is close to 0.08, the existing assumption is validated.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's central claim that the favored hybrid set (ηV=0.452, ηD=0.780) can describe J0952-0607 requires that this star lies inside the non-axisymmetric stability region. The paper uses a single threshold T/W>0.08 for the onset of the l=m=2 f-mode instability, justified only by 'the weak dependence of the instability limit on the star's mass' (Sec. V). This threshold is assumed, not computed, for the 2.35 solar-mass hybrid configuration. Fig. 12 marks the vertical line T/W=0.08, and Fig. 13 uses it to exclude ηV<0.27 while leaving ηV=0.452 allowed. For a star near the maximum mass, the CFS instability boundary is set by the balance between gravitational-wave growth and viscous damping; both the mode eigenfunction and the damping timescale depend on the EoS and on the star's compactness, so the constant threshold is an extrapolation. If the actual threshold for this configuration is below its T/W, then J0952-0607 would be unstable and the favored hybrid sequence could not describe it, undermining the central comparison with the hadronic EoS.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies rapidly rotating hybrid neutron stars built from a hypernuclear relativistic density functional hadronic phase (DD2npY-T) and a color-superconducting confining RDF quark phase (ABPR parametrization), matched by a Maxwell construction. Using the RNS code, the authors compute rotating configurations, analyze the Kepler frequency and the empirical C(M) relation, and compare the fastest millisecond pulsars against stability regions in the mass–frequency plane. The central claim is that the heaviest and fastest galactic pulsar J0952-0607 (M=2.35±0.17 M⊙, f=709 Hz) can be described by hybrid stars with strong vector coupling ηV=0.452, while the purely hadronic hypernuclear DD2npY-T configuration cannot, and that this constitutes an indication for quark matter in neutron stars. The paper also quotes radius limits R1.4≤14.90 km and R0.7<11.49 km derived from the revised C(M) relation.","tokens_in":27443,"tokens_out":7865,"duration_ms":73692,"significance":"The numerical core of the paper is sound: the RNS code is a well-tested tool for rapidly rotating relativistic stars, and the EoSs are taken from prior work, so the reported sequences and stability maps are valuable. The study usefully extends the empirical C(M) relation to hybrid stars with an early deconfinement transition, and the proposed 1000 Hz pulsar forecast is a concrete, falsifiable prediction. The most striking physical conclusion, that J0952-0607 is incompatible with the specific hadronic DD2npY-T EoS but compatible with the favored hybrid set, is interesting and would be important if robust. However, the central comparison relies on a constant T/W=0.08 f-mode instability threshold that is assumed, not computed, for a 2.35 M⊙ near-maximum-mass star, and the comparison against hadronic matter is made with only one hadronic EoS. The radius limits are also calibrated on the same model family that they constrain. These issues are load-bearing for the paper's headline claims, although they appear addressable with additional calculations or a more careful framing.","major_comments":[{"comment":"The exclusion ηV≲0.27 and the viability of the favored ηV=0.452 set both rest on the assumption that the l=m=2 f-mode instability onset is T/W=0.08 for the 2.35 M⊙ J0952-0607 configuration, with the value taken from 1.4 M⊙ stars. The text justifies this by the 'weak dependence of the instability limit on the star's mass', but the CFS instability boundary depends on the mode eigenfunction and on the damping timescale, both of which are sensitive to compactness and to the EoS. No f-mode calculation, or even a bracketing estimate, is provided for the actual hybrid sequences. This is load-bearing: if the true threshold for the favored sequence lies below its T/W value, J0952-0607 would be unstable and the central comparison with the hadronic EoS would collapse. Please either compute the f-mode instability boundary for these configurations or demonstrate that a range of thresholds, including values below 0.08, leaves the conclusion unchanged.","section":"Section V, Figs. 12–13"},{"comment":"The logistic C(M) parametrization is fitted to the same hybrid EoS family that is then used to derive the upper radius limits quoted in the abstract (R1.4≤14.90 km, R0.7<11.49 km). The fit residuals are reported as 'maximal deviations ±8 Hz and relative error below 1%', but no propagation of this uncertainty into the radius limits is given, and the limits are presented as direct constraints. Because the relation is calibrated on the model set, these numbers should be presented as model-dependent estimates with propagated errors, not as independent bounds on R1.4 or R0.7.","section":"Section IV, Eq. (21)"},{"comment":"The statement that J0952-0607 is 'out of reach' for the purely hadronic hypernuclear configuration is based on a single hadronic EoS, DD2npY-T, and no maximum mass or maximum Kepler mass for this EoS is quoted in the paper. A stiffer hadronic EoS, or a different hyperon coupling, could in principle reach 2.35 M⊙ at 709 Hz, so the claim should be qualified as a statement about the DD2npY-T family, or supported by a scan over hadronic EoSs. Without that, the title's suggestion that MSPs are 'indicators for quark matter' goes beyond the presented evidence.","section":"Section V and abstract"}],"minor_comments":[{"comment":"Equation (3) contains a typo: the two derivatives in the expression for Δε are both written as dpq/dμ; the second should be dph/dμ.","section":"Sec. II C, Eq. (3)"},{"comment":"The statement that the logistic function has 'inflection point equal to E(Monset)·M K' is not correct as written; for the form in Eq. (21) the inflection point in M_K is (D+ln5)/E, so the sentence should be rewritten.","section":"Sec. IV, text after Eq. (21)"},{"comment":"The caption and text should clarify whether the cyan hatched region is computed for the fixed rest mass 2.1 M⊙ (gravitational mass ~1.9 M⊙) shown in the right panel or separately for a gravitational mass of 2.35 M⊙; as written, the panel label and the text are confusing.","section":"Fig. 13 caption"},{"comment":"The fitted oblateness coefficient 2√a≈0.772 is quoted without an uncertainty; the lower panels show deviations up to 6%, so an error bar on the fit value should be given.","section":"Sec. III B, Fig. 5"},{"comment":"A stray displayed formula 'dJ/dεc|Mb=const > 0' appears in the text after the paragraph discussing back-bending; this fragment should be removed or integrated into a proper sentence.","section":"Sec. II C, after the back-bending discussion"}],"recommendation":"major_revision","confidential_remarks":"The main risk is the T/W threshold assumption: if the authors can compute or conservatively bracket the f-mode instability for the near-maximum-mass hybrid configurations, the central claim becomes much stronger. The radius-limit claims are secondary and should be softened to model-dependent estimates. The paper is within the journal's scope and the numerical machinery is appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a solid, model-dependent study that deserves a serious referee. The genuinely new pieces are the logistic C(M) parametrization that folds in the deconfinement onset mass, and the use of the T/W instability to exclude weak vector couplings. The oblateness fit (2√a≈0.772) is a nice, clean result. The numerical work is standard and reliable: RNS code, EoSs from previous papers, and astrophysical constraints applied consistently.\n\nThe main soft spot is the non-axisymmetric stability criterion. The paper takes T/W=0.08 as the onset of the l=m=2 f-mode instability for all masses, justified by the weak mass dependence of the limit. For a 2.35 solar-mass star near the maximum mass, that is an extrapolation. The f-mode threshold does depend on compactness and EoS, and if the actual threshold for the favored hybrid configuration is below its T/W, J0952 would be unstable and the preferred parameter set would be excluded. This does not undo the hadronic-vs-hybrid comparison, which is driven by mass-shedding and the hadronic maximum mass, but it does make the specific exclusion ηV<0.27 and the stability of the favored set conditional on an uncomputed threshold. I'd ask the authors to either compute the threshold for their sequences or cite a dedicated study that justifies the extrapolation.\n\nSecond, the C(M) fit is self-referential in the mild sense: the logistic function is fitted to their own hybrid EoS family, and the same family is then used to derive radius upper limits. The limits (R1.4≤14.90 km, etc.) are consistent with earlier bounds, so no harm done, but they are quoted without propagated uncertainties and should be presented as conditional on this EoS family. Third, the data and scripts are not archived; 'available upon request' is not enough for a numerical paper of this type.\n\nReading it, I don't see a fatal flaw. The central argument—that a strong-vector-coupling hybrid EoS can accommodate J0952 while the DD2npY-T hadronic EoS cannot—is internally consistent. But it is a model-based existence proof, not a detection of quark matter. Other hadronic EoSs with higher maximum masses could conceivably do the job, so the conclusion is conditional on the EoS family.\n\nWho this is for: anyone working on rotating neutron stars, millisecond pulsars, or hybrid EoS. It deserves peer review; I'd recommend acceptance after revision, with the T/W threshold and uncertainty propagation addressed.","headline":"A solid, model-dependent study that adds a useful empirical relation and a plausible quark-matter interpretation of J0952-0607, but the f-mode threshold is assumed, not computed.","tokens_in":28059,"tokens_out":3882,"would_cite":true,"duration_ms":39648,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["85A15","83C05","76Y05"],"pacs":["97.60.Gb","26.60.Kp","12.38.Mh"],"model":"deepseek-v4-flash","headline":"Fastest-spinning pulsars may contain quark-matter cores, this paper argues.","keywords":["hybrid stars","millisecond pulsars","quark matter","color superconductivity","Kepler frequency","neutron star equation of state","J0952-0607","deconfinement phase transition"],"falsifier":"A decisive test is the discovery of a pulsar with a spin frequency clearly above 716 Hz: the paper predicts that a hadronic-only equation of state would be excluded for any such object at its measured mass, while the hybrid equations of state with strong vector coupling remain viable up to Kepler frequencies near 2.9 kHz. A clean radius measurement of J0952-0607 (for example by NICER) that is consistent with the hadronic mass-radius relation would similarly weaken the quark-core interpretation.","tokens_in":2256,"feed_emoji":"🌟","tokens_out":2221,"duration_ms":38401,"temperature":0.7,"pith_summary":"The paper studies rapidly rotating hybrid neutron stars, whose interiors contain a color-superconducting quark matter core beneath a hypernuclear hadronic outer core. Its central claim is that the fastest spinning and heaviest galactic millisecond pulsar, J0952-0607, cannot be reproduced by the purely hadronic equation of state considered here, but is naturally explained by hybrid star configurations when the quark matter has a sufficiently strong vector coupling. If correct, this provides an observational hint that deconfinement to quark matter occurs in neutron star cores, and it revises the empirical relation between spin frequency, mass, and radius that is used to constrain dense matter. The paper also derives new upper limits on neutron star radii from the fastest known pulsar, PSR J1748-2446ad, and from hypothetical 1000 Hz pulsars.","feed_headline":"Quark cores may explain the fastest spinning pulsar","feed_subtitle":"A hybrid star model reproduces J0952-0607's mass and spin where purely hadronic matter fails.","key_machinery":"The central object is the hybrid equation of state built by the Maxwell construction: a hadronic DD2npY-T equation of state for hypernuclear matter matched to the confining relativistic density functional (RDF) quark matter model with color superconductivity, parameterized by the ABPR form p = (A4 $μ^{4}$)/($2π^{2}$) + ($Δ^{2}$ $μ^{2}$)/$π^{2}$ − B, where A4, Δ, and B depend on the two free couplings ηV (vector) and ηD (diquark). The ratio of vector to scalar coupling sets the stiffness of quark matter and thereby the maximum hybrid star mass, while the diquark coupling sets the onset density of deconfinement. The rotating star configurations are computed with the RNS code, which solves the Einstein equations for axisymmetric uniformly rotating perfect-fluid stars, yielding the Kepler frequency, the onset of quasi-radial oscillations, and the T/W instability criterion.","core_discovery":"The paper argues that the observed mass and spin of the black-widow pulsar J0952-0607 (gravitational mass 2.35 ± 0.17 solar masses, spin 709 Hz) lie within the stable region for hybrid stars with a strong vector coupling (ηV = 0.452, ηD = 0.775 or 0.780), whereas for the purely hadronic DD2npY-T hypernuclear equation of state this pulsar lies beyond the stability limit at the 1σ level. Rotation shifts the onset of deconfinement to higher masses, and for these parameter sets the quark matter core extends enough to support the high mass at high spin. The same analysis maps the regions of angular velocity and mass where stars are hadronic, hybrid, or unstable, and shows that the gravitational-radiation-driven f-mode instability (T/W ≤ 0.08) excludes parameter values ηV ≲ 0.27 for the heaviest spinning stars. The paper also shows that the empirical C-factor in the Kepler-frequency relation fK = C(M/M⊙)^{1/2}(R/10 km)^{-3/2} rises from C ≈ 1.088 kHz for hadronic stars to C ≈ 1.16 kHz for hybrid stars with quark cores, and provides an analytic parametrization of C as a function of the quark onset mass.","pith_inferences":["If a pulsar spinning near or above 1000 Hz is ever found, the paper's framework predicts that a hadronic-only equation of state would be ruled out for it, making the search for such objects a direct test of quark matter in neutron stars.","The paper's C-factor parametrization could be tested against future NICER-style mass-radius measurements of rapidly rotating pulsars: a measured radius for a known spinning mass that violates the hadronic lower bound would independently support the hybrid interpretation.","The stability window analysis suggests a sharper population-level test: if several heavy millisecond pulsars are found above the T/W threshold for hadronic configurations but below it for hybrid ones, the statistical clustering of spins and masses would separate hadronic from hybrid equations of state.","The constraint ηV ≳ 0.27 from the f-mode instability effectively implies that any allowed quark matter description must be rather stiff, which could be compared with lattice or perturbative QCD estimates of the speed of sound at several times nuclear saturation density."],"forward_implications":["If J0952-0607 has a quark matter core, then the heaviest millisecond pulsars become direct evidence that deconfinement occurs inside neutron stars, not only in mergers or heavy-ion collisions.","The revised C-factor relation gives a simple tool to convert a measured pulsar spin frequency into a lower bound on the mass and radius of a non-rotating star, now applicable to hybrid stars with early deconfinement.","The radius upper limits R1.4 ≤ 14.90 km (from the 716 Hz pulsar J1748-2446ad) and R1.4 ≤ 11.90 km (for a hypothetical 1000 Hz pulsar) sharpen the constraints on the dense matter equation of state, with the 1000 Hz case strongly discriminating between hadronic and hybrid scenarios.","The T/W ≤ 0.08 stability window for J0952-0607 excludes vector couplings ηV ≲ 0.27 in this model family, providing a microphysical constraint on quark matter parameters.","The clustering of observed millisecond pulsars with similar spin frequencies finds a natural explanation in the spin-evolution model where the deconfinement transition slows down the spin-down rate, producing a waiting-time pile-up."],"supporting_citations":[{"why":"Supplies the ABPR parametrization of the RDF EoS, the mapping from coupling constants to A4, Δ, B, and the parameter ranges consistent with astrophysical constraints.","marker":"[47] Gärtlein et al. 2023"},{"why":"Defines the confining RDF quark matter model with color superconductivity and the medium-dependent vector and diquark couplings.","marker":"[46] Ivanytskyi & Blaschke 2022"},{"why":"Provides the ABPR parameterization that the authors fit, fixing the functional form of the quark matter EoS.","marker":"[48] Alford et al. 2005"},{"why":"Supplies the hadronic DD2npY-T equation of state with hyperons, the baseline against which hybrid stars are compared.","marker":"[49] Shahrbaf et al. 2022"},{"why":"Gives the mass measurement of J0952-0607 (2.35 ± 0.17 M⊙), the key observational anchor of the central claim.","marker":"[4] Romani et al. 2022"},{"why":"Provides the 716 Hz spin frequency of PSR J1748-2446ad, used to set the radius limits and the spin cutoff discussion.","marker":"[17] Hessels et al. 2006"},{"why":"Establishes the original empirical C = 1.08 kHz (hadronic) and C = 1.15 kHz (quark star) values that the paper revises for hybrid stars.","marker":"[86] Haensel et al. 2009"},{"why":"Provides the spin-evolution model with magnetic field decay and accretion that explains the pulsar clustering in frequency.","marker":"[94] Poghosyan et al. 2001"}],"fun_headline_variants":["Quark cores spin record pulsar J0952-0607","Hybrid stars with quark matter beat hadronic models","Fast spin points to quark cores in neutron stars","Rotation lifts quark matter cores to match pulsar","Quark phase transition explains fastest pulsar"],"cache_read_input_tokens":30080,"weakest_assumption_plain":"The central claim depends on the choice of quark matter model: if the real deconfined phase has a different stiffness or transition density than the RDF model with ABPR parametrization, then the purely hadronic equation of state might also accommodate J0952-0607, and the conclusion that this pulsar requires quark matter would no longer follow.","fun_headline_variants_meta":{"raw":{"variants":["Quark cores spin record pulsar J0952-0607","Hybrid stars with quark matter beat hadronic models","Fast spin points to quark cores in neutron stars","Rotation lifts quark matter cores to match pulsar","Quark phase transition explains fastest pulsar"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00059,"raw_usage":{"total_tokens":2892,"prompt_tokens":1190,"completion_tokens":1702,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":806,"completion_tokens_details":{"reasoning_tokens":1627}},"tokens_in":806,"tokens_out":1702,"duration_ms":10973,"temperature":1.0,"reasoning_tokens":1627,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T18:31:08.927495+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive test is the discovery of a pulsar with a spin frequency clearly above 716 Hz: the paper predicts that a hadronic-only equation of state would be excluded for any such object at its measured mass, while the hybrid equations of state with strong vector coupling remain viable up to Kepler frequencies near 2.9 kHz. A clean radius measurement of J0952-0607 (for example by NICER) that is consistent with the hadronic mass-radius relation would similarly weaken the quark-core interpretation.","supporting_citations":[],"review_version":1}