{"id":"7030f8f6-bf7d-4d2f-a041-4b754a912fa1","arxiv_id":"1909.00005","paper_version":2,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"Existing PIENU and PIBETA rare pion decay data constrain ALP-pion mixing to sin^2 theta below about 1e-5 for ALP masses from 10 MeV to the pion mass.","lead":"Two old measurements of rare pion decays, the positron spectrum from pi+ to e nu and the photon opening angle from pi+ to pi0 e nu, can be repurposed to constrain how much a light axion-like particle mixes with the neutral pion. The new bounds exclude mixing angles down to around 1e-5 for ALP masses between about 10 and 140 MeV, complementing existing limits from kaon decays.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The mapping from pi+ -> a e nu data to sin^2 theta assumes the pi0-mixing amplitude (Eq. 6) dominates; the residual rho-exchange term in Eq. (5) is not quantified for small sin theta.","rationale":"Both the reader and I identify the same soft spot: the hadronic matrix element relation in Eq. (6) is the step that turns a rare-decay rate into a constraint on ALP-pion mixing. The paper's one-sentence suppression argument is plausible in vector-meson-dominance counting and is likely correct at the O(10%) level for m_a near m_pi, but it is not a derivation, and the danger is concentrated at small sin theta where the mixing amplitude vanishes. This is load-bearing because the headline reach sin^2 theta ~ 1e-5 sits in that regime. No code or independent check of Eq. (6) is provided, so the check I propose would settle the issue. If the check passes, the ACCEPT verdict stands; if it fails, the contours shift or lose their interpretation as bounds on sin theta. I therefore recommend CONDITIONAL rather than unconditional acceptance: the authors should either supply the EFT check or explicitly state that the bounds apply only when direct ALP-hadron operators are subdominant to the mixing contribution. This does not call into question the PIENU/PIBETA statistical treatment, which is validated against sterile-neutrino bounds and clearly labeled as approximate.","tokens_in":11033,"tokens_out":19483,"duration_ms":176812,"concrete_test":"Evaluate the second term of Eq. (5) in the chiral Lagrangian with vector mesons for a representative UV completion, e.g., the derivative ALP-quark coupling (d_mu a/f_a) bar q gamma^mu gamma5 q matched as in Ref. [18]. Compute the virtual-rho exchange contribution to pi+ -> a e nu using the rho decay constant and the ALP-rho-pi vertex, and compare its magnitude with sin theta times the pibeta amplitude from Eq. (6) across 10 MeV < m_a < 140 MeV and sin^2 theta at the Fig. 3 exclusion boundary. If the ratio |A_rho|/|sin theta A_pi0| exceeds about 10% anywhere in the excluded region, the contours in Fig. 3 overstate the sin theta reach and must be revised or presented as model-dependent.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central exclusions in Fig. 3 rely on Eq. (9), which converts a bound on Br(pi+ -> a e nu) into a bound on sin^2 theta. That conversion assumes the hadronic amplitude is dominated by the off-shell pi0 insertion in the first term of Eq. (5), so that A_mu ~ sin theta <pi*0| bar d gamma_mu u |pi+> (Eq. 6). The discarded second term of Eq. (5), which receives its leading chirally-unsuppressed contribution from virtual rho exchange, is dismissed by the parametric statement that it is 'suppressed by m_a^2/m_rho^2'. No bound on the coefficient of that term is given. In a generic UV completion the same ALP-quark operators that produce the mixing also generate the rho-pole contribution. Because the mixing term is proportional to sin theta, which can be arbitrarily small, the ratio (direct term)/(mixing term) grows as sin theta -> 0; for sin^2 theta ~ 1e-5 the direct term can dominate unless its coefficient is numerically bounded. Thus the Fig. 3 contours are not yet demonstrated to be model-independent bounds on sin theta. The paper itself writes 'Other contributions may involve mixing with other hadrons, or other UV operators', but the exclusion curves are presented without this caveat.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proposes new constraints on the mixing angle between a light axion-like particle (ALP) and the neutral pion, using existing rare pion decay data. The authors argue that UV completions of ALPs lighter than m_pi generically induce ALP-pion mixing, and that the semileptonic decay pi+ -> a e nu can be related, via mixing with an off-shell pi0, to the Standard Model decay pi+ -> pi0 e nu. Using the PIENU positron-energy residuals from a previous sterile-neutrino analysis and the PIBETA diphoton opening-angle spectrum, they derive 95% CL exclusions in the sin^2 theta versus m_a plane, reaching sin^2 theta ~ 10^-5 for masses between about 25 and 130 MeV. They also explore a 'pure mixing' scenario in which the ALP-photon coupling is set by sin theta times the pion-photon coupling, and combine their bounds with existing beam-dump and collider limits. The analysis assumes no tree-level ALP-lepton couplings, a dominant diphoton decay mode, and either prompt or invisible ALP lifetimes, with the pure mixing scenario interpolating between the two regimes.","tokens_in":11266,"tokens_out":9722,"duration_ms":93648,"significance":"If the derived bounds are correct, they provide a genuinely new probe of ALP-pion mixing that is complementary to kaon-based limits: the process pi+ -> a e nu arises from a tree-level charged-current amplitude and is not subject to the ALP-top-quark short-distance contributions that complicate K+ -> pi+ a estimates. The central formula, Eq. (9), is cleanly derived and cross-checked against the known pi+ -> pi0 e nu rate to O[(1-r)^8], and the PIENU analysis is validated by reproducing the published pi+ -> e N constraints. The paper is also transparent about many of its approximations, including the uncorrelated-bin assumption and the approximate nature of the PIBETA bound. However, the central claim of model-independent constraints on sin theta is weakened by an unquantified hadronic assumption: the direct ALP-hadron contribution in Eq. (5) is dismissed parametrically, but no bound on its coefficient is given, so the mapping from the branching-ratio limit to sin^2 theta is not established for generic UV completions. Because of this, the paper should be revised to either quantify that contribution or carefully qualify the central claims.","major_comments":[{"comment":"The central step of the paper is the approximation in Eq. (6) that the pi+ -> a e nu amplitude is saturated by the off-shell pi0 insertion from the first term of Eq. (5). The second term, whose leading chirally-unsuppressed contribution comes from virtual rho exchange, is dismissed because it is 'suppressed by m_a^2/m_rho^2'. This is a parametric suppression, but no bound on the coefficient of this term is given. Since sin theta is treated as a free phenomenological parameter, the ratio of the direct rho-exchange contribution to the mixing contribution scales as (1/sin theta) times hadronic matrix elements; for sin^2 theta ~ 10^-5, sin theta ~ 3 x 10^-3, and for m_a near m_pi one has m_a^2/m_rho^2 ~ 0.02, so a coefficient of order a few makes the direct term comparable to the mixing term. The exclusion contours in Fig. 3 and the abstract's claim of 'strong and novel bounds on the ALP-pion mixing angle' are therefore not established as model-independent bounds on sin theta. The authors should either derive a quantitative estimate of the direct term in a concrete EFT (e.g., chiral perturbation theory or vector-meson dominance) or explicitly qualify all central claims and figures as valid only when the direct hadronic contribution is negligible, with a discussion of how generic UV completions would alter the interpretation.","section":"PIENU residuals bound"},{"comment":"The PIENU exclusion limits in the left panel of Fig. 3 rely on treating all residual bins as uncorrelated. The only justification given is that this assumption reproduces the published pi+ -> e N bounds from Ref. [10]. Energy-spectrum residuals from a global fit are typically correlated, and the impact of possible bin correlations on the 95% CL contours is not quantified. The authors should provide a robustness test, for example a prescription with a simple positive-correlation model, or explicitly state that the derived limits are an approximation whose uncertainty from bin correlations is not assessed.","section":"PIENU residuals bound"},{"comment":"The PIBETA exclusion in Fig. 3 is derived from an order-of-magnitude requirement that the integrated ALP contribution in the 160-180 degree range does not exceed the quoted 0.6% branching-ratio uncertainty. This is not a statistical bound based on the full shape of the theta_gamma_gamma distribution, and the approximate treatment is described only in the text. Since the PIBETA bound is presented as one of the two new constraints in the abstract and in Fig. 3, its approximate character should be carried through to the abstract and the figure caption, or the bound should be replaced by a more rigorous likelihood-based estimate using the published spectrum.","section":"PIBETA diphoton bound"}],"minor_comments":[{"comment":"The statement that the theoretical approximations introduce no more than O(10%) uncertainties refers to the form-factor approximation f_+(q^2) ~ 1; it does not cover the uncertainty from the neglected second term in Eq. (5). Please clarify the scope of that uncertainty estimate.","section":"Conclusions"},{"comment":"Several references are incomplete: Refs. [28], [31], and [36] lack author lists and titles. Please complete them.","section":"References"},{"comment":"The y-axis label 'epsilon_a 10^5/Gamma[pi e2] x dGamma/dE_cal' is difficult to parse; consider rewriting it as, for example, '10^5/Gamma(pi->e nu) x dGamma/dE_cal (MeV^-1)'.","section":"Fig. 2"},{"comment":"The phrase 'direct sin theta constraints' in the text overstates the result given the caveats in the following section; a more qualified wording such as 'constraints under the stated assumptions' would be more accurate.","section":"ALP-pion mixing"},{"comment":"The sentence 'The irreducible background from Br[pi+ -> pi0 e nu] is much smaller than the experimental precision' could be expanded to state the numerical value of the pi beta branching ratio (about 1.04 x 10^-8) here rather than only in the PIBETA section, for reader convenience.","section":"Introduction"}],"recommendation":"major_revision","confidential_remarks":"This is a well-written and useful phenomenological paper, and the central idea of using rare pion decays to constrain ALP-pion mixing is sound in principle. My main reservation is that the central claim of model-independent bounds on sin theta is not fully supported because the dominance of the pi0-mixing amplitude over direct hadronic contributions is assumed without a quantitative estimate. The bin-correlation and PIBETA approximations also deserve more prominent caveats. These issues appear addressable within the manuscript's scope, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Punchline: this is a solid phenomenological paper. It shows that existing PIENU and PIBETA rare pion decay data can be reinterpreted as constraints on ALP-pion mixing, reaching sin^2 theta around 1e-5 in parts of the parameter space. The derivation of Eq. (9) is clean, and the validation against the published pi+ -> e N bounds is a good sanity check.\n\nWhat is new is not the mixing formalism (that goes back to Krauss-Wise and Bardeen et al.) nor the idea of using PIENU residuals, but the actual exclusion regions in the sin^2 theta vs m_a plane. The bounds are complementary to kaon-based limits because they come from tree-level charged-current processes, so they are less dependent on assumptions about the ALP-top coupling. The paper is honest about its approximations: the uncorrelated PIENU bins, the approximate PIBETA threshold, and the f+ ~ 1 form factor are all flagged.\n\nThe main soft spot is the hadronic matrix element. Eq. (5) has two terms; the paper drops the second, arguing that the leading virtual-rho contribution is suppressed by m_a^2/m_rho^2. That parametric suppression is not enough on its own: the first term is proportional to sin theta, which can be as small as 1e-5 in the excluded region, and the second term has no such factor. Without a bound on the coefficient of the rho-pole term, the contours in Fig. 3 are not yet demonstrated to be model-independent bounds on sin theta. The authors do include a sentence noting other hadronic mixing or UV operators could contribute, but the abstract and conclusions present the constraints more definitively. The exclusions should be labeled as applying to scenarios where direct ALP-hadron couplings are subdominant.\n\nThe other limitations are minor. No code or data files are shipped, but the analysis is simple enough to reproduce. The PIBETA bound is approximate by their own admission. The citation pattern looks appropriate.\n\nWho is this for? Anyone doing ALP phenomenology in the MeV-GeV range. It gives a new way to leverage existing pion data and will be a useful reference for future rare-pion analyses. It deserves peer review. If I were the editor, I would send it out; my own verdict would be close to accept, with revisions to make the direct-hadronic-coupling caveat more prominent.","headline":"Clean, transparent analysis turning rare pion decay data into new ALP-pion mixing exclusions; main caveat is the unquantified direct hadronic contribution.","tokens_in":11875,"tokens_out":11791,"would_cite":true,"duration_ms":111147,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["13.20.Cz","14.80.Va"],"model":"deepseek-v4-flash","headline":"Existing pion-decay data from PIENU and PIBETA exclude axion-pion mixing down to $\\sin^2\\vartheta \\sim 10^{-5}$ at 95% confidence for light axion-like particles.","keywords":["axion-like particles","ALP-pion mixing","rare pion decays","PIENU","PIBETA","leptonic pion decay","pion beta decay","mixing angle constraints"],"falsifier":"A dedicated stopped-pion measurement of $\\pi^+\\to (a\\to\\gamma\\gamma)e\\nu$ reconstructing the diphoton opening angle for $10\\text{ MeV}\\lesssim m_a\\lesssim m_\\pi$ would settle the central claim: the predicted opening-angle distribution beyond the $\\pi^0$ edge, Eq. (11), should appear at the rate set by Eq. (9) for the claimed $\\sin^2\\vartheta$, and its absence would falsify the off-shell-$\\pi^0$ dominance assumption.","tokens_in":10755,"feed_emoji":"⚛️","tokens_out":9937,"duration_ms":87788,"temperature":0.7,"pith_summary":"Ultraviolet completions that produce axion-like particles (ALPs) lighter than the neutral pion generically induce a small admixture of the pion in the ALP state, parametrized by a mixing angle $\\vartheta$. The paper shows that this mixing can be constrained directly with data that already exist: the positron energy spectrum of $\\pi^+\\to e\\nu$ measured by PIENU, and the diphoton opening-angle spectrum of $\\pi^+\\to\\pi^0 e\\nu$ measured by PIBETA. In the ALP mass range $10\\text{ MeV}\\lesssim m_a<m_\\pi$, the derived 95% CL exclusions reach $\\sin^2\\vartheta$ of order $10^{-5}$, with PIBETA providing the strongest limit for prompt ALPs above roughly 100 MeV. Because these processes are tree-level charged-current decays, the bounds are largely independent of kaon-based ALP limits and apply to any UV model that generates ALP-pion mixing. If correct, they close a substantial part of the $\\sin^2\\vartheta$--$m_a$ plane.","feed_headline":"Pion decay data squeeze axion-pion mixing to 10^-5","feed_subtitle":"PIENU and PIBETA spectra exclude light axions mixed with pions across a wide mass range.","key_machinery":"The carrying object is the ALP-pion mixing angle $\\sin\\vartheta$, defined by the physical ALP state containing an admixture of the neutral pion. Its work is to connect an experimentally accessible process, $\\pi^+\\to a e\\nu$, to existing pion-decay measurements: the amplitude is assumed to be dominated by the off-shell $\\pi^0$ insertion $\\langle a|\\pi^{*0}\\rangle\\langle\\pi^{*0}|\\bar d\\gamma^\\mu u|\\pi^+\\rangle = \\sin\\vartheta\\langle\\pi^{*0}|\\bar d\\gamma^\\mu u|\\pi^+\\rangle$, and the $\\pi^+\\to\\pi^0 e\\nu$ rate is fixed by form factors normalized by the Ademollo-Gatto theorem ($f_+(0)\\simeq 1$). This yields the branching-ratio formula, Eq. (9), that maps PIENU positron-spectrum residuals and PIBETA diphoton-opening-angle data onto the $\\sin^2\\vartheta$--$m_a$ plane, with separate treatment of prompt and invisible ALP lifetimes.","core_discovery":"The central claim is that the rare decay $\\pi^+\\to a e\\nu$ — which inherits the pion's vector current through an off-shell $\\pi^0$ insertion proportional to $\\sin\\vartheta$ — lets one convert precision measurements of $\\pi^+\\to e\\nu$ and $\\pi^+\\to \\pi^0 e\\nu$ into direct constraints on ALP-pion mixing. The paper computes the ratio $\\mathrm{Br}[\\pi a3]/\\mathrm{Br}[\\pi e2]$ via Eq. (9), using the conserved-vector-current form-factor description of $\\pi^+\\to\\pi^0 e\\nu$, treats the ALP lifetime in prompt and invisible regimes, and uses PIENU fit residuals plus the PIBETA $\\theta_{\\gamma\\gamma}$ endpoint to set 95% CL exclusions. In the pure-mixing scenario with no direct ALP-photon coupling, these constraints combine with beam-dump, fixed-target, and collider limits to cover many decades of $\\sin^2\\vartheta$ for $10\\text{ MeV}\\lesssim m_a\\lesssim m_\\pi$.","pith_inferences":["Beyond the paper: the same mixing-template logic could be applied to other semileptonic pseudoscalar decays with well-measured $\\pi^0$ channels, such as kaon or charmed-meson decays, though backgrounds, form-factor uncertainties, and short-distance contributions would need separate evaluation.","Beyond the paper: if future $\\pi^+\\to e\\nu$ experiments release the full correlated residual covariance, the prompt-versus-invisible dichotomy in the PIENU analysis could be replaced by a single likelihood, likely sharpening the reach below $\\sin^2\\vartheta = 10^{-5}$.","Beyond the paper: the bounds assume no tree-level ALP-lepton couplings; in models with such couplings, the same spectra would instead constrain those couplings, so the $\\sin^2\\vartheta$ contours should be read as a benchmark of a restricted but representative scenario."],"forward_implications":["The PIENU positron-spectrum residuals exclude $\\sin^2\\vartheta \\gtrsim 10^{-5}$ at 95% CL for ALP masses between about 10 and 130 MeV, in both the prompt and invisible lifetime regimes.","The PIBETA diphoton opening-angle edge sets the strongest prompt-ALP bound for $m_a\\gtrsim 100$ MeV, and extending the analysis to opening angles below 160 degrees would extend the reach to lower masses.","Because the constraints arise from tree-level charged-current processes, they survive in UV models where kaon-based ALP-pion mixing bounds are spoiled by ALP-top penguin contributions.","In the pure-mixing scenario, the new data plus existing CHARM, E137, E141, and LEP limits nearly cover the $\\sin^2\\vartheta$--$m_a$ plane for $10\\text{ MeV}\\lesssim m_a\\lesssim m_\\pi$.","Improved treatments of the $\\pi^+\\to\\pi^{*0}$ form factors would extend the bounds below 10 MeV toward the massless-ALP limit, with only O(10%) theoretical uncertainty at the current mass range."],"supporting_citations":[{"why":"PIENU stopped-pion sample and detector description defining the experimental configuration.","marker":"[7]"},{"why":"High-precision PiENU measurement of Br[$\\pi^+\\to e\\nu$] that provides the normalization.","marker":"[8]"},{"why":"PDG world-average Br[$\\pi^+\\to e\\nu$] value used as the denominator in the branching-ratio formula.","marker":"[9]"},{"why":"PIENU binned low-energy $E_{\\rm cal}$ residuals from which the $\\sin^2\\vartheta$ exclusions are extracted.","marker":"[10]"},{"why":"PIBETA measurement of Br[$\\pi^+\\to\\pi^0 e\\nu$] and the $\\theta_{\\gamma\\gamma}$ opening-angle spectrum used for the diphoton bound.","marker":"[13]"},{"why":"Establishes the off-shell-$\\pi^0$ mixing amplitude estimate for $\\pi^+\\to a e\\nu$ used in Eq. (6).","marker":"[5]"},{"why":"Together with [5], provides the reference for estimating $\\pi^+\\to a e\\nu$ via pion mixing.","marker":"[6]"},{"why":"Chiral-Lagrangian matching that parametrizes ALP-pion mixing and justifies $\\sin\\vartheta$ as the phenomenological parameter.","marker":"[18]"},{"why":"Standard form-factor treatment of $\\pi^+\\to\\pi^0 e\\nu$ used to derive the rate in Eq. (7).","marker":"[19]"},{"why":"Leutwyler-Roos analysis supporting $f_+(0)\\simeq 1$ via the Ademollo-Gatto theorem, bounding the form-factor uncertainty.","marker":"[23]"}],"fun_headline_variants":["Rare pion decays set tight limits on axion-pion mixing","PIENU and PIBETA data constrain light axions","Axion-pion mixing squeezed by pion decay data","Pion decay spectra exclude mixed axions","New bounds on axion-pion mixing from pion decays"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole chain rests on the assumption that the ALP production in $\\pi^+\\to a e\\nu$ goes almost entirely through the ALP mixing with a virtual neutral pion, so that any other hadronic or ultraviolet contribution is negligible at the quoted level.","fun_headline_variants_meta":{"raw":{"variants":["Rare pion decays set tight limits on axion-pion mixing","PIENU and PIBETA data constrain light axions","Axion-pion mixing squeezed by pion decay data","Pion decay spectra exclude mixed axions","New bounds on axion-pion mixing from pion decays"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000312,"raw_usage":{"total_tokens":1719,"prompt_tokens":830,"completion_tokens":889,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":446,"completion_tokens_details":{"reasoning_tokens":809}},"tokens_in":446,"tokens_out":889,"duration_ms":8123,"temperature":1.0,"reasoning_tokens":809,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T10:05:33.108386+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A dedicated stopped-pion measurement of $\\pi^+\\to (a\\to\\gamma\\gamma)e\\nu$ reconstructing the diphoton opening angle for $10\\text{ MeV}\\lesssim m_a\\lesssim m_\\pi$ would settle the central claim: the predicted opening-angle distribution beyond the $\\pi^0$ edge, Eq. (11), should appear at the rate set by Eq. (9) for the claimed $\\sin^2\\vartheta$, and its absence would falsify the off-shell-$\\pi^0$ dominance assumption.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"High-precision PiENU measurement of Br[$\\pi^+\\to e\\nu$] that provides the normalization."},{"cited_title":"Search for Massive Neutrinos in the Decay pi --> e nu","cited_arxiv_id":"1106.4055","evidence_quote":"PIENU binned low-energy $E_{\\rm cal}$ residuals from which the $\\sin^2\\vartheta$ exclusions are extracted."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Together with [5], provides the reference for estimating $\\pi^+\\to a e\\nu$ via pion mixing."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Standard form-factor treatment of $\\pi^+\\to\\pi^0 e\\nu$ used to derive the rate in Eq. (7)."},{"cited_title":"Cowan, H","cited_arxiv_id":null,"evidence_quote":"Leutwyler-Roos analysis supporting $f_+(0)\\simeq 1$ via the Ademollo-Gatto theorem, bounding the form-factor uncertainty."}],"review_version":1}