Search

Aims. We introduce a new deep-learning approach for the reconstruction of 3D dust density and temperature distributions from multi-wavelength dust emission observations on the scale of individual star-forming cloud cores (<0.2 pc). Methods. We constructed a training data set by processing cloud cores from the Cloud Factory simulations with the POLARIS radiative transfer code to produce synthetic dust emission observations at 23 wavelengths between 12 and 1300 µm. We simplified the task by reconstructing the cloud structure along individual lines of sight (LoSs) and trained a conditional invertible neural network (cINN) for this purpose. The cINN belongs to the group of normalising flow methods and it is able to predict full posterior distributions for the target dust properties. We tested different cINN setups, ranging from a scenario that includes all 23 wavelengths down to a more realistically limited case with observations at only seven wavelengths. We evaluated the predictive performance of these models on synthetic test data. Results. We report an excellent reconstruction performance for the 23-wavelength cINN model, achieving median absolute relative errors of about 1.8% in log(n/m−3) and 1% in log(Tdust/K), respectively. We identify trends towards an overestimation at the low end of the density range and towards an underestimation at the high end of both the density and temperature values, which may be related to a bias in the training data. After limiting our coverage to a combination of only seven wavelengths, we still find a satisfactory performance with average absolute relative errors of about 2.8% and 1.7% in log(n/m−3) and log(Tdust/K). Conclusions. This proof-of-concept study shows that the cINN-based approach for 3D reconstruction of dust density and temperature is very promising and it is even compatible with a more realistically constrained wavelength coverage. [ABSTRACT FROM AUTHOR]

The Central Molecular Zone (CMZ; the central ∼500 pc of the Galaxy) is a kinematically unusual environment relative to the Galactic disc, with high-velocity dispersions and a steep size–linewidth relation of the molecular clouds. In addition, the CMZ region has a significantly lower star formation rate (SFR) than expected by its large amount of dense gas. An important factor in explaining the low SFR is the turbulent state of the star-forming gas, which seems to be dominated by rotational modes. However, the turbulence driving mechanism remains unclear. In this work, we investigate how the Galactic gravitational potential affects the turbulence in CMZ clouds. We focus on the CMZ cloud G0.253+0.016 ('the Brick'), which is very quiescent and unlikely to be kinematically dominated by stellar feedback. We demonstrate that several kinematic properties of the Brick arise naturally in a cloud-scale hydrodynamics simulation, that takes into account the Galactic gravitational potential. These properties include the line-of-sight velocity distribution, the steepened size–linewidth relation, and the predominantly solenoidal nature of the turbulence. Within the simulation, these properties result from the Galactic shear in combination with the cloud's gravitational collapse. This is a strong indication that the Galactic gravitational potential plays a crucial role in shaping the CMZ gas kinematics, and is a major contributor to suppressing the SFR, by inducing predominantly solenoidal turbulent modes. [ABSTRACT FROM AUTHOR]

The formation of supermassive stars is believed to be an essential intermediate step for the formation of the massive black hole seeds that become the supermassive black holes powering the quasars observed in the early Universe. Numerical simulations have shown that supermassive stars can form in atomic-cooling haloes when protostars reach accretion rates higher than ∼10−2 M⊙ yr−1 and fragmentation is suppressed on pc scales. It is, however, still uncertain if a supermassive star still emerges when fragmentation occurs at smaller scales and a cluster of stars is formed instead. In this work, we explore the problem of massive object formation due to the interplay of collisions and accretion in star clusters at low metallicity. We model a small embedded cluster of accreting protostars following subparsec scale fragmentation during the collapse of a primordial gas cloud, and follow its evolution by performing N -body plus hydrodynamical simulations. Our results show that supermassive stars with 103 and 104 M⊙ are always formed due to the interplay of collisions and accretion, and in some cases these objects are part of a binary system. The resulting supermassive star is surrounded by tens of smaller stars with typical masses in the range 1–100 M⊙. [ABSTRACT FROM AUTHOR]

We study the evolution of embedded star clusters with the goal to reproduce 30 Doradus, specifically the compact star cluster known as R136 and its surrounding stellar envelope, which is believed to be part of an earlier star formation event. We employ the high-precision stellar dynamics code Nbody6+ + GPU to calculate the dynamics of the stars embedded in different evolving molecular clouds modelled with the 1D cloud/clusters code warpfield. We explore clouds with initial masses of M cloud = 3.16 × 105 M⊙ that (re)-collapse allowing for the birth of a second generation. We explore different star formation efficiencies to find the best set of parameters that can reproduce the observations. Our best-fit models correspond to a first generation of stars with a total mass M in the range |$1.26 \!-\! 2.85\times \,\,10^4\,$|  M⊙. As the initial stellar feedback is insufficient to unbind the parental cloud, the gas re-collapses after about 2–4 million years and builds up a second generation of stars with M ≈ 6.32 × 104 M⊙. We can match the observed stellar ages, the radius of the shell of swept up cloud material, and the fact that the second generation of stars is more concentrated than the first one. This is independent of the cluster starting out with mass segregation or without. By comparing with recent measurements of mass segregation and density profile in the central region of the cluster we again find close agreement, providing further evidence for a re-collapse scenario building up multiple generations of stars in 30 Doradus. [ABSTRACT FROM AUTHOR]

Many studies concluded that magnetic fields suppress star formation in molecular clouds and Milky Way like galaxies. However, most of these studies are based on fully developed fields that have reached the saturation level, with little work on investigating how an initial weak primordial field affects star formation in low metallicity environments. In this paper, we investigate the impact of a weak initial field on low metallicity dwarf galaxies. We perform high-resolution arepo simulations of five isolated dwarf galaxies. Two models are hydrodynamical, two start with a primordial magnetic field of 10 |$^{-6} \, \mu$| G and different sub-solar metallicities, and one starts with a saturated field of 10 |$^{-2} \, \mu$| G. All models include a non-equilibrium, time-dependent chemical network that includes the effects of gas shielding from the ambient ultraviolet field. Sink particles form directly from the gravitational collapse of gas and are treated as star-forming clumps that can accrete gas. We vary the ambient uniform far ultraviolet field, and cosmic ray ionization rate between 1 per cent and 10 per cent of solar values. We find that the magnetic field has little impact on the global star formation rate (SFR), which is in tension with some previously published results. We further find that the initial field strength has little impact on the global SFR. We show that an increase in the mass fractions of both molecular hydrogen and cold gas, along with changes in the perpendicular gas velocity dispersion and the magnetic field acting in the weak-field model, overcome the expected suppression in star formation. [ABSTRACT FROM AUTHOR]

Collisions between giant molecular clouds (GMCs) are one of the pathways for massive star formation due to the high densities created. However, the enhancement of the star formation rate (SFR) is not well constrained. In this study, we perform a parameter study of cloud–cloud collisions and investigate how the resulting SFR depends on the details of set-up. Our parameter study explores variations in collision speed, magnetic field inclination (with respect to the collisional axis), and resolution, as defined by the number of cells per Jeans length. In all our collision simulations, we find a factor of 2–3 increase in the SFR compared to our no collision simulation, with star formation beginning sooner with (a) high collisional velocities, (b) parallel orientation between the magnetic field and collision axis, (c) and lower resolution. The mean virial parameter of high density (and thus possible star-forming) gas increases with collisional velocity, but has little variation with magnetic field inclination. The alignment of the velocity and magnetic field remains uniform in low-density environments but becomes more perpendicular with increasing density, indicating the compression of the magnetic field by collapsing gas. Comparing the trends in the SFR with other GMC collision studies, we find good agreement with studies that account for the gravitational boundedness of the gas in their star formation algorithm, but not with those that simply form stars above a prescribed density threshold. This suggests that the latter approach should be used with caution when modelling star formation on resolved cloud scales. [ABSTRACT FROM AUTHOR]

During primordial star formation, the main cooling channel is provided by H2 and super-molecules, such as H2 or H2, at sufficiently high densities. When the latter form at |$n_{\rm H}\, \ge \, 10^{14}$|  cm−3, collision-induced emission (CIE) provides efficient gas cooling. We investigate how CIE cooling affects the formation of metal-free binaries comparing simulations with and without this process. Irrespective of the cooling mechanism, we find a typical protostellar mass range between 0.01 and 100 M⊙. However, models with only H2 line cooling produce a greater number of low-mass protostars that exhibit stronger variations in their radial velocities than the high-mass protostars. Similarly, in models with both H2 cooling and CIE cooling, significant variations in the radial velocities are found for protostars in the intermediate-mass range. The initial number of fragments N max decreases with increasing strength of turbulence. Cooling via super-molecules lets the most massive protobinaries (MMPBs) efficiently accrete mass. The maximum mass accretion rate |$\dot{M}_{\rm max}$| for the MMPBs is more than an order of magnitude higher in the presence of CIE cooling than for pure H2 line cooling. As a result, compact binaries with a semimajor axis as small as 3.57 au may form through the H2– H2 cooling channel. Our results indicate that, in addition to the MMPBs, most population III (Pop. III) binaries should be in eccentric i.e. non-circular orbits. This provides an important connection to the eccentric binaries reported in previous studies, which were found to exhibit rich temporal accretion signals during their evolution. [ABSTRACT FROM AUTHOR]

The processes of star formation and feedback, regulating the cycle of matter between gas and stars on the scales of giant molecular clouds (GMCs; ∼100 pc), play a major role in governing galaxy evolution. Measuring the time-scales of GMC evolution is important to identify and characterize the specific physical mechanisms that drive this transition. By applying a robust statistical method to high-resolution CO and narrow-band H α imaging from the PHANGS survey, we systematically measure the evolutionary timeline from molecular clouds to exposed young stellar regions on GMC scales, across the discs of an unprecedented sample of 54 star-forming main-sequence galaxies (excluding their unresolved centres). We find that clouds live for about 1−3 GMC turbulence crossing times (5−30 Myr) and are efficiently dispersed by stellar feedback within 1−5 Myr once the star-forming region becomes partially exposed, resulting in integrated star formation efficiencies of 1−8 per cent. These ranges reflect physical galaxy-to-galaxy variation. In order to evaluate whether galactic environment influences GMC evolution, we correlate our measurements with average properties of the GMCs and their local galactic environment. We find several strong correlations that can be physically understood, revealing a quantitative link between galactic-scale environmental properties and the small-scale GMC evolution. Notably, the measured CO-visible cloud lifetimes become shorter with decreasing galaxy mass, mostly due to the increasing presence of CO-dark molecular gas in such environment. Our results represent a first step towards a comprehensive picture of cloud assembly and dispersal, which requires further extension and refinement with tracers of the atomic gas, dust, and deeply embedded stars. [ABSTRACT FROM AUTHOR]

In recent years, it has been speculated that in extreme low-metallicity galactic environments, stars form in regions that lack H2. In this paper, we investigate how changing the metallicity and ultraviolet (UV) field strength of a galaxy affects the star formation within, and the molecular gas Kennicutt–Schmidt (KS) relation. Using extremely high-resolution arepo simulations of isolated dwarf galaxies, we independently vary the metallicity and UV field to between 1 per cent and 10 per cent solar neighbourhood values. We include a non-equilibrium, time-dependent chemical network to model the molecular composition of the interstellar medium and include the effects of gas shielding from an ambient UV field. Crucially, our simulations directly model the gravitational collapse of gas into star-forming clumps and cores and their subsequent accretion using sink particles. In this first publication, we find that reducing the metallicity and UV field by a factor of 10 has no effect on star formation and minimal effect on the cold, dense star-forming gas. The cold gas depletion times are almost an order of magnitude longer than the molecular gas depletion time due to the presence of star formation in H  i dominated cold gas. We study the H2 KS relationship that arises naturally within the simulations and find a near-linear power-law index of N  = 1.09 ± 0.014 in our fiducial |$10{{\ \rm per\ cent}}$| solar metallicity model. As the metallicity and UV field are reduced, this becomes moderately steeper, with a slope of N  = 1.24 ± 0.022 for our |$1{{\ \rm per\ cent}}$| solar metallicity and |$1{{\ \rm per\ cent}}$| solar UV field model. [ABSTRACT FROM AUTHOR]

To investigate how molecular clouds react to different environmental conditions at a galactic scale, we present a catalogue of giant molecular clouds (GMCs) resolved down to masses of ∼10 M⊙ from a simulation of the entire disc of an interacting M51-like galaxy and a comparable isolated galaxy. Our model includes time-dependent gas chemistry, sink particles for star formation, and supernova feedback, meaning we are not reliant on star formation recipes based on threshold densities and can follow the physics of the cold molecular phase. We extract GMCs from the simulations and analyse their properties. In the disc of our simulated galaxies, spiral arms seem to act merely as snowplows, gathering gas, and clouds without dramatically affecting their properties. In the centre of the galaxy, on the other hand, environmental conditions lead to larger, more massive clouds. While the galaxy interaction has little effect on cloud masses and sizes, it does promote the formation of counter-rotating clouds. We find that the identified clouds seem to be largely gravitationally unbound at first glance, but a closer analysis of the hierarchical structure of the molecular interstellar medium shows that there is a large range of virial parameters with a smooth transition from unbound to mostly bound for the densest structures. The common observation that clouds appear to be virialized entities may therefore be due to CO bright emission highlighting a specific level in this hierarchical binding sequence. The small fraction of gravitationally bound structures found suggests that low galactic star formation efficiencies may be set by the process of cloud formation and initial collapse. [ABSTRACT FROM AUTHOR]

Feedback from massive stars plays a key role in molecular cloud evolution. After the onset of star formation, the young stellar population is exposed by photoionization, winds, supernovae, and radiation pressure from massive stars. Recent observations of nearby galaxies have provided the evolutionary timeline between molecular clouds and exposed young stars, but the duration of the embedded phase of massive star formation is still ill-constrained. We measure how long massive stellar populations remain embedded within their natal cloud, by applying a statistical method to six nearby galaxies at |$20{-}100~\mbox{${\rm ~pc}$}$| resolution, using CO, Spitzer 24 |$\rm \, \mu m$|⁠ , and H α emission as tracers of molecular clouds, embedded star formation, and exposed star formation, respectively. We find that the embedded phase (with CO and 24 |$\rm \, \mu m$|  emission) lasts for 2−7 Myr and constitutes |$17{-}47{{\ \rm per\ cent}}$| of the cloud lifetime. During approximately the first half of this phase, the region is invisible in H α, making it heavily obscured. For the second half of this phase, the region also emits in H α and is partially exposed. Once the cloud has been dispersed by feedback, 24 |$\rm \, \mu m$| emission no longer traces ongoing star formation, but remains detectable for another 2−9 Myr through the emission from ambient CO-dark gas, tracing star formation that recently ended. The short duration of massive star formation suggests that pre-supernova feedback (photoionization and winds) is important in disrupting molecular clouds. The measured time-scales do not show significant correlations with environmental properties (e.g. metallicity). Future JWST observations will enable these measurements routinely across the nearby galaxy population. [ABSTRACT FROM AUTHOR]

We present here the first of a series of papers aimed at better understanding the evolution and properties of giant molecular clouds (GMCs) in a galactic context. We perform high-resolution, three-dimensional arepo simulations of an interacting galaxy inspired by the well-observed M51 galaxy. Our fiducial simulations include a non-equilibrium, time-dependent, chemical network that follows the evolution of atomic and molecular hydrogen as well as carbon and oxygen self-consistently. Our calculations also treat gas self-gravity and subsequent star formation (described by sink particles), and coupled supernova feedback. In the densest parts of the simulated interstellar medium (ISM), we reach sub-parsec resolution, granting us the ability to resolve individual GMCs and their formation and destruction self-consistently throughout the galaxy. In this initial work, we focus on the general properties of the ISM with a particular focus on the cold star-forming gas. We discuss the role of the interaction with the companion galaxy in generating cold molecular gas and controlling stellar birth. We find that while the interaction drives large-scale gas flows and induces spiral arms in the galaxy, it is of secondary importance in determining gas fractions in the different ISM phases and the overall star formation rate. The behaviour of the gas on small GMC scales instead is mostly controlled by the self-regulating property of the ISM driven by coupled feedback. [ABSTRACT FROM AUTHOR]

Star formation is an inefficient process and in general only a small fraction of the gas in a giant molecular cloud (GMC) is turned into stars. This is partly due to the negative effect of stellar feedback from young massive star clusters. Recently, we introduced a novel 1D numerical treatment of the effects of stellar feedback from young massive clusters on their natal clouds, which we named warpfield. Here, we present version 2 of the warpfield code, containing improved treatments of the thermal evolution of the gas and the fragmentation of the feedback-driven shell. As part of this update, we have produced new cooling and heating tables that account for the combined effects of photoionization and collisional ionization on the cooling rate, which we now make publicly available. We employ our updated version of warpfield to investigate the impact of stellar feedback on GMCs with a broad range of masses and surface densities and a variety of density profiles. We show that the minimum star formation efficiency (SFE) εmin, i.e. the SFE above which the cloud is destroyed by feedback, is mainly set by the average cloud surface density. An SFE of 1–6  per cent is generally sufficient to destroy a GMC. We also find an SFE per free-fall time εff ∼ 0.3  per cent, in good agreement with recent observations. Our results imply that feedback alone is sufficient to explain the low observed SFE of GMCs. Finally, we show that very massive clouds with steep density profiles – possible proxies of the giant clumps observed in galaxies at |$z$| ≈ 2 – are more resilient to feedback than typical GMCs, with εmin between 1 and 12  per cent. [ABSTRACT FROM AUTHOR]

Molecular clouds are turbulent structures whose star formation efficiency (SFE) is strongly affected by internal stellar feedback processes. In this paper, we determine how sensitive the SFE of molecular clouds is to randomized inputs in the star formation feedback loop, and to what extent relationships between emergent cloud properties and the SFE can be recovered. We introduce the yule suite of 26 radiative magnetohydrodynamic simulations of a 10 000 solar mass cloud similar to those in the solar neighbourhood. We use the same initial global properties in every simulation but vary the initial mass function sampling and initial cloud velocity structure. The final SFE lies between 6 and 23 per cent when either of these parameters are changed. We use Bayesian mixed-effects models to uncover trends in the SFE. The number of photons emitted early in the cluster's life and the length of the cloud provide the strongest predictors of the SFE. The H ii regions evolve following an analytic model of expansion into a roughly isothermal density field. The more efficient feedback is at evaporating the cloud, the less the star cluster is dispersed. We argue that this is because if the gas is evaporated slowly, the stars are dragged outwards towards surviving gas clumps due to the gravitational attraction between the stars and gas. While star formation and feedback efficiencies are dependent on non-linear processes, statistical models describing cloud-scale processes can be constructed. [ABSTRACT FROM AUTHOR]

We use the polaris radiative transfer code to generate predictions of the two main observables directly sensitive to the magnetic field morphology and strength in filaments: dust polarization and gas Zeeman line splitting. We simulate generic gas filaments with power-law density profiles assuming two density-field strength dependencies, six different filament inclinations, and nine distinct magnetic field morphologies, including helical, toroidal, and warped magnetic field geometries. We present idealized spatially resolved dust polarization and Zeeman-derived field strengths and directions maps. Under the assumption that dust grains are aligned by radiative torques, dust polarization traces the projected plane-of-the-sky magnetic field morphology. Zeeman line splitting delivers simultaneously the intensity-weighted line-of-sight field strength and direction. We show that linear dust polarization alone is unable to uniquely constrain the 3D field morphology. We demonstrate that these ambiguities are ameliorated or resolved with the addition of the Zeeman directional information. Thus, observations of both the dust polarization and Zeeman splitting together provide the most promising means for obtaining constraints of the 3D magnetic field configuration. We find that the Zeeman-derived field strengths are at least a factor of a few below the input field strengths due to line-of-sight averaging through the filament density gradient. Future observations of both dust polarization and Zeeman splitting are essential for gaining insights into the role of magnetic fields in star- and cluster-forming filaments. [ABSTRACT FROM AUTHOR]

We study the dynamical state of massive cores by using a simple analytical model, an observational sample, and numerical simulations of collapsing massive cores. From the analytical model, we find that cores increase their column density and velocity dispersion as they collapse, resulting in a time evolution path in the Larson velocity dispersion-size diagram from large sizes and small velocity dispersions to small sizes and large velocity dispersions, while they tend to equipartition between gravity and kinetic energy. From the observational sample, we find that: (a) cores with substantially different column densities in the sample do not follow a Larson-like linewidth–size relation. Instead, cores with higher column densities tend to be located in the upper-left corner of the Larson velocity dispersion σv,3D–size R diagram, a result explained in the hierarchical and chaotic collapse scenario. (b) Cores appear to have overvirial values. Finally, our numerical simulations reproduce the behaviour predicted by the analytical model and depicted in the observational sample: collapsing cores evolve towards larger velocity dispersions and smaller sizes as they collapse and increase their column density. More importantly, however, they exhibit overvirial states. This apparent excess is due to the assumption that the gravitational energy is given by the energy of an isolated homogeneous sphere. However, such excess disappears when the gravitational energy is correctly calculated from the actual spatial mass distribution. We conclude that the observed energy budget of cores is consistent with their non-thermal motions being driven by their self-gravity and in the process of dynamical collapse. [ABSTRACT FROM AUTHOR]

The escape of ultraviolet photons from the densest regions of the interstellar medium (ISM) - giant molecular clouds (GMCs) - is a poorly constrained parameter which is vital to understanding the ionization of the ISM and the intergalactic medium. We characterize the escape fraction, fesc,GMC, from a suite of individual GMC simulations with masses in the range 104-6 M☉ using the adaptive-mesh refinement code FLASH. We find significantly different fesc,GMC depending on the GMC mass that can reach >90 per cent in the evolution of 5 × 104 and 105 M☉ clouds or remain low at ∼5 per cent for most of the lifetime of more massive GMCs. All clouds show fluctuations over short, sub-Myr time-scales produced by flickering HII regions. We combine our results to calculate the total escape fraction (fesc,tot) from GMC populations in dwarf starburst and spiral galaxies by randomly drawing clouds from a GMC mass distribution (dN/dM α Mα, where α is either -1.5 or -2.5) over fixed time intervals.We find typical fesc,tot values of 8 per cent for both the dwarf and spiral models. The fluctuations of fesc,tot, however, are much larger for the dwarf models with values as high as 90 per cent. The photons escaping from the 5 × 104 and 105 M☉ GMCs are the dominant contributors to fesc,tot in all cases. We also show that the accompanying star formation rates (SFRs) of our model (∼2 × 10-2 and 0.73 M☉yr-1) are consistent with observations of SFRs in dwarf starburst and spiral galaxies, respectively. [ABSTRACT FROM AUTHOR]

Line emission is strongly dependent on the local environmental conditions in which the emitting tracers reside. In this work, we focus on modelling the CO emission from simulated giant molecular clouds (GMCs), and study the variations in the resulting line ratios arising from the emission from the J = 1-0, J = 2-1, and J = 3-2 transitions. We perform a set of smoothed particle hydrodynamics simulations with time-dependent chemistry, in which environmental conditions - including total cloud mass, density, size, velocity dispersion, metallicity, interstellar radiation field (ISRF), and the cosmic ray ionization rate (CRIR) - were systematically varied. The simulations were then post-processed using radiative transfer to produce synthetic emission maps in the three transitions quoted above. We find that the cloud-averaged values of the line ratios can vary by up to±0.3 dex, triggered by changes in the environmental conditions. Changes in the ISRF and/or in the CRIR have the largest impact on line ratios since they directly affect the abundance, temperature, and distribution of CO-rich gas within the clouds. We show that the standard methods used to convert CO emission to H2 column density can underestimate the total H2 molecular gas in GMCs by factors of 2 or 3, depending on the environmental conditions in the clouds. [ABSTRACT FROM AUTHOR]

We investigate how star formation is spatially organized in the grand-design spiral NGC 1566 from deep Hubble Space Telescope photometry with the Legacy ExtraGalactic UV Survey. Our contour-based clustering analysis reveals 890 distinct stellar conglomerations at various levels of significance. These star-forming complexes are organized in a hierarchical fashion with the larger congregations consisting of smaller structures, which themselves fragment into even smaller and more compact stellar groupings. Their size distribution, covering a wide range in length-scales, shows a power law as expected from scale-free processes. We explain this shape with a simple 'fragmentation and enrichment' model. The hierarchical morphology of the complexes is confirmed by their mass-size relation that can be represented by a power law with a fractional exponent, analogous to that determined for fractal molecular clouds. The surface stellar density distribution of the complexes shows a lognormal shape similar to that for supersonic non-gravitating turbulent gas. Between 50 and 65 per cent of the recently formed stars, as well as about 90 per cent of the young star clusters, are found inside the stellar complexes, located along the spiral arms.We find an age difference between young stars inside the complexes and those in their direct vicinity in the arms of at least 10 Myr. This time-scale may relate to the minimum time for stellar evaporation, although we cannot exclude the in situ formation of stars. As expected, star formation preferentially occurs in spiral arms. Our findings reveal turbulent-driven hierarchical star formation along the arms of a grand-design galaxy. [ABSTRACT FROM AUTHOR]

We present three-dimensional radiation-hydrodynamical simulations of the impact of stellar winds, photoelectric heating, photodissociating and photoionizing radiation, and supernovae on the chemical composition and star formation in a stratified disc model. This is followed by a sink-based model for star clusters with populations of individual massive stars. Stellar winds and ionizing radiation regulate the star formation rate at a factor of ~10 below the simulation with only supernova feedback due to their immediate impact on the ambient interstellar medium after star formation. Ionizing radiation (with winds and supernovae) significantly reduces the ambient densities for most supernova explosions to ρ < 10-25 g cm-3, compared to 10-23g cm-3 for the model with only winds and supernovae. Radiation from massive stars reduces the amount of molecular hydrogen and increases the neutral hydrogen mass and volume filling fraction. Only this model results in a molecular gas depletion time-scale of 2 Gyr and shows the best agreement with observations. In the radiative models, the Hα emission is dominated by radiative recombination as opposed to collisional excitation (the dominant emission in non-radiative models), which only contributes ~1-10 per cent to the total Hα emission. Individual massive stars (M ≥ 30M☉) with short lifetimes are responsible for significant fluctuations in the Hα luminosities. The corresponding inferred star formation rates can underestimate the true instantaneous star formation rate by a factor of ~10. [ABSTRACT FROM AUTHOR]

The carbon monoxide (CO) rotational transition lines are the most common tracers of molecular gas within giant molecular clouds (MCs).We study the ratio (R2-1/1-0) between CO's first two emission lines and examine what information it provides about the physical properties of the cloud. To study R2-1/1-0, we perform smooth particle hydrodynamic simulations with timedependent chemistry (using GADGET-2), along with post-process radiative transfer calculations on an adaptive grid (using RADMC-3D) to create synthetic emission maps of a MC. R2-1/1-0 has a bimodal distribution that is a consequence of the excitation properties of each line, given that J = 1 reaches local thermal equilibrium while J = 2 is still sub-thermally excited in the considered clouds. The bimodality of R2-1/1-0 serves as a tracer of the physical properties of different regions of the cloud, and it helps constrain local temperatures, densities and opacities. Additionally, this bimodal structure shows an important portion of the CO emission comes from diffuse regions of the cloud, suggesting that the commonly used conversion factor of R2-1/1-0 ~ 0.7 between both lines may need to be studied further. [ABSTRACT FROM AUTHOR]

We present a method to simultaneously model the dust far-infrared (FIR) spectral energy distribution (SED) and the total infrared - carbon monoxide (CO) integrated intensity (SIR-ICO) relationship. The modelling employs a hierarchical Bayesian (HB) technique to estimate the dust surface density, temperature (Teff), and spectral index at each pixel from the observed FIR maps. Additionally, given the corresponding CO map, the method simultaneously estimates the slope and intercept between the FIR and CO intensities, which are global properties of the observed source. The model accounts for correlated and uncorrelated uncertainties, such as those present in Herschel observations. Using synthetic data sets, we demonstrate the accuracy of the HB method, and contrast the results with common non-hierarchical fitting methods. As an initial application, we model the dust and gas on 100 pc scales in the Magellanic Clouds from Herschel FIR and NANTEN CO observations. The slopes of the logSIR-logICO relationship are similar in both galaxies, falling in the range 1.1-1.7. However, in the Small Magellanic Cloud the intercept is nearly three times higher, which can be explained by its lower metallicity than the Large Magellanic Cloud (LMC), resulting in a larger SIR per unit ICO. The HB modelling evidences an increase in Teff in regions with the highest ICO in the LMC. This may be due to enhanced dust heating in the densest molecular regions from young stars. Such simultaneous dust and gas modelling may reveal variations in the properties of the interstellar medium and its association with other galactic characteristics, such as star formation rates and/or metallicities. [ABSTRACT FROM AUTHOR]

We report on observations of the hydroxyl radical (OH) within The HI, OH, Recombination line survey (THOR) pilot region. The region is bounded approximately between Galactic coordinates l = 29.°2 to 31.°5 and b = -1. ° 0 to +1.°0 and includes the high-mass star-forming region W43. We identify 103 maser sites, including 72 with 1612 MHz masers, 42 showing masers in either of the main-line transitions at 1665 and 1667MHzand four showing 1720MHz masers. Most maser sites with either main-line or 1720 MHz emission are associated with star formation, whereas most of the 1612 MHz masers are associated with evolved stars. We find that nearly all of the main-line maser sites are co-spatial with an infrared source, detected by GLIMPSE. We alsofind diffuse OH emission, as well as OH in absorption towards selected unresolved or partially resolved sites. Extended OH absorption is found towards the wellknown star-forming complex W43 Main. [ABSTRACT FROM AUTHOR]

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, smallscale 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. [ABSTRACT FROM AUTHOR]

We run numerical simulations of molecular clouds, adopting properties similar to those found in the central molecular zone (CMZ) of the Milky Way. For this, we employ the moving mesh code AREPO and perform simulations which account for a simplified treatment of timedependent chemistry and the non-isothermal nature of gas and dust. We perform simulations using an initial density of n0 = 103 cm-3 and a mass of 1.3 × 105M☉. Furthermore, we vary the virial parameter, defined as the ratio of kinetic and potential energy, α = Ekin/ [ABSTRACT FROM AUTHOR]

It has been shown that the behaviour of primordial gas collapsing in a dark matter minihalo can depend on the adopted choice of three-body H2 formation rate. The uncertainties in this rate span two orders of magnitude in the current literature, and so it remains a source of uncertainty in our knowledge of Population III star formation. Here, we investigate how the amount of fragmentation in primordial gas depends on the adopted three-body rate. We present the results of calculations that follow the chemical and thermal evolution of primordial gas as it collapses in two dark matter minihaloes. Our results on the effect of three-body rate on the evolution until the first protostar forms agree well with previous studies. However, our modified version of GADGET-2 smoothed particle hydrodynamics also includes sink particles, which allows us to follow the initial evolution of the accretion disc that builds up on the centre of each halo, and capture the fragmentation in gas as well as its dependence on the adopted three-body H2 formation rate. We find that the fragmentation behaviour of the gas is only marginally affected by the choice of three-body rate co-efficient, and that halo-to-halo differences are of equal importance in affecting the final mass distribution of stars. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Recent observational and theoretical investigations have emphasized the importance of filamentary networks within molecular clouds as sites of star formation. Since such environments are more complex than those of isolated cores, it is essential to understand how the observed line profiles from collapsing cores with non-spherical geometry are affected by filaments. In this study, we investigate line profile asymmetries by performing radiative transfer calculations on hydrodynamic models of three collapsing cores that are embedded in filaments. We compare the results to those that are expected for isolated cores. We model the five lowest rotational transition line (J = 1–0, 2–1, 3–2, 4–3 and 5–4) of both optically thick (HCN, HCO+) as well as optically thin (N2H+, H13CO+) molecules using constant abundance laws. We find that less than 50 per cent of simulated (1–0) transition lines show blue infall asymmetries due to obscuration by the surrounding filament. However, the fraction of collapsing cores that have a blue asymmetric emission line profile rises to 90 per cent when observed in the (4–3) transition. Since the densest gas towards the collapsing core can excite higher rotational states, upper level transitions are more likely to produce blue asymmetric emission profiles. We conclude that even in irregular, embedded cores one can trace infalling gas motions with blue asymmetric line profiles of optically thick lines by observing higher transitions. The best tracer of collapse motions of our sample is the (4–3) transition of HCN, but the (3–2) and (5–4) transitions of both HCN and HCO+ are also good tracers. [ABSTRACT FROM PUBLISHER]

Low-metallicity star formation poses a central problem of cosmology, as it determines the characteristic mass scale and distribution for the first and second generations of stars forming in our Universe. Here, we present a comprehensive investigation assessing the relative impact of metals and magnetic fields, which may both be present during low-metallicity star formation. We show that the presence of magnetic fields generated via the small-scale dynamo stabilizes the protostellar disc and provides some degree of support against fragmentation. In the absence of magnetic fields, the fragmentation time-scale in our model decreases by a factor of ∼10 at the transition from Z = 0 to Z > 0, with subsequently only a weak dependence on metallicity. Similarly, the accretion time-scale of the cluster is set by the large-scale dynamics rather than the local thermodynamics. In the presence of magnetic fields, the primordial disc can become completely stable, therefore forming only one central fragment. At Z > 0, the number of fragments is somewhat reduced in the presence of magnetic fields, though the shape of the mass spectrum is not strongly affected in the limits of the statistical uncertainties. The fragmentation time-scale, however, increases by roughly a factor of 3 in the presence of magnetic fields. Indeed, our results indicate comparable fragmentation time-scales in primordial runs without magnetic fields and Z > 0 runs with magnetic fields. [ABSTRACT FROM PUBLISHER]

Observations of atomic or molecular lines can provide important information about the physical state of star-forming regions. In order to investigate the line profiles from dynamical collapsing massive star-forming regions (MSFRs), we model the emission from hydrodynamic simulations of a collapsing cloud in the absence of outflows. By performing radiative transfer calculations, we compute the optically thick HCO+ and optically thin N2H+ line profiles from two collapsing regions at different epochs. Due to large-scale collapse, theMSFRs have large velocity gradients, reaching up to 20 km s-1 pc-1 across the central core. The optically thin lines typically contain multiple velocity components resulting from the superposition of numerous density peaks along the line of sight. The optically thick lines are only marginally shifted to the blue side of the optically thin line profiles, and frequently do not have a central depression in their profiles due to self-absorption. As the regions evolve, the lines become brighter and the optically thick lines become broader. The lower-order HCO+ (1-0) transitions are better indicators of collapse than the higher-order (4-3) transitions. We also investigate how the beam sizes affect profile shapes. Smaller beams lead to brighter and narrower lines that are more skewed to the blue in HCO+ relative to the true core velocity, but show multiple components in N2H+. High-resolution observations (e.g., with Atacama Large Millimeter Array) can test these predictions and provide insights into the nature of MSFRs. [ABSTRACT FROM AUTHOR]

We analyse the cosmological simulations performed in the recent work of Greif et al., which followed the early growth and merger history of Population III (Pop III) stars while resolving scales as small as 0.05 R⊙. This is the first set of cosmological simulations to self-consistently resolve the rotation and internal structure of Pop III protostars. We find that Pop III stars form under significant rotational support which is maintained for the duration of the simulations. The protostellar surfaces spin from ∼50 per cent to nearly 100 per cent of Keplerian rotational velocity.These rotation rates persist after experiencing multiple stellar merger events. In the brief time period simulated (∼10 yr), the protostars show little indication of convective instability, and their properties furthermore show little correlation with the properties of their host minihaloes. If Pop III protostars within this range of environments generally form with high degrees of rotational support, and if this rotational support is maintained for a sufficient amount of time, this has a number of crucial implications for Pop III evolution and nucleosynthesis, as well as the possibility for Pop III pair-instability supernovae, and the question of whether the first stars produced gamma-ray bursts. [ABSTRACT FROM AUTHOR]

The first stars to form in the universe are believed to have distribution of masses biased toward massive stars. This contrasts with the present-day initial mass function, which has a predominance of stars with masses lower than 1M☉. Therefore, the mode of star formation must have changed as the universe evolved. Such a transition is attributed to a more efficient cooling provided by increasing metallicity. Especially dust cooling can overcome the compressional heating, which lowers the gas temperature thus increasing its instability to fragmentation. The purpose of this paper is to verify if dust cooling can efficiently cool the gas, and enhance the fragmentation of gas clouds at the early stages of the universe. To confirm that, we calculate a set of hydrodynamic simulations that include sink particles, which represent contracting protostars. The thermal evolution of the gas during the collapse is followed by making use of a primordial chemical network and also a recipe for dust cooling. We model four clouds with different amounts of metals (10-4, 10-5, 10-6 Z☉, and 0), and analyze how this property affect the fragmentation of star-forming clouds. We find evidence for fragmentation in all four cases, and hence conclude that there is no critical metallicity below which fragmentation is impossible. Nevertheless, there is a clear change in the behavior of the clouds at Z ≲ 10-5 Z☉, caused by the fact that at this metallicity, fragmentation takes longer to occur than accretion, leading to a flat mass function at lower metallicities. [ABSTRACT FROM AUTHOR]

The first galaxies form due to gravitational collapse of primordial halos. During this collapse, weak magnetic seed fields get amplified exponentially by the small-scale dynamo - a process converting kinetic energy from turbulence into magnetic energy. We use the Kazantsev theory, which describes the small-scale dynamo analytically, to study magnetic field amplification for different turbulent velocity correlation functions. For incompressible turbulence (Kolmogorov turbulence), we find that the growth rate is proportional to the square root of the hydrodynamic Reynolds number, Re1/2. In the case of highly compressible turbulence (Burgers turbulence) the growth rate increases proportional to Re1/3. With a detailed chemical network we are able to follow the chemical evolution and determine the kinetic and magnetic viscosities (due to Ohmic and ambipolar diffusion) during the collapse of the halo. This way, we can calculate the growth rate of the small-scale dynamo quantitatively and predict the evolution of the small-scale magnetic field. As the magnetic energy is transported to larger scales on the local eddy-timescale, we obtain an estimate for the magnetic field on the Jeans scale. Even there, we find that equipartition with the kinetic energy is reached on small timescales. Dynamically relevant field structures can thus be expected already during the formation of the first objects in the Universe. [ABSTRACT FROM PUBLISHER]

ABSTRACT The linewidth-size (σ- R) relationship of the interstellar medium (ISM) has been extensively measured and analysed, in both the local ISM and in nearby normal galaxies. Generally, a power law describes the relationship well with an index ranging from 0.2 to 0.6, and is now referred to as one of the 'Larson's Relationships'. The nature of turbulence and star formation is considered to be intimately related to these relationships, so evaluating the σ- R correlations in various environments is important for developing a comprehensive understanding of the ISM. We measure the linewidth-size relationship in the Central Molecular Zone (CMZ) of the Galactic Centre using spectral line observations of the high-density tracers N2H+, HCN, H13CN and HCO+. We construct dendrograms, which map the hierarchical nature of the position-position-velocity (PPV) data, and compute the linewidths and sizes of the dendrogram-defined structures. The dispersions range from ∼2 to 30 km s−1 in structures spanning sizes 2-40 pc. By performing Bayesian inference, we show that a power law with exponent 0.3-1.1 can reasonably describe the σ- R trend. We demonstrate that the derived σ- R relationship is independent of the locations in the PPV data set where σ and R are measured. The uniformity in the σ- R relationship indicates that turbulence in the CMZ is driven on the large scales beyond ≳30 pc. We compare the CMZ σ- R relationship to that measured in the Galactic molecular cloud Perseus. The exponents between the two systems are similar, suggestive of a connection between the turbulent properties within a cloud to its ambient medium. Yet, the velocity dispersion in the CMZ is systematically higher, resulting in a scaling coefficient that is approximately five times larger. The systematic enhancement of turbulent velocities may be due to the combined effects of increased star formation activity, larger densities and higher pressures relative to the local ISM. [ABSTRACT FROM AUTHOR]

ABSTRACT 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 interstellar medium. We adopt initial conditions characteristic of the warm neutral medium and study two different flow velocities - a slow flow with vflow = 6.8 km s−1 and a fast flow with vflow = 13.6 km s−1. 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 pre-stellar 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. [ABSTRACT FROM AUTHOR]

ABSTRACT We combine the output of hydrodynamical simulations of Population III star cluster formation with stellar evolution models, and calculate the evolution of protostars experiencing variable mass accretion rates due to interactions within a massive disc. We find that the primordial protostars are extended 'fluffy' objects for the bulk of their pre-main-sequence lifetimes. Accretion luminosity feedback from such objects is high, but, as shown in previous work, it has a minimal effect on the star cluster. The extended radii of the protostars, combined with the observation of close encounters in the simulations, suggest that mergers will occur in such systems. Furthermore, mass transfer between close protostellar binaries with extended radii could lead to massive tight binaries, which are a possible progenitor of gamma-ray bursts. [ABSTRACT FROM AUTHOR]

ABSTRACT We investigate the formation of the first stars at the end of the cosmic dark ages with a suite of three-dimensional, moving-mesh simulations that directly resolve the collapse of the gas beyond the formation of the first protostar at the centre of a dark matter minihalo. The simulations cover more than 25 orders of magnitude in density and have a maximum spatial resolution of 0.05 R⊙, which extends well below the radius of individual protostars and captures their interaction with the surrounding gas. In analogy to previous studies that employed sink particles, we find that the Keplerian disc around the primary protostar fragments into a number of secondary protostars, which is facilitated by H2 collisional dissociation cooling and collision-induced emission. The further evolution of the protostellar system is characterized by strong gravitational torques that transfer angular momentum between the secondary protostars formed in the disc and the surrounding gas. This leads to the migration of about half of the secondary protostars to the centre of the cloud in a free-fall time, where they merge with the primary protostar and enhance its growth to about five times the mass of the second most massive protostar. By the same token, a fraction of the protostars obtain angular momentum from other protostars via N-body interactions and migrate to higher orbits. On average, only every third protostar survives until the end of the simulation. However, the number of protostars present at any given time increases monotonically, suggesting that the system will continue to grow beyond the limited period of time simulated here. [ABSTRACT FROM AUTHOR]

ABSTRACT We study the influence of initial conditions on the magnetic field amplification during the collapse of a magnetized gas cloud. We focus on the dependence of the growth and saturation level of the dynamo-generated field on the turbulent properties of the collapsing cloud. In particular, we explore the effect of varying the initial strength and injection scale of turbulence and the initial uniform rotation of the collapsing magnetized cloud. In order to follow the evolution of the magnetic field in both the kinematic and the non-linear regime, we choose an initial field strength of with the magnetic to kinetic energy ratio, Em/ Ek∼ 10−4. Both gravitational compression and the small-scale dynamo initially amplify the magnetic field. Further into the evolution, the dynamo-generated magnetic field saturates but the total magnetic field continues to grow because of compression. The saturation of the small-scale dynamo is marked by a change in the slope of B/ρ2/3 and by a shift in the peak of the magnetic energy spectrum from small scales to larger scales. For the range of initial Mach numbers explored in this study, the dynamo growth rate increases as the Mach number increases from vrms/ cs∼ 0.2 to 0.4 and then starts decreasing from vrms/ cs∼ 1.0. We obtain saturation values of Em/ Ek= 0.2-0.3 for these runs. Simulations with different initial injection scales of turbulence also show saturation at similar levels. For runs with different initial rotation of the cloud, the magnetic energy saturates at Em/ Ek∼ 0.2-0.4 of the equipartition value. The overall saturation level of the magnetic energy, obtained by varying the initial conditions, is in agreement with previous analytical and numerical studies of small-scale dynamo action where turbulence is driven by an external forcing instead of gravitational collapse. [ABSTRACT FROM AUTHOR]

Observations are revealing the ubiquity of filamentary structures in molecular clouds. As cores are often embedded in filaments, it is important to understand how line profiles from such systems differ from those of isolated cores. We perform radiative transfer calculations on a hydrodynamic simulation of a molecular cloud in order to model line emission from collapsing cores embedded in filaments. We model two optically thick lines, CS(2-1) and HCN(1-0), and one optically thin line, N2H+ (1-0), from three embedded cores. In the hydrodynamic simulation, gas self-gravity, turbulence, and bulk flows create filamentary regions within which cores form. Though the filaments have large dispersions, the N2H+ (1-0) lines indicate subsonic velocities within the cores. We find that the observed optically thick line profiles of CS(2-1) and HCN(1-0) vary drastically with viewing angle. In over 50% of viewing angles, there is no sign of a blue asymmetry, an idealized signature of infall motions in an isolated spherical collapsing core. Profiles that primarily trace the cores, with little contribution from the surrounding filament, are characterized by a systematically higher HCN(1-0) peak intensity. The N2H+ (1-0) lines do not follow this trend. We demonstrate that red asymmetric profiles are also feasible in the optically thick lines, due to emission from the filament or one-sided accretion flows onto the core. We conclude that embedded cores may frequently undergo collapse without showing a blue asymmetric profile, and that observational surveys including filamentary regions may underestimate the number of collapsing cores if based solely on profile shapes of optically thick lines. [ABSTRACT FROM AUTHOR]

ABSTRACT The formation of stars is a key process in the early universe with far-reaching consequences for further cosmic evolution. While stars forming from truly primordial gas are thought to be considerably more massive than our Sun, stars in the universe today have typical masses below one solar mass. The physical origin of this transition and the conditions under which it occurs are highly debated. There are two competing models: one based on metal-line cooling as the primary agent and one based on dust cooling. The recent discovery of the extremely metal poor star SDSS J1029151+172927 provides a unique opportunity to distinguish between these two models. Based on simple thermodynamic considerations we argue that SDSS J1029151+172927 was more likely formed as a result of dust continuum cooling rather than cooling by metal lines. We conclude that the masses of extremely metal poor stars are determined by dust-induced fragmentation. [ABSTRACT FROM AUTHOR]

The cosmic dark ages ended a few hundred million years after the big bang, when the first stars began to fill the universe with new light. It has generally been argued that these stars formed in isolation and were extremely massive--perhaps 100 times as massive as the Sun. In a recent study, Clark and collaborators showed that this picture requires revision. They demonstrated that the accretion disks that build up around Population III stars are strongly susceptible to fragmentation and that the first stars should therefore form in clusters rather than in isolation. We here use a series of high-resolution hydrodynamical simulations performed with the moving mesh code AREPO to follow up on this proposal and to study the influence of environmental parameters on the level of fragmentation. We model the collapse of five independent minihalos from cosmological initial conditions, through the runaway condensation of their central gas clouds, to the formation of the first protostar, and beyond for a further 1000 years. During this latter accretion phase, we represent the optically thick regions of protostars by sink particles. Gas accumulates rapidly in the circumstellar disk around the first protostar, fragmenting vigorously to produce a small group of protostars. After an initial burst, gravitational instability recurs periodically, forming additional protostars with masses ranging from ~0.1 to 10 M⊙. Although the shape, multiplicity, and normalization of the protostellar mass function depend on the details of the sink-particle algorithm, fragmentation into protostars with diverse masses occurs in all cases, confirming earlier reports of Population III stars forming in clusters. Depending on the efficiency of later accretion and merging, Population III stars may enter the main sequence in clusters and with much more diverse masses than are commonly assumed. [ABSTRACT FROM AUTHOR]

It has recently been demonstrated that coherent relative streaming velocities of order 30km s-1 between dark matter and gas permeated the universe on scales below a few Mpc directly after recombination. Here we use a series of high-resolution moving-mesh calculations to show that these supersonic motions significantly influence the virialization of the gas in minihalos and delay the formation of the first stars. As the gas streams into minihalos with bulk velocities around 1 km s-1 at z ~ 20, the additional momentum and energy input reduces the gas fractions and central densities of the halos, increasing the typical virial mass required for efficient cooling by a factor of three and delaying Population III star formation by Δz ≃ 4. Since the distribution of the magnitude of the streaming velocities is narrowly peaked around a non-negligible value, this effect is important in most regions of the universe. As a consequence, the increased minimum halo mass implies a reduction of the absolute number of minihalos that can be expected to cool and form Population III stars by up to an order of magnitude. We further find that the streaming velocities increase the turbulent velocity dispersion of the minihalo gas, which could affect its ability to fragment and hence alter the mass function of the first stars. [ABSTRACT FROM AUTHOR]

We use three-dimensional smoothed particle hydrodynamics simulations together with a dynamical ray-tracing scheme to investigate the build-up of the first H ii regions around massive Population III stars in minihaloes. We trace the highly anisotropic breakout of the ionizing radiation into the intergalactic medium, allowing us to predict the resulting recombination radiation with greatly increased realism. Our simulations, together with Press–Schechter type arguments, allow us to predict the Population III contribution to the radio background at via bremsstrahlung and 21-cm emission. We find a global bremsstrahlung signal of around , and a combined 21-cm signature which is an order of magnitude larger. Both might be within the reach of the planned Square Kilometer Array experiment, although detection of the free–free emission is only marginal. The imprint of the first stars on the cosmic radio background might provide us with one of the few diagnostics to test the otherwise elusive minihalo star formation site. [ABSTRACT FROM AUTHOR]

We introduce a simple and efficient algorithm for diffusion in smoothed particle hydrodynamics (SPH) simulations and apply it to the problem of chemical mixing. Based on the concept of turbulent diffusion, we link the diffusivity of a pollutant to the local physical conditions and can thus resolve mixing in space and time. We apply our prescription to the evolution of an idealized supernova remnant and find that we can model the distribution of heavy elements without having to explicitly resolve hydrodynamic instabilities in the post-shock gas. Instead, the dispersal of the pollutant is implicitly modelled through its dependence on the local velocity dispersion. Our method can thus be used in any SPH simulation that investigates chemical mixing but lacks the necessary resolution on small scales. Potential applications include the enrichment of the interstellar medium in present-day galaxies, as well as the intergalactic medium at high redshifts. [ABSTRACT FROM AUTHOR]

We investigate the properties of the first galaxies at with highly resolved numerical simulations, starting from cosmological initial conditions and taking into account all relevant primordial chemistry and cooling. A first galaxy is characterized by the onset of atomic hydrogen cooling, once the virial temperature exceeds , and its ability to retain photoheated gas. We follow the complex accretion and star formation history of a system by means of a detailed merger tree and derive an upper limit on the number of Population III (Pop III) stars formed prior to its assembly. We investigate the thermal and chemical evolution of infalling gas and find that partial ionization at temperatures catalyses the formation of and hydrogen deuteride, allowing the gas to cool to the temperature of the cosmic microwave background. Depending on the strength of radiative and chemical feedback, primordial star formation might be dominated by intermediate-mass Pop III stars formed during the assembly of the first galaxies. Accretion on to the nascent galaxy begins with hot accretion, where gas is accreted directly from the intergalactic medium and shock heated to the virial temperature, but is quickly accompanied by a phase of cold accretion, where the gas cools in filaments before flowing into the parent halo with high velocities. The latter drives supersonic turbulence at the centre of the galaxy and could lead to very efficient chemical mixing. The onset of turbulence in the first galaxies thus likely marks the transition to Pop II star formation. [ABSTRACT FROM AUTHOR]

We investigate how the dynamical state of a turbulently supported, , molecular cloud affects the properties of the cluster it forms, focusing our discussion on the star formation efficiency (SFE) and the initial mass function (IMF). A variety of initial energy states are examined in this paper, ranging from clouds with to clouds with , and for both isothermal and piece-wise polytropic equations of state (similar to that suggested by Larson). It is found that arbitrary SFEs are possible, with strongly unbound clouds yielding very low SFEs. We suggest that the low SFE in the Maddelena cloud may be a consequence of the relatively unbound state of its internal structure. It is also found that competitive accretion results in the observed IMF when the clouds have initial energy states of . We show that under such conditions the shape of the IMF is independent of time in the calculations. This demonstrates that the global accretion process can be terminated at any stage in the cluster's evolution, while still yielding a distribution of stellar masses that is consistent with the observed IMF. As the clouds become progressively more unbound, competitive accretion is less important and the protostellar mass function flattens. These results predict that molecular clouds should be permeated with a distributed population of stars that follow a flatter than Salpeter IMF. [ABSTRACT FROM AUTHOR]

We present results of hydrodynamic simulations of star formation triggered by cloud-cloud collisions. During the early stages of star formation, low-mass objects form by gravitational instabilities in protostellar discs. A number of these low-mass objects are in the sub-stellar mass range, including a few objects of planetary mass. The disc instabilities that lead to the formation of low-mass objects in our simulations are the product of disc-disc interactions and/or interactions between the discs and their surrounding gas. [ABSTRACT FROM PUBLISHER]

Recent studies of dense clumps/cores in a number of regions of low-mass star formation have shown that the mass distribution of these clumps closely resembles the initial mass function (IMF) of field stars. One possible interpretation of these observations is that we are witnessing the fragmentation of the clouds into the IMF, and the observed clumps are bound pre-stellar cores. In this paper, we highlight a potential difficulty in this interpretation, namely that clumps of varying mass are likely to have systematically varying lifetimes. This ‘time-scale’ problem can effectively destroy the similarity between the clump and stellar mass functions, such that a stellar-like clump mass function results in a much steeper stellar IMF. We also discuss some ways in which this problem may be avoided. [ABSTRACT FROM AUTHOR]