Combined theoretical study of the D+→π+ηη\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D^+ \rightarrow \pi ^+ \eta \eta $$\end{document} and D+→π+π0η\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D^+ \rightarrow \pi ^+ \pi ^0 \eta $$\end{document} reactions

We study the D+→π+ηη\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D^+ \rightarrow \pi ^+ \eta \eta $$\end{document} and D+→π+π0η\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D^+ \rightarrow \pi ^+ \pi ^0 \eta $$\end{document} reactions, which are single Cabibbo suppressed and can proceed both through internal and external emission. The primary mechanisms at quark level are considered, followed by hadronization to produce three mesons in the D+\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D^+$$\end{document} decay, and after that the final state interaction of these mesons leads to the production of the a0(980)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a_0(980)$$\end{document} resonance, seen in the π+η\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pi ^+ \eta $$\end{document}, π0η\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pi ^0 \eta $$\end{document} mass distributions. The theory has three unknown parameters to determine the shape of the distributions and the ratio between the D+→π+ηη\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D^+ \rightarrow \pi ^+ \eta \eta $$\end{document} and D+→π+π0η\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D^+ \rightarrow \pi ^+ \pi ^0 \eta $$\end{document} rates. This ratio restricts much the sets of parameters but there is still much freedom leading to different shapes in the mass distributions. We call for a measurement of these mass distributions that will settle the reaction mechanism, while at the same time provide relevant information on the way that the a0(980)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a_0(980)$$\end{document} resonance is produced in the reactions.


Introduction
The D weak decays into three light mesons have proved to be very rich, allowing one to dig into the weak reaction mechanism [1][2][3] as well as providing information on the meson-meson interaction [4][5][6][7][8][9][10][11][12][13][14][15][16][17]. A review on this latter issue can be seen in Ref. [18]. The rich field of mesonmeson interactions gives rise to many mesonic resonances that show up in most of the D decays. Normally, the possible pairs with the tree final mesons lead to some mesonic resonances and it is common to see the effects of several resonances in just one decay. For instance, in the D 0 → K − π + η a e-mail: ikeno@tottori-u.ac.jp (corresponding author) b e-mail: melahat.bayar@kocaeli.edu.tr c e-mail: oset@ific.uv.es reaction [19], one has the contributions of the a 0 (980), κ(K * 0 (700)) and K * (890), among other resonances that only have a minor effect in the mass distributions [20]. Another example would be the D 0 → ηπ + π − reaction measured at BESIII [21], where one expects contribution from the ρ 0 , a 0 (980), f 0 (500) and f 0 (980) resonances. It is interesting to look for reactions which show only effects of a single resonance. This was the case of the D + s → π + π 0 η reaction measured at BESIII [22], which was shown to be dominated by the a 0 (980), seen in the π + η and π 0 η mass distribution [23]. Actually in this reaction, the π + π 0 can also come from the ρ + , but in the experiment [22] this contribution was eliminated with a simple cut, demanding that M(π + π 0 ) > 1 GeV.
In the present paper, we want to study two reactions, also measured in Ref. [21] although without information on mass distributions: the D + → π + π 0 η and D + → π + ηη reactions. In the first one we expect to have contributions of the ρ + and a 0 (980) and in the second one of the a 0 (980). As in Ref. [22], the ρ + contribution in the D + → π + π 0 η reaction can be eliminated with the same cut, requiring that M(π + π 0 ) > 1 GeV, and then the two reactions can be related and will show only the a 0 (980) resonance. Actually, one of these reactions, the D + → π + π 0 η, has been studied theoretically in Ref. [24] and shown to be Cabibbo suppressed, exhibiting clearly the a 0 (980) excitation. We plan to study the two reactions together, along the lines of Ref. [24], using the extra information provided in Ref. [21] about the ratio of branching fractions of these two decay modes. This information puts constraints on the free parameters of the theory and allows us to make predictions on mass distributions which can be tested in future runs of the reactions with more statistics. Although the final state in the D + → π + π 0 η reaction is the same as in D + s → π + π 0 η measured in Ref. [22], the reaction mechanism is different here since it allows contributions from external emission and internal emission [25], while in the D + s → π + π 0 η external emission did not occur and the process proceeded via internal emission [23]. In addition, the D + → π + π 0 η reaction is Cabibbo suppressed and there are two topologically different mechanisms to proceed for each of the different modes. This necessarily introduces more unknowns in the theoretical description, in spite of which it was shown in Ref. [24] that the a 0 (980) signal should show up clearly in the reaction. The data on the branching ratios for D + → π + π 0 η and D + → π + ηη will limit somewhat the theoretical freedom and allow us to make more constraint predictions in both reactions.
Let us begin with external emission. We can have the primary mechanisms at the quark level depicted in Fig. 1.
At the same time we can also have internal emission which is depicted in Fig. 2.
Next we proceed to hadronize those mechanisms introducing aqq pair SU(3) singletūu +dd +ss, which is depicted in Figs. 3 and 4.
In Fig. 3a, we will have the hadronization of the us pair as where M is the qq matrix, which we write in terms of physical mesons as where we have used the standard η, η mixing of Ref. [26] and omitted the η which has a very high mass and does not play a role in the generation of the low energy scalar resonances [27]. By looking at the diagrams in Figs. 3 and 4 we see which is the qq pair that is hadronized, and following the method of Eq. (1) we see the pairs of mesons which are produced in each case. The mesons that appear in the hadronization are Since Fig. 3a, c have the same topology and the same Cabibbo suppressing factor, we can sum them and the same can be said about Fig. 3b, d. For internal emission, we can also sum the contributions of Fig. 4a, c and Fig. 4b, d.
However, there is a subtlety about the diagrams of type Fig. 3a, c. The effective W P P (P pseudoscalar meson) vertex can be evaluated with effective chiral Lagrangians W μ [P, ∂ μ P]T − with W μ the W field and T − a matrix related to the Cabibbo-Kobayashi-Maskawa elements [28,29]. If we wish to get the two pseudoscalar mesons in s-wave, which we need to produce the scalar resonances, we get such a contribution with this Lagrangian with μ = 0, which pro-duces a vertex proportional to p 0 − p 0 in the rest frame of W , and hence vanishes for particles with equal mass. This is discussed in Ref. [30]. This means we can keep the π K and πη terms in Fig. 3a, c but we must omit the K +K 0 contribution in Fig. 3c. These mechanisms were neglected in Ref. [24], but we keep them here and play with the relative weight with respect to the hadronization in Fig. 3b, d in order to get agreement with the branching ratios of D + → π + π 0 η and D + → π + ηη [21]. We should make some comments at this point because this mechanism, giving contribution when the masses of the pseudoscalar mesons are different, produces a violation of SU(3) and one may wonder why not other sources of SU(3) breaking are considered. At this point we must recall that chiral Lagrangians are SU(3) invariant except for a mass term where the physical masses are used, breaking then SU(3). Keeping the right masses is also important when it comes to thresholds and phase space of reactions. It is then quite customary in this kind of problems to assume SU(3) symmetry except in terms which involve directly the masses of the particles. Some tests of SU(3) breaking have been done in other weak decays hadronizing withūu +dd + (1 − α)ss, finding that even a 20 % breaking of SU(3) symmetry, taking α = 0.20, leads only to 5 % changes in the decay widths [31]. For all these reasons we keep these terms, but we also perform a study removing them and discuss the consequence on the mass distributions.
This said, we obtain from hadronization of the upper qq pair of the external emission in Fig. 3a c the hadronic contribution and from hadronization of the lower qq pair (Fig. 3b, d) Similarly, summing the contributions of Fig. 4a, c for internal emission, we find and summing the contribution of Fig. 4b, d In the former equations, H 1 , H 2 , H 1 , H 2 symbolize the trios of mesons that stem from Figs. 3 and 4 upon hadronization of a given qq pair, prior to final state interaction, and can be considered as the tree level terms that come from the weak decays. We can further simplify these expressions. Since K − K + + K 0K 0 has isospin I = 0, the combination (K − K + + K 0K 0 )π + cannot give rise to π 0 ηπ + upon rescattering since π 0 η has I = 1. It could contribute to ηηπ + production but ηη is far away from the narrow f 0 (980) resonance and will also be ineffective in this channel. Similar arguments can be done with the contribution π + π − + π 0 π 0 2 which has I = 0 and hence cannot give π 0 η upon rescattering. It could give ηη but the same argument as above holds. This means that far practical purposes we can write The same arguments can be used concerning H 1 and H 2 and we can then rewrite for effective purposes and H 1 , H 2 are not changed. We should note that eliminating the ηηπ + terms in H 2 , H 1 , H 2 , which do not contribute to the π + π 0 η production studied in Ref. [24], we obtain the same terms as in Ref. [24]. In order to keep a close analogy to the results of Ref. [24], we give weights to the different terms; As discussed in Ref. [24], one expects B ∼ 1 3 A, as it corresponds to the internal emission color suppression, and γ ∼ 1. We shall also assume values around these ratios taking A = 1 with arbitrary normalization.
The next step consists of allowing the interaction of pairs of particles to finally have either the π + π 0 η or π + ηη final state. The possible ways to have π + π 0 η are given in Fig. 5.
Similarly, if we wish to produce π + ηη, we will have the mechanisms depicted in Fig. 6.
Note that in Fig. 6 we are neglecting K 0K 0 → ηη because the f 0 (980) resonance strength already becomes negligible above the ηη threshold. By looking at the diagrams in Figs. 5, 6, we can write the final amplitude for the two reactions. For this purpose we define the weight of the different terms as: Rescattering of terms to produce π + π 0 η And the amplitudes are now written as, Fig. 6 Rescattering of terms to produce π + ηη where in t D + →π + ηη we have taken into account the factor 2 of symmetry in the amplitude because of the two η and included a factor 1 √ 2 such that when squaring t we get the factor 1 2 of symmetry in the width for two identical particles in the final state.
The differential width, up to an arbitrary normalization, is given by Ref. [32], and according to Figs. 5, 6 we have π + (1), η(2), π 0 (3) for π + π 0 η production and η(1), π + (2), η(3) for π + ηη production. The single differential mass distribution d d M inv (12) is obtained integrating over M 23 with the limits of the Dalitz plot that are shown in the PDG [32]. The scattering matrices are calculated using the chiral unitary approach [27] with the coupled channels, K + K − , K 0K 0 , π 0 η using with the potential V i j between the channels given in Ref. [33] and the G function regularized with a cut off q max = 600 MeV (q max for | q | in the d 3 q integration of the loop G function ) as done in Refs. [33,34]. Since with these channels we only get scattering matrices in the neutral states, taking into account the I = 1 character of all these amplitudes which involve πη, and the phase convention of the isospin multiples (K + , K 0 ), (K 0 , −K − ), (−π + , π 0 , π − ), we find We should note that while the t i j scattering matrices appearing in Eqs. (22), (23) are dominated by the a 0 (980) resonance, they contain more than the pole term since the chiral unitary approach produces the full scattering matrix and not only the resonance part of it. Another small technical detail is that since our amplitudes are good up to about M inv = 1200 MeV, we make a smooth extrapolation of Gt at higher energies, as done in Refs. [20,23,35] and the results barely change for different sensible extrapolations.
In the calculations we consider the t i j scattering matrices with πη in the final states in s-wave. As discussed in the study of the pp → d K +K 0 and pp → dπ + η reactions in Ref. [36], one interesting property of the πη system reflected in chiral dynamics is that it does not couple in p-wave to lowest order in the chiral counting [37]. It neither couples to vector mesons as shown in Ref. [38]. It can couple in higher orders but such effects are suppressed by more than one order of magnitude with respect to the dominant s-wave [37]. The π + π 0 pair could come in p-wave through the ρ + η production channel. We shall discuss this point in the next section anticipating that its contribution can be removed with a cut M inv (π + π 0 ) > 1 GeV, as has been done by the BESIII collaboration in the study of the D + s → π + π 0 η reaction [22] 3 Results We have seen that we have four parameters A, β, B, γ . Aβ gives the strength of the hadronization of the upper vertex in external emission (Fig. 3a plus c). A the strength for hadronization of the lower qq components of external emission ( Fig. 3b plus d). B the strength for the hadronization of the upper qq pair in internal emission (Fig. 4a plus c) and Bγ the strength for hadronization of the lower qq component in internal emission (Fig. 4b plus d). One parameter, we take A for this purpose, provides an arbitrary normalization and we take it 1 or −1, since we only evaluate the shapes of the distributions and the ratios of the widths for D + → π + ηη and D + → π + π 0 η. For the second reaction, the PDG [32] provides the branching ratio obtained from CLEO [39] The D + → π + ηη is measured for the first time in BESIII [21], where also the D + → π + π 0 η decay is studied and the following branching ratios, based on improved statistics, are reported We can see that the result for B(D + → π + π 0 η) of CLEO (Eq. (27)) and BESIII (Eq. (29)) are different, and even incompatible counting errors, although they are close.
Our strategy is to find a set of three parameters β, B, γ that provide a ratio of R = B(D + → π + ηη)/B(D + → π + π 0 η) = 1.33 ± 0.16, (30) summing relative errors in quadrature. We shall then search for parameters that give the ratio of Eq. (30) between 1 − 2 if possible. First we start from A = 1 and take values of B, β, γ in the range: The reason for the range of parameters is that we expect B to be suppressed by the number of colors, and then B ∼ 1 3 . The relative sign between external emission and internal emission is favored to be positive in analyses of c → pπ + π − measured at BESIII [40] and B + → J/ψωK + measured by the LHCb collaboration [41]. Yet, the present process is different, and as in Ref. [24] we shall explore the results with opposite sign. We should also warn at this point that the D + → π + π 0 η process will have contribution from ρ + η. This can come from the diagram of Fig 1b when  will cause a reduction of this mechanism, but the absence of hadronization of the ρ + η channel can compensate for that Fig. 7 Solid lines: Results for d /d M inv (π + η) for D + → π + π 0 η (upper curve) and D + → π + ηη (lower curve). The dashed curves stand for the phase space of D + → π + π 0 η (upper curve) and D + → π + ηη (lower curve). The set of parameters used, A = 1, β = 1, B = 1/3, γ = 1.
and we anticipate a clear contribution of the ρ + η channel comparable to the case that we evaluate based on s-wave process that will generate the a 0 (980) resonance with a fairly large strength. However, there is a way to remove the ρ + η contribution, imposing a cut M inv (π + π 0 ) > 1 GeV, which was already used in D + s → π + π 0 η [22], which facilitates the comparison with our results. This was the case when comparing those results of Ref. [22] with theory based upon only s-wave interaction in Ref. [23], where a very good agreement was found. Unlike the work of Ref. [23], based on only internal emission, which had no free parameters up to an arbitrary normalization, the present reactions contain both external and internal emission and we have three free parameters.
Before we do the fit to the data, we find illustrative to see our results with a standard set of parameters: A = 1, β = 1, B = 1 3 , γ = 1. The results of d d M inv for D + → π + π 0 η and D + → π + ηη are shown in Fig. 7. The results are shown in arbitrary units but all figures share the same normalization. We also plot there the shape of the phase space calculated taking only the term 1 in Eqs. (22), (23). What we see is that the difference between the shape of phase space and the results of our model is huge, indicating the important role played by the final state interaction of the meson components originated in the first step. This should be sufficient to show the value of these reactions to provide information on the meson-meson interaction. In Fig. 8, we also show the results with the same parameters for the π 0 η distribution in the D + → π + π 0 η decay. The strength at the peak is similar but not the same, which is clear from Eq. (22), but the features are qualitatively similar. From now on we shall only plot results for the π + η distribution.
With the caveats discussed above, we make a survey of the ratio R with different parameters at the range of Eq. (31) and we observe that it is impossible to get a ratio bigger than  (π 0 η) for D + → π + π 0 η. The dashed curve stands for the phase space of D + → π + π 0 η. Same set of parameters as in Fig. 7. 0.5, about a factor three lower than the ratio of Eq. (30). One of the reasons is that the D + → π + ηη reaction has a more reduced phase space than D + → π + π 0 η. In Fig. 9, we show the results for d /d M inv (π + η) with the sets of parameters, In the range of these parameter R varies between R ∈ [0 : 0.46].
In Fig. 9, we see that the D + → π + ηη reaction has a neat a 0 (980) signal, with the sharp peak corresponding to the cusp-like shape of the a 0 (980). The shapes for D + → π + π 0 η are rather different, with a broad bump at low invariant masses, caused by the tree level contribution, and much strength at high invariant masses, caused again by the tree level and the π 0 η producing the a 0 (980). This broad contribution at higher invariant masses was also visible in the D s → π + π 0 η experiment of Ref. [22]. This large contribution, in a region of phase space not allowed in D + → π + ηη must be seen as the main reason on why it is difficult to get ratio R bigger than 0.5.
The important point concerning the results of Fig. 9, is to note that in all cases and for the two reactions the peak of the a 0 (980) is clearly seen, and in our approach we have not introduced the resonance by hand, but comes automatically from the rescattering of the particles, and not only the final particles observed, but also rescattering of KK pairs produced in a first step. This conclusion is the same one reached in Ref. [24].
In view of the impossibility to obtain a ratio R bigger than 0.5 with A = 1, we take now A = −1 with the other parameters in the same range of Eq. (31). In this case we can obtain bigger ratios than before. We select three parameter sets, In the range of these parameter R varies between R ∈ [0 : 3.44].
The mass distributions for these sets are shown in Fig. 10. The case (a) in Eq. (33) leads to a big ratio R = 2.22 but at the expense of having very small individual rates due to large cancellations. What we observe here is that the shape of the a 0 (980) is much more clearly seen in the D + → π + ηη reaction than in the D + → π + π 0 η one. In this latter case the a 0 (980) is barely seen as a small cusp at the KK threshold.
What we see is that the shape of the π + η mass distributions depends strongly one the set of parameters that we take. The other reading of this finding is that the measurement of the mass distributions in both reactions, which is not done so far, will provide information on the reaction mechanisms for these two decays. The remarks obtained here should be a motivation for further measurement of the reactions, Fig. 11 Results for D + → π + π 0 η with and without the cut M inv (π + π 0 ) > 1 GeV for the set A = 1, β = 3, B = 0.15, γ = 0.33 with larger statistics that allow the mass distributions to be obtained.
To finalize the work, we show in Figs. 11 and 12 the results for two selected distributions of Figs. 9 and 10 when we make a cut demanding that M inv (π + π 0 ) > 1 GeV. Since our variables are M inv (π + η), M inv (π 0 η), we use the relationship, In Fig. 11, we choose the results of Fig. 9 for set (a) of Eq. (32) for D + → π + π 0 η and plot the corresponding results without the cut and with the cut of M inv (π + π 0 ) > 1 GeV. We can see that the effect is a reduction of the contribution at large M inv (π + η), similar to what was found in Ref. [24], and also at low invariant masses.
In Fig. 12, we show the same results for the set (c) of Fig. 10, the parameters of (c) of Eq. (33). The effect of the cut is similar to the former case, with the reduction of the strength at larger M inv (π + η), and also at low invariant masses. Finally, we perform a new calculation eliminating the β term which explicitly broke SU(3) symmetry through a term distinguishing the masses of the pseudoscalar mesons belonging to the SU(3) octet. Now we only have the B and γ parameters to vary. The first thing to note is that we do not find any solution for R 1.35 in the range of Eq. (31). All the values of R obtained are below R = 0.24 . We extend the range a bit to find a solution and we obtain A = −1, B = −0.68. γ = −1.5 for which R = 1.35. We plot now in Fig. 13 the mass distributions corresponding to Fig. 10 with the set of parameters (b) which gives also R = 1.35. As we can see, the results for the πη mass distribution in the case of π + ηη are very similar in strength and shape. The πη mass distributions for the π + π 0 η case differ somewhat, although the features are similar, with the a 0 (980) showing up as a cusp and a larger strength at larger invariant masses. From these pictures, it looks clear that while the peak of the a 0 should be visible in all cases, the π + ηη final state shows the a 0 (980) much clearer.

Conclusions
We have studied the D + → π + ηη and D + → π + π 0 η reactions from a perspective where the possible mechanisms at the quark level are considered with unknown strength. Both reactions are single Cabibbo suppressed and we find that both the internal and external emission mechanisms are possible. After this, hadronization of the qq pairs of the different mechanisms is considered, keeping not only π + ηη or π + π 0 η but also other intermediate states which upon rescattering can produce these final states. The pairs of mesons obtained after the hadronization are let to interact, using the chiral unitary approach to account for that interaction. The reactions selected are relevant because one can only produce the a 0 (980) resonance upon rescattering and the f 0 (500) and f 0 (980) do not show up. This makes the reactions useful to learn about the interaction of pseudoscalar mesons in the scalar sector with isospin 1.
In order to constrain the values of the parameters we used the ratio R of the branching ratios of D + → π + ηη and D + → π + π 0 η. This indeed puts much constraints on the parameters, but we still had freedom to obtain values of R bigger than 1 for different sets, leading to different shapes of the πη mass distributions. We observed that, in all cases the a 0 (980) signal was visible in the π + η invariant mass distributions, but was more neat in the D + → π + ηη reaction. We conclude that, while we clearly expect to see the a 0 (980) signal in the πη mass distributions, there are still uncertainties in the theory concerning the actual shape tied to the details of the reaction mechanism. In this sense, when the actual mass distributions are measured, information will be available that allows us to pin down the reaction mechanisms and come out with more assertive conclusions on the role played by the a 0 (980) resonance in these reactions.