{"id":"6e937036-2144-49b0-a251-8bf7eb7cc7e9","arxiv_id":"2508.21754","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A parameter study of beta-stable neutron star matter with neutron dark decay shows the symmetry energy and dark-sector interactions jointly set radii and tidal deformability, keeping the decay scenario compatible with observed masses.","lead":"This paper calculates how the nuclear symmetry energy, which sets the neutron-proton balance in dense matter, changes neutron star models that include the proposed decay of neutrons into dark matter particles. It finds that dark-matter interaction strengths and the symmetry energy together control the predicted radii and tidal deformability, and that this dark decay cannot be ruled out by current neutron star mass measurements.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The μ_n=μ_χ equilibrium condition in Eq. (17) is load-bearing but not established: if φ escapes and the reverse channel is suppressed, the χ abundance, EoS, and all M-R/Λ predictions change.","rationale":"The paper is a transparent parameter study with a clear numerical recipe, and the authors explicitly acknowledge that the environmental dependence of neutron dark decay is unresolved. Even so, the chemical-equilibrium hypothesis is the load-bearing step: every macroscopic prediction—composition fractions (Fig. 4), M-R curves (Figs. 1–3), and the concluding statement that the decay channel 'cannot be excluded'—follows from imposing μ_n = μ_χ. If the correct steady state is set by forward-decay kinetics plus Pauli blocking rather than by reversible equilibrium with an escaping φ, the resulting χ abundance can differ, and with it the EoS and all derived observables. This is not a disagreement with external consensus; it is an internal gap between the modeled equilibrium and the stated microphysics of an escaping φ. The rate-equation test proposed above would settle the issue without requiring new observational data or a change of nuclear model. I agree with the reader's identification of this as the weakest assumption, and the verdict remains conditional pending that test.","tokens_in":15807,"tokens_out":17251,"duration_ms":239987,"concrete_test":"Analytic rate-equation check: for fixed local baryon densities, write the Boltzmann collision term for n → χ + φ using the vacuum dark-decay width implied by the beam-bottle anomaly, the final-state Pauli factor (1 − f_χ), and no inverse term (n_φ ≈ 0). Integrate over a neutron-star lifetime at representative densities (n_b ≈ 0.3–3 n0) and compare the resulting n_χ with the solution of Eq. (17). If the kinetic n_χ differs by more than ~20% at any density in this range, recompute the TOV and tidal-deformability curves using the kinetic n_χ(r); if the maximum masses or mass-gap branches change by more than the observational bands, the central compatibility claim fails unless the equilibrium assumption is separately justified.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The entire EoS is built on Eq. (17), μ_n = μ_χ, which treats the dark decay n ↔ χ + φ as reversible with μ_φ = 0. The paper explicitly states (Secs. II and III.C) that φ escapes from the star and does not contribute to the energy density or pressure. In that open-system limit the reverse process χ + φ → n is suppressed because the φ abundance is negligible, so detailed balance does not enforce the equality. The forward decay is instead stopped by Pauli blocking of the χ Fermi sea, which at T=0 only requires μ_χ ≥ μ_n, not equality; the final n_χ is then history-dependent and not fixed uniquely by thermodynamic potentials. Consequently the composition fractions in Fig. 4 and the M-R/Λ curves in Figs. 1–3 are conditional on an equilibrium assumption that the paper itself flags as unresolved: the abstract's final sentence and concluding remark (e) state that environmental dependence of the decay is an open problem. If the actual steady-state χ population differs from the Eq. (17) value, the claimed sensitivity to η and the mass-gap objects could shift substantially.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the astrophysical implications of the neutron dark-decay hypothesis n → χ + φ inside neutron stars, focusing on the role of the nuclear symmetry energy. The hadronic sector is described by a low-order Taylor expansion of the energy density around saturation (Eqs. 5–9), parameterized by η = (K0 L^2)^{1/3}. The dark sector is treated as a Fermi gas with repulsive self-interactions and, in some variants, repulsive baryon–dark matter interactions. The composition is determined by the chemical equilibrium conditions μ_n = μ_χ and μ_n = μ_p + μ_e (Eq. 17), with the dark boson φ escaping freely. The resulting EoS is used to compute mass–radius curves and tidal deformabilities. The main claims are that the symmetry-energy parameter strongly affects the EoS and observables, that appropriate dark-sector repulsion can keep the model compatible with 2 M⊙ pulsars, and that baryon–DM interactions alone can produce mass-gap compact objects. The paper concludes that the dark-decay channel cannot be excluded by current astrophysical constraints.","tokens_in":16146,"tokens_out":19216,"duration_ms":235387,"significance":"If the model assumptions are accepted, this is a useful systematic study: it extends earlier pure-neutron-matter treatments to β-stable matter with protons and electrons, separately considers dark self-interactions and baryon–dark interactions, and uses standard TOV and tidal-deformability machinery. The finding that baryon–DM repulsion can produce mass-gap objects is interesting and connects to recent two-fluid models. The paper is also transparent that the environmental dependence of the decay is an open problem. The main limitations are that the composition and all derived observables rest on an unvalidated equilibrium saturation assumption, that the EoS is a low-order expansion extrapolated to very high densities, and that the parameter η entangles the symmetry energy with the symmetric-matter incompressibility.","major_comments":[{"comment":"The equality μ_n = μ_χ fixes nχ and, through Eq. (20), the entire EoS and all M–R/Λ results. The paper states that φ escapes and does not contribute to the EoS; in this open-system limit the reverse process χ + φ → n is suppressed, so Eq. (17) is not a detailed-balance chemical equilibrium. It can be interpreted as the Pauli-blocking endpoint of the one-way decay, but that requires the decay to have run to saturation inside the star and to be unaffected by the medium. The abstract and Concluding remark (e) explicitly leave the environmental dependence open. Since Fig. 4 and Figs. 1–3, 5 are conditional on this saturation assumption, the conclusions should be re-framed as a particular scenario, or the sensitivity to a partially populated χ Fermi sea should be quantified.","section":"§III.C, Eq. (17)"},{"comment":"The hadronic EoS is a Taylor expansion around n0 truncated at second order in (n−n0) and first order in S(n). It is used up to total densities n_t ≈ 1.5 fm⁻³ (Fig. 4), where (n−n0)/n0 ≈ 8. At these densities the omitted Ksym term in Eq. (6) is not necessarily small: with |Ksym| ≈ 100 MeV, its contribution to the energy density is of order several hundred MeV fm⁻³. The quantitative claims about the η dependence of Mmax, R1.4, and Λ1.4 are therefore not robust to the truncation. I recommend either restricting the analysis to densities where the expansion is controlled or benchmarking against a more complete EoS.","section":"§III.A, Eqs. (5)–(9)"},{"comment":"The parameter η = (K0 L^2)^{1/3} is varied by increasing K0 and L simultaneously (Table I). Thus the effects attributed to the nuclear symmetry energy are entangled with the stiffness of symmetric nuclear matter through K0. The abstract's central claim that the symmetry energy critically shapes the EoS is not isolated by this design. To support the title and abstract, the authors should vary L at fixed K0, or otherwise decorrelate the symmetry-energy slope from the incompressibility, at least for representative cases.","section":"Table I and §III.A"}],"minor_comments":[{"comment":"There are numerous typographical errors: 'defomability' in the Fig. 1 caption, 'amnd' in Ref. [19], 'e.t.c.' in Concluding remark (e), a missing initials in Ref. [27], and 'Tanjia Hinderer' in Ref. [57]. Please proofread carefully.","section":"General"},{"comment":"The legends mix 'with DM' and 'without DM' curves for ten values of η, making individual curves difficult to distinguish. The observational shaded regions are described only in the captions. Consider separating the panels or using distinct line styles, and include a table of Mmax, R1.4, and Λ1.4 values.","section":"Figs. 1–3"},{"comment":"The text says both the dark boson and the neutrino escape. While the neutrino condition is standard, the escaping dark boson means the star is not in global thermodynamic equilibrium; the local equilibrium used here should be clearly stated as a steady-state assumption rather than an equilibrium of the closed system.","section":"§III.C"},{"comment":"The mass bounds in Eqs. (3)–(4) are stated allowing for nonzero mφ, but the calculations set mφ = 0. A brief sentence explaining that mφ = 0 is consistent with the quoted bounds would avoid confusion.","section":"§II"}],"recommendation":"major_revision","confidential_remarks":"The paper is a coherent parameter study, but the central claim about the role of the symmetry energy is weakened by the joint variation of K0 and L in η, and the equilibrium condition in Eq. (17) needs a clear physical justification or a sensitivity test. The authors are transparent about the environmental uncertainty, which is to their credit. I believe the paper could be publishable after major revision, but the conclusions need to be re-scoped to the assumptions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear [Colleague],\n\nWhat you should know: this is a careful, systematic parameter study of neutron dark decay in beta-stable neutron star matter, and its main conclusions are conditional on an equilibrium condition that the authors assume rather than derive. It deserves a referee, but any use of its mass-radius or tidal-deformability curves should carry the caveat that μ_n = μ_χ in Eq. (17) is an imposed condition.\n\nThe genuinely new part is the scan over the symmetry-energy combination η = (K0 L^2)^{1/3} together with dark-matter self-interaction and baryon-dark-matter interaction strengths in beta-stable matter. That combination is not in the prior literature. The paper shows how η affects stellar radii and tidal deformability, and how a strong baryon-dark-matter interaction can stiffen the EoS enough to produce mass-gap objects. The TOV and Love-number machinery is standard, and the comparison with observations is qualitative rather than fitted. The authors also explicitly flag the environmental uncertainty in the decay, which is honest.\n\nThe soft spots are real. The load-bearing equation is μ_n = μ_χ. If the dark boson φ escapes, the reverse process χ + φ → n is suppressed, so detailed balance does not enforce equality. At T=0 the forward decay only stops when the χ Fermi sea blocks it, which requires μ_χ ≥ μ_n, not equality. The χ abundance, and hence the whole EoS, becomes history-dependent. The paper acknowledges this as an open problem in the abstract and in remark (e), but then proceeds as if equilibrium were guaranteed. That makes every M-R and Λ prediction conditional. I would ask the authors to either justify the equilibrium (e.g., through a trapped φ background) or present the results explicitly as a bracketing scenario.\n\nTwo smaller issues. The symmetry-energy expansion in Eq. (6) drops Ksym and is extrapolated to densities near 1.5 fm^-3; that is a standard low-order parameterization but could shift the quantitative radii and Λ values. And the notation is confusing: m_φ is used both for the massless decay product and for the massive mediator in the self-interaction. Clean that up.\n\nThe final conclusion—that this decay 'cannot be excluded' by mass constraints—should be softened to 'cannot be excluded within this equilibrium model.' As written, it promises more than the calculation delivers.\n\nWho gets value: nuclear astrophysicists working on dark matter in neutron stars or symmetry-energy constraints. This is not a landmark, but it is a reproducible scan, and the equilibrium caveat is an important cautionary tale.\n\nMy recommendation: send it to peer review, with the explicit request that the authors confront the μ_n = μ_χ issue and tone down the conclusion. With those changes it would be a solid contribution.\n\nBest.","headline":"Systematic, honest sensitivity scan of neutron dark decay in beta-stable stars, but the assumed μ_n=μ_χ equilibrium makes all M-R/Λ predictions conditional on an unproven condition.","tokens_in":16644,"tokens_out":4702,"would_cite":true,"duration_ms":51225,"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 argues that the nuclear symmetry energy and repulsive dark-sector interactions together determine whether neutron dark decay can occur inside neutron stars without violating observed mass constraints.","keywords":["neutron dark decay","neutron star equation of state","nuclear symmetry energy","dark matter self-interaction","beta-stable matter","tidal deformability","mass-gap compact objects","neutron lifetime anomaly"],"falsifier":"A decisive test would be a high-precision joint measurement of the radius and dimensionless tidal deformability Λ of a 1.4-solar-mass neutron star (for example from a loud binary neutron star merger). If the measured point falls outside the envelope of all curves the paper generates across its full range of η and dark-sector coupling strengths, the symmetry-energy-tuned dark-decay scenario is falsified. Equally decisive would be evidence that the neutron dark decay rate is suppressed or enhanced in dense matter, since that would break the chemical-equilibrium assumption on which the entire cal","tokens_in":15666,"feed_emoji":"⭐","tokens_out":7268,"duration_ms":79615,"temperature":0.7,"pith_summary":"This paper asks whether the proposed dark decay of neutrons—an extra decay channel n→χ+φ put forward to explain the beam-versus-bottle neutron lifetime discrepancy—can happen inside neutron stars without contradicting observed masses. The authors find that the answer hinges on the nuclear symmetry energy, the quantity governing how the energy of asymmetric nuclear matter rises with neutron excess. For beta-stable matter with neutrons, protons, electrons, and equilibrated dark fermions, the symmetry-energy parameter η=(K0L^2)^(1/3) sets how stiff the equation of state is. When dark particles repel each other strongly, or when baryon-dark matter repulsion is added, the equation of state can stay stiff enough to support two-solar-mass neutron stars, and in some cases even objects in the mass-gap region. The paper concludes that the dark decay channel cannot be excluded on the basis of current astrophysical mass constraints, although the environmental dependence of the decay inside dense matter remains an open problem.","feed_headline":"Neutron dark decay survives in stars with tuned symmetry energy","feed_subtitle":"Tuning the nuclear symmetry energy and dark-sector repulsion keeps 2-solar-mass neutron stars compatible with the decay channel.","key_machinery":"The key machinery is an equation of state built from a total energy density E_tot(n_n, n_p, n_χ, n_e) that couples a parabolic nuclear-matter ansatz—whose stiffness is set by η=(K0L^2)^(1/3), where K0 is the incompressibility and L the symmetry-energy slope—to a gas of repulsively interacting dark fermions described by Yukawa-type self- and cross-interactions. Chemical equilibrium conditions μ_n=μ_χ and μ_n=μ_p+μ_e, together with charge neutrality, fix the particle fractions; the pressure follows from the Gibbs-Duhem relation. Solving the TOV equations for these equations of state yields mass-radius curves and tidal deformabilities. The η parameter is the knob connecting finite-nucleus input","core_discovery":"The core discovery is that, unlike several earlier studies that considered only pure neutron matter, the neutron dark decay scenario is not automatically lethal for neutron stars. Treating the star as beta-stable npe matter plus an equilibrated gas of dark fermions from n→χ+φ, and parameterizing the nuclear equation of state by η=(K0L^2)^(1/3), the authors show that the symmetry energy (through L and K0) and the repulsive interactions within the dark sector and between dark and baryonic matter act together to set the particle fractions and total pressure. For strong dark self-interactions, the dark particles appear mainly at high density and leave stellar structure nearly unchanged; for weak","pith_inferences":["If the neutron dark decay rate is density- or pressure-dependent—which the paper explicitly leaves open—the chemical equilibrium assumption would break down, and the predicted mass-radius and tidal deformability curves would shift; a measurable cooling anomaly or a dark-fraction signal in neutron stars could probe this directly.","The paper treats the dark and baryonic components as a single fluid; a two-fluid treatment could yield different radii and tidal signatures, so merger waveforms might be able to distinguish the two descriptions.","Direct measurements of the symmetry-energy slope L (for instance from neutron-skin experiments) would narrow the η band and sharpen the predicted relationship between dark-sector coupling and observable neutron-star properties.","The mass-gap objects produced in the pure baryon-dark matter interaction case would have unusually large radii and high tidal deformabilities, offering a way to distinguish them from black holes in gravitational-wave events."],"forward_implications":["If the dark decay channel is real, neutron star equations of state should be evaluated with beta-stable matter rather than pure neutron matter, because the symmetry energy changes the dark-particle fraction and hence the stiffness.","Strongly repulsive dark-matter self-interactions (small zχ) keep the mass-radius curves nearly identical to the no-dark-matter case, so the scenario remains compatible with observations such as GW170817 and the high-mass pulsar measurements.","Adding repulsive baryon-dark matter interactions, without self-interaction, can stiffen the equation of state enough to produce objects in the mass-gap region, giving a concrete signature for future searches.","The tidal deformability Λ at 1.4 solar masses spans up to two orders of magnitude as η and the dark-sector couplings are varied, making it the most sensitive observable for distinguishing or constraining the dark-decay scenario.","The dark decay channel cannot currently be excluded by mass constraints alone; only more precise radius or tidal measurements can tighten the allowed parameter space."],"supporting_citations":[{"why":"Proposes the n→χ+φ dark decay channel as an explanation of the neutron lifetime anomaly; this is the scenario under study.","marker":"[7]"},{"why":"Supplies the Yukawa interaction parametrization and the ranges of dark self- and baryon-dark coupling parameters used here.","marker":"[9]"},{"why":"Shows dark matter-baryon repulsion can support two-solar-mass neutron stars; the direct baseline this paper extends to beta-stable matter.","marker":"[15]"},{"why":"Establishes minimum dark-fermion mass constraints from the existence of low-mass neutron stars; used to set mχ=938 MeV.","marker":"[16]"},{"why":"Introduces η=(K0L^2)^(1/3) as a single stiffness regulator for the nuclear equation of state; this parameter is the paper's symmetry-energy knob.","marker":"[44]"},{"why":"Provides the specific K0-L-η grid and the chemical-equilibrium treatment adopted for the hadronic equation of state.","marker":"[46]"},{"why":"Gives the repulsive self-interaction energy density for dark fermions used to build the dark-matter equation of state.","marker":"[47]"},{"why":"The GW170817 constraint on Λ at 1.4 solar masses used to check whether the dark-decay equations of state are observationally acceptable.","marker":"[59]"}],"fun_headline_variants":["Symmetry energy key to neutron dark decay survival","Dark repulsion lets neutron stars endure dark decay","Beta-stable matter opens door to neutron dark decay","Tuned nuclear forces keep neutron decay from killing stars"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The calculation assumes the dark fermion χ produced by neutron decay reaches chemical equilibrium with neutrons (μ_n = μ_χ) and forms a single equilibrated fluid filling the star; if the decay rate depends on density, pressure, or magnetic field—which the paper explicitly flags as an open problem—the equilibrium composition and all resulting mass-radius curves no longer describe the star.","fun_headline_variants_meta":{"raw":{"variants":["Symmetry energy key to neutron dark decay survival","Dark repulsion lets neutron stars endure dark decay","Beta-stable matter opens door to neutron dark decay","Tuned nuclear forces keep neutron decay from killing stars"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000277,"raw_usage":{"total_tokens":1478,"prompt_tokens":728,"completion_tokens":750,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":472,"completion_tokens_details":{"reasoning_tokens":689}},"tokens_in":472,"tokens_out":750,"duration_ms":9150,"temperature":1.0,"reasoning_tokens":689,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T13:58:26.253916+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive test would be a high-precision joint measurement of the radius and dimensionless tidal deformability Λ of a 1.4-solar-mass neutron star (for example from a loud binary neutron star merger). If the measured point falls outside the envelope of all curves the paper generates across its full range of η and dark-sector coupling strengths, the symmetry-energy-tuned dark-decay scenario is falsified. Equally decisive would be evidence that the neutron dark decay rate is suppressed or enhanced in dense matter, since that would break the chemical-equilibrium assumption on which the entire cal","supporting_citations":[{"cited_title":"Fornal and B","cited_arxiv_id":null,"evidence_quote":"Proposes the n→χ+φ dark decay channel as an explanation of the neutron lifetime anomaly; this is the scenario under study."},{"cited_title":"Bastero-Gil, T","cited_arxiv_id":null,"evidence_quote":"Supplies the Yukawa interaction parametrization and the ranges of dark self- and baryon-dark coupling parameters used here."},{"cited_title":"Grinstein, C","cited_arxiv_id":null,"evidence_quote":"Shows dark matter-baryon repulsion can support two-solar-mass neutron stars; the direct baseline this paper extends to beta-stable matter."},{"cited_title":"McKeen Ann E","cited_arxiv_id":null,"evidence_quote":"Establishes minimum dark-fermion mass constraints from the existence of low-mass neutron stars; used to set mχ=938 MeV."},{"cited_title":"Sotani, K","cited_arxiv_id":null,"evidence_quote":"Introduces η=(K0L^2)^(1/3) as a single stiffness regulator for the nuclear equation of state; this parameter is the paper's symmetry-energy knob."},{"cited_title":"Divaris, A","cited_arxiv_id":null,"evidence_quote":"Provides the specific K0-L-η grid and the chemical-equilibrium treatment adopted for the hadronic equation of state."},{"cited_title":"Nelson, S","cited_arxiv_id":null,"evidence_quote":"Gives the repulsive self-interaction energy density for dark fermions used to build the dark-matter equation of state."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The GW170817 constraint on Λ at 1.4 solar masses used to check whether the dark-decay equations of state are observationally acceptable."}],"review_version":1}