Quark-mass dependence in $\omega\to3\pi$ decays

We study the quark-mass dependence of $\omega\to3\pi$ decays, based on a dispersion-theoretical framework. We rely on the quark-mass-dependent scattering phase shift for the pion-pion $P$-wave extracted from unitarized chiral perturbation theory. The dispersive representation then takes into account the final-state rescattering among all three pions. The described formalism may be used as an extrapolation tool for lattice QCD calculations of three-pion decays, for which $\omega\to3\pi$ can serve as a paradigm case.


Introduction
Despite tremendous progress in simulating Quantum Chromodynamics (QCD) on space-time lattices using physical quark masses, many studies of complicated observables within lattice QCD are still performed with light quarks that are heavier than they are in the real world (see e.g. Refs. [1,2] for reviews). To extrapolate such simulations to the physical point, additional theoretical input is required, which should ideally be based on systematically improvable effective field theories. At low energies, the effective field theory that controls the quark-mass dependence by construction is chiral perturbation theory (ChPT) [3][4][5], which describes the interactions of the pseudo-Goldstone bosons of spontaneous chiral symmetry breaking, the pions (as well as kaons and the η).
However, the vast majority of states in QCD are resonances, and to perform chiral extrapolations for these is less straightforward. A popular tool in this regard has been to employ unitarized versions of ChPT, such as the inverse amplitude method (IAM) [6][7][8][9]: a resummation of higher-order effects obeying S-matrix unitarity (i.e., probability conservation) allows one to generate poles a e-mail: isken@hiskp.uni-bonn.de on unphysical Riemann sheets in the complex-energy plane, the signatures of resonances. The IAM can be justified using dispersion theory; scattering amplitudes constructed via the IAM match smoothly on the ChPT expansion at low energies. In this manner, the properties of elastic resonances such as the f 0 (500) and the ρ(770) in pion-pion as well as the K * 0 (700) and K * (892) in pion-kaon scattering have been investigated with respect to their quark-mass dependence [10,11].
Nevertheless, by far not all hadronic resonances appear in two-body scattering processes. The lightest resonance that decays only into a three-body final state (in QCD in the isospin limit) is the ω(782), with its dominant decay ω → 3π. Clearly, the quark-mass dependence of the ω cannot be assessed within an IAMtype formalism; it could at best be studied within the appropriate partial wave of the 3π → 3π scattering process, and the formalism to study such processes on the lattice is currently under intense investigation [12][13][14][15][16][17][18]. In this article, we suggest an approach to assess the quark-mass dependence of the ω → 3π decay amplitude based on dispersion relations. We employ the so-called Khuri-Treiman equations [19] that require the pionpion two-body phase shift as input, which we extract from the known quark-mass-dependent IAM partial wave. While we still need to rely on effective field theory ideas to describe the variation of the ω mass with the quark masses, the dispersive framework allows us to predict its quark-mass-dependent width. The idea to employ dispersion theory to extend the applicability of IAM-generated phase shifts is not new: it has already been applied to describe the pion vector form factor [20], as well as, in a formalism closely related to what we present here, to the reaction γπ → ππ [21].
The outline of the present article is as follows. We recall the description of pion-pion scattering with the IAM formalism in Sect. 2. The Khuri-Treiman formal-ism for ω → 3π is described in Sect. 3. Supplementary assumptions to describe the quark-mass dependence of the ω width are collected in Sect. 4, before we show results in Sect. 5. We summarize our findings in Sect. 6.
2 Pion-pion scattering and the ρ resonance in one-loop unitarized ChPT Before starting with the discussion of the quark-mass dependence of the ω as a three-pion resonance, we briefly summarize the investigation of the quark-mass dependence in ππ → ππ scattering. The results of this section will be an essential input for the study of ω → 3π. We follow the formalism introduced in Refs. [10,22]. For the later purpose our investigation will focus on the ρ(770) resonance, which appears as a pole in the P -wave scattering amplitude. We will treat the ππ rescattering as elastic in this whole section.
The partial-wave decomposition of the ππ → ππ scattering amplitude T I of definite isospin I is defined by with the partial-wave amplitude t I ℓ of angular momentum ℓ, the Legendre polynomial P ℓ , and the s-channel scattering angle z s = cos θ s . Below any inelastic thresholds the partial-wave amplitude is given in terms of the scattering phase shift δ I ℓ only, where σ(s) = 1 − 4M 2 π /s. Since we are interested in the P -wave only (I = ℓ = 1) we will drop the labels for simplicity from now on.
In ChPT the pion mass is given in terms of the light quark masses as an expansion M 2 π = 2Bm + O(m 2 q ), where the leading term is known as the Gell-Mann-Oakes-Renner (GMOR) relation [23], with the two light quark masses combined inm = 1 2 (m u + m d ). The constant B is related to the scalar quark condensate in the chiral limit, which measures the strength of spontaneous symmetry breakdown in QCD. We will work in the isospin limit, meaning that m u = m d and M π ± = M π 0 ≡ M π . A brief discussion of the isospinbreaking effects in ω → 3π can be found in Sect. 4.3. The GMOR relation implies that studying the quarkmass dependence is equivalent to an investigation of the pion-mass dependence. Hence from now on we will refer to the pion-mass dependence instead. Since we are interested in pion-mass-dependent quantities, it turns out to be useful to define parameters at the physical point as e.g.M π ≡ M phys π .
In the ChPT power counting the P -wave amplitude for ππ → ππ scattering up to next-to-leading order can be written as where t i denotes the contribution of chiral order p i . The P -wave projection of the scattering amplitude T I given in Ref. [4] yields where we made use of the abbreviation The value for pion decay constant in the chiral limit F is taken from the ratio F π /F = 1.064(7) [1,[24][25][26][27][28], where F π = 92.28 (9) MeV [29] is the pion decay constant at the physical point. For our purpose it is beneficial to work with F instead of F π , since F is independent of M π . We treat the combination of low-energy constants (LECs)l =l 2 −l 1 , which occurs in the ChPT expression at next-to-leading order, as a free parameter that will be fixed in the following. Note thatl is also independent of M π , since the individual mass dependences ofl 1 and l 2 cancel [4]. This amplitude however cannot capture the effects of the ρ resonance, which we expect to be the dominant effect in the P -wave above the threshold region, since unitarity is only fulfilled perturbatively (Im t 4 = σ|t 2 | 2 ). Furthermore t 2 and t 4 are polynomials in s (up to cuts encoded in the σ dependence), thus the analytic structure of the standard ChPT expression t ChPT does not allow for any poles on the second Riemann sheet.
In order to include the ρ resonance into our amplitude we will use the IAM. This method allows us to construct an amplitude that fulfills unitarity exactly. Up to next-to-leading order the IAM yields which is equivalent to t ChPT up to corrections of O(p 6 ). Note that crossing symmetry is now only fulfilled perturbatively.

Pole position and residue
The characteristic properties of the ρ resonance are encoded in the pole position and residue of the amplitude on the second Riemann sheet. By analytic continuation the amplitude on the second sheet can be expressed in terms of the amplitude on the first sheet [30]. We employ the specific representation [31] t II (s) = t I (s) where t I and t II denote the amplitudes on the first and second sheets, respectively. Thus the pole position of the amplitude is determined by where s pole corresponds to This allows us to identify the mass M ρ , as well as the decay width Γ ρ of the ρ resonance (that is assumed to be a purely elastic resonance with a single decay channel ππ). The location of the ρ pole at the physical pion mass stemming from two studies of ππ scattering with Roy-type dispersion relations [32,33] is used to constrain the up to now undetermined LEC. This is done in the following way: we minimize the distance of the ρ pole position at the physical point of the IAM amplitude with respect to the most precise extraction of the pole position from the GKPY analysis 1 of Ref. [33]. This yieldsl = 5.73 (8).
The pole position of the IAM at the physical point is then given by √ s pole = 0.7620(15) GeV − i0.0778(11) GeV , which is in good agreement with Refs. [10,22] and the real part of the pole position coincides with Refs. [32,33] within the error bars. Nevertheless we observe a tension in the imaginary part, which is ∼ 4 MeV (∼ 2 − 3 standard deviations) larger than the imaginary parts from the Roy analyses. Thus here we reach the limits of the one-loop IAM description with only one free parameter; for more elaborate studies with an O(p 6 ) IAM amplitude containing several LECs, see Refs. [22,34]. However, for the purposes of our study of ω → 3π, we consider the one-loop IAM a sufficiently reasonable description of ππ scattering. Since the amplitude is not limited to the physical value of the pion mass, we are able to calculate M ρ and Γ ρ as a function of M π . The trajectory of the pole position on the second Riemann sheet is displayed in Fig. 1. As expected from its quark content, the mass of the ρ increases if the pion becomes heavier. This behavior can be described to good approximation as a linear function in M 2 π given by where M ρ (0) and a can be matched to the pion-massdependent pole trajectory extracted from Eq. (8), which yields M ρ (0) = 0.7480(16) GeV , a = 0.719(9) GeV −1 . (12) Similar observations have been made by investigating chiral symmetry constraints [35]. The available phase space for the decay decreases with growing pion mass, since the ρ mass increases much more slowly than the pion mass [10]. Thus the width of the ρ becomes smaller for larger values of M π . The coupling g ρππ of the ρ to the ππ system is defined via the residue where the normalization factors are chosen such that it coincides with the naive expression as obtained from a Lagrangian-based narrow-width approximation or a vector-meson-dominance (VMD) model. Equation (13) yields a numerical value of |g ρππ | = 6.12 (4), which is in fair agreement with other determinations [22,33]. Note that the coupling g ρππ extracted from the IAM is pion-mass independent to very good approximation. Thus the pion-mass dependence of Γ ρ is driven by the phase space factor only, see Eq. (14), as confirmed by lattice QCD calculations [36].
As defined in Eq. (2) the scattering phase shift of the IAM amplitude can be extracted via δ(s) ≡ arg t(s).
The results for different values of the pion mass are shown in Fig. 2. As expected, the slope of the phase shift becomes steeper for heavier pions, while the whole curve moves to the right (decreasing width and increasing mass of the ρ). This behavior is also observed by various lattice QCD calculations carried out at different pion masses [36][37][38][39][40][41][42][43][44][45][46][47]. At the physical pion mass the phase shift is in perfect agreement with the Roy analyses of Refs. [48,49] in the low-energy regime up to √ s ∼ 0.8 GeV. Above this energy the Roy solutions are typically continued to an asymptotic value of π. The IAM amplitude on the other hand behaves like lim s→∞ t IAM (s) = − 3π 6l + 2 + 3πi (15) in the high-energy limit, which depends only on the value of the LEC. Thus the phase of the IAM amplitude will not reach π for all reasonable values ofl. Nevertheless the phase shift gives a reasonable parametrization up to the ρ resonance region and thus will be used as a key ingredient for the dispersive representation of the ω → 3π amplitude, as described in Sect. 3. In order to test the effect of the discrepancy of the IAM compared to the Roy solution phase shifts, we will also use the parametrization to account for the correct asymptotic behavior. The parameters δ a and δ b are fixed by ensuring continuity of the phase shift and its derivative at some highenergy scale Λ 2 ∼ 1.2 GeV 2 , at which we switch from the IAM to the asymptotic phase shift. Note that this high-energy continuation has no relevant influence on the ω → 3π amplitude, since its effects can be fully absorbed by adjusting the unknown normalization of the amplitude, see Eq. (31) below. Similarly, any hypothetical, sub-leading, pion-mass dependence in Λ 2 would be far too small to be of relevance.

Formalism and dispersive representation of the ω → 3π amplitude
The transition amplitude for the ω → 3π decay is defined as In our convention the Mandelstam variables are defined according to which fulfill the relation In the s-channel center-of-mass system, t and u can be expressed as where z s = cos θ s is the scattering angle with the Källén function λ(x, y, z) = x 2 + y 2 + z 2 − 2(xy + xz + yz). Similar expressions hold for the t and u-channels, respectively. The physical thresholds in the three channels are Since the transition ω → 3π is of odd intrinsic parity, the amplitude can be further decomposed into a kinematic prefactor and a scalar function F (s, t, u) containing the dynamical information, Here ǫ µ (P ) denotes the polarization vector of the ω meson. The modulus of the amplitude is given by The expression in the square brackets is also known as the Kibble cubic [50]. Due to Bose symmetry only odd partial waves are allowed to contribute to the process. Thus the partialwave decomposition for the scalar function F (s, t, u) in the s-channel reads where P ′ ℓ (z s ) denotes the differentiated Legendre polynomials. A particular partial wave can be projected out by making use of As the available phase space in the ω → 3π decay is rather small, the dominant contribution will come from the ℓ = 1 partial wave (see Ref. [51] for a discussion of potential F -wave contributions). Neglecting discontinuities from F -and higher partial waves allows us to decompose the scalar function into a sum of singlevariable functions [51][52][53][54][55][56] where F (s) possesses only a right-hand cut. This kind of decomposition is known as a reconstruction theorem [57][58][59][60]. The symmetry in this decomposition reflects the fact that the process is invariant under the exchange of the pions. Combining Eqs. (26) and (27) leads to The right-hand cut of the partial wave f 1 (s) is contained in F (s), while its left-hand cut contributions reside entirely in the projection of the crossed-channel single-variable functionsF (s).
The dispersive framework to describe this decay was already used in previous studies of ω → 3π [51,53] as well as the closely related processes γ ( * ) → 3π [52,[54][55][56]61]. It uses the formalism of the so-called Khuri-Treiman equations [19,62], which is based on analytic continuation of a crossed scattering amplitude in the decay mass [63]. Final-state pion-pion rescattering is assumed to be elastic, and can be described in terms of the phase shift only, see Eq. (2). The unitarity relation for the partial wave f 1 (s) is given by 2 where δ(s) is the ππ P -wave phase shift. Inserting Eq. (28) and noting that disc f 1 (s) = disc F (s) along the right-hand cut, we conclude where the Omnès function Ω(s) is given by [66] The basis function F α (s) can be constructed independently of the numerical value of the subtraction constant α, which therefore can be determined a posteriori. As α serves as an overall normalization of the amplitude, at physical pion masses it is fixed to the total rate Γ (ω → 3π), with the energy dependence of the amplitude or the Dalitz plot distribution then being a theoretical prediction [51].

Decay width and M π dependencies
The process ω → 3π gives the by far dominant contribution to the total ω decay width, B(ω → 3π) = (89.2± 0.7)%. Besides that, the main subleading contributions stem from B(ω → πγ) = (8.4 ± 0.2)% (electromagnetic) and B(ω → ππ) = (1.5 ± 0.1)% (isospin-breaking). Together these contributions account for > 99% of the decay width at the physical point [29]. As long as we restrict ourselves to the isospin limit and strong contributions only, the decay width is fully driven by ω → 3π. Thus the decay width Γ (ω → 3π) ≡ Γ ω is obtained by integrating the squared amplitude over phase space according to This expression has several M π dependencies besides the explicit ones (integration boundaries and Mandelstam variables), which will be discussed in the following.

Pion-mass dependence of the ω mass
In contrast to the case of the ρ, for which we can derive the pion-mass dependence of the complete pole position in the ππ → ππ P -wave amplitude by means of the IAM, we are not in the position to do the same for the ω within some 3π → 3π amplitude of the appropriate quantum numbers. We will discuss the complicated pion-mass dependence of the width of the ω, or the imaginary part of its pole in the complex plane, which is the main focus of this study, in the following; for the pion-mass dependence of its mass, the corresponding real part, we have to resort to symmetry arguments based on effective Lagrangians. These will relate M ω (M 2 π ) to M ρ (M 2 π ), which we have discussed in Sect. 2.1.
We briefly recapitulate the analysis of the leading symmetry-breaking effects in the masses of the vector meson nonet [67][68][69]. Here, the vector mesons are treated as static matter fields; the effective Lagrangian is organized in terms of increasing chiral dimension as well as using the expansion in the inverse number of colors 1/N c . We neglect isospin breaking and electromagnetic effects [68], and ignore deviations from ideal mixing. In this approximation, the symmetry-breaking part of the effective Lagrangian can be written as where contains the (nonrelativistic) vector-meson fields, and χ = diag(M 2 π , M 2 π , 2M 2 K − M 2 π ) breaks SU (3) flavor symmetry due to the different quark masses. Among the terms in Eq. (34), the quark-mass-independent operator ∝ δ is 1/N c suppressed and breaks nonet symmetry; the term ∝ a is chirally suppressed, but the dominant flavor-breaking term in the large-N c limit; and the operators ∝ b and ∝ c are both chirally and 1/N c suppressed. The term ∝ c leads to a common shift in all nonet masses and hence cannot be discerned from the common mass M V using experimental data only; on account of the fact that we can show the operator ∝ b indeed to be strongly suppressed below, we will neglect the former in the following.

Equation (34) then leads to the vector-meson masses
which allows us to extract the coupling constants according to In particular the value for b depends quite sensitively on the precise values inserted for the masses of the broad ρ and K * resonances; if we employ the real parts of their pole positions [33,70], we find a = 0.57(1) GeV −1 , b = −0.045(20) GeV −1 , δ = 20(2) MeV. We hence conclude that the formal 1/N c suppression of b/a translates, in fact, into a numerical suppression by more than an order of magnitude; we will therefore neglect b, too. Furthermore, we observe that the determination of a based on SU (3) symmetry leads to an estimate that is about 20% smaller than the value deduced from the one-loop IAM representation, see Eq. (12). As we expect that SU (3) breaking effects ought to affect the relation between ρ and ω observables less, we use the arguments above to employ a pion-mass dependence of M ω that equals the one of M ρ up to the constant offset δ, hence At higher orders in the chiral expansion, Goldstoneboson loops induce nonanalytic dependencies of the ω π π π ρ Fig. 3 VMD tree-level diagram for ω → 3π given by a ω → ρπ and a subsequent ρ → ππ decay.
vector meson masses of the form O(m 3/2 q ) and O(m 2 q log m q ), which have been studied extensively in the literature [67][68][69][71][72][73]; such terms will obviously break the similarity in M ρ (M 2 π ) and M ω (M 2 π ) due to the different coupling of ρ and ω to pions. We ignore such terms in the present study solely based on the observation in Sect. 2.1 that a linear dependence of M ρ on M 2 π is sufficient to describe the behavior of the ρ pole of the O(p 4 ) IAM amplitude.

Subtraction constant
As derived in Sect. 3, we require one subtraction constant in the dispersive representation of the ω → 3π decay amplitude in order to maintain a convergent integral representation. Since this subtraction constant is not fixed by unitarity (and, for the process at hand, cannot be matched to ChPT as for other processes such as γπ → ππ [21,52]), we need to fix its pion-mass dependence in a different way.
In Sect. 2 we have recounted that the coupling g ρππ at the ρ → ππ vertex is (essentially) pion-mass independent and in good agreement with a narrowwidth formula or a VMD model. In an isobar model of subsequent two-body decays, ω → 3π is typically understood in terms of processes ω → ρπ, followed by ρ → ππ decays, see Fig. 3 (cf. e.g. Refs. [74][75][76]). Reducing the dispersive representation Eq. (31) to such a simplified picture, we find the subtraction constant α in one-to-one correspondence with the product of coupling constants g ωρπ ×g ρππ . We therefore conjecture that, by analogy, it is reasonable to assume g ωρπ , and hence α, to be also pion-mass independent, and we fix the subtraction constant to the total decay width Γ ω at the physical point.

Isospin-breaking effects
Up to now all calculations have been carried out under the assumption of isospin symmetry (M π ± = M π 0 ). We now want to briefly discuss the influence of isospin breaking effects in ω → 3π; a more detailed discussion can be found in Ref. [77]. It is well known that in the context of precision analyses of η → 3π, in particular comparing η → π + π − π 0 and η → 3π 0 , taking into account the pion-mass difference for the available phase space at least is mandatory [78][79][80][81]. In this section, we therefore only investigate isospin breaking in the kinematical contribution, ignoring all dynamical effects (i.e. using |F (s, t, u)| 2 = const.), expecting this to be the dominant change.
The pion-mass difference ∆ π = M 2 π ± − M 2 π 0 = 1.26116(15) × 10 −3 GeV 2 (39) originates from two sources: electromagnetic effects and the light-quark-mass difference (m u − m d ). At leadingorder in ChPT (for three flavors), the latter can be evaluated to hence this effect is very small: the pion-mass difference is dominantly caused by electromagnetism, and we neglect the effect of the difference of the light quark masses completely. Consequently, as long as we neglect higher-order corrections of O(e 2 m q ), the pionmass difference ∆ π stays constant when varying M π . This allows us to relate the neutral pion mass to the charged one according to Obviously, this relation breaks down at a minimal charged pion mass of M π ± = √ ∆ π ≈ 35.5 MeV. It turns out that isospin breaking in the kinematical dependence of Γ ω as given in Eq. (33) when tuning the charged pion mass gives only a tiny correction to Γ ω of less than 2%. Due to the connection of the pion masses in Eq. (41) it is easy to see that the effect of isospin breaking becomes smaller when increasing the mass of the pions; on the other hand, since the ratio of the ω mass to the pion masses becomes larger when approaching the chiral limit, the increasing isospinbreaking effects are lifted by the ω mass. Thus using the isospin limit is entirely justified for our purposes.

Results
In this section we discuss the final results of our dispersive representation of ω → 3π. First of all we want to compare the resulting pion-mass-dependent Omnès function or two-body final-state interaction 3 (2-body  FSI) to the single-variable amplitude F α (3-body FSI) to study the dynamical effects generated by the interaction with the third pion as depicted in Fig. 4.
Besides pion masses close to the chiral limit, the 3body FSI leads to an enhancement of the modulus of the single-variable amplitude compared to the Omnès function. While the peak position in the absolute values (due to the ρ resonance) is essentially identical for 2and 3-body FSI at a given value of M π , the phases behave rather differently. The Omnès function fulfills Watson's theorem and thus its phase is identical to the ππ scattering phase shift (the IAM phase). Due to the dynamical effects stemming from the interaction with the third pion (and the generation of a three-pion cut), this does not hold for the single-variable amplitude (see Sect. 3). Thus the argument of the single-variable function is shifted compared to the input IAM phase. The lower and upper phase-space boundaries are marked by dashed vertical lines, and hence denote the kinematical range directly accessible in the decay. Already here we want to point out that mainly the tails of the ρ resonance will only contribute to the dynamics when increasing the pion mass (cf. Fig. 6).
In order to study the dynamical effects on the decay width we consider three different scenarios: first we consider only kinematic contributions to the decay width (all pion-pion dynamics are disregarded, δ(s) = 0, and thus |F (s, t, u)| 2 = const.), secondly we allow for 2-body rescattering effects (meaning |F (s, t, u)| 2 ∝ |Ω(s) + Ω(t) + Ω(u)| 2 ), and third the full 3-body dynamics are taken into account. A comparison of the different cases is displayed in Fig. 5. Pion-mass-dependent phase space boundaries of ω → 3π (solid lines) and position of the "on-shell" ρ in the respective s-, t-, and u-channels (dashed lines). The color coding is kept identical to Figs. 1 and 2 with the following pion masses: 0 (blue), 1 2M π (cyan),Mπ (green), 3 2M π (yellow), 2Mπ (orange), and 5 2M π (red). For the last two the phase space already vanishes, and thus the ω becomes stable with respect to the investigated decay mode.
First of all we notice that the ω → 3π width decreases with increasing pion mass. This is not surprising since the mass of the three pions is increasing faster than the mass of the ω, cf. Eq. (38). Hence the phase space shrinks for larger pion masses as depicted in Fig. 6. For M π > 1 3 M ω ≈ 1.96M π the masses of the three pions exceed the ω mass, thus the reaction ω → 3π is no longer allowed and the ω becomes stable with respect to the considered decay channel (i.e., in QCD in the isospin limit).
We now study the effects of the 2-and 3-body dynamics. First of all we notice that the kinematical prefactor given in Eq. (24) vanishes at the phase space boundaries in all directions. This leads to a stronger weighting of the inner region compared to the outskirts of the phase space when evaluating the integral Eq. (33). Secondly, the dynamics are mainly governed by the ρ resonance in the respective 2-body channels, leading to a three-band structure (one for each channel) in the function F (s, t, u). The peak position and width of these bands will be determined by the respective 2-and 3-body FSI effects and the information of the scattering phase input from the IAM. We conclude that strong imprints of the ρ as a dynamical effect will only affect the decay width close to the chiral limit, since the ρ is only allowed to go on-shell within the phase space boundaries for M π < 0.15M π . This effect is even reinforced by the increasing mass of the ρ and its decreasing width with growing pion mass, as investigated in Sect. 2. Thus the main contributions to the ω width will come from the tail of the ρ resonances at low pion masses, see Fig. 6. This is consistent with the result given in Fig. 5. We observe that the impact of the 2-and 3-body dynamics on the decay width is very strong close to the chiral limit, where it leads to a reduction of the decay width compared to the one of pure kinematics. The damping due to the 3-body dynamics is four times stronger than the influence of the 2-body effects and results in a width smaller by about one third than the purely kinematic effects. In the chiral limit we find for the widths of the three discussed scenarios, where we tookΓ ω ≡Γ (ω → 3π) = 7.57(9) MeV [29]. Since the subtraction constant of all three curves is fixed at the physical point, the difference between them shrinks when approaching this point. Above the behavior is opposite, here the 2-and 3-body dynamics generate a decay width that is larger than the one obtained from pure kinematic effects. Since the phase space at the physical point and above does not allow for on-shell ρ mesons in the respective 2-body systems, the influence of the dynamical effects above this point is very small. Thus the 2-and 3-body FSI curves do not differ strongly from the curve of pure kinematics. Considering Eq. (9) we are in the position to predict a pion-mass-dependent trajectory of the ω pole position in the complex-energy plane of 3π → 3π scattering similar to the ρ resonance in Fig. 1. A plot of this pole trajectory is given in Fig. 7. Note that, strictly speaking, our formalism does not determine such a complex pole position: rather, we have discussed mass (via effective Lagrangians) and width (via the Khuri-Treiman formalism) individually. Given the smallness of the width of the ω, we regard the error committed thereby as negligible, although this may not obviously be so in the general case of an arbitrary 3π resonance.
Since dynamical effects in the investigated process are limited due to the small phase space, we want to emphasize that the dispersive representation derived in Sect. 3 is valid for general V → 3π decays, with V denoting an arbitrary isoscalar vector meson, at least to the extent that the elastic approximation in the pion-pion rescattering is justifiable. Thus we are able to describe φ → 3π by just replacing the decay mass according to M ω → M φ within the same formalism [51,53]. With a mass of 1019 MeV [29] its decay allows for much richer dynamics due to the larger phase space, in particular the three ρ bands are visible inside the Dalitz plot [82,83]. In the quark-model picture, the φ can be understood as a pure ss state, thus its mass does not depend on the pion mass at leading order. This follows from Eq. (36) when using the GMOR relation to relate M 2 K = B(m + m s ) + O(m 2 q ). The results for the pion-mass-dependent partial decay width Γ (φ → 3π) are displayed in Fig. 8. Since the ρ is allowed to go onshell up to M π ≈ 1.7M π , the dynamical imprints are much stronger compared to the case of the ω → 3π even above the physical point.
Overall, given how subtle rescattering effects beyond two-body rescattering are usually thought to be (and more often than not neglected altogether in experimental Dalitz plot analyses), it is remarkable to see that in both cases studied here, ω → 3π and φ → 3π, the 3-body FSI effects tend to affect the quark-mass dependence of the (partial) widths about as strongly as the 2-body FSI.
We wish to add a few caveats concerning the precision of the predictions shown in this section. Error bands are estimated solely based on the uncertainty within the one-loop IAM phase shift input, as well as the one in M ω (M 2 π ), which is nonetheless based on the large-N c expansion and leading-order symmetry breaking in the quark masses. We do not attempt to quantify corrections due to higher orders in either case, and refer to Ref. [34] for future work concerning an improved pion-mass-dependent phase shift. Furthermore, a highprecision measurement of the φ → 3π Dalitz plot [82] revealed the need to include a second subtraction in the dispersive representation of the decay amplitude [51], which has also not been considered here; experimental data on the ω → 3π Dalitz plot is not conclusive in this respect yet [84].

Summary
We have investigated the pion-mass (and hence the quark-mass) dependence of the ω → 3π decay width, generalizing previous studies of two-pion resonances based on the inverse amplitude method. To this end, we have employed a dispersive formalism, based on Khuri-Treiman equations, that uses inverse-amplitudemethod phase shifts as input. The pion-mass dependence of the ω mass is estimated using a symmetry relation based on chiral perturbation theory for vector mesons. Deviations from phase space behavior alone, induced by the pion-mass-dependent decay amplitude, are clearly visible, although suppressed for larger-thanphysical pion masses due to the smallness of phase space. We have demonstrated for the decay of the heavier φ into three pions that this need not be the case in general. Remarkably, a simple description of the decays in terms of two-body rescattering alone does not yield a good approximation to the full pion-mass dependence seen.
The three-pion decays of the lightest isoscalar vector mesons only serve as a paradigm case for the investigation of three-body resonances; extensions to other, similar decays within the same formalism ought to be tested in the future.