quantum chromodynamics ( qcd ) and the gluon field were first proposed by nambu in 1966 @xcite just after the introduction of color degrees of freedom @xcite , and qcd has been established as the fundamental gauge theory of the strong interaction , through the explanation of asymptotic freedom @xcite , the success of perturbative qcd for high - energy phenomena based on the parton model @xcite , and lattice qcd calculations @xcite . in spite of many successes of qcd , there are still unsolved problems in the low - energy region of qcd , owing to its strong - coupling nature .
indeed , the mechanism of color confinement is regarded as one of the most difficult important problems in theoretical physics @xcite , and spontaneous chiral - symmetry breaking @xcite is also a difficult issue in describing it quantitatively directly from qcd @xcite .
one of the difficulties of qcd lies on the large gap between the fundamental fields ( quarks and gluons ) and the observable particles ( hadrons ) .
in fact , quarks and gluons , which are the building blocks of qcd , have no physical asymptotic states , as a result of color confinement or the gauge invariance .
this is in contrast to the ordinary perturbation theory , where all the observable phenomena can be directly described with the fundamental fields appearing in the lagrangian .
furthermore , the field - theoretical description of the confined particles is an interesting but unsolved difficult subject . in general ,
the green function is one of the most basic quantities to describe the motions and the interactions of particles @xcite , but most green s functions of quarks and gluons are still unknown in the nonperturbative description .    as for the essence of nonperturbative qcd ,
the central issue is gluon dynamics rather than quarks @xcite .
in fact , the strong gluon interaction makes the qcd vacuum highly nontrivial , and color confinement and chiral symmetry breaking are realized even at the quenched level @xcite .
then , the analysis of gluon properties is the key point to clarify the nonperturbative aspects of qcd . in particular ,
the gluon propagator , _
i.e. _ , the two - point green function is one of the most basic quantities in qcd , and has been investigated with much interest in various gauges , such as the landau gauge @xcite , the coulomb gauge @xcite , and the maximally abelian ( ma ) gauge @xcite , in the context of various aspects of qcd .
for example , an infrared - vanishing gluon propagator @xcite is proposed from the mathematical analysis of the gribov horizon @xcite , and the infrared singularity of the gluon propagator has been investigated in terms of color confinement @xcite , and also from the viewpoint of renormalon cancellation @xcite .
dynamical gluon - mass generation @xcite is also an important subject related to the infrared gluon propagation . in spite of massless perturbative gluons ,
the gluon field is conjectured to acquire a large effective mass as the self - energy through the self - interaction of gluons in a nonperturbative manner . for , the glueball states ,
color - singlet bound states of gluons , are considered to be fairly massive _ e.g. _ , about 1.5gev for the lowest @xmath19 and about 2gev for the lowest @xmath20 , as indicated in lattice qcd calculations @xcite .
so far , a large effective gluon mass estimated about @xmath21gev has been proposed in various context of qcd physics : analytical framework based on qcd @xcite , lattice qcd calculations @xcite , pomeron physics @xcite , glueball phenomenology @xcite , and heavy - quark phenomenology @xcite .
since the color @xmath22 symmetry is unbroken in qcd , the effective - mass generation of gluons is quite different from the higgs mechanism , which is a standard mass generation in quantum field theories . in the electro - weak unified theory , the weak bosons , @xmath23 and @xmath24 ,
obtain quite a large mass of about 100gev , as a result of the higgs mechanism , _
i.e. _ , spontaneous breaking of the @xmath25 gauge symmetry .
however , the gluon mass generation is not a result of spontaneous gauge - symmetry breaking , but stems from much more complicated nonperturbative dynamics of gluons , and generally depends on the gauge choice .
the nonperturbative effects originate from the strong - coupling infrared region of qcd .
actually , a recent lattice qcd study clarifies that the relevant energy scale for confinement is the infrared gluon component below 1.5gev @xcite : the string tension is almost unchanged even by cutting off high - momentum gluons above 1.5gev . reflecting the asymptotic freedom or the running coupling @xmath26 , qcd exhibits various different features according to the energy scale
here , we roughly classify three scale regions of qcd as  ultraviolet ( uv ) " ,  infrared ( ir)/intermediate ( i m ) " , and  deep - infrared ( deep - ir ) " regions , in terms of the length @xmath27 .
* we define the uv region as @xmath28 0.1fm , which corresponds to the high - energy region more than a few gev s . in this region ,
perturbative qcd is approximately applicable to the reaction process , and the inter - quark potential is almost coulomb - like @xcite .
* we define the ir / im region as 0.1fm @xmath29 1fm , which ranges from a few hundred mev to a few gev in energy . in this region ,
the system is described by quark - gluon degrees of freedom in a nonperturbative way , which is usually substituted with some effective models @xcite .
* we define the deep - ir region as @xmath30 1fm , which is low - energy below @xmath31 0.2gev . in this region ,
the perturbative running coupling @xmath26 diverges , and the confinement effect is extremely large , so that quark - gluon degrees freedom are hidden and the system is described by hadrons .
so far , the gluon propagator has been studied mainly in the landau gauge both in analytic framework @xcite and in lattice qcd @xcite . in the uv region
, the gluon propagator takes a massless perturbative form , @xmath32 , apart from the tensor factor . in the ir / im and the deep - ir regions , the perturbative approach breaks down , and we need a nonperturbative approach such as lattice qcd calculations @xcite .
basic important ideas were proposed and investigated in the pioneering early - time lattice studies @xcite , and recent lattice studies @xcite followed them and presented high - precision data with the main interests in the deep - ir behavior by using huge - volume lattices @xcite .
nevertheless , the functional form of the gluon propagator is unclear still now .    in this paper
, we study the functional form of the gluon propagator in the landau gauge in su(3 ) lattice qcd monte carlo calculations , especially for _ the ir / im region _ of @xmath33fm , which is considered to be relevant for the quark - hadron physics @xcite , and also aim to describe nonperturbative gluon properties , based on the obtained function form of the gluon propagator .
the organization of this paper is as follows . in sec.ii , we give the formalism of the landau gauge fixing and the gluon propagator both in continuum and in lattice qcd . in sec.iii , we show the lattice qcd result of the gluon propagator in the landau gauge both in the coordinate space and in the momentum space . in sec.iv
, we estimate the effective gluon mass from the gluon propagator and the effective - mass plot in lattice qcd .
in sec.v , we investigate the functional form of the gluon propagator in the landau gauge , by analyzing the lattice qcd data .
we show that the landau - gauge gluon propagator is well described with the yukawa - type function in the ir / im region . in sec.vi , as the applications of the yukawa - type gluon propagator , we derive analytic expressions for the zero - spatial - momentum propagator , the effective mass , and the spectral function of the gluon field .
section vii is devoted to summary and discussions .
in this section , we review the formalism of the landau gauge fixing and the gluon propagator in the euclidean space - time .
the landau gauge is one of the most popular gauges in qcd . as a remarkable feature
, the landau gauge keeps lorentz covariance and global @xmath22 symmetry . owing to these symmetries and the transverse property ,
the color and lorentz structure of the gluon propagator is uniquely determined .    in the euclidean space - time formalism such as lattice qcd
, the landau gauge is usually defined so as to minimize the gauge - field fluctuation .
we then expect that only the minimal fluctuation of the gluon field survives in the landau gauge , and the physical essence of gluon properties can be investigated without suffering from large stochastic fluctuations of gauge degrees of freedom .      to begin with , let us consider the landau gauge fixing and its physical meaning in euclidean qcd . in the su@xmath34 continuum qcd , the gluon field is expressed as @xmath35 with the generator @xmath36 and @xmath37 .
the landau gauge is usually defined by the local condition on the gauge field as @xmath38 in euclidean qcd , the landau gauge has a global definition to minimize the global quantity , @xmath39 by the gauge transformation .
this global definition is more strict , and the local condition ( [ eq : locallandau ] ) is derived from the minimization of @xmath40 .
the global quantity @xmath40 can be regarded as the total amount of the gauge - field fluctuation in the euclidean space - time . in the global definition ,
the landau gauge has a clear physical interpretation that it maximally suppresses the artificial gauge - field fluctuation relating to the gauge degrees of freedom .    in lattice qcd
, the theory is formulated on the discretized space - time @xcite .
the qcd action is constructed from the link - variable @xmath41 , instead of the gauge field @xmath42 .
the link - variable is defined as @xmath43 , with the lattice spacing @xmath44 and the gauge coupling constant @xmath45 .
the gauge transformation of the link - variable is given by @xmath46 with the gauge function @xmath47 .    in lattice qcd
, the landau gauge fixing is also expressed in terms of the link - variable : the landau gauge is defined by the maximization of @xmath48 by the gauge transformation of the link - variable .    for small @xmath44 ,
using the expansion @xmath49 in terms of lattice spacing @xmath44 , @xmath50 is expressed as @xmath51 therefore , the maximization of @xmath50 corresponds to the minimization of the gauge - field fluctuation as well as the continuum theory .
this minimization of the gluon - field fluctuation in turn justifies the expansion in eq.([eq : linkvariableexpansion ] ) .      in this subsection
, we formulate the gluon propagator in the landau gauge . to begin with
, we extract the gluon field from the gauge - fixed link - variable , which is obtained by the gauge transformation of the link - variable to maximize @xmath50 . using the expansion ( [ eq : linkvariableexpansion ] ) ,
we define the bare gluon field @xmath52 as @xmath53   - \frac{1}{2iagn_c } { \rm tr } \left [ u_\mu(x ) - u_\mu^\dagger(x ) \right ] , \label{eq : gluonfield}\ ] ] where the second term is added to make @xmath54 traceless .
the renormalized gluon field @xmath55 is obtained by multiplying a real renormalization factor @xmath56 as @xmath57 we abbreviate @xmath55 as @xmath58 hereafter .
the gluon field @xmath58 defined above is traceless and hermite , and is expressed as @xmath35 with @xmath37 . in the landau gauge which maximizes @xmath50
, the gluon field @xmath58 satisfies the local condition , @xmath59 with the forward or backward derivative @xmath60 on the lattice .
the gluon propagator @xmath0 is defined by the two - point function as @xmath61 note that the time - ordered product is unnecessary in the euclidean metric , and the translational invariance of the vacuum leads to the @xmath62-dependence . from eq.([eq : locallandaul ] ) , @xmath63 satisfies the transverse property , @xmath64    next , we consider the gluon propagator @xmath65 in the momentum space , which is defined by the fourier transformation of the coordinate - space propagator as @xmath66 on the @xmath67 lattice , this fourier transformation is discretized , and the momentum - space gluon propagator is expressed as @xmath68 here , the discretized momentum @xmath69 and the continuum momentum @xmath70 are defined as @xcite @xmath71 with @xmath72 .    in the landau gauge ,
the color and tensor structure of the propagator is uniquely determined as @xmath73 from the @xmath22 global symmetry , the lorentz symmetry , and the transverse property , @xmath74 therefore , we only have to consider the scalar factor of the gluon propagator , @xmath75 , which is a function of the continuum - momentum squared @xmath76 .    in the coordinate space
, we investigate the scalar combination of the gluon propagator @xmath77 as a function of the four - dimensional euclidean distance , @xmath78 here , the denominator factor @xmath79 in eq.([eq : defpropagator ] ) has been introduced considering the tensor factor @xmath80 .    the scalar factor @xmath81 in eq.([eq : gpmom ] )
is expressed by the fourier transformation of @xmath82 as @xmath83 from which one can prove that @xmath82 depends only on @xmath84 near the continuum limit . in this paper , we call @xmath82 and @xmath81  scalar - type propagator " .
we perform @xmath85 lattice qcd monte carlo calculations at the quenched level using the standard plaquette action . here
, we adopt three different lattices with the lattice parameter @xmath86=5.7 , 5.8 , and 6.0 .
the used lattice size is @xmath87 , @xmath88 , and @xmath89 at @xmath1 = 5.7 , 5.8 , and 6.0 , respectively . in this lattice calculation
, we mainly use @xmath90 .
the lattice spacing @xmath44 is found to be @xmath91 , and @xmath92fm , at @xmath1 = 5.7 , 5.8 , and 6.0 , respectively , when the scale is determined so as to reproduce the string tension as @xmath93mev from the static q@xmath94 potential @xcite .
the number of used gauge configurations is 50 , 40 and 30 for @xmath1 = 5.7 , 5.8 , and 6.0 , respectively .
the gauge configurations are picked up every 1,000 sweeps after a thermalization of 20,000 sweeps .
we summarize the parameter and calculation conditions in table  [ tab : latparam ] .
.[tab : latparam ] the lattice parameter @xmath1 , lattice size , and the gauge - configuration number @xmath95 .
the corresponding lattice spacing @xmath44 and the lattice volume in the physical unit are added .
the lattice spacing @xmath44 is determined so as to reproduce the string tension @xmath93mev . [ cols="^,^,^,^,^",options="header " , ]     as a remarkable fact , the yukawa - type function @xmath96 well reproduces the lattice qcd data in the whole region of @xmath7fm . on the other hand , the dipole - type function @xmath97 fails to reproduce the whole region of the lattice data , since it gives too strong reduction at the short distance
. the massive propagator @xmath98 also fails to reproduce the whole region of the lattice data , as was already shown .    to see the difference of the three fit results clearer
, we show in fig.[fig : threefittingratio ] the ratio of the lattice qcd data @xmath99 to the three fit - functions on the scalar - type gluon propagator , _
i.e. _ , @xmath100 , @xmath101 , and @xmath102 .
one finds @xmath103 in the whole region of @xmath104fm , while @xmath100 and @xmath102 differ from the unity for small @xmath84 .     the ratio of the lattice qcd data @xmath99 at @xmath1 = 6.0 to the fit - functions @xmath98 , @xmath96 and @xmath97 on the scalar - type gluon propagator , _
i.e. _ , @xmath100 , @xmath101 , and @xmath102 .
]      we thus find an appropriate functional form of the gluon propagator in the landau gauge .
in fact , the scalar - type gluon propagator @xmath82 in the coordinate space is well described by the yukawa - type function @xmath105 with @xmath106gev and @xmath107 in the whole range of @xmath33fm .    as a summary figure
, we show in fig.[fig : yukawasummary ] the comparison between the obtained yukawa - type function @xmath96 and all the lattice qcd data of the scalar - type gluon propagator @xmath82 at @xmath1=5.7 , 5.8 , and 6.0 in the range of @xmath33fm . for clearer presentation of the yukawa - functional behavior of the lattice gluon propagator @xmath82
, we also show the logarithmic plot of @xmath108 in fig.[fig : yukawasummary](b ) .
all the lattice qcd data at @xmath1=5.7 , 5.8 , and 6.0 are found to be well reproduced with the yukawa - type function in the range of @xmath33fm .
note also that all the lattice data of @xmath108 locate around a straight line in the logarithmic plot of fig.[fig : yukawasummary](b ) .
the yukawa - type function @xmath105 ( solid line ) with @xmath109gev and @xmath110 obtained by the fit analysis at @xmath90 , and the lattice qcd data of the scalar - type gluon propagator @xmath82 in the landau gauge at @xmath1=5.7 , 5.8 , and 6.0 in the range of @xmath33fm .
the lower figure is the logarithmic plot of @xmath108 and @xmath111 .
, title="fig : " ]   the yukawa - type function @xmath105 ( solid line ) with @xmath109gev and @xmath110 obtained by the fit analysis at @xmath90 , and the lattice qcd data of the scalar - type gluon propagator @xmath82 in the landau gauge at @xmath1=5.7 , 5.8 , and 6.0 in the range of @xmath33fm . the lower figure is the logarithmic plot of @xmath108 and @xmath111 .
, title="fig : " ]    next , we investigate the gluon propagator @xmath75 in the momentum space in terms of the yukawa - type function . in fig.[fig : yukawamom ] , we show the scalar - type gluon propagator @xmath75 in the landau gauge obtained in lattice qcd at @xmath112 , and @xmath113 , which is the fourier transformation of the yukawa - type function @xmath96 .
the horizontal axis is @xmath114 . here
, we use the same parameters @xmath115gev and @xmath116 as those used in the best - fit analysis for the coordinate - space gluon propagator . from figs .
[ fig : latresmom ] and [ fig : yukawamom ] , the lattice qcd data of @xmath81 at @xmath1=5.7 , 5.8 , and 6.0 are found to be approximated with @xmath117 in the range of @xmath118 .
we also perform the best - fit analysis of the lattice data of the gluon propagator @xmath81 in the momentum space with @xmath113 . for the lattice data at @xmath1=6.0 in the fit range of @xmath119 ,
the best - fit parameters are found to be @xmath120gev and @xmath121 , which are close to the values obtained from the fit - analysis of the coordinate - space gluon propagator .
the yukawa - type propagator in the momentum space , _ i.e. _ , @xmath113 ( solid line ) with @xmath115gev and @xmath116 , the same values used in fig.[fig : yukawasummary ] .
the horizontal axis is @xmath114 .
the symbols denote the lattice qcd data of the scalar - type gluon propagator @xmath81 in the landau gauge at @xmath122 , where the momentum is defined as @xmath123 . ]    as a caution , in the uv region of @xmath124 , the gluon propagator @xmath75 gradually deviates from the yukawa form @xmath125 and gradually approaches the perturbative form @xmath32 .
however , the coordinate - space gluon propagator @xmath82 is found to be almost unchanged for @xmath33fm by correcting the uv deviation .    also in the deep - ir region of @xmath126 , as is briefly summarized in appendix  b
, there appears some deviation between @xmath75 and @xmath127 .
in fact , such a deep - ir deviation is indicated by recent huge - volume lattice qcd studies @xcite . in the momentum space ,
the true gluon propagator @xmath75 turns out to take a saturated value smaller than the yukawa - type propagator @xmath127 in the deep - ir region of @xmath128gev . in other words , @xmath129gev is the lower bound on the applicability of the yukawa - type propagator @xmath127 to the gluon propagator .
based on the huge - volume lattice data @xcite , we investigate the deep - ir - corrected gluon propagator in the coordinate space in appendix  b. as the conclusion , even after the correction in the deep - ir region , the yukawa - type function is found to work well for @xmath7fm .
( see fig.[fig : hugelatticecorr ] in appendix  b. )    as the main conclusion of this paper , we summarize the functional form of the gluon propagator in the landau gauge obtained in su(3 ) lattice qcd .    1 .   the coordinate - space gluon propagator @xmath82 in the landau gauge is well described by the four - dimensional yukawa - type function as @xmath130 with @xmath131 600mev and @xmath132 0.16 , for the whole region of @xmath133fm .
( this is valid even after the possible correction in the uv and deep - ir regions . )
the gluon propagator @xmath81 in the momentum space is also well described by the corresponding new - type propagator ( four - dimensional fourier transformed yukawa - type function ) as @xmath134 with @xmath131 600mev and @xmath132 0.16 ( same values ) , in the momentum region of @xmath135 .
note here that all the component of the gluon propagator @xmath136 and @xmath65 in the landau gauge can be analytically expressed , starting from the yukawa - type function .
( see sec.ii . )
quantitatively , the yukawa - type propagator @xmath82 exhibits a slower decreasing feature , and @xmath75 exhibits faster decreasing , in comparison with the ordinary massive propagator .
it seems suggestive to rewrite @xmath75 as @xmath137 where @xmath138 corresponds to the wave - function renormalization of the gluon field , in a similar manner to the schwinger - dyson formalism .
near the on - shell - like condition of @xmath139 , @xmath138 tends to diverge as @xmath140 , which leads to anomalous gluon propagation and may mimic the gluon confinement , _ i.e. _ , the absence of on - shell gluon states .    to see the physical meaning of the massive parameter @xmath14 in the yukawa - type function , we compare the yukawa - type propagator form @xmath5 with the massive propagator form @xmath141 in the coordinate space . in spite of a significant difference at the short distance ,
their difference is just the prefactor at the large distance , where @xmath142 and the main reduction factor is @xmath143 .
then , the mass parameter @xmath131 600mev in the yukawa - type gluon propagator directly corresponds to the effective gluons mass in the infrared region .
( see sec.vi c. )    here , we briefly comment on the other functional forms for the gluon propagator . up to now , many functional forms of the gluon propagator @xmath75 have been considered and compared with the lattice qcd result @xcite .
some functional forms well describe the gluon propagator better than the yukawa - type function , but they need four or more fit parameters and take highly non - analytical complicated form . on the other hand ,
the yukawa - type function @xmath96 has _ only two parameters _ @xmath144 and takes an analytical form , which are its advantages . for example , owing to the analyticity , the yukawa - type gluon propagator @xmath96 leads to analytical applications of gluonic nonperturbative quantities , as will be demonstrated in the next section .
in this section , as the applications of the yukawa - type gluon propagator , we derive analytical expressions for the zero - spatial - momentum propagator @xmath11 , the effective mass @xmath12 , and the spectral function @xmath13 of the gluon field .
all the derivations can be analytically performed , starting from the yukawa - type gluon propagator @xmath96 .
second , we consider the zero - spatial - momentum propagator @xmath11 , associated with the yukawa - type propagator @xmath96 .
@xmath11 was introduced in sec.iv-b in the context of the effective mass in lattice qcd . here
, we mainly deal with the continuum formalism with infinite spatial volume . for the simple argument
, we first neglect the temporal periodicity , which is justified for large temporal lattice size .
we start from the yukawa - type gluon propagator , @xmath145 with @xmath146 . like eq.([eq : zmpcorrl ] ) , the zero - spatial - momentum propagator is given by @xmath147 in the continuum formalism .
using the three - dimensional polar coordinate of @xmath148 , we calculate this integral as follows : @xmath149 with @xmath150 . here , we have used eq.([eq : integformmodbessel ] ) on the modified bessel function .
thus , we derive an analytical expression for the zero - spatial - momentum propagator , @xmath151    for the actual comparison with the lattice qcd data , we take account of the temporal periodicity , which is used in lattice calculations . in this case , @xmath11 is given as @xmath152 .
\label{eq : corranalytical}\ ] ] in fig.[fig : corryukawa ] , we show the lattice qcd result of @xmath11 in the landau gauge , and the theoretical curve of eq.([eq : corranalytical ] ) with @xmath14=0.624gev and @xmath153=0.162 , the same values used in the previous section .
the lattice qcd data of @xmath11 are found to be well described by the theoretical curve , associated with the yukawa - type gluon propagator .     the zero - spatial - momentum propagator @xmath11 of gluons in the landau gauge .
the symbols are the lattice qcd data at @xmath112 , and the solid line denotes the theoretical curve of eq.([eq : corranalytical ] ) , derived from the yukawa - type propagator with @xmath14=0.624gev and @xmath153=0.162 , the same values in fig.[fig : yukawasummary ] . ]
we also consider @xmath11 in the lattice formalism , for more direct comparison with the lattice qcd data . here
, we use the expression of @xmath11 with @xmath127 , @xmath154 on the lattice , @xmath155 is given as @xmath156 with @xmath157 ( @xmath158 ) in the lattice unit .
then , we obtain an analytical expression for the zero - spatial - momentum propagator , @xmath159 no significant numerical difference turns out to be found between the continuum and the lattice formulae , under the condition of our lattice qcd calculation .
this means that the integral quantity @xmath11 is not so sensitive to the details of the uv - behavior of the propagator @xmath96 .    in the uv region ,
the yukawa - type propagator @xmath160 deviates from the correct behavior of the perturbative propagator @xmath32 . for the quantitative estimate of the influence from the deviation in the uv region
, we calculate the zero - spatial - momentum propagator @xmath11 using the uv - corrected yukawa propagator @xmath161 , which is @xmath127 for @xmath162gev and @xmath32 in the uv region of @xmath163gev .
the difference of @xmath11 is found to be at most 1% between the cases with @xmath127 and @xmath161 .
thus , the integral quantity @xmath11 is insensitive to the uv behavior of the gluon propagator , so that it is also expected to be insensitive to the discritization error @xcite in the uv region in the lattice calculation .
second , we investigate the effective mass @xmath12 , as the consequence of the yukawa - type propagator @xmath96 . the effective - mass plot is a general useful technique for the mass estimation in lattice qcd , and was actually examined for gluons in sec.iv-b . for simplicity
, we here treat the three - dimensional space as a continuous infinite - volume space , while the temporal variable @xmath164 is discrete and periodic . for the simple notation , we here use the lattice unit for @xmath164 .
when the temporal periodicity can be neglected , the zero - spatial - momentum propagator @xmath11 is expressed by eq.([eq : d0yukawa ] ) , and we obtain an analytical expression of the effective mass , @xmath165 from the asymptotic form of @xmath166 , the effective mass of gluons is approximated as @xmath167 for large @xmath164 .
this functional form indicates that @xmath12 is an increasing function and approaches @xmath14 from below , as @xmath164 increases .
note that the mass parameter @xmath14 in the yukawa - type gluon propagator directly corresponds to the effective mass @xmath12 of gluons for large @xmath164 .
in fact , @xmath131 600mev has a definite physical meaning of the effective gluon mass in the infrared region .
note also that the simple analytical expression of eq.([eq : emgyukawa ] ) or ( [ eq : meffasymptotic ] ) reproduces the anomalous increasing behavior of the effective mass @xmath12 of gluons , as was observed in fig.[fig : effmass ] . as for the increasing behavior ,
there is a general argument : this can occur when the spectral function is not positive definite @xcite , although the concrete form of the gluon spectral function is not yet known .
instead , this framework with the yukawa - type gluon propagator gives an analytical and quantitative method , and is found to well reproduce the lattice result .
( the actual comparison is demonstrated with eq.([eq : meffanaytical ] ) . )
next , we take account of the temporal periodicity , which is used in lattice qcd calculations . in this case , the effective mass @xmath12 is defined by
cosh - type " as eq.([eq : emcosh ] ) , and the zero - spatial - momentum propagator @xmath11 is given by eq.([eq : corranalytical ] ) .
then , the effective mass @xmath12 of gluons is expressed as @xmath168 } { \cosh\left [ m_{\rm eff}(t ) ( n_t/2-t)\right]}= \\ & & \frac{(t+1)k_1(m(t+1 ) ) + ( n_t -(t+1))k_1(m(n_t-(t+1)))}{tk_1(mt)+(n_t - t)k_1(m(n_t - t))}. \nonumber\end{aligned}\ ] ]    in fig.[fig : effmassyukawa ] , we show the theoretical curve obtained by eq.([eq : meffanaytical ] ) together with the lattice result of @xmath12 . here
, we take @xmath14=0.624gev , the same value used in the previous section .
the lattice qcd data of @xmath12 are found to be well described by the theoretical curve , associated with the yukawa - type gluon propagator .
the effective mass @xmath12 of gluons in the landau gauge .
the symbols denote the lattice qcd data at @xmath112 , and the solid line denotes the theoretical curve of eq.([eq : meffanaytical ] ) derived from the yukawa - type propagator with @xmath14=0.624gev , the same value used in fig.[fig : yukawasummary ] . ]      from the analytical expression of the zero - spatial - momentum propagator @xmath169 in eq.([eq : d0yukawa ] ) , we can derive the spectral function @xmath13 of the gluon field , associated with the yukawa - type gluon propagator .
for simplicity , we take continuum formalism with infinite space - time .
the relation between the spectral function @xmath13 and the zero - spatial - momentum propagator @xmath11 is given by the laplace transformation , @xmath170 when the spectral function is given by a @xmath17-function such as @xmath171 , which corresponds to a single mass spectrum , one finds a familiar relation of @xmath172 . for the physical state ,
the spectral function @xmath13 gives a probability factor , and is non - negative definite in the whole region of @xmath173 .
this property is related to the unitarity of the @xmath174-matrix .    in general ,
the laplace transformation is defined by @xmath175 and the inverse laplace transformation is expressed as @xmath176    then , from eq.([eq : corrvssf ] ) , the spectral function @xmath13 is expressed as @xmath177 with @xmath178 , @xmath179 and @xmath180 .    performing the partial integration in eq.([eq : integformmodbessel ] ) on the modified bessel function
, we obtain a formula of the laplace transformation , @xmath181 which leads to the inverse laplace transformation , @xmath182 by differentiating this formula by @xmath173 , we find @xmath183 where the second term includes an infinite factor besides the @xmath17-function .
then , we apply this formula to eq.([eq : sfint ] ) , and obtain the spectral function @xmath13 as @xmath184    for more rigorous derivation , we avoid the singularity at @xmath18 by regularizing eq.([eq : ltk1 ] ) as @xmath185 with a positive infinitesimal @xmath186 , and find the formula of @xmath187 which leads to the regularized spectral function , @xmath188 in the calculation of the laplace transformation from @xmath13 to @xmath11 , we can avoid the divergence of the integral at @xmath18 by using the regularized spectral function @xmath189 .
then , we can perform the integration , and properly obtain @xmath190 , by taking the limit of @xmath191 after the integration .
we have numerically confirmed that the spectral function @xmath192 in eq.([eq : sfyukawar ] ) reproduces @xmath190 by the laplace transformation in the limit of @xmath191 .    in this way , we derive the spectral function @xmath13 of the gluon field , associated with the yukawa - type propagator : @xmath193 which is regularized as @xmath189 in more rigorous derivation . here , @xmath131 600mev is the mass parameter in the yukawa - type function for the landau - gauge gluon propagator .
the first term expresses a negative continuum spectrum , and the second term a @xmath17-functional peak with the residue including an infinite factor , which is positive as @xmath194 at @xmath195 .
we show in fig.[fig : spectral ] the spectral function @xmath13 of the gluon field .
although the appearance of the negative - value region in the gluon spectral function is expected , @xmath13 exhibits two anomalous behaviors : it has a positive @xmath17-functional peak with the residue of @xmath140 at @xmath196 , and it takes negative values for all the region of @xmath197 .
this negative contribution of the spectral function gives unusual behavior of the effective mass @xmath12 of gluons , _
i.e. _ , its increasing behavior on @xmath164 .
the spectral function @xmath13 of the gluon field , associated with the yukawa - type propagator .
the unit is normalized by the mass parameter @xmath131 600mev . as eq.([eq : sfyukawa ] ) indicates , @xmath13 shows anomalous behaviors : it has a positive @xmath17-functional peak with the residue of @xmath140 at @xmath196 , and takes negative values for all the region of @xmath197 . ]    as was discussed in sec.iv , if the spectral function is non - negative definite , like the ordinary hadronic correlator , the effective mass must be a monotonously decreasing function , due to the mathematical nature of the summation of @xmath198 with non - negative coefficients , as can be proven with eq.([eq : dmdt ] ) .
physically , this is due to a larger positive contribution of the excited states to the effective mass at a shorter distance .
however , if the spectral function includes the negative - value region , such a definite tendency is lost @xcite .
this is the case of the gluon field in the landau gauge .
more precisely , from eq.([eq : dmdt ] ) , we find that the increasing property of @xmath12 can be realized , only when there is some suitable coexistence of positive- and negative - value regions in the spectral function @xmath13 . as a remarkable fact ,
the obtained gluon spectral function @xmath13 is negative - definite for all the region of @xmath197 , except for the positive @xmath17-functional peak at @xmath18 .
the negative property of the spectral function in coexistence with the positive peak leads to the anomalous increasing behavior of the effective mass @xmath12 .    in the uv region ,
the true gluon propagator deviates from the yukawa - type function .
then , we estimate the contribution from the uv part of this spectral function @xmath13 to the zero - spatial - momentum propagator @xmath11 . introducing the uv cutoff @xmath199
, we define the integration @xmath200 which reduces to the laplace transformation in the @xmath201 limit . for @xmath202 0.1fm , the @xmath199-dependence of the integral quantity @xmath203
is found to be negligible for @xmath204 gev , _
i.e. _ , @xmath205 . in other words , due to this uv - insensitivity ,
the obtained information on the uv part of the spectral function suffers from a large uncertainty , as is also seen in the maximum entropy method ( mem ) analysis . on the other hand ,
the ir / im part of the spectral function @xmath13 is relatively stable and reliable .
we note that the gluon spectral function @xmath13 is divergent at @xmath195 , and the divergence structure is complicated and consists of two ingredients : a @xmath17-functional peak with a positive infinite residue and a negative wider power - damping peak . on the finite - volume lattice ,
these singularities are to be smeared , and @xmath13 is expected to take a finite value everywhere on @xmath173 .
on the lattice , we conjecture that the spectral function @xmath13 includes a narrow positive peak stemming from the @xmath17-function in the vicinity of @xmath206 and a wider negative peak near @xmath207 in the region of @xmath197 .    in this way
, the yukawa - type gluon propagator suggests an extremely anomalous spectral function of the gluon field in the landau gauge .
note that this framework gives an analytical and concrete expression for the gluon spectral function @xmath13 at the quantitative level .
actually , the resulting effective mass @xmath12 well describes the lattice result , as shown in fig.[fig : effmass ] .
the obtained gluon spectral function @xmath13 is negative almost everywhere , and includes a complicated divergence structure near the  anomalous threshold " , @xmath208 .
these anomalous features of the gluon spectral function may have some relation to the various nonperturbative qcd phenomena , such as the gluon confinement and the gluonic instability of the qcd vacuum , _
e.g. _ , gluon condensation , the savvidy vacuum @xcite , and the copenhagen vacuum @xcite . in any case , the yukawa - type gluon propagator @xmath96 , which models the landau - gauge gluon propagator , would be useful for the analytical and quantitative investigation of nonperturbative qcd .
we have studied the gluon propagator @xmath0 in the landau gauge in the infrared / intermediate region of @xmath209 1.0fm , which is relevant to the quark - hadron physics , in su(3 ) lattice qcd at @xmath1 = 5.7 , 5.8 , and 6.0 at the quenched level . from the gluon propagator analysis and the effective - mass plot ,
we have estimated the effective gluon mass of @xmath2mev in the infrared region of @xmath210fm .
the effective gluon mass exhibits a significant @xmath84-dependence in this region : it takes a smaller value in the smaller @xmath84-region .
we have also studied the functional form of the gluon propagator @xmath0 in lattice qcd . as a remarkable fact ,
the lattice qcd result of the landau - gauge gluon propagator @xmath4 is fairly well described by the yukawa - type form @xmath211 with the mass parameter @xmath6mev in the whole region of @xmath7fm in four - dimensional euclidean space - time .
this yukawa - type propagator corresponds to the new - type propagator @xmath160 in the momentum space , through the fourier transformation , and this also well describes the lattice qcd result of the gluon propagator @xmath212 in the momentum space .
as the application of the yukawa - type gluon propagator , we have derived the analytical expression of the zero - spatial - momentum propagator @xmath11 , and the effective mass @xmath12 .
the obtained analytical functions for @xmath11 and @xmath12 well reproduce the lattice qcd results , in particular the anomalous increasing behavior of @xmath12 .
we have found that the mass parameter @xmath14 of the yukawa - type gluon propagator directly corresponds to the effective gluon mass in the infrared region of @xmath151fm .
we have also derived the analytical expression of the spectral function @xmath13 of the gluon field , associated with the yukawa - type gluon propagator , using the inverse laplace transformation of the temporal propagator @xmath11 . as a remarkable fact ,
the obtained spectral function @xmath13 is negative - definite almost everywhere for @xmath197 , except for a positive @xmath17-functional peak with the residue of @xmath140 at @xmath213 .
the coexistence of negative- and positive - value regions of @xmath13 lead to the anomalous increasing behavior of the effective mass @xmath12 of gluons .
thus , the theoretical analysis with the yukawa - type gluon propagator gives a new analytical and quantitative method for the nonperturbative gluonic phenomena .
the yukawa function @xmath211 in the coordinate space is highly analytic , and this analyticity plays an important role in deriving the analytical expressions for @xmath11 , @xmath12 , and @xmath13 . on the other hand ,
the momentum - space propagator @xmath160 includes a singular cut stemming from the square root .
if this singularity is taken seriously , there would arise a problem , because analyticity of the green function is important in quantum field theories , _
e.g. _ , in the wick rotation converting between the euclidean space and the minkowski space .
of course , this yukawa - type function is an approximate form for the landau - gauge gluon propagator in the region of @xmath214fm .
there would be more regular and more thorough expression for the gluon propagator .    in this paper
, we have mainly considered the gluon propagator in the coordinate space instead of the momentum space , since the coordinate - space variable is more directly obtained in lattice qcd . for the confined particles ,
however , the reason to use the momentum representation would be less clear , compared with ordinary particles . in the ordinary particles , the momentum representation of the green function is clearly useful to express the pole structure , to distinguish the on - shell and off - shell states , and so on .
however , for the confined field , there is no on - shell state , _
i.e. _ , no physical asymptotic state , so that there is no definite reason to use the momentum representation , besides the total momentum conservation .
indeed , it is difficult to image the non - zero momentum - space propagator with no pole , and the coordinate - space representation may be more convenient for some description of the confined field , similar to potential problems in quantum mechanics , where the coordinate - space wave - function is convenient .
the yukawa - type gluon propagator @xmath96 includes @xmath143 as the main reduction factor in the infrared region , and @xmath215gev can be regarded as the infrared effective gluon mass . in terms of the infrared reduction , a simple  constituent gluon picture " may be approximately obtained as @xmath216 for the glueball mass @xmath217 . in general , the glueball mass @xmath217 can be estimated from the infrared reduction of the glueball correlator @xmath218 with the glueball operator , _
e.g. _ , @xmath219 for the scalar glueball , with the field strength tensor @xmath220 @xcite . by the wick contraction , @xmath221 can be expressed by some combination of the gluon green functions . in the framework with the yukawa - type gluon propagator
, the leading reduction term of @xmath221 can be expressed with some derivative of @xmath222 for large @xmath84 .
then , @xmath221 includes @xmath223 as one of the main infrared reduction factors , and the lowest - glueball mass is roughly estimated as @xmath224 1.3gev , by neglecting all the prefactor and higher - order terms , which express the interactions between gluons . in spite of the crude estimate , this value gives the same order of the lowest - glueball mass of about 1.5gev obtained in lattice qcd @xcite .    in this subject
, there remain difficult problems related to the confinement mechanism , the infinite - volume limit , the gribov copies , and the gauge dependence .
for example , it is a highly difficult mathematical problem to find out the precise description of the physical hadronic states like glueballs in terms of the confined particles , quarks and gluons . as other example
, zwanziger s theorem @xcite is derived from the argument of the gribov horizon : all connected gluon correlation functions including the gluon propagator must vanish at zero momentum in the infinite - volume limit .
however , our lattice qcd results and the yukawa - type propagator @xmath75 indicate an infrared non - vanishing property of the gluon propagator as @xmath225
. this may be due to the absence of an infinite - volume effect in the deep - ir region , as is conjectured by analytical studies @xcite .
for this problem , it is desired to clarify the deep - ir behavior of the gluon propagator , and recent huge - volume lattice studies @xcite and a recent analytical study based on the schwinge - dyson equation @xcite also indicate the infrared non - vanishing property of the gluon propagator . as for the gribov - copy problem , it is reported that the gribov - copy effect is quantitatively rather small in the actual lattice qcd calculation for the landau - gauge gluon propagator @xcite , while the ghost propagator slightly suffers from it .
however , this is a fundamental problem in qcd , and it would be serious in the argument of the large - volume limit , so that it is also desired to remove the gribov copies and to extract the fundamental moduli region .    finally , we discuss the yukawa - type gluon propagation and a possible dimensional reduction due to the stochastic behavior of the gluon field in the infrared region . as shown in this paper , the landau - gauge gluon propagator is well described by the yukawa function in _ four_-dimensional euclidean space - time . however , the yukawa function @xmath5 is a natural form in _ three_-dimensional euclidean space - time , since it is obtained by the three - dimensional fourier transformation of the ordinary massive propagator @xmath226 .
in fact , the yukawa - type propagator has a  three - dimensional " property . in this sense , as an interesting possibility , we propose to interpret this yukawa - type behavior of the gluon propagation as an  effective reduction of the space - time dimension " .
such a  dimensional reduction " sometimes occurs in stochastic systems , as parisi and sourlas pointed out for the spin system in a random magnetic field @xcite . on the infrared dominant diagrams ,
the @xmath227-dimensional system coupled to the gaussian - random external field is equivalent to the @xmath228-dimensional system without the external field .
in fact , the space - time dimension of the theory is apparently reduced by two . for the system coupled to the gaussian - random external source ,
the dimensional reduction is universal and is associated with a hidden supersymmetry : in the superspace formalism with @xmath229 , the integration over two grassmann - variables ( @xmath230 ) reduces the space - time coordinates @xmath231 by two @xcite .
we note that the gluon propagation in the qcd vacuum resembles the situation of the system coupled to the stochastic external field .
in fact , as is indicated by a large positive value of the gluon condensate @xmath232 in the minkowski space , the qcd vacuum is filled with a strong color - magnetic field @xcite , which can contribute spontaneous chiral - symmetry breaking @xcite , and the color - magnetic field is considered to be highly random at an infrared scale @xcite . since gluons interact each other , the propagating gluon is violently scattered by the other gluon fields which are randomly condensed in the qcd vacuum at the infrared scale , as schematically shown in fig.[figgluonprogagation ] .     a schematic figure for a propagating gluon in the qcd vacuum .
the qcd vacuum is filled with color - magnetic fields which are stochastic at an infrared scale , and the gluon propagates in the random color - magnetic fields . ]    actually in the infrared region , the gluon field shows a strong randomness due to the strong interaction , and this infrared strong randomness is considered to be responsible for color confinement , as is indicated in strong - coupling lattice qcd @xcite . in the chiral random matrix theory for qcd @xcite , the infrared randomness of gluons is taken account as an essence of qcd in a simplified manner , and the gluon field is replaced by a gaussian - random external field coupled to the quark field .
notice that there are two kinds of randomness in the gluon field : one is a completely - random gauge degree of freedom , which is fake , and the other is a net physical gluonic degree of freedom , which is not completely but highly random at the infrared scale . for the argument of physical randomness , these two concepts should be clearly divided , since the gauge degree of freedom is just a fake . even after the removal of the fake gauge degree of freedom by gauge fixing , the gluon field exhibits a strong randomness @xcite accompanying a quite large fluctuation at the infrared scale .    as a generalization of the parisi - sourlas mechanism
, we conjecture that the infrared structure of a theory in the presence of the quasi - random external field in higher - dimensional space - time has a similarity to the theory without the external field in lower - dimensional space - time . from this point of view
, the yukawa - type behavior of gluon propagation may indicate an  effective reduction of the space - time dimension " by one , due to the stochastic interaction between the propagating gluon and the other gluon fields in the qcd vacuum , of which net physical fluctuation is highly random at the infrared scale .    in this paper , the yukawa - type gluon propagator
is obtained phenomenologically from lattice qcd results , but we expect some deeper theoretical reasons for the yukawa - type propagation , which may be an effective reduction of the space - time dimension , due to the stochastic behavior of the infrared gluon field . in any case , the yukawa - type gluon propagator would provide a new analytical framework for the study of nonperturbative qcd .    h.  s. is grateful to professor j.m .
cornwall for his useful suggestions on the dynamical gluon mass .
h.  s. is supported in part by the grant for scientific research [ ( c ) no .
19540287 ] from the ministry of education , culture , science and technology ( mext ) of japan .
this work is supported by the global coe program ,  the next generation of physics , spun from universality and emergence " at kyoto university .
the lattice qcd calculations have been done on nec - sx8 at osaka university .
in this appendix , we derive several fourier transformations used in sec.[sec : functionalforms ] .          we here calculate the fourier transformation of the massive propagator @xmath226 .
first , we consider the fourier integral of @xmath226 in the four - dimensional euclidean space - time , @xmath243 by rotating the coordinate , we set @xmath244 without loss of generality .
then , the integral is expressed as @xmath245 with @xmath246 . with the three - dimensional polar coordinate of @xmath247 ,
@xmath248 is written as @xmath249 with @xmath250 and @xmath251 . using the integral representation of the modified bessel function , @xmath252 we obtain the fourier transformation formula , @xmath253      next ,
we consider the fourier integral of @xmath226 in the two - dimensional euclidean space - time , @xmath254 by rotating the coordinate , we set @xmath255 without loss of generality , and integrate @xmath256 as @xmath257 using mehler s integral representation of the modified bessel function , @xmath258 we obtain the fourier transformation formula , @xmath259
in the deep - ir region , there has been reported to appear some deviation on the momentum - space gluon propagator @xmath75 between small - size lattice data and huge - volume lattice data @xcite . in this appendix
, we demonstrate that the deviation in the deep - ir region does not affect the yukawa - type behavior of the coordinate - space gluon propagator @xmath82 in the ir / im region of @xmath33fm .
figure [ fig : hugelattice ] shows the scalar - type gluon propagator @xmath81 in the recent lattice - qcd calculation with a huge volume @xmath265 at @xmath266 , taken from ref.@xcite . here ,
a renormalization constant is multiplied for the huge - volume lattice data so as to adjust them to the renormalization condition ( [ renormcond ] ) at @xmath267=4gev . in the momentum space ,
the true gluon propagator @xmath75 turns out to take a saturated value smaller than the yukawa - type propagator @xmath127 in the deep - ir region of @xmath128gev . in other words , @xmath129gev is the lower bound on the applicability of the yukawa - type propagator @xmath127 to the gluon propagator .     the infrared behavior of the gluon propagator @xmath75 .
the triangle symbols denote recent huge - volume lattice data taken from ref.@xcite .
the solid line denotes the yukawa - type propagator @xmath117 , and the dashed line the deep - ir - corrected propagator @xmath268 with @xmath269gev . ]
taking account of the deviation in the deep - ir region , we define the deep - ir - corrected momentum - space propagator as @xmath270 with the ir - saturation momentum of @xmath271 .
this value of @xmath272 is determined so as to consist with the huge - volume lattice result in the deep - ir region . using this deep - ir - corrected propagator @xmath273
, we calculate the deep - ir - corrected coordinate - space propagator by the fourier transformation as @xmath274    in fig.[fig : hugelatticecorr ] , we show the yukawa - type propagator @xmath96 and this deep - ir - corrected propagator @xmath275 , together with the lattice qcd data .
the difference between @xmath96 and @xmath275 is fairly small for @xmath7 fm .
the yukawa - type propagator @xmath96 ( solid line ) , and deep - ir - corrected propagator @xmath275 ( dashed - line ) , together with the lattice data .
the difference between them is fairly small in the ir / im region of @xmath33fm . ]
99 y.  nambu , proc . of _ preludes theoretical
physics , in honor of v.f .
weisskopf _ ( north - holland , 1966 )
han and y.  nambu , phys . rev .
* 139 * , b1006 ( 1965 ) .
gross and f.  wilczek , phys .
lett . * 30 * , 1343 ( 1973 ) ; h.d .
politzer , phys .
lett . * 30 * , 1346 ( 1973 ) .
feynman , proc . of _ high energy collision of hadrons _
( stony brook , n.y . , 1969 ) ;
bjorken and e.a .
paschos , phys . rev . * 185 * , 1975 ( 1969 ) .
w.  greiner , s.  schramm , and e.  stein , _ quantum chromodynamics _ ,
( springer , 2007 ) , and its related references .
wilson , phys .
d*10 * , 2445 ( 1974 ) ; j.b .
kogut and l.  susskind , phys .
rev . d*11 * , 395 ( 1975 ) .
m.  creutz , phys .
rev . lett .
* 43 * , 553 ( 1979 ) ; phys .
d * 21 * , 2308 ( 1980 ) . as a recent textbook , _
e.g. _ , h.j .
rothe , _ lattice gauge theories : an introduction _ , 3rd edition , ( world scientific , 2005 ) , and its related references .
for example ,  yang - mills and mass gap " is nominated in millennium problems by clay mathematics institute .
see the website : http://www.claymath.org/millennium/ y.  nambu , g.  jona - lasinio , phys .
rev . * 122 * , 345 ( 1961 ) ; phys . rev . *
124 * , 246 ( 1961 ) .
miransky , _ dynamical symmetry breaking in quantum field theories _ ( world scientific , 1993 ) ; + k.  higashijima , prog .
. suppl . *
104 * , 1 ( 1991 ) . as an instructive review ,
e.g. _ , j.e .
mandula , phys . rept .
* 315 * , 273 ( 1999 ) , and its references . j.m .
cornwall , phys .
rev . d*26 * , 1453 ( 1982 ) ; phys .
rev . d*76 * , 025012 ( 2007 ) , and references therein .
gribov , nucl . phys . *
b139 * , 1 ( 1978 ) .
d.  zwanziger , nucl . phys . *
b364 * , 127 ( 1991 ) ; phys .
b*257 * , 168 ( 1991 ) ; nucl . phys .
* b378 * , 525 ( 1992 ) ; phys .
rev . d*65 * , 094039 ( 2002 ) ; phys .
d*69 * , 016002 ( 2004 ) .
m.  stingl , phys .
rev . d*34 * , 3863 ( 1986 ) ; _ erratum - ibid .
_ d*36 * , 652 ( 1987 ) ; u.  hbel , r. konning , h.g .
reusch , m. stingl , and s. wigard , z. phys .
a * 336 * , 423 ( 1990 ) .
l. von  smekal , a.  hauck , and r.  alkofer , phys .
* 79 * , 3591 ( 1997 ) ; r.  alkofer and l. von  smekal , phys . rept . * 353 * , 281 ( 2001 ) , and references therein .
aguilar , d.  binosi , and j.  papavassiliou , phys .
rev . d*78 * , 025010 ( 2008 ) .
kondo , phys .
b*678 * , 322 ( 2009 ) .
mandula and m.  ogilvie , phys .
b*185 * , 127 ( 1987 ) .
r. gupta , g. guralnik , g. kilcup , a. patel , s.r .
sharpe , and t.  warnock , phys .
d * 36 * , 2813 ( 1987 ) .
bernard , c.  parrinello , and a.  soni , phys .
* , 1585 ( 1994 ) .
p.  marenzoni , g.  martinelli , n.  stella , and m.  testa , phys .
b*318 * , 511 ( 1993 ) ; p.  marenzoni , g.  martinelli , and n.  stella , nucl . phys . *
b455 * , 339 ( 1995 ) .
a.  cucchieri , nucl .
* b508 * , 353 ( 1997 ) ; nucl . phys . * b521 * , 365 ( 1998 ) . d.b .
leinweber , j.i .
skullerud , a.g .
williams , and c.  parrinello , phys .
rev . d*58 * , 031501 ( 1998 ) ; phys .
rev . d*60 * , 094507 ( 1999 ) , _ erratum - ibid .
_ d*61 * , 079901 ( 2000 ) .
d.  becirevic , ph .
boucaud , j.p .
leroy , j.  micheli , o.  pne , j.  rodriguez - quintero , c.  roiesnel , phys .
d*60 * , 094509 ( 1999 ) ; phys . rev . d*61 * , 114508 ( 2000 ) .
bonnet , p.o .
bowman , d.b .
leinweber , and a.g .
williams , phys .
rev . d*62 * , 051501 ( 2000 ) ; f.d.r .
bonnet , p.o .
bowman , d.b .
leinweber , a.g .
williams , and j.m .
zonotti , phys .
rev . d*64 * , 034501 ( 2001 ) .
k.  langfeld , h.  reinhardt , and j.  gattnar , nucl .
b621 * , 131 ( 2002 ) .
s.  furui and h.  nakajima , phys .
rev . d*69 * , 074505 ( 2004 ) .
bowman , u.m .
heller , d.b .
leinweber , m.b .
parappilly , and a.g .
williams , phys .
rev . d*70 * , 034509 ( 2004 ) ; phys .
rev . d*76 * , 094505 ( 2007 ) .
a.  sternbeck , e .-
ilgenfritz , m.  mueller - preussker , and a.  schiller , phys .
d*72 * , 014507 ( 2005 ) .
bogolubsky , e .-
ilgenfritz , m.  mller - preussker , and a.  sternbeck , phys .
b*676 * , 69 ( 2009 ) ; pos(lat2007 ) , 290 ( 2007 ) .
p.j .  silva and o.  oliveira ,
phys . rev . d*74 * , 034513 ( 2006 ) ; a.  cucchieri , t.  mendes , o.  oliveira , and p.j .
silva , phys .
d*76 * , 114507 ( 2007 ) .
a.  cucchieri and t.  mendes , pos(lat2007 ) , 297 ( 2007 ) ; phys . rev .
lett . * 100 * 241601 ( 2008 ) .
j.  greensite and s.  olejnik , phys .
rev . d*67 * , 094503 ( 2003 ) ; j.  greensite , s.  olejnik , and d.  zwanziger , phys . rev . d*69 * , 074506 ( 2004 ) , and their references .
a.  cucchieri and d.  zwanziger , phys .
b*524 * , 123 ( 2002 ) , and references therein .
k.  amemiya and h.  suganuma , phys .
rev . d*60 * , 114509 ( 1999 ) ; h.  suganuma , k.  amemiya , a.  tanaka , and h.  ichie , nucl .
* a670 * , 40 ( 2000 ) ; h.  suganuma , k.  amemiya , h.  ichie , n.  ishii , h.  matsufuru , t.t .
* b * ( proc . suppl . )
* 106 - 107 * , 679 ( 2002 ) .
kondo , phys .
rev . d*57 * , 7467 ( 1998 ) ; phys .
rev . d*58 * , 105016 ( 1998 ) ; phys .
rev . d*58 * , 105019 ( 1998 ) .
s.  mandelstam , phys .
rev . d*20 * , 3223 ( 1979 ) . m. baker , j. s. ball , and f. zachariasen , nucl .
b186 * , 531 , 560 ( 1981 ) ; phys . rept . * 209 * , 73 ( 1991 ) . h.  suganuma , s.  sasaki and h.  toki , nucl . phys . *
b435 * , 207 ( 1995 ) ; prog .
. phys . * 94 * , 373 ( 1995 ) .
m.  beneke and v.m .
braun , nucl . phys . *
b426 * , 301 ( 1994 ) ; m. beneke , phys . rept . * 317 * , 1 ( 1999 ) .
n.  brambilla , y.  sumino , and a.  vairo , phys .
b*513 * , 381 ( 2001 ) ; y.  sumino , phys .
b*571 * , 173 ( 2003 ) .
n.  ishii , h.  suganuma , and h.  matsufuru , phys .
rev . d*66 * , 094506 ( 2002 ) ; phys .
rev . d*66 * , 014507 ( 2002 ) .
bernard , phys .
b*108 * , 431 ( 1982 ) ; nucl . phys . *
b219 * , 341 ( 1983 ) .
f.  halzen , g.  krein , and a.a .
natale , phys .
rev . d*47 * , 295 ( 1993 ) .
cornwall and a. soni , phys .
b*120 * , 431 ( 1983 ) .
hou , c.s  luo , and g.g .
wong , phys .
* , 014028 ( 2001 ) .
g.  parisi and r.  petronzio , phys .
b*94 * , 51 ( 1980 ) . a.  yamamoto and h.  suganuma , phys .
* 101 * , 241601 ( 2008 ) ; phys . rev . d*79 * , 054504 ( 2009 ) .
h.  suganuma , t.t .
takahashi , and h.  ichie , _ color confinement and hadrons in quantum chromodynamics _ , ( world scientific , 2004 ) 249 ; t.t .
takahashi , h.  suganuma , y.  nemoto , and h.  matsufuru , phys .
rev . d*65 * , 114509 ( 2002 ) ; t.t .
takahashi , h.matsufuru , y.nemoto , and h.  suganuma , phys .
lett .  * 86 * , 18 ( 2001 ) .
h.  iida , m.  oka , and h.  suganuma , eur .
j. * a23 * , 305 ( 2005 ) . c.  itzykson and j.  zuber , _ quantum field theory _
( mcgraw - hill inc .
, 1980 ) . f.  de soto and c.  roiesnel , jhep * 09 * , 007 ( 2007 ) .
savvidy , phys .
b*71 * , 133 ( 1977 ) .
nielsen and p.  olesen , nucl .
b144 * , 376 ( 1978 ) ; j.  ambjorn and p.  olesen , nucl . phys . *
b170 * , 60 ( 1980 ) .
g.  parisi and n.  sourlas , phys .
lett . * 43 * , 744 ( 1979 ) .
h.  suganuma and t.  tatsumi , ann . phys . * 208 * , 470 ( 1991 ) ; prog . theor . phys .
* 90 * , 379 ( 1993 ) .
verbaarschot and t.  wettig , ann .
nucl . part .
sci . * 50 * , 343 ( 2000 ) , and its references .
h.  ichie and h.  suganuma , nucl . phys . *
b548 * , 365 ( 1999 ) ; nucl . phys . *
b574 * , 70 ( 2000 ) .