We present the first numerical simulations that self-consistently follow the formation of dense molecular clouds in colliding flows. Our calculations include a time-dependent model for the H2 and CO chemistry that runs alongside a detailed treatment
of the dominant heating and cooling processes in the ISM. We adopt initial conditions characteristic of the warm neutral medium and study two different flow velocities - a slow flow with v = 6.8 km/s and a fast flow with v = 13.6 km/s. The clouds formed by the collision of these flows form stars, with star formation beginning after 16 Myr in the case of the slower flow, but after only 4.4 Myr in the case of the faster flow. In both flows, the formation of CO-dominated regions occurs only around 2 Myr before the onset of star formation. Prior to this, the clouds produce very little emission in the J = 1 -> 0 transition line of CO, and would probably not be identified as molecular clouds in observational surveys. In contrast, our models show that H2-dominated regions can form much earlier, with the timing depending on the details of the flow. In the case of the slow flow, small pockets of gas become fully molecular around 10 Myr before star formation begins, while in the fast flow, the first H2-dominated regions occur around 3 Myr before the first prestellar cores form. Our results are consistent with models of molecular cloud formation in which the clouds are dominated by dark molecular gas for a considerable proportion of their assembly history.
We present a new systematic way of setting up galactic gas disks based on the assumption of detailed hydrodynamic equilibrium. To do this, we need to specify the density distribution and the velocity field which supports the disk. We first show that
the required circular velocity has no dependence on the height above or below the midplane so long as the gas pressure is a function of density only. The assumption of disks being very thin enables us to decouple the vertical structure from the radial direction. Based on that, the equation of hydrostatic equilibrium together with the reduced Poisson equation leads to two sets of second-order non-linear differential equation, which are easily integrated to set-up a stable disk. We call one approach `density method and the other one `potential method. Gas disks in detailed balance are especially suitable for investigating the onset of the gravitational instability. We revisit the question of global, axisymmetric instability using fully three-dimensional disk simulations. The impact of disk thickness on the disk instability and the formation of spontaneously induced spirals is studied systematically with or without the presence of the stellar potential. In our models, the numerical results show that the threshold value for disk instability is shifted from unity to 0.69 for self-gravitating thick disks and to 0.75 for combined stellar and gas thick disks. The simulations also show that self-induced spirals occur in the correct regions and with the right numbers as predicted by the analytic theory.
Primordial star formation appears to result in stars at least an order of magnitude more massive than modern star formation. It has been proposed that the transition from primordial to modern initial mass functions occurs due to the onset of effectiv
e metal line cooling at a metallicity Z/Z_sun = 10^{-3.5}. However, these simulations neglected molecular hydrogen cooling. We perform simulations using the same initial conditions, but including molecular cooling, using a complex network that follows molecular hydrogen formation and also directly follows carbon monoxide and water. We find that molecular hydrogen cooling allows roughly equivalent fragmentation to proceed even at zero metallicity for these initial conditions. The apparent transition just represents the point where metal line cooling becomes more important than molecular cooling. In all cases, the fragments are massive enough to be consistent with models of primordial stellar masses, suggesting that the transition to the modern initial mass function may be determined by other physics such as dust formation. We conclude that such additional cooling mechanisms, combined with the exact initial conditions produced by cosmological collapse are likely more important than metal line cooling in determining the initial mass function, and thus that there is unlikely to be a sharp transition in the initial mass function at Z/Z_sun = 10^{-3.5}.
Dark stars powered by dark matter annihilation have been proposed as the first luminous sources in the universe. These stars are believed to form in the central dark matter cusp of low-mass minihalos. Recent calculations indicate stellar masses up to
sim1000 solar masses and/or have very long lifetimes. The UV photons from these objects could therefore contribute significantly to cosmic reionization. Here we show that such dark star models would require a somewhat artificial reionization history, based on a double-reionization phase and a late star-burst near redshift $zsim6$, in order to fulfill the WMAP constraint on the optical depth as well as the Gunn-Peterson constraint at $zsim6$. This suggests that, if dark stars were common in the early universe, then models are preferred which predict a number of UV photons similar to conventional Pop. III stars. This excludes dark stars with 100 solar masses that enter a main-sequence phase and other models that lead to a strong increase in the number of UV photons. We also derive constraints for massive as well as light dark matter candidates from the observed X-ray, gamma-ray and neutrino background, considering dark matter profiles which have been steepened during the formation of dark stars. This increases the clumping factor at high redshift and gives rise to a higher dark matter annihilation rate in the early universe. We furthermore estimate the potential contribution from the annihilation products in the remnants of dark stars, which may provide a promising path to constrain such models further, but which is currently still uncertain.
Magnetic fields in the early universe can significantly alter the thermal evolution and the ionization history during the dark ages. This is reflected in the 21 cm line of atomic hydrogen, which is coupled to the gas temperature through collisions at
high redshifts, and through the Wouthuysen-Field effect at low redshifts. We present a semi-analytic model for star formation and the build-up of a Lyman alpha background in the presence of magnetic fields, and calculate the evolution of the mean 21 cm brightness temperature and its frequency gradient as a function of redshift. We further discuss the evolution of linear fluctuations in temperature and ionization in the presence of magnetic fields and calculate the effect on the 21 cm power spectrum. At high redshifts, the signal is increased compared to the non-magnetic case due to the additional heat input into the IGM from ambipolar diffusion and the decay of MHD turbulence. At lower redshifts, the formation of luminous objects and the build-up of a Lyman alpha background can be delayed by a redshift interval of 10 due to the strong increase of the filtering mass scale in the presence of magnetic fields. This tends to decrease the 21 cm signal compared to the zero-field case. In summary, we find that 21 cm observations may become a promising tool to constrain primordial magnetic fields.
We calculate the reionization history for different models of the stellar population and explore the effects of primordial magnetic fields, dark matter decay and dark matter annihilation on reionization. We find that stellar populations based on a Sc
alo-type initial mass function for Population II stars can be ruled out as sole sources for reionization, unless star formation efficiencies of more than 10% or very high photon escape fractions from the parental halo are adopted. When considering primordial magnetic fields, we find that the additional heat injection from ambipolar diffusion and decaying MHD turbulence has significant impact on the thermal evolution and the ionization history of the post-recombination universe and on structure formation. The magnetic Jeans mass changes the typical mass scale of the star forming halos, and depending on the adopted stellar model we derive upper limits to the magnetic field strength between 0.7 and $5 $nG (comoving). For dark matter annihilation, we find an upper limit to the thermally averaged mass-weighted cross section of $10^{-33} mathrm{cm}^3mathrm{/s/eV}$. For dark matter decay, our calculations yield a lower limit to the lifetime of dark matter particles of $3times10^{23}$ s. These limits are in agreement with constraints from recombination and the X-ray background and provide an independent confirmation at a much later epoch.
We present 3D radiation-gasdynamical simulations of an ionization front running into a dense clump. In our setup, a B0 star irradiates an overdensity which is at a distance of 10 pc and modelled as a supercritical 100 M_sol Bonnor-Ebert sphere. The r
adiation from the star heats up the gas and creates a shock front that expands into the interstellar medium. The shock compresses the clump material while the ionizing radiation heats it up. The outcome of this cloud-crushing process is a fully turbulent gas in the wake of the clump. In the end, the clump entirely dissolves. We propose that this mechanism is very efficient in creating short-living supersonic turbulence in the vicinity of massive stars.
Jets and outflows from young stellar objects are proposed candidates to drive supersonic turbulence in molecular clouds. Here, we present the results from multi-dimensional jet simulations where we investigate in detail the energy and momentum deposi
tion from jets into their surrounding environment and quantify the character of the excited turbulence with velocity probability density functions. Our study include jet--clump interaction, transient jets, and magnetised jets. We find that collimated supersonic jets do not excite supersonic motions far from the vicinity of the jet. Supersonic fluctuations are damped quickly and do not spread into the parent cloud. Instead subsonic, non-compressional modes occupy most of the excited volume. This is a generic feature which can not be fully circumvented by overdense jets or magnetic fields. Nevertheless, jets are able to leave strong imprints in their cloud structure and can disrupt dense clumps. Our results question the ability of collimated jets to sustain supersonic turbulence in molecular clouds.
