Extended chiral Khuri-Treiman formalism for $\eta\to 3\pi$ and the role of the $a_0(980)$, $f_0(980)$ resonances

TL;DR

This paper extends the Khuri-Treiman formalism to include inelastic effects and scalar resonances in the analysis of $ ext{eta} o 3 ext{pi}$ decays, improving theoretical predictions with recent experimental data.

Contribution

It introduces an extended dispersive approach incorporating $Kar{K}$ channel effects and scalar resonances, enhancing the accuracy of $ ext{eta} o 3 ext{pi}$ decay amplitude modeling.

Findings

Resonances $f_0(980)$ and $a_0(980)$ significantly influence decay amplitudes.

The extended formalism improves the fit to Dalitz plot data.

Inelastic effects are crucial for accurate low-energy decay descriptions.

Abstract

Recent experiments on $\eta\to 3\pi$ decays have provided an extremely precise knowledge of the amplitudes across the Dalitz region which represent stringent constraints on theoretical descriptions. We reconsider an approach in which the low-energy chiral expansion is assumed to be optimally convergent in an unphysical region surrounding the Adler zero, and the amplitude in the physical region is uniquely deduced by an analyticity-based extrapolation using the Khuri-Treiman dispersive formalism. We present an extension of the usual formalism which implements the leading inelastic effects from the $K\bar{K}$ channel in the final-state $\pi\pi$ interaction as well as in the initial-state $\eta\pi$ interaction. The constructed amplitude has an enlarged region of validity and accounts in a realistic way for the influence of the two light scalar resonances $f_0(980)$ and $a_0(980)$ in the dispersive integrals. It is shown that the effect of these resonances in the low energy region of the $\eta \to 3\pi$ decay is not negligible, in particular for the $3\pi^0$ mode, and improves the description of the energy variation across the Dalitz plot. Some remarks are made on the scale dependence and the value of the double quark mass ratio $Q$.

Tables

Table 1. Table 1: Isospin-conserving (cons.) and isospin-violating (viol.) amplitudes involving π ​ π 𝜋 𝜋 \pi\pi , π ​ η 𝜋 𝜂 \pi\eta and K ​ K ¯ 𝐾 ¯ 𝐾 K\bar{K} channels.

	cons.	viol.
$\pi \pi \to \pi \pi$	$f_j^I$	$-$
\hdashline$\eta \pi \to \eta \pi$	$f_j^{\eta \pi }$	$-$
\hdashline$\eta \pi \to \pi \pi$	$-$	${ℳ}_j^I$
\hdashline$\eta \pi \to K\overline{K}$	$g_{j\left(even\right)}^{\eta \pi }$	${𝒩}_{j\left(even\right)}^0,{𝒩}_{j\left(odd\right)}^1$
\hdashline$K\overline{K}\to \pi \pi$	$g_j^I$	${𝒢}_{j\left(even\right)}^{10},{𝒢}_{j\left(even\right)}^{12},{𝒢}_{j\left(odd\right)}^{01}$
\hdashline$K\overline{K}\to K\overline{K}$	$h_j^I$	${\mathscr{H}}_j^{10}$

Table 2. Table 2: Results for the Dalitz plot parameters of the charged and neutral η → 3 ​ π → 𝜂 3 𝜋 \eta\to 3\pi decays based on the NLO chiral amplitude (in the form given in appendix  A ) and its dispersive extrapolations based on single-channel and coupled-channel Khuri-Treiman equations and a matching procedure described in the text. The last two columns show experimental results from refs.  [ 16 , 15 ] and  [ 54 ] .

$\eta \to {\pi }^{+}{\pi }^{-}{\pi }^0$	$O\left(p^4\right)$	single-ch.	coupled-ch.	KLOE	BESIII
a	-1.328	-1.156	-1.142	-1.095(4)	-1.128(15)
b	0.429	0.200	0.172	0.145(6)	0.153(17)
d	0.090	0.095	0.097	0.081(7)	0.085(16)
f	0.017	0.109	0.122	0.141(10)	0.173(28)
g	-0.081	-0.088	-0.089	-0.044(16)	–
$\eta \to {\pi }^0{\pi }^0{\pi }^0$				PDG
$\alpha$	+0.0142	-0.0268	-0.0319	-0.0318(15)
$\beta$	-0.0007	-0.0046	-0.0056	$-$

Table 3. Table 3: Relative variation of the Dalitz plot parameters from their central values. Second column: the η ​ π 𝜂 𝜋 \eta\pi T 𝑇 T -matrix elements are set to zero, third column: the η ​ π 𝜂 𝜋 \eta\pi T 𝑇 T -matrix has larger phase-shifts: δ η ​ π → K ​ K ¯ ​ ( m a 0 ​ ( 1450 ) ) = 180 ​ ° subscript 𝛿 → 𝜂 𝜋 𝐾 ¯ 𝐾 subscript 𝑚 subscript 𝑎 0 1450 180 ° \delta_{\eta\pi\to K\bar{K}}(m_{a_{0}(1450)})=180\degree (instead of 100 ​ ° 100 ° 100\degree ) and the scattering length is a 0 η ​ π = 21.6 ⋅ 10 − 3 superscript subscript 𝑎 0 𝜂 𝜋 ⋅ 21.6 superscript 10 3 a_{0}^{\eta\pi}=21.6\cdot 10^{-3} (instead of 13.2 ⋅ 10 − 3 ) 13.2\cdot 10^{-3}) ), fourth column L 3 = − 2.65 ⋅ 10 − 3 subscript 𝐿 3 ⋅ 2.65 superscript 10 3 L_{3}=-2.65\cdot 10^{-3} (from ref.  [ 55 ] ) instead of L 3 = − 3.82 ⋅ 10 − 3 subscript 𝐿 3 ⋅ 3.82 superscript 10 3 L_{3}=-3.82\cdot 10^{-3} .

	no $\eta \pi$	large $\eta \pi$	$L_3=-2.65\cdot {10}^{-3}$
$\Delta a/\left|a\right|$	$-0.6\%$	$+0.8\%$	$+3.9\%$
$\Delta b/\left|b\right|$	$+9.0\%$	$-9.6\%$	$-2.4\%$
$\Delta d/\left|d\right|$	$-0.7\%$	$+0.8\%$	$-13.0\%$
$\Delta f/\left|f\right|$	$-6.3\%$	$+6.4\%$	$-11.3\%$
$\Delta g/\left|g\right|$	$-0.2\%$	$+0.3\%$	$+10.8\%$
$\Delta \alpha /\left|\alpha \right|$	$+9.1\%$	$-9.2\%$	$+5.5\%$

Equations

$$Q^{-2}=\frac{m_d^2-m_u^2}{m_s^2-((m_u+m_d)/2{)}^2}.$$

$$\begin{array}[]{l}{\cal T}_{\eta\to\pi^{+}\pi^{-}\pi^{0}}(s,t,u)\equiv A(s,t,u)=-\epsilon_{L}\Big{[}M_{0}(s)-{2\over 3}M_{2}(s)+(s-u)M_{1}(t)\\[5.69046pt] \phantom{{\cal T}_{\eta\pi^{0}\to\pi^{+}\pi^{-}}(s,t,u)=-\epsilon_{L}}+(s-t)M_{1}(u)+M_{2}(t)+M_{2}(u)\Big{]}\ ,\end{array}$$

$${ϵ}_L=Q^{-2}\frac{m_K^2-m_{\pi }^2}{3\sqrt{3}{F}_{\pi }^2}\frac{m_K^2}{m_{\pi }^2}.$$

$$s=(p_{\pi^{+}}+p_{\pi^{-}})^{2},\quad t=(p_{\eta}-p_{\pi^{+}})^{2},\quad u=(p_{\eta}-p_{\pi^{-}})^{2}\$$

$$s+t+u=3s_0,s_0=\frac{1}{3}{m}_{\eta }^2+m_{\pi }^2.$$

$${\cal T}_{\eta\to\pi^{0}\pi^{0}\pi^{0}}(s,t,u)=-\epsilon_{L}\Big{[}M_{0}(s)+M_{0}(t)+M_{0}(u)+\frac{4}{3}\big{(}M_{2}(s)+M_{2}(t)+M_{2}(u)\big{)}\Big{]}\ .$$

$$M_1(w)\to M_1(w)+a_1,M_2(w)\to M_2(w)+a_2+b_{2}w$$

$$M_0(w)\to M_0(w)+a_0+b_{0}w$$

$$a_{0}=-\frac{4}{3}a_{2}+3s_{0}(a_{1}-b_{2}),\quad b_{0}=-3a_{1}+\frac{5}{3}b_{2}\$$

$$M_1(0)=0,M_2(0)=0,M_2^{\prime }(0)=0.$$

$$M_0(s)={\tilde{\alpha }}_0+{\tilde{\beta }}_{0}s+\frac{s^2}{\pi }{\int }_{4m_{\pi }^2}^{\infty }d{s}^{\prime }\frac{disc[M_0(s^{\prime })]}{(s^{\prime }{)}^2(s^{\prime }-s)}$$

$$M_1(s)=\frac{s}{\pi }{\int }_{4m_{\pi }^2}^{\infty }d{s}^{\prime }\frac{disc[M_1(s^{\prime })]}{s^{\prime }(s^{\prime }-s)}$$

$$M_2(s)={\tilde{\beta }}_{2}s+\frac{s^2}{\pi }{\int }_{4m_{\pi }^2}^{\infty }d{s}^{\prime }\frac{disc[M_2(s^{\prime })]}{(s^{\prime }{)}^2(s^{\prime }-s)}$$

$$disc[M_I(s)]\equiv \frac{1}{2i}(M_I(s+iϵ)-M_I(s-iϵ)).$$

$$\left(\begin{array}{c}{\pi }^{+}\\ {\pi }^0\\ {\pi }^{-}\end{array}\right)\sim \left(\begin{array}{c}-|11⟩,\\ |10⟩,\\ |1-1⟩\end{array}\right),\left(\begin{array}{c}{K}^{+}\\ K^0\end{array}\right)\sim \left(\begin{array}{c}{\bar{K}}^0\\ -K^{-}\end{array}\right)\sim \left(\begin{array}{c}|\frac{1}{2}\frac{1}{2}⟩\\ |\frac{1}{2}\frac{-1}{2}⟩\end{array}\right).$$

$${\mathcal{M}}^{I,I_z}=⟨\eta \pi |\hat{T}|\pi \pi ;I{I}_z⟩.$$

$$A(s,t,u)=-\frac{1}{\sqrt{3}}{\mathcal{M}}^{0,0}(s,t,u)-\frac{1}{\sqrt{6}}{\mathcal{M}}^{2,0}(s,t,u)$$

$$\begin{array}[]{l}{\cal T}_{\eta\pi^{-}\to\pi^{0}\pi^{-}}=A(t,s,u)={1\over\sqrt{2}}({\cal M}^{1,-1}(s,t,u)+{\cal M}^{2,-1}(s,t,u))\\[5.69046pt] {\cal T}_{\eta\pi^{+}\to\pi^{+}\pi^{0}}=A(u,t,s)={1\over\sqrt{2}}({\cal M}^{1,1}(s,t,u)+{\cal M}^{2,1}(s,t,u))\ .\\ \end{array}$$

$$⟨j^{\prime }{m}^{\prime }|T_q^k|jm⟩=(-1{)}^{j^{\prime }-m^{\prime }}\left(\begin{array}{ccc}{j}^{\prime }& k& j\\ -m^{\prime }& q& m\end{array}\right)⟨j^{\prime }||T^k||j⟩.$$

$$\begin{array}[]{l}{\cal M}^{1,1}=-{\cal M}^{1,-1}\\ {\cal M}^{2,1}={\cal M}^{2,-1}={\sqrt{3}\over 2}{\cal M}^{2,0}\ .\end{array}$$

$$\begin{array}[]{l}{\cal M}^{0,0}(s,t,u)=-\sqrt{3}\left(A(s,t,u)+{1\over 3}(A(t,s,u)+A(u,t,s))\right)\\[5.69046pt] {\cal M}^{1,1}(s,t,u)={1\over\sqrt{2}}\left(-A(t,s,u)+A(u,t,s)\right)\\[5.69046pt] {\cal M}^{2,1}(s,t,u)={1\over\sqrt{2}}\left(A(t,s,u)+A(u,t,s)\right)\ .\\ \end{array}$$

$${\mathcal{M}}^{0,0}\equiv {\mathcal{M}}^0,{\mathcal{M}}^{1,1}\equiv {\mathcal{M}}^1,{\mathcal{M}}^{2,1}\equiv {\mathcal{M}}^2.$$

$$\begin{array}[]{l@{}r@{}l}{\cal M}^{0}(s,t,u)=&\sqrt{3}\,\epsilon_{L}&\left[M_{0}(s)+{1\over 3}M_{0}(t)+{10\over 9}M_{2}(t)+{2\over 3}(s-u)M_{1}(t)+(t\leftrightarrow u)\right]\\[5.69046pt] {\cal M}^{1}(s,t,u)=&{\sqrt{2}\over 3}\,\epsilon_{L}&\left[3tM_{1}(s)+{3\over 2}M_{0}(t)+{3\over 2}(s-u)M_{1}(t)-{5\over 2}M_{2}(t)-(t\leftrightarrow u\right]\\[5.69046pt] {\cal M}^{2}(s,t,u)=&-\sqrt{2}\,\epsilon_{L}&\left[M_{2}(s)+{1\over 2}M_{0}(t)+{1\over 6}M_{2}(t)-{1\over 2}(s-u)M_{1}(t)+(t\leftrightarrow u)\right]\\[5.69046pt] \end{array}$$

$${\mathcal{M}}^I(s,t,u)=16\pi \sqrt{2}\sum _j(2j+1){\mathcal{M}}_j^I(s)P_j(z)$$

$$t,u=\frac{1}{2}\left(m_{\eta }^2+3m_{\pi }^2-s\pm \kappa (s)z\right),\kappa (s)=\sqrt{(1-4m_{\pi }^2/s){\lambda }_{\eta \pi }(s)}$$

$${\lambda }_{PQ}(s)=(s-(m_{\eta }+m_{\pi }{)}^2)(s-(m_{\eta }-m_{\pi }{)}^2).$$

$$\begin{array}[]{l}{\cal M}^{0}_{0}(s)=\dfrac{\sqrt{6}\,\epsilon_{L}}{32\pi}\,[M_{0}(s)+\hat{M}_{0}(s)]\\[8.5359pt] {\cal M}^{1}_{1}(s)=\dfrac{\epsilon_{L}}{48\pi}\,\kappa(s)\,[M_{1}(s)+\hat{M}_{1}(s)]\\[8.5359pt] {\cal M}^{2}_{0}(s)=-\dfrac{\epsilon_{L}}{16\pi}\,[M_{2}(s)+\hat{M}_{2}(s)]\\ \end{array}$$

$$\begin{array}[]{rl}\hat{M}_{0}=&{2\over 3}{\langle M_{0}\rangle}+{20\over 9}{\langle M_{2}\rangle}+2(s-s_{0}){\langle M_{1}\rangle}+{2\over 3}\kappa{\langle zM_{1}\rangle}\\[5.69046pt] \kappa\hat{M}_{1}=&3{\langle zM_{0}\rangle}-5{\langle zM_{2}\rangle}+{9\over 2}(s-s_{0}){\langle zM_{1}\rangle}+{3\over 2}\kappa{\langle z^{2}M_{1}\rangle}\\[5.69046pt] \hat{M}_{2}=&{\langle M_{0}\rangle}+{1\over 3}{\langle M_{2}\rangle}-{3\over 2}(s-s_{0}){\langle M_{1}\rangle}-{1\over 2}\kappa{\langle zM_{1}\rangle}\ .\\ \end{array}$$

$$⟨z^n{M}_I⟩(s)=\frac{1}{2}{\int }_{-1}^{1}dz{z}^n{M}_I(t(s,z)).$$

$$\mathrm{Im}{\mathcal{M}}_j^I(s)=disc[{\mathcal{M}}_j^I(s)]={\sigma }_{\pi }(s)(f_j^I(s){)}^{\ast }{\mathcal{M}}_j^I(s).$$

$$\sigma_{P}(s)=\sqrt{1-\frac{4m_{P}^{2}}{s}}\,\theta(s-4m_{P}^{2})\$$

Full text

1. 1 Introduction

2. 2 Khuri-Treiman equations in the elastic approximation

3. 2.1 Single-variable amplitudes for $\eta\to 3\pi$

4. 2.2 Isospin amplitudes and crossing relations

5. 2.3 Partial-waves and elastic $\pi\pi$ unitarity relations

6. 2.4 Khuri-Treiman equations in the elastic approximation

7. 3 Beyond the elastic $\pi\pi$ approximation

8. 3.1 $\eta\pi$ contribution to unitarity

9. 3.2 $K\bar{K}$ contributions to unitarity

10. 4 Multichannel Khuri-Treiman equations

11. 4.1 Closed system of unitarity equations

12. 4.2 Coupled Muskhelishvili-Omnès representation

13. 4.3 Matching to the chiral amplitudes

14. 5 Results and comparisons with experiment

15. 5.1 Numerical method

16. 5.2 Input $I=0,1$ T-matrices

17. 5.3 Illustration of the role of the inelastic channels

18. 5.4 Dalitz plot parameters

19. 5.5 The ratio $\Gamma(\eta\to 3\pi^{0})/\Gamma(\eta\to{\pi^{+}}{\pi^{-}}{\pi^{0}})$

20. 5.6 The quark mass ratio $Q$ from the chirally matched dispersive amplitude

21. 5.7 Further experimental constraints on $Q$

22. 6 Conclusions

23. A The amplitude $\eta\to{\pi^{+}}{\pi^{-}}{\pi^{0}}$ at chiral order $p^{4}$ and $e^{2}p^{2}$

24. B Matching equations

25. C Angular integrals and hat functions

26. D Supplementary material

1 Introduction

The physics of QCD in the soft regime is dominated by the phenomenon of spontaneous symmetry breaking because of the presence of three light quarks in the standard model. The low-energy dynamics can then be described accurately through an expansion built from a chiral effective theory (e.g. [1] for a recent review). This approach, which applies in both the Euclidean and Minkowski space-times is, to some extent, complementary to the purely numerical lattice simulation method. In the effective theory, however, part of the information on the non-perturbative QCD dynamics is contained as sets of values of the chiral coupling constants. These are not known, a priori, except for their order of magnitude [2] and must be determined as part of the probing of the effective theory.

In QCD, isospin breaking phenomena are driven by $m_{d}-m_{u}$, the mass difference of the two lightest quarks. For low-energy observables, an isospin breaking ratio which conveniently absorbs some of the next-to-leading (NLO) chiral coupling constants was introduced in ref. [3]

[TABLE]

In general, isospin breaking effects induced by electromagnetism are comparable in size to those proportional to $m_{d}-m_{u}$ and their precise evaluation is made difficult by the poor knowledge of the associated chiral coupling constants [4]. In this respect, the $\eta\to 3\pi$ amplitude plays a special role because these electromagnetic contributions vanish in the $SU(2)$ chiral limit [5] and are thus expected to be suppressed. This has been confirmed in the work of refs. [6, 7] who evaluated the contributions proportional to $e^{2}m_{u}$, $e^{2}m_{d}$.

A number of recent high-statistics experiments have studied the $3\pi^{0}$ decay mode of the $\eta$ [8, 9, 10, 11, 12] as well as the charged mode $\pi^{+}\pi^{-}\pi^{0}$ [13, 14, 15, 16]. An extremely precise knowledge of the energy variation of the amplitudes squared across the Dalitz plot, which are traditionally represented by a set of polynomial parameters, has now become available. These accurate experimental results allow for stringent tests of the theoretical description of the amplitude which must obviously be passed before one attempts to determine $Q$. The Dalitz plot parameters derived directly from the NLO chiral amplitude, which was first computed in [17], are in clear disagreement with experiment. For instance, the prediction for the parameter $\alpha$, involved in the neutral mode, has the wrong sign. The same problems, essentially, are found in the resummed expansion approach discussed recently in ref. [18]. The computation of the amplitude at the next-to-next-to leading (NNLO) chiral order was performed [19]. The comparison with experiment again fails if one assumes a simple naive model for the $O(p^{6})$ couplings $C_{i}$ (classified in ref. [20]) which are involved as six independent combinations. The $\eta$ decay amplitude thus contains crucial information on the true QCD values of these couplings, which are essentially not known at present.

One obvious deficiency of the chiral expansion when calculating scattering or decay amplitudes in physical regions is the lack of exact unitarity (as emphasised e.g. in ref. [21]) which is restored gradually when going to higher orders. In the case of $\eta\to 3\pi$, a large amount of work was devoted to the problem of estimating these unitarity, or final-state rescattering, higher order corrections [22, 23, 24, 25, 26, 27, 28].

In the present paper we reconsider, more specifically, the approach followed in refs. [22, 24, 25]. The main underlying assumption is that the chiral expansion of the $\eta\to{\pi^{+}}{\pi^{-}}{\pi^{0}}$ decay amplitude should be optimally converging in an unphysical region of the Mandelstam plane in the neighbourhood of the Adler zero. The amplitude in the physical region is then deduced from a well-defined extrapolation procedure based on the analyticity properties of amplitudes in QCD, which utilise the set of dispersive equations derived initially by Khuri and Treiman [29] and perfected in the work of refs. [22, 24, 25]. These equations implement crossing-symmetry and unitarity in a more complete way than more naive loop-resummation approaches [26].

In previous work, $\pi\pi$ rescattering was assumed to be elastic. This is essentially exact in the physical $\eta$ decay region. However, the dispersive formalism involves integrals over an energy range extending up to infinity. A property of the $\pi\pi$ scattering amplitude in the isoscalar S-wave is the sharp onset of $K\bar{K}$ inelasticity associated with the $f_{0}(980)$ scalar resonance [30, 31]. Because of isospin violation, the $\eta\pi\to\pi\pi$ amplitude actually exhibits a double resonance effect from both the $f_{0}(980)$ and the $a_{0}(980)$ scalars [32] near the $K\bar{K}$ threshold. Our aim is to propose a generalisation of the Khuri-Treiman formalism which takes into account $K\bar{K}$ inelasticity in the unitarity relations for both $\pi\pi$ scattering and $\eta\pi$ scattering (which may be viewed as an initial-state interaction). Some approximations will be made, which simplify considerably the practical implementation, such that crossing-symmetry is maintained at the level of the $\eta\to 3\pi$ amplitudes but not in the amplitudes involving the $K\bar{K}$ channel. In this multichannel formalism, the double resonance effect is taken into account in the dispersive integrals and, furthermore, the construction of the $\eta\to 3\pi$ amplitude becomes valid in an extended energy region which includes not only the $\eta$ decay region but also a portion of the $\eta\pi\to\pi\pi$ scattering region. This will allow us to study how the energy dependence induced by the two 1 GeV scalar resonances propagate down to the low energy region and quantitatively affects the Dalitz plot parameters. The fact that these resonances could be influential at low energy was pointed out previously in ref. [33].

The plan of the paper is as follows. We first review in sec. 2 the derivation of the one-channel equations in which the amplitudes satisfy elastic $\pi\pi$ unitarity in the $S$ and the $P$-waves. In sec. 3 we write the unitarity relations including the $\eta\pi$ and the $K\bar{K}$ channels. A closed system of unitarity equations involves, besides $\eta\to 3\pi$, isospin violating components of $\eta\pi\to K\bar{K}$, $\pi\pi\to K\bar{K}$ as well as $K\bar{K}\to K\bar{K}$ amplitudes. A multichannel set of Khuri-Treiman integral equations is defined such that the solution amplitudes satisfy these unitarity relations. We then discuss in sec. 4 the matching between the chiral expansion amplitudes and the dispersive ones. We adopt a simple approach which consists in imposing that the difference between the chiral NLO and dispersive amplitudes vanishes at order $p^{4}$. This provides four equations which, in the single-channel case determine completely the dispersive amplitude provided one had introduced exactly four polynomial parameters in the Khuri-Treiman representation. This is generalised to the multichannel situation, in which one introduces 16 polynomial parameters. Finally, in sec. 5, the results on the Dalitz plot parameters are presented and some remarks are made on the determination of the quark mass ratio $Q$.

2 Khuri-Treiman equations in the elastic approximation

We remind below how the dispersion relations-based equations derived by Khuri and Treiman [29] for $K\to 3\pi$ decay can be generalised and applied to $\eta\to 3\pi$, following [22, 24, 25].

2.1 Single-variable amplitudes for $\eta\to 3\pi$

As initially demonstrated in the case of $\pi\pi\to\pi\pi$ (see [34]), amplitudes involving four pseudo-Goldstone bosons satisfy, in a certain kinematical range, an approximate representation in terms of functions of a single variable which have simple analyticity properties (see [35] for a review). In the case of $\eta\to 3\pi$, and neglecting quadratic isospin breaking, there are three functions involved [24, 25]: $M_{0}$, $M_{1}$, $M_{2}$. We will follow the notation of ref. [25] and write the $\eta\to\pi^{+}\pi^{-}\pi^{0}$ amplitude as

[TABLE]

with an overall factor111It is convenient to formally factor out $\epsilon_{L}$ but the amplitude is actually of the form: ${\cal T}_{\eta\to 3\pi}=\epsilon_{L}\,{\cal A}+\Delta{m}^{2}_{K}\,{\cal B}+e^{2}\,{\cal C}$ where $e$ is the electric charge and $\Delta{m}^{2}_{K}$ is the physical $K^{0}-K^{+}$ mass squared difference (see appendix A). $\epsilon_{L}$ which is proportional to the isospin breaking double quark mass ratio $Q^{2}$ given in eq. (1)

[TABLE]

The Mandelstam variables are defined, as usual, as

[TABLE]

and satisfy

[TABLE]

General analyticity properties imply that $\eta\to 3\pi$ decay and $\eta\pi\to\pi\pi$ scattering are described by the same function in different regions of the Mandelstam plane. The corresponding physical regions are illustrated in fig. 1. The analogous representation for $\eta\to 3\pi^{0}$ involves the two functions $M_{0}$ and $M_{2}$ only and reads

[TABLE]

The representations (2) and (6) are accurate in regions of the Mandelstam plane where the imaginary parts of the partial-wave amplitudes with angular momentum $j\geq 2$ in the $s$, $t$ or $u$ channels are negligible (compared to those of the $j=0,1$ partial-waves). In the case of $\eta\to 3\pi$ or $\eta\pi\to\pi\pi$, this condition is satisfied in the range where $s$, $t$, $u$ are sufficiently small compared with the masses squared of the tensor resonances, i.e. |s|,|t|,|u|\mathrel{\raise 2.32501pt\hbox{<}\kern-7.00002pt\lower 2.5187pt\hbox{\sim}}1 $\hbox{GeV}^{2}$. This condition is also satisfied exactly by the amplitude obtained from the chiral expansion up to order $p^{6}$ [19]. This will prove very useful for writing matching conditions.

The functions $M_{I}(w)$ are analytic in $w$ with a cut on the positive real axis: $4m_{\pi}^{2}\leq w<\infty$. Based on Regge theory, we expect that the functions $M_{0}(w)$ and $M_{2}(w)$ should not grow faster than $w$ at infinity, while $M_{1}$ should be bounded by a constant. In the one-channel Khuri-Treiman framework, both $M_{0}$ and $M_{2}$ are usually assumed to behave linearly in $w$ when $w\to\infty$. With this asymptotic behaviour, there is a family of redefinitions of the functions $M_{1}$ and $M_{2}$ which leaves the physical amplitude $A(s,t,u)$ unmodified [25, 36],

[TABLE]

($a_{1}$, $a_{2}$, $b_{1}$ being arbitrary constant parameters) provided one correspondingly redefines $M_{0}$ as

[TABLE]

with

[TABLE]

and using the $s+t+u$ constraint (5). This arbitrariness can be fixed by imposing the three $w=0$ conditions

[TABLE]

In the coupled-channel set-up, to be discussed below, the asymptotic condition on $M_{0}$ will be modified such that $M_{0}(w)$ goes to a constant when $w\to\infty$ instead of behaving linearly. This restricts the allowed redefinitions of $M_{1}$, $M_{2}$ to those which satisfy $a_{1}=5/9\,b_{2}$. In that case, only the two $w=0$ conditions $M_{1}(0)=M_{2}(0)=0$ can be imposed, while $M^{\prime}_{2}(0)$ is determined from the equations.

These properties of the functions $M_{I}$ lead to the following dispersive representations

[TABLE]

defining

[TABLE]

2.2 Isospin amplitudes and crossing relations

We choose the following conventional isospin assignment for the pions and the kaons

[TABLE]

Let us consider the amplitudes which correspond to isospin states of the $\pi\pi$ system,

[TABLE]

One can express ${\cal T}_{\eta\pi^{0}\to\pi^{+}\pi^{-}}=A(s,t,u)$ in terms of $I=0,2$ isospin amplitudes,

[TABLE]

Further relations for the isospin amplitudes can be obtained using crossing symmetries. Under $s-t$ and $s-u$ crossing one obtains,

[TABLE]

Since the isospin breaking operator in QCD, $H_{IB}=-1/2(m_{d}-m_{u})\bar{\psi}\lambda_{3}\psi$, transforms as $I=1$, $I_{z}=0$, one can use the Wigner-Eckart theorem

[TABLE]

which yields the following relations among the ${\cal M}^{I,I_{z}}$ amplitudes

[TABLE]

One can then express the three independent isospin amplitudes in terms of the function $A(s,t,u)$

[TABLE]

In the following we will simply denote

[TABLE]

Inserting the representation (2) we obtain an expression of the three isospin amplitudes in terms of the one-variable functions $M_{I}(w)$

[TABLE]

2.3 Partial-waves and elastic $\pi\pi$ unitarity relations

In order to derive expressions for the discontinuities $\hbox{disc}[M_{I}(s)]$, we must consider partial-waves and their unitarity relations. We can define the partial-wave expansion of the isospin amplitudes as

[TABLE]

where $z$ is the cosine of the scattering angle in the centre-of-mass frame of $\eta\pi\to\pi\pi$ which is related to the Mandelstam variables by

[TABLE]

with

[TABLE]

From the representation of the isospin amplitudes, eq. (21) one easily derives the expression for the partial-waves. The result, for the $j=0,1$ partial-waves, can be written as,

[TABLE]

where the functions $\hat{M}_{I}(s)$ are given by linear combinations of angular integrals of the functions $M_{I}$

[TABLE]

with the notation [25]

[TABLE]

Writing unitarity relations, one must first consider the unphysical situation where the $\eta$ meson is stable, $m_{\eta}\leq 3m_{\pi}$. The physical case is defined using analytic continuation in $m_{\eta}$ as in the classic derivation of generalised unitarity [37]. The contribution from the $\pi\pi$ states to the unitarity relation for the partial-wave ${\cal M}_{j}^{I}$, reads

[TABLE]

with

[TABLE]

and $f_{j}^{I}(s)$ is the $\pi\pi\to\pi\pi$ partial-wave amplitude, which is related to the scattering phase-shift by

[TABLE]

The equality between the imaginary part and the discontinuity in eq. (28) holds when $m_{\eta}<3m_{\pi}$. For the physical value of $m_{\eta}$ the right-hand side of eq. (28) continues to give the discontinuity across the unitarity cut, while the imaginary part must be deduced from the dispersive representation.

2.4 Khuri-Treiman equations in the elastic approximation

In the unphysical situation when $m_{\eta}<3m_{\pi}$, the cuts of the functions $\hat{M}_{I}(w)$ are located in the region $Re[w]<4m_{\pi}^{2}$ such that eqs. (25) correspond to a splitting of the partial wave amplitudes into two functions which have a separated cut structure. When the $\eta$ mass is increased to its physical value, the $m_{\eta}^{2}+i\epsilon$ prescription must be used [37] and this insures that the complex cut of the $\hat{M}_{I}$ functions, which approaches infinitesimally close to the unitarity cut in the region $4m_{\pi}^{2}\leq s\leq(m_{\eta}-m_{\pi})^{2}$, remains well separated from it (see fig. 4 in ref. [38]). Using the fact that $\hat{M}_{I}(s)$ has no discontinuity across the unitarity cut one can deduce from (28) that the discontinuities of the functions $M_{I}(s)$ along $4m_{\pi}^{2}\leq s<\infty$ are given by

[TABLE]

where $j=0$ when $I=0,2$ and $j=1$ when $I=1$. In the sequel, we will drop the $j$ subscript in the $\pi\pi$ phase-shift. Eqs. (31) imply that the functions $M_{I}$ can be written as Muskhelishvili-Omnès (MO) integral representations

[TABLE]

where

[TABLE]

and where the Omnès functions $\Omega_{I}$ are given in terms of the $\pi\pi$ phase-shifts by the usual relation

[TABLE]

The phase-shifts, in the present context, are usually taken to obey the following asymptotic conditions

[TABLE]

which seem rather natural since these conditions are roughly satisfied at the $K\bar{K}$ threshold. In the elastic approximation framework, one can thus take the phases to be constant or quasi-constant in the inelastic region, above 1 GeV.

The polynomial part in the MO representation (2.4) was chosen to have four parameters. This allows one to define a unique dispersive amplitude by implementing four independent matching conditions with the chiral NLO amplitude. Taking into account the asymptotic conditions on the phase-shifts (35) one easily sees that the Khuri-Treiman equations (2.4) implement the following asymptotic behaviour for the functions $M_{I}$

[TABLE]

The functions $M_{I}$ also satisfy the $w=0$ conditions (10) and are thus uniquely defined.

3 Beyond the elastic $\pi\pi$ approximation

Elastic unitarity for $\pi\pi$ scattering is valid exactly below the four pions threshold and approximately up to the $K\bar{K}$ threshold. Close to 1 GeV, the $\pi\pi$ phase-shift increases very sharply under the influence of the $f_{0}(980)$ resonance which also couples strongly to $K\bar{K}$. In order to properly account for the effect of this resonance it is thus necessary to go beyond the elastic unitarity approximation. We discuss in this section how to include both the $\eta\pi$ and the $K\bar{K}$ channels into the unitarity relations and then generate a generalisation of the Khuri-Treiman equations. This will allow us to account for both the $f_{0}(980)$ and the $a_{0}(980)$ resonances in a realistic way.

3.1 $\eta\pi$ contribution to unitarity

Including the $\eta\pi$ channel in addition to $\pi\pi$, the unitarity relation becomes

[TABLE]

where

[TABLE]

and ${f}_{j}^{\eta\pi}(s)$ is the $\eta\pi\to\eta\pi$ partial-wave amplitude. In the energy region where $\eta\pi$ scattering is elastic, which we will assume to extend up to the $K\bar{K}$ threshold, ${f}_{j}^{\eta\pi}(s)$ is related to the scattering phase-shift by

[TABLE]

The $j=1$ partial-wave ${f}_{1}^{\eta\pi}$ corresponds to exotic quantum numbers $j^{PC}=1^{-+}$ and should thus remain rather small up to the 1 GeV region222A resonance possibly exists in this amplitude [39] with a mass $M\simeq 1.3$ GeV. Therefore, $\eta\pi$ rescattering is expected to affect mainly the two $j=0$ amplitudes ${\cal M}_{0}^{0}$ and ${\cal M}_{0}^{2}$. In the elastic regime, using eq. (37) one easily derives that the relation between the amplitudes on both sides of the unitarity cut reads,

[TABLE]

Comparing with the analogous relation in the elastic unitarity case (28) one deduces that including the effect of $\eta\pi$ rescattering in the $\eta\pi\to\pi\pi$ amplitude (which can be viewed as an initial-state interaction) amounts to simply perform the following replacements in the Omnès representations (2.4),

[TABLE]

In practice, $\eta\pi$ rescattering is expected to become significant when the energy approaches the mass of the $a_{0}(980)$ resonance. It becomes necessary, then, to also take into account the $K\bar{K}$ channel.

3.2 $K\bar{K}$ contributions to unitarity

Let us now include the $K\bar{K}$ states into the partial-wave unitarity relations.

• a)

$I=1,j=1$: We are concerned mainly with $j=0$ amplitudes but let us consider the $j=1$ amplitude ${\cal M}_{1}^{1}$ here also for completeness. The $K\bar{K}$ contribution reads

[TABLE]

The amplitude $T_{1}^{\eta\pi^{+}\to{K^{+}}{\bar{K}^{0}}}$ is isospin violating333 The amplitudes $T_{j}^{\eta\pi^{+}\to{K^{+}}{\bar{K}^{0}}}$ are isospin violating (conserving) for odd (even) values of $j$. This can be seen using G-parity: $G|(K\bar{K})^{I}_{j}\rangle=(-1)^{I+j}|(K\bar{K})^{I}_{j}\rangle$. Since $I=1$ for $K^{+}{\bar{K}^{0}}$, $G=+1$ for odd values of $j$ and $-1$ for even values while $G=-1$ for $\eta\pi$. and, at linear order in isospin breaking, one can set $\sigma_{K^{+}K^{0}}(s)=\sigma_{K}(s)$. We will denote the amplitudes appearing above as

[TABLE]

• b)

$I=2,j=0$: For the ${\cal M}_{0}^{2}$ amplitude now, the $K\bar{K}$ contribution reads

[TABLE]

The amplitude $T_{0}^{\eta\pi^{+}\to{K^{+}}{K^{0}}}$ is isospin conserving in this case, since $j$ is even while the amplitude $T_{0}^{{K^{+}}{\bar{K}^{0}}\to\pi^{+}\pi^{0}}$ is isospin violating. We will denote the two amplitudes in eq. (44) as

[TABLE]

• c)

$I=0,j=0$: Finally, for the ${\cal M}_{0}^{0}$ amplitude, the $K\bar{K}$ contributions to the unitarity relations are

[TABLE]

Let us now separate the isospin conserving and the isospin violating contributions. For the kinematical factors, we introduce

[TABLE]

For the $\eta{\pi^{0}}\to{K^{+}}{K^{-}},\,{K^{0}}{\bar{K}^{0}}$ amplitudes, the isospin conserving part, $g^{\eta\pi}_{0}$, was introduced in eq. (45) and we call the isospin violating one ${\cal N}^{0}_{0}$. One has

[TABLE]

For the $K\bar{K}\to(\pi\pi)^{0}$ amplitudes, the isospin conserving and the isospin violating amplitudes are denoted as

[TABLE]

Using this notation, the $K\bar{K}$ contributions in the unitarity relations for the partial-waves ${\cal M}^{I}_{j}$ can be summarised as follows

[TABLE]

These contributions involve new isospin-breaking $K\bar{K}\to\pi\pi$ and $\eta\pi\to K\bar{K}$ amplitudes: ${\cal G}^{10}_{0}$, ${\cal G}^{12}_{0}$, ${\cal N}^{0}_{0}$, ${\cal N}^{1}_{1}$. In order to write a closed set of unitarity equations we must also consider $K\bar{K}\to K\bar{K}$ amplitudes. For these, the isospin conserving/violating components are denoted $h_{j}^{I}$, ${\cal H}^{10}_{j}$

[TABLE]

Table 1 summarises our notation for the various amplitudes involved.

4 Multichannel Khuri-Treiman equations

4.1 Closed system of unitarity equations

We can now write down a closed system of unitarity equations. It will be convenient to introduce a matrix notation for the isospin conserving amplitudes with $I=0$ and $I=1$

[TABLE]

The $I=0$ amplitude ${\cal M}_{0}^{0}$ is now embedded into a system of four coupled unitarity equations

[TABLE]

where

[TABLE]

while the $I=2$ amplitude ${\cal M}_{0}^{2}$ is involved in a system of two unitarity equations,

[TABLE]

Finally, for the $I=1$ amplitude ${\cal M}_{1}^{1}$, the coupled unitarity equations read

[TABLE]

In the following, however, we will disregard the inelasticity effects for ${\cal M}_{1}^{1}$ and continue to use the elastic unitarity equation (28) in this case.

4.2 Coupled Muskhelishvili-Omnès representation

The next step is to write each one of the isospin violating partial-wave amplitudes as a sum of two functions, one having a right-hand cut only and one having a generalised left-hand cut. For physical mass values, the left and right-hand cuts of the various partial-waves appear to be overlapping, but it is possible to separate them unambiguously. In the case of the amplitudes involving the $\eta$ meson, this is done by using the $m_{\eta}^{2}+i\epsilon$ prescription. The amplitude $K\bar{K}\to K\bar{K}$ has a left-hand cut which extends on the real axis in the range $[-\infty,4m_{K}^{2}-4m_{\pi}^{2}]$. Using the $m_{K}^{2}+i\epsilon$ prescription shifts this cut above the real axis.

For the $I=0$ amplitudes one writes

[TABLE]

which generalises eq. (25) while for the $I=2$ amplitudes one can write

[TABLE]

One can now employ the standard Omnès method in order to express the right-cut functions in terms of the left-cut ones. Introducing the matrix notation

[TABLE]

the discontinuity relation for the $\bm{M}_{0}$ functions is deduced from the unitarity relation (53),

[TABLE]

where

[TABLE]

and $\Delta\sigma_{K}$ is given in eq. (47). Eq. (60) generalises the one-channel discontinuity relation (31).

Let us now consider the following matrix

[TABLE]

where $\bm{\Omega}_{I}$ are the $2\times 2$ Muskhelishvili-Omnès matrices corresponding to the $T$-matrices $\bm{T}_{I}$. Making use of the following discontinuity properties of the MO matrices

[TABLE]

one can express the discontinuity of the $\bm{X}$ matrix elements in terms of the $\hat{\bm{M}}_{0}$ functions and $\Delta\sigma_{K}$

[TABLE]

where

[TABLE]

An alternative expression for $\Delta\bm{X}_{a}$ can be derived, using eqs. (4.2))

[TABLE]

which shows that it represents the discontinuity of the following quantity

[TABLE]

across the right-hand cut. The quantity $\Delta\bm{X}_{b}$ is proportional to $\Delta\sigma_{K}$ and it is given by

[TABLE]

We can then write a twice subtracted dispersive representation for $\bm{X}(s)$ and generate, via (62), a MO representation for the $\bm{M}_{0}$ amplitudes

[TABLE]

where $\bm{P}_{0}$ is a $2\times 2$ matrix of polynomial functions,

[TABLE]

and the integral parts are

[TABLE]

One remarks that in the term $\Delta\bm{X}_{b}$ the quark mass ratio $\epsilon_{L}$ appears in the denominator and thus cancels with the overall factor in the complete amplitudes. This part is driven by the physical $K^{+}-K^{0}$ mass difference via the $\Delta\sigma_{K}$ function (see (47)). This mechanism was first studied in ref. [32] who predicted that a large isospin violation should take place at 1 GeV in the $\eta\pi\to\pi\pi$ scattering amplitude, due to the contributions from both the $a_{0}(980)$ and the $f_{0}(980)$ resonances. The set of eqs. (69) account for the other sources of isospin violation as well.

We now consider the case of the $I=2$ amplitudes and we also introduce a matrix notation for the column matrices

[TABLE]

From the unitarity relations including the $\pi\pi$, $\pi\eta$ and $K\bar{K}$ contributions we deduce the discontinuity of the $\bm{M}_{2}$ functions

[TABLE]

As before, we introduce a matrix $\bm{X}_{2}$, multiplying $\bm{M}_{2}$ by inverse MO functions,

[TABLE]

Its discontinuity relation is expressed in terms of the $\hat{\bm{M}}_{2}$ functions and the MO functions

[TABLE]

An alternative useful expression for $\Delta\bm{X}_{2}$ can be derived

[TABLE]

This leads to the following integral MO representation for the matrix $\bm{M}_{2}$,

[TABLE]

with

[TABLE]

Eqs. (69) and  (77) together with the uncoupled equation for $M_{1}$ (2.4) involve 16 polynomial parameters. The polynomial dependence was chosen such that the equations would reduce exactly to the set of elastic equations (2.4) if one switches off the coupling to the $K\bar{K}$ channel as well as $\eta\pi$ rescattering444The one-channel case is recovered by setting $(\bm{T}^{1})_{ij}=0$ and the MO matrices to $[\bm{\Omega}_{0}]_{ij}=\delta_{ij}(1+\delta_{i1}\Omega_{0})$ and $[\bm{\Omega}_{1}]_{ij}=\delta_{ij}$.. One remarks that the asymptotic behaviour in the coupled-channel equations is modified as compared to the one-channel case: since the matrix elements of the MO matrices $\bm{\Omega}_{I}$, decrease as $1/s$ when $s$ goes to $\infty$ the entries of the $\bm{M}_{0}$ matrix (thus the $M_{0}$ function) behave as constants. The asymptotic behaviour of $M_{2}$, in contrast, remains the same as before.

In addition to the polynomial parameters, the right-hand side of eqs. (69), (77) involve a number of “hat functions”. The functions $\hat{M}_{I}$ are determined in terms of the functions $M_{I}$ by the angular integrals (26), but one still needs to determine the other hat functions $\hat{G}_{10}$, $\hat{N}_{0}$, $\hat{H}_{10}$, $\hat{G}_{12}$ which are related to the $K\bar{K}$ amplitudes. For this purpose, one would have to consider all the related crossed channel amplitudes an write similar sets of equations (which would introduce further one-variable functions). Here, since we are mainly interested in the $\eta\pi\to\pi\pi$ amplitude we will content with an approximation for the amplitudes involving the $K\bar{K}$ channel, simply neglecting the integrals which involve the left-cut functions, i.e. we take

[TABLE]

In support of this approximation, one observes that if one were to neglect the left-cut integrals in the $\eta\to 3\pi$ amplitude itself, one would still obtain a qualitatively reasonable description (e.g. [40]). With this approximation, eqs. (69) and  (77) constitute a closed set of equations which form a coupled-channel generalisation of the Khuri-Treiman equations (2.4) for $M_{0}$ and $M_{2}$.

4.3 Matching to the chiral amplitudes

We intend to fix the 16 polynomial parameters by matching to the chiral expansions of the amplitudes involved. For the $\eta\to 3\pi$ amplitude we will use the NLO expansion also including the part of the electromagnetic contributions of order $e^{2}(m_{d}+m_{u})$ from ref. [6]. This makes it possible to display explicitly the term induced by the $K^{+}-K^{0}$ mass difference via unitarity which, in the dispersive representation, is contained in the $\hat{\bm{I}}_{b}$ integrals (see eqs. (68), (71)). The explicit chiral expressions for the functions $M_{I}$ as used here are given in appendix A.

For the isospin violating amplitudes involving the $K\bar{K}$ channel we will use the leading order chiral expansion. At this order, the partial-wave amplitudes have no left-hand cut, which is consistent with the approximation of dropping the left-cut functions in the integral equations. The relevant expressions (including the $O(e^{2})$ contributions) are given below

[TABLE]

with

[TABLE]

and the coupling constant $C$ can be related to the $\pi^{+}-\pi^{0}$ mass difference (see appendix A).

We implement matching conditions for the $\eta\to 3\pi$ amplitude following the simple method of ref. [38] which differs slightly from that of ref. [25]. Let $\bar{M}(s,t,u)$ be the amplitude computed in the chiral expansion at order $p^{4}$. The polynomial parameters of the dispersive amplitude must be determined such that the dispersive and chiral amplitudes coincide for small values of the Mandelstam variables. At order $p^{4}$ one should thus have

[TABLE]

This condition is satisfied automatically for the discontinuities, which implies that one can neglect the discontinuity of the differences of the one-variable functions $M_{I}-\bar{M}_{I}$ and thus expand these differences as polynomials,

[TABLE]

(with $n_{0}=n_{2}=2$, $n_{1}=1$). Inserting these expansions into the amplitude difference $M(s,t,u)-\bar{M}(s,t,u)$ and requiring that the $O(p^{4})$ polynomial part vanishes gives four equations. In the elastic Khuri-Treiman framework these determine the four parameters via the following four equations

[TABLE]

with

[TABLE]

and one must keep in mind that the integrals $\hat{I}_{I}(0)$ which appear on the right-hand sides depend linearly on the four parameters. In ref. [25] the first two of eqs. (4.3) were replaced by two equations related to the position $s_{A}$ of the Adler zero of the chiral amplitude $\bar{M}(s,t,u)$ along the line $u=s$ and the value of its derivative at $s=s_{A}$. We will see below that eqs. (4.3) do actually implement these Adler zero conditions to a rather good approximation. Additionally, approximations were made in ref. [25] in the determination of $\beta_{1}$ and $\gamma_{0}$ (yielding e.g. $\gamma_{0}\simeq 0$), the validity of which depends on the assumed behaviour of the phase $\delta^{0}$ in the inelastic region. Here, the four equations will be solved without approximation. Doing so, one notes that the polynomial parameters get an imaginary part from the contributions of the integrals $\hat{I}_{I}(0)$ which, however, is small (less that 10% of the real part).

In the coupled-channel case, the matching of the $\eta\to 3\pi$ amplitude again gives rise to four equations. In addition, we can match the values at $w=0$ of each one of the four $K\bar{K}$ amplitudes with the chiral ones given in eqs. (4.3) as well as the values of the first and second derivatives (which are vanishing). Altogether, this provides 16 constraints which fix all the polynomial parameters appearing in the coupled-channel Khuri-Treiman equations. The details of these matching equations are given in appendix B.

5 Results and comparisons with experiment

5.1 Numerical method

The main difficulty which is involved in deriving accurate numerical solutions of the Khuri-Treiman equations is tied to the evaluation of the angular integrals ${\langle z^{n}M_{I}\rangle}$ needed to obtain the hat functions $\hat{M}_{I}$ and to the treatment of the singularities of these functions in the computation of the $\hat{I}_{a}$ integrals which are finite. These technical aspects are discussed in detail in ref. [24]. By a change of variables one can rewrite the angular integrals as integrals over $t$

[TABLE]

with

[TABLE]

When $s$ lies in the range $(m_{\eta}-m_{\pi})^{2}<s<(m_{\eta}+m_{\pi})^{2}$, the endpoints $t^{\pm}(s)$ become complex. In fact, using the $m_{\eta}^{2}+i\epsilon$ prescription one sees that $t^{+}$ and $t^{-}$ are placed on opposite sides of the unitarity cut when $s$ gets larger that $(m_{\eta}^{2}-m_{\pi}^{2})/2$. The integral from $t^{-}$ to $t^{+}$ must then be evaluated along a complex contour which circles around the unitarity cut of the functions $M_{I}$. Rather than computing explicitly $M_{I}(w)$ for complex values of $w$, as in ref. [24], we follow here the approach of ref. [38] which consists in inserting the dispersive representations (2.1) of the $M_{I}$ functions into eq. (86) and computing analytically the $t$ integrals. This makes it possible to express the functions ${\langle z^{n}M_{I}\rangle}$ in terms of the discontinuities $\hbox{disc}[M_{I}]$ along the positive real axis. The relevant formulae are recalled in appendix C.

The equations are conveniently solved using an iterative procedure. On the first iteration step, one sets the $\hat{M}_{I}$ functions equal to zero. Then the $\hat{I}$ integrals are also equal to zero and it is straightforward to compute the values of the polynomial parameters from the matching equations and then the functions $M_{I}$, $G_{10}$, $N_{0}$, $H_{10}$, $G_{12}$ as well as the discontinuities $\hbox{disc}[M_{I}]$ (which are given from eqs. (31) in the one-channel case and  (60) (73) in the coupled-channel case). The coupled-channel framework is somewhat more complicated to handle than the single-channel one, essentially because the MO matrices do not obey a simple explicit representation in terms of the $T$-matrix elements and must be solved for numerically, but does not otherwise involve any specific difficulty. Then, on each iteration step, one updates the values of the functions $\hat{M}_{I}$ using $\hbox{disc}[M_{I}]$ from the preceding step, and then compute the $\hat{I}_{a}$ integrals. Then, one updates the values of the polynomial parameters and derive the new values of all the functions $M_{I}$, $G_{10}$, $N_{0}$, $H_{10}$, $G_{12}$ and those of discontinuities $\hbox{disc}[M_{I}]$.

Convergence is found to be reasonably fast. Denoting by $M_{I}^{(n)}$ the result obtained at the $n^{th}$ iteration step we can estimate the rate of convergence from the quantity

[TABLE]

Anticipating on the numerical results to be presented below, with $n=5,6,7$ we find: $\epsilon^{(5)}\simeq 4\cdot 10^{-3}$, $\epsilon^{(6)}\simeq 2\cdot 10^{-4}$, $\epsilon^{(7)}\simeq 4\cdot 10^{-5}$. The $I=2$ amplitude is the one which has the slowest convergence rate.

5.2 Input $I=0,1$ T-matrices

Above the $K\bar{K}$ threshold, the $S$ and $T$ matrices are related by

[TABLE]

Two-channel unitarity implies that all the $T$-matrix elements are determined from three real inputs: a) the phase of $S_{11}$, b) the modulus of $T_{12}$ and c) the phase of $T_{12}$. For $I=0$ scattering, experimental data exist for these quantities up to 2 GeV, approximately. We will use here a determination based on the experimental data of Hyams et al. [41] for the $\pi\pi\to\pi\pi$ phase-shifts and the data of refs. [42] and [43] (above 1.3 GeV) for the $\pi\pi\to K\bar{K}$ amplitude555We note, however, that a recent analysis of the $\pi\pi$ Roy equations in a once-subtracted version[44] shows some tension with the data of Cohen et al. [42] assuming that it saturates $\pi\pi$ inelasticity near the $K\bar{K}$ threshold.. In the higher energy region, a smooth interpolation is performed with the following asymptotic conditions,

[TABLE]

which insure (in general) the existence of a corresponding unique MO matrix [45]. Below the $K\bar{K}$ threshold, a determination of the phase-shift based on the data of ref. [41] together with constraints from the $\pi\pi$ Roy equations (in the twice-subtracted version of ref. [46]) is performed. It is assumed that inelasticity can be neglected in this region, which implies that the phase of the $\pi\pi\to K\bar{K}$ amplitude coincides with the $\pi\pi$ phase-shift (Fermi-Watson theorem). This allows one to determine the modulus of this amplitude below the $K\bar{K}$ threshold [47] where it is also needed. Details on the parametrisation can be found in refs. [48, 49]. Fig. 2 shows the experimental data used and the fitted curves.

In the case of the $I=1$ $T$-matrix, experimental information on $\eta\pi$ and $K\bar{K}$ scattering are indirect and far less detailed than for $I=0$. We will rely here on the chiral $K$-matrix model proposed in ref. [50] which provides a simple interpolation between certain known properties of the prominent $I=1$ scalar resonances $a_{0}(980)$ and $a_{0}(1450)$ and the low energy properties constrained by chiral symmetry. The $T$-matrix reproduces, by construction, the $\eta\pi\to\eta\pi$ and $\eta\pi\to K\bar{K}$ chiral amplitudes at NLO (and more approximately the amplitude $K\bar{K}\to K\bar{K}$). It depends on the values of the $O(p^{4})$ chiral couplings: we use here (as in ref. [50]) a set of values for these taken from a $p^{4}$ fit of ref. [51]. We note that the value of $L_{3}$ in this set is $L_{3}=-3.82\cdot 10^{-3}$. We will use this value also in the computation of the $\eta\to 3\pi$ amplitude, for consistency.

A further constraint can be implemented by computing the $\eta\pi$ and $K\bar{K}$ scalar form-factors from this $T$-matrix. Chiral symmetry relates the $\eta\pi$ and the $K\pi$ scalar radii, which leads to the prediction that ${\langle r^{2}\rangle}_{S}^{\eta\pi}$ should be remarkably small, ${\langle r^{2}\rangle}_{S}^{\eta\pi}\simeq 0.1$ $\hbox{fm}^{2}$. This small value can be reproduced provided the phase $\delta_{\eta\pi\to K\bar{K}}$ raises sufficiently slowly above the $K\bar{K}$ threshold. The phenomenological parameters of the model are also constrained by the properties of the resonances. Fig. 3 shows a typical result for the $\eta\pi$ phase shift and for the phase and modulus of the $\eta\pi\to K\bar{K}$ amplitude666The phase shown is that of $T_{0}(\eta\pi^{+}\to K^{+}{\bar{K}^{0}})=-g_{0}^{\eta\pi}$ (according to the isospin convention of eq. (13)). It satisfies $\delta^{0}_{\eta\pi\to K\bar{K}}(m_{a_{0}(1450)})=100\degree$ which corresponds to ${\langle r^{2}\rangle}^{\eta\pi}_{S}=0.12$ $\hbox{fm}^{2}$, within 20% of the chiral $O(p^{4})$ result. which we will use in the present work.

Finally, the phase shifts which we used for the $I=1$ $\pi\pi$ $P$-wave and the $I=2$ $S$-wave, for which inelasticity is ignored, are shown in fig. 4. In the energy region $\sqrt{s}\leq 0.8$ GeV, these phase shifts are given by the Roy equations solution parametrisation of ref. [46]. They are fitted to experimental data from [41] ($P$-wave) and [52, 53] ($I=2$) in the region $\sqrt{s}\leq 1.5$ GeV and interpolated to $\delta_{1}^{1}(\infty)=\pi$, $\delta_{0}^{2}(\infty)=0$ in the higher energy region.

5.3 Illustration of the role of the inelastic channels

Results for the $\eta\to 3\pi$ amplitude obtained from solving numerically the Khuri-Treiman equations are presented in fig. 5 which shows the real part of the amplitude along the line $t=u$ as a function of $s$. Let us consider the role of the inelastic channels in four energy regions

• a)

In the neighbourhood of $s=1$ $\hbox{GeV}^{2}$ there is a very sharp energy variation, as one could have easily expected, induced by the interference of the $a_{0}$ and $f_{0}$ resonances and the presence of the ${K^{+}}{K^{-}}$ and ${K^{0}}{\bar{K}^{0}}$ thresholds.

• b)

In the region 0.7\mathrel{\raise 2.32501pt\hbox{<}\kern-7.00002pt\lower 2.5187pt\hbox{\sim}}s\mathrel{\raise 2.32501pt\hbox{<}\kern-7.00002pt\lower 2.5187pt\hbox{\sim}}0.97 $\hbox{GeV}^{2}$ the coupled-channel amplitude displays a large enhancement as compared to the single-channel amplitude.

• c)

In the lower energy region, s\mathrel{\raise 2.32501pt\hbox{<}\kern-7.00002pt\lower 2.5187pt\hbox{\sim}}0.7 $\hbox{GeV}^{2}$, on the contrary, the effect of the inelastic channels is to reduce the size of the amplitude. One also observes that in this region the influence of the inelastic channels becomes small.

• d)

In the sub-threshold region, finally, the coupled-channel and single-channel amplitudes are essentially indistinguishable. This is expected since the amplitudes are constrained to satisfy the same chiral matching equations.

It is not difficult to identify the main mechanism which generates the behaviour described above. For this purpose, let us consider the $i=1,j=1$ component of the matrix $\bm{M_{0}}$ and let us absorb the the effect of the integrals $\hat{I}_{a}$, $\hat{I}_{a}$ at $s=0$ into the polynomial matrix, defining

[TABLE]

The real parts of the three components of $\tilde{\bm{P}_{0}}$ which are related to the $K\bar{K}$ channel have the following expressions

[TABLE]

These polynomial coefficients are controlled, we recall, by the matching conditions to the LO chiral $\eta\pi\to K\bar{K}$, $\pi\pi\to K\bar{K}$, and $K\bar{K}\to K\bar{K}$ isospin violating amplitudes (see appendix B). The component ${\rm Re\,}[\bm{P}_{0}]_{21}$ is negative in the region $s>4m_{\pi}^{2}$ and dominates the others in the range s\mathrel{\raise 2.32501pt\hbox{>}\kern-7.00002pt\lower 2.5187pt\hbox{\sim}}0.2 $\hbox{GeV}^{2}$. In fact, the corresponding contribution to $\bm{M}_{0}$ dominates in the whole region $4m_{\pi}^{2}<s<0.95$ $\hbox{GeV}^{2}$. In this region, the contribution from the inelastic channels can thus be written approximately as

[TABLE]

The components of the Omnès matrices which appear in eq. (93) are plotted in fig. 6. The salient feature here is that the real part of the $I=0$ component $(\bm{\Omega}_{0})_{12}$ is positive at low energy and changes sign777The presence of a zero below the $K\bar{K}$ threshold follows from Watson’s theorem which is obeyed by the component $(\bm{\Omega}_{0})_{12}$ and leads to the relation ${\rm Re\,}(\bm{\Omega}_{0})_{12}/{\rm Im\,}(\bm{\Omega}_{0})_{12}=\cot\delta_{0}^{0}$. This implies that ${\rm Re\,}(\bm{\Omega}_{0})_{12}$ vanishes when the phase shift $\delta_{0}^{0}$ goes through $\pi/2$. at $s\simeq 0.73$ $\hbox{GeV}^{2}$. This is the main reason why the inelastic channels decrease the $\eta\to 3\pi$ amplitude below $0.7$ $\hbox{GeV}^{2}$ and increase it above. This behaviour is enhanced by the $I=1$ component $(\bm{\Omega}_{1})_{11}$ which is larger than 1 (reflecting the attractive nature of the $\eta\pi\to\eta\pi$ interaction below 1 GeV).

The behaviour of the amplitude along the $t=s$ (or $u=s$) line, which displays the Adler zero, is shown in fig. 7. In the sub-threshold region, the chiral, the single-channel and the coupled-channel amplitudes are seen to be very close. For the position of the Adler zero, we find

[TABLE]

Finally, the results for the $\eta\to 3\pi^{0}$ amplitude are shown in fig. 8. The influence of the inelastic channels are seen to be quite substantial in this case in the whole low-energy region. One also sees that there is no region in which there is close agreement between the dispersive and the chiral $O(p^{4})$ amplitudes. This, of course, is because of the occurrence of the combination $M_{0}(s)+M_{0}(t)+M_{0}(u)$ and the fact that at least one of the variables $s$, $t$, $u$ must lie above the $\pi\pi$ threshold, thus generating significant $S$-wave rescattering chiral corrections.

5.4 Dalitz plot parameters

One traditionally describes the Dalitz plot in terms of two dimensionless variables $X$, $Y$ such that $|X|,|Y|\leq 1$ and the centre of the Dalitz plot corresponds to $X=Y=0$. In the case of the charged decay amplitude $\eta\to\pi^{+}\pi^{-}\pi^{0}$, the variables $X$, $Y$ are related to the Mandelstam ones by

[TABLE]

with $Q_{c}=m_{\eta}-2m_{\pi^{+}}-m_{\pi^{0}}$. Assuming charge conjugation invariance, the amplitude must be invariant under the transformation $X\to-X$ and a polynomial parametrisation of the amplitude squared can be written as

[TABLE]

In the case of the neutral decay amplitude $\eta\to 3{\pi^{0}}$, $Q_{c}$ must be replaced by $Q_{n}=m_{\eta}-3m_{\pi^{0}}$ in the definition of $X$ and $Y$. Charge conjugation and Bose symmetry imply that the amplitude must be invariant under the two transformations

[TABLE]

with $z=X+iY$. The amplitude squared can thus be represented as

[TABLE]

A direct comparison of the dispersive amplitudes squared with the experimental data from KLOE [16] is shown in fig. 9 and our numerical results for the Dalitz plot parameters are collected in table 2. The numbers quoted in the table are obtained in a way which parallels the experimental determination: a discrete binning of the Dalitz plot is performed (with a few hundred bins) and a global least squares fit of the theoretical (chiral or dispersive) amplitudes squared is performed using the representations (96), (98). The table also shows the two most recent experimental determinations [15, 16].

It is clear, at first, that the amplitudes obtained from solving the Khuri-Treiman equations and constrained to match the chiral NLO amplitudes are in much better agreement with the experimental results in the physical decay region than the NLO amplitude itself. In particular, the parameter $b$ which was too large by a factor of three is reduced by a factor of two and the parameter $\alpha$ which was positive becomes negative. This is in close agreement with the results obtained a long time ago by Kambor et al. [24]. Our main new result is that taking into account the $K\bar{K}$ inelastic channels and the effect of $\eta\pi$ rescattering has a non negligible influence on the Dalitz parameters and tends to further improve the agreement with experiment. The influence of these inelastic channels for the parameters $d$ and $g$ is small (less than 5%) but quite significant for the parameter $b$ which is reduced by 17% and now lies within 15% of the experimental value. This reflects the reduction of the amplitude caused by the inelastic effects at low energy discussed in sec. 5.3. Similarly, the parameter $\alpha$ is modified by approximately 20% by the $K\bar{K}$ inelastic channels and becomes rather close to the experimental result. The parameter $g$ is the only one which shows a mismatch, by a factor of two, with the value measured by KLOE.

Finally, let us consider the sensitivity of the Dalitz plot parameters to the strength of the $\eta\pi$ interaction, which is not precisely known at present. This is illustrated in table 3. The table also shows that varying the $O(p^{4})$ coupling $L_{3}$ has a significant influence, in particular on the parameter $d$.

5.5 The ratio $\Gamma(\eta\to 3\pi^{0})/\Gamma(\eta\to{\pi^{+}}{\pi^{-}}{\pi^{0}})$

Let us quote here the results for the ratio of the $3{\pi^{0}}$ and the ${\pi^{+}}{\pi^{-}}{\pi^{0}}$ decay rates

[TABLE]

We find that the influence of the inelastic channels on this ratio is very small,

[TABLE]

As compared with the chiral $O(p^{4})$ result

[TABLE]

this ratio is thus predicted to increase under the effect of the final-state interactions, by only 2%. The experimental status of this quantity is not completely clear at present, as the PDG[54] quotes two different numbers

[TABLE]

Besides, the CLEO collaboration [56] has performed an experiment dedicated to the determinations of the $\eta$ meson decay branching fractions and they quote

[TABLE]

as the most precise determination of the $3{\pi^{0}}/{\pi^{+}}{\pi^{-}}{\pi^{0}}$ ratio.

5.6 The quark mass ratio $Q$ from the chirally matched dispersive

amplitude

It must first be reminded that, once the electromagnetic interaction is taken into account, the quark mass ratio $Q$ is no longer invariant under the QCD renormalisation group since the quark mass variation with the scale depends on its electric charge

[TABLE]

(with $e_{u}=2/3$, $e_{d}=e_{s}=-1/3$). This implies the following scale variation for the factor $\epsilon_{L}$

[TABLE]

It can then easily be verified, using the equations of appendix A which include the $e^{2}m_{\pi}^{2}$ contributions [6], that the scale invariance of the complete NLO chiral amplitude is restored thanks to the combination of two of the electromagnetic coupling constants [4], $K_{9}^{r}+K_{10}^{r}$. Indeed, as shown in refs. [57, 58] this combination depends not only on the chiral scale $\mu$ but also on the QCD scale $\mu_{0}$ and satisfies

[TABLE]

In practice, this means that in order to determine $Q(\mu_{0})$ from the $\eta\to 3\pi$ amplitude we must specify the values of the electromagnetic chiral couplings $K_{i}$ at the corresponding scale. We choose here $\mu_{0}=0.77$ GeV and estimate the values of the couplings $K_{i}^{r}(\mu,\mu_{0})$ from a resonance saturation model [59]. Such estimates are qualitative at best but it can be shown, based on general order of magnitude arguments, that the uncertainty induced on the amplitude should not exceed a few percent [6, 7].

Having verified that the dispersive amplitude is in qualitative agreement with experiment concerning the Dalitz plot parameters, we can make an estimate of the value of $Q$. In the present approach, the amplitude in the physical decay region is derived as a Khuri-Treiman solution uniquely defined from the chiral NLO amplitude by the set of four matching equations, and thus has a definite dependence on $Q$. We can then estimate the value of $Q$ by the requirement that the dispersive amplitude reproduces the experimental values of the $\eta\to 3\pi$ decay widths

[TABLE]

taken from the PDG [54] (constrained fit). Doing this, we find

[TABLE]

where the quoted errors only reflect the experimental ones on the widths (we comment below on the theoretical error). This shows that the effect of the inelastic channels on the determination of $Q$ is rather small, of the order of 1%, and tends to decrease its value.

The central value of $Q$ is somewhat smaller than the results which are obtained from lattice QCD+QED simulations of hadron masses: $Q=22.9\pm 0.4$ (ref. [60]), $Q=23.4(0.4)(0.3)(0.4)$ (ref. [61]). The error on $Q$ associated with the phase-shifts below 1 GeV is rather small, of the order of 1%. The error associated with the NLO amplitude, essentially related to $L_{3}$, is of the same order. The largest error arises from chiral NNLO contributions to the amplitude which will modify the determination of the polynomial parameters via the matching conditions. Assuming that they induce a 10% relative error in the determination of each one of the four polynomial parameters $\alpha_{0}$, $\beta_{0}$, $\gamma_{0}$, $\beta_{1}$ and assuming the errors to be independent, gives the following theory error on $Q$

[TABLE]

$Q$ being mostly sensitive to the variation of the first two parameters $\alpha_{0}$, $\beta_{0}$. We have also varied the remaining 14 polynomial parameters, assuming a 20% relative error, and found that they have a much smaller influence. Within the error (109), the determination based on $\eta\to 3\pi$ decay is compatible with the lattice QCD results, which confirms that the size of the NNLO corrections to the $\eta$ decay amplitude in the sub-threshold region should not exceed $10\%$.

5.7 Further experimental constraints on $Q$

Our estimate for the error on $Q$ (eq. (109)) was based on a general order of magnitude assumption on the size of the NNLO corrections to the four leading polynomial parameters. More precise information on the size of these corrections can be derived by making use of the precise experimental results on the energy dependence across the Dalitz region, imposing that the dispersive amplitude reproduces these via a least-squares fit. Not all the four leading polynomial parameters can get independently constrained in this way since a ratio of amplitudes is involved. We will make the simple choice to fix one of them, $\alpha_{0}$, from its chiral value and to perform a variation of the three others $\beta_{0}$, $\gamma_{0}$, $\beta_{1}$. We will use the latest KLOE experimental data [16] which consist of a set of amplitudes squared, $|{\cal T}_{exp}(X_{i},Y_{i})|^{2}$, measured over 371 energy bins in the Dalitz region and satisfying the normalisation condition

[TABLE]

We introduce corresponding theory amplitudes ${\cal T}_{th}(X,Y)=M_{c}(X,Y)/M_{c}(0,0.05)$ and define the $\chi^{2}$ as

[TABLE]

allowing for a floating of the normalisation within the experimental error via a parameter $\lambda$.

At first, setting $\lambda=1$ and computing the $\chi^{2}$ with our chirally matched central amplitude with $L_{3}=-3.82\cdot 10^{-3}$ we obtain $\chi^{2}=3079$. If, instead, we use the value recently derived in ref. [55] from $K_{l4}$ decays: $L_{3}=-2.65\cdot 10^{-3}$, one obtains a significantly reduced result: $\chi^{2}=714$. It thus seems reasonable to use this value as a starting point, which fixes the central value of $\alpha_{0}$. We have then searched for a minimum of the $\chi^{2}$ by varying the real parts of the three polynomial parameters $\beta_{0}$, $\gamma_{0}$, $\beta_{1}$ (keeping their imaginary parts fixed) and the normalisation $\lambda$. The other polynomial parameters, in particular $\alpha_{0}$, are kept fixed to their chirally-matched value. In this way, we obtain at the minimum: $\chi^{2}_{min}=387$, which corresponds to a value per degree of freedom $\chi^{2}_{min}/dof=1.055$. The fitted values of the polynomial parameters differ from the chirally-matched values by less than 5% (the numerical values are given in appendix D). This confirms that this difference can be consistently attributed to chiral NNLO effects.

The theory error in this approach is dominated by the variation (by 10%) of the single polynomial parameter $\alpha_{0}$ (we have also added an error associated with the input $T$-matrices and an error from varying the imaginary parts of the polynomial parameters). The error induced by the variation of the parameters $\beta_{0}$, $\gamma_{0}$, $\beta_{1}$ is now computed using the covariance matrix as evaluated by the MINUIT fitting program [62] (this matrix is given in appendix D). This error is added in quadrature with the experimental error on the $\pi^{+}\pi^{-}\pi^{0}$ and $3\pi^{0}$ decay rates. Finally, the result for the quark mass ratio as determined from this fitted amplitude can be written as

[TABLE]

We find it useful to quote the theoretical and experimental errors separately since the former is not necessarily gaussian. Previous determinations of $Q$ which combine chiral constraints and fits to Dalitz plot data have been performed in refs. [63, 28, 64, 65, 66, 67]. In the most sophisticated of these approaches [67], five polynomial parameters are included in the fit and the effect of the ${\pi^{+}}-{\pi^{0}}$ mass difference in the amplitudes is accounted for.

Finally, we quote the values of the $3\pi^{0}$ Dalitz plot parameters which can be predicted from our fitted amplitude

[TABLE]

6 Conclusions

We have proposed an extension of the Khuri-Treiman formalism for the $\eta\to 3\pi$ amplitude which includes the $\eta\pi$ and the $K\bar{K}$ channels in the unitarity equations in addition to the elastic $\pi\pi$ channel. Modulo some approximations (in particular we do not attempt to impose unitarity in the crossed channels involving kaons like $\pi K\to\pi K$ or $\eta K\to\pi K$) the equations for the one-variable functions $M_{0}$ and $M_{2}$ are shown to be simply replaced by $2\times 2$ matrix equations. These are given in eqs. (69) and (77) which involve both the $I=0$ and the $I=1$ Omnès $2\times 2$ matrices. Eq. (69) exhibits explicitly the contribution induced by the physical $K^{0}-K^{+}$ mass difference via unitarity in integral form.

The amplitudes derived from this extended framework should be valid in an energy range which covers the physical decay region and also the physical region of the scattering $\eta\pi\to\pi\pi$ below 1 GeV. Given a fixed number of polynomial parameters, an improved precision at low energy should result from the fact that the effects of the two prominent light scalar resonances $a_{0}(980)$ and $f_{0}(980)$ are taken into account in the dispersive integrals.

Using four polynomial parameters in the $\eta\to 3\pi$ amplitude we have reconsidered the idea of performing a prediction of the amplitude in the physical region as an extrapolation of the $O(p^{4})$ chiral amplitude, uniquely defined by fixing all the polynomial parameters by matching conditions. These are imposed in the form of a set of equations which insure that the differences between the dispersive and the $O(p^{4})$ chiral $\eta\to 3\pi$ amplitude are of order $p^{6}$ or higher. One verifies then, that the chiral and the dispersive $\eta\to\pi^{+}\pi^{-}\pi^{0}$ amplitudes are very close in the neighbourhood of the Adler zero. These conditions also insure, for the charged decay amplitude, that the single and multi-channel dispersive amplitudes are quasi-identical in the whole region $0\leq s\leq 4m_{\pi}^{2}$, $|t-u|\leq(m_{\eta}+m_{\pi})^{2}$. In contrast, for $\eta\to 3\pi^{0}$, one finds that the unitarity induced chiral corrections are significant even in the sub-threshold region.

We have considered the Dalitz plot parameters and we found that the induced influence of the $a_{0}(980)$, $f_{0}(980)$ resonances is not negligible, in particular for the neutral mode. The modifications of the parameters, in the coupled-channel framework, go in the sense of improving the agreement with experiment, in particular for the parameters $a$, $b$, $f$ of the charged mode. The parameter $\alpha$, for the neutral mode is modified by 20% by the effects of the resonances and lies rather close to the experimental value. The remaining differences between the experimental and the dispersive-theoretical amplitude suggest that NNLO contributions are needed in the matching conditions, at the 5-10% level, which seems quite plausible. Some of these NNLO effects could be accounted for in a more general framework which would implement both unitarity and crossing symmetry completely for all the channels involved. This is left for future work.

The $\eta\to 3\pi$ amplitudes constructed in the present approach inherit a well defined dependence on the quark mass ratio $Q$ from that of the chiral NLO amplitude. We can then determine $Q$ such as to reproduce the integrated decay widths. The central value that one obtains is somewhat low compared to the recent determinations from lattice QCD simulations but it is compatible within the uncertainty induced by the NNLO effects in the matching. Some knowledge of these NNLO effects seems necessary in order to improve the precision of the determination of $Q$ by this method.

Acknowledgements

B. M. acknowledges rewarding discussions with Émilie Passemar. M. A. acknowledges financial support from the Spanish MINECO and European FEDER funds under the contracts 27-13-463B-731 (“Juan de la Cierva” program), FIS2014-51948-C2-1-P, FIS2014-57026-REDT, and SEV-2014-0398, and from Generalitat Valenciana under contract PROMETEOII/2014/0068.

Appendix A The amplitude $\eta\to{\pi^{+}}{\pi^{-}}{\pi^{0}}$ at chiral order $p^{4}$ and

$e^{2}p^{2}$

We give below the explicit expressions for the three functions $\bar{M}_{I}$ in the chiral expansion at order $p^{4}$ and including $e^{2}p^{2}$ contributions as used in the present work. They are given in a form which is manifestly scale invariant independently of the values of $m_{\pi}$, $m_{K}$, $m_{\eta}$. The amplitude is identical with the one originally computed in ref. [17] when $e^{2}=0$ and the $m_{\pi}$, $m_{K}$, $m_{\eta}$ values in the NLO part obey the Gell-Mann-Okubo relation $3m_{\eta}^{2}=4m_{K}^{2}-m_{\pi}^{2}$. Following ref. [17] the LEC’s $L^{r}_{5}$, $L^{r}_{7}$, $L^{r}_{8}$ are expressed in terms of the two physical quantities888In eq. (114) $\bar{m}_{P}$ are QCD masses for which we use the values provided by the FLAG review [68] $\bar{m}_{\pi}=0.1348$, $\bar{m}_{K}=0.4942$ and $\bar{m}_{\eta}\simeq m_{\eta}=0.537862$ (all in GeV) which gives $\Delta_{GMO}=0.2068$. Also using $F_{K}/F_{\pi}$ from this review gives $\Delta_{F}=0.2005$. Elsewhere in the chiral formulae we use $F_{\pi}=92.21$ MeV (from the PDG), $m_{\pi}=0.13957$, $m_{K}=m_{K^{+}}=0.493677$ GeV and $L_{3}=-3.82\cdot 10^{-3}$.

[TABLE]

and in terms of the quark mass ratio $(m_{s}+\hat{m})/\hat{m}$ (see ref. [3]). We also include the electromagnetic contributions of order $e^{2}p^{2}$ evaluated in ref. [6] which allows one to identify the piece induced by the physical $K^{+}-K^{0}$ mass difference via unitarity. Further electromagnetic corrections which have been computed in ref. [7] are not included here. The expressions for $\bar{M}_{I}$ given below also implement the $w=0$ conditions (10), which simplifies somewhat the writing of the matching relations.

Using the following notation

[TABLE]

the function $\bar{M}_{0}$ reads:

[TABLE]

The contributions $\Delta\bar{M}_{0}^{a}$, $\Delta\bar{M}_{0}^{b}$ are induced by the electromagnetic interaction and proportional to $e^{2}$,

[TABLE]

The last contribution, $\Delta\bar{M}_{0}^{c}$, is induced by the physical $K^{+}-K^{0}$ mass difference

[TABLE]

The parameters $C$ and $K^{r}_{i}$ in the above expressions are the chiral coupling constants which appear at order $e^{2}$ and $e^{2}\,p^{2}$ respectively [4].

The chiral expression for the function $\bar{M}_{1}$ reads

[TABLE]

and it has no electromagnetic contributions. The function $\bar{M}_{2}$, finally, reads

[TABLE]

with the electromagnetic contribution

[TABLE]

The coupling $C$ can be simply determined from the $\pi^{+}-\pi^{0}$ mass difference,

[TABLE]

The couplings $K_{i}^{r}$ are expected to have an order of magnitude $K_{i}^{r}(\mu=m_{\rho})\sim 1/(16\pi^{2})$ but, otherwise, they are not precisely known. Fortunately, in the $\eta\to 3\pi$ amplitude, they appear multiplied by $e^{2}m_{\pi}^{2}$. This is in contrast to other isospin violating observables like $m^{2}_{K^{+}}-m^{2}_{K^{0}}$ in which they appear multiplied by the larger factor $e^{2}m_{K}^{2}$. A simple resonance saturation estimate [58, 59] gives: $-2K_{3}+K_{4}\simeq-4.0\cdot 10^{-3}$, $K_{5}+K_{6}\simeq 14.4\cdot 10^{-3}$, $K_{9}+K_{10}\simeq 7.5\cdot 10^{-3}$ (with $\mu=\mu_{0}=0.77$ GeV) which suggests that there are cancellations in the combination relevant for $\eta\to 3\pi$

[TABLE]

Appendix B Matching equations

We reproduce below the set of matching relations from which we determine the set of 16 polynomial parameters of the Khuri-Treiman coupled-channel equations (i.e. eqs. (69), (77) and the second one of eqs. (2.4)). In order to simplify the relations it is assumed here that the chiral expressions for the $\eta\to 3\pi$ functions $\bar{M}_{I}$ satisfy, as in appendix A, the $w=0$ relations $\bar{M}_{1}(0)=\bar{M}_{2}(0)=\bar{M}^{\prime}_{2}(0)=0$. Derivatives at $w=0$ are denoted either by dots or by primes and matrix elements of the $I=0,1$ MO matrices are denoted here by $\Omega^{(I)}_{ij}$. The chiral LO expressions for the $K\bar{K}$ amplitudes $\bar{N}_{0}$, $\bar{G}_{10}$, $\bar{H}_{10}$ and $\bar{G}_{12}$ are given in eqs. (4.3). A first set of five relations is

[TABLE]

Then, introducing the notation

[TABLE]

we have

[TABLE]

and

[TABLE]

where $\hat{I}_{ij}$ denote the matrix elements of the matrix sum $\hat{\bm{I}}_{a}+\hat{\bm{I}}_{b}$ (see eqs. (65) (68) (71)). The final four relations read

[TABLE]

Appendix C Angular integrals and hat functions

Using the dispersive representations (2.1) of the functions $M_{I}$, one can express the angular integrals in the following form which displays explicitly their singularity when $s\to(m_{\eta}-m_{\pi})^{2}$. For $I=0,2$ one has

[TABLE]

where $\kappa(s)$ is given in eq. (23). For $I=1$ one has

[TABLE]

The functions $R_{I}^{n}(s)$ which control the singularities arise from the part of the $t$ integration contour which encircle the unitarity cut, they are given by

[TABLE]

(where $t^{\pm}(s)$ is given in eq. (87)) in the $s$ range

[TABLE]

and $R_{I}^{n}(s)=0$ otherwise. In particular, no divergence occurs when $s\to 4m_{\pi}^{2}$ or $s\to(m_{\eta}+m_{\pi})^{2}$.

The kernels which are needed here are given by the following expressions

[TABLE]

and

[TABLE]

The function $L(t^{\prime},s)$ arises from the parts of the $t$ integration contour not taken into account in the functions $R_{I}^{n}$ (the function $L(t^{\prime},s)$ thus vanishes when $s\to(m_{\eta}-m_{\pi})^{2}$, such that the kernels remain finite) it is given by

$4m_{\pi}^{2}\leq s<(m_{\eta}-m_{\pi})^{2}$:

[TABLE]

• 2)

$(m_{\eta}-m_{\pi})^{2}\leq s<m_{\eta}^{2}-5m_{\pi}^{2}$:

[TABLE]

• 3)

$m_{\eta}^{2}-5m_{\pi}^{2}\leq s<\infty$

[TABLE]

Appendix D Supplementary material

The integral equations for the amplitudes $M_{I}$ depend linearly on the 16 polynomial parameters. The matching equations are also linear. One can then express the amplitudes in which, for instance, the four leading parameters $\alpha_{0}$, $\beta_{0}$, $\gamma_{0}$, $\beta_{1}$ are fixed to some arbitrary values and the remaining 12 parameters are fixed from the corresponding matching equations in the form of a linear superposition

[TABLE]

We provide our numerical results for the amplitudes in which the parameters $\alpha_{0}$, $\beta_{0}$, $\gamma_{0}$, $\beta_{1}$ are either 0 or 1 in five data files: MI_0000.dat, MI_1000.dat, MI_0100.dat, MI_0010.dat, MI_0001.dat. In each file, the first column is $s$ (in $\hbox{GeV}^{2}$) and the other columns are the corresponding real and imaginary parts of $M_{0}$, $M_{1}$ and $M_{2}$. We note that these amplitudes depend on the energy value $E_{asy}$ above which the $T$ matrix parameters are set to their asymptotic values. We take here $E_{asy}=10$ GeV. We give below the corresponding values of the four polynomial parameters for several cases considered in this paper.

Chirally matched amplitude, $L_{3}=-3.82\cdot 10^{-3}$:

[TABLE]

• 2)

Chirally matched amplitude, $L_{3}=-2.65\cdot 10^{-3}$:

[TABLE]

• 3)

Fitted amplitude:

[TABLE]

Note that $\alpha_{0}$ and the imaginary parts of $\beta_{0}$, $\gamma_{0}$, $\beta_{1}$ are the same as in 2). The covariance matrix of the three fitted parameters reads

[TABLE]

and, finally, the value of the normalisation parameter $\lambda$ (see eq. (111)) is: $\lambda=1.009\pm 0.001$.