K± → π±a at next-to-leading order in chiral perturbation theory and updated bounds on ALP couplings

The weak decays K- +/- -> pi (+/-) a offer a powerful probe of axion-like particles (ALPs). In this work, we provide a comprehensive analysis of these processes within chiral perturbation theory, extending existing calculations by including complete next-to-leading order (NLO) contributions and isospin-breaking corrections at first order in (m(d) - m(u) ). We show that the consistent incorporation of ALPs in the QCD and weak chiral Lagrangians requires a non-trivial extension of the corresponding operator bases, which we describe in detail. Furthermore, we show that in the presence of an ALP the so-called "weak mass term", which is unobservable in the Standard Model, is non-redundant already at leading order. We find that NLO corrections associated with flavor-violating ALP couplings modify the leading-order result by a few percent, with negligible uncertainties. NLO corrections proportional to flavor-conserving ALP couplings lead to potentially larger corrections, which, however, are accompanied by sizable uncertainties mainly due to the currently limited knowledge of various low-energy constants. We study how these corrections impact bounds on the ALP couplings, first model independently, and then specializing to the case of an ALP with flavor-universal couplings in the UV. Our findings confirm that the decays K- +/- -> pi (+/-) a provide the strongest particle-physics constraints for m a less than or similar to 300 MeV. In addition, we point out that these bounds have interesting implications for the ALP couplings to nucleons, which were so far only constrained by astrophysical measurements and non-accelerator experiments.