Advances in High Energy Physics

Volume 2018, Article ID 3637824, 9 pages

https://doi.org/10.1155/2018/3637824

## Reanalysis of Decay in QCD

Department of Physics, Middle East Technical University, 06800 Ankara, Turkey

Correspondence should be addressed to T. M. Aliev; rt.ude.utem@veilat

Received 9 January 2018; Revised 28 March 2018; Accepted 4 April 2018; Published 10 May 2018

Academic Editor: Enrico Lunghi

Copyright © 2018 T. M. Aliev et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The publication of this article was funded by SCOAP^{3}.

#### Abstract

The strong coupling constants of newly observed baryons with spins and decaying into are estimated within light cone QCD sum rules. The calculations are performed within two different scenarios on quantum numbers of baryons: (a) all newly observed baryons are negative parity baryons; that is, the , , , and have quantum numbers and states, respectively; (b) the states and have quantum numbers and , while the states and have the quantum numbers and , respectively. By using the obtained results on the coupling constants, we calculate the decay widths of the corresponding decay. The results on decay widths are compared with the experimental data of LHC Collaboration. We found out that the predictions on decay widths within these scenarios are considerably different from the experimental data; that is, both considered scenarios are ruled out.

#### 1. Introduction

Lately, in the invariant mass spectrum of , very narrow excited states (, , , , and ) have been observed at LHCb [1]. Quantum numbers of these newly observed states have not been determined in the experiments yet. Hence, various possibilities about the quantum numbers of these states have been speculated in recent works. In [2], the states and are assigned as radial excitation of ground state and baryons with the and , respectively. On the other hand, in [3–7] these new states are assumed as the -wave states with , and , respectively. Moreover the new states are assumed as pentaquarks in [8]. Similar quantum numbers of these new states are assigned in [9]. Analysis of these states is also studied with lattice QCD, and the results indicated that most probably these states have quantum numbers [3]. Another set of quantum number assignments, namely, , , , and , is given in [4]. In [10], it is obtained that the prediction on mass supports assigning as , as or the state with , and as .

In this work, we estimate the strong coupling constants of in the framework of light cone QCD sum rules. In our calculations, two different possibilities on quantum numbers of baryons are explored:(a)All newly observed states have negative parities. More precisely, and have quantum numbers , , and have , and has quantum numbers .(b)Part of newly observed baryons have negative parity, and another part represents a radial excitations of ground state baryons; that is, has , has , and and states have quantum numbers and , respectively.

Note that the strong coupling constants of decay within the same framework are studied in [11] and in chiral quark model [12], respectively. However, the analysis performed in [11] is incomplete. First of all, the contribution of negative parity baryons is neglected entirely. Second, in our opinion the numerical analysis presented in [11] is inconsistent.

The article is organized as follows. In Section 2 the light cone sum rules for the coupling constants of decay are derived. Section 3 is devoted to the analysis of the sum rules obtained in the previous section. In this section, we also estimate the widths of corresponding decay, and comparison with the experimental data is presented.

#### 2. Light Cone Sum Rules for the Strong Coupling Constants of Transitions

For the calculation of the strong coupling constants of transitions we consider the following two correlation functions in both pictures:where is the interpolating current of baryon and is the interpolating current of baryon:where , and are color indices, is the charge conjugation operator, and is arbitrary parameter.

We calculate () employing the light cone QCD sum rules (LCSR). According to the sum rules method approach, the correlation functions in (1) can be calculated in two different ways:(i)In terms of hadron parameters.(ii)In terms of quark-gluons in the deep Euclidean domain.

These two representations are then equated by using the dispersion relation, and we get the desired sum rules for corresponding strong coupling constant. The hadronic representations of the correlation functions can be obtained by saturating (2) and (3) with corresponding baryons.

Here we would like to note that the currents , , and interact with both positive and negative parity baryons. Using this fact for the correlation functions from hadronic part we getThe matrix elements in (5) are determined aswhere is strong coupling constant of the corresponding decay, are the residues of the corresponding baryons, and is the Rarita-Schwinger spinors. Here the sign corresponds to positive (negative) parity baryon. In further discussions, we will denote the mass and residues of ground and excited states of baryons as , , and for scenario (a) and for scenario (b); the same notation is used as in previous case by just replacing to . Moreover, the mass and residues of baryons are denoted as , , and . Using the matrix elements defined in (6) for the correlation functions given in (1) we get (for case (a))whereThe result for scenario (b) can be obtained from (8) and (9) by following replacements:

Note that to derive (9), we used the following formula for performing summation over spins of Rarita-Schwinger spinors:and in principle one can obtain the expression for the hadronic part of the correlation function. At this stage two problems arise. One of them is dictated by the fact that the current interacts not only with spin 3/2 but also with spin 1/2 states. The matrix element of the current with spin state is defined asthat is, the terms in the RHS of (12) and the right end contain contributions from states, which should be removed. The second problem is related to the fact that not all structures appearing in (9) are independent. In order to cure both these problems we need ordering procedure of Dirac matrices. In present work, we use ordering of Dirac matrices as . Under this ordering, only the term contains contributions solely from spin states. For this reason, we will retain only terms in the RHS of (9).

In order to find sum rules for the strong coupling constants of transitions we need to calculate and from QCD side in the deep Euclidean region, , . The correlation from QCD side can be calculated by using the operator product expansion.

Now let us demonstrate steps of calculation of the correlation function from QCD side. As an example let us consider one term of correlation ; that is, consider

By using Wick’s theorem, this term can be written asFrom this formula, it follows that, to obtain the correlation function(s) from QCD side, first of all we need the expressions of light and heavy quark propagators. The expressions of the light quark propagator in the presence of gluonic and electromagnetic background fields are derived in [13]

The heavy quark propagator is given aswhere is the Euler constant.

For calculation of the correlator function(s) we need another ingredient of light cone sum rules, namely, the matrix elements of nonlocal operators and between vacuum and the -meson, that is, and . Here is the any Dirac matrix, and is the gluon field strength tensor, respectively. These matrix elements are defined in terms of -meson distribution amplitudes (DAs). The DAs of meson up to twist- are presented in [12].

From (8) and (9) it follows that the different Lorentz structures can be used for construction of the relevant sum rules. Among of six couplings, we need only and and and for cases (a) and (b), respectively. For determination of these coupling constants, we need to combine sum rules obtained from different Lorentz structures. From (8) and (9) (for case (a)) it follows that the Lorentz structures , , , and , , , and appear. We denote the corresponding invariant functions and , and , respectively. Explicit expressions of the invariant functions and are very lengthy, and therefore we do not present them in the present study.

The sum rules for the corresponding strong coupling constants are obtained by choosing the coefficients aforementioned structures and equating to the corresponding results from hadronic and QCD sides. Performing doubly Borel transformation with respect to variable and in order to suppress the contributions of higher states and continuum we get the following four equations (for each transition):where superscript means Borel transformed quantities,

The masses of the initial and final baryons are close to each other; hence in the next discussions, we set . In order to suppress the contributions of higher states and continuum we need subtraction procedure. It can be performed by using quark-hadron duality; that is, starting some threshold the spectral density of continuum coincides with spectral density of perturbative contribution. The continuum subtraction can be done using formulaFor more details about continuum subtraction in light cone sum rules, we refer readers to work [14].

As we have already noted in case (a) we need to determine two coupling constants () and () for each class of transitions. From (18) it follows that we have six unknown coupling constants but have only four equations. Two extra equations can be obtained by performing derivative over () of the any two equations. In result, we have six equations and six unknowns and the relevant coupling constants and can be determined by solving this system of equations.

The results for scenario (b) can be obtained from the results for scenario (a) with the help of aforementioned replacements.

From (18), it follows that, to estimate strong coupling constants and responsible for the decay of and , we need the residues of and baryons. For calculation of these residues for , we consider the following two point correlation functions:

The interpolating currents and couple not only to ground states, but also to negative (positive) parity excited states; therefore their contributions should be taken into account. In result, for physical parts of the correlation functions we getwhere the dots denote contributions of higher states and continuum. The matrix elements in these expressions are defined as

As we already noted, only the structure describes the contribution coming from baryons. Therefore we retain only this structure.

For the physical parts of the correlation function, we get

Here in the last term, upper (lower) sign corresponds to case (a) (case (b)).

Denoting the coefficients of the Lorentz structures and operators , and , as , , respectively, and performing Borel transformations with respect to , for spin case, we find

The expressions for spin case formally can be obtained from these expressions by replacing , , and . The invariant functions , from QCD side can be calculated straightforwardly by using the operator product expansion. Their expressions are presented in [15] (see also [5]).

Similar to the determination of the strong coupling constant, for obtaining the sum rule for residues we need the continuum subtraction. It can be performed in following way. In terms of the spectral density the Borel transformed can be written asThe continuum subtraction can be done by using the quark-hadron duality and for this aim it is enough to replace

It follows from the sum rules that we have only two equations, but six (three masses and three residues) are unknowns. In order to simplify the calculations, we take the masses of as input parameters. Hence, in this situation, we need only one extra equation, which can be obtained by performing derivatives over () on both sides of the equation. Note that the residues of baryons are calculated in a similar way.

#### 3. Numerical Analysis

In this section we present our numerical results of the sum rules for the strong coupling constants responsible for and decay derived in previous section. The Kaon distribution amplitudes are the key nonperturbative inputs of sum rules whose expressions are presented in [12]. The values of other input parameters are

The sum rules for and contain the continuum threshold , Borel variable , and parameter in interpolating current for spin particles. In order to extract reliable values of these constants from QCD sum rules, we must find the working regions of , , and in such a way that the result is insensitive to the variation of these parameters. The working region of is determined from conditions that the operator product expansion (OPE) series is convergent and higher states and continuum contributions should be suppressed. More accurately, the lower bound of is obtained by demanding the convergence of OPE and dominance of the perturbative contributions over the nonperturbative one. The upper bound of is determined from the condition that the pole contribution should be larger than the continuum and higher states contributions. We obtained that both conditions are satisfied when lies in the rangeThe continuum threshold is not arbitrary and related to the energy of the first excited state; that is, . Analysis of various sum rules shows that varies between and GeV, and in this analysis is chosen. As an example, in Figures 1 and 2 we present the dependence of the residues of and on for the scenario (a) at and several fixed values of , respectively. From these figures, we obtain that when lies between and the residues exhibit good stability with respect to the variation of and the results are practically insensitive to the variation of . And we deduce the following results for the residues: