Triangle mechanism in the decay process $B_0 \to J/\psi K^0 f_0(980)(a_0(980))$

The role of the triangle mechanism in the decay process $B_0\to J/\psi K^0f_0 \to J/\psi K^0\pi^+\pi^-$ and $B_0\to J\psi K^0a_0\to J/\psi K^0 \pi^0\eta$ is probed. In these process, the triangle singularity appears from the decay of $B^0$ into $J/\psi\phi K^0$ then $\phi$ decays into $K^0\bar{K^0}$ and $K^0\bar{K^0}$ merged into $f_0$ or $a_0$ which finally decay into $\pi^+\pi^-$ and $\pi^0\eta$ respectively. We find that this mechanism leads to a triangle singularity around $M_{\rm inv}(K^0f_0(a_0))\approx1520\ {\rm MeV}$, and gives sizable branching fractions ${\rm{Br}}(B_0\to J/\psi K^0f_0\to J/\psi K^0\pi^+\pi^-)=1.38\times10^{-6}$ and $ {\rm{Br}}(B_0\to J/\psi K^0a_0\to J/\psi K^0\pi^0\eta)=2.56\times10^{-7}$. This potential investigation can help us obtain the information of the scalar meson $f_0(980)$ or $a_0(980)$.


I. INTRODUCTION
The dynamics of the strong interaction are described by quantum chromodynamics (QCD), and hadron spectroscopy is a method of studying QCD, meanwhile, hadron spectroscopy is the basic theory of the strong interaction. All the time, understanding the spectrum of hadron resonances [1] and establishing a connection with the QCD, is one of the important goals of hadron physics. The conventional quark models in the low-lying hadron spectrum successfully explain that baryon is a complex of three quarks, and meson is a combination of quark and antiquark [2,3]. However, even if the model provides a great deal of data about the meson and baryon resonances [2][3][4], we also cannot rule out other more exotic components, especially considering that the QCD Lagrangian includes not only quarks but gluons as well. This leads to other configurations of color singlets, such as glueballs made purely of gluons, mixed states made of a quark and gluon excitations, and multiquarks. There is also the possibility of more quark states in a hadron, such as qqqq and qqqqq, which were already mentioned in the Ref. [5]. To understand quantitatively the QCD of quarks and gluons, over the years, numerous related experiments have been conducted in search of evidence for these exotic components in the mesonic and baryonic spectrum [6][7][8][9][10].
The triangle singularity (TS) has been discussed in Ref. [11], and Landau has systematized it in Ref. [12]. The TS was fashionable in the 1960s [13][14][15][16]. In addition to ordinary hadronic, molecular, or multiquarks states, TSs can produce peaks, but these peaks are produced by simple kinematic effects. The Coleman-Norton theorem [17] states that the Feynman amplitude has a singularity on the physical boundary as time moves forward if the decay process can be interpreted as taking place during the conservation of energy and momentum in space-time, and all internal particles really exist on the shell. In the process of particle 1 decays into particles 2 and 3, particle * xuanluo@ahu.edu.cn 1 first decays into particles A and B, then A decays into particles 2 and C, and finally particles B and C fused into an external particle 3. Particles A, B, and C are intermediate particles, and if the momenta of these intermediate particles can take on-shell values, a singularity will occur. A new way to understand this process is proposed in Ref. [18], this method does not compute the entire amplitude of the Feynman diagram including the triangle loop. The condition for producing a TS is q on+ = q a− [18], where q on is the on-shell momentum of particle A or B in the rest system of particle 1, q a− defines one of the two solutions to the momentum of particle B when B and C produce particle 3 on the shell. Since the process of the triangle mechanism involves the fusion of hadrons, the presence of hadronic molecular states plays an important role in having measurable strength. Therefore, the study of singularity is also a useful tool to study the molecular states of hadrons.
The isospin violation in production of the f 0 (980) or a 0 (980) resonances generation and its mixing have long been controversial. We have to give up the idea of trying to establish a "f 0 − a 0 mixing parameter" out of different reactions because it was shown that the isospin violation depends a lot on the reaction [19][20][21][22]. The spark was raised by the puzzle of the anomalously large isospin violation in the η(1450) → π 0 f 0 (980) decay [23], which was due to a TS and explained in Refs. [21,22]. From the point of view of the f 0 (980) and a 0 (980) themselves, the above reaction is also very enlightening. These resonances are not produced directly but are produced by the rescattering of KK, the mechanism for the formation of these resonances in the chiral unitary approach [24][25][26][27]. In addition, there are a large number of processes with the same mechanism, such as the process τ − → ν τ π − f 0 (a 0 ) [28], the process B 0 s → J/ψπ 0 f 0 (a 0 ) [29], the process D + s → π + π 0 f 0 (a 0 ) [30], the process B − → D * 0 π − f 0 (a 0 ) [31] and the same TS was shown in Refs. [32,33] to provide a plausible explanation for the peak observed in the πf 0 (980) final state, it is easy to envisage many reactions of this type [34], it also inspires us to find more processes like this type and search for TS enhanced isospin-violating reactions producing the f 0 (980) The Feynman diagrams of the decay process B0 → J/ψK 0 f0(980) → J/ψK 0 π + π − involving a triangle loop.
In the present work we study the reactions B 0 → J/ψK 0 f 0 (980) and B 0 → J/ψK 0 a 0 (980), both decays modes are allowed. The process followed by the φ decay into K 0K 0 and the K 0K 0 fuse into f 0 (a 0 ) generate a singularity, we show that it develops a TS at an invariant mass M inv (K 0 R) 1520 MeV. Meanwhile, we can obtain which show the shapes of the f 0 (980) and a 0 (980) resonances in the π + π − or π 0 η mass distributions respectively. We can restrict the integral in M inv (π + π − ) or M inv (π 0 η) to this region when calculating the mass distribution Then we integrate over the π + π − or π 0 η invariant masses and obtain which shows a clear peak around 1520 MeV. The further integration over M inv (K 0 R) provides us branching ratios for B 0 → J/ψK 0 π + π − B 0 → J/ψK 0 π 0 η, we find that the mass distribution of these decay processes showed a peak associated with TS. In addition, the corresponding decay branching ratio is obtained, and we find the branching fractions Br(B 0 → J/ψK 0 f 0 (a 0 )) = 1.007 × 10 −5 , Br(B 0 → J/ψK 0 f 0 (980) → J/ψK 0 π + π − ) = 1.38 × 10 −6 and Br(B 0 → J/ψK 0 a 0 (980) → J/ψK 0 π 0 η) = 2.56 × 10 −7 .In any case, the main aim of the present work is to point out the presence of the TS in this reaction. This work provides one more measurable example of a TS, which has been quite sparse up to now. The singularity generated by this process can also play an early warning role for future experiments.
We take Fig. 1 for example to perform the following discussion since Fig. 1 and Fig. 2 have nearly the same amplitude. Now we want to find the position of TS in complex-q plane, instead of evaluating the whole amplitude of a Feynman diagram including triangle loop, analogously to Refs. [18,36,37] we can use where the q on+ is the on-shell three momenta of the K 0 in the center of mass frame of φK 0 , meanwhile, q a− defines one of the two solutions for the momentum of φ when φK 0 are on-shell to produce f 0 and the q a− can be obtained by analyzing the singularity structure of the triangle loop. q on+ and q a− are given by we use s denotes the squared invariant mass of φ and K 0 , and λ(x, y, z) = It is easy to realize that E K 0 and p * K 0 are the energy and momentum of K 0 in the center of mass frame of the φK 0 system, v and γ are the velocity of the f 0 and Lorentz boost factor. In addition, we can easily get The three intermediate particles must be on the shell, in the φK 0 center of mass frame, we let the mass of f 0 slightly larger than the mass sum of K 0 andK 0 , now we can determine that the mass of K 0 is 497.61 MeV, because we're making the mass of f 0 a little bit greater than the sum of the masses of K 0 andK 0 , so now we're making the mass of f 0 991 MeV, we can find a TS at around √ s = 1520 MeV. If we use in Eq. 1 complex masses (M − iΓ/2) of mesons which include widths of K 0 andK 0 , the solution of Eq. 1 is then (1520−10i) MeV and the solution implies that the TS has a "width" 20 MeV.
The differential decay width over the invariant mass distribution φK 0 can be written as where we use the polarization summation formula we have Then one can obtain where the integration is performed from B. The triangle mechanism in the decay B0 → J/ψK 0 f0(980) , f0(980) → π + π − In the previous subsection, we have calculated the transition strength of the decay process B 0 → J/ψK 0 f 0 (980), now we calculate the contribution of the vertex φ → K 0K 0 , we can obtain this V P P vertex from the chiral invariant lagrangian with local hidden symmetry given in Refs. [39][40][41][42].
where the stand for the trace of the flavor SU(3) matrices and g is the coupling in the local hidden gauge , The V and P in Eq. 14 are the vector meson matrix and pseudoscalar meson matrix in the SU(3) group, P and V are given by we get where Now we can write the invariant mass distribution where p J/ψ is the momentum of the J/ψ in the B 0 rest frame , andp K 0 = | k| is the momentum of the K 0 in the K 0 f 0 (980) rest frame , where Then we have The K 0K 0 → π + π − and K 0K 0 → π 0 η scattering has been studied in detail in Refs. [43,44] within the chiral unitary approach, where altogether six channels were taken into account, including π + π − , π 0 π 0 , K + K − , K 0K 0 , ηη and π 0 η. In the present study we use this as input, and we shall see simultaneously both the f 0 (980) (with I = 0) and a 0 (980) (with I = 1) productions. Now we can write down the double differential mass distribution for the decay process B 0 → J/ψK 0 f 0 (980) → J/ψK 0 π + π − . Now for the case of f 0 (980), we only have the decay f 0 → π + π − , and thus where p J/ψ is the momentum of the J/ψ in the B 0 rest frame,p K 0 is the momentum of the K 0 in the K 0 f 0 (980) ret frame, andp π is the momentum of the π in the π + π − rest frame: Fig. 3 the contribution of the triangle loop to the total amplitude and the triangle loop defined in Eq. (17). In order to satisfy the TS condition of the Eq. (1), all intermediate particles must be on the shell so the mass sum of K 0K 0 must be less than that of f 0 (a 0 ). On the other hand, if the mass of f 0 (a 0 ) is too large, the Eq. (1) will no longer satisfy. From the perspective of the above considerations, we plot the real and imaginary parts of t T , as well as the absolute value with M inv (R) fixed at 986, 991, and 996 MeV. It can be observed that Re(t T ) has a peak around 1518 MeV, Im(t T ) has a peak around 1529 MeV, and there is a peak for |t T | around 1520 MeV. As discussed in Refs. [45,46], the peak of the imginary part is related to the TSs, while the one of the real part is related to the K * K threshold. Note that around 1520 MeV and above the TS dominates the reaction.
In the upper panel of Fig. 5 we plot Eq. 23 for the B 0 → J/ψK 0 π + π − decay, and similarly in the lower panel of Fig. 5 for the B 0 → J/ψK 0 π 0 η decay as a function of M inv (R), where in both figures we fix M inv (K 0 R)=1500 MeV, 1520 MeV, and 1540 MeV and vary M inv (R). We can see that the distribution with the largest strength is near M inv (K 0 R)=1520 MeV, we can also observe a strong peak when M inv (π + π − ) around 980 MeV in the upper panel of Fig. 5, we can also get very similar results in the lower panel of Fig. 5. We see that most of the contribution to the width Γ comes from M inv (K 0 R)=M R , and we have strong contributions for M inv (π + π − ) ∈ [500 MeV, 990 MeV] and M inv (π 0 η) ∈ [800 MeV, 990 MeV]. Therefore, when we calculate the mass distribution dΓ dMinv(K 0 R) , we restrict the integral to the limits already mentioned and perform the integration.
We show the Eq. (23) for both B 0 → J/ψK 0 π + π − and B 0 → J/ψK 0 π 0 η. When we integrate over M inv (R) we obtain dΓ dMinv(K 0 R) which we show in Fig. 6. We see a clear peak of the distribution around 1520 MeV, for f 0 and a 0 production. At the same time, we can observe that the peak of M inv (K 0 a 0 ) is significantly lower than the peak of M inv (K 0 f 0 ).

IV. CONCLUSION
We have performed the calculations for the reactions B 0 → J/ψK 0 f 0 (980)(a 0 (980)) and shown that they develop a TS for an invariant mass of 1520 MeV in (K 0 R). This TS shows up as a peak in the invariant mass distribution of these pairs with an apparent width of about 20 MeV. We have applied the experimental data of the branching ratio of the decay B 0 → J/ψK 0 f 0 (980)(a 0 (980)) to determine the coupling strength of the B 0 → J/ψK 0 f 0 (980)(a 0 (980)) vertex.
This peak is a result of the singularity of the triangle and may be misidentified with resonance when the experiment is completed. In this sense, the work done here should serve as a warning not to treat it as resonance when this peak is seen in future experiments. It is important to discover new conditions about TSs, and to allow for this possibility when experimentally observed peaks can avoid associating these peaks with resonance. The value of this work lies in identifying a TS for a suitable reaction and then preparing the results and research to interpret the peak correctly when it is observed.