Low energy spectrum of SU(2) lattice gauge theory

Prepotential formulation of gauge theories on honeycomb lattice yields local loop states, which are exact and orthonormal being free from any spurious loop degrees of freedom. We illustrate that, the dynamics of orthonormal loop states are exactly same in both the square and honeycomb lattices. We further extend this construction to arbitrary dimensions. Utilizing this result, we make a mean field ansatz for loop configurations for SU(2) lattice gauge theory in 2+1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2+1$$\end{document} dimension contributing to the low energy sector of the spectrum. Using variational analysis, we show that, this type of mean loop configurations has two distinct phases in the strong and weak coupling regime and shows a first order transition at g=1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$g=1$$\end{document}. We also propose a reduced Hamiltonian to describe the dynamics of the theory within the mean field ansatz. We further work with the mean loop configuration obtained towards the weak coupling limit and analytically calculate the spectrum of the reduced Hamiltonian. The spectrum matches with that of the existing literature in this regime, establishing our ansatz to be a valid alternate one which is far more easier to handle for computation.


Introduction
Understanding the low energy behaviour of gauge theories is one of the most important problem of particle physics. Formulation of gauge theories on discrete lattice [1] has shown many major investigations in this direction over past few decades. Using Monte Carlo method, many important physical quantities can be computed numerically [2,3]. However, understanding the vacuum as well as excited states in this sector is still open for investigation. Hence, there should always be attempts to make analytic approximations. This present work proposes such an approximation for SU (2) lattice gauge theories in 2 + 1 dimensions, which can as well be generalized to higher dimensions and higher gauge groups. a e-mail: tpir@iacs.res.in We work within the Hamiltonian formulation of lattice gauge theory [4] and use prepotential [5][6][7][8][9] framework, which gives an useful reformulation in terms of gauge invariant loops. Loop formulation of gauge theories is always desired by theorists [10][11][12][13][14][15][16][17][18] as one can get rid of spurious gauge degrees of freedom. However, working with loop does not guarantee to work with only physical degrees of freedom as the loop space itself is highly over-complete [19][20][21]. Working in terms of gauge invariant loops in the weak coupling regime is again particularly difficult, as all possible loops of all shapes and sizes do contribute to the low energy spectrum of weakly coupled gauge field theories. In this scenario, the prepotential formulation [5][6][7][8][9] gives a great advantage over the standard Wilson loop approach as it is possible to extract only physically relevant loop degrees of freedom and study the dynamics of those.
In prepotenntial formulation of lattice gauge theory [5][6][7][8][9], one constructs gauge invariant loop variables, locally at each site. This particular feature is extremely useful for analytic calculations [22][23][24][25][26][27] and also for the recent progress in quantum simulating lattice gauge theory [28,29]. However, the original motivation for formulating prepotential approach was to make better understanding at the weak coupling regime of the theory. Towards this direction, a very important step is to construct the exact and orthonormal loop Hilbert space, containing only physical degrees of freedom. In this work, (also in another parallel recent work by Anishetty et al. [27]), we have proposed a general technique of constructing explicit orthonormal loop states for SU (2) lattice gauge theory in any arbitrary dimension. This is done even without going into the complicated Clebsch-Gordon coeffiecients specific to SU (2) and hence is generalizable to SU(3) in a straightforward way.
The prepotential formulation on square lattice reveals that the physically relevant loops are the non-intersecting ones [22,23]. We propose that, if one virtually splits each site of the lattice into two virtual ones following Fig. 8, the resulting lattice is an hexagonal one in two dimension. It has also been demonstrated that the dynamics of orthogonal loop states on square lattice is exactly equivalent to that of all possible loops on hexagonal one. The prepotential formulation on hexagonal lattice keeps all the important features of this particular formulation intact, such as local loop description by constructing intertwiners and Abelian weaving of those intertwiners leads to standard Wilson loops of the theory [5][6][7][8][9]. Moreover, the extra advantage on hexagonal lattice is, that the local loop space constructed at each site is exact and orthonormal. This is a tremendous advantage for the purpose of any computation be it analytic or numerical as one needs to work within a really small Hilbert space without bothering about complicated Mandelstam constraint. In [27], it has also been shown that there exists, more than one way to split each point, and the resulting lattices turn out to be of different types (i.e hexagon,octagon, square etc) as well as of different translational symmetries. It has also been argued in [27] that, dynamics on these virtual lattices are exactly equivalent to that of the original lattice by gauge fixing on the virtual links connecting splitted lattice sites. One can choose any splitting scheme as per the calculational convenience. For the purpose of present work we fix the splitting scheme given in Fig. 8 at each site, and get a hexagonal lattice to work with. We establish its equivalence with the original lattice by computing the dynamics explicitly for random loop configurations generated on both square and hexagonal lattice.
Within this framework, we propose a mean field configuration describing the low energy spectrum of the theory. It is true, that the mean field analysis cannot take us to the continuum limit, but is worth studying as it provides an excellent analytical tool to understand the vacuum structure of the for g < 1 and to study the low energy dynamics in this regime. In this work, we start from prepotential formulation of pure SU(2) gauge theory on square lattice. Then using the virtual point splitting technique of [27], we move to the virtual hexagonal lattice. On this particular lattice, we make a mean field ansatz that only an average loop configuration contribute to the low energy spectrum of SU(2) lattice gauge theory. We also show that this average loop configuration has two distinct phases at the weak and strong coupling regime of the theory and shows a first order phase transition at g = 1 denoting two distinct vacuum at two regimes of the theory. We also construct a reduced Hamiltonian, starting from the Kogut-Susskind Hamiltonian, which keeps the dynamics of the theory within the mean field ansatz. We further perform a variational calculation, to fix the mean field configuration at different values of coupling. This analysis confirms that the average fluxes flowing across each site shows a distinct jump as one moves from strong to weak coupling regime at g = 1.
We are however interested in the the loop configuration at small values of g. This analysis reveals that large fluxes contribute to the low energy spectrum at weak coupling regime as The orientation denoted in this figure defines the convention to construct the local loop operators and states in (7) opposed to the zero flux at strong coupling regime. For the original Kogut-Susskind Hamiltonian, largest contribution should come from large loops carrying large fluxes. Working with prepotentials makes us free from considering large loops at all, as all loops are now local [5][6][7][8][9]22,23,27]. Hence, our ansatz of the vacuum in weak coupling regime, consists of only large fluxes flowing across the sites. We finally compute the lower lying spectrum of that reduced Hamiltonian at weak coupling regime and show that we get reasonably acceptable results within this approximation.
The organization of the paper is as follows: in Sect. 2 we briefly discuss the loop formulation of SU(2) gauge theory on hexagonal lattice and compare its dynamics with that of the square lattice. In Sect. 3, we illustrate the the origin of hexagonal lattice from square lattice by virtual point splitting in arbitrary dimension. In Sect. 4, we discuss the average loop configuration for prepotential formulation of SU(2) theory in 2+1 dimentions on virtual hexagonal lattice. In Sect. 5, we propose a reduced Hamiltonian and discuss its dynamics within the mean field ansatz and finally compute low energy spectrum within this ansatz. Finally we summarize and discuss future aspects of this study in Sect. 6.

Loop formulation of SU(2) gauge theory on hexagonal lattice
Let us consider SU(2) pure gauge theory, formulated on 2 dimensional spatial lattice consisting of hexagonal plaquettes as shown in Fig. 1. As in the Kogut-Susskind Hamiltonian formulation [4], each of the links (n, i) originating from site Link operators consisting of two parts in prepotential, one of which increases flux along that link by one unit and another one decreases n along direction i 1 carries a link variables U (n, i) and there are SU(2) generators E a (n, i), a = 1, 2, 3 present at each end of the links. The Hamiltonian, in terms of these canonically conjugate variables reads as: Within prepotential framework [5][6][7][8][9], we attach e set of prepotential doublet (as we are working with SU(2)) a † α (L)&a † α (R) at left and right end of each link (n, i) with α = 1, 2 and i = 1, 2, 3. In terms of preotentials [5][6][7][8][9], the electric field is given by, satisfying the SU(2) algebra at each end. The link variable at site n along the direction i takes the form given by: where, N i = a † (L)·a(L) = a † (R)·a(R) counts the number of prepotentials along direction i. From (3), we find that, the link operator is a sum of a creation part, which increases the flux by one unit and another annihilation part, where flux is decreased by same. Hence, pictorially we can think of the link as shown in Fig. 2. For hexagonal plaquette this decomposition yields 2 6 plaquette terms (as opposed to 2 4 terms in the case of square plaquette). Likewise on square lattice, the prepotential formulation on hexagonal lattice, yields a local loop description of the theory which we will show to be free from any loop redundancy. Each of the plaquette operators basically consists of six local loop operators glued together following Abelian Gauss law. Let us first concentrate on the local loop operators and loop states constructed on at each site of the hexagonal plaquette. At each site of the hexagonal lattice, links emerge in three directions, labelled as 1, 2, 3. Each link is associated with a prepotential creation and annihilation operators where, i, j are the direction indices with i = j. These local loop operators on the hexagonal lattice are represented pictorially in Fig. 3. It is clear from the above set of equations that, acting on strong coupling vacuum only the first operator, i.e the one given in (6) will give non-zero contribution and will build up the local loop Hilbert space. For two dimensional hexagonal lattice, we characterize the local loop space by three independent linking numbers l i j denoting the flux flowing along three (i j) directions, namely (12), (23) and (31) as shown in Fig. 4. Hence, the most general loop state at each site is given by: where, N is the normalization factor. Note that, on hexagonal lattice, only three linking numbers are present at each characterizing a complete basis. The local and orthonormal loop states at each site is explicitly obtained as, The |l i j basis on hexagonal lattice is exactly equivalent to the number operator basis, where three number operators n 1 , n 2 , n 3 counts the number of prepotentials along each direction. These two basis are related as: n 1 = l 12 + l 31 , n 2 = l 12 + l 23 , n 3 = l 23 + l 31 (9) or equivalently, l 12 = 1 2 (n 1 + n 2 − n 3 ) , From the above set of equations, note that, the l i j variables are truly independent positive integers ranging from zero to infinity, whereas the transformed n i variables not truly independent but derived from l i j basis following (9). Unlike l i j variables, the ranges of n i 's are restricted to satisfy the triangle inequalities given in (10) keeping the l i j 's positive semidefinite. In this regard, working with the l i j basis is more convenient as it involve no further constraint unlike the n i basis, however one can always move from one to another following (9) or (10). Next is to find out the action of the local loop operators (4,5,6) on orthonormal local loop states (8).
These actions are far more simple than that in case of square lattice, derived in [22,23]. The first and simplest operator (6) involves no annihilation operator, and hence simply gives: Note that, the coefficient is obtained to create another normalized state from the one on which it acts. This loop action is pictorially represented in Fig. 5a showing that the resultant flux along i j direction increases by one unit. The next loop operator in (6) involves only annihilation operator and its action on a general loop state is: The right hand side of (13) is obtained by shifting both the annihilation operators of O −− i j to the right using commutation relations and adjusting the normalization factor as given in (8).
Another loop operator (6) involves one creation as well as one annihilation operator and in the very same way its action on arbitrary loop state is obtained as: Note that, as there are only three links emerging from each site, for a given operator O +− i j , the direction k, on which there will be change in flux configuration, but no change in net flux is always fixed, and hence there is no summation over index k, in the right hand side of (13). The last two loop actions (12,13) are also represented pictorially in Fig. 5b, c respectively.
Let us now concentrate to the complete plaquette, around which the Hamiltonian dynamics evolves. As stated earlier, plaquettes are the smallest Wilson loops and the basic gauge invariant variables to appear in the Kogut-Susskind Hamiltonian. But in terms of prepotentials, it is not the fundamental one but consists of six vertices as shown in Fig. 1 each consisting of the full set of loop operators and states. All possible loop operators at each vertex which are given in Fig. 6, which are of type listed in (4)(5)(6). Their action on loop state as derived above is listed in detail in Table 1 for convenience.
One important thing to note here is that, as per our convention of defining loop operators and loop states, some of the loop actions bring −ve sign in the coefficient as shown in Table 2. But this happens only for the mixed operators (i.e type O +− ) which involves one creation and one annihilation operator. But when we consider a closed loop such as a plaquette, we see that these type of mixed vertices always appear in pairs. Moreover note that, for the full plaquette operator (or any closed loop), the mixed terms can only appear in pairs and hence, the plauette operators (or any closed loop operators) are always positive by our convention.
Further note that, the action of loop operators on any loop state consists of two parts, one contains a numerical coefficient or number operators and another is some shift operators for the linking numbers. The coefficient that appears in Table 1 are calculated follows from the convention that, the shift operators are always at right most position and coefficients (function of number operators) are at left.
Having set the action of the loop operators on arbitrary loop states, one can easily compute the matrix element of the magnetic part of the Hamiltonian (1) within orthonormal loop basis, characterized by l 12 , l 23 , l 31 (or, n 1 , n 2 , n 3 ) basis. The magnetic Hamiltonian consists of 2 6 = 64 terms, each of which is a set of six local loop operator at each of the six vertices of the hexagon, the action of which on respective local loop states are computed following the Table 2.

Dynamics on hexagonal lattice vs. dynamics on square lattice
At this point we compare the Hamiltonian dynamics on hexagonal lattice with that on the square lattice numerically. For this purpose, we generate an arbitrary but valid loop configuration around one particular hexagonal plaquette of the lattice. This is done by specifying a set of three positive semidefinite integers denoting n 1 , n 2 , n 3 at each of the alternate site (say at site a, c and e of the plaquette in Fig. 1). We denote these alternate sites as the even sites of the lattice. The set of three integers when satisfy triangle inequalities are accepted as a valid loop configuration. For the neighbouring odd sites, i.e for sites b,d and f, two numbers are  Fig. 6 Pictorial representation of all possible loop actions around a plaquette 'abcde f ' given in Fig. 1 fixed by even sites a, c, e, and the third one is generated randomly satisfying triangle inequalities. The prefixed ones are: and randomly generate n 2 (b), n 3 (d) & n 1 ( f ) satisfying triangle inequalities at b,d and f sites as well. 2 Having gener- Table 1 Loop actions around the hexagonal plaquette, given in Fig. 6. The coefficients in the loop action are given in the last column Sl no. At vertices Action on |l 12 , l 23 , l 31 Coefficient 3.
c and f

5.
a and d O ++

7.
a and d

11.
b and e Table 2 Few sample results for the comparison between dynamics around square plaquette and hexagonal plaquette. The first and second and column denotes two loop states around a hexagonal plaquette between which there exist a non zero matrix element for plaquette term Tr U plaquette of (1). These configuration has been translated to a valid orthonormal loop configuration around a square plaquette using the pre-scription given in this section and the matrix element M s is calculated following Appendix A. The same M h is also calculated for hexagonal plaquette following (16). Results are shown upto 8 decimal places upto which these two matches exactly. Note, these are only few sample results from numerical simulation

Initial state
Final state Square plaquette Hexagonal plaquette ated a loop configuration denoted by n 1 , n 2 , n 3 at each of the six sites around the plaquette, we readily compute the corresponding l 12 , l 23 , l 31 at each of them following (10). All possible loop operators residing at each of the six vertices around the hexagonal plaquette, changes the loop configurations following a coefficient listed in Table 2. The dynamics of an arbitrary hexagonal plaquette in compact form is given below: Footnote 2 continued Abelian Gauss law (21) is satisfied is another option. It requires a detail study to find out which one is the most efficient one. Here, as we are only interested in dynamics around a chosen plaquette, we do not bother to generate loop configurations throughout the lattice.
where, C ± i for i = 1, 2, ..., 6 's are some algebraic coefficients which are functions of number operators. In the full plaquette operator, they always come in pairs as in the six vertices of the plaquette. Product of two such C ± i s are the vertex coefficient C j for j = 1, 2, .., 12 listed in Table 2. We compute those coefficeints for our configuration to find out  Table 2.
Having calculated numerical values of the matrix elements of the magnetic Hamiltonian for hexagonal plaquette, we will now compare those matrix elements with that for a square plaquette (Matrix elements of magnetic Hamiltonian on a square lattice is given in Appendix A).
An important point to note at this point is, even for a square plaquette, the loop configuration actually changes over a virtual hexagonal plaquette consisting of four real fluxes flowing around a plaquette together with two internal fluxes (coupled angular momenta, see Appendix A for detail). We provide the dictionary of identifying square plaquette to hexagonal plaquette in Fig. 7.
The dictionary for shifting between square (given in Appendix A) and hexagonal plaquette (given in Sect. 2) is as follows: -Each site (s ≡ a, b, c, d) on the square lattice has the following orthonormal angular momentum to have net angular momentum zero at site x: -Identify the flux around a square plaquette abcd as: -At each site (s ≡ a, b, c, d, e, f ) of the hexagonal plaquette, the orthonormal states are characterized by l x 12 , l x 23 , l x 31 or equivalently by following (9). 3 -Identify the flux (marked with 'h') around the hexagonal plaquette 'abcdef' as: -Now, the same flux around the hexagonal plaquette, can be identified as the dynamic flux around square plaquette 'abcd' in the following way: -Moreover, the external links of plaquettes 'abcd' and 'abcdef' can also be identified, but we are not writing them explicitly as we find them to remain unchanged in this particular plaquette dynamics.
Having established the connection between the particular square and hexagonal plaquette of interest we can now compute the matrix elements for magnetic Hamiltonian for both the cases. We now compare the above calculated dynamics around an hexagonal plaquette with that of the square plaquette given in (A.6). For this purpose we identify the fluxes around the hexagonal plaquette with those around a square plaquette as given in Fig. 7. To compare the dynamics on square lattice and that on the hexagonal lattice, we simulate a random loop configuration on hexagonal lattice and compute the matrix element of the magnetic Hamiltonian as discussed above. Next we identify the same loop configuration on square lattice following prescription listed above, and compute the matrix element of Magnetic Hamiltonian on square lattice for this state following (A.6). In this comparison it is easy to observe that the nontrivial delta functions in (16) are exactly same as those arising in evaluating the 6j symbols in (A.6). More importantly, our numerical calculation using random loop configuration reveals that the numerical value of the non-zero matrix elements for each and every cases matches exactly (upto a sign) with each other for the calculations done on square lattice and hexagonal lattice. The discrepancy in sign arises as for the particular convention of defining the loop states on square plaquette, that we have chosen, each and every term becomes positive. We repeat this comparison for 1000 random loop configurations and find this exact matching for each and every case. For the purpose of illustration, we only quote a few sample results in Table 2. Hence, this numerical study proves that the dynamics of loop states on a square plaquette is identical to that on an hexagonal plaquette as long as one is interested only in orthonormal loop states, which are actually relevant for exact physical degrees of freedom.

Point Splitting and virtual hexagonal lattice
In the last section, we have established the equivalence in dynamics of orthonormal loop configurations in a square and hexagonal lattice. In this section, we prescribe a virtual point splitting technique, which translates any square lattice to its hexagonal counter part. As a result of this transition, we gain a theory formulated in terms of only explicitly orthonormal loop degrees of freedom at each site, and pay the price of an extra Abelian Gauss law constraint. This price is actually negligible as the square lattice already had the Abelian constraints to solve and the extra one in hexagonal lattice is on very same footing as those.
The two dimensional hexagonal lattice, that we demonstrated here, is obtained by virtually splitting of each lattice site of a two dimensional square lattice. Elaborating a bit, consider a lattice site 'x' as shown in Fig. 8. From this site on a 2d lattice, clearly 4 links emerge in 1, 2,1,2 directions each carrying n i , (i = 1, 2,1,2) number of prepotentials (or angular momentum fluxes). We split this site 'x' into two sites 'x e & x o '. Links from direction 1 & 2 meet at site x e , where as links from direction1 &2 meet at x o . The splitted sites x 1 -x 2 are connected by a virtual link along direction 3. This same splitting done at each and every sites on the square lattice lifts the lattice to an hexagonal structure as shown in Fig. 12. Now, prepotential formulation on this hexagonal lattice yields a local loop formulation of lattice gauge theory, exactly eqiuivalent to the original square lattice, but contains only orthonormal and physical loop degrees of freedom as this is free from complicated Mandelstam constraints. This makes the analysis on hexagonal lattice simpler for practical purpose of analytical as well as numerical computation.
At this point we must explicitly match the degrees of freedom of these two systems as well. For the site x on the square lattice, there is six linking numbers or six loop degrees of freedom as earlier. Moreover, there was two Abelian Gauss law along two directions of the lattice (21) and one Mandelstam constraint (implying only non-intersecting loops contribute to physical degrees of freedom). Resulting only three phys- ical degrees of freedom. Now, coming back to hexagonal lattice, two sites, say x 1 & x 2 corresponds to actual site x on the square lattice and together should have only three degrees of freedom. Each site of the hexagonal lattice contains three linking numbers or three loop degrees of freedom, hence total six loop degrees of freedom matches with that of the square plaquette case. Unlike square plaquette, here there is no Mandelstam constraint at all to solve as there is no concept of intersecting loops. Together with the two Abelian Gauss law constraints along directions 1 and 2 (same as square plaquette), there exists one more Abelian Gauss law constraint on the link connecting two splitted site (along direction 3) and hence yielding exact degrees of freedom as square plaquette. Hence, counting of degrees of freedom goes as: where, S is the splitting index, which denotes each site has been splitted into S virtual sites; d is the dimension of the lattice, and d v is the number of virtual links connecting the splitted site, which is d v = S−1. It is straightforward to show that this analysis smoothly extends to any higher dimension as well.

The Hamiltonian and Average loop configurations
In this work, we consider pure SU(2) gauge theory defined on a square lattice. The system is described by the Kogut-Susskind Hamiltonian, The virtual point splitting, discussed in last section now enables us to define the theory on hexagonal lattice. On hexagonal lattice, the hexagonal plaquettes are surrounded by six links, out of which four are links of the original square lattice along directions 1 & 2, and remaining two are virtual links, resulted from point splitting along direction 3. The electric fields are defined at each end of the links of the original square lattice. Hence even for hexagonal lattice, only the electric fields for links along directions 1 & 2 contribute to the Hamiltonian given in (20). In prepotential formulation, this electric part of Hamiltonian H E counts the fluxes which are actually related to number of prepotentials sitting at each end of the links of the original lattice [5][6][7][8][9].
The magnetic part of the Hamiltonian H B is anyway more complicated to analyse. Clearly, at weak coupling regime, this part contributes most, and hence it is essential to simplify it as much as possible to make analytic calculations feasible. The part however is altered from the square lattice, as it must contain trace of the product of link operators along the full hexagonal plaquette to make the smallest closed loop.
The magnetic Hamiltonian on hexagonal lattice contains 2 6 different gauge invariant plaquette terms for prepotential formulation of SU(2) gauge theory on hexagonal lattice. Pictorially, these plaquette terms contain n d number of dotted link and 6 − n d number of solid links around a plaquette, for n d = 1, 2, .., 6. These solid and dashed line comes in different combination yielding all of the 64 plaquette terms as illustrated in Fig. 11.
Let us now, concentrate on the loop configurations of the system and proceed towards a mean field ansatz (i.e average loop configuration) for low energy limit of the theory.
We have already characterized the local loop states at a particular site x of a hexagonal lattice in terms of linking numbers as |l 12 , l 23 , l 31 x in (7) and also discussed how these linking numbers are related to the occupational number basis in (9). The original Wilson loops, which are non-local are obtained by weaving these local loops at neighbouring sites along the three directions, using an extra Abelian Gauss law [5][6][7][8][9] constraint at each link direction as given below:

Mean Field Ansatz
In this sub-section, we make an ansatz for the vacuum loop configuration of the SU(2) lattice gauge theory. Strong coupling vacuum of the system is well-known and consists of 0 flux state. Whereas, in the naive continuum limit, as g → 0, all the loop configurations contribute to the low energy spectrum. However, the maximum contribution is expected to come from large loops carrying large fluxes. In prepotential formulation, the size of the loop is not relevant, as all the loops has been made local. Abelian weaving along the links give rise to the standard Wilson loops. In this context, let us make a general ansatz for the vacuum loop configuration of the system irrespective of the coupling regime. Let us assume that the low energy loop configurations are given by same amount of fluxes flowing across each site. Note that, we are working on a virtual hexagonal lattice in two dimensions as discussed before. Note that, on these lattices, there exists two different types of flux, flowing across each site. The standard flux flowing along the direction of original links, i.e between {12} direction is the real flux. Moreover, there exist the fluxes, which flows from a real direction to a virtual direction, namely {23} and {31}. We consider this, and make an ansatz that, where, l i j (x)'s are the linking quantum number (measuring fluxes flowing across {i j} directions) specifying the loop states in (7) at a site x. The electric part of the Hamiltonian (20) in Prepotential formulation reads as: where, n 1 and n 2 counts the number of prepotentials residing on links along 1&2 directions respectively. This occupation numbers are related to the linking numbers at each site x and hence with mean values L&M as, Implying, Next we concentrate on the magnetic part of the Hamiltonian within this mean field approximation. The magnetic Hamiltonian contains 64 plaquette operators given in Fig. 11. Each of these 64 plaquette terms consists of 6 local loop operators given in Fig. 6, which comes with a coefficients listed in Table 1 which are functions of number operators (i.e linking  Fig. 10 The numerical value of the mean field variable L (M) for which the Hamiltonian function reaches a minima , i.e L| H min is plotted against the coupling g. The curve shows a clear jump in the average loop configuration, moving from strong to weak coupling regime of the theory numbers or occupation numbers). However, within the mean field ansatz they can be regarded as C-numbers (i.e functions of constants L and M) as listed in 3. Combining these coefficients we get 2 6 coefficients seating in front of each of the Magnetic term given in Fig. 9, which we call as C i , for i = 1 to 64. Now, in order to fix the numerical values of L and M, in our analysis, we calculate the Hamiltonian for the limiting case, when the all the plaquette operators (except the coefficients) becomes 1. Now, the Hamiltonian function for each plaquette reads as: We now minimize this function for different values of the coupling g. This minimization yields a set of L and M values at different couplings, for which the Hamiltonian function reaches a minima. We plot this set of L values at H min in Fig. 10 to get a clear notion of phase transition between the mean field phases at weak and strong coupling regime occurring exactly at g = 1. A very similar curve is obtained for M as well, which shows the exact same nature except the fact that for each values of g, the numerical value of M is much less than the corresponding L in weak coupling regime, i.e L >> M >> 0. This same result will be again obtained in a more precise calculation for the proposed reduced Hamiltonian in the next section. Finding the spectrum for the full Kogut-Susskind Hamiltonian even in the mean field approximation is still a challenge at this stage. In the next section we construct a reduced Hamiltonian which describes the dynamics within the Hilbert space of the above mentioned average loop configurations.

The reduced Hamiltonian and its spectrum
In the last section we have made a mean field ansatz, in which an arbitrary loop state defined at each site of the virtual hexagonal lattice is given by the average value for the linking number variables is given by: In this section, we will reduce the full Kogut-Susskind Hamiltonian to a sub Hamiltonian, which describes the dynamics within this mean field ansatz. With this reduced system, one can get reasonable physical results with minimal calculational effort and hence establishes this as a valid toy system to understand weak coupling regime of the gauge theory analytically, numerically as well as by quantum simulating the system.

The Sub-Hamiltonian
As we have already explained in the last section, the electric part of the Kogut-Susskind Hamiltonian in the mean field approximation reads as: Next we concentrate on the magnetic part of the Hamiltonian within this mean field approximation. The magnetic Hamiltonian contains 64 plaquette operators, each of which is a combination of six local vertex operators listed in Table 3 within the mean field ansatz. Let us now analyze each and every plaquette term illustrated in Fig. 11: (1) This term together with its Hermitian conjugate pair (2 terms) creates or annihilates flux around a full plaquette. (2) This term along with its rotationally symmetric (6 terms) and hermitian conjugate pairs (6 more terms) increases or decreases the length of the Wilson loops by 5 units. (3) Same as (2), total 12 terms, changes length by 4 units. (4), (5), (6) and the rotationally symmetric 9 terms ( 9 more hermitian conjugate terms) merge (separate) two loops and construct  7) and its hermitian conjugate term (total 2 terms) are rotationally symmetric and merges three loops to construct e bigger loop of their combined length. (8), (9) and their rotationally symmetric 6 terms for each constitute a hermitian conjugate set of 12 terms, each of which merges two loops and construct a bigger one of the same length. Finally (10) and its 6 rotationally symmetric terms are hermitian conjugate set and changes shape of a loop without changing its length. Let us now make the following observations: Among the 64 plaquette terms given in figure , there are certain terms which create or annihilate net fluxes around a plaquette. These plaquette terms can indeed build up the complete loop space, starting from strong coupling vacuum. However, we are interested in a particular state of the system, where each of the loops takes an average value throughout the lattice, hence there really exists no scope for creating or annihilating any net flux around any of the plaquettes. Hence, we choose only a subset of these 64 plaquette terms which keeps the dynamics within the mean field ansatz. More specifically, we choose a sub-magnetic Hamiltonian, which does not create or annihilate any net flux around a plaquette, but rearranges the internal loop configurations. Note that, (vii),(viii),(ix) and (x) type terms in Fig. 11 are loop operators that do not change any net flux flowing across the lattice, but just rearranges the intertwiners across different sites. These plaquette operators form a rotationally symmetric hermitian operatorThose are the terms which has equal number of creation and annihilation of links around a plaquette, pictorially represented by a plaquette with 3 solid link and 3 dotted link. These special terms are explicitly given as below: 1. Plaquette consisting of alternate solid and dashed line.
There are only two options for this, which constitutes the rotationally symmetric and Hermitian operator. We denote this operator by 2. Plaquette consisting of three consecutive solid and dashed lines. There are six options for this, which constitutes the rotationally symmetric and Hermitian operator. We denote this operator by 3. Plaquette consisting of two consecutive solid line, two consecutive dashed line and then a single solid and dashed line along with their Hermitian conjugate plaquette terms. We denote this type of operator by Each of these two types of plaquette terms has 6 rotationally symmetric contribution. Hence this particular type of rotational symmetric Hermitian operator contains total of 12 individual plaquette operators.
Combining the above three, we find that the magnetic part of the reduced Hamiltonian contains 2 + 6 + 12 = 20 plaquette operators and is completely rotationally symmetric and Hermitian.
Let us now construct a particular basis consisting of the state |L , M and the twenty plaquette operators listed above, acting on it. Hence we have a 21 dimensional basis vectors listed as below: where, the braces in the last three states denote the sets consisting of rotationally symmetric and Hermitian conjugate states. Each of these states consists of one of those three types of vertices at each site. We first calculate the local action of E 2 on each vertex, depending on whether it is a junction of two solid line (pp-vertex ≡ |pp ), two dashed line (mmvertex ≡ |mm ) and of one solid and one dashed line (pm vertex ≡ |pp ). Using the 12 vertices described in Table 3, the action of E 2 within mean field approach is given by: However, our basis now consists of 21 states given in (33), each having six of the above vertices. The action of E 2 on these states are given as below: At this point, we import another notation for the basis. Let the basis be denoted as |i , f ori = 0, .., 20, with |0 = |L , M , and |i = H i |0 for i = 1, 2, .., 20, where H i denotes the 20 plaquette operators described before (given as H pmpmpm , H pppmmm , H ppmmpm ). Hence, the reduced Hamiltonian takes the form: The electric part of the Hamiltonian acting on a state of this basis gives: and the magnetic part gives: In this particular basis, the Hamiltonian matrixĤ reduced for each plaquette takes the form given below: Note that, the diagonals have contribution from the electric term as well from the constant term in magnetic Hamiltonian. Other than the diagonal entries, other non zero matrix elements of the Hamiltonian lies along the first row and first column of the matrix, and those are basically given by the coefficients sitting in front of the 20 plaquette terms discussed above. We label these coefficients as C i , for i = 1, .., 20. The advantage of working with such a basis is that the Hamiltonian matrix takes a very special form, namely form of an arrowhead matrix which in turn enables us to solve for the eigenvalues analytically. Let us assume the Hamiltonian matrix to satisfy the following eigenvalue equation: This matrix equation is equivalent to the set of eigenvalue equations mentioned below: Now, from (42), Putting (43) back in (41), we get, where,λ = 1 − g 2 λ + g 4 E i . Solving (44) would yield theλ, which in tern gives all the eigenvalues of the 21×21 Hamilto-nian matrix. To simplify (44), we further make the following observation. Note that, the diagonal elements g 2 E i + 1/g 2 are highly degenerate, as one finds in (36) that, E i 's can take only three possible values, namely 6v, (2v + 2v p + 2v m ) and (4v + v p + v m ). Hence, the sum in (44) reduces to the following sum of three terms: where,C 1 ,C 2 ,C 3 are combinations of the coefficients C i taking into account of the degeneracy. (45) is a transcendental equation and can be solved graphically, however approaching a rough solution is quite easy. We plot the right hand side of this equation and get divergences at three poles precisely at λ = 0, g 4 , 2g 4 respectively. The plot of left hand side, i.e straight line y = x cuts the right hand side curve very close to the position of the poles. The exact value ofλ yields the value of the eigenvalues as: Note that, the smallest eigenvalue λ, actually corresponds to the highestλ from (46). To get the exact numerical value of the solution one needs to consider the coefficientsC i 's appearing in right hand side of (45). Note that, the coeffi-cientsC i 's are functions of the mean fields L&M. Hence, it is necassary to calculate the exact value of these mean fields.
For this purpose, we take the following approach.
-For some fixed g, we numerically calculate (using Mathematica) the value of L and M, for which lowest eigenvalue of the Hamiltonian matrix reaches a minima. In the strong coupling regime, i.e for g ≥ 1, that minima is always at L = M = 0. However, for smaller and smaller values of g, the minima is at larger and larger values of L&M, as listed in Table 4. This calculation establishes the naive analysis done in Sect. 4 of this paper which shows the existence of two different mean field phases of the system in strong and weak coupling regimes. Note that, we are interested in the weak coupling regime of the theory as the continuum limit lies there. Up to this point of this work, we have not used any assumption for weak coupling limit, except taking a mean field ansatz. We now fix the mean field configuration in a way, such that we are in the weak coupling regime of the theory. From the variational study discussed above, we see, that the ground state energy shows a perfect first order phase transition at coupling g = 1, above which the vacuum is the strong coupling vacuum, which in the mean field ansatz gives |L = 0, M = 0 . However, as g → 0, the lowest energy mean field state turns out to be a state comprising of large average flux at each site, implying the weak coupling vacuum to be consisted of loops carrying large fluxes, throughout the lattice. As g → 0, L >> M >> 0. The coefficients, listed in Table 3 also changes to particular limiting values as g → 0.
-We now consider a particlular value of g, in the weak coupling regime, and the corresponding mean field configuration. For each configuration, we exactly diagonalize the Hamiltonian matrix, and calculate the eigenvalues. We list the spacings between the lowest one and first excited one as Δλ 1 and similarly between the first and second excited one as Δλ 2 . We list these two gaps and their ratios at different values of the coupling constant g → 0 in Table 5. It is very much clear from Table 5, that the gaps are scaling as ∼ g 2 and the ratio of lowest two energy gaps converges to the numerical value of 1. Note that, the absolute value of masses are subjected to be renormalized. However, the ratio of the two consecutive mass-gaps are always physical. In our study we have shown it to converge to numerical value of 1, for g → 0 for arbitrarily large lattices, as all the computation was done locally at each site.
The scaling of mass-gaps ∼ g 2 is consistent with both the weak coupling perturbation expansion in [30] 4 , variational analysis [31,32] 5 , using cluster algorithm [34] 6 for SU(2) gauge theory on 2 + 1 dimensional lattice as well as with Monte Carlo study on 3dimensional Euclidean lattice [33]. 7 On a further note, the analytically calculated value of the ratio of energy differences for the reduced Hamiltonian within the mean field ansatz, matches with previous calculations at weak coupling limit for finite lattices. For example, the mass spectrum of 0 ++ sector for SU(2) lattice gauge theory in 2 + 1 dimensions is obtained in [31]  This result is obtained with a lattice consisting of 25 plaquettes and in the weak coupling region 1.8 < 1/g 2 < 3.6. Note that, the results obtained in [31] are consistent with the results obtained with Monte Carlo calculations within Euclidean formalism as well.
(47) establishes the validity of our much simplified system to extract out weak coupling results for practical purposes, which can be exploited to address various problems in future researches.

Summary and future directions
In this work we have proposed and justified an effective mean field description for the low energy spectrum of SU(2) lattice gauge theory in 2 + 1 dimension and have analytically calculated the spectrum at the weak coupling regime. Starting from prepotential formulation on the spatial 2D square lattice, we perform virtual splitting of each lattice site into two and end up with a virtual hexagonal lattice. On this hexagonal lattice, all of the local loop states in prepotential formulations constitutes an exact and orthonormal loop basis, with no further Mandelstam constraints. We have proposed a mean value ansatz for the loop configurations throughout the lattice contributing to the low energy spectrum of the theory. We have shown that such average loop configurations have two distinct phases at the strong and weak coupling regime. Next, we have chosen a reduced Hamiltonian, from the full Kogut-Susskind Hamiltonian, which keeps the dynamics of the loops confined into our ansatz. Variational study shows that this reduced system with mean value-loop configuration shows a clear jump between the weak and strong coupling vacuum. As we are interested to explore weak coupling regime of the theory, we choose the relevant average loop configuration in that regime and calculate the spectrum for the reduced Hamiltonian we choose. In this spectrum we find ΔE ∼ g 2 which is the expected weak coupling behaviour for mass gap of the theory. The spacings of the spectrum obtained in this work is as well consistent with the available literature at weak coupling regime of 2 + 1 dimensional SU(2) lattice gauge theory. We have discussed in detail, how the point splitting lattices are constructed in higher dimensions, which can as well be exploited to extend this work beyond 2 + 1 dimensions in a straight forward way.
In a recent and parallel work [27], the point splitting lattice is constructed and utilized to analytically study the weak coupling limit of SU(2) lattice gauge theory in 2 + 1 dimension as well. In that work, they have used the path integral representation of the phase space to analytically compute the dispersion relation at the lowest order in weak coupling perturbation expansion.
However, the particular study demonstrated in this paper shows that the physical results at the weak coupling regime of SU(2) gauge theory can be extracted from a much simpler mean field approximation made within prepotential formulation of the theory. Being completely gauge invariant, and formulated only in terms of relevant physical degrees of freedom, this technique is suited for both analytic calculations and numerical simulations. From analytic perspective, this study gives a clear notion of the weak coupling vacuum for pure gauge theory and its dynamics. From numerical perspective, this particular formulation is most suited for quantum Monte-Carlo simulation of Hamiltonian lattice gauge theory using a complete gauge invariant basis characterized by only integers. Till date, this aspect has not been studied extensively, but worth investigating in near future. That study will lead to explore some of the very important physics such as calculation of the entanglement entropy of lattice gauge theory. Last but not the least, there is a tremendous progress going on, in the recently developed research interests for quantum simulating gauge theories using both analog [28,29,[35][36][37][38] and digital quantum simulators [39][40][41][42][43]. The prepotential formulation has already been explored to propose quantum simulator for gauge theories [28,29] and is a promising framework to define a whole new prepotential paradigm in quantum simulating QCD [43]. This present work, shows the way to construct quantum simuator to simulate the loop dynamics of non-Abelian lattice gauge theory beyond strong coupling limit. The work in this direction is in progress and will be reported shortly.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: The results presented in this work are obtained by semi-analytic approach, where no explicit numerical simulation has been performed. Numerical results are obtained by performing algebraic calculations on Mathematica. The graphs presented in this work are obtained from explicit mathematical expressions contained in this paper and hence no explicit data set are provided.] Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .

Appendix A: Orthonormal loop basis and its dynamics on square lattice
Within prepotential framework of the theory one is equipped with a set of orthonormal basis states, defined at each lattice site [5][6][7][8][9]22,23]. For SU(2) lattice gauge theory on 2 dimensional spatial lattice, such an orthonormal basis can be easily obtained in terms of angular momentum fluxes. In this section, we briefly quote the result of plaquette (the smallest closed loop) dynamics from [5][6][7][8][9] which we have used to compare with the dynamics of hexagonal plaquette.
At each site of a two dimensional square lattice, four links are attatched in directions 1, 2,1 &2 carrying the angular momentum fluxes j 1 , j 2 , j1 & j2 respectively ( j i = n i /2). The gauge invariant state at a site x must have the net angular momentum, i.e the sum of the four angular momentums along the four directions zero. One can add these four angular momentum according to the following scheme: | j 1 , j 2 , j 12 , j1, j2, j12 = j 12 , j 1212 = 0 x (A.1) As, the abelian Gauss law implies j1| x = j 1 | x−e 1 and j2| x = j 2 | x−e 2 following (21). Hence, each site, one choice of such orthonormal basis at each site x is [5][6][7][8][9] | j 1 , j 2 , j 12 x (A.2) Fig. 12 Dynamics of loops on a square plaquette: orthonormal loop states are actually around the hexagon which fluctuate in plaquette dynamics. The hexagon arises as point splitting at each vertices, i.e a → a 1 ,