Context. Even though turbulent motions are found everywhere in astrophysical systems, the origin of this turbulence is poorly understood. When cosmic structures form, they grow in mass via accretion ...from their surrounding environment. Aims. We propose that accretion is able to drive internal turbulent motions in a wide range of astrophysical objects and study this process in the case of galaxies, molecular clouds, and protoplanetary disks. Methods. We use a combination of numerical simulations and analytical arguments to predict the level of turbulence as a function of the accretion rate, the dissipation scale, and the density contrast, and compare our models with observational data. Results. We find that in Milky Way type galaxies the observed level of turbulence in the interstellar medium can be explained by accretion, provided that the galaxies gain mass at a rate comparable to the rate at which they form stars. This process is particularly relevant in the extended outer disks beyond the star-forming radius. For it to drive turbulence in dwarf galaxies, the accretion rate needs to exceed the star formation rate by a large factor, so we expect other sources to dominate. We also calculate the rate at which molecular clouds grow in mass when they build up from the atomic component of the galactic gas and find that their internal turbulence is likely to be driven by accretion as well. It is the very process of cloud formation that excites turbulent motions on small scales by establishing the turbulent cascade. In the case of T Tauri disks, we show that accretion can drive subsonic turbulence if the rate at which gas falls onto the disk is comparable to the rate at which disk material accretes onto the central star. This also explains the observed relation of accretion rate and stellar mass, $\dot{M}$ $\propto$ $M_\star^{1.8}$. The efficiency required to convert infall motion into turbulence is a few percent in all three cases. Conclusions. We conclude that accretion-driven turbulence is a universal concept with far-reaching implications for a wide range of astrophysical objects.
The SILCC (SImulating the Life-Cycle of molecular Clouds) project aims to self-consistently understand the small-scale structure of the interstellar medium (ISM) and its link to galaxy evolution. We ...simulate the evolution of the multiphase ISM in a (500 pc)2 × ±5 kpc region of a galactic disc, with a gas surface density of
$\Sigma _{_{\rm GAS}} = 10 \;{\rm M}_{\odot }\,{\rm pc}^{-2}$
. The flash 4 simulations include an external potential, self-gravity, magnetic fields, heating and radiative cooling, time-dependent chemistry of H2 and CO considering (self-) shielding, and supernova (SN) feedback but omit shear due to galactic rotation. We explore SN explosions at different rates in high-density regions (peak), in random locations with a Gaussian distribution in the vertical direction (random), in a combination of both (mixed), or clustered in space and time (clus/clus2). Only models with self-gravity and a significant fraction of SNe that explode in low-density gas are in agreement with observations. Without self-gravity and in models with peak driving the formation of H2 is strongly suppressed. For decreasing SN rates, the H2 mass fraction increases significantly from <10 per cent for high SN rates, i.e. 0.5 dex above Kennicutt–Schmidt, to 70–85 per cent for low SN rates, i.e. 0.5 dex below KS. For an intermediate SN rate, clustered driving results in slightly more H2 than random driving due to the more coherent compression of the gas in larger bubbles. Magnetic fields have little impact on the final disc structure but affect the dense gas (n ≳ 10 cm−3) and delay H2 formation. Most of the volume is filled with hot gas (∼80 per cent within ±150 pc). For all but peak driving a vertically expanding warm component of atomic hydrogen indicates a fountain flow. We highlight that individual chemical species populate different ISM phases and cannot be accurately modelled with temperature-/density-based phase cut-offs.
Abstract
We present 3D ‘zoom-in’ simulations of the formation of two molecular clouds out of the galactic interstellar medium. We model the clouds – identified from the SILCC simulations – with a ...resolution of up to 0.06 pc using adaptive mesh refinement in combination with a chemical network to follow heating, cooling and the formation of H2 and CO including (self-) shielding. The two clouds are assembled within a few million years with mass growth rates of up to ∼10−2 M⊙ yr−1 and final masses of ∼50 000 M⊙. A spatial resolution of ≲0.1 pc is required for convergence with respect to the mass, velocity dispersion and chemical abundances of the clouds, although these properties also depend on the cloud definition such as based on density thresholds, H2 or CO mass fraction. To avoid grid artefacts, the progressive increase of resolution has to occur within the free-fall time of the densest structures (1–1.5 Myr) and ≳200 time-steps should be spent on each refinement level before the resolution is progressively increased further. This avoids the formation of spurious, large-scale, rotating clumps from unresolved turbulent flows. While CO is a good tracer for the evolution of dense gas with number densities n ≥ 300 cm−3, H2 is also found for n ≲ 30 cm−3 due to turbulent mixing and becomes dominant at column densities around 30–50 M⊙ pc−2. The CO-to-H2 ratio steadily increases within the first 2 Myr, whereas XCO ≃ 1–4 × 1020 cm−2 (K km s−1)−1 is approximately constant since the CO(1−0) line quickly becomes optically thick.
We study the properties of clumps formed in three-dimensional weakly magnetized magnetohydrodynamic simulations of converging flows in the thermally bistable, warm neutral medium (WNM). We find the ...following. (1) Similarly to the situation in the classical two-phase medium, cold, dense clumps form through dynamically triggered thermal instability in the compressed layer between the convergent flows, and are often characterized by a sharp density jump at their boundaries though not always. (2) However, the clumps are bounded by phase-transition fronts rather than by contact discontinuities, and thus they grow in size and mass mainly by accretion of WNM material through their boundaries. (3) The clump boundaries generally consist of thin layers of thermally unstable gas, but these layers are often widened by the turbulence, and penetrate deep into the clumps. (4) The clumps are approximately in both ram and thermal pressure balance with their surroundings, a condition which causes their internal Mach numbers to be comparable to the bulk Mach number of the colliding WNM flows. (5) The clumps typically have mean temperatures 20 ≲〈T〉≲ 50 K, corresponding to the wide range of densities they contain (20 ≲n≲ 5000 cm−3) under a nearly isothermal equation of state. (6) The turbulent ram pressure fluctuations of the WNM induce density fluctuations that then serve as seeds for local gravitational collapse within the clumps. (7) The velocity and magnetic fields tend to be aligned with each other within the clumps, although both are significantly fluctuating, suggesting that the velocity tends to stretch and align the magnetic field with it. (8) The typical mean field strength in the clumps is a few times larger than that in the WNM. (9) The magnetic field strength in the densest regions within the clumps (n∼ 104 cm−3) has a mean value of B∼ 6 μG but with a large scatter of nearly two orders of magnitude, implying that both sub- and supercritical cores are formed in the simulation. (10) In the final stages of the evolution, the clumps' growth drives them into gravitational instability, at which point star formation sets in, and the pressure in the clumps' centres increases even further.
It is widely accepted that supersonic, magnetized turbulence plays a fundamental role for star formation in molecular clouds. It produces the initial dense gas seeds out of which new stars can form. ...However, the exact relation between gas compression, turbulent Mach number and magnetic field strength is still poorly understood. Here, we introduce and test an analytical prediction for the relation between the density variance and the rms Mach number
in supersonic, isothermal, magnetized turbulent flows. We approximate the density and velocity structure of the interstellar medium as a superposition of shock waves. We obtain the density contrast considering the momentum equation for a single magnetized shock and extrapolate this result to the entire cloud. Depending on the field geometry, we then make three different assumptions based on observational and theoretical constraints: B independent of ρ, B∝ρ1/2 and B∝ρ. We test the analytically derived density variance-Mach number relation with numerical simulations, and find that for B∝ρ1/2, the variance in the logarithmic density contrast,
, fits very well to simulated data with turbulent forcing parameter b= 0.4, when the gas is super-Alfvénic. However, this result breaks down when the turbulence becomes trans-Alfvénic or sub-Alfvénic, because in this regime the turbulence becomes highly anisotropic. Our density variance-Mach number relations simplify to the purely hydrodynamic relation as the ratio of thermal to magnetic pressure β0→∞.
We present collapse simulations of strongly magnetized, turbulent molecular cloud cores with masses ranging from 2.6 to 1000 M in order to study the influence of the initial conditions on the ...turbulence-induced disc formation mechanism proposed recently by Seifried et al. We find that Keplerian discs are formed in all cases independently of the core mass, the strength of turbulence or the presence of global rotation. The discs appear within a few kyr after the formation of the protostar, are 50-150 au in size, and have masses between 0.05 and a few 0.1 M. During the formation of the discs the mass-to-flux ratio stays well below the critical value of 10 for Keplerian disc formation. Hence, flux-loss alone cannot explain the formation of Keplerian discs. The formation of rotationally supported discs at such early phases is rather due to the disordered magnetic field structure and due to turbulent motions in the surroundings of the discs, two effects lowering the classical magnetic braking efficiency. Binary systems occurring in the discs are mainly formed via the disc capturing mechanism rather than via disc fragmentation, which is largely suppressed by the presence of magnetic fields.
We review the current status of knowledge concerning the early phases of star formation during cosmic dawn. This includes the first generations of stars forming in the lowest mass dark matter halos ...in which cooling and condensation of gas with primordial composition is possible at very high redshift (
z
>
20
), namely metal-free Population III stars, and the first generation of massive black holes forming at such early epochs, the so-called black hole seeds. The formation of black hole seeds as end states of the collapse of Population III stars, or via direct collapse scenarios, is discussed. In particular, special emphasis is given to the physics of supermassive stars as potential precursors of direct collapse black holes, in light of recent results of stellar evolution models, and of numerical simulations of the early stages of galaxy formation. Furthermore, we discuss the role of the cosmic radiation produced by the early generation of stars and black holes at high redshift in the process of reionization.
Context. It has been proposed that giant molecular complexes form at the sites of streams of diffuse warm atomic gas that collide at transonic velocities. Aims. We study the global statistics of ...molecular clouds formed by large scale colliding flows of warm neutral atomic interstellar gas under pure hydrodynamic and ideal MHD conditions. The flows deliver material as well as kinetic energy and trigger thermal instability leading eventually to gravitational collapse. Methods. We perform adaptive mesh refinement MHD simulations that, for the first time in this context, treat cooling and self-gravity self-consistently. Results. The clouds formed in the simulations develop a highly inhomogeneous density and temperature structure, with cold dense filaments and clumps condensing from converging flows of warm atomic gas. In the clouds, the column density probability density distribution (PDF) peaks at ~ 2 $\times$ 1021 cm-2 and decays rapidly at higher values; the magnetic intensity correlates weakly with density between n ~ 0.1 and 104 cm-3, and then varies roughly as $n^{1/2}$ for higher densities. Conclusions. The global statistical properties of such molecular clouds are reasonably consistent with observational measurements. Our numerical simulations suggest that molecular clouds form by the moderately supersonic collision of warm atomic gas streams.
ABSTRACT
We present synthetic dust polarization maps of simulated molecular clouds with the goal to systematically explore the origin of the relative orientation of the magnetic field ($\mathbf {B}$) ...with respect to the cloud sub-structure identified in density (n; 3D) and column density (N; 2D). The polarization maps are generated with the radiative transfer code polaris, which includes self-consistently calculated efficiencies for radiative torque alignment. The molecular clouds are formed in two sets of 3D magnetohydrodynamical simulations: (i) in colliding flows (CF), and (ii) in the SILCC-Zoom simulations. In 3D, for the CF simulations with an initial field strength below ∼5 μG, $\mathbf {B}$ is oriented either parallel or randomly with respect to the n-structures. For CF runs with stronger initial fields as well as all SILCC-Zoom simulations, which have an initial field strength of 3 μG, a flip from parallel to perpendicular orientation occurs at high densities of $n_\rm {trans}\, \simeq$ 102–103 cm−3. We suggest that this flip happens if the cloud’s mass-to-flux ratio, μ, is close to or below the critical value of 1. This corresponds to a field strength around 3–5 μG, close to the Galactic average. In 2D, we use the method of Projected Rayleigh Statistics (PRS) to study the relative orientation of $\mathbf {B}$. If present, the flip in orientation occurs in the projected maps at $N_\rm {trans}\, \simeq$ 1021 − 21.5 cm−2. This value is similar to the observed transition value from sub- to supercritical magnetic fields in the interstellar medium. However, projection effects can strongly reduce the predictive power of the PRS method: Depending on the considered cloud or line-of-sight, the projected maps of the SILCC-Zoom simulations do not always show the flip, although it is expected given the 3D morphology. Such projection effects can explain the variety of recently observed field configurations, in particular within a single cloud. Finally, we do not find a correlation between the observed orientation of $\mathbf {B}$ and the N-PDF.
Numerical simulations of outflows formed during the collapse of 100-M⊙ cloud cores are presented. We derive a generalized criterion from magnetohydrodynamical wind theory to analyse the launching ...mechanism of these outflows. The criterion is successfully applied to the whole outflow structure and cases with sub-Keplerian disc rotation. It allows us to decide whether an outflow is driven centrifugally or by the toroidal magnetic pressure. We show that quantities such as the magnetic field line inclination or the ratio of the toroidal to poloidal magnetic field alone are insufficient to determine the driving mechanism of outflows. By performing 12 simulations with variable initial rotational and magnetic energies, we are able to study the influence of the initial conditions on the properties of outflows and jets around massive protostars in detail. Our simulations reveal a strong effect of the magnetic field strength on the morphology of outflows. In runs with weak fields or high rotational energies, well-collimated, fast jets are observed, whereas for strong fields poorly collimated, low-velocity outflows are found. We show that the occurrence of a fast jet is coupled to the existence of a Keplerian protostellar disc. Despite the very different morphologies, all outflows are launched from the discs by centrifugal acceleration with the toroidal magnetic field increasingly contributing to the gas acceleration further away from the discs. The poor collimation of the outflows in runs with strong magnetic fields is a consequence of the weak hoop stresses. This in turn is caused by the slow build-up of a toroidal magnetic field due to strongly sub-Keplerian disc rotation. The mass and momentum outflow rates are of the order of 10−4 M⊙ yr−1 and 10−4 M⊙ km s−1 yr−1, respectively. The mass ejection/accretion ratios scatter around a mean of 0.3 in accordance with observational and analytical results. Based on our results, we suggest an evolutionary scenario for the earliest stage of massive star formation in which initially poorly collimated outflows develop, which successively get better collimated during their evolution due to the generation of fast jets.