the idea behind the sp was put forward by efetov in 1979 in the context of layered superconductor with parallel magnetic field @xcite . it was suggested that the field can suppress the inter - layer josephson coupling so that the low energy properties of this 3d system can be described as being essentially of 2d character .
later , in ref . @xcite it was shown that the frustration due to the magnetic field is not sufficient to fully suppress the coupling .    in the context of quantum 1d chains ( equivalent to 2d classical layers by the virtue of the quantum to classical mapping ) the possibility of the decoupling between chains
has also been explored @xcite . the main argument for such a decoupling
is based on the scaling dimensions of the josephson coupling determined with respect to the luttinger liquid parameter in each chain : if it is larger than 2 , the coupling should become irrelevant @xcite .
these proposals have been criticized in refs.@xcite where it was shown that the inter - chain tunneling is always relevant .    further argumentation in favor of the sp and , actually , the name _ sliding phases _ have been proposed in refs.@xcite where the inter - layer gradient couplings between classical xy variables in each layer have been considered in addition to the josephson one .
such gradient terms can independently control the scaling dimensions of the josephson coupling and of the vortex fugacity in each layer so that the first one can become irrelevant above some temperature @xmath0 ( of the _ dimensional reduction _
@xcite ) which is below the temperature of the berizinskii - kosterlitz - thouless ( bkt ) transition in the layers .
thus , there is a range of temperatures where the sp are supposed to exist .
this approach was also developed for the case of quantum 1d luttinger liquids coupled by both the josephson and the gradient terms @xcite which are the analog of the andreev - bashkin drag effect @xcite .
it is important to note that the proposal of sp is based on applying the rg logic to compact variables characterized by global u(1 ) symmetry . while these early suggestions were more of a purely academic interest ,
expanding capabilities of ultra - cold - atoms techniques in recent years emphasize the importance of these suggestions especially in the context of possible new phases in composite lattices @xcite and in the presence of disorder @xcite . in more general terms
, the question is if it is possible to realize a phase transition from a low- to higher- dimensional behavior .    here
we will analyze a simplest classical xy system characterized by the gradient interactions and the josephson coupling @xmath1 .
the gradient terms are chosen in such a way that the sp is supposed to exist in some range where the renormalized value @xmath2 of @xmath1 scales to zero as layers size @xmath3 grows .
we will present the results of the large scale monte carlo simulations of this system .
our analysis is based on the dual formulation of the model  in terms of the closed loops .
the main result is that no sp state exists in such a system .
instead , the value of @xmath2 is always finite .
this behavior will be compared with the standard asymmetric xy layered model where no sp are expected to occur .
we will also derive the analytical result for @xmath2 in the asymptotic limit when the intra - layer stiffness is much larger than @xmath1 .
the numerical results have been found to be consistent with the analytical solutions for both models .
our paper is organized as follows . in sec.[sec : i ] we introduce the bilayer model and provide the rg solution for sp .
then , we construct the dual representation in sec . [
sec : dual ] . using the duality we have found the asymptotic analytical solution for the renormalized josephson coupling @xmath2 in sec .
[ sec : as ] .
the monte carlo simulations of the bilayer model are presented in sec.[sec : num ] .
then , in sec.[sec : nz ] we present the results on a stack of bilayers along the same lines as for the bilayer . finally , in sec .
[ sec : dis ] we discuss the implications of our analytical and numerical results and also provide an alternative model for the sp .
here we introduce a model of two asymmetric parallel layers , each being a square lattice of linear size @xmath4 ( in terms of the inter - site shortest distance ) characterized by two fields @xmath5 and @xmath6 on the layers @xmath7 , respectively.the action can be written as h&= & - _ ij [ t_1 ( _ ij _ 1 - a_ij ) + t_2 ( _ ij _ 2 -g_2 a_ij ) + & + & a^2_ij ] - _ i u(_2(i)-_1(i ) ) [ 2n ] where @xmath8 and @xmath9 are parameters ; @xmath10 denotes summation over nearest neighbor sites within each layer ; @xmath11 ; @xmath12 is a bond vector field ( that is , @xmath13 ) oriented along the bond @xmath10 .
it is introduced in order to generate the `` current - current '' interaction ( cf .
@xcite ) consistent with the compact nature of the fields @xmath14 .
this action is to be used in the partition function z = da d_1 d_2 ( - h ) [ zz2 ] where the temperature is absorbed into the the parameters @xmath15 .
our focus is on verifying the applicability of the rg analysis to the renormalization of the josephson coupling @xmath1 .
hence , we will not discuss physical origins of the variables and the parameters .      in the approximation ignoring compact nature of the variables ,
the terms @xmath16 and @xmath17 are replaced by @xmath18 and @xmath19 , respectively .
then , the gaussian integration over @xmath12 can be carried out explicitly in eq.([zz2 ] ) , so that ( [ 2n ] ) in terms of the remaining variables becomes h_0= _ ij k_ab _ ij _
b - _ i u(_2-_1 ) , [ 2n0 ] where the matrix @xmath20 is related to the original parameters as k_11&= & , k_22= , + k_12 & = & - .
the stability of @xmath21 is guaranteed by k_11k_22 - k_12 ^ 2= > 0 .
[ stab ] the condition for sp can be obtained along the lines of the logic @xcite which ignores the compactness of @xmath14 .
specifically , introducing the variables @xmath22 and @xmath23 and , then , integrating out @xmath24 , the resulting partition function becomes z_0=d^-h _ , h
_ = d^2 x , [ z_0 ] where the notation k= [ k ] is introduced and the long wave limit is considered  so that the summation along the layers is replaced by the integration @xmath25 .
as long as the compactness of @xmath26 is ignored , eqs .
( [ z_0 ] ) represent the standard sine - gordon model in 2d .
the rg analysis predicts ( see in , e.g. , @xcite ) that at @xmath27 the renormalization renders the josephson coupling @xmath1 irrelevant in the thermo - limit @xmath28 .
more specifically , the renormalized @xmath1 should flow to zero as @xmath29 .
such a behavior is supposed to occur together with the persistence of the algebraic order along the planes .
this requirement imposes further restrictions on the values of @xmath30 .    without loss of generality
let s assume @xmath31 and introduce the notations : @xmath32 as a measure of temperature , and @xmath33 .
then , the condition @xmath34 for sp becomes t > t_d=. [ ssp ]    in order to guarantee the algebraic order in each layer no bkt transition should occur in the layers . in order to determine possible types of vortices responsible for the transition
, we examine the form ( [ 2n0 ] ) by `` reinstating '' the compactness of the variables in the limit @xmath35 ( which is supposed to renormalize to zero ) . in the presence of the gradient coupling a vortex
may have a composite structure @xcite : @xmath36 circulations in the component 1 may be bound to @xmath37 circulations in the component 2 . finding the condition for the proliferation of such composite vortices can be achieved along the line of kosterlitz - thouless argument developed for simple xy - model in 2d .
the free energy of such a complex is f_v= l -2 l. [ vort ] then , @xmath38   - 2 t $ ] , and the stability against the bkt transition is guaranteed by the positivity of @xmath39 or t < t_(q_1,q_2)= [ ( q_1 + xq_2)^2 + ( y - x^2)q^2_2 ] , [ tbkt ] where the minimization with respect to @xmath40 must be performed .
proliferation of simple vortices @xmath41 corresponds to @xmath42 .
it is , however , not always a minimal value as long as @xmath43 .
let s also note that , since @xmath44 by definition , there are no solutions for @xmath45 if @xmath46 .
let s look for a solution when @xmath47 is integer , that is , @xmath48 . then , keeping in mind the stability requirement ( [ stab ] ) , that is , @xmath49 , the minimal vortex corresponding to @xmath50 has free energy lower than that of the simple vortex if @xmath51 . in this case ,
the solution for @xmath52 exists if @xmath53 .
proliferation of the composite vortex pairs corresponds to disordering of the original fields @xmath54 , while the composite field @xmath55 remains ( algebraically ) ordered .
this mechanism constitutes the formation of thermally induced bound phases ( or using the language of superfluidity  _ thermally paired superfluid _ , tps , @xcite ) . in other words
, the system behaves in a such a way that the algebraic order persists only in the composite field @xmath56.[since @xmath57 we call such a composite phase as thermally bound superfluid ( tbs ) by analogy with the tps ] .
this effect does not require that @xmath47 is necessarily integer .
if @xmath47 is non - integer , its closest integer part will determine the power of @xmath58 .
if @xmath59 , the lowest energy vortex is the simple one , and the solution for the inequality @xmath45 also requires @xmath53 . for @xmath60 , @xmath61 while @xmath62 as long as @xmath63 is kept constant .
such a limit corresponds to the largest range of @xmath64 where sp are to be anticipated for the two - layer model .
however , for practical purposes of simulations using too large @xmath47 leads to slower convergence .
thus , we choose @xmath65 corresponding to a reasonably wide range where sp is anticipated to exist . taking into account eq.([k ] ) , the condition for the sp can be written as < t < .
[ spt ] keeping in mind the chosen values @xmath66 , this becomes @xmath67 or @xmath68 .
the results of the simulations will be conducted at the `` temperature '' @xmath64 in the middle of the interval @xmath69 , that is , @xmath70 .
more specifically , @xmath71 . in terms of the original parameters @xmath72 , the relations @xmath73 , @xmath74 , @xmath75 and @xmath76 must be satisfied .
[ such a fine tuning is of no concern , since the focus is on the consistency of the paradigm of sp rather than on physical origin of the values ] .    concluding this section , the presented analysis based on the rg finds that it is possible to find a range of temperatures where the sequence of phases is as presented in fig .
[ pd ] in the panel ( a ) : at @xmath77 the josephson coupling is relevant .
we call this range as josephson phase in fig .
[ pd ] . at @xmath78
there is the sp where the symmetry @xmath79 is promoted to u(1)@xmath80 u(1 ) . in the range @xmath81 the tbs phase
is characterized by the composite field @xmath82 .
thus , the broken symmetry is partially restored through the subgroup @xmath83 , where @xmath84 $ ] . at higher temperatures @xmath85
the composite field @xmath82 becomes disordered .
-8 mm    -5 mm      the partition function @xmath86 can be evaluated by the high - temperature expansion method ( see e.g. in @xcite ) in terms of @xmath87 and the explicit integration over the variables .
this approach allows obtaining @xmath86 in terms of the integer bond variables  powers of the corresponding taylor series .
since the resulting configurational space consists of closed loops of the bond currents , further simulations can be effectively performed by the worm algorithm @xcite .
as will be also shown , the language of loops also allows obtaining analytic expressions for the renormalized josephson coupling @xmath2 which are exact in the asymptotic limit .
we will be utilizing the villain approximation @xcite for the cosines to obtain the so called j - current version @xcite of eqs.([zz2]),([2n ] ) : z=_\{m_a , ij , m_i } dda e^-h_v , [ vill ] h_v&= & _ ij [ ( _ ij _ 1 - a_ij + 2m_1,ij)^2 + & + & ( _ ij _ 2 -g_2 a_ij + 2m_1,ij ) ^2 + a^2_ij ] + & + & _ i ( _ 2(i)-_1(i ) + 2m_i)^2 , [ 2nv ] where @xmath88 ( @xmath89 ) are integer numbers defined along bonds between two nearest sites @xmath90 and @xmath91 along the planes and @xmath92 is an integer assigned to a site @xmath90 and oriented from the layer 1 to the layer 2 .
the villain approximation proves to be very accurate for establishing the transition points as well as in general if the effective constants @xmath93 are properly expressed in terms of the corresponding bare values @xmath87 ( see in ref.@xcite ) .
the `` renormalization '' can be essentially ignored for @xmath94 , so that in what follows we will be using @xmath95 .
similarly , for the josephson coupling @xmath96 one should take @xmath97 and , if @xmath98 , the corresponding relation is @xmath99 @xcite .
after using the poisson identity for each integer and performing the integrations over @xmath100 and @xmath101 , the resulting expression becomes z=_\{j_1,ij},\{j_2,ij } , \{j_z , i } e^-h_j , [ z2 ] h_j= _ ij ; a , b ( k^-1)_ab j_a , ij j_b , ij + _
i j_i^2 , [ h_j ] where @xmath102 labels layers and @xmath103 is the matrix inverse to @xmath30 introduced in eqs.([2n0]),([ktg ] ) ; the summation runs over the integer bond currents @xmath104 ( oriented from site @xmath90 to site @xmath91 so that @xmath105 ) within each corresponding layer 1,2 as well as over the integer currents @xmath106 oriented along the bond connecting the site @xmath90 in the layer 1 to the site @xmath90 in the layer 2 .
all the configurations are restricted by the kirchhoff s current conservation rule
the total of all j - currents incoming to any site must be equal to the total of all outcoming currents from the same site .
it is useful to note that eqs.([z2],[h_j ] ) can also be obtained directly from eq.([2n0 ] ) by reinstating the compact nature of the variables : @xmath107 and @xmath108 so that the matrix @xmath109 is viewed as being independent from the parameters in the action ( [ 2n ] ) .
the system ( [ z2],[h_j ] ) features statistics of closed loops .
if @xmath35 , there are two sorts of loops  one in each layer .
thus , each configuration is characterized by definite values of the windings @xmath110 in the @xmath111th layer along the @xmath112 directions of the planes .
it is straightforward to show that statistics of these windings determine the renormalized values @xmath113 of the matrix @xmath30 along the line of the approach @xcite . more specifically _
ab = _ = , w_a , w_b , .
[ kr ] this expressions are valid for periodic boundary conditions ( pbc ) .
it is important to note that @xmath113 represents an exact linear response ( at zero momentum ) with respect to thouless phase twists .
in other words , if there are externally imposed infinitesimal constant gradients @xmath114 ( violating the pbc ) of the phases @xmath14 , the free energy acquires the contribution @xmath115 . on the other hand , in the presence of the gradient the integrand of the partition function gets the factor @xmath116 .
comparing both expressions leads to the relation ( [ kr ] ) .    as a test of consistency ,
we have checked numerically that in the regime where the sp state is supposed to exist ( that is , @xmath117 ) , the deviations of @xmath113 from the bare values @xmath30 are within the statistical error less than 1% for all tested sizes of the layers @xmath118 .
significant deviations are observed only as the system approaches fully disordered state  that is , where the fields @xmath119 as well as the composite one @xmath82 become disordered . in this case , @xmath113 flow to zero as @xmath3 increases .
the deviations remain small ( about 2 - 3% ) even in the regime where @xmath82 is the only ordered field .
the emergence of the tbs is detected by observing that windings @xmath120 in the layers 1 and 2 are changing exactly by the increment @xmath121 , @xmath122 ( plus or minus ) , respectively .    at finite values of @xmath1
the loops belong to both layers so that no separate windings can be introduced .
however , the sums @xmath123 remain well defined and can be used to evaluate the rigidity @xmath124 of the fields along the layers . in a general case of @xmath125 symmetric ( with respect to the @xmath126 directions ) layers
@xmath127 : & = & _ w^2 _ [ stif22 ] + w _ & = & _ ij , a=1,2, ... n_z j_a , ij , [ ww ] where for a given @xmath112 in ( [ ww ] ) the bond @xmath128 ( as well as @xmath129 ) is oriented along the direction @xmath130 .    our focus here on the renormalized value @xmath2 of the josephson coupling @xmath1 in the sp regime .
if the periodic boundary conditions are also imposed perpendicular to the layers ( along @xmath131-direction ) , the inter - layer response @xmath2 is given by windings @xmath132 along @xmath131-direction : u_r = w^2_z , w_z = _ i j_i , [ wu ] where the summation @xmath133 of the currents @xmath134 ( oriented along @xmath131-direction ) is performed over all sites of all layers .
similarly to the cases ( [ kr ] ) and ( [ stif22 ] ) , eq.([wu ] ) represents the full linear response at zero momentum  that is , the renormalized value @xmath2 of the josephson coupling @xmath1 .    at this point
, we should comment on how to interpret the pbc for two layers , @xmath135 . while in the case @xmath136 it is a natural procedure to link the @xmath137th layer to the first one , @xmath138 , by the josephson term , the case @xmath135 needs an auxiliary construction because the layers 1 and 2 are coupled already directly .
the formal procedure , then , consists of adding a third layer , @xmath139 , with no rigidity along @xmath126 directions and coupled by the josephson term to both layers , @xmath7 .
if the coupling @xmath140 between the layers 1 and 3 and the coupling @xmath141 between the layers 2 and 3 add up as @xmath142 , in the dual action ( [ h_j ] ) the sum in the last term can be extended to the layers @xmath143 , while the first term is still confined to the layers @xmath7 .
the key to this procedure is the kirchhoff s rule : the j - current from a site @xmath144 along @xmath131-direction from the layer 2 to the layer 3 must be exactly the same as the current from the site @xmath144 in the layer 3 to the layer 1 . then , in the form ( [ h_j ] ) the same value @xmath145 can be used for the currents from the layer 1 to the layer 2 directly or through the layer 3 .
as mentioned above , the dual representation allows obtaining analytically the asymptotic values for @xmath2 .
let s begin with the trivial case of zero stiffnesses @xmath30 and arbitrary number of layers , @xmath146 . in this case
, the action ( [ h_j ] ) trivially becomes h_a= _
i j^2_i , j_i=0,1 , 2 , ... [ a ] where the summation runs over all sites @xmath90 of only _
one _ layer , say , @xmath138 . in this expression
the kirchhoff rule dictates that the current @xmath134 at a given site along @xmath131-direction must be the same for all values of @xmath131 .
thus , such a current with @xmath147 constitutes one closed loop characterized by the winding @xmath148 .
this allows constructing the partition function exactly as z_a=^l^2 [ za ] where @xmath149 is the number of sites in one layer .
the stiffness ( [ wu ] ) can be found by taking into account that the total winding along z - direction is @xmath150 , where the summation runs over @xmath149 sites of only one layer .
then , using statistical independence of different sites we find u_r=. [ ur1 ] this expression shows that , as long as @xmath125 is finite , the josephson coupling remains relevant even if there is no in - plane order .    in the limit @xmath151
only the term @xmath152 is important , so that eq.([ur1 ] ) becomes u_r= .
[ urr ]    eq.([urr ] ) can be used even in the case when there is finite stiffness along the layers , with @xmath125 substituted by some effective value @xmath153 , that is , u_r= . [ genm ] this is easily justified in the `` ideal gas '' approximation of rare fluctuations of the j - currents in z - direction . in the case of @xmath154
the asymptotic asymmetric limit corresponds to @xmath155 ( so that the system is well above the disordering transition ) and @xmath156 .
the loop proliferation can be viewed from the perspective of the worm algorithm @xcite where one open end of a string of j - currents walks randomly until it meets another open end so that the loop is formed .
then , the most of the path is residing in a layer with only occasional jumps between neighboring layers .
these jumps are controlled by the exponential smallness of @xmath157 .
thus , the full action can be well approximated by @xmath158 , eq.([a ] ) , with @xmath159 .
this leads to eq.([genm ] ) with @xmath160 : u_r= .
[ urxy ] later we will present the numerical evidence that the stiffness perpendicular to the layers of a simple xy layered model in the asymptotic limit can well be described by the above equation .
below we will show that for the case of the two asymmetric layers , the effective value of @xmath153 is @xmath161 in eq.([genm ] ) .      here
we present the results of monte carlo simulations of the bilayer in the regime of sp in the limit @xmath156 .
the action ( [ h_j ] ) can be represented in the notations @xmath162 as h_j= _
ij+ _ i , [ h12 ] where @xmath163 and @xmath164 refer to the inplane bond currents in the layers 1 and 2 , respectively ; the actual values of the parameters used in the simulations have been discussed at the end of sec.[sec : rg ] : @xmath165 .
the structure of the loops is determined by the energy of creating a j - current element along a given direction .
a typical energy to create an additional j - current element in the plane 2 can be estimated as @xmath166 .
thus , large loops with a typical values @xmath167 can exist in the plane 2 .
in contrast , the energy to create an isolated element in the plane 1 ( with no @xmath168 currents along the same bond in the layer 2 ) costs much larger energy : @xmath169 .
thus , the probability to create such an element is exponentially suppressed as @xmath170 , and no entropy contribution ( due to 4 optional directions along the plane ) can compensate for such a low value .
this implies that no large isolated loops can exist in the layer 1 .
the only option to create a large loop in the layer 1 is if each element @xmath163 is coupled to currents @xmath171 along the same bonds in the layer 2 .
a typical energy of this combined element is @xmath172 .
this strong asymmetry between the layers has immediate implication for the windings along @xmath131-direction  the minimal length @xmath153 of the element @xmath134 must be @xmath161 in eq.([genm ] ) .
thus , the stiffness @xmath2 in the limit @xmath98 becomes u_r = 4e^-1/u_v= u^2 , [ ur3 ] where the asymptotic expression @xmath173 @xcite has been used .
the results of the simulations is shown in fig.[fign2 ] .
-8 mm    -5 mm    the first striking feature to notice is that @xmath2 does not depend on the layers size @xmath3 over 2 orders of magnitude of @xmath3 and over 7 orders of @xmath145 ( which is actually @xmath174 of the bare coupling ) .
second , the numerically found value @xmath2 follows the analytical result ( [ ur3 ] ) with high accuracy  even for values @xmath175 .
both features are in the striking conflict with the rg prediction stating that @xmath2 should decay as @xmath176 in the sp regime ( @xmath177 ) .
it should be also noted that the stiffness along the layers ( [ stif22 ] ) remains finite and much larger than @xmath2 , that is , @xmath178 for all simulated sizes from @xmath179 to @xmath180 .
this justifies the validity of eq.([ur3 ] ) even in the case @xmath175 .
as it became evident from the previous analysis , no sp can occur in the double layer . referring to the sketch of the possibilities , fig .
[ pd ] , the option ( b ) is realized . here
we will address a possibility of sp in a @xmath125-layers setup .
in other words , we will be looking for a behavior where the renormalized stiffness ( the inter - layer josephson coupling ) @xmath2 decays as a function of @xmath125 in the limit @xmath182 , while the stiffness along planes remains finite.[this would be a `` weaker '' version of the sp ] .
we consider the pbc setup : the odd @xmath183 and the even @xmath184 layers are characterized by the inplane stiffnesses @xmath185 and @xmath186 , respectively , with the nearest layers coupled by the current - current term @xmath187 ( the same for all pairs of layers ) as well as by the josephson coupling @xmath188 , where @xmath189 is the xy variable defined on a site @xmath144 belonging to the layer @xmath131 .    in the linearized with respect to @xmath190 approximation analogous to eq.([2n0 ] ) the model becomes h_0&= & _ z=1,3,5, ... \ { _ ij h_z;ij + & - & u _ x , y [ ( _ z+1-_z ) + ( _ z-1-_z ) ) ] } [ nnz ] where the summation runs over odd values of @xmath131 and the notation h_z;ij = ( _ ij _ z)^2 + ( _ ij _
z+1)^2 + + k_12_ij _ z(_ij _ z+1+_ij _ z-1 ) [ hnnz ] is used .
if the compact nature of @xmath191 is ignored , the one - loop rg flow equation for @xmath2 reads = ( 2- ( _ z+1 - _ z)^2 _ s ) u_r [ rgn ] where @xmath192 refers to the rg shell integration over the inplane momenta @xmath193 , with @xmath194 being the cutoff and @xmath195 .
we note that , due to the pbc along @xmath131-direction , the mean @xmath196 does not depend on @xmath131 .    using discrete fourier representation along @xmath131 direction with doubled unit cell containing two layers ( the odd and the even ) with two sorts of phases @xmath197 and @xmath198 along odd and even layers , respectively , the part @xmath199 can be diagonalized and the correlator in eq.([rgn ] ) found .
this gives eq.([rgn ] ) rewritten as = 2(1- ) u_r , [ rgn2 ] where t^-1_d=^(nz/2)-1_m=0 , [ tdnz ] and the wavevectors along @xmath131 take values dictated by the periodic boundary conditions @xmath200 .
here we use the same notations @xmath201 introduced in sec.[sec : rg ] .
thus , at @xmath177 rg predicts irrelevance of @xmath2 .    the upper limit on @xmath64 can be obtained from the requirement of no free 2d vortices in the limit @xmath35 .
there are , actually , two options : i ) looking for a composite vortex characterized by phase windings @xmath36 and @xmath37 in odd and even layers , respectively , forming a string of length @xmath125 perpendicular to the layers ; ii ) considering independent vortices @xmath202 only in odd layers ( characterized by smallest stiffness @xmath185 ) and characterized by the temperature @xmath203 . as the analysis shows , the option when composite vortices form finite strings ( say , @xmath204 in layers @xmath205 and @xmath206 in the layer @xmath207 by analogy with the @xmath135 case ) costs larger energy than in the case ii ) .
the option i ) is characterized by the vortex energy @xmath208\propto n_z$ ] .
the minimum for a @xmath209 integer is achieved at @xmath210 .
thus , in order to compensate for the factor @xmath211 , the system must be very close to the instability @xmath212 .
simulations in this region turn out to be problematic .
thus , we will conduct simulations in the range @xmath213 , provided the condition @xmath214 holds in the limit where @xmath215 remains finite for @xmath216 . specifically , the condition @xmath217 reads ^nz/2 -1_m=0 > 8 .
[ tdt ] it can surely be achieved for large enough @xmath47 in the limit @xmath218 .
replacing the summation by integration in this limit and considering @xmath219 , eq.([tdt ] ) gives @xmath220 .
for the simulations we have chosen @xmath221 and @xmath222 , which gives @xmath223 with @xmath224 chosen in the middle of the interval between @xmath225 and @xmath0 .
the chosen value of @xmath0 corresponds to the limit @xmath226 , and for any finite @xmath125 , the actual @xmath0 from eq.([tdnz ] ) is below this value .
at this point we note that for any finite @xmath1 the system is 3d and , strictly , speaking the notion of the bkt transition becomes inadequate for large enough @xmath125 : even at @xmath227 the odd layers would still have xy order due to the proximity to even layers . here , however , we presume that sp scenario holds and @xmath2 vanishes at large @xmath3 . practically , simulations are performed at finite @xmath1 and we always control that the xy stiffness ( helicity modulus ) along the layers remains finite and independent of @xmath3 .
the dual formulation in terms of the closed loops of integer j - currents ( along bonds in and between the layers ) can be achieved similarly to the case @xmath135 by reinstating the compactness of @xmath191 in eqs.([nnz],[hnnz ] ) through the villain approach : @xmath228 along the planes and @xmath229 for josephson coupling , where @xmath230 refers to an arbitrary ( oriented ) integer defined on the bond @xmath231 belonging to the plane @xmath131 and @xmath232 stands for an integer on a bond connecting site @xmath90 in the plane @xmath131 to the same site in the plane @xmath233 .
the partition function is obtained as a result of integration over all @xmath234 and summations over all bond integers .
the j - currents enter through the poisson identity @xmath235 applied to each bond integer .
this allows explicit integration over all phases @xmath191 as well as over the bond integers @xmath236 .
there are two types of j - currents : inplane @xmath237 within each `` elementary cell '' ( along @xmath131 ) and between the planes @xmath238 .
the label @xmath239 refers to j - current defined on the bond @xmath231 belonging to a plane with odd @xmath131 .
accordingly , @xmath240 stands for the inplane current on the layer with even @xmath131 .
then , @xmath238 denotes the current from the site @xmath90 from the plane @xmath131 to the plane @xmath233 .
the integration over phases @xmath241 generates the kirchhoff constraint  similarly to the bilayer case .
finally the ensemble can be represented as z&= & _ \ { } ( - h_j ) , + h_j&= & _ ij ; z , zv_ab(z - z ) j^(a)_z , ijj^(b)_z,ij + & + & _ i , z j^2_z , i , [ hdual ] where the matrix @xmath242 is defined in terms of the matrix @xmath30 .
it reflects the asymmetry between odd and even layers .
explicitly , @xmath243 , for @xmath244 being even , describes the interaction between odd layers , and @xmath245 is defined between even layers ; @xmath246 $ ] refers to the interaction between odd and even layers ( that is , @xmath131 is odd ) , and v_22(z)= _ q_m , [ vzz22 ] with @xmath247 and the summation running over @xmath248 .
the dual formulation for @xmath125 layers allows obtaining the asymptotic expression for @xmath2 within the same logic used for deriving eq.([ur3 ] ) .
we will repeat it here .
the loop formation can be viewed as a process of random walks of two ends of a broken loop
exactly along the line of the worm algorithm @xcite .
such a walk of each end is controlled by energetics of creating one bond element @xmath249 in a randomly chosen direction  either along a given plane or perpendicular to it . very similar to the case of the two layers , the energy to create such an element alone along the odd layer costs energy @xmath250 , while the same element along an even layer costs energy @xmath251 .
the only option for creating a loop in an odd layer is if its energy is compensated by parallel elements in the even plane .
this feature is caused by the strong current - current interaction @xmath252 .
thus , if the walk occurs along @xmath131-direction from some even layer @xmath131 toward the neighboring odd layer @xmath233 , the subsequent move along the odd layer will be too energetically costly so that the walker would either move further toward @xmath253 layer or will go back to the original layer @xmath131 .
thus , the inter - layer elements are characterized by either @xmath254 or @xmath255 .
the weight of such a process is either @xmath256 or @xmath257 , respectively .
even if the walker makes a step or two along the layer @xmath233 ( which is a highly improbable event ) and then chooses to go toward the layer @xmath253 , the contribution to the partition function will be further reduced exponentially by the energy of the element @xmath258 along the odd plane .
thus , such processes can be ignored , and we arrive at the conclusion that @xmath2 given by eq.([ur3 ] ) must be valid for arbitrary @xmath125 in the asymptotic limit .
the model ( [ hdual ] ) has been simulated by the worm algorithm @xcite .
the renormalized inter - layer stiffness @xmath2 was found for a wide range of layer sizes , @xmath259 and @xmath260 .
the resulting data is presented in figs .
[ fig3],[fig4 ] . as can be seen in fig .
[ fig3 ] , the solution ( [ ur3 ] ) plays the role of the envelop for the family of the curves @xmath2 vs @xmath261 for various @xmath3 and @xmath262 .
we note that the stiffness @xmath263 along the layers ( as determined by eq.([stif22 ] ) ) remains independent of the sizes and much larger ( @xmath264 ) than @xmath2 .
this justifies the applicability of the asymptotic limit for eq.([ur3 ] ) .
we have also controlled that the system is far enough from any possible composite phases @xcite state by measuring the lowest order correlator @xmath265 and observing that it exhibits long - range order for values @xmath266 in the limit @xmath267.[in the composite phase state such a correlator is short ranged ] .
thus , the system is well in the putative sp state .
its behavior , however , is drastically different from the rg prediction .
-8 mm    -5 mm    -8 mm    -5 mm    -8 mm    -5 mm    -8 mm    -5 mm    at this point we should discuss the deviations of the numerical curves from the analytical result seen in fig.[fig3 ] . to some extent this behavior could have been interpreted as the evidence of sp .
there is , however , one important observation : the value of @xmath268 below which such suppression begins _ decreases _ as ( u^*)^-1 = ( l^2/n_z ) , = 1.000.02 [ ugm2 ] for @xmath269 in the main logarithmic approximation .
this behavior is demonstrated in fig.[fig4 ] , where the value @xmath270 corresponds the offset for @xmath2 being taken at 1/10 of the value given by the analytical expression ( [ ur3 ] ) .
the suppression of @xmath271 in the limit @xmath28 indicates that it is a finite size effect . in other words , starting from small @xmath3 at some @xmath272 , the renormalized stiffness @xmath2 can be well below the asymptotic limit ( [ ur3 ] ) .
as @xmath3 increases while @xmath145 is kept fixed , @xmath2 eventually approaches the limit ( [ ur3 ] ) .
this feature is clearly in a stark contrast with the rg prediction of sp where @xmath2 is supposed to flow to zero for fixed @xmath145 as a power of @xmath182 .
the deviation from the limit ( [ ur3 ] ) has a very simple explanation : it is essentially a consequence of the generic exponential suppression of any order in a quasi-1d limit @xmath226 .
[ there is no such suppression in the case of @xmath135 because the loops along z - direction are independent from the inplane ones ] . in the context of our system
it can be interpreted as the effect of zero modes ( in each layer ) fluctuations .
indeed , excitations along the planes renormalize @xmath145 to @xmath2 at short distances
. then , as long as @xmath273 the only remaining lowest energy degrees of freedoms are zero modes , that is , excitations with @xmath274 being independent from the @xmath126 positions along the planes ( and dependent on @xmath131 ) .
then , the effective action becomes s_0= -u_rl^2 _ z[(^(0)_z+1 - ^(0)_z ) ] .
[ zero ] it should be used in the partition function @xmath275 .
its analytical evaluation gives u^(0)_r ~ ( -n_z/(u_rl^2 ) ) [ zero ] for the value of the stiffness along @xmath131-direction  that is , the renormalised @xmath2 on large scale .
keeping in mind that at short scales @xmath276 for @xmath156 , one arrives at the relation ( [ ugm2 ] ) with @xmath277 .
it corresponds to the requirement @xmath278 for some @xmath268 so that eq.([zero ] ) describes the exponential suppression at @xmath279 .
clearly , such a quasi-1d suppression is also present in the standard xy model ( where no sp are anticipated to exist ) . in order to demonstrate this
explicitly we have also simulated a simple xy model given by the system z_xy&=&d _ z ( -h_xy ) , [ hxy ] + h_xy&=&- _ ij , z [ ( _ ij _ z ) + u ( _ z _ z ) ] , with some @xmath280 ( guaranteeing that no bkt transition occurs in each layer for @xmath35 ) , and @xmath281 . in the dual representation this system is described by h_xy _
xy = _ ij , z j_z , ij^2 + _ i , z j_i , z^2 , [ hxy2 ] where @xmath282 and @xmath238 are the same j - currents introduced above for the model ( [ hdual ] ) .
the results of the simulations of this model are presented in fig.[fig1],[fig2 ] . in the asymptotic limit
the inter - layer stiffness is described by eq.([urxy ] ) .
then , according to the above discussion the value @xmath270 determining the start of the deviations is given by @xmath283 , that is , with the slope @xmath284 which should be compared with the numerical value @xmath285 in fig.[fig2 ] .
thus , both models demonstrate essentially the same 3d behavior .
the only difference is the slope of the renormalized josephson coupling @xmath286 vs its bare value @xmath145 .
it is determined by the minimal length of the elementary j - current in the direction perpendicular to the layers .
the rg approach to 2d systems proves to be every effective in many cases including 2d xy model when it can be mapped on the sine - gordon ( sg ) one @xcite .
a successful implementation of the rg analysis to the josephson coupling was demonstrated in ref.@xcite where a single weak link can make one channel luttinger liquid insulating .
the merit of rg , however , should be taken with caution when applied to the dimensional reduction situations in layered systems . in this case
there is no exact mapping between xy and sg representations at finite inter - layer josephson coupling , and the approximation ignoring the compact nature of the variables becomes uncontrolled .
as our analysis of one particular layered system shows , no sp exists in such a system despite the rg prediction : the system shows essentially the 3d behavior of the asymmetric xy model .
while it is not obvious to us what is wrong with treating josephson coupling between 2d layers by rg in the original representation of fields , the dual formulation in terms of the closed loops gives a very important insight .
specifically , the sp means that as layer size @xmath28 , a suppression of the josephson coupling between layers would require that the number of times elements of closed loops fluctuate between layers must scale slower than @xmath149 so that the density of such events is zero in the limit @xmath287 .
the loops statistics , however , is controlled by local energies of creating finite elements and the entropy due to 6 directions in 3d vs 4 along layers .
thus , as long as there is a finite energy to cross between neighboring layers , the entropy will lead to a finite density of crossings for large enough @xmath3 .
similar argument can be applied to quantum wires in terms of the quantum to classical mapping where imaginary time is treated as an extra dimension .
the dual approach and the argumentation along the line of the numerical algorithm @xcite , treating closed loops formation as a process of the worm head wondering around and eventually finding its tail , allowed us to expose the actual stages of the renormalization of the josephson coupling : i ) at small scales josephson coupling is controlled by exponentially suppressed random and independent ( in the asymptotic limit ) events of crossings between layers .
it can be viewed as an ideal gas of j - currents between the layers .
this stage leads to the renormalized coupling , in general , represented by eq.([genm ] ) with @xmath288 .
ii ) if the number of layers @xmath125 increases , with @xmath3 being fixed , quasi 1d fluctuations further suppress the coupling exponentially as demonstrated in eq.([zero ] ) .
here we have discussed a local model characterized by short range interactions between the inter - layer j - current elements .
this feature in combination with the low density of such elements justifies the `` ideal gas '' approximation for them , which in its turn leads to finite values of the renormalized inter - layer josephson coupling .
the question may be raised if a presence of long - range forces between the inter - plane j - currents @xmath134 can change the situation and lead to the sp or its weaker version  where @xmath289 with the growing number of layers @xmath125 in the limit @xmath287 . in this respect we note that in order to realize this , fluctuations of the difference of the j - currents with positive and negative orientations must be macroscopically suppressed . in this case
the fluctuation of the winding numbers in @xmath131-direction @xmath290 will scale slower than @xmath149 so that @xmath291 .
this may be caused by interactions between the inter - layer j - currents decaying not faster than the second power of their separation .
more specifically , the following additional repulsive term h_sp=_i , j u(_i -_j ) j_i j_j [ spx ] in the simple xy j - current model ( [ hxy2 ] ) with @xmath292 having the long range tail @xmath293 with @xmath294 will generate the energy contribution @xmath295 in terms of the windings in z - direction . consequently , the renormalized josephson coupling ( [ wu ] ) would scale as @xmath296 .
in one particular example long - range forces are introduced into the inter - layer josephson in the standard xy model ( [ hxy ] ) by some effective gauge - type term @xmath297 , where @xmath298 refers to the derivatives along the layers of some soft mode @xmath299 , with @xmath300 being a parameter .
the resulting interaction in the dual form ( [ spx ] ) becomes @xmath301 and , thus , it suppresses the inter - layer josephson as @xmath302 in the limit @xmath28 for fixed @xmath125 .
one of us ( abk ) acknowledges helpful discussions of the sliding phases with charles kane and boris svistunov .
he thanks kitpc ( the program `` precision many - body physics of strongly correlated quantum matter '' , may - june 2014 ) and aspen center for physics ( the program `` beyond quasiparticles : new paradigms for quantum fluids '' , august 2015 ) for hospitality when this work was in progress .
this work was supported by the national science foundation under the grant phy1314469 , and by a grant for computer time from the cuny hpcc under nsf grants cns-0855217 and cns - 0958379 .
99 k. b. efetov , sov .
jetp * 49 * , 905 ( 1979 ) .
s. e. korshunov and a. i. larkin , phys .
b * 46 * , 6395 ( 1992 ) .
p. w. anderson , phys .
lett . * 64 * , 1839 ( 1990 ) .
x. g. wen , phys.rev .
* b 42 * , 6623 ( 1990 ) .
c. castellani , c. di castro , and w. metzner , phys .
lett . * 69 * , 1703(1992 ) .
m. fabrizio , and a. parola , phys.rev.lett . * 70*,226(1993 ) . c. s. ohern and t. c. lubensky , phys.rev.lett .
* 80 * , 4345(1998 ) . c. s. ohern , t. c. lubensky , and j. toner , phys.rev.lett . *
83 * , 2745 ( 1999 ) ; s. l. sondhi , kun yang , phys . rev . *
b63*,054430(2001 ) .
r. mukhopadhyay , c. l. kane , and t. c. lubensky , phys .
* b 63 * , 081103(r ) ( 2001 ) ; a. vishwanath and d. carpentier , phys .
* 86 * , 676 ( 2001 ) ; a. f. andreev and e. p. bashkin , sov .
jetp * 42 * , 164(1976 ) .
m. a. cazalilla , a. f. ho , and t. giamarchi , new j. phys .
* 8 * , 158 ( 2006 ) d. pekker , g. refael , and e. demler , phys .
105 * , 085302 ( 2010 ) ; t. c. lubensky , p. m. chaikin , principles of condensed matter physics , cambridge university press ( 2000 ) .
e. babaev , nucl.phys .
* b 686*,397 ( 2004 ) . v. m. kaurov , a. b. kuklov , and a. e. meyerovich phys .
lett . * 95 * , 090403 ( 2005 ) .
e. k. dahl , e. babaev , a. sudbo , phys .
rev . * b78 * , 144510 ( 2008 ) . g. parisi , in _ statistical field theory , frontiers in physics _
, addison - wesley , reading , ma , 1988 .
prokofev , b.v .
svistunov , and i.s .
tupitsyn , phys .
a * 238 * , 253 ( 1998 ) ; jetp * 87 * , 310 ( 1998 ) ; n.v .
prokofev and b.v .
svistunov , phys .
rev . lett .
* 87 * , 160601(2001 ) .
j. villain , j. phys .
( paris ) * 36 * , 581(1975 ) .
m. wallin , e.s .
sorensen , s.m .
girvin , and a. p. young , phys . rev .
* b 49 * , 12115 ( 1994 ) . w. janke and h. kleinert , nucl . phys .
* b270 * , 135(1986 ) .
e. i. pollock and d. m. ceperley , phys . rev . *
b 36 * , 8343 ( 1987 ) .
a.m. polyakov , _ gauge fields and strings _ , in contemporary concepts in physics , vol.3 , harwood academic publishers , chur - melbourne , ch.4 ( 1987 ) .
c.l.kane , m.p.a .
fisher , phys .
. lett . * 68*,1220 ( 1992 ) .


in the last few years , the universality of the stellar initial mass function ( imf ) has been questioned by the results of theoretical studies , which consistently predict that the first stars ( hereafter pop iii stars ) had characteristic masses of 100 - 600 @xmath3 ( omukai & nishi 1998 ; abel , bryan & norman 2002 ; bromm , coppi & larson 2002 ; omukai & palla 2003 ) , about 100 times more massive than those observed today ( larson 2003 ) .    indeed , the physical conditions in primordial star - forming regions appear to systematically favor the formation of very massive stars . in particular , ( i ) the fragmentation scale of metal - free clouds is typically @xmath4 ( abel , bryan & norman 2002 ; bromm , coppi & larson 2002 ) ( ii ) because of the absence of dust grains the radiative feedback from the forming star is not strong enough to halt further gas accretion ( omukai & palla 2003 ) .
( iii ) since the accretion rate is as large as @xmath5-@xmath6 , the star grows up to @xmath7 within its lifetime ( stahler , palla & salpeter 1986 ; omukai & nishi 1998 ; ripamonti et al .
2002 ) .    in spite of the significant progresses made by numerical simulations and semi - analytic models ,
many aspects of primordial star formation remain to be fully understood , the more important of which are the radiative feedback effects which operate during protostellar accretion and which are likely to ultimately set the final stellar mass . by modeling the structure of the accretion flow and the evolution of the protostar , tan & mckee ( 2004 )
have recently shown that radiative feedback becomes dynamically significant at protostellar masses @xmath8 , and are likely to constrain the mass of first stars in the range @xmath9 .
on the other hand , observations of present - day stellar populations show that stars form according to a salpeter imf with a characteristic mass of 1 @xmath3 , below which the imf flattens . thus ,
unless the current picture of primordial star formation is lacking in some fundamental physical process , a transition in the modes of star formation must have occurred during cosmic evolution .
the physical processes responsible for this transition , as well as its cosmological consequences , are currently the subject of systematic investigation . in particular , following the early study by yoshii & sabano ( 1980 ) , it has been shown that a key element driving this transition is the metallicity of the star - forming gas ( omukai 2000 ; bromm et al . 2001 ; schneider et al . 2002 , 2003 ) .
the fragmentation properties of the collapsing clouds change as the mean metallicity of the gas increases above a critical threshold , @xmath10 ( schneider et al .
2002 , 2003 ) .
the fragmentation of clouds with @xmath11 proceeds only down to relatively massive ( @xmath12 ) cores , whereas in clouds with @xmath13 lower mass fragments can be formed . within the critical metallicity range
, low - mass cores can form if a sufficient amount of metals are depleted onto solid dust grains , which provide an efficient cooling channel in the high density regime ( schneider et al .
2003 ; omukai et al .
2005 ) .    on cosmic scales ,
the transition between an epoch dominated by the formation of massive pop iii stars to an epoch dominated by the formation of ordinary pop ii / i stars is controlled by the strength of the so - called _ chemical feedback _ : the explosions of the first massive supernovae starts to metal - enrich the gas out of which subsequent stellar generations form .
in particular , in regions of the universe where the metallicity exceeds @xmath14 , low - mass stars can form .
chemical feedback self - propagates through the metal enrichment process and therefore depends on many poorly constrained parameters , such as the star formation efficiency at high redshifts , the number of supernovae per stellar mass formed , the efficiency of metal ejection and mixing . in particular , the shape of the primordial imf governs the number of pop iii stars with masses in the pair - instability supernova range ( pisn , heger & woosley 2002 ) 140 @xmath3 - 260 @xmath3 ( which contribute to the metal enrichment ) relative to that of pop iii stars which end up as massive black holes ( which instead do not contribute to the metal enrichment as all the heavy elements are captured by the black hole ) .    using a semi - analytic approach , scannapieco , schneider & ferrara ( 2003 )
have modeled the probability that a collapsing halo has been impacted by metal - enriched outflows from neighboring halos .
exploring different chemical feedback strengths , they concluded that the transition between the pop iii and pop ii / i dominated star formation epochs is essentially controlled by the spatial distribution of metals . because metal enrichment is generally inefficient ( the filling factor of metals
is always found to be @xmath15 ) , the transition epoch is extended in time , with coeval pop iii and pop ii / i star formation episodes , and pop iii stars still contributing to the global star formation rate at redshifts @xmath16 . these findings have very important implications for the detection of pop iii stars through their strong ly@xmath17 emission ( scannapieco et al .
2003 ) as well as the corresponding pair - instability supernova rate ( scannapieco et al .
2005 ) .    in this work ,
we investigate the cosmic transition between pop iii and pop ii / i stars following a complementary approach .
we use the code pinocchio developed by monaco , theuns & taffoni ( 2002 ) and taffoni , monaco & theuns ( 2002 ) to produce random catalogues of dark matter halos at different redshifts .
we then propagate metal enrichment along the hierarchy of mergers of star - forming galaxies to predict the number and mass distribution of metal - free objects , which can host pop iii stars , at each redshift .
we then couple these findings to the star formation history ( of pop iii and pop ii / i stars ) and to the history of cosmic reionization and quantify their interplay through radiative feedback effects which can contribute to the extinction of pop iii stars .
a similar approach has been recently undertaken by furlanetto & loeb ( 2005 ) with the aim of assessing the plausibility of the double reionization scenario suggested by cen ( 2003 ) as a way to meet both the constraints on the gas neutral fraction implied by observations of gunn - peterson troughs in the spectra of distant quasars ( fan et al .
2002 , 2003 ) and on the gas ionized fraction derived by the high optical depth to thomson scattering measured by the wmap satellite ( spergel et al . 2003 ; kogut et al .
2003 ) .    here
we instead focus on the extinction epoch of pop iii stars which we believe to be regulated by feedback both of radiative and chemical types .
we investigate the implications that complementary observations on galaxy number counts at redshifts @xmath1 ( bouwens et al . 2005 ) and on the history of cosmic reionization ( kogut et al . 2003 ) set on the transition epoch between pop iii and pop ii / i stars and on the nature of pop iii stars as sources of uv photons contributing to cosmic reionization .
in particular , the detection of galaxies at the highest redshifts is challenging because of the redshifting of uv light into the infrared and the well - known limitations of current ir instruments . in spite of these difficulties ,
the search for high redshift galaxies , based on the dropout technique , has provided the first tentative detection of galaxies at redshifts beyond @xmath18 ( bouwens et al .
2004 ) , stimulating a re - analysis of the small amount of deep ir data .
bouwens et al .
( 2005 ) looked at the prevalence of galaxies at @xmath19 by applying the dropout technique to the wide variety of deep fields that have been imaged with the near infrared camera and multi - object spectrometer ( nicmos ) on board of the hubble space telescope .
using an appropriate selection criterion for high-@xmath20 sources , bouwens et al .
( 2005 ) concluded that the number of @xmath21 sources must be three or fewer . when combined with the high optical depth to thomson scattering inferred from the _ wilkinson microwave anisotropy probe _
( wmap ) satellite ( kogut et al 2003 ) , these data already provide a demanding benchmark for theoretical models .
the paper is organized as follows . in section  2
we describe the different feedback effects that control the transition from a pop  iii to a pop  ii dominated star formation epoch .
we illustrate the adopted chemical feedback model , predicting the number and mass function of halos which can host pop iii stars at different redshifts . in section  3 , we couple these predictions to a star formation and reionization model .
in section  4 we summarize the observations of galaxy counts at redshifts @xmath22 and discuss how these can be used to set important constraints on the model . finally , in section  5 we present the main results and in section  6 we discuss their implications .    throughout the paper ,
we adopt a @xmath23cdm cosmological model with parameters @xmath24 , @xmath25 , @xmath26 , @xmath27 , @xmath28 and @xmath29 , which lie within the experimental errorbars of wmap experiment ( spergel et al .
we also use the ab magnitude system , where @xmath30 is the spectral energy density within a given passband in units of erg s@xmath31 @xmath32 hz@xmath31 ] .
at least three types of feedback effects can be relevant to the extinction of the first massive stars : ( i ) radiative feedback caused by h@xmath34 photo - dissociating radiation emitted by the first stellar sources , which can inhibit gas cooling and star formation inside small halos ( ii ) chemical feedback due to the release of metals and dust in the first pop iii supernova explosions , which pollute the gas out of which subsequent generations of pop ii stars form ( iii ) radiative feedback due to photo - heating filtering caused by cosmic reionization , the increased temperature of the cosmic gas suppresses the formation of galaxies below the jeans mass .
a thorough discussion of each of these processes is beyond the scope of this paper and can be found in the recent extensive review by ciardi & ferrara ( 2005 ) .    as a first step ,
we build up the hierarchy of mergers of dark halos hosting the star formation sites .
we use pinocchio code ( monaco , theuns & taffoni 2002 ; taffoni , monaco & theuns 2002 ) to generate random catalogues of dark matter halos in a hierarchical universe .
our reference run has @xmath35 particles within a volume of comoving side 6 mpc  @xmath36 so that the smallest resolved dark matter halos have total mass @xmath37 ( in the merger trees only halos with at least 10 particles are considered ) .
1 summarizes the details of our reference pinocchio run and shows the redshift evolution of the largest and smallest mass present in the simulation box and of the minimum halo mass that is allowed to form stars .
the latter corresponds , at each redshift , to halos with virial temperatures @xmath38k , m_min = 10 ^ 8 ( ) ^-3/2 [ eq : min ] so that stars do not form in halos which rely only on h@xmath34 cooling ( radiative feedback ) .
the chosen resolution ensures that we can follow the formation of stars in the lowest mass halos at the highest redshift .
due to radiative feedback , the first star - forming halos have a total mass @xmath39 and start to form at @xmath40 .
this resolution is appropriate to follow the transition from a pop  iii to a pop  ii  star formation mode around redshift @xmath1 , which is the main scope of this analysis .    in fig .
2 we compare the resulting mass function of existing halos at redshifts 30 , 20 , 10 and 5 in the simulation box , to the analytic predictions of press & schechter ( 1974 ) and sheth & tormen ( 1999 ) .
the errorbars are poisson errors on the counts in each mass - bin .
the agreement is very good down to redshift @xmath16 showing that the chosen resolution ensures completeness in the simulation volume in the redshift range relevant to our analysis .      at each redshift , existing halos which are allowed to form stars ( @xmath41  k ) are classified as hosting pop ii ( pop iii ) stars depending on whether the halo itself or any of its progenitors have ( have not ) already experienced an episode of star formation , leading to the formation of sn progenitors which then metal - enrich the gas .
note that even a single 200 @xmath3 pisn can enrich the gas within a halo with mass @xmath42 to a metallicity level well above the critical value as ,    z = 0.02 z _ ( ) ^3/2    where , following heger & woosley ( 2002 ) , we have assumed that the mass of metals released in a pisn explosion is 45% of the initial progenitor mass .
this way of taking into account chemical feedback is very conservative as ( i ) we are assuming that all metals synthesized during a star formation episode remain confined ( or are re - accreted ) within the host halo and/or ( ii ) even if a fraction of metals is ejected out of the host halo , the metals that are kept within the potential well are sufficient to pollute the gas within the host system or within any unpolluted halo that the system will later merge with , to a metallicity level larger than @xmath14 . on the other hand , we are neglecting metal pollution from neighboring galaxies through enriched - winds ( scannapieco et al .
2003 ; sigward , ferrara & scannapieco 2005 ) .
3 shows the fraction of halos which can host pop iii stars ( hereafter pop iii halos ) at different redshifts as a function of their mass .
the three different curves correspond to different assumptions on the shape of the primordial imf .
in particular , we randomly classify pop  iii halos as sn - forming or bh - forming depending on whether they host stars with masses in the pisn range or outside this range .
sn - forming halos are metal - enriched and pollute their descendants .
bh - forming halos are not enriched in metals and therefore do not propagate chemical feedback along their subsequent hierarchy of mergers .
the fraction of sn - forming halos , @xmath43 , is therefore a free parameter that controls the strength of chemical feedback .
the three curves shown in each panel of fig .
3 correspond ( from bottom to top ) to @xmath44 .
each curve represents an average over 50 realizations of the random classification procedure and the errors are poisson errors on the counts in each mass bin .
when @xmath45 , all pop iii halos are sn - forming and chemical feedback is strong . the total fraction of pop iii halos decreases with time , being @xmath46 at @xmath47 and @xmath48 at @xmath49 ; at all redshifts , pop iii stars are confined to form in the smallest halo mass bins ( close to the minimum mass to form stars shown as the vertical dotted line in each panel ) . for @xmath50 ( 50% of pop iii halos host sn ) , the number of pop iii halos at each redshift is larger than the previous case ( @xmath51 at @xmath47 and @xmath52 at @xmath53 ) but pop iii halos are still confined to the smallest mass bins ( @xmath54 ) . finally , when @xmath55 the strength of chemical feedback is significantly reduced ( as only 10% of pop iii halos at each redshift host sn ) and therefore the fraction of pop iii halos is dominant at all redshifts , being @xmath56 at @xmath47 and @xmath57 at @xmath49 , not negligible even in halos with the largest mass present in the simulation box at each redshift .
hereafter we only consider @xmath55 and @xmath58 as representative of weak and strong chemical feedback models .    despite the strength of chemical feedback ,
pop iii halos continue to form stars even at redshifts as small as @xmath16 , in agreement with the results of scannapieco et al .
thus , independently of whether it preferentially propagates through galaxy mergers or through metal - enriched winds , chemical feedback always leads to a smooth transition between a pop iii and a pop ii dominated star formation epoch .
note that the above results depend only on the parameter @xmath43 and on the hierarchy of galaxy mergers , and are totally independent of the star formation efficiency or of the value of critical metallicity for the transition in fragmentation scales of collapsing pre - stellar clouds .    [ cols="^,^,^,^,^,^,^,^,^,^,^,^,^,^,^,^,^,^",options="header " , ]      so far , we have taken into account radiative feedback only through the choice of the minimum halo mass that can host stars .
namely , we have assumed that mini - halos with virial temperatures less than @xmath59k do not significantly contribute to the cosmic star formation history .
however , radiative feedback does affect larger halos , particularly when the universe is approaching the epoch of cosmic reionization .
indeed , the increase in temperature of cosmic gas during reionization will suppress the formation of small galaxies with mass below the jeans mass . to quantify this effect ,
we follow the method proposed by gnedin ( 2000 ) and adopt the filtering mass as the characteristic scale below which the gas fraction in star forming objects , hence the fraction of stellar mass , is reduced with respect to the universal value . in particular , before reionization , the average baryon fraction ( the ratio of the mass in baryons over the total halo mass ) has the universal value @xmath60 while after reionization it can be approximated by the fitting formula ( gnedin 2000 ) ,    f_b = , [ eq : fil ]    where the characteristic mass , @xmath61 , is the total mass of objects which , on average , retain 50% of their universal baryon mass , @xmath62 .
this characteristic mass is very well approximated by the linear - theory filtering mass , @xmath63 , with m_f^2/3 = _
0^a da^m_j^2/3(t , a^)[1-()^1/2 ] ,    where the jeans mass is given by    m_j(t ) = 2.56 10 ^ 5 ( _ m h^2)^-1/2(t/)^3/2 ( 1+z)^-3/2 ,    with @xmath64 and @xmath65 the temperature and molecular weight of the cosmic gas .
here we assume the temperature of ionized regions to be @xmath66  k and we modulate the increase in temperature of the cosmic gas due to reionization with the volume filling factor of ionized regions at the corresponding redshift @xmath67 ( see section  [ sec : reio ] ) .
thus , the characteristic mass depends on the thermal history of the gas and it accounts for the finite time required for pressure to influence the gas distribution .
the redshift evolution of the filtering mass is shown in fig .  1 for a model where reionization occurs at redshift @xmath68 ( short - dashed line ) and @xmath69 ( long - dashed line ) .
the filtering mass steeply increases with time and , as expected , its value in the redshift range @xmath70 depends strongly on the reionization history . at redshift @xmath1 ,
the filtering mass is @xmath71 for the late reionization case , and it is a factor ten larger for the early reionization case .    as can be inferred from fig .  1 and  3 , following reionization
, photoionization feedback can significantly suppress the formation of pop iii stars since halos with masses @xmath72 have evaporated at least 50% of their mass .
following barkana & loeb ( 2001 ) , we compute the redshift evolution of the filling factor of hii regions as , q_hii(z ) = _
z^ dz^ ^f(z^,z ) , [ eq : reio ]    where @xmath73 and @xmath74 are the present - day number densities of hydrogen and baryons ( @xmath75 is the hydrogen mass fraction ) , and @xmath76 is the production rate of ionizing photons that we can express as ,    = f_esc n _ _ ,    where @xmath77 is the escape fraction , the fraction of ionizing photons that is able to escape the host galaxy , @xmath78 is the time - averaged number of ionizing photons emitted per unit stellar mass formed , and @xmath79 is the comoving star formation density , _
= f_f_b _ m_min^ dm m ( m , z ) where @xmath80 is the star formation efficiency , @xmath81 is the baryon fraction given by eq .
( [ eq : fil ] ) , @xmath42 is the minimum halo mass to form stars ( see eq .
[ eq : min ] ) and @xmath82 is the number density of halos ( pop ii and pop iii ) at redshift @xmath20 .
the parameter @xmath77 also affects the relative strength of the stellar continuum and the ly@xmath17 line and nebular luminosity of each galaxy .
the latter contributions are particularly important for massive pop  iii stars .
finally , the function @xmath83 takes into account the effect of recombinations . assuming a time - independent volume - averaged clumping factor @xmath84 , common to all hii regions ,
we can write ,    f(z^,z ) = - c [ f(z^)-f(z ) ] ,    and    f(z)= ,    where @xmath85 is the hydrogen recombination rate .    for a given reionization history , it is straightforward to compute the optical depth to electron scattering as ,    _ e(z)=_0^z dz^c _
t n_e(z^ ) ,    where @xmath86 is the thomson cross section and @xmath87 is the number density of free electrons at redshift @xmath20 . in this paper
we only consider the evolution of the filling factor of hii regions ( @xmath88 ) .
however , for the purpose of computing the optical depth to electron scattering , we assume that @xmath89 and we neglect the contribution of electrons coming from heiii regions ( see also haiman & holder 2003 )
. under these assumptions we can write @xmath90 where @xmath67 is given by eq.(6 ) .
bouwens et al . ( 2005 ) looked at the prevalence of galaxies at @xmath19 by applying the dropout technique to the wide variety of deep f110w- and f160w - band ( hereafter @xmath91 and @xmath92 , respectively ) fields that have been imaged with the near infrared camera and multi - object spectrometer ( nicmos ) .
the principal data set is constituted by two 1.3 arcmin@xmath93 deep nicmos parallels taken with the advanced camera for surveys ( acs ) hubble ultra deep field ( udf ; each has about @xmath94 orbits data ) .
complementary shallower fields possessing similar @xmath95 imaging have been analyzed as well .
the 5@xmath96 limiting magnitude for the udfs is @xmath97 in @xmath91 and @xmath98 in @xmath92 ( @xmath99 aperture ) .
since these limiting magnitudes are reached only in the deepest regions of the two nicmos udfs , we adopt the conservative limiting magnitude @xmath100 ( bouwens , private communication ) .    the primary selection criterion for high-@xmath20 sources
is @xmath101 .
using this criterion , bouwens et al .
( 2005 ) found eleven sources .
eight of these are ruled out as credible @xmath21 sources , either as a result of detection ( @xmath102 ) blue - ward of @xmath91 or because of their colors red - ward of the break ( @xmath103 ) .
the nature of the three remaining sources could not be assessed from the data , but this number appears consistent with the expected contamination from low - redshift interlopers .
hence , bouwens et al .
( 2005 ) concluded that the actual number of @xmath21 sources in the nicmos parallel fields must be three or fewer .    adopting the same selection criterion as in bouwens et al .
( 2005 ) , we compute the number of @xmath104 galaxies detectable in the nicmos udfs predicted by our models . assuming that the three candidates selected by bouwens et al .
( 2005 ) are indeed @xmath105 galaxies , we can set interesting constraints on the models . in particular , for a given chemical feedback strength @xmath43 ,
we compute the number density of pop iii and pop ii halos at each redshift .
we then assume a universal value for the star formation efficiency @xmath106 and escape fraction @xmath77 of ionizing photons and we characterize each stellar population ( pop iii and pop ii / i stars ) by assuming a specific template emission spectrum .    in the following , we consider two models , a and b , whose details are given in table  1 .
we assume two representative chemical feedback strengths ( @xmath43 = 0.1 for model a and @xmath43 = 1 for model b ) . in both cases ,
the emission properties of pop  ii stars have been computed using a simple stellar population model taken from the galaxev library ( bruzual & charlot 2003 ) .
pop  ii stars have an initial metallicity of @xmath107 and are assumed to form according to a salpeter imf ( with masses in the range 0.1 @xmath3 - 100 @xmath3 ) .
pop  iii stars are assumed to be very massive ( with masses larger than 100 @xmath108 ) and their emission properties are computed using the stellar spectra of schaerer ( 2002 ) , including the nebular continuum emission found to be very important for stars with strong ionizing fluxes ( bromm , kudritzki & loeb 2001 ; schaerer 2002 ) , we are not forced to make any assumption about the shape of the imf of pop  iii stars as long as these are assumed to be very massive . ] .    for each model
, we then compute the number of sources per unit solid angle with observed flux in the range @xmath109 and @xmath110 as    @xmath111    where @xmath112 is the comoving volume element per unit solid angle and redshift , @xmath113 is the comoving number of objects at redshift @xmath20 with observed flux in the range @xmath114 $ ] , given by    @xmath115    here @xmath116 is the halo formation rate , and the factor @xmath117 converts the number density per unit mass in unit flux interval .
the average flux @xmath109 from a halo of mass @xmath118 at redshift @xmath20 that has formed at @xmath119 is ,    @xmath120    where @xmath121 is the template luminosity per unit solar mass for a population of age @xmath122 ( the time elapsed between the redshift @xmath123 and @xmath20 ) , @xmath124 is the luminosity distance , @xmath125 is the instrumental bandwidth , @xmath126 and @xmath127 are the restframe frequencies corresponding to the observed ones . finally , @xmath128 is the intergalactic medium ( igm ) effective optical depth at @xmath129 between redshift @xmath130 and @xmath20 ( see sect .
2.2 of salvaterra & ferrara 2003 )
we compute the expected number of @xmath1 galaxies detectable in the nicmos udfs using eqs .
[ eq : counts1]-[eq : counts3 ] and applying the same selection criterion of bouwens et al .
( 2005 ) ( @xmath101 , integrated up to the magnitude limit of the survey , @xmath100 ) .
assuming that only the three most luminous sources should be above the detection threshold , we derive the maximum allowed value for the star formation efficiency , @xmath80 , corresponding to a given escape fraction of ionizing photons .
the results obtained for models a and b are shown in fig .  4 in terms of comoving star formation density in pop  iii and pop  ii stars as a function of redshift .
model c will be discussed in section [ sec : reio ] .    in each panel , vertically ( horizontally ) shaded regions represent the range of star formation histories for pop  iii ( pop  ii ) stars allowed by the nicmos - udfs counts when a value @xmath131 in assumed ( see table  1 for the corresponding values of @xmath80 ) . for both models and all the explored values of @xmath77 ,
the three most luminous objects in nicmos udfs are always represented by galaxies hosting pop  iii stars , as their strong ly@xmath17 and nebular emission ( schaerer 2002 ) is dominant . as a consequence , @xmath80
is restricted to relatively small values , which increase with increasing @xmath77 .
model b , being characterized by the strongest chemical feedback , and therefore by the smallest fraction of halos hosting pop  iii stars at redshift @xmath132 ( see also fig .
3 ) , has larger allowed star formation efficiencies than model a.    fig .
4 shows that in model a chemical feedback is very inefficient ( @xmath133 ) and pop  iii stars dominate the cosmic star formation history down to redshift @xmath134 .
conversely , model b shows that when chemical feedback is maximally efficient , the star formation density in pop  iii stars is sub - dominant at all but the highest redshifts , but pop  iii stars still make a significant contribution to the cosmic star formation history down to redshift 5 .    because of the moderate star formation efficiencies implied by the observed nicmos udfs counts , radiative feedback is negligible for the most massive galaxies . to evaluate its impact on the predicted star formation densities shown in fig .  4 , we have assumed a clumping factor @xmath135 . for @xmath136 ,
the corresponding reionization redshifts vary in the range @xmath137 for model a , and in the range @xmath138 for model b.    in model a , radiative feedback effect on the overall pop  iii star formation density manifests itself as a progressive decrement of the upper limit ( which corresponds to the largest @xmath77 and @xmath80 ) of the shaded region for @xmath139 . in model b , where the star formation efficiencies implied by the observed counts are larger and , at each redshift , pop  iii stars are confined to form only in the smallest halos , radiative feedback is more significant ; around redshift 10 , the lines corresponding to the maximum and minimum pop  iii star formation efficiencies cross each other ( the highest the star formation efficiency and escape fraction , the strongest the radiative feedback ) .
overall , the observed galaxy counts at redshift @xmath1 favor a cosmic star formation history where the transition from pop  iii to pop  ii stars is governed by strong chemical feedback , with all pop  iii halos efficiently propagating metal enrichment to their descendants along the hierarchy of galaxy mergers .
in fact , in less efficient chemical feedback scenarios ( @xmath140 ) , the predicted cosmic star formation density is dominated by pop iii stars down to redshift @xmath53 .
stronger constraints on models a and b come from their corresponding reionization histories
. figs .  5 and  6 show isocontours of redshifts of reionization as a function of the clumping factor and the escape fraction ( the star formation efficiency is set by the observed counts at @xmath1 ) .
we define the redshift of reionization , @xmath141 , as the redshift when the filling factor of ionized regions , given by eq .
( [ eq : reio ] ) is @xmath142 .
the shaded regions indicate the parameter space which leads to an optical depth to electron scattering within the range @xmath143 observed by the wmap satellite ( spergel et al .
2003 ; kogut et al .
2003 ) , with the dashed line corresponding to the central value .
as it can be inferred from the figures , both models lead to reionization histories which are in poor agreement with observations .
indeed , large escape fractions and moderate clumping factors are required to reproduce the observed range of optical depths .
independent observations of the escape fraction give , with few exceptions , values in the range @xmath144 ( ciardi & ferrara 2005 and references therein ) whereas numerical simulations ( gnedin & ostriker 1997 ; iliev , scannapieco & shapiro 2005 ) and semi - analytic models ( madau , haardt & rees 1999 ) predict clumping factors in the range @xmath145 .
5 shows that in model a , even assuming an homogeneous igm ( @xmath146 ) , the escape fraction has to be @xmath147 to be marginally consistent with the observations . in model b ,
the agreement improves but still requires extreme assumptions on either the escape fraction or the clumping factor : when @xmath144 , the predicted optical depth lies within -1@xmath96 error from the observed value only if @xmath148 .
our analysis suggests that it is very hard to reconcile the observed number counts at redshift @xmath1 with the reionization history inferred by wmap observations if pop  iii stars are assumed to be all very massive , with masses @xmath149 .
this result is somehow counter - intuitive because metal - free very massive stars have always believed to play an important role in an early epoch of cosmic reionization ( cen 2003 ; haiman & holder 2003 ; sokasian et al .
2004 , 2004 ; venkatesan , tumlinson & shull 2003 ; wyithe & loeb 2003 ; but see also ciardi , ferrara & white 2003 , and ricotti & ostriker 2004 for alternative interpretations ) .
still , due to their strong ly@xmath17 and nebular luminosity , pop  iii stars with masses @xmath150 always dominate the predicted galaxy number counts at @xmath1 , even assuming the strongest chemical feedback model .
the paucity of candidate objects observed in nicmos udfs implies that star formation in halos with @xmath151  k must have occurred very inefficiently . even if all three objects identified by bouwens et al .
( 2005 ) are truly @xmath1 sources , the contribution of pop  iii stars to cosmic reionization is strongly constrained .
the constraints are particularly severe when chemical feedback is less efficient ( model a ) , because pop  iii stars can form at redshift @xmath1 in relatively large , and therefore luminous , galaxies .
even if chemical feedback is assumed to be maximally efficient and pop  iii stars are confined to form in the smallest halos , their contribution to cosmic reionization is in agreement with observations only for moderate ( @xmath152 ) clumping factors .
unless an additional type of feedback comes into play , the transition between pop  iii to pop  ii halos , regulated by chemical and radiative feedback , is very smooth .
thus , to relax the severe upper limits placed on the star formation efficiencies by nicmos udfs counts , the emission properties of pop  iii stars have to be modified .    in particular ,
if pop  iii stars are assumed to form according to a salpeter imf with masses in the range @xmath153 , their ly@xmath17 and nebular emission is strongly suppressed . in this limit , ( see model c in table  1 ) , all pop  iii halos contribute to chemical feedback , @xmath154 , because each burst of star formation leads to the formation of sn progenitors and nicmos udfs counts at @xmath1 are dominated by pop  ii galaxies . as a consequence ,
the upper limit on the star formation efficiency set by the observed counts is @xmath155 , much higher than those found for models a and b and independent of the escape fraction .
the predicted evolution of the star formation density in pop ii and pop iii stars is shown in the lowest panel of fig .
similar to model b , pop ii star dominate the star formation history at all but the highest redshifts .
indeed , in our approach the transition between pop iii and pop ii dominated star formation epochs is controlled only by the strength of chemical feedback ( @xmath156 ) and it is independent of the assumed star formation efficiencies ( see section [ sec : chem ] ) . however , contrary to model b , the star formation density depends on the escape fraction only indirectly , through the strength of radiative feedback . for this reason
, the shaded regions representing the range of allowed star formation densities are not degenerate only after reionization which , for @xmath157 and @xmath158 occurs at redshift @xmath159 ( for pop ii stars , which are preferentially hosted in large halos , this effect is negligible ) .
the implications for the history of reionization are summarized in fig .
7 , which shows that model c is in much better agreement with observations than model a and b. the optical depth observed by wmap can be reproduced with @xmath144 even if the clumping factor is as high as @xmath160 .
thus , when the observed @xmath1 counts are used to set an upper limit to the allowed star formation efficiency , pop  iii stars forming with a conventional imf are more efficient sources of cosmic reionization than very massive stars . note that this conclusion is independent on issues related to chemical and radiative feedback .
indeed , both models b and c assume maximally efficient chemical feedback and radiative feedback plays a minor role for model c because the predicted @xmath1 counts are dominated by massive pop  ii galaxies , which are largely unaffected by photoevaporation .
fig.7 shows that the best - fit value of optical depth implied by wmap observations , @xmath161 , can only be reproduced by model c with @xmath144 if @xmath162 .
if larger values of @xmath84 are to be favored or the observed optical depth will be corrected up - ward in future wmap data releases , a revision of model c is to be considered .    however , unless the nicmos udfs counts have under - estimated by a large factor the number of @xmath1 sources , our results suggest that the only plausible revision of model c is to modify the primordial imf by increasing the characteristic mass but allowing pop  iii stars to form only in the conventional mass range , @xmath163 . in fig .
8 we show the predicted optical depth as a function of the pop  iii characteristic stellar mass , @xmath164 .
we assume that the characteristic stellar mass enters in the definition of the stellar imf as suggested by larson ( 1998 ) ,    ( 1 + ) ^-1.35 , which , for @xmath165 leads to the standard salpeter imf .    for a given @xmath164 , we compute @xmath166 considering the proper emission rate of ionizing photons ( taken from table  1 of ciardi & ferrara 2005 ) , assuming @xmath167 , three different clumping factors @xmath168 and 30 , and that either all 3 or none of the candidates identified by bouwens et al .
( 2005 ) in nicmos udfs represent @xmath1 galaxies ( upper and lower limits to each shaded region ) .
8 shows that the predicted optical depth depends more on the clumping factor than on the assumed value of @xmath164 . if @xmath135 , the predicted @xmath166 is in agreement with the observed values even assuming @xmath165 , consistent with the results shown in fig .
7 . if @xmath169 , wmap and nicmos udfs data can be reconciled only if all three candidates are confirmed to be @xmath1 galaxies and if the characteristic pop  iii stellar mass is @xmath170 .
finally , if the observed @xmath166 will be corrected up - ward in future analyses of wmap data , none of the models considered , independently of the assumed @xmath164 , will be able to reconcile the observed reionization history with the counts at @xmath171 , unless @xmath172 .
the tentative detections of @xmath1 galaxies in nicmos udfs by bouwens et al .
( 2005 ) provides tight constraints on the relative contributions of pop  iii and pop  ii stars to the galaxy luminosity function .
when coupled to the reionization history of the universe implied by the wmap satellite , this collection of data sets represents already a very demanding benchmark for galaxy evolution models .
stimulated by these data , we have followed the assembling history of high - redshift galaxies hosting the first stars in a hierarchical model for structure formation . adopting simple prescriptions for chemical and radiative feedbacks along with available data
, we have constrained the evolution of the imf of the first stars .
in particular we find that :    1 .   two types of feedbacks need to be considered : ( ia ) radiative feedback , which suppresses star formation in h@xmath34-cooling halos and the formation of low - mass galaxies because of photoionization heating following reionization , and ( ib ) chemical feedback due to the progressive enrichment of ( pop  ii ) star forming gas by the metals released by the first ( pop  iii ) stellar explosions .
2 .   if chemical feedback is weak ( @xmath133 ) pop  iii stars can form in relatively large halos and dominate the evolution of the cosmic star formation density down to redshift @xmath16 .
more realistic scenarios are found assuming @xmath154 , where pop  iii stars are confined to form only in the smallest halos at each redshift .
3 .   the combined effect of radiative and chemical feedbacks fails to suppress the formation of massive pop  iii stars at @xmath139 .
if pop  iii stars are assumed to form with a top - heavy imf , that is with masses @xmath173 , pop  iii galaxy are always found to dominate the predicted @xmath1 nicmos udfs counts . because of their strong ly@xmath17 and nebular luminosities , very tight upper limits on the star formation efficiency in @xmath174  k halos are placed by the paucity of observed counts .
the corresponding reionization histories fail to reproduce the observed wmap optical depths unless uncomfortably large values for the escape fraction of ionizing photons and/or too small values for the igm clumping factor are assumed .
reionization constraints from wmap and observed @xmath1 counts in nicmos udfs can be reconciled if pop  iii stars form according to a larson imf with stellar masses in the conventional range @xmath175 and a characteristic mass @xmath176 if the clumping factor is @xmath177 , or @xmath178 if the clumping factor is as large as 30 . in the latter case ,
at least 1 candidate identified by bouwens et al .
( 2005 ) has to be a galaxy at @xmath1 .
our analysis is based on the assumption that chemical feedback propagates along the hierarchy of galaxy mergers from progenitors to descendants rather than through metal - enriched galaxy outflows ( scannapieco et al . 2003 ) . by adopting this simplified scheme
, we are implicitly assuming that ( i ) even if a fraction of metals is ejected out of the host halo , the metals that are kept within the potential well are sufficient to pollute the gas within the host system or within any unpolluted halo that the system will later merge with , to a metallicity level larger than @xmath14 and that ( ii ) metal pollution from neighboring galaxies is negligible . in spite of these simplifying assumptions ,
the resulting pop  iii and pop  ii star formation histories are consistent with the predictions of scannapieco et al .
in particular , even assuming a maximally efficient chemical feedback , the transition is always predicted to be very smooth , with pop  iii stars which continue to form in the smallest halos at each redshift down to @xmath16 .
interestingly enough , our conclusion is independent of the primordial imf and/or star formation efficiency .
we show that models where all pop  iii halos are able to metal pollute their descendants ( @xmath154 ) are favored by observations .
the smoothness of the cosmic imf transition is therefore strictly linked to the inhomogeneous character of the metal enrichment process , independently of how it is implemented in the models .    even if we have adopted a very approximate description of the thermal evolution of the cosmic gas , we find that radiative feedback following reionization can play an important role in the extinction of pop  iii stars . in particular ,
if reionization occurs at redshifts @xmath179 ( @xmath1 ) the baryon fraction in @xmath180 halos is reduced to 70% ( 99% ) of its universal value ( @xmath60 ) at redshift 10 and to 10% ( 30% ) of this value at redshift 5 .
thus , the baryon content in the smallest halos where pop  iii stars can form at low redshift is severely limited by radiative feedback .
this is strongly related to the epoch of cosmic reionization which , in turn , depends on the degree of inhomogeneities which characterize the cosmic gas .
we have parametrized this quantity through a constant clumping factor @xmath84 and we have explored different values of @xmath84 ranging from 1 to 50 .
however , it is likely that the degree of clumpiness of the cosmic gas varies with time , rising from @xmath181 at high redshift to a few tens at @xmath1 as structures go increasingly non - linear .
recently , iliev et al . (
2005 ) have studied the effect of igm clumping at redshifts @xmath182 , providing a fit to the evolution of the clumping factor ( see their fig .  7 ) where @xmath183 at redshift 30 and raises to @xmath184 at redshift 10 .
adopting this evolution for the clumping factor and @xmath185 we find that models a and b are unable to reproduce the optical depths inferred from wmap data whereas in model c we find that reionization is complete at redshift 13.5 and that the resulting optical depth is @xmath186 in very good agreement with the observed value .
thus , our analysis suggests that pop  iii stars forming with masses in the same dynamical range as pop  ii / i stars , but possibly with higher characteristic stellar mass , appear to be favored by observations .
several scenarios that enable the formation of low - mass stars in metal - free environments have been proposed , among which the bi - modal imf ( nakamura & umemura 1999 , 2002 ) , enhanced hd - cooling in relic hii regions ( uehara & inutsuka 2000 ; nakamura & umemura 2002 ) , the influence of a strong uv radiation field due to the vicinity of a very massive star ( omukai & yoshii 2003 ) , fragmentation in shock - compressed shells induced by the first very massive sn explosions ( mackey , bromm & hernquist 2003 ; salvaterra , ferrara & schneider 2004 ) .
however , all these processes either require very special conditions , generally related to the existence of a previous generation of very massive stars , or have failed to be realized in realistic numerical simulations . in spite of this , observations seem to indicate that either radiative feedback effects during the accretion phase on protostellar cores and/or a combination of the above mentioned processes can lead to the formation of pop  iii stars with masses @xmath187 . according to the results of our analysis
, the latter contribution to the pop  iii stellar imf must be dominant at least at redshifts @xmath1 .
if this scenario is correct , relics of the low - mass pop  iii star formation mode should be still present in our galaxy , eventually contributing to the lowermost end of the metallicity distribution function ( salvadori , schneider & ferrara in prep ) .
we acknowledge the use of pinocchio v 1.0 package written by pierluigi monaco , tom theuns & giuliano taffoni ( http://www.daut.univ.trieste.it/pinocchio ) . this work has been partially supported by the research and training network `` the physics of the intergalactic medium '' set up by the european community under the contract hprn - ct-2000 - 00126 .
99 abel t. , bryan g. l. , norman m. l. 2002 , science , 295 , 93 barkana , r. & loeb , a. 2001 , phr , 349 , 125 bouwens r. j. et al .
2004 , , apj , 616 , l79 bouwens r. j. , illingworth g. d. , thomson r. i. , & franx m , 2005 , apj , 624 , l5 bromm , v. , ferrara a. , coppi p. s. , larson r. b. 2001 , mnras , 328 , 969 bromm , v. , kudritzki , r. p. & loeb , a. 2001 , apj , 552 , 464 bromm , v. , coppi p. s. , larson r. b. 2002 , apj , 564 , 23 bruzual , g. & charlot , s. 2003 , 344 , 1000 ciardi , b. & ferrara , a. 2005 , space science reviews , 116 , 625 ciardi , b. , ferrara , a. , white , s. d. m. 2003 , mnras , 344 , l7 cen , r. 2003 , apj , 591 , 12 furlanetto , s. & loeb , a. 2005 , apj in press , ( astro - ph/0409736 ) gnedin , n. , ostriker , j. 1997 , apj , 474 , 581 gnedin , n. 2000 , apj , 542 , 535 haiman , z. & holder , g. p. 2003 , apj , 595 , 1 heger , a. & woosley , s. 2002 , apj , 567 , 532 iliev , i. , scannapieco , e. , shapiro , p. 2005 , apj , 624 , 491 kogut , a. et al .
2003 , apjs , 148 , 161 larson , r. b. 1998 , mnras , 301 , 569 mackey , j. , bromm , v. , hernquist , l. 2003 , apj , 586 , 1 madau , p. , haardt , f. , rees , m.  j. 1999 , apj , 514 , 648 monaco , p. , theuns , t. & taffoni , g. 2002 , mnras , 331 , 587 nakamura , f. , umemura , m. 1999 , apj , 515 , 239 nakamura , f. , umemura , m. 2002 , apj , 569 , 549 omukai , k. & nishi , r. 1998 , apj , 508 , 141 omukai , k. 2000 , apj , 534 , 809 omukai , k. & palla , f. 2003 , apj , 589 , 677 omukai , k. & yoshii , y. 2003 , apj , 599 , 746 omukai , k. , tsuribe , t. , schneider , r. , ferrara , a. 2005 , apj , 626 , 627 ricotti , m. & ostriker , j. p. , 2004
, mnras , 350 , 539 ripamonti , e. , haardt , f. , ferrara , a. , & colpi , m. , 2002 , mnras , 334 , 401 salvaterra , r. , ferrara a. , schneider r. 2004 , new astronomy , 10 , 11 salvaterra , r. & ferrara a. 2003 , mnras , 339 , 973 scannapieco , e. , schneider , r. & ferrara , a. 2003 , apj , 589 , 1 scannapieco , e. , madau , p. , woosley , s. , heger , a. , ferrara , a. 2005 , apj , in press ( astro - ph/0507182 ) schaerer , d. 2002 , a&a , 382 , 28 schneider r. , ferrara a. , natarajan p. , omukai k. 2002 , apj , 571 , 30 schneider r. , ferrara a. , salvaterra r. , omukai k. , bromm v. 2003 , nature , 422 , 869 sigward , f. , ferrara , a. & scannapieco , e. 2005 , mnras , 358 , 755 sokasian , a. , yoshida , n. , abel , t. , hernquist , l. , springel , v. 2004 , mnras 350 , 47 spergel , d. l. et al .
2003 , apjs , 148 , 175 stahler , s. w. , palla , f. , & salpeter , e. e. 1986 , apj , 302 , 590 tan , j. c. & mckee , c. f. 2004 , apj , 603 , 383 taffoni , g. , monaco , p. & theuns , t. 2002 , mnras , 333 , 623 uehara , h. & inutsuka , s. 2000 , apj , 531 , l91 venkatesan , a. , tumlinson , j. , shull , j.  m. 2003 , apj , 584 , 621 wyithe , j. s. b. & loeb , a. 2003 , apj , 588 , l69 yoshii , y. , sabano , y. 1980 , pasj , 32 , 229