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$ 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$ light quark flavors. Moreover, we study the decay channels $\eta/\eta'\rightarrow \pi\pi a$ and we estimate the corresponding decay widths, using the existing bounds on the $U(1)_{PQ}$ breaking scale.


Introduction
It is well known that the QCD Lagrangian L QCD 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], d N ≃ m 2 π m 3 N e|θ| ≃ 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 neutron electric dipole moment and the θ angle and a more detailed discussion can be found in Refs. [3,4,5]; see also Ref. [6] 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) in 1977 [7] and developed by Weinberg and Wilczek in 1978 [8,9]. The key idea (see also Ref. [10] 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 QCD with a dynamical degree of freedom, namely the combination θ − a P Q Sa fa : on the vacuum we get θ − a P Q Sa fa = 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), is also anomalous, with the related U(1) axial current J µ 5 =qγ µ γ 5 q satisfying the relation ∂ µ J µ 5 = 2LQ. 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 [7,8,9] 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. [11]). In order to bypass these experimental bounds, the so-called "invisible axion" models were developed, such as the Kim-Shifman-Vainshtein-Zakharov (KSVZ) model [12,13] and the Dine-Fischler-Srednicki-Zhitnisky (DFSZ) model [14,15], 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, Ref. [16]): 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. [17] 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 axion-photon-photon coupling constant. This last, in general, is the sum of two contributions, g aγγ = g 0 aγγ + g QCD 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 [12,13]), while g QCD 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 QCD 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. [18] and references therein).
The aim of this paper is to compute the axion-photon-photon coupling constant g QCD 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, Di Vecchia, Veneziano, et al. [19,20,21], 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 Sec. 2, for the benefit of the reader, we briefly recall this chiral effective Lagrangian model, as well as its "axionized" version (see Ref. [22]). Using this model, in Sec. 3 we compute the axion-photon-photon coupling constant g QCD aγγ and the result is compared with the one obtained using the Chiral Effective Lagrangian with L = 2 light quark flavors.
Sec. 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 possibile hadronic decays involving also the axion, these are the ones involving the lowest-energy hadrons.
In Sec. 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 Ref. [23] and then elaborated on in Refs. [24,25,26], which can be interpreted as an extension of the model considered in Sec. 2 with the inclusion of a U(1) axial condensate.
In Sec. 6 we report numerical estimates for the axion-photon-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 Sec. 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.  [19,20,21] 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. [20] for a detailed discussion) The mesonic field U is represented by a 3 × 3 complex matrix, which can be written in terms of the quark fields as U ij ∼q jR q iL , up to a multiplicative constant. * Under a general 2) * 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 fieldstrength tensorF µν .
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 ).
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 so correctly reproducing the U(1) axial anomaly of the fundamental theory.
The model can be easily "axionized" (see Ref. [22]), 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 parametrizes 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): 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 θ ij is an orthogonal mixing matrix.
From the astrophysical bounds on the scale f a [16] (or better on f a /a P Q , but a P Q ∼ O(1) for the more realistic axion models [27]) 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 [22]: This expression is in perfect agreement with the well-known relationship (valid at the leading order in b) [13]: m 2 a = 2b 2 F 2 π χ QCD , between the squared mass of the axion and the topological susceptibility of QCD, χ QCD ≡ −i d 4 x T Q(x)Q(0) | QCD , considering the expression of χ QCD which is found using the WDV model (see Refs. [22,28] 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 A µ is the electromagnetic field and Q = diag(2/3, −1/3, −1/3) is the quark electric-charge matrix (in units of e, the absolute value of the electron charge), and (ii) by adding the following term, which reproduces the electromagnetic anomaly of the U(1) and SU(3) axial currents (see Ref. [29]): where G = e 2 N C 32π 2 ε µνρσ F µν F ρσ , F µν being the electromagnetic field-strenght 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 immediatly sees that this term contains an axion-photonphoton interaction of the type (1.4), with the following expression for the axion-photonphoton coupling constant: * g aγγ = α e.m. πF π θ π 3 Sa + 1 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: (3.4) * As we have already said in the Introduction, this is indeed the model-independent contribution g QCD 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.
the axion eigenvector is given by (3.5) 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 Sa | ∆=0 = 0. The eigenvector equations are, in this case: 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 first-order 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: while the corrections to the other mixing parameters are of order O(∆ 2 ). Finally, substituting the expressions (3.7) and (3.10) into Eq. (3.3), we find the following result: 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 (χEL) with L = 2 flavors at LO, i.e.,

Hadronic decays with the axion
This section is devoted to the study of the following processes:          η → π 0 + π 0 + a, η → π + + π − + a, η ′ → π 0 + π 0 + a, η ′ → π + + π − + a. 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. [30]), but never using the WDV Lagrangian (2.7). Moreover, the particular processes (4.1) have never been explicitly investigated before. * From an experimental point of view, * However, in the recent Ref. [31] similar processes, such as a → 3π or a → η(η ′ )ππ, involving QCDscale axionlike particles with masses m π m a 3 GeV, have been investigated, using also (for the case m a 1 GeV) chiral effective Lagrangian techniques.
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 Ref. [23] and then elaborated on in Refs. [24,25,26]: 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 U(1) axial anomaly is implemented as in the WDV model (by properly introducing the auxiliary field Q), so that it correctly satisfies the transformation property (2.4) under the chiral group, but it also includes an extra U(1) axial condensate, in addition to the usual chiral condensate qq . This U(1) axial condensate has the form C U (1) = O U (1) , where, for a theory with L light quark flavors, O U (1) is a 2L-quark local operator that has the chiral transformation properties of O U (1) ∼ det st (q sR q tL ) + det st (q sL q tR ), where s, t = 1, . . . , L are flavor indices. The color indices (not explicitly indicated) are arranged in such a way that (i) O U (1) 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 Ref. [25]. 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 ij ∼q jR q iL , and a new field variable X ∼ det stqsR q tL , associated with the U(1) axial condensate, which under a general SU(3) L ⊗SU(3) R ⊗U(1) A chiral transformation [see Eq. (2.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. [25,28], 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. [23,26]).
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 π χ QCD , considering the expression of χ QCD which is found using the extended model (see Ref. [28]). 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. [28]).
Concerning the axion-photon-photon coupling g aγγ , the following result is found: Also in this case, setting the "natural" value ω 1 = 1 we recover the WDV expression (3.11), independently of the other parameters of the extended model.
. Anyway, if we neglect the isospin violations (∆ = 0), this term does not contribute to the processes (4.1). If, for the reasons explained above, we take ω 1 = 1, which is the "natural" choice (at least at T = 0), we easily derive (proceeding as in Sec. 4 and making use of the results already found in Refs. [24,25]) the following expressions for the coupling constants g ηππa and g η ′ ππa : and Am m s , can be interpreted as the η ′ decay constant (see Refs. [23,24,25]) andφ is the mixing angle between π 8 and S, which turns out to be a bit larger than the value ϕ in Eq. (4.5), being [24,25]: tanφ = F η ′ Fπ tan ϕ. We observe that in the limit F X → 0 we have F η ′ → F π andφ → ϕ, and the expressions (5.7) and (5.8) correctly reduce to the WDV results (4.7) and (4.8).

Numerical results
In this section we report numerical estimates for the axion-photon-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. [33], where the value of f π ≡ √ 2F π is reported).
• 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): 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 axion-photon-photon coupling constant g aγγ , obtained using the expression (3.11) that we have derived in Sec. 3 using the "axionized" WDV model (as we have seen in Sec. 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 (χEL) with L = 2 light quark flavors at LO [O(p 2 )] and NLO [O(p 4 )] (see Ref. [18] and references therein).

Conclusions: summary of the results and prospects
In this paper we have investigated the most interesting decay processes involving axions, photons and the lightest pseudoscalar mesons, making use of the "axionized" version of a chiral effective Lagrangian model proposed by Witten, Di Vecchia, Veneziano, et al. (WDV), 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 particular, in Sec. 3 we have computed the axion-photon-photon coupling constant g QCD 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 Sec. 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 Sec. 6 we have reported the numerical estimate for the axion-photonphoton 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 corrections. . . ). 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.