The First Stars: formation under X-ray feedback
We investigate the impact of a cosmic X-ray background (CXB) on Population III stars forming in a minihalo at . Using the smoothed particle hydrodynamics code gadget-2, we attain sufficient numerical resolution to follow gas collapsing into the centre of the minihalo from cosmological initial conditions up to densities of , at which point we form sink particles. This allows us to study how the presence of a CXB affects the formation of and in the gas prior to becoming fully molecular. Using a suite of simulations for a range of possible CXB models, we follow each simulation for 5000 after the first sink particle forms. The CXB provides two competing effects, with X-rays both heating the gas and increasing the free electron fraction, allowing more H to form. X-ray heating dominates below , while the additional H cooling becomes more important above . The gas becomes optically thick to X-rays as it exits the quasi-hydrostatic ‘loitering phase,’ such that the primary impact of the CXB is to cool the gas at intermediate densities, resulting in an earlier onset of baryonic collapse into the dark matter halo. At the highest densities, self-shielding results in similar thermodynamic behaviour across a wide range of CXB strengths. Consequently, we find that star formation is relatively insensitive to the presence of a CXB; both the number and the characteristic mass of the stars formed remains quite similar even as the strength of the CXB varies by several orders of magnitude.
keywords:stars: formation — stars: Population III — cosmology: theory — early Universe — dark ages, first stars
The formation of the first stars in the Universe marked a watershed moment in cosmic history. It was during this as-yet unobserved epoch that our Universe began its transformation from the relatively simple initial conditions of the Big Bang to the complex tapestry of dark matter (DM), galaxies, stars and planets that we see today (Barkana & Loeb, 2001; Miralda-Escudé, 2003; Bromm et al., 2009; Loeb, 2010). The radiation from these so-called Population III (Pop III) stars swept through the Universe, beginning the process of reionisation (Kitayama et al., 2004; Sokasian et al., 2004; Whalen et al., 2004; Alvarez et al., 2006; Johnson et al., 2007; Robertson et al., 2010), while the heavy elements forged in their cores and released in the violent supernova explosions marking their deaths began the process of chemical enrichment (Madau et al. 2001; Mori et al. 2002; Bromm et al. 2003; Heger et al. 2003; Umeda & Nomoto 2003; Tornatore et al. 2007; Greif et al. 2007, 2010; Wise & Abel 2008; Maio et al. 2011; recently reviewed in Karlsson et al. 2013). These effects are strongly dependent on the characteristic mass of Pop III stars, which determines their total luminosity and ionising radiation output (Schaerer, 2002), and the details of their eventual demise (Heger et al., 2003; Heger & Woosley, 2010; Maeder & Meynet, 2012). As a result, developing a thorough knowledge of how environmental effects influence the properties of these stars is crucial to understanding their impact on the intergalactic medium (IGM) and subsequent stellar generations.
While the complexities of the various physical processes involved have so far prevented a definitive answer to this question, the basic properties of the first stars have been fairly well established, with the consensus that they formed in dark matter ‘minihaloes,’ having on the order of at (Couchman & Rees, 1986; Haiman et al., 1996; Tegmark et al., 1997). Pioneering numerical studies of the collapse of metal-free gas into these haloes, where molecular hydrogen was the only available coolant, suggested that Pop III stars were very massive—on the order of 100—due to the lack of coolants more efficient than (e.g., Bromm et al., 1999, 2002; Abel et al., 2002; Yoshida et al., 2003; Bromm & Larson, 2004; Yoshida et al., 2006; O’Shea & Norman, 2007). More recent simulations, benefiting from increased resolution, have found that significant fragmentation of the protostellar disc occurs during the star formation process, with the protostellar cores ranging in mass from 0.1 to tens of solar masses (Stacy et al., 2010; Clark et al., 2011a, b; Greif et al., 2011a, 2012; Stacy & Bromm, 2013; Susa et al., 2014; Hirano et al., 2014, 2015), with a presumably flat initial mass function (IMF; Dopcke et al. 2013).
One intriguing outcome of these studies is that while the protostellar disc does indeed fragment, it only marginally satisfies the Gammie (2001) criterion for disc instability (Clark et al., 2011b; Greif et al., 2011a, 2012). While protostellar feedback is insufficient to stabilise the disc (Smith et al., 2011; Stacy et al., 2012a), it is possible that an external heating source could serve to stabilise the disc and prevent fragmentation. One promising source of such an external background is far-ultraviolet radiation in the Lyman-Werner (LW) bands (11.2-13.6 ). While lacking sufficient energy to interact with atomic hydrogen, LW photons efficiently dissociate molecules, which serve as the primary coolant in primordial gas. This diminishes the ability of the gas to cool, but studies have found that the critical LW flux required to suppress cooling is far above the expected mean value of such radiation (Dijkstra et al., 2008). Another possible heating source was recently explored by Smith et al. (2012; see also Ripamonti et al. 2009, 2010) who investigated the ability of DM annihilation to suppress fragmentation of the protostellar disc. While such heating is unable to suppress star formation, it does serve to stabilise the disc, suppressing fragmentation within 1000 AU of the central protostar, at least for a while (Stacy et al., 2012b, 2014).
While LW radiation alone is unable to reliably suppress cooling in minihaloes, significant sources of LW photons, i.e., active star forming regions, contain large numbers of massive stars, and possibly mini-quasars as well (Kuhlen & Madau, 2005; Jeon et al., 2012, 2014). Not only do massive stars end their lives as supernovae, leaving behind remnants that are significant sources of X-rays, a significant fraction of these stars are likely to be in tight binaries (e.g., Clark et al., 2011b; Greif et al., 2012; Mirocha, 2014) and produce high-mass X-ray binaries (HMXBs). As the cross-section of neutral hydrogen for X-rays is relatively small, such photons easily escape their host haloes, building up a cosmic X-ray background (CXB; Oh 2001; Haiman et al. 2000; Venkatesan et al. 2001; Glover & Brand 2003; Cen 2003; Kuhlen & Madau 2005; Jeon et al. 2012, 2014). This CXB serves to both heat and increase the ionisation fraction of gas in neighbouring minihaloes, which in turn can serve to increase the fraction of the gas by increasing the number of free electrons available to act as catalysts.
In this paper we consider the effects of such an X-ray background on Pop III star formation. In Section 2 we provide the cosmological context for this study, estimating both the expected intensity of the CXB, and the amount of additional heating required to prevent minihalo collapse. Our numerical methodology is described in Section 3, while our results are found in Section 4. Finally, our conclusions are gathered in Section 5. Throughout this paper we adopt a CDM model of hierarchical structure formation, using the following cosmological parameters: , , , and .
2 Cosmological Context
2.1 Early Cosmic X-ray Background
The predominant source of X-rays at high redshifts was likely HMXBs. Supermassive black holes were not yet common during this era, and while supernovae produce significant X-ray radiation, their transient nature precludes them from efficiently building up an X-ray background. We can calculate the energy density of X-rays produced by HMXBs as follows:
where is the mass fraction of stars that form HMXBs, is the comoving star formation rate density (SFRD) as a function of redshift, and , and are the X-ray luminosity, mass and lifespan of a typical HMXB, respectively. The X-ray background accumulates over the Hubble time , and the factor of accounts for the conversion from a comoving SFRD to a physical energy density.
We employ the SFRD calculated by Greif & Bromm (2006), but see Campisi et al. (2011) for a more recent calculation in the context of Pop III gamma-ray burst observations. The Greif & Bromm (2006) SFRD incorporates both Pop III and Population I/II star formation with self-consistent reionisation and chemical enrichment. Their SFRD history only extends out to ; we extrapolate back to using a simple log-linear fit.
The mass fraction of stars forming HMXBs, , is determined by the mass fraction of stars forming black holes , the fraction of black holes in binary systems and the fraction of binaries close enough for mass transfer to occur . As their IMF is nearly flat with a characteristic mass of a few (Bromm, 2013), we make the plausible assumption that half of all Pop III stars end up forming black holes. Recent work by Stacy & Bromm (2013) found that just over half of Pop III stars end up in binary pairs; we set accordingly. While the orbital distribution of nearby solar-type stars is well-studied (e.g., Duquennoy & Mayor, 1991), that of Pop III stars is still very uncertain (but see Stacy & Bromm, 2013). We therefore assume that the orbital distribution of Pop III stars is flat, i.e, is constant, with a minimum orbital radius of and a maximum of , which is approximately the size of the self-gravitating baryonic core of a minihalo. For mass transfer onto the BH to occur, the companion star must at some point exceed its Roche limit; hence we designate all binaries with orbital distances less than as ‘close.’ Doing so, we find that approximately of Pop III binaries will be close enough for mass transfer to occur. Thus, for Pop III stars, we can estimate as follows:
Assuming HMXBs accrete material at close to the Eddington limit, we can approximate
where is the bolometric luminosity. Following Jeon et al. (2014), we assume that the spectral distribution of the emerging radiation takes the form of a thermal multi-colour disc at frequencies below , and that of a non-thermal power law at higher frequencies (e.g., Mitsuda et al., 1984). Assuming the entire accretion luminosity is emitted between and , the fraction of the total luminosity emitted between 1 and 10 is approximately 30%. We choose 1 for the lower limit as the cross section of neutral hydrogen increases rapidly at lower frequencies. X-rays below thus cannot efficiently contribute to a pervasive background. We therefore approximately set
where and . Following Inayoshi & Omukai (2011), we employ
where is the specific X-ray intensity, , and is the normalisation factor. In addition to this fiducial estimate for , henceforth referred to as model , we consider three additional models with ten, one hundred and one thousand times the intensity of , as shown in Figure 1.
2.2 Jeans Mass Filtering under X-ray Feedback
The presence of an ionising X-ray background will necessarily heat the gas in the IGM. As a result, collapse into a minihalo will be suppressed when the thermal energy of the gas exceeds the baryonic gravitational potential energy of the minihalo in question. The baryonic gravitational potential energy of a minihalo prior to collapse is given by
where is the gravitational constant, is the virial mass, and is the virial radius of the minihalo. The gas in the minihalo will naturally undergo adiabatic heating as it is compressed, but the cooling provided by the formation of molecular hydrogen keeps this process from halting collapse. An additional heating mechanism, provided here by X-rays, is required to prevent dissipational collapse of the gas. We approximate this excess thermal energy as
where is the X-ray heating rate per unit volume and is the freefall time on which gravity draws gas into the minihalo, given by
where is the background density. To first order then, the critical heating rate for suppressing collapse is given by
Utilising the fact that
where is the average matter density at the present epoch, we can solve for as a function of halo mass and redshift:
where we have normalised to typical minihalo values.
Prior to minihalo virialisation, gas densities are relatively low, such that attenuation of the CXB as it penetrates the minihalo is negligible (see Figure 4). Given this, is related to the intensity of the CXB as follows:
where is the gas number density and is the contribution to the heating rate from secondary ionisation. While the contribution of secondary ionisation events to the heating rate has a complicated dependence on the ionisation fraction (Shull & van Steenberg, 1985, see Section 3.4 for details), never enhances by more than a factor of two in the regime considered here. When secondary ionisation heating is negligible, we can similarly estimate the critical X-ray background required to suppress collapse:
When secondary ionisation heating is maximally effective then, will be a factor of two lower. Both limits for are shown in Figure 1 for both a and a minihalo. These can be interpreted as the approximate range in which the CXB will begin to have a significant impact on the collapse of gas into such a minihalo.
3 Numerical Methodology
3.1 Initial Setup
We use the well-tested -body smoothed particle hydrodynamics (SPH) code gadget-2 (Springel, 2005). We initialised our simulations in a periodic box of length 140 kpc (comoving) at in accordance with a CDM model of hierarchical structure formation. An artificially enhanced normalisation of the power spectrum, , was used to accelerate structure formation. See Stacy et al. (2010) for a discussion of the validity of this choice. High resolution in this simulation was achieved using a standard hierarchical zoom-in technique for both DM and SPH particles. Three nested levels of additional refinement at 40, 30 and 20 kpc (comoving) were added, each centred on the point where the first minihalo forms in the simulation. As resolution increases, each ‘parent’ particle is replaced by eight ‘child’ particles, such that at the greatest refinement level, each original particle has been replaced by 512 high-resolution particles. These highest-resolution SPH particles have a mass , such that the mass resolution of the simulation is , where is the number of particles used in the SPH smoothing kernel (Bate & Burkert, 1997).
3.2 Thermodynamics and Chemistry
Our chemistry and cooling network is the same as that described in Greif et al. (2009). We follow the abundances of , , , , , , , , , , and e. All relevant cooling mechanisms, including cooling via and collisional excitation and ionisation, recombination, bremsstrahlung and inverse Compton scattering are accounted for. Of particular importance is cooling via the ro-vibrational modes of , which are excited by collisions with and atoms and other molecules. At high densities, additional processes must also be included in order to properly model the gas evolution. For example, three-body formation and the associated heating becomes important above (Turk et al., 2011). The formation rates for these reactions are uncertain; we employ the intermediate rate from Palla et al. (1983). At densities greater than , the ro-vibrational lines for become optically thick, decreasing the efficiency of such cooling. We employ the Sobolev approximation and an escape probability formalism to account for this (see Yoshida et al. 2006; Greif et al. 2011a for details).
3.3 Optical Depth Estimation
Over the length scale of our cosmological box, the X-ray optical depth of primordial gas is everywhere except approaching the centre of the star-forming minihalo. As the CXB radiation penetrates the minihalo, it will necessarily be attenuated due to the high column density of the intervening gas. In order to estimate the optical depth, we directly calculate the column density approaching the centre of the accretion disc along both polar and equatorial lines of sight in the absence of a CXB. As shown in Figure 2, the column density remains essentially constant for the duration of the simulation, and has a simple power-law dependence on radius over several orders of magnitude up to the resolution limit of our simulation (). This same power-law behaviour can be seen versus density as well, as shown in Figure 3. In addition we notice that for a given gas density, the column density along the pole is roughly a factor of 10 lower than along the equator. Performing an ordinary least squares fit to the combined data from several snapshots for , we find that the column density along the pole and equator is well fit by
respectively, and does not change appreciably with increasing .
To account for this difference in column density, we assume every line of sight within 45 degrees of the pole (Hosokawa et al., 2011) experiences column density and every other line of sight experiences . Radiation within the opening angle is attenuated by while the remainder of the background is attenuated by , allowing us to define an effective optical depth such that
Accounting for this attenuation and incorporating the cosmic X-ray background described in Section 2.1, the resulting ionisation and heating rates for the case at as a function of total number density are shown in Figure 4.
3.4 X-ray Ionisation and Heating
To study the effects of X-ray ionisation and heating on primordial star formation, we implement a uniform CXB as discussed in Section 2.1. For further details, see Jeon et al. (2012, 2014). Accounting for attenuation of the incident X-ray radiation while penetrating the minihalo, the primary ionisation rate coefficient for chemical species can be written as
where and .
We include the effects of secondary ionisation from energetic electrons released by the absorption of X-ray photons by adopting the fitting formulae of Shull & van Steenberg (1985; see also Valdés & Ferrara 2008; Furlanetto & Stoever 2010), who calculated the fractions of the initial electron energy going into heating the surrounding gas, as well as into secondary ionisations of and ( and , respectively). While such secondary ionisation events have a significant impact on the ionisation fraction of and , secondary ionisations of are negligible (Shull & van Steenberg, 1985), and are not included here. The effective ionisation rates are thus given by
Here is the heating rate at which excess energy from the initial X-ray photoionisation is released into the gas, given by
where is the ionisation threshold of the species in question; , 24.6 and 54.4 for hydrogen, neutral helium and singly ionised helium, respectively.
The fraction of the initial electron energy going into secondary ionisations depends on the hydrogen ionisation fraction as follows:
Thus the total heating rate , including contributions from both primary and secondary ionizations, can be written as
3.5 Sink Particles
We employ the sink particle method described in Stacy et al. (2010). When the density of a gas particle exceeds , we replace it and all non-rotationally-supported particles within the accretion radius with a single sink particle. We set equal to the resolution length of the simulation: . Here,
where . The sink thus immediately accretes the majority of the particles within its smoothing kernel, such that its mass is initially close to . Once the sink is formed, additional gas particles and smaller sinks are accreted as they approach within of that sink particle. After each accretion event, the position and momentum of the sink particle is set to the mass-weighted average of the sink and the accreted particle.
Following the creation of a sink particle, its density, temperature and chemical abundances are no longer updated. The sink’s density is held constant at , and its temperature is kept at 650, typical for collapsing gas reaching this density; the pressure of the sink is set correspondingly. Assigning a temperature and pressure to the sink particle in this fashion allows it to behave as an SPH particle. This avoids the creation of an artificial pressure vacuum, which would inflate the accretion rate onto the sink (see Bromm et al., 2002; Martel et al., 2006). The sink’s position and momentum continue to evolve through gravitational and, initially, hydrodynamical interactions with the surrounding particles. As it gains mass and gravity becomes the dominant force, the sink behaves less like an SPH particle and more like a non-gaseous -body particle.
We perform a total of nine simulations, following each for 5000 years after the formation of the first sink particle in the simulation. Beyond our fiducial case of and the standard case of , we examine three additional simulations with , , and times . The results of a simulation were indistinguishable from . Not only is the CXB at high redshifts subject to fairly large uncertainties, but as Pop III star formation is highly biased there is likely a significant amount of cosmic variance. In addition, we also consider the optically-thin limit where . While physically unrealistic, these numerical experiments serve as comparison cases to more fully elucidate the physics in the properly shielded simulations.
4.1 Initial Collapse
Collapse under Enhanced Cooling
Figure 5 shows the final simulation output on various scales for both the and the shielded case. While the expected filamentary structure is visible in all cases, the effects of X-ray heating in the case are readily apparent, with the low-density filamentary gas experiencing significant heating. With the usual definition of the virial radius as , where , and are the average halo density, the density at the point of virialisation and the background density at the time of halo virialisation, we find that our minihalo collapses at with and , typical for the minihalo environment (Bromm, 2013).
After the minihalo virialises, the gas continues to collapse in accordance with the standard picture of Pop III star formation (e.g., Stacy et al., 2010; Greif et al., 2012; Stacy & Bromm, 2013). This picture is modified slightly in the presence of an X-ray background, but in the case the gas heats adiabatically until reaching , attaining temperatures of 1000. Beyond this point, the gas is able to cool via the ro-vibrational modes of , reaching a minimum temperature of 200 at a density of . After exiting the quasi-hydrostatic ‘loitering phase’ (Bromm et al., 2002), the gas enters runaway collapse until , when three-body reactions become important. This process turns the gas fully molecular by , the density at which we form sink particles. As seen in Figure 6, this basic picture holds true even as we vary the strength of the CXB by several orders of magnitude.
We find that X-ray heating dominates at low densities in all cases, while above the additional cooling catalysed by X-ray ionisation exceeds it; the precise density at which this transition occurs depends on the strength of the CXB. This enhanced cooling significantly impacts the subsequent evolution of the gas as it collapses to high densities; the minimum temperature of the gas in the loitering phase approaches the CMB floor as increases, and—in the unshielded simulations—remains cooler than the case in all later stages of the collapse. This allows the gas to more easily fulfil the Jeans criterion and thus collapse sooner, as demonstrated by Figure 7.
When shielding is properly accounted for, the gas becomes optically thick to X-rays as it exits the loitering phase and proceeds to runaway collapse (see Figure 4). In the absence of continued X-ray ionisation, the free electron fraction decays, re-converging with that of the case. Consequently, the thermodynamic behaviour of the gas at high densities in the shielded simulations is remarkably similar, even as we vary the CXB strength by several orders of magnitude as shown in Figure 6. This convergence under a wide range of X-ray backgrounds is similar to the behaviour noted by Stacy et al. (2011) and Greif et al. (2011b) in the presence of dark matter–baryon streaming, though we observe an earlier—rather than a later—collapse.
Minihalo collapse in the case is completely suppressed. While the gas in the very centre of the minihalo does initially begin to cool via , this process is quickly overwhelmed by heating from the increasingly strong X-ray background. Despite reaching densities and temperatures approaching and this cold core is eventually dissipated as the CXB continues to heat the gas. The reason for this suppression is demonstrated clearly in Figure 8, where we have shown the enclosed mass and the Jeans mass as a function of radius for both the and the case. Here, we calculate the Jeans mass using the average density and temperature within a given radius, extending out to the virial radius of the halo. The case is shown just prior to the formation of the first sink particle in that simulation; the case is shown at the maximum density reached over the course of the simulation. While the enclosed mass in the minihalo exceeds the Jeans mass on all scales, and thus is universally collapsing, the cold core in the minihalo never gains sufficient mass to exceed the Jeans criterion.
4.2 Sink Particle Formation and Accretion
Build-up of a Central Disc
Figure 9 shows the growth over time of all sink particles formed in our simulations, from the formation of the first sink to simulation’s end 5000 later. Sink particles are formed when the gas in the centre of the minihalo reaches , and in all cases a central disc forms within of the first sink particle, as seen in Figure 10, where we have shown the density structure of the central AU of each simulation. In the absence of any X-ray irradiation, the accretion disc fragments and forms a stable binary system, in agreement with previous studies (e.g., Stacy et al., 2010; Clark et al., 2011a, b; Greif et al., 2011a, 2012). When self-shielding is properly accounted for, the thermodynamic state of the gas in X-ray irradiated minihaloes is remarkably similar to the case by the time it reaches sink formation densities. The mass distribution, however, is in all cases significantly steeper than in the case, as shown in Figure 11. This is a direct consequence of the earlier collapse discussed in Section 4.1.1: the sooner the gas collapses to high densities, the less time is available for gas to accumulate in the centre of the minihalo.
Once sink formation densities are achieved, there is no clear trend in the subsequent behaviour of the sink particles formed—either in number or accretion rate. As seen in Table 1, both the virial mass of the minihalo and the cold core mass (defined here as ) decrease with increasing , as expected given the earlier collapse induced by a stronger CXB. However, at higher densities the case breaks this trend, suggesting that the final stages of the collapse are somewhat chaotic, and influenced more by small-scale randomness related to turbulence than by the strength of the CXB.
In the absence of shielding, the gas invariably fragments, forming a binary or small multiple within . When shielding is properly accounted for though, excess cooling of the disc (; see Figure 6) is eliminated, and fragmentation is suppressed. In fact, only a single sink particle forms in both the and cases, and while the simulation still fragments, it does so considerably less than in the absence of shielding. We may quantify this using the Toomre Q parameter (Toomre, 1964):
where is the gas sound speed, is the epicyclic frequency of the disc, and is the surface density; we replace with the orbital frequency , as appropriate for Keplerian discs. While the parameter specifically applies to infinitely thin isothermal discs, it is correct to within a factor of order unity when applied to thick discs (Wang et al., 2010), as is the case here. Figure 12 shows the parameter evaluated in mass-weighted spherical shells centred on the accretion disc after the first sink particle has formed. As the mass within these shells is dominated by the disc component, applying this analysis to the disc particles alone would have a negligible impact on the results (e.g., Greif et al., 2012). The and simulations maintain throughout the disc, and are thus susceptible to fragmentation. On the other hand, save for the central few hundred AU, where they approach the resolution limit of the simulation, the and discs stay well above , hence the lack of fragmentation.
While the minihalo also collapses early, in the same manner as the and simulations, it still experiences significant fragmentation. This is primarily due to the specific details of its collapse history, modulated by small-scale turbulence, rather than the precise value of , and can be understood by examining the velocity profile of the minihalo, shown in Figure 13 along with the Shu accretion rate (Shu, 1977) for mass-weighted spherical bins out to . While the inner regions of the accretion disc exhibit approximately Keplerian rotation in line with the other simulations, the gas at large radii experiences significantly less rotational support. This results in a much larger infall velocity, leading in turn to a higher accretion rate. This overwhelms the ability of the accretion disc to support itself against the growth of perturbations, resulting in significant fragmentation.
5 Summary and Conclusions
We have performed a suite of cosmological simulations employing a range of CXB models, focusing on the impact of such a background on Pop III stars forming in a minihalo. As three-body processes turn the gas fully molecular by , following the evolution of the gas up to this density allows us to fully capture the impact of a CXB on and cooling in the gas. After the gas reaches we form sink particles, enabling us to study the subsequent evolution of the system, which we follow for an additional 5000.
X-rays have two competing effects on primordial gas, as ionising the gas serves to both heat it and increase the number of free electrons, which catalyse the formation of molecular hydrogen. As is the main coolant in primordial gas, this actually enhances cooling. We find that heating dominates in low density gas, but is overwhelmed at higher densities by the enhanced cooling X-rays provide. The transition between these two regimes occurs between and 100, depending on the strength of the CXB.
Previous work investigating the impact of a CXB on structure formation in the early universe has generally been found to increase the supply of cold gas in a given halo (Haiman et al., 2000; Venkatesan et al., 2001; Glover & Brand, 2003; Cen, 2003; Kuhlen & Madau, 2005; Jeon et al., 2012). In particular, Kuhlen & Madau (2005) found the increase in cold gas available for star formation was greatest in – minihaloes, where the increase could exceed 1–2 orders of magnitude for moderate X-ray ionisation rates. This comports with our findings—allowing the X-ray irradiated haloes to evolve to the same cosmic time as the minihalo would result in a relatively larger supply of cold dense gas. As X-rays penetrating the minihalo begin to experience significant attenuation above though, the primary impact of a CXB is on gas in the loitering phase and below. While X-ray heating dominates below , minihaloes that overcome this impediment collapse earlier, as the cooler gas in the loitering phase requires a smaller Jeans mass to proceed to high densities.
The X-ray background is severely attenuated by the time it reaches sink formation densities. As a result its impact on the gas temperature is largely neutral, with the subsequent behaviour of the sink particles formed influenced more by small-scale turbulence than the strength of the CXB. Consequently, the characteristic mass of the stars formed is quite stable even as we vary the CXB strength by several orders of magnitude, and does not change dramatically even when the supply of cold gas in the centre of the minihalo is significantly increased, as in the case. Instead, this causes the disc to fragment, forming several protostellar cores.
It should be noted that these findings are somewhat sensitive to the density at which the minihalo becomes opaque to X-rays. While we found no difference in the total column density between simulations, a more robust estimate of the local column density would be beneficial (e.g., Hartwig et al., 2015; Safranek-Shrader et al., 2012). Additionally, there is a possibility that the more rapid minihalo collapse induced by a strong CXB may have an impact on the velocity dispersion of the infalling gas, suppressing fragmentation and possibly increasing the mass of the stars formed. While our simulations lack sufficient resolution to verify this, the findings of Clark et al. (2011a) suggest that stars forming in pre-ionised minihaloes experience more turbulence, resulting in a somewhat larger characteristic mass than in pristine minihaloes. Finally, our findings may have implications for reionisation and the 21-cm signal (Furlanetto et al., 2006; Mirocha, 2014): while the impact on the number of fragments and characteristic mass of Pop III appears to be nearly neutral, low-density gas is still smoothed by X-ray heating, thus resulting in a lower IGM clumping factor.
The authors thank Paul Clark for many insightful comments, and acknowledge the Texas Advanced Computing Center (TACC) at The University of Texas at Austin for providing HPC resources under XSEDE allocation TG-AST120024. This study was supported in part by NSF grants AST-1009928 and AST-1413501, and by the NASA grant NNX09AJ33G. A.S. gratefully acknowledges support through NSF grants AST-0908553 and AST-1211729. J.H. thanks Chalence Safranek-Shrader for many enlightening conversations. This research has made use of NASA’s Astrophysics Data System and Astropy, a community-developed core Python package for Astronomy (Robitaille et al., 2013).
- Abel T., Bryan G. L., Norman M. L., 2002, Science, 295, 93
- Alvarez M. A., Bromm V., Shapiro P. R., 2006, ApJ, 639, 621
- Barkana R., Loeb A., 2001, Phys. Rep., 349, 125
- Bate M. R., Burkert A., 1997, MNRAS, 288, 1060
- Belczynski K., Bulik T., Fryer C. L., 2012, ApJ, submitted (arXiv:1208.2422)
- Bromm V., 2013, Rep. Prog. Phys., 76, 112901
- Bromm V., Coppi P. S., Larson R. B., 1999, ApJ, 527, L5
- Bromm V., Coppi P. S., Larson R. B., 2002, ApJ, 564, 23
- Bromm V., Larson R. B., 2004, ARA&A, 42, 79
- Bromm V., Yoshida N., Hernquist L., 2003, ApJ, 596, L135
- Bromm V., Yoshida N., Hernquist L., McKee C. F., 2009, Nature, 459, 49
- Campisi M. A., Maio U., Salvaterra R., Ciardi B., 2011, MNRAS, 416, 2760
- Cen R., 2003, ApJ, 591, 12
- Clark P. C., Glover S. C. O., 2014, MNRAS, 444, 2396
- Clark P. C., Glover S. C. O., Klessen R. S., Bromm V., 2011a, ApJ, 727, 110
- Clark P. C., Glover S. C. O., Smith R. J., Greif T. H., Klessen R. S., Bromm V., 2011b, Science, 331, 1040
- Couchman H. M. P., Rees M. J., 1986, MNRAS, 221, 53
- Dijkstra M., Haiman Z., Mesinger A., Wyithe J. S. B., 2008, MNRAS, 391, 1961
- Dopcke G., Glover S. C. O., Clark P. C., Klessen R. S., 2013, ApJ, 766, 103
- Duquennoy A., Mayor M., 1991, A&A, 248, 485
- Furlanetto S. R., Peng Oh S., Briggs F. H., 2006, Cosmology at low frequencies: The 21 cm transition and the high-redshift Universe
- Furlanetto S. R., Stoever S. J., 2010, MNRAS, 404, 1869
- Gammie C. F., 2001, ApJ, 553, 174
- Glover S. C. O., Brand P. W. J. L., 2003, MNRAS, 340, 210
- Greif T. H., Bromm V., 2006, MNRAS, 373, 128
- Greif T. H., Bromm V., Clark P. C., Glover S. C. O., Smith R. J., Klessen R. S., Yoshida N., Springel V., 2012, MNRAS, 424, 399
- Greif T. H., Glover S. C. O., Bromm V., Klessen R. S., 2010, ApJ, 716, 510
- Greif T. H., Johnson J. L., Bromm V., Klessen R. S., 2007, ApJ, 670, 1
- Greif T. H., Johnson J. L., Klessen R. S., Bromm V., 2009, MNRAS, 399, 639
- Greif T. H., Springel V., White S. D. M., Glover S. C. O., Clark P. C., Smith R. J., Klessen R. S., Bromm V., 2011a, ApJ, 737, 75
- Greif T. H., White S. D. M., Klessen R. S., Springel V., 2011b, ApJ, 736, 147
- Haiman Z., Abel T., Rees M. J., 2000, ApJ, 534, 11
- Haiman Z., Thoul A. A., Loeb A., 1996, ApJ, 464, 523
- Hartwig T., Clark P. C., Glover S. C. O., Klessen R. S., Sasaki M., 2015, ApJ, 799, 114
- Heger A., Fryer C. L., Woosley S. E., Langer N., Hartmann D. H., 2003, ApJ, 591, 288
- Heger A., Woosley S. E., 2010, ApJ, 724, 341
- Hirano S., Hosokawa T., Yoshida N., Omukai K., Yorke H. W., 2015, MNRAS, 448, 568
- Hirano S., Hosokawa T., Yoshida N., Umeda H., Omukai K., Chiaki G., Yorke H. W., 2014, ApJ, 781, 60
- Hosokawa T., Omukai K., Yoshida N., Yorke H. W., 2011, Science, 334, 1250
- Inayoshi K., Omukai K., 2011, MNRAS, 416, 2748
- Jeon M., Pawlik A. H., Bromm V., Milosavljević M., 2014, MNRAS, 440, 3778
- Jeon M., Pawlik A. H., Greif T. H., Glover S. C. O., Bromm V., Milosavljević M., Klessen R. S., 2012, ApJ, 754, 34
- Johnson J. L., Greif T. H., Bromm V., 2007, ApJ, 665, 85
- Karlsson T., Bromm V., Bland-Hawthorn J., 2013, Rev. Mod. Phys., 85, 809
- Kitayama T., Yoshida N., Susa H., Umemura M., 2004, ApJ, 613, 631
- Kuhlen M., Madau P., 2005, MNRAS, 363, 1069
- Loeb A., 2010, How Did the First Stars and Galaxies Form? Princeton University Press, Princeton, NJ
- Madau P., Ferrara A., Rees M. J., 2001, ApJ, 555, 92
- Maeder A., Meynet G., 2012, Rev. Mod. Phys., 84, 25
- Maio U., Khochfar S., Johnson J. L., Ciardi B., 2011, MNRAS, 414, 1145
- Martel H., Evans II N. J., Shapiro P. R., 2006, ApJS, 163, 122
- Miralda-Escudé J., 2003, Science, 300, 1904
- Mirocha J., 2014, MNRAS, 443, 1211
- Mitsuda K. et al., 1984, Astronomical Society of Japan, 36, 741
- Mori M., Ferrara A., Madau P., 2002, ApJ, 571, 40
- Oh S. P., 2001, ApJ, 553, 499
- O’Shea B. W., Norman M. L., 2007, ApJ, 654, 66
- Osterbrock D. E., Ferland G. J., 2006, Astrophysics of gaseous nebulae and active galactic nuclei. University Science Books, Sausalito, CA
- Palla F., Salpeter E. E., Stahler S. W., 1983, ApJ, 271, 632
- Ripamonti E., Iocco F., Bressan A., Schneider R., Ferrara A., Marigo P., 2009, in Proceedings of Science
- Ripamonti E., Iocco F., Ferrara A., Schneider R., Bressan A., Marigo P., 2010, MNRAS, 406, 2605
- Robertson B. E., Ellis R. S., Dunlop J. S., McLure R. J., Stark D. P., 2010, Nature, 468, 49
- Robitaille T. P. et al., 2013, A&A, 558, A33
- Safranek-Shrader C., Agarwal M., Federrath C., Dubey A., Milosavljević M., Bromm V., 2012, MNRAS, 426, 1159
- Schaerer D., 2002, A&A, 382, 28
- Shu F. H., 1977, ApJ, 214, 488
- Shull J. M., van Steenberg M. E., 1985, ApJ, 298, 268
- Smith R. J., Glover S. C. O., Clark P. C., Greif T. H., Klessen R. S., 2011, MNRAS, 414, 3633
- Smith R. J., Iocco F., Glover S. C. O., Schleicher D. R. G., Klessen R. S., Hirano S., Yoshida N., 2012, ApJ, 761, 154
- Sokasian A., Yoshida N., Abel T., Hernquist L., Springel V., 2004, MNRAS, 350, 47
- Springel V., 2005, MNRAS, 364, 1105
- Stacy A., Bromm V., 2013, MNRAS, 433, 1094
- Stacy A., Bromm V., Loeb A., 2011, ApJ, 730, L1
- Stacy A., Greif T. H., Bromm V., 2010, MNRAS, 403, 45
- Stacy A., Greif T. H., Bromm V., 2012a, MNRAS, 422, 290
- Stacy A., Pawlik A. H., Bromm V., Loeb A., 2012b, MNRAS, 421, 894
- Stacy A., Pawlik A. H., Bromm V., Loeb A., 2014, MNRAS, 441, 822
- Susa H., Hasegawa K., Tominaga N., 2014, ApJ, 792, 32
- Tegmark M., Silk J., Rees M. J., Blanchard A., Abel T., Palla F., 1997, ApJ, 474, 1
- Toomre A., 1964, ApJ, 139, 1217
- Tornatore L., Ferrara A., Schneider R., 2007, MNRAS, 382, 945
- Turk M. J., Clark P. C., Glover S. C. O., Greif T. H., Abel T., Klessen R., Bromm V., 2011, ApJ, 726, 55
- Umeda H., Nomoto K., 2003, Nature, 422, 871
- Valdés M., Ferrara A., 2008, MNRAS, 387, L8
- Venkatesan A., Giroux M. L., Shull J. M., 2001, ApJ, 563, 1
- Wang H.-H., Klessen R. S., Dullemond C. P., Van Den Bosch F. C., Fuchs B., 2010, MNRAS, 407, 705
- Whalen D. J., Abel T., Norman M. L., 2004, ApJ, 610, 14
- Wise J. H., Abel T., 2008, ApJ, 685, 40
- Yoshida N., Abel T., Hernquist L., Sugiyama N., 2003, ApJ, 592, 645
- Yoshida N., Omukai K., Hernquist L., Abel T., 2006, ApJ, 652, 6