Analyticity of ηπ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta \pi $$\end{document} isospin-violating form factors and the τ→ηπν\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tau \rightarrow \eta \pi \nu $$\end{document} second-class decay

We consider the evaluation of the ηπ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta \pi $$\end{document} isospin-violating vector and scalar form factors relying on a systematic application of analyticity and unitarity, combined with chiral expansion results. It is argued that the usual analyticity properties do hold (i.e. no anomalous thresholds are present) in spite of the instability of the η\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta $$\end{document} meson in QCD. Unitarity relates the vector form factor to the ηπ→ππ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta \pi \rightarrow \pi \pi $$\end{document} amplitude: we exploit progress in formulating and solving the Khuri–Treiman equations for η→3π\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta \rightarrow 3\pi $$\end{document} and in experimental measurements of the Dalitz plot parameters to evaluate the shape of the ρ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho $$\end{document}-meson peak. Observing this peak in the energy distribution of the τ→ηπν\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tau \rightarrow \eta \pi \nu $$\end{document} decay would be a background-free signature of a second-class amplitude. The scalar form factor is also estimated from a phase dispersive representation using a plausible model for the ηπ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta \pi $$\end{document} elastic scattering S\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S$$\end{document}-wave phase shift and a sum rule constraint in the inelastic region. We indicate how a possibly exotic nature of the a0(980)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a_0(980)$$\end{document} scalar meson manifests itself in a dispersive approach. A remark is finally made on a second-class amplitude in the τ→ππν\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tau \rightarrow \pi \pi \nu $$\end{document} decay.


Introduction
Isospin-breaking phenomena involving light pseudoscalar mesons are particularly interesting probes of the three flavour chiral expansion as they are driven by the parameter which involves the three light-quark masses. One of the major goals of non-perturbative approaches to QCD is to arrive at an accurate determination of the light-quark masses. One issue in the direct determination of from the very precisely a e-mail: moussall@ipno.in2p3.fr known difference between the masses of the charged and the neutral kaon is to properly evaluate the electromagnetic contribution to this difference. At leading chiral order, it is given by Dashen's low-energy theorem [1]. There are, however, possible substantial corrections from next-to-leading O(e 2 m s ) effects which suggests to explore, in parallel, other isospin-violating processes.
In this respect, the η → 3π decay amplitude is of particular interest since electromagnetic contributions are absent at leading order [2] and found to be rather small at nextto-leading order (NLO) [3][4][5]. Furthermore, there has been considerable progress, on the experimental side, in the precision of the measurements of the Dalitz plot parameters for both the η → 3π 0 amplitude [6][7][8][9] and the η → π + π − π 0 amplitude [10]. There is a price to pay, however, in that the chiral expansion has an inherently slow convergence in the treatment of final-state interactions, e.g. in the NLO expression of the η → 3π amplitude [11], ππ rescattering is treated only at leading chiral order. The amplitude has now been computed to NNLO in chiral perturbation theory (ChPT) [12]. A partly analytic representation of the ππ rescattering part at NNLO, accounting for some effects of higher order, was obtained in Ref. [13]. A treatment of rescattering in the framework of non-relativistic effective field theory has been discussed in Ref. [14]. An alternative approach is to combine the chiral expansion with a more general representation which encodes exact unitarity, analyticity and crossing symmetry [15]. A rigorous framework was proposed by Khuri and Treiman (KT) [16] who derived a set of integral equations for the analogous K → 3π problem. Application to the η → 3π amplitude was first discussed in Ref. [17]. The KT equations were more recently generalised to account for both S and P-wave elastic rescattering [18,19] and numerical solutions were constructed. Updates of these analyses, which take into account the recent experimental data, have been presented [20][21][22]. This progress have motivated us to reconsider the problem of evaluating the ηπ isospin-violating vector and scalar form factors exploiting, as systematically as possible, their analyticity properties and matching with chiral NLO calculations [23,24]. The vector form factor, in particular, probes the ηπ → ππ amplitude, via unitarity, in a kinematical region different from that of the decay, but where the KT equations should still be applicable. These two form factors are measurable, in principle, from the τ → ηπν decay mode. This mode, being forbidden in the isospin-symmetric limit, is a clean example of the "second-class currents" as introduced by Weinberg [25], which are yet to be discovered experimentally. An upper bound on the branching fraction, B ηπ < 1.4 × 10 −4 , was obtained by the CLEO collaboration [26] which was slightly improved to B ηπ < 9.9 × 10 −5 by Babar [27]. The Belle collaboration has quoted B ηπ < 7.3 × 10 −5 as a preliminary result [28] which, however, was not subsequently confirmed. Theoretical estimates for this branching fraction [23, [29][30][31][32][33][34] yield values in the range 10 −6 to 10 −5 , which do not seem so small as compared to the number of τ pairs available at Babar: N τ 4.3 × 10 8 or Belle N τ 9.0 × 10 8 . One difficulty faced by the B factories was that a substantial number of ηπ pairs were produced from background modes, like τ − → ηπ − π 0 ν, which have to be subtracted. A drastic reduction of this background should be possible at τ -charm factories [35,36], which could make detailed measurements of the ηπ mode possible. An increase in the luminosity by a factor of 50 is expected at Belle II [37].
Here, we will consider not only the integrated τ → ηπν τ branching fraction but also the detailed dependence as a function of the ηπ invariant mass. This dependence carries nontrivial dynamical information. Close to s = 0 it can be related, via ChPT, to the isospin-breaking quark mass ratio while, at higher energy, the shape of the ρ resonance peak in the vector form factor can be related to η → 3π decay properties. From an experimental point of view, the observation of the ρ peak would be a background-free signal of a second-class amplitude. A peak at the a 0 (980) mass is also expected from the scalar form factor. As was noted long ago [38], the mode τ → ηπν probes the "nature" of the scalar meson a 0 (980) in a clean way via its coupling to theūd operator.
The ηπ scalar form factor has a potential for constraining extensions of the Standard Model which contain charged Higgs bosons. For illustration, in the two-Higgs model proposed in Ref. [39], 1 the energy dependence of the form factor is modified as follows: where the ζ are coupling constants. The influence of the charged Higgs in this form factor could be enhanced because of the m u − m d denominator, depending on the relative sign and size of ζ u , ζ d . The constraints which are already available (specifically from B → τ ν) are not so stringent: |ζ l ζ d /m 2 H + | < 0.1 GeV −2 [40]. In order to derive a similar level of constraint, one should be able to evaluate f ηπ 0 in the Standard Model (and also be able to measure it, of course) with a precision of 20 % at s = 1 GeV 2 .
The plan of the paper is as follows. After introducing some basic formulae and notation, we list the contributions from the light two-meson states to the unitarity relations of the two form factors. We also discuss contributions with one photon. Then we recall the main results from the NLO ChPT calculations [23,24]: the values of the form factors at s = 0 and their first derivatives will be used as input in the dispersive representations. In order to derive these, it is important to check the possible presence of anomalous thresholds, since the η meson is unstable: we present arguments that they are actually absent. We then discuss the dispersive evaluation of the vector form factor using as input πη → ππ amplitudes satisfying the KT equations and constrained by experimental data in the physical decay region. Finally, we estimate the scalar form factor from a phase dispersive representation, using a modelling of ηπ → ηπ elastic scattering borrowed from Ref. [41].

Definitions and basic unitarity relations
The semi-leptonic weak decay amplitudes τ → ηπν and η → lπν (with l = e, μ) are induced by the usual Fermi Lagrangian The ηπ matrix element of the charged vector current is expressed in terms of two form factors (we follow the same notation as Ref. [23] except that we call the ηπ invariant mass squared s instead of t), with When writing unitarity relations it is convenient to introduce the scalar form factor f ηπ The expression for the differential decay width of the τ lepton which derives from the Fermi Lagrangian (3) and the definition of the form factors (4) then reads where S EW is the logarithmically enhanced universal radiative correction factor [42] (S EW = 1.0201 [43]) and λ being the Källén function λ(x, y, z) = x 2 + y 2 + z 2 − 2(x y + yz + xz). For the physical τ decay, the variable s lies in the range (m η + m π ) 2 ≤ s ≤ m 2 τ . The analogous formula for the differential decay of the η, η → l ± π ∓ ν l with l = e or l = μ reads In this case, the variable s is restricted to the range m 2 l ≤ s ≤ (m η − m π ) 2 .

Unitarity relations for
We consider the ηπ + centre-of-mass system and choose the z-axis along the three-momentum of the η meson. The vector form factor is easily seen to be proportional to the matrix element of the third component of the vector current in this frame out η( p η )π + ( p π )| j ud where q ηπ (s) is the centre-of-mass momentum. The form factor f ηπ + (s) can be defined as an analytic function of s with a cut along the positive real axis starting at s th = 4m 2 π (see the discussion of the absence of anomalous thresholds in Sect. 4 below). The discontinuity across the cut has the form of a generalised unitarity relation and is given as a sum over a complete set of states, with The lightest state contributing to the unitarity relation is n = π 0 π + with angular momentum l = 1. The next-tolightest contribution is from four pion states n = π 0 π + ππ. However, we expect such contributions not to be effectively relevant below 1 GeV because of phase-space suppression and we will ignore them here. Let us consider successively the contributions from the lightest two-body states n = π 0 π + , n = ηπ + and n =K 0 K + (a) n = π 0 π + : The π 0 π + matrix element of the vector current, involves two form factors since m π + = m π 0 . The unitarity relation for f ηπ + involves only the vector form factor, F π V (s). Using Eq. (11), we can write the ππ contribution in the unitarity relation as follows: where z = cos θ , θ being the scattering angle in the centre-of-mass system. We expect this contribution to be important below 1 GeV, because of the presence of the ρ(770) resonance. (b) n = ηπ + Next, the contribution from the ηπ state to the unitarity relation reads This contribution involves the ηπ → ηπ amplitude projected on the P-wave. The quantum numbers of the state (ηπ ) l=1 are exotic: J PC = 1 −+ . We expect the (ηπ → ηπ) l=1 amplitude to be very small below 1 GeV. This is borne out by the ηπ scattering model proposed in Ref. [41], which predicts that the P-wave phase shift is of the order of −1 • at 1 GeV. (c) n =K 0 K + Finally, let us write the contribution of theK 0 K + state in the unitarity relation, which is useful for comparing with the chiral calculation. In this case, the kaon vector form factor appears, defined from The corresponding contribution in the unitarity relation reads In the above expression, isospin breaking is contained in the amplitudeK 0 K + → ηπ + projected on the P-wave. Let us recall the reason: since G|K + = |K 0 , G|K 0 = −|K + one has G|K 0 K + = (−1) l+1 |K 0 K + . Since ηπ has G-parity −1, this implies that the partial-wave amplitudes (K 0 K + → ηπ + ) l with odd angular momentum l vanish in the isospin limit.

Unitarity relations for
Unitarity relations for the scalar form factor can be derived in exactly the same way as above noticing that, in the centreof-mass frame, the matrix element of the zeroth component of the vector current is proportional to One can then derive a relation for the discontinuity along the cut, analogous to Eq. (11), As before, let us consider the contributions from the lightest two-particle states ππ, πη and KK .
The contribution from the ηπ + states to the unitarity relation reads It has a form similar to Eq. (15) for the vector form factor except that it involves the ηπ → ηπ amplitude projected on the S-wave instead of the P-wave. This contribution is enhanced by the presence of the a 0 (980) resonance and thus must be the dominating one below 1 GeV. (c) n =K 0 K + : Finally, the contribution fromK 0 K + involves the corresponding scalar form factor and one has the following expression: disc f ηπ 0 As compared to the analogous contribution for f ηπ + , the relation (24) involves the ηπ + →K 0 K + amplitude projected on the S-wave, which is isospin conserving. Isospin breaking is contained in the mass difference factor K 0 K + .

Some electromagnetic contributions to the unitarity relations
In the unitarity equations discussed above, we have considered only hadronic states in the sums over n. Since we are studying isospin-breaking form factors, electromagnetic contributions are present and at order e 2 one should also consider states involving one photon. Note that EM contributions have already appeared, e.g. in Eqs. (21), (24), which are proportional to the mass differences m 2 π 0 − m 2 π + (which is mainly electromagnetic) and m 2K 0 − m 2 K + (which is partly of electromagnetic origin). These contributions are dominant in the chiral counting; they are included in the NLO chiral expressions. We will not discuss EM contributions in their full generality here and simply mention the contributions of the two lightest states n = γ π and n = γ ππ in the unitarity relations: The γ π matrix element of the vector current can be expressed in terms of one form factor where e γ is the polarisation vector of the photon. 3 At leading order in the chiral expansion, the value of this form factor at s = 0 is given by the anomaly Going to the centre-of-mass frame, one sees that the matrix element (25) vanishes for j ηπ 0 . The unitarity contribution from γ π thus concerns only the vector form factor. One finds the following expression for the discontinuity: Evaluating this contribution precisely would require some modelling of the amplitude ηπ → γ π. It is likely that this amplitude should be small below 1 GeV since no resonant contribution from the ρ-meson is allowed in the isospin limit.
In principle, the states n = γ ππ can contribute to the unitarity relations for both f ηπ + and f ηπ 0 . We will consider here only the latter one, the evaluation of which is simplified by using the relation between f ηπ 0 and the matrix element of the divergence, out η( p η )π( p π )|i∂ μ j ud together with the Ward identity for the divergence, Equation (29) makes it easy to evaluate the matrix element involving γ ππ in terms of the pion vector form factor out γ (λ)π 0 π + |∂ μ j ud with s ππ = ( p π 0 + p π + ) 2 . One can then write the unitarity relation in the form (where dLips 3 is the three-body Lorentz invariant phasespace measure). A resonant contribution from the a 0 (980) to the amplitude ηπ → γ ππ which appears in Eq. (31) is possible. However, a suppression of this contribution in the region below 1 GeV is expected because of the three-body phase space.
In summary, below the KK threshold, the dominant contribution is from n = π 0 π + [enhanced by the ρ(770) resonance] for the vector form factor and from n = ηπ + (enhanced by the a 0 (980)) for the scalar form factor. We will use this result in the sequel in order to evaluate the two form factors with the help of dispersion relations. In order to suppress the sensitivity of the integrals to the region s ≥ 1 GeV 2 it is necessary to introduce subtractions. We now recall the results which have been obtained in ChPT at NLO from which we will be able to estimate the subtraction constants.

Results from ChPT at order p 4
The ηπ form factors have been computed in ChPT at nextto-leading order in Refs. [23,24]. We collect below some of the results which are relevant to our study.

Form factors at s = 0
Consider first the form factors at s = 0. At leading order in ChPT they are simply equal to the π 0 -η mixing angle, where is given in (1). At the same order, can be determined using the experimental values of the pseudoscalar meson masses m π + , m π 0 , m K + , m K 0 , together with Dashen's lowenergy theorem [1], which gives The corrections of order p 4 , including also the electromagnetic e 2 p 2 piece, were written in Ref.
[23] (see also [44]) in the following form: with (an additional small electromagnetic contribution, proportional to m 2 π / ηπ has been neglected). Here, L r 7 , L r 8 are the standard low-energy coupling constants of the strong chiral Lagrangian at O( p 4 ) [45], while Z , S r 2 , S r 3 are O(e 2 ) and O(e 2 p 2 ) electromagnetic couplings [46,47]. One can also express f ηπ + (0) in terms of the two π 0 -η mixing angles 1 , 2 introduced in Ref. [45], (the electromagnetic contributions to 1 , 2 can be found in [48]). In using Eq. (35), one must also update the determination of the mixing angle , e.g. from the pseudoscalar meson masses, including O( p 4 ) and O(e 2 p 2 ) contributions. At present, however, the O(e 2 p 2 ) LECs, which play an important numerical role for the meson masses, are not known in a model independent way. Fortunately, it was observed in Ref.
[23] that a very simple relation holds between f ηπ + (0) and an isospin-breaking difference involving the K 0 π ± and K ± π 0 vector form factors. In its updated form, it can be written in terms of the form factorsf K π + defined in Ref. [44] to include the isospin-breaking effects from QCD and part of the radiative corrections as This relation was not used in earlier work on the secondclass ηπ amplitudes because the precision of the experimental results on the K l3 form factors was insufficient. The situation has considerably improved in recent years and, using the averaged experimental result quoted in the review [49], one obtains Clearly, this is a very significant enhancement as compared to the leading order result. For comparison, the following range of values was quoted in Ref.
[23]: f ηπ + (0) = [1.22-1.37] × 10 −2 based on the chiral expression (35), using the LO result for and order of magnitude estimates for the electromagnetic coupling constants. A result completely compatible with (40) can be obtained if the value of the quark mass ratio is enhanced from its LO result by 20-30 % due to NLO effects. There are indications that this could indeed be the case from model estimates as well as lattice QCD calculations of the chiral corrections to Dashen's low-energy theorem [see e.g. the FLAG review [50]]. We will return to the question of the quark mass ratio in Sect. 5.4 in connection with η → 3π decay. The fact that compatible results are obtained using either the low-energy theorem relation (38) or the chiral expression (35), which have different O( p 6 ) corrections, is an indication that O( p 6 ) corrections should be of a natural size (5-10 %, say) in spite of the large size of the NLO correction.

Vector form factor
We reproduce below the expression of the vector form factor f ηπ + (s) from Ref.
[23], in a slightly re-expressed form, which involves the scalar loop functionsJ P Q (s) [45], The result for f Equation (42) allows one to deduce the value of the derivative of the vector form factor at s = 0, which will serve us, together with f ηπ + (0), to normalise the dispersive construction of the form factor. One findṡ In the last equality, we have used the value of the coupling L r 9 given by [51], L r 9 (m ρ ) = (5.93 ± 0.43) 10 −3 .

Scalar form factor
Next, the scalar form factor f ηπ 0 from Ref.
[23] can be expressed as follows: with One expects that the scalar form factor, evaluated at the point s = ηπ = m 2 η − m 2 π should satisfy a Callan-Treiman relation [52,53] where F (3) P is defined from the matrix element of the axial current j 3 Indeed, using the expressions for F π 0 computed in Ref. [47] at chiral order p 4 and e 2 p 2 one obtains that CT is proportional to m 2 π and reads The Callan-Treiman relation (46) could be used to evaluate f ηπ 0 ( ηπ ) rather precisely if the ratio F η /F π were known accurately from lattice QCD. For now, we must rely only on the chiral expansion and, as in the case of the vector form factor, we will use the value of the derivative of the scalar form factor at s = 0 as an input to the dispersive calculation. Using the chiral expression (45) with the value of the LEC L 5 : L r 5 (m ρ ) = (1.20 ± 0.05) × 10 −3 deduced from the ratio F K /F π = 1.192 ± 0.005 [50] one obtainṡ f ηπ 4 Absence of anomalous thresholds in ηπ form factors in a toy model A usually accepted property of form factors involving two stable particles (like the pion or the nucleon electromagnetic form factors) is that they can be defined as analytic functions of the energy variable s = ( p 1 − p 2 ) 2 , with a cut along the positive real axis, the discontinuity along this cut being given by unitarity relations (e.g. [55]).
Here, we wish to consider the ηπ form factors f ηπ + , f ηπ 0 and, since the η meson is unstable, one should be concerned about the presence of anomalous thresholds. We illustrate in Fig. 1 the two types of diagrams which involve the η → 3π decay amplitude at one vertex. We will consider here contributions of the form of Fig. 1a in which the normal threshold is s th = 4m 2 π . The contributions of the second type, as shown in Fig. 1b have a much higher normal threshold s th = 16m 2 π and the discontinuity function is expected to be very much suppressed in the region s 1 GeV 2 because of the four-body phase space. Let us discuss here the question of the anomalous threshold in the toy model case where Fig. 1a represents a Feynman diagram with local vertices. The form factor can then be represented as an integral [e.g. [56]] where K ηπ (t , s) corresponds to the simple triangle Feynman diagram with external momenta p 2 1 =m 2 η , p 2 2 = m 2 π , p 2 3 = s and internal masses m 2 1 = m 2 2 = m 2 π , m 2 3 ≡ t (see Fig. 2), and the weight function is In Eq. (50), additional terms (including the UV divergent ones) have been omitted since they are not concerned with the possibility of an anomalous threshold. We will vary the value ofm 2 η and, since it can be considered as an energy variable, the amplitude may be defined by appending an infinitesimal positive imaginary part to it, i.e.
We first takem 2 η to be sufficiently small (e.g.m 2 η = m 2 π ) such that an ordinary dispersion relation (DR) holds and then increasem 2 η until it reaches the physical value m 2 η . The ordinary DR for the triangle graph reads, where L ηπ (t , s ) is the discontinuity function of the triangle graph, We note that the discontinuity of the form factor is then given as an integral over L ηπ (t , s) As discussed in Refs. [57,58] the presence of anomalous thresholds can be inferred from studying the motion of the singularities of the function L ηπ upon varyingm 2 η : if one of the singularities crosses the unitarity cut, it is then necessary to deform the contour in the dispersion representation (53), in order to properly define its analytical continuation as a function ofm 2 η , thereby introducing an anomalous threshold. In the present case (54), the singularities of the function L ηπ (t , s) are given by the solutions of the equation Let us consider three cases, depending on the value of the mass variable t  . 3 Illustration of the motion of the two singularities of the discontinuity function L ηπ (t , s) in the s plane for a given value of t when varying the external mass squaredm 2 η ≡ xm 2 η . The plot corresponds to t = 5m 2 π , = 10 −3 m 2 π 1. t = 4m 2 π : In this case, the two singularities coincide and are given by which is above the unitarity cut.
: This is the most interesting situation, the motion of the two singularities s ± (m 2 η ) as a function of x =m 2 η /m 2 η is illustrated in Fig. 3. The figure shows that while s + (m 2 η ) remains above the real axis, the other singularity does cross the real axis very close to 4m 2 π . It is easy to see that the crossing occurs whenm 2 η = m 2 π + 2t and the value of the crossing point is which is located strictly below the normal threshold. 2 : In this case s + (m 2 η ) remains above the real axis and s − (m 2 η ) remains below it.
For larger values of t , it is easily verified that s ± (m 2 η ) do not come close to the unitarity cut. The conclusion of this discussion is that, for a Feynman diagram, the amplitude as given from Fig. 1a does not involve any anomalous threshold. We will argue in the next section that the same conclusion holds in a more realistic approach where the η → 3π amplitude is given from Khuri-Treiman equations solutions. The fact that the η meson is unstable, i.e. m η > 3m π , manifests itself in the violation of real analyticity, the discontinuity function is complex and does not coincide with the imaginary part of the diagram (the latter, indeed, is given by the sum of two Cutkosky contributions, corresponding to two physically allowed ways of cutting the diagram; see Fig. 1a).

Dispersive evaluation of f ηπ +
We have argued in Sect. 2 that the dominant contribution to the discontinuity of f ηπ + (s), with s 1 GeV 2 , was from the n = π 0 π + state. Then disc[ f ηπ + ] was found to be proportional to the pion vector form factor F π V , which is well known from experiment, and to the l = 1 projection of the ηπ + → π 0 π + amplitude. In order to evaluate this amplitude, we will make use of the work of Refs. [18,19], who developed and solved a set of Khuri-Treiman [16] equations and applied the results in the physical region of the η → 3π decay. We briefly recall the main features of this formalism below (a very detailed account can be found in the thesis [59]). As it encodes the analyticity properties of the ηπ + → π 0 π + amplitude, this formalism should be suitable for evaluating the amplitude in the partly unphysical region needed for computing disc[ f ηπ + ].

Brief review of the Khuri-Treiman formalism
The Khuri-Treiman (KT) equations implement dispersion relations, crossing symmetry and unitarity in the approximation where a single state, n = ππ, is retained in the unitarity relations. This approximation is acceptable when the Mandelstam variables s, t, u are smaller than 1 GeV 2 in magnitude. It was noted in Refs. [18,19] that, in the same region, the contributions from the discontinuities of the partial-waves J ≥ 2 can also be neglected in the dispersion relations such that the decomposition theorem [60] may be applied to the η → 3π amplitude. As a result, it can be expressed in terms of three functions of a single variable, M I (w), Based on the usual Regge phenomenology for estimating the asymptotic behaviour, it was concluded in Ref. [19] that M 0 , M 2 should obey converging DRs with two subtractions and M 1 a converging DR with a single subtraction, In writing these DRs one has further made use of freedom to redefine M I by linear functions without modifying the ηπ → ππ amplitude in Eq. (59) because of the constraint s + t + u = m 2 η + 3m 2 π ≡ 3s 0 . The functions M I (w) are analytic in the complex w plane except for a cut along the real axis along [4m 2 π , ∞]. The discontinuity along this cut is obtained from the unitarity relations of the J = 0, 1 partialwave projections of the ηπ → ππ amplitudes and they read, In Eq. (62) the functionsM I are linear combinations of the angular integrals, with t (s, z) = 1 2 (3s 0 −s+κ(s)z), κ(s) = 1 − 4m 2 π /s λ ηπ (s).
The explicit expressions ofM I in terms of the angular integrals read [19] Equations (61) For completeness, let us mention the analogous relation for the scalar form factor, Adapting the discussion of Sect. 4 about the presence of anomalous thresholds to the present, more realistic situation, requires one to investigate the singularities of the functionsM I (w) i.e. of the angular integrals given in (63). This can be done by inserting the dispersive representations of the functions M I [Eqs. (61)] into the angular integrals [Eqs. (63)] from which one obtains an expression of the functions z n M I as integrals over kernels K (n) (t , w), The kernels which are needed here are given explicitly in Appendix A. They involve the logarithmic function which controls their singularity structure. When performing the angular integration in Eqs. (63), (69) one must keep in mind that the path of integration from z = −1 to z = 1 must eventually be deformed in order not to intersect with the cut of the functions M I , i.e. 4m 2 π ≤ t (w, z) < ∞, as is explained in detail in Ref. [18]. The following expression of the logarithmic function exactly encodes these prescriptions on the integration path in a simple way: where displaying explicitly the i prescription. From the form of the logarithmic function (70), (71) we can infer the following consequences: 1. Absence of anomalous thresholds: Inserting the representation ofM 1 ,M 0 in terms of the kernels (68) one obtains an integral representation of the form factor discontinu- in terms of the logarithm function (70) which is analogous to Eq. (55). Furthermore, the singularities of the logarithms are exactly the same. Therefore, the conclusion about the absence of anomalous thresholds, as discussed in Sect. 4, applies also in the present realistic situation.
2. Cuts C of theM I functions: They are given, using the integral representations in terms of kernels, as the ensemble of the singularities of the logarithmic function (70) (see [61]) i.e. the points w which satisfy t ± (w) = t . This relation can be recast as a quadratic equation in w identical to (56). Therefore, The curve (72) includes the negative real axis; it extends into the complex plane and approaches infinitesimally close the unitarity cut on the positive real axis in the region 4m 2 π ≤ w ≤ (m η −m π ) 2 without, however, crossing it: this is shown on Fig. 4. As a consequence, the functionsM I (s ) are unambiguously defined on the real axis in the integration range 4m 2 π ≤ s < ∞.
The point w = m 2 − ≡ (m η − m π ) 2 requires special attention. At this point the difference between t + (w) and t − (w) is infinitesimal but their imaginary parts have different signs which implies that the logarithmic function L(t , w) diverges when goes to zero and t lies in the range 4m 2 π ≤ t < m π (m η + m π ),

Matching with ChPT
In the form given by Eqs. (61), (62), (65) the KT equations are linear integral equations with a singular Cauchy kernel. The most general solutions of such integral equations involve an arbitrary number of polynomial parameters [62,63]. Physically, the polynomial growth is limited by asymptotic conditions on the amplitudes. In practice, however, the system of equations is valid only in the elastic scattering region, while the integrals run up to infinity. The polynomial part must thus be considered as a parametrisation of the corrections to the effects of the integration regions s , t ≥ 1 GeV 2 over the solutions in the elastic scattering region. For our purposes, we will consider here a four-parameter family of solutions. The polynomial dependence is exhibited by introducing the Omnès functions, where δ I is equal to the S or P-wave ππ phase shift with isospin I in the elastic region. The functions M I are then expressed as and they ensure that the discontinuities of Eqs. (62) are obeyed. The two subtraction constants which appear in the dispersive representation (61) are simply related to the polynomial parameters:α 0 = α 0 ,β 0 = β 0 + α 0˙ (0). Plugging this representation into Eq. (65) one obtains a set of linear integral equations for the functionsM I . As these are not singular equations anymore, one may thus expect that, for given values of the parameters α 0 , β 0 , γ 0 , β 1 , if a solution exists forM I , it should be unique [19]. The most natural idea for determining the polynomial parameters is by matching the KT amplitude with the amplitude computed in the chiral expansion [17][18][19] in the region where the variables s, t, u are small. More precisely, if M(s, t, u) is the amplitude computed to chiral order p N , then the parameters α 0 , β 0 , β 1 , γ 0 should be such that Considering the N = 4 case, a first observation is that the differences of the discontinuities in each of the component functions M I are of chiral order p 6 , independently of the values of the polynomial parameters. This implies that the O( p 4 ) parts of the differences M I (w) − M I (w) must be polynomial. Imposing the requirement that the O( p 4 ) parts of polynomial expansions of the differences M I − M I vanish yields the following four matching equations [19]: where the functions M a , a ,Î a and their derivatives are all to be taken at w = 0. In order to solve the set of Eqs. (81), one must keep in mind that the integralsÎ a carry an implicit linear dependence on the four polynomial parameters. This dependence must be determined by using four independent KT solutions in which one of the polynomial parameters is set to 1 and the others to 0. The chiral NLO amplitude M(s, t, u) was first computed in Ref. [11]; it is given in Ref. [19] in a form which involves a single O( p 4 ) coupling constant, L 3 . We will use here the value L 3 = (−3.04 ± 0.43)10 −3 [64].

Numerical solutions and comparisons with the η → 3π data
We have constructed numerical solutions of the set of KT equations by iteration. The main differences with earlier work [18,19] is that: (a) we have used the kernel representations (68) for performing the angular integrations, which should be somewhat faster than the integration over a complex contour method used previously, (b) The matching with ChPT was done via the four matching Eqs. (81), which were solved with no approximations.
In Ref. [19] only the last two matching equations were implemented while the first two were replaced by imposing the requirement that the amplitude M(s, t, u) along the line t = s has an Adler zero at the same position and with the same slope  [59] and the dashed line is the phase used in Ref. [18]. Also shown are some experimental ππ data from Refs. [65,66] as the chiral NLO amplitude. We will see below (Fig. 6) that the first two matching relations ensure essentially equivalent constraints on the Adler zero.
Concerning the phases δ I (s) for I = 1 and I = 2, for which inelasticity sets in rather smoothly, we take the phases to be equal to the corresponding ππ scattering phase shifts up to √ s cut = 1.4 GeV and, for s > s cut interpolate to the asymptotic values δ 1 (∞) = π and δ 2 (∞) = 0.
In the case of the I = 0 S-wave, inelasticity sets in sharply around the KK threshold. We have employed two different phase choices in the inelastic region: a) the one used in Ref. [18] (which we call cutoff (1)): for s larger than s cut = (0.865) 2 GeV 2 , the phase is interpolated rapidly to δ 0 (∞) = 0 and b) a condition similar to that used in Ref. [59]) (which we call cutoff (2)): for s larger than 4m 2 K , the phase is interpolated slowly to δ 0 (∞) = π . These phases are illustrated in Fig. 5. When used in the matching relations (81) these different conditions lead to rather different values for some of the polynomial parameters, 4 reflecting differences in the values of the integralÎ 0 (0) as well as in the values of the Omnès function and its derivatives at s = 0 (see Table 1). However, the matching conditions ensure that the complete KT amplitude depends only moderately on the cutoff conditions. Figure 6 illustrates some results, showing amplitudes along the line t = s where an Adler zero is present in the chiral amplitude.
In order to assess the reliability of the ηπ → ππ amplitude resulting from KT solutions with ChPT matching, let us compare with experimental results. From the integrated decay rate of the charged mode: ex p η→π 0 π − π + = 300 ± 11 eV [67], one obtains for the central value of the double  Precise measurements of the differential decay distributions across the Dalitz plot have been performed recently for both the charged [10] and the neutral decay modes [6][7][8][9]. It is customary to represent these differential distributions in terms of a polynomial of two independent energy variables X , Y (defined such that X 2 + Y 2 ≤ 1 and X = Y = 0 corresponds to the centre of the Dalitz plot, where the three pions have equal kinetic energies; see Appendix C); we have Comparison between the experimental values of the Dalitz plot parameters for η → π 0 π + π − and η → 3π 0 and the predictions of the KT solution amplitudes with NLO matching and two different cutoff conditions (see Table 1). for the charged mode, and for the neutral mode. The experimental values of the Dalitz plot parameters a, b, d, f and α are shown in Table 2 together with results corresponding to KT solutions amplitudes. Implementation of rescattering effects via the KT equations leads to significantly improved results with respect to the simple use of the chiral NLO amplitude (see [11]) but the two Dalitz parameters b and d are still predicted to be too large. This discrepancy between the theoretical amplitude and experiment indicates that further effects need to be taken into account. These could be either chiral O( p 6 ) effects at the level of the matching equations or further rescattering effects. In this respect, one sees clearly from Table 2 that the ππ phase choice which does include the f 0 (980) resonance leads to better results for the Dalitz plot parameters than the choice which does not. It is then not unlikely that the a 0 (980) resonance should be taken into account as well. It is also worth noting that preliminary results from analysis of new data sets by KLOE and WASA have been presented (see [68], p.16) which go in the direction of improving the agreement with the theoretical predictions.
For our present purposes, in addition to the KT amplitudes which obey the ChPT matching relations, we construct an amplitude which reproduces more closely the experimental results on the Dalitz plot. In order to do so, we allow the two polynomial parameters β 1 , γ 0 to vary freely (still assuming their imaginary parts to be negligible) and use them as fit parameters. The last two polynomial parameters α 0 , β 0 are then fixed by the two conditions: 1) that the amplitude reproduces the position of the Adler zero s A of the NLO amplitude and 2) that the central value of the quark mass double ratio Q from lattice QCD (the recent FLAG review [50] gives Q latt = 22.6(7)(6) from N f = 2 + 1 simulations) is reproduced. The results from this amplitude for Dalitz plot parameters are displayed in the last column of Table 2 and the corresponding polynomial parameters are shown on the last line of  Table 1.

Results for f ηπ +
In the elastic scattering region, the discontinuity of f ηπ + (s) was given in Eq. (14) in terms of the pion form factor F π V and the l = 1 projection of the ηπ + → π 0 π + amplitude. From the decomposition theorem, this projected amplitude may be expressed in terms of M 1 +M 1 (Eq. (66)). We can now calculate this quantity from our solutions of the KT equations. Figure 7 shows the numerical results for the modulus of M 1 +M 1 . It illustrates the strong sensitivity to the choice of the polynomial parameters. The solution corresponding to fitted parameters (last line in Table 1) has a significantly smaller resonance peak, which is related to the smaller size of the parameter β 1 .
Concerning the pion form factor, we used the experimental measurements of |F π V | from τ decays (which provide exactly the same form factor as needed here) by the Belle collaboration [69]. The measurement covers the energy range 0.297 ≤ √ s ≤ 1.255 GeV which is essentially adequate for our purposes. We relied on the fit performed in Ref. [69] in terms of Gounaris-Sakurai (GS) functions [70] for performing an extrapolation of |F π V | in the small energy region down to the threshold and performing the numerical integrations. For the phase δ π V , we assumed elastic unitarity to hold below 1 GeV and thus took δ π V to be equal to the l = 1 ππ phase shift, in accordance with Watson's theorem. Above 1 GeV, we use the phase as predicted by the GS parametrisation. 5 Asymptotically, from the usual QCD-based arguments [71,72], one expects the form factor f ηπ + (s) to behave as 1/s log(s), it should thus obey a convergent dispersion relation. We will actually use DRs for f ηπ + (s)/s N in order to suppress the contribution from the integration region above 1 GeV. In practice, we will use N = 1 or N = 2 and check the stability of the result. For example, with N = 2, the DR reads 5 Below 1 GeV, the phase δ π V produced by the GS parametrisation and the ππ phase shift are quite close, differing by 5 − 10 %, except at low energy, √ s < 0.5 GeV, where the difference is more significant.
KT+match 1-subtr. KT+match 2-subtr. KT+fit 1-subtr. KT+fit 2-subtr. The results obtained using these dispersion relations together with ηπ → ππ amplitudes from KT equations and also using the parametrisation given by the Belle collaboration [69] for the pion form factor are shown in Fig. 8. The figure shows that the DRs with N = 1 or N = 2 yield very similar results in the region √ s ≤ 1 GeV. As one can expect from the behaviour of M 1 +M 1 (see Fig. 7), the KT amplitude with matched polynomial parameters gives rise to a larger ρ peak than the one with the fitted parameters and the peak is more similar in shape to that of the pion vector form factor.

Dispersive estimate of f ηπ 0
In this section we provide a qualitative estimate of the scalar form factor, assuming that there should be analogies between the πη scattering amplitude and the well-known ππ, π K scattering amplitudes and also between f ηπ 0 and the π K scalar form factor.

Phase dispersive representation
The unitarity relations obeyed by the scalar form factor f ηπ 0 , associated with the two-body channels ππ, πη and KK were written in Sect. 2.2. Below the KK threshold, the contribution from πη largely dominates since ππ is comparatively suppressed by isospin symmetry. It is then convenient to write the form factor as a phase dispersive representation, 6 where ζ =ḟ ηπ 0 (0)/ f ηπ 0 (0) and φ ηπ is the phase of the form factor. The representation (85) uses two subtractions in order to suppress, as much as possible, the contributions from the higher-energy regions. The values of the form factor and its derivative at s = 0 can be taken from ChPT at O( p 4 ); see Sect. 3. Below the KK threshold, ηπ scattering can be assumed to be essentially elastic such that, in this region, the form factor phase can be identified with the l = 0 elastic scattering phase shift δ ηπ 0 by Watson's theorem. The S-wave ηπ scattering phase shift has not yet been directly measured, but detailed experimental information exists on the properties of the scalar resonances a 0 (980) and a 0 (1450). Furthermore, chiral symmetry constrains the phase shift to be very small at low energy [73]. For definiteness, we will make use of the ηπ-scattering model proposed in Ref. [41], which encodes these various pieces of information in a simple way. It uses constraints on ηπ scattering derived from the η → ηππ decay amplitude [74]. The ηπ amplitude is written in the following form: where the first term is the constant current algebra contribution, (accounting for η-η mixing, δ being the corresponding octetsinglet angle 7 ). The function A 0 represents a sum over treelevel amplitudes associated with the a 0 (980), a 0 (1450) resonances and F 0 is a similar sum involving the two isoscalar resonances σ (600) and f 0 (980). These amplitudes are computed from a resonance chiral Lagrangian and thus behave as O( p 4 ) at small values of the Mandelstam variables. 8 We used the set of resonance parameters from Eqs. (4.7), (A17), (A1) of Ref. [41]. The partial-wave amplitudes being given by we define the phase shift from the ansatz which applies also in the inelastic scattering region. The l = 0 phase shift from this model is shown in Fig. 9. It is qualitatively in good agreement with other approaches which have been proposed, like the chiral unitary approach [77]. It is also in agreement with one of the models used in Ref. [78] and probed against the recent high-statistics measurements of γ γ → ηπ scattering by the Belle collaboration [79]. Above the KK threshold, ηπ scattering becomes inelastic under the effect, mainly, of the two-body channels KK and η π . The phase φ ηπ must then differ from δ ηπ 0 . A global constraint on φ ηπ arises from imposing the requirement that the form factor, as given from Eq. (85), exhibits no exponential divergence asymptotically. This gives rise to a sum rule, Using the chiral expansion results for ζ (see Sect. 3.3), the sum rule indicates that one should have φ ηπ << δ ηπ 0 in the inelastic region. In order to estimate more precisely how φ ηπ behaves, one may rely on an analogy with the phase of the 7 A different convention for F π and for the mixing angle was used in Ref. [41]. 8 This model predicts a 0 3.2 × 10 −2 for the l = 0 scattering length. This value is somewhat larger than the NLO ChPT low-energy theorem result [75] a 0 = (−0.02 ± 0.77) × 10 −2 . It can possibly be accommodated in schemes where higher order effects associated with OZI violations are included [76]. K π scalar form factor. In that case, there is enough experimental information on the elastic as well as the leading twobody inelastic T -matrix elements, such that the form factor can be deduced from solving a set of Muskhelishvili-Omnès equations. The analysis performed in Ref. [80] shows that the form factor phase displays a sharp drop shortly after the onset of the leading inelastic channel. 9 A similar behaviour has also been observed in the case of the ππ scalar form factor associated with theūu +dd operator [82]. The authors argue that the presence of this phase drop is necessary in order to reproduce the correct value of the pion scalar radius via a sum rule analogous to Eq. (90). We then propose the following simple model for the phase φ ηπ assuming a fast decrease at s = s 1 and a constant value for s > s 2 with 4m 2 K < s 1 < s 2 : (a slight smoothing of the θ function is implemented in practice). For each value of s 1 , the value of s 2 is determined such that the sum rule (90) is exactly satisfied. Figure 9 illustrates the behaviour of φ ηπ with √ s 1 = 1.05, 1.15 GeV. In this model, s 1 is bounded from above: √ s 1 < 1.2 GeV, for otherwise it is not possible to satisfy the sum rule (90).

Results for f
ηπ 0 Using the phase φ ηπ as described above, one can compute the scalar form factor from the dispersive representation (85). Results are shown in Fig. 10 for several values of s 1 . The point s 1 , where the phase drops, corresponds to a dip in the modulus Obviously, if the dip is located very close to the KK threshold, the peak of the a 0 (980) resonance is strongly reduced. This corresponds, in this approach, to a reduced coupling of the resonance to theūd scalar operator and thus to an exotic nature of the a 0 (980). One may formally define a coupling constant from the matrix element of the vector current involving the a 0 state [83], In a dispersive approach, it is in principle possible to identify such a coupling constant from the residue of the a 0 resonance pole of the scalar form factor on the second Riemann sheet. In a simpler way, one may obtain an estimate by matching the shape of the dispersive form factor with a Breit-Wigner-type shape. For this purpose, we have used the chiral resonance approach of Ref. [23] in which one can vary the value of F a 0 (via that of an O( p 6 ) chiral coupling constant, d r ) while the value of the form factor at the origin is kept fixed. In this manner, from the phase dispersive representation, with the largest allowed value of the dip parameter, we find  [83,84].
In the same framework, it was found in Ref. [85] that, on the contrary, the a 0 coupling should be strongly suppressed, while estimates based on the bag model give a range of [0.2-2.0] MeV [86]. A result from a lattice QCD simulation with N f = 2 dynamical quarks has been given [87], which corresponds to F a 0 = [0.8-0.9] MeV.

Application to the τ → ηπν and η l3 decays
We can now compute the contributions from the vector and scalar form factors to the differential decay width of the τ into ηπν [the relevant formula was given in Eq. (7)]. Figure 11 shows the contribution from the vector form factor calculated from a KT amplitude with fitted polynomial parameters (lower curve) and that associated with parameters determined from ChPT matching (upper curve). The vector contribution is somewhat suppressed here by the kinematics but our calculations suggest that it should lead to a clearly visible ρ-meson peak. The contribution to the differential decay width from the scalar form factor is shown for three different values of the dip s 1 .
The corresponding values for the integrated branching fraction of the τ → ηπν mode are given in Table 3    to fitted polynomial parameters for the vector form factor. It is difficult to precisely evaluate the error, in particular in the case of the scalar form factor, because of the various assumptions and model dependence involved. A plausible guess however, in our approach, is that the scalar contribution to the branching fraction should lie in the range 10 0.1 × 10 −5 ≤ B S ≤ 0.6 × 10 −5 . This tends to be smaller than most previous estimates which are often based on simple scalar-dominance models. Conversely, in the dispersive approach, it seems difficult to accommodate a value for B S as small as quoted in Ref. [34] even if the position of the dip s 1 is very close to the KK threshold, because of the contribution from below the resonance region.
In the case of the vector contribution, in view of some possible shifts in the experimental value of the η → 3π Dalitz plot parameters (see [68], p.16), a plausible range for the branching fraction should be 0.
Finally, for the η l3 decays, we find the following central values for the branching fractions (adding the two charge modes): The results, in this case, are practically identical to those computed with the chiral NLO form factors [23,24]. An experimental upper bound on the η e3 mode branching fraction has been obtained recently by the BESIII [88] collaboration

Conclusions
In this work, we have reconsidered the ηπ isospin-violating vector and scalar form factors and the related energy distribution in the second-class τ → ηπν decay which should be measurable at future B or τ -charm factories. We have started from the NLO ChPT results for f ηπ + (0) and for the derivativesḟ ηπ + (0),ḟ ηπ 0 (0). In particular, for f ηπ + (0) a relation was established [23] with the K + l3 and K 0 l3 decays which can now be used thanks to recent experimental progress [49]. These results at s = 0 could be checked, in principle, in lattice QCD simulations including isospin violation.
In order to evaluate the form factors in the resonance regions we further relied on their analyticity properties. We argued that these should be the same as in the more familiar cases of the ππ or the π K form factors, i.e., no anomalous threshold should be present despite the fact that the η meson is unstable. Its instability only generates a few technical complications in the case of f ηπ + : the discontinuity along 10 Varying only the position of the dip parameter s 1 yields a range B S = [0.17-0.30] × 10 −5 . the unitarity cut is complex and, furthermore, it displays a divergence at the pseudothreshold s = (m η − m π ) 2 . Below 1 GeV, the essential contribution to the discontinuity is proportional to the ηπ → ππ amplitude, projected on the P-wave. We constructed a four-parameter family of solutions of the Khuri-Treiman equations which we use in the dispersion relation for f ηπ + . The shape of the vector form factor, in particular the ρ-meson peak, is then correlated with the Dalitz plot parameters of the η → 3π amplitude (in particular, the parameter d). Upon using the recent experimental constraints for the Dalitz plot, we find the ρ-meson peak to be suppressed as compared with earlier evaluations and that its shape differs from the naive vector dominance model.
In the case of the scalar form factor f ηπ 0 , we used a phase dispersive representation. For the ηπ scattering phase shift in the elastic region, we relied on the ηπ scattering model proposed in Ref. [41]. This model should be reasonable, at the qualitative level, but it is clear that there is much to be improved on our knowledge of ηπ scattering. Again in this case, lattice QCD simulations, which are making steady progress in evaluating meson-meson interactions (see [89] for recent work), could provide unique information, e.g., on the value of the scattering length.
At energies above the KK threshold, we argued that a plausible behaviour for the phase is that it should display a sharp fall-off (which corresponds to a dip in the modulus), by analogy with the cases of the ππ or the π K scalar form factors, where the corresponding phase can be generated from dynamical models. In this approach the exotic (non-exotic) nature of the a 0 (980) resonance corresponds to the dip being situated close (far) from the resonance position. This feature of the phase can then be used in association with a global constraint from a sum rule, which relates the integral over the phase to the logarithmic derivative of the form factor at the origin. Varying the position of the dip generates the main source of uncertainty in this approach. The sum rule restricts the range of variation of the dip but the uncertainties still remain much larger than the 20 % level, required to make the τ → ηπν process competitive for constraining the parameters of particle physics models involving charged Higgs bosons.
The Dalitz plot coefficients parameterise the variation of the square of the amplitude from the centre of the plot The parametrisation accounts for the invariance of the amplitude under the transformation X → −X which results from charge conjugation invariance.
In the case of the decays into three neutral pions one similarly introduces two variables with Q n = m η − 3m π 0 . The amplitude is invariant under Bose symmetry transformations π 0 i ↔ π 0 j . Using Eq. (106), one deduces that it must be invariant under the following transformations of the X , Y variables: The expansion of the amplitude squared around the centre of the Dalitz plot thus has the form, (108) Appendix D: Second-class amplitude in τ → π 0 π + ν decay We remark here that the ππ scalar form factor f ππ 0 [see Eq. (20)], while involving no resonance contribution (to first order in isospin breaking) can be estimated in an essentially model independent way in the low-energy region. The π 0 π + system with l = 0 must be in an isospin I = 2 state. From Watson's theorem, the phase of the scalar form factor φ 0 , must coincide with the l = 0, I = 2 ππ scattering phase shift δ 2 0 in an energy range s < s in 1 GeV 2 where ππ scattering is elastic to a good approximation. We can then express f ππ 0 as a phase dispersive representation, (109) (using f ππ 0 (0) = F π V (0) = 1). At energies s << s in we can neglect the effect of the second exponential in Eq. (109) and thus obtain an approximation of the form factor in terms of the known I = 2 phase shift.
(111) Figure 12 shows that the asymmetry is very small, except, however, in the energy region √ s ≤ 300 MeV, where it is positive and larger that 10 %.