{"id":"b74b168f-9eaa-4b3e-861e-5e3bf1aa6a34","arxiv_id":"2507.18916","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Radial oscillation frequencies of neutron stars are computed in Starobinsky and Gauss-Bonnet extended gravity, revealing a dynamical exterior and a low-density plateau in the fundamental mode.","lead":"This paper models how neutron stars oscillate in Starobinsky gravity, a modified theory of gravity with extra curvature terms. It finds that strong modifications make the star's exterior react to the oscillations and can shift the density at which stars become unstable.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Stability criterion selects only exponentially decaying exterior modes; without a complex-frequency search, the sign of the fundamental ω² may not determine stability in the non-Sturm-Liouville problem.","rationale":"The reader identified the same weakest assumption: stability is judged solely by the sign of the fundamental-mode squared frequency, with overtone modes discarded because their perturbation functions oscillate rather than decay at infinity. My reading of the manuscript confirms that this is the central unprotected step. The paper explicitly contrasts the GR Sturm-Liouville case with the modified-gravity case, acknowledges that Birkhoff's theorem fails, and then imposes an exponential-decay requirement at infinity that is stronger than the stated asymptotic condition (45). Because radiative modes with 1/r falloff also satisfy the limit, the paper's exclusion of oscillatory modes needs a physical justification. Moreover, for large α the effective scalar mass is small, so several reported positive-ω² modes are plausibly above the massive-field threshold, making the boundary between normal modes and radiative modes practically important. A complex-frequency search or an energy-based stability argument would settle whether the fundamental-mode sign is sufficient. I do not see a demonstrated fatal error; the interior equations, the α→0 limit, and the Jordan/Einstein-frame consistency check are useful internal controls. Therefore the conditional verdict is appropriate, with the mode-selection/stability criterion as the condition to be resolved.","tokens_in":22952,"tokens_out":10313,"duration_ms":121349,"concrete_test":"Recompute the α=100α⋆ low-density model in Table I (ρ0=1.94×10⁻³ρ⋆) with outgoing-wave boundary conditions at infinity, δf,δφ ∼ r⁻¹ e^{+ikr}, k=√(ω²−m²_eff), and solve for complex ω by complex shooting. If an unstable mode with Im(ω)<0 appears, sign-based stability fails; if none appears and Re(ω) matches the reported fundamental frequency, the normal-mode criterion is supported. Separately, check whether the discarded oscillatory overtones satisfy Eq. (45) as written.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The load-bearing step is the mode-selection and stability criterion. After Eq. (45) and before Fig. 3, the paper restricts physical modes to those whose exterior perturbations decay exponentially at infinity and discards all overtone modes because they 'exhibit oscillatory behavior' rather than decay. This is not a harmless technicality: once Birkhoff's theorem is abandoned, the exterior is a dynamical radiative system, and the perturbation problem is no longer a Sturm-Liouville eigenvalue problem. Stability is then not guaranteed to be decided by the sign of the lowest real ω². Oscillatory solutions with 1/r falloff do satisfy the stated boundary condition (45), which only requires δf,δφ → 0 at infinity; the additional exponential-decay requirement is an extra condition not derived from the field equations. For the large-α cases emphasized in the paper, the scalar mass is small, so positive-ω² modes can lie above the massive-scalar threshold, where only radiative quasinormal modes exist. The paper neither solves for complex frequencies nor proves that no growing radiative mode exists. Without such an argument, the reported fundamental-mode sign may not settle stability, even if the interior equations and the α→0 and Jordan/Einstein-frame checks are correct.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies adiabatic radial oscillations of neutron stars in Starobinsky gravity and in a Gauss-Bonnet extension, working mainly in the Jordan-frame scalar-tensor representation with the SLy equation of state. It derives the modified TOV equations and the linear radial perturbation equations, solves them by a shooting method, and reports three main results: the exterior spacetime responds dynamically to fluid oscillations because Birkhoff's theorem fails; for large coupling α the fundamental-mode squared frequency becomes nearly independent of central density at low densities; and near the maximum-mass configuration the stability transition still approximates the GR behavior. The paper also verifies the α→0 GR limit and checks Jordan/Einstein frame consistency using a polytropic equation of state in Appendix A.","tokens_in":23164,"tokens_out":5920,"duration_ms":71554,"significance":"If the stability conclusions are correct, these are interesting results for higher-curvature gravity phenomenology: they identify a qualitative new feature (near-density-independence of the fundamental mode for large α) and a regime in which the standard dM/dρ0 stability criterion remains approximately valid. The paper is careful to state the GR limit and to provide a frame-consistency check, and the use of a realistic SLy equation of state increases the astrophysical relevance. However, the central claim depends on a mode-selection and stability criterion that is not established for the non-Sturm-Liouville problem considered here, and part of the Gauss-Bonnet analysis is deferred to unavailable supplemental material.","major_comments":[{"comment":"The explicit equations for the Gauss-Bonnet extension are not available in the manuscript. The functions Fi in Eq. (24), the coefficients Aij in Eq. (34), the vacuum coefficients Ãij in Eq. (35), the boundary coefficients δλ2 and δf2 in Eqs. (38)-(39), and the Lagrangian pressure perturbation for β≠0 are all referred to the Supplemental Material, but no such material is included with the manuscript as reviewed. Since the β≠0 results in Figs. 1-5 and Tables I-II are a stated part of the paper's central claims, these expressions must either be included in the paper or the corresponding numerical results cannot be independently verified.","section":"Sec. II.B, II.C, II.D"}],"minor_comments":[{"comment":"The caption states that entries marked with a dagger denote squared-frequency values, but the column header reads ω/(2π) [kHz]; the notation should be made uniform so that the reader can tell which entries are ω/(2π) and which are ω²/(2π)².","section":"Table III and its caption"},{"comment":"Set (5) is written as β = −10β, but Table I and the surrounding text indicate β = −10β⋆; the missing subscript should be restored.","section":"Sec. III, list of coupling sets"},{"comment":"The phrase 'the asymptotic asymptotic behavior' contains a duplicated word and should read 'the asymptotic behavior.'","section":"Appendix A, after Eq. (A21)"},{"comment":"The paper does not report numerical convergence tests or error estimates for the shooting method, although the claim that dM/dρ0 > 0 holds 'up to the second decimal place in mass' in Sec. III would benefit from such an estimate.","section":"General numerical analysis"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the paper is worth reading. It re-derives radial oscillation formalism for Starobinsky gravity in the Jordan-frame scalar-tensor representation, treating the exterior as dynamical because Birkhoff's theorem doesn't hold, and extends to a Gauss-Bonnet version. The new physical results—exterior perturbation response and a low-density plateau in ω² for large α—are clearly presented and internally consistent. The authors verify the α→0 GR limit analytically and show Jordan/Einstein frame agreement numerically for a polytropic EOS. That's solid work.\n\nWhere I get off the bus is the stability criterion. After Eq. (45) they keep only exterior solutions that decay exponentially and discard overtone modes as oscillatory. But the stated boundary condition is only that δf and δφ vanish at infinity; exponential decay is an extra restriction. Once the exterior participates, the eigenvalue problem is not Sturm-Liouville, and stability isn't decided by the sign of the lowest real ω². You need to check for complex frequencies—radiative quasinormal modes with growing amplitude. The paper doesn't. For the large-α models they emphasize, the scalar mass is small, so the fundamental mode sits close to the radiative threshold and overtones sit above it. Without a complex-frequency search, the claim that a model is stable when ω²>0 is not established.\n\nThe other soft spot is verifiability of the Gauss-Bonnet part. The coefficient functions for β≠0 are all in 'supplemental material' that isn't in the arXiv submission, and there's no code or data release, no convergence or error analysis. A referee can't check the GB numbers from the manuscript.\n\nThese are not demonstrated fatal errors. The β=0 equations are mostly in the text, and the internal checks are meaningful. But the stability conclusions are conditional on the mode-selection assumption, so I'd treat them as promising rather than settled.\n\nThis paper deserves serious refereeing. It opens a real question—whether the usual fundamental-mode criterion survives when the exterior is radiative—and the authors are honest about the non-Sturm-Liouville structure. I'd ask for the supplemental material, a complex-frequency search, and convergence tests. I'd also bring it to the reading group; the mode-selection issue is a good discussion topic.","headline":"Solid derivation and new physics in the exterior response, but the stability conclusions rest on an unjustified mode-selection rule and the Gauss-Bonnet part is unverifiable from the manuscript.","tokens_in":23691,"tokens_out":5957,"would_cite":false,"duration_ms":62250,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83D05","83C55"],"pacs":[],"model":"deepseek-v4-flash","headline":"Birkhoff's theorem fails in Starobinsky gravity: a neutron star's exterior spacetime responds to its radial oscillations, flattening the fundamental frequency against central density at low densities while preserving the maximum-mass…","keywords":["Starobinsky gravity","higher-curvature gravity","radial oscillations","neutron star stability","Birkhoff theorem","scalar-tensor gravity","Gauss-Bonnet extension","SLy equation of state"],"falsifier":"Compute the full complex-frequency spectrum of the linearized perturbation system, or evolve it numerically in the time domain, for a model the paper labels stable using the SLy equation of state and large $\\alpha$ (for example $\\rho_0 \\approx 2.5 \\times 10^{15}\\,\\mathrm{g/cm^3}$ with $\\alpha = 1000\\alpha_\\star$); if any decaying-in-time (growing) mode with nonzero imaginary part of $\\omega$ appears, the sign-of-$\\omega^2$ criterion is incomplete and the stability conclusions change. Alternatively, recompute the nearly flat $\\omega^2$ versus $\\rho_0$ curve with a stiff polytropic equation of state: if the plateau disappears, the effect is an artifact of the SLy EOS rather than a generic property of the gravity theory.","tokens_in":22745,"feed_emoji":"🌌","tokens_out":8014,"duration_ms":78432,"temperature":0.7,"pith_summary":"The paper investigates whether neutron stars remain stable against radial collapse when the gravitational action contains a curvature-squared term, as in Starobinsky gravity, or a further Gauss-Bonnet extension. It argues that the higher-derivative character of these theories breaks Birkhoff's theorem, so the spacetime outside a pulsating star is not static but responds dynamically to the fluid's oscillations. Working in the Jordan-frame scalar-tensor representation and using the SLy equation of state, the authors compute fundamental-mode spectra and find two effects: at low central densities with a large curvature-squared coupling the squared fundamental frequency becomes nearly independent of the central density, while near the maximum-mass configuration the stability-to-instability transition still occurs approximately where $dM/d\\rho_0$ changes sign, as in general relativity. If correct, the results imply that modified gravity can leave observable signatures in neutron-star pulsation frequencies and that the standard GR stability criterion is only partially preserved.","feed_headline":"Curvature-squared gravity flattens neutron star pulsation frequencies","feed_subtitle":"In Starobinsky gravity, low-density stars pulse at a nearly fixed rate, but the maximum-mass stability rule survives.","key_machinery":"The load-bearing object is the scalar-tensor representation of the theory with Lagrangian $L = R + \\Phi R - \\tfrac{1}{2}\\mu^2\\Phi^2 + U(\\Phi)L_{GB}$, in which Starobinsky gravity corresponds to $\\beta = 0$ and the Gauss-Bonnet extension to $\\beta \\ne 0$; the massive scalar $\\Phi$ carries the higher-curvature physics. The argument runs through the modified Tolman-Oppenheimer-Volkoff equilibrium equations and a four-field radial perturbation system (fluid displacement $\\xi$, metric perturbations $\\delta\\lambda$ and $\\delta f$, scalar perturbation $\\delta\\phi$), solved by shooting with boundary conditions at the center, at the surface (vanishing Lagrangian pressure perturbation), and at infinity (exponential decay of $\\delta f$ and $\\delta\\phi$). The key criterion is the sign of the lowest eigenvalue $\\omega^2$ of this system together with the requirement that the perturbation functions decay, not oscillate, at infinity.","core_discovery":"On its own terms, the paper's central claim is that in Starobinsky gravity and its Gauss-Bonnet extension, radial stability of neutron stars is governed by a dynamical exterior spacetime rather than the static vacuum of general relativity. Because of the added massive scalar degree of freedom, Birkhoff's theorem fails; the exterior metric and scalar-field perturbations obey their own differential equations, and only the fundamental mode produces perturbation functions that decay exponentially at infinity, while overtones oscillate and are discarded. For the SLy equation of state, large values of the Starobinsky coupling $\\alpha$ make the fundamental squared frequency $\\omega^2$ nearly independent of central density across low-density stellar models, and for high central densities the sign change of $\\omega^2$ still coincides approximately with the maximum-mass configuration. The authors verify that the same frequencies are obtained in the Jordan and Einstein frames of the equivalent scalar-tensor theory and extend the analysis to a Gauss-Bonnet coupling, noting that when $\\alpha$ is large the Gauss-Bonnet coupling's effect on the frequency becomes negligible.","pith_inferences":["Because the harmonic ansatz fixes real $\\omega$, the paper does not compute damping timescales; if the exterior responds dynamically, the fundamental mode likely becomes quasi-normal and radiates scalar radiation, which a time-domain evolution could reveal.","The exclusion of overtone modes on asymptotic grounds suggests the perturbation problem is not a standard Sturm-Liouville system; a full spectral analysis allowing complex frequencies could uncover radiative or growing modes that would alter the stability verdict for some configurations.","The near-flatness of the frequency-density curve at large $\\alpha$ suggests a possible observational discriminant: a population of low-mass neutron stars with nearly identical pulsation frequencies would be a signature of such curvature-squared gravity, though nuclear-physics EOS effects would have to be controlled first.","Whether the density-independence plateau and the persistence of the maximum-mass stability transition hold for other $f(R)$ or higher-derivative theories is left open by the paper; recomputing the same spectrum, for example in cubic quasi-topological gravity, would test its generality."],"forward_implications":["In these theories the exterior of a radially pulsating neutron star moves: the metric and scalar field outside the star carry time-dependent perturbations, so GR's assumption of a static Schwarzschild exterior during pulsation does not apply.","For sufficiently large $\\alpha$, the fundamental-mode squared frequency of low-density neutron stars becomes nearly constant with central density, meaning the oscillation frequency carries almost no information about the star's central density or mass in that regime.","Near the maximum-mass configuration, the criterion $dM/d\\rho_0 > 0$ still marks the transition from stability to instability to within about a percent in mass, extending the GR static-stability result to these higher-curvature models.","Only the fundamental mode is admitted by the asymptotically flat boundary conditions; overtone radial modes, which oscillate rather than decay at infinity, cannot be used to characterize stellar oscillations in this framework.","The Jordan-frame and Einstein-frame computations agree for the same star, so the frame choice does not affect the predicted radial oscillation frequencies in Starobinsky gravity."],"supporting_citations":[{"why":"Defines the curvature-squared (R^2) Starobinsky gravity that is the paper's primary subject.","marker":"[11]"},{"why":"Supplies the non-perturbative self-consistent neutron-star framework in R^2 gravity and the Jordan/Einstein frame transformation used in the consistency check.","marker":"[17]"},{"why":"Earlier radial-stability study using an effective curvature fluid that the paper redoes with the standard GR formalism.","marker":"[23]"},{"why":"Derives the equilibrium structure and boundary conditions for neutron stars in Gauss-Bonnet extended Starobinsky gravity, which the present perturbation analysis builds on.","marker":"[37]"},{"why":"Provides the Gauss-Bonnet extended Starobinsky gravity Lagrangian and the asymptotic Yukawa falloffs of the scalarized vacuum solutions.","marker":"[42]"},{"why":"Supplies the canonical general-relativity formalism for radial oscillations, Lagrangian perturbations, and the Sturm-Liouville spectral structure.","marker":"[48]"},{"why":"Provides GR radial-oscillation reference data and the polytropic equation of state used to cross-check the Jordan and Einstein frame results.","marker":"[50]"},{"why":"Supplies the unified SLy baryonic equation of state of dense matter used in the main stability analysis.","marker":"[52]"},{"why":"Provides the analytical representation of the unified SLy equation of state used in the numerical calculations.","marker":"[53]"},{"why":"States the static stability criterion $dM/d\\rho_0 > 0$ that the paper verifies approximately near the maximum-mass configuration.","marker":"[54]"}],"fun_headline_variants":["Starobinsky gravity flattens low-density neutron star frequencies","Dynamical exterior controls neutron star radial stability","Neutron star stability transition survives in Starobinsky gravity","Birkhoff's theorem fails for pulsating neutron stars in Starobinsky gravity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper assumes that the sign of the fundamental-mode squared frequency, computed with a real frequency and with overtones discarded because they oscillate at infinity, fully determines stability; it does not rule out growing complex-frequency modes of the same system.","fun_headline_variants_meta":{"raw":{"variants":["Starobinsky gravity flattens low-density neutron star frequencies","Dynamical exterior controls neutron star radial stability","Neutron star stability transition survives in Starobinsky gravity","Birkhoff's theorem fails for pulsating neutron stars in Starobinsky gravity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00019,"raw_usage":{"total_tokens":1344,"prompt_tokens":956,"completion_tokens":388,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":572,"completion_tokens_details":{"reasoning_tokens":316}},"tokens_in":572,"tokens_out":388,"duration_ms":4439,"temperature":1.0,"reasoning_tokens":316,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T18:05:41.730684+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the full complex-frequency spectrum of the linearized perturbation system, or evolve it numerically in the time domain, for a model the paper labels stable using the SLy equation of state and large $\\alpha$ (for example $\\rho_0 \\approx 2.5 \\times 10^{15}\\,\\mathrm{g/cm^3}$ with $\\alpha = 1000\\alpha_\\star$); if any decaying-in-time (growing) mode with nonzero imaginary part of $\\omega$ appears, the sign-of-$\\omega^2$ criterion is incomplete and the stability conclusions change. Alternatively, recompute the nearly flat $\\omega^2$ versus $\\rho_0$ curve with a stiff polytropic equation of state: if the plateau disappears, the effect is an artifact of the SLy EOS rather than a generic property of the gravity theory.","supporting_citations":[{"cited_title":"Strongest constraint in $f(R) = R+ \\alpha R^2$ gravity: stellar stability","cited_arxiv_id":"2008.00536","evidence_quote":"Earlier radial-stability study using an effective curvature fluid that the paper redoes with the standard GR formalism."},{"cited_title":"Neutron stars 1: Equation of state and structure,","cited_arxiv_id":null,"evidence_quote":"States the static stability criterion $dM/d\\rho_0 > 0$ that the paper verifies approximately near the maximum-mass configuration."}],"review_version":2}