skip to content

ERC STARTING GRANT RADFEEDBACK (GRANT NO. 679852)

STAR FORMATION AND FEEDBACK IN MOLECULAR CLOUDS

Prof. Stefanie Walch 

Within the SILCC-Zoom project, we study the early impact of ionizing radiation on forming molecular clouds. In Haid et al. (2019) we present our first sub-parsec resolution radiation-hydrodynamic simulations of two molecular clouds self-consistently forming from a turbulent, multiphase ISM. The clouds have similar initial masses of few 104 M⊙, escape velocities of ∼5 km s−1, and a similar initial energy budget. We follow the formation of star clusters with a sink-based model and the impact of radiation from individual massive stars with the tree-based radiation transfer module TREERAY (Wünsch et al., in prep., and Wünsch et al., 2018, MNRAS, 475, 3393). Photoionizing radiation is coupled to a chemical network to follow gas heating, cooling, and molecule formation and dissociation. For the first 3 Myr of cloud evolution, we find that the overall star formation efficiency is considerably reduced by a factor of ˜4 to global cloud values of <10 percent as the mass accretion of sinks that host massive stars is terminated after ≤1 Myr. Despite the low efficiency, star formation is triggered across the clouds. Therefore, a much larger region of the cloud is affected by radiation and the clouds begin to disperse. The time-scale on which the clouds are dispersed sensitively depends on the cloud sub-structure and in particular on the amount of gas at high visual extinction. The damage of radiation done to the highly shielded cloud is delayed. We also show that the radiation input can sustain the thermal and kinetic energy of the clouds at a constant level. Our results strongly support the importance of ionizing radiation from massive stars for explaining the low-observed star formation efficiency of molecular clouds. 

expand: march_2019.png
Image Credit: Stefan Reissl. Image created with the help of the POLARIS code (Reissl et al., 2018, 2018ascl.soft07001R) based on simulation data from Haid et al. (2019) MNRAS, 482, 4062. The colour composite shows red: Hα, green: C+, blue: CO (1-0).

Image Credit: Stefan Reissl. Image created with the help of the POLARIS code (Reissl et al., 2018, 2018ascl.soft07001R) based on simulation data from Haid et al. (2019) MNRAS, 482, 4062. The colour composite shows red: Hα, green: C+ , blue: CO (1-0).

NON EQUILIBRIUM CHEMISTRY AND DESTRUCTION OF CO BY X-RAY FLARES

Sources of X-rays such as active galaxies and X-ray binaries are often variable by orders of magnitude in luminosity over time scales of years. During these flares and for some time afterwards, the surrounding gas is out of chemical and thermal equilibrium. We introduce a new implementation of X-ray radiative transfer coupled to a time dependent chemical network for use in 3D magnetohydrodynamic simulations. ( J. Mackey, S. Walch et al., 2018)

Column density of H2 (left), CO (middle), and the ratio of N(CO)/N(H2) (right) for a fractal molecular cloud irradiated by an external X-ray radiation field of different strengths (see left labels) for 4 Myr. CO is more effectively destroyed than H2, leading to a decreasing CO-to-H2 ratio with increasing X-ray flux.

BIPOLAR OUTFLOW FROM PROTOSTAR FORMATION

We study the influence of episodic outflow feedback on a turbulent, over critical Bonnor-Ebert sphere with a mass of MBES=2.7M⊙ and a radius of rBES=0.056 pc. The mass resolution is 4.5×10−6M⊙ per SPH particle. 
We follow the collapse of the Bonnor-Ebert sphere until tend=180 kyr. A single protostar with a final mass of M⋆=0.45M⊙ forms at t=90 kyr. This protostar launches an S-shaped chain of outflow bullets, called Herbig-Haro objects. This S-shaped form arises due to the varying orientation of the angular momentum axis of the inner accretion disk (IAD), LIAD, caused by anisotropic accretion of the sink particle from the turbulent envelope and outer disk. (Rohde et al., in prep)

Column density plot of a turbulent Bonnor-Ebert sphere at t120 kyr. The core collapsed and formed a protostar that launches a bipolar outflow. Individual outflow bullets can be seen, they are labelled from A to E, according to the ejection history. Bullet A is the oldest one at low velocity; Bullet B and C have higher velocity and overrun bullet A; Bullet D hits the leading shock front and decelerates; Bullet E is about to hit the leading shock front.

VELOCITY CHANNEL MAP OF C18O(1-0) OBSERVATIONS OF A FRAGMENTING FILAMENT

Recent observations of star-forming filaments show that they are split into smaller velocity-coherent filaments, termed fibres. However, these fibres are defined in position-position-velocity (PPV) space, and so from observations it is unclear exactly how they relate to physically cohesive structures in position-position-position (PPP) space. Using simulations of fragmenting non-equilibrium filaments, we produce synthetic observations of the C18O (1-0) line and analysis the resulting PPV cube in the same manner as observers have previously done. As a result we find numerous fibres. Having access to the PPP data from the simulation we show that the fibres do not come from physically cohesive and distinct structures. Rather we show that the mapping from PPV to PPP (and vice-versa) is complicated by the presence of internal shocks within the filament, driven by the filaments accretion from its surrounding. (Clarke et al., 2018)

A velocity channel map movie of synthetic C18O (1-0) observations of a fragmenting non-equilibrium filament. Numerous sub-structures are apparent, these are similar to the fibres seen in real observations.

SYNTHETIC COLUMN DENSITY MAP BASED ON DUST

We compare synthetic observations of dust emission and several molecular line tracers (for example 12CO, 13CO and HI) using zoom-in simulations of molecular cloud formation. We use the SILCC zoom-in simulations presented in Seifried et al. (2017) and new simulations including magnetic fields presented in Seifried et al. (2018). These simulations connect, for the first time simultaneously, the larger scale environment, chemical network, self-gravity and high spatial resolution. The formation of molecular clouds is simulated together with their surroundings, providing more realistic initial conditions for cloud formation than just using an isolated box. The chemical components of the gas are coupled together and simulated simultaneously with the cloud dynamics. The zoom-in simulations reach a resolution of approximately 0.1 pc, showing also the evolution of filamentary structures within the clouds. These simulations give a realistic starting point for making synthetic observations of cloud evolution from large scale structures to molecular cloud and filament scales, and comparing to real observations using several tracers. This enables us to compare, for example, the differences between cloud structures as revealed by different tracers, and their correlation. As an example, the above figure shows a synthetic column density map based on dust emission. We also study the velocity structure of the cloud based on the synthetic line data.

Synthetic column density map based on dust emission for SILCC zoom-in simulations.