Strong Coupling Constants of the Doubly Heavy $ \Xi_{QQ} $ Baryons with $ \pi $ Meson

The doubly charmed $\Xi_{cc}^{++} (ccu)$ state is the only listed baryon in PDG, which was discovered in the experiment. The LHCb collaboration gets closer to discovering the second doubly charmed baryon $\Xi_{cc}^{+} (ccd)$, hence the investigation of the doubly charmed/bottom baryons from many aspects is of great importance that may help us not only get valuable knowledge on the nature of the newly discovered states, but also in the search for other members of the doubly heavy baryons predicted by the quark model. In this context, we investigate the strong coupling constants among the $\Xi_{cc}^{+(+)}$ baryons and $\pi^{0(\pm)}$ mesons by means of light cone QCD sum rule. Using the general forms of the interpolating currents of the $\Xi_{cc}^{+(+)}$ baryons and the distribution amplitudes (DAs) of the $\pi$ meson, we extract the values of the coupling constants $g_{\Xi_{cc} \Xi_{cc} \pi}$. We extend our analyses to calculate the strong coupling constants among the b-partner baryons with $\pi$ mesons, as well, and extract the values of the strong couplings $g_{\Xi_{bb} \Xi_{bb} \pi}$. It is observed that the values of the couplings under study in bottom channels are about 4 times greater than those of the charmed channels. The results of this study may help experimental groups in the analyses of the data related to the strong coupling constants among the hadronic multiplets.


INTRODUCTION
The search for doubly heavy baryons and determination of their properties constitute one of the main directions of the research in the experimental and theoretical high energy physics. There is only one doubly charmed baryon, Ξ ++ cc , listed in the PDG. The searches for other members of the doubly heavy baryons in the experiments, as the natural outcomes of the quark model, are in progress. Theoretical investigations on properties of the doubly heavy baryons, are necessary as their results can help us better understand their structure and the dynamics of the QCD as the theory of the strong interaction.
The search for doubly heavy baryons is a long-standing issue. First evidence was reported by the SELEX experiment for Ξ ++ cc decaying into Λ + c K − π + and pD + K − in final states [1,2]. The mass measured by SELEX, averaged over the two decay modes, was found to be 3518.7 ± 1.7 MeV/c 2 . However, this has not been confirmed by any other experiments so far. The FOCUS [3], BaBar [4], LHCb [5] and Belle [6] experiments did not find any evidence up to 2017. In 2017, the doubly charmed baryon Ξ ++ cc was observed by the LHCb collaboration via the decay channel Ξ ++ cc → Λ + c K − π + π [7], and confirmed via measuring another decay channel Ξ ++ cc → Ξ + c π + [8]. The weighted average of its mass for the two decay modes was determined to be 3621.24±0.65(stat.)±0.31(syst.) MeV/c 2 . Recently, with a data sample corresponding to an integrated luminosity of 9 fb −1 at the centre-of-mass energies of 7, 8 and 13 TeV, the LHCb Collaboration published the results of a search for the doubly charmed baryon Ξ + cc [9]. The upper limit of the ratio of the production cross-sections between the Ξ + cc and Λ + c baryons times the branching fraction of the Ξ ++ cc → Λ + c K − π + decay, was improved by an order of magnitude than the previous search. However, still no significant signal is observed in the mass range from 3.4 to 3.8 GeV/c 2 . Future LHCb searches with further improved trigger conditions, additional Ξ + cc decay modes, and larger data samples should significantly increase the Ξ + cc signal sensitivity. Theoretical studies on the properties and nature of the doubly heavy baryons can play an important role in searching for other members and help us get useful knowledge on the internal structures of the observed resonances. There have been many theoretical efforts aimed at understanding the properties of the doubly-heavy baryon states, see e.g. Refs. . However, most researches are focused on the mass and weak decays of the doubly heavy baryons and the number of studies dedicated to their strong decays and the strong couplings of these baryons with other hadrons is very limited. In this context, we investigate the strong coupling constants among the doubly heavy Ξ cc /Ξ bb baryons and π mesons. We use the well established non-perturbtive method of Light Cone QCD sum rule (LCSR) (for more about this method see, e.g., [37][38][39][40][41] and references therein) as a powerful theoretical tool to calculate the coupling constants under study. In the calculations, we use the general forms of the interpolating currents for Ξ cc and Ξ bb baryons and DAs of the pions.
The outline of the paper is as follows. In Sec. 2, the light cone sum rules for the coupling constants of the doubly heavy baryons with π mesons are obtained. In Sec. 3, we present the numerical results and discussions, and Sec. 4 is reserved for our conclusions.

THEORETICAL FRAMEWORK
The starting point is to choose a suitable correlation function (CF) in terms of hadronic currents sandwiched between the QCD vacuum and the on-shell pseudoscalar meson (here the pion). The QCD vacuum interacts via vacuum fluctuations with the initial and final states and leads to non-perturbative contributions to the final results via the quark-quark, quark-gluon and gluon-gluon condensates. The DAs of π meson also contain non-perturbative information.
To calculate the physical observables like strong coupling constants, we have to calculate the CF in the large timelike momenta by inserting the full set of hadronic states in the CF and isolating the ground state from the continuum and excited states. On the other hand, in QCD or theoretical side, we need to study the CF in the large space-like momenta via the operator product expansion (OPE) that separates the perturbative and non-perturbative contributions in terms of different operators and distribution amplitudes of the particles under consideration. Thus, in QCD side, by contracting the hadronic currents in terms of quark fields we write the CF in terms of the light and heavy quarks propagators as well as wavefunctions of the considered meson in x-space. To proceed, we perform the Fourier transformation to transfer the calculations to the momentum space. To suppress the unwanted contributions coming from the higher states and continuum, we apply the Borel transformation as well as continuum subtraction, supplied by the quark-hadron duality assumption. The two representations of the same CF are then connected to each other via dispersion integrals. These procedures introduce some auxiliary parameters to the calculations, which are then fixed based on the standard prescriptions of the method.

Correlation Function
In order to calculate the strong coupling constants among doubly heavy baryons and light pseudoscalar π meson, we start our discussion by considering the following lightcone correlation function: where P(q) denotes the pseudoscalar mesons of momentum q, T is the time ordering operator, and η represents the interpolating currents of the doubly heavy baryons. The general expressions of the interpolating currents for the spin1/2 doubly heavy baryons in their symmetric and antisymmetric forms can be written as where t is an arbitrary auxiliary parameter and the case, t = −1 corresponds to the Ioffe current. Here Q ( ′ ) and q stand for the heavy and light quarks respectively; a, b, and c are the color indices, C stands for the charge conjugation operator and T denotes the transposition. For the doubly heavy baryons with two identical heavy quarks, the antisymmetric form of the interpolating current is zero and we just need to employ the symmetric form, η S . In the following, we calculate the correlation function in Eq. 1 in two different windows.

Physical Side
To obtain the physical side of correlation function, we insert complete sets of hadronic states with the same quantum numbers as the interpolating currents and isolate the ground state. After performing the integration over four-x, we get where dots in the above equation stand for the contribution of the higher states and continuum. To proceed we introduce the matrix elements for spin-1/2 baryons as where g B 1 B 2 P , representing the strong decay B 1 → B 2 P, is the strong coupling constant among the baryons B 1 and B 2 as well as the P meson, λ B 1 and λ B 2 are the residues of the corresponding baryons and u(q, s) is Dirac spinor with spin s. Putting the above equations all together, and performing summation over spins, we get the following representation of the correlator for the phenomenological side: where dots represent other structures come from spin summation as well as the contributions of higher states and continuum. We will use the explicitly presented structure to extract the value of the strong coupling constant, g B 1 B 2 P . We apply the double Borel transformation with respect to the variables p 2 1 = (p + q) 2 and p 2 2 = p 2 : where

QCD Side
In QCD side, the correlation function is calculated in deep Euclidean region with the help of OPE. To proceed, we need to determine the correlation function using the quark propagators and distribution amplitudes of the π meson. The Π QCD (p, q) can be written in the following general form: where the Π (p + q) 2 , p 2 is an invariant function that should be calculated in terms of QCD degrees of freedom as well as the parameters inside the DAs. Inserting, for instance, η S in the CF and using the Wick theorem to contract all the heavy quark fields we get the following expression in terms of the heavy quark propagators and π meson matrix elements: are the matrix elements for the light quark contents of the doubly heavy baryons. To proceed, we need to know the explicit expression for the heavy quark propagator that is where K 1 and K 2 are the modified Bessel functions of the second kind, and G µν ab ≡ G µν A t A ab with A = 1, 2 . . . 8, and t A = λ A /2, where λ A are the Gell-Mann matrices. The first two terms correspond to perturbative or free part and the rest belong to the interacting parts. The next step is to use the heavy quark propagator and the matrix elements π(q)|q b ′ α (0)q b β (x)|0 in Eq. 8. This leads to different kinds of contributions to the CF. Figs. 1 and 2 are the Feynman diagrams correspond to the leading and next-to-leading order contributions, respectively which are considered in this work. Only matrix elements corresponding to these diagrams are available. To calculate the leading order contribution, the heavy quark propagators are replaced by just their perturbative parts. This contribution can be computed using π meson two-particle DAs of twist two and higher.
The next-to-leading order contributions can also be calculated by choosing the gluonic parts in Eq. 9 for one of the heavy quark propagators and leaving the other with its perturbative term. They can be expressed in terms of pion three particles DAs. The terms involving more than one gluon field that proportional to four-particle DAs or more are neglected as they are not available Now, we concentrate on the strong decay Ξ ++ cc → Ξ ++ cc π 0 with the aim of calculating the corresponding strong coupling constant g Ξ ++ cc Ξ ++ cc π 0 . The other channels have similar procedures. To proceed, we replace the heavy quark propagators in 8 by their explicit expression and perform the summation over the color indices by applying the replacement Figure 1: The leading order diagram contributing to QCD side. Figure 2: The one-gluon exchange diagrams corresponding to the next-to-leading contributions.
Now, using the expression one can relate the CF to the DAs of the pion with different twists. Here the summation over J runs as Γ J = 1, γ 5 , γ µ , iγ 5 γ µ , σ µν / √ 2.
Following the similar way one can calculate the contributions involving the gluon field. As a result, the CF is found in terms of the QCD parameters as well as the matrix elements whose expressions in terms of the wave functions of the pion with different twists are given in the Appendix. Inserting the expression of the above-mentioned matrix elements in term of wave functions of different twists we get the CF in x space. This is followed by the Fourier and Borel transformations as well as continuum subtraction. To proceed we need to perform the Fourier transformation of the following kind: where the expressions in the brackets denote different possibilities arise in the calculations, the blank subscript in the left hand side indicates no indices regarding no x α in the configuration, G(α i ) and f (u) represent wave functions coming from the three and two-particle matrix elements and n is a positive integer. The measure is used in the calculations. To start the Fourier transformation, we use for positive integer n and We also use the following representation of the Bessel functions K ν (see also Ref. [42]): As an example let us consider the following generic form: We substitute Eqs. 15, 16 and 17 into 18. Then to perform the x-integration we go to the Euclidean space by Wick rotation and get where P = p + q(vα g + α q ) . Changing variables from y 1 and y 2 to ρ and z as and taking derivative with respect to P α and P β , we get Now we perform the double Borel transformation using where u 0 = M 2 1 /(M 2 1 + M 2 2 ). After integrating over ρ we have The next step is to perform the continuum subtraction in order to more suppress the contribution of the higher states and continuum. This is done by using the replacement, for n > 0. Finally, using some variable changing we set the lower limit of integration over s as (m Q 1 + m Q 2 ) 2 . By calculation of all the Fourier integrals and applying the Borel transformation and continuum subtraction, for QCD side of the calculations in Borel Scheme, we get, where we only present the explicit expression of the function Π (0) (M 2 , s 0 , t) and ignore to represent the Π (GG) (M 2 , s 0 , t) function, corresponding to the gluon condensate contributions, because of its very lengthy expression. The Π (0) (M 2 , s 0 , t) function is found as The sum rule for the coupling constant under study is found by matching the coefficients of the structure / q/ pγ 5 from both the physical and QCD sides. As a result, we get: As is seen, the sum rules for coupling constants contain the residues of doubly heavy baryons, which are borrowed from Ref. [10]. Similarly, we obtain the sum rules for other strong coupling constants under consideration. f π 130.2 ± 1.2 MeV Table 1: Some input values used in the calculations. They are mainly taken from [43], except the ones that the references are given next to the numbers.

NUMERICAL ANALYSIS
In this section, we numerically analyze the sum rules for the strong coupling constants of the π mesons with Ξ cc and Ξ bb baryons and discuss the results. The sum rules for the couplings g ΞccΞccπ , g Ξ bb Ξ bb π contain some input parameters like, quark masses the mass and decay constant of the π meson and the masses and residues of doubly heavy baryons. They were extracted from experimental data or calculated from nonperturbative methods. The values of some of these parameters together with quark masses are given in Tables 1. As we previously mentioned the values of the residues of baryons are used from Ref. [10]. Another set of important input parameters are the π meson wavefunctions of different twists, entering the DAs. These wavefunctions are given as [44,45]: , where C k n (x) are the Gegenbauer polynomials and h 10 = 7 4 η 4 w 4 + 3 20 a π 2 , g 0 = 1, The constants inside the wavefunctions are calculated at the renormalization scale of µ = 1 GeV 2 and they are given as a π 1 = 0, a π 2 = 0.44, η 3 = 0.015, η 4 = 10, w 3 = −3 and w 4 = 0.2 [44,45].
Finally, the sum rules for the coupling constants contain three auxiliary parameters: Borel mass parameter M 2 , continuum threshold s 0 and the general parameter t entered the general spin1/2 currents. We should find working regions of these parameters, at which the results of coupling constants have relatively small variations with respect to the changes of these parameters. To restrict these parameters, we employ the standard prescriptions of the method such as the pole dominance, convergence of the OPE and mild variations of the physical quantities with respect to the auxiliary parameters. As a result, we obtain the working region of the Borel parameter for the Ξ cc channel as The continuum threshold for this channel is chosen in the interval, For Ξ bb baryon these parameters lie in the intervals: The working window for the parameter t is obtained in Ref. [10] as |t| ≤ 2, which we also use in our analyses. As examples, we display the dependence of the strong coupling constant g Ξ ++ cc Ξ ++ cc π 0 , which is obtained from the sum rule for the strong coupling form factor at q 2 = m 2 π , with respect to M 2 and s 0 in Figs. 3 and 4 at t = −2.

SUMMARY AND CONCLUSIONS
The doubly charmed Ξ ++ cc (ccu) baryon is the only listed doubly heavy baryon in PDG discovered in the experiment so far. The LHCb collaboration gets closer to observing other member Ξ + cc (ccd), as well. Therefore, the investigation of the doubly charmed/bottom baryons from many aspects is of great importance that may help us in the course of search for new members of the doubly heavy baryons predicted by the quark model. The strong coupling constants among the hadronic multiplets are fundamental objects that can help us to explore the nature and structure of the participating particles as well as the properties of QCD as the theory of strong interaction.
We calculated the strong coupling constants g ΞccqΞccqπ 0,± and g Ξ bbq Ξ bbq π 0,± , with q being either u or d quark, in the framework of the light cone QCD sum rule and using the general form of the interpolating currents for the doubly heavy baryons and the π meson's DAs. Based on the standard prescriptions of the method, we fixed the auxiliary parameters entering the calculations. We extracted the values of the strong coupling constants at different channels. Our results may be checked via different theoretical models and approaches. The obtained results may help us in constructing the strong interaction potential among the doubly heavy baryons and the pseudoscalar mesons.
Our results may also help experimental groups in analyses of the obtained related data in hadron colliders.