Carbon monoxide (CO) is widely used as a tracer of molecular hydrogen (H2) in metal-rich galaxies, but is known to become ineffective in low-metallicity dwarf galaxies. Atomic carbon has been ...suggested as a superior tracer of H2 in these metal-poor systems, but its suitability remains unproven. To help us to assess how well atomic carbon traces H2 at low metallicity, we have performed a series of numerical simulations of turbulent molecular clouds that cover a wide range of different metallicities. Our simulations demonstrate that in star-forming clouds, the conversion factor between C i emission and H2 mass, X
CI, scales approximately as X
CI ∝ Z
−1. We recover a similar scaling for the CO-to-H2 conversion factor, X
CO, but find that at this point in the evolution of the clouds, X
CO is consistently smaller than X
CI, by a factor of a few or more. We have also examined how X
CI and X
CO evolve with time. We find that X
CI does not vary strongly with time, demonstrating that atomic carbon remains a good tracer of H2 in metal-poor systems even at times significantly before the onset of star formation. On the other hand, X
CO varies very strongly with time in metal-poor clouds, showing that CO does not trace H2 well in starless clouds at low metallicity.
The very first stars to form in the universe heralded an end to the cosmic dark ages and introduced new physical processes that shaped early cosmic evolution. Until now, it was thought that these ...stars lived short, solitary lives, with only one extremely massive star, or possibly a very wide binary system, forming in each dark-matter minihalo. Here we describe numerical simulations that show that these stars were, to the contrary, often members of tight multiple systems. Our results show that the disks that formed around the first young stars were unstable to gravitational fragmentation, possibly producing small binary and higher-order systems that had separations as small as the distance between Earth and the Sun.
The direct collapse model for the formation of massive seed black holes in the early Universe attempts to explain the observed number density of supermassive black holes (SMBHs) at z ∼ 6 by assuming ...that they grow from seeds with masses M > 104 M⊙ that form by the direct collapse of metal-free gas in atomic cooling haloes in which H2 cooling is suppressed by a strong extragalactic radiation field. The viability of this model depends on the strength of the radiation field required to suppress H2 cooling, J
crit: if this is too large, then too few seeds will form to explain the observed number density of SMBHs. In order to determine J
crit reliably, we need to be able to accurately model the formation and destruction of H2 in gas illuminated by an extremely strong radiation field. In this paper, we use a reaction-based reduction technique to analyse the chemistry of H2 in these conditions, allowing us to identify the key chemical reactions that are responsible for determining the value of J
crit. We construct a reduced network of 26 reactions that allows us to determine J
crit accurately, and compare it with previous treatments in the literature. We show that previous studies have often omitted one or more important chemical reactions, and that these omissions introduce an uncertainty of up to a factor of 3 into previous determinations of J
crit.
The SILCC (SImulating the Life-Cycle of molecular Clouds) project aims to self-consistently understand the small-scale structure of the interstellar medium (ISM) and its link to galaxy evolution. We ...simulate the evolution of the multiphase ISM in a (500 pc)2 × ±5 kpc region of a galactic disc, with a gas surface density of
$\Sigma _{_{\rm GAS}} = 10 \;{\rm M}_{\odot }\,{\rm pc}^{-2}$
. The flash 4 simulations include an external potential, self-gravity, magnetic fields, heating and radiative cooling, time-dependent chemistry of H2 and CO considering (self-) shielding, and supernova (SN) feedback but omit shear due to galactic rotation. We explore SN explosions at different rates in high-density regions (peak), in random locations with a Gaussian distribution in the vertical direction (random), in a combination of both (mixed), or clustered in space and time (clus/clus2). Only models with self-gravity and a significant fraction of SNe that explode in low-density gas are in agreement with observations. Without self-gravity and in models with peak driving the formation of H2 is strongly suppressed. For decreasing SN rates, the H2 mass fraction increases significantly from <10 per cent for high SN rates, i.e. 0.5 dex above Kennicutt–Schmidt, to 70–85 per cent for low SN rates, i.e. 0.5 dex below KS. For an intermediate SN rate, clustered driving results in slightly more H2 than random driving due to the more coherent compression of the gas in larger bubbles. Magnetic fields have little impact on the final disc structure but affect the dense gas (n ≳ 10 cm−3) and delay H2 formation. Most of the volume is filled with hot gas (∼80 per cent within ±150 pc). For all but peak driving a vertically expanding warm component of atomic hydrogen indicates a fountain flow. We highlight that individual chemical species populate different ISM phases and cannot be accurately modelled with temperature-/density-based phase cut-offs.
The SILCC project (SImulating the Life-Cycle of molecular Clouds) aims at a more self-consistent understanding of the interstellar medium (ISM) on small scales and its link to galaxy evolution. We ...present three-dimensional (magneto)hydrodynamic simulations of the ISM in a vertically stratified box including self-gravity, an external potential due to the stellar component of the galactic disc, and stellar feedback in the form of an interstellar radiation field and supernovae (SNe). The cooling of the gas is based on a chemical network that follows the abundances of H+, H, H2, C+, and CO and takes shielding into account consistently. We vary the SN feedback by comparing different SN rates, clustering and different positioning, in particular SNe in density peaks and at random positions, which has a major impact on the dynamics. Only for random SN positions the energy is injected in sufficiently low-density environments to reduce energy losses and enhance the effective kinetic coupling of the SNe with the gas. This leads to more realistic velocity dispersions (
$\sigma _\mathrm{H\,{\small {I}}}\approx 0.8\sigma _{300\rm{-}8000\,\mathrm{K}}\sim 10\hbox{-}20\,\mathrm{km}\,\mathrm{s}^{-1}$
,
$\sigma _\mathrm{H\,\alpha }\approx 0.6\sigma _{8000-3\times 10^5\,\mathrm{K}}\sim 20\hbox{-}30\,\mathrm{km}\,\mathrm{s}^{-1}$
), and strong outflows with mass loading factors (ratio of outflow to star formation rate) of up to 10 even for solar neighbourhood conditions. Clustered SNe abet the onset of outflows compared to individual SNe but do not influence the net outflow rate. The outflows do not contain any molecular gas and are mainly composed of atomic hydrogen. The bulk of the outflowing mass is dense (ρ ∼ 10−25–10−24 g cm−3) and slow (v ∼ 20–40 km s−1) but there is a high-velocity tail of up to v ∼ 500 km s−1 with ρ ∼ 10−28–10−27 g cm−3.
We investigate how uncertainties in the chemical and cooling rate coefficients relevant for a metal-free gas influence our ability to determine the critical ultraviolet field strength required to ...suppress H2 cooling in high-redshift atomic cooling haloes. The suppression of H2 cooling is a necessary prerequisite for the gas to undergo direct collapse and form an intermediate mass black hole. These black holes can then act as seeds for the growth of the supermassive black holes (SMBHs) observed at redshifts z ∼ 6. The viability of this model for SMBH formation depends on the critical ultraviolet field strength, J
crit: if this is too large, then too few seeds will form to explain the observed number density of SMBHs. We show in this paper that there are five key chemical reactions whose rate coefficients are uncertain enough to significantly affect J
crit. The most important of these is the collisional ionization of hydrogen by collisions with other hydrogen atoms, as the rate for this process is very poorly constrained at the low energies relevant for direct collapse. The total uncertainty introduced into J
crit by this and the other four reactions could in the worst case approach a factor of five. We also show that the use of outdated or inappropriate values for the rates of some chemical reactions in previous studies of the direct collapse mechanism may have significantly affected the values of J
crit determined by these studies.
We use the moving mesh code arepo coupled to a time-dependent chemical network to investigate molecular gas in simulated spiral galaxies that is not traced by CO emission. We calculate H2 and CO ...column densities, and estimate the CO emission and CO–H2 conversion factor. We find that in conditions akin to those in the local interstellar medium, around 42 per cent of the total molecular mass should be in CO-dark regions, in reasonable agreement with observational estimates. This fraction is almost insensitive to the CO-integrated intensity threshold used to discriminate between CO-bright and CO-dark gas. The CO-dark molecular gas primarily resides in extremely long (>100 pc) filaments that are stretched between spiral arms by galactic shear. Only the centres of these filaments are bright in CO, suggesting that filamentary molecular clouds observed in the Milky Way may only be small parts of much larger structures. The CO-dark molecular gas mainly exists in a partially molecular phase which accounts for a significant fraction of the total disc mass budget. The dark gas fraction is higher in simulations with higher ambient UV fields or lower surface densities, implying that external galaxies with these conditions might have a greater proportion of dark gas.
Abstract
We present 3D ‘zoom-in’ simulations of the formation of two molecular clouds out of the galactic interstellar medium. We model the clouds – identified from the SILCC simulations – with a ...resolution of up to 0.06 pc using adaptive mesh refinement in combination with a chemical network to follow heating, cooling and the formation of H2 and CO including (self-) shielding. The two clouds are assembled within a few million years with mass growth rates of up to ∼10−2 M⊙ yr−1 and final masses of ∼50 000 M⊙. A spatial resolution of ≲0.1 pc is required for convergence with respect to the mass, velocity dispersion and chemical abundances of the clouds, although these properties also depend on the cloud definition such as based on density thresholds, H2 or CO mass fraction. To avoid grid artefacts, the progressive increase of resolution has to occur within the free-fall time of the densest structures (1–1.5 Myr) and ≳200 time-steps should be spent on each refinement level before the resolution is progressively increased further. This avoids the formation of spurious, large-scale, rotating clumps from unresolved turbulent flows. While CO is a good tracer for the evolution of dense gas with number densities n ≥ 300 cm−3, H2 is also found for n ≲ 30 cm−3 due to turbulent mixing and becomes dominant at column densities around 30–50 M⊙ pc−2. The CO-to-H2 ratio steadily increases within the first 2 Myr, whereas XCO ≃ 1–4 × 1020 cm−2 (K km s−1)−1 is approximately constant since the CO(1−0) line quickly becomes optically thick.
We use a suite of high-resolution molecular cloud simulations carried out with the moving mesh code arepo to explore the nature of star-forming filaments. The simulated filaments are identified and ...categorized from column density maps in the same manner as for recent Herschel observations. When fitted with a Plummer-like profile, the filaments are in excellent agreement with observations, and have shallow power-law profiles of p ∼ 2.2 without the need for magnetic support. When data within 1 pc of the filament centre are fitted with a Gaussian function, the average full width at half-maximum (FWHM) is ∼0.3 pc, in agreement with predictions for accreting filaments. However, if the fit is constructed using only the inner regions, as in Herschel observations, the resulting FWHM is only ∼0.2 pc. This value is larger than that measured in IC 5146 and Taurus, but is similar to that found in the Planck Galactic cold cores and in Cygnus X. The simulated filaments have a range of widths rather than a constant value. When the column density maps are compared to the 3D gas densities, the filaments seen in column density do not belong to a single structure. Instead, they are made up of a network of short ribbon-like sub-filaments reminiscent of those seen in Taurus. The sub-filaments are pre-existing within the simulated clouds, have radii similar to their Jeans radius, and are not primarily formed through fragmentation of the larger filament seen in column density. Instead, small filamentary clumps are swept together into a single column density structure by the large-scale collapse of the cloud.
We show that hydrodynamic turbulent cloud simulations naturally produce large filaments made up of a network of smaller and coherent subfilaments. Such simulations resemble observations of filaments ...and fibres in nearby molecular clouds. The subfilaments are dynamical features formed at the stagnation points of the turbulent velocity field where shocks dissipate the turbulent energy. They are a ubiquitous feature of the simulated clouds, which appear from the beginning of the simulation and are not formed by gradual fragmentation of larger filaments. Most of the subfilaments are gravitationally subcritical and do not fragment into cores, however, there is also a significant fraction of supercritical subfilaments which break up into star-forming cores. The subfilaments are coherent along their length, and the residual velocities along their spine show that they are subsonically contracting without any ordered rotation on scales of ∼0.1 pc. Accretion flows along the subfilaments can feed material into star-forming cores embedded within the network. The overall mass in subfilaments and the number of subfilaments increases as the cloud evolves. We propose that the formation of filaments and subfilaments is a natural consequence of the turbulent cascade in the complex multiphase interstellar medium. Subfilaments are formed by the high wavenumber, small-scale modes in the turbulent velocity field. These are then stretched by local shear motions and gathered together by a combination of low wavenumber modes and gravitational contraction on larger scales, and by doing so build up the extended filaments observed in column density maps.