The $\pi\eta$ interaction and $a_0$ resonances in photon-photon scattering

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\to \pi\eta$ and $\gamma\gamma\to 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\Gamma/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\to\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.


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 twoparticle) threshold (see the review [1]). The a 0 (980) resonance was discovered a long time ago and seen in both the KK and the πη channels [2,3] but its properties are still imprecisely known. This is partly because πη production experiments using an η beam, analogous to those which have allowed to determine the ππ or π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 πη − KK 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. They are free of initial-state interactions and satisfy dispersion relations which can be constrained, in the case of ππ or πη by both soft photon [6,7] and soft pion [8] low-energy theorems. As was illustrated in the seminal papers [9,10] 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 γγ → ηπ 0 cross-sections were first measured by the Crystal Ball collaboration [11]. Measurements with much higher statistics were recently performed by the Belle collaboration [12]. There are also experimental results on the η → π 0 γγ decay amplitude [13,14]. Recent experimental data exist also for γγ → K 0K0 [15] and γγ → K + K − [16] 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, Γ = 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 [17] based on a specific meson-meson T -matrix model [18] applied to the πη S-wave, from which the corresponding γγ amplitude is deduced from a Muskhelishvili-Omnès (MO) construction [19,20]. 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 [21] who implemented a coupled-channel generalisation of Lüscher's single-channel method [22]. This result, however, corresponds to a pion mass m π = 391 MeV and it is not known how the pole location evolves upon varying m π (see, however, ref. [23]). Fits to the Belle data with a conventional a 0 (1450) were performed in ref. [24] 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. [25] focusing mostly on the I = 0, 2 channels. In the I = 1 sector, they considered the KK channel but not πη.
We perform here a global fit which takes into account both πη and KK 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. [26] 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. [17], we differ only by using SU (2) chiral symmetry, which allows to fix one of the subtraction constants through a soft pion theorem. 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 γγ → (KK) I=1 amplitudes are then combined with I = 0 amplitudes taken from a previous work [27] which considered γγ → (ππ) I=0,2 , (KK) I=0 in order to reconstruct the physical K + K − , K 0K0 amplitudes. A global fit of photon-photon to meson-meson data was performed some time ago based on unitarised chiral amplitudes [28]. 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 γγ → πη amplitudes is that they can be probed experimentally both in the scattering regime: s γγ ≥ (m η + m π ) 2 and in the decay regime: 0 ≤ s γγ ≤ (m η − m π ) 2 . In the decay region the amplitudes are largely dominated by the light vector meson exchanges in the cross-channels: γπ → ρ, ω → γη, which were first computed in ref. [29]. In this low-energy region the rescattering contributions can be estimated in the SU (3) chiral expansion [30]. They proceed essentially via γγ → K + K − → π 0 η but, at low energy, the π + π − contribution (which is isospin violating) is not negligible and must be included. We will present a dispersive calculation of this contribution which gives rise to a cusp in the energy distribution dΓ η→π 0 γγ /ds γγ .
The plan of the paper is as follows. After recalling some general properties of photonphoton amplitudes and introducing the notation in sec. 2, we write the unitarity relations for the partial-waves in sec. 3 and present the modelling of the left-hand cut of these. In sec. 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 sec. 5 and the information that can be deduced on the a 0 resonances are discussed.

Basic ingredients 2.1 Kinematics
We consider two-photon to two-meson scattering amplitudes with M 1 M 2 = πη or KK. The Mandelstam variables are defined as usual by The various physical regions in the s, t − u Mandelstam plane for the γγ → πη, γπ → γη and η → γγπ processes are shown on fig. 1. In the centre-of-mass frame of the two photons the momenta are expressed as where Takingẑ to a unit vector along the z axis andv to be a unit vector with polar angles θ, φ such thatẑ ·v = cos θ, we can express t, u in terms of cos θ as, The polarisation vectors of the two photons 1 (q 1 , λ), 2 (q 2 , λ ) in this frame read and they satisfy 1 · q 1 = 2 · q 2 = 1 · q 2 = 2 · q 1 = 0 .

Amplitudes
We denote the γγ → ηπ helicity amplitude as L λλ (we have factored out the electric charge e and the dependence on the azimuthal angle φ) while for the γγ → KK amplitudes we use the same notation as in previous work: K c λλ (s, t) for charged kaons and K n λλ (s, t) for neutral kaons. Helicity amplitudes are convenient for performing the partial wave expansion [31]. They can be expressed in terms of tensor amplitudes using the reduction formulas, e.g.
By gauge invariance, the tensor amplitude W µν must satisfy the two Ward identities One can form two independent tensors which satisfy (9) which we take as The tensor amplitude W µν can then be expressed in terms of two scalar amplitudes A, B W µν (q 1 , q 2 , p 1 , p 2 ) = A(s, t, u)T µν 1 + B(s, t, u)T µν 2 (12) which satisfy dispersion relations. Using eqs. (5) (6) one can easily express the helicity amplitudes L λλ in terms of the two scalar amplitudes, Assuming unpolarised photon beams the differential cross-section reads, Concerning the η → πγγ 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

Right-hand cut and unitarity relations
It is convenient to collect the three I = 1 scattering amplitudes πη → πη, πη → KK and KK → KK into a 2 × 2 matrix and we can define the partial wave expansions as where θ is the scattering angle in the centre-of-mass system. The unitarity relation for the partial waves are easily derived and read 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 Im l 0++ (s) which also give the discontinuities of the partial-waves (extended to complex values of s) across the right-hand cut.

ππ (isospin violating) contribution in S-wave unitarity
The unitarity relation written above (25) collects the contributions from the πη and the KK states. Let us also consider the contribution from the π + π − state, which has the form where dΦ is the phase-space integration measure. This contribution, being proportional to the ππ → ηπ amplitude, is isospin violating but it is enhanced at low energy due the large size of the γγ → π + π − amplitude [30]. The ChPT evaluation performed in ref. [30] amounts to using the O(p 2 ) tree-level amplitudes for both ππ → ηπ and γγ → π + π − 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 GeV 2 . In such a region, the π + π − → ηπ 0 amplitude can be approximated in terms of three one-variable functions M 0 , M 1 , M 2 (see [32,33]), An isospin violating parameter L has been factorised which may be taken as [34] The three M I functions obey a set of coupled Khuri-Treiman integral equations. We refer to ref. [34] for a complete review of work on this subject and an updated evaluation of the K 0 −K + mass difference based on experimental data on η → 3π decays: m 2 K 0 − m 2 K + QCD = (6.24 ± 0.38) · 10 −3 GeV 2 , which gives The amplitudes corresponding to a given ππ isospin state I, I z M I Iz ≡ ηπ|T |ππ; I I z (30) are easily expressed using crossing symmetry and the Wigner-Eckart theorem. In the unitarity relation (26) the amplitudes with I = 0, 2, M 00 are needed, M 20 , 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 ππ contribution to the unitarity relation is finally expressed as follows, where disc [l 0,++ (s)] ππ ≡ l 0, and h I 0,++ (s) are the two S-wave γγ → (ππ) I amplitudes with I = 0, 2. This ππ discontinuity of l 0++ can be estimated from eq. (33) using inputs from ref. [27] for γγ → ππ and from [26] for ππ → ηπ. 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 η − m π ) 2 induced by the endpoint of the left-hand cut in the functionsM 0 ,M 2 which overlaps with the right-hand cut as a result of the instability of the η. As a further consequence, the phase of the partial-waves M I +M I violate Watson's theorem and do not cancel with the phases of h I 0,++ such that the discontinuity has both a real and an imaginary part.

Left-hand cut: Born amplitudes
A left hand cut in the γγ partial-waves is generated by singularities in the cross-channels γP 1 → γP 2 . In the case when P 1 = P 2 = K + or π + the leading singularity is the kaon or pion pole and the corresponding (so-called Born) amplitudes read  Figure 2: Discontinuity of the γγ → πη S-wave amplitude (real part) across the π + π − cut computed using dispersive results for the γγ → ππ and ππ → ηπ amplitudes, compared to the chiral O(p 4 ) result.
with P = K + , π + . We will need the I = 1 component of the K + Born amplitude projected on the S-wave which reads, with β P (s) = 1 − 4m 2 P /s. We also recall the expressions of the J = 2 Born helicity amplitudes

Left-hand cut: vector meson exchanges
Leading contributions to the left-hand cut in the π 0 η amplitudes are modelled from the ρ, ω, φ vector meson exchanges. We define the V P γ coupling constants G V P 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 γγ → 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 2 V . 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 = ±1 when s = s V , which thus corresponds to a branch point of L V (s). Another endpoint singularity occurs when s = ∞. 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: −∞ < 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 γγ → (KK) 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 V P are collected in table 1 below.

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 → ∞ 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 ∼ u ∼ −s ∼ −∞. In this regime, L ++ (s, t) can be estimated from QCD-based methods [35,36] 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 α V 0.5, α 0.9 GeV −1 and β V (t) is a smooth function when t < 0. Assuming only that the integrand evolves smoothly between these two regimes when 0 ≤ z ≤ 1, one deduces that the l 0++ (s) should not grow faster than √ s when s → ∞. Furthermore, the J = 0 partial-waves obey soft-photon theorems [6,7] which imply that the ratios l 0++ (s)/s and (k 1 0++ (s) − k 1,Born 0++ (s))/s remain finite when s → 0. Therefore, they can be expressed as unsubtracted dispersion relations in terms of the left-hand and right-hand cuts discontinuities, with m + = m π + m η . 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 KK threshold and that two-channel unitarity remains a sufficiently good approximation up to √ 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 matrix 1 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 Multiplying the amplitudes (l 0++ , k 0++ ) by the inverse of the matrix Ω 0 removes the righthand cuts. We can use this property in writing once-subtracted dispersion relations for the two functions This provides MO-type dispersive representations for the γγ amplitudes l 0++ (s), k 0++ (s) in terms of their imaginary parts on the left-cut and two parameters, 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 A) 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, The dispersive representations (59), (60) are equivalent to those considered previously in ref. [17]. 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 appendixA) 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 ± 0.08) · 10 −3 taken from ref [37]. 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 KK amplitude k 1 0++ is shown in appendix B (see fig. 12).

Dispersive construction of the isospin-violating S-wave
In sec. 3.2 we have considered the unitarity contribution to the γγ → πη S-wave amplitude induced by the π + π − intermediate state, which is isospin violating. Let us callL ++ the isospin-violating part of the γγ → πη helicity amplitude andl 0++ the corresponding S-wave.
L ++ can be defined as a matrix element, We will attempt here to estimate the amplitudel 0++ at low energy only and we write a dispersive representation keeping only the contribution from the ππ 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λ in eq. (66) by using the soft pion limit. This limit provides a relation between the amplitudel 0++ (s = s A ) and a matrix element of the pseudo-scalar operator p 0 = i(ūγ 5 u +dγ 5 d) which, in turn, can be estimated using ChPT. This is detailed in appendix A.3. Using eqs. (121) ,(122) from this appendix we obtain the following result forλλ where the loop functions G π , G K are given in eq. (100). The discontinuity disc [l 0++ (s)] ππ (given in eq. (33)) has a singularity at the pseudo-threshold s = (m η − m π ) 2 induced by the ππ → ηπ partial-wave (see fig. 2). The integral in eq. (66), however, is finite. It is defined by using the m 2 η + i limiting prescription exactly in the same way as those which appear in the Khuri-Treiman equations for the η → 3π amplitude (see e.g. [32,38]).
The result, in the low energy region relevant for η → πγγ is shown in fig. 4 and compared to the corresponding chiral O(p 4 ) result from eq. (105). The chiral and dispersive real parts agree at s = m 2 η , as a result of the soft pion relation, but the two amplitudes differ substantially at lower energy: • The cusp at the ππ 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 2 π .
• The p 4 amplitude has a zero at s = 4/3m 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.

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 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 ηπ channel by a constant C a 2 ηπ defined from the Lagrangian The couplings to the γγ 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 ηπ (s) = λ ηπ (s)/4s. The experimental values of the branching fractions of the main hadronic decay modes are [5] B ηπ = (14.5 ± 1.2)%, B KK = (4.9 ± 0.8)%, and the 2γ width given by the PDG is Γ a 2 γγ = 1.00 ± 0.06 keV .
From this, one obtains the following values for the coupling constants choosing C a 2 ηπ to have a positive sign. Upon performing fits to the differential cross-sections (see below) the two couplings C a 2 γγ , D a 2 γγ will get separately determined. Defining a coupling constant C a 2 KK in the same way as C a 2 ηπ 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 Γ T (s) in the denominator, for which we use the same modelling 2 as in ref. [12] Γ T (s) = Γ ηπ (s) + Γ KK (s) + Γ ρπ (s) + Γ ω2π (s) and a Blatt-Weisskopf related function (with W 2 (x) = 9 + 3x 2 + x 4 ). Finally, the J = 2 amplitudes l 2λλ 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, 2 For two-body decay modes we take Γ12(s) = Γ 0 12 × (q12(s)/q12(s0)) 2l+1 W l (q12(s)R)/W l (q12(s0)R) with s0 = m 2 a 2 and l = 2 as in [12] while for ρπ we take l = 1 and W1 = 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. q12(s) = 0 if s < (m1 + m2) 2 .
The corresponding J = 2 (KK) 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 (KK) I=0 amplitudes in order to be able to construct the K + K − and K 0K0 amplitudes separately and compare with the experimental results. From the value of the KK 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 KK and C a 2 KK to have the same sign.

Experimental inputs
We will compare our model for the γγ amplitudes with precise experimental data on γγ → πη from the Belle collaboration [12], as was done recently in ref. [17]. In addition, we consider also here γγ → KK data in order to provide further constraints on the coupledchannel 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 [15]. Experimental data for the charged kaons channel K + K − are also available in this low energy range [39] but they are older and have much less statistics. We will restrict ourselves to the energy range E ≤ 1.4 GeV: we can use 448 differential cross-section points for πη (0.85 ≤ E ≤ 1.39 GeV) and 240 differential cross-section points for K S K S (1.105 ≤ E ≤ 1.395 GeV).

Parameters of the T -matrix
Concerning the S-wave, firstly, we employ the T -matrix model of ref. [26] which involves six parameters. It uses a chiral K-matrix type representation, which ensures one-channel unitarity below the KK threshold and two-channel unitarity above, together with the chiral expansion. This kind of approach was initiated in [40], see ref. [41] for a review. The T -matrix is written as where the matrix Φ reads it involves four phenomenological polynomial parameters α i , β i and the one-loop functions J P 1 P 2 are given in (131). 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 ηπ and KK 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 λ, The form of g 1 , g 2 is derived from a resonance chiral Lagrangian [42] such that K (6) has chiral order p 6 provided λ 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 α i , β i , λ 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 [43] 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. [44] (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 1.4 GeV. The T -matrix is computed from eq. (83) for E ≤ 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 δ 11 (∞) = 2π, δ 22 (∞) = 0 and η(∞) = 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.

Fits results
In addition to the six S-wave T -matrix parameters listed above, further parameters must be introduced which describe couplings to the γγ 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γ , D a 2 2γ , C f 2 2γ , D f 2 2γ in the fitting as well as the mass and width of the a 2 (1320) resonance: m a 2 , Γ 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. [26] (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 γγ data in this manner. Relaxing the T -matrix parameters, reasonably good fits become possible and we actually found two distinct minimums of the total χ 2 combining the πη and the K S K S data, In the case of πη the first number corresponds to the region E < 1.1 GeV. In this region fit II is better than fit I while the overall χ 2 is slightly smaller in fit I (χ 2 = 428) than in fit II (χ 2 = 439). The χ 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 [45].  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. Fig. 5 shows the behaviour of the two phase-shifts δ 11 , δ 22 and of the inelasticity parameters η as a function of energy, corresponding to the two fits. At low energies, the πη phase-shift changes sign and becomes negative. This low-energy behaviour is not anticipated in simple hadronic models of the πη amplitude (e.g. [46]) but it was also observed to emerge from fitting the γγ data in ref. [24]. In fit II the πη scattering length (defined as in [47]) is found to have the following value For comparison, the first two terms in the chiral expansion of the scattering length give m π a πη 0 = 6.17 × 10 −3 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 πη phase-shift in fit I around E 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 3 , This resonance is lighter and narrower than the standard a 0 (1450). Since the phase-shift increases from π/2 to π (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 [12] which displays a resonance (called a 0 (Y ) in that reference) which has very similar features while using a parametrisation rather different from ours 4 . 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, s a 0 = 1421(5) − i 175(4) MeV (fit II) .
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. Fig. 6 and fig. 8 show a sample of γγ → πη, K S K 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 λλ = 0 ++ and the 2 ++ amplitudes. We have seen already that the fast  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 (θ) − 1 = − √ 6 d 2 20 (θ)) and no specific 0 ++ resonance effect remains. We also note that the a 0 in fit I essentially decouples from the KK 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 πη 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 θ. This discrepancy could be caused by an isospin breaking effect close to the resonance peak. Our model assumes isospin symmetry and uses m K ≡ 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  Figure 7: Cross-sections for γγ → πη, K S K S , K + K − integrated in the range | cos θ| < 0.8. The data are from refs. [12,15,39], they are compared with the two fits results. resolution should be smaller than 2m K + − 2m K 0 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γ coupling constants, a 2 mass and width) are collected in table 3 below. The couplings D T 2γ are found, as expected, to be smaller in magnitude than the C T 2γ 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 shifted 5 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,λλ . The presence of this term is essential for obtaining a correct description of the amplitude in the η decay region s < (m η − m π ) 2 . In the 1 GeV energy region, one   Table 3: Values of the coupling constants, of the subtraction constant b k and of the mass and width of the a 2 (1320) resonance resulting from the two fits.
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.

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 C. One may define coupling constants from the residues of the resonance poles (e.g. [49] in the context of photon-photon amplitudes), the couplings to the πη and KK channels are thus defined as 16πT (II) and similarly for the third Riemann sheet. The coupling to the γγ 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 √ s a 0 = 1000.7 ( 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 2 η but its exact value is not known, we will vary it in the range s A = m 2 η ± 3m 2 π . 2) Set of the O(p 4 ) parameters L i : in the BE14 determination [44], 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. [44]).
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 mesons 6 and b) we varied the products of the couplings of the vector mesons by ±20%.

5.5
The η → π 0 γγ decay amplitude We consider now the modelled amplitudes in the kinematical region relevant for the decay η → π 0 γγ. Including the J = 0 and J = 2 partial-waves 7 , the distribution as a function of the γγ invariant mass squared, s, reads also accounting for the J = 0 isospin violating amplitudel 0,++ . Fig.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 π + π − 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 Γ = 0.237 eV which is on the low side of the most recent experimental determinations Γ exp = 0.285 ± 0.031 ± 0.061 eV (Crystal Ball at the AGS [13]), Γ exp = 0.33 ± 0.03 eV (MAMI [14]). A smaller value was reported by the KLOE-2 collaboration [51] but it has not been confirmed and the data is currently being reanalysed. The amplitudes in the η 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 2 η − 3m 2 π , m 2 η + 3m 2 π ] and comparing with the experimental results. Given 7 Contributions from J ≥ 4 partial-waves can be checked to be completely negligible.  s (GeV 2 ) s A =m 2 η +3m 2 π s A =m 2 η s A =m 2 η -3m 2 π AGS(08) MAMI(07-09) Figure 10: Experimental data on the η → π 0 γγ energy distribution [13,14] compared with predictions from our amplitudes showing the influence of the position of the Adler zero.
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 Γ η→π 0 γγ = 0.237 +0.060 −0.043 eV .

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 (πη, KK) 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 crosschannel 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 πη and K S K 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 [39].
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 πη and the KK channels is perfectly compatible with both the πη and the K S K 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. [17] in which the mass and width are found to be both significantly larger than the PDG values. Data with a smaller energy resolution would be useful to resolve these remaining ambiguities. The γγ → K + K − , K s K s cross-sections close to the KK 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. [28] and with the estimates made in ref. [52] but not with those from ref. [25]. Experimental data in this near-thresold region would obviously be very constraining. Finally, it will be quite interesting to see how the pole position determined in a lattice QCD simulation [21] evolves when the value of m π is decreased. work is supported in part by the European Union's Horizon2020 research and innovation programme (HADRON-2020) under the Grant Agreement n • 824093.
A Chiral expansion results for γγ → P 1 P 2 amplitudes A.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 π + π − and K + K − which were given in eq. (35). The leading contributions to the amplitudes which involve two neutral mesons: π 0 π 0 , K 0K0 and π 0 η appear at order p 4 . They were computed in refs. [30,[53][54][55]. The basic one-loop function which occurs in these p 4 amplitudes can be written as, for small z values, |z| << m 2 P , it behaves as The expression for the γγ → π 0 π 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 oneloop functions G P while the B amplitudes vanish at this order. Next, the expression for γγ → K 0K0 at O(p 4 ) reads The NLO expression for the γγ → π 0 η amplitude, which is of particular interest here, was worked out in ref. [30]. They considered a chiral framework which also includes the η meson as a light meson in addition to the pseudo-Goldstone octet. The relation between the η and η mesons and the octet and singlet chiral fields involves a mixing angle θ, denoting c θ = cos θ, s θ = sin θ. The expression for the γγ → π 0 η amplitude as a function of θ reads, with t θ = tan θ. The isospin conserving part in eq. (105) agrees with ref. [56]. It reproduces the result from ref. [30] upon setting s θ = −1/3, c θ = √ 8/3 and neglecting the terms proportional to t θ (which are not really negligible). In the standard ChPT one must set c θ = 1, s θ = 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 γγ → ηη Finally, the order p 4 contributions to the γγ → π + π − and K + K − amplitudes read and Combining eqs. (103) and (110) we get the NLO part of the γγ → (KK) I=1 , (KK) I=0 amplitudes A.2 Soft pion limit (isospin conserving part) Let us consider a limit when the pion in the γγ → π 0 η 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 2 π ) corrections. In this limit, firstly, one has s = m 2 η , t = u = 0 (soft pion limit) , and using standard soft pion methods one shows that the amplitude γγ → π 0 η vanishes exactly at this point. When p 1 = 0 the two tensors T µν 1 , T µν 2 become degenerate such that the soft limit amplitude, in tensorial form, reads lim p 1 =0 W µν (q 1 , q 2 ; p 1 , p 2 ) = (A(m 2 η , 0) + 2m 2 η B(m 2 η , 0)) T µν 1 = 0 .
The combination which appears in eq. (114) 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, L ++ (s A , t = u) = 0 (115) with s A = m 2 η + O(m 2 π ). Let us consider the implication of this result for the j = 0 partial wave. The expansion of the helicity amplitude with cos θ = 0 (which corresponds to t−u = 0) reads L ++ (s, cos θ = 0) = j even (2j + 1) (−1) j/2 (j − 1)!! j!! l j,++ (s) .

A.3 Soft pion limit (isospin violating part)
In the case of the isospin violating piece of the γγ → π 0 η amplitude, the soft pion limit does not lead to an Adler zero. Let us define this amplitude in terms of the following matrix element,L ++ ≡ π 0 (p 1 )η(p 2 )| 1 2 (m d − m u )(ūu −dd)|γ(q 1 , +)γ ( q 2 , +) , 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 3 5 is the axial charge The commutator is easily worked out and gives lim p 1 =0L ++ = (m d − m u ) 2F π η(p 2 )|i(ūγ 5 u +dγ 5 d)|γ(q 1 , +)γ ( q 2 , +) .
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(ūγ 5 u +dγ 5 d) can couple either to a single η meson or to ηπ + π − . Representative one-loop diagrams are shown in fig. 11. Computing them gives η(p 2 )|p 0 (0)|γ(q 1 , +)γ(q 2 , As one might expect, using this result in the soft pion relation (121) reproduces the isospin violating amplitude at chiral order p 4 when s = m 2 η (see eq. (106)), up to O(m 4 π /m 4 η ) 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 expressl 0++ (m 2 η ) in terms of L . The chiral calculation (122) 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.
B Dispersive γγ → ππ, (KK) I=0 S-waves Input for the γγ → (KK) I=0 S-wave is needed in order to reconstruct the amplitudes for the physical K + K − and K S K S states. We briefly recall the dispersive coupled-channel representation for the I = 0 S-wave taken from ref. [27]. It has the following form   111)). The three parameters b 0 , b 0 , b 0 K where then determined from fits to γγ → π + π − , π 0 π 0 data. Having done so the KK amplitude k 0 0++ was generated, which we have used in the present work together with the I = 1 amplitude. Fig. 12 shows the dispersive results for these two KK amplitudes. At low energies the amplitude k 0 0++ is significantly larger in magnitude than k 1 0++ because it gets isoscalar ππ rescattering contributions (which contain the broad σ 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 0K0 cross-section as compared to the K + K − one close to the KK threshold.

C 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,