Confinement of fermions in tachyon matter

In this paper we develop a phenomenological model inspired by QCD that mimics QCD theory. We use gauge theory in color dielectric medium ($G(\phi)$) coupled with fermion fields to produce scalar and vector confinement in chromoelectric flux tube scenario. Abelian theory will be used to approximate the non-Abelian QCD theory in a consistent manner. We will calculate vector and scalar glueballs and compare the result to the existing simulation and experimental results and projections. The QCD-like vacuum associated with the model will be calculated and its behavior studied relative to changing quark masses. We will also comment on the relationship between tachyon condensation, dual Higgs mechanism, QCD monopole condensation and their association with confinement. The behavior of the QCD string tension obtained from the vector potential of the model will be studied to establish vector dominance in confinement theories.


I. INTRODUCTION
Scalar and vector confinement [1] in 3 + 1 dimensional world has been predicted by hadron spectroscopy [2], confinement in string picture and by QCD lattice simulation but no success has been made in solving it analytically from 'first principle' of QCD. It has been shown in quarkonia phenomenology in 3 + 1 dimensional world that the best fit for meson spectroscopy are found for a convenient mixture of vector and scalar potentials. The combined vector and scalar potential has also been studied in many perspectives. Dirac equation [3] and Schwinger model [4] are among common examples employed to achieve both linear and coulomb-like potentials [5][6][7][8][9]. The attempt to confirm this result predicted by QCD lattice spectroscopy has led to the use of mass gap equation to generate mass as a dependent parameter even in quark systems with no mass current to create dynamic quarks [10]. Pure scalar and vector potentials have been dealt with in Refs. [11][12][13][14][15] as a point in focus.
In this work, we use the Lagrangian density for confinement of electric field in tachyon matter [16,17] coupled with fermion fields [17] to produce a Lagrangian density for fermionic tachyons through transformations. In this approach, the dielectric function coupled to the gauge field and the fermion mass produce the needed strong interaction between (anti-) quark pairs and the scalar field φ describes the dynamics of the tachyons. We will use a dilaton in gauge theories [18] to determine the coupling constants of the colored particles by transforming the exponential dilaton potential to conform with our chosen tachyon potential considered in this work.
We show that both scalar and vector confinements coincide with the tachyon condensation. This phenomenon is not completely new because, tachyon fields play a role similar to Higgs fields where Higgs mechanism proceeds via tachyon condensation [19]. Again, both Higgs and tachyon fields have some properties in common, they are both associated with instability or fast decaying with negative mass squared. The tachyons are expected to condense to a value of order of the string scale. Tachyons with large string couplings are considered charged, in that regard, their condensation leads to dual Higgs mechanism. This naturally translates into confinement in line with the QCD monopole condensation scenario [20]. Moreover, we will compute both vector and scalar glueball masses associated with the model. Glueballs are simply bound states of pure gluons, mixture of quark and gluon(hybrid) and multi quark bound states. The glueball spectrum has been a subject of interest for some decades now, with the focus on unraveling its states in the context of QCD theory [21].
The success of this type of confinement is based on Nambu-Goto string or the chromoelectric flux tube picture. The flux tube scenario is generally observed in static quark frame, such that there is no chromomagnetic field to induce the spin-orbit interaction of the quarks. The only interaction that is present, in this case, is the kinematics Thomas spin-orbit interaction, which is a relativistic correction. This is true for both vector and scalar confinement potentials and it is confirmed by lattice simulation results [22]. 'Thomas precession' is considered scalar, hence scalar potential is expected but it does not guarantee its connection with QCD theory. They only relate at long range spin-orbit interactions of QCD, this is precisely the infrared (IR) regime of the theory. The vector potential is achieved for short range Thomas spin-orbit interaction and the scalar potential is consistent with long range Thomas spin-orbit interaction. Both models do not depend on the quarks' spin-orbit interactions in agreement with QCD theory [23]. Though, scalar potential models have been phenomenologically accepted and used in many hadron models, it is still struggling to plant it root firmly in fundamental QCD as highlighted in [24,25]. Eventually, it was established that in a slowly moving quark frame, QCD predicts both spin dependent [26] and spin independent [27] relativistic corrections.
We will use an Abelian QED throughout the computations, but the color dielectric function G(r) modifies the gauge field and by extension the QCD vacuum of the model [28,29]. We will establish a relationship between the tachyon potential and the color dielectric function in a suitable manner. It is known that the Abelian part of the non-Abelian QCD string tension is 92%, this represents the linear part of the net potential. Thus, we can make an approximation of the non-Abelian field using an Abelian approach [30,31] -for recent development on this, see [32]. Also, it is established that if Abelian projection is followed as suggested by 't Hooft, non-Abelian QCD is reduced to an Abelian theory with charges and monopoles occurring, when such monopoles condense, confinement results. The idea of monopole condensation is a very useful one. It has been numerically shown that monopole condensation actually occurs in the confinement phase of QCD [33]. Again, the color dielectric function coupled to the gauge field F µν F µν , contains only low momentum components. Therefore, the use of Abelian approximation is justified in phenomenological QCD theory. In sum, these properties enable us to apply phenomenological field theory for QCD in this investigation to establish confinement of quarks and gluons in the infrared regime [6,[34][35][36].
The color dielectric function G is responsible for the long distance dynamics that bring about confinement in the IR regime of the model. It also facilitates the strong interactions between quarks and gluons. The scalar field φ(r) is responsible for the dynamics of the self-interacting gluon fields and the color dielectric function. We will use a Lagrangian density that describes the dynamics of the gauge, the scalar field associated with the tachyon and gluon dynamics and fermion fields coupled with mass at zero temperature [36]. The motivation for using this approach are: Firstly, we are able to compute both the scalar and vector potentials by considering a heavy anti-quark source surrounded by a relatively light and slowly moving quarks. This enables us compute the effect of short range spin dependent 'Thomas precession' (vector potential) and the effect of long range spin independent 'Thomas precession' (scalar potential). Secondly, we are able to apply phenomenological effective field theory to identify the color dielectric function with the tachyon potential in a simple form. Also, this approach makes it easy to observe how the QCD vacuum is modified by the dielectric function resulting in gluon condensation. Finally, vector and scalar glueball masses are easily calculated from the tachyon potential and the Lagrangian respectively.
The paper is organized as follows. In Sec. II we review the Maxwell's Lagrangian with source in color dielectric medium. In Sec. III we introduce the Lagrangian density for the model. In Sec. III A we derive the QCD-like vacuum of the model. In Sec. III B we choose a suitable tachyon potential for the study. In Sec. III C we compute the scalar and vector potentials for confinement, vector and scalar glueball masses and explore the physics involved. In Sec. III D, we compare our parameters with those from well-known phenomenological models. In Sec. III E we discuss the dual Higgs mechanism associated with the tachyon (and QCD magnetic monopole) condensation and confinement. In Sec. IV we present the results and analysis. Our final comments are contained in Sec. V.

II. MAXWELL'S EQUATIONS MODIFIED BY DIELECTRIC FUNCTION
In this section we will review electromagnetic theory in color dielectric medium. Beginning with Maxwell's Lagrangian with source its equations of motion are The equations of motion of the electromagnetic field create spherical symmetric solutions similar to the well-known point charge scenario associated with Coulomb's field [37]. Therefore, for electromagnetic field immersed in a dielectric medium G(φ), where φ(r) is the scalar field representing the dynamics of the medium, we can construct a Lagrangian density with its equations of motion given as We choose the indices to run as µ = 1, 2, 3 and ν = 0. These choices are done carefully to obtain chromoelectric flux confinement and to eliminate chromomagnetic fields in the rest frame of the particles. As a result, the equation of motion in Eq.(4) can be reduced to By this equation, we have eliminated the effect of magnetic field contributions leaving only the electric field contributions that will run through this paper. From Eq.(5), E is the electric field coupled to the dielectric function G(φ). Rewriting this equation in spherical coordinates, following the assumption that E(r) and φ(r) are only functions of r, we get Integrating this differential equation yields where q = 4 3 πR 3 ρ and E =| E |= E r . We observe that the dielectric function modifies the magnitude of the electric field function.

III. LAGRANGIAN DENSITY OF THE MODEL
In this paper we will focus on a model described by a Lagrangian density given by where m qq is the 'bare' quark mass or current quark mass when the particles under consideration are light, and sometimes referred to as the running quark mass if the particles involved are heavy. It should be seen to be the same as the mass term that appear in the QCD Lagrangian [39]. While M (φ) = m qq G(φ) is the constituent quark mass function which will latter be identified with the scalar potential S(r) in subsequent sections. The 'bare' mass (m qq ) leads to explicit chiral symmetry breaking while the M (r) permits both explicit and dynamical chiral symmetry breaking similar to renormalized mass in QCD Lagrangian [40,41], we will discuss it in details subsequently. It should be noted that due to the mass term together with gauge term, F µν F µν , the Lagrangian is certain to produce both scalar and vector potential contributions to the fermions. The equations of motion of this Lagrangian density are Again the indices are ν = 0 and µ = j = 1, 2, 3 as defined in the previous section. These choices are made deliberately to avoid the creation of chromomagnetic field so we can focus on the chromoelectric field which creates the flux tube picture, relevant for our analysis. Therefore, Eq.(11) becomes where we have substituted F j0 = −E. Expressing the above equation in spherical coordinates and integrating the results for electric field solution, yields the same result as Eq. (7).
Expanding Eq. (9) in radial coordinates to ease our analysis, we get Here we replace F µν F µν = −2E 2 and zero otherwise. This follows from the choice of indices defined above. For simplicity, we will also substitute (ε 0 = 1) thus, which implies It has already been established that G(φ) = V (φ) -see [43,44] and references therein -for slowly varying tachyons. This result will be used throughout this paper. In the above equation, if we consider a relatively large distance of particle separation from the charge q source, we can ignore the term with λ 2 , hence the equation reduces to We have used the general definitionψψ ρ(r) = qδ( r) in the above equation. Switching on the perturbation around the vacuum, φ 0 = 1/α, i.e. φ(r) → φ 0 + η(r), here η(r) is a small fluctuation about the true vacuum of the potential, then, Eq.(17) becomes where in the last step we have used Eq. (17) and anticipated the property of the scalar potential that we shall define shortly.

A. Determining QCD-like Vacuum and Gluon Condensation
This section will be a continuation of the review under 'Gluodynamics and QCD-like Vacuum' contained in Ref. [44]. We know that the Lagrangian for gluodynamics is symmetric under conformal transformation, when treated classically i.e. | v | → 0. Without any quantum corrections, its energy-momentum tensor trace is zero, θ µ µ = 0. As a result, it produces vanishing gluon condensate F µν F µν = 0, in the classical limit and non-vanishing gluon condensate F µν F µν = 0, with quantum corrections. Thus, quantum effects distorts the scale invariance [45] and brings about QCD energy-momentum tensor (θ µν ) trace anomaly a phenomenon well known in QCD theory. Here, β(g) is the QCD beta-function of the strong coupling g, with a leading term This model produces vacuum expectation value We will now compute the trace of the energy-momentum tensor from the Lagrangian Eq. (8) and compare the result with the result obtained above. This comparison is possible because, the third and the forth terms in Eq.(8) clearly breaks the scale invariant making it possible for comparison with Eq. (21). We calculate the energy-momentum tensor trace (θ µ µ ) by substituting Eq.(9) into the expression and this yields, Here, we have redefined Comparing Eq.(23) to Eq.(21) we get, We rescaleṼ ef f (φ) to include the vacuum energy density −| v |, i.e., consequently, This equation follows the classical limit, where the gluon condensate vanishes when | v | → 0 [46]. Again, we will demonstrate in the subsequent sections that, the gluon condensate increases with mass and remains non vanishing at m qq = 0, i.e. when the quark mass is 'removed' after confinement. This is a consequence of chromoelectric flux tube confinement.

B. Choosing The Appropriate Tachyon Potential
We select a suitable tachyon potential which produces tachyon condensation at low energies. This potential follows the restriction These restrictions are necessary to stabilize an asymptotically free system as well as its vacuum [47]. To proceed with the computations, we will require a suitable definition of the three dimensional Dirac delta function (δ( r)) that appears in Eq. (18). Consequently, we define δ( r) in the limit of step function as [48] δ( r) = where R is the radius of the hadron and r is the inter-particle separation distance. Solving Eq. (18) in the region r ≤ R, i.e. we are considering the presence of the particles inside the hadron, we get where K = 2 3qm qq 4πR 3 − 1 α 2 . This equation has solutions given as where Thus, Eq.(31) has two different solutions each corresponding to different physical regimes. The hyperbolic cosine represents a solution in the low energy regime (IR regime) whilst the sine function presents solution in the high energy regime where the particles are asymptotically free (UV regime) and some contribution from the IR regime where the particles are in a confined state. We will consider these solutions separately as we proceed.
Outside the hadron, r > R, δ( r) = 0, we have where K 0 = −2α 2 , which has a solution A quark which is kicked out of the pointlike region this way, at relatively high energy, behaves like a massless particle but remains confined to the hadron as would be shown latter. For now, we will focus our attention on the solution obtained in the region r ≤ R (inside the hadron) represented by Eq. (32). The color dielectric function for this solution takes the form Tachyon fields are generally unstable with negative mass squared, therefore it is a common phenomenon to find their vacuum states being unstable. In this case, we choose a potential whose true vacuum is at V (φ)| φ0 = 0. Switching on perturbation expansion about φ 0 will results in a generation of square mass proportional to V (φ)| φ0 . This stabilizes the tachyon fields and reduce their velocities significantly, making them viable for analysis in the infrared (IR) regime of the model. Simply put, tachyon fields are naturally associated with instability in quantum field theory hence, we perturb the fields about its true vacuum where the potential has its minima. As a result, our perturbation should be understood as a mechanism to stabilize the tachyon fields [49]. Following from the above equation, the color dielectric function in the IR regime becomes Substituting the above result for G(η) into the electric field equation in Eq. (7) we get The same tachyon potential and the procedure adopted here was used in [44] to confine light quarks at a finite temperature. Consequently, the tachyons are expected to generate mass at V (φ)| φ0 leading to a particle-like state resulting in tachyon condensation [50]. It follows that, the vacuum stability is independent of the constituent quark mass appearing in the Lagrangian density. These properties observed from this type of potential makes it more efficient for computing confinement potentials and glueball masses for heavy or light (anti-)quark systems. From the afore analysis, we can associate f α = 1/α to the decay constant of the tachyons, the higher its value the faster the tachyons decay and by extension, the more unstable the QCD-like vacuum created in this process and the vice versa.

C. Potentials, Glueball Masses and Constituent Quark Masses
1. The potential of the particles inside the hadron Using the well-known relation for calculating electromagnetic potentials to determine the confinement potential V c (r, m qq ). Now we apply the hyperbolic cosine function at the left side of Eq. (32) representing the solution in the IR regime to get here, β is a constant representing the depth of the delta function well, i.e. β = δ(R → 0) = 3/(4πR 3 ). The above equation gives the net potential observed by massive (anti-)particle pairs ((anti-)quark pairs). In effect, particles of the same kind (two particles or two anti particles) repel each other whilst particles of different kinds (particle and anti-particle) attract; for example, color (blue) attracts an anti-color (anti-blue) of the same kind or colors of diffident kinds (blue and green) attracts each other and the vice versa [11]. We choose the negative part of the potential which correspond to the potential of an anti-particle. It is appropriate to choose q = −1 also corresponding to an anti-charge. We will latter find that q = −1 is associated with an anticharge which is identifiable with an anti-color charge carried by the gluons. Considering that, we are dealing with the potential of an anti-particle, this choice is justifiable. Also, we will choose R = (3/4π) 1/3 i.e. β = 1 and c = 0 (c is the integration constant). Consequently, the net static potential (a combined scalar and vector potentials) seen by an antiparticle confinement in the region r ≤ R, in the rest frame of a heavy source is given by and its string tension is given as 1 The scalar potential contribution here is due to the heavy point-like source situated at the center of the hadron. In this regime or at relatively large inter particle separations r, particles are always confined with hadron degrees of freedom. This result does not present UV effects, i.e., Coulombic potential effect at r → 0. If we choose σ c ∼ 1GeV/fm and α = 0.99 such that f α 1 we obtain m qq = 10.10MeV, which lies within the range of u and d-quark masses. It is important to add that this potential is good for investigating light quarks such as u, d and s-quarks, since the potential already lies in the stable regime of the theory where all the quarks are confined into hadrons. Heavy quarks such as c-quarks starts experiencing the degeneracy (pair production) as shown in Fig. 1 with dashed lines. Therefore, this model is recommended for confining light quarks such as u, d or s only.
This result accounts for both spin-dependent and spin-independent relativistic corrections as QCD predicts for slowly moving quark systems. Because there are spinless quarks (comparatively light) in a slow motion around the [1] The parameters combined as follows are in general related according to the form mqq + 1/(2πα ) 1/2 α 2 , whereα 2 = α 2 /(2πα ) 1/2 . This is because the precise relationship between tachyon potential and color dielectric function is G(φ) = (2πα ) 2 V (φ), where α is the Regge slope which has dimension of length squared. However, we have assumed 2πα = 1 along the paper. heavy point-like anti-quark source at the origin. The co-motion of the light quarks relative to the central massive static anti-quark source creates spin-dependent corrections at short ranges (vector-type interaction) and spin-independent correction at long distances (scalar-type interaction) [51] giving rise to the net potential. To this end, when the surrounding co-moving quarks are far away from the heavy source, the string tension that binds them to the massive anti-quark source will break leading to pair production (a (anti-)quark pair creation), a phenomenon well known in QCD theory. However, in the presence of a heavy anti-quark source, the new quark created in the process remain attractive towards the original source whereas the new anti-quark created in the process remain attractive towards the 'break-away' quark. In this regard, the (anti-)quark pairs will always remain confined once there exist, at least, a heavy source. The process remains the same if we consider a heavy quark source as well. This characteristic is depicted by the linear nature of the string tension changing with m qq and constant when m qq = 0.
We will now analyse the behaviour of the particles in the UV regime where the particles are relatively close i.e. small inter particle separation r, with quark and gluon degrees of freedom. We will use the solution of the sine function at the right side of Eq.(32) for this analysis, Substituting this result into the color dielectric function Eq.(36) results in, Substituting the above equation into Eq. (38), with string tension Where we have chosen the positive part of the potential corresponding to λ = q = 1 andc = 0. It is known in perturbation QCD that the dominant interaction at small distances, r → 0, is Coulombic one-gluon exchange (OGE). We can estimate the running quark mass, (m qq ), using this result, if we consider a typical hadron of mass 1GeV with radius 1fm as determined in electron scattering with estimated string tension σ s ∼ 1GeV/fm and a decay of f α = 1. We can estimate the running quark mass as m qq = 4GeV, a typical mass for heavy quarks. The mass is slightly less than the mass of a b-quark, m b = 4.18 ± 0.03GeV on mass-independent subtraction scheme (M S) at a scale of µ = 2GeV as reported in the Review of Particle Physics of the Particle Data Group [39]. In addition, light quark masses (such as u, d and s-quark masses) are difficult to determine due to their small masses relative to hadron scale, so it is sometimes difficult to classify and significantly identify their influence within hadrons. Thus, this model is vaible for investigating confinement of quarks with masses m qq ≥ 4GeV, bellow this threshold the quarks are expected to be asyptotically free in the UV regime. It is therefore, convenient to study the behaviour of light quarks using the IR potential developed in Eq. (40), instead of the Cornell-like potential in Eq.(45) known for confining heavy quarks [52].

String Tension
Generally, the Cornell potential for confining heavy quarks is giving by where e and σ are the Cornell fit parameters, but σ is related to the lattice spacing a as σ ∼ 1/a 2 and V 0 is the self energy of the static source. The separation distance, r, between the particles determine the magnitude of a, i.e. large distances imply wide a and vice versa [52]. Determination of e is strictly phenomenological, its larger value correspond to smaller mass and the vice versa. In relation to our model, we normalized it at V 0 =c = 0, this is to remove the divergence known to be caused by self-energy contribution in the continuum limit and e = 0.25 [53]. Unlike the quenched approximation where the particle and antiparticle numbers are independently conserved, in this model the conserved quantity is the baryon number due to the presence of the current quark mass and a constituent quark mass function. As would be shown explicitly later, light quark-antiquark pairs are created in the vacuum resulting from the transition in the string tension connecting the two static sources leading to hadronization. When the energy carried by the string tension connecting the two static sources exceed its critical value at some separation distance r = r * or some critical mass m qq = m c , the string will 'break' and decay into light static mesons. Thus, the Coulomb effect e is rather pronounced for higher masses and weaker if the mass involved is slightly weaker than the infrared mass -see the results published in Ref. [55]. Thus in the limit r → r * the potentials stop rising with distance r and the static source quark becomes screened by the light quarks formed in the vacuum. Similar behaviour is observed when we keep the distance fixed and increase the mass i.e. m qq → m c . At the ground state σ c > σ s . To obtain confinement, the σ c can contain current mass within the range 0 ≤ m qq ≤ 10.10MeV whereas σ s can contain current mass within the range 4 ≤ m qq ≤ m c GeV due to the UV contributions. The running masses for the IR and the UV regimes can be compared with the results of the running masses calculated using various QCD sum rules within the size of a hadron. The masses are giving asm u (1GeV) = 5.2 ± 0.5MeV,m d (1GeV) = 9.2 ± 0.5MeV,m s (1GeV) = 159.5 ± 8.8MeV and m b (1GeV) = 5.8 ± 0.06GeV [54,56].

Vector Potential
Using the solution at the left side of Eq.(34) i.e. outside the regime of the pointlike source, we will have two separate potential contributions; vector potential due to the gluonic sector and the massless quark and a scalar potential energy contribution from the system of the pointlike particles (hadron). Here, the hadron serves as a massive point-like source that confines the massless quark. We will compute the vector potential in this section while we reserve the calculations of the scalar potential energy and the net potential energy for the next section. Using Eqs. (37) and (38), the vector potential becomes, with string tension This corresponds to the short range Thomas spin-orbit interactions as highlighted in the introduction. Again, faster tachyon decay means weak string tension while slow decay means strong string tension and stronger confinement. Similar results is obtained if we set m qq = 0 in Eq.(39).

Scalar Potential
To determine the resulting scalar potential due to the hadron, we compare our results from Eq.
Here, we have used P µ = ( E c , P ), A µ = (Φ(r), A), assuming that the quarks are scalar (spinless) and static, we have E = A = Φ(r) = 0. In spherical wall potential, we can impose the restriction V (r) = S(r), for r ≤ R, 0, for r > R where r is the inter-quark separations and R is the radius of the hadron. Here, V (r) = S(r) represents the scalar potential in the Dirac equation. Thus, Eq.( 50) can be rewritten as (cαp +βm 0 c 2 )ψ + S(r)ψ = 0, whereα andβ are Dirac matrices. Also,p → −i∇ ( = 1) is the momentum operator, m 0 is the rest mass of the fermions and c is the speed of light in vacuum. The vector potential is normally introduced by minimal substitution of the momentum P µ → P µ − gA µ whilst the scalar potential is introduced by the mass term m → m 0 + S, g is a real coupling constant. It should be noted that the vector and scalar potentials are coupled differently in the Dirac equation [57]. Comparing Eq.(10) and Eq.(52), we can identify the first and second terms of both equations as the interaction terms and the third terms as the scalar potentials [58,59]. We leave out the vector potential because it has already been calculated in Sec. III C 3 so, it is of no further interest. The resulting scalar potential seen by the fermions is This result represents Thomas spin-orbit interactions at long ranges, where the short range Thomas spin-orbit interactions are partially dominated giving way to scalar interactions only. As a result, the interactions are thought of as being concentrated on the various quark coordinates. Interestingly, the scalar potential energy is simply a product of m qq and the color dielectric function. This gives an indication that confinement in this scenario has a direct relation with tachyon condensation and m qq of the system. Meanwhile, some authors have predicted vector dominance over scalar for inter quark potentials, suggesting that scalar dominance will imply that all quark combinations are confined. Contrary to that, only the energetic combinations are preferred phenomenologically [11]. We will attempt to analyse this assertion by computing and comparing the magnitudes of their coupling constants as they are expected to appear in the Dirac equation and its significance. The net confining potential energy of a quark outside the hadron will be the sum of Eqs. (48) and (53) yielding Here, we have multiply the vector potential in Eq.(48) with an anti-particle charge −q to obtain a vector potential energy that goes into the net potential energy.
To end this section, we will present the color dielectric function and the scalar potential for the sine function used in investigating the particles inside the hadron in the UV regime. Consequently, the color dielectric function will be and the scalar potential of the hadron reads for λ = q = 1.

Glueball Masses
This model is certain to produce both vector and scalar glueballs just as vector and scalar potentials calculated above. Scalar glueball mass has been estimated to have a value within the range 1.5 to 1.7GeV. This value has been affirmed by data and calculations [60]. Specific value was reached in Ref. [61] to be 1.7GeV for fitness, through unquenched calculation. On the other hand, the existence of vector glueball mass has been predicted with an extimated value of 3.8GeV by quenched Lattice QCD [62]. The vector glueball masses are expected to be observed at the ongoing Beijing Spectrometer Experiment (BESIII) and hopefully future PANDA experiment at the FAIR Lab [63]. Recent findings published by BESIII facility [64] points to vector particles of mass 3.77 to 4.60GeV with precision.
From the tachyon potential in Eq.(28) we can directly calculate the 'vector glueball' mass (m gb ) as The scalar glueball mass is then calculated directly from the Lagrangian in Eq.(8) as where, we have substituted ψψ = qδ( r) → qβ, q = −1 and β = 1. Following the analysis in Sec. III A, we find that F µν F µν | φ0 = 0. Therefore, the model produces vanishing gluon condensate in its vacuum (φ 0 ). We can now rewrite Eq. (58) in terms of the string tension σ c and the vector glueball mass as, From the afore analyses, we find that the vector and scalar glueballs depend on the tachyon decay constant, α = 1/f α . Whilst the vector glueball depends directly on f α the scalar glueball only depends on f α through the string tension σ c . Knowing the estimated values of scalar and the vector glueball masses and the value of the string tension σ c ∼ 1GeV/fm we can fix the tachyon decay constant within the range, 0 < f α ≤ 1. If we consider that the tachyons are decaying at half of its maximum value, i.e. f α = 0.5, we will obtain m gbφ = 2GeV and m gb = 4GeV. This decay regime gives a good agreement with the vector glueball mass within the range of the experimental values. That notwithstanding, vector and scalar glueballs are more likely to be observed for slowly decaying tachyons.

Constituent quark masses
We would like to make brief comments on the constituent quark mass function M (φ) that was mentioned bellow Eq. (8). This function is the same as the scalar potential energy S(r) determined from the Dirac's equation and discussed under Sec. III C 4. Here, we will seek to determine its value for r → 0 and r → r * i.e. the excited state and the ground state respectively. Now expanding the solution for the hyperbolic cosine in Eq.(32) for 2|K|r 1, we get Using the expression for G(η) in Eq.(36) and the above solution, we can write Thus, the constituent quark mass for the ground state energy will be M (r → r * ) = 4m qq [53], where r * represents the distance within which the quarks are confined and beyond it we have degeneracy. Therefore, the constituent quark mass (M r * ) for the potential in Eq.(40) is 40.4MeV, and degeneracy should be expected beyond this mass limit.
We can now proceed to equivalently calculate the constituent quark mass (M 0 ) for the highest excited state r → 0 by using the sine function solution in Eq. (32), this yields Therefore, M (r → 0) = 4m qq is the constituent quark mass in this regime, if we choose a running quark mass of m qq = 4GeV as calculated above, we get, M 0 = 16GeV. So, the constituent mass (M ) in this model framework lies between 0 ≤ M ≤ 16GeV, where 0 ≤ M r * ≤ 40.4MeV gives the dynamics of the IR regime while the UV regime can be studied within 4 ≤ M 0 ≤ 16GeV. The quark masses that appear in phenomenological models are generally the constituent quark mass [39] which dynamically breaks the chiral symmetry. In non perturbative theories the range for dynamical chiral symmetry breaking, Λ χ is about 1GeV [42]. It is therefore conventional to say a quark is heavy if M > Λ χ , in this case, we have explicit chiral symmetry breaking for c, b and t-quarks whilst quarks are classified as light if M < Λ χ , leading to spontaneous chiral symmetry breaking dominance. Quarks in this category include u, d and s. In nonrelativistic quark models, one of the useful parameters is the constituent quark mass which is determined to be M u = M d = 350MeV for light quarks in single gluon exchange (OGE) interaction, while the predicted mass for a c-quark is also M c = 1.6GeV. Constituent quark mass models are used to study the effect of dynamical chiral symmetry breaking and they are independent of the current quark mass. The constituent quark mass changes depending on how measurements are made in a particular model framework. The constituent quark masses are greater compared to current quark mass and are usually free parameters to fit in potential models. The masses vary depending on whether you are using a meson fits, baryon fits, hadron fits or other phenomenological models.

D. Identification of Coupling Constants
In this section, we will compare the result from equation (16) with the well-known phenomenological models [11,65], using dilaton in gauge theory to confine quarks and gluons with N c colors. For emphasis, since our model is phenomenological, it lacks all the degrees of freedom required to fully represent pure QCD theory, but our results agree with the QCD spectrum for scalar and vector potentials.
The results presented above still keeps the electromagnetic charge q and do not contain the color numbers N c . This comparison is intended to help fill in these gaps. Furthermore, the scalar and the vector potentials in Dirac equation are of the same weight but coupled differently. This comparison will enable us to determine their couplings and relate them to their individual strengths in the QCD-like model. The same method was used in [43] -and references herein -to determine couplings as well. The dilaton model is given as Transforming the exponential potentials of the above equation to conform with the tachyon potential used, i.e, , hence equation (63) can be rewritten as It should be noted that the purpose of this section is to determine the coupling constants, so the equations have been simplified to achieve that goal. Comparing equations (16) and (64), we can find ξ = 2, and α 2 = 1/f φ = 1/f 2 α . Here ξ represents the coupling strength of the fermions, f φ is the mass decay constant of the dilaton and it is related to the decay constant of the tachyons f α , g is the gluon charge and it is related to the electromagnetic charge q. We can rewrite equations (14), (40) and (41) as and for g → g/(2 √ 2) respectively. The scalar potential energy observed by the heavy pointlike anti-quark source is also given as S(r, m qq ) = 2m qq α 2 η(r) 2 = m qq 1 r 2 (1 + m qq )α cosh 2 (2 (1 + m qq )αr) .
We find from the above equation (36) for the dielectric function G that, the larger the quark mass the more stable the QCD-like vacuum and the faster the tachyons condense [66]. Since tachyon condensation implies confinement, then heavier quarks are more likely to be confined than light quarks [67]. This potential follow the same analysis as the dielectric function with the tachyons condensing faster in a multiple of m qq showing a rather stronger confinement but vanishes at m qq = 0. This result represents the energy flow between a quark and an anti-quark pairs at long ranges.
On the other hand, Eq.(66) represents the net potential for confinement observed for a heavy anti-quark source in the origin surrounded by relatively light spinless quarks in a slow motion.
It is important to state that, the net confinement potential seen in Eq.(66) corresponds to a mixture of ξ = 0 (no mass coupling) and ξ = 1 (Kaluza-Klein type or mass coupling) [68]. The combined ξ = 1 and ξ = 0 gives the net confinement of the quarks at all masses inside the hadron. Also, Eq.(68) gives the scalar potential energy, i.e., the flow of energy between two or more quarks. The coupling of the scalar potential, ξ = 2, exposes its weakness relative to the vector potential in Eq.(48) by a ratio of 0 : 2. By simple interpretation, the scalar potential energy must be coupled strongly in order to coexist with the vector potential energy which needs no coupling as seen in our model framework and presented in expression Eq. (54). It is noticeable that the scalar potential energy vanishes at m qq = 0 while the vector potential energy does not as evident in Eq. (54). Consequently, quarks are always in a confined state [69] under the net confining potential Eq.(66) and at m qq = 0 we retrieve the vector potential Eq. (48). This indicates that the chromoelectric flux generated by the colored particles remain confined even if we 'remove' the quarks after confinement under the net confining potential (V c ) [2]. It should be noted also that, the net confinement potential of the color particles is the sum of the vector and the scalar potentials as presented in (66) and the net potential energy is also the sum of the vector and scalar potential energies as expressed in Eq. (54).
After knowing all the coupling constants and the nature of the Dirac delta function, few comments on Sec. III A will be necessary. Noting that the Dirac delta function is well defined inside the hadron r R, with depth β = 3/(4πR 3 ) = 1, G(φ) = V (φ) and an anti-gluon color charge q → g = −1 (N c 1), Eq.(24) takes the form Substituting Eq. (35), recalling that V (η(r)) = G(η(r)) = 2α 2 η 2 for φ(r) → φ 0 + η(r), into the above expression yields Since tachyon condensation is faster whenṼ ef f → 0, see Fig. 3, tachyons condense with increasing mass and attain its minimum condensation at m qq = 0 corresponding to the maximum value ofṼ . As a results, Eq.(27) yields, By this result, tachyon condensation increases with increasing mass m qq , and attain its minimum when m qq = 0. There is a QCD monopole condensation associated with both instances (i.e. m qq = 0 and m qq > 0 ), since the gluon condensate do not vanish in any of the two instances. The condensate is higher with increasing mass and relatively low at no mass. Naturally, confinement results in both instances.

E. Tachyon Condensation and Dual Abelian Higgs Mechanism
For the sake of simplicity, but without lost of generality, in the previous sections we have postponed the discussion concerning the dual description of the confinement in terms of the Higgs mechanism. We shall now complete the discussion by extending the the Lagrangian in Eq.(8) by imposing gauge invariance on the scalar sector. As such, the degrees of freedom should be carefully changed. Firstly we have to consider a charged scalar field, i.e., a complex scalar field We will leave out the fermion (spinors) coupling from the original Lagrangian in Eq. (8) for this analyses because it has no direct influence on the outcome of the intended result. Consequently, where F µν = ∂ µ A ν −∂ ν A µ andF µν = ∂ µÃν −∂ νÃµ are two independent Abelian field strengths and D µ = ∂ µ −iqÃ µ is the Abelian covariant derivative associated with the dual gauge fieldÃ µ which is responsible for the magnetic monopole description in the dual Higgs mechanism. Besides its original U (1) gauge invariance, the Lagrangian becomes invariant under the followingŨ (1) gauge transformation We can now analyze the dual Higgs mechanism of the model. Going forward, the potential of the model in Eq.(28) will also take the form The minima of this potential function is given as The non-zero vacuum expectation value breaks theŨ (1) gauge symmetry spontaneously.
Let us now proceed as in our previous analyses by expanding about the vacuum of the potential so, for ease of calculations, we choose physical vacuum configuration φ 1min = φ 0 and φ 2min = 0.
Now, we choose two real scalar fields η(r) and ζ(r) [70,71] to represent small fluctuations about the vacuum of the potential. So we can represent the shifted vacuum as As such, we can conveniently parametrize the φ field as The Lagrangian for the small fluctuations about the vacuum up to quadratic order becomes Notice that the scalar field η plays the role of the usual massive Abelian Higgs field and the scalar field ζ is related to the massless Goldstone boson. Now choosing the gaugẽ we can write Thus, the Lagrangian takes the form This is precisely the Lagrangian (8) evaluated around the vacuum plus the Lagrangian of the dual gauge field. Notice that G(η) = V (η) = 2α 2 η 2 is consistent with (35). Furthermore, as we have well discussed in the previous sections, in the limit η → 0 we have tachyon condensation and confinement. Then at such a limit, we are left with the Lagrangiañ which describes the dynamics of a massive dual gauge field -a consequence of the dual Higgs mechanism.
The equations of motion of the dual gauge field are and that for the static fields give where the magnetic monopole charge and current densities ρ m , j m , respectively, are defined by the current density j µ = q 2 φ 2 0Ã µ . So these equations govern a dual superconductor with the dual London equation ∇ × j m = (1/λ 2 ) Ẽ , where λ = (q 2 φ 2 0 ) −1/2 is the penetration depth [72]. This scenario leads to the description of condensed magnetic monopoles and confined electric field, as we previously assumed.

IV. RESULTS AND ANALYSIS
Plotting the results of equations (66), (67), (36), (68), (45), (46), (56) and (55) in Fig. 1, Fig. 2, Fig. 3, Fig. 4, Fig. 5, Fig. 6, Fig. 7 and Fig. 8, we assume that f α = β = 1, g/4π = λ = ±1 and N c 1. Fig. 1 shows graphically the relationship between the inter-quark potential V c (r, m qq ) with (r, m qq ) for a heavy anti-quark source. The graph shows a steady increase in the gradient from m qq1 to m qq3 representing an increase in the strength of confinement from m qq1 to m qq3 . Fig. 2 represents the graph of σ c (m qq ) against m qq , it shows a linear increase in σ c (m qq ) against m qq with its foot intersecting the σ c (m qq ) axis at 1 indicating the strength of interaction even at m qq = 0. The linearity in σ c (m qq ) depict the confinement that exist between the (anti-) quark pairs. Fig. 3 shows the graph of the QCD-like vacuum which is equivalent to the tachyon potential. Tachyon condensation is related to monopole condensation which translates into confinement. The deeper the depth of the curve the more condensed the tachyons and the stronger the strength of the confinement and the vice versa. Again, Fig. 4 shows the relationship between the scalar potential S(r, m qq ) with (r, m qq ). We notice that the greater the quark mass the deeper the depth of the curve and the smaller the minima of the curves representing faster gluon condensation. The above discussions so far are related to the IR regime. But, same is true for Fig. 5, Fig. 6, Fig. 7 and Fig. 8 which incorporate some UV characteristics, the only difference is that heavy quarks (m qq ≥ 4GeV) are required to achieve confinement in this regime due to the UV characteristics. Considering all the factors discussed above, we can deduce that the vector potentials are dominant over the scalar potentials for light quarks in the IR regime and heavy quarks UV regime respectively.
By way of analysis, we will compare the vector and the scalar confinement as captured in literature. The vector potential takes into account the angular momentum of the quarks whilst the scalar potentials do not. The scalar potential actually thrives on the basis that the angular momentum is partially eliminated. In the vector confinement, the energy of the quarks are carried by the gluons (string tension) whereas in the scalar confinement, the energy is centered on the quark coordinates. The scalar potential is most useful in a non rotating quark frames whilst the vector potential is a viable option for a rotating quark frames. Obviously, the 'Thomas precession' that gives rise to vector confinement differs significantly from the one that gives rise to the scalar confinement [3,4,22,23].  The graph in the left panel is for different values of mqq. It shows the strength of qq confinement for different values of mqq as their separation distance r vary. The graph shows a steady increase in its gradient as mqq is increased from 0 to 0.7GeV. By inference, the particles are strongly confined with increasing mass. That notwithstanding, a system of fermions remain confined even if the mass of the fermions is 'removed' (mqq = 0) after confinement due to the chromoelectric flux tube confinement, in consistence with the QCD theory. The graph in the right panel is for infinite limit of quark masses, it rises sharply, indicating confinement and flatten up shortly indicating hadronization. Λχ = 1GeV is the threshold mass bellow which quarks are classified as light and beyond which they are classified as heavy. The linear nature of the graph indicates that the qq is confined for increasing mqq. The foot of the graph at 1 shows that the particles remain confined even at mqq = 0 in consistency with QCD spectrum for heavy quarks.   The potential vanishes at mqq = 0, for both graphs and increases in depth as mqq increases as shown in the left panel. The graph in the right panel has a deeper depth and represents stronger confinement. The S(r, mqq) deceases with mass from the left until it attain its minimum at S(r, mqq) → 0 and start rising from the minimum towards the right, this is true for both graphs. The decrease implies stable confinement while the rise towards the right signifies degeneracy and screening of the static antiquark source. These are characteristics demonstrated by light mesons [53].  The graph in the left panel rises steadily with increasing mass indicating strong confinement with increasing mass whilst the graph in the right panel shows stronger confinement with infinite increase in quark mass.
FIG. 6: A graph of string tension σ s (m qq ) resulting from the Cornell-like potential against mass, m qq , for a heavy quark source.  The potential vanishes at mqq = 0. The scalar potential increases with increasing mass. Thus, the higher the quark mass the more likely the quarks will be confined even at higher energy regimes (small r) for finite quark masses as shown in the left panel. For infinite quark masses (right) the scalar potential also increases as mqq → ∞, consequently confining the quarks at high energy regimes.  The graph in the left panel shows slight increase in tachyon condensation with increasing quark mass, so increase in the quark mass goes into confining the quarks at higher energies r → 0. Likewise, in the right panel, the tachyon condensation increases steadily as mqq → ∞ resulting into confinement at higher energy regimes.

V. CONCLUSIONS
The model gives an insightful details about QCD theory in both the IR and the UV regimes making it more efficient for consideration. We have calculated the net potential of the slowly moving quarks inside the hadron, vector potential, scalar potential energy, net potential energy and the string tensions associated with them in both regimes. Vector and scalar glueball masses which are consequences of the IR regime were also studied and their values for f α = 0.5 (consequently, φ 0 = 500MeV) and σ c ∼ 1GeV/fm were calculated and compared with the existing QCD lattice results for quenched (no quark fluctuation) and unquenched (hybrid fluctuation) approximations respectively.
The dominance of the vector potential expected in QCD theory was clearly demonstrated by computing the coupling strengths ξ for both the vector and the scalar potentials. We also find that the scalar potential (energy) must be strongly coupled in order to coexist with the vector potential (energy) which needs no coupling at all in this model framework. The vector and the scalar potentials should be understood as resulting from relativistic spin-orbit corrections at short and long ranges respectively. The vector potential (energy) shows confinement of chromoelectric flux at zero mass but the scalar potential (energy) vanishes at the same mass. Thus, the vector potential (energy) is more suitable for confining light quarks whilst at least a heavy source is required to obtain scalar potential (energy). Furthermore, we established the relation between tachyon condensation, dual Higgs mechanism, QCD monopole condensation [73] and confinement. Finally, we intend to continue our series on this subject by studying confinement of fermions at a finite temperature.