{"id":"fdc2006e-3fb9-4235-970b-fb5b799b5332","arxiv_id":"2509.07102","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"NLO chiral perturbation theory for non-degenerate two-flavor real and pseudoreal gauge theories is derived, fitted to Sp(4) lattice data, and applied to pion dark matter self-interactions.","lead":"This paper derives next-to-leading-order formulas for pion properties in QCD-like theories with two unequal quark masses, for gauge groups with real or pseudoreal representations. Fitting the leftover constants to lattice data for an Sp(4) theory, it shows the corrections alter dark-matter self-interaction predictions by tens to hundreds of percent and eliminate the naive SIMP dark-matter mass window.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The SIMP exclusion rests on the NLO cross-section at M_pi/F_pi ~ 6.5-9, where the chiral expansion parameter is ~0.27-0.5 and NNLO-sized effects (decay-constant splitting) are already visible in the fitted data; with no truncation-error estimate, the 'window entirely closed' claim in Fig.","rationale":"The paper is a serious, careful ChPT study: the NLO formulas for non-degenerate masses are checked against known limits [25,26], the Mathematica files are released (Ref. [41]), the fit marginalizes over latent variables, and L̃1–L̃3 are consistent across β. I found no fatal error; the verdict CONDITIONAL stands. The single most load-bearing concern is convergence control of the quantitative headline. The exclusion ('At NLO, the viable mass window is entirely closed') is a prediction of Eq. (6.5) in a regime where the chiral expansion parameter is 0.27–0.5, the NLO term is comparable to or larger than LO, and the fitted data already display an NNLO effect — the decay-constant splitting — that the NLO formulas cannot describe and that is removed by averaging (Sec. 5.1). The post hoc r = 5.5 cut and the β-dependence of L̃4 are further indications that the regime is marginal, even though L̃4 does not enter the degenerate cross-section. Because L̃1, the dominant combination in (6.5), is fixed by the mass and decay-constant fits at exactly these high-mass points, the cross-section is not an independent test of the EFT in this region but an interpolation whose NNLO uncertainty is unquantified; the authors themselves concede (footnote 19) that fitting NNLO formulas would shift the NLO LECs. My proposed check — refitting on a reduced-r window or adding an NNLO-sized theory error — settles whether the exclusion survives a conservative truncation-error estimate. I agree with the reader that the load-bearing premise is the quantitative validity of the NLO EFT at the fitted lattice points; my refinement concerns the mechanism (L̃4 as diagnostic, L̃1 and the F-averaging as the direct channel). Since the concern is the one the reader identified and no new failure mode emerged, the verdict remains CONDITIONAL.","tokens_in":34783,"tokens_out":21691,"duration_ms":171586,"concrete_test":"Refit the Sec. 5.2 β=7.2 data with the points closest to the vector-meson regime removed — restrict the fit to r ≤ 3.5 instead of the paper's r ≤ 5.5, and drop the highest-M_pi/F_pi scattering-length point — then recompute the Fig. 11 NLO/LO crossing and the Fig. 13 exclusion. If the crossing point moves by more than ~0.5 in M_pi/F_pi, or if a viable window reopens below 4π, the 'entirely closed' claim is an artifact of fitting NLO formulas where they are not converging. As a complementary analytic check, add an NNLO-sized theory error of order (M_pi/4πF_pi)^2 times the NLO term to the fit likelihood; if the window reopens under this error, the central claim requires NNLO control that the paper does not provide.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's sharpest claim is that at NLO the SIMP mass window is entirely closed (Fig. 13). The exclusion is driven by the NLO 2→2 cross-section, Eq. (6.5) and Fig. 11, evaluated at M_pi/F_pi ≈ 6.5 and beyond, using LECs fitted to lattice data at M_pi/F_pi ≈ 5–7. For the claim to be controlled, the NLO cross-section near the crossing must be accurate at the ~10–30% level; three facts make this assumption fragile. (1) Convergence: the expansion parameter (M_pi/4πF_pi)^2 ≈ 0.27–0.5 in the quoted region, and the NLO correction is already 100% of LO at M_pi/F_pi ≈ 4.0 (Sec. 6.1), so an NNLO term of expected size shifts the crossing point by order one in M_pi/F_pi. (2) The fitted data exhibit NNLO-sized physics that the NLO formulas cannot describe: the lattice decay constants split while Eq. (3.17) does not, so the authors average F_pi,1 and F_pi,3 (Sec. 5.1), and the fit requires a post hoc removal of r > 5.5 points because 'in this regime our EFT is no longer applicable' (Sec. 5.1). Because the LECs are extracted with the same NLO formulas at these same high masses, any NNLO or discretization contamination is absorbed into the LECs and re-emitted in the cross-section; the dominant combination in Eq. (6.5) is L̃1 (coefficient 456 vs −112 and 105), pinned by the mass and decay-constant fits (5.7)–(5.9) in which the F_pi splitting is averaged away. (3) The β-dependence of L̃4 (Fig. 5) signals that the fitted data are not continuum-controlled; L̃4 itself drops out of the degenerate cross-section (6.5), so its role is diagnostic. The authors label the SIMP application an 'illustration' because vector mesons alter the freeze-out curve (App. C), yet the Fig. 13 caption states the exclusion without that caveat, and Sec. 6.1 calls the fit-region prediction 'fully trustable'. No algebraic inconsistency is apparent; the issue is that the quantitative claim 'entirely closed' carries no convergence estimate.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper constructs next-to-leading-order (NLO) chiral perturbation theory for QCD-like gauge theories with N_F=2 Dirac fermions in real (SU(4)/SO(4)) and pseudoreal (SU(4)/Sp(4)) representations with non-degenerate quark masses. It derives NLO formulas for pion masses, quark condensates, decay constants, and 2→2 scattering amplitudes, and shows that the degenerate-mass limits reproduce earlier results while correcting two coefficient errors in the literature. The authors then fit the scale-independent combinations L̃1,...,L̃4 of NLO low-energy constants to Sp(4) lattice data from Refs. [29,30] using a Bayesian MCMC procedure with latent variables, and use the fitted LECs to compute the non-relativistic 2→2 pion self-interaction cross-section. Applied to SIMP dark matter, the NLO cross-section exceeds the Bullet Cluster bound at Mπ/Fπ≈6.5 instead of ≈9 at LO, and the authors conclude that the previously viable SIMP mass window is entirely closed at NLO.","tokens_in":35280,"tokens_out":8886,"duration_ms":76072,"significance":"If the analysis is correct, this is a valuable step: it provides NLO formulas for non-degenerate masses in both real and pseudoreal two-flavour theories, a first LEC determination for Sp(4) with N_F=2 from lattice spectroscopy and scattering, and a concrete demonstration that NLO corrections shift pion self-interaction predictions in a way that matters for SIMP models. The paper ships machine-readable Mathematica code and a data release, reproduces known degenerate limits, and is careful to separate statistical fit uncertainties from the missing continuum extrapolation. The LEC extraction is not circular: the lattice data are external measurements, and the NLO formulas are independent of the dark-matter cross-section claim. The main weakness is not the derivation itself but the quantitative control of the EFT at the large values of Mπ/Fπ used in the fit and in the application; unless a truncation-error estimate is added, the 'window entirely closed' conclusion remains conditional.","major_comments":[{"comment":"The LEC fit is performed at Mπ/Fπ ≈5–7, where the chiral expansion parameter (Mπ/4πFπ)² is 0.16–0.31, and Sec. 6.1 reports that the NLO correction is already 100% of the LO cross-section at Mπ/Fπ≈4.0. The NLO formulas (5.7)–(5.10) are therefore not a small correction at the fitted points, and no NNLO truncation estimate is given. Because the fit uses these same NLO formulas, any NNLO or discretization contamination in the lattice data is absorbed into L̃1–L̃4 and is re-emitted in the cross-section (Eq. (6.5)), whose dominant LEC term is 456L̃1. Please add a quantitative truncation-error estimate, or explicitly frame the extraction and the subsequent SIMP exclusion as conditional on NLO convergence at Mπ/Fπ≳5.","section":"Secs. 5.1, 5.2; Fig. 11"},{"comment":"The lattice data of Ref. [29] show a splitting between Fπ,1 and Fπ,3 that is absent at NLO in Eq. (3.17). The authors construct a single averaged dataset with a systematic error equal to half the channel difference (Sec. 5.1, Fig. 9). This prescription does not propagate the fact that the NLO model cannot describe the splitting, and it neglects correlations between the two channels; because L̃1 carries the dominant coefficient (456) in Eq. (6.5) and is pinned by the mass and decay-constant fits (5.7)–(5.9), the effect on the predicted cross-section needs to be assessed explicitly.","section":"Sec. 5.1, Eqs. (3.17), (5.9)"},{"comment":"The fitted L̃4 shows strong β-dependence, which the authors attribute to discretization artifacts; the continuum limit is not performed. The analysis adopts the β=7.2 values as central and does not add a lattice-spacing systematic to L̃1–L̃3. Although L̃1–L̃3 are reported as mutually consistent across β, a quantitative statement of the residual β-sensitivity of the cross-section combination (456L̃1−112L̃2+105L̃3)/(245760π³) in Eq. (6.5) is necessary to support the precision of the Bullet-Cluster crossing. Please provide this propagation or a conservative systematic band.","section":"Sec. 5.2, Fig. 5"},{"comment":"The caption 'At NLO, the viable mass window is entirely closed' is stronger than the caveats in the text allow. Appendix C states that the freeze-out curve in Fig. 13 is an illustration and that vector mesons with Mρ<2Mπ for Mπ/Fπ≳4 will alter the number-changing processes that fix Mπ/Fπ for a given Mπ. The Bullet-Cluster exclusion itself is more robust, but the closing of the window combines this exclusion with the freeze-out curve; the latter is not computed with the same NLO precision. Please soften the conclusion or present the window closing as conditional on the NLO EFT and on neglect of vector-meson effects in freeze-out.","section":"Sec. 6.1, Fig. 13, Appendix C"}],"minor_comments":[{"comment":"The terms 'M^2_{4,6}' (in the L^r_5 term) and 'M^4_{4,1}' (in the log term) appear to be typos for M^2_{π,6} and M^4_{π,1}; please correct.","section":"Eq. (4.18)"},{"comment":"The statement that the fit is robust to moderate variations of the prior width is not quantified; a short sentence or table giving the range of prior widths tested would be helpful.","section":"Sec. 5.2"},{"comment":"Since the β-dependence is a central systematic, showing all three β datasets in the same plot, or a small table of M^2_{π,1}/M^2_{π,3} at each r, would make the discretization effect easier to judge.","section":"Fig. 6"},{"comment":"The factor-of-two difference from Ref. [25] in the NLO condensate would benefit from a one-sentence explanation of the normalization convention, as it is otherwise easy to misread as an error.","section":"Footnote 8"},{"comment":"The statement that the 3σ band would shrink to the line width at 1σ is useful; please also state the central value and 1σ range of the crossing point Mπ/Fπ where σ2→2/Mπ equals the Bullet Cluster bound.","section":"Sec. 6.1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript's main quantitative claim in Fig. 13 goes somewhat beyond what the text's own caveats support; in revision, the authors should make the 'illustrative' status explicit in the figure caption. The self-citation via Ref. [29] is not a concern because the lattice data are external measurements."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things about this one. First, the genuinely new deliverable is the set of NLO chiral formulas for non-degenerate quark masses in SU(4)/Sp(4) and SU(4)/SO(4), with a first extraction of LEC combinations from Sp(4) lattice data. Second, the headline dark-matter claim—that NLO closes the SIMP window entirely—is not backed by a convergence estimate, so treat it as conditional.\n\nThe algebra is the paper's strength. Sections 3 and 4 systematically extend the degenerate results of Bijnens–Lu to r ≠ 1, they reproduce the degenerate limits, and they flag and fix two sign/transposition errors in earlier literature. That is real and useful. The authors also ship the formulas in a Mathematica repository, which lets anyone cross-check the long expressions. The LEC fit is a reasonable first attempt: they marginalize over latent variables, check prior sensitivity, and show beta-dependence of one combination. The LECs come from external lattice data, so the central calculation is not circular. Honest limitations are stated in Section 7.\n\nThe soft spots are concentrated where the paper goes from fitting to predicting. The lattice data are finite-spacing, not continuum-extrapolated; a post hoc cut at r = 5.5 is applied because the EFT is no longer applicable; the NLO decay constant does not split, so the lattice F_pi values are averaged with a half-difference systematic. The fit region is M_pi/F_pi ~ 5–7, with expansion parameter ~0.27, and the NLO correction is already 100% of LO at M_pi/F_pi ~ 4. Against that, the claim that the Bullet Cluster bound is crossed at 6.5 rather than 9, and the Fig. 13 caption \"entirely closed,\" carry no NNLO estimate. Given that NNLO-size physics (decay-constant splitting) is already visible in the fit data, the LECs themselves may absorb NNLO and discretization contamination, then re-emit it into the cross section. L̃1 dominates the coefficient in Eq. (6.5) and is pinned by fits that average away the F_pi splitting, so this is not a marginal worry. The authors do label the SIMP application an illustration because vector mesons change the freeze-out curve, but the abstract and Fig. 13 caption are stronger than that caveat.\n\nMy read: the paper deserves a serious referee. The NLO formulas are a solid reference result with reproducible code; the fit is useful even if the systematics need tightening; the SIMP exclusion should be framed with an explicit truncation-error band, not as a closure. I would send it to review, and I would cite the formulas.","headline":"Non-degenerate NLO chiral formulas and a first Sp(4) LEC fit are real contributions; the SIMP 'window entirely closed' claim needs a truncation-error band before it can be taken at face value.","tokens_in":35969,"tokens_out":2636,"would_cite":true,"duration_ms":22640,"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":"At next-to-leading order, the viable SIMP pion dark-matter mass window is entirely closed.","keywords":["chiral perturbation theory","next-to-leading order","low-energy constants","Sp(4) gauge theory","pion dark matter","strongly interacting massive particles","self-interacting dark matter","lattice QCD-like theories"],"falsifier":"Measure the threshold $2\\to2$ pion cross section directly on the lattice for the mass-degenerate $Sp(4)$ theory at $M_\\pi/F_\\pi\\approx 6.5$, using the scattering length in the 14-dimensional channel at several volumes with a continuum extrapolation; if the lattice cross section lies below the Bullet Cluster bound while the NLO prediction lies above it, the claimed closure of the SIMP window is wrong.","tokens_in":34517,"feed_emoji":"🌌","tokens_out":15413,"duration_ms":121314,"temperature":0.7,"pith_summary":"The paper establishes next-to-leading-order chiral perturbation theory expressions for QCD-like gauge theories with two non-degenerate fermion flavors in real and pseudoreal representations, covering pion masses, condensates, decay constants, and $2\\to2$ scattering amplitudes. It then fits the NLO low-energy constants to existing lattice spectroscopy and scattering data for the $Sp(4)$ gauge theory with two fundamental flavors, and applies the fitted result to strongly interacting massive particle (SIMP) dark matter. The central conclusion is that NLO corrections substantially raise the non-relativistic pion self-scattering cross section: the Bullet Cluster constraint is crossed at $M_\\pi/F_\\pi\\approx 6.5$ instead of $\\approx 9$ at leading order, and in the degenerate $Sp(4)$ SIMP benchmark the viable mass window is entirely closed. If this is right, leading-order-only treatments of pion dark matter self-interactions are quantitatively inadequate precisely where SIMP models would have to live.","feed_headline":"Corrections beyond leading order close the pion dark matter window","feed_subtitle":"Fitted corrections push the self-interaction bound to pion mass ratio 6.5, ruling out the benchmark model.","key_machinery":"The load-bearing machinery is the $\\mathcal{L}_2+\\mathcal{L}_4$ chiral Lagrangian for the cosets $SU(4)/Sp(4)$ and $SU(4)/SO(4)$, with the $\\mathcal{L}_4$ term containing the low-energy constants $L_0,\\ldots,L_8,H_2$. The calculation combines one-loop diagrams from $\\mathcal{L}_2$ with tree-level $\\mathcal{L}_4$ vertices to produce NLO expressions for masses, condensates, decay constants, and the $2\\to2$ scattering amplitudes; these are reduced to the non-relativistic threshold cross section $\\sigma_{2\\to2}$ and to the scattering length $a_0^{\\mathrm{MS}}$ in the 14-dimensional channel. The fit uses scale-independent combinations $\\tilde L_1,\\ldots,\\tilde L_4$ of the renormalized LECs so that chiral logarithms drop out of the fitted observables. The same LEC combination controls both the fitted scattering length and the dark-matter cross section, which is why the lattice fit can be converted directly into a prediction for self-interactions.","core_discovery":"On the paper's own terms, the claim is that NLO chiral EFT with low-energy constants taken from lattice data gives the quantitative description of dark-pion self-interactions in $SU(4)/Sp(4)$-type theories, and that this description rules out the simplest $Sp(4)$ SIMP model. The paper derives NLO formulas for pion masses, condensates, decay constants, and scattering amplitudes with two non-degenerate quark masses, shows that pion mass splitting first appears at NLO, fits the four accessible combinations of NLO LECs ($\\tilde L_1,\\ldots,\\tilde L_4$) to the $Sp(4)$ lattice data, and evaluates the threshold cross section. The fitted NLO cross section exceeds the Bullet Cluster bound for $M_\\pi/F_\\pi\\gtrsim 6.5$, compared with $\\approx 9$ at LO; for the benchmark $M_\\pi=0.2$ GeV the whole viable SIMP window disappears.","pith_inferences":["If the exclusion is correct, realizing SIMP dark matter with this gauge group requires one of the loopholes the paper lists: working where vector mesons invalidate the pion-only EFT, adding a portal that destabilizes the singlet, or moving to $M_\\pi/F_\\pi$ below the lattice-fit region where the LECs are untested.","The same fitted LECs can be used to correct early-Universe phenomenology beyond the WZW term, since the NLO $2\\to2$ amplitudes enter thermal averaging near freeze-out and are likely to shift the relic-density curve shown in the SIMP parameter space.","A direct lattice determination of the threshold $2\\to2$ cross section in the $SU(4)/Sp(4)$ theory at $M_\\pi/F_\\pi\\approx 6.5$ would test the NLO prediction without any dark-matter assumptions.","The fitting pipeline transfers to other real and pseudoreal theories: once lattice data with split masses appear for $SU(4)/SO(4)$, the same procedure yields LECs and a NLO self-interaction prediction for composite-Higgs-related pion phenomenology."],"forward_implications":["A leading-order-only analysis of SIMP dark matter in this theory is not conservative: NLO moves the self-interaction bound by about $\\Delta(M_\\pi/F_\\pi)\\approx 2.5$ and closes the window.","In the degenerate $Sp(4)$ case, the fitted LECs turn the earlier random-LEC estimate of the NLO cross section into a definite prediction with quantified $3\\sigma$ uncertainty.","For non-degenerate masses the singlet pion is lighter than the four-plet; if the singlet is stable, singlet-only dark matter has a much smaller self-scattering cross section, broadening the viable parameter space.","For $SU(4)/SO(4)$ theories the NLO mass hierarchy typically keeps the $\\pi_6,\\pi_7$ doublet lightest, and its self-interaction is six times larger than the single-pion estimate used previously.","The NLO formulas provide the missing interpolation between the chiral limit used for composite Higgs models and the finite-mass regime where lattice simulations are performed."],"supporting_citations":[{"why":"Supplies the lattice spectroscopy data for Sp(4) with non-degenerate quark masses that anchor the pion-mass and decay-constant part of the LEC fit.","marker":"[29]"},{"why":"Supplies the lattice scattering length in the 14-dimensional channel that anchors the scattering part of the fit.","marker":"[30]"},{"why":"Provides the NLO chiral Lagrangian and renormalization coefficients for QCD-like theories that the paper extends to non-degenerate masses.","marker":"[25]"},{"why":"Provides the meson-scattering amplitudes and loop function that the paper adapts, corrects, and reduces to threshold cross sections.","marker":"[26]"},{"why":"Established that NLO corrections can reshape SIMP viability and supplied the random-LEC treatment that the fitted LECs replace.","marker":"[24]"},{"why":"Defines the SIMP scenario and gives the LO self-scattering cross section and WZW 3-to-2 rate that set the benchmark the paper re-evaluates.","marker":"[18]"},{"why":"Supplies the scale-independent LEC redefinition used to eliminate chiral logarithms from the fitted observables.","marker":"[27]"},{"why":"Shows from lattice data that vector mesons become light at large $M_\\pi/F_\\pi$, motivating the $r=5.5$ cutoff and the vector-meson caveat on freeze-out.","marker":"[43]"}],"fun_headline_variants":["NLO corrections rule out benchmark pion dark matter model","Fitted NLO corrections close pion dark matter window","Pion dark matter: NLO fit kills benchmark model","NLO corrections close SIMP dark matter window","Lattice-constrained NLO rules out pion dark matter"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole exclusion rests on believing that the next-to-leading-order chiral expansion still works at the fairly heavy pions of the lattice data, even though the fit needed a cut at $r=5.5$ because the vector meson was getting close to the heavier pion.","fun_headline_variants_meta":{"raw":{"variants":["NLO corrections rule out benchmark pion dark matter model","Fitted NLO corrections close pion dark matter window","Pion dark matter: NLO fit kills benchmark model","NLO corrections close SIMP dark matter window","Lattice-constrained NLO rules out pion dark matter"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000612,"raw_usage":{"total_tokens":2842,"prompt_tokens":936,"completion_tokens":1906,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":552,"completion_tokens_details":{"reasoning_tokens":1829}},"tokens_in":552,"tokens_out":1906,"duration_ms":14340,"temperature":1.0,"reasoning_tokens":1829,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T16:14:38.126293+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the threshold $2\\to2$ pion cross section directly on the lattice for the mass-degenerate $Sp(4)$ theory at $M_\\pi/F_\\pi\\approx 6.5$, using the scattering length in the 14-dimensional channel at several volumes with a continuum extrapolation; if the lattice cross section lies below the Bullet Cluster bound while the NLO prediction lies above it, the claimed closure of the SIMP window is wrong.","supporting_citations":[{"cited_title":"Technicolor and other QCD-like theories at next-to-next-to-leading order","cited_arxiv_id":"0910.5424","evidence_quote":"Provides the NLO chiral Lagrangian and renormalization coefficients for QCD-like theories that the paper extends to non-degenerate masses."}],"review_version":1}