Abstract

The chromopolarizability of a quarkonium describes the quarkonium’s interaction with soft gluonic fields and can be measured in the heavy quarkonium decay. Within the framework of dispersion theory which considers the final state interaction (FSI) model-independently, we analyze the transition and obtain the chromopolarizability and the parameter . It is found that the FSI plays an important role in extracting the chromopolarizability from the experimental data. The obtained chromopolarizability with the FSI is reduced to about 1/2 of that without the FSI. With the FSI, we determine the chromopolarizability and the parameter Our results could be useful in studying the interactions of charmonium with light hadrons.

1. Introduction

The chromopolarizability of a quarkonium parametrizes the quarkonium’s effective interaction with soft gluons, and it is an important quantity in the heavy quark effective theory. Within the multipole expansion in QCD in terms of the chromopolarizability, many processes can be described, including the hadronic transitions between quarkonium resonances [1, 2] and the interaction of slow quarkonium with a nuclear medium [3]. A recent interest of the chromopolarizabilities of and comes from the hadrocharmonium [38] interpretation of the and observed by the LHCb Collaboration, and it is found that the can be interpreted as a -nucleon bound state if [9].

There are a few studies of the chromopolarizabilities of and , some of which are not in line with each other. Calculated in the large- limit in the heavy quark approximation, the values of the chromopolarizabilities of the and are obtained: and [6, 7, 10, 11]. Within a quarkonium-nucleon effective field theory, the chromopolarizability of the is determined through fitting the lattice QCD data [12] of the -nucleon potential, and the result is [13, 14]. Based on an effective potential formalism given in [15] and a recent lattice QCD calculation [16], the chromopolarizabilities of are extracted to be [9]. On the other hand, the determination of the transitional chromopolarizability is of importance since it acts a reference benchmark for either of the diagonal terms due to the Schwartz inequality: [4]. The perturbative prediction in the large limit is [6, 7, 10, 11]. While being extracted from the process of , the result is [4, 17]. Taking account of the FSI in a chiral unitary approach, it is found that the value of may be reduced to about of that without the FSI [18].

Since the FSI plays an important role in the heavy quarkonium transitions and modifies the value of significantly, it is thus necessary to account for the FSI properly. In this work we will use the dispersion theory to take into account of the FSI and extract the value of . Instead of the chiral unitary approach [18, 19], in which the scalar mesons (, , and ) are dynamically generated, in the dispersion theory the FSI is treated in a model-independent way consistent with scattering data. Another update of our calculation is that we consider the FSIs of separate partial waves, namely, the - and -waves, instead of only accounting for the -wave as in the parametrization in [17, 18].

The theoretical framework is described in detail in Section 2. In Section 3, we fit the decay amplitudes to the data for the transition and determine the chromopolarizability and the parameter . A brief summary will be presented in Section 4.

2. Theoretical Framework

First we define the Mandelstam variables of the decay process

The amplitude for the transition between -wave states and of heavy quarkonium can be written as [4, 13]where the factor appears due to the relativistic normalization of the decay amplitude, is the chromopolarizability, and denotes the chromoelectric field. is the first coefficient of the QCD beta function, , where and are the number of colors and of light flavors, respectively. , and , where is a parameter that can be determined from the data.

The above result of the QCD multipole expansion together with the soft-pion theorem can be reproduced by constructing a chiral effective Lagrangian for the transition. Since the spin-dependent interactions are suppressed by the charm mass, the heavy quarkonia can be expressed in terms of spin multiplets, and one has , where contains the Pauli matrices and and annihilates the and states, respectively [20]. The effective Lagrangian, at the leading order in the chiral as well as the heavy quark nonrelativistic expansion, reads [2123]where is the velocity of the heavy quark. The Goldstone bosons of the spontaneous breaking of chiral symmetry can be parametrized according to where denotes the pion decay constant.

The amplitude obtained by using the effective Lagrangians in (3) isMatching the amplitude in (2) to that in (5), we can express the low-energy couplings in the chiral effective Lagrangian in terms of the chromopolarizability and the parameter

The partial-wave decomposition of can be easily performed by using the relationwhere is the 3-momentum of the final vector meson in the rest frame of the initial state with , , and is the angle between the 3-momentum of the in the rest frame of the system and that of the system in the rest frame of the initial .

Parity and -parity conservations require the pion pair to have even relative angular momentum . We only consider the - and -wave components in this study, neglecting the effects of higher partial waves. Explicitly, the - and -wave components of the amplitude read

There are strong FSI in the system especially in the isospin-0 S-wave, which can be taken into account model-independently using dispersion theory [2231]. We will use the Omnès solution to obtain the amplitude including FSI. In the region of elastic rescattering, the partial-wave unitarity conditions readBelow the inelastic threshold, the phases of the partial-wave amplitudes of isospin and angular momentum coincide with the elastic phase shifts modulo , as required by Watson’s theorem [32, 33]. It is known that the standard Omnès solution of (9) is as follows: where the is a polynomial, and the Omnès function is defined as [34]

At low energies, and can be matched to the chiral representation. Namely, in the limit of switching off the FSI, i.e., , the polynomials can be identified exactly with the expressions given in (8). Therefore, the amplitudes including the FSI take the form

Now we discuss the phase shifts used in the calculation of the Omnès functions. For the -wave, we use the phase of the nonstrange pion scalar form factor as determined in [35], which yields a good description below the onset of the threshold. For the -wave, we employ the parametrization for given by the Madrid–Kraków collaboration [36]. Both phases are guided smoothly to for .

It is then straightforward to calculate the invariant mass spectrum and helicity angular distribution for usingwhere the Legendre polynomial

3. Phenomenological Discussion

The unknown parameters are the low-energy constants and in the chiral Lagrangian (3), which can be expressed in terms of the chromopolarizability and the parameter as in (6). In order to determine and , we fit the theoretical results to the experimental invariant mass spectra and the helicity angular distribution from the BES decay data [37] and the corresponding decay width [38]. The fit results are plotted in Figure 1, where the red solid and blue dashed curves represent the results with or without the FSI, respectively. The fit parameters as well as the are shown in Table 1. One observes that the experimental data can be well described regardless of whether the FSI is included. This is due to the simple shapes of the invariant mass distribution and the helicity angular distribution in this process and does not mean the FSI is not important. Since the dipion mass invariant mass reaches about 600 MeV in such a decay, the FSI is known to be strong in this energy range and needs to be considered. On the other hand, one can readily see from (6) and (8) that while the chromopolarizability determines the overall decay rate, the parameter characterizes the -wave contribution, and we do not find significant correlation between and .

We observe that the FSI modifies the value of the chromopolarizability significantly, and resultant value with the FSI is almost 1/2 of that without the FSI. The obtained value with the FSI, , coincides with the suspicion [3] with the value from the calculation [9] based on the recent lattice QCD data of -nucleon potential [16]. It should be mentioned that the value of with the FSI obtained here is different from the one in [18], , and also our result without the FSI slightly differs from those in [17, 18]. The reasons are that the chiral unitary approach instead of dispersion theory is used to account for the FSI in [18], and we use the updated experimental data [37, 38] and a general theoretical amplitude rather than the one only containing the -wave as employed in [17, 18].

For the parameter , as shown in Table 1 its value is affected little by the FSI. One notes that a detailed study of the process using the Novikov-Shifman model [2] has been performed by BES [39], and based on the joint - distribution this parameter was determined as . We have tried fitting the same old BES data [39] and our changes slightly and is still much smaller than the BES one. In the Novikov-Shifman model, the amplitude reads [2]If we make the same approximation, namely, neglect the terms except the ones, as in [2] and set in the expression of 3-momentum , our amplitude without the FSI agrees with (15). While numerically we find that some neglected terms are at the same order as the term in (15), this may account for the difference of between ours and that in [39]. On the other hand, we have checked that the contribution of the -wave, which is characterized by the parameter , to the total rate is less than two percent, and the same observation has been made in [39].

4. Conclusions

We have used dispersion theory to study the FSI in the decay . Through fitting the data of the mass spectra and the angular distributions, the values of the chromopolarizability and the parameter are determined. It is found that the effect of the FSI is quite sizeable in the chromopolarizability , and the one with FSI is almost of that without FSI. The parameter , which accounts for the -wave contribution, is affected little by the FSI. The results obtained in this work would be valuable to understand the chromopolarizability of charmonia and will have applications for the studies of the nucleon-charmonia interaction.

Data Availability

All the data used in this work are from [37, 38].

Conflicts of Interest

The author declares that they have no conflicts of interest.

Acknowledgments

We acknowledge Feng-Kun Guo for the proposal for this work and for the useful comments on the manuscript. This research is supported in part by the Fundamental Research Funds for the Central Universities under Grant no. 06500077.