Strong coupling constants and radiative decays of the heavy tensor mesons

TL;DR

This paper calculates the strong coupling constants and radiative decay widths of heavy tensor mesons using three-point QCD sum rules, providing insights into their decay mechanisms and interactions.

Contribution

It introduces a comprehensive analysis of tensor-vector-pseudoscalar vertices and computes their coupling constants and decay widths within the QCD sum rules framework.

Findings

Calculated strong coupling constants for various heavy tensor mesons.

Determined radiative decay widths based on the coupling constants.

Provided analytical functions for strong form factors as a function of Q^2.

Abstract

In this article, we analyze tensor-vector-pseudoscalar(TVP) type of vertices $D_{2}^{*+}D^{+}\rho$, $D_{2}^{*0}D^{0}\rho$, $D_{2}^{*+}D^{+}\omega$, $D_{2}^{*0}D^{0}\omega$, $B_{2}^{*+}B^{+}\rho$, $B_{2}^{*0}B^{0}\rho$, $B_{2}^{*+}B^{+}\omega$, $B_{2}^{*0}B^{0}\omega$, $B_{s2}^{*}B_{s}\phi$ and $D_{s2}^{*}D_{s}\phi$, in the frame work of three point QCD sum rules. According to these analysis, we calculate their strong form factors which are used to fit into analytical functions of $Q^{2}$. Then, we obtain the strong coupling constants by extrapolating these strong form factors into deep time-like regions. As an application of this work, the coupling constants for radiative decays of these heavy tensor mesons are also calculated at the point of $Q^{2}=0$. With these coupling constants, we finally calculate the radiative decay widths of these tensor mesons.

Tables

Table 1. Table 1: Input parameters used in this analysis.

Parameters	Values($MeV$ )	Parameters	Values	Parameters	Values
$m_{D_2^{+\ast }}$	$2465.4\pm 1.3$ Tanabashi	$m_{B_{s2}^{\ast }}$	$5839.85\pm 0.17$ Tanabashi	$⟨\frac{{\alpha }_s{G}^2}{\pi }⟩$	$\left(0.012\pm 0.004\right)$ $Ge{V}^4$ Narison1 ; Narison2 ; Narison3
$m_{D_2^{0\ast }}$	$2460.7\pm 0.4$ Tanabashi	$m_{B_s}$	$5366.89\pm 0.19$ Tanabashi	$⟨g_s^{3}GGG⟩$	$\left(0.045\pm 0.002\right)Ge{V}^6$ Narison1 ; Narison2 ; Narison3
$m_{D^{+}}$	$1869.65\pm 0.05$ Tanabashi	$m_{\rho }$	$775.26\pm 0.25$ Tanabashi	$f_{\varphi }$	$\left(229\pm 3\right)MeV$ Tanabashi
$m_{D^0}$	$1864.83\pm 0.05$ Tanabashi	$m_{\omega }$	$782.65\pm 0.12$ Tanabashi	$f_{\rho }$	$\left(210\pm 4\right)MeV$ Tanabashi
$m_{B_2^{+\ast }}$	$5737.2\pm 0.7$ Tanabashi	$m_{\varphi }$	$1019.46\pm 0.016$ Tanabashi	$f_{\omega }$	$197\pm 8MeV$ Bharucha
$m_{B_2^{0\ast }}$	$5739.5\pm 0.7$ Tanabashi	$m_u$	${2.2}_{-0.4}^{+0.5}$ Tanabashi	$f_D$	$203.7\pm 4.7$ $MeV$ Tanabashi
$m_{B^{+}}$	$5279.32\pm 0.14$ Tanabashi	$m_d$	${4.7}_{-0.3}^{+0.5}$ Tanabashi	$f_{B_2^{\ast }}$	$110\pm 11$ $MeV$ Wang22
$m_{B^0}$	$5279.63\pm 0.15$ Tanabashi	$m_b$	${4180}_{-35}^{+25}$ Tanabashi	$f_{D_2^{\ast }}$	($182\pm 20$) $MeV$ Wang22
$m_{D_{s2}^{\ast }}$	$2569.1\pm 0.8$ Tanabashi	$m_c$	$1275\pm 25$ Tanabashi	$f_B$	$188\pm 25MeV$Tanabashi
$m_{D_s}$	$1968.3\pm 0.07$ Tanabashi	$m_s$	${95}_{-3}^{+9}$ Tanabashi	$f_{B_{s2}^{\ast }}$	$134\pm 11MeV$Wang22
$f_{D_{s2}^{\ast }}$	$222\pm 21$ Wang22	$f_{B_s}$	$231\pm 16$ Wang22	$f_{D_s}$	$257.8\pm 4.1MeV$ Tanabashi

Table 2. Table 2: Strong coupling constants and fitted parameters about the central value, upper bound and lower bound.

Mode	g($Ge{V}^{-2}$)	$A_c$	$B_c$	$A_u$	$B_u$	$A_l$	$B_l$
$D_2^{\ast +}$$\to$ $D^{+}\rho$	${4.89}_{-0.41}^{+0.80}$	3.062	-0.781	3.596	-0.764	2.659	-0.801
$D_2^{\ast 0}$$\to$ $D^0\rho$	${4.71}_{-0.41}^{+0.80}$	3.127	-0.684	3.698	-0.665	2.698	-0.706
$B_2^{\ast +}$$\to$ $B^{+}\rho$	${6.11}_{-0.97}^{+1.34}$	5.814	-0.083	7.169	-0.064	4.844	-0.100
$B_2^{\ast 0}$$\to$ $B^0\rho$	${6.03}_{-0.94}^{+1.30}$	5.733	-0.084	7.049	-0.065	4.785	-0.102
$D_2^{\ast 0}$$\to$ $D^0\omega$	${5.07}_{-0.74}^{+1.00}$	3.333	-0.686	4.034	-0.667	2.810	-0.707
$D_2^{\ast +}$$\to$ $D^{+}\omega$	${5.28}_{-0.75}^{+0.98}$	3.262	-0.785	3.915	-0.757	2.769	-0.804
$B_2^{\ast +}$$\to$ $B^{+}\omega$	${6.53}_{-0.74}^{+0.88}$	6.196	-0.084	7.107	-0.068	5.426	-0.105
$B_2^{\ast 0}$$\to$ $B^0\omega$	${6.44}_{-0.69}^{+0.91}$	6.105	-0.086	7.056	-0.066	5.397	-0.104
$B_{s2}^{\ast }$$\to$ $B_s\varphi$	${3.60}_{-0.42}^{+0.54}$	3.162	-0.126	3.702	-0.107	2.741	-0.143
$D_{s2}^{\ast }$$\to$ $D_s\varphi$	${6.93}_{-0.02}^{+0.40}$	2.102	-1.149	2.254	-1.136	1.907	-1.241

Table 3. Table 3: The decay widths for different radiative decays.

Radiative decay	$\Gamma$($keV$)
$D_2^{\ast +}$$\to$ $D^{+}\gamma$	${0.464}_{-0.079}^{+0.126}$
$D_2^{\ast 0}$$\to$ $D^0\gamma$	${2.730}_{-0.491}^{+0.776}$
$D_{s2}^{\ast }$$\to$ $D_s\gamma$	${1.658}_{-0.010}^{+0.198}$
$B_2^{\ast +}$$\to$ $B^{+}\gamma$	${1108}_{-202}^{+315}$
$B_2^{\ast 0}$$\to$ $B^0\gamma$	${275}_{-48.0}^{+78.1}$
$B_{s2}^{\ast }$$\to$ $B_s\gamma$	${37.6}_{-8.2}^{+12.2}$

Equations

$$\displaystyle\Pi_{\mu\nu\tau}(p,p^{\prime})=i^{2}\int d^{4}x\int d^{4}ye^{i(p-p^{\prime}).x+ip^{\prime}.y}\Big{\langle}0|\mathcal{T}\Big{(}J_{\tau}(x)J_{\mathbb{P}}(y)J^{\dagger}_{\mu\nu}(0)\Big{)}|0\Big{\rangle},$$

$$\displaystyle J_{\mu\nu}(z)=i\overline{Q}(z)\Big{(}\gamma_{\mu}\overleftrightarrow{D}_{\nu}+\gamma_{\nu}\overleftrightarrow{D}_{\mu}-\frac{2}{3}\widetilde{g}_{\mu\nu}\overleftrightarrow{D\!\!\!\!/}\Big{)}q(z)$$

$$J_{\mathbb{P}}(y)=\overline{Q}(y)i{\gamma }_{5}q(y)$$

$$J_{\tau }(x)=\overline{q}(x)i{\gamma }_{\tau }q(x)$$

$${\overleftrightarrow{D}}_{\mu }=({\overrightarrow{\partial }}_{\mu }-i{g}_s{G}_{\mu })-({\overleftarrow{\partial }}_{\mu }+i{g}_s{G}_{\mu })$$

$${\tilde{g}}_{\mu \nu }=g_{\mu \nu }-\frac{p_{\mu }{p}_{\nu }}{p^2}$$

$${\Pi }_{\mu \nu \tau }^{hadr}(p,p^{\prime },q)=$$

$$⟨\mathbb{T}(p)|J_{\mu \nu }^{†}(0)|0⟩=f_{\mathbb{T}}{M}_{\mathbb{T}}^2{\xi }_{\mu \nu }^{\ast }(p),$$

$$⟨0|J_{\tau }(0)|\mathbb{V}(q)⟩=f_{\tau }{M}_{\tau }{\zeta }_{\tau }(q),$$

$$⟨0|J_{\mathbb{P}}(0)|\mathbb{P}(p^{\prime })⟩=\frac{f_{\mathbb{P}}{M}_{\mathbb{P}}^2}{m_Q+m_q},$$

$$⟨\mathbb{P}(p^{\prime })\mathbb{V}(q)|L_{\mathbb{VPT}}|\mathbb{T}(p)⟩=g{\epsilon }^{\alpha \beta \lambda \rho }{p}_{\alpha }{\xi }_{\beta \eta }{p}^{\prime \eta }{q}_{\lambda }{\zeta }_{\rho }^{\ast }.$$

$$\displaystyle\xi^{*}_{\mu\nu}\xi_{\beta\eta}=\frac{1}{2}\big{[}\widetilde{g}_{\mu\beta}\widetilde{g}_{\nu\eta}+\widetilde{g}_{\mu\eta}\widetilde{g}_{\nu\beta}\big{]}-\frac{1}{3}\widetilde{g}_{\mu\nu}\widetilde{g}_{\beta\eta}$$

$${\zeta }_{\rho }^{\ast }{\zeta }_{\tau }=g_{\rho \tau }-\frac{p_{\rho }{p}_{\tau }}{p^2}$$

$${\Pi }_{\mu \nu \tau }^{hadr}(p,p^{\prime },q)=\frac{g{f}_{\mathbb{P}}{M}_{\mathbb{P}}^2{f}_{\tau }{M}_{\tau }{f}_{\mathbb{T}}{M}_{\mathbb{T}}^2}{(m_Q+m_q)(M_{\tau }^2-q^2)(M_{\mathbb{P}}^2-p^{\prime 2})(M_{\mathbb{T}}^2-p^2)}$$

$${\Pi }_{\mu \nu \tau }^{OPE}(p,p^{\prime },q)=-i^4\int d^{4}x\int d^{4}y{e}^{i(p-p^{\prime }).x+i{p}^{\prime }.y}tr[S_{nm}^q(y-x){\gamma }_{\tau }{S}_{mk}^q(x-z){\Gamma }_{\mu \nu }{S}_{kn}^Q(z-y){\gamma }_5{|}_{z=0}]$$

$${\Gamma }_{\mu \nu }={\gamma }_{\mu }{\overleftrightarrow{D}}_{\nu }+{\gamma }_{\nu }{\overleftrightarrow{D}}_{\mu }-\frac{2}{3}{\tilde{g}}_{\mu \nu }\overleftrightarrow{D/}$$

$$S_{nm}^q(x)=$$

$$S_{kn}^Q(x)$$

$$f^{\lambda \alpha \beta }=(k/+m_Q){\gamma }^{\lambda }(k/+m_Q){\gamma }^{\alpha }(k/+m_Q){\gamma }^{\beta }(k/+m_Q)$$

$$f^{\alpha \beta \mu \nu }=(k/+m_Q){\gamma }^{\alpha }(k/+m_Q){\gamma }^{\beta }(k/+m_Q){\gamma }^{\mu }(k/+m_Q){\gamma }^{\nu }(k/+m_Q)$$

$${\Pi }_{\mu \nu \tau }^{pert}(p,p^{\prime },q)$$

$$\displaystyle\Gamma_{\mu\nu}=\gamma_{\mu}(2k_{\nu}-2p^{\prime}_{\nu}+p_{\nu})+\gamma_{\nu}(2k_{\mu}-2p^{\prime}_{\mu}+p_{\mu})-\frac{2}{3}\Big{[}(g_{\mu\nu}-\frac{p_{\mu}p_{\nu}}{p^{2}})(2k\!\!\!/-2p^{\prime}\!\!\!\!\!/+p\!\!\!/)\Big{]}$$

$${\rho }_{\mu \nu \tau }^{pert}(s,u,q^2)=$$

$$A=\frac{(u+m_q^2-m_Q^2)(s+u-q^2)-2u(u-q^2+m_q^2-m_Q^2)}{\lambda (s,u,q^2)}$$

$$B=\frac{(u-q^2+m_q^2-m_Q^2)(s+u-q^2)-2s(u+m_q^2-m_Q^2)}{\lambda (s,u,q^2)}$$

$$\lambda (s,u,q^2)=(s+u-q^2{)}^2-4su$$

$${\Pi }_{\mu \nu \tau }^{pert}(p,p^{\prime })={\int }_{s_1}^{s_0}{\int }_{u_1}^{u_0}\frac{{\rho }_{\mu \nu \tau }^{pert}(s,u,q^2)}{(s-p^2)(u-p^{\prime 2})}dsdu{|}_{-1\le \frac{2s(m_q^2-m_Q^2+u)+(m_Q^2-m_q^2-u+q^2)(s+u-q^2)}{\sqrt{(m_Q^2-m_q^2-u+q^2{)}^2-4s{m}_q^2}\sqrt{\lambda (s,u,q^2)}}\le 1}$$

$${\Pi }_{\mu \nu \tau }^{⟨GG⟩}=$$

$${\Pi }_{\mu \nu \tau }^{⟨GGG⟩}=$$

$$I_0^{abc}$$

$$I_{\alpha 01}^{abc}$$

$$I_{\alpha 10}^{abc}$$

$$\displaystyle\frac{gf_{\mathbb{P}}M^{2}_{\mathbb{P}}f_{\tau}M_{\tau}f_{\mathbb{T}}M_{\mathbb{T}}^{2}}{(m_{q}+m_{Q})(M_{\tau}^{2}+Q^{2})}\Big{(}-\frac{M_{\mathbb{T}}^{2}+M_{\mathbb{P}}^{2}-Q^{2}}{4M_{\mathbb{T}}^{2}}\Big{)}\frac{1}{M_{1}^{2}M_{2}^{2}}exp[-\frac{M_{\mathbb{T}}^{2}}{M_{1}^{2}}]exp[-\frac{M_{\mathbb{P}}^{2}}{M_{2}^{2}}]$$

$$=$$

$$g(Q^2)=Aexp[B{Q}^2]$$

$$\Gamma =\frac{1}{2J+1}\sum \frac{|p|}{8\pi M_i^2}|T{|}^2$$

Full text

**1 Introduction **

With rapid developments of high-energy physics experiments, more and more new states of mesons have been confirmed by D0, CDF and LHCb collaborationsAbazov1 ; Aaltonen1 ; Abazov2 ; Aaltonen2 ; Beringer ; Aaij ; Aaltonen3 . The heavy-light mesons, which are composed of a heavy quark and a light quark, can be classified into the spin doublets in the heavy quark limit. For example, the $1S(0^{-},1^{-})$ doublets $(B,B^{*})$, $(D,D^{*})$, $(B_{s},B_{s}^{*})$, $(D_{s},D_{s}^{*})$ and the $1P(1^{+},2^{+})$ doublets $(B_{1},B_{2}^{*})$, $(B_{s1},B_{s2}^{*})$, $(D_{1},D_{2}^{*})$, $(D_{s1},D_{s2}^{*})$ have also been confirmed in experimentsBeringer . That is to say, the quantum numbers $I(J^{P})$ for heavy tensor mesons $D_{2}^{*}$, $B_{2}^{*}$, $D_{s2}^{*}$ and $B_{s2}^{*}$ are $\frac{1}{2}(2^{+})$, $\frac{1}{2}(2^{+})$, $0(2^{+})$ and $0(2^{+})$ respectively.

Compared with the $1S(0^{-},1^{-})$ and $1P(0^{+},1^{+})$ states of the heavy mesons, the $1P(1^{+},2^{+})$ doublets have been drawn little attentionSwanson ; Klempt . The strong decay processes $D_{2}^{*}\rightarrow D^{*}\pi$, $D\pi$Abazov1 ; Aubert ; Sanchez ; Aaij2 , $D_{s2}^{*}\rightarrow DK$Abazov1 , $B_{2}^{*}\rightarrow B^{*}\pi$, $B\pi$Abazov1 ; Aaltonen2 , $B_{s2}^{*}\rightarrow BK,B^{*}K$Aaltonen1 ; Abazov2 ; Aaij have been observed in experiments. In our previous work, we have studied these strong decay processes, obtained their strong coupling constants and strong decay widthsWang1 ; Li ; GuoLiang . As a continuation of these work, we study the strong vertices $D_{2}^{*+}D^{+}\rho$, $D_{2}^{*0}D^{0}\rho$, $D_{2}^{*+}D^{+}\omega$, $D_{2}^{*0}D^{0}\omega$, $B_{2}^{*+}B^{+}\rho$, $B_{2}^{*0}B^{0}\rho$, $B_{2}^{*+}B^{+}\omega$, $B_{2}^{*0}B^{0}\omega$, $B_{s2}^{*}B_{s}\phi$ and $D_{s2}^{*}D_{s}\phi$ and obtain its strong coupling constants. These strong coupling constants not only play an essential role for understanding the inner structure of these mesons but also can help us to know about its decay behaviors. Besides, the strong coupling constants about the heavy-light mesons can also help us understanding the final-state interactions in the heavy quarkonium (or meson) decaysCasalbuoni ; LiuX23 ; GuoFK . With a fitted function about the strong form factors in Section 3, we can also obtain the coupling constants for the radiative decays with intermediate momentum $Q^{2}=0$, which will be used to calculate the radiative decay widths of these mesons.

To study the decay behaviors of the mesons, we can adopt several theoretical models including perturbative and non-perturbative methods. The QCD sum rules, proposed by Shifman, Vainshtein, and ZakharovShifman , connects hadron properties and QCD parametersReind . It has been widely used to study the properties of the hadronsIoffe ; Belyaev ; Brac1 ; Brac2 ; Alie3 ; Alie4 ; Doi ; Altm ; Wzg3 ; Cerq ; Rodr ; Yazi ; Khos1 ; Khos2 ; Rein ; Pasc ; Wzg5 ; Wang2 ; Navar ; Khodj ; GuoLY2 ; GuoLY4 ; Azizi4 ; Azi1 ; Azi2 ; LiuX ; HeJ ; LiuX2 ; ChenHX1 ; ChenHX2 ; Chun ; ChenW ; ChenW2 ; Zhjr1 ; Zhjr2 ; Zhjr3 ; Wzg . In this work, we analyze the tensor-vector-pseudoscalar(TVP) type of vertices and the radiative decays using the three-point QCD sum rules. This paper is organized as follows. After the Introduction, we study the tensor-vector-pseudoscalar(TVP) type of strong vertices using the three point QCD sum rules with vector mesons being off-shell. In Sec.3, we present the numerical results and discussions. Finally, the paper ends with the conclusions.

**2 QCD sum rules for hadronic coupling constants **

For tensor-vector-pseudoscalar(TVP) type of vertices, its three-point correlation function is written as,

[TABLE]

where $J_{\mu\nu}$, $J_{\tau}$, and $J_{\mathbb{P}}$ denote interpolating currents of heavy tensor mesons, vector mesons and pseudoscalar mesons. These interpolating currents have the same quantum numbers with studied mesonsIoffe ; Cola ,

[TABLE]

where

[TABLE]

**2.1 The hadronic side **

To obtain hadronic representation, we insert a complete set of intermediate hadronic states into the correlation $\Pi_{\mu\nu\tau}(p,p^{\prime})$. These intermediate states have the same quantum numbers with the current operators $J_{\mu\nu}$, $J_{\tau}$, and $J_{\mathbb{P}}$. After isolating ground-state contributions of these mesonsShifman ; Reind , the correlation function is expressed as,

[TABLE]

The matrix elements appearing in this equation are substituted with the following parameterized equations,

[TABLE]

with $q=p-p^{\prime}$. Here, $f_{\mathbb{T}}$, $f_{\mathbb{P}}$ and $f_{\tau}$ are decay constants of the tensor mesons, pseudoscalar mesons and vector mesons, and $g$ is the strong form factor of tensor-vector-pseudoscalar(TVP) type of vertices. Besides, $\xi_{\mu\nu}$, $\zeta_{\rho}$ are polarization vectors of the tensor mesons and vector mesons with the following properties,

[TABLE]

With these above equations, the correlation function $\Pi_{\mu\nu\tau}(p,p^{\prime},q)$ can be expressed as follows,

[TABLE]

**2.2 The OPE side **

In this part, we will briefly outline the operator product expansion(OPE) for the correlation function $\Pi_{\mu\nu\tau}(p,p^{\prime},q)$ in perturbative QCD. Firstly, we contract all of the quark fields with Wick’s theorem, and rewrite the correlation function as follows,

[TABLE]

where

[TABLE]

and $S^{q}$($S^{Q}$) denote light(heavy) quark propagators which can be expressed asPasc ; Rein .

[TABLE]

[TABLE]

[TABLE]

where $t^{a}=\frac{\lambda^{a}}{2}$, the $\lambda^{a}$ is the Gell-Mann matrix, and $n$, $m$, $k$ are color indicesReind . In the covariant derivative, the gluon $G_{\mu}(z)$ in Eq.(4) has no contributions as $G_{\mu}(z)=\frac{1}{2}z^{\lambda}G_{\lambda\mu}(0)+\cdots=0$. Using equations (4),(5), (6) and (7), the perturbative contribution of the correlation function is written as

[TABLE]

where

[TABLE]

Putting all the quark lines on mass-shell by the Cutkosky s rules, we compute the integrals both in coordinate and momentum spaces. Then, we can obtain the spectral density by taking the imaginary parts of the correlation function,

[TABLE]

where

[TABLE]

During these derivations, we set $s=p^{2}$, $u=p^{\prime 2}$ and $q=p-p^{\prime}$ in the spectral densities. As a result, we can see that there are several different structures on hadronic side and OPE side. In general, we can choose either structure to study the hadronic coupling constant. In our calculations, we observe that the structure $\varepsilon^{\nu\tau pp^{\prime}}p_{\mu}$ can lead to pertinent result. Using dispersion relation, the perturbative term can be written as,

[TABLE]

For non-perturbative terms, we take into account the contributions of $\langle q\overline{q}\rangle$, $\langle\overline{q}g\sigma.Gq\rangle$, $\langle g^{2}G^{2}\rangle$ and $\langle f^{3}G^{3}\rangle$. After performing double Borel transformation, we find that contributions of non-perturbative terms come only from condensate terms $\langle g^{2}G^{2}\rangle$, $\langle f^{3}G^{3}\rangle$. The expressions of these condensate terms are written as,

[TABLE]

[TABLE]

where

[TABLE]

[TABLE]

[TABLE]

**3 The results and discussions **

Estimating the parameters of the lowest-lying hadronic state are in general plagued by the presence of unknown subtraction terms, the spectral function of excited and continuum states. This situation can be substantially improved by applying to both OPE side and phenomenological side the Borel transformationCola . Thus, we perform the double Borel transform with respect to the variables $P^{2}=-p^{2}$, $P^{\prime 2}=-p^{\prime 2}$ and match OPE side with the hadronic representation Eq.(3), invoking the quark-hadron duality. Finally, we obtain the QCDSR as follows,

[TABLE]

Here, $Q^{2}=-q^{2}$, parameters $s_{0}$ and $u_{0}$ are used to further reduce the contributions from excited and continuum states. Its values are employed as $s_{0}=(m_{i}+\Delta_{i})^{2}$ and $u_{0}=(m_{o}+\Delta_{o})^{2}$, where $m_{i}$ and $m_{o}$ are ground state masses of the in-coming and out-coming hadron. In general, $\Delta_{i}$ and $\Delta_{o}$ are expected to be $0.3GeV\sim 0.5GeV$, which can guarantee the values of $s_{0}$ and $u_{0}$ be close to the mass squared of the first excited state of these in-coming and out-coming hadronsWang1 . Parameters $M_{1}^{2}$ and $M_{2}^{2}$ in Eq.(14) are Borel parameters. In order to determine optimal values about these above parameters, two criteria should be considered. First, pole contribution should be as large as possible comparing with contributions of higher and continuum states. Secondly, we should also ensure OPE convergence and the stability of our results. That is to say, the results which are extracted from sum rules, should be independent of the Borel parameters. One can consult Ref.Li ; GuoLiang for technical details of these processes. As for the other parameters in Eq.(14), their values are all listed in Table 1.

The strong form factor $g$ from Eq.(14) are obtained in deep space-like region $q^{2}\rightarrow-\infty$, where the intermediate mesons are off-shell. In order to obtain strong coupling constants, we must extrapolate these results into deep time-like region. It is noticed that there are no exact expressions for the dependence of the strong form factors on $Q^{2}$. In this work, we find that the results can be fitted into the exponential function,

[TABLE]

In Figs.1-10, we show the values of the strong form factors on $Q^{2}$ which are obtained from Eq.(14) and the fitting curve, in which it is marked as Central value and Fitted curve of Central value. Thus, we can obtain the strong coupling constants by taking $Q^{2}=-m_{on-shell}^{2}$ for intermediate mesons in the fitting function Eq.(15). The values of fitted parameters $A$ and $B$ in Eq.(15) and the strong coupling constants are all listed in Table II.

The uncertainties of strong form factors in Eq.(14) mainly come from input parameters $m_{D_{2}^{*+}}$,$m_{D_{2}^{0*}}$,$m_{u}$,$m_{b}$,$f_{\rho}$, $f_{\phi}$, $\langle\overline{q}q\rangle$, $\cdots$ Theoretically, we can calculate its values with uncertainty transfer formula $\delta=\sqrt{\Sigma_{i}(\frac{\partial f}{\partial x_{i}})^{2}(x_{i}-\overline{x}_{i})^{2}}$, where $f$ denotes the strong form factor in Eq.(14), and $x_{i}$ denotes input parameters. For simplicity, the upper and lower limits of the results are estimated by taking $f^{upper(lower)}=f(\overline{x}_{i}\pm\Delta x_{i})$, which are marked as Upper bound and Lower bound in Figs.1-10. After these approximations, they are also fitted into exponential functions and are also extrapolated into the physical regions in order to get the uncertainties of the strong coupling constants. These results are all listed in Table II.

Finally, we give an analysis of the radiative decays of the heavy tensor mesons $\mathbb{T}\rightarrow\mathbb{P}\gamma$. The coupling constants of these radiative decays $g_{\mathbb{T}\mathbb{P}\gamma}$ can be easily obtained by setting $Q^{2}=0$ in Eq.(15). The radiative decay width can be expressed as the following representation,

[TABLE]

[TABLE]

where $i$ and $f$ denote the initial and final state mesons, $J$ is the total angular momentum of the initial meson, $\sum$ denotes the summation of all the polarization vectors, and $T$ denotes the scattering amplitudes. The radiative decays $\mathbb{T}\rightarrow\mathbb{P}\gamma$ can be described by the following electromagnetic lagrangian $\pounds$

[TABLE]

From this lagrangian, the decay amplitude can be written as,

[TABLE]

Here, $p_{\alpha}$, $p^{\prime\eta}$ and $q_{\lambda}$ are the four momenta of the tensor meson, pseudoscalar meson and $\gamma$, $\xi$, $\zeta$ and $\varepsilon$ are their polarization vectors, respectively. With Eqs.(16) and (17), we can obtain the radiative decay width of $\mathbb{T}\rightarrow\mathbb{P}\gamma$,

[TABLE]

where $\alpha=\frac{1}{137}$,$Q_{u}=\frac{2}{3}$,$Q_{d}=Q_{s}=-\frac{1}{3}$. Considering different decay channels, we obtain the widths of different radiative decays which are listed in Table III. From referenceTanabashi , we can see the decay widths of the tensor mesons, $\Gamma(D_{2}^{*0})=47.5\pm 1.1MeV$, $\Gamma(D_{2}^{*\pm})=46.7\pm 1.2MeV$, $\Gamma(D_{s2}^{*})=16.9\pm 0.8MeV$, $\Gamma(B_{2}^{*0})=24.2\pm 1.7MeV$, $\Gamma(B_{2}^{*+})=20\pm 5MeV$, $\Gamma(B_{s2}^{*})=1.47\pm 0.33MeV$. From these experimental data, we observe that the branching ratios of the calculated radiative decays are of the order of $10^{-2}\sim 10^{-5}$, which are measurable in the future by LHCb. In referenceAliev10 , the radiative decays of the heavy tensor mesons were also analyzed in the framework of the light cone QCD sum rules method. We observe that our results for mesons $D_{2}^{*}$ and $D_{s2}^{*}$ are comparable with its results. For mesons $B_{2}^{*}$ and $B_{s2}^{*}$, the results from QCD sum rules and light cone QCD sum rules vary widely, which need to be further studied by other theoretical methods or in experiments.

**4 Conclusion **

In this paper, we analyze the tensor-vector-pseudoscalar(TVP) type of vertices in the cases of light vector mesons $\rho$, $\omega$ and $\phi$ being off-shell. We firstly calculate its strong form factors in space-like regions($q^{2}<0$). Then, we fit the form factors into exponential functions which are used to extrapolate into time-like regions($q^{2}>0$) to obtain strong coupling constants. These strong coupling constants are important parameters in studying the strong decay behaviors of the tensor mesons in the future. Setting intermediate momentum $Q^{2}=0$ in the fitted analytical functions about strong form factors, we also obtained the coupling constants of the radiative decays of the tensor mesons. With these coupling constants, we calculate the radiative decay widths of these tensor mesons and compare our results with experimental data and those of other research groups.

**Acknowledgment **

This work has been supported by the Fundamental Research Funds for the Central Universities, Grant Number $2016MS133$, Natural Science Foundation of HeBei Province, Grant Number $A2018502124$.