Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userStriations, integrals, hourglasses and collapse -- thermal instability driven magnetic simulations of molecular clouds
56
0
0.0
(
0
)
Added by
Christopher Wareing
Publication date
2020
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
The MHD version of the adaptive mesh refinement (AMR) code, MG, has been employed to study the interaction of thermal instability, magnetic fields and gravity through 3D simulations of the formation of collapsing cold clumps on the scale of a few parsecs, inside a larger molecular cloud. The diffuse atomic initial condition consists of a stationary, thermally unstable, spherical cloud in pressure equilibrium with lower density surroundings and threaded by a uniform magnetic field. This cloud was seeded with 10% density perturbations at the finest initial grid level around n=1.1 cm^{-3} and evolved with self-gravity included from the outset. Several cloud diameters were considered (100 pc, 200 pc and 400 pc) equating to several cloud masses (17,000 Msun, 136,000 Msun and 1.1x10^6 Msun). Low-density magnetic-field-aligned striations were observed as the clouds collapse along the field lines into disc-like structures. The induced flow along field lines leads to oscillations of the sheet about the gravitational minimum and an integral-shaped appearance. When magnetically supercritical, the clouds then collapse and generate hourglass magnetic field configurations with strongly intensified magnetic fields, reproducing observational behaviour. Resimulation of a region of the highest mass cloud at higher resolution forms gravitationally-bound collapsing clumps within the sheet that contain clump-frame supersonic (M~5) and super-Alfvenic (M_A~4) velocities. Observationally realistic density and velocity power spectra of the cloud and densest clump are obtained. Future work will use these realistic initial conditions to study individual star and cluster feedback.
rate research
Read More
This paper describes 3D simulations of the formation of collapsing cold clumps via thermal instability inside a larger cloud complex. The initial condition was a diffuse atomic, stationary, thermally unstable, 200pc diameter spherical cloud in pressure equilibrium with low density surroundings. This was seeded with 10% density perturbations at the finest initial grid level (0.29pc) around n_H = 1.1cm^{-3} and evolved with self-gravity included. No magnetic field was imposed. Resimulations at a higher resolution of a region extracted from this simulation (down to 0.039pc), show that the thermal instability forms sheets, then filaments and finally clumps. The width of the filaments increases over time, in one particular case from 0.26 to 0.56pc. Thereafter clumps with sizes of around 5pc grow at the intersections of filaments. 21 distinct clumps, with properties similar to those observed in molecular clouds, are found by using the FellWalker algorithm to find minima in the gravitational potential. Not all of these are gravitationally bound, but the convergent nature of the flow and increasing central density suggest they are likely to form stars. Further simulation of the most massive clump shows the gravitational collapse to a density >10^6 cm^{-3}. These results provide realistic initial conditions that can be used to study feedback in individual clumps, interacting clumps and the entire molecular cloud complex.
Dust continuum and molecular observations of the low column density parts of molecular clouds have revealed the presence of elongated structures which appear to be well aligned with the magnetic field. These so-called striations are usually assumed to be streams that flow towards or away from denser regions. We perform ideal magnetohydrodynamic (MHD) simulations adopting four models that could account for the formation of such structures. In the first two models striations are created by velocity gradients between ambient, parallel streamlines along magnetic field lines. In the third model striations are formed as a result of a Kelvin-Helmholtz instability perpendicular to field lines. Finally, in the fourth model striations are formed from the nonlinear coupling of MHD waves due to density inhomogeneities. We assess the validity of each scenario by comparing the results from our simulations with previous observational studies and results obtained from the analysis of CO (J = 1 - 0) observations from the Taurus molecular cloud. We find that the first three models cannot reproduce the density contrast and the properties of the spatial power spectrum of a perpendicular cut to the long axes of striations. We conclude that the nonlinear coupling of MHD waves is the most probable formation mechanism of striations.
We have used the AMR hydrodynamic code, MG, to perform 3D hydrodynamic simulations with self-gravity of stellar feedback in a spherical clumpy molecular cloud formed through the action of thermal instability. We simulate the interaction of the mechanical energy input from 15 Msun, 40 Msun, 60 Msun and 120 Msun stars into a 100 pc-diameter 16,500 Msun cloud with a roughly spherical morphology with randomly distributed high density condensations. The stellar winds are introduced using appropriate non-rotating Geneva stellar evolution models. In the 15 Msun star case, the wind has very little effect, spreading around a few neighbouring clumps before becoming overwhelmed by the cloud collapse. In contrast, in the 40 Msun, 60 Msun and 120 Msun star cases, the more powerful stellar winds create large cavities and carve channels through the cloud, breaking out into the surrounding tenuous medium during the wind phase and considerably altering the cloud structure. After 4.97 Myrs, 3.97 Myrs and 3.01 Myrs respectively, the massive stars explode as supernovae (SNe). The wind-sculpted surroundings considerably affect the evolution of these SN events as they both escape the cloud along wind-carved channels and sweep up remaining clumps of cloud/wind material. The `cloud as a coherent structure does not survive the SN from any of these stars, but only in the 120 Msun case is the cold molecular material completely destabilised and returned to the unstable thermal phase. In the 40 Msun and 60 Msun cases, coherent clumps of cold material are ejected from the cloud by the SN, potentially capable of further star formation.
We examine the proposal that the HI high-velocity clouds (HVCs) surrounding the Milky Way and other disc galaxies form by condensation of the hot galactic corona via thermal instability. Under the assumption that the galactic corona is well represented by a non-rotating, stratified atmosphere, we find that for this formation mechanism to work the corona must have an almost perfectly flat entropy profile. In all other cases the growth of thermal perturbations is suppressed by a combination of buoyancy and thermal conduction. Even if the entropy profile were nearly flat, cold clouds with sizes smaller than 10 kpc could form in the corona of the Milky Way only at radii larger than 100 kpc, in contradiction with the determined distances of the largest HVC complexes. Clouds with sizes of a few kpc can form in the inner halo only in low-mass systems. We conclude that unless even slow rotation qualitatively changes the dynamics of a corona, thermal instability is unlikely to be a viable mechanism for formation of cold clouds around disc galaxies.
The magnetic field of molecular clouds (MCs) plays an important role in the process of star formation: it determins the statistical properties of supersonic turbulence that controls the fragmentation of MCs, controls the angular momentum transport during the protostellar collapse, and affects the stability of circumstellar disks. In this work, we focus on the problem of the determination of the magnetic field strength. We review the idea that the MC turbulence is super-Alfv{e}nic, and we argue that MCs are bound to be born super-Alfv{e}nic. We show that this scenario is supported by results from a recent simulation of supernova-driven turbulence on a scale of 250 pc, where the turbulent cascade is resolved on a wide range of scales, including the interior of MCs.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
