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.
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.
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.
ABSTRACT
Our understanding of how molecular clouds form in the interstellar medium (ISM) would be greatly helped if we had a reliable observational tracer of the gas flows responsible for forming the ...clouds. Fine structure emission from singly ionized and neutral carbon (C ii, C i) and rotational line emission from CO are all observed to be associated with molecular clouds. However, it remains unclear whether any of these tracers can be used to study the inflow of gas on to an assembling cloud, or whether they primarily trace the cloud once it has already assembled. In this paper, we address this issue with the help of high-resolution simulations of molecular cloud formation that include a sophisticated treatment of the chemistry and thermal physics of the ISM. Our simulations demonstrate that both C i and CO emission trace gas that is predominantly molecular, with a density n ∼ 500–1000 cm−3, much larger than the density of the inflowing gas. C ii traces lower density material (n ∼ 100 cm−3) that is mainly atomic at early times. A large fraction of the C ii emission traces the same structures as the C i or CO emission, but some arises in the inflowing gas. Unfortunately, this emission is very faint and will be difficult to detect with current observational facilities, even for clouds situated in regions with an elevated interstellar radiation field.
We present TreeCol, a new and efficient tree-based scheme to calculate column densities in numerical simulations. Knowing the column density in any direction at any location in space is a ...prerequisite for modelling the propagation of radiation through the computational domain. TreeCol therefore forms the basis for a fast, approximate method for modelling the attenuation of radiation within large numerical simulations. It constructs a healpix sphere at any desired location and accumulates the column density by walking the tree and by adding up the contributions from all tree nodes whose line of sight contributes to the pixel under consideration. In particular, when combined with widely-used tree-based gravity solvers, the new scheme requires little additional computational cost. In a simulation with N resolution elements, the computational cost of TreeCol scales as NlogN, instead of the N
5/3 scaling of most other radiative transfer schemes. TreeCol is naturally adaptable to arbitrary density distributions and is easy to implement and to parallelize, particularly if a tree structure is already in place for calculating the gravitational forces. We describe our new method and its implementation into the smoothed particle hydrodynamics (SPH) code gadget2 (although note that the scheme is not limited to particle-based fluid dynamics). We discuss its accuracy and performance characteristics for the examples of a spherical protostellar core and for the turbulent interstellar medium. We find that the column density estimates provided by TreeCol are on average accurate to better than 10 per cent. In another application, we compute the dust temperatures for solar neighbourhood conditions and compare with the result of a full-fledged Monte Carlo radiation-transfer calculation. We find that both methods give similar answers. We conclude that TreeCol provides a fast, easy to use and sufficiently accurate method of calculating column densities that comes with little additional computational cost when combined with an existing tree-based gravity solver.
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.
The first generation of stars, often called Population III (or Pop III), form from metal-free primordial gas at redshifts
z
∼ 30 and below. They dominate the cosmic star-formation history until
z
∼ ...15-20, at which point the formation of metal-enriched Population II stars takes over. We review current theoretical models for the formation, properties, and impact of Pop III stars and discuss existing and future observational constraints. Key takeaways from this review include the following:
Primordial gas is highly susceptible to fragmentation and Pop III stars form as members of small clusters with a logarithmically flat mass function.
Feedback from massive Pop III stars plays a central role in regulating subsequent star formation, but major uncertainties remain regarding its immediate impact.
In extreme conditions, supermassive Pop III stars can form, reaching masses of several 10
5
M
. Their remnants may be the seeds of the supermassive black holes observed in high-redshift quasars.
Direct observations of Pop III stars in the early Universe remain extremely challenging. Indirect constraints from the global 21-cm signal or gravitational waves are more promising.
Stellar archeological surveys allow us to constrain both the low-mass and the high-mass ends of the Pop III mass distribution. Observations suggest that most massive Pop III stars end their lives as core-collapse supernovae rather than as pair-instability supernovae.
We study the impact of stellar winds and supernovae on the multiphase interstellar medium using three-dimensional hydrodynamical simulations carried out with flash. The selected galactic disc region ...has a size of (500 pc)2 x plus or minus 5 kpc and a gas surface density of 10 M... pc super( -2). The simulations include an external stellar potential and gas self-gravity, radiative cooling and diffuse heating, sink particles representing star clusters, stellar winds from these clusters that combine the winds from individual massive stars by following their evolution tracks, and subsequent supernova explosions. Dust and gas (self-) shielding is followed to compute the chemical state of the gas with a chemical network. We find that stellar winds can regulate star (cluster) formation. Since the winds suppress the accretion of fresh gas soon after the cluster has formed, they lead to clusters that have lower average masses (10 super( 2)-10 super( 4.3) M...) and form on shorter time-scales (10 super( -3)-10 Myr). In particular, we find an anticorrelation of cluster mass and accretion time-scale. Without winds, the star clusters easily grow to larger masses for ~5 Myr until the first supernova explodes. Overall, the most massive stars provide the most wind energy input, while objects beginning their evolution as B-type stars contribute most of the supernova energy input. A significant outflow from the disc (mass loading ...1 at 1 kpc) can be launched by thermal gas pressure if more than 50 per cent of the volume near the disc mid-plane can be heated to T > 3 x 10 super( 5) K. Stellar winds alone cannot create a hot volume-filling phase. The models that are in best agreement with observed star formation rates drive either no outflows or weak outflows. (ProQuest: ... denotes formulae/symbols omitted.)
ABSTRACT
We use hydrodynamical simulations to study the Milky Way’s central molecular zone (CMZ). The simulations include a non-equilibrium chemical network, the gas self-gravity, star formation, and ...supernova feedback. We resolve the structure of the interstellar medium at sub-parsec resolution while also capturing the interaction between the CMZ and the bar-driven large-scale flow out to $R\sim 5\, {\rm kpc}$. Our main findings are as follows: (1) The distinction between inner (R ≲ 120 pc) and outer (120 ≲ R ≲ 450 pc) CMZ that is sometimes proposed in the literature is unnecessary. Instead, the CMZ is best described as single structure, namely a star-forming ring with outer radius R ≃ 200 pc which includes the 1.3° complex and which is directly interacting with the dust lanes that mediate the bar-driven inflow. (2) This accretion can induce a significant tilt of the CMZ out of the plane. A tilted CMZ might provide an alternative explanation to the ∞-shaped structure identified in Herschel data by Molinari et al. (3) The bar in our simulation efficiently drives an inflow from the Galactic disc (R ≃ 3 kpc) down to the CMZ (R ≃ 200 pc) of the order of $1\rm \, M_\odot \, yr^{-1}$, consistent with observational determinations. (4) Supernova feedback can drive an inflow from the CMZ inwards towards the circumnuclear disc of the order of ${\sim}0.03\, \rm M_\odot \, yr^{-1}$. (5) We give a new interpretation for the 3D placement of the 20 and 50 km s−1 clouds, according to which they are close (R ≲ 30 pc) to the Galactic Centre, but are also connected to the larger scale streams at R ≳ 100 pc.