{"id":"f0e79ca6-4215-40f3-9a33-6ef643c1fad5","arxiv_id":"2608.12653","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A modified Polyakov-loop model that stays finite at zero temperature produces hybrid neutron-star equations of state with quarkyonic or deconfined quark cores and maximum masses above 2 solar masses, for selected parameter sets.","lead":"This paper changes a quark-matter model so its 'freezing and melting' order parameter still works at zero temperature, and uses the result to build neutron-star interiors. With carefully chosen parameters, it finds stars with quark cores can weigh more than two suns.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Stability of the claimed stable hybrid-star branches is not established: the paper explicitly omits a radial oscillation analysis (Sec. III.A.1), and the M-R slope/turning-point criterion is not sufficient for stars with a sharp phase interface.","rationale":"The reader's verdict CONDITIONAL and rationale already note the missing radial oscillation analysis, but the reader's weakest_assumption focuses on the Polyakov-loop regulator in Eq. (20). I regard the stability gap as the single most load-bearing concern because the abstract's headline is 'Stable massive cold hybrid stars...' and the paper explicitly disclaims the analysis needed to support that word. The Polyakov regulator is a model choice whose physical status is debatable, but the paper frames its results as a qualitative model study; even if the regulator is purely phenomenological, the model's equations could still be internally consistent and produce the stated EOS. The stability issue, by contrast, is a direct gap between the claim and the evidence within the paper's own framework. It is also a concrete, checkable issue rather than an interpretive dispute. I therefore agree with the CONDITIONAL verdict (no change) but partially disagree with the reader's identification of the weakest assumption; the stability analysis is the more decisive missing piece.","tokens_in":21494,"tokens_out":6774,"duration_ms":68945,"concrete_test":"For the maximum-mass hybrid-star configurations of Figs. 8c (deconfined core), 11c (quarkyonic core), and 15c (three-phase), solve the Chandrasekhar radial oscillation equations with the correct junction conditions at the phase interfaces (continuity of the Lagrangian displacement and the Eulerian pressure perturbation across each sharp interface). Compute the squared fundamental eigenfrequency ω² and the first interface-coupled mode. If any ω² < 0 for a configuration claimed stable, the central claim fails; if all are positive, the stability concern is resolved and the conditional verdict can be upgraded.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that stable massive hybrid stars (M_max > 2 M_sun) exist with quarkyonic and/or deconfined quark cores. In the model realizations of Sec. III.B, stability is inferred only from the positive slope of the mass-radius branch (e.g., Figs. 8c, 11c, 15c). The turning-point criterion dM/dε_c = 0 locates where an eigenmode changes sign for ordinary one-phase stars, but for a star containing a first-order interface (the Maxwell-constructed hadron-quark discontinuity and the quarkyonic-deconfined discontinuity), stability must be checked with the radial oscillation equations or an equivalent interface-stability condition; a locally positive dM/dR branch can still have an unstable interface mode. The paper explicitly acknowledges the gap: 'However, we emphasize that a more detailed radial oscillations analysis is not performed in this study' (Sec. III.A.1, after Fig. 1c). This admission covers all the claimed stable configurations, since the method used to classify branches as stable/unstable is the same throughout. If the fundamental radial mode or the interface-coupled mode were negative for the maximum-mass configurations in Figs. 8c, 11c, and 15c, the headline claim of 'stable massive cold hybrid stars' would be false. The claim is therefore currently unverified in exactly the part that the abstract emphasizes.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript introduces a modified Polyakov-loop NJL (mPNJL) model in which the Polyakov potential is made finite at zero temperature by replacing T^4 with the quark-gluon Stefan-Boltzmann pressure and by substituting T0/T with T0/sqrt(T^2 + g(μ_f)) (Eqs. (16) and (20)). The quark sector is matched to several hadronic EOSs (SFHo, DD2, DD2hyp, NL3ωρ) through Maxwell constructions, and the resulting EOSs are integrated in the TOV equations. The authors systematically vary T0, α0, η2, GV, Gvv, and a bag constant B0, and identify parameter regions in which hybrid stars with quarkyonic (Φ=0) and/or deconfined (Φ>0) quark cores reach M_max > 2M_sun. The paper states that no radial-oscillation analysis is performed and that stability conclusions are based on the behavior of the mass-radius branches.","tokens_in":21866,"tokens_out":9119,"duration_ms":91290,"significance":"If the stability claim were fully established, the paper would provide a useful qualitative map of how Polyakov-potential parameters, vector couplings, and a bag constant control the cold dense-matter phase structure and hybrid-star properties. The study is transparent: the parameter scans are clearly described, the Polyakov-loop equation has analytical solutions, and the use of several hadronic EOSs and beta-equilibrated matter makes the setup concrete. Its main contribution would be the explicit realization of a T=0 deconfinement transition with a quarkyonic window and the identification of repulsive vector interactions as necessary for stiff quark cores. However, the headline statement that 'stable massive cold hybrid stars' are obtained is currently not supported by the analysis presented, because stability is inferred from M-R branch morphology rather than from a radial-mode or interface-stability calculation. The paper is therefore best viewed, at this stage, as a model-building and parameter-dependence study with suggestive but unverified astrophysical conclusions.","major_comments":[{"comment":"The paper explicitly concedes in Sec. III.A.1, immediately after Fig. 1c, that 'a more detailed radial oscillations analysis is not performed in this study.' Despite this, Sec. III.B repeatedly refers to configurations as stable (e.g., 'the stable maximum mass hybrid star' for the DD2 case of Fig. 8c, the quarkyonic-core stars of Fig. 11c, and the three-phase star of Fig. 15c), and the abstract asserts that 'Stable massive cold hybrid stars ... are obtained.' The stability classification appears to rely entirely on the shape of the mass-radius branch. For EOSs with Maxwell-constructed first-order phase transitions, the standard turning-point (dM/dε_c = 0) argument is not sufficient: the fundamental radial mode can change sign before the turning point, and interface or two-phase eigenmodes can be unstable even when the one-phase M-R slope is positive. Since this stability check is absent, the central claim is unverified. I recommend adding a radial-oscillation analysis with proper junction conditions at the hadron-quark and quarkyonic-deconfined interfaces, or explicitly revising all stability claims in the abstract and conclusions to state that only branch monotonicity has been checked.","section":"Sec. III.A.1 and Sec. III.B (Figs. 1c, 8c, 11c, 15c)"},{"comment":"The replacement T0/T -> T0/sqrt(T^2 + η2 μ_f^2) is introduced as a regulator for the divergent expression in Eq. (17) at T=0, and η2 is fixed by matching the Taylor expansion of Eqs. (17) and (20) to second order at T_dec^lat = 170 MeV. This is an ad hoc construction, and the paper offers no independent evidence that the traced Polyakov loop retains its confinement-deconfinement interpretation at zero temperature after this substitution. The existence of the quarkyonic (Φ=0) branch and its Maxwell transition to the deconfined (Φ>0) branch—the basis for the quarkyonic-core stars in Sec. III.B.2—is directly determined by this regulator. Because the series in Eq. (21) is truncated at quadratic order, the phase structure could depend on the omitted higher-order terms. I suggest two concrete checks: (i) extend g(μ_f) to higher orders and verify that the two-transition structure and the M-R stability classification are robust; (ii) compare the predicted Φ(μ_B) and transition chemical potentials with available functional-QCD or lattice-based constraints at finite density. Without such checks, the quarkyonic phase may be an artifact of the truncation.","section":"Sec. II.B, Eqs. (20)-(23)"},{"comment":"The concluding section states that, for the NL3ωρ case, 'we find a genuine, mechanically stable three-phase structure,' in which the maximum-mass configuration contains a deconfined core, a quarkyonic shell, and a hadronic envelope. This statement is not supported by the analysis in Sec. III.B.3: the three-phase star contains two first-order interfaces, and the maximum-mass configuration in Fig. 15c is again classified as stable solely from the slope of the M-R curve. Interface-coupled modes can be unstable even when the central density is below the one-phase turning point, so the term 'mechanically stable' overstates what has been demonstrated. This is the same type of gap as Major Comment 1, but it is worth flagging separately because the three-phase configuration is the most novel result of the paper.","section":"Sec. III.B.3 and Conclusion (Fig. 15)"}],"minor_comments":[{"comment":"The eight-quark interaction is written in Eq. (1) with \\bar{\\psi}\\gamma^\\mu\\lambda^0\\psi while the rest of the Lagrangian uses q fields; the mean-field reduction leading to the 4/3 G_vv (Σ_f ρ_vf)^4 term in Eq. (4) should be shown explicitly or the notation should be unified.","section":"Eqs. (1), (4), (30)"},{"comment":"The text should state the mass dimension of g(μ_f) and verify that η2 as defined in Eq. (23) indeed gives g units of energy squared; this is easy to check but is not currently spelled out.","section":"Sec. II.B, Eq. (23)"},{"comment":"The legends use 'Hadronic, Quarkyonic(Φ=0), Deconfined Quark(Φ≠0)' with color coding, but many curves are difficult to distinguish in printed grayscale; distinct line styles would help the reader follow the parameter shifts.","section":"Figs. 1-6"},{"comment":"The caveat that radial-oscillation stability has not been computed appears only in Sec. III.A.1 and is not restated in Sec. III.B or the abstract; given the strong wording of the abstract, the caveat should be prominently repeated wherever the word 'stable' is used.","section":"Sec. III.A.1 and Abstract"},{"comment":"The series for g(μ_f) starts at n=1, so odd powers of μ_f are formally present; the paper assumes their coefficients vanish by symmetry, but μ_f is not symmetric about zero and this assumption should be justified explicitly.","section":"Sec. II.B, Eq. (21)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a competent and readable parameter study, but the abstract's central stability claim is stronger than the analysis supports. The missing radial-oscillation analysis is a fixable gap, so I do not recommend rejection; however, the authors should be required either to add the stability calculation or to downgrade all 'stable' claims to 'configurations lying on rising M-R branches that have not been tested for radial stability.' The regulator dependence of the quarkyonic phase (Eq. (20)) is a second concern that should be addressed with sensitivity checks before the model is presented as a physical description of cold dense matter."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the paper is a transparent parameter map for a modified PNJL model, and that map is useful. But the abstract's headline claim—stable massive cold hybrid stars—is not established, because stability is read off the slope of mass-radius curves and the authors explicitly say (Sec. III.A.1) that no radial-oscillation analysis is done. For ordinary one-phase stars the turning-point criterion is fine; for stars with a Maxwell-constructed hadron-quark interface it is not sufficient, since an interface mode can be unstable even when the branch slope is positive. That caveat covers all the 'stable' branches later in the paper, because the same classification is used throughout.\n\nWhat is actually new: the specific mPNJL potential with T0→T0/sqrt(T²+η₂ μ²) (Eq. 20), combined with the 8-quark vector interaction, beta-equilibrated quark matter, and double Maxwell constructions. The thermodynamics looks internally consistent, and the paper is honest about what comes from Refs. [73,74] and what is new. The citation pattern looks fair; earlier PNJL0 papers are cited and the difference is stated. The scans of T0, α0, η2, GV, Gvv, and B0 give a clean qualitative picture—for instance, repulsive vector interactions are what allow a quark core without immediately destroying the star. The authors also state plainly that with the soft SFHo EOS no quark-core configuration reaches 2 M_sun.\n\nThe soft spots are, in proportion: (1) The T=0 Polyakov loop is used as a deconfinement order parameter without independent support; Eq. (20) is a regulator chosen for finiteness, and calling the Φ=0 branch 'quarkyonic' is an interpretation, not a derivation. (2) B0 is introduced ad hoc to shift the transition, so the >2 M_sun examples are properties of chosen parameters, not predictions. The paper says it is a qualitative study, and that is the right framing. (3) The speed-of-sound comments are fine as observations, not constraints.\n\nWho is this for? People building hybrid-star EOS libraries and planning Bayesian fits; they will find the parameter dependencies useful. It deserves peer review, because the framework is coherent and the gap is fixable in revision—the referee should ask for a radial-oscillation analysis of the maximum-mass configurations, or at minimum a clear softening of the abstract and the word 'stable.' My recommendation: send it to review. With the stability caveat resolved, it would be a solid contribution.","headline":"Useful parameter map for a modified PNJL hybrid-star EOS, but the 'stable massive hybrid stars' headline is not supported until a radial-oscillation analysis is done.","tokens_in":22372,"tokens_out":3899,"would_cite":false,"duration_ms":43175,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper proposes a modified Polyakov-Nambu-Jona-Lasinio model that keeps the deconfinement order parameter alive at zero temperature, and shows that the resulting hybrid-star equations of state can support quark cores above two solar…","keywords":["modified Polyakov-Nambu-Jona-Lasinio model","hybrid stars","quarkyonic matter","confinement-deconfinement transition","neutron star equation of state","Polyakov loop","Maxwell construction","vector quark interactions"],"falsifier":"A nonperturbative QCD calculation of the Polyakov loop at zero temperature as a function of baryon chemical potential, or an imaginary-chemical-potential lattice measurement, could test whether a sharp jump from zero to nonzero actually occurs at the model's quarkyonic-to-deconfined transition; if no such jump exists, the transition is a regulator artifact. A purely empirical cross-check is that if precision mass-radius and gravitational-wave data exclude the stiff low-density equations of state used here, the predicted above-two-solar-mass quark cores cannot exist.","tokens_in":1899,"feed_emoji":"⭐","tokens_out":2514,"duration_ms":83939,"temperature":0.7,"pith_summary":"The paper proposes a version of the Polyakov-Nambu-Jona-Lasinio quark model in which the Polyakov loop potential depends explicitly on the quark chemical potential, so that the confinement-deconfinement order parameter does not vanish at zero temperature. Combining this quark sector with hadronic equations of state through a Maxwell construction, the authors find cold neutron stars can contain a core of confined (quarkyonic) quark matter, deconfined quark matter, or both. They report that such hybrid stars can reach maximum masses above two solar masses, provided the low-density hadronic equation of state is stiff and repulsive vector interactions among quarks are strong. The paper also maps how each model parameter shifts the phase transitions, which matters for interpreting future mass-radius and gravitational-wave constraints on dense matter.","feed_headline":"Cold hybrid stars with quark cores can exceed two solar masses","feed_subtitle":"A chemical-potential-aware Polyakov loop keeps deconfinement alive at zero temperature and enables massive hybrid stars.","key_machinery":"The central object is the modified Polyakov loop potential $U(\\Phi,T,\\mu_f)$, obtained by replacing $T^4$ with the Stefan-Boltzmann pressure on the left of the standard potential and by substituting $T_0/T \\to T_0/\\sqrt{T^2+\\eta_2\\mu_f^2}$ on the right. The coefficient $\\eta_2$ is fixed by matching the low-density expansion at the lattice deconfinement temperature, so the potential remains finite at $T=0$ and the coefficients $a(T,\\mu_f)$ and $b(T,\\mu_f)$ carry the chemical-potential dependence. This object supplies the confinement-deconfinement transition in cold matter: its stationary condition gives analytic branches for $\\Phi$, and Maxwell constructions using those branches produce the hadronic, quarkyonic, and deconfined segments of the hybrid-star equation of state that are then integrated in the Tolman-Oppenheimer-Volkoff equations.","core_discovery":"The central claim is that a modified Polyakov-loop potential, built by replacing $T_0/T$ with $T_0/\\sqrt{T^2+\\eta_2\\mu_f^2}$, remains finite at $T=0$ and gives the Polyakov loop $\\Phi$ a nontrivial, chemical-potential-driven behavior in cold dense matter. In $\\beta$-equilibrated, charge-neutral quark matter, minimizing the thermodynamic potential yields analytic branches $\\Phi=0$ (confined quarkyonic matter), $\\Phi=1$, and $0<\\Phi<1$ (deconfined quark matter). Maxwell-constructed hybrid equations of state built from this quark sector and the hadronic SFHo, DD2, DD2hyp, and NL3$\\omega\\rho$ models can produce stable hybrid stars with a deconfined quark core, a quarkyonic core, or a three-phase deconfined-quarkyonic-hadronic structure, with maximum masses above $2M_\\odot$ in the stiff-equation-of-state cases. The paper finds that repulsive vector interactions, $G_V$ and $G_{vv}$, are essential for a stable quark core, and that in quarkyonic-core maximum-mass stars the central squared speed of sound exceeds the conformal value $c_s^2 = 1/3$.","pith_inferences":["An inference the paper leaves implicit: if the zero-temperature Polyakov loop is not a genuine order parameter for deconfinement, the quarkyonic-versus-deconfined distinction reduces to a parameter choice of the regulator.","A natural next check, not performed here, is a radial-oscillation analysis of the three-phase branch, since the paper's stability criterion is the slope of the mass-radius curve.","A testable extension: the Maxwell flat segments in the equation of state should produce characteristic plateaus in mass-radius curves that future precision mass-radius and tidal-deformability measurements could distinguish."],"forward_implications":["Cold neutron-star interiors can plausibly contain quarkyonic or deconfined quark matter without requiring high temperatures.","Maximum masses above two solar masses are achievable with quark cores only when the low-density hadronic equation of state is stiff and quark vector repulsion is sufficiently strong.","Quarkyonic-core maximum-mass stars require the central speed of sound to exceed the conformal limit, whereas deconfined-core stars can remain below it.","Depending on the parameters, the model predicts either one or two phase transitions, leading to qualitatively different mass-radius and tidal-deformability signatures.","With the soft SFHo hadronic equation of state, all quark-core configurations stay below two solar masses, so the existence of massive hybrid stars is tied to low-density stiffness."],"supporting_citations":[{"why":"Supplies the method of making the Polyakov potential chemical-potential dependent, including the regularization that keeps it finite at zero temperature.","marker":"[73]"},{"why":"Provides the analytical solution for the Polyakov loop branches used to identify confined and deconfined phases.","marker":"[74]"},{"why":"Motivates replacing $T^4$ with the Stefan-Boltzmann pressure through a nonperturbative analysis of the Polyakov potential.","marker":"[69]"},{"why":"Provides the perturbative $T_0(\\mu_f)/T$ substitution whose low-temperature divergence the model regularizes.","marker":"[70]"},{"why":"Gives the standard Polyakov potential parametrization and lattice-fitted coefficients that the modified potential generalizes.","marker":"[43]"},{"why":"Supplies the lattice deconfinement temperature used to fix the coefficient $\\eta_2$.","marker":"[75]"},{"why":"Introduces the eight-quark vector interaction whose cubic density dependence stiffens the quark equation of state.","marker":"[54]"},{"why":"Provides the SFHo hadronic equation of state used as the soft low-density baseline.","marker":"[83]"},{"why":"Provides the DD2 hadronic equation of state whose stiffness allows hybrid stars to exceed two solar masses.","marker":"[84]"},{"why":"Supplies the Maxwell construction procedure used to join the hadronic and quark phases.","marker":"[87]"}],"fun_headline_variants":["Quarkyonic cores push hybrid stars past 2 solar masses","Modified Polyakov loop stabilizes heavy hybrid stars","Cold dense stars with quark cores top two solar masses","Hybrid stars with quarkyonic matter surpass 2 solar masses","Chemical-potential-aware Polyakov loop yields massive hybrid stars"],"cache_read_input_tokens":24448,"weakest_assumption_plain":"The load-bearing premise is that replacing $T_0/T$ with $T_0/\\sqrt{T^2+\\eta_2\\mu^2}$ keeps the Polyakov loop a genuine order parameter for deconfinement at zero temperature; the paper offers no independent evidence for that identification, and if it fails the quarkyonic-versus-deconfined distinction is an artifact of the regulator.","fun_headline_variants_meta":{"raw":{"variants":["Quarkyonic cores push hybrid stars past 2 solar masses","Modified Polyakov loop stabilizes heavy hybrid stars","Cold dense stars with quark cores top two solar masses","Hybrid stars with quarkyonic matter surpass 2 solar masses","Chemical-potential-aware Polyakov loop yields massive hybrid stars"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001241,"raw_usage":{"total_tokens":5164,"prompt_tokens":1090,"completion_tokens":4074,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":706,"completion_tokens_details":{"reasoning_tokens":3992}},"tokens_in":706,"tokens_out":4074,"duration_ms":28061,"temperature":1.0,"reasoning_tokens":3992,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:03:38.214575+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A nonperturbative QCD calculation of the Polyakov loop at zero temperature as a function of baryon chemical potential, or an imaginary-chemical-potential lattice measurement, could test whether a sharp jump from zero to nonzero actually occurs at the model's quarkyonic-to-deconfined transition; if no such jump exists, the transition is a regulator artifact. A purely empirical cross-check is that if precision mass-radius and gravitational-wave data exclude the stiff low-density equations of state used here, the predicted above-two-solar-mass quark cores cannot exist.","supporting_citations":[{"cited_title":"Dexheimer, R","cited_arxiv_id":null,"evidence_quote":"Provides the analytical solution for the Polyakov loop branches used to identify confined and deconfined phases."},{"cited_title":"R¨ oßner, T","cited_arxiv_id":null,"evidence_quote":"Motivates replacing $T^4$ with the Stefan-Boltzmann pressure through a nonperturbative analysis of the Polyakov potential."},{"cited_title":"Dexheimer and S","cited_arxiv_id":null,"evidence_quote":"Provides the perturbative $T_0(\\mu_f)/T$ substitution whose low-temperature divergence the model regularizes."},{"cited_title":"Costa, M","cited_arxiv_id":null,"evidence_quote":"Supplies the lattice deconfinement temperature used to fix the coefficient $\\eta_2$."},{"cited_title":"Ratti, S","cited_arxiv_id":null,"evidence_quote":"Introduces the eight-quark vector interaction whose cubic density dependence stiffens the quark equation of state."},{"cited_title":"Xin, S.-x","cited_arxiv_id":null,"evidence_quote":"Provides the SFHo hadronic equation of state used as the soft low-density baseline."},{"cited_title":"Pereira, Chiral Transition and Deconfinement in Hybrid Stars, Master’s thesis, University of Coimbra (2016)","cited_arxiv_id":null,"evidence_quote":"Provides the DD2 hadronic equation of state whose stiffness allows hybrid stars to exceed two solar masses."},{"cited_title":"Olive, Chinese Physics C38, 090001 (2014)","cited_arxiv_id":null,"evidence_quote":"Supplies the Maxwell construction procedure used to join the hadronic and quark phases."}],"review_version":1}