the standard thin accretion disk model @xcite has been widely applied to x - ray binaries and active galactic nuclei . due to the basic assumption of the energy balance between the viscous heating and the radiative cooling ,
such a model was known to be invalid for super - eddington accretion case , where the advective cooling is probably significant . instead
, the slim disk model @xcite was introduced to describe super - eddington accretion disks .
however , there exists some conflict between the theory and the observation .
the theory predicts that the advection is negligible for @xmath1 ( e.g. , * ? ? ?
* ; * ? ? ?
* ) , where @xmath2 is the eddington luminosity , which indicates that the standard disk model should be valid up to @xmath2 . on the contrary
, some works on the black hole spin measurement showed that the standard disk model is likely to be inaccurate for @xmath3 ( e.g. , * ? ? ? * ) .
moreover , even the recent general model for optically thick disks ( e.g. , * ? ? ?
* ; * ? ? ?
* ) , which unifies the standard thin disk and the slim disk , could not help to obtain a self - consistent spin parameter for @xmath3 ( e.g. , * ? ? ?
* ) . in our opinion
, the above conflict may be resolved if the vertical structure is well incorporated .
most previous works on accretion disks focused on the radial structure in cylindrical coordinates ( @xmath4 , @xmath5 , @xmath6 ) . for the vertical structure , however , a simple well - known relationship
 @xmath7 " or  @xmath8 " was widely adopted , where @xmath9 is the half - height of the disk , @xmath10 is the sound speed , and @xmath11 is the keplerian angular velocity .
such a relationship comes from the vertical hydrostatic equilibrium with two additional assumptions .
one is the approximation of gravitational potential : @xmath12 , and the other is a one - zone approximation or a polytropic relation @xmath13 in the vertical direction ( e.g. , * ? ? ?
* ) , where @xmath14 is the total ( gas plus radiation ) pressure and @xmath15 is the density .
obviously , the above assumptions work well for geometrically thin disks , but may be inaccurate for the mass accretion rate @xmath16 approaching the eddington one @xmath17 , for which the disk is probably not thin .
consequently , the relationship  @xmath8 " may be invalid for @xmath18 .    without the potential approximation , our two previous works investigated the geometrical thickness of accretion disks and the validity of the relationship  @xmath8 " .
@xcite adopted the explicit gravitational potential in cylindrical coordinates and found that the above relationship is inaccurate for @xmath18 , and therefore the disk can be geometrically thick .
@xcite took spherical coordinates to avoid the approximation of gravitational potential , and found that an advection - dominated accretion disk is likely to be quite thick .
in these two works , however , the polytropic relation is still adopted in the vertical direction , which takes the place of the energy balance per unit volume between the viscous heating and the advective cooling plus the radiative cooling .
the validity of such a polytropic relation , however , remains questionable , in particular for large @xmath16 due to dominant radiation pressure .
the purpose of this paper is to revisit the vertical structure of radiation pressure - supported disks by taking into account the local energy balance and to study the variation of energy advection with mass accretion rates .
the paper is organized as follows .
equations and boundary conditions are derived in section  2 .
a global view of the solutions in the @xmath19-@xmath20 diagram is presented in section  3 . for a typical radius @xmath21 , the vertical structure and the energy advection
are investigated in section  4 .
the two - dimensional solutions and the convective stability are studied in section  5 .
summary and discussion are made in section  6 .
we consider a steady state axisymmetric accretion disk in spherical coordinates ( @xmath20 , @xmath22 , @xmath5 ) and use the newtonian potential , @xmath23 , where @xmath24 is the black hole mass .
following @xcite , we assume @xmath25 for simplicity , which means a hydrostatic equilibrium in the @xmath22 direction .
simulations , ( e.g. , * ? ? ?
* figure  3 ) , however , revealed that @xmath26 will be significant for extremely high accretion rates such as @xmath27 .
as shown in the following sections , our solutions mainly correspond to @xmath16 around @xmath17 . for such accretion rates ,
the validity of @xmath25 remains a question .
the basic equations of continuity and momentum take the forms ( e.g. , * ? ? ?
* ) : @xmath28 @xmath29 @xmath30 @xmath31 where @xmath32 and @xmath33 are respectively the radial and azimuthal velocity , @xmath34 and @xmath35 are respectively the radial and vertical radiation flux , @xmath36 is the gas pressure , @xmath37 is the opacity of electron scattering , and @xmath38 is the @xmath39 component of the viscous stress tensor , @xmath40 . following the spirit of @xmath41 stress prescription , we assume the kinematic viscosity coefficient @xmath42 , where @xmath10 is the sound speed defined below ( equation  ( 7 ) ) , and @xmath43 is the keplerian velocity .    we would stress that , even though the @xmath41 stress prescription has been widely adopted for theoretical studies , simulations of magnetorotational turbulence have shown that the stress does not well scale locally with the pressure .
for instance , the simulations on thin disks by a shearing box showed that , the time and box - averaged results are likely to support that the stress is proportional to the thermal ( gas plus radiation ) pressure ( e.g. , * ? ? ?
* figure  3 ) .
however , figure  16 of @xcite shows that the maximal thermal pressure is located on the equatorial plane , whereas figure  11 shows that the maximal stress is obviously not at the same place .
these two figures reveal that the stress is not proportional to the pressure locally . in the present study , for simplicity
, we will keep the local @xmath41 stress prescription for numerical calculation , which is a weak point of this work .
the energy equation including gas and radiation is written as ( e.g. , * ? ? ?
* ) @xmath44 = -p { { \mbox{\boldmath $ \nabla$ } } } \cdot { { \mbox{\boldmath $ v$ } } } - { { \mbox{\boldmath $ \nabla$ } } } { { \mbox{\boldmath $ v$ } } } { { \mbox{\boldmath $ : $ } } } { { \mbox{\boldmath $ p$ } } } - { { \mbox{\boldmath $ \nabla$ } } } \cdot { { \mbox{\boldmath $ f$ } } } + \phi_{\rm vis } \ , \ ] ] where @xmath45 and @xmath46 are the internal energy density of the gas and the radiation , respectively .
@xmath47 is the radiation pressure tensor , @xmath48 is the viscous dissipative function , and the radiation flux @xmath49 is expressed as @xmath50 in this work , we focus on the region inside the photosphere , so we can take the well - known eddington approximation , i.e. , @xmath51 and the eddington tensor @xmath52 .
since we only study the radiation pressure - supported disks , the gas pressure @xmath36 and the gas internal energy density @xmath45 will be dropped in our calculation . in order to avoid directly solving the partial differential equations , some assumptions on the radial derivatives ( @xmath53 ) are required . following the spirit of self - similar assumptions ( e.g. , * ? ? ?
* ; * ? ? ?
* ) , we adopt the following radial derivatives for @xmath10 and @xmath46 : @xmath54 where the sound speed @xmath10 is defined as @xmath55 based on the above two radial derivatives , the following four derivatives can be inferred from equations  ( 1)-(7 ) : @xmath56    with all the above derivatives , we can remove the
@xmath53 " terms in equations ( 2 ) and ( 4)-(6 ) , and the following equations are then obtained from equations ( 2)-(6 ) : @xmath57 @xmath58 @xmath59 @xmath60 @xmath61 @xmath62    the seven equations , equations  ( 7)-(13 ) , enable us to solve for the seven variables : @xmath32 , @xmath33 , @xmath10 , @xmath15 , @xmath46 , @xmath34 , and @xmath35 . there are two differential equations in this system .
in addition , the position of the surface is unknown .
thus , totally three boundary conditions are required to determine a unique solution .
an obvious boundary condition on the equatorial plane is @xmath63 .
however , this condition is not applicable for numerical calculation since it is automatically matched as indicated by equation  ( 9 ) .
combining equations ( 9 ) and ( 11 ) we can derive the following equation : @xmath64 an alternative boundary condition on the equatorial plane is then obtained from the above equation ( the left - hand side is zero thus the right - hand side should also be zero ) : @xmath65 the second boundary condition is a definition of the surface .
we define the photosphere as the position above which the optical depth is around unity .
the condition can be written as @xmath66 where @xmath67 ( @xmath68 ) is the polar angle of the photosphere .
the third condition is related to the mass accretion rate : @xmath69
in our calculation we set @xmath70 , @xmath71 @xmath72g@xmath73 , and @xmath74 , where the value of @xmath41 is taken from recent simulations @xcite .
the eddington accretion rate is expressed as @xmath75 , where @xmath76 is a radiative efficiency of the flow .
we choose @xmath77 since it is comparable to the schwarzschild black hole efficiency of 0.057 .
the dimensionless accretion rate is defined as @xmath78 .    with the equations and boundary conditions in section 2
, we can numerically derive the @xmath22-direction distribution of physical quantities for a given @xmath19 at a certain radius @xmath20 .
the radiation pressure - supported disk solutions in the @xmath19-@xmath20 diagram are shown in figure  1 , where @xmath79 is the schwarzschild radius .
the parameter space is divided into three regions by two parallel solid lines , roughly with @xmath80 .
the region above the upper solid line is denoted by  outflow " , where we can not find solutions .
no solution exists probably due to the assumption of @xmath81 in advance . in our view , a real flow located in this region may have @xmath82 and the inflow accretion rate may decrease inward . the physical understanding could be that , for high accretion rates and particularly for the inner radii , the viscous dissipation may be sufficiently large such that the radiation pressure is too strong to be balanced by the gravitational force .
thus , outflows may be driven by the radiation pressure and the inflow @xmath19 drops inward .
on the other hand , simulations of supercritical accretion flows ( e.g. , * ? ? ? * figure  6 ) showed that the inflow accretion rate roughly follows the @xmath80 relationship for @xmath83 at @xmath84 ( corresponding to @xmath85 due to the definition of @xmath17 with @xmath77 ) .
the slope of the upper solid line in figure  1 , which may be regarded as maximal accretion rates due to our calculation , agrees well with the slope in the above simulations .
the region under the lower solid line is denoted by  gas pressure " , where no solution is found either . in our understanding , it is probably because the gas pressure can not be ignored in this region , which may be in conflict with the radiation pressure - supported assumption . we would point out that , the lower solid line in this diagram is higher than the well - known line which separates the inner and middle regions of standard thin disks @xcite .
the reason is that , the gas and radiation pressure are comparable for the latter , whereas the radiation pressure - supported disk may require the accretion rate to be higher such that the radiation pressure sufficiently dominate over the gas pressure , and therefore the effect of gas pressure on the vertical structure can be completely ignored .
the region between the two solid line , denoted by  radiation pressure " , which means that the radiation pressure is completely dominated , corresponds to the solutions of our main interest in this work . in section  4
, we will focus on the vertical structure and the energy advection at a typical radius , @xmath86 , as indicated by the vertical dashed line in figure  1 . in section  5
, we will study the two - dimensional solutions for a typical accretion rate @xmath87 in the range @xmath88 and @xmath89 , as indicated by the horizontal dot - dashed line .
in addition , we would point out that for inner radii such as @xmath90 , the two solid lines in figure  1 may deviate from a real black hole accretion system due to the newtonian potential used in this work .
in this section , we will focus on the solutions at a typical radius @xmath91 .
figure  2 shows the vertical structure of the disk with @xmath87 . in figure  2(a ) , the dot - dashed , dotted , solid , and dashed lines show the vertical distribution of the dimensionless density ( @xmath92 ) , radial velocity ( @xmath93 ) , azimuthal velocity ( @xmath94 ) , and sound speed ( @xmath95 ) , respectively , where @xmath96 is the density on the equatorial plane .
it is seen that @xmath15 significantly decreases , whereas @xmath10 and @xmath97 increases , from the equatorial plane to the surface . in figure  2(b ) , the solid line shows the variation of @xmath98 ( defined in equation  ( 15 ) ) , where the photosphere ( @xmath99 ) is located at @xmath100 , quite close to the polar axis .
the disk seems to be extremely thick according to the position of the photosphere .
however , the profile of @xmath15 implies that most of the accreted matter exists in a moderate range around the equatorial plane , such as @xmath101 , which is more clear in figures  4 and 5 ( discussed below ) .
the dashed line shows the variation of @xmath102 .
there exists @xmath103 for the whole solution , which indicates that the eddington approximation is valid and the solution is therefore self - consistent .    with a more general viscosity
, @xcite studied a geometrically thick , radiation pressure - supported model for supercritical accretion disks .
they showed that there exists a narrow empty funnel along the rotation axis with a half - opening angle @xmath104 .
as seen in our figure  2(a ) , the density drops sharply close to the photosphere , thus a nearly empty funnel also seems to exist in our model .
the difference is that , the disk surface in @xcite is the position where some physical quantities such as @xmath32 diverges , whereas our surface is defined as the position where equation  ( 15 ) is matched , and no divergence appears in our solutions .    for a real disk with @xmath87 , the photosphere may exist between @xmath105 and the present result ( @xmath106 ) .
our argument is as follows .
there are two possible reasons that may cause the present photosphere quite close to the polar axis .
one is that we have ignored the radiation force from one side ( e.g. , @xmath107 and @xmath108 ) to the other ( e.g. , @xmath107 and @xmath109 ) .
the other reason is that we consider only the @xmath39 component of the stress tensor , which may cause inaccurate results for small @xmath22 , such as strong shearing of the angular velocity @xmath110 in the vertical direction , where @xmath111 .
nevertheless , the @xmath112 component assumption may work well for moderate @xmath22 : @xmath101 .
the fact that the surface condition could not be matched in the range @xmath101 indicates that the half - opening angle of the disk ( @xmath113 ) is likely to be larger than @xmath114 .
the profile of @xmath10 in figure  2(a ) is quite different from that in the previous works with a vertical polytropic assumption ( e.g. , * ? ? ?
* ; * ? ? ?
. under the polytropic relation @xmath13 ( normally @xmath115 , and for radiation pressure - dominated case , @xmath14 can be replaced by @xmath116 ) , @xmath10 will decrease continuously from the equatorial plane to the surface .
the reason for the opposite behavior of @xmath10 , as implied in figure  2(a ) , is that @xmath15 drops faster than @xmath46 from the equatorial plane to the surface .
figure  3 shows the variation of the quantity @xmath117 with @xmath22 for @xmath118 ( dashed line ) , @xmath87 ( solid line ) , and @xmath119 ( dotted line ) .
if the polytropic relation works well , @xmath117 should be a constant of @xmath120 .
it is clearly shown in figure  3 that , however , @xmath117 varies significantly with @xmath22 rather than being a constant .
more importantly , @xmath121 indicates that @xmath122 is negative thus unacceptable .
we therefore argue that the polytropic relation should be unsuitable for describing the vertical structure of radiation pressure - supported disks .
moreover , since the energy advection is relevant to @xmath10 , we may expect essentially different results on the strength of advection .
figure  4 shows the variation of the vertically averaged advection factor @xmath123 with the mass accretion rate @xmath19 , where @xmath123 is defined as @xmath124 .
the quantities @xmath125 and @xmath126 are expressed as follows : @xmath127 @xmath128 where @xmath129 and @xmath130 are respectively the advective cooling rate and the viscous heating rate per unit volume , as implied by the left - hand side of equation  ( 11 ) .
the solid line in figure  4 corresponds to the total accretion rate integrating from @xmath67 to @xmath131 , as shown by equation  ( 16 ) , whereas the dashed line corresponds to the specific accretion rate integrating from @xmath132 to @xmath133 , i.e. , @xmath134 the reason why we calculate for @xmath135 is that the @xmath112 stress assumption may work well for @xmath101 . as shown by the horizontal range of the solid and dashed lines , most of the accreted matter exists in this specific range , e.g. , @xmath136 corresponding to @xmath87 .
the figure also shows that @xmath123 rapidly increases with increasing @xmath19 in the range @xmath137 .
more importantly , the value of @xmath123 ( @xmath138 ) indicates that the energy advection is significant even for sub - eddington accretion disks .
such a result is quite different from the previous one , where advection was found to be significant only for super - eddington accretion case .
@xcite introduced an elegant formula to describe the @xmath139 relationship based on their numerical solutions under the well - known paczyski - wiita potential @xcite .
their equations  ( 15)-(19 ) imply that , for the position @xmath91 , advection is negligible for @xmath140 ( @xmath141 ) .
for the whole disk , advection is negligible for @xmath142 ( @xmath143 ) .
such a critical @xmath16 for the whole disk was confirmed by some recent global solutions under the general relativity .
figure  4.11 of @xcite shows that advection is negligible for @xmath1 for any spin parameter @xmath144 . for @xmath145 ,
i.e. , the schwarzschild black hole , the critical @xmath16 is just around @xmath17 . in our opinion ,
the different results on the advection between the above two works and ours are related to the different approach to describing the vertical structure .
a significant difference is that @xcite and @xcite chose the cylindrical coordinates whereas we adopt the spherical coordinates . of course
, the final results should not depend on the coordinates used .
however , as pointed out by @xcite , there are some interesting differences between the equations written in cylindrical and spherical coordinates .
there is no centrifugal force in the @xmath6 direction in cylindrical coordinates , whereas there is no gravitational force in the @xmath22 direction in spherical coordinates .
@xcite claimed that it is exactly this property that makes the spherical coordinates much better adapted for describing the flow near the black hole horizon .
here , we argue that the spherical coordinates should be more suitable for describing geometrically thick disks as follows . in cylindrical coordinates , in the @xmath6 direction , whether with a polytropic relation between the pressure and the density @xcite , or with the local energy balance @xcite , an approximation for the gravitational force , i.e. , @xmath146 , was adopted for describing the vertical structure .
such an approximation will probably be invalid for @xmath147 .
in particular for @xmath148 , the approximate force goes to infinity whereas the real force ought to vanish .
thus , the cylindrical coordinates seem unsuitable for studying geometrically thick disks . in other words , a geometrically thick disk solution in cylindrical coordinates may not be self - consistent .
@xcite limited their solutions by @xmath149 probably due to this reason .
as shown by their figure  10 , the maximal value of @xmath150 for @xmath151 is @xmath152 .
for higher @xmath16 , the value of @xmath150 will be even larger thus the solution based on the approximate force may be inaccurate . on the contrary , in spherical coordinates , there is no need to make approximation for the gravitational force .
the centrifugal force in the @xmath22 direction , which takes the place of the @xmath6-direction gravitational force in cylindrical coordinates , is derived in this work by solving the vertical differential equations .
thus , our approach to the vertical structure seems to be more reasonable .
we would agree that the solutions in @xcite are likely to be self - consistent since their @xmath150 is significantly less than unity , in particular for the solutions with @xmath153 .
then what are the reasons for the quantitative difference in the advection for @xmath153 between their solutions and ours ? in our understanding , there exist three possible reasons as follows .
first , as mentioned in section  3 of @xcite for their numerical methods , the vertical structure is derived by a given advection factor @xmath123 in advance .
the value of @xmath123 is probably obtained by solving the radial structure on the equatorial plane . moreover ,
their @xmath123 is assumed to be uniform in the @xmath6 direction . on the contrary , we obtain a varying @xmath123 by solving the vertical equations . as shown by our figure  8 , @xmath123 increases significantly with @xmath6 .
we can therefore expect that the vertically averaged @xmath123 at a cylindrical radius will also be significantly larger than that at @xmath154 , which may explain why our @xmath123 is larger than that in @xcite .
second , @xcite assumed a uniform @xmath155 and @xmath33 in the @xmath6 direction , and used the keplerian strain to calculate the viscous dissipation . in our method , however , we include varying @xmath32 and @xmath33 in the vertical direction , and the viscous dissipation is calculated based on @xmath33 instead of @xmath156 .
third , @xcite assumed @xmath157 whereas we have @xmath25 .
as stressed by @xcite , since the stationary accretion flows resemble quasi - spherical flows ( @xmath158 constant ) much more than quasi - horizontal flows ( @xmath159 constant ) , @xmath25 may be a more reasonable approximation than @xmath157 .
moreover , the above three reasons may also be responsible for the different results in the convective stability , as will be discussed in section  5.3 .
as mentioned in section  1 , some works on the black hole spin measurement showed that the standard thin disk model is likely to be inaccurate for @xmath3 ( e.g. , * ? ? ?
* ; * ? ? ?
one explanation is that the inner disk edge is still located at the innermost stable circular orbit ( isco ) , but its emission is shaded by the outer disk .
thus , the inner disk radius obtained from the spectral fitting is not true .
however , @xcite showed that the disks in black hole and neutron star x - ray binaries trace the same evolutionary pattern for @xmath3 .
in addition , for the neutron star system xte j1701 - 462 , the boundary emission area maintains nearly constant despite the varying luminosity of the disk ( * ? ? ?
* figure  17 ) , which indicates that the neutron star s surface is not shaded .
@xcite therefore argued that the inner disk of the black hole system should not be shaded either due to the similar phenomenon .
they suggested that the inner disk radius moves outward because of the increasing radiation pressure .    in our opinion , from the energy advection , it is easy to understand that the standard disk model seems to be inaccurate above @xmath160 .
as revealed by the lines in figure  4 , @xmath123 is likely to be non - negligible ( probably @xmath161 ) for @xmath162 at @xmath21 .
we would point out that , compared with the paczyski - wiita potential , the newtonian potential in the present work may magnify the viscous heating rate at small radii such as @xmath163 , thus the real @xmath123 at @xmath163 may be smaller than the values showed in figure  4 . on the other hand , for the same @xmath19 ,
since the viscous heating rate at a smaller radius such as @xmath164 will probably be larger than that at @xmath91 , so does the advection factor .
we can therefore expect that , even for the paczyski - wiita potential , the advection at the position close to the isco should be non - negligible for @xmath162 .
consequently , the standard thin disk model , based on the energy balance between the viscous heating and the radiative cooling with the advective cooling being ignored , may be inaccurate .
figure  2 shows that @xmath15 decreases significantly with decreasing @xmath22 , and figure  4 implies that most of the accreted matter exists in the range @xmath101 . in order to have a more clear view ,
we define an averaged dimensionless height as @xmath165 , where the surface density @xmath166 takes the form : @xmath167 figure  5 shows the variation of @xmath123 with @xmath168 ( solid line ) . even though the photosphere is close to the polar axis ,
the averaged height @xmath168 is geometrically slim with @xmath169 .
furthermore , the figure shows that @xmath123 increases with increasing @xmath168 or @xmath19 , which agrees with the classic picture . for quantitative comparison ,
we plot the function @xmath170 ( dashed line ) in figure  5 due to the relationship @xmath171 introduced by @xcite , which is equivalent to @xmath172 here .
it is seen that @xmath123 is not well proportional to @xmath173 . in the range @xmath174 or @xmath175 ,
however , we may regard the formula @xmath170 as a rough approximation .
in section  4 , we focus on the solutions at a typical radius @xmath21 . in this section we will study the disk solutions for various radii .
since the vertical solutions are based on the assumptions of partial derivatives in the radial direction ( presented in section  2.1 ) , it is necessary to derive vertical solutions for various radii to check whether these assumptions are self - consistent .
following the example solution in figure  2 , we study the two - dimensional solutions for @xmath87 in the range @xmath176 and @xmath89 .
figure  6 shows the radial variations of @xmath10 and @xmath46 ( solid lines ) for five polar angles , i.e. , @xmath177 , @xmath178 , @xmath179 , @xmath180 , and @xmath181 . for comparison , the radial profile of @xmath156 , which is proportional to @xmath182 ,
is shown in figure  6(a ) , and an example slope of @xmath183 is shown in figure  6(b ) .
the figure shows that @xmath10 and @xmath46 behave roughly as @xmath184 and @xmath183 , respectively , which agrees with the original assumptions of radial derivatives .
as mentioned in section 2.1 , once the radial derivatives of @xmath10 and @xmath46 are given , the other ones can be inferred from equations ( 1)-(7 ) .
thus , we can expect that the radial derivatives of @xmath32 , @xmath33 , @xmath15 , and @xmath34 in the two - dimensional solutions should also be in agreement with the assumptions .
our solutions are therefore likely to be self - consistent .    in our calculation ,
the location of the photosphere does not vary much with the radius , i.e. , @xmath185 for various radii . as discussed in section  4.1 , the real position of the photosphere is likely to be located in the range @xmath186 .
we will make some comparison with simulations for the photosphere in section  6 .
as shown in figure  6 , the derivatives of @xmath10 and @xmath46 deviate a little for @xmath187 and @xmath177 .
such a divergence may be well understood from figure  1 , which shows that the solution for @xmath87 and @xmath187 is quite close to the lower solid line , which indicates that the gas pressure may not be negligible .
in particular for the equatorial plane , the mass density has the maximal value there , thus the gas pressure may be most significant at this position .      in this section
, we will study the convective stability of the radiation pressure - supported disks .
the well - known solberg - hiland conditions in cylindrical coordinates ( @xmath4 , @xmath5 , @xmath6 ) take the forms ( e.g. , * ? ? ?
* ) : @xmath188 @xmath189 where @xmath190 is the specific angular momentum per unit mass , @xmath191 is the total pressure , @xmath192 is the specific heat at constant pressure , and @xmath193 is the entropy expressed as @xmath194 where @xmath195 is the adiabatic index .
the @xmath4 and @xmath6 components of the well - known brunt - visl frequency are written as @xmath196 @xmath197 and the epicyclic frequency takes the form : @xmath198 thus , the first solberg - hiland condition , equation  ( 21 ) , can be simplified as @xmath199 where @xmath200 is defined as an effective frequency . for accretion disks
, there usually exists @xmath201 ( as shown in figure  8) , so the second solberg - hiland condition , equation  ( 22 ) , reduces to @xmath202 in numerical calculation , we adopt @xmath203 and @xmath204 according to the radiation pressure - supported assumption .
based on the two - dimensional solutions for @xmath87 in section  5.1 , we can obtain the variations of physical quantities in cylindrical coordinates and therefore investigate the convective stability by equations  ( 24)-(25 ) .
we take the cylindrical radius @xmath205 as a typical position to study the convective stability .
figure  7 shows the @xmath6-direction variations of @xmath206 , @xmath207 , @xmath208 , @xmath209 , and @xmath210 , where the former four quantities are normalized by @xmath211 , and the last one is normalized by @xmath212 .
the positive values for both @xmath209 and @xmath210 indicate that the disk should be convectively stable . for the equatorial plane ,
the result of @xmath213 can be inferred by equation  ( 15 ) of @xcite , which revealed that the disk will always be convectively stable for @xmath204 at @xmath154 .    for further understanding the convectively stable results for @xmath214 , we plot figure  8 to show the @xmath6-direction variations of @xmath15 , @xmath46 , @xmath190 , @xmath215 , and the advection factor @xmath123 .
it is seen that @xmath15 drops faster than @xmath46 with increasing @xmath6 , so the quantity @xmath215 increases with @xmath6 . from equation
( 23 ) we immediately have @xmath216 which is known as the schwarzschild criterion for a constant angular velocity at a cylindrical surface ( @xmath217 ) , corresponding to the so - called barytropic flows where the pressure depends only on the density . as the dashed line in figure  8 shows , the angular momentum @xmath190 ( or equivalently the angular velocity @xmath110 ) does not vary significantly with @xmath6 .
thus , the convectively stable results are easy to understand from @xmath218 .    in a similar study ( vertical structure based on the local energy balance ) of the general model for optically thick disks ,
however , the disk was found to be convectively unstable ( e.g. , * ? ? ?
* ; * ? ? ?
as interpreted by the three possible reasons in section  4.2 , the significant difference in the results between their works and ours is probably related to the different approach to describing the vertical structure .
in addition , here we would mention two more details which may help to understand the difference in convective stability .
first , the profiles of @xmath15 and @xmath46 in figure  8 reveal that the absolute value of radial velocity @xmath97 increases with @xmath6 ( @xmath219 inferred from equations  ( 7 ) and ( 10 ) ) . compared with the uniform @xmath155 in @xcite , our increasing @xmath97 with @xmath6
may result in faster drop of @xmath15 in the @xmath6 direction ( steeper slope of @xmath15 ) if we simply assume the mass supply to be comparable .
second , compared with @xcite , the advection in our solutions is significantly stronger , which means that for the same @xmath19 thus comparable viscous heating rate , the vertical radiation flux @xmath220 will be less in our results .
thus , @xmath46 may decrease slower in the @xmath6 direction ( flatter slope of @xmath46 ) due to the less @xmath220 and lower @xmath15 ( inferred from equation  ( 6 ) ) . as indicated by equation
( 26 ) , the flatter slope of @xmath46 and the steeper slope of @xmath15 will both make contribution to @xmath218 , and the disk is therefore likely to be convectively stable .
we would stress that , our solutions are limited by the radiation pressure - supported case .
thus , the present results can not directly show the convective stability of disks for either the gas pressure - supported case or the case of comparable gas and radiation pressure .
actually , we have made some additional calculation for thin disks in cylindrical coordinates to check the convective stability , following the method of @xcite but without considering advection .
we found that the disk is convectively stable for the gas pressure - supported case , whereas the disk is convectively unstable for the case of radiation pressure being significant .
thus , we would agree with @xcite on the convectively unstable disks for moderate accretion rates such as @xmath221 , corresponding to significant radiation pressure and non - negligible gas pressure . moreover , for @xmath222 , the disk will be geometrically thin , thus there is no difference between the assumptions @xmath25 and @xmath157 , and @xmath123 is probably negligible . as a consequence , the solutions of @xcite ought to be accurate .
furthermore , as revealed by the profiles of @xmath46 and @xmath15 in figure  8 , @xmath117 is less than unity in the @xmath6 direction . following the argument in section  4.1
, the polytropic relation seems unsuitable either in the @xmath6 direction .
in addition , as shown in figures  7 and 8 , our example solution for @xmath87 at @xmath205 is terminated at @xmath223 .
the reason is that the solutions in spherical coordinates is limited by @xmath224 ( as shown by the horizontal dot - dashed line in figure  1 ) , which corresponds to @xmath223 at @xmath205 in cylindrical coordinates .      for the region close to the equatorial plane
, we can make some analysis of the convective stability by the taylor expansion method . obviously , we have @xmath226 at @xmath154 from symmetric conditions .
thus , for @xmath225 , the value of @xmath227 can be estimated by the second - order derivative at @xmath228 : @xmath229    based on equations  ( 7)-(12 ) , we can eliminate @xmath32 , @xmath46 , and @xmath35 and therefore obtain a set of three equations for the three quantities @xmath33 , @xmath10 , and @xmath15 . by using the taylor expansion method , together with the boundary condition of equation  ( 14 )
, we derive the following three relationships for the second - order derivatives of @xmath33 , @xmath10 , and @xmath15 ( the ram - pressure term in equation  ( 8) is ignored ) : @xmath230 @xmath231 @xmath232 where @xmath233 is defined as @xmath234 , which is a small value in the analysis .
in our solutions we have @xmath235 at @xmath154 , e.g. , @xmath236 for @xmath87 at @xmath237 .
for the simple case with @xmath238 , the above three relationships provide @xmath239 the second - order derivative of entropy at @xmath240 can therefore be derived by the following coordinate transformation : @xmath241 thus , equations ( 27)-(28 ) indicate @xmath218 for the region close to the equatorial plane .
the disk in this region is therefore likely to be convectively stable .
in this paper , we have studied the vertical structure , energy advection , and convective stability of radiation pressure - supported disks in spherical coordinates . in the @xmath22 direction , we replaced the pressure - density polytropic relation by the local energy balance per unit volume between the viscous heating and the advective cooling plus the radiative cooling , and obtained the distribution of physical quantities such as @xmath15 , @xmath32 , @xmath33 , @xmath10 , @xmath46 , and @xmath35 .
the photosphere was found close to the polar axis and therefore the disk seems to be extremely thick .
however , most of the accreted matter exists in a moderate range around the equatorial plane such as @xmath242 .
we showed that the polytropic relation is unsuitable for describing the vertical structure of radiation pressure - supported disks .
more importantly , we found that the energy advection is significant even for slightly sub - eddington accretion disks , which is quite different from the previous result that the advection is of importance only for super - eddington accretion disks .
we argued that , the non - negligible advection may help to understand why the standard thin disk model is likely to be inaccurate for @xmath3 .
in addition , we studied the two - dimensional solutions to check our basic assumptions of radial derivatives , which indicates our solutions to be self - consistent .
furthermore , we investigated the convective stability of the disks in cylindrical coordinates by the two - dimensional solutions derived in spherical coordinates .
the disk solutions satisfy the solberg - hiland conditions , which reveals that the disk ought to be convectively stable .    to our knowledge , there are mainly two series of simulations on optically thick accretion flows .
one is three - dimensional , radiation magnetohydrodynamic ( rmhd ) simulations by a shearing box , which focused on thin disks ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
the other is two - dimensional , either radiation hydrodynamic ( rhd ) or rmhd simulations for global flows , which included thin disks and super - eddington accretion flows ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
we would like to compare our numerical results with the latter owing to the global , radiation pressure - supported , and geometrically not thin case .
@xcite studied rhd simulations in spherical coordinates with the @xmath41 stress assumption for only the @xmath112 component , for extremely high accretion rates @xmath243 , @xmath244 , and @xmath245 . moreover ,
@xcite studied rmhd simulations in cylindrical coordinates for @xmath246 . in our results ,
the position of the photosphere is quite close to the polar axis with @xmath185 .
as discussed in section  4.1 , the real position may exist in the range @xmath186 .
in the simulations , @xcite showed that the photosphere is around @xmath247 corresponding to @xmath248 .
@xcite did not mention the position of photosphere , but their figure  4 for the density distribution reveals that the photosphere should exist at a certain @xmath67 significantly less than @xmath180 .
thus , we would express that , for accretion rates around and beyond @xmath17 , the photosphere may exist far from the equatorial plane with @xmath249 .
in addition , we discuss the possible link between our disk model and ultraluminous x - ray sources .
as shown by the upper solid line in figure  1 , there exists a maximal accretion rate @xmath250 varying with the radius .
we argue that the possible upper limit of the accretion rate may help to understand why most ultraluminous x - ray sources ( ulxs ) are not in thermal dominant state .
as mentioned in the review paper of @xcite , there may exist three classes of black holes in ulxs : normal stellar mass black holes ( @xmath251 ) , massive stellar black holes ( @xmath252 ) , and intermediate mass black holes ( @xmath253 ) .
a massive stellar black hole with a moderate super - eddington accretion rate seems to account for most sources up to luminosities @xmath254 a few @xmath255 erg s@xmath73 .
the slim disk model , which is the classic model for super - eddington accretion disks , predicts a dominant thermal radiation from the disk .
however , observations have shown that ulxs are not in the thermal dominant state except only a few sources such as m82 x-1 @xcite and hlx-1 @xcite . in our understanding
, the radiation of ulxs may be interpreted by an optically thick disk with @xmath256 plus strong outflows .
the disk will provide a thermal radiation , which is normally not dominant because of the moderate @xmath250 . on the other hand
, the outflows may make contribution to the non - thermal radiation by the bulk motion comptonization @xcite or through the jet of the radiation - pressure driven and magnetically collimated outflow @xcite .
the author is particularly grateful to ramesh narayan for constructive suggestions and beneficial discussions .
the author also thanks ivan hubeny , da - bin lin , aleksander sadowski , lijun gou , yucong zhu , and shan - shan weng for helpful discussions , and the referee for providing useful comments to improve the paper .
this work was supported by the national natural science foundation of china under grants 11073015 and 10833002 , the national basic research program ( 973 program ) of china under grant 2009cb824800 , and the scholarship from china scholarship council under grant 2009835057 .
abramowicz , m. a. , chen , x. , kato , s. , lasota , j .-
p . , & regev , o. 1995 , , 438 , l37 abramowicz , m. a. , czerny , b. , lasota , j .-
p . , & szuszkiewicz , e. 1988 , , 332 , 646 abramowicz , m. a. , lanza , a. , & percival , m. j. 1997 , , 479 , 179 begelman , m. c. , & meier , d. l. 1982 , , 253 , 873 davis , s. w. , narayan , r. , zhu , y. , et al .
2011 , , 734 , 111 feng , h. , & kaaret , p. 2010
, , 712 , l169 feng , h. , & soria , r. 2011 , new astron . rev .
, 55 , 166 gu , w .- m . , & lu , j .- f .
2007 , , 660 , 541 gu , w .- m . , xue , l. , liu , t. , & lu , j .- f .
2009 , , 61 , 1313 hirose , s. , blaes , o. , & krolik , j. h. 2009a , , 704 , 781 hirose , s. , krolik , j. h. , & blaes , o. 2009b , , 691 , 16 hirose , s. , krolik , j. h. , & stone , j. m. 2006 , , 640 , 901 hshi , r. 1977 , prog .
phys . , 58 , 1191 kato , s. , fukue , j. , & mineshige , s. 2008 , black - hole accretion disks : towards a new paradigm ( kyoto : kyoto univ . press ) krolik , j. h. , hirose , s. , & blaes , o. 2007 , , 664 , 1045 lin , d. , remillard , r. a. , & homan , j. 2009 , , 696 , 1257 mcclintock , j. f. , shafee , r. , narayan , r. , et al .
2006 , , 652 , 518 narayan , r. , & yi , i. 1994 , , 428 , l13 narayan , r. , & yi , i. 1995 , , 444 , 231 ohsuga , k. , & mineshige , s. 2011 , , 736 , 2 ohsuga , k. , mineshige , s. , mori , m. , & kato , y. 2009 , , 61 , l7 ohsuga , k. , mori , m. , nakamoto , t. , & mineshige , s. 2005 , , 628 , 368 paczyski , b. , & wiita , p. j. 1980 , , 88 , 23 sadowski , a. 2011 , phd thesis , arxiv:1108.0396 sadowski , a. , abramowicz , m. , bursa , m. , et al .
2011 , , 527 , a17 sadowski , a. , abramowicz , m. , bursa , m. , et al .
2009 , , 502 , 7 shakura , n. i. , & sunyaev , r. a. 1973 , , 24 , 337 straub , o. , bursa , m. , sadowski , a. , et al . 2011 , , 533 , a67 tassoul , j .- l . 2000 , stellar rotation ( cambridge : cambridge univ .
press ) titarchuk , l. , & zannias , t. 1998 , , 493 , 863 watarai , k .- y . , fukue , j. , takeuchi , m. , & mineshige , s. 2000 , , 52 , 133 weng , s .- s .
, & zhang , s .- n .
2011 , , 739 , 42


the understanding of confinement and , more generally , of the principles that govern the low energy dynamics in qcd is , perhaps , the most interesting unanswered question in the theory of strong interactions .
the strongly coupled dynamics of the four dimensional non abelian gauge theories has so far proven to be forbiddingly complicated .
a fruitful approach to cope with such situation can be to turn to simpler , more tractable models , that nevertheless capture at least some of the important features of the real theory .    from this point of view , three dimensional non abelian gauge theories are of particular interest since they are , indeed , confining , and have a strong coupling dynamics .
on the other hand they can be continuously deformed into the weak coupling region without encountering any phase transition .
the prototypical case is the @xmath1 gauge theory with an adjoint higgs field - the georgi - glashow model .
these weakly coupled theories are confining , just like their strongly coupled counterparts . here , however , confinement can be studied using semiclassical methods , as first done by polyakov for the @xmath1 case  @xcite , and later on generalized to @xmath2 theories  @xcite .
a related feature of these models is that their spectra has a large gap between the light ( light ` photons ' ) and heavy ( heavy massive gauge bosons ) states .
this makes the identification of the relevant low energy degrees of freedom relatively easy , and also allows one to construct the low energy effective lagrangian explicitly .
it turns out that this effective lagrangian exhibits the confinement phenomenon in a very simple and straightforward way _ on the classical level_. in a nutshell , the natural degrees of freedom turn out to be magnetic vortex fields first introduced by thooft  @xcite , and the confinement mechanism is realized as confinement of topological solitons for those vortex fields .
the structure of the effective lagrangian is robustly determined by the spontaneously broken magnetic @xmath0 symmetry  @xcite .
this approach has been motivated by earlier work  @xcite , it was advanced and formulated in  @xcite , and further clarified in  @xcite .
it has also been argued in  @xcite , that the same type of effective lagrangian remains valid in the strongly coupled limit of pure gluodynamics .
in this situation , the components of the vortex field are naturally identified with the scalar and pseudoscalar glueballs , which according to lattice results are indeed the two lowest states in the spectrum of @xmath3 dimensional gluodynamics  @xcite .
the confinement mechanism for the strongly coupled regime remains , in this picture , essentially the same as for the weakly coupled phase .
all previous work dealt exclusively with theories having no matter fields in fundamental representation .
the purpose of this note is to address the question of what effect dynamical fundamental ` quarks ' have on this low energy structure . that that effect exists and must be nontrivial can be seen from one basic observation : in the presence of fundamental charges , the vortex field ( which is the effective low energy degree of freedom ) is no longer a local field .
the effective lagrangian can not , therefore , be a simple local scalar theory as in  @xcite .
we will show below that , in the presence of fundamental matter , the magnetic @xmath0 in the effective lagrangian is _ gauged _ , the value of the gauge coupling being inversely proportional to the mass of the ( lightest ) dynamical quark .
we will construct such an effective lagrangian , which moreover in the limit of the infinitely heavy quarks reduces to the effective theory of  @xcite .
the presence of fundamental quarks implies the existence of a conserved baryon number charge .
we will see how this charge is represented in the low energy theory , and will study the structure of finite energy baryon states .
quite surprisingly , we find that these states naturally have a bag - like structure .
the value of the vortex field in the region of space which contains the quarks , is much lower than its value in the vacuum .
the potential energy associated with this difference in the expectation values can also be naturally identified with the bag constant .
the plan of this paper is the following . in section  [ glue ] , we give a lightning review of the results of  @xcite .
we discuss the magnetic @xmath0 symmetry , its order parameter ( the vortex operator ) , and the structure of the effective lagrangian in pure gluodynamics . in section  [ fund ]
we explain the basic physics of the changes in this structure in the presence of the fundamental charges . in section  [ low ]
we show how to define the relevant low energy degrees of freedom , and how to construct the effective lagrangian in this case . in section  [ baryon ]
, we study the structure of the baryon states .
finally we conclude in section  [ discus ] with a short discussion .
consider the pure su(n ) gauge theory in @xmath3 dimensions @xmath4 as was argued by thooft in  @xcite , this theory possesses a global @xmath0 symmetry which is spontaneously broken in the vacuum .
the explicit realization of this symmetry was found in @xcite .
the generator of the @xmath0 group is ( up to a multiplicative constant factor ) the fundamental wilson loop along the spatial boundary of the system @xmath5 in  @xcite , the conservation of @xmath6 was proved taking the pure gluodynamics limit of the georgi - glashow model .
this , however , can also be shown directly by calculating the commutator of @xmath6 with the yang - mills hamiltonian @xmath7 : @xmath8=\lim_{c\to\infty}\oint_c dx_i { \rm tr } \left [ p e_i(x)\exp\{i\oint_{c(x , x ) } dy_ia^i(y)\ } \right ] \to_{c\to\infty}0 \;. \label{commutator}\ ] ] here , the integral in the exponential on the right hand side starts and ends at the point of insertion of the electric field .
the vanishing of the commutator follows from the fact that it only involves electric fields at spatial infinity , and in a theory with a finite mass gap those fields should vanish at infinity .
to dispel any doubt about the vanishing of this commutator we note that the situation is completely analogous to the commutator of any ` conserved charge ' which is defined as an integral of a local charge density @xmath9 the commutator of such a charge with a hamiltonian contains also a surface term , since the charge density @xmath10 never commutes with the hamiltonian density but rather gives a total derivative in the commutator . for a conserved charge , due to the continuity equation , this surface term is equal to the circulation of the spatial component of the current @xmath11=i\oint_c dx^ij_i\;.\ ] ] the vanishing of this term again is the consequence of the vanishing of the physical fields at infinity in a theory with a mass gap .
when the charge is not conserved , the commutator contains also a bulk term in addition to the surface contribution .
it is the absence of the bulk terms that is the distinctive property of a conserved charge .
the same conclusion is reached if , rather than considering the generator of the algebra , one considers the commutator of the group element for either continuous or discreet symmetry groups .
we see , therefore , that the commutator in ( [ commutator ] ) indeed tells us that @xmath6 is a conserved operator .
next , a simple argument establishes that the magnetic @xmath0 symmetry is spontaneously broken in the vacuum if the wilson loop @xmath12 has an area law behaviour  @xcite .
the argument goes like this .
the vev of @xmath12 is the overlap of the vacuum state @xmath13 with the state @xmath14 which is obtained from the vacuum by acting on it with @xmath12 .
when acting on the vacuum state , @xmath12 performs the @xmath0 transformation at all points within the area @xmath15 bounded by the loop .
if the vacuum wavefunction depends on the configuration of the @xmath0 non invariant degrees of freedom ( the state in question is not @xmath0 invariant ) the action of @xmath12 affects the state everywhere inside the loop . in the local theory with finite correlation length
the overlap between the two states approximately factorizes into the product of the overlaps taken over the region of space of linear dimension of order of the correlation length @xmath16 @xmath17 where the label @xmath18 is the coordinate of the point in the center of a given small region of space . for @xmath18 outside the area @xmath15
the two states @xmath19 and @xmath20 are identical and therefore the overlap is unity .
however for @xmath18 inside @xmath15 the states are different and the overlap is therefore some number @xmath21 smaller than unity .
the number of such regions inside the area is obviously of order @xmath22 and thus @xmath23 the vev of @xmath12 then falls off as an area .
conversely , if the vacuum is @xmath0 invariant , the wavefunction does not depend on the configuration of the non invariant degrees of freedom .
the action of @xmath12 then alters the state only along the perimeter and @xmath12 has the perimeter law behaviour in the unbroken phase .
thus the area law behaviour of the fundamental wilson loop is tantamount to the breaking of the magnetic @xmath0 symmetry in the vacuum of pure gluodynamics .
the only requirement that this argument presupposes is the existence of * local * degrees of freedom which are non invariant under @xmath0 , or in other words the existence of a * local * order parameter .
such an operator can indeed be constructed explicitly .
this is the magnetic vortex creation operator @xmath24  @xcite,@xcite,@xcite .
the defining property of @xmath24 is that it satisfies the following commutation relation with the spatial fundamental wilson loop @xmath25 @xmath26 where @xmath27 is the linking number of the curve @xmath28 and the point @xmath18 , or in other words , @xmath29 if @xmath30 , and @xmath31 if @xmath32 , where @xmath15 is the region of the plane bounded by @xmath28 . an explicit representation for @xmath33
is given by writing it as the operator that performs a ` singular gauge transformation ' @xmath34 here , the hypercharge generator @xmath35 is defined as @xmath36 the electric field is written in the matrix notation @xmath37 , with @xmath38 denoting the @xmath2 generator matrices in the fundamental representation , and @xmath39 is the hypercharge density due to gluons @xmath40 $ ] .
the function @xmath41 is the planar angle of the vector @xmath18 . using gauss law @xmath42 this operator ( on the states that satisfy gauss law )
may be written in an alternative representation @xmath43 the integration here is along the branch cut in the definition of the planar angle @xmath44 , which is an infinite line that starts at the point @xmath18 and goes to infinity .
there are two crucial properties of the definition eq.([v2 ] ) that insure that the operator @xmath24 is a local physical field :    1
.   gauge invariance .
the definition is not explicitly gauge invariant .
nevertheless , one can show  @xcite,@xcite that @xmath24 transforms a physical state into another physical state .
that is to say that the matrix elements of @xmath24 between a physical and an unphysical state vanish . by physical state we mean a state which satisfies the non abelian gauss law . by virtue of this property the definition eq.([v2
] ) can be used as long as the matrix elements of @xmath24 ( or any power of @xmath24 with insertions of an arbitrary number of gauge invariant operators ) are calculated between physical states .
2 .   locality .
again , the definition is not explicitly local , since it contains the integral of the electric field along an infinite line .
however in a theory that does not contain fundamental charges , @xmath24 in fact does not depend on the curve @xmath28 , but only on its endpoint @xmath18 .
physically , this is due to the fact that the bohm - aharonov phase between @xmath24 and any charged state in the theory vanishes .
more formally , consider changing the position of the curve @xmath28 to @xmath45 .
this adds to the phase in the definition eq.([v2 ] ) a contribution @xmath46 , where @xmath15 is the area bounded by @xmath47 . in our normalization ,
the hypercharge of gluons is @xmath48 or @xmath49 and no particles with smaller value of the hypercharge are present in the theory .
therefore the hypercharge within any closed area is an integer multiple of the gauge coupling @xmath50 , and the extra phase factor is always unity .
one can also show directly  @xcite , using canonical commutation relations , that @xmath33 commutes with all the local gauge invariant operators @xmath51 , unless @xmath52 , which formally establishes the locality of @xmath24 .
the universal properties of the low energy dynamics of the pure yang mills theory are determined by the spontaneous breaking of the magnetic @xmath0 symmetry .
the effective low energy lagrangian is constructed in terms of the vortex field @xmath24 in principle we have to include in this lagrangian higher order potential terms which stabilize the potential energy at large values of @xmath24 .
although we will not indicate them explicitly , it is assumed in the following that they are indeed present whenever necessary . ]
@xmath53 this lagrangian can be derived in a controllable way in the weakly coupled phase@xcite,@xcite .
the coupling constants in eq.([ldualgg ] ) are determined in the weakly coupled region from perturbation theory and dilute monopole gas approximation .
it was argued in @xcite that the same low energy lagrangian is also valid in the strongly coupled regime ( pure gluodynamics ) .
several features of this lagrangian are worth noting :    * first , @xmath54 , and therefore the @xmath0 symmetry is spontaneously broken . *
second , the field components of @xmath24 can be reinterpreted in terms of the physical glueball fields . using a polar decomposition for @xmath24 @xmath55
we find that the fluctuation field @xmath10 is the charge conjugation even scalar , while the phase @xmath56 is the conjugation odd pseudoscalar . indeed , according to lattice results  @xcite , the two lightest glueballs in the spectrum of pure gluodynamics are the scalar @xmath57 and pseudoscalar @xmath58 states .
the lagrangian ( [ ldualgg ] ) , is thus the natural effective lagrangian that one would write to describe the dynamics of the lightest particles . in this respect
the identification of the glueball fields with the modulus and the phase of @xmath24 constrains the possible terms in the effective lagrangian by imposing the @xmath0 symmetry on possible interaction terms .
* finally , this effective lagrangian exhibits in a simple way confining properties of the theory .
the colour charged ` gluon ' states ( or massive vector bosons in the weakly coupled regime ) are represented in the low energy theory by topological solitons carrying unit winding of the field @xmath24 . if not for the @xmath59 symmetry breaking term in the potential in ( [ ldualgg ] )
, the lowest lying state in the topologically nontrivial sector would have a rotationally invariant hedgehog - like structure .
its energy would be logarithmically divergent in the infrared , which corresponds to the logarithmic coulomb potential in 2 + 1 dimensions .
however since the symmetry of the potential is @xmath0 rather than @xmath59 , the energy of the hedgehog configuration diverges quadratically with the volume .
this is so since the number of vacuum states is discreet and in the hedgehog configuration the phase of the field @xmath24 is far from its vacuum value everywhere in space .
the configuration that actually minimizes the energy in the one soliton sector breaks rotational invariance by confining the winding of the phase of @xmath24 to some quasi one - dimensional strip , like in figure 1 .
it is clear that , dynamically , the width of this strip is of the order of the inverse glueball mass and the strip itself is nothing but the confining electric flux string .
the energy of such an isolated soliton now diverges _ linearly _ in the infrared , due to the finite energy density per unit length of the string .
[ fig1 ]   in the state of unit charge in the presence of the symmetry breaking terms in the effective lagrangian , for the @xmath60 model.,title="fig : " ]    the soliton we are discussing is a dynamical particle in the effective theory - it thus represents a particle with an adjoint color charge , and the string is the adjoint string .
this string is of course not absolutely stable since it can break into a soliton - antisoliton pair when it is too long .
the fundamental string also appears naturally in the low energy description .
it is the stable domain wall which separates the regions of space with different expectation values of the vortex field @xmath24 @xcite , see figure 2 .
in this section we wish to discuss the modifications introduced by the presence of fundamental quarks in the structure just described . for simplicity
we will consider scalar quarks .
the lagrangian of the theory we are interested in is @xmath61 where the scalar field @xmath62 transforms according to the fundamental representation of the @xmath2 color group .
first thing to note is that the fundamental wilson loop still commutes with the hamiltonian .
this is obvious , since the extra terms in the hamiltonian in the presence of the quarks do not involve electric field operator , but only vector potential .
the wilson loop commutes with the vector potential , and thus with the additional terms in the hamiltonian .
we thus conclude that the theory still has the magnetic @xmath0 symmetry .
this may seem somewhat surprising at first sight .
we have seen in the previous section that spontaneous breaking of the magnetic @xmath0 implies the area law for the wilson loop , and conversely the perimeter law of @xmath25 implies unbroken @xmath0 . in the theory with fundamental charges
the wilson loop is known to have perimeter law due to breaking of the confining string at any finite value of the fundamental mass @xmath63 .
we might then conclude that the magnetic @xmath0 is restored at any , arbitrarily large value of @xmath63 but is broken at @xmath64 .
the common lore is that the @xmath0 breaking phase transition for @xmath65 is first order , and this then implies a discontinuous behaviour of the theory in the infinite mass limit .
this of course is completely counterintuitive and in fact plain wrong .
the caveat in this line of reasoning is the following .
the relation between the behaviour of the wilson loop and the mode of the realization of the magnetic symmetry hinges crucially on the existence of a local order parameter of the magnetic @xmath0 . in the absence of such an order parameter
it is not true that the wilson loop locally changes the quantum state inside the loop and this invalidates the whole argument .
in particular in the absence of a local order parameter , the @xmath0 symmetry can be spontaneously broken but the wilson loop can have a perimeter law .    in fact it is easy to see that for any finite @xmath63 the magnetic @xmath0 does not have a local order parameter .
the only candidate for such an order parameter is the vortex operator @xmath33 defined in eq.([v2 ] ) , since it has to be local also relative to the purely gluonic operators by any explicitly local gauge invariant operator .
such modifications are however irrelevant as far as the locality properties of @xmath24 are concerned , and we will not consider them in the following . ] .
however in the presence of fundamental quarks the operator @xmath33 is not local anymore . to see this
consider the dependence of @xmath66 on the curve @xmath28 which enters its definition .
as before the operators @xmath67 and @xmath68 are related by @xmath69 where @xmath15 is the area bounded by @xmath47 .
as before , due to the gauss law the integral in the exponential is equal to the total hypercharge in the area @xmath15 .
however the hypercharge of fundamental quarks has eigenvalues @xmath70 .
the extra phase factor is therefore not unity anymore but can rather take values @xmath71 depending on the state and the choice of the contour @xmath28 .
the status of the magnetic @xmath0 is thus quite different in a theory with fundamental quarks - it does not have a local order parameter .
nevertheless it is clear that at least as long as the mass @xmath63 is large @xmath72 , the relevant degrees of freedom for the effective infrared dynamics should still be the vortices @xmath24 and the main factor which determines their dynamics should still be the magnetic @xmath0 . at large @xmath63
the dependence of @xmath24 on the curve @xmath28 is weak , since the probability of the vacuum fluctuations which involve fundamental charges is small .
the probability of appearance of a virtual @xmath73 pair separated by a distance @xmath16 is suppressed by the exponential factor @xmath74 .
the typical distance scale for the `` glueball '' physics is @xmath75 .
thus at these distances such fluctuations are unimportant and should not affect much the dynamics .
things are different however if one is interested also in the baryonic sector of the theory .
the baryons are necessarily heavy and in order to be able to discuss their structure we must understand the main dynamical effects also at shorter distances .
therefore our aim now is to understand what is the main effect of the non locality discussed above on the dynamics of magnetic vortices .
the situation we have just described - a symmetry without a local order parameter - is not exceptional in quantum field theory .
this is precisely the property of the global part of any abelian gauge group .
consider for example quantum electrodynamics .
the global electric charge is of course a physical gauge invariant charge with the corresponding gauge invariant local charge density .
@xmath76 nevertheless there is no local operator that carries this charge .
this is a direct consequence of the gauss law @xmath77 any physical , gauge invariant operator that carries @xmath78 must also carry the long range electric field , which can not fall off faster than a power of the distance .
the gauge invariant qed lagrangian is written in terms of `` local '' charged fields @xmath79 .
but appearances are deceptive : these fields are not gauge invariant , and therefore not physical .
a gauge invariant charged field can be constructed from @xmath79 by multiplying it by a phase factor @xmath80 with the @xmath81-number field @xmath82 satisfying @xmath83 for any @xmath82 satisfying this condition the operator @xmath84 is gauge invariant , and therefore physical .
it is however necessarily nonlocal .
different choices of @xmath82 define different gauge invariant operators and correspond to different gauge fixings .
thus for @xmath85 the field @xmath84 is the field @xmath79 in the coulomb gauge , while @xmath86 corresponds to the axial gauge @xmath87 , and so on .
different definitions of @xmath84 differ from each other by a phase factor , which is precisely the gauge ambiguity of the original field @xmath79 .
the @xmath59 gauge group is the most natural abelian gauge symmetry to consider in continuum field theory .
one can however also consider discrete groups like @xmath0@xcite . in this case again no local operator that carries the global @xmath0 charge exists .
the various gauge invariant charged operators are nonlocal and differ from each other by a local @xmath0 - valued phase .
these different operators again correspond to different gauge fixings of the local @xmath0 group .
this is precisely the structure that emerged in our discussion in the earlier part of this section .
we have    \1 . global magnetic @xmath0 symmetry generated by the fundamental wilson loop .
the set of nonlocal vortex operators @xmath66 , which all carry the @xmath0 charge and differ from each other by a @xmath0 valued phase factors .
it is very suggestive therefore to think about @xmath66 as of different gauge fixed versions of a field charged under local @xmath0 .
this leads us to expect that the low energy theory we are after should be a @xmath0 gauge theory of the magnetic vortex field @xmath24 .
in fact coming back to the discussion in the beginning of this section , we see that from this vantage point it is obvious why the wilson loop has a perimeter law , even if the global @xmath0 is broken spontaneously .
the action of a wilson loop of a finite size inside the contour is a gauge transformation .
thus in physical terms locally inside the contour the new state is _ the same _ as the old one and so the overlap between the two locally is unity .
the only nontrivial contributions to the overlap come from the region close to the contour , thus giving the perimeter law .    in the next section
we shall construct the effective gauge theory and discuss its relation to the original qcd lagrangian .
although it is possible to give a more ` analytic ' derivation of the effective lagrangian , we prefer to write it down directly guided by the previous discussion .
we will then explain the physical meaning of the various fields that appear in it .
the easiest way to construct a @xmath0 gauge theory in the continuum is to consider a @xmath59 gauge theory with the higgs field of charge @xmath88 which has a large expectation value  @xcite .
consider therefore the following lagrangian : @xmath89 here @xmath90 .
we take the parameters such that , @xmath91 and @xmath92 .
the ` higgs ' field @xmath93 has a large expectation value and breaks the @xmath59 gauge symmetry down to its @xmath0 subgroup @xmath94 .
it should be noted that this gauge subgroup does not correspond to ` singular ' gauge transformations , since partial derivatives commute on the gauge parameter , despite the fact that there are @xmath95-like singularities .    below the scale determined by the expectation value @xmath96 ,
the field @xmath93 is practically frozen and its fluctuations are unimportant
. in this regime the model indeed describes the locally @xmath0 invariant theory . the global part of the gauge group is our @xmath0 magnetic symmetry generated by the wilson loop .
the larger @xmath59 gauge structure at this point is just an auxiliary trick which enables us to write down a discrete gauge theory in continuous notations .
we will however see later that it does in fact has a real physical meaning of its own and arises naturally in the effective theory .
let us first discuss how this lagrangian reduces to the effective lagrangian of the pure yang - mills theory eq.([ldualgg ] ) in the limit of the large quark mass . in this limit
not only the field @xmath93 must decouple , but also the gauge interactions of the field @xmath24 must vanish . there is another consistency requirement . in the limit of zero gauge coupling
@xmath97 the gauge @xmath59 symmetry becomes global @xmath59 and is broken due to non vanishing expectation value of @xmath93 .
the spectrum therefore contains a massless goldstone boson .
this goldstone boson is of course the longitudinal component of the gauge field @xmath98 .
there is however no such massless particle in the pure gluodynamics nor in the effective lagrangian eq.([ldualgg ] ) .
this means that the couplings in eq.([efflfund ] ) should depend on the quark mass in such a way that the vector particles remains heavy for any finite @xmath63 and its mass goes to zero very sharply only in the limit when it is completely decoupled . in terms of the goldstone boson
couplings it means that @xmath99 must be larger than any scale relevant to the dynamics of the vortex field @xmath24 .
all these conditions can be met by choosing for example @xmath100 with this choice the mass of the vector boson @xmath101 stays finite as @xmath102 and can be arbitrarily large .
the goldstone boson in the decoupling limit has an infinite @xmath99 and is completely decoupled just as the `` invisible axion '' .
we now want to understand how the basic fields present in the effective lagrangian arise in the fundamental theory eq.([lqcd ] ) .
we start with considering the symmetries .
the lagrangian eq.([efflfund ] ) has two global @xmath59 symmetries .
one is the global part of the local @xmath59 with the conserved current @xmath103 .
note that due to the direct coupling between the @xmath24 and @xmath93 fields the @xmath59 rotations of @xmath24 and @xmath93 separately are not symmetries .
the other current , conserved by virtue of the homogeneous maxwell equation is @xmath104 .
the two charges have quite different nature .
the second one , the dual magnetic flux @xmath105 has a local order parameter
. it can be constructed in a way similar to the vortex field in qed @xcite,@xcite . the global @xmath59 gauge charge ,
on the other hand does not have a local order parameter , as discussed earlier .
both these conserved currents should also exist in the fundamental theory .
it is fairly straightforward to identify them .
the qcd lagrangian has one obvious global @xmath59 charge - the baryon number .
this charge has local order parameters - gauge invariant baryon fields of qcd , and is therefore identified with the dual magnetic flux @xmath106    the second charge can be expressed in terms of the spatial current components of the first one @xmath107 \label{qt}\ ] ] it is thus the vorticity associated with the baryon number current .
the vortex operator @xmath24 which appears in eq.([efflfund ] ) is not gauge invariant and is only physical after complete gauge fixing of the @xmath59 gauge group .
after such a gauge fixing ( which amounts to multiplying @xmath24 by an operator valued phase ) the gauss law requires that on the physical states the physical operator @xmath24 carry the charge @xmath108 . due to the identification eqs.([f],[qt ] ) this leads to a surprising conclusion that any physical operator @xmath24 in the effective theory creates a vortex of the original baryon number current .
this somewhat unexpected conclusion is in fact quite natural for an eigenoperator of the magnetic @xmath0 . consider the vortex operator @xmath24 defined by eq.([v2 ] ) .
in qcd just like in pure gluodynamics , it has an alternative representation of an operator of the singular gauge transformation of the form eq.([v1 ] ) @xmath109 the only difference is that now @xmath110 is the hypercharge operator due to both , gluons and fundamental quarks @xmath111+y_{\alpha\beta}(\phi^*_\alpha\pi_\beta- \phi_\alpha\pi^*_\beta)\big]\ ] ] consider the action of this operator on a quark field @xmath112 .
the transformed quark field @xmath113 is a gauge transform of @xmath112 everywhere except along the branch cut of the function @xmath114 . across this
cut the phase of @xmath115 is discontinuous - it jumps by @xmath116 for all color components of @xmath115 . due to this discontinuity the baryon current - the global @xmath59 current of @xmath112 - does not vanish at points along the cut . calculating explicitly the action of @xmath24 on the baryon number current we find @xmath117 where @xmath118 is a unit vector normal to the branch cut @xmath28 at the point @xmath18 .
it is natural to define the local vorticity associated with the baryon number as the vorticity defined in this way is precisely the circulation of the phase of @xmath112 . ]
@xmath119 \label{vort}\ ] ] the vortex operator @xmath24 therefore creates a vortex of baryon number current with fractional vorticity @xmath116 @xmath120    the operator @xmath93 due to the gauss law also carries baryon vorticity . since its gauge coupling is @xmath88 times the coupling of @xmath24 , it creates one unit of vorticity .    in fact
this simple exercise also helps us to identify the value of the gauge coupling constant @xmath121 in the effective theory . comparing eq.([vort ] ) with eq.([qt ] ) we find @xmath122 the same relation is obtained by comparing the current algebra in the fundamental and the effective theories .
the commutator of the baryon charge density with the baryon current density in the fundamental theory is @xmath123=i\phi^*\phi\partial_i\delta^2(x - y)\ ] ] in the effective theory using the canonical commutators that follow from eq.([efflfund ] ) and the identification eq.([f ] ) we find @xmath123=i{e^2\over 4\pi^2}\partial_i\delta^2(x - y)\ ] ] again we deduce eq.([e2 ] ) .
the operator on the right hand side of this equation in the effective theory is indeed a constant .
recall , that our effective theory should be valid at long distances . in this regime in the leading order in the derivative expansion the operator @xmath124 should be approximated by its expectation value .
taking into account fluctuations of @xmath124 is tantamount to including higher derivative terms in the effective lagrangian eq.([efflfund ] ) .
at this , higher order in derivative expansion the gauge coupling constant would become a dynamical field . while this is perfectly legitimate ,
it is certainly beyond our present framework .
when taking the expectation value of @xmath124 we should remember that it has to be calculated with the cutoff @xmath125 appropriate for the effective theory .
this cutoff must be above the characteristic scale of the pure gluodynamics ( determined by the string tension ) but below the heavy quark mass . with this in mind
we get @xmath126 so that finally @xmath127 which is consistent with the expected scaling in the large mass limit eq.([largem ] ) .
having understood the origin of the fields and the symmetries in the effective theory , we would like to see how it encodes the qualitative features of the low energy qcd physics . since we are considering a heavy quark theory , below the fundamental mass scale the spectrum should be the same as in pure gluodynamics .
we have already seen that , in the infinite mass limit , the effective theory reduces to that of ( [ ldualgg ] ) .
indeed , even at finite but large @xmath63 , this is the case at low energies .
formally , this can be seen as follows .
since the vev of the field @xmath93 is large , we can impose the unitary gauge condition on it . in this unitary gauge the phase of @xmath93 disappears .
the modulus of @xmath93 is very heavy , and so is the vector field @xmath98 .
thus at low energies we recover the effective theory of pure glue sector .
there is however a set of configurations , on which the unitary gauge can not be imposed .
those are configurations in which @xmath93 vanishes at some points in space .
indeed it is these configurations that are important for the baryonic sector .
recall that @xmath93 is a vortex of baryon number current .
thus one expects that the baryon charge is associated with the vortex configuration of the field @xmath93 . in the core of the vortex
the field @xmath93 of course has to vanish , and so the unitary gauge is not admissible .
thus in the baryon sector we can not think of @xmath93 as frozen at its expectation value and instead have to treat it as a dynamical field .    that the baryon does indeed carry vorticity of the field @xmath93 ,
can be seen by the following simple argument .
the baryon number is represented in the effective theory by the dual magnetic flux
. the baryonic state must therefore be the dual magnetic vortex .
such a vortex in a nonsingular gauge has a vector potential of the form @xmath128 . to have a finite energy
it must be accompanied by the winding of the phases of both @xmath24 and @xmath93 .
in fact since @xmath93 carries @xmath88 times the charge of @xmath24 , the only states that are allowed energetically are those that carry @xmath88 vortices of @xmath93 .
this is natural in view of the ` dual ' relation between the effective and the fundamental theories .
the field @xmath93 is the ` vortex ' dual to the fundamental quark .
we thus expect that its elementary vortex would represent the fundamental quark itself , and so finite energy states must contain @xmath88 such elementary vortices .
a single vortex must be confined . in the same sense , @xmath24 is ` dual ' to the adjoint gluon field .
the elementary vortex of @xmath24 is then , in a sense , the ` constituent ' gluon .
in fact , such single gluon should not exist as a finite energy state either , and we expect it also to be confined .
the single vortex of @xmath24 should therefore bind either with the anti - vortex of the same type , or with @xmath88 vortices of @xmath93 .
the former type of state is a glueball , and exists in the pure glue theory  @xcite , while the latter type is a baryon .
interestingly enough , this line of reasoning leads us to expect that the baryon must have a bag - like structure .
namely , the quarks are bound to the @xmath24 field vortex . inside this vortex , the value of @xmath24 is small - in fact it exactly vanishes in the middle , and than rises quite slowly ( relative to the scale of @xmath129 ) towards the edges . recall that @xmath130 corresponds to a non confining state  @xcite .
the quarks are therefore sitting in the ` perturbative ' region of space - where there are no confining forces .
only when they separate far from each other - into the region with non vanishing @xmath24 , the linear potential pulls them inside again .
let us look at this more carefully .
consider for simplicity the @xmath60 symmetric case @xmath131 .
the baryon is the dual vortex with dual magnetic flux @xmath132 .
far from the vortex core , the field configuration is pure gauge , with the phases of @xmath24 and @xmath93 following the vector potential : @xmath133 the parameters of the model are such that the field @xmath24 is much lighter than both @xmath93 and @xmath134 .
thus the size of the vortex core of @xmath24 is large - of the order of the inverse glueball mass .
the two @xmath93-vortices which have a very small size core , sit inside this core .
since the length associated with the dual magnetic field is much smaller than the core size of the @xmath24 - vortex , the dual flux is concentrated on the @xmath93 - vortices . from the low energy point of view
, the picture is that two point - like magnetic vortices sit inside a soft core of a @xmath24 field vortex .
the field configuration looks roughly as depicted on figure 3 .
the magnetic flux is concentrated in the points @xmath135 and @xmath136 .
the phase of the field @xmath93 follows the variation of the vector potential very closely .
the phase of @xmath24 is also trying to do that , but it can not quite follow it all the way , since on the line between the two vortices it would have to be discontinuous . the most important energy contribution ( apart form the core energy of small vortices ) comes therefore from the vicinity of this line .
the phase of @xmath24 obviously has to interpolate across this line between the values @xmath48 and @xmath137 .
since the modulus of @xmath24 is not extremely rigid , it will be smaller along this line than in the immediate neighborhood .
both , the variation of the modulus and the phase of @xmath24 along the line connecting the small vortices contribute to the energy which is clearly linear in the distance @xmath138 .    [ fig3 ] .
arrows represent the direction ( phase ) of @xmath24 , and the big dots at @xmath135 and @xmath136 correspond to the positions of the @xmath93 field vortices.,title="fig : " ]      to study the structure of the baryon more quantitatively let us fix the gauge such that @xmath139\ ] ] @xmath140 this is a valid gauge choice for the configuration with two vortices of unit winding at the positions @xmath141 and @xmath142 .
due to the conditions on the parameters of our model , the vector field @xmath98 follows the phase of @xmath93 in the whole space .
thus @xmath143\ ] ] where : @xmath144 , @xmath145 , and @xmath16 , @xmath146 are the positions of the left and right vortices , respectively : @xmath147 , @xmath148 .
there are two interesting limiting situations .
the first is when the distance @xmath149 is larger than the dynamical distance scales @xmath150 and @xmath151 which determine the masses of the glueballs .
the second interesting situation is the reverse , that is when the two vortices are sitting well inside the glueball correlation length .
let us look at them in turn .
when the distance between the vortices is large we expect the potential between them to be linear with the string tension calculated in the theory without the dynamical @xmath93 field .
we will study the interquark potential in our low energy theory in the semiclassical approximation .    to find the minimal energy configuration
, we have to solve the classical equations of motion for the field @xmath24 at a fixed configuration of @xmath93 and @xmath134 given by eqs.([eq : uconf],[eq : bconf ] ) .
let us concentrate on the points which are close to the @xmath18 - axis , with @xmath152 .
the main contribution to the energy comes from this region of space . in this region
the vector potential @xmath134 vanishes , and the phase of the field @xmath93 is zero .
thus the equations of motion for the field @xmath24 are the same as in the pure gluodynamics .
also as long as we stick to this region , the configuration of @xmath24 is translationally invariant in the @xmath18 direction .
what determines the energy then are the boundary conditions on the field @xmath24 . in this configuration clearly
the phase of the field @xmath24 is @xmath153 far above the @xmath18 axis and @xmath154 far below the axis , fig.3 .
thus both , the equations and the boundary conditions are precisely the same as for the domain wall separating the two degenerate vacua in the effective theory of pure gluodynamics ( @xmath155 ) .
there is an extra contribution to the energy that comes from the region of space close to the points @xmath135 and @xmath136 . but
this energy does not depend on @xmath149 and is subleading for large @xmath149 .
the rest of the space does not contribute to the energy , since the field configuration there is pure gauge .
thus as expected the energy in this regime is @xmath156 , where @xmath157 is the domain wall tension ( fundamental string tension ) calculated in pure gluodynamics  @xcite .
it is perhaps more interesting to consider the opposite situation , that is when the distance between the quarks is smaller than the glueball correlation length .
this is the regime in which we do not expect to see any stringy structure .
instead we can ask whether the lowest energy configuration has any resemblance to a bag . to study this question we take the limit @xmath158 . in this case clearly the phase of the field @xmath24 will follow the phase of @xmath93 in the whole space @xmath159 given this condition , only the variation of the radial component of @xmath24 has to be determined .
since the problem has rotational symmetry , the equation of motion for the modulus @xmath10 becomes @xmath160 with @xmath161 , which is assumed to be positive throughout this paper .
the relevant boundary condition in this case is that at infinity @xmath10 approaches its vacuum value @xmath162 while in the vortex core it vanishes @xmath163 . with these boundary conditions , eq .
( [ eq : eqrho1 ] ) has the familiar form of the @xmath164 static kink equation in @xmath165 dimensions with the solution @xmath166 the energy of this solution is @xmath167_{a=0 } \;=\ ; \lambda \int_0^\infty \ , 2 \pi dr r \ , ( v^4 - \rho^4 ) \;=\ ; \frac{\pi^3}{12 } \ , \frac{\eta}{\lambda } \;=\ ; \frac{\pi^3}{12 } \ , ( \mu^2 - \frac{\xi u}{\lambda } ) \;.\ ] ]    the picture is thus just as described in the previous subsection .
since the two quark state is accompanied by the vortex of the field @xmath24 , the quarks effectively ` dig a hole ' in the vacuum . in their immediate vicinity
the modulus @xmath10 vanishes , and therefore there is a ` bag ' of the nonconfining state .
the radius of this bag is given by the mass of the scalar glueball @xmath168 .
it is interesting to note that , although for large separation @xmath149 the energy of the string gets contributions from both , the scalar and the pseudoscalar glueballs ( the modulus and the phase of @xmath24 ) , the ` bag constant ' is determined solely by the scalar glueball . for small @xmath149 the phase of @xmath24 is not excited and only @xmath10 deviates from the vacuum state inside the ` bag ' .
this is consistent with the common lore that the inside of the bag is distinguished from the outside by the value of the @xmath169 condensate .
in fact since @xmath10 has vacuum quantum numbers and interpolates in our effective theory the scalar glueball , it is naturally associated with the operator @xmath169 which has a large overlap with the scalar glueball in qcd .
in this note we have considered the effect of heavy fundamental quarks on the effective low energy theory description of nonabelian gauge theories in 2 + 1 dimensions .
we found that the status of the magnetic @xmath0 symmetry changes .
it becomes a local symmetry with the ` gauge coupling ' inversely proportional to the quark mass .
thus the low energy theory becomes a discrete gauge theory .
we have also studied the structure of the baryon in the framework of the effective theory .
the picture that emerges is very reminiscent to the ` bag ' .
the quarks in the baryon sit in the middle of the @xmath24-vortex , where the energy density differs from the vacuum energy density .
the ` bag constant ' is equal to the difference of these two energy densities , and parametrically in the effective theory is given by @xmath170 , where @xmath171 is the mass of the scalar glueball and @xmath172 is the expectation value of the vortex radial field @xmath10 .
an important thing to note is that the bag we are talking about here arises in a very different situation than in the usual bag model  @xcite .
there the bag describes the structure of the baryon containing _ light _ quarks .
the radius of the bag in this situation is determined by the balance of the vacuum pressure and the pressure due to the free motion of the light quarks inside , and in fact depends on the quark wave function . in our case
the inside of the bag contains heavy quarks .
their kinetic energy is small , and we have treated them here as static .
the radius of the bag thus is determined purely by the dynamics of the scalar glueball field and is not sensitive to the state of the heavy quarks .
this is true for low lying excitations for which the radius of the quark state is smaller than the inverse glueball mass .
when these two scales are comparable presumably the quark pressure will also be important and will play a role in the energy balance . thus in this intermediate regime we expect the @xmath24-vortex to be similar to the bag in the usual bag model . for states of even larger size
the potential between quarks is linear with the fundamental string tension .
the bag picture should therefore go smoothly into the string picture
. it would be interesting to study these questions further in the framework of the effective theory discussed in this note .
lastly we note that although we have so far considered scalar quarks , extending this discussion to spin @xmath173 quarks is fairly straightforward .
it is well known @xcite that in the abelian higgs model the spin of the vortex is controlled by the coefficient of the chern - simons term of the vector field . thus adding the chern simons term for the vector field @xmath98 @xmath174 with @xmath175 , will endow the vortex of the @xmath93-field with a half integer spin - vortex and not the @xmath24-vortex , since in the dynamics of our effective model the vector potential @xmath98 follows the phase of the field @xmath93 . ] .
our effective model with this extra term is relevant to the description of qcd with one flavour of heavy spinor quarks . just like the fermion mass term in 2 + 1 dimensions , the chern simons term breaks the parity symmetry in the baryonic sector of the theory .
the extra term does not affect the classical analysis we have preformed here , but is clearly relevant for determination of the spin of the baryon as well as the spectrum of the excitations .
this work was part of a collaboration supported by a fundacin antorchas - the british council grant .
c.  d.  f.  is a member of conicet , and is partially supported by anpcyt ( pict 97/1040 ) and fundacin antorchas grants .
a.  k.  is supported by pparc .
99 a.  m.  polyakov , _ nucl .
_ _ b120 _ , 429 ( 1977 ) . s.  r.  das and s.  r.  wadia , _ phys .
* b106 * , 386 ( 1981 ) , erratum - ibid . *
b108 * 435 ( 1982 ) ; n.  j.  snyderman _ nucl .
_ * b218 * , 381 ( 1983 ) .
g.  thooft , _ nucl .
phys . _ * b138 * , 1 ( 1978 ) .
s.  samuel , _ nucl .
phys . _ * b154 * , 62,(1979 ) ; a. kovner , b. rosenstein and d. eliezer , _ mod .
lett . _ * a5 * , 2661 ( 1990 ) ; _ nucl .
phys . _ * b350 * , 325 ( 1991 ) ; a. kovner and b. rosenstein , _ phys .
lett . _ * b266 * , 443 ( 1991 ) .
a.  kovner and b.  rosenstein , _ int .
j.  mod .
phys . _ * a7 * , 7419 ( 1992 ) .
a.  kovner and b.  rosenstein , _
j.  high  energy  phys . _ * 003 * 9809 ( 1998 ) . c.  p.
korthals - altes and a.  kovner , ` magnetic z(n ) symmetry in hot qcd and the spatial wilson loop , ' hep - ph/0004052 .
m.  j.  teper , _ phys .
_ * d59 * 014512 ( 1999 ) . l. krauss and f. wilczek _ phys .
lett . _ * 62 * , 1221 ( 1989 ) ; j. preskill and l. krauss _ nucl .
phys . _ * b341 * 50 ( 1990 ) .
a.  chodos , r.  l.  jaffe , k.  johnson and c.  b.  thorn , _ phys .
_ * d10 * ( 1974 ) 2599 ; a.  chodos , r.  l.  jaffe , k.  johnson , c.  b.  thorn and v.  f.  weisskopf , phys . rev .
* d9 * ( 1974 ) 3471 .
see , for example : a. lerda , _ anyons _ , lecture notes in physics , springer - verlag , berlin , 1992 .