Analysis of the strong vertices of ΣcΔD∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Sigma _{c}\Delta D^{*}$$\end{document} and ΣbΔB∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Sigma _{b}\Delta B^{*}$$\end{document} in QCD sum rules

In this work, we analyze the strong vertices ΣcΔD∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Sigma _{c}\Delta D^{*}$$\end{document} and ΣbΔB∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Sigma _{b}\Delta B^{*}$$\end{document} using the three-point QCD sum rules under the tensor structures iϵρταβpαpβ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i\epsilon ^{\rho \tau \alpha \beta }p_{\alpha }p_{\beta }$$\end{document}, pρp′τ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p^{\rho }p'^{\tau }$$\end{document} and pρpτ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p^{\rho }p^{\tau }$$\end{document}. We firstly calculate the momentum dependent strong coupling constants g(Q2)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g(Q^{2})$$\end{document} by considering contributions of the perturbative part and the condensate terms ⟨q¯q⟩\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\langle {\overline{q}}q\rangle $$\end{document}, ⟨gs2GG⟩\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\langle g_{s}^{2}GG \rangle $$\end{document}, ⟨q¯gsσGq⟩\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\langle {\overline{q}}g_{s}\sigma Gq\rangle $$\end{document} and ⟨q¯q⟩2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\langle {\overline{q}}q\rangle ^{2}$$\end{document}. By fitting these coupling constants into analytical functions and extrapolating them into time-like regions, we then obtain the on-shell values of strong coupling constants for these vertices. The results are g1ΣcΔD∗=5.13-0.49+0.39GeV-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g_{1\Sigma _{c}\Delta D^{*}}=5.13^{+0.39}_{-0.49}\,\hbox {GeV}^{-1}$$\end{document}, g2ΣcΔD∗=-3.03-0.35+0.27GeV-2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g_{2\Sigma _{c}\Delta D^{*}}=-3.03^{+0.27}_{-0.35}\,\hbox {GeV}^{-2}$$\end{document}, g3ΣcΔD∗=17.64-1.95+1.51GeV-2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g_{3\Sigma _{c}\Delta D^{*}}=17.64^{+1.51}_{-1.95}\,\hbox {GeV}^{-2}$$\end{document}, g1ΣbΔB∗=20.97-2.39+2.15GeV-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g_{1\Sigma _{b}\Delta B^{*}}=20.97^{+2.15}_{-2.39}\,\hbox {GeV}^{-1}$$\end{document}, g2ΣbΔB∗=-11.42-1.28+1.17GeV-2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g_{2\Sigma _{b}\Delta B^{*}}=-11.42^{+1.17}_{-1.28}\,\hbox {GeV}^{-2}$$\end{document} and g3ΣbΔB∗=24.87-2.82+2.57GeV-2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g_{3\Sigma _{b}\Delta B^{*}}=24.87^{+2.57}_{-2.82}\,\hbox {GeV}^{-2}$$\end{document}. These strong coupling constants are important parameters which can help us to understand the strong decay behaviors of hadrons.


I. INTRODUCTION
The physics of charmed hadrons became an interesting subjects since the observations of J/ψ meson [1, 2]and charmed baryons (Λ c , Σ c ) [3].Up to now, lots of charmed baryons have been discovered by different experimental collaborations [4].Moreover, many bottom baryons such as Λ b , Ξ b , Σ b , Σ * b and Ω b have also been confirmed in experiments by CFD and LHCb collaborations [5][6][7][8][9][10].Although scientists have devoted much of their energy to this field, but the details of some charmed and bottom baryons are still less known.Thus, many experimental plans for the research of charmed and bottom baryons have been proposed by PANDA [11], J-PARC [12] and many other facilities.Under this circumstance, theoretical research on production of the baryons is very interesting and important.The strong coupling constants of baryons is an important input parameter which can help us to understand their production and decay processes [13].This is the first motivation for us to carry out the present work.
Since the observation of X(3872) by Belle collaboration in 2003 [14], exotic hadrons which are beyond the usual quarkmodel emerged like bamboo shoots after a spring rain [15][16][17][18][19][20][21][22][23][24][25][26][27][28].Some exotic states were interpreted as hadronic molecular states because their masses are close to the known twohadrons thresholds [29].However, the study of mass spectra is insufficient to understand the inner structure of these exotic states.We need to further study their strong decay behaviours, where the strong coupling constants are particularly important.For examples, in Ref [30], the authors predicted two pentaquark molecular states D * Σ c and D * Σ * c with the QCD sum rules.These two states were named as P c (4470) and P c (4620) which have the isospin I = 3  2 .If we studied their two-body strong decay P c (4470/4620) → J/ψ∆, this process can be described by the triangle diagram in  The strong interaction between the hadrons is nonperturbative in the low energy region, which can not be studied from the QCD first principle.But, as an important parameter, the strong coupling constant is urgently needed in studying the production and strong decay process of hadrons.Thus, some phenomenological methods are employed to analyze the strong vertices [31][32][33][34][35][36][37][38][39][40][41][42][43][44][45].The QCD sum rules (QCDSR) [46] and the light-cone sum rules (LCSR) are powerful phenomenological methods to study the strong interaction.In recent years, some coupling constants have been analyzed with LCSR by considering the higher-order QCD corrections and subleading power contributions [47,48].These studies show that considering the higher-order QCD corrections and subleading power contributions is very important for the accuracy of the results.In our previous work, we have analyzed the strong vertices Σ c ND, Σ b NB, Σ * c ND, Σ * b NB, Σ c ND * and Σ b NB * in the frame work of QCDSR basing on threepoint correlation function [38,40,41], where the higher-order perturbative corrections were neglected.As a continuation of these works, we analyze the strong vertices Σ c ∆D * and Σ b ∆B * using the three-point QCDSR under the tensor structure iǫ ρταβ p α p ′ β , p ρ p ′τ and p ρ p τ .According to our previous work, it showed that the subleading power contributions are really important for the final results.Considering higherorder corrections should make the final results more accurate, however it will also make the calculations of the three-point QCDSR very complicated.Thus, we neglect contributions from these corrections in the present work.
The layout of this paper is as follows.After the introduc-tion in Sec.I, the strong coupling constants of the vertices Σ c ∆D * and Σ b ∆B * are analyzed by QCD sum rules in Sec.II.In these analyses, the off-shell cases of the vector mesons are considered.In the QCD side, the perturbative contribution and vacuum condensate terms qq , g 2 s GG , qg s σGq and qq 2 are also considered.In Sec.III, we present the numerical results and discussions.Sec.IV is reserved for our conclusions.Some calculation details and important formulas are shown in Appendix A and B. The first step to analyze strong coupling constants with QCD sum rules is to write the following three-point correlation, where T is the time ordered product, and respectively.These interpolating currents can be ex-pressed as [49], where i, j and k represent the color indices and C denotes the charge conjugation operator.
The correlation function can be handled at both hadron and quark level in the framework of QCD sum rules, where the former is called the phenomenological side and the later is called the QCD side.Matching the calculation of these two sides by quark hadron duality, the sum rules for the strong coupling constants can be obtained.

A. The phenomenological side
In the phenomenological side, a complete sets of hadron states with the same quantum numbers as the hadronic interpolating currents are inserted into the correlation function.After isolating the contributions of ground and excited states, the expression of the correlation function can be written as [50], where h.c.denotes the contributions of higher resonances and continuum states.From this above equation, we can see that the current J µ ∆ (0) couples not only with the baryon J P = 3 2 + but also with the state of 1 2 + .Similarly, the meson current J ν D * [B * ] (0) couples with both the vector meson with J P = 1 − and the pseudoscalar meson with J P = 0 − .Therefore, there will be some redundant terms, that is the second, third and fourth term in Eq. ( 3)).They will disturb the items that we are interested in(the first term in Eq. ( 3)).These redundant matrix elements can be parameterized by the following equations, where N represents baryon with spin parity 1  4)), the projection operators (g µρ − p ′µ p ′ρ p ′2 ) and (g ντ − q ν q τ q 2 ) are employed in Eq. ( 3).The matrix elements about the vertex Σ c [Σ b ]∆D * [B * ] can be written as follows, where P = p + p ′ .The matrix elements appearing in Eq. ( 3) are substituted with Eqs. ( 4) and ( 5).Then, the correlation function in the phenomenological side can be written as the following form, From Eq. ( 6), we can see that the correlation function will have so complex tensor structure, e.g./ p / p ′ γ ρ γ τ γ 5 , / pγ ρ γ τ γ 5 , / p / p ′ g ρτ γ 5 , p ρ p τ γ 5 , • • • that the calculation become tedious and lengthy.Theoretically, if all the criteria of QCD sum rules are satisfied, each tensor structure can lead to the same results.For simplicity, we choose the tensor structure in the following ways, Πphy 1 , Πphy 2 and Πphy 3 are named as scalar invariant amplitudes which can be obtained by using Eq. ( 6) and ( 7), where, B. The QCD side In the QCD side, we firstly contract all of the quark fields in the correlation function with Wick's theorem, Here, S mn u[d] (x) and S mn c[b] (x) are light and heavy quark full propagators which can be written as [51,52], where 2 , λ a (a = 1, .., 8) are the Gell-Mann matrixes, i and j are color indices, σ αβ = i 2 [γ α , γ β ] and f λαβ , f αβµν have the following forms, Taking the same way as the phenomenological side, the correlation function in QCD side can also be written as, After conducting operator product expansion(OPE) and taking their imaginary part, we can obtain the spectral density of correlation function.Finally, the correlation function can be written as following form by using the dispersion relation, where s = p 2 , u = p ′2 , q = p − p ′ , s 1 and u 1 are the kinematic limits which are taken as (2m q + m Q ) 2 and 9m 2 q respectively.The QCD spectral density ρi (s, u, q 2 ) can be obtained by Cutkosky's rules [53][54][55][56][57][58], and their calculation details are briefly discussed in Appendix A. Full expressions of the QCD spectral density for different tensor structures are shown in Appendix B. The contributions of perturbative part and the vacuum condensation terms including qq , g 2 s GG , qg s σGq and qq 2 are all considered, where their Feynman diagrams are shown in Fig. 2.

C. The strong coupling constants
We take the change of variables p 2 → −P 2 , p ′2 → −P ′2 and q 2 → −Q 2 and perform double Borel transformation [59,60] to both the phenomenological and QCD sides.The variables P 2 and P ′2 are replaced by T 2 1 and T 2 2 which are called the Borel parameters.Then we take . Finally, we can obtain the following equations about the strong coupling constants g i (i = 1, 2, 3) using the quark-hadron duality condition, FIG. 2: Feynman diagrams for the perturbative part and vacuum condensate terms.
The momentum dependent coupling constants can be expressed as, where s 0 and u 0 are the threshold parameters which are introduced to eliminate the h.c.terms in Eq. ( 6).They satisfy the relations, , where m and m ′ are the masses of the ground and first excited states of the baryons.
In the framework of QCD sum rules, two conditions should also be satisfied, which are the pole dominance and convergence of OPE.To analyze the pole contribution, we write down, Then, the pole contribution can be defined as [50], The convergence of OPE is quantified via the contributions of the vacuum condensates of dimension n, which is defined as, where ρQCD i (s, u, Q 2 ) and ρQCD,n i (s, u, Q 2 ) represent the spectral densities of total and the nth dimension vacuum condensates, respectively.Fixing Q 2 = 3 GeV 2 in Eqs. ( 18) and ( 19), we plot the pole contributions with variation of the Borel parameter for different tensor structures in Fig. 3. To satisfy the convergence of OPE, we should also find a good plateau which is generally called 'Borel window'.Then, an appropriate Borel parameter in the Borel window is selected to make pole contributions larger then 40%.Considering these above requirements, the Borel windows are selected as 5(24) GeV 2 ≤ T 2 ≤ 7( 26 As for the gluon condensate g 2 s GG , it plays a less important role since |D(4)| < 1%.Therefore, the convergence of OPE is well satisfied.
By taking different values of Q 2 , we finally obtain the momentum dependent coupling constants g(Q 2 ) whose values are  shown in Fig. 6.In order to obtain the on-shell values of these coupling constants, it is necessary to extrapolate these results into the time-like regions (Q 2 < 0).This process is realized by fitting g(Q 2 ) with appropriate analytical functions and setting the vector meson ).To our knowledge, there are no specific expressions for the momentum dependent strong coupling constants which describe the interactions between hadrons.We only know that the value of running coupling constant α s (Q) decreases with the increment of square of momentum.Commonly, when we choose appropriate fitting functions, two conditions should be considered.The first is that the coupling constants should be well fitted by the fitting functions in the space-like regions (Q 2 > 0).Secondly, the on-shell values of the strong coupling constants, which are obtained by extrapolating the fitting functions into deep time-like regions, should converge.Based on our previous work, the combination of exponential and polynomial functions usually satisfies these conditions.In this work, the coupling constants of vertex Σ c ∆D * are well fitted by the combination of exponential and polynomial functions.For vertex of bottom baryon, the exponential function is not well convergent in Q 2 = −m 2 B * because the square mass of the vector bottom meson is much larger than that of charmed meson.Thus, the polynomial function is employed to fit the coupling constants of the vertex Σ b ∆B * .Finally, the momentum dependent strong coupling constants can be fitted into the following analytical functions, where a, b, c, d, e and f are the fitted parameters whose values are show in Tables I and II.The fitting curves for vertices Σ c ∆D * and Σ b ∆B * are also shown in Fig. 6.Finally, the on-shell values of strong coupling constants are obtained by setting 21),

Fig. 1 .
From this figure, we can see that analysis of strong vertices P c Σ c D * , P c Σ * c D * , DD * J/ψ, D * D * J/ψ, Σ c ∆D, Σ * c ∆D, Σ c ∆D * and Σ * c ∆D * is essential for us to study the strong decay behaviors of these two exotic states.This constituents the second motivation of our present work.
II.THE QCD SUM RULES FOR VERTICES Σ c ∆D * AND Σ b ∆B *

2 + 2 + and 3 2 +
, D[B] is the pseudoscalar charmed(bottom) meson, U(p, s) and U α (p, s) are the spinor wave functions of the baryon with spin parity 1 , respectively, ε β is the polarization vector of the vector meson D * [B * ], λ N is the pole residues, f D[B] is the decay constant.To extract the contributions of Σ c [Σ b ], D * [B * ] and ∆, and eliminate the contaminations of the redundant terms(see Eq. ( . Their values are m Σ c = 2.45 GeV, m Σ b = 5.81 GeV, m ∆ = 1.23 GeV, m D * = 2.01 GeV and m B * = 5.33 GeV, m u(d) = 0.006 ± 0.001 GeV, m c = 1.275 ± 0.025 GeV and m b = 4.18 ± 0.03

TABLE I :
The parameters of the analytical function for the coupling constants of vertex Σ c ∆D * .

TABLE II :
The parameters of the analytical function for the coupling constants of vertex Σ b ∆B * .