{"id":"71a9f828-fc26-4afd-baf4-12f7d85229ba","arxiv_id":"2411.17828","paper_version":3,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":8,"one_line_summary":"Using a 3-flavor NJL model with axions, the authors find that a CP-violating angle theta=pi reduces the chiral transition density and, with a vector repulsion, yields stable hybrid neutron stars with enhanced f-mode oscillation frequencies.","lead":"This paper calculates how a CP-violating angle in the QCD axion sector changes the quark matter equation of state and then predicts changes in neutron star mass-radius curves and oscillation frequencies. The key claimed result is that axions can stabilize hybrid neutron stars with quark matter cores and boost the fundamental oscillation mode frequency.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Axion results use θ=π, which the paper's own effective potential identifies as a maximum; a dynamical QCD axion would relax to θ=0 and the reported effects would disappear.","rationale":"The most load-bearing assumption is the physical relevance of θ=π, exactly as the reader's weakest_assumption states. The paper's own Figure 7 and Section VI A show that the axion effective potential has a maximum at θ=π and minima at θ=2nπ; therefore, a dynamical axion treated as a field (as implied by θ=⟨a⟩/fa in Section II) would not reside at θ=π without an external pinning mechanism. The paper provides none. This is not merely a disagreement with 'outside current consensus'—it is an internal inconsistency between the model's effective potential and the parameter value used for all headline results. The authors also acknowledge in the introduction that experimental limits constrain θ to be extremely small, further underscoring that θ=π is not a justified representative case. The concrete test of minimizing the effective potential self-consistently would decisively settle the concern: if the minimum is always θ=0, the reported axion-driven stabilization and f-mode enhancement disappear, confirming the reader's REJECT verdict. I see no independent support (e.g., code release, parameter-free derivation, or falsifiable prediction) that rescues the central claim as stated. The paper may be useful as a study of θ-vacuum effects in hybrid stars, but not as a physical claim about QCD axions in neutron stars. Therefore the reader's verdict of REJECT is unchanged by this stress-test pass.","tokens_in":31239,"tokens_out":3256,"duration_ms":34357,"concrete_test":"Recompute the hybrid-star EOS, M-R relations, and f-mode frequencies with θ determined dynamically at each baryon chemical potential by imposing ∂Ω/∂θ=0 and ∂²Ω/∂θ²>0 on the zero-temperature thermodynamic potential of Eq. (3), solved simultaneously with the gap equations (5)–(10) and the beta-equilibrium/charge-neutrality conditions. If, as Fig. 7 indicates, the global minimum is always θ=0 mod 2π (or if no stationary point exists near θ=π), the θ=π-based results in Figs. 9–12 collapse to the θ=0 curves, confirming that the claimed axion stabilization and 110 Hz f-mode enhancement are artifacts of fixing θ at a maximum rather than minimizing the axion potential.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim (abstract; Section VI C, Figs. 9 and 12) is that axions stabilize hybrid neutron stars and enhance f-mode frequencies, with the representative case θ=π, Gv=0.1Gs yielding Mmax=2.05 M⊙ and a 110 Hz enhancement. However, the paper's own axion effective potential, shown in Fig. 7 and discussed in Section VI A, has degenerate minima at θ=2nπ and maxima at θ=(2n+1)π for all chemical potentials studied, including the zero-temperature finite-density cases relevant for neutron stars. The text explicitly states that at θ=(2n+1)π 'the thermodynamic potential is still a maximum.' For a dynamical QCD axion, which the paper invokes through θ=⟨a⟩/fa and which the PQ mechanism is designed to relax to CP-conserving values, the field should settle at a minimum, θ=0 mod 2π, not at the maximum θ=π. No mechanism—pinning, domain walls, boundary conditions, or a density-induced shift of the minimum—is provided to justify θ=π. Moreover, Section I cites the experimental bound θ<0.7×10−11 for the CP-violating parameter, making θ=π even more remote as a physical input. Since all axion-dependent equation-of-state, mass-radius, and f-mode results are computed at θ=π, the central claim is conditional on an unphysical and internally inconsistent choice of the axion field value. The axion-specific effects vanish if the field is instead placed at the global minimum, reducing the results to the θ=0 hybrid-star case, which has no pure quark core for Gv=0.1Gs and a smaller f-mode enhancement.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies the effects of a CP-violating axion angle theta on cold dense quark matter within a three-flavor Nambu-Jona-Lasinio model, and uses the resulting equation of state, combined with a relativistic mean-field hadronic EOS via a Gibbs construction, to compute hybrid neutron star mass-radius relations and quadrupolar f-mode frequencies in the Cowling approximation. The main claims are that theta=pi lowers the chiral transition density, allows a pure quark matter core in stable hybrid stars when the vector coupling Gv=0.1 Gs, increases the maximum mass to 2.05 solar masses, and enhances f-mode frequencies by about 110 Hz relative to nucleonic stars.","tokens_in":31625,"tokens_out":5115,"duration_ms":45567,"significance":"If the theta=pi input were physically justified, the paper would provide a plausible model exploration of how axion-induced CP violation could alter hybrid star structure and oscillation signatures. The gap equations and thermodynamic potential treatment are standard and carefully reduced to zero temperature, and the Gibbs construction and Cowling f-mode formalism are applied consistently. The paper also makes a clear comparison with NICER and GW170817 mass-radius constraints. However, the central physical claim depends entirely on the axion field sitting at a maximum of its own effective potential, which no mechanism in the paper stabilizes. Since a dynamical QCD axion would relax to theta=0, where the reported effects largely disappear, the significance for axion physics is currently not established.","major_comments":[{"comment":"The paper's own effective potential, shown in Fig. 7 and discussed in Sec. VI A, has degenerate minima at theta=2n pi and maxima at theta=(2n+1) pi; the text explicitly states that at theta=(2n+1) pi 'the thermodynamic potential is still a maximum.' All axion-dependent results—the EOS in Fig. 8, the mass-radius curves and core profiles in Figs. 9 and 10, and the f-mode frequencies in Fig. 12—are computed at theta=pi. For the QCD axion, which the paper identifies with theta=<a>/f_a and which is designed to relax to CP-conserving values, the field would settle at a minimum, theta=0 mod 2 pi, not at the maximum theta=pi. No pinning mechanism, domain-wall stabilization, or density-induced shift of the minimum is provided. Given that the introduction itself cites the experimental bound theta<0.7 x 10^-11, the representative choice theta=pi is both dynamically and phenomenologically unjustified, and the central claim of axion-induced stabilization of hybrid stars with pure quark cores collapses unless a mechanism to maintain theta=pi is supplied.","section":"Sec. VI A, Fig. 7; Sec. VI C"},{"comment":"The claimed axion-specific enhancements are not robust under a dynamically selected theta. At theta=0 with Gv=0.1 Gs, the model does not produce a pure quark matter core (Fig. 3(b)), and the f-mode enhancement relative to a nucleonic star is 90 Hz rather than 110 Hz (Fig. 12). The difference between theta=0 and theta=pi is therefore not only quantitatively modest but also rests on the unphysical theta=pi choice. The abstract's statement that 'with the presence of axions, it is possible to have stable hybrid neutron stars having an inner core of quark matter both in pure quark matter phase as well as in a mixed phase with hyperonic matter' is thus conditional on an input that the paper's own dynamics rule out; as it stands, it overstates what the calculation establishes.","section":"Abstract; Sec. VI C, Figs. 3(b), 12"}],"minor_comments":[{"comment":"The text contains 'we we have introduced'; please correct to 'we have introduced'.","section":"Sec. II, Eq. (14)"},{"comment":"The in-text references to 'Fig 5 b' and 'Fig 5 c' in the discussion of the strange quark condensates are mislabeled; they should refer to Figs. 6(c) and 6(d).","section":"Sec. VI A, figures captions"},{"comment":"The sentence 'the maximum mass becomes 2.05 M_sun ans satisfies the maximum mass constraint' contains a typo: 'ans' should be 'and'.","section":"Sec. VI C"},{"comment":"The caption describes 'The dark gray and light gray region here correspond to 50% and 90% confidence interval', while the figure legend states '90% (solid) and 50% (dashed)'; please align the caption with the legend.","section":"Fig. 9 caption"},{"comment":"The claim that the inner-crust polytrope treatment yields radius differences of less than 0.5 km compared to a unified description is stated without a quantitative reference; please cite the specific figure or table from Ref. [69].","section":"Sec. VI C, inner crust"}],"recommendation":"reject","confidential_remarks":"The central flaw is the arbitrary placement of the axion at theta=pi despite the paper's own effective potential showing that this is a maximum, not a minimum. The dynamical axion would relax to theta=0, where the reported axion-specific effects largely disappear. I do not see a way to repair this within the current framework without introducing a new stabilization mechanism, which is beyond the scope of the manuscript. The paper's most striking conclusions are therefore conditional on an input that its own calculation rules out."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is a clean, internally consistent model calculation: three-flavor NJL with axion coupling, hyperonic RMF hadronic EOS, Gibbs mixed phase, and Cowling-approximation f-modes. The new bit is combining the axionic 3-flavor NJL with hyperonic RMF in a Gibbs construction and pushing it through to oscillation frequencies. That combination is genuinely new, the numerics look careful, and the parameters are fitted to external low-energy data rather than to the target results, so there is no circular fitting. Credit where due: the gap equations, the Gibbs treatment, and the f-mode machinery are standard but competently executed, and the paper is upfront about the Cowling approximation's ~15-20% frequency overestimate. The soft spot is load-bearing. The central physical claim—that axions stabilize hybrid stars and boost f-modes—is computed at θ=π, yet the paper's own Fig. 7 shows the thermodynamic potential has degenerate minima at θ=2nπ and maxima at θ=(2n+1)π, for every chemical potential studied. The text explicitly says the potential is still a maximum at θ=π. For a dynamical QCD axion, the field relaxes to the minimum, θ=0 mod 2π, where the reported effects vanish. The paper cites the experimental bound θ<0.7×10⁻¹¹ in the introduction, which makes θ=π an even harder sell. No mechanism—pinning, domain walls, boundary conditions, or a density-induced shift of the minimum—is proposed. So the abstract's claim about axions producing stable hybrid stars with pure quark cores is conditional on a hand-imposed, non-equilibrium angle. The calculation is fine as a parametric study of θ-vacuum effects; it is not fine as a statement about physical axions in neutron stars. Minor issues: no code or data shipped, the error assessment is a literature estimate, and the 15-20% Cowling systematic could affect the magnitude of the 110 Hz enhancement, though less so the model-to-model difference. Who gets value from this: people working on θ-vacuum effects in dense QCD and on hybrid-star asteroseismology as a model exercise. It deserves a serious referee, but the referee should demand either a dynamical story for θ=π or a reframing of the results as θ-vacuum rather than axion effects. My own verdict is skeptical: the physical claim is not supported as stated, but the model machinery is usable and the paper is honest about its limitations.","headline":"A technically sound NJL+RMF hybrid-star calculation whose headline axion effect rests entirely on θ=π, a point the paper's own effective potential identifies as a maximum.","tokens_in":713,"tokens_out":738,"would_cite":false,"duration_ms":61475,"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":"Axions could make neutron star cores ring measurably louder, this paper argues, by changing the quark matter equation of state.","keywords":["axion","strong CP violation","Nambu-Jona-Lasinio model","hybrid neutron stars","hadron-quark phase transition","Gibbs construction","f-mode oscillations","hyperonic matter"],"falsifier":"Compute the minimum of the zero-temperature finite-density axion effective potential without fixing $\\theta$; if the global minimum remains at $\\theta=2n\\pi$ for all neutron-star densities, the predicted equation-of-state softening and f-mode enhancements disappear. Observationally, a precise measurement of the f-mode frequency of a $\\sim2\\,M_\\odot$ neutron star would test the predicted $\\sim110$ Hz offset from a nucleonic star of the same mass, since the offset is the paper's concrete, mass-matched signature.","tokens_in":31009,"feed_emoji":"🌌","tokens_out":9655,"duration_ms":78050,"temperature":0.7,"pith_summary":"The paper tries to establish that axion-induced CP violation in cold dense quark matter can change the hadron-quark phase transition enough to matter for neutron star structure and gravitational-wave asteroseismology. Within a three-flavor Nambu-Jona-Lasinio (NJL) model with the axion angle $\\theta=\\langle a\\rangle/f_a$ coupled through the flavor-mixing determinant term, a nonzero $\\theta$ lowers the critical chemical potential for chiral restoration, and at $\\theta=\\pi$ it shifts the onset of the hadron-quark mixed phase to lower baryon densities. Combining this quark-matter equation of state with a hyperonic relativistic mean-field hadronic equation of state through a Gibbs construction, the paper finds stable hybrid neutron stars with a pure quark matter core and a hyperonic outer core that satisfy modern NICER and GW170817 constraints, provided quark vector repulsion is included. The paper further claims that the quadrupolar fundamental oscillation mode ($f$-mode) of such hybrid stars is enhanced by roughly 110 Hz relative to a nucleonic star of the same mass, making the composition potentially audible in gravitational waves.","feed_headline":"Axions could add 110 Hz to neutron star f-mode signals","feed_subtitle":"A 2.05-solar-mass hybrid star with a quark core stays inside NICER and GW170817 bounds—if the axion angle sits at π.","key_machinery":"The load-bearing object is the zero-temperature effective potential $\\Omega(I_s^i,I_p^i,\\theta,\\mu)$ of the three-flavor Nambu-Jona-Lasinio model, whose flavor-mixing determinant interaction carries the axion phase $e^{\\pm i\\theta}$. It generates coexisting scalar and pseudoscalar quark condensates, whose jumps mark the first-order chiral transition, and its dependence on $\\theta$ gives the axion potential: degenerate minima at $\\theta=2n\\pi$ and maxima at $\\theta=(2n+1)\\pi$. This potential supplies the quark-matter equation of state that is joined to a hyperonic relativistic mean-field hadronic equation of state by a Gibbs construction enforcing $\\beta$ equilibrium and global charge neutrality; the resulting speed-of-sound profile feeds the Cowling-approximation equations for the non-radial $f$-mode oscillations.","core_discovery":"The central claim is that at $\\theta=\\pi$, where CP violation is maximal, the quark matter equation of state in the three-flavor NJL model with vector coupling $G_v=0.1\\,G_s$ supports hybrid neutron stars with a maximum mass of $2.05\\,M_\\odot$, in agreement with the NICER measurement of PSR J0740+6620 and the GW170817 tidal constraints. The star has a pure quark matter core of radius $\\sim 1.6$ km, quark matter in a mixed phase out to $\\sim 6.5$ km, and a hyperonic outer core; without the vector term the same $\\theta=\\pi$ equation of state gives only $1.84\\,M_\\odot$ and fails the two-solar-mass constraint. The $\\theta=\\pi$ case also lowers the onset of the mixed phase from about $0.44$ fm$^{-3}$ (at $\\theta=0$) to $0.35$ fm$^{-3}$ for $G_v=0.1\\,G_s$. For the maximum-mass hybrid star, the paper computes an $f$-mode frequency near $2.1$ kHz, about $110$ Hz above a nucleonic star of the same mass, with hyperons and the enlarged quark core contributing separately to the shift.","pith_inferences":["The paper's representative case $\\theta=\\pi$ is a maximum, not a minimum, of its own axion effective potential; a dynamical axion would relax to $\\theta=2n\\pi$, where the CP-violating effects vanish, so the reported signals require an unspecified mechanism to pin the angle.","Even if the 110 Hz enhancement were observed, it would not identify axions uniquely: the paper itself shows hyperons already raise $f$-mode frequencies by about 260 Hz for a $2.35\\,M_\\odot$ star, so mass-radius-tidal joint measurements would be needed to separate axionic quark cores from merely hyperonic stars.","One testable extension is to give the axion a finite relaxation time: if the dense-matter vacuum is metastable at $\\theta=\\pi$ on neutron-star timescales, the $f$-mode shift should depend on stellar age or formation history, which could be searched for in populations of young versus old pulsars."],"forward_implications":["If $\\theta=\\pi$ with $G_v=0.1\\,G_s$, hybrid neutron stars with a pure quark matter core satisfy the $2\\,M_\\odot$ mass constraint; if $G_v=0$, they do not.","Axion presence lowers the onset density of the hadron-quark mixed phase, so quark matter can appear well below the hyperon threshold or coexist with hyperons in the same star.","The $f$-mode frequency shift (about 110 Hz for the $\\theta=\\pi$, $G_v=0.1\\,G_s$ maximum-mass star) turns the composition of the core into a potentially observable gravitational-wave signature.","A nonmonotonic speed of sound with sharp drops at phase boundaries is a generic prediction of these Gibbs-constructed hybrid equations of state, linking the model to the broad behavior inferred from neutron star observations."],"supporting_citations":[{"why":"Prior study of axion effects on hybrid star stability with a Maxwell construction; the present paper extends this to a Gibbs construction and f-modes.","marker":"[37]"},{"why":"Supplies the CP-violating NJL effective potential and the variational treatment of scalar/pseudoscalar condensates that this paper adapts to three flavors at finite density.","marker":"[29]"},{"why":"Provides the ground-state variational construct for the thermodynamic potential with CP violation that underlies Eq. (3).","marker":"[30]"},{"why":"Gives the mixed-phase Gibbs construction and the Cowling-approximation f-mode equations used here.","marker":"[41]"},{"why":"Supplies the NL3$\\omega\\rho$ parameter set for the hadronic relativistic mean-field equation of state.","marker":"[59]"},{"why":"Supplies the SU(3) NJL parameter set ($\\Lambda$, $G_s$, $K$, quark masses) used for the quark matter equation of state.","marker":"[56]"},{"why":"NICER mass measurement of PSR J0740+6620 that sets the $2\\,M_\\odot$ maximum-mass constraint the hybrid stars must satisfy.","marker":"[5]"},{"why":"GW170817 tidal-deformability constraints used to check consistency of the hybrid equation of state.","marker":"[2]"},{"why":"Establishes the two degenerate CP-conjugate vacua at $\\theta=\\pi$ that explain the pseudoscalar condensate discontinuity.","marker":"[7]"}],"fun_headline_variants":["Axions boost neutron star f-modes by 110 Hz","Quark cores with axions yield 2.05-solar-mass neutron stars","CP violation at θ=π enables quark-core hybrid stars","Maximal CP violation stabilizes quark cores in neutron stars","Hybrid stars with axion quark cores reach 2.05 solar masses"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the axion angle can be held at $\\theta=\\pi$, where the paper's own effective potential has a maximum rather than a minimum; a dynamical axion would normally settle at $\\theta=2n\\pi$ and erase the reported equation-of-state and f-mode effects.","fun_headline_variants_meta":{"raw":{"variants":["Axions boost neutron star f-modes by 110 Hz","Quark cores with axions yield 2.05-solar-mass neutron stars","CP violation at θ=π enables quark-core hybrid stars","Maximal CP violation stabilizes quark cores in neutron stars","Hybrid stars with axion quark cores reach 2.05 solar masses"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001116,"raw_usage":{"total_tokens":4730,"prompt_tokens":1114,"completion_tokens":3616,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":730,"completion_tokens_details":{"reasoning_tokens":3525}},"tokens_in":730,"tokens_out":3616,"duration_ms":25027,"temperature":1.0,"reasoning_tokens":3525,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:47:53.009218+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the minimum of the zero-temperature finite-density axion effective potential without fixing $\\theta$; if the global minimum remains at $\\theta=2n\\pi$ for all neutron-star densities, the predicted equation-of-state softening and f-mode enhancements disappear. Observationally, a precise measurement of the f-mode frequency of a $\\sim2\\,M_\\odot$ neutron star would test the predicted $\\sim110$ Hz offset from a nucleonic star of the same mass, since the offset is the paper's concrete, mass-matched signature.","supporting_citations":[],"review_version":1}