The $\pi f_0(500)$ decay of the $a_1(1260)$

We evaluate the $a_1(1260) \to \pi \sigma (f_0(500))$ decay width from the perspective that the $a_1(1260)$ resonance is dynamically generated from the pseudoscalar-vector interaction and the $\sigma$ arises from the pseudoscalar-pseudoscalar interaction. A triangle mechanism with $a_1(1260) \to \rho \pi$ followed by $\rho \to \pi \pi$ and a fusion of two pions within the loop to produce the $\sigma$ provides the mechanism for this decay under these assumptions for the nature of the two resonances. We obtain widths of the order of $13-22$ MeV. Present experimental results differ substantially from each other, suggesting that extra efforts should be devoted to the precise extraction of this important partial decay width, which should provide valuable information on the nature of the axial vector and scalar meson resonances and help clarify the role of the $\pi\sigma$ channel in recent lattice QCD calculations of the $a_1$.


I. INTRODUCTION
The nature of the axial vector meson resonances is a subject of continuous debate. It is well known that quark models, very successful in correlating a vast amount of data concerning meson spectra, have difficulties reproducing the properties of axial vector mesons [1,2]. On the other hand, the advent of the chiral unitary approach, combining chiral dynamics with a unitary approach in coupled channels, showed that the lowest-energy axial vector mesons can be fairly well reproduced from the interaction of vector mesons with pseudoscalars [3][4][5][6]. The interaction used between vector-pseudoscalar (V P ) channels is taken from the chiral Lagrangian of Ref. [7]. In Ref. [8], for the particular case of the ρπ interaction, and in Ref. [9] for the general V P interaction, it was shown that this Lagrangian can be obtained from the exchange of vector mesons within the local hidden gauge approach [8,[10][11][12]. This conclusion is in line with the findings in the pseudoscalar-pseudoscalar (P P ) interaction [13]. However,the P V interaction shows a peculiarity since one can make the transition P V → V P mediated by a pseudoscalar exchange of one P from the decay of V → P P . Fortunately, such contributions are very small compared to those from vector exchange, as recently shown in Ref. [9], something corroborated in the study of D decays into three mesons looking at meson final state interactions [14]. Thus, one can rely on the information for the axial vector mesons obtained within the chiral unitary approach with the chiral Lagrangian of Ref. [7]. Lattice QCD has also shown interest in describing the axial vector resonances. Concretely, the a 1 (1260) was already studied in Ref. [15] (see Ref. [16] in connection to this work). More detailed work on the a 1 (1260) and b 1 (1235) is shown in Ref. [17], and the b 1 (1235) is also investigated along this line considering the coupled channels πω, πφ in Ref. [18]. The first extraction of a three-body resonance from lattice QCD, using up to three pion operators and a suitable three-body finite-volume formalism was achieved recently for the example of the a 1 (1260) [19]. The finite-volume mapping (FVU) relies on unitarity for three particles [20] leading to the aforementioned pseudoscalar exchange which develops singularities in the finite volume [21]. The largest source of systematics in the a 1 (1260) pole extraction from lattice QCD is the role of the πσ channel [19]. The prediction of the partial decay width a 1 → πσ in the present paper is, therefore, essential to assess the role of such sub-dominant channels in upcoming lattice QCD calculations.
As these examples show, combined work using effective field theories, quark models, lattice QCD simulations and experiment will ultimately unveil the nature of the axial vector meson resonances, and work goes along in these direction. One of the successful strategies pursuing this aim is to look for particular decay modes of these states which reveal the coupling of the resonances to some bound meson-meson components. Let us take as an example the f 1 (1285), which appears as a K * K + c.c. single-channel state, bound by about 100 MeV [3,4,6]. The f 1 (1285) can naturally decay to K * K (due to the K * width), then the K * decays to Kπ and KK fuses into a 0 (980), which couples strongly to KK. This process constitutes a triangle mechanism that leads to the f 1 (1285) → πa 0 decay which is tied to the coupling of the f 1 (1285) to the bound K * K component. This work is done in Ref. [22] and a good agreement with experiment is found. In a similar way one can study the f 1 (1285) → πKK decay, and a good agreement with experiment is found in Ref. [23].
The triangle mechanisms have been widely used to test molecular components, with the information encoded in the resonance decay vertex of the first step, or on the fusion vertex of the last step. In Refs. [24,25] a triangle mechanism was used to explain the unexpectedly large isospin breaking of the η(1405) → π 0 f 0 (980) decay compared to η(1405) → π 0 a 0 (980) observed at BES [26] (see also following works of Refs. [27,28]). Recent work on triangle mechanisms to describe different processes can be seen in Refs. [29][30][31][32][33][34][35].
One particular case of a triangle mechanisms is the one that develops a triangle singularity (TS). This was early discussed in Refs. [36,37], which showed that a singularity in the triangle diagram emerged when the three particles in the loop are placed simultaneously on shell, while being collinear and the decay of the first step and fusion of the second one are possible at a classical level (Coleman-Norton theorem [38]). The topic has experienced a rebirth recently since one finds nowadays many physical examples of reactions which are interpreted in terms of a TS. A modern reformulation of the problem, intuitive and practical, is given in Ref. [39] and a recent review on the topic is given in Ref. [40].
Many examples of TS are discussed in Ref. [40]. Here we just mention two examples which have particular relevance. One of them is the TS developed in the a 1 (1260) decay into πf 0 (980), related to the topic that we face in this paper (the a 1 (1260) → πσ(f 0 (500)) decay) which was originally associated to a new resonance "a 1 (1420)" by the COMPASS Collaboration [41], but which was soon interpreted in terms of a TS [42][43][44] and acknowledged as such by the COMPASS Collaboration [45]. The other recent example is the case of the pp → π + d fusion reaction, which, as shown in Ref. [46], has an unexpectedly large strength due to a TS. This work, together with the one of Ref. [47], present an alternative explanation of the peak observed in the pn → π 0 π 0 d (pn → π + π − d) reaction, so far associated to a dibaryon "d * (2380)" [48].
It is interesting to recall that not all triangle diagrams develop a TS, actually these are exceptional cases, but it is easy to test if a triangle diagram develops a TS applying a simple equation (see Eq. (18) of Ref. [39]). We shall take advantage of this equation too, because it not only serves to see if there is a TS, but also to have a feeling of the strength of the loop by observing how far this equation is from being fulfilled.
In the present work we address the issue of the partial decay width of the a 1 (1260) into πσ(f 0 (500)). This can be seen for instance in the a + 1 (1260) → π + π 0 π 0 decay, where the π 0 π 0 will come from the σ meson, hence, this is part of the three pion decay of the a 1 (1260).
This latter problem has been extensively studied from different points of view and it is not our purpose to go through it again. A thorough account and discussions on the different approaches can be seen in the works of Refs. [49][50][51], but in none of those works the πσ decay mode has been isolated.
Experimentally there are only a few experiments tabulated in the PDG [52] (although "they are not taken for averages, fits, limts, etc") and the branching ratios are and in Ref. [57] this decay mode is reported as "seen". As we can see, the decay ratios are very different, and with large uncertainties. The data of Ref. [53] come from the τ − → ν τ π − π 0 π 0 decay measured at CLEO II. These data should be taken with certain caution since m σ = 860 MeV and Γ σ = 880 MeV are used in the substructure fits, which are very different from present values [52,58,59]. 1 Dalitz plots with a more realistic σ resonance 1 From Ref. [58] we find and from Ref. [59] We take values close to these as m σ = 458 MeV and Γ σ = 464 MeV, coming from the use of Eqs. (32), (33).
were predicted by the three-body framework of Ref. [60]. That work and Ref. [49] use two-body amplitudes for the sub-channels in contrast to earlier work [61] in this direction.
However, none of these frameworks or related theoretical studies [62,63] have performed a refit of the CLEO Dalitz distributions [53] which are among the most precise sources of information on the a 1 decay modes. The ALEPH data on the lineshape [64], in contrast, served for the determination of the pole position of the a 1 [19, 43, 61] but does not contain information on branching ratios.
In order to evaluate the a 1 (1260) → πσ decay width, we follow the procedure of Ref. [50] in the study of a 1 (1260) → π + π + π − decay, where the tree level a 1 (1260) → ρπ together with final state interaction of two pions via a triangle loop mechanism are considered. The formalism follows closely the one of Refs. [22,44] and is based on the chiral unitary approach for the axial vector mesons [4] and for scalar mesons [65]. The present work has to be considered as complementary to the one of Ref. [50], where the π + π − mass distribution is obtained, that is largely dominated by the a 1 → ρπ tree level mechanism, and that is shown to be consistent with the experimental data of ARGUS [66]. In the low energy part of the π + π − invariant mass spectrum, some extra strength appears which should be attributed to the π + σ production but the πσ decay mode was not isolated. The precise evaluation of this decay mode is the purpose of the present work.

II. FORMALISM
We follow closely the approach of Ref. [22] for the f 1 (1285) → πa 0 (980) decay. In that case the f 1 (1285) decays to K * K , the K * to Kπ and then the KK fuse to give the a 0 (980).
In the present case, we have the a 1 (1260) coupled to ρπ andK * K − K * K , with the largest coupling to ρπ [4]. In view of this, and similar to Ref. [50], the mechanisms to produce the σ are given by the diagrams of Fig. 1. One can see that none of these triangle mechanisms produces a TS [39], yet the diagrams Fig. 1 (b) and Fig. 1 (c) are highly suppressed with respect to diagram Fig. 1 (a). Indeed, in diagram Fig. 1 (a) one can place the ρ, π, π intermediate states on shell, yet the condition of the Coleman-Norton theorem [38] is not satisfied, since, in a collinear decay, the π coming from the ρ decay has a smaller velocity than the one from a 1 decay, hence, it cannot catch up with the latter pion to produce the σ. The processes in Fig. 1 (b) and Fig. 1 (c) are much weaker since, first, the a 1 →K * K through a 1 → K * K in the first step; (c) through a 1 →K * K in the first step. Fig 1 showing explicitly the ρ + π 0 and ρ 0 π + decays of the a 1 .
Particle momenta are shown in parentheses.
coupling is smaller than the a 1 → ρπ one [4], second, because the σ couples strongly to ππ but very weakly to KK [65] and, third, because the KK state is highly off shell in the loop if we wish to produce the σ. Hence, we shall not consider diagrams Fig. 1 (b) and Fig. 1 (c) in our evaluation.
It is curious to note that, conversely, the diagrams Fig. 1 (b) and Fig. 1 (c) are those responsible for the a 1 (1260) → πf 0 (980) decay, replacing the σ by the f 0 (980) [44]. This is because in this case the diagrams develop a TS, and also the f 0 (980) coupling to KK is much stronger than the coupling of the σ to KK [67]. This decay mode was the one originally assigned to the a 1 (1420) resonance, now accepted as a TS [45]. Yet, the contribution of this mode to the a 1 (1260) width is small, of the order of 1 MeV [44] which gives us an idea of the smallness of the diagrams in Fig. 1 (b) and Fig. 1 (c) in the much less favored actual situation of σ production.
We, thus, proceed with the diagram of Fig. 1(a) to produce πσ and we show in detail in Fig. 2 the diagrams to consider. We need the states of the σ and a + 1 in terms of π and ρ multiplets (−π + , π 0 , π − ), (−ρ + , ρ 0 , ρ − ), |ρπ, I = 1, We take the coupling g σ,ππ of the σ to ππ from the chiral unitary approach [65] with a single channel ππ in I = 0, but we also consider the ππ, I = 0 amplitudes in coupled channels.
Further details will be given at the beginning of Sec. III and Appendix A. In the chiral unitary approach, for convenience, one works with a unitary normalization where Eq. (3) has a coefficient 1/ √ 6 instead of 1/ √ 3. This √ 2 factor has to be restored at the end for the external pions so that Similarly, from Ref. [4], in the a 1 rest frame, the decay vertex reads where a 1 , ρ are the polarization vectors of the a 1 , ρ, with g a 1 ,ρπ = (−3795 + i 2330) MeV . Thus, Analogously, we have (see Appendix A) We also need the ρ → ππ Lagrangian which is given by with V, P the SU(3) matrices for the vector mesons and pseudoscalars, respectively [4], which gives where we keep only the spatial components of the ρ polarization, given the coupling of Eq. (6) and the fact that the 0 component provides only a very small contribution even for ρ mesons with the momenta involved in the a 1 → ρπ decay (see Appendix A of Ref. [68]).
With these ingredients we obtain the amplitudes corresponding to the diagrams of Fig. 2 (a1) and (a2): where Π is the product of propagators with ω 2 ρ = m 2 ρ + q 2 . In Eq. (14) we have kept the positive-energy part of the ρ propagator, corresponding to a heavy particle that will be mostly on shell in the loop, but keep the full propagators for the two pions in the loop. Separating the pion propagators into their positive and negative energy parts, where ω 2 ( q ) = q 2 + m 2 π , we can readily perform the q 0 integration of the amplitudes t (a1) and t (a 2 ) of Eqs. (12) and (13). We obtain where which, taking into account that can be cast as In Ref. [50] the integral of Eq. (19) is regularized in terms of a cut off in q of 630 MeV and an additional form factor. We regularize it with the intrinsic cut offs that stem from the chiral unitary approach. Indeed, as seen in Ref. [69], the scattering amplitudes appear with the factorization implying that the vertices are factorized as and q is the boosted value of the momentum q to the frame where the σ is at rest (see formula in Appendix B). Thus, finally, and, following Ref. [4], we take q max = 800 MeV and q max = 750 MeV (as discussed in Sec. III), by means of which we obtain the width for a 1 → πσ, averaging over the a 1 initial polarizations, where M I is the invariant mass of the σ and p π ≡ | p π | = λ 1 /2 M 2 a 1 , m 2 π , M 2 I / (2M a 1 ) here and in the following. Eq. (23) assumes a σ with fixed mass M I , but the σ has, in fact, a large mass distribution. We can formally take this into account (it will done more rigorously below) assuming a Breit-Wigner mass distribution and, thus, But now we see that formally |g σ,ππ | 2 (M 2 I −m 2 σ +i m σ Γ σ ) −1 is the Breit-Wigner representation of the ππ amplitude. Hence, we replace, which allows us to use the realistic scattering amplitude t I=0 ππ,ππ from Ref. [65], or from Ref. [70] where the coupled channels π + π − , π 0 π 0 , K + K − , K 0K 0 and ηη are used. Then, as found in Hence, we can write where the dependence on M a 1 is explicitly highlighted. Since the a 1 (1260) also has a large width we can do a convolution of the width of Eq. (28) with the mass distribution of the a 1 (1260) in analogy to Eq. (25). We use the most accurate results for the a 1 (1260) mass and width from the COMPASS collaboration [71], and obtain with N = To calculate the a 1 → πσ width we use the model of Ref. [65] with a single channel |ππ, I = 0 amplitude using unitary normalization (extra 1/ √ 2 in Eq. (3)), with and G the ππ loop function regularized with a cut off q max = 750 MeV, such as to get the I = 0, ππ phase shifts. The coupling that we obtain is given in Eq. (5). Thus, and we have, We have also carried out the calculations using the model of Ref. [70], considering explicitly the π + π − , π 0 π 0 , K + K − , K 0K 0 , and ηη channels, and using Eq. (27), and the results are practically the same except that including the extra channels demands the use of a different cut off, q max = 600 MeV. The results that we show are those based on Eqs. (32)- (35). We shall estimate uncertainties at the end using results of different models for ππ scattering.
In Fig. 3 we show first the results obtained with Eq. (23) as a function of M I . As expected for lower σ masses Γ grows since p π is bigger, and we find a cusp like peak at the two-pion threshold as it should be. We already see from Fig. 3 that around M I 500 MeV the width is of the order of 33 MeV.
Next we use Eq. (28) in terms of t I=0 ππ,ππ of Eqs. (32)- (35). We show the results in terms of M a 1 which we plot in Fig. 4. We see that Γ has been reduced, since around M a 1 = 1255 MeV it has a value about 14 MeV, around one half what we estimated from Fig. 3. We should explain this difference, which we do in Fig. 5. There we plot the imaginary part of the Breit-Wigner form of the σ, versus that of the realistic t I=0 ππ,ππ amplitude, Im t I=0 ππ,ππ . What we observe is that these two magnitudes are very different, both in strength and shape. In particular, the imaginary part of the full ππ amplitude goes to zero at threshold due to unitarity, as it should be, while the Breit-Wigner approximation does not. At higher M I the situation is opposite. This discrepancy should be stressed since in conventional experimental partial-wave analyses it is customary to use Breit-Wigner forms for the σ, although in some cases more sophisticated forms are taken, like in Ref. [72] based upon the formulation of Ref. [73]. Altogether we find about a factor of two difference in the evaluation of the width by both methods, the realistic considering the errors of Eq. (29).

A. Uncertainty analysis
The branching ratio of Eq. (38) reflects only the uncertainties in the width of Eq. (30). To estimate uncertainties from the use of other models we resort to a simple method by looking at the couplings obtained by different groups which are tabulated in Table 5 of Ref. [59]. We simply scale the result of Eq. (37) by the squared of the ratio of the couplings shown in that in Eqs.
(1), (2) that might come from treating the σ as an ordinary resonance as discussed before. The formulation presented here should help perform more accurate analyses of data in future analyses. Actually, with this work finished, there has been a recent experiment [72] for the D + s → π + π − π + η decay with emphasis on different partial wave contributions, and from where one can induce the branching ratio that we have evaluated. From Table I of that work one finds a weight for a 1 (1260) + η, with a + 1 observed in ρ 0 π + , W 1 = 55.4 ± 3.9 ± 2.0. Then a + 1 η, with a + 1 observed in f 0 (500)π + , has a weight of W 2 = 8.1 ± 1.9 ± 2.1. But this only accounts for f 0 (500) → π + π − . The π 0 π 0 decay has a strength of 1 /2 of that of π + π − . Furthermore, since the weight for a 1 (1260) decay into ρ + π 0 is the same as for ρ 0 π + , and the total branching ratio of a 1 → ρπ is 60% [52], we have where we have summed errors in quadrature. This number is in good agreement with our results. It is interesting to see that with our a 1 → ρπ coupling of Eq. (7), we obtain a branching fraction of a 1 to ρπ of 51%, close to the one reported in the PDG [52]. If, instead of the rather uncertain experimental a 1 → ρπ 60% branching ratio, we use our theoretical value of 51%, the numbers in Eq. (39) change to (5.6 ± 2)%. Hence, a range of Γ(a 1 → σπ)/Γ a 1 ≈ (3.6 − 9)% is a safe estimate from the experiment of Ref. [72].
We have evaluated the decay width of the a 1 (1260) resonance into πσ(f 0 (500)) taking into account that the a 1 (1260) is a state dynamically generated by the vector-pseudoscalar interaction, concretely with the channels ρπ and K * K , and the σ arises from the unitarized interaction of pseudoscalar-pseudoscalar mesons, mostly of two pions. Given this starting point, the mechanism that leads to the a 1 → πσ decay is given by a triangle mechanism in which the a 1 (1260) decays to ρπ, the ρ decays to ππ and one of these pions, together with the one coming from the a 1 → ρπ decay, fuse into the σ meson. The mechanism allows to have the intermediate particles on shell, but does not develop a triangle singularity, yet it provides a relatively large strength for that decay. We find decay widths in this channel of the order of 13 − 22 MeV, which amount to a branching fraction of (3.3 − 6.6)%. This is a small fraction and several experiments, so far largely contradictory, have been devoted to estimate that width. On the other hand, the recent BESIII experiment [72] allows to extract this branching ratio to be in the range (3.6 − 9)%, which overlaps with the range that we obtain.
One of the findings is that the width obtained depends much on how the mass distribution of the σ is considered. We found that treating the σ as an ordinary Breit-Wigner resonance, using realistic coupling, mass and width obtained in theoretical works, gives rise to results which differ substantially from those obtained by using the imaginary part of the I = 0, ππ scattering amplitude as the mass distribution. This should serve as a warning for future experimental analyses of this partial decay width. The work of Ref. [72] uses a more elaborate and realistic picture for the σ resonance extracted from Ref. [73].
Any information on the πσ and other subdominant branching ratios of the a 1 (1260) helps assess systematic effects in recent efforts to calculate the properties of three-body resonances in lattice QCD [19]. In particular, the smallness of that branching ratio found in this study strengthens the truncation to πρ channels made in the finite-volume frameworks of Refs. [16,19].
The results of this study are tied to the assumptions about the nature of both the a 1 (1260) and the σ resonances as being dynamically generated from the pseudoscalar-vector and pseudoscalar-pseudoscalar interaction. Given the ongoing discussion on the nature of the axial vector and scalar meson resonances, the precise determination of this width should provide valid information concerning this relevant issue, and it is our hope that the present work stimulates further efforts in this direction.