Processing math: 100%

Predictions of Υ(4S)hb(1P,2P)π+π transitions

  • We study the contributions of intermediate bottomonium-like Zb states and the bottom meson loops in the heavy quark spin flip transitions Υ(4S)hb(1P,2P)π+π. Depending on the constructive or destructive interferences between the Zb-exchange and the bottom meson loops mechanisms, we predict two possible branching ratios for each process: BRΥ(4S)hb(1P)π+π(1.2+0.80.4×106) or (0.5+0.50.2×106), and BRΥ(4S)hb(2P)π+π(7.1+1.71.1×1010) or (2.4+0.20.1×1010). The contribution of the bottom meson loops is found to be considerably larger than that of the Zb-exchange in the Υ(4S)hb(1P)ππ transitions, while its decay rates are not comparable to those of heavy quark spin conserved Υ(4S)Υ(1S,2S)ππ processes. We also predict the contribution of the charm meson loops in the branch fractions of ψ(3S,4S)hc(1P)ππ.
  • Hadronic transitions Υ(mS)Υ(lS)ππ and Υ(mS)hb(nP)ππ are important processes for understanding heavy-quarkonium dynamics and low-energy quantum chromodynamics (QCD). Because bottomonia are expected to be compact and nonrelativistic, the QCD multipole expansion (QCDME) method [14] is often used to analyze these transitions, where pions are emitted because of the hadronization of soft gluons. The decay rates of Υ(2S,3S)Υ(1S,2S)ππ can be well described by QCDME [5]. Because the total spin of bˉb system in Υ(mS) and hb(nP) are 1 and 0, respectively, in general, the heavy quark spin flip Υ(mS)hb(nP)ππ processes are expected to be suppressed, compared with the heavy quark spin conserved Υ(mS)Υ(nS)ππ processes. Within the framework of QCDME, studies [57] predicted that the branching fraction of Υ(3S)Υ(1P)ππ is suppressed by two orders of magnitude, relative to that of Υ(3S)Υ(1S)ππ, while Ref. [8] predicted a suppression of at least three orders of magnitude. The prediction of Ref. [8] is supported by experimental data [9]. In the decay processes Υ(5S)Υ(lS)π+π (l=1,2,3) and Υ(5S)hb(nP)π+π (n=1,2) where the two charged bottomoniumlike resonances Zb(10610)± and Zb(10650)± were observed, Υ(5S)hb(nP)π+π proceeds at a rate comparable to the Υ(5S)Υ(lS)π+π processes [10, 11]. The mechanism that mitigates this expected suppression has remained controversial. In Refs. [12, 13], the Υ(5S)hb(nP)π+π processes were interpreted via bottom meson loops mechanism, while genuine S-matrix Zb poles are required as in Refs. [1417]. The meson loops mechanism has been investigated by many previous works [1824] to study the dipion and η transitions of higher charmonia and bottomonia because the branch ratios and dipion invariant mass spectra cannot be described by QCDME.

    In this work, we will study whether the bottom meson loops mechanism can produce the Υ(4S)hb(nP)π+π transitions at decay ratios comparable to Υ(4S)Υ(lS)π+π. Because in the dipion emission processes of the Υ(4S) the crossed-channel exchanged Zb cannot be on-shell, these transitions are expected to be good channels to study the bottom meson loops' effect. In our previous works [25, 26], by using the nonrelativistic effective field theory (NREFT), we calculated the effects of the bottom meson loops as well as the Zb-exchange in the Υ(4S)Υ(1S,2S)ππ processes, and we found that the experimental data can be described well. Here, within the same theoretical scheme, we will calculate the contributions of the bottom meson loops and Zb-exchange in the Υ(4S)hb(nP)π+π processes, and give theoretical predictions for the decay branching ratios. We find that the contribution of the bottom meson loops is considerably larger than that of the Zb-exchange in the Υ(4S)hb(1P)π+π process, and it cannot produce a rate comparable to that of Υ(4S)Υ(1S,2S)π+π.

    The remainder of this paper is organized as follows. In Sec. 2, the theoretical framework is described in detail. In Sec. 3, we provide the theoretical predictions for the decay branching fractions of Υ(4S)hb(1P,2P)π+π, and discuss the contributions of different mechanisms. The study is concluded in Sec. 4.

    To calculate the contribution of the mechanism Υ(mS)Zbπhb(nP)ππ, we need the effective Lagrangians for the ZbΥπ and Zbhbπ interactions [27],

    LZbΥπ=j=1,2CZbjΥ(mS)πΥi(mS)Zibjuμvμ+h.c.,

    (1)

    LZbhbπ=j=1,2gZbjhb(nP)πϵijkZibjujhkb+h.c.,

    (2)

    where Zb1 and Zb2 denote Zb(10610) and Zb(10650), respectively, and vμ=(1,0) is the velocity of the heavy quark. The Zb states are given in the matrix as

    Zibj=(12Z0ibjZ+ibjZibj12Z0ibj).

    (3)

    The pions can be parametrized as Goldstone bosons of the spontaneous breaking of the chiral symmetry:

    uμ=i(uμuuμu),u=exp(iΦ2Fπ),Φ=(12π0π+π12π0),

    (4)

    where Fπ=92.2MeV is the pion decay constant.

    To calculate the box diagrams, we need the Lagrangian for the coupling of the Υ to the bottom mesons and the coupling of the hb to the bottom mesons [27, 28],

    LΥHH=igJHH2JHaσˉHa+h.c.,

    (5)

    LhbHH=ig12hibHaσiˉHa+h.c.,

    (6)

    where JΥσ+ηb denotes the heavy quarkonia spin multiplet, Ha=Vaσ+Pa with Pa(Va)=(B(),ˉB()0) collects the bottom mesons, and ABA(B)(A)B. We also need the Lagrangian for the axial coupling of the pion fields to the bottom and antibottom mesons, which at the lowest order in heavy-flavor chiral perturbation theory is given by [2933]

    LHHΦ=gπ2ˉHaσuabˉHbgπ2HaHbσuba,

    (7)

    where ui=2iΦ/F+O(Φ3) denotes the three-vector components of uμ, as defined in Eq. (4). Here, we use gπ=0.492±0.029 from a recent lattice QCD calculation [34].

    Because the Υ(4S) meson is above the BˉB threshold and decays predominantly into BˉB pairs, the loop mechanism with intermediate bottom mesons may be important in the transitions Υ(4S)hb(nP)π+π. By following the formalism set-up based on NREFT [28, 35, 36], we will analyze the power counting of different types of loops. In NREFT, the expansion parameter is the velocity of the intermediate heavy meson, namely νX=|mXmB()mB()|/mB(), which is small because the bottomonia X are close to the B()ˉB() thresholds. In this power counting, each nonrelativistic propagator scale as 1/ν2, and the measure of one-loop integration scales as d4lν5.

    There are five different kinds of loop contributions, namely the box diagrams displayed in Fig. 1(b), (c), triangle diagrams displayed in Fig. 2(a)(c), and the bubble loop in Fig. 2(d). We analyze them one by one as follows:

    Figure 1

    Figure 1.  Feynman diagrams considered for Υ(mS)hb(nP)ππ processes. Crossed diagrams of (a) and (b) are not shown explicitly.

    Figure 2

    Figure 2.  Loop diagrams not considered in the calculations. The corresponding power counting arguments are given in the main text.

    First, we analyze the power counting of the box diagrams, namely Fig. 1(b), (c). As shown in Eq. (7), the vertex of B()B()π is proportional to the external momentum of the pion qπ. The ΥB()ˉB() vertex is in a P-wave, and the hbB()ˉB() vertex is in an S-wave; therefore, the loop momentum must contract with the external pion momentum, and hence the P-wave vertex scales as O(qπ). Thus, the box diagrams scales as ν5q3π/ν8=q3π/ν3.

    For the triangle diagram Fig. 2(a), the leading ΥB()ˉB()π vertex given by gJHHπJˉHaHbu0ab [37] is proportional to the energy of the pion, Eπqπ. Therefore, Fig. 2(a) is counted as mBν5q2π/ν6=mBq2π/ν, where the factor mB is introduced to match the dimension with the scaling for the box diagrams.

    In Fig. 2(b), the leading hbB()ˉB()π vertex given by ghbHHπhibHaσjˉHbϵijkukab [29] is proportional to the momentum of the pion qπ. The loop momentum due to the ΥB()ˉB() coupling has to contract with the external pion momentum. Thus, Fig. 2(b) scales as ν5q3π/ν6=q3π/ν.

    The leading B()B()ππ vertex comes from the chiral derivative term Ha(iD0)baHb=Ha(i0iV0)baHb [38, 39], in which the pion pair produced by the vector current, Vμ=12(uμu+uμu), cannot form a positive-parity and C-parity state; therefore, this leading vertex does not contribute to the Υ(mS)hb(nP)ππ processes. Isoscalar, PC=++ pion pairs only enter in the next order O(q2π) from point vertices. Therefore, Fig. 2(c) scales as ν5q3π/ν6=q3π/ν.

    In Fig. 2(d), both the initial and final vertices are proportional to qπ; therefore, the bubble loop scales as mBν5q2π/ν4=mBq2πν.

    Therefore, we expect that the ratios of the contributions of the box diagrams, triangle diagram Fig. 2(a)(c), and the bubble loop Fig. 2(d) are

    q3πν3:mBq2πν:q3πν:q3πν:mBq2πν=1:mBν2qπ:ν2:ν2:mBν4qπ,

    (8)

    where qπ(mΥ(4S)mhb(nP))/2 and ν=(νΥ(4S)+νhb(nP))/2, with νΥ(4S)0.06, νhb(1P)0.35, and νhb(2P)0.24. Thus, for the Υ(4S)hb(1P)π+π transition, the ratios in Eq. (8) are 1:0.67:0.04:0.04:0.03. For the Υ(4S)hb(2P)π+π transition, the ratios are 1:0.75:0.02:0.02:0.02. Therefore, according to the power counting, the box diagrams and triangle diagram in Fig. 2(a) are dominant among the loop contributions, and they are of the same order. The Υ(4S) is below the B()ˉB()π threshold, and the coupling gJHHπ in the triangle diagram Fig. 2(a) is unknown. Thus, for a rough estimation of the loop contributions, we will only calculate the box diagrams in the present study. All the box and triangle loop contributions discussed here are ultraviolet-finite and do not require additional introduction of counterterms.

    The decay amplitude for

    Υ(mS)(pa)hb(nP)(pb)π(pc)π(pd)

    (9)

    is described in terms of the Mandelstam variables

    s=(pc+pd)2,t=(papc)2,u=(papd)2.

    (10)

    By using the effective Lagrangians in Eqs. (1) and (2), the tree amplitude of Υ(mS)Zbπhb(nP)ππ can be obtained as

    MZb=2mΥ(mS)mhb(nP)F2πϵabjϵaΥ(mS)ϵbhb(nP)i=1,2mZbiCZbiΥ(mS)πgZbihb(nP)π{p0cpjd1tm2Zbi+p0dpjc1um2Zbi}.

    (11)

    The nonrelativistic normalization factor mY has been multiplied with the amplitude for every heavy particle, with Y=Υ(mS),hb(nP),Zbi. The widths of the Zb states are neglected in the present study, because they are of the order of 10MeV and are considerably smaller than the difference between the Zb masses and the Υ(mS)π/hb(nP)π threshold.

    Now, we discuss the calculation of the box diagrams. In the box diagrams Fig. 1(b) and (c), we denote the top left intermediate bottom meson as M1, and the other intermediate bottom mesons as M2, M3, and M4, in counterclockwise order. For the pseudoscalar or vector content of [M1,M2,M3,M4], there are twelve possible patterns and we number them in order: 1, [PPPV]; 2, [PPVV]; 3, [PVPV]; 4, [PVVP]; 5, [VVPP]; 6, [VPVP]; 7, [VPPV]; 8, [PVVV]; 9, [VPVV]; 10, [VVPV]; 11, [VVVP]; 12, [VVVV]. For each pattern, we also need to consider six possibilities of different flavor of the intermediate bottom mesons: [B()+,B(),B()+,B()0], [B()+,B(),ˉB()0,B()0], [B()0,ˉB()0,B()0,B()+], [B(),B()+,B(),ˉB()0], [ˉB()0,B()0,B()+,B()], and [ˉB()0,B()0,ˉB()0,B()]. The full amplitude contains the sum of all possible amplitudes.

    For the tensor reduction of the loop integrals, it is convenient to define q=pb and thperpendicular momentum q=pcq(qpc)/q2, which satisfy qq=0. The result of the amplitude of the box diagrams can be written as

    Mloop=ϵaΥ(mS)ϵbhb(nP){ϵabiqiA1+ϵabiqiA2+ϵbijqiqjqaA3+ϵbijqiqjqaA4+ϵaijqiqjqbA5+ϵaijqiqjqbA6}.

    (12)

    Details on the analytic calculation of the box diagrams and explicit expressions of Ai (i=1,2,...,6) are given in Appendix A.

    The decay width for Υ(mS)hb(nP)ππ is given by

    Γ=s+st+t|MZb+Mloop|2dsdt768π3m3Υ(mS),

    (13)

    where the lower and upper limits are given as

    s=4m2π,s+=(mΥ(mS)mhb(nP))2,t±=14s{(m2Υ(mS)m2hb(nP))2[λ12(s,m2π,m2π)λ12(m2Υ(mS),s,m2hb(nP))]2},λ(a,b,c)=a2+b2+c22(ab+ac+bc).

    (14)

    To estimate the contribution of the Zb-exchange mechanism, we need to know the coupling strengths of ZbΥ(4S)π and Zbhb(nP)π. The mass difference between Zb(10610) and Zb(10650) is considerably smaller than the difference between their masses and the Υ(mS)π/hb(nP)π threshold, and they have the same quantum numbers, and thus the same coupling structures as given by Eqs. (1) and (2). Therefore, it is very difficult to distinguish their effects from each other in the dipion transitions of Υ(4S), so we only used one Zb, the Zb(10610), which approximately combines both Zb states' effects. In Ref. [26], we studied the Υ(4S)Υ(mS)ππ processes to extract the coupling constant |CZbΥ(4S)π|=(3.3±0.1)×103, which combine the effects from both Zb states. For the couplings of Zbhb(nP)π, in principle, they can be extracted from the partial widths of the Zb states decay into hb(nP)π(n=1,2)

    |gZbhbπ|={6πF2πmZbΓZbhbπ|pf|3mhb}12,

    (15)

    where |pf|λ1/2(m2Zb,m2hb,m2π)/(2mZb). The branching fractions of the decays of both Zb states into hb(nP)π(n=1,2) has been given in [40], where the Zb line shapes were described using Breit-Wigner forms. If we naively use these branching fractions, we would obtain

    |gnaiveZb1hb(1P)π|=0.019±0.003,|gnaiveZb2hb(1P)π|=0.021±0.003,|gnaiveZb1hb(2P)π|=0.068±0.011,|gnaiveZb2hb(2P)π|=0.077±0.010.

    (16)

    Here, all the Zbhbπ couplings are labeled by a superscript "naive" because this is not the appropriate way to extract coupling strengths in this case; the Zb states are very close to the B()ˉB thresholds. Therefore, the Flatté parametrization for the Zb spectral functions should be used, which will lead to much larger partial widths into (bˉbπ) channels, and thus the relevant coupling strengths. As analyzed in Ref. [41], in the the Flatté parametrization the sum of the partial widths of the Zb(10610) other than that for the BˉB channel should be greater than the nominal width, which is approximately 20MeV. While summing over all the Υ(nS)π(n=1,2,3) and hb(mP)π(m=1,2), and the branching fractions in Ref. [40] is approximately 14% or 3MeV in terms of partial widths. Therefore, for a rough estimation, we use three times the results of Eq. (16):

    |gZbhb(1P)π|0.057,|gZbhb(2P)π|0.204.

    (17)

    We find that even after considering the enlarging factor of three for the couplings |gZbhb(nP)π|, the Zb-exchange contribution is still considerably smaller than the bottom meson loops contributions.

    In the calculation of the box diagrams, the coupling strength gJHH(4S) can be extracted from the measured open-bottom decay widths of the Υ(4S), and we have gJHH(4S)=1.43±0.01GeV3/2. For the hbBˉB coupling g1, we can use the results from Ref. [27]. In [27], the Zb-exchange mechanism in the Υ(5S)hb(1P,2P)ππ processes was studied assuming the Zb states are B()ˉB bound states and the physical coupling of the Zb states to the bottom and anti-bottom mesons, z1, as well as the product g1z1 was determined. By using their results z1=0.75+0.080.11GeV1/2 and g1z1=0.40±0.06GeV1, we can extract that g1=0.53+0.190.13GeV1/2. In [27], in order to reduce the number of free parameters, the couplings of hb(1P)BˉB and hb(2P)BˉB are assumed to be the same.

    By using the coupling strengths above, we can predict the decay branching fractions of Υ(4S) hb(1P,2P)π+π. Depending on the sign of the couplings in Eq. (17), the interferences can be constructive or destructive between the Zb-exchange and box graph mechanisms; therefore, there are two possible results for each process

    BRΥ(4S)hb(1P)π+π(1.2+0.80.4×106)or(0.5+0.50.2×106),BRΥ(4S)hb(2P)π+π(7.1+1.71.1×1010)or(2.4+0.20.1×1010).

    (18)

    We find that the BRΥ(4S)hb(1P)π+π is at least one order of magnitude smaller than the branching fractions BRΥ(4S)Υ(1S,2S)π+π, which are approximately 8×105 given in PDG [42], and the BRΥ(4S)hb(2P)π+π is small owing to the very small phase space. We discuss the Υ(4S)hb(1P)π+π transition in further detail. To illustrate the effects of the Zb-exchange and box graph mechanisms in Υ(4S)hb(1P)ππ, we give the predictions only including the Zb-exchange terms or only including the box diagrams

    BRZbΥ(4S)hb(1P)π+π=0.6+0.10.1×107,BRBoxΥ(4S)hb(1P)π+π=0.8+0.70.3×106.

    (19)

    We can observe that the bottom meson loops contribution is considerably larger than the Zb-exchange contribution, while it is two orders of magnitude smaller than the Υ(4S)Υ(1S,2S)π+π transitions. Note that the direct gluon hadronization mechanism contribution within QCDME for the Υ(4S)hb(1P)π+π process has not been calculated thus far. In references [58], QCDME predict that the branching fraction of Υ(3S)hb(1P)ππ is 2-3 orders of magnitude suppressed compared to of Υ(3S)Υ(1S)ππ. The three orders of magnitude suppression is supported by experiment [9]. Because the mass difference between the Υ(4S) and hb(1P) is approximately 0.68 GeV, the pions in the Υ(4S)hb(1P)π+π process can also be considered to be in the soft region. If one approximates that the gluon hadronization mechanism within QCDME in Υ(4S)hb(1P)π+π is also 2-3 orders of magnitudes suppressed compared to that in the Υ(4S)Υ(1S)π+π process, as in the Υ(3S) decay cases. Then, the gluon hadronization mechanism contribution is at most at the same order of the bottom meson loops contribution. Owing to the lack of exact information on the gluon hadronization within QCDME and the neglecting of the triangle diagram Fig. 2(a) as discussed in section 2.3, it is important to note that the results presented in this paper are order-of-magnitude estimates.

    In Fig. 3, we plot the distributions of the ππ and hbπ invariant mass spectra, and the distribution of cosθ, where θ is defined as the angle between the initial Υ(mS) and the π+ in the rest frame of the ππ system. To illustrate the effects of different mechanisms, the contributions of the box diagrams, Zb-exchange, their sum with the constructive interference, and the sum with destructive interference are indicated by dark green dashed, magenta dotted, red solid, and blue dot-dashed lines, respectively. There is a broad bump at approximately 0.5 GeV in the dipion invariant mass distribution. The ππ invariant mass spectra with unknown normalization predicted within QCDME in Ref. [5] showed a peak at low ππ masses. Thus, the ππ invariant mass spectra can be useful to identify the effects of the bottom meson loops and the gluon hadronization mechanism with future experimental data. Further, the angular distribution is far from flat. In the Υ(mS)hb(nP)ππ process, the isospin conservation combined with Bose symmetry requires the pions to have even relative angular momentum. Therefore, there is a large D-wave component from the box diagrams if higher partial waves are neglected.

    Figure 3

    Figure 3.  (color online) Theoretical predictions of the distributions of the ππ and hbπ invariant mass spectra, and helicity angular distributions in the Υ(4S)hb(1P)ππ process. Dark green dashed, magenta dotted, red solid, and blue dot-dashed lines represent the contributions of the box diagrams, Zb-exchange, their sum with constructive interference, and the sum with destructive interference, respectively.

    The Υ(mS)hb(nP)ππ are heavy quark spin flip processes, and they are forbidden in the heavy quark limit. We checked that in the heavy quark limit, i.e. mB=mB, all the box diagrams were cancelled with each other; therefore, the bottomed loops did not contribute to the Υ(mS)hb(nP)ππ transitions. With the small mass splitting of B and B in the real world, as shown in Eqs. (18) and (19), the bottomed meson loops contribution does not produce Υ(4S)hb(1P)ππ at a rate comparable to the heavy quark spin conserved Υ(4S)Υ(1S,2S)ππ transitions. Note that the datasets collected at Υ(4S) by BABAR and Belle II collaborations are 471×106 and 772×106 [43], respectively. Thus, they should contain several hundreds of Υ(4S)hb(1P)ππ events according to our calculation. We hope future experimental analysis by BABAR and Belle can test our predictions. As stated in the introduction, the observed Υ(5S)hb(nP)π+π proceed at a rate comparable to the Υ(5S)Υ(lS)π+π processes [10, 11]. The enhancements may be caused by the effects of the on-shell Zb exchange and the two-cut condition complexity of the bottom meson loops in the Υ(5S) decays. A detailed analysis of the Υ(5S) hb(nP)π+π processes is beyond the scope of this study.

    Owing to the similarity between the bottomonium and charmonium families, we can extend the box diagrams calculation to give a rough estimation of the branch fractions of the ψ(3S)(ψ(4040))hc(1P)π+π and ψ(4S)(ψ(4415))hc(1P)π+π transitions. The relevant Feynman diagrams can be obtained by replacing the external Υ(mS) and hb(nP) by ψ(mS) and hc(nP), respectively, and replacing the intermediate B() by D() in Fig. 1(b) and (c). The experimental decay widths of the ψ(3S,4S)D()ˉD() transitions are not given in PDG, and we will use the theoretical predictions of the decay widths in Ref. [44] to estimate the coupling strengths gJHH(ψ(3S)) and gJHH(ψ(4S)). Because among the different decay modes D()ˉD(), the DˉD and DˉD modes are dominant for ψ(3S) and ψ(4S), respectively, we will use the corresponding coupling constants in the calculation, namely gJHH(ψ(3S))=gψ(3S)DˉD=0.97GeV3/2 and gJHH(ψ(4S))=gψ(4S)DˉD=0.25GeV3/2. For the hcDˉD coupling, we use the result from Ref. [45], ghcDˉD=(mχc0/3)/(fχc0)=(3.415/3)/(0.297)GeV1/2=3.59GeV1/2. The predictions of the box diagrams contributions to the branch fractions of ψ(3S,4S)hc(1P)π+π are

    BRBoxψ(3S)hc(1P)π+π=2.9×105,BRBoxψ(4S)hc(1P)π+π=4.5×103.

    (20)

    The prediction of BRBoxψ(3S)hc(1P)π+π is below the upper limit given in PDG [42]. As expected the branch fractions BRBoxψ(3S,4S)hc(1P)π+π are considerably larger than BRBoxΥ(4S)hb(1P)π+π, because the mass splitting of D and D is considerably larger than that of B and B. Note that this is just a preliminary rough estimation, owing to the lack of sufficient information concerning the ψ(3S,4S)D()ˉD() coupling constants and the neglecting of the loop diagrams with intermediate D1 state in the present calculation. A detailed theoretical study of the ψ(3S,4S) hc(1P)π+π transitions will be pursued in the future.

    In this study, we investigate the effects of Zb exchange and bottom meson loops in the heavy quark spin flip transitions Υ(4S)hb(nP)ππ(n=1,2). The bottom meson loops are treated in the NREFT scheme, in which the dominant box diagrams are considered. We find that the bottom meson loops contribution is considerably larger than the Zb-exchange contribution in the Υ(4S) hb(1P)ππ transition, while it can not produce decay rates comparable to the heavy quark spin conserved Υ(4S)Υ(1S,2S)ππ processes. The theoretical prediction of the decay rate and the dipion invariant mass spectra of Υ(4S)hb(1P)ππ in this work may be useful for identifying the effect of the bottom meson loops with future experimental analysis. We also predict the branch fractions of ψ(3S,4S)hc(1P)ππ contributed from the charm meson loops.

    In this appendix, first, we discuss the parametrization and simplification of the scalar four-point integrals in the box diagrams. Then, we introduce a tensor reduction scheme to deal with higher-rank loop integrals. Finally, we give the amplitude of the box diagrams for the Υ(mS)hb(nP)ππ process.

    For the first topology, as shown in Fig. A1, the scalar integral evaluated for the initial bottomonium at rest (p=(M,0)) reads

    Figure A1

    Figure A1.  Kinematics used for calculating four-point integrals.

    We are grateful to the referees' useful suggestions and constructive remarks, which helped formulate the present version of this manuscript. We are grateful to Martin Cleven for the collaboration in the early stages of this study. We acknowledge Guo-Ying Chen, Meng-Lin Du, and Qian Wang for their helpful discussions, and Feng-Kun Guo for a careful reading of the manuscript and valuable comments.

    J(0)1id4l(2π)41[l2m21+iϵ][(pl)2m22+iϵ][(lq1q2)2m23+iϵ][(lq1)2m24+iϵ]i16m1m2m3m4×d4l(2π)41[l0l22m1m1+iϵ][l0M+l22m2+m2iϵ]1[l0q01q02(l+q)22m3m3+iϵ][l0q01(lq1)22m4m4+iϵ].

    By performing the contour integration, we find

    μ12μ23μ242m1m2m3m4d3l(2π)31[l2+c12iϵ][l2+2μ23m3lq+c23iϵ][l22μ24m4lq1+c24iϵ],

    where we defined

    c122μ12(m1+m2M),c232μ23(m2+m3M+q01+q02+q22m3),c242μ24(m2+m4M+q01+q212m4),μij=mimjmi+mj.

    The second topology in Fig. A1 is just the crossed diagram of the first topology with q1q2, so the scalar integral reads

    J(0)2=μ12μ23μ242m1m2m3m4d3l(2π)3×1[l2+c12iϵ][l2+2μ23m3lq+c23iϵ][l22μ24m4lq2+c24iϵ],

    where

    c242μ24(m2+m4M+q02+q222m4).

    For the third topology, we have

    J(0)3id4l(2π)41[l2m21+iϵ][(pl)2m22+iϵ][(pq2l)2m23+iϵ][(lq1)2m24+iϵ]i16m1m2m3m4d4l(2π)41[l0l22m1m1+iϵ][l0M+l22m2+m2iϵ]×1[l0+q02M+(l+q2)22m3+m3iϵ][l0q01(lq1)22m4m4+iϵ].

    By the contour integration, we find

    μ12μ342m1m2m3m4d3l(2π)31[l2+d12iϵ][l22μ34m4lq12μ34m3lq2+d34iϵ]×[μ24[l22μ24m4lq1+d24iϵ]+μ13[l2+2μ13m3lq2+d13iϵ]],

    where we defined

    d122μ12(m1+m2M),d342μ34(m3+m4q0+q212m4+q222m3),d242μ24(m2+m4M+q01+q212m4),d132μ13(m1+m3M+q02+q222m3).

    In all the three cases, the remaining three-dimensional momentum integration will be carried out numerically.

    Because the ΥB()ˉB() vertex scales with the momentum of the bottom meson pair, for topology I we have to deal with

    μ12μ23μ242m1m2m3m4d3l(2π)3f(l)[l2+c12iϵ][l2+2μ23m3lq+c23iϵ][l22μ24m4lq1+c24iϵ],

    where f(l)={1,li} for the fundamental scalar and vector integrals, respectively. A convenient parametrization of the tensor reduction is

    J(1)i1=μ12μ23μ242m1m2m3m4d3l(2π)3li[l2+c1iϵ][l22μ23m3lq+c2iϵ][l22μ24m4lq1+c3iϵ]qiJ(1)1+qi1J(2)1,

    where q1=q1q(qq1)/q2. The expressions of the scalar integrals J(r)1 can easily be disentangled and should be evaluated numerically. The corresponding expressions for topology II and III can be obtained by changing the denominators accordingly.

    We define the scalar integrals J1(i,r,k) based on the J(r)1 in the tensor reduction of vector integral in Eq. (A10), where i=1,2,3 denotes the three topologies of the box diagrams shown in Fig. A1, r=1,2 refers to the two components J(r)1, and k=1,2,...,12 represents the twelve patterns with different pseudoscalar or vector content of the intermediate bottom mesons in [M1,M2,M3,M4] as displayed in Sec. 2.3.

    We give the amplitude of the box diagrams for the Υ(mS)hb(nP)ππ process, namely the Al(l=1,2,...,6) in the Eq. (12).

    A1=8g1gJHHg2πF2πq2{q2{pcpd[J1(1,1,3)+J1(2,1,3)+J1(3,1,8)]+pcq[J1(1,1,9)+J1(1,1,11)J1(2,1,12)+J1(3,1,9)+J1(3,1,11)]+q2[J1(1,2,9)+J1(1,2,11)J1(2,2,12)+J1(3,2,9)+J1(3,2,11)]}+pcq{pcq[J1(1,1,9)+J1(1,1,11)J1(2,1,12)+J1(3,1,9)+J1(3,1,11)]+pdq[J1(1,1,12)J1(2,1,9)J1(2,1,11)+J1(3,1,10)]+q2[J1(1,2,9)+J1(1,2,11)J1(1,2,12)+J1(2,2,9)+J1(2,2,11)J1(2,2,12)+J1(3,2,9)J1(3,2,10)+J1(3,2,11)]}},

    A2=8g1gJHHg2πF2π{pcpd[J1(1,2,3)+J1(2,2,3)+J1(3,2,8)]+pcq[J1(1,1,9)+J1(1,1,11)J1(2,1,12)+J1(3,1,9)+J1(3,1,11)]+pdq[J1(1,1,12)J1(2,1,9)J1(2,1,11)+J1(3,1,10)]+q2[J1(1,2,9)+J1(1,2,11)J1(1,2,12)+J1(2,2,9)+J1(2,2,11)J1(2,2,12)+J1(3,2,9)J1(3,2,10)+J1(3,2,11)]},

    A3=8g1gJHHg2πF2πq2{q2[J1(1,1,9)+J1(1,1,11)J1(1,1,12)J1(1,2,2)+J1(1,2,9)+J1(1,2,10)+J1(1,2,12)J1(2,1,9)+J1(2,1,11)J1(2,1,12)+J1(2,2,2)J1(2,2,10)J1(2,2,11)J1(3,1,9)J1(3,1,10)+J1(3,1,11)J1(3,2,2)+J1(3,2,9)+J1(3,2,10)+J1(3,2,12)]+pcq[J1(1,2,9)J1(1,2,11)+J1(1,2,12)+J1(2,2,9)J1(2,2,11)+J1(2,2,12)+J1(3,2,9)+J1(3,2,10)J1(3,2,11)]},

    A4=8g1gJHHg2πF2πq4{q4[J1(1,1,2)J1(1,1,10)J1(1,1,11)J1(2,1,2)+J1(2,1,9)+J1(2,1,10)+J1(2,1,12)+J1(3,1,2)J1(3,1,11)J1(3,1,12)]+q2pcq[J1(1,1,9)J1(1,1,11)+J1(1,1,12)J1(1,2,9)+J1(1,2,11)J1(1,2,12)+J1(2,1,9)J1(2,1,11)+J1(2,1,12)J1(2,2,9)+J1(2,2,11)J1(2,2,12)+J1(3,1,9)+J1(3,1,10)J1(3,1,11)J1(3,2,9)J1(3,2,10)+J1(3,2,11))(pcq)2[J1(1,2,9)J1(1,2,11)+J1(1,2,12)+J1(2,2,9)J1(2,2,11)+J1(2,2,12)+J1(3,2,9)+J1(3,2,10)J1(3,2,11)]},

    A5=8g1gJHHg2πF2πq2{q2[J1(1,1,6)+J1(1,2,8)+J1(1,2,10)J1(2,1,6)+J1(2,2,6)J1(2,2,8)J1(2,2,10)J1(3,1,3)+J1(3,1,6)J1(3,1,8)+J1(3,2,3)+J1(3,2,12)]+pcq[J1(1,2,6)+J1(2,2,6)+J1(3,2,3)J1(3,2,6)+J1(3,2,8)]},

    A6=8g1gJHHg2πF2πq4{q4[J1(1,1,6)J1(1,1,8)J1(1,1,10)+J1(2,1,8)+J1(2,1,10)J1(3,1,6)+J1(3,1,8)J1(3,1,12)]+q2pcq[J1(1,1,6)J1(1,2,6)+J1(2,1,6)J1(2,2,6)+J1(3,1,3)J1(3,1,6)+J1(3,1,8)J1(3,2,3)+J1(3,2,6)J1(3,2,8)](pcq)2[J1(1,2,6)+J1(2,2,6)+J1(3,2,3)J1(3,2,6)+J1(3,2,8)]}.

    [1] M. B. Voloshin and V. I. Zakharov, Phys. Rev. Lett., 45: 688 (1980) doi: 10.1103/PhysRevLett.45.688
    [2] V. A. Novikov and M. A. Shifman, Z. Phys. C, 8: 43 (1981) doi: 10.1007/BF01429829
    [3] Y. P. Kuang and T. M. Yan, Phys. Rev. D, 24: 2874 (1981)
    [4] Y. P. Kuang, Front. Phys. China, 1: 19 (2006) doi: 10.1007/s11467-005-0012-6
    [5] Y.-P. Kuang, S. F. Tuan, and T. M. Yan, 37: 1210 (1988)
    [6] S. F. Tuan, A7, 3527 (1992); 42: 3207 (1990)
    [7] Y.-P. Kuang and T. M. Yan, 24: 2874 (1981); 41: 155 (1990)
    [8] M. B. Voloshin, 43: 1011(1986)
    [9] J. P. Lees et al (BaBar Collaboration), Phys. Rev. D, 84: 011104 (2011) doi: 10.1103/PhysRevD.84.011104
    [10] I. Adachi (Belle Collaboration), arXiv: 1105.4583[hep-ex]
    [11] A. Bondar et al (Belle Collaboration), Phys. Rev. Lett., 108: 122001 (2012) doi: 10.1103/PhysRevLett.108.122001
    [12] D. Y. Chen and X. Liu, Phys. Rev. D, 84: 094003 (2011) doi: 10.1103/PhysRevD.84.094003
    [13] E. S. Swanson, Phys. Rev. D, 91(3): 034009 (2015) doi: 10.1103/PhysRevD.91.034009
    [14] F. K. Guo, C. Hanhart, Q. Wang et al, Phys. Rev. D, 91(5): 051504 (2015) doi: 10.1103/PhysRevD.91.051504
    [15] W. S. Huo and G. Y. Chen, Eur. Phys. J. C, 76(3): 172 (2016) doi: 10.1140/epjc/s10052-016-4013-0
    [16] F.-K. Guo, C. Hanhart, Y. S. Kalashnikova et al, Phys. Rev. D, 93(7): 074031 (2016) doi: 10.1103/PhysRevD.93.074031
    [17] Q. Wang, V. Baru, A. A. Filin et al, Phys. Rev. D, 98(7): 074023 (2018) doi: 10.1103/PhysRevD.98.074023
    [18] D. Y. Chen, J. He, X. Q. Li et al, Phys. Rev. D, 84: 074006 (2011) doi: 10.1103/PhysRevD.84.074006
    [19] D. Y. Chen, X. Liu and S. L. Zhu, Phys. Rev. D, 84: 074016 (2011) doi: 10.1103/PhysRevD.84.074016
    [20] C. Meng and K. T. Chao, Phys. Rev. D, 77: 074003 (2008) doi: 10.1103/PhysRevD.77.074003
    [21] C. Meng and K. T. Chao, Phys. Rev. D, 78: 074001 (2008) doi: 10.1103/PhysRevD.78.074001
    [22] B. Wang, X. Liu, and D. Y. Chen, Phys. Rev. D, 94(9): 094039 (2016) doi: 10.1103/PhysRevD.94.094039
    [23] Y. A. Simonov and A. I. Veselov, Phys. Rev. D, 79: 034024 (2009) doi: 10.1103/PhysRevD.79.034024
    [24] B. Wang, H. Xu, X. Liu et al, Front. Phys. (Beijing), 11: 111402 (2016) doi: 10.1007/s11467-016-0564-7
    [25] Y. H. Chen, M. Cleven, J. T. Daub et al, Phys. Rev. D, 95(3): 034022 (2017) doi: 10.1103/PhysRevD.95.034022
    [26] Y. H. Chen and F. K. Guo, Phys. Rev. D, 100(5): 054035 (2019) doi: 10.1103/PhysRevD.100.054035
    [27] M. Cleven, F.-K. Guo, C. Hanhart et al, Eur. Phys. J. A, 47: 120 (2011) doi: 10.1140/epja/i2011-11120-6
    [28] F.-K. Guo, C. Hanhart, and U.-G. Meißner, Phys. Rev. Lett., 103: 082003 (2009); 104: 109901(E) (2010)
    [29] S. Fleming and T. Mehen, Phys. Rev. D, 78: 094019 (2008) doi: 10.1103/PhysRevD.78.094019
    [30] G. Burdman and J. F. Donoghue, Phys. Lett. B, 280: 287 (1992) doi: 10.1016/0370-2693(92)90068-F
    [31] M. B. Wise, Phys. Rev. D, 45: R2188 (1992) doi: 10.1103/PhysRevD.45.R2188
    [32] T. M. Yan, H. Y. Cheng, C. Y. Cheung et al, Phys. Rev. D, 46: 1148 (1992); 55: 5851(E) (1997)
    [33] R. Casalbuoni, A. Deandrea, N. Di Bartolomeo et al, Phys. Rept., 281: 145 (1997) doi: 10.1016/S0370-1573(96)00027-0
    [34] F. Bernardoni et al (ALPHA Collaboration), Phys. Lett. B, 740: 278 (2015) doi: 10.1016/j.physletb.2014.11.051
    [35] F.-K. Guo, C. Hanhart, G. Li et al, Phys. Rev. D, 83: 034013 (2011) doi: 10.1103/PhysRevD.83.034013
    [36] M. Cleven, Q. Wang, F.-K. Guo et al, Phys. Rev. D, 87: 074006 (2013) doi: 10.1103/PhysRevD.87.074006
    [37] T. Mehen and J. W. Powell, Phys. Rev. D, 88: 034017 (2013) doi: 10.1103/PhysRevD.88.034017
    [38] I. W. Stewart, Nucl. Phys. B, 529: 62 (1998) doi: 10.1016/S0550-3213(98)00374-5
    [39] J. Hu and T. Mehen, Phys. Rev. D, 73: 054003 (2006) doi: 10.1103/PhysRevD.73.054003
    [40] A. Garmash et al (Belle Collaboration), Phys. Rev. Lett., 116(21): 212001 (2016) doi: 10.1103/PhysRevLett.116.212001
    [41] Y.-H. Chen, J. T. Daub, F.-K. Guo et al, Phys. Rev. D, 93: 034030 (2016) doi: 10.1103/PhysRevD.93.034030
    [42] M. Tanabashi et al (Particle Data Group Collaboration), Phys. Rev. D, 98: 030001 (2018)
    [43] E. Kou et al (Belle-II Collaboration), arXiv: 1808.10567[hep-ex]
    [44] T. Barnes, S. Godfrey, and E. S. Swanson, Phys. Rev. D, 72: 054026 (2005) doi: 10.1103/PhysRevD.72.054026
    [45] D. Y. Chen and X. Liu, Phys. Rev. D, 84: 034032 (2011) doi: 10.1103/PhysRevD.84.034032
  • [1] M. B. Voloshin and V. I. Zakharov, Phys. Rev. Lett., 45: 688 (1980) doi: 10.1103/PhysRevLett.45.688
    [2] V. A. Novikov and M. A. Shifman, Z. Phys. C, 8: 43 (1981) doi: 10.1007/BF01429829
    [3] Y. P. Kuang and T. M. Yan, Phys. Rev. D, 24: 2874 (1981)
    [4] Y. P. Kuang, Front. Phys. China, 1: 19 (2006) doi: 10.1007/s11467-005-0012-6
    [5] Y.-P. Kuang, S. F. Tuan, and T. M. Yan, 37: 1210 (1988)
    [6] S. F. Tuan, A7, 3527 (1992); 42: 3207 (1990)
    [7] Y.-P. Kuang and T. M. Yan, 24: 2874 (1981); 41: 155 (1990)
    [8] M. B. Voloshin, 43: 1011(1986)
    [9] J. P. Lees et al (BaBar Collaboration), Phys. Rev. D, 84: 011104 (2011) doi: 10.1103/PhysRevD.84.011104
    [10] I. Adachi (Belle Collaboration), arXiv: 1105.4583[hep-ex]
    [11] A. Bondar et al (Belle Collaboration), Phys. Rev. Lett., 108: 122001 (2012) doi: 10.1103/PhysRevLett.108.122001
    [12] D. Y. Chen and X. Liu, Phys. Rev. D, 84: 094003 (2011) doi: 10.1103/PhysRevD.84.094003
    [13] E. S. Swanson, Phys. Rev. D, 91(3): 034009 (2015) doi: 10.1103/PhysRevD.91.034009
    [14] F. K. Guo, C. Hanhart, Q. Wang et al, Phys. Rev. D, 91(5): 051504 (2015) doi: 10.1103/PhysRevD.91.051504
    [15] W. S. Huo and G. Y. Chen, Eur. Phys. J. C, 76(3): 172 (2016) doi: 10.1140/epjc/s10052-016-4013-0
    [16] F.-K. Guo, C. Hanhart, Y. S. Kalashnikova et al, Phys. Rev. D, 93(7): 074031 (2016) doi: 10.1103/PhysRevD.93.074031
    [17] Q. Wang, V. Baru, A. A. Filin et al, Phys. Rev. D, 98(7): 074023 (2018) doi: 10.1103/PhysRevD.98.074023
    [18] D. Y. Chen, J. He, X. Q. Li et al, Phys. Rev. D, 84: 074006 (2011) doi: 10.1103/PhysRevD.84.074006
    [19] D. Y. Chen, X. Liu and S. L. Zhu, Phys. Rev. D, 84: 074016 (2011) doi: 10.1103/PhysRevD.84.074016
    [20] C. Meng and K. T. Chao, Phys. Rev. D, 77: 074003 (2008) doi: 10.1103/PhysRevD.77.074003
    [21] C. Meng and K. T. Chao, Phys. Rev. D, 78: 074001 (2008) doi: 10.1103/PhysRevD.78.074001
    [22] B. Wang, X. Liu, and D. Y. Chen, Phys. Rev. D, 94(9): 094039 (2016) doi: 10.1103/PhysRevD.94.094039
    [23] Y. A. Simonov and A. I. Veselov, Phys. Rev. D, 79: 034024 (2009) doi: 10.1103/PhysRevD.79.034024
    [24] B. Wang, H. Xu, X. Liu et al, Front. Phys. (Beijing), 11: 111402 (2016) doi: 10.1007/s11467-016-0564-7
    [25] Y. H. Chen, M. Cleven, J. T. Daub et al, Phys. Rev. D, 95(3): 034022 (2017) doi: 10.1103/PhysRevD.95.034022
    [26] Y. H. Chen and F. K. Guo, Phys. Rev. D, 100(5): 054035 (2019) doi: 10.1103/PhysRevD.100.054035
    [27] M. Cleven, F.-K. Guo, C. Hanhart et al, Eur. Phys. J. A, 47: 120 (2011) doi: 10.1140/epja/i2011-11120-6
    [28] F.-K. Guo, C. Hanhart, and U.-G. Meißner, Phys. Rev. Lett., 103: 082003 (2009); 104: 109901(E) (2010)
    [29] S. Fleming and T. Mehen, Phys. Rev. D, 78: 094019 (2008) doi: 10.1103/PhysRevD.78.094019
    [30] G. Burdman and J. F. Donoghue, Phys. Lett. B, 280: 287 (1992) doi: 10.1016/0370-2693(92)90068-F
    [31] M. B. Wise, Phys. Rev. D, 45: R2188 (1992) doi: 10.1103/PhysRevD.45.R2188
    [32] T. M. Yan, H. Y. Cheng, C. Y. Cheung et al, Phys. Rev. D, 46: 1148 (1992); 55: 5851(E) (1997)
    [33] R. Casalbuoni, A. Deandrea, N. Di Bartolomeo et al, Phys. Rept., 281: 145 (1997) doi: 10.1016/S0370-1573(96)00027-0
    [34] F. Bernardoni et al (ALPHA Collaboration), Phys. Lett. B, 740: 278 (2015) doi: 10.1016/j.physletb.2014.11.051
    [35] F.-K. Guo, C. Hanhart, G. Li et al, Phys. Rev. D, 83: 034013 (2011) doi: 10.1103/PhysRevD.83.034013
    [36] M. Cleven, Q. Wang, F.-K. Guo et al, Phys. Rev. D, 87: 074006 (2013) doi: 10.1103/PhysRevD.87.074006
    [37] T. Mehen and J. W. Powell, Phys. Rev. D, 88: 034017 (2013) doi: 10.1103/PhysRevD.88.034017
    [38] I. W. Stewart, Nucl. Phys. B, 529: 62 (1998) doi: 10.1016/S0550-3213(98)00374-5
    [39] J. Hu and T. Mehen, Phys. Rev. D, 73: 054003 (2006) doi: 10.1103/PhysRevD.73.054003
    [40] A. Garmash et al (Belle Collaboration), Phys. Rev. Lett., 116(21): 212001 (2016) doi: 10.1103/PhysRevLett.116.212001
    [41] Y.-H. Chen, J. T. Daub, F.-K. Guo et al, Phys. Rev. D, 93: 034030 (2016) doi: 10.1103/PhysRevD.93.034030
    [42] M. Tanabashi et al (Particle Data Group Collaboration), Phys. Rev. D, 98: 030001 (2018)
    [43] E. Kou et al (Belle-II Collaboration), arXiv: 1808.10567[hep-ex]
    [44] T. Barnes, S. Godfrey, and E. S. Swanson, Phys. Rev. D, 72: 054026 (2005) doi: 10.1103/PhysRevD.72.054026
    [45] D. Y. Chen and X. Liu, Phys. Rev. D, 84: 034032 (2011) doi: 10.1103/PhysRevD.84.034032
  • 加载中

Cited by

1. Liu, S.-D., Wang, F., Jia, Z.-S. et al. Pionic transitions of the spin-2 partner of X (3872) to χcJ[J]. Physical Review D, 2024, 110(5): 054048. doi: 10.1103/PhysRevD.110.054048
2. Jia, Z.-S., Zhang, Z.-H., Qin, W.-H. et al. Hunting for Xb via hidden bottomonium decays Xb →ππχbJ[J]. Physical Review D, 2024, 109(3): 034017. doi: 10.1103/PhysRevD.109.034017

Figures(4)

Get Citation
Yun-Hua Chen. Predictions of Υ(4S)hb(1P,2P)π+π transitions[J]. Chinese Physics C, 2020, 44(2): 023103. doi: 10.1088/1674-1137/44/2/023103
Yun-Hua Chen. Predictions of Υ(4S)hb(1P,2P)π+π transitions[J]. Chinese Physics C, 2020, 44(2): 023103.  doi: 10.1088/1674-1137/44/2/023103 shu
Milestone
Received: 2019-09-11
Article Metric

Article Views(3099)
PDF Downloads(24)
Cited by(2)
Policy on re-use
To reuse of Open Access content published by CPC, for content published under the terms of the Creative Commons Attribution 3.0 license (“CC CY”), the users don’t need to request permission to copy, distribute and display the final published version of the article and to create derivative works, subject to appropriate attribution.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Predictions of Υ(4S)hb(1P,2P)π+π transitions

    Corresponding author: Yun-Hua Chen, yhchen@ustb.edu.cn
  • School of Mathematics and Physics, University of Science and Technology Beijing, Beijing 100083, China

Abstract: We study the contributions of intermediate bottomonium-like Zb states and the bottom meson loops in the heavy quark spin flip transitions Υ(4S)hb(1P,2P)π+π. Depending on the constructive or destructive interferences between the Zb-exchange and the bottom meson loops mechanisms, we predict two possible branching ratios for each process: BRΥ(4S)hb(1P)π+π(1.2+0.80.4×106) or (0.5+0.50.2×106), and BRΥ(4S)hb(2P)π+π(7.1+1.71.1×1010) or (2.4+0.20.1×1010). The contribution of the bottom meson loops is found to be considerably larger than that of the Zb-exchange in the Υ(4S)hb(1P)ππ transitions, while its decay rates are not comparable to those of heavy quark spin conserved Υ(4S)Υ(1S,2S)ππ processes. We also predict the contribution of the charm meson loops in the branch fractions of ψ(3S,4S)hc(1P)ππ.

    HTML

    1.   Introduction
    • Hadronic transitions Υ(mS)Υ(lS)ππ and Υ(mS)hb(nP)ππ are important processes for understanding heavy-quarkonium dynamics and low-energy quantum chromodynamics (QCD). Because bottomonia are expected to be compact and nonrelativistic, the QCD multipole expansion (QCDME) method [14] is often used to analyze these transitions, where pions are emitted because of the hadronization of soft gluons. The decay rates of Υ(2S,3S)Υ(1S,2S)ππ can be well described by QCDME [5]. Because the total spin of bˉb system in Υ(mS) and hb(nP) are 1 and 0, respectively, in general, the heavy quark spin flip Υ(mS)hb(nP)ππ processes are expected to be suppressed, compared with the heavy quark spin conserved Υ(mS)Υ(nS)ππ processes. Within the framework of QCDME, studies [57] predicted that the branching fraction of Υ(3S)Υ(1P)ππ is suppressed by two orders of magnitude, relative to that of Υ(3S)Υ(1S)ππ, while Ref. [8] predicted a suppression of at least three orders of magnitude. The prediction of Ref. [8] is supported by experimental data [9]. In the decay processes Υ(5S)Υ(lS)π+π (l=1,2,3) and Υ(5S)hb(nP)π+π (n=1,2) where the two charged bottomoniumlike resonances Zb(10610)± and Zb(10650)± were observed, Υ(5S)hb(nP)π+π proceeds at a rate comparable to the Υ(5S)Υ(lS)π+π processes [10, 11]. The mechanism that mitigates this expected suppression has remained controversial. In Refs. [12, 13], the Υ(5S)hb(nP)π+π processes were interpreted via bottom meson loops mechanism, while genuine S-matrix Zb poles are required as in Refs. [1417]. The meson loops mechanism has been investigated by many previous works [1824] to study the dipion and η transitions of higher charmonia and bottomonia because the branch ratios and dipion invariant mass spectra cannot be described by QCDME.

      In this work, we will study whether the bottom meson loops mechanism can produce the Υ(4S)hb(nP)π+π transitions at decay ratios comparable to Υ(4S)Υ(lS)π+π. Because in the dipion emission processes of the Υ(4S) the crossed-channel exchanged Zb cannot be on-shell, these transitions are expected to be good channels to study the bottom meson loops' effect. In our previous works [25, 26], by using the nonrelativistic effective field theory (NREFT), we calculated the effects of the bottom meson loops as well as the Zb-exchange in the Υ(4S)Υ(1S,2S)ππ processes, and we found that the experimental data can be described well. Here, within the same theoretical scheme, we will calculate the contributions of the bottom meson loops and Zb-exchange in the Υ(4S)hb(nP)π+π processes, and give theoretical predictions for the decay branching ratios. We find that the contribution of the bottom meson loops is considerably larger than that of the Zb-exchange in the Υ(4S)hb(1P)π+π process, and it cannot produce a rate comparable to that of Υ(4S)Υ(1S,2S)π+π.

      The remainder of this paper is organized as follows. In Sec. 2, the theoretical framework is described in detail. In Sec. 3, we provide the theoretical predictions for the decay branching fractions of Υ(4S)hb(1P,2P)π+π, and discuss the contributions of different mechanisms. The study is concluded in Sec. 4.

    2.   Theoretical framework

      2.1.   Lagrangians

    • To calculate the contribution of the mechanism Υ(mS)Zbπhb(nP)ππ, we need the effective Lagrangians for the ZbΥπ and Zbhbπ interactions [27],

      LZbΥπ=j=1,2CZbjΥ(mS)πΥi(mS)Zibjuμvμ+h.c.,

      (1)

      LZbhbπ=j=1,2gZbjhb(nP)πϵijkZibjujhkb+h.c.,

      (2)

      where Zb1 and Zb2 denote Zb(10610) and Zb(10650), respectively, and vμ=(1,0) is the velocity of the heavy quark. The Zb states are given in the matrix as

      Zibj=(12Z0ibjZ+ibjZibj12Z0ibj).

      (3)

      The pions can be parametrized as Goldstone bosons of the spontaneous breaking of the chiral symmetry:

      uμ=i(uμuuμu),u=exp(iΦ2Fπ),Φ=(12π0π+π12π0),

      (4)

      where Fπ=92.2MeV is the pion decay constant.

      To calculate the box diagrams, we need the Lagrangian for the coupling of the Υ to the bottom mesons and the coupling of the hb to the bottom mesons [27, 28],

      LΥHH=igJHH2JHaσˉHa+h.c.,

      (5)

      LhbHH=ig12hibHaσiˉHa+h.c.,

      (6)

      where JΥσ+ηb denotes the heavy quarkonia spin multiplet, Ha=Vaσ+Pa with Pa(Va)=(B(),ˉB()0) collects the bottom mesons, and ABA(B)(A)B. We also need the Lagrangian for the axial coupling of the pion fields to the bottom and antibottom mesons, which at the lowest order in heavy-flavor chiral perturbation theory is given by [2933]

      LHHΦ=gπ2ˉHaσuabˉHbgπ2HaHbσuba,

      (7)

      where ui=2iΦ/F+O(Φ3) denotes the three-vector components of uμ, as defined in Eq. (4). Here, we use gπ=0.492±0.029 from a recent lattice QCD calculation [34].

    • 2.2.   Power counting of loops

    • Because the Υ(4S) meson is above the BˉB threshold and decays predominantly into BˉB pairs, the loop mechanism with intermediate bottom mesons may be important in the transitions Υ(4S)hb(nP)π+π. By following the formalism set-up based on NREFT [28, 35, 36], we will analyze the power counting of different types of loops. In NREFT, the expansion parameter is the velocity of the intermediate heavy meson, namely νX=|mXmB()mB()|/mB(), which is small because the bottomonia X are close to the B()ˉB() thresholds. In this power counting, each nonrelativistic propagator scale as 1/ν2, and the measure of one-loop integration scales as d4lν5.

      There are five different kinds of loop contributions, namely the box diagrams displayed in Fig. 1(b), (c), triangle diagrams displayed in Fig. 2(a)(c), and the bubble loop in Fig. 2(d). We analyze them one by one as follows:

      Figure 1.  Feynman diagrams considered for Υ(mS)hb(nP)ππ processes. Crossed diagrams of (a) and (b) are not shown explicitly.

      Figure 2.  Loop diagrams not considered in the calculations. The corresponding power counting arguments are given in the main text.

      First, we analyze the power counting of the box diagrams, namely Fig. 1(b), (c). As shown in Eq. (7), the vertex of B()B()π is proportional to the external momentum of the pion qπ. The ΥB()ˉB() vertex is in a P-wave, and the hbB()ˉB() vertex is in an S-wave; therefore, the loop momentum must contract with the external pion momentum, and hence the P-wave vertex scales as O(qπ). Thus, the box diagrams scales as ν5q3π/ν8=q3π/ν3.

      For the triangle diagram Fig. 2(a), the leading ΥB()ˉB()π vertex given by gJHHπJˉHaHbu0ab [37] is proportional to the energy of the pion, Eπqπ. Therefore, Fig. 2(a) is counted as mBν5q2π/ν6=mBq2π/ν, where the factor mB is introduced to match the dimension with the scaling for the box diagrams.

      In Fig. 2(b), the leading hbB()ˉB()π vertex given by ghbHHπhibHaσjˉHbϵijkukab [29] is proportional to the momentum of the pion qπ. The loop momentum due to the ΥB()ˉB() coupling has to contract with the external pion momentum. Thus, Fig. 2(b) scales as ν5q3π/ν6=q3π/ν.

      The leading B()B()ππ vertex comes from the chiral derivative term Ha(iD0)baHb=Ha(i0iV0)baHb [38, 39], in which the pion pair produced by the vector current, Vμ=12(uμu+uμu), cannot form a positive-parity and C-parity state; therefore, this leading vertex does not contribute to the Υ(mS)hb(nP)ππ processes. Isoscalar, PC=++ pion pairs only enter in the next order O(q2π) from point vertices. Therefore, Fig. 2(c) scales as ν5q3π/ν6=q3π/ν.

      In Fig. 2(d), both the initial and final vertices are proportional to qπ; therefore, the bubble loop scales as mBν5q2π/ν4=mBq2πν.

      Therefore, we expect that the ratios of the contributions of the box diagrams, triangle diagram Fig. 2(a)(c), and the bubble loop Fig. 2(d) are

      q3πν3:mBq2πν:q3πν:q3πν:mBq2πν=1:mBν2qπ:ν2:ν2:mBν4qπ,

      (8)

      where qπ(mΥ(4S)mhb(nP))/2 and ν=(νΥ(4S)+νhb(nP))/2, with νΥ(4S)0.06, νhb(1P)0.35, and νhb(2P)0.24. Thus, for the Υ(4S)hb(1P)π+π transition, the ratios in Eq. (8) are 1:0.67:0.04:0.04:0.03. For the Υ(4S)hb(2P)π+π transition, the ratios are 1:0.75:0.02:0.02:0.02. Therefore, according to the power counting, the box diagrams and triangle diagram in Fig. 2(a) are dominant among the loop contributions, and they are of the same order. The Υ(4S) is below the B()ˉB()π threshold, and the coupling gJHHπ in the triangle diagram Fig. 2(a) is unknown. Thus, for a rough estimation of the loop contributions, we will only calculate the box diagrams in the present study. All the box and triangle loop contributions discussed here are ultraviolet-finite and do not require additional introduction of counterterms.

    • 2.3.   Tree-level amplitudes and box diagram calculation

    • The decay amplitude for

      Υ(mS)(pa)hb(nP)(pb)π(pc)π(pd)

      (9)

      is described in terms of the Mandelstam variables

      s=(pc+pd)2,t=(papc)2,u=(papd)2.

      (10)

      By using the effective Lagrangians in Eqs. (1) and (2), the tree amplitude of Υ(mS)Zbπhb(nP)ππ can be obtained as

      MZb=2mΥ(mS)mhb(nP)F2πϵabjϵaΥ(mS)ϵbhb(nP)i=1,2mZbiCZbiΥ(mS)πgZbihb(nP)π{p0cpjd1tm2Zbi+p0dpjc1um2Zbi}.

      (11)

      The nonrelativistic normalization factor mY has been multiplied with the amplitude for every heavy particle, with Y=Υ(mS),hb(nP),Zbi. The widths of the Zb states are neglected in the present study, because they are of the order of 10MeV and are considerably smaller than the difference between the Zb masses and the Υ(mS)π/hb(nP)π threshold.

      Now, we discuss the calculation of the box diagrams. In the box diagrams Fig. 1(b) and (c), we denote the top left intermediate bottom meson as M1, and the other intermediate bottom mesons as M2, M3, and M4, in counterclockwise order. For the pseudoscalar or vector content of [M1,M2,M3,M4], there are twelve possible patterns and we number them in order: 1, [PPPV]; 2, [PPVV]; 3, [PVPV]; 4, [PVVP]; 5, [VVPP]; 6, [VPVP]; 7, [VPPV]; 8, [PVVV]; 9, [VPVV]; 10, [VVPV]; 11, [VVVP]; 12, [VVVV]. For each pattern, we also need to consider six possibilities of different flavor of the intermediate bottom mesons: [B()+,B(),B()+,B()0], [B()+,B(),ˉB()0,B()0], [B()0,ˉB()0,B()0,B()+], [B(),B()+,B(),ˉB()0], [ˉB()0,B()0,B()+,B()], and [ˉB()0,B()0,ˉB()0,B()]. The full amplitude contains the sum of all possible amplitudes.

      For the tensor reduction of the loop integrals, it is convenient to define q=pb and thperpendicular momentum q=pcq(qpc)/q2, which satisfy qq=0. The result of the amplitude of the box diagrams can be written as

      Mloop=ϵaΥ(mS)ϵbhb(nP){ϵabiqiA1+ϵabiqiA2+ϵbijqiqjqaA3+ϵbijqiqjqaA4+ϵaijqiqjqbA5+ϵaijqiqjqbA6}.

      (12)

      Details on the analytic calculation of the box diagrams and explicit expressions of Ai (i=1,2,...,6) are given in Appendix A.

      The decay width for Υ(mS)hb(nP)ππ is given by

      Γ=s+st+t|MZb+Mloop|2dsdt768π3m3Υ(mS),

      (13)

      where the lower and upper limits are given as

      s=4m2π,s+=(mΥ(mS)mhb(nP))2,t±=14s{(m2Υ(mS)m2hb(nP))2[λ12(s,m2π,m2π)λ12(m2Υ(mS),s,m2hb(nP))]2},λ(a,b,c)=a2+b2+c22(ab+ac+bc).

      (14)
    3.   Phenomenological discussion
    • To estimate the contribution of the Zb-exchange mechanism, we need to know the coupling strengths of ZbΥ(4S)π and Zbhb(nP)π. The mass difference between Zb(10610) and Zb(10650) is considerably smaller than the difference between their masses and the Υ(mS)π/hb(nP)π threshold, and they have the same quantum numbers, and thus the same coupling structures as given by Eqs. (1) and (2). Therefore, it is very difficult to distinguish their effects from each other in the dipion transitions of Υ(4S), so we only used one Zb, the Zb(10610), which approximately combines both Zb states' effects. In Ref. [26], we studied the Υ(4S)Υ(mS)ππ processes to extract the coupling constant |CZbΥ(4S)π|=(3.3±0.1)×103, which combine the effects from both Zb states. For the couplings of Zbhb(nP)π, in principle, they can be extracted from the partial widths of the Zb states decay into hb(nP)π(n=1,2)

      |gZbhbπ|={6πF2πmZbΓZbhbπ|pf|3mhb}12,

      (15)

      where |pf|λ1/2(m2Zb,m2hb,m2π)/(2mZb). The branching fractions of the decays of both Zb states into hb(nP)π(n=1,2) has been given in [40], where the Zb line shapes were described using Breit-Wigner forms. If we naively use these branching fractions, we would obtain

      |gnaiveZb1hb(1P)π|=0.019±0.003,|gnaiveZb2hb(1P)π|=0.021±0.003,|gnaiveZb1hb(2P)π|=0.068±0.011,|gnaiveZb2hb(2P)π|=0.077±0.010.

      (16)

      Here, all the Zbhbπ couplings are labeled by a superscript "naive" because this is not the appropriate way to extract coupling strengths in this case; the Zb states are very close to the B()ˉB thresholds. Therefore, the Flatté parametrization for the Zb spectral functions should be used, which will lead to much larger partial widths into (bˉbπ) channels, and thus the relevant coupling strengths. As analyzed in Ref. [41], in the the Flatté parametrization the sum of the partial widths of the Zb(10610) other than that for the BˉB channel should be greater than the nominal width, which is approximately 20MeV. While summing over all the Υ(nS)π(n=1,2,3) and hb(mP)π(m=1,2), and the branching fractions in Ref. [40] is approximately 14% or 3MeV in terms of partial widths. Therefore, for a rough estimation, we use three times the results of Eq. (16):

      |gZbhb(1P)π|0.057,|gZbhb(2P)π|0.204.

      (17)

      We find that even after considering the enlarging factor of three for the couplings |gZbhb(nP)π|, the Zb-exchange contribution is still considerably smaller than the bottom meson loops contributions.

      In the calculation of the box diagrams, the coupling strength gJHH(4S) can be extracted from the measured open-bottom decay widths of the Υ(4S), and we have gJHH(4S)=1.43±0.01GeV3/2. For the hbBˉB coupling g1, we can use the results from Ref. [27]. In [27], the Zb-exchange mechanism in the Υ(5S)hb(1P,2P)ππ processes was studied assuming the Zb states are B()ˉB bound states and the physical coupling of the Zb states to the bottom and anti-bottom mesons, z1, as well as the product g1z1 was determined. By using their results z1=0.75+0.080.11GeV1/2 and g1z1=0.40±0.06GeV1, we can extract that g1=0.53+0.190.13GeV1/2. In [27], in order to reduce the number of free parameters, the couplings of hb(1P)BˉB and hb(2P)BˉB are assumed to be the same.

      By using the coupling strengths above, we can predict the decay branching fractions of Υ(4S) hb(1P,2P)π+π. Depending on the sign of the couplings in Eq. (17), the interferences can be constructive or destructive between the Zb-exchange and box graph mechanisms; therefore, there are two possible results for each process

      BRΥ(4S)hb(1P)π+π(1.2+0.80.4×106)or(0.5+0.50.2×106),BRΥ(4S)hb(2P)π+π(7.1+1.71.1×1010)or(2.4+0.20.1×1010).

      (18)

      We find that the BRΥ(4S)hb(1P)π+π is at least one order of magnitude smaller than the branching fractions BRΥ(4S)Υ(1S,2S)π+π, which are approximately 8×105 given in PDG [42], and the BRΥ(4S)hb(2P)π+π is small owing to the very small phase space. We discuss the Υ(4S)hb(1P)π+π transition in further detail. To illustrate the effects of the Zb-exchange and box graph mechanisms in Υ(4S)hb(1P)ππ, we give the predictions only including the Zb-exchange terms or only including the box diagrams

      BRZbΥ(4S)hb(1P)π+π=0.6+0.10.1×107,BRBoxΥ(4S)hb(1P)π+π=0.8+0.70.3×106.

      (19)

      We can observe that the bottom meson loops contribution is considerably larger than the Zb-exchange contribution, while it is two orders of magnitude smaller than the Υ(4S)Υ(1S,2S)π+π transitions. Note that the direct gluon hadronization mechanism contribution within QCDME for the Υ(4S)hb(1P)π+π process has not been calculated thus far. In references [58], QCDME predict that the branching fraction of Υ(3S)hb(1P)ππ is 2-3 orders of magnitude suppressed compared to of Υ(3S)Υ(1S)ππ. The three orders of magnitude suppression is supported by experiment [9]. Because the mass difference between the Υ(4S) and hb(1P) is approximately 0.68 GeV, the pions in the Υ(4S)hb(1P)π+π process can also be considered to be in the soft region. If one approximates that the gluon hadronization mechanism within QCDME in Υ(4S)hb(1P)π+π is also 2-3 orders of magnitudes suppressed compared to that in the Υ(4S)Υ(1S)π+π process, as in the Υ(3S) decay cases. Then, the gluon hadronization mechanism contribution is at most at the same order of the bottom meson loops contribution. Owing to the lack of exact information on the gluon hadronization within QCDME and the neglecting of the triangle diagram Fig. 2(a) as discussed in section 2.3, it is important to note that the results presented in this paper are order-of-magnitude estimates.

      In Fig. 3, we plot the distributions of the ππ and hbπ invariant mass spectra, and the distribution of cosθ, where θ is defined as the angle between the initial Υ(mS) and the π+ in the rest frame of the ππ system. To illustrate the effects of different mechanisms, the contributions of the box diagrams, Zb-exchange, their sum with the constructive interference, and the sum with destructive interference are indicated by dark green dashed, magenta dotted, red solid, and blue dot-dashed lines, respectively. There is a broad bump at approximately 0.5 GeV in the dipion invariant mass distribution. The ππ invariant mass spectra with unknown normalization predicted within QCDME in Ref. [5] showed a peak at low ππ masses. Thus, the ππ invariant mass spectra can be useful to identify the effects of the bottom meson loops and the gluon hadronization mechanism with future experimental data. Further, the angular distribution is far from flat. In the Υ(mS)hb(nP)ππ process, the isospin conservation combined with Bose symmetry requires the pions to have even relative angular momentum. Therefore, there is a large D-wave component from the box diagrams if higher partial waves are neglected.

      Figure 3.  (color online) Theoretical predictions of the distributions of the ππ and hbπ invariant mass spectra, and helicity angular distributions in the Υ(4S)hb(1P)ππ process. Dark green dashed, magenta dotted, red solid, and blue dot-dashed lines represent the contributions of the box diagrams, Zb-exchange, their sum with constructive interference, and the sum with destructive interference, respectively.

      The Υ(mS)hb(nP)ππ are heavy quark spin flip processes, and they are forbidden in the heavy quark limit. We checked that in the heavy quark limit, i.e. mB=mB, all the box diagrams were cancelled with each other; therefore, the bottomed loops did not contribute to the Υ(mS)hb(nP)ππ transitions. With the small mass splitting of B and B in the real world, as shown in Eqs. (18) and (19), the bottomed meson loops contribution does not produce Υ(4S)hb(1P)ππ at a rate comparable to the heavy quark spin conserved Υ(4S)Υ(1S,2S)ππ transitions. Note that the datasets collected at Υ(4S) by BABAR and Belle II collaborations are 471×106 and 772×106 [43], respectively. Thus, they should contain several hundreds of Υ(4S)hb(1P)ππ events according to our calculation. We hope future experimental analysis by BABAR and Belle can test our predictions. As stated in the introduction, the observed Υ(5S)hb(nP)π+π proceed at a rate comparable to the Υ(5S)Υ(lS)π+π processes [10, 11]. The enhancements may be caused by the effects of the on-shell Zb exchange and the two-cut condition complexity of the bottom meson loops in the Υ(5S) decays. A detailed analysis of the Υ(5S) hb(nP)π+π processes is beyond the scope of this study.

      Owing to the similarity between the bottomonium and charmonium families, we can extend the box diagrams calculation to give a rough estimation of the branch fractions of the ψ(3S)(ψ(4040))hc(1P)π+π and ψ(4S)(ψ(4415))hc(1P)π+π transitions. The relevant Feynman diagrams can be obtained by replacing the external Υ(mS) and hb(nP) by ψ(mS) and hc(nP), respectively, and replacing the intermediate B() by D() in Fig. 1(b) and (c). The experimental decay widths of the ψ(3S,4S)D()ˉD() transitions are not given in PDG, and we will use the theoretical predictions of the decay widths in Ref. [44] to estimate the coupling strengths gJHH(ψ(3S)) and gJHH(ψ(4S)). Because among the different decay modes D()ˉD(), the DˉD and DˉD modes are dominant for ψ(3S) and ψ(4S), respectively, we will use the corresponding coupling constants in the calculation, namely gJHH(ψ(3S))=gψ(3S)DˉD=0.97GeV3/2 and gJHH(ψ(4S))=gψ(4S)DˉD=0.25GeV3/2. For the hcDˉD coupling, we use the result from Ref. [45], ghcDˉD=(mχc0/3)/(fχc0)=(3.415/3)/(0.297)GeV1/2=3.59GeV1/2. The predictions of the box diagrams contributions to the branch fractions of ψ(3S,4S)hc(1P)π+π are

      BRBoxψ(3S)hc(1P)π+π=2.9×105,BRBoxψ(4S)hc(1P)π+π=4.5×103.

      (20)

      The prediction of BRBoxψ(3S)hc(1P)π+π is below the upper limit given in PDG [42]. As expected the branch fractions BRBoxψ(3S,4S)hc(1P)π+π are considerably larger than BRBoxΥ(4S)hb(1P)π+π, because the mass splitting of D and D is considerably larger than that of B and B. Note that this is just a preliminary rough estimation, owing to the lack of sufficient information concerning the ψ(3S,4S)D()ˉD() coupling constants and the neglecting of the loop diagrams with intermediate D1 state in the present calculation. A detailed theoretical study of the ψ(3S,4S) hc(1P)π+π transitions will be pursued in the future.

    4.   Conclusions
    • In this study, we investigate the effects of Zb exchange and bottom meson loops in the heavy quark spin flip transitions Υ(4S)hb(nP)ππ(n=1,2). The bottom meson loops are treated in the NREFT scheme, in which the dominant box diagrams are considered. We find that the bottom meson loops contribution is considerably larger than the Zb-exchange contribution in the Υ(4S) hb(1P)ππ transition, while it can not produce decay rates comparable to the heavy quark spin conserved Υ(4S)Υ(1S,2S)ππ processes. The theoretical prediction of the decay rate and the dipion invariant mass spectra of Υ(4S)hb(1P)ππ in this work may be useful for identifying the effect of the bottom meson loops with future experimental analysis. We also predict the branch fractions of ψ(3S,4S)hc(1P)ππ contributed from the charm meson loops.

    Appendix A: Remarks on box diagrams and four-point integrals
    • In this appendix, first, we discuss the parametrization and simplification of the scalar four-point integrals in the box diagrams. Then, we introduce a tensor reduction scheme to deal with higher-rank loop integrals. Finally, we give the amplitude of the box diagrams for the Υ(mS)hb(nP)ππ process.

    • A.1.   Scalar four-point integrals

    • For the first topology, as shown in Fig. A1, the scalar integral evaluated for the initial bottomonium at rest (p=(M,0)) reads

      Figure A1.  Kinematics used for calculating four-point integrals.

      We are grateful to the referees' useful suggestions and constructive remarks, which helped formulate the present version of this manuscript. We are grateful to Martin Cleven for the collaboration in the early stages of this study. We acknowledge Guo-Ying Chen, Meng-Lin Du, and Qian Wang for their helpful discussions, and Feng-Kun Guo for a careful reading of the manuscript and valuable comments.

      J(0)1id4l(2π)41[l2m21+iϵ][(pl)2m22+iϵ][(lq1q2)2m23+iϵ][(lq1)2m24+iϵ]i16m1m2m3m4×d4l(2π)41[l0l22m1m1+iϵ][l0M+l22m2+m2iϵ]1[l0q01q02(l+q)22m3m3+iϵ][l0q01(lq1)22m4m4+iϵ].

      By performing the contour integration, we find

      μ12μ23μ242m1m2m3m4d3l(2π)31[l2+c12iϵ][l2+2μ23m3lq+c23iϵ][l22μ24m4lq1+c24iϵ],

      where we defined

      c122μ12(m1+m2M),c232μ23(m2+m3M+q01+q02+q22m3),c242μ24(m2+m4M+q01+q212m4),μij=mimjmi+mj.

      The second topology in Fig. A1 is just the crossed diagram of the first topology with q1q2, so the scalar integral reads

      J(0)2=μ12μ23μ242m1m2m3m4d3l(2π)3×1[l2+c12iϵ][l2+2μ23m3lq+c23iϵ][l22μ24m4lq2+c24iϵ],

      where

      c242μ24(m2+m4M+q02+q222m4).

      For the third topology, we have

      J(0)3id4l(2π)41[l2m21+iϵ][(pl)2m22+iϵ][(pq2l)2m23+iϵ][(lq1)2m24+iϵ]i16m1m2m3m4d4l(2π)41[l0l22m1m1+iϵ][l0M+l22m2+m2iϵ]×1[l0+q02M+(l+q2)22m3+m3iϵ][l0q01(lq1)22m4m4+iϵ].

      By the contour integration, we find

      μ12μ342m1m2m3m4d3l(2π)31[l2+d12iϵ][l22μ34m4lq12μ34m3lq2+d34iϵ]×[μ24[l22μ24m4lq1+d24iϵ]+μ13[l2+2μ13m3lq2+d13iϵ]],

      where we defined

      d122μ12(m1+m2M),d342μ34(m3+m4q0+q212m4+q222m3),d242μ24(m2+m4M+q01+q212m4),d132μ13(m1+m3M+q02+q222m3).

      In all the three cases, the remaining three-dimensional momentum integration will be carried out numerically.

    • A.2.   Tensor reduction

    • Because the ΥB()ˉB() vertex scales with the momentum of the bottom meson pair, for topology I we have to deal with

      μ12μ23μ242m1m2m3m4d3l(2π)3f(l)[l2+c12iϵ][l2+2μ23m3lq+c23iϵ][l22μ24m4lq1+c24iϵ],

      where f(l)={1,li} for the fundamental scalar and vector integrals, respectively. A convenient parametrization of the tensor reduction is

      J(1)i1=μ12μ23μ242m1m2m3m4d3l(2π)3li[l2+c1iϵ][l22μ23m3lq+c2iϵ][l22μ24m4lq1+c3iϵ]qiJ(1)1+qi1J(2)1,

      where q1=q1q(qq1)/q2. The expressions of the scalar integrals J(r)1 can easily be disentangled and should be evaluated numerically. The corresponding expressions for topology II and III can be obtained by changing the denominators accordingly.

    • A.3.   Amplitudes

    • We define the scalar integrals J1(i,r,k) based on the J(r)1 in the tensor reduction of vector integral in Eq. (A10), where i=1,2,3 denotes the three topologies of the box diagrams shown in Fig. A1, r=1,2 refers to the two components J(r)1, and k=1,2,...,12 represents the twelve patterns with different pseudoscalar or vector content of the intermediate bottom mesons in [M1,M2,M3,M4] as displayed in Sec. 2.3.

      We give the amplitude of the box diagrams for the Υ(mS)hb(nP)ππ process, namely the Al(l=1,2,...,6) in the Eq. (12).

      A1=8g1gJHHg2πF2πq2{q2{pcpd[J1(1,1,3)+J1(2,1,3)+J1(3,1,8)]+pcq[J1(1,1,9)+J1(1,1,11)J1(2,1,12)+J1(3,1,9)+J1(3,1,11)]+q2[J1(1,2,9)+J1(1,2,11)J1(2,2,12)+J1(3,2,9)+J1(3,2,11)]}+pcq{pcq[J1(1,1,9)+J1(1,1,11)J1(2,1,12)+J1(3,1,9)+J1(3,1,11)]+pdq[J1(1,1,12)J1(2,1,9)J1(2,1,11)+J1(3,1,10)]+q2[J1(1,2,9)+J1(1,2,11)J1(1,2,12)+J1(2,2,9)+J1(2,2,11)J1(2,2,12)+J1(3,2,9)J1(3,2,10)+J1(3,2,11)]}},

      A2=8g1gJHHg2πF2π{pcpd[J1(1,2,3)+J1(2,2,3)+J1(3,2,8)]+pcq[J1(1,1,9)+J1(1,1,11)J1(2,1,12)+J1(3,1,9)+J1(3,1,11)]+pdq[J1(1,1,12)J1(2,1,9)J1(2,1,11)+J1(3,1,10)]+q2[J1(1,2,9)+J1(1,2,11)J1(1,2,12)+J1(2,2,9)+J1(2,2,11)J1(2,2,12)+J1(3,2,9)J1(3,2,10)+J1(3,2,11)]},

      A3=8g1gJHHg2πF2πq2{q2[J1(1,1,9)+J1(1,1,11)J1(1,1,12)J1(1,2,2)+J1(1,2,9)+J1(1,2,10)+J1(1,2,12)J1(2,1,9)+J1(2,1,11)J1(2,1,12)+J1(2,2,2)J1(2,2,10)J1(2,2,11)J1(3,1,9)J1(3,1,10)+J1(3,1,11)J1(3,2,2)+J1(3,2,9)+J1(3,2,10)+J1(3,2,12)]+pcq[J1(1,2,9)J1(1,2,11)+J1(1,2,12)+J1(2,2,9)J1(2,2,11)+J1(2,2,12)+J1(3,2,9)+J1(3,2,10)J1(3,2,11)]},

      A4=8g1gJHHg2πF2πq4{q4[J1(1,1,2)J1(1,1,10)J1(1,1,11)J1(2,1,2)+J1(2,1,9)+J1(2,1,10)+J1(2,1,12)+J1(3,1,2)J1(3,1,11)J1(3,1,12)]+q2pcq[J1(1,1,9)J1(1,1,11)+J1(1,1,12)J1(1,2,9)+J1(1,2,11)J1(1,2,12)+J1(2,1,9)J1(2,1,11)+J1(2,1,12)J1(2,2,9)+J1(2,2,11)J1(2,2,12)+J1(3,1,9)+J1(3,1,10)J1(3,1,11)J1(3,2,9)J1(3,2,10)+J1(3,2,11))(pcq)2[J1(1,2,9)J1(1,2,11)+J1(1,2,12)+J1(2,2,9)J1(2,2,11)+J1(2,2,12)+J1(3,2,9)+J1(3,2,10)J1(3,2,11)]},

      A5=8g1gJHHg2πF2πq2{q2[J1(1,1,6)+J1(1,2,8)+J1(1,2,10)J1(2,1,6)+J1(2,2,6)J1(2,2,8)J1(2,2,10)J1(3,1,3)+J1(3,1,6)J1(3,1,8)+J1(3,2,3)+J1(3,2,12)]+pcq[J1(1,2,6)+J1(2,2,6)+J1(3,2,3)J1(3,2,6)+J1(3,2,8)]},

      A6=8g1gJHHg2πF2πq4{q4[J1(1,1,6)J1(1,1,8)J1(1,1,10)+J1(2,1,8)+J1(2,1,10)J1(3,1,6)+J1(3,1,8)J1(3,1,12)]+q2pcq[J1(1,1,6)J1(1,2,6)+J1(2,1,6)J1(2,2,6)+J1(3,1,3)J1(3,1,6)+J1(3,1,8)J1(3,2,3)+J1(3,2,6)J1(3,2,8)](pcq)2[J1(1,2,6)+J1(2,2,6)+J1(3,2,3)J1(3,2,6)+J1(3,2,8)]}.

Reference (45)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return