Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userHydrodynamic simulations of mechanical stellar feedback in a molecular cloud formed by thermal instability
70
0
0.0
(
0
)
Added by
Christopher Wareing
Publication date
2017
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
We have used the AMR hydrodynamic code, MG, to perform 3D magnetohydrodynamic simulations with self-gravity of stellar feedback in a sheet-like molecular cloud formed through the action of the thermal instability. We simulate the interaction of the mechanical energy input from a 15 solar mass star and a 40 solar mass star into a 100 pc-diameter 17000 solar mass cloud with a corrugated sheet morphology that in projection appears filamentary. The stellar winds are introduced using appropriate Geneva stellar evolution models. In the 15 solar mass star case, the wind forms a narrow bipolar cavity with minimal effect on the parent cloud. In the 40 solar mass star case, the more powerful stellar wind creates a large cylindrical cavity through the centre of the cloud. After 12.5 Myrs and 4.97 Myrs respectively, the massive stars explode as supernovae (SNe). In the 15 solar mass star case, the SN material and energy is primarily deposited into the molecular cloud surroundings over ~10^5 years before the SN remnant escapes the cloud. In the 40 solar mass star case, a significant fraction of the SN material and energy rapidly escapes the molecular cloud along the wind cavity in a few tens of kiloyears. Both SN events compress the molecular cloud material around them to higher densities (so may trigger further star formation), and strengthen the magnetic field, typically by factors of 2-3 but up to a factor of 10. Our simulations are relevant to observations of bubbles in flattened ring-like molecular clouds and bipolar HII regions.
We present a numerical study of the evolution of molecular clouds, from their formation by converging flows in the warm ISM, to their destruction by the ionizing feedback of the massive stars they form. We improve with respect to our previous simulations by including a different stellar-particle formation algorithm, which allows them to have masses corresponding to single stars rather than to small clusters, and with a mass distribution following a near-Salpeter stellar IMF. We also employ a simplified radiative-transfer algorithm that allows the stellar particles to feed back on the medium at a rate that depends on their mass and the local density. Our results are as follows: a) Contrary to the results from our previous study, where all stellar particles injected energy at a rate corresponding to a star of ~ 10 Msun, the dense gas is now completely evacuated from 10-pc regions around the stars within 10-20 Myr, suggesting that this feat is accomplished essentially by the most massive stars. b) At the scale of the whole numerical simulations, the dense gas mass is reduced by up to an order of magnitude, although star formation (SF) never shuts off completely, indicating that the feedback terminates SF locally, but new SF events continue to occur elesewhere in the clouds. c) The SF efficiency (SFE) is maintained globally at the ~ 10% level, although locally, the cloud with largest degree of focusing of its accretion flow reaches SFE ~ 30%. d) The virial parameter of the clouds approaches unity before the stellar feedback begins to dominate the dynamics, becoming much larger once feedback dominates, suggesting that clouds become unbound as a consequence of the stellar feedback. e) The erosion of the filaments that feed the star-forming clumps produces chains of isolated dense blobs reminiscent of those observed in the vicinity of the dark globule B68.
We present Herschel SPIRE Fourier Transform Spectrometer (FTS) observations of N159W, an active star-forming region in the Large Magellanic Cloud (LMC). In our observations, a number of far-infrared cooling lines including CO(4-3) to CO(12-11), [CI] 609 and 370 micron, and [NII] 205 micron are clearly detected. With an aim of investigating the physical conditions and excitation processes of molecular gas, we first construct CO spectral line energy distributions (SLEDs) on 10 pc scales by combining the FTS CO transitions with ground-based low-J CO data and analyze the observed CO SLEDs using non-LTE radiative transfer models. We find that the CO-traced molecular gas in N159W is warm (kinetic temperature of 153-754 K) and moderately dense (H2 number density of (1.1-4.5)e3 cm-3). To assess the impact of the energetic processes in the interstellar medium on the physical conditions of the CO-emitting gas, we then compare the observed CO line intensities with the models of photodissociation regions (PDRs) and shocks. We first constrain the properties of PDRs by modelling Herschel observations of [OI] 145, [CII] 158, and [CI] 370 micron fine-structure lines and find that the constrained PDR components emit very weak CO emission. X-rays and cosmic-rays are also found to provide a negligible contribution to the CO emission, essentially ruling out ionizing sources (ultraviolet photons, X-rays, and cosmic-rays) as the dominant heating source for CO in N159W. On the other hand, mechanical heating by low-velocity C-type shocks with ~10 km/s appears sufficient enough to reproduce the observed warm CO.
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.
We report on the filaments that develop self-consistently in a new numerical simulation of cloud formation by colliding flows. As in previous studies, the forming cloud begins to undergo gravitational collapse because it rapidly acquires a mass much larger than the average Jeans mass. Thus, the collapse soon becomes nearly pressureless, proceeding along its shortest dimension first. This naturally produces filaments in the cloud, and clumps within the filaments. The filaments are not in equilibrium at any time, but instead are long-lived flow features, through which the gas flows from the cloud to the clumps. The filaments are long-lived because they accrete from their environment while simultaneously accreting onto the clumps within them; they are essentially the locus where the flow changes from accreting in two dimensions to accreting in one dimension. Moreover, the clumps also exhibit a hierarchical nature: the gas in a filament flows onto a main, central clump, but other, smaller-scale clumps form along the infalling gas. Correspondingly, the velocity along the filament exhibits a hierarchy of jumps at the locations of the clumps. Two prominent filaments in the simulation have lengths ~15 pc, and masses ~600 Msun above density n ~ 10^3 cm-3 (~2x10^3 Msun at n > 50 cm-3). The density profile exhibits a central flattened core of size ~0.3 pc and an envelope that decays as r^-2.5, in reasonable agreement with observations. Accretion onto the filament reaches a maximum linear density rate of ~30 Msun Myr^-1 pc^-1.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
