Confinement of Fermions in Tachyon Matter at Finite Temperature

We study a phenomenological model that mimics the characteristics of QCD theory at finite temperature. The model involves fermions coupled with a modified Abelian gauge field in a tachyon matter. It reproduces some important QCD features such as confinement, deconfinement, chiral symmetry, and quark-gluon-plasma (QGP) phase transitions. The study may shed light on both light and heavy quark potentials and their string tensions. Flux tube and Cornell potentials are developed depending on the regime under consideration. Other confining properties such as scalar glueball mass, gluon mass, glueball-meson mixing states, and gluon and chiral condensates are exploited as well. The study is focused on two possible regimes, the ultraviolet (UV) and the infrared (IR) regimes.


Confinement of heavy quark states
QQ is an important subject in both theoretical and experimental studies of hightemperature QCD matter and quark-gluon-plasma phase (QGP) [1]. The production of heavy quarkonia such as the fundamental state of cc in the Relativistic Heavy Iron Collider (RHIC) [2] and the Large Hadron Collider (LHC) [3] provides basics for the study of QGP. Lattice QCD simulations of quarkonium at finite temperature indicate that J/ψ may persist even at T = 1:5T c [4], i.e., a temperature beyond the deconfinement temperature. However, in simple models such as the one under consideration, confinement is obtained at T < T c , deconfinement and QGP phase at T ≥ T c . Several approaches of finite temperature QCD have been studied over some period now, but the subject still remains opened because there is no proper understanding of how to consolidate the various different approaches [5][6][7][8][9][10]. Ever since, the Debye screening of heavy QQ potential was put forward [11], leading to the suppression of quarkonia states such as J/ψ. Attention has been directed to investigations to under-stand the behavior of QQ at the deconfined phase by nonrelativistic calculations [12][13][14][15][16] of the effective potential or lattice QCD calculations [4][5][6][7][8]17].
We will study the behavior of light quarks such as up (u), down (d), and strange (s) quarks and also heavy quarks such as charm (c) and other heavier ones with temperature (T). We will elaborate on the net confining potential, vector potential, scalar potential energy, string tension, scalar glueball masses, glueball-meson mixing states, and gluon condensate with temperature. Additionally, confinement of quarks at a finite temperature is an important phenomenon for understanding the QGP phase. This phase of matter is believed to have formed few milliseconds after the Big Bang before binding to form protons and neutrons. Consequently, the study of this phase of matter is necessary for understanding the early universe. However, creating it in the laboratory poses a great challenge to physicists, because immensely high energy is needed to break the bond between hadrons to form free particles as it existed at the time. Also, the plasma, when formed, has a short lifetime, so they decay quickly to elude detection and analyses. Some theoretical work has been done in tracing the signature of plasma formation in a series of heavy ion collisions. Nonetheless, a phenomenological model that throws light on the possible ways of creating this matter is necessary. In this model, the QGP is expected to be created at extremely high energy and density by regulating temperature, T > T c . Besides these specific motivations for this paper, all others pointed out in [18,19] are also relevant to this continuation. This paper is meant to complement the first two papers in this series [18,19]. Particularly, it is a continuation of [18]. We will try as much as possible to maintain the notations in [18] in order to make it easy to connect the two. That notwithstanding, we will redefine some of the terms for better clarity to this paper when necessary. Generally, we will investigate how thermal fluctuations affect confinement of the fermions at various temperature distribution regimes. The Lagrangian adopted here has the same structure as the one used in [18], so we will not repeat the discussions.
Even though strong interactions are known to be born out of non-Abelian gauge field theory, no one has been able to compute confining potentials through this formalism. Therefore, several phenomenological models have been relied upon to describe this phenomenon. Common among which are linear and Cornell potential models. Bound states of heavy quarks form part of the most relevant parameters for understanding QCD in high-energy hadronic collisions together with the properties of QGP. As a result, several experimental groups such as CLEO, LEP, CDF, D0, Belle, and NA5 have produced data, and currently, ongoing experiments at BaBar, ATLAS, CMS, and LHCb are producing and are expected to produce more precise data in the near future [20][21][22][23][24][25][26].
The paper is organized as follows: In Section 2, we present the model and further calculate the associated potentials and string tensions in the ultraviolet (UV) and the infrared (IR) regimes. In Section 3, we present the technique for introducing temperature into the model. We divide Section 4 into three subsections where we elaborate on vector potential in Section 4.1, scalar potential energy in Section 4.2, the quark and gluon condensates in Section 4.3, and chiral condensate in Section 4.4. In Section 5, we present our findings and conclusions.

The Model
We start with the Lagrangian density: with a potential given as The Lagrangian together with the potential and tachyon condensation governed by the color dielectric function GðϕÞ and dynamics of the scalar field produces the required QCD characteristics.
The equations of motion are where μ, ν = 0, 1, 2, 3. We are interested in studying static confining potentials that confine the electric flux generated in the model framework, leaving out the magnetic flux which will not be necessary for the current study. With this in mind, we choose magnetic field components F ij = 0, such that F μν F μν = −F j0 F j0 = 2E 2 and the current density j ν = ψqγ ν ψ restricted to charge density, i.e., j ν = ðρ, 0 ! Þ, where ρ = q/ðð4/ 3ÞπR 3 Þ. Thus, in spherical coordinates, Equation (5) becomes where λ = q/4πε 0 . Expanding Equation (3) in spherical coordinates, we obtain 1 we set GðϕÞ = ð2πα′Þ 2 VðϕÞ and ignore the term in the order λ 2 . Shifting the vacuum of the potential such that ϕðrÞ ⟶ ϕ 0 + η, where ηðrÞ is a small perturbation about the vacuum ϕ 0 = 1/α. We can express Equation (7) as We can also determine the bosonic mass from the relation Also, we define the Dirac delta function such that where R is the radius of the hadron and r is the interparticle separation distance between the quarks. Now, developing the Laplacian in Equation (8) for particles inside the hadron gives Advances in High Energy Physics here, we have set A similar expression can be seen in [27,28], where the mass of the tachyonic mode is determined to be inversely related to the Regge slope, i.e., M 2 ∝ 1/α′ and R ∝ ffiffiffiffi ffi α′ p . Accordingly, m 2 ϕ can be identified as the scalar glueball mass obtained from the bound state of gluons (quenched approximation) through the tachyon potential VðϕÞ. Thus, the tachyonic field ϕ determines the dynamics of the gluons. Expanding Equation (11) leads to where K = qm 2 q − m 2 ϕ . This equation has two separate solutions representing two different regimes, i.e., the infrared (IR) and the ultraviolet (UV) regimes: for r ≤ R,respectively. We can now determine the confinement potential from the electromagnetic potential where E r is the electric field modified by the color dielectric function GðrÞ. From Equation (6), we can write With the shift in the vacuum, ϕ ⟶ η + ϕ 0 , the tachyon potential and the color dielectric function take the form and the net confining potential becomes Here, we choose the negative part of Equation (15) corre-sponding to the potential of antiparticle; consequently, we choose q = −1 corresponding to an anticolor gluon charge. Besides, confinement in this regime increases with interquark separation r until a certain critical distance r = r * ; beyond r * , hadronization sets in. Since λ depends on the charge q, we can choose λ = −1. Noting that ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi 2ðm 2 ϕ + m 2 q Þ q r ≪ 1, the string tension can be determined from Equation (18) as where m 2 L can be identified as the glueball-meson mixing state obtained from unquenched approximation [29][30][31][32]. Now, using σ c~1 GeV/fm as adopted in [18,19] and m 2 ϕ = 0:99 GeV 2 corresponding to m ϕ = 995 MeV, while r * = 1/ ffiffiffiffiffiffiffi 2σ c p [33]. This glueball mass is closer to the scalar glueball mass determined for isoscalar resonance of f 0 ð980Þ reported in the PDG as ð990 ± 20Þ MeV [34]. This choice of glueball mass was considered appropriate because we expect m q to be small but not negative or zero. This choice leads to m q = 100 MeV; this is also within the range of strange (s) quark mass. In the PDG report, m s = ð92:03 ± 0:88Þ MeV, but much higher values up to 400 MeV are obtained using various fit approaches [35][36][37], and m s ð1 GeVÞ = ð159:50 ± 8:8Þ MeV are obtained using various sum rules [38,39]. In this model, the hadronization is expected to start when r * = 1/ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi 2ðm 2 ϕ + m 2 q Þ q = 0:71 fm, but for string models such as this, r ≫ r * [33].
On the other hand, the potential in the UV regime becomes Here, we chose the positive part of the potential, i.e., a potential of a particle, so we set λ = q = 1. Also, with the string tension in both cases, we have set the integration constants (c) to zero (c = 0) for simplicity. Assuming that σ s~1 GeV/fm as used in [18], we can determine m s = 1:73 GeV, precisely the same as the lightest scalar glueball mass of quantum number J PC = 0 ++ [34,[40][41][42][43][44], and the critical distance is r * s = 1/ ffiffiffiffi ffi σ s p . The string tension σ is related to the glueball mass by mð0 ++ Þ/ ffiffiffi σ p = const:, where the magnitude of the constant is dependent on the approach adopted for the particular study. In this model framework, we have two different ratios, one in the IR regime and the other in the UV regime, i.e., m L / ffiffiffiffi ffi σ c p = 1 and m s / ffiffiffiffi ffi σ s p = 1:73, respectively.
Hence, the ratio is dependent on the regime of interest. The string tension and the dimensionless ratio are precisely 3 Advances in High Energy Physics known physical quantities in QCD, because their corrections are known to be the leading order of Oð1/N 2 c Þ in SUð∞Þ limit. These quantities have been determined by lattice calculations, QCD sum rules, flux tube models, and constituent glue models to be exact [45][46][47].
From the afore analyses, we can determine m q = 2 GeV, while a bound state of quark-antiquark pair would be m 2 q = 4 GeV 2 . This is compared with the mass of a charm quark m c ð m c Þ = ð1:280 ± 0:025Þ GeV as reported by PDG [34] and a slightly higher values up to m c = 1:49 GeV in various fits [35][36][37]. Of cause, the model can take any quark mass greater than one (m q > 1) but less than some critical mass m. Hence, it could be used to successfully study chamonia and bottomonium properties. The only restriction here is that m ϕ < m q ≤ m in the UV regime for the model to work efficiently. In fact, the major results in this section, i.e., the potentials and the string tensions, were reached in [18] and thoroughly discussed at T = 0. So, an interested reader can refer to it for step-by-step computations and justifications.
In the model framework, we developed two different potentials, linear confining potential in the IR and a Cornell-like potential in the UV regimes. The potentials obtained from the UV and the IR regimes have confining strengths, m 2 s and m 2 L , respectively. Linear confining models have been studied by several groups [48][49][50][51][52][53][54][55][56][57][58][59], and its outcome establishes a good agreement with experimental data for quarkonium spectroscopy together with its decay properties. This potential is motivated by the string model for hadron, where quark and an antiquark pair are seen to be connected by a string that keeps them together. In spite of the successes of this model, it fails to explain the partonic structure of the color particles as observed from deep-inelastic scattering experiments. Consequently, since the QCD theory has generally been accepted as the fundamental theory governing strong interactions, the linear confining model gives a comprehensive description of some sectors of the theory, identified as the IR regime [60,61]. The Cornell potential, on the other hand, calculated in the UV regime of the model is motivated by lattice QCD (LQCD) calculations [62][63][64][65][66][67][68][69][70][71][72]. Generally, it is represented as where e and σ are positive constants which are determined through fits. It consists of Coulomb-like part with strength e = 0:5 due to single gluon exchange and the linear confining part with strength σ = m 2 s /3 in the model framework [33,73,74]. It is a useful potential in QCD theory because it accounts for the two most important features of the theory, i.e., color confinement and asymptotic free nature of the theory. It reproduces the quarkonium spectrum as well. Also, e and σ can be determined in the framework of heavy quark systems through fitting to experimental data [48][49][50]75].
Moreover, the IR model takes in quark masses in the range of 0 to 1 GeV, while the UV model takes in masses greater than ffiffi ffi 3 p ; consequently, there is intermediate mass that can be exploited between 1 and ffiffi ffi 3 p GeV. This mass range is expected to show hadronization in the IR regime and an asymptotic free behavior in the UV regime without a stable linear confining contribution (QGP). Likewise, the energy suitable for investigating the IR properties is 1:4 GeV and below, deducing from the critical distance r * , while in the UV regime the maximum energy for confinement is 1 GeV and below judging from the magnitude of r * s , but the process takes place at enormously high energy r ⟶ 0, where Coulomb potential contribution is dominant. As a result, the model does not account for the behavior of the particles at the intermediate energy regime, apart from knowing that the particles are "asymptotically free" at r ⟶ 0 and confined at r ⟶ r * s and beyond. How the particles behave at the intermediate energy regimes is beyond the scope of this paper. Additionally, perturbative QCD (PQCD) involves hard scattering contribution which requires higher energy and momentum transfer in order to take place. This is referred to as the "asymptotic limit"; thus, the initial and final states of the processes are clearly distinct. Therefore, it has been widely argued that the current energy regime for exploring the QCD theory describes the "subasymptotic" regime, calling for a revision in pure perturbation treatment of the theory. In this case, new mechanisms or models are required to investigate the dynamics of elastic scattering in the intermediate energy region [76][77][78][79][80]. One way to investigate this energy region ("soft" behavior) is to use dispersion relation together with operator product expansion (quark-hadron duality) which are collectively referred to as the QCD sum rule [81][82][83]. Thus, we can study the IR behavior at the lower end of the intermediate energy region while the UV regime lies within the upper end of the intermediate energy region to infinity. The CEBAF (Continuous Electron Beam Accelerator Facility) at Jefferson LAB was initially built to experimentally investigate these intermediate energy regions before it was recently upgraded to cover the high-energy regime.

Thermal Fluctuations
In this section, we discuss and propose how temperature can be introduced into the model. Since we have determined in the previous section that the string tension which keeps the quark and antiquark pair in a confined state is dependent on the glueball-meson mixing states m L and m s [29,32], we can determine the confining potentials at a finite temperature if we know exactly how the glueball-meson states fluctuate with temperature. We calculate the temperature fluctuating glueball-meson states directly from Equation (1). Here, we will write the equations in terms of the glueball field η, redefined in dimensionless form such that χ = η/ϕ 0 for mathematical convenience. The difference between the gluon field and the glueball fields η is that the glueball fields are massive while the gluon fields are not. In this paper, the field ϕ that describes the dynamics of the gluons contains a tachyonic mode; when the tachyons condense, they transform into glueballs η with mass m ϕ . Accordingly, 4 Advances in High Energy Physics [84], where the factor ð2πα′Þ 2 was introduced to absorb the dimensionality of VðηÞ and render GðηÞ dimensionless as explained above. Besides, the theory of gluodynamics is classically invariant under conformal transformation giving rise to vanishing gluon condensation hF μν F μν i = 0 [18,19]. However, when the scale invariance is broken by introducing quantum correction, say −jε v j, the gluon condensate becomes nonvanishing, i.e., hF μν F μν i ≠ 0. This phenomenon is referred to as the QCD energy momentum tensor trace (θ μ μ ) anomaly [85]. In this paper, if one disregards the fermions and considers only the gluon dynamics in Equation (1), the potential VðϕÞ breaks the scale symmetry and brings about gluon condensation. This phenomenon has been elaborated below in Section 4.3.
To determine the thermal fluctuations, we need to define the field quanta distribution of the mean gauge field hF μν F μν i, hΔ 2 i with χ = χ + Δ and the spinor fields h ψψi. We find that the fluctuating scalar glueball mass m * 2 ϕ does not directly depend on the glueball field χ; hence, no correction to the scalar field will be required. The field quanta distribution of the fields is 2 where n B ðE B Þ = ðe βE A − 1Þ −1 is the Bose-Einstein distribution function and n F ðE ψ Þ = ð1 + e βE ψ Þ −1 is the Fermi-Dirac distribution function. Also, ν and ν q are the degeneracies of the gluons and the quarks, respectively, and β = 1/T. We can now analytically solve these integrals by imposing some few restrictions. We assume a high-energy limit, such that E ≈ p for c = 1 = ℏ corresponding to a single particle energy E 2 = m * 2 + p 2 ; therefore, we find [86,87] which results in 3 and we find the quark condensate 4 Here, we used the transformation x = E/T and the critical temperatures However, from Equation (23), we can define the dimensionless quantity α 2 ϕ 2 0 ð2πα′Þ 2 =g 2 ; as a result, Now, substituting the thermal averages in Equations (26) and (28) of the fields into the above equation while we absorb the second term derived from the potential into the definition of the first and the third terms at their ground states, To proceed, we use the standard definition for determining the QCD vacuum: Using this expression and Equation (2), we obtain maintaining the term with dependence on ϕ and discarding the constant, the vacuum becomes B 0 = m 2 ϕ ϕ 2 0 /2 [84]. Noting that hF μν F μν i 0 = B 0 represents the vacuum gluon 5 Advances in High Energy Physics condensate, we can express here, we setg 2 = 2, hF μν F μν i 0 /ϕ 2 0 ≈ m 2 ϕ /2 and h ψψi 0 /ϕ 2 0 ≈ m q /2. The string tension is determined by QCD lattice calculations to reduce sharply with T, vanishing at T = T c representing melting of hadrons [88][89][90]. At T = 0, we retrieve the result for glueball-meson mixing state m 2 L in Equation (19), i.e., We can also retrieve the result in the UV regime m 2 s if we set m 2 ϕ ⟶ −m 2 ϕ in Equation (34) at T = 0. The nonperturbative feature of QGP is accompanied by a change in the characteristics of the scalar or the isoscalar glueballs and the gluon condensate [91]. This is evident in lattice QCD calculation of pure SUð3Þ c theory, pointing to sign changes [86,87]. Using Equations (17) and (34), the glueball potential can be expressed as Accordingly, we can express Equations (18), (19), (20), and (21) in terms of temperature as with string tension with string tension As the temperature is increasing, the confining part shows some saturation near T ≈ T cψ ≈ T cA ≈ T c and vanishes completely at T = T c indicating the commencement of deconfinement and the initiation of the QGP phase. This weakens the interaction of the particles such that the string tension that binds the Q and Q reduces and eventually becomes asymptotically free at T = T c . Beyond the critical temperature T > T c , we have the QGP state where the quarks behave freely and in a disorderly manner [92,93]. The possibility of studying the QGP state in detail with this model exists, but that is beyond the scope of this paper. When we set σ c ðm q , TÞ = 0 and calculate for a common critical temperature by assuming T cψ = T cA = 1 and m 2 ϕ = σ c − m 2 q , we will have where in the last step we substitute m 2 ϕ = σ c ð0Þ − m 2 q , and in solving the equation, we bear in mind that σ c ≃ 1 GeV/fm. Therefore, consequently, On the other hand, the common critical temperature in the UV regime can be calculated from Substituting σ s ≃ 1 GeV/fm and solving the equation accordingly, It has been found that the bound state of the charmanticharm state dissolves at T = 1:1T c [94]. However, the question as to whether heavier quark bound states dissolve at T = T c [95][96][97] or temperatures higher than 6 Advances in High Energy Physics deconfinement temperature T > T c [4,6,10] still remains open. These two separate pictures have informed different phenomenological models based on confinement and deconfinement transitions to QGP states to explain the observed suppression of J/ψ produced in the RHIC. We present a simple model based on the projection that the QQ bound state melt at T = T c .

Vector and Scalar Potentials and Gluon and
Chiral Condensates 4.1. Vector Potential. To determine the vector potential, we solve Equation (11) outside the hadron, i.e., δð r ! Þ = 0, so for IR and UV regimes, respectively; these solutions are equivalent to the solutions in Equation (14) at m q = 0. Now, substituting the solution at the left side of Equation (48) into (16) and into (15), we can determine the vector potential to be with string tension In terms of temperature fluctuations, the potential can be expressed as bearing in mind that the corresponding temperature fluctuating string tension [98] becomes These results can also be derived from Equation (34) for m q = 0.

Scalar Potential Energy.
The scalar potential energy [99][100][101] for confinement is calculated by comparing Equation (4) with the Dirac equation where b α and b β are Dirac matrices and SðrÞ, which is well defined inside and on the surface of the hadron, and zero otherwise is the scalar potential (for detailed explanations, see [18]). Hence, the scalar potential obtained by comparing Equations (4) and (54) can be expressed as For an antiparticle, we set q = −1, so we can rewrite the scalar potential in terms of the temperature as S r, m q , T À Á = m q m * 2 ϕ r 2 cosh 2 ffiffi ffi We can now write the net potential energy by adding the vector potential energy and the scalar potential energy; hence, [102,103], and in terms of temperature, we can express 4.3. Gluon Condensates. We calculate the energy momentum tensor trace, θ μ μ , from the relation Moreover, the energy momentum tensor trace θ μ μ for the classical QCD chiral effective Lagrangian [93,104] is given as where βðgÞ = −bα s /ð4πÞ (b = 11 depicts pure gluodynamics) is the QCD β-function, m f is the current quark mass matrix, and q f is the quark field. SimplifyingṼ eff in Equation (60), we find so Accordingly, comparing this with (62), we can identify ∑ f m f = m 2 L (glueball-meson mixing mass), ∑ f q f q f = η 2 (glueball field), and bα s /ð4πÞ = GðηÞ = −βð1/r 2 Þ QCD βfunction. To determine the strong running coupling, we note that GðηÞ = ηG ′ ðηÞ/2, so using the renormalization group theory [105] β we can deduce consequently, the strong running coupling becomes α s ð1/r 2 Þ = GðηÞ/2, where we can relate Q ⟶ 1/r. Using the solution at the right side of Equation (14), we can express Noting that Q 2 ≡ −q 2 , where Q 2 is the space-like momentum and q 2 is the four-vector momentum. To eliminate the Landau ghost pole that occurs at q 2 ⟶ 0, we assume that there is dynamically generated "gluon mass" at q 2 ⟶ 0, i.e., q 2 ≅ m 2 A [106,107]. Now, relating the results in Equation (60) with the standard vacuum expectation value of QCD energy momentum tensor trace, [18,19], we can writẽ Now, we rescaleṼ eff with the QCD vacuum energy density −jε v j such thatṼ eff ⟶ −jε v jṼ eff , and using the potential of the glueball fields defined in Equation (17), we obtain here, we usedṼ eff = 1/4ðm 2 ϕ + qm 2 q Þη 2 = ðm 2 L /4Þη 2 for q = 1. Also, at the classical limit ε v ⟶ 0, we have vanishing gluon condensate hF μν F μν i = 0. The vacuum expectation value (VEV) of the gluon condensate was determined in [108][109][110] using Yang Mills theory with an auxiliary field ϕ, where the fluctuations around ϕ give rise to the glueball mass, mð0 ++ Þ. It was observed that in the nonvanishing VEV, the gluon acquire a common mass, m A , and the ratio mð0 ++ Þ/m A in the leading order is The nonvanishing condensate was attributed to the pair condensate of transverse gluons. Juxtaposing that to the results in Equation (70), we can identify Other calculations based on QCD theory [111] also show that the nonvanishing gluon condensate in the 8 Advances in High Energy Physics absence of fermions can be expressed in terms of the gluon mass as where N is the number of colors and α s is the strong coupling constant. This relation also show a proportionality between gluon mass and the nonvanishing condensate-more discussions relating to this subject can also be found in [112,113]. Again, from the left side of Equation (70), we can deduce the expression hðGðηÞ/2ÞF μν F μν i = hα s ð1/r 2 ÞF μν F μν i; hence, the expression in the angle brackets at the right side of Equation (70) is equivalent to Equation (67). Following the discussions at the latter part of Section 2, m L = 1 GeV so we can determine the gluon mass as m A = 500 MeV. This can be compared to the result determined from QCD lattice simulation, projecting the gluon mass to be m A = 600~700 MeV [114][115][116]. Some heavier gluon masses have also been determined closed to~1 GeV using phenomenological analysis [117,118] and some QCD lattice studies [119]. Quark and gluon condensates are responsible for confinement, glueball formation, and hadron mass formation [80,104,120,121].
The QCD vacuum at the ground state enables us to study the characteristics of the QGP, dynamics of phase transitions, and hadronization. All these properties are a result of the nonperturbative nature of QCD theory in some regimes (IR) and cannot be studied using the usual perturbative QCD theory.
To determine the temperature fluctuations in the gluon mass m A and the gluon condensate, we need to correct the temperature in the glueball field χ. The gluon mass is formed due to screening of the gluons in the vacuum or at a temperature where the constituent quarks dissolve into gluons [122]. The temperature can be introduced by defining the glueball field such that χ = χ + Δ, so η = ϕ 0 ð χ + ΔÞ, with the restriction hΔi = 0 to avoid the occurrence of cross-terms in the thermal average. Again, χ represents the mean glueball field, and the angle brackets represent the thermal average. We express the fluctuation hΔ 2 i in terms of field quanta distribution as defined for the gauge and the fermion fields in Equations (26) and (27), respectively, i.e., Reverting to Equation (3) and the definition η = ϕ 0 ð χ + ΔÞ, we can determine the thermal average directly from the equation of motion to be clearly, η = χ = 0 is a solution; hence, hη 2 i = ϕ 2 0 hΔ 2 i.

Moreover, substituting Equation (74) for temperature fluctuations in η 2 into Equation (70), the gluon condensate becomes
The QCD vacuum at an extremely high temperature and density has an important significance in both elementary particle physics and cosmology. In cosmology, it can be used to explain the evolution of the universe at extremely higher temperatures when matter becomes super dense and hadrons dissolve into a "soup" of their constituents-quarks and gluons. At such superdense matter regimes, the quarks become very close and asymptotically free, so the quarks are no longer confined into hadrons. It is important to note that the energy density is related to temperature as ρ ∝ T 4 , so the density increases with temperature. This phase of matter was proposed [123,124] after the asymptotic free [125,126] nature of the QCD theory [127] was determined. Currently, almost all the accepted models for the cooling of the universe are based on phase transitions, either first or second order of spontaneous symmetry breaking of fundamental interactions [128][129][130]. On the other hand, the standard model (SM) of elementary particle physics suggests two of such phase transitions [130]. It shows an electroweak symmetry (EW) breaking at temperatures in the order of T~100 GeV which generates mass to elementary particles such as quarks. It is also associated with the observed baryon-number-violation of the universe [131]. This transition is determined to be an analytical crossover in lattice simulations [132].
Again, spontaneous chiral symmetry breaking is next; it is known to occur at temperatures in the order of T~200 GeV. Nevertheless, we have hadronic matter below this temperature, and above it, we have an expected transit into the QGP phase. This phase is also important in understanding the evolution of the early universe. Some references on this subject using nonperturbative lattice QCD models can be found in [133,134]. Particularly, knowing that the baryon chemical potential μ B is expected to be smaller than the usual hadron mass, μ B ≈ 45 MeV at ffiffiffiffiffiffiffi ffi s NN p = 200 GeV [135], and almost vanishing μ B ≈ 0 GeV in the early universe. It is suitable to model inflation of the early universe at high temperatures and low baryon densities. There is strong evidence that confinement of quarks into hadrons is a low-energy phenomenon [133,136,137], but it strongly suggests that QCD phase transition is a crossover. Numerically, it has been established that at μ B = 0 GeV, the two-phase transitions possibly coincide, i.e., deconfinement and chiral symmetry restoration rendering the chiral effective theory invalid [138]. From our model framework, Equation (76), we obtain confinement at T < T c , deconfinement and restoration of chiral symmetry at T = T c , and QGP phase at T > T c .

Advances in High Energy Physics
Deducing from Equation (76), the thermal fluctuating gluon mass becomes where g 2 = m 2 L /48ϕ 2 0 is a dimensionless coupling constant [122,139]. This results look similar to the Debye mass obtained from the leading order of QCD coupling expansion at higher temperatures and zero chemical potential, m D = gT [140]. This mass is obtained from the lowest order perturbation QCD theory [141], but it is known to be resulting from the IR behavior of the theory [142]. This mass occurs at temperatures greater than the deconfinement temperature T c due to the melting of QQ bound states [11,140]; generally, it is referred to as Debye screening mass of color charges.
Furthermore, the experimental evidence of QCD vacuum condensates dates back from the pre-QCD era 1950-1973 to post-QCD after 1974, and it still remains an issue of interest. The QCD vacuum condensate informs color confinement and chiral symmetry breaking which are fundamentals of strong interactions. The condensate is a very dense nonperturbative state of matter consisting nonvanishing gluon and quark condensates, in theories containing effective quarks, interacting in a haphazard manner [80]. Since quarks and gluons are not observed directly in nature, these characteristics are difficult to observe experimentally; only colorless (baryons) or color neutral (mesons) hadrons are observed. However, much is known about the QCD vacuum, credit to the fast-evolving QCD sum rule which explores the gauge invariant formalism to explain the nonperturbative nature of the condensates [143,144]. That notwithstanding, the QCD sum rule gives good results at the intermediate energy region and poses some short falls at large distance regimes where color confinement and chiral symmetry breaking are more pronounced. So, a resort is made to the Vacuum Correlator Method (VCM) to give a comprehensive description of all the possible QCD phenomena involved ( [145][146][147][148]. The VCM is based on gauge invariant Green's function of colorless or color neutral objects expressed in path integral formalism using field correlators instead of propagators.

Chiral
Condensate. Chiral symmetry breaking is one of the known physical properties of the nonperturbation regime of QCD theory besides the famous color confinement. Mass splitting of the chiral partner observed in the hadron spectrum and Goldstone bosons (π ± and π 0 ) which appear due to spontaneous symmetry breaking leading to color confinement are strong evidence of chiral symmetry breaking [149][150][151] in QCD vacuum [152]. The presence of quark mass breaks the chiral symmetry explicitly. With these background and other theoretical considerations, it is believed that there is chiral condensate in the vacuum which is proportional to the expectation value of the fermionic operator or the quark condensate h ψψi. However, increasing temperature augments the thermal excitations of the hadrons due to an increase in the density (ρ ∝ T 4 ) of their quark constitu-ents. This decreases the vacuum condensate until it eventually vanishes at a critical temperature T cq resulting into restoration of the chiral symmetry. At this temperature, the hadrons undergo phase transition to deconfinement and quark-gluon plasma (QGP) phase. Increasing baryon density has the same consequences on the QGP phase and chiral condensate similar to temperature increase [153,154]. To begin the calculation of the chiral condensate, one needs to take into account the dependence of the hadron masses on the current quark mass m q . Some known first principle approaches to this subject that give a consistent account of hadrons and their quark mass dependent are Chiral Perturbation Theory (ChPT) [155], lattice QCD (lQCD) [156], and Dyson-Schwinger Equation (DSE) [157,158].
However, the dominant degrees of freedom for QCD at low energies are hadrons [159,160]. Particularly, pions and kaons interact weakly according to Goldstone theory; hence, they can be treated as free particles. Consequently, the standard relation for calculating the chiral condensate, h qqi, can be derived from Equation (1) as [155], where h qqi is significant for spontaneous chiral symmetry breaking (SχSB) and the angle brackets represent thermal average. Subtracting the fundamental vacuum condensate h qqi 0 to ensure that the condensate has its maximum T = 0 [161], we can express the above relation as here, details of the vacuum condensate h qqi 0 [162,163] are of no importance to the analyses. All the information needed to study the contributions of the hadrons is contained in the medium dependent term. From Equation (1) and using the expression in Equation (17) derived from ϕðrÞ ⟶ ϕ 0 + η, we obtain Thus, the chiral condensate is proportional to the glueball potential GðηÞ which has its maximum condensate at GðηÞ ⟶ 0 [18,19]. The quark density h ψψi is a measure of the strength of the condensate and SχSB. At h ψψi = 0, we have an exact chiral symmetry, deconfinement, and QGP phase, while nonvanishing quark condensate, h ψψi ≠ 0, is the regime with SχSB and color confinement [164,165]. Therefore, h ψψi serves as an order parameter that determines the phase transitions. Since the explicit pion degrees of freedom is significant for studying chiral symmetry restoration at low temperature and quark densities, h ψψi is constituted with up (u) and down (d) quarks. Eventually, protons and neutrons are also composed of 10 Advances in High Energy Physics

Advances in High Energy Physics
Additionally, in terms of temperature where in the first step, we used the results in Equations (27) and (74), bearing in mind that hη 2 i ≃ ϕ 2 0 hΔ 2 i and Thus, at T = T cq , the chiral symmetry gets restored in the model framework. If one intends to investigate the behavior of the chiral condensation with varying quark mass m q , we can define the critical temperature as corresponding to a condensate Also, an evidence from lQCD confirms that at the confinement phase the chiral symmetry is spontaneously broken down [165,166] to the flavor group, i.e., with associated three Goldstone bosons (π ± and π 0 ) which spontaneously break the chiral symmetry. For three quark flavors, we have here, there are eight Goldstone bosons (π ± , π 0 , K ± , K 0 , K 0 , and η) involved (see References [167,168]).
Again, a highly excited state in high-energy hadronic collisions leads to the formation of disoriented chiral condensate which can later decay into ordinary vacuum through coherent emission of low momentum pions. This process is theoretically synonymous to the Higgs mechanism that leads to the release of Goldstone bosons. This leaves a signature of color confinement eventually [169]. Even though the idea of chiral condensate was speculative at its inception in the 1990's [170][171][172][173][174], it has attracted several theoretical and experimental attention subsequently. Besides its simplicity, it is also motivated by the discovery of the "so-called" Centauro events in cosmic ray [175][176][177] where clusters consisting charged and neutral pions were observed. The chiral condensate has also been studied theoretically in the light of high-energy heavy ion collisions [178,179] due to the high energies involved at the collision zone; a hot chirally symmetric state (QGP) is formed in the process. Because of the fast expanding nature of the system at the early stages, it is quenched down to a low temperature where chiral symmetry is spontaneously broken down. Several experiments have since been set up to investigate this phenomena, keys among them are the cosmic ray experiment [175,176], nucleon-nucleon collisions at CERN [180][181][182], Fermi LAB [183], particularly, MiniMAX experiment [184], nucleon-nucleon collisions at CERN SPS ( [185][186][187][188], and the RHIC [189,190]. It also forms part of the heavy ion collision programme carried out with the multipurpose detector ALICE, at LHC [191,192].

Conclusion
Following the discussions in [18] for the constituent quark masses, we can deduce that the constituent quark masses of this model are Mðr ⟶ r * Þ = 2m q = 200 MeV and Mðr ⟶ 0Þ = 2m q = 4 GeV for the IR and the UV regimes, respectively. Thus, the potential in the IR and the UV regimes can take masses within the ranges 0 ≤ m q ≤ 200 MeV and 2 ≤ m q ≤ 4 GeV, respectively. Hadronization is expected to set in, in the IR and the UV regimes for m q > 200 and m q > 4, respectively, since we have adopted the lattice simulation 13 Advances in High Energy Physics results for the string tension σ~1 GeV/fm, and noticing that σ = 1/ð2πα ′ Þ as shown in many confining string models [28,193]. The choice ð2πα′Þ 2 = 1 used throughout the paper is in order. In the framework of the model, the glueball field χ does not contribute to the fluctuations in the scalar glueball mass m * 2 ϕ , while the only candidate that contributes to the fluctuating gluon mass m * 2 A is the glueball field. On the other hand, the gauge fields, the spinor fields, and the glueball fields all contribute to the gluon condensate hðm 2 ϕ η 2 /4ÞF μν F μν i and chiral condensate. Thus, the condensates are important in understanding QCD theory but difficult to study experimentally due to the haphazard nature of interactions among these fields in the vacuum. Also, ν q and ν are the degeneracies of quarks and gluons, and they take in the values, 6 and 16 for SUð2Þ and SUð3Þ representations, respectively. These degeneracies are important in determining the critical temperatures of the model. The critical temperatures are small when the degeneracies are infinitely large, and when there is no degeneracies at all (ν, ν q ⟶ 0), the critical temperature becomes infinitely large; same is true for quarks and gluons as presented in Equations (29) and (86).
The model produces two forms of temperature corrections to the string tension, −T 2 coming from the spinors (quarks) and −T 4 from the gauge fields (gluons). The −T 2 correction to the string tension has been corroborated by some QCD lattice calculations [89,194,195] and some phenomenological models [19]. That notwithstanding, −T 4 T c1 T c2 14 Advances in High Energy Physics correction to the string tension has also been proposed by some phenomenological models [28]. Both corrections give the correct behavior of the string tension, i.e., it should reduce sharply with temperature and break or vanish at T = T c for simple models such as the one discussed here. For simplicity, we used T cA = T cψ for some of the analyses-particularly, the potentials-but there is no evidence that these two critical temperatures have the same magnitude. In any case, such assumption is informed and does not affect the results or the analyses. However, there is a discussion in [196] suggesting that T cψ > T cA (T cψ~2 70 MeV and T cA~1 70 MeV) or at least a discrepancy of about 3% reported in [197]. Using the magnitude of the string tension and the scalar glueball mass m ϕ calculated above, we obtained two different glueballmeson states corresponding to m L = 1 GeV and m s = 1:73 GeVfor the IR and UV regimes, respectively. The gluon mass was also determined as m A = 500 MeV. The critical distance for confinement in the IR regime has been determined to be r * = 0:71 fm, and its corresponding value in the UV regime is r * s = 1/ ffiffiffiffi ffi σ s p = 1 fm. Similarly, r * and r * s can be expressed as a function of temperature like the string tensions.
Some of the major results obtained are displayed on a graph to make it easy to see their behavior. Since we have extensively studied and discussed the corresponding results for T = 0 in [18], we will concentrate on the results with temperature fluctuations. We plot the glueball potential Equation (36) and its behavior with temperature in Figure 1. The confining potential Equation (37) in the IR regime for finite and infinite m q is plotted in Figure 2 and its string tension Equation (38) plotted in Figure 3. The potential in the UV regime Equation (39) is plotted in Figure 4 displaying how Cornelllike potential obtained varies with temperature for finite and infinite quark mass limits and their string tension in Equation (40) displayed in Figure 5. The gluon mass Equation (77) which possesses all the characteristics of Debye mass is displayed in Figure 6. The vector potential Equation (51) which represents chromoelectric flux confinement is displayed in Figure 7. The scalar potential energy Equation (55) and the corresponding net potential energy Equation (58) for finite and infinite m q are plotted in Figures 8 and 9, respectively. The gluon condensates calculated in Equations (70) and (76) are also displayed in Figures 10 and 11, respectively. The color dielectric function in Equation (17) represents the glueball potential. We have higher glueball condensation when GðrÞ ⟶ 0, so it follows the same discussions as Figures 8 and 9. Also, an increase in quark mass m q will lead to an increase in glueball condensation [198,199]. We explored a phase transition from the low-energy IR regime to the high-energy UV regime by studying the characteristics of T c1 and T c2 in Equations (42) and (45) displayed in Figure 12. We find that the critical temperatures decrease with an increase in quark mass thereby increasing the strength of confinement. We observed that the light quarks that are confined in the IR regime 0 ≤ m q ≤ 1 are relegated to the QGP phase in the UV regime which confines quarks with m q > ffiffi ffi 3 p . Finally, the chiral condensate h qqi is calculated in Equation (85) and displayed in Figure 13.
The glueball condensation increases with an increase in depth of the curve. Hence, the condensate increases from T = 0:9T c (black), T = 0:5 (red), to T = 0 (blue).
As it is shown in Figure 2(a), the gradient of the graph increases with decreasing temperature from T = 0:9T c (black), T = 0:5 (red), to T = 0 (blue) [26,93]. Within this temperature range, there is confinement and chiral symmetry breaking. In the right panel, we show how an increase in m q ⟶ ∞ affects the behavior of the potential and confinement of the particles. As m q is increasing, the potential increases and the confinement becomes stronger.
Here, we show the behavior of the string tension in the IR regime with varying temperatures. The string breaks quickly for light quarks while the heavier quarks have relatively longer life time even though they all vanish at T = T c1 . A regime is where all the bond states are expected to dissolve into a "soup" of their constituents. It also gives an insight into the behavior of the string tension for chiral limit m q = 0 and chiral effective m q ≠ 0 regimes.
From Figure 4(a), the gradient of the curves decreases with increasing T from T = 0 (blue), T = 0:5T c (red), to T = 0:9T c (black), indicating a decrease in binding of the quarks as the energy of the system is increasing. At T = T c (dashed), we have 15 Advances in High Energy Physics a deconfinement phase and restoration of the chiral symmetry [200,201]. Figure 5 follows the same behavior as discussed in Figure 3. However, in this regime, the string breaking is explicit as the curves intercept the T/T c2 axis at different points. That notwithstanding, the bond states of the light and the heavy quarks will dissolve at T = T c2 . Also, the fluctuating glueball mass m 2 s ðTÞ in this regime is related to the string tension by m s ðTÞ = ffiffiffiffiffiffiffiffiffiffiffiffiffi ffi 3σ s ðTÞ p . Figure 7 is similar to Figure 2 for m q = 0, T = 0 (blue), T = 0:5T c (red), and T = 0:9T c (black). This means that the gluons remain confined even if the quark mass is "removed" (m q = 0) after confinement. This is known as chromoelectric flux confinement.
Confinement is stronger with increasing depth of the curves [18,33]. Thus, the smaller the minima of the curves, the more condensed the glueballs and stronger the confinement. The minima increases from T = 0 (blue), T = 0:5T c (red), to T = 0:9T c (black).
The condensation has its maximum value at η = 0 and reduces steadily with increasing η until it vanishes.
The condensate reduces sharply with increasing temperature until it vanishes at T = T cη .
We show a transition between the IR regime corresponding to lighter quark masses to the UV regime corresponding to heavier quark masses. The nonphysical behavior of T c1 at masses m q ≤ 1 is the regime where confinement of light quarks is observed in the IR regime. The curve becomes constant near T c1 ≃ 1, the point where the string tension in the IR regime gets saturated and begins to degenerate. Consequently, the IR regime can take masses within 0 ≤ m q ≤ 1; beyond this threshold, hadronization sets in and the string tension starts decaying. On the other hand, the UV regime takes in heavier quark masses (green curve), m q > ffiffi ffi 3 p . The nonphysical behavior observed for m q ≤ ffiffi ffi 3 p corresponds to the QGP regime with negative and constantly decaying σ s . As m q increases, T c2 decreases and the confinement becomes stronger.
An increase in quark mass courses the condensate to reduce sharply and eventually vanish at T/T cq [153,162]. This behavior is opposite to that of the string tension in Figures 3 and 5 as expected in QCD lattice simulations [202]. Consequently, confinement increases with increasing quark mass while the chiral condensate decreases with increasing quark mass.

1
It is important to state at this point that the color dielectric function and the potential are related as GðϕÞ = ð2πα ′ Þ 2 V ðϕÞ [203], where α ′ is the Regge slop with dimension of length squared. So, it absorbs the dimension of the potential rendering the GðϕÞ dimensionless. But we set ð2πα ′ Þ = 1 throughout the paper for simplicity. 2 This follows the Matsubara formalism of field theory at finite temperature where space time becomes topological. Here, one makes use of the Euclidean imaginary time and solve the path integral by imposing a periodic condition on the gauge field, ω n = 2nπT, and an antiperiodic condition, ω n = ð2n + 1ÞπT, for the fermion fields for n ∈ ℤ. The imaginary time [204,205] transforms as β = 1/T, where T is temperature. 3 We used the standard integrals Here, we used

Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
The authors declare that they have no conflicts of interest.