understanding the deconfining phase transition of qcd is a long - standing problem in high energy physics .
the situation is somewhat clearer in pure su(n ) yang - mills theory , where it is understood that on length scales larger than @xmath2 the relevant effective degrees of freedom are su(n ) spin variables ( thermal wilson loops ) in three dimensions , obtained by compactifying the euclidean time direction . in this picture
, the deconfining phase transition manifests as spontaneous breaking of the z(n ) center symmetry @xcite .
numerical simulations show that the interaction measure , @xmath3 , times @xmath4 , is approximately flat from @xmath5 to @xmath6 , for up to six colors @xcite .
similar behavior is observed in three space - time dimensions , for @xmath7 times @xmath8 @xcite .
this shows that the leading corrections to an `` ideal '' gas term , @xmath9 in four dimensions , and @xmath10 in three dimensions , are terms @xmath11 , in both four and three dimensions .
this differs from corrections due to a bag pressure as in the mit bag model , which would be independent of temperature .
this behavior can be reproduced within matrix models of the thermal wilson loop  @xcite with coupling constants which are rather simple functions of @xmath12 .    to understand this phenomenon better , it is useful to study the dynamics of the eigenvalues of the wilson loops in different representations of the gauge group .
the purpose of this work is to study the effective potential of bare wilson loops and of their eigenvalues , in the fundamental and adjoint representations of su(2 ) , through ab - initio simulations in discretized space - time .
our goal is to shed some light on the validity of the perturbative ansatzes for the effective action , which are used in the construction of effective theories .
effective potentials for variables related to the phases of the eigenvalues of thermal wilson loops have also been obtained from functional methods such as the functional renormalization group or dyson - schwinger equations typically using background landau - dewitt gauges @xcite , and our results may be useful to benchmark these non - perturbative continuum quantum field - theoretic computations . also , we wish to verify that the ansatz used to parametrize the results obtained in refs .
@xcite for the classical potential of a three dimensional matrix model of wilson loops in the fundamental representation is valid also for the full gauge theory .
we simulate pure su(2 ) gauge theory using the standard wilson action s= _
( 1- re tr  * u * _ )  , [ slatt ] where the sum runs over all plaquettes of the four dimensional lattice .
lattice configurations are generated using an exact heatbath algorithm @xcite which generates link variables according to the local boltzmann distribution given by the staple matrix .
su(2 ) matrices are stored in the quaternionic representation u = a_0 * 1 * + i _ j=1 ^ 3 a_j _ j , _
j=0 ^ 3 a_i^2=1  , [ eq : quater1 ] where @xmath13 and @xmath14 are the pauli matrices .
we construct thermal wilson loops @xmath15 in the fundamental representation in the usual way by multiplying temporal links : * l*(*x*)^f=_t=0^n_t-1u_4(*x*,t )  .[eq : ploop1 ] loops in higher representations can be generated from the fundamental loop @xcite .
we currently only also consider the adjoint representation which is constructed as * l*_ij^a=2tr(t_i * l*^f t_j * l*^f )  , [ eq : adj1 ] where @xmath16 are the three generators in the fundamental representation .
we are interested in the per - site effective potential @xmath17 for the polyakov loop @xmath18 and for the eigenvalues of the wilson line @xmath19 , in the fundamental or adjoint representations . in other words
, @xmath20 describes fluctuations about the volume averaged mean field @xmath21 .
the latter , in turn , determines the couplings @xmath22 , @xmath23 , @xmath24 in @xmath20 , see below .    in our work ,
`` sites '' are understood as the lattice - points of the three dimensional spin - system of the thermal wilson loops .
our approach differs from the computation of a potential for volume averaged fields as in refs .
.    we also comment on the difference of our approach to that of ref .
@xcite . in this work
they take the wilson line to be constant in time .
this is allowed , and in the continuum , corresponds to a static gauge , @xmath25 .
there is still a residual gauge freedom to fix spatial gauge transformations at a given time .
for example , one could fix to coulomb gauge , @xmath26 . instead , ref .
@xcite computes an effective potential by assuming that the thermal wilson line is constant in space . in the continuum this corresponds to setting @xmath27 for all three spatial directions , and
is not an allowable gauge condition .
because of this , their results are very different from ours or from refs .  @xcite .
in particular , @xcite finds a potential which is completely flat in the confined phase .    to obtain @xmath17
we first compute the per - site probability distribution @xmath28 by histogramming each observable . from these distributions
we obtain its classical or `` constraint '' effective potential @xmath29 .
the effective potential is then obtained from the legendre transform of the moment - generating function @xmath30 : w(h)= dx(-v_0(x)+hx )  .
[ eq : momgen ] here , @xmath31 stands for the respective observable and the integral runs over the entire range of @xmath31 .
we obtain v_()= - ( )  , where ( ) = w(h())-h ( )  ,   =  .
[ eq : legendre1 ] in most cases considered here it is simpler to carry out the legendre transformation via ( ) = _ h ( h - w(h ) )  .[eq : legendre2 ] this is our method of choice .
, for each data - set we take a proper parametrization of @xmath32 as input . for a given @xmath33 we change @xmath34 from @xmath35 to @xmath36 in steps of @xmath37 . for each @xmath34 , @xmath30 is computed from eq .
( [ eq : momgen ] ) using a @xmath38 point gaussian quadrature with @xmath39 .
we obtain the maximum @xmath40 with respect to @xmath34 .
we repeat this process for @xmath41 equidistant values of @xmath33 in the domain of @xmath42 .
we find that this procedure yields results which are stable under further increase of the resolution .
] we use eq .
( [ eq : legendre1 ] ) only in some special cases where analytical approximations to @xmath30 can be easily obtained .
we use the string tension to fix the lattice spacing and thus the temperature @xmath43 in physical units .
the exact procedure is discussed at length in appendix [ sec : physunits ] ( in all figures , the quoted @xmath8 are understood to have a relative uncertainty of less than @xmath44 ) .
we change the temperature by varying @xmath45 .
the advantage of this fixed scale approach is that it allows us to clearly disentangle the effects of temperature changes and renormalization . in the following
we present a few formulae which shall be used throughout this work . using eqs .
( [ eq : quater1 ] ) and ( [ eq : adj1 ] ) with @xmath46 one can easily obtain the compact expression * l*_ij^a=2(a_i a_j + a_0 _ k=1 ^ 3 _ ijk a_k + _ ij(a_0 ^ 2- ) )  , [ eq : funadjoint_entries ] which relates the elements of the adjoint matrices to the quaternionic parameters @xmath47 . from eq .
( [ eq : funadjoint_entries ] ) it immediately follows that tr*l*^a=|tr*l*^f|^2 - 1  .[eq : funadjoint ] from ( [ eq : quater1 ] ) it follows that the eigenvalues @xmath48 of the fundamental loop are given by _ 1,2=a_0  .
the eigenvalues form a pair of complex conjugates which lie on the unit circle .
also , they are related to the polyakov loop @xmath49 by ^f = tr*l*^f = a_0 = ( _ 1 + _ 2 )  . likewise , for the adjoint wilson loop ( which is a real 3-by-3 matrix ) one can easily show that the eigenvalues @xmath50 are _ 1=1  ,  _
2,3=(tr*l*^a -1 )  .
the non - unit eigenvalues @xmath51 again form a pair of complex conjugates with @xmath52 .
thus for both representations a single phase @xmath53 uniquely fixes all eigenvalues .
the @xmath19 are related by a similarity transform to @xmath54 and @xmath55 respectively , so that @xmath56 also fixes the trace of @xmath19 through the relations tr*l*^f= 2 ( ^f )  , tr*l*^a= 1 + 2 ( ^a )  .[eq : wilsontrafo ]
we simulate @xmath57 gauge theory at @xmath58 and @xmath59 , which corresponds to @xmath0 at @xmath60 .
the lattice spacing is @xmath61 and @xmath62 , respectively .
we use time - like lattice sizes of @xmath63 . for every lattice size and @xmath64 value
several hundreds of independent lattice configurations were sampled from the equilibrium distribution .
bin sizes for the histograms were chosen such that the results appear as smooth lines at the resolutions of the figures presented here .
all distributions @xmath28 presented here are normalized .
there is a subtle issue regarding ergodicity of simulations in the deconfined phase where the @xmath65 center symmetry is spontaneously broken : in principle , at any finite lattice volume , an ergodic algorithm which runs forever will tunnel infinitely often between the two groundstates .
hence , all probability densities @xmath28 will display a mixing of these groundstates and show no signs of a broken symmetry .
the physics presented here however is concerned with the thermodynamic limit , for which the potential barrier between the groundstates becomes infinitely large .
tunneling events should therefore not occur , which strictly speaking violates ergodicity . from a practical point of view
it is justified to ignore this problem as long as the lattice volumes considered are sufficiently large so that tunneling is strongly suppressed .
one should take care however to ensure that this is truly the case for a given temperature , since such tunneling events occur more frequently as one approaches @xmath0 from above .
one can achieve this , for instance , by monitoring the volume average of the polyakov loop and looking for sign flips .
the results presented here are unaffected by this issue .
the lowest temperatures in the deconfined phase which we consider here are well above @xmath0 , so that even on a @xmath66 lattice tunneling is practically impossible and the infinite volume limit is well approximated . for each dataset discussed in this section simulations
were done on lattices with @xmath67 and @xmath68 .
we find that both cases give identical results ( indistinguishable histograms ) , which confirms that finite - volume effects play no role .
likewise , this confirms that statistical errors are negligible given the resolutions of the presented figures . without restricting generality we chose to represent the symmetry broken phase as the state with @xmath69 .
we begin with the probability density @xmath70 of the polyakov loop @xmath71 in the fundamental representation .
[ fig : su2fun1 ] shows @xmath70 for @xmath58 and @xmath59 at @xmath63 , which corresponds to temperatures ranging from slightly below @xmath0 to roughly @xmath72 ( smaller @xmath64 corresponds to lower @xmath12 for fixed @xmath45 ) .
( top ) and @xmath59 ( bottom ) .
the different temperatures correspond to @xmath63 respectively .
[ fig : su2fun1],title="fig : " ]   ( top ) and @xmath59 ( bottom ) .
the different temperatures correspond to @xmath63 respectively .
[ fig : su2fun1],title="fig : " ]    one can clearly see the effect of spontaneous breaking of the @xmath65 ( @xmath73 ) symmetry : right at @xmath0 the distributions are symmetric around @xmath74 .
in fact , we have confirmed that they are identical for any @xmath75 $ ] at the respective @xmath76 . at higher temperature @xmath70
is skewed and it develops a peak at non - zero @xmath49 which moves towards @xmath77 as one raises the temperature . in fig .
[ fig : su2fun2 ] one can see that the corresponding constraint effective potential @xmath78 gets tilted , which is exactly what one expects at a second order phase transition .
( top ) and @xmath59 ( bottom ) .
the different temperatures correspond to @xmath63 respectively .
[ fig : su2fun2],title="fig : " ]   ( top ) and @xmath59 ( bottom ) .
the different temperatures correspond to @xmath63 respectively .
[ fig : su2fun2],title="fig : " ]    what is striking is that the lattice spacing affects the distribution . comparing the simulations at the two different values of the coupling , @xmath58 ( corresponding to a larger lattice spacing ) clearly approaches the perturbative vacuum ( @xmath79 ) much faster with increasing temperature .
it is well known that the expectation value of the bare polyakov loop is scale dependent and vanishes in the continuum limit ( hence the need to construct renormalized loops @xcite ) . however , here we show that in fact the entire distribution is affected in a profound way . in the continuum @xmath80 limit
@xmath28 is deformed into the distribution at @xmath81 .
one way of thinking of this is that the @xmath81 distribution is an attractive ultra - violet fixed point of renormalization group transformations .
what is also noteworthy is that the distribution appears not to change at @xmath82 . which we are not sensitive to .
such would be consistent with the findings of other authors @xcite . ]    to make these statements quantitative we model the distributions using appropriate parametrizations .
it turns out that the @xmath81 distribution is exactly p^(t_c)()=  , [ eq : vand1 ] which corresponds to a random walk on the su(2 ) group manifold with the appropriate measure .
the corresponding constraint potential is the vandermonde potential @xcite v^(t_c)_0()=-(1-^2 ) - ( )  .[eq : vand2 ] at @xmath83 , for the entire range of temperatures considered here , the following ansatz reproduces the simulated curves accurately : v_0()=v^(t_c)_0 ( ) + a(t)-b(t)+c(t)^2  .[eq : vand3 ] this @xmath84 corresponds to the potential for a single `` local spin '' @xmath85 given that @xmath86 .
if instead one adopts the convention that @xmath87 by performing a z(2 ) transformation then @xmath84 is obtained from the above by letting @xmath88 .
the potential  ( [ eq : vand3 ] ) corresponds to the distribution p()=(-a(t)+b(t)-c(t)^2 )  .
[ eq : vand4 ] we obtain the parameters @xmath89 from a @xmath90 fit .
we do not plot the fits as they are indistinguishable from the simulation results .
the temperature dependence of the parameters for both values of @xmath64 are shown in fig .
[ fig : su2evfun5 ] .
( top ) and @xmath59 ( bottom ) . @xmath22 and @xmath24 are scaled up by factors of @xmath91 .
the connecting lines are there to guide the eye only but do not represent actual models .
[ fig : su2evfun5],title="fig : " ]   ( top ) and @xmath59 ( bottom ) . @xmath22 and @xmath24 are scaled up by factors of @xmath91 .
the connecting lines are there to guide the eye only but do not represent actual models .
[ fig : su2evfun5],title="fig : " ]    it appears that the term linear in @xmath49 is dominant for the temperature range investigated here .
it is also precisely this term which drives breaking of the @xmath65 symmetry ( it in some sense describes the interaction of a particular site with the external field generated by the symmetry broken state of the lattice as a whole ) .
the ansatz ( [ eq : vand3 ] ) is exactly of the form which was used in refs .
@xcite to parametrize the constraint effective potential for @xmath49 in a @xmath92 effective matrix theory of @xmath57 wilson lines , with a similar behavior of the parameters @xmath89 .
( top ) and @xmath59 ( bottom ) .
the quoted temperatures correspond to @xmath63 respectively .
[ fig : su2fun2eff],title="fig : " ]   ( top ) and @xmath59 ( bottom ) .
the quoted temperatures correspond to @xmath63 respectively .
[ fig : su2fun2eff],title="fig : " ]    we obtain the effective potential @xmath20 by numerically solving eq .
( [ eq : legendre2 ] ) , using the parametrization  ( [ eq : vand3 ] ) of @xmath84 as input .
obtaining the supremum numerically is straightforward .
the integral in eq .
( [ eq : momgen ] ) is computed numerically via gaussian quadrature .
the results are shown in fig .
[ fig : su2fun2eff ] .
we have confirmed that the expectation value = _
-1 ^ 1d p ( )  , coincides with the point at which @xmath20 is minimal , as it should .
( top ) and @xmath59 ( bottom ) .
the quoted temperatures correspond to @xmath63 respectively .
[ fig : su2evfun3],title="fig : " ]   ( top ) and @xmath59 ( bottom ) .
the quoted temperatures correspond to @xmath63 respectively .
[ fig : su2evfun3],title="fig : " ]     ( top ) and @xmath59 ( bottom ) .
the quoted temperatures correspond to @xmath63 respectively .
[ fig : su2evfun4],title="fig : " ]   ( top ) and @xmath59 ( bottom ) .
the quoted temperatures correspond to @xmath63 respectively .
[ fig : su2evfun4],title="fig : " ]    from the same ensembles we obtain the per - site distribution of the phase @xmath56 of the eigenvalues of the wilson loop @xmath19 .
[ fig : su2evfun3 ] and [ fig : su2evfun4 ] show the distributions @xmath93 and the potentials @xmath94 respectively .
one can see that right at @xmath0 the distribution is peaked around @xmath95 , which is the state of maximal eigenvalue repulsion .
increasing the temperature drives the peak towards @xmath96 though eigenvalue repulsion persists even above @xmath0 .
such an effect is known to occur also in matrix models of the wilson line @xcite .
we can transform the parametrization of @xmath70 given in eq .
( [ eq : vand4 ] ) to @xmath93 .
( [ eq : wilsontrafo ] ) provides a one - to - one mapping from @xmath97 $ ] to @xmath98 $ ] ; we can thus use the standard integral transformation law to obtain @xmath93 from our models ( [ eq : vand1 ] ) and ( [ eq : vand4 ] ) .
@xmath99 $ ] maps into the same distribution @xmath70 , thus the restriction to positive @xmath56 is legitimate . using _
-1 ^ 1dp ( ) = _ ( = -1)^(=1 ) p ( ( ) ) d=_^0 p()d  , we arrive at p^(t_c)()=^2= ( 1+(2 + ) )  , [ eq : vand5 ] for the distribution right at @xmath0 and p ( ) & = & ^2 [ eq : vand6 ] + & & ( -a(t)+b(t)()-c(t)^2 ( ) )  , for the general case .
the parameters @xmath89 are unchanged .
( [ eq : vand5 ] ) and ( [ eq : vand6 ] ) again describe the simulated results to such accuracy that the lines are indistinguishable and thus the fit - curves are left out of the figures .
we again use our model function v_0()= & & [ eq : vand7 ] + & & -1.2 cm -(^2 ) + a(t)-b(t)()+c(t)^2 ( )  , to obtain @xmath100 . the legendre transformation ( [ eq : legendre2 ] ) is carried out numerically , as before .
the results are shown in fig .
[ fig : su2evfun4eff ] .
( top ) and @xmath59 ( bottom ) .
the different temperatures correspond to @xmath63 respectively .
[ fig : su2evfun4eff],title="fig : " ]   ( top ) and @xmath59 ( bottom ) .
the different temperatures correspond to @xmath63 respectively .
[ fig : su2evfun4eff],title="fig : " ]    in order to make our results portable ( to compare them with other calculations or use them as input for effective theories ) we would like to obtain analytical expressions for @xmath100 and @xmath20 .
one possible approach is to again use polynomial models ( together with the vandermonde contribution ) and adjust the parameters to fit the data .
it turns out , however , that this approach is not satisfactory as it requires high - order polynomials and many parameters in both cases .
given how well our models ( [ eq : vand4 ] ) and ( [ eq : vand6 ] ) described the simulated results for @xmath70 and @xmath93 , it would therefore be nice if the legendre transformation ( [ eq : legendre1 ] ) could be carried out exactly .
unfortunately this is not possible in closed form .
nevertheless , an approximate solution can still be obtained using the saddle point approximation .
the next - to - leading order result is easily obtained , and we will compare this result to our data in the following .
consider a stationary point @xmath101 of the exponential appearing in eq .
( [ eq : momgen ] ) such that v_0(x_0)=h  .[eq : saddle1 ] we can expand the exponent in eq .
( [ eq : momgen ] ) around this point . neglecting terms containing derivatives of higher than second order we can carry out the integral : w(h ) & = & _ -^dx ( -v_0(x)+hx ) + & & -1.4 cm _ -^ dx ( -v_0(x_0)+hx_0 -v_0(x_0)(x_0-x)^2 ) + & = & -v_0(x_0)+hx_0 + ( )  .
it holds that = x_0 + o ( v_0 )  , and thus the nlo effective potential is given by v_(x)v_0(x)+ ( )  .
the leading order is just the classical potential .
this formula can be applied to the potentials for both @xmath49 and @xmath56 . from eqs .
( [ eq : vand7 ] ) and  ( [ eq : vand3 ] ) we obtain the second derivatives v_0()= 2c(t)+  , and v_0()= + b(t)()+2 c(t)(^2()-^2 ( ) )  . in figs .
[ fig : su2fun2eff_spcompare ] and [ fig : su2evfun4eff_spcompare ] we compare the nlo saddle point approximation to the exact results for both @xmath20 and @xmath100 .
the approximate curves are corrected by an additive shift to account for the loss of normalization of @xmath28 which occurs .
the shift is such that @xmath17 is moved to zero at its minimum . in both cases @xmath102
is chosen as an example .
it appears that the approximation is quite good .
( top ) and @xmath59 ( bottom ) .
@xmath102 simulation result compared to nlo saddle point expansion .
the sp curve has been shifted vertically to bring the minimum to zero . [
fig : su2fun2eff_spcompare],title="fig : " ]   ( top ) and @xmath59 ( bottom ) .
@xmath102 simulation result compared to nlo saddle point expansion .
the sp curve has been shifted vertically to bring the minimum to zero .
[ fig : su2fun2eff_spcompare],title="fig : " ]     ( top ) and @xmath59 ( bottom ) .
@xmath102 result compared to nlo saddle point expansion .
the sp curve has been shifted vertically to bring the minimum to zero .
[ fig : su2evfun4eff_spcompare],title="fig : " ]   ( top ) and @xmath59 ( bottom ) .
@xmath102 result compared to nlo saddle point expansion .
the sp curve has been shifted vertically to bring the minimum to zero .
[ fig : su2evfun4eff_spcompare],title="fig : " ]
lastly , we investigate how well the classical potentials @xmath84 and @xmath94 can be approximated by taylor polynomials around the confined ( @xmath103 ) or deconfined ( @xmath104 ) vacua .
expanding about @xmath105 we obtain v_0,1 ( ) & & - ( ) + a(t ) + b(t ) ( - ) + & & + ( c(t)+1)(- ) ^2 - ( - ) ^3 + & & + ( - ) ( - ) ^4  .[eq : vand8 ] likewise , about @xmath106 v_0,2 ( ) & & -(^2 ( ) ) + a(t ) -b(t ) + c(t ) + & & + ( -c(t ) ) ^2 + & & + ( - ) ^4  .[eq : vand9 ] the logarithmic term has not been expanded since it diverges at @xmath106 .    finally , an expansion around @xmath107 yields v_0,3 ( ) & & -(^2 ( ) ) + a(t ) + b(t ) + c(t ) + & & + ( -c(t ) ) ( -)^2 + & & + ( - ) ( -)^4  .[eq : vand10 ] note that  ( [ eq : vand10 ] ) is nothing but the z(2 ) transform of  ( [ eq : vand9 ] ) corresponding to @xmath108 and @xmath88 .
[ fig : su2evfun_taylor ] compares these expansions to the simulation results for @xmath102 , which corresponds to the highest temperature for each respective @xmath64 .
it can be seen that eq .
( [ eq : vand8 ] ) works well , even if one terminates the expansion at @xmath109 .
the approximation smoothly approaches the simulation results when at higher orders , and @xmath110 reproduces the exact result to high accuracy . convergence is slower for eq .
( [ eq : vand9 ] ) and very slow for ( [ eq : vand10 ] ) .
the latter is not extremely suprising since @xmath107 is the wrong groundstate .
( top ) and @xmath59 ( bottom ) .
@xmath102 results compared to taylor expansions .
[ fig : su2evfun_taylor],title="fig : " ]   ( top ) and @xmath59 ( bottom ) .
@xmath102 results compared to taylor expansions .
[ fig : su2evfun_taylor],title="fig : " ]    a similar expansion can be constructed for the logarithmic contribution to @xmath84 .
expanding eq .
( [ eq : vand3 ] ) to @xmath111 about @xmath74 we get v_0,1 ( ) & & - ( ) + a(t)- b(t ) + & & + ( c(t)+)^2 + ^4+^6  .[eq : vand11 ] fig .
[ fig : su2fun_taylor ] shows the convergence of this expansion towards the simulated curves for @xmath102 .
high orders are required to accurately reproduce the data .
( top ) and @xmath59 ( bottom ) .
@xmath102 results compared to taylor expansions .
[ fig : su2fun_taylor],title="fig : " ]    ( top ) and @xmath59 ( bottom ) .
@xmath102 results compared to taylor expansions .
[ fig : su2fun_taylor],title="fig : " ]      from the same ensembles , we can also construct distributions @xmath112 ( here without any normalization factors before the trace ) for the adjoint loop as well as for its eigenvalues .
the results are shown in figs .
[ fig : su2adj1 ] and [ fig : su2evadj1 ] , respectively . figs .
[ fig : su2adj2 ] and [ fig : su2evadj2 ] show the corresponding potentials @xmath29 .
what is striking is that temperature effects appear to be small .
nevertheless , a slight suppression of @xmath93 around @xmath107 can be observed with increasing temperature .
( top ) and @xmath59 ( bottom ) .
the different temperatures correspond to @xmath63 respectively .
[ fig : su2adj1],title="fig : " ]   ( top ) and @xmath59 ( bottom ) .
the different temperatures correspond to @xmath63 respectively .
[ fig : su2adj1],title="fig : " ]     ( top ) and @xmath59 ( bottom ) .
the different temperatures correspond to @xmath63 respectively .
[ fig : su2adj2],title="fig : " ]   ( top ) and @xmath59 ( bottom ) .
the different temperatures correspond to @xmath63 respectively .
[ fig : su2adj2],title="fig : " ]     ( top ) and @xmath59 ( bottom ) .
the different temperatures correspond to @xmath63 respectively .
[ fig : su2evadj1],title="fig : " ]   ( top ) and @xmath59 ( bottom ) .
the different temperatures correspond to @xmath63 respectively .
[ fig : su2evadj1],title="fig : " ]     ( top ) and @xmath59 ( bottom ) .
the different temperatures correspond to @xmath63 respectively .
[ fig : su2evadj2],title="fig : " ]   ( top ) and @xmath59 ( bottom ) .
the different temperatures correspond to @xmath63 respectively .
[ fig : su2evadj2],title="fig : " ]    we parametrize these results in a way similar to those for the fundamental representation .
we begin by transforming the parametrization of @xmath93 given in eq .
( [ eq : vand6 ] ) to the adjoint representation .
once this is achieved , @xmath113 can be easily obtained . to approach this problem we require a one - to - one mapping @xmath114 from the fundamental to the adjoint phase .
to construct such a mapping a few subtleties arise .
( [ eq : funadjoint ] ) implies that ^2(^f ) = ( 1+(^a ) )  .[eq : funadjointphase ] completing one full rotation in the fundamental phase @xmath115 ( from @xmath116 to @xmath117 ) generates two full cycles of the square cosine on the left - hand side of eq .
( [ eq : funadjointphase ] ) , whereas one rotation of the adjoint phase @xmath118 is only one cycle of the right - hand side .
both sides are symmetric around @xmath106 , thus it is sufficient to consider only @xmath119 $ ] which can be mapped to one full rotation of the adjoint phase .
second , the left - hand side possesses an inflection point at @xmath120 .
thus the mapping must be constructed piece - wise in the intervals @xmath121 $ ] and @xmath122 $ ] .
taking positive roots in eq .
( [ eq : funadjointphase ] ) and considering that @xmath123\\         -\cos\left(\frac{x}{2}\right ) & : x \in [ \pi,2\pi ]       \end{array }     \right.\ ] ] we obtain @xmath124\\         -\cos\left(\frac{\phi^a}{2}\right ) & : \phi^f \in [ \pi/2,\ldots,\pi ]       \end{array }     \right.\ ] ] and thus @xmath125\\         \pi - \frac{\phi^a}{2 } & : \phi^f \in [ \pi/2,\ldots,\pi ]       \end{array }     \right.\ ] ] this is a one - to - one mapping , covering the entire range of both @xmath126 and @xmath127 and which bijectively maps @xmath128 .
we can now carry out the transformation according to _ 0^d^fp(^f)&=&_0^d^fp(^f ) + _ ^d^fp(^f ) + & = & _ 0^ p(^a)d^a  .
we obtain for the adjoint phase @xmath118 p(^a)&=&p^(t_c)(^a ) ( -a(t)-c(t)^2 ( ) ) + & & ( b(t ) ( ) )  , [ eq : phia_distrib1 ] with p^(t_c)(^a)=(^a)=^2 ( )  .
likewise , for the adjoint loop @xmath127 we obtain p(tr*l*^a)&=&p^(t_c)(tr*l*^a)(-a(t)-(1+tr*l*^a ) ) + & & ( )  , with p^(t_c)(tr*l*^a)=  .
the simulation results ( figs .  [
fig : su2adj1 ] and [ fig : su2evadj1 ] ) are again reproduced to such accuracy that the curves are virtually indistinguishable .
we carry out the legendre transformations  ( [ eq : legendre2 ] ) to obtain @xmath129 and @xmath130 .
the results are shown in figs .
[ fig : su2adj2eff ] and [ fig : su2evadj2eff ] .
unfortunately , approximating these results via the saddle - point method is not possible due to the lack of a proper expansion point .
( top ) and @xmath59 ( bottom ) .
the different temperatures correspond to @xmath63 respectively .
[ fig : su2adj2eff],title="fig : " ]    ( top ) and @xmath59 ( bottom ) .
the different temperatures correspond to @xmath63 respectively .
[ fig : su2adj2eff],title="fig : " ]     ( top ) and @xmath59 ( bottom ) .
the different temperatures correspond to @xmath63 respectively .
[ fig : su2evadj2eff],title="fig : " ]   ( top ) and @xmath59 ( bottom ) .
the different temperatures correspond to @xmath63 respectively .
[ fig : su2evadj2eff],title="fig : " ]    we again obtain a taylor expansion of the classical potential @xmath131 . expanding ( [ eq : phia_distrib1 ] ) around @xmath132 to @xmath110 yields v_0(^a)&&-()+a(t ) + & & + ( 1+c(t)-- ) ( ) ^2 + & & + ( + + ) ( ) ^4 [ eq : phia_distrib2 ] in fig .  [ fig : su2evadj_taylor ] a comparison of  ( [ eq : phia_distrib2 ] ) and the simulation for @xmath102 is shown .
a fourth order expansion reproduces the data with high accuracy .
an expansion around @xmath106 is not possible due to the logarithmic term in the potential .
( top ) and @xmath59 ( bottom ) .
@xmath102 results compared to taylor expansions .
[ fig : su2evadj_taylor],title="fig : " ]   ( top ) and @xmath59 ( bottom ) .
@xmath102 results compared to taylor expansions .
[ fig : su2evadj_taylor],title="fig : " ]
in this work we simulated pure su(2 ) gauge theory at finite temperature , for two different fixed values of @xmath64 and hence of the lattice spacing .
we changed the temperature by changing the time - like extent @xmath45 of the lattice .
for each temperature we computed the classical potential of the bare polyakov loop and of the phase of the eigenvalues of the thermal wilson loop ( wrapped around the periodic boundary of the time direction ) in the fundamental and adjoint representation , using a histogram method . from these potentials
, we obtained the effective potential via legendre transformation .
we obtained model functions for the classical potentials for each case ( three parameters are required ) and investigated how the parameters depend on temperature and lattice spacing .
we discussed how the models can be approximated by taylor expansions and how they can be converted into one another .
we showed how the effective potential of the fundamental loop can be approximated to high precision by analytical expressions , when a next - to - leading order saddle expansion is used to compute the legendre transformation .
several conclusions may be drawn from our work : the most obvious conclusion is that both the classical and effective potentials of the bare loop and its eigenvalues are well described by simple models , reminiscent of landau - ginzburg theories but with a symmetry breaking term , where the whole dependence on temperature and lattice - spacing is absorbed into three parameters .
we have thus confirmed that the ansatz used in refs .
@xcite to parametrize the classical potential in a su(2 ) matrix model of wilson lines is valid also in the full gauge theory , with the same qualitative behavior of the parameters .
also , we have shown that a next - to - leading order approximation of the effective potential of the fundamental loop is reasonable .
we stress that what we have computed in this work is the effective potential for the _ bare _ loop , which vanishes in the continuum limit . in a given irreducible representation of the gauge group ,
a renormalized loop is extracted from the bare loop by dividing by a renormalization constant , which depends upon the representation @xcite .
one of the motivations of this work was the idea that the eigenvalues of the wilson line might be insensitive to renormalization .
our results contradict this .
instead , the eigenvalues of the bare loop are dominated by the effects of the vandermonde determinant . on the lattice ,
the vandermonde determinant arises naturally at each point in spacetime . in the action ,
its contribution is proportional to @xmath133 , where @xmath134 is the lattice spacing , and so it diverges in the continuum limit .
in retrospect , it is natural that the lattice effective potential for a site potential is dominated by a vandermonde term .    for the future
, it will be interesting to investigate how the effective potential of renormalized loops can be obtained from our results .
also , one should extract the non - perturbative contribution to the effective potential from our data and compare it to model calculations .
furthermore , our investigation should be extended to higher representations of su(2 ) and also to other groups , such as su(3 ) or exceptional groups with a trivial center , such as g(2 ) .
investigating the effects of dynamical fermions on the results discussed here is currently under way , using standard staggered fermions
. it would be interesting to extend this investigation to non - zero chemical potential , which is possible for two colors , where there is no sign problem .
we thank david scheffler for helpful discussions , for proof - reading and for a script to create histograms .    this work was supported by the deutsche forschungsgemeinschaft within sfb 634 , by the helmholtz international center for fair within the loewe initiative of the state of hesse , by the u.s .
department of energy under contracts # de - fg02 - 09er41620 and # de - ac02 - 98ch10886 , by the city university of new york through the psc - cuny research award program , grant 66514 - 00  44 , and by the european commission , fp7-people-2009-rg , no .
249203 .
all results presented here were obtained using nvidia geforce gtx 580 graphics cards .
to fix the physical units for the results discussed in this work we use a standard method which determines the lattice spacing @xmath134 through measurements of the string tension in lattice units at zero temperature ( which in practice means lattices of large volume which are either hypercubic or have @xmath135 ) . in this appendix
we give a short summary of the method ( which is derived in detail in several textbooks ) and present numerical results for @xmath136 and related quantities .
of course the running of the coupling with the cut - off scale has been investigated both analytically and numerically for both two- and three - color gauge theory , in many previous works ( see e.g. @xcite and references therein ) .
our main motivation for repeating such an investigation is to make use of the advanced computing power of the gpus to measure at larger @xmath137 and @xmath45 and with larger statistics as was previously possible and thus reduce systematic and statistical errors .
all physical units quoted in this work implicitly refer to the results presented in this appendix .    in pure gauge theory the static quark - antiquark potential
is known @xcite to be well described by v(r)=a+(1/r)+r + o(1/ r^2).[eq : quarkpot3 ] the linear coefficient @xmath138 is the string tension . on the lattice we can measure the potential in lattice units , _
i.e. _ @xmath139 , through the use of rectangular wilson loops @xmath140 in the @xmath141 plane ( our definition of @xmath142 implies that the trace is taken ) .
these are related to the potential by w(n_x , n_t ) = e^-n_t a v(a n_x ) ( 1+o ( e^-n_t a e ) )  , [ eq : quarkpot1 ] where @xmath143 is the energy gap between the ground state and the first excited state of a quark - antiquark pair .
since the second term is exponentially suppressed with @xmath144 , in the limit of large @xmath144 the potential can be directly obtained through the use of _ creutz ratios _ a v(a n_x )  .[eq : quarkpot2 ] since eq .
( [ eq : quarkpot3 ] ) is written in lattice units as a v(a n_x)=a a+(1/n_x)+ a^2 n_x + o(1/ ( a n_x ) ^2 )  , [ eq : quarkpot4 ] it is clear that what we actually obtain from a @xmath90 fit of the data is the dimensionless product @xmath145 .
fixing the string tension @xmath138 to its physical @xmath146 value ( or some other definite value @xmath147 ) then determines @xmath134 for a given @xmath64 .
the lattice spacing @xmath134 depends on the choice of @xmath138 of course , and for theories where this is not directly related to observable physics ( as in su(2 ) gauge theory ) one in fact often leaves @xmath138 undefined , which implies that @xmath134 and all derived quantities are expressed in units of an external parameter .
the crucial point here is that the @xmath64 dependence of the lattice spacing is in fact temperature independent .
this implies that once @xmath148 is known at @xmath146 , this fixes the temperature ( through @xmath149 ) and the physical units of all measured quantities for any given @xmath45 and @xmath64 .    to improve the signal - to - noise ratio for our string tension measurements
we employ a method known as ape smearing @xcite .
this method consists of repeatedly replacing all gauge links @xmath150 with averages of the form u_(n)=(1-)u_(n)+ _
c_(n )  , [ eq : apesmear1 ] where @xmath151 represents the staple matrix c_(n)&=&u_(n ) u_(n+ ) u_^(n+ ) + & + & u_^(n- ) u_(n- ) u_(n-+ )  , [ eq : apesmear2 ] and re - unitarizing the link ( projecting back to @xmath152 ) after each iteration
. the crucial effect of this procedure is to enhance the ground state overlap .
ape smearing is applied to a copy of the lattice ( the monte - carlo thus remains unaffected ) on which the wilson loops are subsequently measured .
a good choice for the parameter @xmath153 is known @xcite to be = 1 -  ;   w = 0.2   = 0.5454   .
we adopt this choice as well and apply the smearing @xmath154 times before measuring wilson loops on a given lattice configuration .
we measure @xmath155 on lattices of sizes @xmath156 , @xmath157 and @xmath158 , for @xmath64 values ranging from @xmath159 to @xmath160 using ( [ eq : quarkpot2 ] ) and ( [ eq : quarkpot4 ] ) .
sample sizes range from measurements on each site of several hundreds ( for the larger @xmath64 values ) to several tens of thousands ( for the smaller @xmath64 values ) of independent ( smeared ) gauge field configurations . from these results
we obtain @xmath161 , which is shown in fig .
[ fig : runcoupl1 ] .
wilson loops are measured up to spatial and timelike extents of @xmath162 and @xmath163 respectively .
we find consistency with the data sets quoted in ref .
@xcite .
one must consider potential sources of systematic errors .
firstly there are finite volume and discretization errors , for which we check by comparing the results obtained from lattices of different dimensions .
we find that our different choices give consistent results in @xmath64 ranges where they overlap ( at least two different lattices overlap in any given region ) , therefore we consider these errors to be under control .
second , one must take care to choose the timelike extent @xmath144 of the wilson loops large enough , such that the contribution from the exited states to the creutz ratios ( [ eq : quarkpot2 ] ) can be neglected .
this is tricky since the signal diminishes with rising @xmath144 .
we find that we can obtain , with reasonably small statistical errors , consistency between @xmath166 and @xmath167 on the two smaller lattices for each value of @xmath64 .
on the largest lattice consistent results are obtained on @xmath168 .
[ fig : runcoupl1 ] plots both respective choices for each lattice as the same dataset .
we consider this source of error to be under control as well .
the last , and potentially most severe , source of systematic errors are the higher order contributions to the potential ( [ eq : quarkpot3 ] ) , which one must account for .
we find that an excellent fit of the data can be obtained by discarding higher order terms and restricting the fit to @xmath169 .
alternatively one could introduce an additional term @xmath170 , or throw away more of the short - distance results .
although we find that the resulting string - tension measurements are much noisier if either of the latter two options is chosen , we find evidence that using @xmath169 introduces a systematic error such that the string tension is over - estimated ( and thus the temperature is under - estimated ) for the region @xmath171 if higher - order terms are excluded .
we choose to combine two sets of string tension data as our final results .
one with @xmath169 where the next term of ( [ eq : quarkpot3 ] ) is included and one without such a term but fitting only to @xmath172 .
we find that they are consistent .
we take this as evidence that the systematic errors from higher - order terms are small for our final results .
the data points shown in the figures in this section represent the measurements with the smallest errorbars . to avoid cluttering of the figures ,
we refrain from showing additional points which were obtained but come with huge errorbars .
all points were used however to determine the parameters of the model functions discussed in the following .
our main goal is to obtain the @xmath64 dependence of @xmath8 and its inverse for different @xmath45 .
this can be achieved if both @xmath175 for a given @xmath45 and @xmath176 are known .
furthermore , it is desirable to obtain not only a finite set of values but model - functions for each of these , valid over a large as possible range . a reasonable approach is to attempt to fit the data with known perturbative formulae , possibly with additional correction terms .    from @xmath161
we directly obtain @xmath177 .
the results are plotted in fig .
[ fig : runcoupl2 ] for different @xmath45 .
we find that it is possible to model our results by taking the inverse of the leading logarithmic term of the perturbative running - coupling formula ( see e.g. @xcite ) .
defining r()= ( )  , [ eq : runcouplfit2 ] and using a =  , n_t = r ( )  .[eq : runcouplfit1 ] we find that our data shown in figs .
[ fig : runcoupl1 ] and [ fig : runcoupl2 ] are well described for the entire range of @xmath64 considered with the parameters d= 1.98(1 ) ,  b= 0.305(6 )  .[eq : runcouplfit3 ] this is similar to what was done in ref .
@xcite although the @xmath64 range considered here is larger . it is in fact somewhat amazing that a leading - order formula works so well over such a large range .
it is also very convenient , since eq .
( [ eq : runcouplfit2 ] ) can be easily inverted .
if we restrict our fit - window to the same as was used in @xcite , we obtain a consistent value for the parameter @xmath178 . if we know a @xmath179 for a given @xmath45 we are done , for we can then obtain = =( )  , [ eq : totc ] and likewise ( ) b+_c=  .
[ eq : runcouplfit4 ] our numerical results for @xmath180 can analogously be converted to @xmath8 .
the parameter @xmath138 drops out entirely .
we show @xmath64 as a function of @xmath8 in fig .
[ fig : runcoupl3 ] , both as raw data and model - function .
we collected @xmath76 values for @xmath182 from refs .
@xcite and fit them to interpolate for other @xmath45 values , using _
c(n_t)= a_0 + b_0 ( n_t ) - ( ( n_t))^c_0  , [ eq : nlo_bc ] we find that the next - to - leading logarithm is necessary to appropriately describe such a range of @xmath45 .
our fit - parameters converge to a_0 = 1.1579(6 )  ,  b_0=0.9398(1 )  ,  c_0=1.627(2 )  , with @xmath183 .
( [ eq : nlo_bc ] ) is used in eq .
( [ eq : totc ] ) to determine @xmath76 for all @xmath45 instead of the raw data .
one last thing should be mentioned : we realize that the results presented in this appendix go beyond what is necessary for the study presented in the main sections .
in particular , to obtain @xmath8 , when simulating on a lattice of time - like extent @xmath45 at a given @xmath64 which is known to be @xmath76 for a different @xmath184 , one can simply use the relation =  .[eq : totc_ntontc ] we conducted the investigation discussed here hoping that it will be useful in a more general setting and we have chosen to present @xmath8 values throughout this paper as our general method provides them . equation ( [ eq : totc_ntontc ] ) however is very useful to us to estimate errors .
it turns out that the @xmath8 values quoted throughout this paper differ by no more than @xmath185 from what ( [ eq : totc_ntontc ] ) predicts , which is quite satisfactory .
g.  t hooft , nucl .
b * 138 * , 1 ( 1978 ) .
a.  m.  polyakov , phys .
b * 72 * , 477 ( 1978 ) .
l.  susskind , phys .
d * 20 * , 2610 ( 1979 ) .
l.  d.  mclerran and b.  svetitsky , phys .
d * 24 * , 450 ( 1981 ) .
a.  dumitru , j.  lenaghan and r.  d.  pisarski , phys .  rev .
d * 71 * , 074004 ( 2005 ) [ hep - ph/0410294 ] .
a.  dumitru , y.  hatta , j.  lenaghan , k.  orginos and r.  d.  pisarski , phys .
d * 70 * , 034511 ( 2004 ) [ hep - th/0311223 ] .
n.  meisinger , t.  r.  miller and m.  c.  ogilvie , phys .
d * 65 * , 034009 ( 2002 ) [ hep - ph/0108009 ] .
p.  n.  meisinger and m.  c.  ogilvie , phys .
d * 65 * , 056013 ( 2002 ) [ hep - ph/0108026 ] .
r.  d.  pisarski , phys .
d * 74 * , 121703 ( 2006 ) [ hep - ph/0608242 ] .
a.  dumitru , y.  guo , y.  hidaka , c.  p.  k.  altes and r.  d.  pisarski , phys
.  rev .
d * 83 * , 034022 ( 2011 ) [ arxiv:1011.3820 [ hep - ph ] ] .
a.  dumitru , y.  guo , y.  hidaka , c.  p.  k.  altes and r.  d.  pisarski , phys .  rev .
d * 86 * , 105017 ( 2012 ) [ arxiv:1205.0137 [ hep - ph ] ] .
t.  umeda , s.  ejiri , s.  aoki , t.  hatsuda , k.  kanaya , y.  maezawa and h.  ohno , phys .
d * 79 * , 051501 ( 2009 ) [ arxiv:0809.2842 [ hep - lat ] ] .
m.  panero , phys .
lett .   * 103 * , 232001 ( 2009 ) [ arxiv:0907.3719 [ hep - lat ] ] . s.  datta and s.  gupta , phys .
d * 82 * , 114505 ( 2010 ) [ arxiv:1006.0938 [ hep - lat ] ] .
s.  borsanyi , g.  endrodi , z.  fodor , s.  d.  katz and k.  k.  szabo , jhep * 1207 * , 056 ( 2012 ) [ arxiv:1204.6184 [ hep - lat ] ] .
m.  caselle , l.  castagnini , a.  feo , f.  gliozzi , u.  gursoy , m.  panero and a.  schafer , jhep * 1205 * , 135 ( 2012 ) [ arxiv:1111.0580 [ hep - th ] ] . b.  lucini and m.  panero , phys .
rept .   * 526 * , 93 ( 2013 ) [ arxiv:1210.4997 [ hep - th ] ] .
l.  fister and j.  m.  pawlowski , arxiv:1301.4163 [ hep - ph ] .
j.  braun , h.  gies and j.  m.  pawlowski , phys .
b * 684 * , 262 ( 2010 ) [ arxiv:0708.2413 [ hep - th ] ] .
d.  smith , phys .
d * 82 * , 034503 ( 2010 ) [ arxiv:0911.4037 [ hep - lat ] ] .
d.  smith ( 2010 ) , _ lattice simulation of a center symmetric three dimensional effective theory for su(2 ) yang - mills _
thesis , goethe university of frankfurt , germany .
a.  d.  kennedy and b.  j.  pendleton , phys .
b * 156 * , 393 ( 1985 ) .
s.  gupta , k.  huebner and o.  kaczmarek , phys .
d * 77 * , 034503 ( 2008 ) [ arxiv:0711.2251 [ hep - lat ] ] . c.  s.  fischer , l.  fister , j.  luecker and j.  m.  pawlowski , arxiv:1306.6022 [ hep - ph ] . k.  langfeld and j.  m.  pawlowski , arxiv:1307.0455 [ hep - lat ] .
d.  diakonov , c.  gattringer and h.  -p .
schadler , jhep * 1208 * , 128 ( 2012 ) [ arxiv:1205.4768 [ hep - lat ] ] .
b.  a.  berg , _ markov chain monte carlo simulations and their statistical analysis _
, world scientific publishing , singapore , 2004    j.  greensite , phys .
d * 86 * , 114507 ( 2012 ) [ arxiv:1209.5697 [ hep - lat ] ] . j.  greensite and k.  langfeld , arxiv:1301.4977 [ hep - lat ] .
a.  dumitru and d.  smith , phys .
d * 77 * , 094022 ( 2008 ) [ arxiv:0711.0868 [ hep - lat ] ] .
d.  smith , nucl .
a * 820 * , 227 ( 2009 ) [ arxiv:0810.1129 [ hep - lat ] ] .
g.  s.  bali and k.  schilling , phys .
d * 47 * , 661 ( 1993 ) [ hep - lat/9208028 ] .
j.  c.  r.  bloch , a.  cucchieri , k.  langfeld and t.  mendes , nucl .
b * 687 * , 76 ( 2004 ) [ hep - lat/0312036 ] . c.  michael , phys .
b * 283 * , 103 ( 1992 ) [ hep - lat/9205010 ] .
a.  cucchieri , hep - lat/0209076 .
m.  luscher , r.  sommer , p.  weisz and u.  wolff , nucl .
b * 413 * , 481 ( 1994 ) [ hep - lat/9309005 ] .
s.  necco and r.  sommer , phys .
b * 523 * , 135 ( 2001 ) [ hep - ph/0109093 ] .
m.  luscher and p.  weisz , jhep * 0207 * , 049 ( 2002 ) [ hep - lat/0207003 ] .
m.  albanese _ et al . _
[ ape collaboration ] , phys .
b * 192 * , 163 ( 1987 ) .
n.  cardoso and p.  bicudo , j.  comput .
phys .   * 230 * , 3998 ( 2011 ) [ arxiv:1010.4834 [ hep - lat ] ] .
f.  geles and c.  b.  lang , arxiv:1103.5368 [ hep - lat ] . c.  morningstar and m.  j.  peardon , phys .
d * 69 * , 054501 ( 2004 ) [ hep - lat/0311018 ] .
j.  fingberg , u.  m.  heller and f.  karsch , nucl .
b * 392 * , 493 ( 1993 ) [ hep - lat/9208012 ] .
p.  de forcrand and l.  von smekal , phys .
d * 66 * , 011504 ( 2002 ) [ hep - lat/0107018 ] .
a.  velytsky , int .
j.  mod .
c * 19 * , 1079 ( 2008 ) [ arxiv:0711.0748 [ hep - lat ] ] .
x.  cheng and e.  t.  tomboulis , phys .
d * 86 * , 074507 ( 2012 ) [ arxiv:1206.3816 [ hep - lat ] ] .