arxiv: 2604.25673 · v1 · submitted 2026-04-28 · ✦ hep-ph

Recognition: unknown

Axion-like particle-meson production in semileptonic τ decays

Jin Hao, Yu-Xuan Bai, Zhi-Hui Guo

Authors on Pith no claims yet

Pith reviewed 2026-05-07 16:04 UTC · model grok-4.3

classification ✦ hep-ph

keywords axion-like particlessemileptonic tau decayschiral effective field theorymeson mixingbranching ratiosinvariant mass distributionsforward-backward asymmetryisospin breaking

0 comments

The pith

Chiral effective theory supplies concrete predictions for axion-like particle production in semileptonic tau decays after fitting resonance parameters to standard spectra.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper calculates rates and kinematic distributions for tau decays that produce an axion-like particle together with a pion or kaon. It builds a next-to-leading-order mixing matrix among the neutral mesons and the axion-like particle that includes linear isospin breaking, then uses this matrix to obtain the hadronic form factors that enter the decay amplitudes. Resonance parameters in those form factors are fixed by fits to three well-measured tau decay channels that do not involve the axion-like particle. The resulting numerical predictions for branching ratios, invariant-mass spectra, and forward-backward asymmetries give definite targets that experiments can search for. A reader would care because these predictions turn a general theoretical possibility into a set of observable numbers in a clean leptonic environment.

Core claim

In this work we explore the semileptonic τ decays into the axion-like particle (a)-meson final states within chiral effective field theory. The next-to-leading-order mixing matrix for the π⁰-η-η′-a system with the linear isospin-breaking effects is exploited and then implemented to calculate the hadronic form factors relevant to the τ decays. The resonance parameters entering the form factors are determined from fits to the experimental spectra of τ⁻ → π⁻ π⁰ ν_τ, τ⁻ → K_S π⁻ ν_τ, and τ⁻ → K⁻ η ν_τ. We then focus on the predictions to the branching ratios, invariant-mass distributions, and forward-backward asymmetries from the τ⁻ → P a ν_τ processes, with P = π⁻ and K⁻.

What carries the argument

The next-to-leading-order mixing matrix among π⁰, η, η′, and the axion-like particle that incorporates linear isospin-breaking effects and determines the hadronic form factors used in the decay amplitudes.

If this is right

Numerical values for the branching ratios of τ⁻ → π⁻ a ν_τ and τ⁻ → K⁻ a ν_τ follow directly from the fitted form factors.
The invariant-mass distributions of the π a and K a systems acquire specific shapes set by the resonance content of the form factors.
Forward-backward asymmetries in the same decays become additional observables that experiments can measure to isolate signals.
The same mixing matrix and form-factor framework can be reused for any other semileptonic tau channel involving the axion-like particle.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

These predictions supply ready-made signal templates that future high-statistics tau data sets can use to set limits or discover an axion-like particle.
If the measured distributions deviate from the calculated ones, the discrepancy could point either to a breakdown of the effective-theory assumptions or to additional new-physics contributions not included in the mixing matrix.
The method of fitting resonance parameters in the absence of the new particle and then inserting it into the mixing matrix offers a template that could be applied to other light pseudoscalar searches in meson decays.

Load-bearing premise

Resonance parameters extracted from ordinary tau decays without axion-like particles continue to apply once the axion-like particle is added to the meson mixing matrix.

What would settle it

A measured branching ratio or invariant-mass distribution for τ⁻ → π⁻ a ν_τ that lies well outside the range computed from the fitted form factors would falsify the prediction.

Figures

Figures reproduced from arXiv: 2604.25673 by Jin Hao, Yu-Xuan Bai, Zhi-Hui Guo.

**Figure 1.** Figure 1: Differential decay widths for τ − → π−aντ (left) and τ− → K−aντ (right) as functions of √ s, the invariant mass of the ALP-meson system. Results are shown for ALP masses with ma = 0.05, 0.1, 0.15, 0.2, and 0.3 GeV, the QCD-axion benchmark with ma = 0 is also displayed for comparison. metries AF B, defined in Eq. (12). Here we focus on the distributions of AF B as functions of the invariant mass of the ALP-… view at source ↗

**Figure 2.** Figure 2: Forward-backward asymmetries AF B for τ− → π−aντ (left) and τ− → K−aντ (right) as functions of √ s, the invariant mass of the ALP-meson system. Results are shown for ALP masses with ma = 0.05, 0.1, 0.15, 0.2, and 0.3 GeV, the QCD-axion benchmark with ma = 0 is also displayed for comparison. the forward-backward asymmetries for the τ − → π −/K−aντ processes, by taking several different values for the ALP ma… view at source ↗

read the original abstract

In this work we explore the semileptonic $\tau$ decays into the axion-like particle ($a$)-meson final states within chiral effective field theory. The next-to-leading-order mixing matrix for the $\pi^0$-$\eta$-$\eta'$-$a$ system with the linear isospin-breaking effects, is exploited and then implemented to calculate the hadronic form factors relevant to the $\tau$ decays. The resonance parameters entering the form factors are determined from fits to the experimental spectra of $\tau^- \to \pi^- \pi^0 \nu_\tau$, $\tau^- \to K_S \pi^- \nu_\tau$, and $\tau^- \to K^- \eta \nu_\tau$. We then focus on the predictions to the branching ratios, invariant-mass distributions, and forward-backward asymmetries from the $\tau^-\to P a \nu_\tau$ processes, with $P=\pi^-$ and $K^-$. Our results provide a quantitative basis for future searches of the axion-like particle signals in semileptonic $\tau$ decays.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper supplies explicit predictions for ALP branching ratios and spectra in tau to pi/K plus ALP decays by fitting resonances to standard channels and inserting them into an NLO mixing matrix.

read the letter

The main takeaway is that this work produces concrete numbers for branching ratios, invariant-mass distributions, and forward-backward asymmetries in τ− → π− a ντ and τ− → K− a ντ. The authors extend the usual chiral EFT treatment of semileptonic tau decays by adding the ALP to the π0-η-η' mixing matrix at next-to-leading order with linear isospin breaking, then reuse resonance parameters fitted to the three measured channels τ− → π−π0 ντ, τ− → KS π− ντ, and τ− → K− η ντ to generate the new form factors and predictions. Those specific ALP-channel results are not already in the cited literature, so the quantitative basis for searches is genuinely new material. The approach is standard and the abstract shows they follow the usual workflow of fitting first, then predicting. The fits to real data give the predictions some anchoring, and the focus on measurable distributions makes the output directly usable by experimental groups. The soft spot is the untested assumption that the resonance parameters stay fixed once the ALP enters the mixing. The paper fits those parameters in the absence of the ALP and then inserts the ALP without a re-fit or sensitivity check, so any back-reaction on the effective form factors is not quantified. At the same time, NLO chiral EFT with only linear isospin breaking is applied near the tau mass, where resonance effects and higher-order terms can be sizable; the linear approximation may not fully capture the dynamics. These issues are not fatal but mean the quoted predictions carry an unstated systematic uncertainty. The paper is aimed at ALP phenomenologists and tau-decay experimentalists who need ready-to-use numbers for search strategies. A reader working on ALP constraints or on improving hadronic form factors in tau decays will get value from the explicit results. The calculation is competent and grounded enough to deserve a serious referee who can verify the mixing-matrix implementation and the fit stability.

Referee Report

2 major / 2 minor

Summary. The manuscript calculates axion-like particle (ALP) production in association with a pseudoscalar meson in semileptonic tau decays (τ⁻ → P a ν_τ with P = π⁻, K⁻) within next-to-leading-order chiral effective field theory. It constructs the NLO mixing matrix for the π⁰-η-η'-a system including linear isospin-breaking effects, derives the relevant hadronic form factors, determines resonance parameters by fitting to the measured spectra of three standard (ALP-free) tau decays, and then predicts branching ratios, invariant-mass distributions, and forward-backward asymmetries for the ALP channels.

Significance. If the central assumptions hold, the work supplies concrete, data-constrained predictions that could serve as a quantitative reference for experimental searches of ALPs in tau decays at Belle II or future facilities. The use of resonance parameters extracted directly from experiment and the systematic inclusion of NLO mixing are methodological strengths that align with standard practice in hadronic tau physics.

major comments (2)

[Description of resonance-parameter fits and subsequent ALP predictions] The resonance parameters (masses, widths, and couplings) are extracted from fits to the experimental spectra of τ⁻ → π⁻π⁰ν_τ, τ⁻ → K_S π⁻ν_τ, and τ⁻ → K⁻η ν_τ performed in the absence of the ALP. These parameters are then inserted unchanged into the form factors that incorporate the ALP through the NLO π⁰-η-η'-a mixing matrix. No re-fit, stability test, or sensitivity analysis is presented to verify that the parameters remain appropriate once the ALP modifies the mixing and the effective form-factor expressions. Because the predicted branching ratios, mass spectra, and asymmetries rest directly on these parameters, this assumption is load-bearing for the central quantitative results.
[Chiral Lagrangian and mixing-matrix construction] The NLO chiral EFT with linear isospin breaking is applied at energies up to m_τ ≈ 1.78 GeV. The manuscript does not quantify the expected truncation error or discuss whether resonance saturation and higher-order terms remain adequately captured when the ALP is included in the mixing matrix. This truncation directly affects the form factors used for the ALP channels.

minor comments (2)

[Numerical results and figures] The figures showing invariant-mass distributions and forward-backward asymmetries would be strengthened by including uncertainty bands propagated from the fitted resonance parameters.
[Throughout the text] Notation for the ALP mass and coupling constant should be defined once and used consistently in all equations and figure labels.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive assessment of the manuscript and for the constructive major comments. We address each point below and will incorporate revisions to strengthen the analysis and presentation.

read point-by-point responses

Referee: The resonance parameters (masses, widths, and couplings) are extracted from fits to the experimental spectra of τ⁻ → π⁻π⁰ν_τ, τ⁻ → K_S π⁻ν_τ, and τ⁻ → K⁻η ν_τ performed in the absence of the ALP. These parameters are then inserted unchanged into the form factors that incorporate the ALP through the NLO π⁰-η-η'-a mixing matrix. No re-fit, stability test, or sensitivity analysis is presented to verify that the parameters remain appropriate once the ALP modifies the mixing and the effective form-factor expressions. Because the predicted branching ratios, mass spectra, and asymmetries rest directly on these parameters, this assumption is load-bearing for the central quantitative results.

Authors: We appreciate the referee identifying this key assumption. The resonance parameters are fixed by fits to the three standard tau decay channels, which directly constrain the vector and axial-vector resonance contributions. The ALP enters the calculation exclusively via the NLO pseudoscalar mixing matrix, which rescales the effective hadronic currents without modifying the resonance propagators or their intrinsic parameters. This separation follows from the chiral power counting. To quantify the robustness of this procedure, we will add a dedicated sensitivity analysis in the revised manuscript: the resonance parameters will be varied within their experimental fit uncertainties, and the resulting bands will be shown on the predicted ALP branching ratios, invariant-mass distributions, and forward-backward asymmetries. revision: yes
Referee: The NLO chiral EFT with linear isospin breaking is applied at energies up to m_τ ≈ 1.78 GeV. The manuscript does not quantify the expected truncation error or discuss whether resonance saturation and higher-order terms remain adequately captured when the ALP is included in the mixing matrix. This truncation directly affects the form factors used for the ALP channels.

Authors: The referee correctly notes the energy reach of the calculation. Our framework employs resonance saturation of higher-order counterterms, a standard and validated approach in hadronic tau physics that successfully reproduces the reference decay spectra. Because the ALP is treated as an additional light pseudoscalar with mixing generated at NLO, the truncation error is expected to be comparable to that in the ALP-free channels. In the revised manuscript we will add an explicit discussion of truncation uncertainties, including a qualitative estimate based on the relative magnitude of NLO corrections observed in the fits and the size of resonance contributions. revision: yes

Circularity Check

0 steps flagged

No significant circularity: resonance parameters fitted to external ALP-free data then transferred to distinct ALP channels

full rationale

The paper determines resonance parameters by fitting to experimental spectra of the standard channels τ⁻→π⁻π⁰ν_τ, τ⁻→K_S π⁻ν_τ and τ⁻→K⁻η ν_τ (explicitly without ALP), then inserts those fixed parameters into the NLO π⁰-η-η'-a mixing matrix and form factors to generate predictions for the separate processes τ⁻→P a ν_τ (P=π⁻,K⁻) where ALP mass and coupling remain free parameters. This is a conventional transfer of external constraints to new-physics observables; the output distributions and branching ratios are not forced by construction to reproduce the input spectra, nor does any equation reduce the ALP predictions to the fitted inputs. No load-bearing self-citation, ansatz smuggling, or uniqueness theorem is invoked in the provided derivation chain.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 1 invented entities

The central claim rests on the validity of NLO chiral EFT for the mixed meson-ALP system, the applicability of resonance parameters fitted to standard decays, and the linear isospin-breaking approximation in the mixing matrix. No machine-checked proofs or independent data releases are provided.

free parameters (1)

resonance parameters (masses, widths, couplings)
Determined from fits to the experimental spectra of τ−→π−π0ντ, τ−→KSπ−ντ, and τ−→K−ηντ as stated in the abstract.

axioms (2)

domain assumption Next-to-leading-order mixing matrix for the π0-η-η'-a system with linear isospin-breaking effects
Exploited to calculate the hadronic form factors relevant to the τ decays.
domain assumption Chiral effective field theory framework accurately describes the relevant hadronic interactions
Used throughout for form factor construction and decay rate calculations.

invented entities (1)

axion-like particle a no independent evidence
purpose: To model its production in semileptonic tau decays
Hypothetical particle introduced into the mixing matrix and form factors; no independent evidence or falsifiable prediction outside the model is supplied.

pith-pipeline@v0.9.0 · 5501 in / 1630 out tokens · 59234 ms · 2026-05-07T16:04:32.352968+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

31 extracted references · 30 canonical work pages · 1 internal anchor

[1]

R. D. Peccei and H. R. Quinn,Phys. Rev. Lett.38, 1440 (1977), doi:10.1103/ PhysRevLett.38.1440

1977
[2]

R. D. Peccei and H. R. Quinn,Phys. Rev. D16, 1791 (1977), doi:10.1103/PhysRevD. 16.1791

work page doi:10.1103/physrevd 1977
[3]

Weinberg, Phys

S. Weinberg,Phys. Rev. Lett.40, 223 (1978), doi:10.1103/PhysRevLett.40.223

work page doi:10.1103/physrevlett.40.223 1978
[4]

Wilczek, Phys

F. Wilczek,Phys. Rev. Lett.40, 279 (1978), doi:10.1103/PhysRevLett.40.279

work page doi:10.1103/physrevlett.40.279 1978
[5]

The landscape of QCD axion models

L. Di Luzio, M. Giannotti, E. Nardi and L. Visinelli,Phys. Rept.870, 1 (2020), arXiv:2003.01100 [hep-ph], doi:10.1016/j.physrep.2020.06.002

work page doi:10.1016/j.physrep.2020.06.002 2020
[6]

J. E. Kim and G. Carosi,Rev. Mod. Phys.82, 557 (2010),arXiv:0807.3125 [hep-ph], doi:10.1103/RevModPhys.82.557, [Erratum: Rev.Mod.Phys. 91, 049902 (2019)]

work page doi:10.1103/revmodphys.82.557 2010
[7]

P. W. Graham, I. G. Irastorza, S. K. Lamoreaux, A. Lindner and K. A. van Bibber, Ann. Rev. Nucl. Part. Sci.65, 485 (2015),arXiv:1602.00039 [hep-ex], doi:10.1146/ annurev-nucl-102014-022120

work page arXiv 2015
[8]

I. G. Irastorza and J. Redondo,Prog. Part. Nucl. Phys.102, 89 (2018), arXiv:1801.08127 [hep-ph], doi:10.1016/j.ppnp.2018.05.003

work page doi:10.1016/j.ppnp.2018.05.003 2018
[9]

K. Choi, S. H. Im and C. Sub Shin,Ann. Rev. Nucl. Part. Sci.71, 225 (2021), arXiv:2012.05029 [hep-ph], doi:10.1146/annurev-nucl-120720-031147

work page doi:10.1146/annurev-nucl-120720-031147 2021
[10]

Achasov, et al., STCF conceptual design report: V olume I – Physics & detector, Front

M. Achasovet al.,Front. Phys. (Beijing)19, 14701 (2024),arXiv:2303.15790 [hep-ex], doi:10.1007/s11467-023-1333-z. April 29, 2026 1:2 ws-ijmpa Axion-like particle-meson production in semileptonicτdecays9

work page doi:10.1007/s11467-023-1333-z 2024
[11]

Chenget al.(2 2025),arXiv:2502.08907 [hep-ex]

H.-Y. Chenget al.(2 2025),arXiv:2502.08907 [hep-ex]

work page arXiv 2025
[12]

Hao, C.-G

J. Hao, C.-G. Duan and Z.-H. Guo,Front. Phys. (Beijing)21, 076201 (2026), arXiv:2507.00383 [hep-ph], doi:10.15302/frontphys.2026.076201

work page doi:10.15302/frontphys.2026.076201 2026
[13]

Z.-X. Li, A. Li, J. Hao, C.-G. Duan and Z.-H. Guo,Phys. Rev. D113, 014005 (2026), arXiv:2510.04172 [hep-ph], doi:10.1103/j9wc-wgd9

work page doi:10.1103/j9wc-wgd9 2026
[14]

Escribano, S

R. Escribano, S. Gonzalez-Solis and P. Roig,Phys. Rev. D94, 034008 (2016), arXiv:1601.03989 [hep-ph], doi:10.1103/PhysRevD.94.034008

work page doi:10.1103/physrevd.94.034008 2016
[15]

D. A. L. Aguilar, J. Rend´ on and P. Roig,JHEP01, 105 (2025),arXiv:2409.05588 [hep-ph], doi:10.1007/JHEP01(2025)105

work page doi:10.1007/jhep01(2025)105 2025
[16]

Gao, Z.-H

R. Gao, Z.-H. Guo, J. A. Oller and H.-Q. Zhou,JHEP04, 022 (2023), arXiv:2211.02867 [hep-ph], doi:10.1007/JHEP04(2023)022

work page doi:10.1007/jhep04(2023)022 2023
[17]

R. Gao, J. Hao, C.-G. Duan, Z.-H. Guo, J. A. Oller and H.-Q. Zhou,Eur. Phys. J. C 85, 97 (2025),arXiv:2411.06737 [hep-ph], doi:10.1140/epjc/s10052-025-13807-9

work page doi:10.1140/epjc/s10052-025-13807-9 2025
[18]

Ecker, J

G. Ecker, J. Gasser, A. Pich and E. de Rafael,Nucl. Phys. B321, 311 (1989), doi: 10.1016/0550-3213(89)90346-5

work page doi:10.1016/0550-3213(89)90346-5 1989
[19]

di Cortona, E

G. Grilli di Cortona, E. Hardy, J. Pardo Vega and G. Villadoro,JHEP01, 034 (2016), arXiv:1511.02867 [hep-ph], doi:10.1007/JHEP01(2016)034

work page doi:10.1007/jhep01(2016)034 2016
[20]

Kaiser and H

R. Kaiser and H. Leutwyler,Eur. Phys. J. C17, 623 (2000),arXiv:hep-ph/0007101, doi:10.1007/s100520000499

work page doi:10.1007/s100520000499 2000
[21]

Guo, Z.-H

X.-K. Guo, Z.-H. Guo, J. A. Oller and J. J. Sanz-Cillero,JHEP06, 175 (2015), arXiv:1503.02248 [hep-ph], doi:10.1007/JHEP06(2015)175

work page doi:10.1007/jhep06(2015)175 2015
[22]

Erler,Electroweak radiative corrections to semileptonic tau decays, Rev

J. Erler,Rev. Mex. Fis.50, 200 (2004),arXiv:hep-ph/0211345

work page internal anchor Pith review arXiv 2004
[23]

Navas et al

Particle Data Group Collaboration (S. Navaset al.),Phys. Rev. D110, 030001 (2024), doi:10.1103/PhysRevD.110.030001

work page doi:10.1103/physrevd.110.030001 2024
[24]

Fujikawaet al.),Phys

Belle Collaboration (M. Fujikawaet al.),Phys. Rev. D78, 072006 (2008), arXiv:0805.3773 [hep-ex], doi:10.1103/PhysRevD.78.072006

work page doi:10.1103/physrevd.78.072006 2008
[25]

Epifanovet al.),Phys

Belle Collaboration (D. Epifanovet al.),Phys. Lett. B654, 65 (2007), arXiv:0706.2231 [hep-ex], doi:10.1016/j.physletb.2007.08.045

work page doi:10.1016/j.physletb.2007.08.045 2007
[26]

Inamiet al.),Phys

Belle Collaboration (K. Inamiet al.),Phys. Lett. B672, 209 (2009),arXiv:0811.0088 [hep-ex], doi:10.1016/j.physletb.2009.01.047

work page doi:10.1016/j.physletb.2009.01.047 2009
[27]

Wang, Z.-H

J.-B. Wang, Z.-H. Guo and H.-Q. Zhou,Phys. Rev. D109, 075030 (2024), arXiv:2312.15240 [hep-ph], doi:10.1103/PhysRevD.109.075030

work page doi:10.1103/physrevd.109.075030 2024
[28]

Wang, Z.-H

J.-B. Wang, Z.-H. Guo and H.-Q. Zhou,Phys. Rev. D112, L031701 (2025), arXiv:2505.21276 [hep-ph], doi:10.1103/xsvj-g811

work page doi:10.1103/xsvj-g811 2025
[29]

Wang, Z.-H

J.-B. Wang, Z.-H. Guo, Z. Lu and H.-Q. Zhou,JHEP11, 029 (2024), arXiv:2403.16064 [hep-ph], doi:10.1007/JHEP11(2024)029

work page doi:10.1007/jhep11(2024)029 2024
[30]

Cao and Z.-H

X.-H. Cao and Z.-H. Guo,Phys. Rev. D110, 095025 (2024),arXiv:2408.15825 [hep-ph], doi:10.1103/PhysRevD.110.095025

work page doi:10.1103/physrevd.110.095025 2024
[31]

Guo, X.-H

Q.-Q. Guo, X.-H. Cao, Z.-H. Guo and H.-Q. Zhou,Phys. Rev. D113, 034030 (2026), arXiv:2504.06786 [hep-ph], doi:10.1103/pw8q-xknh

work page doi:10.1103/pw8q-xknh 2026