Abstract

The glueballs lead to gluon and QCD monopole condensations as by-products of color confinement. A color dielectric function coupled with a Abelian gauge field is properly defined to mediate the glueball interactions at confining regime after spontaneous symmetry breaking (SSB) of the gauge symmetry. The particles are expected to form through the quark-gluon plasma (QGP) hadronization phase where the free quarks and gluons start clamping together to form hadrons. The QCD-like vacuum , confining potential , string tension , penetration depth , superconducting and normal monopole densities (), and the effective masses ( and ) will be investigated at finite temperature . We also calculate the strong “running” coupling and subsequently the QCD -function. The dual superconducting nature of the QCD vacuum will be investigated based on monopole condensation.

1. Introduction

The Large Hadron Collider (LHC) is currently attracting the attention of particle physicists to investigate the dynamics of the standard model (SM) at the TeV scale and to probe for a possible new physics. One of the major expectations at the LHC is the observation of jet substructures. Generally, jets are collimated bundles of hadrons constituted by quarks and gluons at short distances or high energies [13]. Jets played a significant role in the discovery of gluons () [47] and top quark () [8, 9]. These observations are featured prominently in classifying Quantum Chromodynamics (QCD) as a theory for strong interaction within the standard model. The most abundant secondary particles produced in heavy ion collisions are pions which emerge from hot and dense hadron gas states. These pions can be studied from the projected quark-gluon plasma (QGP) and the deconfinement phase through emitted photons and dileptons produced at this phase. Since the photon and the dilepton do not interact with the hadron matter through strong force, they decay quickly [1013] making it easy to study the pions. Due to the large production of jets and pions at the LHC and Relativistic Heavy Ion Collider (RHIC), they are suitable candidates for the search for new particles and the testing of the SM properties. Quarks and gluons behave as quasi-free particles at higher energies (short distances) due to asymptotically free nature of the QCD theory. When the energy among these color particles is reduced to about 1 GeV and below or at a separation of order 1 fm or higher, color confinement sets in, such that the quarks and the gluons coexist as hadrons. At lower enough energies, hadronization is highly favoured leading to the formation of light pions. Hadron jet was first produced by annihilation of electron-positron to produce a two-jet event, . This two-jet structure [14] was first observed in 1975 at SPEAR (SLAC), and the spin nature of the quarks was also established [15]. The jets are produced when fly apart; more pairs are produced which recombine with the existing pairs to form mesons and baryons as a by-product [16]. The first three-jet structure from the electron-positron annihilation process was observed by analysing 40 hadronic events at the center-of-mass energy 27.4 GeV from SPEAR [16]. These jets were subsequently observed clearly and separately [17, 18]. The discovery of the gluon through the three-jet process played a major role in the discovery of the Higgs particle [1924] by ATLAS collaboration [25] and CMS collaboration [26] using the LHC at CERN. On the other hand, pion-pion annihilation serves as the principal source of dileptons ( and ) produced from the hadron matter [27, 28]. Thus, the proper study of the main source of dilepton spectra observed experimentally proposes significant observables which help to understand the pion dynamics in a dense nuclear matter that exist at the beginning of the collision [29].

We consider a complex scalar field model coupled to Abelian gauge fields in two different ways. At relatively low energies, the particles undergo transformation into glueballs. The emission of gluon which mediates the strong interactions is a non-Abelian feature observed in the nonperturbative regime of QCD theory. However, in the model framework, the gluon emission is a result of the modification of the Abelian gauge field by the color dielectric function. In this approach, we approximate the non-Abelian gauge field responsible for color confinement with an Abelian gauge field coupled with a dimensionless color dielectric function . The color dielectric function is responsible for regulating the long-distance dynamics of the photon propagator, so it does not decouple at higher wavelengths where glueballs are expected to be in a confined state. The color dielectric function is defined in terms of the glueball field after SSB, to give meaning to color confinement and its associated properties. While initiates the annihilation and production processes, is the glueball field which brings about confinement and gluon condensation. The Abelian approach to QCD theory was first proposed by ‘t Hooft to justify magnetic monopole condensation in QCD vacuum leading to color superconductivity [30, 31] that also suggested “infrared Abelian dominance.” Fast forward, it has also been shown that about 92% of the QCD string tension is Abelian [32], so an Abelian approximation will not be out of place. Further studies of Abelian dominance for color confinement in and lattice QCD can be found in Ref. [3337]. An interested reader can also refer to Ref. [38] for more recent results on Abelian dominance in QCD theory. This approach has also been adopted in other confinement models—see [3941] and references therein for more justifications.

The annihilation of through “scalar QED” where the particles are relatively free can be used to address hadronization due to the presence of high-density deconfined mater. The spontaneous symmetry breaking of the symmetry to an energy regime is suitable for describing color confinement with a glueball field . Because the running coupling, , of the strong interaction decreases with momentum increases (short distances), annihilation can be studied from perturbation theory () similar to annihilation and creation of pions () and Bhabha scattering () [42]. A low energy analyses of annihilation will be made to show color confinement, bound states of the gluons (glueball masses) [4345], and the QCD vacuum [46, 47] responsible for the gluon condensate. We will shed light on monopole condensation in the QCD vacuum and its role in screening the QCD vacuum from external electric and magnetic field penetrations similar to superconductors. There is a long-standing belief that quark masses are obtained through the spontaneous symmetry breaking mechanism facilitated by the nonzero expectation value of the Higgs field. That notwithstanding, quark masses continue to be an input parameter in the standard model [48]. As a free parameter in the SM, the magnitude of the mass is determined phenomenologically and the result compared with sum rules, lattice simulation results, other theoretical techniques, and experimental data. The Dyson-Schwinger integral equation for the quark mass function has two known solutions, i.e., trivial and nontrivial solutions. The latter is obtained under low energy and nonperturbative conditions while the former leads to unphysical, massless quarks. Thus, the nontrivial solution leads to the creation of quark masses resulting in “dynamical chiral symmetry breaking,” which is a consequence of confinement [4954]. Monopoles can be said to be a product of grand unification theory (GUT). The mixing of the strong and the electroweak interactions due to the higher gauge symmetries in GUT breaks spontaneously at higher energies or extremely short distances. So, physical features of monopoles such as size and mass are studied through the energy of spontaneous symmetry breaking [55]. Monopoles play important role in color confinement. Monopole condensation produces a dual superconductor which squeezes the uniform electric field at the confinement phase into a thin flux tube picture. The flux tube is formed between quark and an antiquark to keep them confined. This scenario has been extensively investigated in lattice gauge theory [5658] and Polyakov loop [59] to establish its confinement properties. The involvement of monopoles in chiral symmetry has recently been investigated as well [60, 61]. The impact of monopoles in quark-gluon plasma (QGP) has also been studied in [6264].

The paper is organised as follows: In Section 2, we introduce the Lagrangian density that will be the basics for this study. In Section 3, we present spontaneous symmetry breaking of the model presented in the previous sections; this section is divided into four subsections; in Section 3.1, we investigate confinement of glueballs; in Section 3.2, we study the effective masses; we dedicate Section 3.3 to gluon condensation, and in Section 4, we investigate the monopole condensation. We proceed to present the analysis in Section 5.1 and the final findings in Section 5.2. We have adopted the natural units , except otherwise stated.

2. The Model

We start with a Lagrangian density developed by exploring gauge-invariant properties (see [39]): where , , and are the covariant derivative, strength of the dual gauge field, and the gauge field, respectively. This model can be used to study the annihilation and creation of identical particles and its transformation into glueballs at low energies. The complex scalar fields and its conjugate describe particle and its antiparticle with the same mass but different charges, respectively. The electromagnetic interactions among the scalar fields with an exchange of single photon produced from the dual gauge field dynamics are the suitable scenario to address a high energy limit. However, we shall focus on the strong interactions, where the scalar fields modify to glueball fields at relatively low energies, which will lead to color confinement and its associated properties mediated by the modified gauge field dynamics . The equations of motion are

We will adopt a complex scalar field potential of the form where

We can isolate the electromagnetic interactions among the fields [65] for study using Equation (3). This enables us to visualize how the fields annihilate and create identical particles [66, 67] through single-photon exchange, , in high energy regime [68].

3. Confinement, Effective Masses, Gluon Condensation, and Dual Superconductivity

To discuss the color confinement and its consequences, we need to follow spontaneous symmetry breaking of the symmetry in Lagrangian Equation (1) to give mass to the resulting glueball fields and their associated products such as the monopole condensation and gluon condensation which form the basics for color confinement. We will proceed with the symmetry breaking from the transformations:

The vacuum expectation value of the potential Equation (5), breaks the symmetry spontaneously. Considering two real scalar fields and representing small fluctuations about the vacuum, a shift in the vacuum can be expressed as where . Hence, the potential takes the form

Consequently, we can suitably parametrize the scalar field such that,

The Lagrangian in Equation (1) becomes where is associated with the glueball fields and is also associated with the massless Goldstone bosons which will be swallowed by the massive gauge fields through the gauge transformation [39]:

Pions are strongly interacting elementary particles with an integer spin. Therefore, they are bosons which are not governed by Pauli’s exclusion principle and can exist as relativistic or nonrelativistic particles. As a result, they can be represented by spin-0 scalar fields. In the nonrelativistic regime, they exist as Goldstone bosons , signifying the breakdown of chiral symmetry [69]. Pions are generally produced through matter annihilations such as , , and , which undergo transition from baryon structures to mesons. This is an interesting phenomenon in low energy hadron physics [70, 71]. Pions have a flavour structure and other quantum numbers that permit us to classify them as bound states of quark and an antiquark. However, the valence quarks which characterize the flavour structure are surrounded by gluons and quark-antiquark pairs. Physically, glue-rich components mix with pions () forming an enriched spectrum of isospin-zero states and hybrid states [72]. The glueball spectrum and their corresponding quantum numbers are known in lattice gauge theory predictions [43, 73]. Additionally, the existence of glueballs has not been decisive because of the fear of possible mixing with quark degrees of freedom [74]. In the model framework, the glueballs coexist with the Goldstone bosons which are subsequently swallowed by the massive gauge fields leaving out the glueballs () for analysis. The glueball decay occurs when the separation distance between the valence gluons exceeds some thresholds (>1 fm) leading to hadronization.

Consequently, the Lagrangian can be simplified as where . The equations of motion for this Lagrangian are

The Feynman propagator for the glueball field in Equation (12) reads where .

3.1. Confinement

Here, we are interested in calculating static confining potential, such that only chromoelectric flux that results in confining electric field is present, while chromomagnetic flux influence is completely eliminated. Thus, in spherical coordinates, Equations (15) and (17) become, respectively, where is the integration constant. We are also now setting (similar process was adopted in Ref. [3941]), where is a dimensionful constant associated with the Regge slope that absorbs the dimensionality of , so remains dimensionless and Equation (19) reads

In the last step, we have assumed that the particles are far away from the charge source , and in such a limit, terms of the order can be disregarded. This equation has a solution

Substituting this result into the electric field Equation (20), and using the well-known expression for determining the electrodynamic potential, to determine the confining potential , we get so, setting leads to with string tension

In the last step of Equation (25), we substituted and for simplicity.

3.2. Effective Masses

Glueballs are simply bound states of gluons, and they are “white” or colorless in nature. They are flavour blind and are capable of decaying. Scalar glueball with quantum number is observed by lattice QCD simulations as the lightest with a mass of about 1.7 GeV [4345].

The dispersion relation for the glueball and the gluon excitations including thermal fluctuations can be expressed as

We can calculate the effective thermal fluctuating glueball mass by taking the second derivative of the Lagrangian density Equation (14) with respect to :

We redefine the glueball field to include the thermal fluctuations as , with restriction, . The angle brackets represent thermal average, and is the mean field of the glueballs. This equation can be solved by determining explicitly the nature of and using field quantum distributions. By the standard approach, we can express Here, represents the Bose-Einstein distribution function and , where is the temperature. We can analytically solve Equation (31) at a high energy regime, where , and is the gluon/screening mass. Therefore, where we have substituted and used the standard integral

In the last step, we have substituted being the critical temperature. Following the same analysis, where we have assumed , with being the glueball mass, and considered the standard integral,

Thus, Equation (30) becomes

When we take a thermal average of Equation (15), we get ; therefore, and are exact solutions of the glueball fields [75]. We further demonstrate in Figures 1 and 2 that has solution at and increases steady as increases to its maximum at , corresponding to no glueball fields and the melting of the glueballs.

3.3. Gluon Condensation

Classical gluodynamics is described by the Lagrangian density: which is invariant under scale and conformal symmetries, , and also produces vanishing gluon condensate . Meanwhile, the symmetries are broken when quantum correction is added to the Lagrangian, i.e., resulting in nonvanishing gluon condensate, . This correction is the consequence of the well-known scale anomaly observed in QCD energy-momentum () trace where is the “so-called” beta-function of the strong coupling , with a leading term

Accordingly, we have nonzero vacuum expectation, expressed as

The second term in Equation (14) acts similar to the quantum correction and explicitly breaks the scale and the conformal symmetries, so its energy-momentum tensor satisfies Equation (42) [40, 7678]. We will compute the trace of the energy-momentum tensor of Equation (14) using the relation

Substituting the equation of motion Equation (15) into Equation (43) yields where in the last step we have substituted with and representing the derivative of and with respect to , respectively. In order to make the expression easy to be compared with Equation (42), we rescale as , consequently,

We recover the classical expression for when we set . Following the potential defined below Equation (14), we can express thus, Equation (45) becomes

The QCD vacuum is seen as a very dense state of matter comprising gauge fields and quarks that interact in a haphazard manner. These characteristics are hard to be seen experimentally because quark and gluon fields cannot be observed directly; only the color neutral hadrons are observable. The appearance of the mass term in the condensate is also significant because it leads to chiral symmetry breakdown in the vacuum. Additionally, the gluon mass can be derived from the vacuum as a function of the glueball mass as

Generally, scalar glueball mass is related to the gluon mass as considering the leading order [79, 80] from Yang-Mills theory with an auxiliary field , here, the scalar glueball mass represents fluctuations around . The gluon mass is determined to be  MeV as obtained from lattice simulations [8183]. Heavier gluon masses have also been observed in the range of ~1 GeV by phenomenological analysis [84, 85], lattice simulation [86], and analytical investigations [87]. The glueball mass is responsible for all the confinement properties and the chiral symmetry breaking in the vacuum. Observing that the string tension is  GeV/fm, we can identify , corresponding to a gluon mass , from the model framework.

In terms of temperature, we use Equation (35), and the gluon condensate becomes

Corresponding to a temperature fluctuating gluon mass where is a dimensionless coupling constant. This expression looks like the Debye mass derived from the leading order of QCD coupling expansion [8890]. It carries nonperturbative property of the theory. Comparing Equations (40) and (44), we identify

Additionally, we can determine the strong running coupling from the renormalization group theory [91]

therefore,

Comparatively, the strong running coupling can be identified as , with , as the space-like momentum associated with the four-vector momentum as . Thus, the color dielectric function is associated with the QCD -function and the strong running coupling. We can write the Feynman propagator for the interaction by considering the left-hand side of Equation (45): where we have added the gauge fixing term and made the Fourier transform

In the last step, the term in the square brackets is precisely the photon propagator normally associated with Abelian gauge fields, but the multiplicative factor is defined such that the photon propagator does not decouple at higher wavelengths. This influences the photons to behave like gluons.

Also, one can retrieve all the major results obtained in [40], if we consider temperature fluctuations in the scalar field . In that case, we move away from the known form of field theory, i.e., perturbation about the vacuum, , to account for the density of particles and their interactions about the vacuum and the thermal distributions in the scalar field. So, we can substitute into the potential in Equation (5), where plays a similar role as introduced above. Thus, here, we have used the high-temperature approximation , and the potential becomes where ,

The glueball mass can be calculated from the fluctuations in Equation (58) around the vacuum, such that where is the vacuum of the potential. We observe that at , similar to Equation (37). The differences in the degree of the temperature in both equations arise because Equation (37) has temperature correction to the gauge field , while in Equation (37), we have temperature correction to the scalar field . Following the same procedure as the one used in deriving Equation (47), we can express

Here, we have replaced with from Equation (47). The negative sign and differences between this result and Ref. [40] might be due to the differences in the methods adopted for calculating the energy-momentum trace tensor . Also, if we consider the IR part of the potential in Equation (25) as studied in Ref. [40], the potential becomes with string tension which is similar to the result in Ref. [40]. Therefore, the approach adopted in this paper accounts for the temperature corrections to both the glueball field in Equation (35) and the gauge field , while in Ref. [40], only the thermal correction to the scalar field was considered.

4. Dual Superconductivity

In this section, we will discuss briefly color superconductivity and proceed to elaborate on dual superconductivity in detail. The fields annihilate and create identical fields through a decay process, . The dual gauge field which mediates the interaction of the scalar fields during the annihilation process is also responsible for the monopole condensation. In the model framework, the scalar field can be seen as point-like, asymptotically free and degenerate at the high energy regime [92]. In the context of this work, the high energy regime is in reference to the phase at which the scalar fields are annihilating while low energy regime is after the symmetry breaking or the glueball regime. In effect, multiparticle states are formed in the high energy regime. In such a high particle density regime, the charged scalar fields form Bose-Einstein condensation [93, 94] similar to induced isospin imbalance systems leading to color superconductivity [95, 96]. With this type of condensation, there is a higher occupation number at the ground state than the excited states. So, temperature increase goes into reducing the occupation numbers [97]. There are several works on nonzero isospin chemical potential () and baryon chemical potential () in pion condensation [98100], and they are related to the two flavour quarks that constitute the pions as and , respectively. It has also been established [101] analytically at the quark level that the critical isospin chemical potential [102, 103] for pion superfluidity is precisely the pion mass in vacuum. This behaviour and other related isospin chemical potentials in pion condensate have been investigated using the Nambu-Jona-Lasinio (NJL) model, ladder-QCD [98].

In the model framework, the color superconducting property can be studied from the propagator, . The high density region where the particles are asymptotically free forms deconfined matter leading to Bose-Einstein condensation resulting in the color superconducting phase [96, 104, 105]. On the other hand, at the low energy region, the fields modify into glueballs with real mass, , leading to chiral symmetry breaking and color confinement. Since quark-quark combinations do not form color singlet, Cooper pair condensate will break the local color symmetry giving rise to color superconductivity [106]. This phenomenon was first studied by Bardeen, Cooper, and Schrieffer (BCS) [107]. The BCS mechanism for pairing in general appears to be more robust in dense quark matter than superconducting metals due to the nature of interactions among quarks. The virtual fermions in the vacuum make it color diamagnetic preventing both the electric and magnetic fields from penetrating [2].

We turn attention to the dual superconductivity in color confinement whose property we can derive from Equation (16). This subject was first proposed by Nambu, ‘t Hooft, and Mandelstam [56]. While the color superconducting phase takes place at the deconfined matter regime, the dual superconductivity is observed at the confining phase [108]. In the model framework, dual superconductivity is observed at the high glueball condensation region, , i.e., a region of strong confinement. We will use the rest of this section to discuss dual superconductivity in detail within the model framework. In this phenomenon, the QCD vacuum is seen as a color magnetic monopole condensate resulting in one-dimensional squeezing of the uniform electric flux that form on the surface of the quark and antiquark pairs by dual Meissner effect. This results in the formation of a color flux tube between quark and antiquark used in describing string pictures of hadrons [38]. In this regard, magnetic monopole condensation is crucial to color confinement and dual superconductivity. Monopoles also play a role in chiral symmetry breaking [109] and strongly coupled QGP. They are involved in the decay of plasma formed immediately after relativistic heavy ion collision [6264]. Even though monopoles have not been seen or proven experimentally, there are many theoretical bases for which one can believe their existence. The theoretical evidence for its existence is as strong as any undiscovered theoretical particle. Polyakov [110] and ‘t Hooft [111] discovered that monopoles are the consequences of the general idea of unification of fundamental interactions at short distances (high energies). Dirac on the other hand demonstrated the presence of monopoles in QED. Monopoles in general play an insightful role in understanding color confinement. They also provide some useful explanations to the features of superstring theory and supersymmetric quantum theory where the use of duality is a common phenomenon. The discovery of these monopoles some day will be interesting and possibly revolutionary since almost all existing technologies are based on electricity and magnetism.

Monopole condensation can now be exploited from the condensation of particles in the region , where . Hence, from Equation (16), we can express

The equations of motion for the static fields are here, and are the dual versions of the static electric and magnetic fields, while and are the magnetic current and charge densities, respectively. The Lorentz force associated with these fields can be expressed as where is the speed of the particle within the fields and is the monopole charge. The homogeneous equation represents the uniform electric field present at the confining phase. Combining the dual version of Ampere’s law on the right side of Equation (65) with the dual London equation responsible for persistent current generated by the monopoles we obtain the expression

Noting that , we can identify as the London penetration depth [112, 113]. Developing the Laplacian in Equation (68) in spherical coordinates yields this equation has a solution given as where is a constant. Using dimensional analysis, we can fix the constant as ; therefore, where . This means that the electric field is exponentially screened from the interior of the vacuum with penetration depth , a phenomenon known as the color Meissner effect. Also, using London’s accelerating current relation, and writing the continuity equation in the form of London equations, together with the equation at the left side of Equation (65), we obtain

This equation has a solution similar to Equation (69), i.e.,

Thus, the magnetic field is equally screened exponentially from the interior of the vacuum by a penetration depth . In sum, both time-varying electric and magnetic fields are screened with the same depth .

Now, comparing with the results in Refs. [112, 114], where we rearrange the expression obtained for in the model framework, to make it suitable for comparison. Juxtaposing these two equations leads to , and , where are mass, charge, and electron number density of a superconducting material. Hence, the electron charge is equivalent to the monopole charge [115, 116], and the electron number density and mass are related to the monopole density and mass, respectively. Combining Equation (67) with the electric field expression at the right side of Equation (65), we obtain a fluxoid quantization relation: where and is the quantum of the electric flux [112]. The penetration depth could also be expressed as a function of temperature

We retrieve the initial penetration depth below Equation (68), if we set . Accordingly, the penetration depth increases with temperature until it becomes infinite at , where deconfinement and restoration of chiral symmetry coincide [117, 118]. At this point, the field lines are expected to penetrate the vacuum causing it to lose its superconducting properties. This expression is precisely the same as the one obtained for metallic superconductors [114]. At , the magnetic and electric field lines go through the vacuum. This phenomenon is referred to as the Meissner effect. The number density of the superconducting monopoles can also be expressed in terms of temperature as

The number density decreases sharply with temperature and vanishes at . On the other hand, at , similar to the two-fluid model where normal monopole fluid mixes with superconducting ones, the density of the monopoles can be calculated from the relation , where is the normal monopole density [119]. Additionally, at , , and in the limit , . Consequently,

Thus, the color superconducting monopoles at conduct with dissipationless flow; on the other hand, the normal monopoles conduct with finite resistance at the same temperature range. Consequently, decreases sharply with increasing while increases steadily with increasing .

5. Analysis and Conclusion

5.1. Analysis

At a high energy regime, there is high particle density which forms deconfined matter leading to hadronization [120, 121], [68, 122125]. In this regime, the process takes place at low temperatures but in high momentum and particle density. This leads to an asymptotic freedom behaviour among the particles because color confinement depends on how far or close the particles are together. Thus, higher density means the particles are closer together leading to deconfinement. Additionally, Equation (55) represents the gluon propagator which mediates the interaction of the glueballs in the IR regime of the model. Hadronization [126] comes into play in this regime when the momenta of the particles are reduced below a particular threshold set by the QCD scale,  MeV, or the separation distance between particle and antiparticle pairs is greater than . For the sake of clarity, there are two forms of hadronization mentioned here, one arising from high energy QGP formation and the other from QCD string decay into hadrons in the low energy regime. The QGP hadronization is believed to have been formed immediately after the Big Bang when the QGP starts cooling down to Hagedorn temperature, , or higher baryon densities above where the free quarks and gluons start clamping together into hadrons. This is also observed at the initial stages of heavy ion collisions [121, 127, 128]. On the other hand, the string decay hadronization occurs due to nonperturbative effects. Under this picture, new hadrons are formed out of quark-antiquark pairs or through a gluon cascade [126].

Theoretically, Light-Front (LF) is one of the suitable theories for describing relativistic interactions due to its natural association with the light cone. Hence, under this theory, one can perform Fock expansion where the valence and the nonvalence quark components are fully covered [129, 130]. Hence, there is much information on the partonic structure of hadrons stored in the pions produced during hadronization that can be further exploited. These processes give an insight into the transition to the hypothetical QGP phase of matter where deconfinement and chiral symmetry restoration are expected. These studies gained more attention when ultrarelativistic heavy ion collision experiments in BNL Brookhaven and CERN Geneva were able to create matter under extreme conditions of temperature and densities necessary for phase transition [131]. Pion-pion annihilation is particularly necessary since pions are the second most abundant particles produced during the heavy ion collision from hot and dense hadron gas. Such annihilation also results in the production of photon and dilepton which form the basics for probing quark-gluon plasma phase and further hadronizations. These particles are known to leave the dense matter phase quickly without undergoing strong interactions [1012]. However, the leptons are capable of annihilating into quarks and gluons (partons) in a process popularly referred to as gluon bremsstralhaung [16, 17]. The free parameters of the model are , , , , , , and , where is the degeneracy of the gluons and its acceptable value for massless gluons in representation is while its value under representation is . It is known from QCD lattice simulation that ; therefore, we can determine from Equation (28) that [4345, 132, 133] which is precisely the scalar glueball mass. For zero degeneracy, i.e., , the critical temperature Equation (34) goes to infinity and all the glueballs get melted leaving pure gluons [75]. We determined the gluon mass to be in the model framework. This value lies within the values obtained from lattice simulation as referred to in Section 3.3. Also, is the magnitude of the QCD vacuum at the ground state. It has different values depending on the model under consideration. Its value estimated from the sum rule for gluon condensate is [134], and the Bag model is also quoted as . From the model , recalling that is related to and as , has the dimension of energy and is dimensionless. Again, when we combine Equations (28) and (37), we can express the string tension as a function of temperature :

Here, we retrieve Equation (28) at , , and at , we move into the quark-gluon plasma (QGP) phase where the particles interact in a disorderly manner [46, 135, 136]. Thus, the potential in Equation (27) can be expressed in terms of temperature as where we recover Equation (27) at . Thermal deconfinement and the restoration of chiral symmetry occur at [117, 118] due to the dissolution of the glueball mass, and we have the QGP phase at [46, 135, 136] where the glueball mass becomes unstable. As established in group representation, it is known from the hadron spectrum that the chiral symmetry is spontaneously broken due to the nonperturbative dynamics of the theory, leading to color confinement. Considering quark representation in , the color charges may be screened by the gauge degrees of freedom. Under this representation, we can possibly have quark-gluon color singlet thermionic states, quark-antiquark color singlets, or gluon-gluon color singlets. Generally, in this paper, we investigate glueballs which will fall under the gluon-gluon color singlet states.

5.2. Conclusion

The process of high energy annihilation of to the production of through the interaction during hadronization in the high particle density region is briefly addressed. Thus, we focus on the discussion to cover the low energy regime where symmetry is spontaneously broken through the Abelian Higgs mechanism to give mass to the resulting glueballs. The scalar field in this case plays the role of the Higgs field, which undergoes modification into glueballs with mass . This leads to color confinement of glueballs. We explored the behaviour of the scalar field potential and the glueball potential with the decay constant , and the results are presented in Figures 3 and 4. The model explains other confining (IR) properties such as gluon condensation, glueball mass, gluon mass, string tension, and dual superconductivity through monopole condensation. The results for color confinement, screened electric (magnetic) field due to the monopoles, and the gluon condensate are presented in Figures 15, respectively. The glueball fields acquire their mass through SSB, and their mass remains the most relevant parameter throughout the confinement properties. The magnitude of the glueball mass is precisely the same as the observed lightest scalar glueball mass . This mass is responsible for the string tension that keeps the particles in a confined state and appears in the monopole condensate as well. It also appears in the gluon condensate , capable of hadronizing to form light pions in the low energy regime. In this paper, two forms of hadronization were discussed; the first one takes place at the high energy regime where the individual quarks and gluons recombine to form hadrons/pions. The second is at the low energy regime where the quark and antiquark pairs/valence gluons that form pions/glueballs get separated due to long separation distances, and the string tension that holds them together breaks leading to hadronization. We also investigated the effect of temperature on the effective scalar glueball mass , the gluon mass , string tension , confining potential , gluon condensate , the penetration depth , the superconducting monopole densities , and the normal monopole density , and their results are presented in Figures 611, respectively. We calculated the QCD -function and the strong “running” coupling through renormalization group theory to enhance the discussions on gluon mass generation. Additionally, the glueball mass and the gluon masses were calculated, and the outcome was compared with lattice simulation results, analytical study, or phenomenological analyses to ascertain their reasonability. Finally, the model produces the correct behaviour of confining potential at , where the potential has linear growth with , Equation (27), consistent with the Cornell potential model. The potential keeps growing linearly for with thermal deconfinement phase at . However, most of the results on QCD lattice calculations point to a temperature correction of order to the string tension [137139]; meanwhile, the model proposes correction in order of in Equation (83). This can be regularized at to obtain the correct order . Indeed, we did not find such regularization necessary because there are some proposed phenomenological models that suggest correction [140] to the string tension as well. But we deem it necessary to provide some explanations based on the model. In the previous paper Ref. [40], we obtained a correction of order as elaborated in Equation (63). We observed that in Ref. [40], the temperature correction to the string tension came from the scalar field as elaborated in the latter part of Section 3.3. On the other hand, the temperature correction in the gauge field leads to a correction in the order of as discussed in Section 3.2 and elaborated in Section 5.1 in Equations (83) and (84).

Data Availability

The data used to support the findings of this study have been made already available in the literature.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

Acknowledgments

We would like to thank CNPq, CAPES, and CNPq/PRONEX/FAPESQ-PB (Grant no. 165/2018), for partial financial support. FAB also acknowledges support from CNPq (Grant no. 312104/2018-9).