Study of the interactions of the axion with mesons and photons using a chiral effective Lagrangian model

We investigate the most interesting decay processes involving axions, photons and the lightest pseudoscalar mesons, making use of a chiral effective Lagrangian model with L=3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L=3$$\end{document} light quark flavors, which also includes the flavor-singlet pseudoscalar meson and implements the U(1) axial anomaly of the fundamental theory. In particular, we compute the electromagnetic coupling of the axion to photons and we compare our result with the prediction of the Chiral Effective Lagrangian with L=2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L=2$$\end{document} light quark flavors. Moreover, we study the decay channels η/η′→ππa\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta /\eta '\rightarrow \pi \pi a$$\end{document} and we estimate the corresponding decay widths, using the existing bounds on the U(1)PQ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U(1)_{PQ}$$\end{document} breaking scale.


Introduction
It is well known that the QCD Lagrangian L QC D can be extended by adding the term L θ = θ Q, where Q = g 2 64π 2 ε μνρσ G a μν G a ρσ is the so-called topological charge density and θ is a free parameter, which can assume any value in [0, 2π). This θ -term (or topological term) introduces an explicit breaking of the CP symmetry in the strong sector (referred to as strong-CP violation). Despite the fact that Q = ∂ μ K μ , where K μ is the so-called Chern-Simons current, its contribution is nonzero because of topologically nontrivial configurations of gauge fields, such as instantons. So far, however, no CP violation in the strong sector has been observed experimentally, constraining θ to be zero or extremely small. In particular, one can find a relationship between θ and the neutron electric-dipole moment [1], 10 −16 |θ | e · cm, where m N is the neutron mass, whereas m π is the pion mass. From experimental data [2] we know that d N < 10 −26 e · cm, which leads to the upper bound |θ | < 10 −10 . (More refined relations among the a e-mail: giacomo.landini@phd.unipi.it b e-mail: enrico.meggiolaro@unipi.it (corresponding author) neutron electric dipole moment and the θ angle and a more detailed discussion can be found in Refs. [3][4][5][6]; see also Ref. [7] for a recent lattice determination). This "fine-tuning" problem, known in the literature as the "strong-CP problem", is one of the open issues of the Standard Model. Among the several possible solutions, the most appealing is surely the one proposed by Peccei and Quinn (PQ) [8,9] and developed by Weinberg and Wilczek [10,11].
The key idea (see also Ref. [12] for a recent review) is to extend the Standard Model by adding a new pseudoscalar particle, called "axion", in such a way that there is a new U (1) global symmetry, referred to as U (1) P Q , which is both spontaneously broken at a scale f a and anomalous (i.e., broken by quantum effects), with the related current satisfying the relation ∂ μ J μ P Q = a P Q Q, where a P Q is the so-called color anomaly parameter. The most general Lagrangian describing the QCD degrees of freedom and the axion has the following form: where S a is the axion field, which under U (1) P Q transforms nonlinearly as The term L int [∂ μ S a , ] describes the interactions between the axion and the quark fields and it is strongly model dependent. The effect of this extension of the Lagrangian is to replace the static θ parameter of L QC D with a dynamical degree of freedom, namely the combination θ − a P Q S a f a : on the vacuum we get θ − a P Q S a f a = 0. Performing a U (1) P Q transformation (1.2) with γ = θ a P Q , we can rotate away the θ term, so obtaining a manifestly CP-conserving theory (with S a = 0). Moreover, it is well known that the U (1) axial symmetry of QCD with L light quark flavors (taken to be massless in the ideal chiral limit; the physically relevant cases are L = 2, with the quarks up and down, and L = 3, including also the strange quark), U (1)A : q i → q i = e iβγ 5 q i , i = 1, . . . , L , (1.3) is also anomalous, with the related U (1) axial current J μ 5 = qγ μ γ 5 q satisfying the relation ∂ μ J μ 5 = 2L Q. Therefore, we find that the U (1) A ⊗ U (1) P Q transformations with the parameters β and γ satisfying the constraint 2Lβ + a P Q γ = 0, form a U (1) subgroup which is spontaneously broken but anomaly-free (in the chiral limit): as a consequence, a new (pseudo-)Nambu-Goldstone boson appears in the spectrum, the axion.
In the original Peccei-Quinn-Weinberg-Wilczek (PQWW) model [8,10,11] the scale f a was identified with the electroweak breaking scale v ≈ 250 GeV, but this leads to large couplings between the axion and the Standard Model fields, which have been ruled out by experiments (see, for example, Ref. [13]). In order to bypass these experimental bounds, the so-called "invisible axion" models were developed, such as the Kim-Shifman-Vainshtein-Zakharov (KSVZ) model [14,15] and the Dine-Fischler-Srednicki-Zhitnisky (DFSZ) model [16,17], in which new heavy quarks or scalar fields, charged under U (1) P Q but neutral with respect to the Standard Model gauge group, are introduced. In these models, the U (1) P Q breaking scale f a is a free parameter of the theory and, assuming f a v, a very light axion with small couplings to the Standard Model fields is predicted, a scenario which is still compatible with the experimental bounds. At present, the more precise bounds on the U (1) P Q breaking scale come from astrophysical and cosmological considerations (see, for example, Refs. [18,19]): 10 9 GeV f a 10 17 GeV. All these models predict an axion-photon-photon coupling and therefore the electromagnetic decay of the axion in two photons: most of the experimental research concerning the axion is focused on this process (see, for example, Ref. [20] for an exhaustive review on both the theoretical aspects and the experimental research of axions and axion-like particles). The electromagnetic interaction of the axion is usually parametrized as where a is the "physical" axion [as we will see, the field S a , which appears in the Lagrangian (1.1), has nonzero mixings with the QCD degrees of freedom, such as the pseudoscalar meson fields], F μν is the electromagnetic field-strength tensor,F μν = 1 2 ε μνρσ F ρσ is its dual, and g aγ γ is the axionphoton-photon coupling constant. This last, in general, is the sum of two contributions, g aγ γ = g 0 aγ γ +g QC D aγ γ , where g 0 aγ γ is the model-dependent contribution proportional to the electromagnetic anomaly of the U (1) P Q symmetry (which can also be simply zero, as it happens in the original KSVZ model [14,15]), while g QC D aγ γ is the model-independent contribution coming from the minimal coupling to QCD (i.e., the mixing of the axion with the pseudoscalar mesons π 0 , η, and η ). The coupling constant g QC D aγ γ has been computed using the Chiral Effective Lagrangian with L = 2 light quark flavors both at the leading order (LO) in the momentum expansion [O( p 2 )] and at the next-to-leading order (NLO) [O( p 4 ] (see Ref. [21] and references therein).
The aim of this paper is to compute the axion-photonphoton coupling constant g QC D aγ γ and, moreover, to study the decay processes involving the axion and the lightest pseudoscalar mesons, making use of a chiral effective Lagrangian model proposed by Witten et al. [22][23][24][25][26][27][28], which describes the Nambu-Goldstone bosons originated by the spontaneous breaking of the SU (3) L ⊗ SU (3) R chiral symmetry (with L = 3 light quark flavors) and the flavor-singlet pseudoscalar meson, implementing the U (1) axial anomaly of the fundamental theory.
In Sect. 2, for the benefit of the reader, we briefly recall this chiral effective Lagrangian model, as well as its "axionized" version (see Ref. [29]).
Using this model, in Sect. 3 we compute the axion-photonphoton coupling constant g QC D aγ γ and the result is compared with the one obtained using the Chiral Effective Lagrangian with L = 2 light quark flavors. Section 4 is devoted to the study of the hadronic decays η/η → ππa (which, of course, cannot be studied using the Chiral Effective Lagrangian with L = 2 light quark flavors, since the η and η degrees of freedom are integrated out). Among all the possible hadronic decays involving also the axion, these are the ones involving the lowest-energy hadrons.
In Sect. 5 we study the effects of a possible U (1) axial condensate on the various quantities that we have evaluated in the previous sections: we do this by using a chiral effective Lagrangian model proposed in Refs. [30][31][32] and then elaborated on in Refs. [33][34][35][36][37], which can be interpreted as an extension of the model considered in Sect. 2 with the inclusion of a U (1) axial condensate.
In Sect. 6 we report numerical estimates for the axionphoton-photon coupling constant (making also a comparison with the prediction of the Chiral Effective Lagrangian with L = 2 light quark flavors) and for the widths of the hadronic decays η/η → ππa.
Finally, in Sect. 7 we summarize and critically comment on the results that we have obtained in the previous sections for the electromagnetic and the hadronic processes involving the axion (considering the existing bounds on the U (1) P Q breaking scale) and we also give some prospects for further theoretical and experimental studies of the hadronic decays η/η → ππa. The effective Lagrangian model proposed by Witten et al. [22][23][24][25][26][27][28] describes the Nambu-Goldstone bosons originated by the spontaneous breaking of the SU (3) L ⊗ SU (3) R chiral symmetry and the flavor-singlet pseudoscalar meson, implementing the U (1) axial anomaly of the fundamental theory.
We will refer to it as the "WDV model". Even if this model was derived and fully justified in the large-N c limit (N c being the number of colors), the numerical results obtained for the physical value N c = 3 turn out to be quite consistent with experimental data. The Lagrangian is given by (see Ref. [23] for a detailed discussion) (2.1) The mesonic field U is represented by a 3×3 complex matrix, which can be written in terms of the quark fields as U i j ∼ q j R q i L , up to a multiplicative constant. 1 Under a general At zero temperature (after integrating out the scalar meson fields) we can adopt the usual nonlinear parametrization: where λ a (a = 1, . . . , 8) are the usual generators of SU (3) (Gell-Mann matrices), normalized so as Tr [λ a λ b ] = 2δ ab , and π a (x) are the nonsinglet pseudoscalar-meson fields, while S(x) is the flavor-singlet pseudoscalar-meson field. Moreover: • F π is the pion decay constant.
• B is a constant (with the dimension of a mass) which relates the squared masses of the pseudoscalar mesons and the quark masses: for example, m 2 π = B(m u + m d ). 1 Throughout this paper, we shall use the following notations for the left-handed and right-handed quark fields: q L ,R ≡ 1 2 (1 ± γ 5 )q, with γ 5 ≡ −iγ 0 γ 1 γ 2 γ 3 . Moreover, we shall adopt the convention ε 0123 = −ε 0123 = 1 for the (Minkowskian) completely antisymmetric tensor ε μνρσ (= −ε μνρσ ) which appears in the expressions of the topological charge density Q and of the dual electromagnetic field-strength tensor F μν .
The topological charge density Q is introduced as an auxiliary field, whereas A is a parameter which (at least in the large-N c limit) can be identified with the topological susceptibility in the pure Yang-Mills theory It is easy to see that the anoma- so correctly reproducing the U (1) axial anomaly of the fundamental theory. The model can be easily "axionized" (see Ref. [29]), essentially by promoting the parameter θ to the axion field S a (apart from a multiplicative constant a P Q / f a ) and adding a kinetic term for it, i.e., where N = f a e i Sa fa parameterizes the axion field in the standard notation for Nambu-Goldstone bosons. It is convenient to integrate out the auxiliary field Q using its equations of motion: The resulting Lagrangian is given by (2.7) Expanding the Lagrangian up to the second order in the fields, we get the following squared-mass matrix for the fields π 3 , π 8 , S, S a (the mass term for the fields π 1 , π 2 , π 4 , π 5 , π 6 , and π 7 being diagonal): where we have defined The fields π 3 , π 8 , S, S a can be written in terms of the "physical" fields π 0 , η, η , a, associated with the mass eigenstates of Eq. (2.8), as follows: where θ i j is an orthogonal mixing matrix. From the astrophysical bounds on the scale f a [18,19] (or better on f a /a P Q , but a P Q ∼ O(1) for the more realistic axion models [38]) we have: 10 −18 b 10 −10 . As a consequence, it is surely legitimate to perform the computations only at the leading order in b. In particular, diagonalizing the squared-mass matrix, we can derive the following expression (at the leading order in b) for the squared mass of the axion [29]: This expression is in perfect agreement with the well-known relationship (valid at the leading order in b) [15]: between the squared mass of the axion and the topological susceptibility of QCD, considering the expression of χ QC D which is found using the WDV model (see Refs. [29,39] and references therein).

Electromagnetic decay of the axion
In order to investigate the electromagnetic decay of the axion, we have to introduce the electromagnetic interactions into the Lagrangian (2.7). This is done by (i) replacing the derivative of the field U with the corresponding covariant derivative where G = e 2 N C 32π 2 ε μνρσ F μν F ρσ , F μν being the electromagnetic field-strength tensor. Using Eq. (2.3), this term can be rewritten explicitly in terms of the meson fields, as follows: Making use of Eq. (2.10), one immediately sees that this term contains an axion-photon-photon interaction of the type (1.4), with the following expression for the axion-photonphoton coupling constant: 2 where α e.m. = e 2 4π 1 137 is the fine-structure constant. To find the mixing parameters in Eq. (2.10), we have to solve the equations for the eigenvectors of the matrix (2.8).
In particular, using the following notation: the axion eigenvector is given by First we shall derive our expressions for = 0, i.e., neglecting the experimentally small violations of the SU (2) V isospin symmetry. For = 0 the mass matrix becomes diagonal with respect to π 3 , which can thus be identified with π 0 : therefore θ π 3 S a | =0 = 0. The eigenvector equations are, in this case: 2 As we have already said in the Introduction, this is indeed the modelindependent contribution g QC D aγ γ coming from the minimal coupling to QCD: for simplicity, in the rest of the paper we will refer to it simply as g aγ γ , ignoring the model-dependent contribution g 0 aγ γ proportional to the electromagnetic anomaly of the U (1) P Q symmetry.
where m 2 a | =0 is given by the expression (2.11) with = 0, i.e., with m u = m d =m, and the third equation is the normalization condition. At the leading order in b, the following results are found: Let's now consider the realistic case = 0. If we write the squared-mass matrix (2.8) as we can evaluate the eigenvalues and the eigenstates of the matrix M 2 at the first order in the parameter , by treating the term δM 2 as a small perturbation. In particular, using firstorder perturbation theory, we obtain for the axion eigenstate: where |π 0 =0 = |π 3 and m 2 π 0 | =0 = 2Bm. Therefore, at the leading order in and b: (3.11) We observe that, if we take the formal limits m s → ∞ and A → ∞, this result correctly reduces to the corresponding expression derived with the Chiral Effective Lagrangian (χ E L) with L = 2 flavors at LO, i.e., (3.12)

Hadronic decays with the axion
This section is devoted to the study of the following processes: Among all the possible hadronic decays involving also the axion, these are the ones involving the lowest-energy hadrons. (Since, as we shall see below, every axion in the final or initial state implies a factor b in the decay amplitude, multiaxion processes are extremely suppressed and we disregard them.) The couplings of the axion with hadrons in general (and with the lightest mesons in particular) have already been investigated in the past literature, in many cases using also chiral effective Lagrangian techniques (see, e.g., Refs. [41][42][43]), but never using the WDV Lagrangian (2.7). Moreover, the particular processes (4.1) have never been explicitly investigated before. 3 From an experimental point of view, there are well-known bounds on the decay widths of η/η → π 0 π 0 and η/η → π + π − : they will be compared to our predictions in Sect. 7.
In order to compute the amplitudes of the processes (4.1), we must derive the interaction vertices between the axion and the pseudoscalar mesons. This can be achieved by expanding the WDV Lagrangian (2.7) up to the fourth order in the fields. We thus obtain the following quartic Lagrangian: where f abc are the SU (3) structure constants (defined as [λ a , λ b ] = 2i f abc λ c ). In particular, only the following term of the quartic Lagrangian is relevant for studying the processes (4.1): where, as usual, π ± = π 1 ∓iπ 2 √ 2 are the charged pion fields. As in the previous section, we shall work at the leading order in the parameter b. Moreover, considering also the explorative nature of this study, we shall neglect (for simplicity) isospin violations ( = 0). 4 With these approximations, the following (relevant) mixing parameters are found diagonalizing the squared-mass matrix (2.8): where ϕ is the mixing angle between π 8 and S, given by [45]: . (4.5) In particular, being = 0, π 3 can be simply identified with π 0 . Making use of Eq. (2.10) and of the expressions (4.4), the following quartic interaction terms are found from Eq.

Effects of an extra U(1) axial condensate
In this section, we will try to answer the following question: considering the relevance of the U (1) axial symmetry in defining the physical aspects of a hypothetical axion (i.e., its mass and its interactions), could a (no matter how small) hypothetical U (1) axial condensate significantly modify these expectations? More precisely, we will study the effects of a possible U (1) axial condensate on the various quantities that we have evaluated in the previous sections, by using a chiral effective Lagrangian model proposed in Refs. [30][31][32] and then elaborated on in Refs. [33][34][35][36][37]: it can be interpreted as an extension of the WDV model with the inclusion of a U (1) axial condensate and therefore we will refer to it as the "extended model". In this model the is a color singlet and (ii) C U (1) is a genuine 2L-quark condensate, i.e., it has no disconnected part proportional to some power of the quark-antiquark chiral condensate qq . The explicit form of this condensate has been discussed in detail in Refs. [35,36]. In what follows we shall consider the case L = 3. The Lagrangian of the extended model is thus written in terms of the topological charge density Q, the usual mesonic field U i j ∼q j R q i L , and a new field variable X ∼ det stqs R q t L , associated with the U (1) axial condensate, which under a general SU 2)] transforms as: In the usual nonlinear parametrization, the field X can be written as where F X is essentially the vacuum expectation value of X ( X = F X √ 2 ), i.e., the U (1) axial condensate, and S X is an exotic flavor-singlet pseudoscalar field. The model can be "axionized" in the same way as the WDV model. The Lagrangian of the "axionized extended model" is written as: Integrating out the auxiliary field Q, one obtains: The model is characterized, with respect to the WDV model, by three new parameters: ω 1 , κ 1 , and F X . As already observed in Refs. [35,36,39], the Lagrangian of the extended model reduces to that of the WDV model by first choosing ω 1 = 1 and then letting F X → 0. Therefore, ω 1 = 1 seems to be the most "natural" choice, at least at low temperatures, near T = 0, where minimal deviations from the results of the WDV model are expected (on the other side, ω 1 must necessarily vanish above the chiral transition temperature in order to avoid a singular behaviour of the anomalous term: see Refs. [30][31][32]37]). Expanding the Lagrangian up to the second order in the fields, one finds the following squared-mass matrix for the fields π 3 , π 8 , S, S X , S a : The eigenstates of this matrix are the usual pseudoscalar mesons π 0 , η and η , plus another exotic pseudoscalar state, called η X , and the axion. At the leading order in b, the following value for the squared mass of the axion is found: Also in this case [see the discussion after Eq. (2.11)], the expression for m 2 a turns out to be in agreement with the relation m 2 a = 2b 2 F 2 π χ QC D , considering the expression of χ QC D which is found using the extended model (see Ref. [39]). Moreover, we notice that for ω 1 = 1 the mass of the axion in the extended model is smaller than the one obtained in the WDV model, due to the positive corrective factor in the denominator. If, instead, we consider ω 1 = 1, the result coincides precisely with the result (2.11) of the WDV model, independently of the other parameters (κ 1 and F X ) of the extended model. This is not a totally unexpected result since in this particular case the potential coincides with that of the WDV model, apart from a term independent of the axion field (a more detailed explanation of this can be found in Ref. [39]).

Numerical results
In this section we report numerical estimates for the axionphoton-photon coupling constant and for the decay widths of the hadronic processes (4.1) with the axion. For the numerical computations, we have used the following values of the known parameters: • F π = (92.1 ± 1.2) MeV (see Ref. [46], where the value of f π ≡ √ 2F π is reported). • A = (180 ± 5 MeV) 4 (see Ref. [6] and references therein). • For what concerns the parameter B and the quark masses m u , m d , m s , we can make use of the well-known relations (see, e.g., Ref. [1]) between Bm u , Bm d , Bm s and the squared pseudoscalar-meson masses, derived using leading-order chiral perturbation theory (and ignoring small corrections due to the mixing with the axion): (6.1) The pseudoscalar-mesons masses are given by [46] ⎧ We also need m η = 547.862 (17) MeV and m η = 957.78(6) MeV.

Axion-photon-photon coupling constant g aγ γ
In Table 1 we report the numerical estimate for the axionphoton-photon coupling constant g aγ γ , obtained using the expression (3.11) that we have derived in Sect. 3 using the "axionized" WDV model (as we have seen in Sect. 5, this expression is not modified using, in place of the WDV model, a "natural" extension of it which also includes an extra U(1) axial condensate): for comparison, we also report the corresponding estimates derived using the Chiral Effective Lagrangian (χ E L) with L = 2 light quark flavors at LO [O( p 2 )] and NLO [O( p 4 )] (see Ref. [21] and references therein). Table 1 Numerical values of the axion-photon-photon coupling constant g aγ γ , obtained using Eq. χ E L (L = 2) at LO [21] 3 .59 ± 0.05 χ E L (L = 2) at NLO [21] 3 .42 ± 0.07 3.29 ± 0.06

Hadronic decay widths with the axion
The decay widths for the processes (4.1) are given by where the amplitudes g ηππa and g η ππa are given by Eqs.
(4.7) and (4.8) respectively and (3) (M|m 1 , m 2 , m 3 ) is the phase space (with the usual "relativistic" normalization) for three particles of masses m 1 , m 2 , m 3 with total energy M in the center-of-mass system. The exact expression is rather complicated (see Eq. (3.18) in Refs. [35,36], and also Ref. [47] and references therein), but it is surely a good approximation to take m a 0, considering the experimental upper bound on the axion mass m a 10 −2 eV [18,20]. The expression for the phase space for two particles of mass m and one massless particle turns out to be Inserting the numerical values of the constants (6.1)-(6.2) and F π , we obtain the following results: In particular, in Sect. 3 we have computed the axionphoton-photon coupling constant g QC D aγ γ and the result is given by the expression (3.11), that we have compared with the one obtained using the Chiral Effective Lagrangian with L = 2 light quark flavors. As we have verified in Sect. 5, this expression (as well as the expression for the mass of the axion) is not modified using, in place of the WDV model, a "natural" extension of it which also includes an extra U(1) axial condensate.
In Table 1 of Sect. 6 we have reported the numerical estimate for the axion-photon-photon coupling constant (3.11): for comparison, we have also reported the corresponding estimates derived using the Chiral Effective Lagrangian with L = 2 light quark flavors at LO and NLO. Comparing our result with the estimate found using the L = 2 Chiral Effective Lagrangian at LO, we get a value which is about 9% smaller, and it is also a bit smaller than (but almost compatible within the errors with) the value obtained using the L = 2 Chiral Effective Lagrangian at NLO. Of course, in the hypothesis that this type of process will be observed in the future, it will be important to know the level of accuracy of a given theoretical estimate, when comparing it with the experimental result, and in this perspective our result will be surely relevant. Looking at the values reported in Table 1, one could optimistically consider our result as a more "precise" determination, with respect to the result obtained using the L = 2 Chiral Effective Lagrangian at LO and NLO (because, maybe, our effective model, already at tree level, is able to reproduce results with an accuracy comparable to the one which is obtained using the L = 2 Chiral Effective Lagrangian after including many higher-order correc-tions…). Adopting, instead, a more conservative approach, one could simply consider our result as an alternative determination using a chiral effective Lagrangian model, which (when compared with other similar determinations) allows to estimate a sort of "systematic uncertainty" for this kind of theoretical predictions.
Then, in Sect. 4 we have performed an explorative study of the hadronic decays η/η → ππa (which, among all the possible hadronic decays involving also the axion, are the ones involving the lowest-energy hadrons): the expressions for the amplitudes g ηππa and g η ππa are given by Eqs. (4.7) and (4.8) respectively. In Sect. 6, Eq. (6.5), we have reported the numerical estimates for the corresponding decay widths: these are the main original results obtained in this paper. (Moreover, as we have found in Sect. 5, the addition of a possible U(1) axial condensate, while not modifying the η → ππa decays, makes the η → ππa decay widths a bit smaller by a factor 0.88 (F π /F η ) 2 ≤ 1.) Considering the existing experimental bounds on b (based on astrophysical and cosmological considerations) [18,19], 10 −18 b 10 −10 , the decay widths (6.5) turn out to very small, smaller than about 10 −22 MeV. 5 As far as we know, no experimental search for these processes has been attempted up to know. However, even if the electromagnetic decay of the axion (a → γ γ ) certainly remains the most promising process which might provide some experimental signature of the axion, we believe that it would be worthwhile to look also for these possible decay processes η/η → ππa in future η-factory experiments.
We conclude by observing that our estimates (6.5) for the widths of the η/η → ππa decays are based on the expressions (4.7) and (4.8) for the amplitudes g ηππa and g η ππa , which have been obtained directly (i.e., at leading order) from our chiral effective Lagrangian model (described in Sect. 2). It is plausible that these LO estimates will receive large contributions from chiral loop corrections at NLO and NNLO, and from strong final-state rescattering (as it happens, for example, in the η/η → 3π decays). Alternatively, one could consider the approach described in Refs. [44,50], in which one takes into account an extended chiral effective Lagrangian model, which also includes the lowestlying nonet of scalar mesons (i.e., a linearized version of the [nonlinear] chiral effective Lagrangian model described in Sect. 2): in this alternative approach, also contributions to the η/η → ππa amplitudes coming from scalar-meson exchanges are taken into account.
We believe that it would be worthwhile to go beyond the explorative study undertaken in this paper and to better investigate the η/η → ππa decays following the "guidelines" mentioned above: some progress in these directions is expected in the near future.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: Data sharing not applicable to this article as no datasets were generated or analysed during the current study.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copy- 5 Just for comparison, we recall here the experimental bounds on the CP-violating η/η decays in two pions [46]: ex p (η → π 0 π 0 ) < 4.6 × 10 −7 MeV, ex p (η → π + π − ) < 1.7 × 10 −8 MeV, ex p (η → π 0 π 0 ) < 7.8×10 −5 MeV, ex p (η → π + π − ) < 3.5×10 −6 MeV. We also observe that, even with the largest value of b allowed by the abovementioned astrophysical bounds, i.e, b 10 −10 , the decay widths (6.5) turn out to be about a factor 10 −3 (for η → ππa) and 10 −2 (for η → ππa) smaller than the model-independent bounds on the rates of the rare (CP-violating) decays η(η ) → ππ, which have been derived in Refs. [48,49], using the experimental limits on the neutron electric dipole moment. right holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .