Steps towards generalization of tensionless string theory with contact interactions as Wilson loop of non-Abelian Yang–Mills theory

We propose a possible modification to the tensionless string model with contact interactions. The proposed model aims to reproduce the expectation value of a non-Abelian Wilson loop in Yang–Mills theory by integrating out string degrees of freedom with a fixed worldsheet boundary. To reproduce path-ordering along the worldsheet boundary, we introduce Lie algebra-valued fields on the string worldsheet, whose dynamics are determined by the topological BF action. Without bulk contributions, we show that the model describes the non-Abelian Wilson loop, neglecting the effects of self-interactions. Finally, we test the reproduction of the Wilson loop with three-point interaction in the case of SU(2).


Introduction
Yang-Mills theories serve as important building blocks of the standard model.They are considered one of the most successful theories of fundamental particles ever tested experimentally.There are several viewpoints suggesting possible connections between Yang-Mills and string theories.Dating back to the seventies, it was known that string theory reproduces tree diagrams of Yang-Mills theory in the limit of low energies or infinite string tension [1].This perspective proves useful for understanding the structure of amplitude relations for both string theory and field theory.Many structural aspects in the amplitudes of both theories were widely studied [2][3][4][5][6][7][8][9].Another connection has been noticed by 't Hooft [10] that in the large N limit, a pure SU (N ) Yang-Mills theory diagram matches a perturbative expansion of a string theory with the string coupling constant 1/N .
Quite recently, the formulation of Abelian Yang-Mills theories as the tensionless limit of a spinning string with contact interactions was formulated in [11] and [12].In addition, by using the worldline approach, one can relate the theory of interacting strings to Wilson lines in spinor quantum electrodynamics in the tensionless limit [13].The idea of this formalism was initiated by Mansfield [14] which expresses an electromagnetic field strength tensor of two moving charges as a string-like object supported on a worldsheet Σ bounded by particle worldlines expressed in the form where dΣ µν is an infinitesimal area element on the surface defined as X µ (ξ) is a spacetime vector with the worldsheet coordinates ξ a , a = 1, 2 lying in the domain Σ with boundary ∂Σ on which X µ ∂Σ = ω µ and ∂ a is the derivative with respect to world-sheet co-ordinates ξ a .The expression (1) satisfies the Gauss' law, ∂ µ F µν = J ν with a current density supported on the worldsheet boundaries Substituting (1) in the pure electromagnetic action, S = 1 4 d 4 xF µν F µν , gives Due to the appearance of the delta function, the action is non-vanishing when the worldsheet coordinates coincide, ξ = ξ, or when any two points coincide, i.e.X(ξ) = X( ξ).This splits (4) into two pieces as S EM = q 2 4 δ 2 (0)Area(Σ) + q 2 4 Σ dΣ µν (X(ξ))δ 4 (X(ξ) − X( ξ))dΣ µν (X( ξ)) The first piece contains the area of the worldsheet corresponding to the Nambu-Goto action of string theory, albeit with a divergent coefficient.The latter piece provides a more interesting interpretation.It implies a contact interaction that occurs when the worldsheet intersects with itself.This is unusual from string perspectives in which the standard interactions in string theory are caused by joining and splitting worldsheets.Similar interactions have previously been discussed by Kalb and Ramond [15].From now on we denote the second term of (5) by S I [X].
By replacing the Nambu-Goto action with the classically equivalent Polyakov action it can be shown perturbatively [16] that the partition function of a tensionless four-dimensional string with the contact interaction S I whose worldsheet Σ spans the closed loop ∂Σ is similar to the Wilson loop for Abelian gauge theory associated with the closed curve ∂Σ in flat Euclidean space at the first leading order.
To illustrate this, let us determine an expectation value of the contact interaction S I Σ where we defined Ω Σ , the worldsheet average of a quantity Ω over all surfaces Σ spanning ∂Σ, as with a normalization constant Z P .According to the expression of S I in (5), we apply the Fourier decomposition to the delta function to obtain where V µν k (ξ) is the vertex operator defined as The fact that the theory lives in non-critical dimensions may concern some readers as it usually leads to the notorious Weyl anomaly.Luckily, this issue is cured by the presence of Dirac-delta function.It was shown in [12,14] that the expectation value of the delta-function decouples from the scale of the worldsheet metric.Therefore, this allows us not to integrate over g ab but rather choose a value for it.
The expression (8) can be separated into two terms by introducing a projection operator P k which is defined as The operator projects any 4-vectors onto their transverse directions compared to the vector k.Therefore, the vertex operator now takes the form where Ṽ µν k (ξ) is a projected vertex operator defined as Note that as k µ P k (X) µ = 0, the vertex operator must satisfy k µ V µν k = 0. Consequently, when inserting (11) into (8), we find To obtain the above expression, Stoke's theorem was used.It turns out that averaging over the worldsheet using the standard string action will suppress the first term.To see this, we use Wick's theorem to evaluate the expectation of products of fields.According to Wick's theorem for the bosonic string, G(ξ, ξ ′ ) is the Green function for the worldsheet Laplacian.The colons indicate normal ordering and since the field X can be expanded around the classical field X c , so the expectation value of the normal ordered part is It is not hard to find that the expression for the projected vertex operator is The Green's function G(ξ, ξ ′ ) can be evaluated from the heat kernel G via where the Laplacian △ is given by The general form of the heat kernel can be written using the Seeley-DeWitt expansion [17] which can be modified to a manifold with boundary [18][19].If σ r (ξ, ξ ′ ) is twice the square of the length of the geodesic path connecting between ξ and ξ ′ with r reflections at the boundary, then the heat kernel is obtained by writing The function Ω r can be expanded as a power series of τ which is where a r n (ξ, ξ ′ ) are called the Seeley-DeWitt coefficients.According to [18], for ξ = ξ ′ , the coefficients of the first few orders are evaluated as where R is the Ricci scalar.It is not difficult to see from the expression (20) that the heat kernel is divergent at a small τ , thus the value of Green's function at co-incident points should be regulated with a short-distance cut-off, ǫ, via an integral In this limit as τ → 0 and for ξ ≈ ξ ′ , it is sufficient to obtain the asymptotic version of (20) by including only zero and one reflection terms as At the co-incident points, σ 0 = 0. Therefore, G ǫ (ξ, ξ) reads As a result, the Green's function at co-incident points G(ξ, ξ) is zero on the worldsheet boundary and diverges as ξ moves away from the boundary into the interior of the worldsheet.Since the theory was in the Euclidean signature, k 2 > 0, the projected vertex operator ( 16) is suppressed inside the worldsheet for which α ′ k 2 is finite.This suppression gets further amplified when taking the tensionless limit α ′ → ∞ into consideration.What remains in the expectation of S I is that of the second term in ( 13) To discard the bulk term in (13), one needs to take good care of some possible divergences appearing when the vertex operators are placed close to each other either in the bulk or near the boundary of the worldsheet.A detailed analysis can be found in [12].The expression ( 26) is exactly the integral over ∂Σ of the Abelian gauge field propagator in the Lorenz gauge The subscript YM notifies that the expectation was made in the Abelian Yang-Mills theory.This is the first non-trivial term in the perturbative expansion of the Wilson loop exp (−q This suggests that the expectation value of the Wilson loop could be expressed as the worldsheet average of the exponential of S I .However, a difficulty arises as divergences appear when exponentiating S I which potentially spoils the suppression.Fortunately, no such terms are produced in the supersymmetric generalization of the model.It appears that the expectation value of the super Wilson loop for (non-supersymmetric) Abelian gauge theory can be expressed as the worldsheet average of the spinning string with a contact interaction [11], [12].
In this paper, we would like to construct a suitable modification to the string model to reproduce the expectation value of a non-Abelian Wilson loop in the Yang-Mills theory.The difficulty which impedes the non-Abelian generalization is incorporating Lie algebras into the theory.This can be done by simply introducing Lie algebra-valued degrees of freedom into the worldsheet.However, one needs to find suitable dynamics to describe them.These new degrees of freedom have to generate the interaction vertices of Yang-Mills theory and reduce to path-ordered generators on the boundary.If successfully formulated, the theory would provide an alternative channel to investigate a non-Abelian Yang-Mills theory.
The outline of this paper is as follows.In section two, we develop the tensionless string model to incorporate a boundary path-ordering of Lie generators and at the end of the section a partition function for the non-Abelian generalization of the string model with contact interactions is proposed.In section three, the evaluation of the partition function is calculated on the worldsheet boundary.We review an effective BF theory in section four as well as compute a worldsheet gauge propagator.In section five, we evaluate a functional average of the square of string contact interactions which is partial to the partition function and compare the result with the expectation value in the SU (2) gauge theory of the Wilson loop.Finally, we summarise our results in section six.

Production of Boundary Path-ordering of the String Model
To generalize the string model to describe the non-Abelian Yang-Mills theory, at the very least we need to introduce Lie algebra-valued world-sheet degrees of freedom φ R (ξ) to try to reproduce the Lie algebra structure of Yang-Mills propagators Consequently, we need a Lagrangian to describe the dynamics of φ R .This has to be gauge-invariant to preserve the spacetime gauge invariance of the contact interaction (29) and Weyl invariant to satisfy the usual organizing principle of string theories.In addition, the newly introduced degrees of freedom are expected to reduce to path-ordered generators on the boundary to represent a non-Abelian gauge theory.We propose that the action which meets the minimum criteria is the topological BF action [20,21], i.e.
with the field strength tensor The action (30) has a close connection to the Yang-Mills action in two dimensions as they are equivalent in the zero coupling constant limit [22,23].Remember that both fields φ and A are Lie algebra-valued fields.They can be written in terms of a set of generators {T R } as φ = φ R T R and A = A R T R . 1 The worldsheet 1-form A is intrinsic to the worldsheet and it differs from the actual spacetime gauge field A in the gauge theory whose dynamics we wish to reformulate.It can be shown that when we define the partition function corresponding to the BF action as one can obtain a boundary correlation function of scalar fields as path-ordering of Lie generators.A Wilson loop is inserted along the boundary of Σ.To remove all the gauge redundancy, we apply the axial gauge-fixing condition via the insertion with a fixed vector n.Therefore, the partition function (31) takes the form where the covariant derivative D µ is defined as ∂ µ + A µ .To obtain this, we used the fact that the integrand and the measures are gauge invariant which can then be renamed from (φ Λ , A Λ ) to (φ, A).By introducing a source term for the scalar field, one can construct a generating functional as Integrating out the field φ generates the constraint via δ-function as For simplicity, we will content ourselves to consider the worldsheet with the topology of a disk D 2 which is topologically equivalent to the upper halfplane.Accordingly, it is convenient to work in Cartesian coordinates, hence, (35) taking the form where A y ] and we chose the reference vector n to be a unit vector pointing in y-direction.
After integrating out A y , the generating functional (36) becomes This requires us to solve the constraint, ∂ y A x = − 1 2 J.By setting where I is an index that represents different insertions of the field φ(x I ) on the worldsheet and λ I is a functional dependent on x I , one can obtain the solution as where θ(x) is a Heaviside step function.
On the boundary, the generating function takes the form Therefore, an expectation value of a product of φ on the boundary is Although we are able to incorporate the boundary path-ordered generators into the string model, the object on the left-hand side of (41) is still not what we desire as it is not gauge invariant which implies that the result depends on what gauge conditions we utilized.Accordingly, to make the observable gauge independent, we insert Wilson lines W C defined as where C is an arbitrary curve whose endpoints are x 1 and x 2 .Therefore, one can modify the observables (41) to be The appearance of trace and Wilson lines is to make the product gauge invariant.
To evaluate (43), let first consider just two insertions of φ at the points x 1 and x 2 with a Wilson line joining both points, i.e.
in the axial gauge condition where we set A y = 0 at each point.Using the expression (39), the Wilson line W C (x 1 , x 2 ) is simply an ordered product of group elements when the contour C intersects with the vertical lines S 1 and S 2 as illustrate in the figure 1. Surely, the value of the Wilson line depends on the path C we are considering.For the given path shown in the figure, it Fig. 1: A Wilson line connecting from the point x 1 to x 2 reads exp (−λ(x 1 )) exp (−λ(x 2 )).Note that the sign of the exponents rely on the orientation of the path C when cutting through the vertical lines.
If the insertion points x 1 and x 2 are moved to the boundary, the Wilson line will reduce to one as no path can go below the inserted points.As a result, the boundary expectation value of (43) is Remember that the value (44) does not depend on which gauge choice and the contour paths we chose.Consequently, we propose a modification to the string contact interaction S I (5) for incorporating the boundary path-ordered generators as This form of action enjoys gauge symmetry at the worldsheet level.To avoid picking a particular path of the Wilson lines, we average over all possible paths using the approach [16].The expectation value over all the paths C is given by with where h rs is the induced metric along the path.
We then construct the functional integral to be with a suitable normalization factor Z. This functional integral is simply an expectation value of exp(−S φ,A I ) whose dynamics are described by the string action S P and the topological BF action S BF with a Wilson loop inserted along the boundary.
In this paper, we aim to show that the expectation of the non-Abelian Wilson loop in gauge theory W [∂Σ] is equivalent to the functional integral Ψ.However, we will not be able to give a full demonstration of the equality between the two models due to complications of potential divergences but rather show that the model (48) contains the required ingredients to reproduce the non-Abelian gauge theory Wilson loop.A more complete proof needs to wait until the supersymmetric model is developed.

Evaluation of Ψ: Boundary Contributions
In this section, we would like to evaluate the partition function Ψ along the worldsheet boundary ignoring any contributions from the worldsheet interior.
To evaluate Ψ, we first expand the exponential of the contact interaction S φ,A I into the power series of S φ,A I as Let first consider the first non-trivial term of (49), i.e.
To evaluate this, we rewrite S φ,A I (29) in terms of vertex operators V µν k (ξ) as Similar to the Abelian calculation, the projection operator P k in (10) was used to write with a generic function F (ξ).If we give where Note that for the above expression we denote the primed and un-primed variables to be the objects determined on the point ξ ′ and ξ respectively.However, the terms from the third to the eighth lines of (53) can be discarded as the contraction between U and V contains k µ V µν which is zero.For simplification, we introduce the vertex operators as follows: where we use the parameter τ to parameterize the worldsheet boundary on which the field X µ is ω µ and denote ω as ∂ω/∂τ .With these definitions and the expression (53), we can rewrite the contact interaction (51) as Notice that the terms which involve the vertex operators C and D are not present in the Abelian model.When averaging (56), it is likely that the bulk terms will get suppressed due to an appearance of e ±ikX in the interior of the worldsheet.This is because the Wick theorem of the exponential vanishes in the worldsheet's interior at the tensionless limit.Remember that the Green's function at coincident points G(ξ, ξ) is divergent in the worldsheet interior as discussed earlier.However, we do not expect that all the bulk terms get completely suppressed by this manner as there is possible that some contractions among the bulk terms might provide singularity to negate the suppression.This issue will be explored in the section five where a contraction from the worldsheet gauge fields provides a singularity necessary to reproduce a self-interaction in the Yang-Mills theory.
Considering only the boundary terms, the boundary expectation of the contact interaction is where we give ω ≡ ω(τ ) and ω ′ ≡ ω(τ ′ ).According to (44), one finds Consequently, the expectation of S φ,A I takes the form where S I Σ is the expectation of the Abelian contact interaction expressed in (26).The expression (60) is the exact prescription to reformulate the expectation value of the non-Abelian Wilson loop omitting self-interaction terms.This can be seen by evaluating the boundary contributions of the expectation value of e −S φ,A I which is At the order q 2n , the expression (61) describes a Wilson loop with n pairs of gauge propagators which freely propagate between the boundary.The expression of the non-Abelian Wilson loop is presented in Appendix A. Remember that the calculation was evaluated on the worldsheet boundary ignoring the effect of bulk terms in the contact interaction on which potential divergences may arise.In the Abelian case, such terms are absent by introducing supersymmetry on the worldsheet.However, a similar generalization of the non-Abelian case has not yet been achieved.

Effective BF Theory and Worldsheet Gauge Propagator
To verify if the proposed string model provides a valid description of the non-Abelian Yang-Mills theory or not, we need to find out whether the model can reproduce self-interactions of the gauge fields.According to the previous section, it is obvious that no such terms appear when evaluating solely on the worldsheet boundary.We expect that such terms can be obtained from the worldsheet gauge fields A in the worldsheet interior.To demonstrate this, it requires the computation of a worldsheet gauge field propagator.This can be obtained by inspecting the generating functional We continue the calculation by integrating out the gauge field A. To do this, we expand the fields φ and A i in the Cartan-Weyl basis as where H a and E α are the Cartan and Weyl generators respectively satisfying the following algebra: and where H α is defined as H α = α a H a .The Cartan generators H a are diagonal traceless matrices in the adjoint representation.Note that these bases are ξdependent.The Cartan generators were chosen such that the field φ lies within their subalgebra at every point on the surface.
To relate Lie indices A with the Cartan and Weyl indices a and α, we introduce unit vectors Ĥa A and Êα A in Lie vector space which are defined as δ a A and δ α A respectively.As a result, the inner products among the vectors are and the completeness relation is It is not hard to write the field φ and A i in terms of the unit vectors as Using the relations (67), one can find the topological BF action (30) taking the form Notice that there is no contribution from diagonal components of A A i to the first term as the Cartan subalgebra is commutative.The last term only contributes to the boundary of the worldsheet.Consequently, this gives the generating function (62) as with To obtain the above expression, we consider the case that the worldsheet Σ has the topology of a disk, accordingly, one can parametrize the worldsheet coordinates such that ξ 1 is zero at the boundary.In this way, the coordinates ξ 1 and ξ 2 can be seen as radial and angular-like coordinates.
The path integration of the second line gives a constraint on the theory in the form of the Dirac delta function with m the number of Cartan generators.Without the source J , the constraint implies that the square of the field φ is constant throughout the worldsheet interior with a jump at the boundary, i.e.
We then proceed the calculation by changing the worldsheet coordinates (ξ 1 , ξ 2 ) into the complex coordinates (z, z) where z = ξ 1 + iξ 2 and z = ξ 1 − iξ 2 .In these new coordinates, the field a α i becomes complex fields b α where The partition function now resembles a Gaussian path integral with respect to the fields b and b which is where φA is a unit vector in Lie indices space of the field φ whose value is ĤA a H a .We can then use the Gaussian integration formula to integrate out the complex field b, To apply (76) into (74), one can use a general expression for an inverse matrix ( M −1 ) α β which is where δ j1j2...jn i1i2...in is a generalized Kronecker delta that is related to an antisymmetrization of ordinary Kronecker deltas as The integer n is the number of Weyl generators.In the case of SU (N ), n is equal to N 2 − N .This yields where the matrix component M αβ is given by and the matrix Θ AB relates to the adjugate matrix adj( M ) by the relation Turning back to the (ξ 1 , ξ 2 ) coordinates, the effective action takes the form When J = 0, the expressions (81) and (84) are the effective actions of the BF theory provided earlier in [24] in which the general expression for the effective Lagrangian was presented in a diagrammatic representation.One can include a boundary Wilson loop W ∂Σ into the generating functional (62) by introducing auxiliary either commuting or anti-commuting fields onto a worldline to generate the path-ordering of the Lie algebra generators [25][26][27].According to [28,29], a trace together with a path-ordering operator can be replaced by a functional integral over the Grassmannian field ψ as where the loop ∂Σ is now parameterized by τ .Therefore, a suitable choice of the source term J to regain the notion of the Wilson loop W ∂Σ is Hence, to align with our model (48), we modify the generating functional (62) to be where D(ψ † , ψ, φ) = Dψ † DψDφ.We can then obtain a two-point function of the fields A by It turns out that only the gauge propagator in the Weyl directions survives because when taking the functional derivative of Z BF [J ] with respect to J am with a Cartan index a, it would appear terms containing Θ aB J B in the integrand which is equal to zero due to the commutative property of the Cartan generators.Remind that a derivative of the Dirac delta functionals with respect to J gives the path integral zero as they contain only the source terms in the Cartan directions.For this reason, the double functional derivative of (87) takes the form where G i ( J ) is the constraint defined as tr(( Taking the term in the square bracket into consideration, we obtain If we consider the gauge propagator in the interior of the worldsheet, the two-point function becomes To obtain the above relation, (91) and the fact that the source term inside the worldsheet is zero were used.For the case of SU (2), one can find that the matrix Θ AB = −ǫ ACB φ C and the matrix determinant is 2|φ| 2 .More detail can be found in [24].Thus, (92) becomes In this case, the effective action for SU(2) BF theory is Note that to get the second line in (93), the constraint G i (0) = 0 was implemented which implies that |φ| is constant throughout the worldsheet Σ.

Evaluation of Ψ: Three-point Self-interactions
In this section, we would like to determine whether or not our string model can re-create self-interaction vertices of the gauge theory.To demonstrate this, we provide a perturbative computation of Ψ up to a few orders of q where the effects of self-interactions can be observed.
For simplicity, we content ourselves to work in the simplest case of the non-Abelian model, i.e.SU (2).Since a contribution of the three gauge-boson vertex in the Wilson loop is first observed at O(q 4 ), we expect that the expectation Fig. 2: The simplest Wilson loop with the three gauge boson vertex of (S φ,A I ) 2 contains the expression for the Wilson loop with three-point interaction as illustrated in the figure 2. The question is which terms can generate such a structure.As mentioned earlier that it is highly probable that the bulk terms will get suppressed in the tensionless limit unless there exist any singular terms which potentially spoil the suppression.According to (92), we got a hint that Wick contraction of the worldsheet gauge fields A generates singularity at a coincident point.This could result in forming a new vertex operator at the coincident point.By this argument, we speculate that the three-point vertex interaction to appear in the expectation of where the vertex operators were defined in (55).To evaluate the right-hand side of (95), we begin by considering the products One can show that these products contain the following terms: whose values are singular at the point ξ = ξ ′ through the Dirac delta function when Wick contracted using (93).We can then apply Wick's theorem for the worldsheet gauge field A to determine the singular behavior of the above products.This turns (96a), (96b) and (96c) to be respectively.The above expressions are merely the singular parts among other non-singular terms resulting from the Wick's theorem.Note that there are still gauge fields left uncontracted in the Wilson lines which we will cope with them later.In addition, we have omitted the contribution of the second line of (93) as it will be suppressed by the matter fields X µ when functional averaged in the tensionless limit.Consequently, to obtain the expression for the right-hand side of (95), we need to evaluate the following integrands which are and Substitute (97a) into the first integrand, one obtains To obtain the last expression, we used the fact that k • dP k (X) = 0, hence Similarly, one can find the two remaining integrands, i.e. (98b) and (98c), as and To proceed the calculation, we need to evaluate Wick's contraction of the Wilson lines.For simplicity, we are not going to provide a full expression for that contraction but rather evaluate Wick's contraction of their Taylor expansions keeping up to the first non-trivial terms which are (102) We will show that keeping the series expansion in this way is enough to provide a structure to reproduce the three-point interaction of the Wilson loop.Accordingly, let's consider the double integral on the right-hand side.Using (93), one obtains where which counts the number of times the two curves intersect in an oriented way.
To obtain (103), we choose to work in the gauge choice such that the unit vector φ is constant everywhere in the worldsheet interior Σ but its value varies along the worldsheet boundary ∂Σ.If the two curves C 1 and C 2 lie themselves on the upper half plane with the boundaries on the x-axis whose endpoints are b 1 and b 2 respectively, according to [16], the average over with proportional constant α.Remind that the average over curve was defined in (46).The Wick's contraction of the curve-averaged Wilson loops reads Remember that |φ| 2 is constant throughout the worldsheet.As a consequence, applying (106), one finds the expectation value of the integrand (98a) as where and Similarly, we can find the expectation values of the integrand (98b) and (98c) as and where and As we have everything set up, we shall proceed the calculation of the gauge field expectation of (95), 1  2 (S φ,A I ) 2

A
, by substituting (107), (110), and (111) into (95).This yields where k ≡ (2π) 4 .To obtain the above equation, we utilized the fact that the integral d 4 k is invariant under a sign flip, i.e. k → −k, thus one can rename the variables k and k ′ in ( 110) and (111) to have the same integrand regarding the matter field X µ as in (107).It is not difficult to see that the terms in the squared bracket in the third line sum up to zero as tr(T A T B T C ) is all anti-symmetric in the case of SU (2).The same situation also happens for the terms Θ 1 + 2Θ 2 + Θ 3 whose value sum up to zero.This can be seen using the expression for a trace of a product of four generators [30], i.e.
Therefore, what remains is to compute the integrand (118) To evaluate the quantities Ω, we use a formula for trace of five SU (2) generators A more detailed derivation of the above relation can be found in Appendix B. As a consequence, the integrand (118) now takes the form where We then further simplify (120) using the property of the Levi-Civita symbol which is Therefore, the integrand (120) becomes Substitute (123) back to (116), it yields For the next step, we would like to apply the functional average for the field X to the right-hand side of (124).To do so, we will change the form of the bulk integral into a new boundary vertex operator as we learned in the first section that the expectation inside the worldsheet gets greatly suppressed by an exponential of the Green's function at the co-incident points G(ξ, ξ).This can be done by expanding the integral using the projection of X (10) along (k +k ′ ).Therefore, there exists the term in which we can turn it into an integral along the boundary using integration by parts which is where ωµ (λ) = dω µ dλ (λ).Substituting (127) to (124) and relabeling some dummy parameters, one obtains To relate the right-hand side of the above expression to the expectation of the Wilson loop with three-point self-interaction at the lowest order (A10), we expect that the expectation of the product of φ would provide a notion of path-ordering of the trace of generators.As a consequence, what we need to do next is to determine the expectation of 1 |φ| 2 φ A φ B φ C concerning the remaining fields, i.e.
According to (87), a functional average of an object Ξ with respect to the fields ψ † , ψ, φ is defined as where and G i ( J ) is a boundary constraint defined in (90).The partition function Z[φ, ψ † , ψ] was placed to normalize the expectation which is defined as If we consider a worldsheet of which the topology is equivalent to that of a disk, we can map the worldsheet into an upper-half plane with a boundary on the x-axis.The expectation (129) becomes where and x i = x(τ i ).Remember that this path integral was taken along the worldsheet boundary on which the effective action S SU(2) eff vanishes.We found that if we imposed a boundary condition such that a derivative of the field φ in the normal direction equals zero, i.e. ∂ y φ = 0, the expectation (132) gives path-ordering of the trace of the product of three generators.To illustrate this, let's consider the integral Dφ δ(tr((2∂ x φ) φ))δ(tr((2φδ(y) + J y ) φ)) where the boundary condition was already imposed.Taking a closer look at the constraints, the first constraint δ(tr((2∂ x φ) φ)) implies that a modulus of the scalar field φ is constant throughout the boundary and the second one δ(tr((2φδ(y)+ J y ) φ)) tells us that φ(x) is required to equal to 1 2|φ| ψ † (x)T ψ(x).Note that we are not interested in the solution where φ = 0.
Accordingly, we proceed with the calculation by changing the variables of the measure Dφ and the Dirac deltas to be with the Jacobian factors J 1 and J 2 .It is not difficult to check that both J 1 and J 2 do not depend on either |φ| or φ.Exploiting this fact, it is not necessary to find the exact forms of the Jacobian factors as they will eventually cancel out with ones in Z[φ, ψ † , ψ].Consequently, (132) takes the form 1 where the Jacobian factors J 1 and J 2 have been omitted as discussed and Z[φ, ψ † , ψ] takes the form By integrating out the unit scalar field φ, (138) becomes 1 Notice that the product in the square bracket is evaluated at every point on the boundary except for x = x 3 while the integrand in the second line is at x = x 3 .It is obvious that the integral over |φ| in the first line will get canceled out by the same integrals in Z[φ, ψ † , ψ] for every point x = x 3 .However, at the point x 3 , the integral over |φ| reads This is exactly the integral in the parenthesis of (139) when relabelling c → 1 c so they also cancel out.In consequence, this turns the expectation (129) to be which can be interpreted as according to [28,29].
All things considered, we find that the expectation value of (S φ,A I ) 2 in the proposed string model contains a term which is the form we need to relate to the expectation value in the gauge theory of the Wilson loop at the O(q 4 ) (A10).Nevertheless, we may leave the constant α undetermined as it is probable that similar expressions can also be obtained when including higher order of corrections in Wick's contraction of Wilson loops.Remember that to obtain the expression (144), we only keep Wick's contraction of the Wilson loops up to the first non-trivial term.
It is undeniable that the computation presented in section five relied heavily on the properties of SU (2) group generators.To generalize to SU (N ), one needs to handle more complicated relations among SU (N ) group generators, which involve both structure constants f ABC and totally symmetric third-rank tensors d ABC .We speculate that the three-point vertex could be reproduced from the expectation of the same term in equation (95).
However, the model is still incomplete without a reproduction of fourpoint self-interactions.We speculate that the four-point vertex term could be obtained from the expectation value of (S φ,A I ) 3 as the four-gluon vertex is first observed at the order q 6 .Needless to say, the computation would be significantly more difficult and will need to be left for future investigation.

Concluding Remarks
To conclude, we proposed a possible non-Abelian generalization of the [11,12] via the introduction of Lie algebra-valued worldsheet fields into the string model.The dynamics of the new degrees of freedom are described by the topological BF action.In our model, the string theory has unusual contact interactions defined in (45) which enjoy local gauge symmetries.By inserting a Wilson loop along the worldsheet boundary, we showed that a boundary path-ordering of Lie generators can be reformulated through the expectation in equation (43).
The partition function of the modified tensionless string model with contact interactions was presented in (48).When neglecting the effects of the bulk terms, the partition function is equal to the expectation of a non-Abelian Wilson loop in the usual gauge theory when self-interactions have been turned off.
We then verified that our string model correctly reproduces the three-point self-interaction of SU (2) Yang-Mills theory.We performed this calculation by perturbatively expanding the exponential of the contact interaction up to the second order.We found that at that order, the expectation of (S φ,A I ) 2 contains a precise expression for the expectation value of the SU (2) Wilson loop with three-point interaction at the lowest order.To obtain this result, we used a worldsheet gauge propagator calculated in the section four.
Despite encouraging results, we cannot rush to the conclusion that our string model provides a true description of non-Abelian Yang-Mills theory.There are still a few issues to be further considered.Firstly, we have not yet investigated whether the model can reproduce a four-point interaction, and if it does, we expect that such a vertex can be found from the expectation of (S φ,A I ) 3 .More importantly, it is still unexplored whether the model reconciles with general Lie algebras.Secondly, in this work, we ignore potential divergences that may occur when performing functional integration over the matter field X when any two string vertex operators are placed very close to each other inside the worldsheet.In the Abelian case, this issue was addressed by including supersymmetry on the worldsheet [11,12].Developing a supersymmetric generalization of the tensionless string model to describe non-Abelian super Yang-Mills is a subject for future investigation.Finally, it is also possible that our model will generate additional interacting vertices not found from the Wilson loop in the non-Abelian theory, which could imply potential deviation between the two models.For a small coupling constant q, we can evaluate (A3) perturbatively which requires the calculation of A n .The first non-trivial contribution to W is dξ µ 1 dξ ν 2 A µ (ξ 1 )A ν (ξ 2 ) .(A4) The above term can be written as The main difference between the Abelian and non-Abelian theories is the existence of self-interactions.There are three and four gauge boson vertices which first appear in the expectation of the Wilson loop at q 4 and q 6 respectively.For the three-point vertex, we can calculate its contribution to the From (B12), one finds In the case of SU (2), if we identify T A = 1 2 σ A where σ A are the Pauli matrices, the metric η AB is simply δ AB , thus there is no distinction between upper and lower indices.Moreover, in this scenario, f ABC = ǫ ABC and d ABC = 0.As a consequence, one obtains (B16) Remember that the generator T A is traceless.Readers may consult [30,32] for more detailed discussion on relations among the SU (N ) generators.
Appendix A Non-Abelian Wilson loop of the Yang-Mills TheoryThe Wilson loop can be defined as the trace of the path-ordered exponential of a line integral of the gauge field A along a closed loop C,W [C] = tr P e −q C A•dξ (A1)where P is a path-ordering operator which orders a product of operators along the path C asP(O 1 (ξ 1 )O 1 (ξ 2 ) . . .O N (ξ N ) ≡ O P1 (ξ P1 )O P2 (ξ P2 ) . . .O PN (ξ PN ) (A2)where on the right hand side the operators are ordered according to their positions along the path, i.e. ξ P1 ≥ ξ P2 ≥ . . .≥ ξ PN .The trace in (A1) is computed over colour indices.By Taylor expanding the exponential (A1), the expectation value of the Wilson loop isW [C] = tr P