Magnetic dipole moments of the heavy tensor mesons in QCD

The magnetic dipole moments of the D2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal{D}_2$$\end{document}, and DS2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal{D}_{S_2}$$\end{document}, B2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal{B}_2$$\end{document}, and BS2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal{B}_{S_2}$$\end{document} heavy tensor mesons are estimated in framework of the light cone QCD sum rules. It is observed that the magnetic dipole moments for the charged mesons are larger than that of its neutral counterpart. It is found that the SU(3) flavor symmetry violation is about 10 % in both b and c sectors.

The masses and decay constants of the heavy tensor meson D * 2 (2740) and D * s 2 (2573) states were first studied within the framework of the QCD sum rules method in [13][14][15]. These mesons, as well as the B * 2 (5747) and B * s 2 (5840) tensor mesons have recently been studied within the same approach a e-mail: taliev@metu.edu.tr b e-mail: tbarakat@KSU.EDU.SA c e-mail: savci@metu.edu.tr in [15]. Note that light tensor mesons without, and with strange quark have also been analyzed in framework of the QCD sum rules method in [16,17], respectively.
One of the most promising approaches in investigating the properties of mesons and hadrons is the study of the electromagnetic form factors and multipole moments. These studies can provide useful information about their internal structures. The electromagnetic properties of usual mesons, as well as photons and neutrons have comprehensively been studied from theoretical and experimental sides, and at present it is the subject of the growing interest from both sides. However the study of the electromagnetic properties of the tensor mesons has received less interest, and therefore more effort is needed in this respect. The magnetic moments of the light tensor mesons were investigated in framework of the light cone QCD sum rules method (LCSR) in [18]. In the present work we calculate the magnetic dipole moments of the recently discovered heavy tensor mesons in LCSR.
The paper is organized as follows. In Sect. 2, the light cone QCD sum rules are constructed for the magnetic dipole moments of the heavy tensor mesons. In Sect. 3, the numerical analysis is performed for the obtained sum rules.

Theoretical framework
Before presenting the light cone sum rules for the magnetic dipole moments of the heavy tensor mesons, let us first introduce the matrix element which corresponds to the transition of the heavy tensor meson with momentum p + q to the heavy tensor meson with momentum p in presence of the electromagnetic field, i.e., where F i (Q 2 ) are the form factors, and T Q ( p, ε) means heavy tensor meson with momentum p and polarization tensor ε αβ . Since we are interested with the magnetic dipole moments of heavy tensor mesons, the values of the form factors at Q 2 = −q 2 = 0 need to be calculated. The transition under consideration can be described by the following correlation function, where j μν is the interpolating current of the heavy tensor meson, and j el ρ is the electromagnetic current given as, with the electric charges e q and e Q of the light and heavy quarks, respectively. The coupling of the tensor meson current j μν to the tensor state is defined as, where m T Q is the mass, and g T Q is the coupling constant of the tensor meson. It is more convenient to rewrite the correlator (2) by introducing the electromagnetic background field strength tensor of the plane wave, in the following form where F means that all vacuum expectation values are calculated in the background electromagnetic field. The correlation function (2) can be obtained by expanding the correlation function (4) in powers of the background field, and taking into account the terms that are linear in F μν which corresponds to the single photon radiation. The advantage of using the background field is that it separates the hard and soft photon radiation in an explicitly gauge invariant way (for more about the details of the background field method, see [19,20]). After these preliminary remarks, we can now proceed deriving the light cone QCD sum rules for the magnetic dipole moments of the heavy tensor mesons.
These sum rules can be obtained by calculating the correlator function in terms of mesons (physical part) from one side; and calculating the same correlation function in terms of quark-gluon degrees of freedom by using the operator product expansion (OPE) in deep Eucledian region from theoretical side. Matching these two representations of the same correlation function, the sum rules for the magnetic dipole moments of the heavy tensor mesons are obtained.
Calculation of the correlation function from the physical side is performed by inserting the complete set of tensor meson states having the same quantum number as that of the interpolating current, and isolating he ground state, as the result of which we get, where dots mean the contribution of the higher states and continuum.
From experimental point of view the multipole form factors are more useful than the form factors given in Eq. (1). The relations between the two sets of the form factors for any arbitrary q 2 are derived in [21]. For the real photon (q 2 = 0) these relations are given as: where G E (0) and G M (0) are the electric and magnetic multipoles. Calculation of the correlator function from the physical side is performed by inserting the complete set of tensor meson states in Eqs. (1)(2)(3)(4), from which for the hadronic part we get, Summation over the polarizations of the heavy tensor mesons can be performed by using the relation, where Using Eqs. (6) and (8) in Eq. (7), for the physical part of the correlation function we get, It follows from Eq. (9) that the correlation function μνραβ ε ρ contains various independent structures, each of which can be used in determination of the multipole moments of heavy tensor mesons. In the present work we restrict ourselves to the calculation of the magnetic dipole moment G M 1 , and for this goal we choose the structure (ε β q ν − ε ν q β )g μα . This structure has advantage over the others since it does not contain contributions coming from the contact terms (for more detail about the contact terms, see for example [22]). As the result we get the following sum rules for the magnetic dipole moment of the heavy tensor mesons, Denoting the coefficient of the (ε β q ν −ε ν q β )g μα structure as the invariant function , we get the following final result for the physical part of the correlator function for the magnetic dipole moment G M 1 , In order to construct the corresponding sum rules the correlation function μνραβ ε ρ needs to be calculated from the QCD side in terms of quark and gluon degrees of freedom using the operator product expansion, for which we need the interpolating current j μν . The interpolating current for the ground state heavy tensor meson with the quantum numbers 2 + can be chosen as, where q and Q denote the light and heavy quark fields, and the derivative operator ↔ Dμ with respect to x μ can be written as, which acts on the right and left sides, and the covariant derivatives in it are defined as where λ a are the Gell-Mann matrices, and A a μ (x) is the gluon field. In the present work we use the Fock-Schwinger gauge in which the external field A a μ (x) satisfies the condition x μ A a μ (x) = 0. The correlation function is calculated from the QCD side in deep Eucledian region p 2 → −∞, and ( p + q) 2 → −∞ after contracting the corresponding heavy and light quark fields, as a result of which we get, where we set y = 0 after performing derivative with respect to y. We see from Eq. (14) that, the expressions of the propagators of the light and heavy quarks are needed in order to calculate the correlation function. The light quark propagator in presence of the external field is calculated in [23], whose expression is given as: where is the scale parameter separating the perturbative and nonperturbative domains. This parameter is estimated in [24] to have the value = (0.5 ÷ 1.0) GeV; and S free (x − y) is the free quark operator whose expression is given as: The propagator for the heavy quark have the following form in coordinate space: where K i (m Q √ −x 2 ) are the modified Bessel functions, and Few words about the expression of the quark propagator are in order. The complete light cone expansion of the light quark propagator in presence of the external field is calculated in [23], which includes the contributions coming from nonlocal threeqGq, and four-particleqqqq,qG 2 q operators. Using the expansion in conformal spin, one can show that aforementioned contributions are small (for more detail see [25]), therefore we shall neglect them in further analysis.
There are three type of contributions to the correlation function: (1) perturbative part, when photon interacts perturbatively with the quark propagator (light or heavy). (2) "Mixed" contributions, which take place when heavy quark propagator interacts with the photon field perturbatively, and light quark fields form quark condensate. (3) "Long distance" contribution. It takes place when photon is radiated at long distance.
The perturbative contribution is calculated from Eq. (14) by replacing heavy or light propagator with, where ρ are the full set of Dirac matrices. As has already been noted, when a photon interacts with the light quark fields matrix elements of the nonlocal operators such as q(x) q(y) andq(x)G μν q(y) appear between vacuum and photon states. These matrix elements are parametrized in terms of the photon distribution amplitudes (DAs), which are the key nonperturbative parameters in light cone sum rules, whose explicit expressions are given below, where ϕ γ (u) is the leading twist-2, ψ v (u), ψ a (u), A and V are the twist-3, and h γ (u), A, T i (i = 1, 2, 3, 4) are the twist-4 photon DAs, and χ is the magnetic susceptibility. The photon DAs are calculated in [20] and their explicit expressions are given in Appendix A. The measure Dα i is defined as Separating out the coefficient of the structure (ε β q ν − ε ν q β )g μα from the QCD and the phenomenological parts of the correlation function and equating them, we get the magnetic moments of the heavy tensor mesons. In order to suppress the contributions of the higher states and continuum, we apply double Borel transformation with respect to the variables p 2 and ( p + q) 2 . As the result of these calculations we obtain the following sum rules for the magnetic dipole moment of the heavy tensor mesons, where The functions i n (n = 1, 2), and j 1 ( f (u)) are defined as: where s 0 is the continuum threshold.
Since we have the same heavy tensor mesons in the initial and final states, we can set M 2 1 = M 2 2 = 2M 2 , as the result of which we have, Physically this means that each quark and antiquark carries the half the photon momentum.

Numerical analysis
This section is devoted to the numerical analysis of the sum rules for the magnetic dipole moments of the heavy tensor mesons obtained in the previous section. The values of the input parameters entering into sum rules are, 2) GeV 2 which are obtained from the mass sum rule analysis for the light baryons [26,27], and B, B * mesons [28]. For the heavy quark masses we have used their M S values, which are given as:m b (m b ) = (4.16 ± 0.03) GeV andm c (m c ) = (1.28 ± 0.03) GeV [29][30][31][32][33][34][35]. The magnetic susceptibility of quarks was estimated in [36][37][38] in framework of the QCD sum rules. As we have noted earlier, the residues and masses of the heavy tensor mesons were calculated within the QCD sum rules method in [13][14][15], and their values are g D 2 = 0.0228 ± 0.0068, g D S 2 = 0.023 ± 0.011, g B 2 = 0.0050 ± 0.0005, g B S 2 = 0.0060 ± 0.0005.
Having decided the values of the input parameters, we are ready now to perform numerical analysis of the sum rules for the magnetic dipole moment of the heavy tensor mesons.The sum rule contains two unphysical parameters: a) The Borel mass parameter M 2 , and continuum threshold s 0 . It is known that the physical results should be independent of the these parameters. Therefore, our primary goal is to find such domain of these parameters for which the magnetic dipole moment is practically independent of them. The "working region" of M 2 is determined as follows: The upper bound of M 2 can be found by requiring that the contributions coming from higher states constitutes about 40 % of the perturbative part. The lower bound of M 2 could be fixed by demanding that the higher twist contributions are less than that of the leading twist contributions. In other words, the light cone expansion with increasing twist should be convergent. These requirements leads us to the following domains for the Borel mass parameter: