Phenomenology of strong interactions: towards an effective theory for low energy QCD

In this paper, we develop models applicable to phenomenological particle physics by using the string analogy of particles. These theories can be used to investigate the phenomenology of confinement, deconfinement, chiral condensate, QGP phase transitions, and even the evolution of the early universe. Other confining properties such as scalar glueball mass, gluon mass, glueball-meson mixing states, QCD vacuum, and color superconductivity can also be investigated in these model frameworks. We use one of the models to describe the phenomenon of color confinement among glueballs at the end of the paper. The models are built based on the Dirac–Born–Infeld (DBI) action modified for open strings with their endpoints on a Dp-brane or brane–anti-brane at a tachyonic vacuum.


I. INTRODUCTION
String theory was conceived in the late 1960s to provide an explanation for the behavior of nuclear matter such as protons and neutrons [1,2].Even though the theory was not successful in explaining the characteristics of quarks and gluons at its inception, it promised an interesting intervention in other areas of physics [3][4][5][6][7].It showed capability in giving an insight into cosmology, astrophysics, and unification of the fundamental forces of nature which has been on the table of physicists for some time now.Consequently, the theory of Quantum Chromodynamics (QCD) was developed in the early 1970s to give a comprehensive explanation of nuclear matter [8,9].The QCD theory is now accepted as the standard theory for strong interactions.Interestingly, recent development shows that string theory and QCD describe the same physics.
The formulation of Quantum Electrodynamics (QED) and QCD are almost the same.They are both formulated from field theory based on gauge theory which forms the foundation of the highly acceptable Standard Model (SM).The fundamental particles of these theories are studied based on the gauge boson they interchange: Photons for electrodynamic force, W ± and Z 0 bosons for electroweak force and gluon for strong nuclear force [10][11][12].Subsequently, gravitational force does not fall under this category, so they are investigated under Einstein's general theory of relativity.Physicists have bunged their hope on string theory to unify all the fundamental forces, this will fall under 'physics beyond the SM' [13][14][15][16][17]. Upon the similarities, quarks and gluons are color particles classified under three conventional colors (red, blue, and green) whilst photons are color neural bosons that mediate electrically charged leptons.Also, gluons self-interact due to their color charges but photons do not [10,12].
Additionally, QCD falls short in explaining color-neutral particles such as bound states of gluons (glueballs) and quarks (hadrons and mesons), so string theory can be resorted for further description.The string-like description of hadrons arising from the quark model [18] is an important phenomenon in applying string theory.Under this picture, when a quark and an antiquark are pulled apart, they behave as if they are connected by a rubber band (gluon) which becomes increasingly difficult to separate when the separation distance between them keeps increasing.This analogy gradually fails when the particles are brought closer and closer together.In this regime string theory fails and QCD theory becomes viable.Now, looking at the particles in terms of fundamental strings, string theory describes hadrons quite well and provides the background for the unification of the fundamental forces.In string theory, the strings are treated as particles, where different particles are associated with different string oscillations.The masses of the particles are also associated with the energy of the oscillating string.The intrinsic spin of the particle is associated with the two perpendicular oscillations of the string with its endpoints fixed on D-branes, similar to the direction of electric and magnetic fields in a photon.Hence, photons and gluons can be identified in terms of the open strings as spin-1 particles.The clockwise and anti-clockwise movement of the closed strings makes it possible to classify it among spin-2 bosons such as graviton.Since, QCD involves color fields, to study it under string theory we assume that the endpoints of the open strings on a D-brane serve as the source and sink of the color charge [20].
It has been conjectured that open bosonic strings studied at a tachyonic vacuum behave as if they are closed strings with no D-branes.However, soliton solutions in this region point to the presence of lower dimensional branes [19,21].This projection is also corroborated in superstring theories [22,23] and evident in the first [24][25][26] and second [27][28][29][30][31][32][33][34] quantizations in string theories [35].At the tachyonic vacuum, the negative energy density V (T 0 ) of the tachyons in the vacuum exactly cancels [36] the energy density of the D-branes ε p i.e.V (T 0 ) + ε p = 0.In non-BPS D-branes ε p = τ p where τ p is the D-brane tension [37][38][39].It should be noted that at |T | = T 0 the total energy density of the tachyons and the brane tension vanishes identically, signifying the absence of open strings at that point.In this view, at a tachyonic vacuum, there will be no physical open string excitations because there are no D-branes.On the contrary, the 'usual' field theory gives an alternative explanation, because shifting the vacuum and expanding around its 'true' minimum can change a negative square mass of the tachyons to a positive one even though it does not completely remove all the states.Since it will not cost any energy to adjust the fundamental strings on the worldvolme of the D-branes at a tachyonic vacuum, it will be difficult to notice their presence.Hence, fluctuations around the vacuum represent the presence of lower dimensional D-branes [19,21].These phenomena have been investigated in references [27,28,30] using open string field theory [40].
In this paper, we modify the Dirac-Born-Infeld (DBI) action, so that the associated open strings with their endpoints on the remaining lower-dimensional branes at the tachyonic vacuum can be analyzed.The objective is to develop models that can mimic QCD theory both in the UV and IR regions with UV safety.The models can be applied in developing potential models such as the linear confining and Cornell potential models.In the analysis, we consider that the string worldsheet falls inside the D-brane worldvolume making it easy for the endpoints of the strings to be connected by the flux line on the worldvolume to form a close string suitable for modeling color confinement in a flux-tube like picture.The dynamics of the strings tangential to the D-brane worldvolume is represented by the gauge field F µν and the component transverse to the worldvolume is represented by a massless scalar field X a .The net flux involved is determined by the source and sink of the flux carried by the endpoints of the string on the lower dimensional remnants of the original D-brane in the vacuum.Also, the condition for minimum energy warrants that the flux does not spread because the source and sink of the flux emanate from a pointlike object on the D-brane worldvolume.
The paper is organized as follows: In Sec.II we review Tachyons and D-branes, divided into two subsections, under Sec.II A we review Tachyons, and in Sec.II B we review D-branes.We present the Modification of the Dirac-Born-Infeld Action which serves as the bases for the study in Sec.III.Also, in subsection III A we review the Dimensional Reduced U(1) Yang-Mills theory and its intended consequences in relation to the Standard Model (SM) of Particle Physics.Section IV and subsequent sections contain our original contributions to the subject.We present the Bosonic String at Tachyonic Vacuum in Sec.IV, where we present the details of dimensional reduction from 10D Dirac-Born-Infeld action to 4D conducive for describing SM particles.We present Gauge Theories Modified with G(φ), in Sec.V, this section was divided into two.We studied Fermions Coupled to Fundamental Strings at Tachyonic Vacuum in Sec.V A and Non-Abelian Gauge Theory in Sec.V B. Section VI contains the Phenomenon of Gluon Confinement, divided into four subsections.We present The Model in Sec.VI A, Confining Potentials in Sec.VI B, Gluon Condensation in Sec.VI C and strong Running Coupling and QCD β-function in Sec.VI D. we present our findings and Conclusions in Sec.VII.

A. Tachyons
Generally, tachyons are classified as particles that travel faster than light or weakly interacting superluminal particles.Relativistically, single particle energy is expressed as E 2 = p 2 c 2 + m 2 c 4 , where p is the spatial momentum and m is the mass of the particle.For a particle to be faster than light, the relativistic velocity β = pc/E > 1.Hence, for tachyons with real p, m must necessarily be imaginary [41].However, this analogy does not make a strong and convincing case for the tachyons.Rather, Quantum Field Theory (QFT) provides some satisfactory explanation to the dynamics of tachyons.QFT suggests that particles that travel faster than light does not exist in nature.So tachyons are simply unstable particles that decay.Based on this understanding, we consider a scalar field, say φ, with the usual kinetic term and a potential V (φ) whose extremum is at the origin.If one carries out perturbation quantization about φ = 0 up to the quadratic term and ignores the higher order terms in the action, we obtain a particle-like state at V (0).These results have two interesting interpretations; if V (0) > 0 we have a particle of mass m 2 φ but for V (0) < 0 we have a tachyonic state with a negative m 2 φ .In this case, tachyons can be given a physical meaning.So far, we know that the tachyons have negative m 2 φ and potential whose maximum is at the origin, thus a small displacement of φ at the origin will cause it to grow exponentially with time towards the true minimum.By the above description, tachyons can be represented by a potential such as where c 3 and c 4 are constants.Hence, tachyonic fields are associated with the 'usual' Higgs fields where the Higgs particles acquire negative square mass at the true minimum of their potential.Accordingly, tachyons in QFT are related to instability that breaks down the perturbation theory in normal field theory.The usual quantization where the cubic and higher order terms are considered small corrections to the quadratic term is no longer tenable.Since V (φ) has its maximum at φ = 0, it renders it classically unstable at that point.So one cannot guarantee that the fluctuation at that point is small.This behavior comes with an inbuilt solution i.e. one can expand the potential about the true minima φ 0 up to the quadratic term and proceed with the perturbation quantization about that point.Thus, the cubic and higher-order terms in the expansion can be discarded.This process will lead to the creation of a particle with a positive mass m 2 φ in the spectrum.This process removes the tachyonic modes in the spectrum.Additionally, in Dp-brane systems, the theory must be invariant under Z 2 symmetry, i.e. φ → −φ and in the presence of brane-anti-brane, the theory must necessarily be invariant under phase symmetry φ → e iα φ.Indeed, there are some benefits in working in tachyonic modes otherwise we can define a new field, η = φ − φ 0 , and express the potential in terms of the new field, Working with this potential from the onset will remove the tachyonic mode from the spectrum because V (η = 0) will be positive.However, there are some benefits to working with tachyonic fields, for instance, they possess high symmetry.This symmetry might not be explicit in V (η) as it is in V (φ).The high symmetry in tachyonic fields leads to the phenomenon of spontaneous symmetry breaking, where the potential has more than one minima i.e.V (φ) = V (−φ) corresponding to ±φ 0 .This phenomenon is well known in elementary particle physics [39].

B. Dp-branes
Recent advances in string theory have provided some explanations regarding nonperturbative features of QCD theory, thanks to the discovery of Dp-branes.The study of Dp-branes gives insight into physically relevant systems such as black holes, supersymmetric gauge theories, and the connection between Yang-Mills theories and quantum gravity.Upon the numerous progress, a consistent nonperturbative background-independent construction of the theory is yet to be developed.This poses a challenge in directly addressing cosmological problems.There are five known ways for which supersymmetric closed strings can be quantized, i.e.IIA, IIB, I, heterotic SO (32) and heterotic E 8 × E 8 superstring theories.Individually, they give a perturbative description of quantum gravity.However, they are connected through duality symmetries [42,43].Indeed, some of these dualities are nonperturbative in nature because the string coupling g s in one theory may have an inverse relation 1/g s with the other.The superstring theories together with M-theory (a theory that unifies superstring theories) consist of an extended object of higher dimensions, called Dp-branes.Each of these theories contains branes but in different complements.Particularly, IIA/IIB superstring theories consist of even/odd Dp-brane dimensional configurations.Using the appropriate duality transformations, one brane can be mapped onto the other even with the strings connecting them.Thus, none of the branes are seen to be more fundamental to others.This process is commonly referred to as 'brane democracy'.
Before moving into the specific kind of strings set out for this study, we will shed light on the features of open and closed strings.Closed strings are topologically the same as a circle S 1 , i.e. [0, 2π].They give rise to a massless set of spacetime fields, identified as graviton g µν , dilaton ϕ, antisymmetric two-form B µν and an infinite set of massive fields upon quantization.Supersymmetric closed strings are also related to massless fields identified with graviton supermultiplets such as Ramond-Ramond p-form fields A  [44,45].The open strings have their endpoints on the Dp-branes whilst close strings have no endpoints.The degree of freedom of the open strings is associated with standing wave modes on the fields while the close strings correspond to left-moving and right-moving waves.The boundary conditions for open strings on the bosonic field X M are Neumann (freely moving endpoints) and Dirichlet boundary conditions (fixed endpoints).The close string on the other hand corresponds to periodic and anti-periodic boundary conditions.Some Dp-branes have unstable configurations both in the supersymmetric or the nonsupersymmetric string theories.The instability is attributed to tachyonic modes with negative square mass m 2 φ < 0 in the open string spectrum, we are interested in investigating the effects of the tachyonic modes [27,46].Some Dp-brane configurations with open strings containing tachyons are: Brane-antibrane; it is a type IIA or IIB string theory with parallel Dp-branes separated by d < l s (l s is the string length scale).They also carry tachyons in their open string spectrum.The difference in their orientation leads to opposite Ramond-Ramond (R-R) charges [47].So the brane and anti-brane pair can annihilate leaving a neutral vacuum because the net R-R charge will be zero.
Wrong-dimension branes; the Dp-brane with wrong dimension for type IIA/IIB with odd/even spatial dimension for p instead of even/odd dimensions carry no charges under classical IIA/IIB supergravity fields.They have tachyons in their open string spectrum.Such branes can annihilate to form a vacuum without violating charge conservation.
Bosonic Dp-branes; just as the wrong-dimension branes of type IIA/IIB string theory, the Dp-brane of any particular dimension in the bosonic string theory have no conserve charge and has tachyons in their open string spectrum.Also, they can annihilate to form a neutral vacuum without violating charge conservation.
Again, even though the non-BPS Dp-branes of type IIA/IIB string theory are unstable owing to the presence of tachyonic modes on their worldvolume [23].We can obtain stable non-BPS branes by taking orientifolds/orbifolds of the theory which projects out the tachyonic modes [23,24,87].

III. MODIFICATION OF DIRAC-BORN-INFELD ACTION
We begin the study with the Born-Infeld action (BI) [48] sometimes referred to as Dirac-Born-Infeld action (DBI) [49,50].The focus will be on Dp-branes [51,52] which are nonperturbative states on which open strings live.They are equally coupled with closed strings, Ramond-Ramond states and other massive fields.The nonperturbative nature of the action makes it possible to describe low energy degrees of freedom of the Dp-branes [53] making it possible for application in low energy QCD.The distinction between the DBI action, other p-branes and supermembrane theories [47,54,55] is the presence of gauge field in the worldvolume of DBI.The gauge field is associated with virtual open string states.Consequently, we obtain confinement of the fundamental strings with their endpoints fixed on the branes.Generally, the action can be expressed as where G µν , B µν and ϕ are the induced metric tensor, antisymmetric tensor, and dilaton field to the Dp-brane worldvolume respectively.Also, is the worldvolume electromagnetic field strength of A µ and S CS is a set of Chern-Simon terms while is the brane tension and g s = e ϕ is the string coupling, with associated string tension In type IIA/IIB string theory with p even/odd are associated with quantum open strings containing massless A µ , µ = 0, 1, • • • , p and X M , M = p + 1, • • • , 9 fields.These fields are the consequence of the gauge field living on the hypersurface and X µ (ξ) transverse excitations.The geometry of Dp-brane is not flat, so we generally define the embedding X µ (ξ), where ξ α represent p + 1 coordinates on the Dp-brane worldvolume (p+1) , and X µ is the ten functions mapping from (p+1) onto the spacetime manifold R 9,1 .Introducing the scalar field into the action, it becomes invariant under diffeomorphism and Abelian gauge transformations.So, a way of fixing the freedom of the former is to adopt a 'static gauge' such that The remaining fields are which are the transverse coordinates to the worldvolume [58][59][60] where d is a spatial dimension.Thus, one can choose d = 9 for superstring theory and d = 25 for bosonic string theory.Under this gauge, the originally d + 1-dimensional global poincaré symmetry spontaneously breaks down to a product of p + 1 dimensional poincaré group with d − p dimensional rotational symmetry group i.e.SO(1, d) → SO(1, p) × SO(d − p).For a Dp-brane possessing p < d extends over p-dimensional subspace of d-dimensional space.The focus will be on a Dp-branes with p-dimensional hyperplanes in d-dimensional space.Again, according to the static gauge fixed above, there are two possible consistent truncations; • X M = 0; corresponds to pure BI theory [48] in E p,1 and • F µν = 0; also corresponds to Dirac's theory [50] of minimal timelike submanifolds of E d,1 .
We can introduce the transverse scalar fluctuations by defining the induced metric as that will approximate the Dp-brane worldvolume to near flat.In this study, we are interested in understanding the dynamics of bosonic open strings with tachyonic modes in their spectrum, so we will proceed systematically toward that objective.The worldvolume theory of non-BPS Dp-brane in IIA/IIB string theory corresponds to a massless U(1) vector field, transverse oscillating scalar field X M and a tachyonic field T [61].Accordingly, the leading order of the action results in dimensional reduction 10-dimensional U(1) Yang-Mills theory.Besides, higher-order corrections, α = l 2 s , to the order of the string scale are also possible.Since we are proceeding with the assumption that the massless fields are slowly varying compared to the string length l s , i.e. we discard the higher order derivatives and write the DBI action [53,56,62] in a simple form without the explicit presence of the tachyons.Therefore, Eq.( 3) takes the form Now we will express the extended form of the DBI by including the dynamics of the tachyons as studied in [56].Accordingly, where V (T ) is the tachyon potential and ∂ µ T ∂ ν T is, as usual, the kinetic energy of the tachyons [63].Under this conjuncture, the action vanishes at the minimum of the tachyon potential V (T 0 ) = 0 [19,22,25].We will now continue the discussion by approximating that the Dp-branes are nearly flat, with constant dilaton, and vanishing antisymmetric two-form term B µν to remove the close string quanta in the system.

A. Dimensional Reduction U(1) Yang-Mills Theory
Studying Dp-branes under Yang-Mills theory enables us to understand the physics of Dp-branes without necessarily applying any complex string theory artifacts.Detailed analyses show there is enough evidence that super Yang-Mills theory carries a lot of information concerning string theory than one may possibly imagine [64].Besides, recent developments in high-energy physics investigations have shown that string theory gives insight into low-energy field theories in the nonperturbative region [59].This has been conjectured to be equivalent to QCD where confinement can be realized in color fluxtube picture.From Eq.(9), in the limit of vanishing B µν and further assuming that the Dp-branes are almost flat, close to the hypersurface i.e.X M = 0, M > p. Again, we suppose that the fields are slowly varying such that ∂ µ X M ∂ µ X M and 2πα F µν are in the same order.Therefore, the action can be expanded as Here, V p is the p-brane worldvolume and g Y M is the Yang-Mills coupling given by, The second term in Eq.( 11) corresponds to U(1) gauge theory in p + 1 dimension coupled with 9 − p scalar fields.Introducing fermions fields, as mentioned below Eq.( 3) into the action, we recover the supersymmetric U(1) Yang-Mills theory in 10 d N = 1 Super Yang-Mills action, In the case of N parallel Dp-branes, the p-dimensional branes must be distinctively labeled from 1 to N. Subsequently, the massless scalar fields living on the individual Dp-brane worldvolume are related to U(1) N gauge group.The fields arising from the strings stretching from one brane to the other are labeled as A µ i,j , where i, j specifies the individual branes that carry the endpoints of the strings.The strings are oriented such that they consist of N (N − 1) fields, corresponding to A α = (A µ , X M ) individual fields.The mass of the strings is proportional to the separation distances between the branes.So the strings become massless [36,52] when the Dp-branes get very close to each other.Hence, the open strings can transform under the adjoint representation U(N).Thus, the corresponding fields can be decomposed similarly to adjoint supersymmetric gauge field U(N) = SU(N ) × U(1) in p + 1-dimensions.With the stacks of Dp-branes crossing each other at some angles, one can also break the U(N) symmetry into SM particle physics i.e.SU(3) × SU(2) × U(1) [65].In sum, the DBI action can be generalized into a non-Abelian gauge group by considering stacks of Dp-branes instead of a Dp-brane.
A major motivation for string dual to QCD is derived from the 't Hooft large N c -limit [66,67].Though, SU(N c ) and U(N c ) have different representations, in the limit of N c → ∞, the difference can be overlooked.Thus, the number of gluons can be approximated as N 2 c .This is more than the quark degree of freedom N f N c , therefore, we expect the dynamics of the gluons to dominate in this regime [68].

IV. BOSONIC STRINGS AT TACHYONIC VACUUM
Generally, there are two known boundary conditions associated with open strings.The Dirichlet (fixed) boundary condition is where the coordinates are normal to the brane, and under Neumann (free) boundary condition is where the coordinates are parallel to the brane [53,69,70].It has been established [71] that a small disturbance normal to the string and the brane are likely to reflect back with a phase shift of π, corresponding to a Dirichlet boundary condition.Some study in this regard has been carried out in [72,73] using Nambu-Goto action for strings with their endpoints on a supergravity background of D3-branes.Since the strings attached to the 3-branes manifest themselves as electric charges, the Neumann boundary condition where the endpoints of the strings are freely oscillating will lead to the production of electromagnetic dipole radiation at the asymptotic outer area of the brane.
In this section, we will consider Eq.( 10) which includes the dynamics of tachyons on open strings.Accordingly, we will keep all the arguments made for the 10-dimensional DBI action, but we will reduce the spacetime dimension from 10 to 4. One of the major differences between superstring theory and particle theory is that the former lives on 10-dimensional spacetime and the latter on 4-dimensional spacetime.Nonetheless, this discrepancy can be dealt with using compactification scheme [74][75][76], where the spacetime is divided into an external non-compact spacetime M 10−d and an internal compact spacetime M d .We can combine these two phenomena into a single expression as We adopt a physically realistic phenomenon where d = 6.Additionally, we can set a compact scale M c = 1/R, where R is the radius of the internally compact spacetime, smaller than the string mass M s = 1/l s [15,77].As a result, the energy E required for this work must lie in the range E M c M s .So, the worldvolume coordinate becomes ξ = (x µ , x 6 ), where µ = 0, 1, 2, 3. Consequently, we will reduce the dimension to 1 + 3-dimension with η µν = diag(+, −, −, −) metric signature.We will also decouple the 6 available transverse fluctuating scalar fields and the antisymmetric tensor i.e., X M = B µν = 0. Keeping the dilaton field constant, we get In the above expression we calculate the determinant up to the second-order derivative discarding higher-order corrections.In view of field theory at tachyonic vacuum, the D-branes do not vanish completely rather, there are lower dimensional D-branes present.So, the endpoints of the strings sitting on the low dimensional D-branes behave like point-like particles which serve as the source and sink of the flux carrying the color particles.For instance, expanding around the 'true' minimum of the potential gets rid of the tachyons leading to a particle with a positive square mass.We consider a field configuration such that T (r) = f (φ(r)), in this case, the potential can be expressed as As a result, the Lagrangian of the system becomes, here, we have introduced a dimensionless quantity G(φ) = (2πα ) 2 V (φ) which will be referred to as a color dielectric function, subsequently.If we set (2πα ) = 1 the string tension becomes T string = 1 and G(φ) = V (φ) [78][79][80][81][82].It is important to mention at this point that the potential of the field, φ, follows the same discussion as contained in Sec.II A. It has its minimum at V (φ = φ 0 ) = 0, where φ 0 is the 'true' vacuum of the potential.To apply this theory to asymptotically free systems, the potential must satisfy additional conditions, necessary for stabilizing its vacuum [81,83,84].These restrictions make it possible to apply Eq.( 18) in investigating asymptotically free particles such as gluons in a QCD-like fashion.Here, the modified Abelian gauge mimics the role of the non-Abelian gauge in the 'usual' QCD theory.We treat the endpoints of the open string on the lower dimensional branes as the source of quark and an anti-quark (if fermions are present) or valence gluons and the string connecting them as gluons that mediate the interactions.Also, the potential V (φ) plays a similar role as a quantum correction in gluodynamics theory that breaks the conformal symmetry to bring about gluon condensation [80,84,86].This model has been used to study the color confinement of glueballs at a finite temperature in Ref. [79].
Considering brane-anti-brane systems, we note that the model must be invariant under the global rotation φ → e iα φ.This warrants an introduction of a complex scalar field φ with potential V (|φ|) [88].So, we can redefine the scalar field in a form, to incorporate the Dp-brane and anti-Dp-brane dynamics.Supposing that the original gauge field F µν is on the worldvolume of the Dp-brane represented by a complex scalar field φ, we will have a dual gauge field Fµν also on the wordvolume of the anti-Dp-brane represented by the conjugate field φ * .Therefore, the string here has its endpoints on remnants of Dp-brane and the anti-Dp-brane parallel to each other.Imposing gauge invariance on the scalar sector, we can define an Abelian covariant derivative where Ãµ is a gauge field which is dual to A µ .Accordingly, the Lagrangian in Eq.( 18), can be extended as where Ãµ are two independent Abelian field strengths.Aside from the original U(1) gauge invariance of the Lagrangian, it is also invariant under the Ũ(1) gauge transformation, This Lagrangian can undergo a spontaneous symmetry breaking (SSB) process (when we choose an appropriate potential) similar to the usual Abelian Higgs mechanism [10,12].In this case, φ plays a role similar to the 'usual' Higgs field in the standard model of particle physics.This process will also lead to the observation of Goldston boson (most likely π 0 due to the Ũ(1)) corresponding to the number of unbroken symmetries signaling confinement as well [79,91].This model has been exploited in detail to study glueballs and color superconductivity in [82].It is important to mention that ( 22) is valid when we consider Dp-brane and an anti-Dp-brane system in the framework of DBI action [89,90].It is suitable for describing massless particles such as glueballs or gluon confinement.
V. GAUGE THEORIES MODIFIED WITH G(φ)

A. Fermions Coupled To Fundamental Strings at Tachyonic Vacuum
In this section, we introduce Dirac's Lagrangian for free particles adopted by Maxwell in the unification of electric and magnetic field interactions.We introduce the dynamics of the fermions which were dropped in the previous discussions as stated below Sec.III.However, it will be introduced here through Dirac's equation modified by G(φ) while taking into consideration all gauge invariance properties.We start with the well-known Dirac's Lagrangian for free particles Even though this Lagrangian is well known in QED, it is also used to describe free nucleons in terms of their composite fermions; protons, and neutrons in strong interaction.It is invariant under global phase rotation, Noticing that the fields in the Lagrangian are color neutral in nature, to apply them in studying color particles such as the ones considered here, we need to modify the fields.Hence, we will modify the bispinors with the G(φ) in order to give them some color features.Thus, we perform the transformations where G is a function therefore, the gradient in Dirac's equation will transform as and Lagrangian (24) becomes The local gauge invariance is violated by G(φ) and the gradient function (∂ µ G 1/2 ).To ensure local gauge invariance is satisfied, we modify all variables involving derivatives, and color neutral fields such as the Dirac γ-matrix with G(φ) and also introduce the electromagnetic field A µ (x) in other to make the equation gauge invariant.Consequently, we will adopt the transformations, where D µ is the gauge-covariant derivative and q is the electric charge.This type of derivative corresponds to momentum transformation, p µ → p µ − qA µ .By this transformation, the gauge field also gets modified and enters the Lagrangian as GA µ .It also introduces a coupling between the electromagnetic field and matter in a form D µ ψ.So Eq.( 28) becomes consequently, the electromagnetic field transforms as Equation ( 30) looks similar to Dirac's Lagrangian with interaction term that couples the gauge field A µ (x) to the conserve current j µ = q ψγ µ ψ with a modified mass term M (r) = mG(φ).Also, ψD µ ψ becomes invariant under local phase rotation, granting interaction between ψ, ψ and A µ with momentum p µ → i∂ µ .In this way, the Lagrangian has been modified, but we ensure that all conservation laws are duly respected.Now, to derive a complete Lagrangian that can mimic QCD theory, we include the kinetic energy term of the gauge field which has been derived in Eq.( 18), thus, This expression looks similar to the usual QED Lagrangian with φ as intermediate field/particle and a mass term modified by color dielectric function, G(φ).As a result, this expression approximates the non-Abelian QCD theory with an Abelian one.This was motivated by the projection that the confining regime of QCD is mostly Abelian dominated [92][93][94][95][96].The spinor fields ψ and ψ represent the quarks and antiquarks and the modified gauge field A µ also describes gluons.All the long-distance behavior of the gluons is absorbed in G(φ).The observed mass of the system M (r) = mG(φ) is expected to have a fixed value at the beginning of the interaction r → r 0 and at the end of the interaction r → r * , so it can be measured precisely.That is, the dielectric function should be such that G(φ(r → r 0 )) = G(φ(r → r * )) = constant, where 1/r 0 is the energy at the beginning of the interaction and 1/r * is the energy at the end of the interaction.Thus, M (r) = mG(φ) is the constituent quark mass function of the system.
We have presented a detailed study of this theory in Refs.[80,81] Since the renormalization theory remains the systematic approach for which UV divergences can be resolved [97], we will compare the result with the renormalized QED Lagrangian to properly identify the nature of the dielectric function in the context of renormalization factor Z(µ), where µ is a scale that comes from the dimensional regularization scheme [98,99].The objective is to ensure that the result obtained in Eq.( 32) does not pose any UV divergences.When the real QED Lagrangian is written in terms of the renormalized factors, it takes the form, where the two equations bear some resemblance, so we can compare them.With Z 3 the gluon propagator renormalization factor, Z 1 the quark-quark-gluon vertex renormalization, and Z 2 the quark self-energy renormalization factor.Additionally, the covariant derivative can be expressed in terms of the renormalized factors as We have substituted the 'conventional' representation of the renormalization factors Z 2 and Z 3 with Z ψ and Z A respectively, and to make the distinction more obvious relative to fermions and the gauge field.The renormalized fields are without subscript '0' [100][101][102][103][104][105][106].Comparing the results in Eqs.( 32) and ( 33), we identify Z A = Z ψ = G and the gauge invariance warrant that, Z A = 1.Thus, in addition to the color properties carried by the color dielectric function, it also absorbs the UV divergences.Consequently, the dielectric function follows the restrictions, From the gauge conditions adopted above, a gluon mass term, will not be invariant under the local gauge transformation in Eq.( 31) because Consequently, local gauge invariance accounts for the existence of massless photons [10].

B. Non-Abelian Gauge Theory
In this section, we will focus on constructing a non-Abelian gauge theory traditionally used for describing the strong nuclear force.Additionally, we buy into the projection that proton and neutron have the same mass and are charge independent due to the strong nuclear force.Therefore, there is an agreement that isospin is conserved in strong interactions.As in the case of QED theory, we will base the discussions on SU(2)-isospin gauge theory introduced by Yang and Mills [107], elaborated by Shaw [11].Here, we will consider Eq.( 24) as the Lagrangian for free nucleons with composite fermions, protons (p), and neutrons (n) The Lagrangian is invariant under global spin rotation also, isospin current, j µ = ψγ µ (τ /2)ψ, is conserved therefore, proton and neutron can be treated symmetrically in the absence of electromagnetic interactions.In that case, the distinction between proton and neutron is arbitrary and conventional under this representation.
To maintain the differences between the Abelian theory treated in Sec.V A and the non-Abelian theory intended for this section, we will replace A µ from the Abelian covariant derivative, D µ , expressed in Eq.( 29) with its non-Abelian counterpart B µ i.e., We have also replaced q with g, the strong coupling constant to make the analysis distinct from Sec. V A and befitting for describing strong interactions.Despite the similar global gauge invariance satisfied by both theories, they also exhibit some major differences; the algebra of the non-Abelian group is more complex and the associated gauge bosons self-interact due to the structure of the non-Abelian group.Accordingly, the field can be expressed as, where  24) will then transform as, Again, following similar procedure as adopted for Eq.( 30) using the necessary transformations, we obtain with gauge invariant transformation similar to Eq.( 31), By comparing Eq.( 43) and Eq.( 45) we can deduce, The SU(2)-isospin generators can be expressed in terms of commutator relation, Conventionally, τ i do not commute with others in different spatial directions, so only one component can be measured at a time and this is taken to be the third component T 3 = 1/2τ 3 , generally in z-direction.
We will now develop the field strength tensor that will form the kinetic term of the gauge field.Starting with the electromagnetism gauge, we can construct, From these transformations, we can develop a gauge invariant field, where the field strength tensor transforms under the dielectric function as We know from the QED relation constructed in Sec.V A that the field strength can be expressed as We adopted the expression for B µ defined in Eq.( 43), regrouping the terms in the above equation yields, higher derivative terms were discarded.We can cast this result in a slightly symmetric form by using the identity Using this identity appropriately leads to, This equation shows the additional terms that come from the non-vanishing commutators due to the non-Abelian group structure.By this expression, we deduce that a term can be added to ∂ ν B µ − ∂ µ B ν to modify F µν to achieve the desired transformation property we require.With this inspiration, the observed electromagnetic field strength tensor can be modified to read, Substituting the definition for D µ in Eq.( 29) into the above expression, we obtain The commutator vanishes for Abelian theories.Applying this transformation to Eq.( 50) yields, Expanding the commutator using the transformation Eq.( 43) to enable us to compare the outcome to the nonvanishing commutator relations in Eq.( 54) leads to, The commutator relations that come after the first term in the last step are the exact terms required to cancel the extra terms in Eq. (54).Therefore, the field strength tensor expressed in Eq.( 57) has the required structure under local gauge transformation.Combining Eqs.( 44) and ( 49) gives rise to modified Yang-Mills Lagrangian [10,107], where M (r) = mG(φ) is the modified nucleon mass.This Lagrangian is also invariant under local gauge transformations and does not permit the existence of mass term M 2 (r)B µ B µ .It should also be noted that introducing the non-Abelian gauge leads to the automatic cancellation of the color dielectric function.Hence, we can infer that the color dielectric function attached to the Abelian gauge induces strong interaction properties.Using the expression in Eq.( 41) and the commutator relation in Eq.( 47), we can rewrite the transformed version of Eq.( 57) as in the last step we have dropped the three isospin generators τ l of the gauge field because they are linearly independent.Generally, non-Abelian gauge groups that do not fall under SU(2), the Levi-Civitá symbol ε jkl is replaced with the antisymmetric structure constant f jkl .The introduction of the SU(2)-isospin symmetry by Heisenberg [118] preceded the development of the quark model.That notwithstanding, the only known fundamental components of the nucleon describing strong interactions are up-quark (u) and down-quark (d).While the proton is composed of two u-quarks and a d-quark, the neutron is also composed of two d-quarks and an u-quark.Indeed, these particles remain the only known constituents of proton and neutron with almost the same mass and coupling force.Hence, Eq.( 59) can be used to study the behavior of quarks and gluons inside hadrons.In that case, the up/down quarks are treated as having the same mass and 0-charge, so they are seen as similar particles with different isospin states, I 3u/d = ±1/2, same as the nucleon.Interestingly, the spin addition of the quark constituents of proton and neutron agree with I 3 = 1/2 and I 3 = −1/2 respectively.Similar to the isospin representation of the nucleon field ψ in Eq.( 38) the quark field can be represented as Nevertheless, the mathematical structure of this theory is the same as QCD, the theory that describes the characteristics of quarks and gluons inside the hadrons.Again, the strong interaction is well known for its invariance under quark color permutations, so we can express the quark wave function for color triplet representation as where R ≡ red, G ≡ green and B ≡ blue.This is invariant under SU(3) transformation of the form where α j , j = 1, 2, 3, • • • , 8 are the eight phase angles and λ j is a 3×3 matrix representing eight independent traceless Hermitian generators of the group.Ignoring the differences in quark masses, it can also be applied in studying u, d and s quark systems, The generators are fundamentally equivalent to Pauli's matrices for the SU(2) representation and they satisfy the Lie algebra, similar to Eq.( 47), f ijk is the structure constant.Consequently, to generalize the relation in Eq.( 60) for QCD under the SU(3) color symmetry, we substitute the antisymmetric tensor ε ijk for the antisymmetric structure constant f ijk [12].Finally, the models exhibit the expected asymptotic free [8,9] behavior at high energy regions while at low energies color confinement and hadronization sets in.Here, M (r) represents constituent quark mass function [119] while m is the bare quark mass.As analyzed in Sec.V A, the constituent mass of the quarks in the asymptotically free region can be determined as M (r → r 0 ) = M 0 , while the constituent mass at the nonperturbative (low energy) region where color confinement is expected, will be M (r → r * ) = M * .Hence, it is possible to have the same constituent quark mass in both the UV (r → r 0 ) and the IR (r → r * ) regimes, if indeed the bare quark mass m in both regimes are the same [81,120], because the G(φ) behaves as, G(φ(r → r 0 )) = G(φ(r → r * )) = constant.However, there is evidence that m is small in the IR regime and large in the UV regime [80,81].Accordingly, if we require color confinement in both regions, M 0 > M * because higher bare quark mass m is required to obtain confinement in the UV region than it is required in the IR region.

VI. PHENOMENOLOGY OF GLUEBALL CONFINEMENT
In this section we will use one of the models built in Sec.IV specifically, Eq.( 18) to build a model that describes the features of confining glueballs.Apart from the insight into the behavior of the glueballs in both the IR and the UV regions, it will serve as a test to the models developed earlier.The model will be based on electric field confinement, commonly referred to as the chromoelectric flux confinement.In that light, we will define the indices of the gauge field such that the chromomagnetic flux is eliminated from the system (i.e.j ν = (ρ, 0)) and only the static sector of the scalar field, i.e. µ = j, is available for analysis.This ensures that the color particles that generate the gluons are static.This section will enable us to see how the Abelian gauge can be used to approximate a non-Abelian theory.Additionally, the color dielectric function G(φ) absorbs the long-distance dynamics of the gluons such that the photon propagator emanating from the Abelian gauge field does not decouple at longer wavelengths.We will further demonstrate that the G(φ) is directly related to the QCD β-function and the strong running coupling constantly.

A. The Model
The equations of motion for Eq.( 18) are, and Expressing the above equations in spherical coordinates, where, F µν F µν = F 0j F 0j + F i0 F i0 = −2E 2 (only electric field components are considered) and V = G as discussed in Sec.IV.Accordingly, where and λ is an integration constant, we also substituted ε 0 = 1 to achieve the desired objective.Now we choose a potential that satisfies all the conditions expressed in Sec.IV thus, where ρ and α = 1/f α are dimensionless constants, f α is the tachyon decay constant.This potential contains tachyonic mode at V (0), so the fields cannot be quantized around this point -see Sec.II A for detailed discussions.However, we can remove the tachyonic modes by shifting the vacuum φ → φ 0 + η and quantizing around the true minimum, φ 0 = ±a/α.Consequently, and the potential can be expanded for the small perturbation η, which will be referred to as a glueball field with a real square mass, m 2 φ , hence, Considering that the particles are sufficiently separated such that color confinement can be observed, we ignore the 1/r 4 term in Eq.( 69) and simplify it as leading to This equation has several solutions but we choose two of such solutions suitable for the analysis, Each solution corresponds to the characteristics of the particle in a particular regime i.e., IR and UV regimes respectively.The solution at the IR regime will give rise to a linear confining potential and the UV solution will lead to a Cornell-like potential.

B. Confining Potentials
In this section we will present the confining potentials derived from the model by considering the electrodynamic potential, Substituting the equation at the left side of the solution Eq.( 77) and Eq.( 70) into Eq.(78) leads to where c is an integration constant that is set to zero in the last step.Considering that m φ r 1, c = 0 and λ = 1 corresponding to the positive part of the potential V , we can deduce the QCD string tension σ to be, Here, we can infer that confinement is occasioned by the magnitude of the glueball mass, in the limit of vanishing glueball mass there will be no confinement.Furthermore, this potential leads to linear growth in r, and at some critical distance, r * = 1/ √ σ L [108] the potential begins to flatten up leading to hadronization.It is known from the flux tube models for confining color particles that, r r * [108].On the other hand, taking the solution at the right side of Eq.( 77) and following the same process as followed above, we get in the last step, we set the integration constant c = 0, also choosing the positive part of the potential corresponding to λ = 1 and ρa 4 = 1, we arrive at the Cornell-like potential for confining heavy quarks i.e., corresponding to a string tension It is important to recognize that the critical distance, r * s = 1/ √ σ s , in this regime marks the transition from the asymptotic freedom region to the confining region.An increase in the distance also leads to an increase in the strength of confinement.
We know from Sec. III that the string tension T string ∼ σ L ∼ σ s ∼ 1 GeV/fm, hence the glueball mass in the IR regime becomes m L ≈ 1 GeV corresponding to glueball mass of isoscalar resonance f 0 (980) [109].The commonly known ratio of m(0 ++ )/ √ σ L in QCD theory in SU(∞) limit [110][111][112], at this regime, can be determined as Likewise in the UV regime, we have a glueball of mass m s ≈ 1.73 GeV corresponding to the lightest scalar glueball mass of resonance f 0 (1710).The result obtained here is precisely the same as the results obtained from QCD lattice calculations [109,[113][114][115][116][117]] also, m s / √ σ s ≈ 1.73.The critical distances become r * s = r * = 1 fm for both the IR and the UV regimes.While critical distance in the IR regime r * refers to the transition from confinement to hadronization regions, critical distance in the UV regime r * s refers to the transition from the asymptotically free region to the confining region.

C. Gluon Condensation
Classical theory for gluodynamics is invariant under the scale transformation x → λx, this leads to a scale current s µ which is related to the energy momentum tensor trace θ µ µ (x) as In the absence of quantum corrections, θ µ µ = 0, the theory remains conformally invariant.This will lead to a vanishing gluon condensation F a µν F aµν = 0. On the other hand, when quantum correction, −|ε v |, is introduced, the conformal symmetry is broken leading to non-vanishing gluon condensate F a µν F aµν = 0 and energy-momentum trace anomaly comes to play with vacuum expectation, The leading term of the QCD β-function is known to be, Now, calculating the energy-momentum tensor trace of Eq.( 18) for the glueball field η using the relation, we get, where G and V represent first derivative with respect to η and Ṽ (η) = V (η) − ηV (η)/4.Also, rescaling Ṽ (η) with the energy density −|ε v | i.e.Ṽ (η) → −|ε v | Ṽ (η) together with the vacuum expectation value in Eq.( 86) we get, with this equation, we recover the classical result in the limit |ε v | → 0 i.e.F µν F µν = 0. Using the potential expressed in Eq.( 74) we can determine, as a result, Eq.( 90) can be expressed as we can identify the gluon mass [82], Furthermore, taking the expectation value of Eq.( 66) in terms of the glueball field η, we can express consequently, the mean glueball field η has two possible solution i.e. η = 0 and η = 1.Higher glueball condensate corresponds to η = 0 whilst lower glueball condensate corresponds to η = 1 [121].Comparing Eqs.( 85) and ( 89), we can relate We can also extract the strong running coupling α s using the renormalization group theory [122], comparatively; Therefore, the strong coupling can be identified as α s (η) = G(η) = α s (1/r 2 ).QCD β-function is naturally a negative quantity showing the asymptomatic freedom nature of the strong coupling.It also reveals the anti-screening behavior of the theory at higher energies.The strong running coupling, on the other hand, gives an insight into the growing precision of hadron scattering experiments at high energy limits.And at low energy limits, within the scale of hadron mass, it enhances understanding of hadron structure, color confinement, and hadronization.Now, substituting the solution of the glueball field η and expanding it for r → 0 we obtain, Also, we can associate the spacelike momentum Q with Q ≡ 1/r, then In terms of the four-vector momentum i.e.Q 2 ≡ −q 2 , the strong coupling becomes and the β-function becomes, We observe that in the limit q 2 → 0 the strong coupling shows a singularity, generally referred to as the Landau singularity.It marks the failure of perturbative QCD.The singularity is attributed to the self-interacting gluons with hadron degrees of freedom in the IR regime leading to color confinement [123].In that case, the gluons dynamically acquire mass at q 2 → 0 i.e. q 2 ∼ = m 2 A [124,125] which increases the coupling infinitely.Hence, the singularity can be removed by fixing a freezing point [124,126] to the strong running coupling at q 2 ∼ = m 2 A i.e., and Thus, the 'so called' gluon mass is more pronounced at q 2 → 0 and its effect gradually fades off in the limit where q 2 → ∞ [127].A recent analysis of this subject is contained in [122].
FIG. 4: Strong Running Couple α s (left) and β-Function (right) against q, with a Landau Ghost Pole The graphs show an unphysical behavior at q → 0, this is due to the presence of dynamically generated gluon mass.The self-interacting gluons and the strong force that exist between them are capable of creating bound states with hadron degrees of freedom.These graphs depict the behavior of αs and β observed from pQCD.
FIG. 5: Strong Running Coupling α s (left) and β-Function (right) against q with a freezing point at q → 0 (a) Left Panel Here, the Landau singularity has been fixed by introducing the gluon mass m A .The presence of the gluon mass is more pronounced at q → 0 and gradually vanishes in the limit q → ∞.Consequently, α(0) 2.3 and β(0) −4.7 from the graphs.

VII. CONCLUSION
We modified the DBI action to develop models that are capable of mimicking the phenomenon of QCD theory using an Abelian gauge field.The models were based on the behavior of opened string with its endpoints on the Dp-brane.In studying color particles, the endpoints of the string serve as the source and sink of the color charges.Additionally, the models are efficient in investigating glueballs when the tachyons condense and transform into glueballs with real square masses that keep them confined.Without fermions, the models are suitable for studying the bound states of gluons and the dynamics of glueballs.To study the dynamics of quarks, we showed how the model can be coupled with standard model fermions systematically.Here, the particles involved are glueball-fermion-mix in a confined state.Moreover, we demonstrated that the color dielectric function coupled with the Abelian gauge capable of causing color confinement vanishes automatically when we introduced the non-Abelian gauge field in Sec.V B. Consequently, the presence of G(φ) coupled with the Abelian gauge field was to induce non-Abelian characteristics.
We also developed one of the models to demonstrate its ability to explain some basic characteristics of strong interactions.We derived the linear and Cornell-like potentials that are used to describe color particles in phenomenological QCD.The linear potential is motivated by the string model of hadrons whilst the Cornell potential is motivated by Lattice QCD calculations.The Cornell potential is particularly important in QCD because it shows both the asymptotic freedom and color-confining behavior exhibited by the model.We also calculated the strong running α s coupling and the QCD β-function and compared their behavior with the traditional QCD theory.However, in the model framework, we are able to fix the non-physical Landau ghost pole that occurs at the low energy region of the model by assuming the existence of gluon mass m A at low momentum region, q 2 ∼ m 2 A .The model leads to the determination of gluon condensate and how the glueball fields contribute to the condensate.
Furthermore, the models can be discretized using path integral formalism and investigated under lattice field theory with the availability of the required computational artifacts.As observed in the model developed in Sec.VI, the linear and the Cornell-like potentials can be used to study hadron and quarkonia spectra.Other hadron properties can also be studied from these models when the appropriate spin contributions to the potential are added.Extending the models to study the characteristics of particles at a finite temperature will pave way for understanding, chiral symmetry breaking and restoration, confinement/deconfinement, and quark-gluon-plasma phase transitions.Finally, the models can be applied in investigating physical systems such as pions.
•• ,µp and gravitino ψ µα .Consequently, the quantum theory of closed strings is naturally related to the theory of gravity in spacetime.On the other hand, open strings are topologically similar to the interval [0, π].They give rise to massless gauge field A µ in spacetime upon quantization.Supersymmetric open strings are also associated with massless gaugino ψ α .Accordingly, open strings have their ends on Dirichlet p-branes (Dp-brane) and the gauge field lives on the worldvolume of the Dp-brane.Upon these differences, the physics of closed and open strings are related at the quantum level.Closed strings were first observed through the one-loop process of an open string.Under this process, close strings appeared as poles in nonplanar one-loop open string diagrams the three non-Abelian gauge fields are b µ = (b 1 µ , b 2 µ , b 3 µ ).The generators τ a (a = 1, 2 and 3) are associated with Pauli's matrices, b ± µ = (b 1 µ ∓ ib 2 µ )/ √ 2 are the charged gauge bosons and the isospin step-up and step-down operators 1/2(τ 1 ± iτ 2 ) exchange p ↔ n following the absorption of b ± µ boson.The gradient in Eq.(