{"id":"bd4d6418-3282-4ffc-bb6c-787ee1990406","arxiv_id":"2412.16636","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Elastic quark cores in hybrid stars increase maximum mass and compactness, and a new fitting model reproduces the resulting pressure anisotropy.","lead":"This paper models hybrid neutron stars with a solid, elastic quark core and finds that the core's shear stiffness increases the maximum mass by a few percent. The result matters because it changes which equations of state can satisfy current pulsar mass measurements and because ultra-compact elastic stars could mimic black holes.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claim that stable elastic hybrid stars can exceed compactness 1/3 is not established, because the turning-point stability criterion does not apply to stars with a sharp phase transition; the paper's own footnote 10 concedes this, and a radial perturbation analysis is required.","rationale":"The reader's weakest assumption identifies exactly the same load-bearing concern, so I agree. The paper's abstract and Section III B describe the ultra-compact models as stable based only on their position before the mass-radius turning point, but the turning-point theorem is not valid for stars with a sharp phase-transition discontinuity; footnote 10 and the conclusion explicitly acknowledge this. A full radial perturbation analysis would settle whether any radially stable configuration actually reaches compactness greater than 1/3. The rest of the paper, including the maximum-mass increase of several percent and the new anisotropy fit, appears internally consistent and is supported by the cross-check between the structural equations and the anisotropic TOV equations in Appendix D, so the overall verdict remains conditional rather than reject. The conclusion should be reflected in the abstract by either performing the stability computation or qualifying the 'stable stars' phrasing.","tokens_in":24713,"tokens_out":7446,"duration_ms":62162,"concrete_test":"For the brown HS model in Fig. 7 and HS-A as a control, solve the relativistic radial pulsation equations for the quasi-Hookean elastic matter derived in Karlovini, Samuelsson and Zarroug (2004), imposing junction conditions at the hadron-quark interface for a sharp first-order phase transition in both slow and rapid phase-conversion limits. Compute the squared frequency of the fundamental radial mode as a function of central pressure across the turning point. If the squared frequency is negative for any configuration with compactness greater than 1/3, the claim that stable stars can exceed that compactness is false; if it is non-negative on the entire high-compactness branch, the concern is resolved.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Section III B and Fig. 7 label all models with central pressure below the mass-radius turning point as the 'stable branch' and use this to claim that stable elastic HSs can exceed compactness 1/3. This relies on the classical turning-point theorem for perfect-fluid stars. For hybrid stars with a first-order hadron-quark transition and a density discontinuity, the theorem need not hold: the number of unstable radial modes can change without a turning point, and a turning point can occur without a change in stability, depending on the phase-conversion timescale (Karlovini, Samuelsson and Zarroug 2004; Pereira, Flores and Lugones 2018). Footnote 10 concedes exactly this, and the conclusion states that the turning-point criterion 'may not apply to a star with a density discontinuity.' Consequently, the abstract's assertion that 'the compactness of stable stars can exceed 1/3' is not established by the paper's calculations: the ultra-compact brown model in Fig. 7 may be radially unstable before the turning point, or stable after it, and no radial perturbation spectrum is computed. The black-hole-mimicker interpretation is therefore conditional on an unverified stability assumption. This is the most load-bearing concern because the novelty of the work, as presented, includes stable ultra-compact objects, not merely the few-percent maximum-mass increase, which is more robust and less controversial.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper models static, spherically symmetric hybrid stars whose quark core is treated as a nonlinearly elastic solid described by the quasi-Hookean equation of state of Karlovini and Samuelsson, with a CSS-type unsheared EOS and an APR fluid envelope. The authors find that the background shear (nonzero even in spherical symmetry) induces tangential pressure exceeding radial pressure, increasing the maximum mass by a few percent and enhancing compactness. They further show that standard phenomenological anisotropy models (Bowers-Liang and Horvat et al.) fail to reproduce the numerically obtained anisotropy profile, and they propose a new parametrized anisotropy fitting formula claimed to be accurate to about 10%. The paper also discusses causality constraints, junction conditions at the hadron-quark interface, and revisits a discrepancy with Karlovini and Samuelsson on the sign of the elastic maximum-mass effect.","tokens_in":25078,"tokens_out":4036,"duration_ms":39513,"significance":"If the central claims are established, the paper provides a useful bridge between relativistic elasticity theory and the phenomenological anisotropic-star literature, and it identifies a potentially observable effect of a crystalline quark core on the M-R relation. The numerical work has genuine strengths: the structural equations are solved with two independent codes, and the consistency between the elasticity equations and the anisotropic TOV equations is checked to about 1e-6; a Mathematica notebook and fitting code are made publicly available. However, the most striking advertised result, that stably supported elastic hybrid stars can exceed compactness 1/3 and serve as black hole mimickers, is not actually established by the calculations, because the stability boundary is identified with the turning point of the M-R relation even though the paper itself concedes that this criterion may fail for stars with a density discontinuity. The proposed anisotropy model is also fitted and evaluated on the same numerical data, so its 10% accuracy is in-sample. These issues materially weaken the headline claims, though the maximum-mass enhancement by several percent appears more robust.","major_comments":[{"comment":"The claim that \"the compactness of stable stars can exceed 1/3\" is not supported by the analysis. Figure 7 labels the portion of each M-R curve before the turning point as the \"stable branch,\" and the abstract and Section III.B use this to infer stable ultra-compact configurations. However, as Footnote 10 and the Conclusion state, the turning-point stability criterion need not apply to stars with a sharp hadron-quark discontinuity, and the number of unstable radial modes can change independently of a turning point depending on the phase-conversion timescale. Since no radial perturbation calculation is performed for these elastic hybrid stars, the compactness-exceeding-1/3 claim and the black-hole-mimicker interpretation are unverified. The paper should either compute the radial mode spectrum (e.g., using the Karlovini-Samuelsson-Zarroug formalism) or explicitly remove the word \"stable\" and the black-hole-mimicker inference from the abstract and Section III.B.","section":"§III.B, Fig. 7; abstract; conclusion"},{"comment":"The qualitative disagreement with Karlovini and Samuelsson (2003) is a load-bearing unresolved issue. The paper reports the opposite trend for the maximum mass as a function of the shear parameter k and merely states that the increasing trend is consistent with earlier studies of anisotropic stars with pt > pr. Since the entire construction is based on Equations (14)-(16) from that reference, an unexplained sign discrepancy in a benchmark case leaves open the possibility of an error in the equations, boundary conditions, or the junction condition. The 1e-6 consistency check between the structural equations and the anisotropic TOV equations is valuable but only checks internal consistency of the present implementation, not correctness against an independent source. The authors should identify the specific source of the difference, such as a sign convention for σ, a different junction condition, or a different stability criterion, or provide a more detailed comparison with the Karlovini-Samuelsson calculation.","section":"Appendix D and §III.B"},{"comment":"The claimed 10% accuracy of the new parametrized anisotropy model is an in-sample result and is therefore overstated. The fitting function in Eqs. (25)-(26) is calibrated on the same numerical solutions whose fractional residuals are plotted in Fig. 8; there is no out-of-sample test at parameter values excluded from the fit or a comparison with an independent family of elastic EOSs. With the many coefficients in Tables II and III (four polynomial degrees in both pc and κ plus inverse powers), the fit has enough flexibility that 10% agreement on training data is not strong evidence that the model \"can accurately capture physically-motivated profiles across a wide parameter space.\" A holdout test, a cross-validation statement, or a more explicit statement that the error is a training-set fitting error would be needed.","section":"§IV, Eqs. (25)-(26), Fig. 8"}],"minor_comments":[{"comment":"The conclusion states that \"the maximum mass increases by around 10%,\" but Section III.B reports maximum-mass increases of 5.242%, 3.263%, 4.648%, and 1.844% for HS-A through HS-D at the largest κ. Please reconcile these numbers; the abstract's \"several percent\" is consistent with the latter, but the conclusion is not.","section":"Section V vs. Section III.B"},{"comment":"The bottom panel of Fig. 2 is difficult to read because the fractional-difference axis uses a very wide logarithmic range and the H and BL model curves cross the numerical curve in multiple places; labeling the vertical axis more clearly or using a more conventional relative-error scale would help.","section":"Fig. 2"},{"comment":"The sentence \"The model is derived by fitting the numerical anisotropy profiles\" should be moved earlier in the section so that the reader immediately understands that Eq. (1) is a fitting formula and not a physics-based EOS; the current wording in the Introduction could be misread as implying a first-principles derivation.","section":"Section IV"},{"comment":"The junction-condition discussion is clear but would benefit from a sentence stating explicitly which of the two cases (zero surface energy vs. thin shell) is used for the main results of Section III, since only the zero-surface-energy case is actually adopted in the body of the paper.","section":"Appendix C"}],"recommendation":"major_revision","confidential_remarks":"The paper has solid numerical machinery and a useful fitting tool, but the headline 'stable compactness > 1/3' claim is not established by the submitted analysis, and the unresolved discrepancy with Karlovini and Samuelsson in Appendix D should be resolved or at least discussed much more carefully before the work is published. The in-sample nature of the anisotropy fit should also be disclosed in the abstract or introduction, otherwise readers may over-interpret the 10% accuracy statement."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing you should know: the structural modeling is solid and new, but the abstract's claim about stable compactness exceeding 1/3 is not backed by the paper's own stability analysis. The authors apply the Karlovini-Samuelsson quasi-Hookean elasticity framework, solve the full background equations with shear, and use it for hybrid stars with APR+CSS EOSs. That is genuinely new. Maximum mass increases by a few percent (1.8–5.2% in the body), which can help some soft EOSs satisfy the PSR J0740+6620 constraint, and they correctly note that current observations cannot distinguish elastic from fluid cores. Their new parametrized anisotropy model, fitting numerical profiles to ~10% across a wide parameter space, is a useful bridge between elastic EOSs and phenomenological anisotropic star models.\n\nThe internal checks are strong: two independent codes agree, and the quasi-Hookean structural equations reproduce the anisotropic TOV equations to about 1e-6. That gives me confidence the solutions are right. The wave-speed causality analysis in Appendix B is a nice addition, and the authors are candid about limitations—footnote 10 and the conclusion explicitly say the turning-point stability criterion may not hold for stars with a density discontinuity. But that candor makes the abstract's wording worse. The paper shows configurations before the M–R turning point can have compactness >1/3, not that they are stable. For hybrid stars with a sharp phase transition, the turning point does not guarantee a change in radial stability; that is known from Karlovini, Samuelsson and Zarroug 2004 and Pereira, Flores and Lugones 2018. Without a radial perturbation spectrum, the black-hole-mimicker claim is conditional. This is the main thing to fix, and it is fixable by rewriting the abstract and the relevant section.\n\nTwo minor issues. The conclusion says the maximum mass increases by 'around 10%', but the body quotes at most 5.2% for the models shown; either an error or an undocumented model. And Appendix D reports an opposite trend to Karlovini and Samuelsson without explaining why; the reproducibility check does not resolve the discrepancy. The new anisotropy fit is in-sample, so the 10% error is descriptive, not predictive—fine as an interpolation tool, but not a validation.\n\nBottom line: this deserves a serious referee and likely publication after revision. The stability framing needs correction and the max-mass inconsistency needs resolution. I would cite it for the anisotropic modeling and elastic-HS framework, and I would bring it to the reading group.","headline":"Solid new elastic-HS modeling, but the abstract's stable-compactness>1/3 claim outruns the stability analysis.","tokens_in":25562,"tokens_out":2824,"would_cite":true,"duration_ms":24974,"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 paper claims that a nonlinearly elastic quark core, modeled with the quasi-Hookean equation of state, raises the maximum mass of hybrid stars by several percent and can push turning-point-stable compactness above 1/3, making such…","keywords":["hybrid stars","elastic quark core","quasi-Hookean equation of state","pressure anisotropy","maximum mass","compactness","black hole mimickers","relativistic elasticity"],"falsifier":"Compute the fundamental radial oscillation eigenfrequency for the maximum-mass elastic hybrid star with $\\kappa = 7 \\times 10^{26} \\, \\mathrm{cm}^{1/2} \\, \\mathrm{g}^{1/2} \\, \\mathrm{s}^{-2}$ in the quasi-Hookean model; if the mode is unstable at or before the turning point, or if a nonlinear evolution drives the ultra-compact configuration to collapse or shed its photon sphere, the claim that these stars are stable black hole mimickers would be refuted.","tokens_in":1715,"feed_emoji":"⭐","tokens_out":5489,"duration_ms":82732,"temperature":0.7,"pith_summary":"This paper asks what happens to the static structure of a hybrid star when its quark core is treated as a real elastic solid, not a perfect fluid. It finds that the shear stress the core supports under its own gravity makes the pressure anisotropic, with tangential pressure exceeding radial pressure, and this extra support raises the maximum mass by several percent (1.8% to 5.2% in the four main models). Some previously disfavored soft equations of state then satisfy the observed mass of PSR J0740+6620, and the stable branch can reach compactness above 1/3, which would make these objects potential black hole mimickers. The paper also proposes a new fitted formula for the pressure-anisotropy profile that matches the elastic-core calculation to within 10%, because standard phenomenological anisotropy models fail to capture it.","feed_headline":"Solid quark cores lift hybrid stars past the 1/3 compactness limit","feed_subtitle":"Elasticity adds pressure support, raising mass limits and letting softer equations of state pass pulsar bounds.","key_machinery":"The load-bearing object is the quasi-Hookean equation of state for a relativistic solid, $\\rho = \\tilde{\\rho} + \\tilde{\\mu} S^2$, together with the structural equations obtained from the general relativistic elasticity formalism. The shear scalar $S^2 = (1/6)(z^{-1} - z)^2$ is controlled by $z = n_r/n_t$, the ratio of radial to tangential linear particle density; the shear modulus is taken as $\\tilde{\\mu} = \\kappa \\sqrt{\\tilde{\\rho}}$ with $\\kappa$ set by crystalline color-superconducting estimates. The anisotropy follows as $\\sigma = -(\\tilde{\\mu}/2)(z^{-2} - z^2)$, and the central equations integrate for the mass $m$, the unsheared pressure $\\tilde{p}$, and $z$, with a junction condition at the quark–nuclear interface. This machinery converts the microscopic rigidity of the core into macroscopic changes in mass, radius, and compactness, and it generates the numerical anisotropy profiles that the new fitting formula reproduces.","core_discovery":"The central discovery is that relaxing the usual assumption of an unsheared background changes hybrid star structure in a specific, calculable way. In the fully relativistic elasticity framework used here, the quasi-Hookean quark core develops a shear strain characterized by $z = n_r/n_t < 1$, so the tangential pressure exceeds the radial pressure throughout the core. That anisotropy supplies extra pressure support: the maximum mass rises monotonically with the shear-modulus coefficient $\\kappa$, and the compactness at the maximum-mass configuration increases by up to about 5.8%. For some equation-of-state parameters the stable, turning-point branch crosses compactness 1/3 before the maximum mass, so these elastic hybrid stars would be horizonless objects with a photon sphere, i.e., black hole mimickers. The paper explicitly flags that this last conclusion depends on the turning point signaling radial instability, which is not guaranteed for stars with a sharp phase transition.","pith_inferences":["Beyond the paper's claims: if the elastic core changes the high-mass structure this much, the tidal deformability and related universal relations of hybrid stars should shift by comparable percentages, so gravitational-wave measurements of high-mass binaries could indirectly test the solid-core hypothesis.","The black-hole-mimicker conclusion is more fragile than the mass enhancement: it inherits the turning-point stability assumption, and the paper's own light-ring discussion suggests nonlinear instabilities may destroy these ultra-compact configurations even if linear radial stability holds.","The 10% fitting formula could be tested against microscopic calculations of the crystalline color-superconducting shear modulus at lower densities, where the $\\kappa \\sqrt{\\rho}$ scaling is known to overestimate rigidity.","A targeted observation of a high-mass pulsar near 2.1 to 2.3 solar masses with a precise radius could decide between elastic and fluid cores, since the elastic models are stiffer precisely in that regime."],"forward_implications":["Some soft equations of state that fail the PSR J0740+6620 mass constraint in fluid models can pass it once elasticity is included.","The maximum mass and the compactness of hybrid stars increase monotonically with the shear-modulus coefficient, so stiffer solid cores make heavier, more compact stars.","Turning-point-stable configurations with compactness above 1/3 exist in this model, placing elastic hybrid stars among candidate black hole mimickers.","The new anisotropy fit gives a 10%-accurate phenomenological bridge between microphysical elastic equations of state and parametrized anisotropic-star models, usable for future tidal and oscillation studies.","Current mass–radius observations cannot distinguish an elastic quark core from a fluid one, because the elastic effects concentrate near the maximum mass."],"supporting_citations":[{"why":"Supplies the fully relativistic elasticity theory in which the background shear of the star is computed.","marker":"[23]"},{"why":"Provides the quasi-Hookean equation of state and the structural equations for m, p̃, and z that the paper integrates for elastic hybrid stars.","marker":"[24]"},{"why":"Showed that elastic stars can be radially stable with radii inside the unstable light ring, the basis for the black-hole-mimicker interpretation.","marker":"[28]"},{"why":"Gives the shear modulus of crystalline color superconducting quark matter and the scaling used to set the range of κ.","marker":"[37]"},{"why":"Defines one of the standard phenomenological anisotropy models that the paper shows fails to fit the elastic-core profile.","marker":"[49]"},{"why":"Defines the other standard anisotropy model that fails against the numerical profile.","marker":"[54]"},{"why":"Provides the high-precision mass measurement of PSR J0740+6620 that soft equations of state must satisfy once elasticity is included.","marker":"[61]"},{"why":"Supplies the constant-speed-of-sound template used for the unsheared quark matter equation of state.","marker":"[64]"},{"why":"Provides the nuclear matter equation of state used for the hadronic envelope.","marker":"[63]"}],"fun_headline_variants":["Elastic quark cores boost hybrid star mass and compactness","Shear pressure lifts hybrid stars toward black hole mimickers","Quark core elasticity lets softer stars pass pulsar constraints","Elastic hybrid stars can exceed the 1/3 compactness barrier","Modeling elastic quark cores: higher mass, higher compactness"],"cache_read_input_tokens":27648,"weakest_assumption_plain":"The black-hole-mimicker claim rests on assuming that the turning point of the mass-radius curve marks the onset of radial instability for hybrid stars with a sharp density discontinuity, an assumption the paper itself notes may not hold.","fun_headline_variants_meta":{"raw":{"variants":["Elastic quark cores boost hybrid star mass and compactness","Shear pressure lifts hybrid stars toward black hole mimickers","Quark core elasticity lets softer stars pass pulsar constraints","Elastic hybrid stars can exceed the 1/3 compactness barrier","Modeling elastic quark cores: higher mass, higher compactness"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000316,"raw_usage":{"total_tokens":1840,"prompt_tokens":1044,"completion_tokens":796,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":660,"completion_tokens_details":{"reasoning_tokens":711}},"tokens_in":660,"tokens_out":796,"duration_ms":7165,"temperature":1.0,"reasoning_tokens":711,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T10:23:09.765734+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the fundamental radial oscillation eigenfrequency for the maximum-mass elastic hybrid star with $\\kappa = 7 \\times 10^{26} \\, \\mathrm{cm}^{1/2} \\, \\mathrm{g}^{1/2} \\, \\mathrm{s}^{-2}$ in the quasi-Hookean model; if the mode is unstable at or before the turning point, or if a nonlinear evolution drives the ultra-compact configuration to collapse or shed its photon sphere, the claim that these stars are stable black hole mimickers would be refuted.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the fully relativistic elasticity theory in which the background shear of the star is computed."},{"cited_title":"Haskell, N","cited_arxiv_id":null,"evidence_quote":"Provides the quasi-Hookean equation of state and the structural equations for m, p̃, and z that the paper integrates for elastic hybrid stars."},{"cited_title":"Anglani, R","cited_arxiv_id":null,"evidence_quote":"Showed that elastic stars can be radially stable with radii inside the unstable light ring, the basis for the black-hole-mimicker interpretation."},{"cited_title":"Cardoso, L","cited_arxiv_id":null,"evidence_quote":"Gives the shear modulus of crystalline color superconducting quark matter and the scaling used to set the range of κ."},{"cited_title":"Choudhury, T","cited_arxiv_id":null,"evidence_quote":"Defines the other standard anisotropy model that fails against the numerical profile."},{"cited_title":"Dev and M","cited_arxiv_id":null,"evidence_quote":"Provides the high-precision mass measurement of PSR J0740+6620 that soft equations of state must satisfy once elasticity is included."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the constant-speed-of-sound template used for the unsheared quark matter equation of state."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the nuclear matter equation of state used for the hadronic envelope."}],"review_version":1}