Revisiting $\gamma^*\to\gamma\pi^0\eta$ near the $\phi(1020)$ using analyticity and the left-cut structure

Amplitudes of the form $\gamma^*(q^2)\to \gamma P_1P_2$ appear as sub-processes in the computation of the muon $g-2$. We test a proposed theoretical modelling against very precise experimental measurements by the KLOE collaboration at $q^2=m^2_\phi$. Starting from an exact parameter-free dispersive representation for the $S$-wave satisfying QCD asymptotic constraints and Low's soft photon theorem we derive, in an effective theory spirit, a two-channel Omn\`es integral representation which involves two subtraction parameters. The discontinuities along the left-hand cuts which, for timelike virtualities, extend both on the real axis and into the complex plane are saturated by the contributions from the light vector mesons. In the case of $P_1P_2=\pi\eta$, we show that a very good fit of the KLOE data can be achieved with two real parameters, using a $T$-matrix previously determined from $\gamma\gamma$ scattering data. This indicates a good compatibility between the two data sets and confirms the validity of the $T$-matrix. The resulting amplitude is also found to be compatible with the chiral soft pion theorem. Applications to the $I=1$ scalar form factors and to the $a_0(980)$ resonance complex pole are presented.


Introduction
The discrepancy between the experimental determinations of the muon g − 2 and its calculation within the Standard Model (e.g. [1]) has recently been confirmed by a new measurement [2]. Since the soft hadronic contribution to the g − 2 has the largest error in the calculation it seems necessary to further investigate the theoretical descriptions of amplitudes involving light mesons together with real or virtual photons and to also further improve our knowledge of the interactions of these mesons among themselves and the properties of the associated light resonances.
We reconsider here the description of amplitudes of the form γ * (q 2 ) → γP 1 P 2 involving one real and one virtual photon. In the context of the g − 2, the channel P 1 P 2 = π + π − is the most relevant one. It was considered in refs. [3,4], also in the case of two virtual photons, using a method which applies when the energy of the P 1 P 2 system is much smaller than 1 GeV. Recently, a generalisation was performed [5], from which the contribution of the f 0 (980) resonance to the muon g − 2 was estimated. We develop here a similar approach, in the case of one virtual photon, which can be applied in an energy region of the P 1 P 2 system slightly larger than 1 GeV. We focus on the case where P 1 P 2 = π 0 η, (KK) I=1 and will test the method against the very precise measurements performed by the KLOE collaboration close to the φ(1020) meson peak [6] (earlier results can be found in [7,8,9,10]). The amplitudes φ → γπ 0 η, φ → γπ 0 π 0 have large resonant contributions induced by the scalar mesons a 0 (980), f 0 (980) respectively. Initial interest in measuring such amplitudes was stimulated by the claim [11] that this would allow to clearly discriminate between models of these presumably exotic mesons (the q 2q2 [12] versus the KK molecule model [13,14]) but some disagreements with this claim have been expressed (e.g. [15]). We refer to [16] for a review on this topic.
Our objective here is rather to probe both the quality of the description of the radiative amplitude which can be achieved and the properties of the πη S-wave T -matrix globally in the low/medium energy range. While the existence of a sharp I = 1 scalar resonance coupling to both πη and KK channels was established long ago [17,18] the detailed behaviour of the phase shifts is still controversial. The pattern which was established in lattice QCD simulations [19] at m π = 391 MeV may or may not extend down to the physical pion mass value, depending on the extrapolation model [20]. A very similar phase shift pattern was found to emerge in the meson-meson scattering model developed in ref. [21] which was applied to γγ → πη scattering in ref. [22] and found to describe reasonably well the data by the Belle collaboration [23]. However, the mass and width properties of the a 0 (980) resonance in this model are not compatible with the PDG values. Moreover, a determination of the πη − KK scattering matrix performed in ref. [24] using the same set of γγ data, together with data on γγ → K S K S from [25] and γγ → K + K − from [26] obtained different results for the phase shifts and for the properties of the a 0 resonance.
Our approach is based on the general analyticity structure of the γ * → γπη, KK partial waves and the use of the Omnès method [27] (extended to several channels [28]) which ensures, by construction, that two-channel unitarity is exactly satisfied. In some previous work [11,29,30] the emphasis has been on a quick determination of the a 0 (980) parameters, but the amplitude models eventually violate unitarity, which may introduce a bias in this determination. The amplitude models used in refs. [31,32] based on a unitarisation of the leading order chiral expansion (uχPT) do fulfil unitarity (the relation to our formalism is worked out in appendix A) but the underlying T -matrix is possibly somewhat unrealistic (having e.g. no left-hand cuts).
We develop an Omnès-type representation for γ * → γπη, γKK S-waves which involves explicit integrations over the left-cut. When the virtuality q 2 is negative or vanishing, this left-cut has a simple structure, lying on the negative real axis, but this changes for timelike virtualities. In that case, the left-cut has several components extending in the complex plane [33]. The integrations have to be done carefully and we explain this in detail. One component of this left-cut eventually overlaps with part of the unitarity cut and turns around the threshold point. In addition the Born amplitude has a pole singularity which also lies on top of the unitarity cut. These singularities induce a deviation of the phase of the γ * → γπη partial wave from the the πη elastic phase shift such that Watson's theorem does not apply [34]. A further consequence concerns the Adler zero which, in the partial wave, is moved to an unphysical Riemann sheet.
This coupled-channel dispersive representation involves two real subtraction parameters. The presence of subtraction parameters is a necessary and unavoidable consequence of the effectivetheory nature of dispersion relations and should not be omitted, they absorb required corrections to the higher energy parts of the integrations. The minimal number of subtractions is determined from the asymptotic bounds on the partial waves and the dynamical constraints which must be satisfied, like the exact soft photon zero in the present case [35]. Because of the very small number of undetermined parameters which are involved, the amplitudes γ * → γπη, similarly to γγ → πη, are very good probes of the final-state interaction T -matrix. We will show that a rather good fit of the KLOE data [6] on φ → γπη can be achieved with two fit parameters, using a T -matrix previously determined from γγ data [24]. We will also study how combining the two sets of data can help improve the determination of the T -matrix and the properties of the a 0 (980) resonance.
The plan of the paper is as follows. After recalling some general properties related to gauge invariance we discuss the origin of the cuts in the partial waves and the approximations to be used for evaluating the discontinuities, based on vector meson exchanges in crossed channels. Starting form a general, exact, dispersive representation for the partial waves of interest we derive a coupled-channel Omnès representation valid in the considered energy region. This representation involves two subtraction parameters and a number of integrals over both the unitarity cut and the left cuts. We then discuss the evaluation of the vector meson coupling constants which are needed (including their relative signs) using experimental inputs together with flavour, chiral and asymptotic constraints. After explaining how to accurately compute all the integrals over the complex contours we perform a detailed comparison with the experimental data and discuss some consequences.

.1 Tensor and helicity amplitudes
We first recall that amplitudes involving two real or virtual photons and two pseudo-scalar mesons γ * (q 2 ) → γ * (q 1 )P 1 (p 1 )P 2 (p 2 ) can be derived from a correlation function involving two electromagnetic currents W µν (p 1 , p 2 ; q 1 , q 2 ) = i d 4 xd 4 y exp(−iq 1 x + iq 2 y) P 1 P 2 |T [j em µ (x)j em ν (y)]|0 (1) Current conservation ∂ α j em α = 0 implies that the tensor W µν must satisfy the two Ward identities One can then expand W µν over a basis of five independent tensors (e.g. [36]) satisfying (2) formed with the three independent momenta q 1 , q 2 and ∆ = p 1 − p 2 , We consider here more specifically the situation where one photon is real and the other is virtual, In this case, the tensors T 4µν , T 5µν are not physically relevant.
The dependence as a function of the virtuality q 2 is expected to display large Breit-Wigner peaks associated with the light vector resonances ρ,ω,φ, such that where T µν V ′ describes the amplitude for the vector meson V ′ to decay into γP 1 P 2 . We will consider here the situation where q 2 is close to the peak of the φ(1020) resonance i.e. q 2 ≃ m 2 φ , and focus on the construction of the amplitude T µν φ , assuming that the background contributions from the the other vector resonances can be neglected or have been subtracted. In the sequel, the index φ will be dropped. This amplitude can be expanded in terms of the independent tensors as where the functions A, B, C are Lorentz scalars depending on the external masses and on the Mandelstam invariants s, t, u, which satisfy s + t + u = q 2 + m 2 1 + m 2 2 . We can derive the helicity amplitudes by contracting the tensorial amplitude with the polarisation vectors of the photon and the φ meson (with a conventional minus sign). The scattering angle θ is defined as the angle between p 1 and q 1 in the centre-of-mass system of the two pseudoscalar mesons (see fig. 1). The Mandelstam variables t, u are expressed as follows as a function of θ where λ 12 (s) = s 2 − 2sΣ 12 + ∆ 2 12 (10) The three independent helicity amplitudes, finally, can be written in terms of the three functions A, B, C as γ( q 1 ) γ * ( q 2 ) Figure 1: Centre-of-mass system for the γ * (q 2 ) → γ(q 1 )P 1 (p 1 )P 2 (p 2 ) amplitude.

Partial waves and their singularities
We will designate the helicity amplitudes of interest as The isospin I = 0, 1 KK amplitudes are related to the K + K − , K 0K 0 ones by and we consider here the coupled-channel system of the two S-waves with I = 1 The angular integration can be expressed in terms of the t-variable, e.g. for the πη amplitude Let us first examine the singularities of the partial wave when s → 0 and s → (m η ± m π ) 2 : 1) Singularity when s → 0: When s → 0, the lower integration boundary in eq. (16) goes to ∞ since one has which leads to a singular behaviour of l 0++ (s) (see [37]). The behaviour of the amplitude L ++ (s, t) in the regime when s → 0 and t → ∞ is driven by the a 2 Regge trajectory which gives, with α 0 ≃ 0.5, α ′ ≃ 0.9 GeV −1 . Using this in eq. (16) one obtains for the S-wave near 2) Singularities when s → (m η ± m π ) 2 : When performing the t-integration (16) with a fixed value of the energy s one must pay attention to the cuts of the function L ++ (s, t) in the t variable. In the t-channel, which we can write φη → γπ the lightest unitarity contributions are from isoscalar 3π states which generates a cut in the t variable lying on the real axis from t = t 0 = 9m 2 π to t = ∞. Similarly, in the u channel φπ → γη the lightest unitarity contributions are from isovector ππ states. This generates a cut in the t variable extending from −∞ to t(u 0 , s) = q 2 + m 2 η + m 2 π − u 0 − s with u 0 = 4m 2 π . This cut can be shifted away from the real axis by appending an infinitesimal imaginary part to q 2 When s is in the physical region, λ πη (s) is real and one easily verifies that The integration path in eq. (16) therefore lies in between the two cuts without touching them for finite values of s. The situation changes when s is in the unphysical region (m η − m π ) 2 < s < (m η + m π ) 2 where λ πη (s) is imaginary. The path of integration must then be distorted to turn around the cut as illustrated in fig. 2. As a consequence, the partial wave l 0++ (s) diverges when s approaches the threshold from below since the integral in eq. (16) remains finite while λ πη (s) in the denominator vanishes. A similar divergence occurs when s approaches the pseudo-threshold (m η − m π ) 2 from above. The discontinuities of l 0++ (s) when s moves across the points (m η ± m π ) 2 reflect the fact that s crosses the cut of the partial-wave amplitude at these points. This cut will be described in more detail below.
The singularities of the partial waves when s crosses the points (m η ± m π ) 2 affects the position of the Adler zero. The existence of an Adler zero in the chiral limit when the pion is soft, i.e. p 1 = 0, implies that the physical helicity amplitude L ++ (s, t) should have an Adler of s in the range (m η − m π ) 2 < s < (m η + m π ) 2 (here s = 0.25 GeV 2 , q 2 = m 2 φ ) and of the integration path for performing the partial-wave projection.
The Adler zero is also present in the J = 0 partial wave for small enough q 2 virtualities 1 . When q 2 increases the zero is no longer present in the first Riemann sheet.

Partial-wave dispersion relations
One can estimate the asymptotic behaviour of the amplitudes L λλ ′ (s, t) or K λλ ′ (s, t) when s → ∞ while t (or u) remains finite, based on Regge theory. It is also possible to make estimates in the regime where all three Mandelstam variables go to infinity based on quark counting rules [38]. These arguments suggest that the J = 0 partial waves l 0++ (s) and k 0++ (s) should remain bounded by √ s asymptotically and should thus satisfy once-subtracted dispersion relations.
where k 1,Born 0++ is the I = 1, J = 0 projection of the K ± pole contribution to K c ++ (s, t). As we will review below (in sec. 3.4) the S-wave amplitude k 1,Born 0++ has the following form where, for real values of s, the function I K (s) is given by The function I K is an analytic function of s with a cut along the negative real axis, The amplitude k 1,Born 0++ (s) is thus also an analytic function of s except for a cut on the negative axis and a pole at s = q 2 (which corresponds to the soft photon limit). It satisfies the following dispersive representation with β(q 2 ) = 8m 2 The φ → γKK amplitudes are symmetric in the t, u variables due to C-invariance. Singularities in the crossed channels φK → γK, apart from the charged kaon pole φK + → K + → γK + , are induced by unitarity contributions from I = 1/2 states: πK, ππK, · · · . We will approximate these by the K * (892) resonance contribution: φK → K * → γK. These cross-channel singularities generate a cut in the partial-wave amplitude, denoted as C KK illustrated in fig. 5 below. This cut includes the negative real axis and also has a complex component.
The amplitudes φ → γπη are not symmetric in the variables t,u. This leads to two different cuts in the partial waves which we will denote as C πη and C ηπ (see fig. 5). The cut C πη corresponds to unitarity contributions in the t-channel φη → γπ which has isospin I = 0 i.e. they are generated from 3π, 5π, KK, . . . intermediate states. We will approximate these by the ω and φ resonance contributions: φη → ω, φ → πγ. In the u-channel φπ → γη the unitarity contributions have I = 1 states ( 2π, 4π, KK, . . .) and will be approximated by the ρ-meson contribution: φπ → ρ → ηγ. In addition to these left-hand (and complex) cuts the partial wave amplitudes have a unitarity right-hand cut on the positive real axis s > (m π + m η ) 2 .
On must also check for the possible presence of an anomalous threshold. Anomalous thresholds occur when one of the endpoints of the left-hand cut requires a deformation of the unitarity cut [39] (see [40] for a review). This happens for amplitudes involving two virtual photons (e.g. γ * γ * → ππ, see [41,4]). In our case, only one photon is virtual. With the t, u-channel exchanges which we have considered and the two-body channels included in the unitarity relation we have checked that no anomalous threshold gets generated when varying the virtuality q 2 .
Finally, taking the soft-photon constraints and the asymptotic bounds into account one can write the following dispersive integral representations for the J = 0 partial waves l 0++ , k 1 0++ , and The discontinuity along the real axis is defined as Along the complex contour, the corresponding definition is given in terms of the parametric representation in eq. (105) below. The dispersive integral representations given in eqs. (29) (30) can be considered as exact. In practice, one intends to estimate the partial-wave amplitudes in the energy region relevant for the φ decay i.e. (m η + m π ) 2 ≤ s ≤ m 2 φ . The dispersive equations (29) (30) must be used in an effective theory sense: one tries to perform accurate approximations to the integrands in a limited energy region e.g. |s ′ | < ∼ 1.5 GeV 2 while absorbing the remaining higher energy contributions into a finite number of low energy parameters. In this work, we implement the following approximations: a) Along the left-hand and complex cuts we approximate the discontinuities by those of the light vector mesons (ρ, ω, φ, K * ) and by the K + exchange contributions to the amplitudes, b) Along the right-hand cut we evaluate the discontinuity from the unitarity relations, including two channels disc l 0++ (s) where T (s) is the two-channel πη − KK S-wave T -matrix with and Using these approximations for the discontinuities in eqs. (29) (30) generates a closed set of coupled-channel Muskhelishvili equations [42]. The solutions of these equations will include subtraction parameters which account for the contributions of the higher energy regions in the dispersive integrals where the approximations made above are no longer valid.

Coupled-channel Omnès-Muskhelishvili representation
The fundamental ingredient for expressing the solutions of these equations is the Ω matrix [42]. The matrix elements Ω ij (s) are defined to be real-analytic functions of s, having a cut on the positive real axis s ∈ [m 2 + , ∞] and the discontinuity across this cut being given by Im In order to write dispersion relations for the matrix elements Ω ij (s) we make the assumption that Ω ij (s) ∼ 1/s when s → ∞. Other choices can obviously be made 2 but this one is convenient because it reproduces the asymptotic behaviour of scalar form-factors in QCD. The I = 1 form factors are then related to the matrix elements Ω ij by simple linear relations (see [44] and ref. [45] for the analogous I = 0 case). Writing unsubtracted dispersion representations for the Ω ij functions provides a relation between Ω(s) and the S-wave T -matrix in the form of a singular integral matrix equation In general, this equation must be solved numerically but the determinant of Ω can be expressed analytically in terms of the two phase shifts δ πη , δ KK , The phases are defined as continuous functions and the following asymptotic condition is imposed where N is the so-called Noether index [46,42]. Taking N = 2 ensures that eq. (38) has a unique solution with Ω(0) = 1. One also sees from eq. (39) that the determinant of the Omnès matrix does not vanish (except at infinity) such that the inverse of Ω is defined for all s. We will use the (conventional [47]) notation for this inverse As a simple consequence of the unitarity relations satisfied by the amplitudes l 0++ (s), k 1 0++ (s) and by the matrix Ω(s) eqs. (34), (37), the two functions obtained by multiplying the amplitudes by the inverse of the Omnès matrix, have no discontinuity across [m 2 Based on this property, one can write various types of Omnès representations involving either only right-cut integrations [27] or only left-cut integrations [48] or both, which are exactly equivalent. We consider here a representation starting from the following two functions, The terms subtracted from k 1 0++ (s) remove the soft-photon singularity at s = q 2 such that the functions φ 1 , φ 2 vanish at the soft photon point. These functions have both left and right-hand cuts.
When s → ∞ the matrix elements D ij (s) ∼ s while the amplitudes l 0++ (s), k 1 0++ (s) are expected to grow no faster than √ s (see [24]). One can thus express φ 1 (s), φ 2 (s) as twicesubtracted dispersion relations. It is natural to take s = q 2 as one of the subtraction points while the second one, s 0 , can be chosen arbitrarily. A variation of s 0 is compensated by a variation of the subtraction constants a 1 , a 2 . It is convenient to choose s 0 to lie on the real axis and away from the cuts, in the subthreshold region i.e. in the range [m 2 − , m 2 + ]. In this case we can expect the subtraction constants to be essentially real. The dispersive representations of φ 1 , φ 2 involve both left-cut and right-cut integrations. They can be written as follows, taking the soft photon constraints into account where a 1 , a 2 are subtraction constants and the integrals along the cuts have the following form Let us now express these integrals in more detail.

1) Right-cut integrals:
The discontinuities of the functions φ i along the right-hand cut are driven by the singular parts of the Born amplitude. Their expressions are deduced from (43) using (42) disc Using the analyticity properties of the matrix elements D ij (s) these right-cut integrals I RC i (s, q 2 ) can be calculated explicitly and one obtains with 2) Left-cut integrals: Along the left-hand cuts now, the discontinuities of the functions φ i (s) are proportional to those of the amplitudes l 0++ , k 1 0++ . Using the model defined in eq. (33), one can write the integrals I LC i in eq. (43) as a sum over five terms The integrals I LC,Born i are associated with the left-cut discontinuity induced by the kaon pole, they read One remarks that the double pole at s ′ = q 2 in the integrand does not cause any numerical difficulty when q 2 > 0 and that the integrals are real when s is in the physical region (also assuming that the subtraction point s 0 is real and positive).
The transitions φ → ωη, φ → φη or φ → K * K ,K * K are all kinematically forbidden. In those cases, ignoring the width of the vector resonances is a good approximation. In the zero-width limit, the expressions for the integrals I LC,V i appearing in (49) are as follows, using the results from the next section on the discontinuities of the vector-exchange amplitudes, where and sgn(z) is a sign factor, sgn(z) = ±1 (see sec. 4.2). The parameters α ω , α φ , α K * are related to resonance chiral Lagrangian coupling constants introduced in the next section, In the case of the ρ meson, the effect of the width must be taken into account because the transition φ → ρπ is kinematically allowed. This can be done in a simple way, consistently with the analyticity properties, by replacing the u-channel pole 1/(u − m 2 ρ ) by a Källen-Lehmann dispersive representation [49,50]. Denoting the zero-width limit of the ρ integrals as I LC,ρ i (s, q 2 ; m ρ ) with where one can express the finite-width result as an integral involving the corresponding Källen-Lehmann spectral function ρ V (µ 2 ; m ρ , Γ ρ ), We will show in sec. 4.3 how to simplify this expression using integration by parts.
The radiative decay amplitudes V → γP 1 P 2 have often been expressed in a form which involves the charged kaon (or charged pion) one-loop triangle function(e.g. [11,51,52,53,54,55,31,56]), which seems different from the Omnès representations presented above. We show in appendix A that expressions in terms of the triangle loop function can indeed be derived from our dispersive representation for certain specific models of T -matrices which, in particular, have no left-hand cut. Physically, of course, all the T -matrix elements for πη − KK scattering do have left-hand cuts (and a complex cut in the case of πη → πη) such that the Kaon-loop representation is not exactly valid.

Born and vector-exchange amplitudes
In this section we review the evaluation of the Born amplitude and that of the vector-exchange amplitudes based on a resonance chiral Lagrangian. The Born amplitude is proportional to the φ → K + K − coupling constant, the magnitude of which can be evaluated from experiment. The vector meson exchange contributions to φ → γKK and φ → γπ 0 η involve the product of a radiative decay V → γP coupling, the magnitude of which can be determined from experiment, and a hadronic V 1 → V 2 P coupling. Among those, only the ω → ρπ and φ → ρπ couplings can be estimated from experiment. We can use flavour symmetry and large N c arguments in order to derive estimates for the other V 1 → V 2 P couplings which are needed. This will also allow us to relate the sign of the Born amplitude to those of the vector exchange amplitudes.

Resonance chiral Lagrangian and mixing angle
We start from the following resonance chiral Lagrangian [57] We differ from ref. [57] in considering a nonet (rather than an octet) of vector mesons, which are encoded in a 3 × 3 matrix with The Lagrangian (57) is of leading order in the large N c and in the chiral expansions i.e. O(N c ), O(m 0 q ) . We consider only the couplings which are relevant to the amplitudes of interest here. In the terms proportional to ǫ µναβ we use the same notation as [58]. The pseudoscalar fields P a are encoded, as usual, in a unitary matrix U = exp(iλ a P a /F π ) and u µ ≡ u † D µ U u † with u = √ U . The external vector and axial-vector sources can be set to v µ = eQA µ , a µ = 0, so that where A µ is the photon field. With these definitions we can express the matrix element of the electromagnetic current, In order to get the correct pattern for the masses of the vector mesons we must add mass terms which break both the flavour and the nonet symmetries, where M q is the quark mass matrix and the parameter ǫ V is sub-leading in the large N c expansion i.e. O(N 0 c ). These flavour and nonet symmetry-breaking terms induce a mixing between the ω and the φ mesons, with a mixing angle ϕ V , such that the λ matrices attached to the physical ω and φ mesons can be written as where ϕ I = arctan(1/ √ 2) ≃ 35.26 • is the "ideal" mixing angle. We can first derive two relations between the Lagrangian parameters and the ρ and K * meson masses and then, the masses of the ω and φ mesons are obtained by diagonalising the singlet-octet mass matrix. After a rotation by the ideal mixing angle this matrix reads In this form, it is easy to see that the experimental fact that m ρ ≃ m ω implies that ǫ V << m 2 ρ , m 2 K * . To leading order in ǫ V , one obtains for the angle ϕ 1 This formula shows that ϕ 1 , even though numerically small, is actually chirally enhanced because the denominator is O(m s − m ud ). This provides some justification for neglecting further Lagrangian terms which break the nonet symmetry. The masses of the ω and φ mesons (i.e. the eigenvalues ofM) are given by to leading order in ǫ V . This implies the mass relation which is verified within 4%. This small discrepancy, however, gives rise to a significant relative uncertainty in the evaluation of ϕ 1 . Using the ω mass gives ϕ 1 = −1.2 • while using the φ mass gives ϕ 1 = −7.3 • .

Signs of the coupling constants
One can derive useful information on the signs of the coupling constants by using asymptotic constraints on the V → P γ * form factors. Let us define the form-factor F V P as follows where q = p V − p P . We have factored out the coupling C V P γ such that the form factor must satisfy F V P (0) = 1, which ensures that the amplitude V → P γ * becomes equal to V → P γ when q 2 = 0. Computing the form factor in the flavour symmetry limit from the Lagrangian (57) one obtains [58] Requiring that the form factor goes to zero when q 2 goes to infinity gives Combining this result with the relation which can be derived from a similar asymptotic constraint on the γπ axial form factor [57] one obtains This relation fixes the relative signs between the vector-exchange amplitudes which are proportional to σ V h V and the Born amplitude which is proportional to g V . Without loss of generality we can choose g V , h V and σ V to be positive.

V P γ and V V P coupling constants
We will need a set of V P γ and V V P couplings which can be defined from the Lagrangians The magnitude of the V P γ couplings which are needed can all be determined from the experimental values of the radiative decay widths, The results are collected in table 1 which also shows the expressions of these couplings using the resonance chiral chiral Lagrangian (57) with the leading order nonet symmetry breaking . terms (61). For the amplitudes which involve an η meson we use a simple η−η ′ mixing description such that the λ matrix attached to the η meson is and take ϕ P = −20 • . These expressions allow us to determine the signs of the C V P γ couplings as shown in the table.
Let us now consider the g V V P couplings. As is well known, one can estimate the two couplings g ωρπ , g φρπ from the experimental values of the ω, φ → 3π decay widths. The V → 3π decay amplitude derived from the resonance chiral Lagrangian (57) can be written in the following form with where P ρ (z) is the ρ-meson propagator taking the width into account (see sec. 4.2). Imposing the following asymptotic constraint 3 on F and using that P ρ (z) ∼ 1/z when z → ∞ yields the following relation which determines the value of The amplitude (76) then becomes identical to that of the original GSW model [60] and to the one obtained from effective Lagrangians implementing a hidden-gauge symmetry [61]. From the expression of the V → π + π − π 0 width which reads with one obtains the following values for the g V ρπ couplings The value of g ωρπ is compatible with the result derived from the experimental measurements of the form factor F ωπ (s) in the region s > (m ω + m π ) 2 in ref [62]. The errors quoted in eq. (82) do not contain the uncertainties induced by the modelling. Varying θ V from the value given by eq. (79) by 20% induces a variation of g ωρπ by 14% and a variation of g φρπ by 8%. It must also be kept in mind that these simple modellings of the V → 3π amplitudes and of the F V π form factors do not correctly account for the unitarity relations and the related dispersive representations [63,64,65,66]. We finally ascribe a 20% uncertainty to the values of g ωρπ and g φρπ . Using these two inputs, together with the C V P γ couplings as shown in the table one can determine the values of the parameters σ V , h V and θ 1 from a least-squares fit. This yields One can then deduce estimates for the values of the couplings g φK * K , g φωη and g φφη which are needed in order to evaluate all the relevant vector-exchange amplitudes, The corresponding numerical values of the effective couplings α V which appear in front of the I = 1 partial-wave amplitudes (see (53) (55)) are given by

Born amplitude
The Born amplitude in φ → γK + K − (see fig. 3) can now be computed from the Lagrangian (57). Including the contributions proportional to the two couplings g V and f V gives with This contribution is suppressed because of the approximate relation f V ≃ 2g V . Actually, our Omnès representations use only the singular parts of the tree amplitudes, other contributions being absorbed into the subtraction constants. Terms like (86) which make a constant contribution to the S-wave can thus be omitted. We can replace g V in eq. (85) using the relation with the φK + K − coupling The coupling g φK + K − can be determined from the φ → K + K − decay width, which gives, Figure 4: Illustration of the vector meson exchange contributions to φ → γKK and φ → γπη amplitudes considered here.
recalling that g V was chosen to have a positive sign. The three independent Born helicity amplitudes deduced from eq. (85) read The J = 0 partial-wave amplitude is then easily obtained where I K + (s) was given in eq. (25).

Vector-exchange amplitudes
We can now express the vector-exchange amplitudes (see fig. 4) in terms of the couplings listed in table 1. These amplitudes involve the following two combinations of the three independent tensors T µν The K * exchange contributions to φ → γK + K − , γK 0K 0 read and the contributions from ω, φ and ρ exchanges to φ → γπ 0 η read These expressions are derived in the zero-width limit. Finite width effects will be discussed in sec. 4.2.
The helicity amplitudes corresponding to eqs. (93), (94) are deduced using the relations (12) and the partial-wave amplitudes are then easily computed. The J = 0 partial-wave projections can be expressed in terms of the following generic function F V with and The expressions of the φ → γKK and φ → γπη S-wave amplitudes associated with vector meson exchanges are given below in terms of the F V functions, and 4 Integrations along the complex cuts

Parametric representations of the cuts
It is easy to check that the function F V vanishes in the limit s → q 2 , as it should. Concerning the singularities, F V has a pole at s = 0 when the two masses m 1 , m 2 are different. The function λ 12 (s) has square-root singularities at s = (m 1 ± m 2 ) 2 but these are not present in F V because it is an even function of λ 12 . The remaining singularity of F V is a cut generated by the log functions. This cut is given by the values of s such that the argument of one of the log's is a negative real number, that is, Solving for s as a function of t provides the following parametric representation of the cut, An important point is that the shape of the cut does not depend on the mass of the exchanged vector meson. Only the lower bound on the variable t depends on this mass, see (100). The cut of the φ → γπη amplitudes has two components, denoted as C πη and C ηπ corresponding to eq. (101) with m 1 , m 2 = m π , m η and m η , m π respectively. The left-cut of the φ → γKK amplitudes, denoted as C KK is given by eq. (101) with m 1 = m 2 = m K . The cuts lie on the real axis when t < (q − m 2 ) 2 or t > (q + m 2 ) 2 otherwise they are complex. One sees that s + goes to −∞ when t → ∞. The s − function can be re-written as which shows that s − goes to 0 when t → ∞. We also note that in the kinematical configuration where q ≥ m V + m 2 the endpoints of the cut s ± (m V , q 2 , m 2 1 , m 2 2 ) are located on top of the unitarity cut. This corresponds to the triangle diagram singularities interpreted in ref. [67].

Testing the dispersive representations on the vector-exchange amplitudes
We first consider the J = 0 partial-wave projections of the vector-exchange amplitudes and discuss their dispersive representations as integrals over the complex cuts. This allows us to check the correctness and the numerical accuracy of such integrals before extending them to the complete Omnès representation of the S-wave. In the figures, the lower bound on the parameter t in the representation (101) has been taken as t 0 = (m K +m π ) 2 for the K * exchange, t 0 = 9m 2 π for ω exchange and t 0 = 4m 2 π for ρ exchange. The effect of the q 2 + iǫ prescription is illustrated in the right bottom plot.
The generic function F V which appears in the expressions of the S-waves (see (95)) must satisfy the following Cauchy representation as an integral over the cut C 12 , with taking into account the asymptotic behaviour, the pole 4 at s = 0 and the cut. The discontinuity across the cut is defined as follows (whereṡ ± (t) ≡ ds ± (t)/dt). Using the explicit form of F V (eq. (95)) the discontinuity reads where sgn(z) = ±1, depending on which one of the two log functions in L V generates the discontinuity. The integrations along the cut are expressed as follows, using the parametric representations Fig . 5 illustrates the three different cuts C πη , C ηπ , C KK which are involved and shows the regions where the sign factor in the discontinuity (106) sgn(z) is +1 or −1 in different colour. For C KK the sign change occurs when the parameter t = q 2 − m 2 K and for C πη when t = q 2 − m 2 η . Along the cut C ηπ there are several sign changes which are indicated in table 2. Table 2: Values of the sign factor along the contour C ηπ as a function of the parameter t, with t 1 = m η (q 2 /(m η + m π ) − m π ), t 2 = m η (q 2 /(m η − m π ) + m π ) .
These three cuts include the negative real axis [−∞, 0], as we have seen, and a complex component. In addition, the πη cuts extend along the positive real axis. These cuts can be shifted away from the unitarity cut by appending an infinitesimally small positive imaginary part to q 2 . Nevertheless, the cut C ηπ crosses the real axis at the two points m 2 ± = (m η ± m π ) 2 . When the energy s crosses one of these points the discontinuity of the πη partial wave exhibits a divergence induced by the factor 1/λ ηπ (s). This is in accordance with the discussion in sec. 2.2 on the path of integration for performing the partial-wave projection when s is in the range m 2 − < s < m 2 + . Vector mesons are resonances and thus have a finite width. The width can be simply taken into account, in a way compatible with analyticity properties, by replacing the pole in the t or u variable in the tensorial vector-exchange amplitudes (93) (94) by a dispersive Källen-Lehmann [49,50] representation where t 0 = 4m 2 π in the case of the ρ. We will assume that the spectral function ρ V satisfies the normalisation condition ∞ t 0 which ensures that the propagator goes as 1/u when u goes to infinity. In the case of the ρ meson, for instance, we will use the following model for the spectral function 5 where γ V is proportional to the physical width Γ ρ and N V is adjusted such that eq. (109) is satisfied. It is easy to see that the model (110) satisfies the condition that ρ V should tend to a δ function when the width Γ V goes to zero, Using the finite width propagators (108), the vector-exchange partial-wave amplitudes, get expressed as integrals over ρ V , e.g.
The corresponding dispersive representations have the same form as eq. (103) in which the functions I 12 (s) are expressed as double integrals, (114) These expressions can be simplified using integration by parts. Let us introduce the following two integrals of ρ V , which can be computed analytically using the model (110) (see appendix B). Integrating by parts, the functions I 12 get expressed in the following way in terms of ρ V 1 and ρ V 2 Figure 6: Deformed contour used for the numerical computation of the functionsΦ (n) which involves a single integration. The two terms involving Φ (n) (s, q 2 ) are boundary contributions which will be detailed below. A remark is in order here: one sees from eq. (101) that s + (t) − q 2 , s − (t) − q 2 both vanish when t = m 2 1 . This generates a potentially highly singular term in the integral of eq. (116). This problem only affects the contour C ηπ since, in this case, m 2 1 = m 2 η which lies within the range of integration. The functions ρ V 1 (t), ρ V 2 (t) were chosen such as to vanish when t = m 2 1 (see (115)) and this removes completely these singularities in the t-integral (116). The singular integrations are now contained in the two functions Φ (1) They can be evaluated analytically in terms of the function L V (s, q 2 , m 2 (118)

Omnès integrations with a finite resonance width
Let us now return to the Omnès representations and reconsider the left-cut integrals I LC,ρ i , which were given in eq. (56) in the form of double integrals. The discussion in the preceding section can be easily adapted to derive expressions in the form of single integrals which are easy and fast to compute. One obtains, where s ± (t) parametrise the contour C ηπ and the functionsΦ These integrals must be computed numerically. When s is in the physical range of the φ decay, (m η + m π ) 2 ≤ s ≤ q 2 , the numerical integration can be performed in a way which avoids the singularities at z = q 2 and z = s by deforming the contour. Taking into account the exact cancellations in the regions of overlapping branches which have a different sign factor (see fig. 5) the contour can be deformed into a) the segment [s − (t 0 ), s + (t 0 )] on the real axis (note that q 2 < s − (t 0 ) ≃ 1.266 GeV 2 ), b) the four segments starting from (m η ± m π ) 2 and joining s ± (q 2 − m 2 π ) and c) the negative real axis. This deformed contour, showing the directions of integrations, is illustrated in fig. 6.

Recovering the Adler zero
The component C ηπ of the complex cut in the l 0++ partial wave which is induced by the ρ meson exchange (and more generally by u-channel exchanges) is responsible for the disappearance of the Adler zero, as was mentioned above. Indeed, this cut crosses the real axis at the two points s = m 2 + and s = m 2 − and this causes a divergence of the partial wave when s approaches m 2 + from below or m 2 − from above. It is possible to recover the Adler zero by considering an analytical continuation of the amplitude to an unphysical Riemann sheet which we will call the B Riemann sheet. Let us illustrate this in the case of the simple ρ-exchange amplitude l ρ 0++ (s). The B Riemann sheet extension l ρ,B 0++ (s) is defined such as to match continuously with l ρ 0++ (s) upon crossing the portion of the cut C ηπ which extends from m 2 + to s + (t 0 ). In the zero-width limit the two functions are simply related by This formula is easily generalised to the case of a finite width using the integral representation of the amplitude and one finds where theρ V 1 (t),ρ V 2 (t) are integrals of the spectral function which vanish at t = t 0 (see (160)). These ρ-exchange partial-wave amplitudes are illustrated in fig. 7 which shows that the B-sheet extension displays an Adler zero close to s = m 2 η .

Comparison with experiment
Given a T -matrix satisfying two-channel unitarity and appropriate asymptotic conditions the Ω-matrix (and its inverse D) can be determined numerically in a unique way. It is then only necessary to compute the left-cut integrals which appear in eqs. (44) and one obtains an expression for the two S-wave φ decay amplitudes l 0++ , k 1 0++ in terms of two parameters a 1 , a 2 . An important check of the numerical implementation (which we have performed) is to verify that the discontinuities of the obtained solutions for l 0++ (s), k 1 0++ (s) across the unitarity cut are exactly given by the unitarity equations (34).
In the numerical results presented below we use the two-channel T -matrix model introduced in ref. [44]. In this model, two-channel unitarity is implemented using a K-matrix approach. The form of this K-matrix is recalled in appendix C: it is constrained such that the low-energy expansion of the corresponding T -matrix matches with the chiral expansion (χPT) up to NLO and it involves six phenomenological parameters. Determinations of these parameters were obtained in ref. [24] by performing fits to experimental measurements of γγ → πη, γγ → K S K S differential cross sections in an energy range E < ∼ 1.4 GeV (a set of 688 data points was used).
These fits are based on Omnès representations of the γγ S-wave amplitudes, exactly analogous to eqs. (43), (44) which depend on two subtraction parameters in addition to the T -matrix parameters. Two different sets of T -matrix parameters were actually found in the γγ fits, giving approximately similar χ 2 minimum values. One of these sets gives rise to a "narrow" a 0 (1450) resonance and has some unphysical features, it will not be considered here. Let us now re-write the expressions of the S-wave amplitudes in the Omnès-dispersive approach l 0++ (s)

Left-cut and right-cut integrals
in which we have definedk 1 0++ (s) to be the part of the the I = 1 KK amplitude which vanishes at the soft photon point i.e. to be dominant over the left-cut ones since they are induced by the diverging part of the Born amplitude. This turns out to be true in the case of the I 1 integrals as can be seen from fig. 8 (top). In the case of the I 2 integrals, however, one sees from the figure that I LC,K * 2 is larger than the real part of I RC 2 .

Detailed comparison with KLOE results
In the physical decay region, the φ → γπ 0 η helicity amplitudes are dominated by the S-wave.
The J ≥ 1 amplitudes are approximated by the tree-level vector-exchange contributions. We have also made an estimate of the J = 2 amplitude induced by the a 2 (1320) tensor meson exchange in the s-channel (φ → γa 2 (1320) → γπη, see appendix D). This contribution turns out to be essentially negligible in the physical region. Using these amplitudes, the differential decay width of the φ → γπ 0 η mode as a function of the πη invariant mass squared can be expressed as follows in terms of the three independent helicity amplitudes We will compare the results from this theoretical model with the experimental measurements of the single differential decay width performed by the KLOE collaboration [6]. They are presented in table 5 of ref. [6] in the form of 49 data points which have been corrected from the background and the energy resolution effects. We will not include the last two points in the χ 2 (see the discussion below). The χ 2 is then computed from the following formula, where σ i is the error on the i th data point, ∆E = 6.35 MeV is the size of the energy bins, Γ φ is the total width of the φ (we take Γ φ = 4.26 MeV and the value of the photon virtuality q = 1020 MeV as in ref. [6]). The determination of the T -matrix parameters performed in ref. [24] is based on a dispersive Omnès description of the γγ → πη, KK amplitudes involving two subtractions parameters. One of these was fixed by assuming a given position of the Adler zero s γγ A in the γγ → πη S-wave amplitude. Using the T -matrix which corresponds to s γγ A = m 2 η in the computation of the φ decay amplitudes and fitting the two parameters a 1 , a 2 to the KLOE data on gets χ 2 (φ)/N dof = 58/45. Taking a slightly larger value for the γγ amplitude Adler zero: 6 s γγ A = m 2 η +3m 2 π gives a T -matrix which provides an even better fit to the φ decay data with χ 2 (φ)/N dof = 40/45 and the values of a 1 , a 2 are a 1 = 0.24 ± 0.01, a 2 = −1.74 ± 0.03 (GeV −1 ).  Figure 9: Normalised differential energy distribution of the φ → γπ 0 η decay width. The dashed blue curve is the result from the dispersive S-wave amplitude and the solid red curve is the result including the J ≥ 1 partial waves. The experimental points are results from the KLOE collaboration: black squares are from [6], brown filled circles and green empty circles are from ref. [7] (the η being detected via the η → 3π and the η → 2γ modes respectively).
The ability to reproduce the φ decay data in the whole physical energy region with only two fit parameters 7 indicates a good compatibility between these data and the γγ scattering data concerning the final-state rescattering T -matrix.
The differential decay width corresponding to this fit is displayed in fig. 9 showing separately the contribution of the S-wave. One sees a small but visible contribution from the higher waves in the energy region below 900 MeV. This contribution is dominated by the J = 1 amplitude which is allowed for γ * → γπη when q 2 = 0. The two S-wave amplitudes l 0++ andk 1 0++ obtained from this fit are plotted in fig. 10 and compared with the tree-level vector-exchange amplitudes. The amplitudek 1 0++ is seen to qualitatively agree in size with the sum of the K * exchange and the regular part of the Born amplitude for energies below 0.8 GeV 2 . In contrast, the behaviour of the πη amplitude l 0++ is dominated by the rescattering mechanism even at low energy. In the energy region below the threshold s ≤ m 2 + the B Riemann sheet extension l (B) 0++ is shown (see sec. 4.4), which is expected to display an Adler zero s φ A = m 2 η + O(m 2 π ). The figure shows that this is indeed the case. The central value of this zero is located at s φ A = m 2 η + 6.2m 2 π and its position varies in the range [m 2 η + 3.7m 2 π , m 2 η + 8.1m 2 π ] when varying the input values of the coupling constants which control the left-hand cut discontinuities. These features indicate that the fitted values of the two parameters a 1 , a 2 have a physically reasonable size.  Below the πη threshold the extension to the B-Riemann sheet is shown (see sec 4.4). The figure on the right displays the result for thek 1 0++ amplitude. Also shown for comparison are the tree-level vectorexchange amplitudes and the singular part of the Born amplitude. The cyan band shows the variation of l 0++ when the coupling constants C V P γ and g φV P are varied.

Combined γγ and φ decay fits
We have seen that the T -matrix determined from γγ data is compatible with the φ decay data from KLOE. It is interesting to study how the T -matrix parameters would be modified if the two data sets were combined. We have performed such a fit in which we have increased the weight of the φ decay data by a factor of two. In this way, the following total χ 2 values are obtained, showing a significant improvement in the φ data χ 2 . At the same time the χ 2 value of the γγ data remains acceptable. For comparison, the fit including only γγ data (with s γγ A = m 2 η + 3m 2 π ) gives χ 2 [φ] = 39.8, χ 2 [γγ] = 421.8 (γγ fit) .
The φ differential decay width in the two fits differ essentially in the higher energy region E πη > ∼ 940 MeV. This is illustrated in fig. 11 which also shows the location of the K + K − , K 0K 0 thresholds. One observes, in particular, a significant improvement on the last two points below the K + K − threshold (points 46 and 47) in the combined fit. However, the two points which are above (points 48 and 49) are not well described in either one of the fits. One observes that the energy of the point 48 is located in between the K + K − and the K 0K 0 thresholds. It is possible that the discrepancy could be explained by isospin breaking effects which were predicted to be enhanced in this energy region [69]. In our model, isospin symmetry is assumed, and we have taken m K = m K + for the kaon mass in order to have the correct mass in the Born amplitude. Because of the importance of the Born amplitude in the rescattering dynamics it seems plausible that the cusp at the K + K − threshold should be significantly more pronounced than the one  Figure 11: Results for dΓ φ→γπη /dE πη when the T matrix is determined from γγ data only (solid red curve) and when it is determined by combining the γγ and the φ data (dotted blue line).
at the K 0K 0 threshold. A complete treatment of isospin breaking effects is difficult since they induce couplings between the I = 1 channels considered here and I = 0 as well as I = 2 channels. One must then use a six-channel T -matrix and it is clearly not possible to determine all the T -matrix elements in a model independent way. An example of a plausible modelling has been proposed e.g. in ref. [70].
On the experimental side, a clear observation of a two-particle threshold effect requires a very good energy resolution (for an example, see ref. [71]). The data points tabulated in ref. [6] are obtained after implementing a procedure for unfolding the energy resolution effects and subtracting the background. It is possible that this procedure becomes somewhat inaccurate near the tip of the Dalitz plot region.
Details of the numerical values of the T -matrix parameters corresponding to different fits are given in appendix C. Fig. 12 shows the phase shifts and inelasticity parameter of the S-matrix corresponding to two different fits. We observe that in these determinations the phase shift δ πη increases above the KK threshold while δ KK decreases. This pattern differs from the one found in the lattice QCD simulation [19].

a 0 (980) complex pole and couplings
In the determinations of the πη − KK two-channel T -matrix performed in ref. [24] based on analysing γγ scattering data, the a 0 (980) resonance was always found to correspond to a pole of the T -matrix located on the second Riemann sheet. No pole was found on the third or fourth Riemann sheet with a real part close to 1 GeV 2 . We recall the formulae which define the second sheet extension of the T -matrix elements (see e.g. [40]) with S 11 (z) = 1 − 2σ πη (z)T 11 (z),σ πη (z) = (m 2  It is indeed easy to verify that the matrix elements T II ij (s) satisfy the continuity equation across the cut Simple definitions of the a 0 coupling constants can be associated with residues of the T -matrix poles, A formal definition of the width Γ a 0 →πη can be associated with the coupling g a 0 πη , The a 0 coupling to K + K − , which is often quoted, is related to g a 0 KK by  Table 3: Mass and width of the a 0 (980) resonance associated with the complex pole position ( √ s a0 = m a0 − iΓ a0 /2) and the related coupling constants in several fits.
The magnitude of the couplings g a 0 πη , g a 0 K + K − as defined above should be approximately equal to those defined in refs. [11,30].
The analytic extensions of the φ radiative decay S-wave amplitudes to the 2nd Riemann sheet are defined as which shows that they should display an a 0 pole. A coupling constant g φa 0 γ can be defined from the residue This definition (which takes into account the presence of the nearby soft photon zero) matches to the one introduced in ref. [30].
The central values of the a 0 (980) mass, width and couplings, corresponding to several fits are collected in table 3. The last column corresponds to a fit which combines the γγ data and the φ data. In this fit, the a 0 mass is somewhat smaller as well as the coupling g φa 0 γ . As discussed in [24] an important source of uncertainty is generated by the "fixed" parameters in the T -matrix, mainly the couplings L i and the ratio c m /c d (see eqs. (171) (172)), as well as the value of the Adler zero in the γγ → πη amplitude. We have performed fits in which these parameters are varied around their central values in the same way as in [24] except for s γγ A which we vary here in the range [m 2 η , m 2 η + 6m 2 π ]. As a result of this, we would quote the following determination of the mass/width from these data and for the a 0 couplings g a 0 ηπ = 2.2 ± 0.2 (GeV) g a 0 K + K − = 2.8 ± 0.2 (GeV)

Scalar I = 1 form factors
As already mentioned, the I = 1 πη and KK scalar form factors are linearly related to the Omnès matrix elements which have been probed by the γγ → πη and the φ → γπη processes. In view of their applications to τ decays, it is convenient to define the scalar form factors from the matrix elements of the charged scalar currentūd(x) as where s = (p π + p η ) 2 = (p K + + p K 0 ) 2 and B 0 is the O(p 2 ) chiral coupling proportional to the quark condensate [72]. In a minimal construction, the form factors are related to the Omnès matrix simply by (the minus sign being caused by |π + η = −|I = 1, I z = 1 ). At leading chiral order, the values at s = 0 are given by Expressions for the NLO corrections can be found e.g. in [44,73]. Using the determination BE14 [74] for the chiral coupling constants L i (as used also here for the T -matrix) one finds for the central values Other chiral constraints that one might consider concern the scalar radii which are not difficult to compute at order p 4 (see [44]). Using the BE14 value for the coupling L 5 , L r 5 = 1.01 · 10 −3 , one gets The minimal dispersive model gives a result which is too small by 50% for the πη radius and too large by 40% for the KK one as compared to the O(p 4 ) values. While these differences can be ascribed, in part, to higher order corrections, it seems interesting to consider a nonminimal dispersive model in which one enforces exactly the O(p 4 ) values for both F P Q S (0) anḋ F P Q S (0). This is easily achieved by replacing F P Q S (0) in eq. (141) by linear functions of s, F P Q S (0)(1 + λ P Q s), and adjusting the slope parameters λ P Q . Fig. 13 shows results for the real and imaginary parts of the πη form factor in the minimal dispersive model using either LO or NLO values for F P Q S (0) and the result from the non-minimal model. Differences between the models remain rather small in the region of the a 0 (980) resonance which is important for the τ decay. Our results are rather similar to those presented recently in ref. [73] but differ somewhat from ref. [75].

Summary and conclusions
We have re-examined the φ radiative decay amplitudes φ → γπη, γKK based on coupledchannel Omnès-type dispersive representations for the S-waves, which seems not to have been done previously. By construction, these amplitudes have the correct analytic structure, with a left-hand cut having several components in addition to the unitarity cut, and the representation involves explicit integrations along these cuts. The discontinuities across the left-hand cut components were assumed to be dominated by the contributions from the light vector-meson resonance exchanges. The coupling constants which are needed can be determined approximately (including the relative signs) by combining experimental inputs with flavour and chiral symmetry arguments. This approach, which uses only the discontinuities of the vector-exchange amplitudes has the advantage of being insensitive to the off-shell behaviour of the resonance propagators (polynomial ambiguities, form factors) and allows to take into account the effects of the resonance widths in a numerically fast way.
The kaon exchange Born contribution to the φ → γK + K − amplitude plays an important role in the dynamics, as is well known, because it carries a pole when s = m 2 φ . It proved convenient to split this amplitude into a singular part and a regular part (which vanishes in the soft photon limit). The integrals involving the first part can be computed analytically in terms of matrix elements of the D matrix. In practice, we describe the πη/KK scattering using a two-channel unitary T -matrix model proposed in [44], which has both left and right-hand cuts, and involves six phenomenological parameters. The corresponding Ω (and D) matrices are computed numerically by solving a set of Muskhelishvili integral equations. In some previous work an apparently different representation was used (e.g. [11,51,52,53,54,55,31] ) which involves a one-loop triangle function and the T -matrix itself. As shown in appendix A this representation holds only for restricted classes of T -matrices which, in particular, have no lefthand cuts.
We have shown that a rather good description of the experimental data on φ → γπη can be obtained fitting only two subtraction parameters, while using in the T -matrix a set of parameters determined previously from photon-photon scattering data. The fitted values of the two subtraction parameters were shown to be physically reasonable, such that the φ → γπη amplitude displays an Adler zero and the φ → γKK amplitude is dominated by the sum of the K * -exchange and the regular part of the K-exchange amplitudes at low energy.
These results seem to confirm the validity of the two-channel T -matrix model and the corresponding Ω matrix. As an application, we give results for the πη and the (KK) I=1 scalar form factors (which are simply related to the Ω matrix), they were shown to be well determined in the region of the a 0 (980) resonance using O(p 4 ) chiral constraints at s = 0. The πη form factor is measurable, in principle, from the τ → πην decay mode. Finally, we give some results on the a 0 (980) complex pole, performing combined fits of γγ data and φ radiative decay data.
Let us illustrate this in the case of the T -matrices based on unitarised χPT used in refs. [55,31]. In these models the T -matrices have no left-hand cut and the D and Ω matrices can be expressed explicitly as follows where T (2) (s) is the chiral O(p 2 ) T -matrix and J P 1 P 2 (s) are one-loop functions which satisfy In the Omnès representations we had separated the Born amplitude into a singular and a regular part we can first express the right-cut integrals as using the unitarity relation for the matrix elements D i2 (s) where the first equality is general and the second one follows from (147). The fact that T i2 (s ′ ) are linear functions of s ′ is the key to the simplification. Indeed, we can express the Born leftcut integrals using a dispersive representation for the productk 1,B 0++ (s)D i2 (s)/(s − q 2 ) (which remains finite when s = q 2 ) obtaining Then, replacing Im [D i2 (s ′ )] with (151) and adding the two integrals one gets We can then choose the subtraction point s 0 such that T i2 (s ′ )/(s ′ − s 0 ) is a constant, which effectively pulls out T (2) i2 (s ′ ) from the integral. The remaining integral is equal to the kaon loop function I loop K (s, q 2 ). Finally, one arrives at the following form for the two γγ * amplitudes, l 0++ (s) using that Ω(s)T (2) (s) = T (s) in this model. The integrals associated with meson exchanges e.g. I LC,V (s) in the Omnès representation can be similarly re-expressed in terms of simple triangle loop functions using eq. (151) which allows one to recover the formulae of ref. [32].

B ρ-meson spectral function
Integrals involving the spectral function ρ V which describe the width of the ρ meson, in particular ρ V 1 , ρ V 2 (eqs. (115)), are easily computed in the model given by eq. (110) by writing the denominator of ρ V , which is a cubic polynomial, as a product of three factors, where z R is a real root and z ± are complex conjugated. These roots must be computed numerically. The integrals of interest involve the following functions For instance, one has Taking the limit t → ∞, one can determine N V from the normalisation condition (109) which yields We can express in a similar way the following two integrals, with with The integrals ρ V 1 , ρ V 2 introduced in eq. (115) are simply related toρ V 1 ,ρ V 2 , Finally, using this model, the propagator of the ρ meson (eq. (108)) is expressed as a sum of four terms C K-matrix parametrisation We briefly recall below the parametrisation of the two-channel (πη, KK) S-wave scattering amplitudes introduced in ref. [44]. The T matrix is expressed in terms of a 2 × 2 symmetric K-matrix, T (s) = (1 − K(s)Φ(s)) −1 K(s) Two-channel unitarity is satisfied, as is well known, provided that the K-matrix elements are real in the physical region s ≥ (m η + m π ) 2 and the Φ-matrix satisfies We use a Φ-matrix which involves four phenomenological parameters Φ(s) = α 1 + β 1 s + 16πJ πη (s) 0 0 α 2 + β 2 s + 16πJ KK (s) whereJ P 1 P 2 (s) are one-loop functions as defined in ref. [72], which satisfyJ P 1 P 2 (0) = 0 and Im [16πJ P 1 P 2 (s)] = σ P 1 P 2 (s). The parameters α i , β i are real and assumed to be of chiral order p 0 . The K-matrix is written as a sum of three terms of increasing chiral order K(s) = K (2) (s) + K (4) (s) + K (6) (s) in which the terms of order p 2 and p 4 are determined by matching to the chiral expansion of the T -matrix i.e.
where Φ (0) is the O(p 0 ) part of Φ. The chiral expansions of T ij up to NLO were first derived in [76] and in [77,78]. The matching can be performed exactly for the matrix elements T 11 , T 12 but only approximately for T 22 in order to have K 22 real. The matrix elements K (4) ij depend on the Gasser-Leutwyler [72] coupling constants L r i . In this work, their values were taken from table 3 of ref. [74] (BE14 set). The third term, K (6) , implements a K-matrix pole  where λ and m 8 are phenomenological parameters while the form of g i is derived from a resonance chiral Lagrangian [79] The central values of c m , c d were taken as The form described above is used in a limited energy region s ≤ s cut with √ s cut = 1.5 GeV. The T -matrix must also be defined for s > s cut in order to solve for the Ω matrix. In this region, a simple interpolation of the two phase shifts (such that the sum goes to 2π at infinity) and of the inelasticity (going to 1 at infinity) is implemented. We collect the numerical values of the sets of parameters corresponding to the fits described in sec. 5.4 in table 4 below. D Estimate of the a 2 (1320) contribution The contributions to the φ → γπη helicity amplitudes induced by the a 2 (1320) tensor resonance, i.e. φ → γa 2 (1320) → γπη have the following form, where the constant C T is proportional to the product of the a 2 πη and a 2 φγ couplings. Unfortunately, there is no experimental information on a 2 → φγ decay. It is possible to get a qualitative estimate of the a 2 φγ coupling, relating it to the known a 2 → 2γ coupling, using vector meson dominance ideas, following ref. [80]. This gives where the coupling constants C a 2 γγ , C a 2 πη , defined in [24], have the values C a 2 γγ = (0.115 ± 0.005) GeV −1 C a 2 πη = (10.8 ± 0.5)GeV −1 .