{"id":"85d31a6f-7972-4b4d-b4e0-4d76f082dbef","arxiv_id":"2505.03142","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Hybrid neutron stars with axion-influenced quark matter are predicted to have enhanced quadrupolar f-mode oscillation frequencies.","lead":"This paper calculates how axions, hypothetical particles, could change the interior composition and oscillation frequencies of neutron stars with quark cores. It finds that including axionic effects in a quark matter model can increase the predicted f-mode frequencies compared to ordinary neutron stars.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The axion effect is modeled as an arbitrary vacuum angle θ up to π in Eq. (4), with no axion potential or relaxation dynamics; the claimed f-mode enhancement therefore rests on an undefended assumption that large θ can exist inside a neutron star.","rationale":"I read the paper in good faith. The authors set out to show that axions modify the quark matter EOS and thereby enhance the f-mode frequencies of hybrid neutron stars. For that claim to hold, the axion must be able to produce a large effective θ in the dense stellar medium. The paper's calculation, however, only scans a constant θ up to π in the NJL determinant (Eq. 4) and never includes the axion field dynamics or the axion potential. This is precisely the weakest point the reader identified. I agree with that diagnosis. The concern is load-bearing because every quantitative output—the phase transition density, the mixed-phase boundaries, the mass-radius curves, and the f-mode frequencies—depends on the prescribed θ. If θ cannot be systematically large inside neutron stars, the central claim about axion-induced enhancement does not follow. The concrete test I propose is to add the axion potential and determine θ_eq as a function of density. Such a computation would settle whether the θ = π regime is physically accessible or merely a parameter-space choice. The paper has independent support only in the sense that it uses established NJL and RMF frameworks, but no numerical data or code is provided, and no machine-checked verification exists. Given that the reader already conditioned acceptance on justifying the θ range, my stress-test supports that conditional verdict rather than moving to rejection: a relatively short, well-defined calculation could resolve the issue.","tokens_in":3975,"tokens_out":5723,"duration_ms":64729,"concrete_test":"Extend the calculation by adding the axion potential V(θ) = m_a^2 f_a^2 (1 − cos θ) to the NJL grand potential, or a more complete medium-dependent axion effective potential, and solve ∂Ω/∂θ = 0 together with beta equilibrium and charge neutrality at each density. Recompute the EOS, the Gibbs construction, and the f-mode frequencies using the resulting θ_eq(n_B). If θ_eq remains far below O(1) throughout the star, the claimed f-mode enhancement is an artifact of the arbitrary θ scan; if θ_eq reaches O(1), report the required axion parameters and verify that they are compatible with astrophysical and laboratory axion bounds.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The most load-bearing assumption is the treatment of axions as a fixed θ parameter. Equation (4) inserts θ only in the 't Hooft determinant, and the paper scans θ up to π, described in Section 2 as the 'maximum CP violation effect at θ = π'. For a QCD axion, θ = a/f_a is a dynamical field, not a free constant: its expectation value is set by minimizing the total thermodynamic potential including the axion potential, and in vacuum the strong CP bound forces |θ| ≲ 10^{-10}. The dense-matter axion potential can in principle have a non-trivial minimum, but the paper never computes it; the EOS, the mixed-phase boundaries, and the f-mode shifts in Figs. 1–2 are all computed as functions of an externally prescribed θ. Hence the headline statement that the enhancement is 'particularly large in the presence of axions' is not supported: the enhancement is due to an assumed large θ, not to a demonstrated axion mechanism. Reference [1] may contain a full derivation, but it is not provided here, and the present paper's self-contained argument relies entirely on the scan in θ.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the effect of axions on the equation of state (EOS) of quark matter in a three-flavor Nambu–Jona-Lasinio (NJL) model with a theta-term in the 't Hooft determinant, matched to a relativistic mean field (RMF) hadronic EOS via the Gibbs construction. The resulting hybrid EOS is used to solve the TOV equations and the linearized non-radial oscillation equations to obtain neutron star mass-radius relations and quadrupolar f-mode frequencies. The central claim is that a non-zero theta parameter (representing axions) lowers the hadron-quark phase transition density, stabilizes hybrid neutron stars, and substantially enhances f-mode frequencies relative to canonical nucleonic stars, especially when combined with hyperons.","tokens_in":4295,"tokens_out":4354,"duration_ms":44239,"significance":"If the claimed enhancement is physically realized, the result would offer an observable gravitational-wave signature of non-nucleonic degrees of freedom (hyperons and quark matter) and of a CP-violating axion background inside neutron star cores. The paper builds on established machinery (NJL and RMF EOS, Gibbs construction, TOV equations, fluid perturbation theory) and presents parameter-sensitivity plots rather than a first-principles derivation. The significance is strongly tempered by the fact that the axion is modeled as a fixed, externally prescribed vacuum angle without computing the axion field dynamics; the headline prediction is therefore a parameter-sensitivity statement rather than a demonstrated consequence of axions. The manuscript also provides no numerical tables and its central oscillation equations are garbled as typeset, hindering reproducibility and verification.","major_comments":[{"comment":"The axion effect is implemented solely as a static angle theta inserted in the 't Hooft determinant, scanned up to theta = pi ('maximum CP violation effect at theta = pi'). For a QCD axion, theta = a/f_a is a dynamical field whose expectation value minimizes the total thermodynamic potential including the axion potential; in vacuum the strong CP bound requires |theta| less than about 10^-10, and at neutron star core densities the paper does not compute the effective theta. Consequently, the statement that the f-mode enhancement is 'particularly large in the presence of axions' is not established: the enhancement is produced by an externally imposed large theta, not by a demonstrated axion mechanism. The authors should either compute the density-dependent effective theta from an axion potential or explicitly reframe the paper as a study of an arbitrary CP-violating angle and discuss the conditions under which axions could realize such a value.","section":"Section 2, Eq. (4)"},{"comment":"The linearized pulsation equations are garbled as typeset, with misplaced superscripts, exponents, and fractions (e.g., 'omega^2 r^2 e^{Lambda - 2 Phi} Z + Phi' Q' and the expression containing 'Phi' e^{-Lambda + 2 Phi} Q / (omega^2 r^2)'). As printed they are neither readable nor verifiable, and these equations are central to the computed f-mode frequencies. The paper should display the oscillation equations in a correct, standard form and define all symbols (Phi, Lambda, Q, Z, omega_BV, c_s^2) with consistent notation.","section":"Eqs. (2)-(3)"},{"comment":"The claimed f-mode enhancement is presented as a single family of curves for one choice of hadronic and quark parameters. The paper does not demonstrate robustness of the enhancement to variations of other NJL couplings (G_s, K), the RMF coupling set, or the matching procedure. Given that the EOS and hence the f-mode frequency depend on several parameters that are only partially varied (theta and G_v), the paper should include at least a sensitivity check with respect to another parameter (e.g., G_s or K) or provide a quantitative discussion of the parameter dependence, otherwise the enhancement may be a fine-tuned result rather than a robust prediction.","section":"Figure 2 and Section 2"}],"minor_comments":[{"comment":"There are numerous typographical errors, including 'on on' in the abstract, 'frquencies', 'e ffects', 'V olkoff', 'Nambu--Jona-Lasino' (inconsistent spelling), 'charge neural and beta equlibriataed', and 'frequncies' in the conclusions; the manuscript should be carefully proofread.","section":"Abstract and throughout"},{"comment":"The TOV equations are written in a single line with 'dp/dr = ... , dm/dr = ...'; these should be displayed as two separate equations for clarity.","section":"Eq. (1)"},{"comment":"The symbols Phi, Lambda, and c_s^2 that appear in the oscillation equations are not defined in the text; their definitions (relating to the equilibrium metric and the adiabatic sound speed) should be provided.","section":"Section 1, after Eq. (3)"},{"comment":"The text says 'when theta = pi and G_v = 0.1 gives larger enhancement' but the sentence is grammatically incomplete; it should specify that this case gives the largest enhancement among the models that satisfy the astrophysical constraints.","section":"Section 2, Fig. 2"},{"comment":"The figure caption refers to 'axion parameter theta' but the figure legend and text do not state the values of theta and G_v clearly for each curve; a table or explicit legend would improve readability.","section":"Section 2, Fig. 1"},{"comment":"The phrase 'second densest object in the universe after a black hole' is imprecise, as black holes are not characterized by a density in the same way as neutron stars; the sentence should be rephrased.","section":"Introduction"}],"recommendation":"major_revision","confidential_remarks":"The core computation follows established methods and the paper may contain an interesting exploration of the role of a CP-violating angle in hybrid neutron star oscillations. However, the identification of the theta-scan with 'axions' is the central interpretative step and is not justified by any dynamical calculation; this is a load-bearing issue that must be addressed. The garbled typesetting of Eqs. (2)-(3) also strongly impedes verification. The paper is likely to be suitable after a major revision that either adds a proper treatment of the axion potential or explicitly restricts the claims to a phenomenological theta parameter."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: the central qualitative claim—axions enhance f-modes in hybrid stars—was already in the authors' own earlier preprint (arXiv:2411.17828). What's actually new here is the hyperonic hadronic sector and a more complete hybrid star description with a scan over the vector coupling Gv. That's a legitimate but incremental step, not a new phenomenon.\n\nThe paper does some things well. The numerical chain is standard: RMF for hadrons, NJL for quark matter, Gibbs construction for the mixed phase, TOV for structure, and the linearized oscillation equations for f-modes. The mass-radius curves and f-mode frequencies behave as expected—hyperons soften, Gv stiffens, and the θ=π, Gv=0.1Gs case satisfies the maximum-mass constraint. The problem is the presentation: Eqs. (2)–(3) are garbled, there are typos throughout, and there are no numerical tables, derivation details, or code, so the numbers can't be checked directly.\n\nThe soft spot is the axion sector, and it is load-bearing. The paper treats θ as a free parameter and scans up to π, calling π the 'maximum CP violation effect.' But for a QCD axion, θ is a dynamical field, not a constant. Its expectation value minimizes the full thermodynamic potential, and in vacuum the strong CP bound forces |θ|≲10^-10. Dense matter can in principle shift that minimum, but the paper never computes the axion potential—it just varies θ by hand. So the f-mode enhancement is a sensitivity to an assumed large CP-violating angle, not a demonstrated consequence of axions. The stress-test note is right about this. The authors cite their previous paper for the full derivation, but this manuscript is not self-contained on the point that matters most. If the dense-matter axion potential really does relax to near θ=π, the claim would be interesting; without computing it, the prediction is conditional on an unexamined premise.\n\nThe novelty question also needs addressing. The title is nearly identical to the earlier preprint, and the abstract doesn't say what is new beyond that work. A referee should ask for a clear statement of the new contribution, not just the hyperon extension.\n\nWho is this for: people working on hybrid star phenomenology who want a parameter study of how f-modes shift with θ, hyperons, and vector coupling. It is a useful data point, but not a first-principles axion prediction.\n\nRecommendation: send it to peer review, but conditional on revisions. Require a justification of the θ range—ideally a computation of the axion effective potential in dense matter—plus a statement of novelty and either clean equations or a public dataset. I wouldn't desk-reject, but I also wouldn't accept as-is.","headline":"The hyperonic extension is a real but incremental step; the axion effect is an assumed θ up to π, not a computed axion potential, so the central prediction is conditional.","tokens_in":4748,"tokens_out":3745,"would_cite":false,"duration_ms":33688,"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":"The paper claims that axionic quark matter enhances the f-mode frequencies of hybrid neutron stars beyond canonical nucleonic predictions, making the mode a probe of non-nucleonic degrees of freedom.","keywords":["axion","neutron star oscillations","f-mode","hybrid star","Nambu–Jona-Lasinio model","hadron-quark phase transition","equation of state","CP violation"],"falsifier":"A single neutron star with a well-measured mass, radius, and f-mode frequency whose measured frequency falls on the canonical nucleonic curve rather than the $\\theta>0$ hybrid curves would falsify the claimed enhancement. Equivalently, a tight laboratory or astrophysical bound showing that the effective $\\theta$ inside dense matter must be far below $\\pi$ would remove the mechanism's quantitative basis.","tokens_in":3812,"feed_emoji":"🌊","tokens_out":4018,"duration_ms":42066,"temperature":0.7,"pith_summary":"This paper argues that axions can leave a measurable imprint in neutron star oscillations. Combining a hadronic equation of state with quark matter described by a three-flavor NJL model that includes a CP-violating axion angle $\\theta$, the authors construct hybrid stars with a mixed hadron-quark phase. They claim that the axion lowers the density at which quark matter appears, stabilizes hybrid stars against astrophysical mass constraints, and enhances the quadrupolar f-mode frequencies compared with canonical nucleonic stars. If correct, a detected f-mode frequency noticeably above nucleonic predictions would be evidence for non-nucleonic degrees of freedom in the star's core.","feed_headline":"Axions could raise neutron-star f-mode frequencies","feed_subtitle":"A quark-matter core with axion effects would make hybrid stars oscillate faster than nucleonic predictions.","key_machinery":"The load-bearing object is the three-flavor NJL Lagrangian with the Kobayashi-Maskawa-'t Hooft determinant term multiplied by $e^{\\pm i\\theta}$, which lets an axion-induced CP-violating angle shift the quark-matter equation of state. Around it sits a two-layer construction: the relativistic mean-field model supplies the hadronic, hyperon-rich equation of state, and the Gibbs condition joins the hadronic and quark phases into a hybrid-star equation of state. The non-radial oscillation frequencies follow from solving the perturbed Einstein equations for quadrupolar modes, with the Brunt-Väisälä frequency entering through the sound-speed profile. The mechanism that produces the claimed signature is that increasing $\\theta$ lowers the density threshold for quark-matter appearance, enlarging the quark core and changing the f-mode frequency.","core_discovery":"The central claim is that in hybrid neutron stars, the fundamental quadrupolar oscillation (f-mode) frequency is enhanced both by hyperons and by quark matter, and the enhancement is largest when axionic effects are present. The paper models axion effects as a vacuum angle $\\theta$ in the flavor-determinant interaction of the three-flavor NJL Lagrangian; at maximum CP violation, $\\theta=\\pi$, the mixed phase begins at lower density (about $1.9n_0$) than without axions, producing a larger quark core and, once vector repulsion is included, a stiffer overall equation of state. The mass-radius curves satisfy modern astrophysical constraints except for the $\\theta=\\pi$, zero-vector-coupling case, which fails the maximum-mass bound. The authors conclude that a measured f-mode frequency noticeably above nucleonic-only predictions would indicate quark or hyperonic content in neutron star matter.","pith_inferences":["A dynamical-axion treatment, where $\\theta$ evolves with density rather than staying fixed, would be the natural next step; the present constant-$\\theta$ scan likely overstates the effect if the axion relaxes toward zero inside the star.","The same $\\theta$-enhanced quark matter should also alter tidal deformability and cooling; those independent observables could be cross-checked against the f-mode claim.","The abrupt sound-speed discontinuities at the mixed-phase boundaries would imprint on the gravitational-wave ringdown spectrum, potentially letting detectors hear the phase-transition structure alongside the frequency shift."],"forward_implications":["Measured f-mode frequencies higher than nucleonic-only predictions become a direct diagnostic for non-nucleonic degrees of freedom rather than just an equation-of-state curiosity.","Gravitational-wave observations of f-modes, for example from post-merger remnants, could place constraints on the axion parameter $\\theta$ if the equation of state is otherwise known.","Vector repulsion in quark matter is necessary for the axion-enhanced hybrid stars to satisfy maximum-mass constraints; without it, the $\\theta=\\pi$ case is ruled out by observations.","Because axions lower the onset density of the hadron-quark mixed phase, precision oscillation measurements could map where the phase transition happens inside the star."],"supporting_citations":[{"why":"Predecessor work that derives the CP-violating quark matter equation of state and motivates the theta-dependent NJL treatment.","marker":"[1]"},{"why":"Supplies the non-radial oscillation equations and the hybrid-star mixed-phase formalism used to compute f-modes.","marker":"[2]"},{"why":"Source of the theta-vacuum NJL Lagrangian with the 't Hooft determinant interaction.","marker":"[3]"},{"why":"Basis for CP violation and chiral symmetry breaking in hot and dense quark matter within the NJL model.","marker":"[4]"},{"why":"Provides the relativistic mean-field hadronic Lagrangian with hyperons used for the outer layer of the hybrid star.","marker":"[5]"}],"fun_headline_variants":["Axion-driven quark cores boost f-mode frequencies","Axions amplify f-modes in hybrid neutron stars","Axion-rich cores raise neutron star oscillation frequencies","Axions in quark matter enhance neutron star f-modes","Axions make hybrid stars ring with higher f-modes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole enhancement rests on treating the axion effect as a fixed vacuum angle that can be as large as $\\theta=\\pi$ inside the star; if the strong CP bound or the dynamics of the axion field keep the effective angle minuscule at stellar densities, the predicted f-mode boost disappears.","fun_headline_variants_meta":{"raw":{"variants":["Axion-driven quark cores boost f-mode frequencies","Axions amplify f-modes in hybrid neutron stars","Axion-rich cores raise neutron star oscillation frequencies","Axions in quark matter enhance neutron star f-modes","Axions make hybrid stars ring with higher f-modes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.0009,"raw_usage":{"total_tokens":3846,"prompt_tokens":886,"completion_tokens":2960,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":502,"completion_tokens_details":{"reasoning_tokens":2884}},"tokens_in":502,"tokens_out":2960,"duration_ms":22983,"temperature":1.0,"reasoning_tokens":2884,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:57:59.969807+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A single neutron star with a well-measured mass, radius, and f-mode frequency whose measured frequency falls on the canonical nucleonic curve rather than the $\\theta>0$ hybrid curves would falsify the claimed enhancement. Equivalently, a tight laboratory or astrophysical bound showing that the effective $\\theta$ inside dense matter must be far below $\\pi$ would remove the mechanism's quantitative basis.","supporting_citations":[],"review_version":1}