Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: the study of the formation of population iii stars in a cosmological context via high - resolution simulations is becoming a mature discipline , with multiple groups finding similar results ( bromm , coppi , & larson 2002 ; bromm & loeb 2004 ; abel , bryan , & norman 2002 ; oshea & norman 2007 ; yoshida et  al . 2006b ; gao et  al . however , all of these calculations operate under a common fundamental assumption  namely , the absence of an ultraviolet background . given the abundance of halos in which population iii stars would form and the general consensus that these stars must be massive , one expects that , though only a small fraction of the volume of the universe would be ionized , a significant background of ultraviolet radiation in the lyman - werner ( lw ) band ( 11.18 - 13.6 ev ) , which is capable of photodissociating molecular hydrogen , would be present for the formation of the bulk of population iii stars ( omukai 2001 ; yoshida et  al . 2003 ; machacek , bryan , & abel 2001 ; wise & abel 2005 ) . given that atomic hydrogen is optically thin to this radiation and that it easily destroys molecular hydrogen , this could have a significant impact on the formation of primordial stars . ( 2001 ) used cosmological amr simulations to study the formation of primordial stars in the presence of a soft uv background , and later including an x - ray background ( machacek , bryan , & abel 2003 ) . they found that the lw background delays the formation of population iii stars and shifts halo formation to higher masses . however , their calculations were not of sufficiently high resolution to give any significant information about the formation of the protostellar cloud at the core of each collapsing cosmological halo . ( 2003 ) performed a suite of large sph simulations of high redshift structure formation , including two calculations with a soft uv background with @xmath6 and @xmath7 , where @xmath8 is the mean intensity of the uv background in the lw band in units of @xmath9  erg s@xmath10  @xmath11  hz@xmath10  sr@xmath10 . they found that gas cooling is suppressed in the higher lw background case , and predicted that star formation would not occur . on the other hand , & haiman ( 2002 ) showed analytically that primordial gas in halos with t@xmath12  k can still collapse to high densities in the presence of a strong molecular hydrogen - dissociating background via atomic hydrogen line radiation , and can eventually form h@xmath13 . they referred to these objects as  second generation objects " in the sense that although still of primordial composition , their formation pathway is different from the first stars . an issue that needs to be addressed concerns the mode of primordial star formation in halos with t@xmath5  k and above . it is apparent from recent studies using large numbers of simulations that smaller halos , with masses of @xmath14  m@xmath15 , appear to always form primordial stars in isolation  ( gao et  al . 2006 ; oshea & norman 2007 ) . is this still the case in halos whose masses are an order of magnitude or more larger ? does a single primordial star form , or several ? this is particularly relevant given that many semi - analytic models of structure formation that follow the chemical evolution of structures at high redshift ( scannapieco et  al . 2005 ; schneider et  al . 2006 ) suggest that halos of this scale which have not been enriched by metals ( and thus will form primordial stars ) exist up to at least @xmath16 . if this is true , there is expected to be a strong photodissociating background at that time  ( yoshida et  al . 2003 ; wise & abel 2005 ) . this paper is the second in a series . in & norman ( 2007 ) ( hereafter referred to as paper i ) , we examined several aspects of population iii star formation in a @xmath17cdm universe , in the absence of an ultraviolet background which photodissociates molecular hydrogen . in this paper , we study population iii star formation in a single cosmological realization varying the strength of the lw background , with the goal of investigating the effect that this background has on the evolution and properties of gas in the halo core ( at radii @xmath18  pc ) and the ultimate fate of the gas at the center of t@xmath5  k halos . we use a single cosmological realization , varying the strength of the photodissociating background over the range suggested by   & abel ( 2005 ) , and look for the cutoff in population iii star formation suggested by   ( 2003 ) .    in agreement with the findings of ( machacek et  al . 2001 ) , we find that the addition of a photodissociating background results in a delay of the collapse of high density gas at the center of the most massive halo in the simulation and , as a result , an increase in the virial mass of this halo at the onset of baryon collapse . we find that , contrary to the results of   ( 2003 ) , star formation is not suppressed for j@xmath0 , but occurs even with backgrounds as high as j@xmath1 . we find that 2 cooling leads to collapse despite the depressed halo core molecular hydrogen fractions f@xmath19 by two multiplicative effects : ( 1 ) the elevated 2 cooling rates per molecule at @xmath2 k , and ( 2 ) time . we find that halo core collapse occurs in the usual way once the gas in the halo core has become dense enough that the cooling time becomes much less than the hubble time . we also observe a relationship between the strength of the photodissociating background and the rate of accretion onto the evolving protostellar cloud core , with higher lw background fluxes resulting in higher accretion rates . this is a simple consequence of the suppression of molecular hydrogen formation ( and thus suppression of cooling ) by the photodissociating background , as well as the higher virial temperatures of the more massive halos at the epoch of collapse . this may have implications for the range of population iii stellar masses . finally , we find that the collapsing halo cores in our simulations do not fragment at densities below @xmath3  @xmath4 regardless of the strength of the soft uv background , suggesting that population iii stars forming in halos with t@xmath5  k may still form in isolation . the organization of this paper is as follows . in section  [ methodology ] we provide a description of enzo , the code used to perform the calculations in this paper and of the simulation setup . the results from our simulations are presented in sections  [ results.meanprop ] through [ results.rephalos ] : section  [ results.meanprop ] discusses some of the mean halo properties observed in the calculations , section  [ results.evolution ] discusses the evolution of the halo core prior to collapse , section  [ results.collapse ] compares spherically - averaged halo properties for all simulations at the epoch of collapse , section  [ results.fixred ] discusses a variety of halo properties at a fixed redshift , and section  [ results.rephalos ] compares the evolution of two representative simulations taken from our ensemble . in section  [ sect.issues ] we discuss neglected physics and possible numerical issues , and in section  [ discuss ] we discuss some of the results presented in this work and their implications . finally , we present a summary of the main results in section  [ summary ] . ` enzo ' is a publicly available , extensively tested adaptive mesh refinement ( amr ) cosmology code developed by greg bryan and others ( bryan & norman 1997a , 1997b ; norman & bryan 1999 ; oshea et  al . 2004 , 2005b ) . the specifics of the enzo code are described in detail in these papers ( and references therein ) , but we present a brief description here for clarity .    the enzo code couples an n - body particle - mesh ( pm ) solver ( efstathiou et  al . 1985 ; hockney & eastwood 1988 ) used to follow the evolution of a collisionless dark matter component with an eulerian amr method for ideal gas dynamics by & colella ( 1989 ) , which allows high dynamic range in gravitational physics and hydrodynamics in an expanding universe . this amr method ( referred to as _ structured _ amr ) utilizes an adaptive hierarchy of grid patches at varying levels of resolution . each rectangular grid patch ( referred to as a `` grid '' ) covers some region of space in its _ parent grid _ which requires higher resolution , and can itself become the parent grid to an even more highly resolved _ child grid_. enzo s implementation of structured amr places no fundamental restrictions on the number of grids at a given level of refinement , or on the number of levels of refinement . however , owing to limited computational resources it is practical to institute a maximum level of refinement , @xmath20 . additionally , the enzo amr implementation allows arbitrary integer ratios of parent and child grid resolution , though in general for cosmological simulations ( including the work described in this paper ) a refinement ratio of 2 is used . since the addition of more highly refined grids is adaptive , the conditions for refinement must be specified . in enzo , the criteria for refinement can be set by the user to be a combination of any or all of the following : baryon or dark matter overdensity threshold , minimum resolution of the local jeans length , local density gradients , local pressure gradients , local energy gradients , shocks , and cooling time . a cell reaching any or all of the user - specified criteria will then be flagged for refinement . once all cells of a given level have been examined , rectangular solid boundaries are determined which minimally encompass the flagged cells on that level . a refined grid patch is then introduced within each such bounding volume , and the results are interpolated to a higher level of resolution .    in enzo , resolution of the equations being solved is adaptive in time as well as in space . the timestep in enzo is satisfied on a level - by - level basis by finding the largest timestep such that the courant condition ( and an analogous condition for the dark matter particles ) is satisfied by every cell on that level . all cells on a given level are advanced using the same timestep . once a level @xmath21 has been advanced in time @xmath22 , all grids at level @xmath23 are advanced , using the same criteria for timestep calculations described above , until they reach the same physical time as the grids at level @xmath21 . at this point grids at level @xmath23 exchange baryon flux information with their parent grids , providing a more accurate solution on level @xmath21 . cells at level @xmath23 are then examined to see if they should be refined or de - refined , and the entire grid hierarchy is rebuilt at that level ( including all more highly refined levels ) . the timestepping and hierarchy rebuilding processes are repeated recursively on every level to the maximum existing grid level in the simulation . two different hydrodynamic methods are implemented in enzo : the piecewise parabolic method ( ppm ) ( woodward & colella 1984 ) , which was extended to cosmology by ( 1995 ) , and the hydrodynamic method used in the zeus magnetohydrodynamics code ( stone & norman 1992a , 1992b ) . we direct the interested reader to the papers describing both of these methods for more information , and note that ppm is the preferred choice of hydro method since it is higher - order - accurate and is based on a technique that does not require artificial viscosity , which smoothes shocks and can smear out features in the hydrodynamic flow .    the chemical and cooling properties of primordial ( metal - free ) gas are followed using the method of ( 1997 ) and ( 1997 ) . this method follows the non - equilibrium evolution of a gas of primordial composition with 9 total species : @xmath24 , @xmath25 , @xmath26 , @xmath27 , @xmath28 , @xmath29 , @xmath30 , @xmath31 , and @xmath32 . the code also calculates radiative heating and cooling following atomic line excitation , recombination , collisional excitation , free - free transitions , molecular line cooling , and compton scattering of the cosmic microwave background , as well as any of approximately a dozen different models for a metagalactic ultraviolet background that heat the gas via photoionization and/or photodissociation . we model the cooling processes detailed in   ( 1997 ) , but use the   & palla ( 1998 ) molecular hydrogen cooling function . the multispecies rate equations are solved out of equilibrium to properly model situations where , e.g. , the cooling time of the gas is much shorter than the hydrogen recombination time . a total of 9 kinetic equations are solved , including 29 kinetic and radiative processes , for the 9 species mentioned above . the chemical reaction equation network is technically challenging to solve due to the huge range of reaction time scales involved ; the characteristic creation and destruction time scales of the various species and reactions can differ by many orders of magnitude . as a result , the set of rate equations is extremely stiff , and an explicit scheme for integration of the rate equations can be costly if small enough timestep are taken to keep the network stable . this makes an implicit scheme preferable for such a set of equations , and enzo solves the rate equations using a method based on a backwards differencing formula ( bdf ) in order to provide a stable and accurate solution . it is important to note the regime in which this chemistry model is valid . according to ( 1997 ) and ( 1997 ) , the reaction network is valid for temperatures between @xmath33 k. the original model discussed in these two references is only valid up to n@xmath34  @xmath4 . however , addition of the 3-body h@xmath13 formation process allows correct modeling of the gas chemistry up until the point where collisionally induced emission from molecular hydrogen becomes an important cooling processes , which occurs at @xmath35  @xmath4 . we do not include heating by molecular hydrogen formation , which will be significant at densities of @xmath36  @xmath4 and above , and may affect temperature evolution at these high densities . a further concern is that the optically thin approximation for radiative cooling breaks down beginning at n@xmath37  @xmath4 . beyond this point , modifications to the cooling function that take into account the non - negligible opacity of the gas to line radiation from molecular hydrogen must be made , as discussed by & abel ( 2004 ) . even with these modifications , a more correct description of the cooling of gas of primordial composition at high densities will require some form of radiation transport , which will greatly increase the cost of the simulations . the simulations discussed in this paper are set up in a similar way to those in paper i. a dark matter - only calculation with @xmath38 particles in a three - dimensional simulation volume which is @xmath39  h@xmath10  mpc ( comoving ) on a side is set up at @xmath40 assuming a `` concordance '' cosmological model with no baryons : @xmath41 , @xmath42 , @xmath43 , @xmath44 ( in units of 100 km / s / mpc ) , @xmath45 , and using an eisenstein & hu power spectrum ( eisenstein & hu 1999 ) with a spectral index of @xmath46 . the cold dark matter cosmological model is assumed . this calculation is then evolved to @xmath47 using a maximum of four levels of adaptive mesh refinement , refining on a dark matter overdensity of 8.0 . at @xmath47 , the hop halo finding algorithm ( eisenstein & hut 1998 ) is used to find the most massive halo in the simulation .    at this point , we generate a new set of initial conditions which contain the same large - scale power as the dark matter - only calculation , but include both dark matter and baryons such that the lagrangian volume in which the halo in the dark matter - only calculation formed is resolved at high spatial and mass resolution using a series of static nested grids , with a @xmath38 root grid and three static nested grids , for an overall effective root grid size of @xmath48 cells . the highest resolution grid is @xmath49 grid cells , and corresponds to a volume @xmath50  h@xmath10 comoving kpc on a side . the dark matter particles in the highest resolution grid are 14.48  h@xmath10  m@xmath15 and the spatial resolution of the highest resolution grid is 586  h@xmath10 parsecs ( comoving ) . previous work shows that this particle mass resolution is adequate to fully resolve the collapse of the halo ( abel et  al . 2002 ; oshea & norman 2007 ) . all simulations are performed using the adaptive mesh cosmology code enzo , which is described in section  [ enzocode ] . the simulations are started at @xmath40 and allowed to evolve until the collapse of the gas within the center of the most massive halo , assuming the presence of an unevolving soft uv background with intensities in the lyman - werner band of j@xmath51 and  @xmath9  erg s@xmath10  @xmath11  hz@xmath10  sr@xmath10 . this range covers a much wider range of parameter space than the results described in ( 2001 ) , and completely encompasses the range of possible j@xmath52 values suggested by & abel ( 2005 ) , for a wide range of mean population iii stellar masses . note that many publications use f@xmath52 instead of j@xmath52 : f@xmath52 has units of erg  s@xmath10  @xmath11  hz@xmath10 , and thus f@xmath53  j@xmath52 . we will use j@xmath52 throughout this paper , and for convenience express values in units of j@xmath54 , where @xmath8 is the mean intensity of the uv background in the lw band in units of @xmath9  erg s@xmath10  @xmath11  hz@xmath10  sr@xmath10 . the equations of hydrodynamics are solved using the ppm method with a dual energy formulation , which is required to adequately resolve the thermal properties of gas in high - mach flows . the nonequilibrium chemical evolution and optically thin radiative cooling of the primordial gas is modeled as described in section  [ enzocode ] , following 9 separate species including molecular hydrogen ( but excluding deuterium ) , with an initial electron fraction of @xmath55 ( which is roughly consistent with ( 1968 ) for an @xmath56 , @xmath57 universe ) . note that the initial electron fraction in the calculation is relatively unimportant to the molecular hydrogen formation rates in halo cores , as the electron fraction at the center of a given halo is controlled primarily by mergers and the shock formed by accretion of gas onto the halo . adaptive mesh refinement is used throughout the innermost high resolution region such that cells are refined by factors of two along each axis , with a maximum of 22 total levels of refinement . this corresponds to a maximum spatial resolution of 115  h@xmath10 au ( comoving ) at the finest level of resolution , with an overall spatial dynamical range of @xmath58 . to avoid effects due to the finite size of the dark matter particles , the dark matter density is smoothed on a comoving scale of @xmath59  pc ( which corresponds to @xmath60 0.03 proper pc at @xmath61 ) . this is reasonable because at that scale in all of our calculations the gravitational potential in the halo of interest is completely dominated by the baryons . grid cells are adaptively refined based upon several criteria . we refine on baryon and dark matter overdensities in cells of 4.0 and 8.0 , respectively . this corresponds to a maximum mass of gas or dark matter per cell ( on the most highly refined static grid ) of m@xmath62 and @xmath63  m@xmath15 , respectively . in addition , the _ minimummassforrefinementlevelexponent _ parameter is set to @xmath64 for both the dark matter and baryon overdensity refinement criteria , meaning that the mass required to refine to a higher level decreases as a function of increasing level , as :    @xmath65    where l is the current level of refinement . the negative exponent means that the mass resolution in the calculations is super - lagrangian  for example , m@xmath66  m@xmath67 . in addition to refining on baryon and dark matter overdensity , these simulations include refinement criteria which ensure that shocks are always well - resolved , that the cooling time in a given cell is always longer than the sound crossing time of that cell , and that the jeans length is always resolved by at least 16 cells . this last criterion guarantees that the truelove criterion ( truelove et  al . 1997 ) is always resolved by a factor of four more cells in each dimension than is strictly necessary , ensuring that no artificial fragmentation will take place . figures  [ fig.image_jlw_0 ] and  [ fig.image_jlw_1em21 ] show projections of baryon density and temperature for the j@xmath68 and j@xmath1 calculations , respectively , at the epoch at which each calculation collapses , defined as the redshift at which the baryon number density reaches @xmath69  @xmath4 . note that due to the rapid evolution of gas at high density , the `` collapse redshift '' depends very weakly on the exact choice of density threshold . the calculations are started from the same set of initial conditions , and the j@xmath1 calculation is clearly a later stage in the evolution of the j@xmath68 run . the satellite halos surrounding the primary halo which collapses at z=24.12 in the j@xmath68 run have merged by z=17.32 when the halo in the j@xmath1 case collapses ( cf . [ fig.merger ] ) . with virial masses of @xmath70 and @xmath71 , respectively , the latter halo is more than 20 times as massive as the former . though quite different in mass , both halos exhibit similar morphologies  they are extremely centrally - concentrated , and only a single condensed object ( that is to say , a primordial protostellar cloud core ) is visible in the highest - resolution panel in each image . the existence of a single cloud core is common to all calculations discussed in this paper . simulation ( @xmath72 ; m@xmath73  m@xmath15 ) . top row : projected log baryon density . bottom row : projected , mass - weighted log baryon temperature . left column : region 2.13 kpc ( proper ) across and deep . middle column : region 133.3 pc ( proper ) across and deep . right column : region 0.130 pc ( proper ) across and deep . the middle and right column are zoomed images focusing on the high baryon density core which forms in the center of the halo . , scaledwidth=90.0% ]     simulation ( @xmath74 ; m@xmath75  m@xmath15 ) . top row : projected log baryon density . bottom row : projected , mass - weighted log baryon temperature . left column : region 2.13 kpc ( proper ) across and deep . middle column : region 533 pc ( proper ) across and deep . right column : region 0.260 pc ( proper ) across and deep . the middle and right column are zoomed images focusing on the high baryon density core which forms in the center of the halo . , scaledwidth=90.0% ]    figure  [ fig.halovals1 ] shows several mean halo quantities for all simulations discussed in this paper , including halo collapse redshift as a function of j@xmath52 , virial mass as a function of j@xmath52 , virial temperature as a function of j@xmath52 , and the virial mass as a function of collapse redshift . simulations where the soft uv background is turned on are shown by solid squares , and the `` control '' j@xmath68 calculation is shown by an open square . in plots where j@xmath52 is shown on the x - axis , the j@xmath68 simulation is placed at log@xmath76 . figure  [ fig.halovals1 ] shows that there is a clear relationship between the lw intensity and the collapse redshift and virial mass of the halo . a larger lw intensity results in a later collapse time and larger virial mass because the halo must be hotter in order to have a cooling time which is less than a hubble time despite the depressed 2 abundance in the high - density gas at the halo s center . this is discussed further in section  [ results.evolution ] the final mass of the halos in the simulations with j@xmath1 is approximately a factor of 20 higher than that in the `` control '' simulation . this is in qualitative agreement with ( 2001 ) , who suggest that there is a `` minimum halo mass '' as a function of the strength of the lw background of the form    @xmath77    this threshold mass is plotted in panel ( c ) of figure  [ fig.halovals1 ] . [ eqn - mach - mthresh ] is only strictly valid over the range @xmath78 , because this was the range simulated by ( 2001 ) . we see that our points for @xmath79 parallel the threshold curve but at a mass approximately four times higher . this difference can be ascribed to two factors . first , the threshold mass is the lowest possible halo mass that can collapse , derived from a statistical sample . since we study only one realization which focuses on the most massive halo in the box , its mass is bound to be higher than the statistical minimum . second , ( 2001 ) s criterion for cooling catches halos at an earlier stage of evolution compared to our data points , which give the halos virial masses at the time of central baryon collapse ( the `` collapse redshift '' ) . examination of panels ( b ) and ( c ) of figure  [ fig.halovals1 ] show that there is some sort of `` phase change '' between j@xmath80 and @xmath81 causing the halo mass at the time of collapse to increase steeply and non - monotonically . the collapse redshift steadily decreases as a function of uv increasing background strength ( i.e. collapse of gas at the center of the halo is delayed ) . however , panel ( c ) shows that the virial mass steadily increases with increasing uv background strength until j@xmath80 , at which point there is a jump in mass by more than a factor of four , above which the halo mass increases only slightly with increasing uv background strength . that something interesting should happen at these lw intensities is consistent with the findings of   ( 2003 ) who found that 2 cooling is strongly suppressed at j@xmath82 . they predicted that star formation would be inhibited since the equilibrium 2 fraction is below the critical fraction for cooling derived by ( 1997 ) . contrary to these predictions , we find collapse not only at j@xmath82 , but also at values as large as j@xmath83 . we note that highest value considered by ( 2001 ) was j@xmath84 , and the highest simulated by   ( 2003 ) was j@xmath82 . therefore our cases j@xmath85 have not been examined before , and certainly not at the resolution of our simulations . how is the gas in the center of these halos able to cool and collapse in such high radiation backgrounds ? this is analyzed in some detail in the next section . a hint is provided in ( d ) of figure  [ fig.halovals1 ] , which plots the halo virial temperature as a function of j@xmath52 . at background strengths above j@xmath80 , the virial temperature is consistently approximately @xmath86  k , which is roughly the temperature at which atomic hydrogen cooling is effective . however , an examination of the radial temperature profiles shows that high - density gas in the center of the halo never reaches these temperatures , but is more typically 2000 k. atomic line cooling is unimportant at such low temperatures ; however the 2 cooling rate per particle is roughly 100 times as large at 2000k as at 500k . 2 cooling still operates in the halo centers despite low 2 abundance due to the higher cooling rates and long evolutionary timescales . in the rest of this section we merely present additional mean properties of the halo and its central region at collapse , and defer discussion of the relevant timescales to section  [ results.evolution ] .     vs. halo collapse redshift . panel ( b ) : halo virial mass vs. halo collapse redshift panel ( c ) : halo virial mass vs. j@xmath52 panel ( d ) : halo virial temperature vs. j@xmath52 the j@xmath68 `` control '' result are shown as an open square ( and is at log j@xmath87 in the panels which are a function of j@xmath52 ) . in the bottom left panel , the dashed line corresponds to the fitting function for threshold mass from ( 2001 ) , eqtn . 8.,scaledwidth=90.0% ]    figure  [ fig.halovals2 ] shows several properties of the halo core at the epoch of collapse , including the core temperature , molecular hydrogen fraction , and spherically - averaged accretion rate as a function of ultraviolet background strength , and the spherically - averaged accretion rate as a function of the molecular hydrogen fraction . all `` core '' values are spherically - averaged and measured at the mass shell where @xmath88  m@xmath15 of gas is enclosed . panel ( a ) shows that the core h@xmath13 fraction decreases significantly with increasing fuv flux , with a corresponding increase in the core temperature ( panel ( b ) ) . this relationship is similar to that noted in   & norman ( 2007 ) , where the amount of molecular hydrogen at densities of @xmath89  @xmath4 varies between simulations , and correlates strongly with baryon temperature . given that the accretion of gas onto the protostellar cloud is subsonic , this results in a strong relationship between the soft uv background flux and accretion rate onto the protostellar cloud , as shown by panel ( c ) , with accretion rates varying by more than a factor of 30 between the j@xmath68 and @xmath90 cases . panel ( d ) shows the strength of the correlation between core h@xmath13 fraction and accretion rate . note that the values discussed above are insensitive to the exact definition of the halo `` core . '' core for several simulations with the same cosmic realization but a range of lw background intensities . panel ( a ) : j@xmath52 vs. baryon core temperature . panel ( b ) : j@xmath52 vs. baryon core @xmath31 fraction . panel ( c ) : j@xmath52 vs. instantaneous accretion rate . panel ( d ) : baryon core @xmath31 fraction vs. accretion rate . the j@xmath68 `` control '' result are shown as an open square ( and is at log j@xmath87 in the panels which are a function of j@xmath52 ) . all `` core '' values are spherically - averaged and measured at the mass shell where @xmath88  m@xmath15 of gas is enclosed . , scaledwidth=90.0% ]    figure  [ fig.halovals3 ] shows the halo core mass and mass fraction at the epoch of collapse ( note that the definition of `` core '' is somewhat different than in the previous figure ) . the halo `` core mass '' is defined as being all gas at a density of n@xmath91  @xmath4 or above , and the core mass fraction is defined as the core mass divided by the virial mass of the halo at the epoch of collapse . this choice of minimum density ensures that we have captured the entirety of the quasistatically contracting analog of a galactic molecular cloud core , and for practical purposes is comparable to the gas fraction of cold , dense gas discussed in   ( 2001 ) . the size of the halo core increases in absolute terms as the strength of the photodissociating background is increased , from @xmath92  m@xmath15 until j@xmath93 , where it reaches @xmath94  m@xmath15 . above this value of j@xmath54 , the core mass then begins to decrease , reaching @xmath95  m@xmath15 at j@xmath1 . the core mass fraction generally decreases ( though with some noise in the relationship ) from @xmath96 at j@xmath68 to @xmath97 at j@xmath1 , though the total overall change is less than a factor of 2 . @xmath4 or above . the `` core mass fraction '' is defined as the `` core mass '' divided by the virial mass of the halo at the epoch of collapse . the j@xmath98 `` control '' result are shown as an open square ( and is at log j@xmath87 in the panels which are a function of j@xmath52 ) . , scaledwidth=90.0% ] we now examine in more detail the evolution of the conditions in the cloud center leading to collapse . in figure  [ fig.evolution ] we plot the central density , temperature , entropy @xmath99 , and 2 fraction versus time for the calculations with @xmath100 ( solid line ) , @xmath101 ( short - dashed line ) , @xmath81 ( long - dashed line ) and @xmath102 ( dot - dashed line ) . in these plots , central values of temperature and density are defined as being the values of these quantities in the cell with the highest baryon density . it is most instructive to compare the two extreme cases , j@xmath103 , which very closely resembles the j@xmath104 case , and the j@xmath83 case , which according to   ( 2003 ) should not collapse at all . in panel ( a ) we see the central density increase modestly from @xmath105 2 @xmath4 to @xmath105 20 @xmath4 over 40 myr , and then increase rapidly thereafter , reaching @xmath106 @xmath4 in a scant additional 15 myr . this evolution is driven by 2 cooling , as can be seen from the entropy evolution in panel ( c ) . as expected , runaway cooling occurs when the 2 fraction reaches @xmath107 ( panel d ) , consistent with the ( 1997 ) analysis . by contrast the j@xmath83 case requires 120 myr for the central density to increases from @xmath105 2 @xmath4 to @xmath105 20 @xmath4 . as shown in panel ( c ) , the first 50 myr of this contraction is adiabatic , followed by an increase in entropy due to some heating event ( mergers ) . the fact that central entropy is not decreasing for 90 myr is indicative of the fact that 2 cooling is unimportant over this interval due to extremely low equilibrium 2 fractions ( f@xmath108 ; ( panel d ) ) . however , something interesting happens in the j@xmath1 calculation at @xmath109 myr . prior to that time the central temperature has crept up to 5000 k due to the increased virial mass . the cooling rate per 2 molecule is about 3 orders of magnitude higher at 5000k than at 500k  a typical temperature at the halo center in low uv background evolutions . we believe this high temperature is due to a combination of adiabatic heating as the potential well deepens , and merger - induced shock heating . as shown in figure [ fig.merger ] , three large halos merge between @xmath110 and @xmath111 , quadrupling the halo mass ( figure [ fig.halovals1 ] , panel ( b ) ) . the elevated cooling rates cool the gas at the halo center to @xmath112 , allowing the gas density to increase slightly . the two effects reduce the central entropy from @xmath113 700 k @xmath114 to @xmath115 250 k @xmath114 . this is followed by another heating event at t=190 myr , presumably due to another merger , followed by a second more catastrophic cooling event . at t=200 myr , the halo core collapses , driving the central density higher , temperature and entropy lower , and 2 fraction higher . why does this collapse occur ? this is analyzed in figure [ fig.timescales ] , where we plot the evolution of the important timescales in the center of the halo . in the j@xmath103 case , the 2-body 2 formation timescale is always roughly an order of magnitude smaller than the photodissociation timescale . consequently , the 2 abundance is out of equilibrium and grows steadily with time ( figure [ fig.evolution ] panel d ) . the 2 cooling time in the central zone ( @xmath116 myr ) is far less than the hubble time ( @xmath117 myr ) , and before too long runaway collapse occurs . at higher values of j@xmath52 , the 2-body 2 formation time always hovers around the photodissociation timescale , regardless of how long that timescale is . this tells us that 2 is in equilibrium : the 2-body formations balance the photo - destructions by the lw background . the equilibrium 2 fraction is given by    @xmath118    where k@xmath119 is the rate coefficient for the formation of h@xmath120 ( the limiting reaction for h@xmath13 formation at the temperature and density range considered here ) , @xmath121 is the electron density , and k@xmath122 is the 2 photodissociation rate due to the solomon process  ( abel et  al . 1997 ; yoshida et  al . 2003 ) . using t @xmath123  k and @xmath124 , we get @xmath125 , which is in good agreement with the 2 fractions we see in figure [ fig.evolution ] . since @xmath126/ n@xmath127 , the smaller @xmath128 is , the shorter its formation time . now let us consider the 2 cooling time evolution ( short - long - dashed line ) . it is far lass than a hubble time for j@xmath129 , ensuring that these halo cores will eventually cool and collapse . in the j@xmath83 case , the cooling time drops below the hubble time for the first time at t=170 myr , and again at t=190 myr where it begins a steady decrease toward @xmath130 yr . despite the low equilibrium 2 fraction , if we wait long enough the inexorable press of time eventually establishes the condition @xmath131 and we get runaway collapse , even for j@xmath83 . note that we get catastrophic cooling because we are in the density regime ( @xmath132 @xmath4 ) where 2 cooling is proportional to the square of the gas density , not linear in the gas density  ( abel , bryan , & norman 2000 ; abel et  al . 2002 ) . ( solid ) , @xmath101 ( short - dashed ) , @xmath81 ( long - dashed ) and @xmath90 ( dot - dashed ) . , scaledwidth=90.0% ]     and @xmath90 ( panels ( a)(d ) , respectively ) : 2-body 2 formation time ( solid ) ; 3-body 2 formation time ( short - dashed ) ; 2 collisional dissociate time ( long - dashed ) ; 2 photodissociation time ( dot - dashed ) ; 2 cooling time ( short - long - dashed ) . , scaledwidth=90.0% ]    , t@xmath133  k halo that collapses at z=17.32 with a lyman - werner background flux @xmath134 . in comic strip order , the redshifts are z=22 , 21 , 20 , 19 , 18 , and 17.32 . field of view is 53.571 comoving kpc ( 2.678 proper kpc at z = 19 ) . logarithm of the projected baryon density ( column density ) is displayed , with color table scaled to the maximum and minimum values in each image . , scaledwidth=90.0% ] figure  [ fig.radprof_zcoll.1 ] shows radial profiles of baryon number density , temperature , enclosed mass , circular velocity , and rms mach number as a function of radius , as well as specific angular momentum as a function of enclosed mass , for all simulations discussed in this paper . to facilitate comparison , results from the output of each simulation where the peak baryon number density at the center of the halo is approximately @xmath135 @xmath4 are shown , in order to capture each halo at a similar evolutionary stage , rather than at a fixed point in time . panel ( a ) shows that all of the simulations have similar density profiles . the scatter in number density at a given radius is readily explained by variation in halo mass  halos which collapse at later times ( due to higher uv background strength ) are more massive , and thus have higher overall baryon densities . this is shown in another way in panel ( c ) , in a plot of enclosed mass as a function of radius which shows that all halos have very similar profiles at @xmath136 pc , but a significant variation at larger radii which is related to halo mass . the plot of temperature as a function of radius in panel ( b ) shows an interesting trend  as the lw background is increased , the overall halo temperature as well as the halo core ( where `` core '' is roughly defined as gas within @xmath137 pc of the halo center ) temperature go up . in the outskirts of the halo , where baryon densities are low and thus cooling times are long , this is due primarily to the increase in halo mass , with the peak temperature corresponding approximately to the halo virial temperature . in the halo core , however , this temperature is correlated more strongly with the h@xmath13 fraction , and thus the uv background , as suggested by figure  [ fig.halovals2 ] . the plots of specific angular momentum as a function of enclosed mass and keplerian velocity fraction ( defined as the local circular velocity divided by the keplerian orbital velocity due to mass within that radius ) as a function of radius in panels ( d ) and ( e ) show some mild trends trends . it appears that halos in simulations with stronger uv backgrounds tend to have less specific angular momentum at a given mass shell , and also tend to have a lower keplerian velocity fraction at a given radius than gas in halos which form in the presence of a lower uv background  indeed , the halo which forms in the j@xmath1 calculation has the least angular momentum of all of the simulations , and the lowest keplerian velocity fraction out to @xmath137  pc . we speculate that this is due to transport and segregation of angular momentum by turbulence  halos whose collapse is delayed have more time for these processes to act within the halo core , resulting in less angular momentum overall when the halo ultimately collapses ( note that this will be discussed more fully in a later paper ) . the plot of rms mach number as a function of radius in panel ( f ) also shows no obvious trend , though at radii below @xmath138 pc , simulations with a smaller uv background generally have higher rms mach numbers . this is due primarily to the gas being colder , and not the gas velocities being higher .    . panel ( a ) : @xmath139 as a function of radius . panel ( b ) : temperature as a function of radius . panel ( c ) : enclosed gas mass as a function of radius . panel ( d ) : specific angular momentum as a function of enclosed mass . panel ( e ) : keplerian velocity fraction as a function of radius . panel ( f ) : gas rms mach number as a function of radius . all quantities except enclosed gas mass are mass - weighted , and all simulations are shown at the point where the maximum number density is approximately @xmath135  @xmath4 . line types and weights correspond to simulations , as follows . black solid line : j@xmath68 . black short - dashed line : j@xmath140 . black long - dashed line : j@xmath141 . black dot short - dashed line : j@xmath142 . black dot long - dashed line : j@xmath80 . black short dashed - long dashed line : j@xmath143 . red solid line : j@xmath144 . red short - dashed line : j@xmath145 . red long - dashed line : j@xmath1 . , title="fig:",scaledwidth=30.0% ] . panel ( a ) : @xmath139 as a function of radius . panel ( b ) : temperature as a function of radius . panel ( c ) : enclosed gas mass as a function of radius . panel ( d ) : specific angular momentum as a function of enclosed mass . panel ( e ) : keplerian velocity fraction as a function of radius . panel ( f ) : gas rms mach number as a function of radius . all quantities except enclosed gas mass are mass - weighted , and all simulations are shown at the point where the maximum number density is approximately @xmath135  @xmath4 . line types and weights correspond to simulations , as follows . black solid line : j@xmath68 . black short - dashed line : j@xmath140 . black long - dashed line : j@xmath141 . black dot short - dashed line : j@xmath142 . black dot long - dashed line : j@xmath80 . black short dashed - long dashed line : j@xmath143 . red solid line : j@xmath144 . red short - dashed line : j@xmath145 . red long - dashed line : j@xmath1 . , title="fig:",scaledwidth=30.0% ] . panel ( a ) : @xmath139 as a function of radius . panel ( b ) : temperature as a function of radius . panel ( c ) : enclosed gas mass as a function of radius . panel ( d ) : specific angular momentum as a function of enclosed mass . panel ( e ) : keplerian velocity fraction as a function of radius . panel ( f ) : gas rms mach number as a function of radius . all quantities except enclosed gas mass are mass - weighted , and all simulations are shown at the point where the maximum number density is approximately @xmath135  @xmath4 . line types and weights correspond to simulations , as follows . black solid line : j@xmath68 . black short - dashed line : j@xmath140 . black long - dashed line : j@xmath141 . black dot short - dashed line : j@xmath142 . black dot long - dashed line : j@xmath80 . black short dashed - long dashed line : j@xmath143 . red solid line : j@xmath144 . red short - dashed line : j@xmath145 . red long - dashed line : j@xmath1 . , title="fig:",scaledwidth=30.0% ] . panel ( a ) : @xmath139 as a function of radius . panel ( b ) : temperature as a function of radius . panel ( c ) : enclosed gas mass as a function of radius . panel ( d ) : specific angular momentum as a function of enclosed mass . panel ( e ) : keplerian velocity fraction as a function of radius . panel ( f ) : gas rms mach number as a function of radius . all quantities except enclosed gas mass are mass - weighted , and all simulations are shown at the point where the maximum number density is approximately @xmath135  @xmath4 . line types and weights correspond to simulations , as follows . black solid line : j@xmath68 . black short - dashed line : j@xmath140 . black long - dashed line : j@xmath141 . black dot short - dashed line : j@xmath142 . black dot long - dashed line : j@xmath80 . black short dashed - long dashed line : j@xmath143 . red solid line : j@xmath144 . red short - dashed line : j@xmath145 . red long - dashed line : j@xmath1 . , title="fig:",scaledwidth=30.0% ] . panel ( a ) : @xmath139 as a function of radius . panel ( b ) : temperature as a function of radius . panel ( c ) : enclosed gas mass as a function of radius . panel ( d ) : specific angular momentum as a function of enclosed mass . panel ( e ) : keplerian velocity fraction as a function of radius . panel ( f ) : gas rms mach number as a function of radius . all quantities except enclosed gas mass are mass - weighted , and all simulations are shown at the point where the maximum number density is approximately @xmath135  @xmath4 . line types and weights correspond to simulations , as follows . black solid line : j@xmath68 . black short - dashed line : j@xmath140 . black long - dashed line : j@xmath141 . black dot short - dashed line : j@xmath142 . black dot long - dashed line : j@xmath80 . black short dashed - long dashed line : j@xmath143 . red solid line : j@xmath144 . red short - dashed line : j@xmath145 . red long - dashed line : j@xmath1 . , title="fig:",scaledwidth=30.0% ] . panel ( a ) : @xmath139 as a function of radius . panel ( b ) : temperature as a function of radius . panel ( c ) : enclosed gas mass as a function of radius . panel ( d ) : specific angular momentum as a function of enclosed mass . panel ( e ) : keplerian velocity fraction as a function of radius . panel ( f ) : gas rms mach number as a function of radius . all quantities except enclosed gas mass are mass - weighted , and all simulations are shown at the point where the maximum number density is approximately @xmath135  @xmath4 . line types and weights correspond to simulations , as follows . black solid line : j@xmath68 . black short - dashed line : j@xmath140 . black long - dashed line : j@xmath141 . black dot short - dashed line : j@xmath142 . black dot long - dashed line : j@xmath80 . black short dashed - long dashed line : j@xmath143 . red solid line : j@xmath144 . red short - dashed line : j@xmath145 . red long - dashed line : j@xmath1 . , title="fig:",scaledwidth=30.0% ]    figure  [ fig.radprof_zcoll.2 ] shows the h@xmath13 fraction , electron fraction , h@xmath120 fraction , ratio of cooling time to sound crossing time , ratio of cooling time to dynamical time , and ratio of sound crossing time to dynamical time as a function of radius , for all simulations . outputs and line types correspond to those in figure  [ fig.radprof_zcoll.1 ] . there is a clear relationship at all radii between h@xmath13 fraction and uv background strength in panel ( a )  as the fuv background is increased , the overall h@xmath13 fraction decreases . this difference is most noticeable at radii of @xmath146 pc , but is maintained at larger and smaller radii . there is a `` kink '' at approximately @xmath101 pc where the values are all quite similar . this corresponds to a baryon number density of @xmath36  @xmath4 in all simulations , which is where 3-body h@xmath13 formation begins to occur . panels ( b ) and ( c ) show that e@xmath120 and @xmath29 fractions track each other , which is to be expected  the local electron fraction controls the amount of h@xmath120 which can be produced , which is the limiting reaction in the dominant mode of h@xmath13 formation for n @xmath147 @xmath4  and there is a general trend with increasing electron and h@xmath120 fractions with increasing uv background flux . the plots of the ratio of cooling time to sound crossing time , ratio of cooling time to dynamical time , and ratio of sound crossing time to dynamical time as a function of radius shown in panels ( d ) , ( e ) and ( f ) display some interesting trends . the three simulations with the lowest uv background strengths have a much longer cooling time than the other calculations , which are all grouped roughly together with no discernible pattern . this agrees well with the plot of the ratio of cooling time to dynamical time , where a similar trend is observed . the difference between the three simulations with the lowest uv background strengths and the others is due to the somewhat higher h@xmath13 fraction in the halo core of these calculations . the gas temperature is @xmath148  k at the center of the halo cores in these calculations , and the gas can not cool further , which results in a very long cooling time . the sound crossing time is also increased , but not as significantly . all simulations have approximately the same ratio of sound crossing to dynamical times at all radii . given that the cooling time is typically longer than both the sound crossing time and dynamical time , one can infer that the collapse of the halo is occurring quasi - statically for all simulations .    . panel ( a ) : molecular hydrogen fraction as a function of radius . panel ( b ) : electron fraction as a function of radius . panel ( c ) : @xmath29 fraction as a function of radius . panel ( d ) : ratio of gas cooling time to sound crossing time as a function of radius . panel ( e ) : ratio of gas cooling time to dynamical time as a function of radius . panel ( f ) : ratio of gas sound crossing time to dynamical time as a function of radius . all quantities except enclosed gas mass are mass - weighted , and all simulations are shown at the point where the maximum number density is approximately @xmath135  @xmath4 . line types and weights correspond to those in figure  [ fig.radprof_zcoll.1 ] . , title="fig:",scaledwidth=30.0% ] . panel ( a ) : molecular hydrogen fraction as a function of radius . panel ( b ) : electron fraction as a function of radius . panel ( c ) : @xmath29 fraction as a function of radius . panel ( d ) : ratio of gas cooling time to sound crossing time as a function of radius . panel ( e ) : ratio of gas cooling time to dynamical time as a function of radius . panel ( f ) : ratio of gas sound crossing time to dynamical time as a function of radius . all quantities except enclosed gas mass are mass - weighted , and all simulations are shown at the point where the maximum number density is approximately @xmath135  @xmath4 . line types and weights correspond to those in figure  [ fig.radprof_zcoll.1 ] . , title="fig:",scaledwidth=30.0% ] . panel ( a ) : molecular hydrogen fraction as a function of radius . panel ( b ) : electron fraction as a function of radius . panel ( c ) : @xmath29 fraction as a function of radius . panel ( d ) : ratio of gas cooling time to sound crossing time as a function of radius . panel ( e ) : ratio of gas cooling time to dynamical time as a function of radius . panel ( f ) : ratio of gas sound crossing time to dynamical time as a function of radius . all quantities except enclosed gas mass are mass - weighted , and all simulations are shown at the point where the maximum number density is approximately @xmath135  @xmath4 . line types and weights correspond to those in figure  [ fig.radprof_zcoll.1 ] . , title="fig:",scaledwidth=30.0% ] . panel ( a ) : molecular hydrogen fraction as a function of radius . panel ( b ) : electron fraction as a function of radius . panel ( c ) : @xmath29 fraction as a function of radius . panel ( d ) : ratio of gas cooling time to sound crossing time as a function of radius . panel ( e ) : ratio of gas cooling time to dynamical time as a function of radius . panel ( f ) : ratio of gas sound crossing time to dynamical time as a function of radius . all quantities except enclosed gas mass are mass - weighted , and all simulations are shown at the point where the maximum number density is approximately @xmath135  @xmath4 . line types and weights correspond to those in figure  [ fig.radprof_zcoll.1 ] . , title="fig:",scaledwidth=30.0% ] . panel ( a ) : molecular hydrogen fraction as a function of radius . panel ( b ) : electron fraction as a function of radius . panel ( c ) : @xmath29 fraction as a function of radius . panel ( d ) : ratio of gas cooling time to sound crossing time as a function of radius . panel ( e ) : ratio of gas cooling time to dynamical time as a function of radius . panel ( f ) : ratio of gas sound crossing time to dynamical time as a function of radius . all quantities except enclosed gas mass are mass - weighted , and all simulations are shown at the point where the maximum number density is approximately @xmath135  @xmath4 . line types and weights correspond to those in figure  [ fig.radprof_zcoll.1 ] . , title="fig:",scaledwidth=30.0% ] . panel ( a ) : molecular hydrogen fraction as a function of radius . panel ( b ) : electron fraction as a function of radius . panel ( c ) : @xmath29 fraction as a function of radius . panel ( d ) : ratio of gas cooling time to sound crossing time as a function of radius . panel ( e ) : ratio of gas cooling time to dynamical time as a function of radius . panel ( f ) : ratio of gas sound crossing time to dynamical time as a function of radius . all quantities except enclosed gas mass are mass - weighted , and all simulations are shown at the point where the maximum number density is approximately @xmath135  @xmath4 . line types and weights correspond to those in figure  [ fig.radprof_zcoll.1 ] . , title="fig:",scaledwidth=30.0% ]    figure  [ fig.radprof_zcoll.3 ] shows radial velocity as a function of radius , accretion time ( defined as @xmath149 ) as a function of enclosed baryon mass @xmath150 as a function of enclosed baryon mass and @xmath150 as a function of time for all simulations . output time and lines correspond to those in figures  [ fig.radprof_zcoll.1 ] and  [ fig.radprof_zcoll.2 ] . panel ( a ) shows a clear relationship between the strength of the soft uv background and the radial velocity , with the calculations that have higher uv backgrounds also having higher infall velocities . this is easily understood by the quasistatic contraction of the gas , which takes place at or below the local sound speed in the gas , as discussed by ( 2002 ) and & norman ( 2007 ) . given that the halo temperatures are systematically higher in calculations where the uv background is stronger , the sound speed is higher , and thus the rate at which the gas at the center of the halo contracts is higher . this leads to lower accretion times ( panel ( b ) ) and higher accretion rates ( panel ( c ) ) , with the accretion rates for the simulations with the strongest uv backgrounds being higher than those with a low ( or no ) uv background . panel ( d ) shows that the time evolution of the accretion rate in all simulations is qualitatively similar , but varies strongly in absolute magnitude and in the time at which the accretion rate onto the evolving protostellar cloud peaks , suggesting that the stars forming in these simulations may have very different evolutionary histories .    . panel ( a ) : radial velocity as a function of radius . panel ( b ) : accretion time as a function of enclosed baryon mass . panel ( c ) : instantaneous accretion rate as a function of enclosed baryon mass . panel ( d ) : estimated accretion rate as a function of time . all quantities except enclosed gas mass are mass - weighted , and all simulations are shown at the point where the maximum number density is approximately @xmath135  @xmath4 . line types and weights correspond to those in figure  [ fig.radprof_zcoll.1 ] . the upper and lower light short - long - dashed curves which extend from the upper left corner of panel ( b ) correspond to the main sequence lifetime of a massive population iii star of that mass and the kelvin - helmholtz timescale of a population iii timescale with a given luminosity and radius . all values are taken from   ( 2002 ) . the three light diagonal short - dashed lines which extend from bottom left to top right in panel ( b ) correspond to masses accreted using constant accretion rates of ( from top to bottom ) @xmath151 , @xmath101 , and @xmath81 m@xmath152yr . the light horizontal short - dashed line in panels ( c ) and ( d ) correspond to the `` critical '' accretion rate of   & palla ( 2003 ) , @xmath153  m@xmath152yr . , title="fig:",scaledwidth=45.0% ] . panel ( a ) : radial velocity as a function of radius . panel ( b ) : accretion time as a function of enclosed baryon mass . panel ( c ) : instantaneous accretion rate as a function of enclosed baryon mass . panel ( d ) : estimated accretion rate as a function of time . all quantities except enclosed gas mass are mass - weighted , and all simulations are shown at the point where the maximum number density is approximately @xmath135  @xmath4 . line types and weights correspond to those in figure  [ fig.radprof_zcoll.1 ] . the upper and lower light short - long - dashed curves which extend from the upper left corner of panel ( b ) correspond to the main sequence lifetime of a massive population iii star of that mass and the kelvin - helmholtz timescale of a population iii timescale with a given luminosity and radius . all values are taken from   ( 2002 ) . the three light diagonal short - dashed lines which extend from bottom left to top right in panel ( b ) correspond to masses accreted using constant accretion rates of ( from top to bottom ) @xmath151 , @xmath101 , and @xmath81 m@xmath152yr . the light horizontal short - dashed line in panels ( c ) and ( d ) correspond to the `` critical '' accretion rate of   & palla ( 2003 ) , @xmath153  m@xmath152yr . , title="fig:",scaledwidth=45.0% ] . panel ( a ) : radial velocity as a function of radius . panel ( b ) : accretion time as a function of enclosed baryon mass . panel ( c ) : instantaneous accretion rate as a function of enclosed baryon mass . panel ( d ) : estimated accretion rate as a function of time . all quantities except enclosed gas mass are mass - weighted , and all simulations are shown at the point where the maximum number density is approximately @xmath135  @xmath4 . line types and weights correspond to those in figure  [ fig.radprof_zcoll.1 ] . the upper and lower light short - long - dashed curves which extend from the upper left corner of panel ( b ) correspond to the main sequence lifetime of a massive population iii star of that mass and the kelvin - helmholtz timescale of a population iii timescale with a given luminosity and radius . all values are taken from   ( 2002 ) . the three light diagonal short - dashed lines which extend from bottom left to top right in panel ( b ) correspond to masses accreted using constant accretion rates of ( from top to bottom ) @xmath151 , @xmath101 , and @xmath81 m@xmath152yr . the light horizontal short - dashed line in panels ( c ) and ( d ) correspond to the `` critical '' accretion rate of   & palla ( 2003 ) , @xmath153  m@xmath152yr . , title="fig:",scaledwidth=45.0% ] . panel ( a ) : radial velocity as a function of radius . panel ( b ) : accretion time as a function of enclosed baryon mass . panel ( c ) : instantaneous accretion rate as a function of enclosed baryon mass . panel ( d ) : estimated accretion rate as a function of time . all quantities except enclosed gas mass are mass - weighted , and all simulations are shown at the point where the maximum number density is approximately @xmath135  @xmath4 . line types and weights correspond to those in figure  [ fig.radprof_zcoll.1 ] . the upper and lower light short - long - dashed curves which extend from the upper left corner of panel ( b ) correspond to the main sequence lifetime of a massive population iii star of that mass and the kelvin - helmholtz timescale of a population iii timescale with a given luminosity and radius . all values are taken from   ( 2002 ) . the three light diagonal short - dashed lines which extend from bottom left to top right in panel ( b ) correspond to masses accreted using constant accretion rates of ( from top to bottom ) @xmath151 , @xmath101 , and @xmath81 m@xmath152yr . the light horizontal short - dashed line in panels ( c ) and ( d ) correspond to the `` critical '' accretion rate of   & palla ( 2003 ) , @xmath153  m@xmath152yr . , title="fig:",scaledwidth=45.0% ] figure  [ fig.radprof_fixred.1 ] shows the status of several spherically - averaged quantities for halos from all simulations at a fixed redshift , @xmath154 , which is shortly before the baryons in the most massive halo in the j@xmath68 simulation collapses to high density . this is at an early point in the evolution of the majority of these halos , and several conclusions can be drawn from the plots of number density , temperature , h@xmath13 fraction , and radial velocity as a function of radius . the molecular hydrogen fraction declines monotonically with increasing uv background strength , in agreement with an estimate for equilibrium h@xmath13 fraction at low ( @xmath155  @xmath4 ) densities :    @xmath156    where k@xmath119 is the rate coefficient for the formation of h@xmath120 ( the limiting reaction for h@xmath13 formation at the temperature and density range considered here ) and k@xmath52 is the rate coefficient for h@xmath13 in the lw band  ( abel et  al . 1997 ; yoshida et  al . halo cores at higher central densities are departing from the approximations used in the above estimate and thus have values of the h@xmath13 fraction that are not quite the equilibrium values . halo cores which have reached core h@xmath13 fractions higher than @xmath157 have been able to cool to significantly below the virial temperature ( since the cooling time is much less than the hubble time at that value of the h@xmath13 fraction ) , allowing the baryon density to increase . the effect of cooling can also be seen in the plot of radial velocity as a function of radius , where the core regions of halos that have cooled are contracting ( albeit very slowly ) , while those which are still at the virial temperature are not showing significant signs of contraction within @xmath158 pc . all halos have some inflow at large radii , due to the infall of gas from adjoining filaments .    . panel ( a ) : number density as a function of radius . panel ( b ) : temperature as a function of radius . panel ( c ) : @xmath31 fraction as a function of radius . panel ( d ) : radial velocity as a function of radius . all quantities except enclosed gas mass are mass - weighted , and all simulations are shown at @xmath154 ( shortly before the j@xmath68 simulation collapses to high density . line types and weights correspond to simulations , as follows . black solid line : j@xmath159 . black short - dashed line : j@xmath140 . black long - dashed line : j@xmath141 . black dot short - dashed line : j@xmath142 . black dot long - dashed line : j@xmath80 . black short dashed - long dashed line : j@xmath143 . red solid line : j@xmath144 . red short - dashed line : j@xmath145 . red long - dashed line : j@xmath1 . , title="fig:",scaledwidth=45.0% ] . panel ( a ) : number density as a function of radius . panel ( b ) : temperature as a function of radius . panel ( c ) : @xmath31 fraction as a function of radius . panel ( d ) : radial velocity as a function of radius . all quantities except enclosed gas mass are mass - weighted , and all simulations are shown at @xmath154 ( shortly before the j@xmath68 simulation collapses to high density . line types and weights correspond to simulations , as follows . black solid line : j@xmath159 . black short - dashed line : j@xmath140 . black long - dashed line : j@xmath141 . black dot short - dashed line : j@xmath142 . black dot long - dashed line : j@xmath80 . black short dashed - long dashed line : j@xmath143 . red solid line : j@xmath144 . red short - dashed line : j@xmath145 . red long - dashed line : j@xmath1 . , title="fig:",scaledwidth=45.0% ] . panel ( a ) : number density as a function of radius . panel ( b ) : temperature as a function of radius . panel ( c ) : @xmath31 fraction as a function of radius . panel ( d ) : radial velocity as a function of radius . all quantities except enclosed gas mass are mass - weighted , and all simulations are shown at @xmath154 ( shortly before the j@xmath68 simulation collapses to high density . line types and weights correspond to simulations , as follows . black solid line : j@xmath159 . black short - dashed line : j@xmath140 . black long - dashed line : j@xmath141 . black dot short - dashed line : j@xmath142 . black dot long - dashed line : j@xmath80 . black short dashed - long dashed line : j@xmath143 . red solid line : j@xmath144 . red short - dashed line : j@xmath145 . red long - dashed line : j@xmath1 . , title="fig:",scaledwidth=45.0% ] . panel ( a ) : number density as a function of radius . panel ( b ) : temperature as a function of radius . panel ( c ) : @xmath31 fraction as a function of radius . panel ( d ) : radial velocity as a function of radius . all quantities except enclosed gas mass are mass - weighted , and all simulations are shown at @xmath154 ( shortly before the j@xmath68 simulation collapses to high density . line types and weights correspond to simulations , as follows . black solid line : j@xmath159 . black short - dashed line : j@xmath140 . black long - dashed line : j@xmath141 . black dot short - dashed line : j@xmath142 . black dot long - dashed line : j@xmath80 . black short dashed - long dashed line : j@xmath143 . red solid line : j@xmath144 . red short - dashed line : j@xmath145 . red long - dashed line : j@xmath1 . , title="fig:",scaledwidth=45.0% ] in this section we compare the evolution of two representative simulations . we choose the calculations which correspond to j@xmath68 and @xmath90 . these two simulations represent the extremes in halo evolution  the j@xmath68 calculation is an example of the evolution of a `` standard '' halo , as discussed in previous literature , and the j@xmath1 simulation is the most extreme example of a population iii protostar forming in a t@xmath5  k halo available from our suite of calculations . as described in section  [ sect.simsetup ] , both calculations start from the same initial conditions , with only the strength of the lw background being different . projections of the log baryon density and temperature for the j@xmath68 and @xmath90 simulations are shown in figures  [ fig.image_jlw_0 ] and  [ fig.image_jlw_1em21 ] , respectively . these projections are at the epoch of collapse , which is @xmath160  @xmath161 for the j@xmath68  @xmath162 simulation . the halo masses are significantly different : @xmath163 and @xmath164  m@xmath15 for the two calculations , reflecting the rapid pace of mergers at that redshift , and the amount of structure evolution that takes place over that relatively short time period ( @xmath165 years ) . this can clearly be seen in the left column of both figures , when the large amount of structure apparent in figure  [ fig.image_jlw_0 ] has merged into the main halo by the time of protostellar cloud formation in the j@xmath1 case , as seen in figure  [ fig.image_jlw_1em21 ] .    despite the differences in evolutionary states , the halos themselves are quite similar in figures  [ fig.image_jlw_0 ] and  [ fig.image_jlw_1em21 ] . both calculations show halos that are strongly centrally concentrated , as shown by the center column in both figures . this column shows projections of the halos which are centered on the baryon peak and are approximately one virial radius across ( and thus are scaled differently in both images ) . no fragmentation of the halo core is visible . the rightmost panel in both calculations shows the center of the halo , where the protostellar cloud is evolving . there is still no evidence for fragmentation up to a central baryon density of @xmath3  @xmath4 . figures  [ fig.radprof_jlw0.1 ] and  [ fig.radprof_jlw1em21.1 ] show the temporal evolution of several spherically - averaged baryon quantities for the j@xmath68 and @xmath90 calculations , respectively . these plots show baryon number density , enclosed baryon mass , temperature , h@xmath13 fraction , ratio of cooling time to sound crossing time , and ratio of cooling time to dynamical time , as a function of radius , evolving from a number density of @xmath166  @xmath4 to @xmath69  @xmath4 in both calculations . some qualitative commonalities are obvious between the two calculations . both halos steadily grow in central density ( panel ( a ) ) , which is coupled to the growing molecular hydrogen fraction in the central region of the halo ( panel ( d ) ) . both halos collapse quasi - statically , as shown in panels ( e ) and ( f ) . also , the temperature profiles ( as shown in panel ( c ) ) are broadly similar , in that the gas in the halo central region cools to some minimum value at r @xmath167 pc ( corresponding to n @xmath168  @xmath4 in both simulations ) and then creeps steadily upward as gas collapses to higher densities . though there are qualitative similarities between the evolution of the two calculations , detailed examination shows some significant quantitative differences . the time that the halo in the j@xmath68 simulation takes to evolve from a central density of @xmath169  @xmath4 to @xmath135  @xmath4 is @xmath170  myr , with 21 myr of that being the time required to evolve from @xmath169 to n @xmath171 @xmath4 . the time required for the j@xmath1 calculation to evolve from @xmath169  @xmath4 to @xmath135  @xmath4 is slightly shorter , @xmath60 17 myrs , with 8 myrs required to evolve to n @xmath172  @xmath4 . the temperature evolution is also somewhat different  the j@xmath68 simulation has a minimum temperature at r @xmath137 pc of 200 k , with n @xmath171@xmath173 and a molecular hydrogen fraction of @xmath174 , and roughly 1000 m@xmath15 of enclosed gas . as the gas within this radius evolves to higher densities , the temperature increases by a factor of a few , peaking at t @xmath175 k at n  @xmath176  @xmath4 , with f@xmath177 . the gas cooling time within r @xmath137 pc in the j@xmath68 simulation is always larger than the sound crossing and dynamical times by at least a factor of two , implying a quasistatic gas contraction at all times . the j@xmath1 simulation evolves somewhat differently . the collapsing gas reaches a minimum temperature of @xmath178  k at r @xmath179 pc , with a baryon density of @xmath180  @xmath4 , an enclosed gas mass of roughly @xmath86 m@xmath15 , and a molecular hydrogen fraction of @xmath181 , two orders of magnitude below that in the j@xmath68 calculation at the equivalent temperature minimum . the low h@xmath13 fraction is due to the high lw background radiation flux . as the halo evolves to higher densities , the gas temperature and h@xmath13 fraction also creep upwards . the temperature of the gas when the peak reaches roughly @xmath135  @xmath4 is @xmath182 k , with a h@xmath13 fraction of @xmath183 . the cooling time at r @xmath179 pc is an order of magnitude higher than the sound crossing or dynamical times , implying an extremely slow contraction of the gas at that radius . as the central density increases , however , the cooling time becomes lower , with t@xmath184/t@xmath185 and t@xmath184/t@xmath186 dipping to @xmath187 at its minimum , but the latter increasing toward the central density peak to a ratio of @xmath188 . this implies a more rapid contraction of gas in the halo center in the simulation with a higher uv background , despite the lack of molecular hydrogen . it is also worth noticing that the j@xmath1 calculation has h@xmath13 fractions at n @xmath168 and @xmath106  @xmath4 of @xmath189 and @xmath190 , respectively , as compared to f@xmath191 and @xmath192 at the same densities in the j@xmath68 calculation . panel ( a ) : number density as a function of radius . panel ( b ) : enclosed gas mass as a function of radius . panel ( c ) : temperature as a function of radius . panel ( d ) : @xmath31 fraction as a function of radius . panel ( e ) : ratio of cooling time to sound crossing time as a function of radius . panel ( f ) : ratio of cooling time to dynamical time as a function of radius . all quantities except enclosed gas mass are mass - weighted . line types and weights correspond to different times , as follows . red short dashed - long dashed line : @xmath193 , @xmath194 years . red dot long - dashed line : @xmath195 , @xmath196 years . red dot short - dashed line : @xmath197 , @xmath196 years . red long - dashed line : @xmath198 , @xmath199 years . red short - dashed line : @xmath200 , @xmath201 years . red solid line : @xmath202 , @xmath203 years . black short dashed - long dashed line : @xmath204 , @xmath205 years . black dot long - dashed line : @xmath206 , @xmath207 years . black dot short - dashed line : @xmath208 , @xmath209 years . black long - dashed line : @xmath210 , @xmath211 years . black short - dashed line : @xmath212 , @xmath213 years . black solid line : @xmath214 , @xmath215 years . , title="fig:",scaledwidth=30.0% ]   simulation . panel ( a ) : number density as a function of radius . panel ( b ) : enclosed gas mass as a function of radius . panel ( c ) : temperature as a function of radius . panel ( d ) : @xmath31 fraction as a function of radius . panel ( e ) : ratio of cooling time to sound crossing time as a function of radius . panel ( f ) : ratio of cooling time to dynamical time as a function of radius . all quantities except enclosed gas mass are mass - weighted . line types and weights correspond to different times , as follows . red short dashed - long dashed line : @xmath193 , @xmath194 years . red dot long - dashed line : @xmath195 , @xmath196 years . red dot short - dashed line : @xmath197 , @xmath196 years . red long - dashed line : @xmath198 , @xmath199 years . red short - dashed line : @xmath200 , @xmath201 years . red solid line : @xmath202 , @xmath203 years . black short dashed - long dashed line : @xmath204 , @xmath205 years . black dot long - dashed line : @xmath206 , @xmath207 years . black dot short - dashed line : @xmath208 , @xmath209 years . black long - dashed line : @xmath210 , @xmath211 years . black short - dashed line : @xmath212 , @xmath213 years . black solid line : @xmath214 , @xmath215 years . , title="fig:",scaledwidth=30.0% ]   simulation . panel ( a ) : number density as a function of radius . panel ( b ) : enclosed gas mass as a function of radius . panel ( c ) : temperature as a function of radius . panel ( d ) : @xmath31 fraction as a function of radius . panel ( e ) : ratio of cooling time to sound crossing time as a function of radius . panel ( f ) : ratio of cooling time to dynamical time as a function of radius . all quantities except enclosed gas mass are mass - weighted . line types and weights correspond to different times , as follows . red short dashed - long dashed line : @xmath193 , @xmath194 years . red dot long - dashed line : @xmath195 , @xmath196 years . red dot short - dashed line : @xmath197 , @xmath196 years . red long - dashed line : @xmath198 , @xmath199 years . red short - dashed line : @xmath200 , @xmath201 years . red solid line : @xmath202 , @xmath203 years . black short dashed - long dashed line : @xmath204 , @xmath205 years . black dot long - dashed line : @xmath206 , @xmath207 years . black dot short - dashed line : @xmath208 , @xmath209 years . black long - dashed line : @xmath210 , @xmath211 years . black short - dashed line : @xmath212 , @xmath213 years . black solid line : @xmath214 , @xmath215 years . , title="fig:",scaledwidth=30.0% ]   simulation . panel ( a ) : number density as a function of radius . panel ( b ) : enclosed gas mass as a function of radius . panel ( c ) : temperature as a function of radius . panel ( d ) : @xmath31 fraction as a function of radius . panel ( e ) : ratio of cooling time to sound crossing time as a function of radius . panel ( f ) : ratio of cooling time to dynamical time as a function of radius . all quantities except enclosed gas mass are mass - weighted . line types and weights correspond to different times , as follows . red short dashed - long dashed line : @xmath193 , @xmath194 years . red dot long - dashed line : @xmath195 , @xmath196 years . red dot short - dashed line : @xmath197 , @xmath196 years . red long - dashed line : @xmath198 , @xmath199 years . red short - dashed line : @xmath200 , @xmath201 years . red solid line : @xmath202 , @xmath203 years . black short dashed - long dashed line : @xmath204 , @xmath205 years . black dot long - dashed line : @xmath206 , @xmath207 years . black dot short - dashed line : @xmath208 , @xmath209 years . black long - dashed line : @xmath210 , @xmath211 years . black short - dashed line : @xmath212 , @xmath213 years . black solid line : @xmath214 , @xmath215 years . , title="fig:",scaledwidth=30.0% ]   simulation . panel ( a ) : number density as a function of radius . panel ( b ) : enclosed gas mass as a function of radius . panel ( c ) : temperature as a function of radius . panel ( d ) : @xmath31 fraction as a function of radius . panel ( e ) : ratio of cooling time to sound crossing time as a function of radius . panel ( f ) : ratio of cooling time to dynamical time as a function of radius . all quantities except enclosed gas mass are mass - weighted . line types and weights correspond to different times , as follows . red short dashed - long dashed line : @xmath193 , @xmath194 years . red dot long - dashed line : @xmath195 , @xmath196 years . red dot short - dashed line : @xmath197 , @xmath196 years . red long - dashed line : @xmath198 , @xmath199 years . red short - dashed line : @xmath200 , @xmath201 years . red solid line : @xmath202 , @xmath203 years . black short dashed - long dashed line : @xmath204 , @xmath205 years . black dot long - dashed line : @xmath206 , @xmath207 years . black dot short - dashed line : @xmath208 , @xmath209 years . black long - dashed line : @xmath210 , @xmath211 years . black short - dashed line : @xmath212 , @xmath213 years . black solid line : @xmath214 , @xmath215 years . , title="fig:",scaledwidth=30.0% ]   simulation . panel ( a ) : number density as a function of radius . panel ( b ) : enclosed gas mass as a function of radius . panel ( c ) : temperature as a function of radius . panel ( d ) : @xmath31 fraction as a function of radius . panel ( e ) : ratio of cooling time to sound crossing time as a function of radius . panel ( f ) : ratio of cooling time to dynamical time as a function of radius . all quantities except enclosed gas mass are mass - weighted . line types and weights correspond to different times , as follows . red short dashed - long dashed line : @xmath193 , @xmath194 years . red dot long - dashed line : @xmath195 , @xmath196 years . red dot short - dashed line : @xmath197 , @xmath196 years . red long - dashed line : @xmath198 , @xmath199 years . red short - dashed line : @xmath200 , @xmath201 years . red solid line : @xmath202 , @xmath203 years . black short dashed - long dashed line : @xmath204 , @xmath205 years . black dot long - dashed line : @xmath206 , @xmath207 years . black dot short - dashed line : @xmath208 , @xmath209 years . black long - dashed line : @xmath210 , @xmath211 years . black short - dashed line : @xmath212 , @xmath213 years . black solid line : @xmath214 , @xmath215 years . , title="fig:",scaledwidth=30.0% ]     simulation . panel ( a ) : number density as a function of radius . panel ( b ) : enclosed gas mass as a function of radius . panel ( c ) : temperature as a function of radius . panel ( d ) : @xmath31 fraction as a function of radius . panel ( e ) : ratio of cooling time to sound crossing time as a function of radius . panel ( f ) : ratio of cooling time to dynamical time as a function of radius . all quantities except enclosed gas mass are mass - weighted . line types and weights correspond to different times , as follows . red short dashed - long dashed line : @xmath216 , @xmath217 years . red dot long - dashed line : @xmath218 , @xmath219 years . red dot short - dashed line : @xmath220 , @xmath221 years . red long - dashed line : @xmath222 , @xmath223 years . red short - dashed line : @xmath224 , @xmath225 years . red solid line : @xmath226 , @xmath227 years . black short dashed - long dashed line : @xmath228 , @xmath229 years . black dot long - dashed line : @xmath230 , @xmath231 years . black dot short - dashed line : @xmath232 , @xmath233 years . black long - dashed line : @xmath234 , @xmath235 years . black short - dashed line : @xmath236 , @xmath237 years . black solid line : @xmath238 , @xmath239 years . , title="fig:",scaledwidth=30.0% ]   simulation . panel ( a ) : number density as a function of radius . panel ( b ) : enclosed gas mass as a function of radius . panel ( c ) : temperature as a function of radius . panel ( d ) : @xmath31 fraction as a function of radius . panel ( e ) : ratio of cooling time to sound crossing time as a function of radius . panel ( f ) : ratio of cooling time to dynamical time as a function of radius . all quantities except enclosed gas mass are mass - weighted . line types and weights correspond to different times , as follows . red short dashed - long dashed line : @xmath216 , @xmath217 years . red dot long - dashed line : @xmath218 , @xmath219 years . red dot short - dashed line : @xmath220 , @xmath221 years . red long - dashed line : @xmath222 , @xmath223 years . red short - dashed line : @xmath224 , @xmath225 years . red solid line : @xmath226 , @xmath227 years . black short dashed - long dashed line : @xmath228 , @xmath229 years . black dot long - dashed line : @xmath230 , @xmath231 years . black dot short - dashed line : @xmath232 , @xmath233 years . black long - dashed line : @xmath234 , @xmath235 years . black short - dashed line : @xmath236 , @xmath237 years . black solid line : @xmath238 , @xmath239 years . , title="fig:",scaledwidth=30.0% ]   simulation . panel ( a ) : number density as a function of radius . panel ( b ) : enclosed gas mass as a function of radius . panel ( c ) : temperature as a function of radius . panel ( d ) : @xmath31 fraction as a function of radius . panel ( e ) : ratio of cooling time to sound crossing time as a function of radius . panel ( f ) : ratio of cooling time to dynamical time as a function of radius . all quantities except enclosed gas mass are mass - weighted . line types and weights correspond to different times , as follows . red short dashed - long dashed line : @xmath216 , @xmath217 years . red dot long - dashed line : @xmath218 , @xmath219 years . red dot short - dashed line : @xmath220 , @xmath221 years . red long - dashed line : @xmath222 , @xmath223 years . red short - dashed line : @xmath224 , @xmath225 years . red solid line : @xmath226 , @xmath227 years . black short dashed - long dashed line : @xmath228 , @xmath229 years . black dot long - dashed line : @xmath230 , @xmath231 years . black dot short - dashed line : @xmath232 , @xmath233 years . black long - dashed line : @xmath234 , @xmath235 years . black short - dashed line : @xmath236 , @xmath237 years . black solid line : @xmath238 , @xmath239 years . , title="fig:",scaledwidth=30.0% ]   simulation . panel ( a ) : number density as a function of radius . panel ( b ) : enclosed gas mass as a function of radius . panel ( c ) : temperature as a function of radius . panel ( d ) : @xmath31 fraction as a function of radius . panel ( e ) : ratio of cooling time to sound crossing time as a function of radius . panel ( f ) : ratio of cooling time to dynamical time as a function of radius . all quantities except enclosed gas mass are mass - weighted . line types and weights correspond to different times , as follows . red short dashed - long dashed line : @xmath216 , @xmath217 years . red dot long - dashed line : @xmath218 , @xmath219 years . red dot short - dashed line : @xmath220 , @xmath221 years . red long - dashed line : @xmath222 , @xmath223 years . red short - dashed line : @xmath224 , @xmath225 years . red solid line : @xmath226 , @xmath227 years . black short dashed - long dashed line : @xmath228 , @xmath229 years . black dot long - dashed line : @xmath230 , @xmath231 years . black dot short - dashed line : @xmath232 , @xmath233 years . black long - dashed line : @xmath234 , @xmath235 years . black short - dashed line : @xmath236 , @xmath237 years . black solid line : @xmath238 , @xmath239 years . , title="fig:",scaledwidth=30.0% ]   simulation . panel ( a ) : number density as a function of radius . panel ( b ) : enclosed gas mass as a function of radius . panel ( c ) : temperature as a function of radius . panel ( d ) : @xmath31 fraction as a function of radius . panel ( e ) : ratio of cooling time to sound crossing time as a function of radius . panel ( f ) : ratio of cooling time to dynamical time as a function of radius . all quantities except enclosed gas mass are mass - weighted . line types and weights correspond to different times , as follows . red short dashed - long dashed line : @xmath216 , @xmath217 years . red dot long - dashed line : @xmath218 , @xmath219 years . red dot short - dashed line : @xmath220 , @xmath221 years . red long - dashed line : @xmath222 , @xmath223 years . red short - dashed line : @xmath224 , @xmath225 years . red solid line : @xmath226 , @xmath227 years . black short dashed - long dashed line : @xmath228 , @xmath229 years . black dot long - dashed line : @xmath230 , @xmath231 years . black dot short - dashed line : @xmath232 , @xmath233 years . black long - dashed line : @xmath234 , @xmath235 years . black short - dashed line : @xmath236 , @xmath237 years . black solid line : @xmath238 , @xmath239 years . , title="fig:",scaledwidth=30.0% ]   simulation . panel ( a ) : number density as a function of radius . panel ( b ) : enclosed gas mass as a function of radius . panel ( c ) : temperature as a function of radius . panel ( d ) : @xmath31 fraction as a function of radius . panel ( e ) : ratio of cooling time to sound crossing time as a function of radius . panel ( f ) : ratio of cooling time to dynamical time as a function of radius . all quantities except enclosed gas mass are mass - weighted . line types and weights correspond to different times , as follows . red short dashed - long dashed line : @xmath216 , @xmath217 years . red dot long - dashed line : @xmath218 , @xmath219 years . red dot short - dashed line : @xmath220 , @xmath221 years . red long - dashed line : @xmath222 , @xmath223 years . red short - dashed line : @xmath224 , @xmath225 years . red solid line : @xmath226 , @xmath227 years . black short dashed - long dashed line : @xmath228 , @xmath229 years . black dot long - dashed line : @xmath230 , @xmath231 years . black dot short - dashed line : @xmath232 , @xmath233 years . black long - dashed line : @xmath234 , @xmath235 years . black short - dashed line : @xmath236 , @xmath237 years . black solid line : @xmath238 , @xmath239 years . , title="fig:",scaledwidth=30.0% ] in this paper we have examined aspects of the formation of population iii stars in the presence of a soft ultraviolet background , using an adaptive mesh refinement cosmological structure formation code . given the nature of the tool and current limitations in computing power , some physics were neglected , and possible numerical issues may arise . we discuss this here . the primordial chemistry model used in these calculations ignores the effects of deuterium , lithium , and the various molecules that form between these elements and ordinary hydrogen . deuterium and lithium have been shown to be unimportant in the temperature and density regimes that we have examined in this paper ( galli & palla 1998 , 2002 ; lipovka , nez - lpez , & avila - reese 2005 ) . however , it is possible that they may be relevant in other situations of importance to population iii star formation  in particular , regions which have been ionized to very high electron fractions may experience significant cooling from the hd molecule , which due to its permanent dipole moment makes it more than 100 times more effective as a cooling agent than molecular hydrogen ( per molecule ) , and has the potential to cool gas down to approximately the temperature of the cosmic microwave background , which scales with redshift as @xmath240  k  ( flower et  al . 2000 ; galli & palla 2002 ; lipovka et  al . this gives a minimum baryon temperature of approximately @xmath241 kelvin at @xmath110 and could reduce the minimum accretion rate onto a primordial protostar significantly . lithium , while in principle an effective coolant as lih , is safely ignored , since only a tiny fraction of lithium is converted into lih  ( mizusawa , omukai , & nishi 2005 ) . a final chemical process that is omitted from our calculation is heating caused by molecular hydrogen formation at high densities ( n @xmath242  @xmath4 ) , the inclusion of which may result in differences in the temperature evolution of the gas in the highest density regimes considered here . this does not significantly affect the conclusions of our paper , as the significant differences we observe between simulations occur at densities much lower than @xmath36  @xmath4 . self - shielding of the photodissociating background by molecular hydrogen in the high - density gas is ignored in these calculations . though this effect could in principle be important , the actual column densities of molecular hydrogen are typically far too small to actually block the soft uv flux . according to ( 2001 ) ( and references therein ) , a column density of @xmath243  @xmath11 is enough for shielding to become important . however , this was derived for a static distribution of h@xmath13 , while the lw band consists of hundreds of individual lines whose width in this case is dominated by doppler broadening . it is useful to note in this case that the average line width is @xmath179  km / s and the rms baryon velocity in our calculations are @xmath244  km / s . in order for self - shielding to be important in the case of a turbulent medium , the column density must be much higher . typical maximum h@xmath13 column densities in our calculations are on the order of @xmath245  @xmath11 , but these occur late in the collapse of the halo core , and in the highest density regions the cooling and h@xmath13 production times are much shorter than the photodissociation time scale , at which point self - shielding becomes unimportant . it is worth noting that there are some regimes where self - shileding can be critical . ( 2007 ) finds that self - shielding can strongly affect the evolution of a collapsing primordial halo which is being illuminated with photodissociating flux by a neighboring star , though the situation is somewhat idealized . additionally , we do not consider the more complicated effects relating to population iii stars which form in halos in the cosmic neighborhoods where previous generations of stars have existed . this allows us to ignore complex radiative , chemical and dynamical effects that would vastly complicate our calculations  ( oshea et  al . 2005a ; johnson , greif , & bromm 2006 ; yoshida et  al . 2006a ; greif et  al . 2007 ; abel , wise , & bryan 2007 ; ahn & shapiro 2007 ) .    a further effect that is ignored in this paper is @xmath29 photodetachment . this could in principle be a significant effect  @xmath29 is the catalyst for molecular hydrogen formation in the dominant 2 formation channel at densities @xmath147  @xmath4 , and it can be detached by photons with energies @xmath246  ev . photons with this energy would be produced in great numbers by the same massive population iii stars that we assume are producing the molecular hydrogen photodissociating background . however , as shown by ( 2001 ) , the rate of photodetachment is orders of magnitude smaller than the rate of @xmath29 formation at densities comparable to that found in the centers of the halos examined in this work , and thus the process of @xmath29 photodetachment can be safely ignored .    in this paper , we examine the effects of molecular hydrogen dissociating backgrounds which range in strength from j@xmath247 to @xmath9  ergs@xmath10  @xmath11  hz@xmath10  sr@xmath10 in the lw band ( @xmath248 ev ) . this range is in good agreement with the expected range of photodissociating backgrounds predicted by   & abel ( 2005 ) , but it is still finite . examination of our simulations show that the case with the lowest uv background ( j@xmath140 ) is almost identical to the `` control '' ( j@xmath68 ) calculation , justifying our choice of minimum value . our upper value is consistent with wise & abel , and examination of figures  [ fig.halovals1 ] and [ fig.radprof_zcoll.1][fig.radprof_zcoll.3 ] suggests that increases in the strength of the uv background ( within reasonable values ) will result in a further delay in halo collapse , but no major change in the mode of star formation observed . this is an effect of the halo properties  a small amount of molecular hydrogen will always exist in these halos , allowing the gas at the center to cool and contract quasi - statically to higher densities . once the gas has collapsed to very high densities , rapid h@xmath13 formation via the 3-body process will occur , essentially independent of the strength of the soft uv background , and the gas will be able to cool down to @xmath249 k very quickly . the effects of magnetic fields are completely ignored in the simulations discussed in this work . magnetic fields are discussed in detail in paper i , but , to summarize , a fairly high seed magnetic field is needed to be dynamically significant at the relatively low densities we explore in this work . the possible importance of magnetic fields has been explored in analytic and semi - analytic work ( tan & mckee 2004 ; silk & langer 2006 ; maki & susa 2007 ) . we will examine the possible evolution of magnetic fields within the context of cosmological amr simulations in a later paper ( oshea & turk 2007 , in preparation ) . the simulations presented in this paper are generated assuming a cosmology that is somewhat different than the currently - favored wmap year iii `` best - fit '' model  ( spergel et  al . most importantly , our value of @xmath250 is 0.9 , which is significantly higher than the wmap value of 0.761 . the general effect of a higher @xmath250 is to cause structure formation to take place earlier , and thus at a given redshift one would expect significantly more halos in our chosen cosmological model than in the wmap year iii model . however , the evolution of any single cosmological halo , such as the one examined in this work , is not particularly affected by this parameter in the sense that we are not examining halo stastical properties . in addition to @xmath250 , the ratio of @xmath251 in our simulations is 0.1337 , while it is @xmath252 in the wmap year iii cosmology . this may result in some small quantitative differences in , e.g. , the redshift of halo collapse , but should not significantly affect our results . we direct the reader to paper i for a detailed discussion of other possible numerical issues , such as the generation of cosmological initial conditions , the assumption that the cold dark matter model is correct , the choice of halo in our simulations , and the size of the simulation volumes used . this paper explores the formation of population iii stars in simulations with a constant soft uv background . our results agree well qualitatively with that of ( 2001 ) ; we both find that a soft uv background can delay the formation of molecular hydogen , and thus the cooling and collapse , of small ( @xmath253  m@xmath15 ) cosmological halos in which population iii stars form . we also find that increasing the soft uv background increases the minimum halo mass required for a halo to collapse ( in a way similar to that of smoothing the dark matter power spectrum at small scales  see & norman ( 2006 ) for a discussion ) . machacek et al . derived a mass threshold for collapse as a function of the lw background flux that agrees well with our simulations , though the halo masses in our calculation are significantly higher . this is due to the halo that we examine being a `` typical '' halo rather than at the threshold mass for star formation . presumably , if we performed these calculations using many halos in a range of cosmological realizations , we would find a minimum halo mass that is somewhat lower than the masses seen in our calculations . our work is a significant improvement upon that of machacek et al . in some respects , as our simulations are much more highly resolved and we examine the evolution of a single halo over a much wider range of soft uv background fluxes . one drawback of our work compared to machacek et al . is that we examine the evolution of a single halo , albeit with a broad range of uv backgrounds . our work is similar to the calculations in   ( 2003 ) which examine population iii star formation in the presence of a photodissociating ultraviolet background . they find that for values of j@xmath52 above @xmath254  erg s@xmath10  @xmath11  hz@xmath10  sr@xmath10 , hydrogen molecules are rapidly dissociated and gas cooling is inefficient , implying that the collapse of gas in the halo center may be delayed until atomic line cooling dominates , in halos with t@xmath12  k. while we see delays in halo collapse , we do see that a tiny amount of 2 can form even at high fuv background strengths , eventually allowing gas to cool and collapse to high densities without ever reaching temperatures at which atomic line cooling would be effective ( t @xmath255  k ) . yoshida et al . make this statement based on a set of simulations using two different values of j@xmath52 ( @xmath254 and @xmath256 ) , which bracket the `` break '' in properties seen in our results . the simulations used to obtain this result had sph particles with m@xmath257  h@xmath10  m@xmath15 and thus an effective mass resolution comparable to the mass of gas within the quasistatically - collapsing central region of a halo . our mass resolution is more than two orders of magnitude higher in the center of the cosmological halos that we examine . the differences between our result and theirs is likely due to the significant differences in resolution and may in fact be consistent once resolution is taken into account . we see a strong relationship between accretion rates onto the protostellar cloud and the strength of the photodissociating background . calculations with higher ultraviolet background strengths typically have a larger spherically - averaged baryon accretion rate onto the evolving protostellar cloud . this is clearly due to variation in halo central temperature relating to the amount of molecular hydrogen existing in the center of the halo at the epoch of collapse . accretion onto the evolving protostellar cloud is subsonic , and thus regulated by the local sound speed ( which scales as t@xmath258 ) . this implies some relationship between the strength of the photodissociating background and the final stellar mass of the primordial star  however , the details of this relationship are unclear , and depends on many factors . for example ,   & palla ( 2003 ) suggest that an increase in accretion rate above @xmath153  m@xmath15/yr will actually result in a _ decrease _ in the final mass of the star due to radiative feedback from the evolving protostar . however , their results use one - dimenstional , fairly idealized models , and geometrical effects may be important . this is explored in   & mckee ( 2004 ) , who use a combination of analytic and semi - analytic models of the evolving system . they suggest that accretion onto the protostar is highly non - spherical , and is in fact mediated by an accretion disk . they argue that a larger accretion rate onto the disk will lead to a larger star , though with some limitations . finally , ( 2001 ) and & yoshii ( 2003 ) study the imf of stars that form in t@xmath12  k metal - free protogalaxies in the presence of an h@xmath13 photodissociating background . they conclude that photodissociation actually decreases the jeans mass of the gas at high densities , and thus reduces the fragmentation mass scale of the clouds and presumably the stellar mass . they use one - dimensional , spherically symmetric simulations , however . we can explore this issue in more detail in later calculations , when appropriate models for shielding , cooling via @xmath29 , and more advanced chemistry have been implemented into enzo . this will be examined in a later paper . it is worth noting that the range of accretion rates observed in this situation are within the range of rates seen in the @xmath17cdm simulations discussed in paper i. this implies that the mass ranges of the resulting stars will not be significantly different than in the `` standard '' population iii star formation scenario . we see no fragmentation of the high - density baryon core , up to a baryon number density of @xmath259  @xmath4 , for the entire range of simulations explored in this paper . this is due to a combination of effects , but predominantly the poor cooling properties of molecular hydrogen . the gas is relatively hot ( @xmath260 kelvin ) and thus has a high sound speed , which helps to damp out perturbations in the halo center which would otherwise result in multiple fragments . this effect is exacerbated in simulations with high ultraviolet background strengths , as the temperatures ( and thus the sound speeds ) are generally higher . this result implies that there is no fundamental change in the mode of population iii star formation as halos grow , and that more massive `` proto - galactic '' halos , with t@xmath261  k , will continue to form a single massive star per halo . these halos will have a much larger binding energy than the smaller , m @xmath14  m@xmath15 halos which have traditionally been examined by numerical simulations , and implies that multiple generations of star formation may be able to take place in a single halo . this leads to the possibility of `` self - enrichment , '' where a single population iii star enriches the high - density center of the halo to metallicities high enough to change the cooling properties of the gas . this could cause a strong change in the imf ( santoro & shull 2006 ; smith & sigurdsson 2007 ) , and these objects could make a significant contribution to the reionization history of the universe  ( mackey , bromm , & hernquist 2003 ; greif & bromm 2006 ) . we will examine whether the high - density gas in the centers of these more massive halos fragments in a later work . the values of accretion rates onto the evolving protostellar cloud which are observed with increasing photodissociating background flux ( which are always comparable to those observed in paper i ) , in addition to the apparent lack of fragmentation of gas in the halo center up to densities of n@xmath262  @xmath4 , implies that there is no major change in the mass range of population iii stars as halo masses increase to t@xmath5  k. this is in mild disagreement with previous semianalytic results by & haiman ( 2002 ) and ( 2003 ) . this is apparently because previous authors assumed that 2 cooling in `` second generation '' ( t@xmath12  k ) objects is unimportant , while we observe that a small 2 fraction , which forms in equilibrium with the uv background , can cool gas effectively to temperatures of @xmath263  k. we do agree with these authors that the majority of primordial star formation will take place in objects with t@xmath5  k , the so - called `` second generation '' objects  ( macintyre , santoro , & thomas 2006 )    our work appears to contradict the results of & bromm ( 2006 ) , who use analytic and semianalytic techniques to follow the collapse of gas within halos of t@xmath264 . they find that , in their model , gas cools to @xmath265  k via atomic line cooling , and then contracts nearly isothermally to high densities , allowing a molecular hydrogen fraction of f@xmath266 to build up independent of local density and temperature . the gas in the center of the halo then cools and allows the gas to fragment on scales of @xmath158  m@xmath15 . we do not see this mode of star formation , possibly because we never generate the significantly ionized halos with t@xmath267  k upon which their scenario depends . rather , we simulate halos that in the cases of j@xmath0 are built up to t@xmath5  k via mergers , and the gas is never ionized to a significant degree . one possible reason for the observed differences is that we ignore the formation of deuterium hydride and its effects on the cooling properties of primordial gas . it is possible that , in some contexts , the inclusion of hd may result in enhanced fragmentation . this will be examined in future studies . the scenario described by oh & haiman may still in fact be able to occur , and deserves further detailed numerical study to determine the fate of the collapsing gas and the implications of the ( presumably ) extremely high accretion rates onto the forming protostellar clouds . our results are similar to those shown by & abel ( 2007 ) , who perform similar simulations examining the collapse of gas in cosmological halos in the presence of a photodissociating background . they also find that , regardless of the strength of the uv background , the collapse of gas in the centers of the cosmological halos in question still occurs via h@xmath13 cooling . this holds even for quite extreme examples , such as when the original electron fraction is set to unphysically low levels ( which should suppress h@xmath13 formation , at least temporarily ) . in this paper we have performed a suite of high dynamical range ( l@xmath268 ) 3d adaptive mesh cosmological simulations of the formation of population iii stars in a @xmath17cdm universe in the presence of a molecular hydrogen photodissocating ( `` lyman - werner '' ) ultraviolet background . the purpose of these calculations is to understand the effect that the soft ultraviolet background has on the evolution of the gas in collapsing cosmological halos and to determine possible effects on the forming protostellar cloud . these simulations are all of a single cosmological realization , but with varied ultraviolet background strengths . our principal results are as follows :    \1 . our calculations show that , as the flux of ultraviolet radiation in the lyman - werner band is increased , population iii star formation in a given cosmological halo is delayed to later times and , as a result , an increase in the virial mass of this halo at the onset of baryon collapse . this is in good agreement with previous work by   ( 2001 ) and   ( 2003 ) . we find that , contrary to previous work , the formation of primordial stars is never completely suppressed , regardless of the strength of the uv background . a small amount of molecular hydrogen always exists in cosmological halos , and allows cooling and halo collapse in gas which is bathed in a strong photodissociating background but has not reached @xmath86  k. the previously - suggested mode of star formation in `` second generation '' halos , where the collapse of gas to high density is completely suppressed until cooling can take place via atomic hydrogen line transitions , is never observed in our calculations , which have photodissociating background strengths up to j@xmath1 .    \3 . though the molecular hydrogen fraction in the centers of halos which are bathed in strong ( j@xmath0 ) lyman - werner radiation is very small at early times ( f@xmath269 ) , the gas at the center of the halo can still collapse due to the efficient cooling of molecular hydrogen at @xmath270  k , and to the extended time these halos require to collapse compared to halos in the presence of much smaller uv backgrounds . we observe that the estimated accretion rate onto the forming protostellar cloud varies strongly as a function of j@xmath52 , with simulations that have a stronger ultraviolet background having higher accretion rates . this is a function of the suppression of molecular hydrogen formation ( and thus suppression of effective cooling ) by the photodissociating background and of the higher virial temperatures of these halos when the gas finally collapses . only a single clump is formed at the center of each collapsing cosmological halo , regardless of the strength of the photodissociating background , up to a baryon number density of @xmath3  @xmath4 . this implies that , as in the more commonly - studied mode of star formation , we will find only a single star per halo even in objects which are massive enough that t@xmath12  k.
Abstract:

Answer: we examine aspects of primordial star formation in the presence of a molecular hydrogen - dissociating ultraviolet background . we compare a set of amr hydrodynamic cosmological simulations using a single cosmological realization but with a range of ultraviolet background strengths in the lyman - werner band . this allows us to study the effects of lyman - werner radiation on suppressing 2 cooling at low densities as well as the high - density evolution of the collapsing cloud core in a self - consistent cosmological framework . we find that the addition of a photodissociating background results in a delay of the collapse of high density gas at the center of the most massive halo in the simulation and , as a result , an increase in the virial mass of this halo at the onset of baryon collapse . we find that , contrary to previous results , population iii star formation is not suppressed for j@xmath0 , but occurs even with backgrounds as high as j@xmath1 . we find that 2 cooling leads to collapse despite the depressed core molecular hydrogen fractions due to the elevated 2 cooling rates at @xmath2 k. we observe a relationship between the strength of the photodissociating background and the rate of accretion onto the evolving protostellar cloud core , with higher lw background fluxes resulting in higher accretion rates . finally , we find that the collapsing cloud cores in our simulations do not fragment at densities below @xmath3  @xmath4 regardless of the strength of the lw background , suggesting that population iii stars forming in halos with t@xmath5  k may still form in isolation .