{"id":"3cb51ea4-77ba-4a1a-b312-d0460689f3b0","arxiv_id":"2608.10141","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Scalar baryons made inside neutron stars require non-perturbatively large repulsive self-couplings, lambda_4 greater than 1000, to permit two-solar-mass stars, and attractive self-interactions make even scalars heavier than the core chemical potential testable.","lead":"Neutron stars can convert neutrons into hypothetical scalar particles that carry baryon number, because the core's chemical potential can exceed the particle masses. This paper maps the scalar masses and self-couplings that are still allowed by the observed two-solar-mass neutron stars.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central exclusion curves assume n→φ chemical equilibrium without any rate calculation; if conversion is slower than the neutron-star lifetime, the m_φ>m_n constraints and the λ_4>1000 claim are not established.","rationale":"The reader and I identify the same weakest link: the chemical-equilibration assumption. I do not find an internal inconsistency in the EoS construction: the monomial and polynomial EoS are derived consistently from the Lagrangian, the TOV setup follows standard enthalpy variables, and the paper provides data tables. The attractive Q-matter section is unusual but internally consistent. The central problem is that the equilibration assumption is both unquantified and consequential. It is stated at the outset and in Sec. IV, so the paper is not misleading, but the abstract's broad claim that non-perturbatively large repulsive self-couplings are required is only valid for models with fast n→φ conversion. A missing rate calculation is the difference between a constraint on all perturbative scalar baryons and a constraint on a subclass that equilibrates. The proposed rate/kinetic test would settle whether this concern actually lands, so the verdict should remain CONDITIONAL as the reader already concluded.","tokens_in":12304,"tokens_out":23579,"duration_ms":237926,"concrete_test":"Choose a benchmark portal y nνφ* and compute the n→φ rate in degenerate neutron-star matter at T≈0 with Fermi-surface phase space and Pauli blocking, using μ_B and Eq. (17); compare the equilibration time to 10^9 yr. If the produced φ abundance stays below the equilibrium abundance, replace chemical equilibrium by a kinetically fixed φ density in the TOV integration; if M_max remains above 2M_⊙, the exclusion regions in Figs. 2-5 do not apply.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The quantitative claims rest on Eqs. (28)-(29), which add the φ EoS to nuclear matter under the condition μ_φ=μ_B. This condition is imposed, not derived. In Sec. IV the authors state that the portal couplings 'are assumed to be sufficiently efficient to induce n→φ transitions over the lifetime of the neutron star and equilibrate baryons,' and the same assumption appears for the dark fermion in Sec. III. No estimate of the conversion rate is given. For m_φ>m_n the vacuum decay n→φν is kinematically forbidden, so the in-medium rate and the approach to equilibrium depend on the unspecified portal coupling, Fermi-surface phase space, and Pauli blocking. A slow rate would mean no condensate forms, the EoS returns to the nuclear one, and Figs. 2-5 place no constraint for m_φ>m_n. The m_φ<m_n region is independently excluded by nucleon decay bounds, so the meaningfully new constraints are for m_φ>m_n; those are exactly the ones that disappear if equilibration is slow. The paper is explicit that this is an assumption, so this is a limitation rather than an internal inconsistency, but it is the load-bearing step for the claimed exclusion of scalar baryons up to roughly 1.5 GeV.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper investigates the impact of a GeV-scale complex scalar field φ carrying baryon number B=1 on the equation of state (EoS) of neutron stars. Under the assumption of chemical equilibrium between the scalar and neutrons, the authors derive the condensate EoS for repulsive monomial self-interactions (λ_N|φ|^N), for polynomial potentials, and for attractive quartic-plus-sextic potentials. Solving the Tolman-Oppenheimer-Volkoff equations with three CompOSE nuclear EoS, they find that supporting a 2M_⊙ neutron star requires λ_4 ≳ 1000 for quartic interactions when m_φ ≲ m_n, and even larger effective couplings for sextic interactions. For attractive potentials, scalar production can occur for m_φ > μ_B^max, yielding constraints near the Q-matter threshold. The paper benchmarks its TOV integration against published CompOSE mass-radius curves and provides data tables as ancillary files.","tokens_in":12473,"tokens_out":7888,"duration_ms":74481,"significance":"If the chemical-equilibrium assumption holds, the paper delivers a clean, analytic, and falsifiable set of constraints on scalar baryons that complement nucleon-decay searches. The derivation of the condensate EoS is transparent, the TOV solver is benchmarked, and the dependence on three nuclear EoS gives an honest spread of results. The conclusion that perturbative scalar baryons are excluded up to ~1.5 GeV is striking and would be an important addition to the dark-baryon literature. However, the new constraints for m_φ > m_n — precisely the region not already excluded by nucleon decays — rest entirely on the unquantified assumption of fast n↔φ equilibration.","major_comments":[{"comment":"The central exclusion for m_φ > m_n in Figs. 2–5 assumes that the portal couplings nνφ* and peφ* are sufficiently efficient to keep φ in chemical equilibrium over the neutron-star lifetime, but no estimate of the rate is provided. For m_φ > m_n the vacuum decay is kinematically forbidden, so the in-medium rate is controlled by the unspecified portal coupling and by Pauli blocking/phase-space factors. If the rate is slower than the star's age, the φ condensate never forms and the EoS returns to the nuclear one, in which case the m_φ > m_n exclusion curves disappear. Since the m_φ < m_n region is independently excluded by nucleon-decay searches, the genuinely new constraints are exactly the ones at risk. Please quantify the minimum rate required for equilibration (e.g., relative to the inverse neutron-star age) and give a representative order-of-magnitude estimate for the portal coupling that would achieve it, or explicitly state how the exclusions should be interpreted if equilibration is not guaranteed.","section":"Sec. IV, after Eq. (14) and before Eq. (15)"},{"comment":"The constraints for m_φ > μ_B^max in the fine-tuned region λ_4^2 ≃ 4λ_6 rely on the homogeneous Q-matter approximation, with the statement that 'a more careful treatment would keep the surface energy that was neglected here.' Because this is the only region that probes arbitrarily large scalar masses, the surface-energy correction could affect the EoS at low pressure and hence the maximum mass. Please estimate the magnitude of the surface-energy contribution (e.g., the surface tension of the Q-ball) and show that it does not shift the exclusion boundary, or qualify the constraint in this region.","section":"Sec. IVC, around Eq. (42)"}],"minor_comments":[{"comment":"The captions display 'Vint = λ_4 |ϕ 4' and 'λ_6/m_ϕ^2 |ϕ 6'; these should be '|ϕ|^4' and '|ϕ|^6' respectively, with the missing absolute-value bars and exponents.","section":"Fig. 2 and Fig. 3 captions"},{"comment":"The notation for the maximal baryon chemical potential is inconsistent: μ_B^max appears in Sec. IVA while μ_max^B is used in Sec. IVC; please unify.","section":"Secs. IVA and IVC"},{"comment":"The abstract states that 'non-perturbatively large repulsive self-couplings are required,' but this statement applies primarily to m_φ ≲ m_n; for m_φ > m_n the required couplings are smaller and the μ_φ threshold is mass-dependent. Please specify the mass range in the abstract to avoid overgeneralization.","section":"Abstract and Sec. IVA"},{"comment":"The paper cites only PSR J1614-2230 and PSR J0348+0432 as 2M_⊙ pulsars. More recent precision measurements (e.g., PSR J0740+6620) would strengthen the constraints; at minimum, a comment on the robustness of the 2M_⊙ threshold would be useful.","section":"Sec. I and Conclusions"},{"comment":"The analytic estimates 'λ_3 ∼200, λ_4 ∼4000, λ_6 ∼4×10^6' are a factor of a few above the numerical λ_4 threshold of 1000. Please state explicitly that these are order-of-magnitude estimates and not exact thresholds.","section":"Sec. IVA, text before Eq. (28)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a well-written letter that makes a strong claim. The main concern is the unquantified equilibration assumption, which is explicitly acknowledged but is load-bearing for the new constraints. I would ask the authors to add a quantitative rate discussion or to reframe the conclusions as conditional. The paper fits the journal's scope and is likely to be of significant interest."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe short version: this paper gives the first clean treatment of an equilibrated baryon-number-carrying scalar condensate in neutron stars, with a tidy analytic EoS and honest benchmarked TOV integrations. But the exclusion curves in Figs. 2-5 only bite if n->phi conversion is fast enough to reach chemical equilibrium over the star's lifetime, and the paper never estimates that rate. Treat the headline lambda_4 > 1000 claim as conditional.\n\nWhat's actually new: the cited fermionic dark-baryon literature doesn't cover scalars, and the Bose-condensate EoS with self-interactions is a natural, well-executed extension. The analytic expressions for the monomial and polynomial potentials are clean, and the attractive-interaction section identifying the Q-matter regime is a nice touch, showing that even m_phi above the core chemical potential can be constrained near the fine-tuned boundary. The TOV solver is benchmarked against CompOSE curves, and the paper provides data tables, which is reproducible enough for a first pass.\n\nThe soft spot is the equilibration assumption. For m_phi < m_n the region is already dead by nucleon-decay bounds, so the genuinely new constraints are for m_phi > m_n. There the vacuum decay is kinematically forbidden, and the in-medium rate depends on the unspecified portal coupling, phase-space, and Pauli blocking. If the rate is slower than ~1e-9 per year, no condensate forms and the EoS is just nuclear; the exclusions evaporate. The authors state this explicitly, so it's a clearly flagged limitation rather than an internal inconsistency, but it is load-bearing.\n\nA minor second issue: the Thomas-Fermi approximation neglects gradients and surface energy, which matters most in the Q-matter regime where they claim constraints on arbitrarily heavy scalars. That's not central to the quartic result but worth a sentence if the paper is revised.\n\nWho's it for: anyone working on dark baryons, neutron decays, or neutron-star probes of new physics. It's a competent letter-format contribution with a clean formalism; the main open question is whether the equilibrium assumption can be justified. I'd send it to a serious referee rather than desk-reject, because an expert can quickly judge whether a rate estimate is feasible, and the formalism itself is worth publishing. If I were refereeing, I'd ask for a concrete portal model and a rate estimate for a representative coupling, even if order-of-magnitude.","headline":"Clean, well-executed extension of dark-baryon neutron-star constraints to scalar condensates, but the central exclusion curves are conditional on an unevaluated n->phi equilibration rate.","tokens_in":13129,"tokens_out":2247,"would_cite":true,"duration_ms":21007,"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":"A baryon-number-carrying scalar produced inside neutron stars needs self-coupling $\\lambda_4 \\gtrsim 1000$ to allow two-solar-mass stars, excluding perturbative scalar baryons up to $\\sim 1.5$ GeV.","keywords":["scalar baryons","neutron stars","baryon chemical potential","equation of state","self-interacting scalar","Q-balls","Tolman-Oppenheimer-Volkoff","dark baryons"],"falsifier":"Compute the in-medium $n\\to\\phi\\,\\bar{\\nu}$ conversion rate from a specified coupling; if the resulting timescale exceeds the roughly $10^9$-year age of the observed neutron stars (rate below about $10^{-9}\\,\\mathrm{yr}^{-1}$), chemical equilibrium is never reached and the $\\lambda_4\\gtrsim 1000$ exclusion curves for $m_\\phi>m_n$ cease to apply. A positive test would be a neutron-star observation whose mass-radius curve requires the stiff $\\lambda_4=1000$ equation of state predicted here.","tokens_in":11995,"feed_emoji":"⭐","tokens_out":13896,"duration_ms":115316,"temperature":0.7,"pith_summary":"This paper asks whether a GeV-scale complex scalar field carrying baryon number could hide inside neutron stars, produced by converting neutrons into scalars in the dense core where the baryon chemical potential exceeds the neutron mass. The authors show that adding such a scalar to the star's equation of state removes the neutrons' Fermi pressure and softens the star; to still reach the observed two-solar-mass maximum mass, the scalar's repulsive quartic self-coupling must be non-perturbatively large, $\\lambda_4 \\gtrsim 1000$ for scalar masses below the neutron mass. For sextic self-interactions the required effective coupling is even larger, above $10^6$. Attractive self-interactions of the Q-ball type can instead produce scalars even when their vacuum mass exceeds the chemical potential, constraining scalar baryons up to arbitrarily large masses in a fine-tuned region. If correct, neutron-star observations exclude perturbative scalar-baryon models with masses up to the roughly 1.5 GeV core chemical potential.","feed_headline":"To survive neutron stars, scalar baryons need λ ≳ 1000","feed_subtitle":"The two-solar-mass pulsars demand repulsive self-couplings far beyond the perturbative range.","key_machinery":"The load-bearing object is a zero-temperature complex scalar condensate in chemical equilibrium with neutrons, treated in the local-density (Thomas-Fermi) approximation. The equations $\\mu_\\phi^2 = m_\\phi^2 + V_{\\rm int}'(X)$, $n_\\phi = 2\\mu_\\phi X$, $\\varepsilon_\\phi = (\\mu_\\phi^2 + m_\\phi^2)X + V_{\\rm int}(X)$, and $P_\\phi = (\\mu_\\phi^2 - m_\\phi^2)X - V_{\\rm int}(X)$, with $\\mu_\\phi = \\mu_B$, give the scalar's pressure and energy density as functions of the baryon chemical potential. These are added to a nuclear equation of state, and the combined equation of state is integrated through the Tolman-Oppenheimer-Volkoff equations using the enthalpy variable of Ref. [47] to produce mass-radius curves whose maximum mass must reach the observed $2\\,M_\\odot$. Monomial, polynomial, and attractive potentials are treated in turn; the attractive case is the homogeneous limit of the non-topological solitons known as Q-balls, and its $\\varepsilon(P=0)>0$ behaviour is the hallmark of self-bound Q-matter.","core_discovery":"The central claim is that chemical equilibrium between neutrons and a complex scalar baryon $\\phi$ with $B(\\phi)=1$ sets in inside neutron stars once the baryon chemical potential $\\mu_B$ reaches the scalar's effective production threshold: $\\mu_B > m_\\phi$ for repulsive self-interactions, or $\\mu_B > \\mu_0 < m_\\phi$ when a negative $|\\phi|^4$ term provides binding energy. The scalar then forms a homogeneous zero-temperature condensate described by $\\mu_\\phi^2 = m_\\phi^2 + V_{\\rm int}'(X)$ and $n_\\phi = 2\\mu_\\phi X$, with the thermodynamic identity $\\varepsilon_\\phi + P_\\phi = \\mu_\\phi n_\\phi$. Because the scalar is bosonic, it contributes no Fermi pressure; without repulsive self-interactions it would condense into an almost pressureless state and drive the maximal neutron-star mass below one solar mass. Demanding $M_{\\rm max} \\ge 2\\,M_\\odot$ therefore leaves only non-perturbative repulsive couplings: $\\lambda_4 \\gtrsim 1000$ for a quartic $|\\phi|^4$ interaction when $m_\\phi < m_n$, and $\\lambda_6/m_\\phi^2 \\gtrsim 10^6$ to $10^8\\,\\mathrm{GeV}^{-2}$ for a sextic $|\\phi|^6$ interaction, depending on the nuclear equation of state. An attractive quartic term shifts production to $\\mu_0 < m_\\phi$, so even scalars heavier than the maximum core chemical potential are constrained in the fine-tuned region $\\lambda_4^2 \\simeq 4\\lambda_6$.","pith_inferences":["The equilibrium assumption is quoted but not derived; a first-principles computation of $n\\to\\phi$ conversion in dense matter is the single most decisive check, and the exclusion figures should be read as conditional on it.","The same condensate framework implies radial density profiles and cooling signatures beyond mass-radius; those are not computed here and could be probed with independent neutron-star observations.","Stronger empirical statements, such as a measured gravitational-wave tidal deformability, could turn the '$\\lambda_4\\gtrsim 1000$' requirement into a direct test of the scalar potential shape rather than just its overall scale."],"forward_implications":["Perturbative scalar-baryon models with $B=1$ and $m_n < m_\\phi < \\mu_B^{\\rm max}\\sim 1.5$ GeV are excluded unless $\\lambda_4 \\gtrsim 1000$ (quartic) or $\\lambda_6/m_\\phi^2 \\gtrsim 10^6$ GeV$^{-2}$ (sextic).","The region $m_\\phi < m_n$ is independently eliminated by fast neutron-decay bounds once conversion is fast enough to equilibrate, so the new bounds bite in the window $m_n < m_\\phi < \\mu_B^{\\rm max}$.","A scalar-baryon explanation of the neutron-lifetime anomaly through $n\\to\\phi\\bar{\\nu}$ requires additional long-range repulsion beyond the quartic self-interaction to survive neutron-star constraints.","Attractive $|\\phi|^4+\\lambda_6|\\phi|^6$ potentials allow neutron-star constraints to reach scalar masses above $\\mu_B^{\\rm max}$ near $\\lambda_4^2 \\simeq 4\\lambda_6$, where the condensate is Q-matter and contributes less than about 1% of the star's mass.","For a $B=2$ scalar $\\xi$ coupled to two neutrons, the same production logic applies with threshold $2\\mu_B^{\\rm max}\\sim 3$ GeV."],"supporting_citations":[{"why":"Observed two-solar-mass neutron star J1614-2230, the mass anchor the computed curves must exceed.","marker":"[15]"},{"why":"Observed two-solar-mass neutron star J0348+0432, the second mass anchor.","marker":"[16]"},{"why":"Establishes the dark-fermion analogue: a baryon-number-carrying particle produced in the core relieves Fermi pressure and lowers the maximal neutron-star mass.","marker":"[13]"},{"why":"Gives the parallel neutron-star bound for dark baryons that this scalar analysis extends.","marker":"[14]"},{"why":"Provides the enthalpy reparametrization of the TOV equations used in the numerical integration.","marker":"[47]"},{"why":"Supplies the homogeneous complex-scalar condensate equations for energy density and pressure from which the scalar EoS is built.","marker":"[48]"},{"why":"Introduces Q-balls and Q-matter, the soliton picture underlying the attractive-potential condensate.","marker":"[53]"}],"fun_headline_variants":["Scalar baryons need λ ≳ 1000 to survive neutron stars","Neutron stars force scalar baryons to self-repel strongly","Two-solar-mass stars demand repulsive scalar baryon interactions","Scalar baryons in neutron stars: repulsion required for 2 Msun"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire exclusion picture assumes that neutrons and $\\phi$ reach chemical equilibrium inside the star; if the $n\\to\\phi$ conversion rate is slower than the neutron star's lifetime, the condensate never forms and the bounds for $m_\\phi>m_n$ disappear.","fun_headline_variants_meta":{"raw":{"variants":["Scalar baryons need λ ≳ 1000 to survive neutron stars","Neutron stars force scalar baryons to self-repel strongly","Two-solar-mass stars demand repulsive scalar baryon interactions","Scalar baryons in neutron stars: repulsion required for 2 Msun"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00087,"raw_usage":{"total_tokens":3792,"prompt_tokens":993,"completion_tokens":2799,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":609,"completion_tokens_details":{"reasoning_tokens":2720}},"tokens_in":609,"tokens_out":2799,"duration_ms":19407,"temperature":1.0,"reasoning_tokens":2720,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:11:31.638425+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the in-medium $n\\to\\phi\\,\\bar{\\nu}$ conversion rate from a specified coupling; if the resulting timescale exceeds the roughly $10^9$-year age of the observed neutron stars (rate below about $10^{-9}\\,\\mathrm{yr}^{-1}$), chemical equilibrium is never reached and the $\\lambda_4\\gtrsim 1000$ exclusion curves for $m_\\phi>m_n$ cease to apply. A positive test would be a neutron-star observation whose mass-radius curve requires the stiff $\\lambda_4=1000$ equation of state predicted here.","supporting_citations":[{"cited_title":"Phase Transitions and the Mass-Radius Curves of Relativistic Stars","cited_arxiv_id":"gr-qc/9802072","evidence_quote":"Provides the enthalpy reparametrization of the TOV equations used in the numerical integration."}],"review_version":1}