Abstract
We revisit the information on the two lightest \(a_0\) resonances and S-wave \(\pi \eta \) scattering that can be extracted from photon–photon scattering experiments. For this purpose we construct a model for the S-wave photon–photon amplitudes which satisfies analyticity properties, two-channel unitarity and obeys the soft photon as well as the soft pion constraints. The underlying I=1 hadronic T-matrix involves six phenomenological parameters and is able to account for two resonances below 1.5 GeV. We perform a combined fit of the \(\gamma \gamma \rightarrow \pi \eta \) and \(\gamma \gamma \rightarrow K_SK_S\) high statistics experimental data from the Belle collaboration. Minimisation of the \(\chi ^2\) is found to have two distinct solutions with approximately equal \(\chi ^2\). One of these exhibits a light and narrow excited \(a_0\) resonance analogous to the one found in the Belle analysis. This however requires a peculiar coincidence between the \(J=0\) and \(J=2\) resonance effects which is likely to be unphysical. In both solutions the \(a_0(980)\) resonance appears as a pole on the second Riemann sheet. The location of this pole in the physical solution is determined to be \(m-i\varGamma /2=1000.7^{+12.9}_{-0.7} -i\,36.6^{+12.7}_{-2.6}\) MeV. The solutions are also compared to experimental data in the kinematical region of the decay \(\eta \rightarrow \pi ^0\gamma \gamma \). In this region an isospin violating contribution associated with \({\pi ^+}{\pi ^-}\) rescattering must be added for which we provide a dispersive evaluation.
Similar content being viewed by others
1 Introduction
The \(a_0(980)\) and \(f_0(980)\) scalar mesons were the first observed members of a family of exotic resonances in QCD which are located very close to an inelastic two-particle (or quasi two-particle) threshold (see the review [1]). The \(a_0(980)\) resonance was discovered a long time ago and seen in both the \(K\bar{K}\) and the \(\pi \eta \) channels [2, 3] but its properties are still imprecisely known. This is partly because \(\pi \eta \) production experiments using an \(\eta \) beam, analogous to those which have allowed to determine the \(\pi \pi \) or \(\pi K\) phase shifts, are not feasible. The \(a_0(980)\) properties have to be determined solely from final-state rescattering effects. As shown by Flatté [4], who proposed a simple replacement for the Breit-Wigner resonance formula accounting for the \(\pi \eta -K\bar{K}\) coupled-channel dynamics, ambiguous results for the \(a_0\) width can be obtained. In the PDG [5], indeed, its width is simply quoted as lying in a range from 50 to 100 MeV. Beyond the value of the width, one would like to determine the positions of the poles on the unphysical Riemann sheets. For resonance states close to an inelastic threshold, there are three sheets (II, III, IV in the case of the \(a_0\)) which are physically relevant (i.e. a pole in one of these can be close to the physical region). It is also clear that a better determination of the resonance properties is closely tied to a better knowledge of the physical scattering T-matrix.
A theoretically motivated treatment of final-state interactions becomes prohibitively difficult for multiparticle final states. In this regard, photon–photon to meson-meson scattering amplitudes are very favourable processes.Footnote 1 They are free of initial-state interactions and satisfy dispersion relations which can be constrained, in the case of \(\pi \pi \) or \(\pi \eta \) by both soft photon [7, 8] and soft pion [9] low-energy theorems. As was illustrated in the seminal papers [10, 11] a predictive representation can be implemented at the level of the partial-waves with a simple modelling of the left-hand cut. From an experimental point of view \(\gamma \gamma \rightarrow \eta \pi ^0\) cross-sections were first measured by the Crystal Ball collaboration [12]. Measurements with much higher statistics were recently performed by the Belle collaboration [13]. There are also experimental results on the \(\eta \rightarrow \pi ^0\gamma \gamma \) decay amplitude [14, 15]. Recent experimental data exist also for \(\gamma \gamma \rightarrow {K^0}{\bar{K}^0}\) [16] and \(\gamma \gamma \rightarrow {K^+}{K^-}\) [17] which will be considered in our study.
There are two puzzling aspects in the data analysis performed by the Belle collaboration which we wish to reconsider: a) they find that the \(a_0(980)\) peak seems to be best described by an ordinary (i.e. essentially elastic) Breit-Wigner function and b) they find an excited \(a_0\), which could correspond to the \(a_0(1450)\), but has a width, \(\varGamma =65^{+2.1}_{-5.4}\) MeV, much smaller than the PDG average as well as a significantly smaller mass. These data have been re-analysed recently [18] based on a specific meson-meson T-matrix model [19] applied to the \(\pi \eta \) S-wave, from which the corresponding \(\gamma \gamma \) amplitude is deduced from a Muskhelishvili-Omnès (MO) construction [20, 21]. Including also the \(a_2(1320)\) resonance but no \(a_0(1450)\) resonance a good description of the data up to 1.1 GeV and a qualitative description up to 1.4 GeV has been achieved. The \(a_0(980)\) pole, in this model, lies on the fourth Riemann sheet. Interestingly, an analogous pole location was found in a lattice QCD calculation of the T-matrix [22] who implemented a coupled-channel generalisation of Lüscher’s single-channel method [23]. This result, however, corresponds to a pion mass \({m_{\pi }}=391\) MeV and it is not known how the pole location evolves upon varying \({m_{\pi }}\) (see, however, Ref. [24]). Fits to the Belle data with a conventional \(a_0(1450)\) were performed in Ref. [25] but only the integrated cross-sections were considered in that work. Detailed global descriptions of photon–photon scattering to two mesons were proposed also in Ref. [26] focusing mostly on the \(I=0,2\) channels. In the \(I=1\) sector, they considered the \(K\bar{K}\) channel but not \(\pi \eta \).
We perform here a global fit which takes into account both \(\pi \eta \) and \(K\bar{K}\) photon–photon data including all the differential cross-sections up to 1.4 GeV. In the \(I=1\) sector, we use the S-wave coupled-channel T-matrix model developed in Ref. [27] which satisfies unitarity, proper analyticity properties and matches to the chiral expansion up to the next-to-leading order at low energy. The S-wave photon–photon amplitudes are then deduced from a general MO representation involving two subtraction constants and implementing a simple description of the left-hand cut from cross-channel vector-meson exchanges. This is quite similar to Ref. [18], we differ mainly by using SU(2) chiral symmetry, which allows to fix one of the subtraction constants through a soft pion theorem.Footnote 2 The \(J=2\) partial-waves are described more phenomenologically as a sum of cross-channel resonance exchange and a direct \(a_2(1320)\) Breit-Wigner amplitude. The \(\gamma \gamma \rightarrow (K\bar{K})_{I=1}\) amplitudes are then combined with \(I=0\) amplitudes taken from a previous work [28] which considered \(\gamma \gamma \rightarrow \) \((\pi \pi )_{I=0,2}\), \((K\bar{K})_{I=0}\) in order to reconstruct the physical \({K^+}{K^-}\), \({K^0}{\bar{K}^0}\) amplitudes. A global fit of photon–photon to meson-meson data was performed some time ago based on unitarised chiral amplitudes [29], this model was also applied to the \(\eta \rightarrow {\pi ^0}\gamma \gamma \) decay [30, 31]. A remarkable qualitative agreement with the data available at the time was achieved using a single arbitrary parameter in the S-wave. Today, a much larger data set is available and the precision has increased significantly. We use here a model for the T-matrix involving six parameters which will be determined by performing fits to these data.
An interesting aspect of the \(\gamma \gamma \rightarrow \pi \eta \) amplitudes is that they can be probed experimentally both in the scattering regime: \(s_{\gamma \gamma } \ge (m_\eta +{m_{\pi }})^2\) and in the decay regime: \(0 \le s_{\gamma \gamma } \le (m_\eta -{m_{\pi }})^2 \). In the decay region the amplitudes are largely dominated by the light vector meson exchanges in the cross-channels: \(\gamma \pi \rightarrow \rho ,\omega \rightarrow \gamma \eta \), which were first computed in Ref. [32]. In this low-energy region the rescattering contributions can be estimated in the SU(3) chiral expansion [33]. They proceed essentially via \(\gamma \gamma \rightarrow {K^+}{K^-}\rightarrow {\pi ^0}\eta \) but, at low energy, the \({\pi ^+}{\pi ^-}\) contribution (which is isospin violating) is not negligible and must also be included [33]. We will present a dispersive calculation of this contribution which gives rise to a cusp in the energy distribution \(d\varGamma ^{\eta \rightarrow {\pi ^0}\gamma \gamma }/ds_{\gamma \gamma }\).
The plan of the paper is as follows. After recalling some general properties of photon–photon amplitudes and introducing the notation in Sect. 2, we write the unitarity relations for the partial-waves in Sect. 3 and present the modelling of the left-hand cut of these. In Sect. 4 we recall the derivation of a MO representation for the S-waves. A dispersive representation for the isospin violating S-wave, valid at low energy, is also derived. The comparison with the experimental data on photon–photon scattering is performed in Sect. 5 and the information that can be deduced on the \(a_0\) resonances are discussed.
2 Basic ingredients
2.1 Kinematics
We consider two-photon to two-meson scattering amplitudes
with \(M_1M_2= \pi \eta \) or \(K\bar{K}\). The Mandelstam variables are defined as usual by
The various physical regions in the s, \(t-u\) Mandelstam plane for the \(\gamma \gamma \rightarrow \pi \eta \), \(\gamma \pi \rightarrow \gamma \eta \) and \(\eta \rightarrow \gamma \gamma \pi \) processes are shown on Fig. 1. In the centre-of-mass frame of the two photons the momenta are expressed as
where
Taking \(\hat{z}\) to a unit vector along the z axis and \(\hat{v}\) to be a unit vector with polar angles \(\theta \), \(\phi \) such that \(\hat{z}\cdot \hat{v}=\cos \theta \), we can express t, u in terms of \(\cos \theta \) as,
The polarisation vectors of the two photons \(\epsilon _1(q_1,\lambda )\), \(\epsilon _2(q_2,\lambda ')\) in this frame read
and they satisfy \(\epsilon _1\cdot q_1 = \epsilon _2\cdot q_2 = \epsilon _1\cdot q_2= \epsilon _2\cdot q_1=0\ . \)
2.2 Amplitudes
We denote the \(\gamma \gamma \rightarrow \eta \pi \) helicity amplitude as \(L_{\lambda \lambda '}\)
(we have factored out the electric charge e and the dependence on the azimuthal angle \(\phi \)) while for the \(\gamma \gamma \rightarrow K\bar{K}\) amplitudes we use the same notation as in previous work: \(K^c_{\lambda \lambda '}(s,t)\) for charged kaons and \(K^n_{\lambda \lambda '}(s,t)\) for neutral kaons. Helicity amplitudes are convenient for performing the partial wave expansion [34]. They can be expressed in terms of tensor amplitudes using the reduction formulas, e.g.
By gauge invariance, the tensor amplitude \(W_{\mu \nu }\) must satisfy the two Ward identities
One can form two independent tensors which satisfy (9) which we take as
with
The tensor amplitude \(W^{\mu \nu }\) can then be expressed in terms of two scalar amplitudes A, B
which satisfy dispersion relations. Using Eqs. (5), (6) one can easily express the helicity amplitudes \(L_{\lambda \lambda '}\) in terms of the two scalar amplitudes,
Assuming unpolarised photon beams the differential cross-section reads,
Concerning the \(\eta \rightarrow \pi \gamma \gamma \) decay amplitude, the double differential distribution in the Dalitz plot reads,
and the distribution as a function of s only (which is the one available experimentally) is given by
2.3 Isospin
Using the following (usual) isospin assignments for the pions and the kaons
while \(\eta \sim {\vert 0,0 \rangle }\), the relations between the amplitudes \(\gamma \gamma \rightarrow {K^+}{K^-}\), \({K^0}{\bar{K}^0}\) and the isospin amplitudes \(\gamma \gamma \rightarrow (K \bar{K})_{I=0,1}\) read,
and the analogous relations between the amplitudes \(\gamma \gamma \rightarrow {\pi ^+}{\pi ^-}\), \({\pi ^0}{\pi ^0}\) and the corresponding isospin \(I=0,2\) amplitudes (which will also be needed) is
3 Partial waves: unitarity, analyticity
3.1 Right-hand cut and unitarity relations
It is convenient to collect the three \(I=1\) scattering amplitudes \(\pi \eta \rightarrow \pi \eta \), \(\pi \eta \rightarrow K\bar{K}\) and \(K\bar{K}\rightarrow K\bar{K}\) into a \(2\times 2\) matrix
and we can define the partial wave expansions as
where \(\theta \) is the scattering angle in the centre-of-mass system. The unitarity relation for the partial waves are easily derived and reads
with
Concerning the \(I=1\) photon–photon amplitudes we define the partial-wave expansion as
The unitarity relations for the S-waves, as will be implemented below, read
which also give the discontinuities of the partial-waves (extended to complex values of s) across the right-hand cut.
3.2 \(\pi \pi \) (isospin violating) contribution in S-wave unitarity
The unitarity relation written above (25) collects the contributions from the \(\pi \eta \) and the \(K\bar{K}\) states. Let us also consider the contribution from the \(\pi ^+\pi ^-\) state, which has the form
where \(d\varPhi \) is the phase-space integration measure. This contribution, being proportional to the \(\pi \pi \rightarrow \eta \pi \) amplitude, is isospin violating but it is enhanced at low energy due the large size of the \(\gamma \gamma \rightarrow \pi ^+\pi ^-\) amplitude [33]. The ChPT evaluation performed in Ref. [33] amounts to using the \(O(p^2)\) tree-level amplitudes for both \(\pi \pi \rightarrow \eta \pi \) and \(\gamma \gamma \rightarrow \pi ^+\pi ^-\) in Eq. (26). An evaluation which goes beyond the chiral expansion can be performed which we now discuss.
We consider only the S-wave contribution in Eq. (26) and a restricted kinematical region such that \(s,t, u < 1\) \(\hbox {GeV}^2\). In such a region, the \({\pi ^+}{\pi ^-}\rightarrow \eta {\pi ^0}\) amplitude can be approximated in terms of three one-variable functions \(M_0\), \(M_1\), \(M_2\) (see [36, 37]),
An isospin violating parameter \(\epsilon _L\) has been factorised which may be taken as [38]
The three \(M_I\) functions obey a set of coupled Khuri-Treiman integral equations, see Ref. [38] for a complete review of work on this subject. We will use here the evaluation of the \({K^0}-{K^+}\) QCD mass difference from Ref. [39] (updated in [38]) based on experimental data on \(\eta \rightarrow 3\pi \) decays: \(\left( {m}^2_{K^0}-{m}^2_{K^+}\right) _{QCD}=(6.24\pm 0.38)\times 10^{-3}\) \(\hbox {GeV}^2\), which gives
The amplitudes corresponding to a given \(\pi \pi \) isospin state \(I,I_z\)
are easily expressed using crossing symmetry and the Wigner-Eckart theorem. In the unitarity relation (26) the amplitudes with \(I=0,2\), \({{\mathcal {M}}}^{00}\) and \({{\mathcal {M}}}^{20}\), are needed which have the following expressions
in terms of the \(M_I\) functions. We denote the angular integrals of these amplitudes as
With this notation, the \(\pi \pi \) contribution to the unitarity relation is finally expressed as follows,
where
and \(h^I_{0,++}(s)\) are the two S-wave \(\gamma \gamma \rightarrow (\pi \pi )^I\) amplitudes with \(I=0,2\). This \(\pi \pi \) discontinuity of \({l}_{0++}\) can be estimated from Eq. (33) using inputs from Ref. [28] for \(\gamma \gamma \rightarrow \pi \pi \) and from [27] for \(\pi \pi \rightarrow \eta \pi \). The result of this estimate is illustrated on Fig. 2 and compared with the chiral calculation at NLO. The dispersive evaluation displays a square-root singularity at \(s=(m_\eta -{m_{\pi }})^2\) induced by the endpoint of the left-hand cut in the functions \(\hat{M}_0\), \(\hat{M}_2\) which overlaps with the right-hand cut as a result of the instability of the \(\eta \). As a further consequence, the phase of the partial-waves \(M_I+\hat{M}_I\) violate Watson’s theorem and do not cancel with the phases of \(h_{0,++}^I\) such that the discontinuity has both a real and an imaginary part.
3.3 Left-hand cut: born amplitudes
A left hand cut in the \(\gamma \gamma \) partial-waves is generated by singularities in the cross-channels \(\gamma P_1 \rightarrow \gamma P_2\). In the case when \(P_1=P_2={K^+}\) or \({\pi ^+}\) the leading singularity is the kaon or pion pole and the corresponding (so-called Born) amplitudes read
with \(P={K^+},\ {\pi ^+}\). We will need the \(I=1\) component of the \(K^+\) Born amplitude projected on the S-wave which reads,
with \(\beta _P(s)=\sqrt{1-4m^2_P/s}\). We also recall the expressions of the \(J=2\) Born helicity amplitudes
3.4 Left-hand cut: vector meson exchanges
Leading contributions to the left-hand cut in the \({\pi ^0}\eta \) amplitudes are modelled from the \(\rho \), \(\omega \), \(\phi \) vector meson exchanges. We define the \(VP\gamma \) coupling constants \(G_{VP}\) through simple Lagrangians
from which one can relate their values to the decay widths of the vector mesons via
The vector-exchange contributions to the \(\gamma \gamma \rightarrow P_1 P_2\) amplitudes are easily computed,
They give rise to poles in the zero-width approximation, which is adequate in the kinematical regions of interest here. The helicity amplitudes are
and the corresponding partial-waves with \(J=0,2\) have the following form
with
which is the cosine of the scattering angle when \(t=m_V^2\). The function \(L_V(s)\) is given by the angular integral
with
it can be expressed in terms of \(X_V\) as
We note that the partial-wave \(l^V_{0++}(s)\) has a soft pion Adler zero at \(s=s_A\) which can be approximated as
From the integral representation (44) one sees that the function \(L_V\) has endpoint singularities at \(z=\pm 1\) when \(s=s_V\), which thus corresponds to a branch point of \(L_V(s)\). Another endpoint singularity occurs when \(s=\infty \). When \(s > s_V\), the denominator remains strictly positive. Therefore, \(L_V\) is an analytic function of s with a cut on the negative real axis: \( -\infty< s < s_V\). The discontinuity of \(L_V\) along the cut is easily determined
from which one deduces the left-cut discontinuities of the vector-exchange partial-waves
We will use these discontinuities in the Muskhelishvili-Omnès representations below. We must consider also the vector exchange contributions for the \(\gamma \gamma \rightarrow (K\bar{K})_{I=1}\) amplitudes, we denote the relevant combination of coupling constants as
The imaginary parts of the \(J=0,2\) partial-wave amplitudes along the left-hand cut read
The updated values of the couplings \(C_{VP}\) are collected in Table 1 below. (a)
4 Representations of the \(J=0,2\) partial-waves
4.1 Muskhelishvili-Omnès representations for the S-waves
In order to write a dispersive representation for \(l_{0++}\) some knowledge concerning its asymptotic behaviour is needed. Let us then consider the angular integral,
in the \(s\rightarrow \infty \) limit. There are two regions of the angular variable z for which the behaviour of the integrand is known: (a) when z is close to 0, then \(t\sim u\sim -s\sim -\infty \). In this regime, \(L_{++}(s,t)\) can be estimated from QCD-based methods [40, 41] according to which one has
(b) When z is close to 1, then \(|t|<< s\) which is the region where Regge theory applies, the leading contribution from the vector meson trajectory gives
with \(\alpha _V\simeq 0.5\), \(\alpha '\simeq 0.9\) \(\hbox {GeV}^{-1}\) and \(\beta _V(t)\) is a smooth function when \(t<0\). Assuming only that the integrand evolves smoothly between these two regimes when \(0\le z\le 1\), one deduces that the \(l_{0++}(s)\) should not grow faster than \(\sqrt{s}\) when \(s\rightarrow \infty \).
Furthermore, the \(J=0\) partial-waves obey soft-photon theorems [7, 8] which imply that the ratios \(l_{0++}(s)/s\) and \((k^1_{0++}(s)-k^{1,Born}_{0++}(s))/s\) remain finite when \(s\rightarrow 0\). Therefore, they can be expressed as unsubtracted dispersion relations in terms of the left-hand and right-hand cuts discontinuities,
with \(m_+=m_\pi +m_\eta \). The isospin symmetry limit has been assumed here, implying that the \(\eta \) meson is stable and the discontinuities are real. At low energies the left-cut discontinuities are dominated by the vector meson exchanges. We will introduce subtractions, in order to reduce the influence of the higher energy regions where the discontinuities are not known. The right-cut discontinuities are given by the unitarity relations. Unitarity is known to be saturated with two channels to a very good approximation in the region of the \(a_0(980)\) resonance. We will assume that elastic unitarity holds below the \(K\bar{K}\) threshold and that two-channel unitarity remains a sufficiently good approximation up to \(\sqrt{s}=1.4\) GeV. Under the assumption of two-channel unitarity the dispersion relations (55) form a set of coupled inhomogeneous Muskhelishvili equations. They can be solved in terms of a two-channel MO matrix which satisfies a homogeneous set of coupled equations in terms of the T matrixFootnote 3
Asymptotic conditions on the phase-shifts are imposed which ensure that the set of equations (56) has a unique solution once initial conditions at \(s=0\) are given (see Sect. 3.2 in Ref. [27] and references therein). At \(s=0\) one can take
Multiplying the amplitudes \((l_{0++},k^1_{0++})\) by the inverse of the matrix \(\varvec{\varOmega }_0\) removes the right-hand cuts. We can use this property in writing once-subtracted dispersion relations for the two functions
This provides MO-type dispersive representations for the \(\gamma \gamma \) amplitudes \(l_{0++}(s)\), \(k^1_{0++}(s)\) in terms of their imaginary parts on the left-cut and two parameters, \(b_l\), \(b_k\)
The functions \(L_i(s)\) are dispersive integrals over the left-hand cuts. We express them in a way which allows to easily implement the presence of an Adler zero at \(s=s_A\) in the amplitude \(l_{0++}\) (see “Appendix B”)
where the functions \(D_{ij}(s)\) are the matrix elements of the inverse of the Omnès matrix,
The functions \(R_i(s)\), secondly, are the dispersive integrals over the right-hand cut,
A relation between the parameters \(b_l\) and \(b_k\) can be derived from imposing that the amplitude \(l_{0++}\) has an Adler zero at \(s=s_A\),
The parameter \(b_k\) will eventually be fitted to the experimental data but we can estimate its order of magnitude by matching the amplitude \(k^1_{0++}\) with the SU(3) chiral expansion. Including the order \(p^4\) contributions (see “Appendix B”) and an estimate of the order \(p^6\) from the vector-meson exchange amplitudes, we obtain
using the determination \(L_9+L_{10}= (1.44\pm 0.08)\cdot 10^{-3}\) taken from Ref. [42]. Below, the value of \(b_k\) will be fitted to the experimental data, the resulting combination \(b_k +L_2(0)+R_2(0)\) turns out to have a sign compatible with (64) and a magnitude smaller by a factor of two. The result of the dispersive construction of the S-wave amplitude \(l_{0++}\) is illustrated in Fig. 3. The corresponding result for the \(K\bar{K}\) amplitude \(k^1_{0++}\) is shown in “Appendix C” (see Fig. 13).
4.2 Dispersive construction of the isospin-violating S-wave
In Sect. 3.2 we have considered the unitarity contribution to the \(\gamma \gamma \rightarrow \pi \eta \) S-wave amplitude induced by the \({\pi ^+}{\pi ^-}\) intermediate state, which is isospin violating. Let us call \(\tilde{L}_{++}\) the isospin-violating part of the \(\gamma \gamma \rightarrow \pi \eta \) helicity amplitude and \(\tilde{l}_{0++}\) the corresponding S-wave. \(\tilde{L}_{++}\) can be defined as a matrix element,
We will attempt here to estimate the amplitude \(\tilde{l}_{0++}\) at low energy only and we write a dispersive representation keeping only the contribution from the \(\pi \pi \) cut,
We have used two subtractions in Eq. (66) in order to strongly reduce the influence of the energy region above 1 GeV. One subtraction constant is fixed from imposing the soft photon zero. We can then estimate the parameter \(\tilde{\lambda }\) in Eq. (66) by using the soft pion limit. This limit provides a relation between the amplitude \(\tilde{l}_{0++}(s=s_A)\) and a matrix element of the pseudo-scalar operator \(p_0=i(\bar{u}\gamma ^5u+\bar{d}\gamma ^5d)\) which, in turn, can be estimated using ChPT. This is detailed in “Appendix B.3”. Using Eqs. (B.22), (B.23) from this appendix we obtain the following result for \(\tilde{\lambda }\)
where the loop functions \(G_\pi \), \(G_K\) are given in Eq. (B.1). The discontinuity \(\mathrm{disc}\,[{l}_{0++}(s)]_{\pi \pi }\) (given in Eq. (33)) has a singularity at the pseudo-threshold \(s=(m_\eta -{m_{\pi }})^2\) induced by the \(\pi \pi \rightarrow \eta \pi \) partial-wave (see Fig. 2). The integral in Eq. (66), however, is finite. It is defined by using the \(m_\eta ^2+i\epsilon \) limiting prescription exactly in the same way as those which appear in the Khuri-Treiman equations for the \(\eta \rightarrow 3\pi \) amplitude (see e.g. [36, 43]).
The result, in the low energy region relevant for \(\eta \rightarrow \pi \gamma \gamma \) is shown in Fig. 4 and compared to the corresponding chiral \(O(p^4)\) result from Eq. (B.6). The chiral and dispersive real parts agree at \(s=m_\eta ^2\), as a result of the soft pion relation, but the two amplitudes differ substantially at lower energy:
The cusp at the \(\pi \pi \) threshold is much more pronounced in the dispersive amplitude which is approximately five times larger in magnitude than the \(p^4\) amplitude at \(s=4m_\pi ^2\).
The \(p^4\) amplitude has a zero at \(s=4/3m_\pi ^2\) in contrast to the dispersive amplitude which has no zero in this region. This is because this zero is unrelated to a soft photon constraint and is an accidental feature of the \(p^4\) amplitude.
4.3 D-wave amplitudes modelling
For the \(J=2\) partial-waves one can write unsubtracted dispersion representations analogous to those for \(J=0\), e.g.,
displaying the kinematical zeros at \(s=m_\pm ^2\) and \(s=0\). In the case of \(J=2\), however, it seems difficult to derive useful constraints from unitarity as in the case of \(J=0\) because in the important energy region of the \(a_2(1320)\) resonance there are too many contributing channels. We will therefore content with a simple Breit–Wigner estimate of the right-hand cut integral in Eq. (68).
We parametrise the coupling of the \(a_2\) resonance to the \(\eta \pi \) channel by a constant \(C^{a_2}_{\eta \pi }\) defined from the Lagrangian
The couplings to the \(\gamma \gamma \) channel involve two constants
The first term in (70) contributes to the \(+-\) helicity state, which is expected to be dominant, and the second term to the \(++\) helicity state. The expressions for the decay widths read,
where \(q_{\eta \pi }(s)=\sqrt{\lambda _{\eta \pi }(s)/4s}\). The experimental values of the branching fractions of the main hadronic decay modes are [5]
and the \(2\gamma \) width given by the PDG is
From this, one obtains the following values for the coupling constants
choosing \(C^{a_2}_{\eta \pi }\) to have a positive sign. Upon performing fits to the differential cross-sections (see below) the two couplings \(C^{a_2}_{\gamma \gamma }\), \(D^{a_2}_{\gamma \gamma }\) will get separately determined. Defining a coupling constant \(C^{a_2}_{K\bar{K}}\) in the same way as \(C^{a_2}_{\eta \pi }\) we get,
where the negative sign derives from flavour symmetry. The Breit–Wigner model for the \(a_2(1320)\) contributions is then taken as
This form is obtained by first computing the amplitudes from the Lagrangians (69) and (70) and then modifying them by introducing a width function \(\varGamma _T(s)\) in the denominator, for which we use the same modellingFootnote 4 as in Ref. [13]
and a Blatt-Weisskopf related function (with \(W_2(x)=9+3x^2+x^4\)). Finally, the \(J=2\) amplitudes \(l_{2\lambda \lambda '}\) are approximated by adding the Breit-Wigner \(a_2(1320)\) contribution in the s-channel to the vector-exchange contributions in the t, u channels,
The corresponding \(J=2\) \((K\bar{K})_{I=1}\) amplitudes are similarly described by a sum of three terms,
where the Breit–Wigner amplitudes are given by
It will be necessary to consider also the \((K\bar{K})_{I=0}\) amplitudes in order to be able to construct the \({K^+}{K^-}\) and \({K^0}{\bar{K}^0}\) amplitudes separately and compare with the experimental results. From the value of the \(K\bar{K}\) branching fraction of the \(f_2(1270)\) resonance [5]:
one derives the following values for the corresponding coupling constants
Based on nonet symmetry we have taken \(C^{f_2}_{K\bar{K}}\) and \(C^{a_2}_{K\bar{K}}\) to have the same sign.
5 Comparison with experiment
5.1 Experimental inputs
We will compare our model for the \(\gamma \gamma \) amplitudes with precise experimental data on \(\gamma \gamma \rightarrow \pi \eta \) from the Belle collaboration [13], as was done recently in Ref. [18]. In addition, we consider also here \(\gamma \gamma \rightarrow K\bar{K}\) data in order to provide further constraints on the coupled-channel dynamics which is believed to be important for the \(a_0(980)\) resonance. Recently, high statistics experimental data have been obtained by the Belle collaboration for the \(K_S K_S\) channel [16]. Experimental data for the charged kaons channel \(K^+K^-\) are also available in this low energy range [44] but they are older and have much less statistics. We will restrict ourselves to the energy range \(E\le 1.4\) GeV: we can use 448 differential cross-section points for \(\pi \eta \) (\(0.85 \le E\le 1.39\) GeV) and 240 differential cross-section points for \(K_SK_S\) (\(1.105\le E\le 1.395\) GeV).
5.2 Parameters of the T-matrix
Concerning the S-wave, firstly, we employ the T-matrix model of Ref. [27] which involves six parameters. It uses a chiral K-matrix type representation, which ensures one-channel unitarity below the \(K\bar{K}\) threshold and two-channel unitarity above, together with the chiral expansion. This kind of approach was initiated in [45], see Ref. [46] for a review. The T-matrix is written as
where the matrix \(\varvec{\varPhi }\) reads
it involves four phenomenological polynomial parameters \(\alpha _i\), \(\beta _i\) and the one-loop functions \(\bar{J}_{P_1P_2}\) are given in (D.8). The K-matrix has the following form
in which the subscript refers to the chiral order. The first two terms are to be computed from the chiral expansion of the scattering amplitudes involving the \(\eta \pi \) and \(K\bar{K}\) channels at order \(p^2\) and \(p^4\) respectively. The last term in Eq. (85) allows for a pole in s and involves two phenomenological parameters \(m_8\) and \(\lambda \),
The form of \(g_1\), \(g_2\) is derived from a resonance chiral Lagrangian [47]
such that \(\varvec{K}_{(6)}\) has chiral order \(p^6\) provided \(\lambda \) is O(1). In this model, the T-matrix has good analyticity properties and coincides with the chiral expansion at low energy up to \(O(p^4)\) provided that the parameters \(\alpha _i\), \(\beta _i\), \(\lambda \) are of chiral order O(1). In addition to these six phenomenological parameters, the T-matrix depends on the values of the \(O(p^4)\) chiral parameters \(L_i\) [48] and on the ratio \(c'_m/c'_d\). We will use here the set of \(L_i\) values from the \(p^6\) fit of Ref. [49] (labelled as BE14 in that reference). The ratio \(c'_m/c'_d\) is expected to be of order \(1-2\), we will use \(c'_m/c'_d=2\) as central value and include the variation as a source of error.
Obviously, such a model which implements two-channel unitarity is mostly justified in the \(a_0(980)\) region and below. We will assume that it remains qualitatively acceptable up to \(E\simeq 1.4\) GeV. The T-matrix is computed from Eq. (83) for \(E \le E_1\), \(E_1=1.5\) GeV. In the higher energy region \(E>E_1\), the T-matrix is described through a simple interpolation of the phase-shifts and the inelasticity such that \(\delta _{11}(\infty )=2\pi \), \(\delta _{22}(\infty )=0\) and \(\eta (\infty )=1\). These conditions introduce a smooth cutoff in the integral equations satisfied by the matrix elements of the MO matrix and ensure the existence of a unique solution.
5.3 Fits results
In addition to the six S-wave T-matrix parameters listed above, further parameters must be introduced which describe couplings to the \(\gamma \gamma \) channel. In the S-wave, two parameters \(b_l\), \(b_k\) were introduced as subtraction constants. Implementing the Adler zero condition we keep only \(b_k\) as an independent parameter. In the D-wave sector, we include the values of the tensor resonance couplings \(C^{a_2}_{2\gamma }\), \(D^{a_2}_{2\gamma }\), \(C^{f_2}_{2\gamma }\), \(D^{f_2}_{2\gamma }\) in the fitting as well as the mass and width of the \(a_2(1320)\) resonance: \(m_{a_2}\), \(\varGamma _{a_2}\). In total, we thus have 6+7 parameters to be fitted.
At first, we have kept the T-matrix parameters fixed to one of the sets of values determined previously in Ref. [27] (which, in particular, use assumed values for the pole positions of the two \(a_0\) resonances). It was not possible to obtain a good fit of the \(\gamma \gamma \) data in this manner: using these sets of parameter values one finds that the \(\pi \eta \) cross-section at the \(a_0(980)\) peak tends to be too large and the energy of the peak tends to be somewhat displaced as compared to experiment. Relaxing the T-matrix parameters, reasonably good fits become possible and we actually found two distinct minimums of the total \(\chi ^2\) combining the \(\pi \eta \) and the \(K_SK_S\) data,
In the case of \(\pi \eta \) the first number corresponds to the region \(E<1.1\) GeV. In this region fit II is better than fit I while the overall \(\chi ^2\) is slightly smaller in fit I (\(\chi ^2=428\)) than in fit II (\(\chi ^2=439\)). The \(\chi ^2\) is defined in a simple and naive way: the correlation matrix is assumed to be diagonal and the statistical and systematic errors provided by the Belle collaboration are added in quadrature. The searches for minimums were performed with the help of the computer code MINUIT [50].
The numerical values of the set of T-matrix parameters resulting from these two fits are shown in Table 2. One notices, in particular, that the value of the pole parameter \(m_8\) differs significantly in the two fits. Figure 5 shows the behaviour of the two phase-shifts \(\delta _{11}\), \(\delta _{22}\) and of the inelasticity parameters \(\eta \) as a function of energy, corresponding to the two fits. At low energies, the \(\pi \eta \) phase-shift changes sign and becomes negative. This low-energy behaviour is not anticipated in simple hadronic models of the \(\pi \eta \) amplitude (e.g. [51]) but it was also observed to emerge from fitting the \(\gamma \gamma \) data in Ref. [25].
In fit II the \(\pi \eta \) scattering length (defined as in [52]) is found to have the following value
For comparison, the first two terms in the chiral expansion of the scattering length give
using the same set of chiral parameters as in the unitary T-matrix.
At higher energies a clear difference between the two fits is the sharp increase of the \(\pi \eta \) phase-shift in fit I around \(E\simeq 1.32\) GeV, typical of a narrow resonance, which we will call \(a'_0\). Indeed, one finds a resonance pole in the T-matrix located on the third Riemann sheet with the value,Footnote 5
This resonance is lighter and narrower than the standard \(a_0(1450)\). Since the phase-shift increases from \(\pi /2\) to \(\pi \) (approximately) it gives rise to a sharp dip (instead of a peak) in the S-wave cross-section. Clearly then, our fit I is quite analogous to the best fit by the Belle collaboration [13] which displays a resonance (called \(a_0(Y)\) in that reference) which has very similar features while using a parametrisation rather different from ours.Footnote 6 The S-wave amplitude in fit II also has an \(a'_0\) resonance pole but it is heavier and much broader than that of fit I,
In this case, the width of the resonance is larger than the experimental one. One eventually does not expect a very accurate determination since the resonance lies close to the cutoff of the model. Finally, the effects of the two \(a_0\) resonances on the \(\pi \eta \) and \(K\bar{K}\) scalar form factors are shown in Fig. 11 in “Appendix A”.
Figures 6 and 8 show a sample of \(\gamma \gamma \rightarrow \pi \eta , K_SK_S\) differential cross sections comparing the experimental results with those from the two fits. The difference between the two fits is remarkably small, which is somewhat puzzling: why is the narrow \(a'_0\) resonance present in fit I not seen much more clearly? The reason for this arises from a specific interference effect between the \(J^{\lambda \lambda '}=0^{++}\) and the \(2^{++}\) amplitudes. We have seen already that the fast energy variation induced by the narrow \(a'_0\) coincides with that induced by the \(a_2(1320)\). More specifically, when \(\sqrt{s}\simeq 1.32\) GeV, the following relations hold approximately between the \(l_{0++}\) and the \(l_{2++}\) amplitudes in the case of fit I
The sum of the \(0^{++}\) and the \(2^{++}\) partial-wave amplitudes can effectively be absorbed into the \(2^{+-}\) resonance amplitude (since the angular functions satisfy the relation \(d^2_{00}(\theta )-1 = -\sqrt{6}\, d^2_{20}(\theta )\)) and no specific \(0^{++}\) resonance effect remains. We also note that the \(a'_0\) in fit I essentially decouples from the \(K\bar{K}\) channel such that no \(0^{++}\) resonance effect is seen in this channel either. Because of these fine-tuned relations between the scalar and the tensor resonances the fit I solution is most likely to be unphysical.
One can see from Fig. 6 that the \(\pi \eta \) differential cross-sections are well described except, however, in one energy bin: \(E=1.01\) MeV. At this energy, the cross-sections are underestimated by our model, see also Fig. 7 which shows the cross-sections integrated over \(\theta \). This discrepancy could be caused by an isospin breaking effect close to the resonance peak. Our model assumes isospin symmetry and uses \(m_K\equiv m_{K^+}\) (in order to correctly implement the kaon Born amplitudes) and the peak of the cross-section occurs exactly at \(E=2m_{K^+}\). Physically, however, the \(K^+\) and the \(K^0\) have slightly different masses which should lead to two cusps in the shape of the integrated cross-section. On the experimental side, the energy resolution should be smaller than \(2m_{K^+}-2m_{K^0}\simeq 8\) MeV in order to clearly observe this effect.
The values of the remaining seven parameters included in the fit (subtraction parameter \(b_k\), tensor mesons to \(2\gamma \) coupling constants, \(a_2\) mass and width) are collected in Table 3 below. The couplings \(D^T_{2\gamma }\) are found, as expected, to be smaller in magnitude than the \(C^T_{2\gamma }\) although this suppression is only by a factor of two in the case of fit I. The values are in qualitative agreement with the PDG expectations except, however, for the \(a_2\) mass which is shiftedFootnote 7 by approximately 10 MeV. This shift is easily seen to be caused by the presence of non-resonant contributions to the \(J=2\) amplitudes modelled here by the vector-meson exchanges, \(l^V_{2,\lambda \lambda '}\). The presence of this term is essential for obtaining a correct description of the amplitude in the \(\eta \) decay region \(s < (m_\eta -{m_{\pi }})^2\). In the 1 GeV energy region, one expects some modifications induced by higher mass exchanges, but we see no reason that it should be completely cancelled. We will consider a variation of this term as a source of error below.
5.4 Properties of the \(a_0\) resonances
In this section we consider in more detail the properties of the two \(a_0\) resonances which can be deduced from our analysis of the photon–photon data. We will focus on the results from fit II since fit I was argued not to be physically relevant. The formulas needed to define the T-matrix elements on the unphysical Riemann sheets are given in “Appendix D”. One may define coupling constants from the residues of the resonance poles (e.g. [54] in the context of photon–photon amplitudes), the couplings to the \(\pi \eta \) and \(K\bar{K}\) channels are thus defined as
and similarly for the third Riemann sheet. The coupling to the \(\gamma \gamma \) channel and the associated width are defined as
The following numerical results are found on sheet II for the position of the \(a_0(980)\) pole and the corresponding coupling constants
From the sheet III resonance, the position of the pole was given in Eq. (92) and the corresponding coupling constants have the following values
The errors quoted in the above formulas are those arising from the fitted parameters (for which MINUIT provides a correlation matrix). In addition to those, one must consider errors associated with further parameters on which the T-matrix depends, which we have assumed to be fixed up to now:
- 1.
Adler zero: the value of Adler zero was fixed to \(s_A=m_\eta ^2\) but its exact value is not known, we will vary it in the range \(s_A=m_\eta ^2\pm 3m_\pi ^2\).
- 2.
Set of the \(O(p^4)\) parameters \(L_i\): in the BE14 determination [49], their values are given with errors but there are strong correlations. In order to estimate an error from this source we used two different sets of \(L_i's\) which correspond to two best fits made with different assumptions (second and third column in Table 3 of Ref. [49]).
- 3.
The ratio of the scalar couplings \(c'_m/c'_d\) (see Eqs. (86), (87)) was varied between 1 and 3 to estimate an uncertainty from this source.
- 4.
Left-hand cut: (a) we added a contribution from a set of axial-vector mesonsFootnote 8 and (b) we varied the products of the couplings of the vector mesons by \(\pm 20\%\).
Table 4 gives the detailed list of the errors generated from these variations on the properties of the two \(a_0\) resonances.
5.5 The \(\eta \rightarrow {\pi ^0}\gamma \gamma \) decay amplitude
We consider now the modelled amplitudes in the kinematical region relevant for the decay \(\eta \rightarrow {\pi ^0}\gamma \gamma \). Including the \(J=0\) and \(J=2\) partial-waves,Footnote 9 the distribution as a function of the \(\gamma \gamma \) invariant mass squared, s, reads (Fig. 8)
also accounting for the \(J=0\) isospin violating amplitude \(\tilde{l}_{0,++}\).
Figure 9 illustrates the contributions of the various amplitudes corresponding to central values of the parameters fitted in the scattering region. The \(2^{+-}\) partial-wave dominates over the other ones near \(s=0\) because the \(J^{++}\) amplitudes are suppressed by the soft-photon zero. The relative role of the S-wave increases with the energy and starts to be dominating above the \(\pi ^+\pi ^-\) threshold. The figure also shows that the isospin-violating S-wave generates a visible cusp at this threshold. The central value of the decay width generated by our amplitudes is \(\varGamma =0.237\) eV which is on the low side of the most recent experimental determinations \(\varGamma _{exp}=0.285\pm 0.031\pm 0.061\) eV (Crystal Ball at the AGS [14]), \(\varGamma _{exp}=0.33\pm 0.03\) eV (MAMI [15]). A smaller value was reported by the KLOE-2 collaboration [56] but it has not been confirmed and the data is currently being reanalysed. The amplitudes in the \(\eta \) decay region are very sensitive to the precise position of the Adler zero. This is illustrated in Fig. 10 showing the effect of varying \(s_A\) in the range \([m_\eta ^2-3m_\pi ^2,m_\eta ^2+3m_\pi ^2]\) and comparing with the experimental results. Given somewhat more precise data the value of \(s_A\) could be included in the fitting. The amplitudes in the decay region are also sensitive to the vector meson coupling constants, the value of which dominate the energy distribution near \(s=0\). Accounting for these main errors we would predict
6 Conclusions
In this work we have reconsidered the properties of the light isovector scalar resonances as can be determined from photon–photon scattering experimental results. For this purpose, we have implemented a standard Muskhelishvili–Omnès integral representation for the \(J=0\) amplitude in which the left-cut is modelled from light vector meson exchanges. The underlying T-matrix satisfies unitarity with two channels (\(\pi \eta ,K\bar{K}\)) and involves six phenomenological parameters. In the case of the \(J=2\) amplitudes the constraints from unitarity are more difficult to implement, a cruder description is used which consists in simply adding the cross-channel vector-exchange and the direct channel tensor resonance amplitudes. In order to constrain the free parameters as unambiguously as possible we performed fits to both \(\pi \eta \) and \(K_SK_S\) data, for which high-statistics data below 1.4 GeV are available at present, and we found two different acceptable solutions to the minimisation. Both solutions are also compatible with the available \(K^+K^-\) data from the ARGUS collaboration [44].
In one of our fits the S-wave amplitude displays a light and narrow \(a'_0\) resonance exactly similar to the one found in the Belle analysis. While this is mathematically allowed we have argued that the fit which displays a broad \(a'_0\) is likely to be more physical. Concerning the \(a_0(980)\) resonance, we find that a rather conventional picture i.e. a pole on the second sheet with a mass and width compatible with the PDG and coupling to both the \(\pi \eta \) and the \(K\bar{K}\) channels is perfectly compatible with both the \(\pi \eta \) and the \(K_SK_S\) data. This is in contrast with the Belle analysis which uses an elastic Breit-Wigner description and also with the recent analysis of Ref. [18] in which the mass and width are found to be both significantly larger than the PDG values. Data with a better energy resolution would be useful to resolve these remaining ambiguities. The \(\gamma \gamma \rightarrow {K^+}{K^-},K_SK_S\) cross-sections close to the \(K\bar{K}\) thresholds are also very sensitive to the position of the \(a_0(980)\). Our results in this energy region are in qualitative agreement with the chiral-unitary calculations from Ref. [29] and with the estimates made in Ref. [57] but not with those from Ref. [26]. Experimental data in this near-threshold region would obviously be very constraining. Finally, it will be quite interesting to see how the pole position determined in a lattice QCD simulation [22] evolves when the value of \(m_\pi \) is decreased.
Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors’ comment: This work is concerned about extracting information from experimental data already available from data repositories using theoretical modelling. The relevant numerical results of this work are presented in tables in the paper.]
Notes
A renewed interest in such processes, both theoretical and experimental, is motivated, more generally, by the problem of evaluating the so-called hadronic light-by-light contribution to the muon \(g-2\) (see e.g. [6]).
A further difference is that the soft photon constraints at \(s=0\) are not imposed in the dispersion relations used in [18].
The \(\pi \eta \) and \(K\bar{K}\) scalar form factors are expected to go like 1/s at infinity (up to log’s) in QCD and must be proportional to the matrix elements of the MO matrix (up to a polynomial). Consequently, the MO matrix elements must vanish when \(s\rightarrow \infty \) like 1/s at least, which is assumed in Eq. (56).
For two-body decay modes we take \(\varGamma _{12}(s)=\varGamma ^0_{12}\times (q_{12}(s)/q_{12}(s_0) )^{2l+1} W_l( q_{12}(s)R)/W_l( q_{12}(s_0)R)\) with \(s_0=m_{a_2}^2\) and \(l=2\) as in [13] while for \(\rho \pi \) we take \(l=1\) and \(W_1=1\). Such Breit–Wigner forms for tensor resonances are widely used but do not have good analyticity properties. In order to define the BW amplitudes below the thresholds we set the corresponding momenta to zero i.e. \(q_{12}(s)=0\) if \(s < (m_1+m_2)^2\).
The errors quoted here are the statistical ones as evaluated by MINUIT. Further errors introduced by varying intrinsic parameters of the model will be considered below.
The Belle collaboration parametrise the S-wave amplitude \(l_{0++}\) as a sum of two Breit–Wigner functions plus a polynomial quadratic in the energy E. This representation involves 11 free parameters.
A recent determination (from the complex pole) of the \(a_2\) mass using data from the COMPASS experiment also finds a shift compared to the PDG value, however in the opposite direction [53].
The coupling constants of the C-odd axial-vectors are not very precisely known (see [55]). We made a simple estimate taking a mass \(M_A=1.2\) GeV and an effective coupling \(G_A=0.30(C_{\rho \pi }C_{\rho \eta }+C_{\omega \pi }C_{\omega \eta })\).
Contributions from \(J\ge 4\) partial-waves can be checked to be completely negligible.
References
E. Klempt, A. Zaitsev, Phys. Rept. 454, 1 (2007). arXiv:0708.4016
A. Astier, L. Montanet, M. Baubillier, J. Duboc, Phys. Lett. B 25, 294 (1967)
R. Ammar, R. Davis, W. Kropac, J. Mott, D. Slate, B. Werner, M. Derrick, T. Fields, F. Schweingruber, Phys. Rev. Lett. 21, 1832 (1968)
S.M. Flatté, Phys. Lett. B 63, 224 (1976)
M. Tanabashi et al. (Particle Data Group), Phys. Rev. D 98(3), 030001 (2018)
F. Jegerlehner, A. Nyffeler, Phys. Rept. 477, 1 (2009). arXiv:0902.3360
F. Low, Phys. Rev. 110, 974 (1958)
H.D.I. Abarbanel, M.L. Goldberger, Phys. Rev. 165, 1594 (1968)
S.L. Adler, Phys. Rev. 139, B1638 (1965)
D. Morgan, M. Pennington, Phys. Lett. B 192, 207 (1987)
D. Morgan, M.R. Pennington, Phys. Lett. B 272, 134 (1991) [Erratum: Nucl. Phys. B 376, 444 (1992)]
D. Antreasyan et al. (Crystal Ball)., Phys. Rev. D 33, 1847 (1986)
S. Uehara et al. (Belle)., Phys. Rev. D 80, 032001 (2009). arXiv:0906.1464
S. Prakhov et al., Phys. Rev. C 78, 015206 (2008)
B.M.K. Nefkens et al. (A2 at MAMI), Phys. Rev. C 90(2), 025206 (2014). arXiv:1405.4904
S. Uehara et al. (Belle), PTEP 2013(12), 123C01 (2013). arXiv:1307.7457
K. Abe et al. (Belle), Eur. Phys. J. C 32, 323 (2003). arXiv:hep-ex/0309077
I. Danilkin, O. Deineka, M. Vanderhaeghen, Phys. Rev. D 96(11), 114018 (2017). arXiv:1709.08595
I.V. Danilkin, L.I.R. Gil, M.F.M. Lutz, Phys. Lett. B 703, 504 (2011). arXiv:1106.2230
N.L. Muskhelishvili, Singular Integral Equations (P. Noordhof, Groningen, 1953)
R. Omnès, Nuovo Cim. 8, 316 (1958)
J.J. Dudek, R.G. Edwards, D.J. Wilson (Hadron Spectrum), Phys. Rev. D 93(9), 094506 (2016). arXiv:1602.05122
M. Lüscher, Nucl. Phys. B 354, 531 (1991)
Z.H. Guo, L. Liu, U.G. Meißner, J.A. Oller, A. Rusetsky, Phys. Rev. D 95(5), 054004 (2017). arXiv:1609.08096
N. Achasov, G. Shestakov, Phys. Rev. D 81, 094029 (2010). arXiv:1003.5054
L.Y. Dai, M.R. Pennington, Phys. Rev. D 90(3), 036004 (2014). arXiv:1404.7524
M. Albaladejo, B. Moussallam, Eur. Phys. J. C 75(10), 488 (2015). arXiv:1507.04526
R. García-Martín, B. Moussallam, Eur. Phys. J. C 70, 155 (2010). arXiv:1006.5373
J.A. Oller, E. Oset, Nucl. Phys. A 629, 739 (1998). arXiv:hep-ph/9706487
E. Oset, J.R. Pelaez, L. Roca, Phys. Rev. D 67, 073013 (2003). arXiv:hep-ph/0210282
E. Oset, J.R. Pelaez, L. Roca, Phys. Rev. D 77, 073001 (2008). arXiv:0801.2633
J.N. Ng, D.J. Peters, Phys. Rev. D 46, 5034 (1992)
L. Ametller, J. Bijnens, A. Bramon, F. Cornet, Phys. Lett. B 276, 185 (1992)
M. Jacob, G.C. Wick, Ann. Phys. 7, 404 (1959)
M. Jacob, G.C. Wick, Ann. Phys. 281, 774 (2000)
J. Kambor, C. Wiesendanger, D. Wyler, Nucl. Phys. B 465, 215 (1996). arXiv:hep-ph/9509374
A. Anisovich, H. Leutwyler, Phys. Lett. B 375, 335 (1996). arXiv:hep-ph/9601237
G. Colangelo, S. Lanz, H. Leutwyler, E. Passemar, Eur. Phys. J. C 78(11), 947 (2018). arXiv:1807.11937
G. Colangelo, S. Lanz, H. Leutwyler, E. Passemar, Phys. Rev. Lett. 118(2), 022001 (2017). arXiv:1610.03494
S.J. Brodsky, G.P. Lepage, Phys. Rev. D 24, 1808 (1981)
M. Diehl, P. Kroll, C. Vogt, Phys. Lett. B 532, 99 (2002). arXiv:hep-ph/0112274
R. Unterdorfer, H. Pichl, Eur. Phys. J. C 55, 273 (2008). arXiv:0801.2482
J. Gasser, A. Rusetsky, Eur. Phys. J. C 78(11), 906 (2018). arXiv:1809.06399
H. Albrecht et al. (ARGUS), Z. Phys. C 48, 183 (1990)
A. Dobado, M.J. Herrero, T.N. Truong, Phys. Lett. B 235, 134 (1990)
J. Oller, E. Oset, A. Ramos, Prog. Part. Nucl. Phys. 45, 157 (2000). arXiv:hep-ph/0002193
G. Ecker, J. Gasser, A. Pich, E. de Rafael, Nucl. Phys. B 321, 311 (1989)
J. Gasser, H. Leutwyler, Nucl. Phys. B 250, 465 (1985)
J. Bijnens, G. Ecker, Ann. Rev. Nucl. Part. Sci. 64, 149 (2014). arXiv:1405.6488
F. James, M. Roos, Comput. Phys. Commun. 10, 343 (1975)
D. Black, A.H. Fariborz, J. Schechter, Phys. Rev. D 61, 074030 (2000). arXiv:hep-ph/9910351
V. Bernard, N. Kaiser, U.G. Meißner, Phys. Rev. D 44, 3698 (1991)
A. Jackura et al. (JPAC, COMPASS), Phys. Lett. B 779, 464 (2018). arXiv:1707.02848
D. Morgan, M.R. Pennington, Z. Phys. C 37, 431 (1988) [Erratum: Z. Phys. C 39, 590 (1988)]
P. Ko, Phys. Rev. D 47, 3933 (1993)
B. Di Micco et al. (KLOE), Acta Phys. Slov. 56, 403 (2006)
N.N. Achasov, G.N. Shestakov, JETP Lett. 96, 493 (2012). arXiv:1210.0739
J. Bijnens, F. Cornet, Nucl. Phys. B 296, 557 (1988)
J.F. Donoghue, B.R. Holstein, Y.C. Lin, Phys. Rev. D 37, 2423 (1988)
F. Guerrero, J. Prades, Phys. Lett. B 405, 341 (1997). arXiv:hep-ph/9702303
R. Escribano, S. Gonzàlez-Solís, R. Jora, E. Royo, (2018). arXiv:1812.08454
Acknowledgements
We would like to thank prof. S. Uehara for many clarifying explanations about the Belle experiments. JXL thanks Prof. Lisheng Geng for useful discussions. JXL is partly supported by the National Natural Science Foundation of China under Grant Nos. 11735003, 11975041, 11961141004 and by the Fundamental Research Funds for the Central Universities. He also gratefully acknowledges the financial support from the China Scholarship Council. This work is supported in part by the European Union’s Horizon2020 research and innovation programme (HADRON-2020) under the Grant Agreement \(\hbox {n}^\circ \) 824093.
Author information
Authors and Affiliations
Corresponding author
Appendices
Appendix A: Scalar form factors
In this section we consider how the two isovector scalar form factors get modified as compared to the results of Ref. [27] when using the set of T-matrix parameters as determined here from the \(\gamma \gamma \) data. The form factors were defined as
where \(B_0\) is the chiral coupling proportional to the quark condensate [48]
Under the assumption of two-channel unitarity these form factors, we recall, are simply related to the MO matrix,
The values at \(s=0\) are estimated from the chiral expansion at order \(p^4\). In Ref. [27] the set of \(L_i\) values from Ref. [49] based on a \(p^4\) fit was used while here we preferred to use a set from a \(p^6\) fit (BE14 set), which gives
Interesting quantities which are related to the scalar form factors are the scalar radii which are proportional to their derivatives at \(s=0\)
At chiral order \(p^4\) the scalar radii depend on a single low-energy coupling, \(L_5\). Using its value from the BE14 set, \(L^r_5=(1.01\pm 0.06)\times 10^{-3} \), gives
Figure 11 shows the absolute values of the form factors using the BE14 chiral couplings and the set of T-matrix parameters fitted to the \(\gamma \gamma \) data. The \(a_0(980)\) resonance displays a smaller peak than in Ref. [27] and its position is slightly shifted, reflecting the modified position of the pole. The \(a'_0\) resonance also differs: in the case of the “narrow \(a'_0\)” fit solution it shows up as a dip in the \(\eta \pi \) form factor while in the case of the “broad \(a'_0\)” solution its effect is hardly visible. Concerning the scalar radii, the dispersive amplitudes give
These results correspond to fit solutions with a “broad” \(a'_0\). As compared with the the chiral \(O(p^4)\) result, the central value is somewhat too small in the case of \(\eta \pi \) and too large in the case of \(K\bar{K}\).
Appendix B: Chiral expansion results for \(\gamma \gamma \rightarrow P_1P_2\) amplitudes
1.1 Appendix B.1: Order \(p^4\)
We collect below the expressions of the photon–photon amplitudes in the chiral expansion at order \(p^4\). The chiral order \(p^2\) coincides with scalar QED and gives rise to the Born amplitudes for \({\pi ^+}{\pi ^-}\) and \({K^+}{K^-}\) which were given in Eq. (35). The leading contributions to the amplitudes which involve two neutral mesons: \(\pi ^0\pi ^0\), \({K^0}{\bar{K}^0}\) and \({\pi ^0}\eta \) appear at order \(p^4\). They were computed in Refs. [33, 58,59,60]. The basic one-loop function which occurs in these \(p^4\) amplitudes can be written as,
for small z values, \(|z| \ll m_P^2\), it behaves as
The expression for the \(\gamma \gamma \rightarrow \pi ^0\pi ^0\) scalar amplitudes at \(O(p^4)\) read,
These formulas illustrate general features: the \(O(p^4)\) part of the A amplitudes have a simple structure involving a sum of products of meson-meson amplitudes at order \(p^2\) with one-loop functions \(G_P\) while the B amplitudes vanish at this order. Next, the expression for \(\gamma \gamma \rightarrow {K^0}{\bar{K}^0}\) at \(O(p^4)\) reads
The NLO expression for the \(\gamma \gamma \rightarrow {\pi ^0}\eta \) amplitude, which is of particular interest here, was worked out in Ref. [33]. They considered a chiral framework which also includes the \(\eta '\) meson as a light meson in addition to the pseudo-Goldstone octet. The relation between the \(\eta \) and \(\eta '\) mesons and the octet and singlet chiral fields involves a mixing angle \(\theta \),
denoting \(c_\theta =\cos \theta \), \(s_\theta =\sin \theta \). The expression for the \(\gamma \gamma \rightarrow {\pi ^0}\eta \) amplitude as a function of \(\theta \) reads,
with \(t_\theta =\tan \theta \). The isospin conserving part in Eq. (B.6) agrees with Ref. [61]. It reproduces the result from Ref. [33] upon setting \(s_\theta =-1/3\), \(c_\theta =\sqrt{8}/3\) and neglecting the terms proportional to \(t_\theta \) (which are not really negligible). In the standard ChPT one must set \(c_\theta =1\), \(s_\theta =0\) at leading order. Doing so, the amplitude simplifies significantly and reads
where the Gell–Mann–Okubo mass relation was used and the kaon loop contribution to the isospin violating part was included for completeness. The quantity \(B_0(m_d-m_u)\) is given at leading chiral order by
Let us quote also the formula for \(\gamma \gamma \rightarrow \eta \eta \)
Finally, the order \(p^4\) contributions to the \(\gamma \gamma \rightarrow {\pi ^+}{\pi ^-}\) and \({K^+}{K^-}\) amplitudes read
and
Combining Eqs. (B.4) and (B.11) we get the NLO part of the \(\gamma \gamma \rightarrow (K\bar{K})^{I=1},\,(K\bar{K})^{I=0}\) amplitudes
1.2 Appendix B.2: Soft pion limit (isospin conserving part)
Let us consider a limit when the pion in the \(\gamma \gamma \rightarrow {\pi ^0}\eta \) amplitude becomes “soft” i.e. \(p_1=0\). This limit is unphysical but it allows to obtain exact results which should hold for the physical amplitude modulo \(O(m_\pi ^2)\) corrections. In this limit, firstly, one has
and using standard soft pion methods one shows that the amplitude \(\gamma \gamma \rightarrow {\pi ^0}\eta \) vanishes exactly at this point. When \(p_1=0\) the two tensors \(T^{\mu \nu }_1\), \(T^{\mu \nu }_2\) become degenerate
such that the soft limit amplitude, in tensorial form, reads
The combination which appears in Eq. (B.15) is the helicity amplitude \(L_{++}\) in the soft pion limit. This implies that the physical helicity amplitude \(L_{++}\) with \(t=u\) should also have an Adler zero as a function of s,
with \(s_A=m_\eta ^2+O(m_\pi ^2)\). Let us consider the implication of this result for the \(j=0\) partial wave. The expansion of the helicity amplitude with \(\cos \theta =0\) (which corresponds to \(t-u=0\)) reads
When \(s=m_\eta ^2\), the amplitudes with \(j\ge 2\) are suppressed by the angular momentum barrier factor: \(l_{2,++}(m_\eta ^2)=O(m_\pi ^2)\), \(l_{4,++}(m_\eta ^2)=O(m_\pi ^4), \ldots \). The soft pion condition therefore implies that the \(j=0\) partial-wave should satisfy \(l_{0,++}(m_\eta ^2)=O(m_\pi ^2)\) such that an Adler zero should be present in this partial-wave amplitude,
1.3 Appendix B.3: Soft pion limit (isospin violating part)
In the case of the isospin violating piece of the \(\gamma \gamma \rightarrow {\pi ^0}\eta \) amplitude, the soft pion limit does not lead to an Adler zero. Let us define this amplitude in terms of the following matrix element,
which involves the isospin violating part of the QCD Lagrangian. Using the usual soft pion techniques, the soft pion limit of this amplitude is expressed in terms of the commutator
where \(Q_5^3\) is the axial charge
The commutator is easily worked out and gives
The matrix element which appears on the right-hand side is a non-perturbative quantity. We can estimate it using the chiral expansion at NLO. The pseudoscalar operator \(p_0=i(\bar{u}\gamma ^5{u}+\bar{d}\gamma ^5{d}) \) can couple either to a single \(\eta \) meson or to \(\eta \pi ^+\pi ^-\). Representative one-loop diagrams are shown in Fig. 12. Computing them gives
As one might expect, using this result in the soft pion relation (B.22) reproduces the isospin violating amplitude at chiral order \(p^4\) when \(s=m_\eta ^2\) (see Eq. (B.7)), up to \(O(m_\pi ^4/m^4_\eta )\) terms. In practice, one can relate \((m_d-m_u)B_0\) to the QCD kaon mass difference \((m^2_{K^0}-m^2_{K^0})_{QCD}\) and then use Eq. (28) in order to express \(\tilde{l}_{0++}(m_\eta ^2)\) in terms of \(\epsilon _L\). The chiral calculation (B.23) of the matrix element should be reliable at the \(20{-}30\%\) level since there are no potentially large \(p^6\) contributions from light vector meson exchanges in this case.
Appendix C: Dispersive Dispersive \(\gamma \gamma \rightarrow \pi \pi ,\ (K\bar{K})_{I=0}\) S-waves
Input for the \(\gamma \gamma \rightarrow (K\bar{K})_{I=0}\) S-wave is needed in order to reconstruct the amplitudes for the physical \(K^+K^-\) and \(K_SK_S\) states. We briefly recall the dispersive coupled-channel representation for the \(I=0\) S-wave taken from Ref. [28]. It has the following form
which is analogous to eqs. (59) for the \(I=1\) amplitudes but more subtraction parameters had to be introduced. In Eq. (C.1) \(\varvec{\varOmega }^{(0)}\) is the \(I=0\) Omnès matrix and the functions \(L^{(0)}_i(s)\) (\(R^{(0)}_i(s)\)) are dispersive integrals on the left (right)-hand cuts which are analogous to the corresponding integrals in Eq****. (59). The parameter \(b^0_K\) was simply fixed by matching to the \(O(p^4)\) amplitude at \(s=0\) (see Eq. (B.12)). The three parameters \(b^0\), \(b^{'0}\), \(b^{'0}_K \) where then determined from fits to \(\gamma \gamma \rightarrow {\pi ^+}{\pi ^-}, {\pi ^0}{\pi ^0}\) data. Having done so the \(K\bar{K}\) amplitude \(k^0_{0++}\) was generated, which we have used in the present work together with the \(I=1\) amplitude. Figure 13 shows the dispersive results for these two \(K\bar{K}\) amplitudes. At low energies the amplitude \(k^0_{0++}\) is significantly larger in magnitude than \(k^1_{0++}\) because it gets isoscalar \(\pi \pi \) rescattering contributions (which contain the broad \(\sigma \) resonance). In the low-energy region, furthermore, the Born amplitude is much larger than the rescattering contributions, in accordance with the chiral counting. In the 1 GeV region the \(I=0,1\) rescattering amplitudes are rather similar which leads to a strong suppression of the \({K^0}{\bar{K}^0}\) cross-section as compared to the \({K^+}{K^-}\) one close to the \(K\bar{K}\) threshold.
Appendix D: The T-matrix on the four Riemann sheets
The second sheet extension of a matrix elements \(T_{ij}\) is defined such as to continue \(T_{ij}(z)\) analytically across the cut, below the first inelastic threshold,
Using the elastic unitarity relation it easy to find an explicit expression for \(T^{(II)}_{ij}\),
with
The second sheet extensions of the \(\gamma \gamma \) amplitudes \(l_{0++}\) and \(k^1_{0++}\) are expressed in a similar way,
The third sheet extension of the T-matrix elements are defined such as to obey continuity equations across the unitarity cut above \(4m_K^2\)
They are easily expressed in matrix form
where \(\tilde{\sigma }_{KK}(z)=\sqrt{4m_K^2/z-1}\). Finally, the fourth sheet extensions of the T-matrix elements are given by
they satisfy continuity equations with \(T^{(II)}_{ij}(z)\) when \(s > 4m_K^2\) and with \(T^{(III)}_{ij}(z)\) when \((m_\eta +{m_{\pi }})^2< s < 4m_K^2\).
In our model, the right-cut structure of T-matrix is generated by the pair of loop functions \(\bar{J}_{\pi \eta }(s)\), \(\bar{J}_{KK}(s)\) which can be written as
This expression defines \(\bar{J}_{12}\) on the first Riemann sheet and shows that it is analytic except for a cut on \([(m_1+m_2)^2,\infty ]\). The second sheet extension of \(\bar{J}_{12}\) is given by
which satisfies the continuity equation \(\bar{J}^{II}_{12}(s-i\epsilon )= \bar{J}_{12}(s+i\epsilon )\) when \(s \ge (m_1+m_2)^2\). The T matrix on the four Riemann sheets can then also be defined as follows in terms of the pair of loop functions
Using the K-matrix type representation (83) of the T-matrix one easily verifies that these definitions satisfy the relevant equations (D.2), (D.6), (D.7) written above.
Rights and permissions
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 copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Funded by SCOAP3
About this article
Cite this article
Lu, J., Moussallam, B. The \(\pi \eta \) interaction and \(a_0\) resonances in photon–photon scattering. Eur. Phys. J. C 80, 436 (2020). https://doi.org/10.1140/epjc/s10052-020-7969-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epjc/s10052-020-7969-8