Abstract

In this paper, we study within the structure of Symplectic Quantum Mechanics a bidimensional nonrelativistic strong interaction system which represent the bound state of heavy quark-antiquark, where we consider a Cornell potential which consists of Coulomb-type plus linear potentials. First, we solve the Schrödinger equation in the phase space with the linear potential. The solution (ground state) is obtained and analyzed by means of the Wigner function related to Airy function for the meson. In the second case, to treat the Schrödinger-like equation in the phase space, a procedure based on the Bohlin transformation is presented and applied to the Cornell potential. In this case, the system is separated into two parts, one analogous to the oscillator and the other we treat using perturbation method. Then, we quantized the Hamiltonian with the aid of stars operators in the phase space representation so that we can determine through the algebraic method the eigenfunctions of the undisturbed Hamiltonian (oscillator solution), and the other part of the Hamiltonian was the perturbation method. The eigenfunctions found (undisturbed plus disturbed) are associated with the Wigner function via Weyl product using the representation theory of Galilei group in the phase space. The Wigner function is analyzed, and the nonclassicality of ground state and first excited state is studied by the nonclassicality indicator or negativity parameter of the Wigner function for this system. In some aspects, we observe that the Wigner function offers an easier way to visualize the nonclassic nature of meson system than the wavefunction does phase space.

1. Introduction

With the discovery of the meson in 1974 [1], the search to describe such heavy quark systems was based on the approach of potential models. The , which has a mass around 3.5 times that of the proton, is the lowest bound state of a charm and anticharm quark [2]. The meson appears to be an excellent device for testing QCD predictions. In quantum chromodynamics (QCD), the computation of meson properties is completely nonperturbative. As a result, direct calculation utilizing lattice QCD (LQCD) techniques is the only general method accessible. However, other strategies are effective for heavy quarkonia. The light quarks in a meson, the mass of the bound state is much larger than the mass of the quark, move at relativistic speeds. On the other hand, the speed of the charm and bottom quarks in their respective quarkonia is tiny enough that relativistic effects in these states are greatly diminished. The velocity, , is believed to be around for charmonia. An extension in powers of can then be used to approximate the computation. This method is known as nonrelativistic QCD (NRQCD).

The masses of heavy quarkonium states can be calculated using effective potential models. The nonrelativistic motion of the quarks that make up the quarkonium state is exploited in this technique to imply that they move in a static potential, comparable to nonrelativistic hydrogen atom models. One of the most prominent potential models is the Cornell (or funnel) potential [3]. Specifically, the system is characterized by a linear combination of the Coulomb and linear potentials. The study of these potentials should take into account the two important features of the Quantum Chromodynamics (QCD), namely, asymptotic freedom and quark confinement [47]. Another interest is that the Cornell potential can be used to analyze the transition between the confined and unconfined phases of matter [810]. This potential is of considerable importance in various branches of physics such as propagation of gravitational waves, particle and nuclear physics [1115], mathematical modeling of the parton vibrations inside hadronic system [1620], quantum chromodynamics, and atomic physics [5, 9, 11, 15]. It takes the form: where the first term is responsible for the interaction between quark and its antiquark by gluon exchange at short distances. The second term is responsible for quark confinement at large distances. It is known that the potential given in Eq. (1) reproduces two specific features of the strong interaction. The Coulomb-like term displays a property known as asymptotic freedom. This phenomenon tells us that strong interaction coupling constant is a function of momentum transfer. At short distances, the momentum transfer increases in a quark-antiquark collision. In this regard, the coupling constant becomes so small that quarks, and gluons can be considered as approximately free, and their interactions can be treated by a perturbation theory [2126]. Over large distances, momentum transfer decreases as the strong interaction coupling constant becomes larger. On the other extreme, it features a mechanism called confinement that keeps quarks and antiquarks permanently within hadrons at 1 Fermi separations [22, 25, 27, 28]. Confinement is explained consistently in the flux tube model. The gluon field between a pair of color charges develops a string-like configuration as the distance between quarks and antiquarks is extended (flux tube). As a result, across extended distances, the energy density in the tube between the quarks maintaining the gluon field remains constant [22, 25, 29]. The energy stored in the field is then proportional to the quark separation, giving us the linear term in the potential of Eq. (1). In that way, at relatively larger separation, the potential energy can create new quarks pairs in colorless forms instead of a free quark. The confinement has been a mystery until now since no one has been able to confirm it; although, there has been recent progress using the techniques of lattice QCD [22, 25, 28].

There are several studies with the Cornell potential in literature. For instance, Flores [9, 30] studied the Schrödinger equation with the Cornell potential using variational method and supersymmetric quantum mechanics [9, 10, 3032]. Bruni et al. [8] calculated the energy configuration for a quark-antiquark pair from the Nambu-Goto action. It was shown that this configuration energy has the shape of a Cornell potential [8]. Khoka et al. [31, 33] utilized the analytical exact iterative method to solve -dimensional Schrödinger equation with an extended Cornell potential [31]. They found the energy and mass spectrum of heavy quarkonia [31]. Omugbe [31] applied the Nikiforov-Uvarov method to obtain the eigensolutions of the radial Schrödinger equation with Cornell potential plus an inversely quadratic potential [11, 31, 34]. Until now, an analysis of such states in quantum phase space via the Wigner function is still lacking. In this work, we investigate confinement term of strong interaction system described by heavy quarks-antiquarks in the phase space.

The first successful formalism of quantum mechanics in phase space was introduced by Wigner in 1932 [3557]. He was motivated by the problem of finding a way to improve the quantum statistical mechanics [4247, 5054, 5870]. In the Wigner formalism, an operator defined in Hilbert space , is associated to a function in a phase space, says . The product of two operators, , both defined in , is associated to a Weyl product or star product of their phase space correspondents, i.e. [6169], where . Such result can be seen as the action of a star operator, , in the function . This fact is important for define the star operators , which are used to build up representations of symmetry groups in a symplectic manifold [5052, 64]. This gives rise, for example, to the Klein-Gordon and Dirac equations written in phase space [44, 61, 62, 64]. The connection with Wigner function is derived, providing a physical interpretation for the representation. These symplectic representations provide a way to consider a perturbative approach for Wigner function on the bases of symmetry groups.

The aim of this work is to analyze the behavior of the Wigner function given in terms of Airy function for the ground state of meson in the formalism of the symplectic quantum mechanics. For this, we consider only the linear potential that is a model more simple of the Cornell potential. In sequence, we studied the Schrodinger equation in the phase space to obtain the energy eigenvalues and corresponding wave functions by using the Bohlin transformation and perturbation theory for the Cornell potential. Then, we calculate the Wigner function for the heavy meson (). Consequently, we calculate the negativity parameter of the Wigner function associated to the ground state and first excited state of meson that present nonclassical behavior. The study of the Wigner function for such a system is crucial in order to learn more about the system’s chaotic nature.

The paper is organized as follows. In Section 2, we write the Schrödinger equation in phase space, and we perform the relation between phase space amplitude and Wigner function. In Section 3, we solve the Schrödinger equation in phase space for the quark-antiquark potential and calculate the Wigner function. In Section 4, the perturbative method is presented in order to solve the Schrödinger equation in phase space with the Cornell potential, and the Wigner function is calculated. The conclusion is presented in Section 5.

2. Symplectic Quantum Mechanics: Outline and Notation

In this section, we construct a formalism to quantum mechanics in phase space. For this purpose, we introduce an Hilbert space associated to phase space, denoted by . In this sense, the association of the Hilbert space with the phase space is given by where . Then, we write , with the aid such that is the dual vector of . This symplectic Hilbert space is denoted by . To construct a symplectic representation of quantum mechanics, we define the momentum and position operator

which satisfy the Heisenberg commutation relation .

Introduced the following operators, We obtain the following set of commutation rules: with all the other commutations being zero. This is the Galilei-Lie algebra with a central extension characterized by , where , and are the generators of translation, Galilean boosts, rotation, and time translation, respectively.

The evolution of the wave function in phase space is derivating with respect to time, and we obtain This is the Schrodinger-like equation in phase space [54].

The association of with a function is given by [54]. which has all the properties of the Wigner function [60]. In the next section, we solve the representation Schrödinger equation in phase space and establish the relation of amplitudes in phase space and the Wigner function for the meson system.

3. Quark Confinement and Schrödinger Equation in Phase Space

In this section, we analyze the quark confinement in the context of phase space. For this purpose, we consider the second term of Cornell potential given in Eq. (1). Therefore, we find a solution for the Scrödinger equation in phase space with linear potential and construct the associated Wigner function. We write the Schrödinger equation in phase space with the linear potential as follows:

And using equations (3) and (4)in (10), we get the Schrodinger equation in the form where we used . By using the transformation , such that With this, we obtain where . Then, the Schrodinger equation is reduced to the Airy equation (49)

where . The solution of this equation is where and are homogeneous Airy functions. However, goes to infinity for , and this solution is not relevant; so, [71]. In this case, the solution of the Eq. (14) is reduced to [71] where is a real function. As consequence, by virtue of associativity, , we write the solution of Wigner function in terms of the Airy function

where is a normalization constant. By condition , we can determine the energy levels

So, the energy of system are determined by where is the zero of the Airy function. In our case, we consider only the energy for ground state

In consequence, the normalization factor, see Ref. [71], is determined by Then, we get

where and are the reduced mass of the constituent quark and antiquark. In the sequence, we analyze the solution above of the Wigner function for meson. The behavior of the Wigner function for the strong interaction of the heavy quarks is shown in (Figures 1 and 2).

For the meson, we are taking into account Eq. (22). The meson is shown in Figure 1 represented by Wigner function for ground state. The experimental parameter, linear confinement parameter, reduced mass, and maximum relative distance between quark and antiquark are all in units of from meson which is shown in Table 1. Figure 1 is observed for ground state of charmonium meson which the Wigner function provides values negative in the phase space. This behavior is related to the character quantum. Then, it was make cut in graph as illustrated in Figure 2. Figure 2 shows that when , the Wigner function presents negative values, curve (a). This is because of quantum interference, also visualized in the curves (b) and (c). In particular, corresponds a limit. It is worth noting that when , the energy is . When varying the kinetic energy to , the Wigner function is displaced to the left, curve (b). We can keep increasing the kinetic energy for , which further evidences the distance for the left of the graph of the Wigner function, curve (c). By varying the kinetic energy for , it is noticed that the graph of the Wigner function is moved away visibly to the left side, curve (d). Thus, it is observed that there is a limit to the right for the existence of the Wigner function not to be zero. This limit approaches the experimental value that refers maximum relative distance (Table 1). Then, a boundary condition is placed at zero. Consequently, the meson only exists between zero and where Wigner function decay. Note that the analysis in phase space is revealing that for the variation of the kinetic energy there is an upper limit of existence of the charm-anticharm meson (Figure 2). Although it is simple to describe the quark case through a linear potential, the confinement is observed, which is not seen with the configuration space solution.

4. Cornell Potential and Schrödinger Equation in Phase Space

In this section, we analyze the quark-antiquark interaction in phase space context. For this purpose, we consider the Cornell potential given in Eq. (1). In this way, we present a solution for the Schrödinger equation with Cornell potential in phase space and the associated Wigner function. The two-dimensional Hamiltonian for the meson is given as where is reduced mass in units of the . In order to solve the Schrödinger equation for this Hamiltonian, the Bohlin transformation is utilized. Then, the Bohlin mapping is defined by [52, 58, 59] by defining which leads to

Substituting Eqs. (25), (26), (28), and (29) in Eq. (23) leads to the Hamiltonian

Consequently, the hypersurface in phase space defined by leads to where and are the canonical momenta conjugate. Eq. (31) can be written as

Observe that Eq. (32) is obtained from the classical Hamiltonian through of the star product. Therefore, the Bohlin mapping that leads to Eq. (31) is a classical transformation. For this reason, the perturbation theory is used to analyze the Eq. (32). So, the equation in phase space is given by where and .

The equation for has the form where corresponds the eigenfunction the unperturbed Hamiltonian.

Defining the operators

where and the star operators , are given by So, the Hamiltonian is

The unperturbed Hamiltonian is defined as and the perturbed part is

The equation that is to be analyzed is given by

The unperturbed equation is

The unperturbed part has solutions given by where and are solutions. The acting operators in the states vectors is

Utilizing the relations , , the ground state solution is where and are Laguerre polynomials, and is a normalisation constant. The eigenvalue solutions obtained of Eq. (41) are

The states excited are obtained from Eq. (47) using operators given in Eqs. (43)–(46). So, the solution for the first-order perturbed is given by

Note that it is required to solve the following expression before a solution for Eq. (49).

Using the orthogonality relations, The ground state of meson is

And for excited states of the meson, the wave functions are

The Wigner function for the meson is given by

To obtain the corrections of the first order of energy that is requiring to solve the matrix, where . Hence, we get

The eigenvalue obtained is and . Then, the corrections of the first order of eigenvalue of Eq. (40) are

Figures 3(a) and 3(b) show plots of the Wigner function associated to the ground state and first excited state for meson. It should be noted that all plots consider the coordinates and constant in order to show a three-dimensional figure, while , . In Figure 3(b), one clearly sees the presence of negative values for the Wigner function at the level excited of the meson. In Figure 3(a), the Wigner function looks strictly positive to the ground state.

With the Wigner function determined for the meson, we can calculate the negativity parameter for this system, and this parameter is correlated to the nonclassicality of the system [73]. The results of this calculation are shown in (Table 2). As can be seen in (Table 2), the negativity parameter increases when and grow. The negativity parameter is a relevant tool in the context of Wigner formalism because one is related to nonclassicality of physical systems as soon as is also related to chaotic behavior. In this way, by examining the negativity parameter, it is possible to investigate elements that the traditional wavefunction formalism does not permit. Then, the study in phase space is crucial in order to understand more about the chaotic nature.

5. Concluding Remarks

In this paper, we have studied heavy quarkonium bound states using the Schrödinger-like equation with an approach the Cornell potential in the Symplectic Quantum Mechanics framework. In particular, this sector is characterized by a linear potential which leads to the quark confinement plus one Coulomb-type potential. This is called of Cornell potential that feature experimentally this sector of quantum chromodynamics QCD. Initially, the Wigner function was obtained for the meson system. It was observed that analyzing in the phase space allows for simple model to observe the mesons confinement and why it is not seen into configuration space solution. We check out through the analysis in the phase space that due to the variation of kinetic energy, it has a maximum limit of the existence to the meson systems given by graphics obtained. After that, we have obtained the energy eigenvalues and the respective wave functions using perturbation theory. Such that the states, denoted a quasiamplitude of probability, are associated with Wigner function by the Weyl product. Using the Wigner functions, the negativity parameter for the ground state and first excited state of the system was derived. Finally, we observe through the analysis in the phase space that due to the change of energy levels, the negativity parameter increases Table 2. Which evidenced an increase in the quantum character of the system. It was not seen in the literature recent works.

Data Availability

Data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors would like to thank CNPq and CAPES of Brazil for their support in this study.