{"id":"ba35a234-abe5-43a6-a0e6-4159cd676327","arxiv_id":"2508.21121","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"A QCD-like dark sector with a theta vacuum and a dark photon portal can simultaneously reproduce the dark matter relic abundance and produce the sharply velocity-dependent self-interactions inferred from small-scale structure.","lead":"This paper proposes a dark matter model made of pions from a QCD-like hidden sector, where a non-zero topological angle and a dark photon portal let the pions both freeze out to the observed cosmic abundance and scatter off each other at galaxy-scale speeds. It is a concrete realization of an idea that could explain the dark matter density and small-scale structure puzzles, including the strong-lensing system SDSS J0946+1006, with a single mechanism.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The benchmark's resonant velocity v_R ~ 85 km/s is a tuned input: it requires m_η ≈ 2m_π0 at θ=0 to a fractional accuracy of ~10^-8, and the paper does not explain its origin; both the SIDM velocity dependence and resonant relic enhancement depend on this near-degeneracy.","rationale":"The reader's weakest_assumption correctly identifies the v_R tuning as the load-bearing point. The paper's central claim combines thermal relic production (§3) with velocity-dependent SIDM (§5). Both rely on the η resonance being near the two-pion threshold: the relic density requires resonant 3→2 processes, and the sharp velocity dependence in Eq. (5.2) is governed by v_R. If m_η is not tuned to 2m_π0 at the 10^-8 level, the Breit-Wigner peak moves out of the galactic velocity range and the cluster bound is no longer evaded; the relic density may still work for some broader resonances, but the specific 'natural explanation' of systems like SDSS J0946+1006 fails. The paper acknowledges this limitation explicitly ('we remain agnostic of the origin of values of v_R'), so this is not an internal inconsistency but a UV-completeness gap. My proposed NLO check is concrete and would quantify whether the tuning is even radiatively stable within the effective theory, making it more severe than a mere choice of parameters. Since this concern aligns with the reader's assessment and does not overturn the paper's technical derivations, the verdict remains CONDITIONAL; no adjustment is needed.","tokens_in":26118,"tokens_out":21550,"duration_ms":225689,"concrete_test":"Compute the one-loop (NLO) chiral perturbation theory correction to m_η − 2m_π0 for the benchmark parameters (m_π = 20 MeV, f_π = 34 MeV, r12 = 0.33, θ = 0.0005). If the correction exceeds ~0.4 eV, then the tree-level v_R condition is not radiatively stable, and the quark mass parameters must be retuned at the level of 10^5 or more. Alternatively, survey the UV models cited in Refs. [51–53] and verify whether any of them enforces m_η = 2m_π0 at θ=0 through a symmetry or other natural mechanism; if not, the tuned v_R is an unresolved gap in the model's UV completion.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central mechanism relies on the η meson being nearly degenerate with the two-pion threshold: v_R is defined in Eq. (2.21) as v_R = 2 sqrt((m_η − 2m_π0)/m_π0). The benchmark uses v_R = 85 km/s, which for m_π0 = 20 MeV gives m_η − 2m_π0 ≈ 0.4 eV, i.e., a fractional mass degeneracy of ~10^-8. The paper explicitly states in Sec. 2.2 that it 'remain[s] agnostic of the origin of values of v_R of the order of non-relativistic velocities' and imposes v_R|θ=0 = 0. This is a severe fine-tuning of the quark mass parameters (m1, m2, m3) and θ, not a derived consequence. Moreover, at one-loop order in chiral perturbation theory, the mass difference m_η − 2m_π0 receives corrections of order m_π^3/(16π^2 f_π^2) ~ 40 keV for the benchmark (m_π = 20 MeV, f_π = 34 MeV), roughly five orders of magnitude larger than the 0.4 eV needed for v_R = 85 km/s. Unless the underlying quark masses are adjusted to cancel these loop contributions, the resonant enhancement for both the 3→2 freeze-out (Sec. 3) and the Breit-Wigner self-interaction cross section in Eq. (5.2) is absent. The paper's central claim is therefore only established for a tuned benchmark, not for the generic θ-vacuum framework.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies dark pion dark matter in a QCD-like theory with a non-zero theta angle, adding a dark photon portal to the Standard Model. It solves coupled Boltzmann equations for π0, π± and η, including semi-annihilation, decays and coannihilations, and finds two families of parameters reproducing Ωh^2 ≈ 0.12. It then derives the DM decay lifetime through an O(p^6) operator, discusses thermalization and constraints from dark photon searches and indirect detection, and shows that the same η resonance that drives freeze-out produces a sharply velocity-dependent Breit-Wigner self-interaction cross section. The benchmark (mπ = 20 MeV, fπ = 34 MeV, r12 = 0.33, θ = 0.0005, vR = 85 km/s) simultaneously satisfies cluster bounds on σ/m at high velocities and has σ/m ≈ 150 cm^2/g at ~60 km/s, as invoked for SDSS J0946+1006. The central claim is that this framework unifies thermal relic production and SIDM phenomenology.","tokens_in":26551,"tokens_out":15247,"duration_ms":163621,"significance":"If the central mechanism holds, the paper provides a concrete, falsifiable model that connects the dark matter relic abundance to a sharp velocity-dependent self-interaction through the same η resonance. The Boltzmann treatment is explicit and internally consistent: rate densities, chemical-potential relations and a new solution family are given in detail, and the dark-photon portal is worked out with thermalization conditions and lifetime estimates. The benchmark yields a specific prediction for σ/m(v) and for the DM lifetime, which is a useful target for future simulations and indirect searches. The paper also correctly flags the absence of N-body simulations for Breit-Wigner self-interactions. However, the two main quantitative outputs—the resonant freeze-out and the resonant SIDM signal—depend on a very finely tuned η–2π0 degeneracy, and the DM lifetime estimate rests on an uncontrolled O(p^6) coefficient. Those issues currently limit the strength of the conclusions.","major_comments":[{"comment":"The central mechanism relies on the tuned relation v_R|θ=0=0. For v_R=85 km/s, Eq. (2.21) gives (m_η−2m_π0)/m_π0=(v_R/2)^2≈2×10^-8, a splitting of ~0.4 eV for m_π0=20 MeV; Sec. 2.2 remains agnostic about this. The resonance is load-bearing: both the 3→2 freeze-out in Sec. 3 and the Breit-Wigner peak in Eq. (5.2) disappear if the degeneracy is absent. In chiral perturbation theory the mass difference receives one-loop corrections ~m_π^3/(16π^2 f_π^2)≈40 keV, five orders above the required splitting, so maintaining v_R=85 km/s needs a cancellation to 10^-5 unless a symmetry protects it. Please provide an origin or counterterm arrangement for this coincidence, or quantify how the predictions degrade as v_R is varied; as written, the headline result is established only for a fine-tuned benchmark.","section":"Sec. 2.2 and Eqs. (2.21), (5.2)"},{"comment":"Viability depends on τ_DM being above the Voyager bound, τ≳3×10^26 s. Equation (4.4) is an estimate 'up to O(1) factors' based on a representative O(p^6) operator with coefficient set to unity. The benchmark star in Fig. 6 has τ_DM≈5×10^26 s, only about twice the bound. An O(1) coefficient of order 4π, or the presence of a second comparable operator, can easily invert the conclusion. Please provide a more controlled estimate of the operator coefficient (e.g., from resonance saturation or lattice input), or demonstrate that the same phenomenological conclusions hold with a conservative coefficient.","section":"Sec. 4.1, Eq. (4.4), Fig. 6"},{"comment":"The SIDM calculation assumes that 'all DM today consists of π0'. In charge assignment 2, however, the π± are stable under U(1)_d (Sec. 4.1), and they are converted to π0 only via π+π−→π0π0, a process that can itself freeze out. The paper does not report the residual π± abundance. A percent-level residual would change the DM mass normalization and the σ/m prediction in Eq. (5.2). Please quantify the late-time Yπ± for the benchmark and for the new diagonal solution branch, and confirm that the π0-only assumption is quantitatively accurate.","section":"Sec. 5.1 and Fig. 7"}],"minor_comments":[{"comment":"The caption contains a typo: 'T able 1' should read 'Table 1'. Also, 'times sin 2 θπη' is ambiguous; it should be written as sin^2 θπη (or otherwise clarified).","section":"Table 1"},{"comment":"The thermally averaged cross sections are given at leading order in the chiral and θ expansions. For the diagonal solution family, which reaches θ ~ 0.1 in Fig. 3, please state the numerical accuracy expected and the range of validity of these truncated expressions.","section":"Sec. 3.1, Eqs. (3.13)–(3.15)"},{"comment":"The text refers repeatedly to the 'vertical line' and 'diagonal line' in Fig. 3, but the figure itself does not label these branches. Adding labels would improve readability.","section":"Fig. 3"},{"comment":"The phrase in Sec. 5.2 and the Abstract that SDSS J0946+1006 'can be naturally explained' is stronger than the caveats in the same section, which note the lack of N-body simulations for Breit-Wigner self-interactions and the ongoing debate about the SIDM interpretation (Refs. [103,104]). Please soften the wording to reflect these uncertainties.","section":"Sec. 5.2"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope and technically solid in its Boltzmann treatment, but the main phenomenological claims are contingent on a very fine-tuned η–2π0 degeneracy and on an uncontrolled O(p^6) lifetime estimate. These are addressable in principle—by treating v_R as an input and scanning over it, or by providing a model-building justification, and by better estimating or marginalizing over the decay operator coefficient. I therefore recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nRead the arXiv:2508.21121 paper on pion DM in a theta vacuum. The headline: it's a technically serious extension of Paper I, but the genuinely new pieces are the dark photon portal, the stability analysis, and a new family of relic-density solutions—not the core mechanism or the velocity-dependent resonance, both of which were already in place. The paper deserves a serious referee, but with the SIDM claim and the v_R tuning stated more honestly.\n\nWhat it does well: the coupled Boltzmann treatment is careful and internally consistent; the new diagonal family in Fig. 3 is a real result, and the discussion of chemical potentials clarifies the dynamics. The anomaly-free charge assignments and the two-step thermalization via charged kaons/pions are sensible, and the dark photon constraints in Fig. 6 are handled competently. The lifetime estimate is at least honest about being 'up to O(1) factors.'\n\nThe soft spots, in order:\n\n1. v_R is a tuned input, and the stress-test arithmetic is right. For m_pi=20 MeV, v_R=85 km/s implies m_eta−2m_pi≈0.4 eV. The one-loop chiral correction is ~40 keV. The paper admits it is agnostic about the origin of v_R and cites possible UV constructions, but the benchmark is still a chosen point, not an outcome of the framework. Both the 3→2 freeze-out and the sharp scattering peak rest on this near-degeneracy. This should be front and center, not tucked into Sec. 2.2.\n\n2. The lifetime in Eq. (4.4) is order-of-magnitude only. The unknown O(p^6) coefficient and 'up to O(1) factors' matter for the benchmark: tau ~5e26 s sits right at the Voyager bound. The yellow exclusion region in Fig. 6 is correspondingly fuzzy.\n\n3. The abstract says the framework 'can be naturally explained' for SDSS J0946+1006, but the paper explicitly says no N-body/gravothermal simulation for Breit-Wigner cross sections exists. The 150 cm²/g at 60 km/s is a benchmark reproduction, not a demonstrated explanation of the lensing system. That sentence should be softened.\n\n4. The benchmark is tuned to hit Omega h²=0.12 and the SIDM cross sections. The relic calculation is a genuine solution of the Boltzmann equations, so it's not circular, but the parameter selection is fitting, not prediction.\n\nThe central mechanism holds up as a coherent EFT framework, but the fine-tuning burden is real. The paper is honest about many limitations in the text; it just doesn't always carry those caveats into the abstract/conclusions.\n\nWho it's for: dark-QCD and SIDM practitioners, and anyone working on MeV-scale thermal relics. I'd bring it to reading group and would cite it. Recommendation: send to peer review; accept after revision, with the astrophysical claim softened and the v_R/lifetime uncertainties given the prominence they deserve.","headline":"Serious and worth refereeing, but the benchmark carries a real tuned-input burden and the SIDM claim outruns the simulations.","tokens_in":27094,"tokens_out":3350,"would_cite":true,"duration_ms":34829,"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 nonzero theta angle lets a dark pion be a thermal dark matter relic, with the same η resonance producing the sharp, velocity-dependent self-interactions that halo data seem to require.","keywords":["dark matter","theta vacuum","QCD-like dark sector","dark pion","self-interacting dark matter","3-to-2 freeze-out","Breit-Wigner resonance","dark photon portal"],"falsifier":"A lattice calculation of the dark gauge theory at the three light quark masses used here could measure m_η−2m_π0 at θ=0; if the degeneracy fails by more than one part in about 10⁸, the resonant catalysis that sets both the relic density and the sharp velocity dependence disappears. Separately, a gravothermal N-body simulation using the paper's Breit-Wigner σ(v) could test whether the benchmark still satisfies the cluster bound σ/m≲0.5 cm²/g at roughly 2000 km/s while producing a lensing subhalo resembling J0946+1006.","tokens_in":25937,"feed_emoji":"⚛️","tokens_out":8839,"duration_ms":91937,"temperature":0.7,"pith_summary":"Dark pions are attractive dark matter candidates, but a simple pion-only strong sector struggles to produce the observed relic density and velocity-dependent self-interactions at the same time. This paper argues that a nonzero topological angle in a QCD-like dark sector solves both problems at once. The angle generates odd-number meson couplings, so the η meson can catalyze 3-to-2 number-changing processes that freeze out to the observed relic abundance, provided the η is nearly degenerate with two π0s. The same η mediates elastic pion scattering through a Breit-Wigner resonance, producing self-interactions that are large near a characteristic velocity around 85 km/s and small in galaxy clusters, exactly the pattern invoked to explain the lensing system J0946+1006. With a dark photon portal providing thermal contact with the Standard Model, the paper shows the scenario passes current constraints and maps out the viable parameter space.","feed_headline":"A dark theta angle turns pions into self-interacting dark matter","feed_subtitle":"The same η resonance sets the relic abundance and gives halos the sharp velocity-dependent collisions they need.","key_machinery":"The load-bearing object is the θ-dependent part of the chiral Lagrangian: for nonzero θ, an odd-parity trilinear ηππ vertex and a five-pion vertex appear, supplying number-changing interactions. The η exchange then provides a resonant s-channel for π0π0 scattering, described by the non-relativistic Breit-Wigner formula with resonance velocity v_R = 2√((m_η−2m_π0)/m_π0). Setting v_R at θ=0 to zero makes both the catalysed 3-to-2 freeze-out and the sharply peaked self-interaction cross-section work; the dark photon portal supplies thermalization, fixes the mediator width, and controls the DM lifetime.","core_discovery":"The paper claims that a nonzero θ angle in a QCD-like dark sector provides a self-contained thermal-relic framework: the same θ-induced interactions that set the pion abundance also control halo-scale self-scattering. The effective chiral Lagrangian acquires odd-meson vertices, and the η meson acts as a catalysing resonance, so that π0π0→η followed by ηπ0→π0π0 realizes resonant 3-to-2 freeze-out. After freeze-out, elastic π0π0→η→π0π0 scattering is a non-relativistic Breit-Wigner resonance with a peak velocity v_R determined by m_η−2m_π0. For a benchmark with m_π0=20 MeV and f_π=34 MeV, the paper obtains Ωh²≈0.12 together with σ/m ≲ 0.5 cm²/g at cluster velocities and ≳ 150 cm²/g near 60 km/s","pith_inferences":["If a UV realization or lattice computation independently fixes m_η≈2m_π0 at θ=0, then v_R becomes a derived number rather than an input, making the framework genuinely predictive from the dark quark masses alone.","The same θ-induced odd couplings should appear in SO and Sp dark gauge groups; if the analog of the η-resonance condition holds there, this mechanism could produce a family of SIDM models with distinct meson spectra and peak velocities.","A resonance velocity near 85 km/s implies a specific ordering of halo evolution: halos with internal velocities near v_R should undergo gravothermal collapse earlier than both slower and faster halos, creating a distinctive mass dependence in the abundance of core-collapsed subhalos.","The framework ties the relic abundance to the SIDM cross-section because both are controlled by the same ηππ coupling; a measurement of the peak velocity and height of σ/m at dwarf-galaxy scales would therefore also constrain the couplings that set the thermal relic density."],"forward_implications":["A single benchmark point with m_π0≈20 MeV, f_π≈34 MeV, and v_R≈85 km/s simultaneously reproduces the relic density and the self-interaction constraints from clusters and galaxies.","The framework predicts a Breit-Wigner peak in σ/m near roughly 100 km/s, so the strongest astrophysical signatures should appear in halos with characteristic velocities near that scale, including dwarf galaxies and lensing substructures.","Because the dark photon controls both thermalization and DM decay lifetime, indirect searches, including Voyager-type electron measurements and future 21-cm power-spectrum observations, can probe much of the allowed parameter space.","Kinetic decoupling after freeze-out and dark-photon mass corrections change the relic density by tens of percent, so quantitative predictions require the full Boltzmann treatment with the dark temperature evolving as 1/a² after decoupling.","A dedicated gravothermal simulation using the Breit-Wigner cross-section is needed before core-collapse predictions, including the interpretation of J0946+1006, can be confirmed or rejected."],"supporting_citations":[{"why":"Defines the benchmark model BM1 and the θ-induced number-changing dynamics that this paper extends with an explicit Standard Model portal.","marker":"[2]"},{"why":"Establishes the strongly interacting massive particle 3-to-2 thermal relic mechanism whose parameter space this scenario must reproduce.","marker":"[1]"},{"why":"Supplies the non-relativistic Breit-Wigner formula and the resonance-velocity parametrization used for the self-interaction predictions.","marker":"[50]"},{"why":"Introduces the dark photon portal used to thermalize the dark sector with the Standard Model.","marker":"[33]"},{"why":"Reports the strong-lensing system J0946+1006, the observational anchor for large σ/m at approximately 60 km/s.","marker":"[42]"},{"why":"Describes the catalysed-annihilation picture underlying resonant 3-to-2 freeze-out via an intermediate state such as η.","marker":"[54]"},{"why":"Provides the measured dark matter relic abundance used to fix Ωh²=0.12 in the parameter scan.","marker":"[55]"}],"fun_headline_variants":["Dark theta angle: pions as self-interacting thermal relic","Pion dark matter from a nonzero theta angle","Theta vacuum makes pions a thermal relic with velocity-dependent collisions","Same eta resonance sets relic abundance and halo self-interactions","Nonzero theta angle yields pions as sharp velocity-dependent DM"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The scenario rests on the dark η being almost exactly twice as massive as the dark π0 at θ=0, a resonance condition accurate to roughly one part in 10⁸ that is put in by hand and not derived from a UV theory; if that mass coincidence is absent, both the relic-density mechanism and the sharp velocity dependence vanish.","fun_headline_variants_meta":{"raw":{"variants":["Dark theta angle: pions as self-interacting thermal relic","Pion dark matter from a nonzero theta angle","Theta vacuum makes pions a thermal relic with velocity-dependent collisions","Same eta resonance sets relic abundance and halo self-interactions","Nonzero theta angle yields pions as sharp velocity-dependent DM"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000712,"raw_usage":{"total_tokens":3029,"prompt_tokens":724,"completion_tokens":2305,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":468,"completion_tokens_details":{"reasoning_tokens":2223}},"tokens_in":468,"tokens_out":2305,"duration_ms":16231,"temperature":1.0,"reasoning_tokens":2223,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T14:34:50.805810+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A lattice calculation of the dark gauge theory at the three light quark masses used here could measure m_η−2m_π0 at θ=0; if the degeneracy fails by more than one part in about 10⁸, the resonant catalysis that sets both the relic density and the sharp velocity dependence disappears. Separately, a gravothermal N-body simulation using the paper's Breit-Wigner σ(v) could test whether the benchmark still satisfies the cluster bound σ/m≲0.5 cm²/g at roughly 2000 km/s while producing a lensing subhalo resembling J0946+1006.","supporting_citations":[],"review_version":1}