Abstract

Considering the formalism of symplectic quantum mechanics, we investigate a two-dimensional nonrelativistic strong interacting system, describing a bound heavy quark-antiquark state. The potential has a linear component that is analyzed in the context of generalized fractional derivatives. For this purpose, the Schrödinger equation in phase space is solved with the linear potential. The ground state solution is obtained and analyzed through the Wigner function for the meson . One basic and fundamental result is that the fractional quantum phase-space analysis gives rise to the confinement of quarks in the meson, consistent with experimental results.

1. Introduction

Over the last decades, strong interaction has been analyzed by different approaches, including quantum chromodynamics (QCD) sum rules and lattice QCD, providing quantitative and qualitative characteristics of the hadronic matter [1]. In particular, systems such as quark-antiquark lead to interesting descriptions and a quantitative test for QCD and for both the particle-physics standard model [215]. In the case of a quarkonium, a popular approach considers the interaction between a quark-antiquark in a meson through a spatial (Euclidean) potential, such as , where is the distance between the quarks. The quantum nature of the state and the mass spectrum are studied by considering as a model for the state time-evolution of the Schrödinger equation. This corresponds to a specific sector of the strong interaction, which is called nonrelativistic QCD.

From the spectral analysis, it appears that the interaction in such systems as a heavy quark and an antiquark (charmonium ) is successfully modeled by the Cornell potential, which is defined by where is the string tension and is the strong coupling constant. The linear term is associated with the confinement, while the Coulomb-like term is a consequence of the asymptotic freedom. This potential has been used to investigate, as an instance, the confinement/deconfinement phase transition in hadronic matter [1620]. In addition, the Schrödinger equation with the Cornell potential as a model has been extensively used to explore quarkonium systems in the configuration space. This is the case of investigations of the heavy quarkonia mass spectroscopy and the bound state properties of and mesons [17, 2128].

Considering theoretical aspects, the heavy quarkonium characteristics have been analyzed by the Schrödinger equation with the Cornell potential through variational method in the framework of supersymmetric quantum mechanics [18, 2931]. The mass spectra of heavy quarks , , and within the framework of the Schrödinger equation with a general polynomial potential were also addressed [1]. In this formulation, the Nikiforov-Uvarov (NU) method [18] was used to calculate energy eigenvalues. The radial Schrödinger equation is extended to higher dimensions, and the NU method is applied to a Cornell-type potential. As a consequence, in order to obtain the heavy quarkonium masses, the energy eigenvalues and the associated wave functions are determined [2]. The eigen-solutions and an inverted polynomial potential were obtained by using the NU procedure [3033].

It is important to notice that although the Cornell potential presents theoretical and experimental consistency with the standard model, aspects of the hadronic matter as confinement are not obviously derived from the model-dependent Schrödinger equation. Indeed, this is the case of solutions of the Schrödinger equation in the Euclidean space representation. However, the recent analysis of the Wigner function of such a system as mesons provides an interesting description of the confinement, by using the characteristics of the phase-space quantum mechanics. These achievements were carried out by considering the symplectic quantum mechanics, in which the Schrödinger equation is written in a phase-space representation [33]. Nevertheless, accounting for the physical richness of phase space, many aspects remain to be explored, such as the fractional structure of the symplectic Schrödinger describing a quark-antiquark system.

The use of fractional calculus has attracted attention in a variety of areas in physics [3439]. For heavy quarkonium systems, methods as the NU formulation and analytical iteration have been explored to provide analytical solutions of the -dimensional radial Schrödinger equation in the framework of fractional space [35, 40]. A category of potentials including the oscillator potential, the Woods-Saxon potential, and the Hulthen potential have also been studied analytically with fractional radial Schrödinger equation by NU method [40, 41]. In order to investigate the binding energy and temperature dissociation, the conformable fractional formulation was extended to a finite temperature context [40]. A fundamental goal of the present work is to survey the applicability of the fractional approach to the study of quark dynamics in phase space.

Then, the behavior of the Wigner function for the ground state of meson is analyzed from the perspective of fractional calculus. For this purpose, the symplectic Schrödinger equation is rewritten in the fractional form with the linear term of the Cornell potential for the heavy meson. Beyond physical aspects, the analysis provides a simpler procedure to study this type of systems.

The work is organized as follows. In Section 2, some aspects of the Schrödinger equation represented in phase space are reviewed in particular to fix the notation. In Section 3, the concept of fractional derivative is implemented in the symplectic Schrödinger equation for the linear part of the Cornell potential. Section 4 is devoted the discussion of outcomes. In Section 5, summary and final concluding remarks are presented.

2. Symplectic Quantum Mechanics: Notation and Wigner Function

Considering a phase-space manifold , where a point is specified by a set of real coordinates , the complex valued square-integrable functions, , such that, , is equipped with a Hilbert space structure, . Here, stands by a vector in the Euclidean manifold, and stands for points in the dual . The point is a vector in the cotangent bundle of , equipped with a symplectic two-form [22]. In this way, can be used to introduce a basis in , denoted by with completeness relation given by . It follows that , where is the dual vector of . The symplectic Hilbert space can be used as the representation space of symmetries. For the nonrelativistic Galilei group, position and momentum operators are written as

A symplectic structure of quantum mechanics is constructed in the following way. The Heisenberg commutation relation is fulfilled. And then, using the following operators: the set of commutation rules are obtained being zero for all the other commutations. It is known as Galilei-Lie algebra and is a central extension. The Galilei symmetries are defined by the operators , and , which stands, respectively, by the generators of spatial translations, Galilean boosts, rotations, and time translations.

The time-translation generator, , leads to the time evolution of a symplectic wave function, i.e.,

The infinitesimal version of this equation reads as the Schrodinger-type equation in [42].

The physical interpretation of this formalism is obtained by the association of with a function , i.e., [4345].

In the next section, the representation symplectic Schrödinger equation in fractional context for the heavy quark system is obtained.

3. Fractional Symplectic Schrödinger Equation for the Confinement Potential

In this section, the symplectic Schrödinger equation is generalized to a fractional-space Schrödinger equation describing two particles interacting to each other by the linear part of the Cornell potential. Using the results of the previous section, the symplectic Schrödinger equation takes the form [33]

Using Equations (2) and (3) in (9), it leads to where natural units are used, such that . By using the transformation , this equations reads where . Writing Equation (11) in fractional from [34], it follows that where and is a scalar factor, and . Thus,

Therefore, Equation (11) in the fractional form is written as where with being a scale factor. It worth noting that if , one obtains the original Equation (11).

4. Discussion of Results

In this section, in order to obtain an analytical function of Equation (15), the fractional parameter is taken as . (For other values, perturbative methods can be used. This will be not addressed in the present paper.) This leads to the following form:

The solution of this equation is given by where and are constants and and are the Kummer functions. One can regard as a physical solution since it is the only one that is square integrable. As a result, one can impose that . Additionally, if , the series becomes a polynomial in of degree not exceeding , where This circumstance allows us to write where , and

Notice that the energy does not depend explicitly on the kinetic energy; thus, the initial condition should be .

For the ground state, making the substitution into and again, one have

Using the fact that is real, the normalized Wigner function of the ground state is given by where the constant depends on the value of .

In Figure 1, the behavior of the and is described. In Figure 2, the difference is plotted as a function of the parameter , considering . This difference has for reached the order of value of experimental measurements [46]. It is worth emphasizing here that the linear part of the Cornell potential only does not provide a spectrum in agreement with experimental measurement [46, 47]. Here, since we have the parameters of the fractional derivatives, those results can be improved for values of and . The next point is to explore the behavior of the Wigner function in order to detail the behavior of the confinement of quark-antiquark.

The Wigner function for the fractional parameter for and different values of are presented in Figure 3. The figure compares fractional Wigner functions to the original one, which is . We observe that the peaks diminish by lowering .

The curves (a) to (c) of Figure 3 show that with the increase of, the peaks of the Wigner functions increase. Additionally, we see that the peaks decrease to zero as goes to zero. Does functions as a fitting parameter for the fractional Wigner function of the meson? The curve (d) is the original Wigner function without the fractional parameters for meson [34]. Figure 4 shows that with the increase of, the distancedecreases. The maximum value of is the best fit for the case of . When compared to the experimental evidence, for comparison, the experimental value for the maximum distance is  MeV-1 [1].

For the case of general Cornell potential (Equation (1)), one can linearize to get an approximation form. In the first approximation, where , , and are constants. Then, the Hamiltonian is

This equation leads to

It is worth noting that Equation (25) is the same as Equation (9) with and . Therefore, the same analysis applies here. The energy is given by where and . It is noteworthy that, when , in the general Cornell potential has the that is responsible by interaction at short distances and corresponding to one gluon exchange. In addition, Table 1 presents the theoretical results from the fractional model for and , calculated from Equation (20), and the respective experimental values. Comparisons were established only for 1S states, as our theoretical model is applicable only to such states. We did not include spin in our theoretical model. We noticed that there is good accuracy between the theoretical and experimental results [48], better than those obtained by other theoretical models [46].

5. Remarks

We have studied the symplectic Schrödinger-like equation in the presence of a linear potential using the formalism of generalized fractional derivatives for nonrelativistic heavy quarkonium bound state. For this purpose, we have investigated the behavior of the Wigner function for the ground-state meson considering the symplectic quantum mechanics and the generalized fractional derivative constructed in [33, 34]. The Wigner function has been obtained for the charmonium state in the fractional form using the generalized fractional derivative as in Ref. [36], where we obtained the classical case at . To obtain an analytical solution, we analyse the case . For this value of , it was observed that the peaks of the Wigner function are lowered by decreasing fractional parameter ; therefore, this parameter can be used as a fitting parameter. For the case of , the value is the best fit considering the experimental evidence. Therefore, the present analysis seems to indicate the relevance of such a generalized fractional model based on the symplectic Schrödinger equation with linear term (Cornell potential) as far as quarkonium dynamics in phase space is concerned. To further the study of a quarkonium system within the fractional and phase-space approaches, we will include a quadratic term (or correction term) at the Cornell potential and other values for the fractional parameter . We also intend to study spinorial systems.

Data Availability

This is not applicable to our paper.

Disclosure

The present research can be found in the arXiv repository with the following identification number arXiv: 2209.12083 [hep-th].

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the CNPq and CAPES (Brazilian Agencies).