Showing 1–50 of 68 results for author: Seifried, D

HYACINTH: HYdrogen And Carbon chemistry in the INTerstellar medium in Hydro simulations

Authors: Prachi Khatri, Cristiano Porciani, Emilio Romano-Díaz, Daniel Seifried, Alexander Schäbe

Abstract: Aims. We present a new sub-grid model, HYACINTH -- HYdrogen And Carbon chemistry in the INTerstellar medium in Hydro simulations -- for computing the non-equilibrium abundances of ${\rm H_2}$ and its carbon-based tracers, namely ${\rm CO}$, ${\rm C}$, and ${\rm C^+}$, in cosmological simulations of galaxy formation. Methods. The model accounts for the unresolved density structure in simulations us… ▽ More Aims. We present a new sub-grid model, HYACINTH -- HYdrogen And Carbon chemistry in the INTerstellar medium in Hydro simulations -- for computing the non-equilibrium abundances of ${\rm H_2}$ and its carbon-based tracers, namely ${\rm CO}$, ${\rm C}$, and ${\rm C^+}$, in cosmological simulations of galaxy formation. Methods. The model accounts for the unresolved density structure in simulations using a variable probability distribution function of sub-grid densities and a temperature-density relation. Included is a simplified chemical network that has been tailored for hydrogen and carbon chemistry within molecular clouds and easily integrated into large-scale simulations with minimal computational overhead. As an example, we applied HYACINTH to a simulated galaxy at redshift $z \sim 2.5$ in post-processing and compared the resulting abundances with observations. Results. The chemical predictions from HYACINTH are in reasonable agreement with high-resolution molecular-cloud simulations at different metallicities. By post-processing a galaxy simulation with HYACINTH, we reproduced the $\rm H\,I-{\rm H_2}$ transition as a function of the hydrogen column density $N_{\rm H}$ for both Milky-Way-like and Large-Magellanic-Cloud-like conditions. Column density maps reveal that ${\rm CO}$ is concentrated in the peaks of the ${\rm H_2}$ distribution, while atomic carbon more broadly traces the bulk of ${\rm H_2}$ in our post-processed galaxy. Based on both the column density maps and the surface density profiles of the different gas species in the post-processed galaxy, we find that ${\rm C^+}$ maintains a substantially high surface density out to $\sim 10 \, \rm kpc$ as opposed to other components that exhibit a higher central concentration. This is similar to the extended $[\rm C\,II]$ emission found in some recent observations at high redshifts. △ Less

Submitted 1 July, 2024; v1 submitted 16 February, 2024; originally announced February 2024.

Comments: 21 pages, 11 figures. Accepted for publication in A&A, version matches accepted version

On the evolution of the observed Mass-to-Length relationship for star-forming filaments

Authors: Jiancheng Feng, Rowan J. Smith, Alvaro Hacar, Susan E. Clark, Daniel Seifried

Abstract: The interstellar medium is threaded by a hierarchy of filaments from large scales (~ 100 pc) to small scales (~ 0.1pc). The masses and lengths of these nested structures may reveal important constraints for cloud formation and evolution, but it is difficult to investigate from an evolutionary perspective using single observations. In this work, we extract simulated molecular clouds from the Cloud… ▽ More The interstellar medium is threaded by a hierarchy of filaments from large scales (~ 100 pc) to small scales (~ 0.1pc). The masses and lengths of these nested structures may reveal important constraints for cloud formation and evolution, but it is difficult to investigate from an evolutionary perspective using single observations. In this work, we extract simulated molecular clouds from the Cloud Factory galactic-scale ISM suite in combination with 3D Monte Carlo radiative transfer code POLARIS to investigate how filamentary structure evolves over time. We produce synthetic dust continuum observations in three regions with a series of snapshots and use the Filfinder algorithm to identify filaments in the dust derived column density maps. When the synthetic filaments mass and length are plotted on an M-L plot, we see a scaling relation of $L\propto M^{0.45}$ similar to that seen in observations, and find that the filaments are thermally supercritical. Projection effects systematically affect the masses and lengths measured for the filaments, and are particularly severe in crowded regions. In the filament Mass-Length (M-L) diagram we identify three main evolutionary mechanisms: accretion, segmentation, and dispersal. In particular we find that the filaments typically evolve from smaller to larger masses in the observational M-L plane, indicating the dominant role of accretion in filament evolution. Moreover, we find a potential correlation between line mass and filament growth rate. Once filaments are actively star forming they then segment into smaller sections, or are dispersed by internal or external forces. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: 19 pages, 18 figures

The Molecular Cloud Lifecycle II: Formation and Destruction of Molecular Clouds Diagnosed via H$_2$ Fluorescent Emission Emission

Authors: Blakesley Burkhart, Shmuel Bialy, Daniel Seifried, Stefanie Walch, Erika Hamden, Thomas J. Haworth, Keri Hoadley, Shuo Kong, Madisen Johnson, Sarah Jeffreson, Mark R. Krumholz, Min-Young Lee, Amiel Sternberg, Neal J. Turner

Abstract: Molecular hydrogen (H$_2$) formation and dissociation are key processes that drive the gas lifecycle in galaxies. Using the SImulating the LifeCycle of Molecular Clouds (SILCC) zoom-in simulation suite, we explore the utility of future observations of H$_2$ dissociation and formation for tracking the lifecycle of molecular clouds. The simulations used in this work include non-equilibrium H$_2$ for… ▽ More Molecular hydrogen (H$_2$) formation and dissociation are key processes that drive the gas lifecycle in galaxies. Using the SImulating the LifeCycle of Molecular Clouds (SILCC) zoom-in simulation suite, we explore the utility of future observations of H$_2$ dissociation and formation for tracking the lifecycle of molecular clouds. The simulations used in this work include non-equilibrium H$_2$ formation, stellar radiation, sink particles, and turbulence. We find that, at early times in the cloud evolution, H$_2$ formation rapidly outpaces dissociation and molecular clouds build their mass from the atomic reservoir in their environment. Rapid H$_2$ formation is also associated with a higher early star formation rate. For the clouds studied here, H$_2$ is strongly out of chemical equilibrium during the early stages of cloud formation but settles into a bursty chemical steady-state about 2 Myrs after the first stars form. At the latest stage of cloud evolution, dissociation outweighs formation and the clouds enter a dispersal phase. We discuss how theories for the molecular cloud lifecycle and the star formation efficiency may be distinguished with observational measurements of H$_2$ fluorescence with a space-based high-resolution FUV spectrometer, such as the proposed Hyperion and Eos NASA Explorer missions. Such missions would enable measurements of the H$_2$ dissociation and formation rates, which we demonstrate can be connected to different phases in a molecular cloud's star-forming life, including cloud building, rapidly star-forming, H$_2$ chemical equilibrium, and cloud destruction. △ Less

Submitted 5 February, 2024; v1 submitted 2 February, 2024; originally announced February 2024.

Comments: Submitted to ApJ, comments welcome

On the origin of V-shaped polarisation spectra in molecular clouds

Authors: Daniel Seifried, Stefanie Walch, Thorsten Balduin

Abstract: In this work we extend previous theoretical works to gain a better understanding of the origin of recently observed polarisation degree spectra of molecular clouds, which show a so-called V-shape, i.e. a pronounced minimum around 350 $μ$m. For this purpose, we present results of semi-analytical dust polarisation models. We benchmark our model against dust polarisation radiative transfer calculatio… ▽ More In this work we extend previous theoretical works to gain a better understanding of the origin of recently observed polarisation degree spectra of molecular clouds, which show a so-called V-shape, i.e. a pronounced minimum around 350 $μ$m. For this purpose, we present results of semi-analytical dust polarisation models. We benchmark our model against dust polarisation radiative transfer calculations performed with POLARIS. We show that V-shaped polarisation spectra can only be obtained if two dust phases, one dense and cold and one warm and dilute phase, are present along the line of sight. In contrast to previous results, no correlation between the alignment efficiency of silicate grains and the dust temperature is required; carbon grains are assumed to be not aligned with the magnetic field. We find that the V-shape is the stronger pronounced the larger the density and temperature contrast between both phases is. Moreover, the destruction of carbon grains by UV radiation in the warm and dilute phase leads to a significantly more pronounced V-shape in the polarisation spectrum. Reducing the alignment efficiency in the cold and dense phase also results in a more pronounced V-shape, its effect, however, is smaller than that of the UV-induced carbon grain destruction. Furthermore, we present a first, self-consistent polarisation spectrum obtained from a 3D, magneto-hydrodynamical molecular cloud simulation. The spectrum matches well with our semi-analytical prediction demonstrating the potential of such complex 3D simulations to study polarisation spectra. Comparing our model results with actual observations indicates that carbon grain destruction in illuminated regions might be required to match these observations. Reducing the alignment efficiency of silicate grains in the cold and dense phase would further improve the match between both data, however, it appears to not be a necessity. △ Less

Submitted 26 October, 2023; originally announced October 2023.

Comments: 15 pages including appendix, 11 figures, submitted to A&A, comments are welcome

Unravelling the structure of magnetised molecular clouds with SILCC-Zoom: sheets, filaments and fragmentation

Authors: S. Ganguly, S. Walch, D. Seifried, S. D. Clarke, M. Weis

Abstract: To what extent magnetic fields affect how molecular clouds (MCs) fragment and create dense structures is an open question. We present a numerical study of cloud fragmentation using the SILCC-Zoom simulations. These simulations follow the self-consistent formation of MCs in a few hundred parsec sized region of a stratified galactic disc; and include magnetic fields, self-gravity, supernova-driven t… ▽ More To what extent magnetic fields affect how molecular clouds (MCs) fragment and create dense structures is an open question. We present a numerical study of cloud fragmentation using the SILCC-Zoom simulations. These simulations follow the self-consistent formation of MCs in a few hundred parsec sized region of a stratified galactic disc; and include magnetic fields, self-gravity, supernova-driven turbulence, as well as a non-equilibrium chemical network. To discern the role of magnetic fields in the evolution of MCs, we study seven simulated clouds, five with magnetic fields, and two without, with a maximum resolution of 0.1 parsec. Using a dendrogram we identify hierarchical structures which form within the clouds. Overall, the magnetised clouds have more mass in a diffuse envelope with a number density between 1-100 cm$^{-3}$. We find that six out of seven clouds are sheet-like on the largest scales, as also found in recent observations, and with filamentary structures embedded within, consistent with the bubble-driven MC formation mechanism. Hydrodynamic simulations tend to produce more sheet-like structures also on smaller scales, while the presence of magnetic fields promotes filament formation. Analysing cloud energetics, we find that magnetic fields are dynamically important for less dense, mostly but not exclusively atomic structures (typically up to $\sim 100 - 1000$~cm$^{-3}$), while the denser, potentially star-forming structures are energetically dominated by self-gravity and turbulence. In addition, we compute the magnetic surface term and demonstrate that it is generally confining, and some atomic structures are even magnetically held together. In general, magnetic fields delay the cloud evolution and fragmentation by $\sim$ 1 Myr. △ Less

Submitted 17 July, 2023; originally announced July 2023.

Comments: 21 pages, 13 figures, accepted for publication in MNRAS

How do supernova remnants cool? -- I. Morphology, optical emission lines, and shocks

Authors: Ekaterina I. Makarenko, Stefanie Walch, Seamus D. Clarke, Daniel Seifried, Thorsten Naab, Pierre C. Nürnberger, Tim-Eric Rathjen

Abstract: Supernovae (SNe) inject $\sim 10^{51}$ erg in the interstellar medium, thereby shocking and heating the gas. A substantial fraction of this energy is later lost via radiative cooling. We present a post-processing module for the FLASH code to calculate the cooling radiation from shock-heated gas using collisional excitation data from MAPPINGS V. When applying this tool to a simulated SN remnant (SN… ▽ More Supernovae (SNe) inject $\sim 10^{51}$ erg in the interstellar medium, thereby shocking and heating the gas. A substantial fraction of this energy is later lost via radiative cooling. We present a post-processing module for the FLASH code to calculate the cooling radiation from shock-heated gas using collisional excitation data from MAPPINGS V. When applying this tool to a simulated SN remnant (SNR), we find that most energy is emitted in the EUV. However, optical emission lines ($[$O III$]$, $[$N II$]$, $[$S II$]$, H$α$, H$β$) are usually best observable. Our shock detection scheme shows that [S II] and [N II] emissions arise from the thin shell surrounding the SNR, while [O III], H$\rm α$, and H$\rm β$ originate from the volume-filling hot gas inside the SNR bubble. We find that the optical emission lines are affected by the SNR's complex structure and its projection onto the plane of the sky because the escaping line luminosity can be reduced by 10 -- 80\% due to absorption along the line-of-sight. Additionally, the subtraction of contaminating background radiation is required for the correct classification of an SNR on the oxygen or sulphur BPT diagrams. The electron temperature and density obtained from our synthetic observations match well with the simulation but are very sensitive to the assumed metallicity. △ Less

Submitted 12 May, 2023; originally announced May 2023.

Comments: 18 pages, 13 figures, accepted for publication in MNRAS

HyGAL: Characterizing the Galactic ISM with observations of hydrides and other small molecules II. The absorption line survey with the IRAM 30 m telescope

Authors: W. -J. Kim, P. Schilke, D. A. Neufeld, A. M. Jacob, Á. Sánchez-Monge, D. Seifried, B. Godard, K. M. Menten, S. Walch, E. Falgarone, V. S. Veena, S. Bialy, T. Möller, F. Wyrowski

Abstract: As a complement to the HyGAL Stratospheric Observatory for Infrared Astronomy Legacy Program, we report the results of a ground-based absorption line survey of simple molecules in diffuse and translucent Galactic clouds. Using the Institut de Radioastronomie Millimétrique (IRAM) 30 m telescope, we surveyed molecular lines in the 2 mm and 3 mm wavelength ranges toward 15 millimeter continuum source… ▽ More As a complement to the HyGAL Stratospheric Observatory for Infrared Astronomy Legacy Program, we report the results of a ground-based absorption line survey of simple molecules in diffuse and translucent Galactic clouds. Using the Institut de Radioastronomie Millimétrique (IRAM) 30 m telescope, we surveyed molecular lines in the 2 mm and 3 mm wavelength ranges toward 15 millimeter continuum sources. These sources, which are all massive star-forming regions located mainly in the first and second quadrants of the Milky Way, form the subset of the HyGAL sample that can be observed by the IRAM 30 m telescope. We detected HCO$^+$ absorption lines toward 14 sightlines, toward which we identified 78 foreground cloud components, as well as lines from HCN, HNC, C$_2$H, and c-C$_3$H$_2$ toward most sightlines. In addition, CS and H$_2$S absorption lines are found toward at least half of the continuum sources. Static Meudon photodissociation region (PDR) isobaric models that consider ultraviolet-dominated chemistry were unable to reproduce the column densities of all seven molecular species by just a factor of a few, except for H$_2$S. The inclusion of other formation routes driven by turbulent dissipation could possibly explain the observed high column densities of these species in diffuse clouds. There is a tentative trend for H$_2$S and CS abundances relative to H$_2$ to be larger in diffuse clouds ($X$(H$_2$S) and $X$(CS) $\sim 10^{-8} - 10^{-7}$) than in translucent clouds ($X$(H$_2$S) and $X$(CS) $\sim 10^{-9} - 10^{-8}$) toward a small sample; however, a larger sample is required in order to confirm this trend. The derived H$_2$S column densities are higher than the values predicted from the isobaric PDR models, suggesting that chemical desorption of H$_2$S from sulfur-containing ice mantles may play a role in increasing the H$_2$S abundance. △ Less

Submitted 19 December, 2022; originally announced December 2022.

Comments: 41 pages, 25 figures, 10 tables, and 5 appendices. Accepted for publication in A&A

Journal ref: A&A 670, A111 (2023)

Hyperion: The origin of the stars A far-UV space telescope for high-resolution spectroscopy over wide fields

Authors: Erika Hamden, David Schiminovich, Shouleh Nikzad, Neal J. Turner, Blakesley Burkhart, Thomas J. Haworth, Keri Hoadley, Jinyoung Serena Kim, Shmuel Bialyh, Geoff Bryden, Haeun Chung, Nia Imara, Rob Kennicutt, Jorge Pineda, Shuo Konga, Yasuhiro Hasegawa, Ilaria Pascucci, Benjamin Godard, Mark Krumholz, Min-Young Lee, Daniel Seifried, Amiel Sternberg, Stefanie Walch, Miles Smith, Stephen C. Unwin , et al. (8 additional authors not shown)

Abstract: We present Hyperion, a mission concept recently proposed to the December 2021 NASA Medium Explorer announcement of opportunity. Hyperion explores the formation and destruction of molecular clouds and planet-forming disks in nearby star-forming regions of the Milky Way. It does this using long-slit, high-resolution spectroscopy of emission from fluorescing molecular hydrogen, which is a powerful fa… ▽ More We present Hyperion, a mission concept recently proposed to the December 2021 NASA Medium Explorer announcement of opportunity. Hyperion explores the formation and destruction of molecular clouds and planet-forming disks in nearby star-forming regions of the Milky Way. It does this using long-slit, high-resolution spectroscopy of emission from fluorescing molecular hydrogen, which is a powerful far-ultraviolet (FUV) diagnostic. Molecular hydrogen (H2) is the most abundant molecule in the universe and a key ingredient for star and planet formation, but is typically not observed directly because its symmetric atomic structure and lack of a dipole moment mean there are no spectral lines at visible wavelengths and few in the infrared. Hyperion uses molecular hydrogen's wealth of FUV emission lines to achieve three science objectives: (1) determining how star formation is related to molecular hydrogen formation and destruction at the boundaries of molecular clouds; (2) determining how quickly and by what process massive star feedback disperses molecular clouds; and (3) determining the mechanism driving the evolution of planet-forming disks around young solar-analog stars. Hyperion conducts this science using a straightforward, highly-efficient, single-channel instrument design. Hyperion's instrument consists of a 48 cm primary mirror, with an f/5 focal ratio. The spectrometer has two modes, both covering 138.5-161.5 nm bandpasses. A low resolution mode has a spectral resolution of R>10,000 with a slit length of 65 arcmin, while the high resolution mode has a spectral resolution of R>50,000 over a slit length of 5 armin. Hyperion occupies a 2 week long, high-earth, Lunar resonance TESS-like orbit, and conducts 2 weeks of planned observations per orbit, with time for downlinks and calibrations. Hyperion was reviewed as Category I, which is the highest rating possible, but was not selected. △ Less

Submitted 13 December, 2022; originally announced December 2022.

Comments: Accepted to JATIS, 9 Figures

SILCC VII -- Gas kinematics and multiphase outflows of the simulated ISM at high gas surface densities

Authors: Tim-Eric Rathjen, Thorsten Naab, Stefanie Walch, Daniel Seifried, Philipp Girichidis, Richard Wünsch

Abstract: We present magnetohydrodynamic (MHD) simulations of the star-forming multiphase interstellar medium (ISM) in stratified galactic patches with gas surface densities $Σ_\mathrm{gas} =$ 10, 30, 50, and 100 $\mathrm{M_\odot\,pc^{-2}}$. The SILCC project simulation framework accounts for non-equilibrium thermal and chemical processes in the warm and cold ISM. The sink-based star formation and feedback… ▽ More We present magnetohydrodynamic (MHD) simulations of the star-forming multiphase interstellar medium (ISM) in stratified galactic patches with gas surface densities $Σ_\mathrm{gas} =$ 10, 30, 50, and 100 $\mathrm{M_\odot\,pc^{-2}}$. The SILCC project simulation framework accounts for non-equilibrium thermal and chemical processes in the warm and cold ISM. The sink-based star formation and feedback model includes stellar winds, hydrogen-ionising UV radiation, core-collapse supernovae, and cosmic ray (CR) injection and diffusion. The simulations follow the observed relation between $Σ_\mathrm{gas}$ and the star formation rate surface density $Σ_\mathrm{SFR}$. CRs qualitatively change the outflow phase structure. Without CRs, the outflows transition from a two-phase (warm and hot at 1 kpc) to a single-phase (hot at 2 kpc) structure. With CRs, the outflow always has three phases (cold, warm, and hot), dominated in mass by the warm phase. The impact of CRs on mass loading decreases for higher $Σ_\mathrm{gas}$ and the mass loading factors of the CR-supported outflows are of order unity independent of $Σ_\mathrm{SFR}$. Similar to observations, vertical velocity dispersions of the warm ionised medium (WIM) and the cold neutral medium (CNM) correlate with the star formation rate as $σ_\mathrm{z} \propto Σ_\mathrm{SFR}^a$, with $a \sim 0.20$. In the absence of stellar feedback, we find no correlation. The velocity dispersion of the WIM is a factor $\sim 2.2$ higher than that of the CNM, in agreement with local observations. For $Σ_\mathrm{SFR} \gtrsim 1.5 \times 10^{-2}\,\mathrm{M}_\odot\,\mathrm{yr}^{-1}\,\mathrm{kpc}^{-2}$ the WIM motions become supersonic. △ Less

Submitted 13 April, 2023; v1 submitted 28 November, 2022; originally announced November 2022.

Comments: 20 pages, 12 figures, accepted by MNRAS

The evolution of HCO$^{+}$ in molecular clouds using a novel chemical post-processing algorithm

Authors: M. Panessa, D. Seifried, S. Walch, B. Gaches, A. T. Barnes, F. Bigiel, L. Neumann

Abstract: Modeling the chemistry of molecular clouds is critical to accurately simulating their evolution. To reduce computational cost, 3D simulations generally restrict their chemistry to species with strong heating and cooling effects. Time-dependent information about the evolution of other species is therefore often neglected. We address this gap by post-processing tracer particles in the SILCC-Zoom mol… ▽ More Modeling the chemistry of molecular clouds is critical to accurately simulating their evolution. To reduce computational cost, 3D simulations generally restrict their chemistry to species with strong heating and cooling effects. Time-dependent information about the evolution of other species is therefore often neglected. We address this gap by post-processing tracer particles in the SILCC-Zoom molecular cloud simulations. Using a chemical network of 39 species and 301 reactions (including freeze-out of CO and H$_{2}$O), and a novel algorithm to reconstruct a density grid from sparse tracer particle data, we produce time-dependent density distributions for various species. We focus upon the evolution of HCO$^{+}$, which is a critical formation reactant of CO but is not typically modeled on-the-fly. We find that $\sim90$% of the HCO$^{+}$ content of the cold molecular gas forms in situ around $n_\textrm{HCO$^+$}\simeq10^3$-$10^4$ cm$^{-3}$, over a time-scale of approximately 1 Myr. The remaining $\sim10$% forms at high extinction sites, with minimal turbulent mixing out into the less dense gas. We further show that the dominant HCO$^{+}$ formation pathway is dependent on the visual extinction, with the reaction H$_{3}^{+}$ + CO contributing 90% of the total HCO$^{+}$ production above $A_\textrm{V,3D}=3$. We produce the very first maps of the HCO$^{+}$ column density, $N$(HCO$^{+}$), and show that it reaches values as high as $10^{15}$ cm$^{-2}$. We find that 50% of the HCO$^+$ mass is located within $A_\textrm{V}\sim10$-30, in a density range of $10^{3.5}$-$10^{4.5}$ cm$^{-3}$. Our maps of $N$(HCO$^{+}$) are shown to be in good agreement with recent observations of the W49A star-forming region. △ Less

Submitted 21 July, 2023; v1 submitted 12 October, 2022; originally announced October 2022.

Comments: 25 pages including appendix, 22 figures, accepted for publication in MNRAS

The Virial Balance of Molecular Clumps and Cores in Colliding Magnetized Flows

Authors: Michael Weis, Stefanie Walch, Daniel Seifried, Shashwata Ganguly

Abstract: We simulate the formation of molecular clouds in colliding flows of warm neutral medium with the adaptive mesh refinement code {\sc Flash}. We include a chemical network to treat heating and cooling and to follow the formation of molecular gas. For the forming molecular clumps and cores in four different simulations with varying initial magnetic field strength between 0.01 - 5$\,μ$G, we carry out… ▽ More We simulate the formation of molecular clouds in colliding flows of warm neutral medium with the adaptive mesh refinement code {\sc Flash}. We include a chemical network to treat heating and cooling and to follow the formation of molecular gas. For the forming molecular clumps and cores in four different simulations with varying initial magnetic field strength between 0.01 - 5$\,μ$G, we carry out a full virial analysis including all time-independent surface and volume terms as well as the time-dependent term. The initial magnetic field strength influences the fragmentation properties of the forming cloud because it prohibits motions perpendicular to the field direction and hence alters, or even suppresses, the formation of filamentary substructures. Molecular clump and core formation occurs anyhow. As a result, with increasing field strength, we find more fragments with a smaller average mass; yet the initial field strength is dynamically not relevant for the fragments which constitute our molecular clumps and cores. %yet the magnetic field overall is dynamically negligible for the fragments which constitute our molecular clumps and cores. The molecular clumps are mostly unbound, probably transient objects, which seem to be weakly confined by ram pressure or thermal pressure, indicating that they are swept up by the turbulent flow. They experience significant fluctuations in the mass flux through their surface, indicating that the Eulerian reference frame gives rise to a dominant time-dependent term due to their ill-defined nature. We define the cores to encompass molecular gas, which is additionally highly shielded. Most cores are in gravitational-kinetic equipartition and are already well described by the common virial parameter $α_\mathrm{vir}$ (as can be seen from the Heyer relation), while some undergo minor dispersion by kinetic surface effects. △ Less

Submitted 24 August, 2022; originally announced August 2022.

Comments: 19 pages, 35 figures, submitted to MNRAS

Understanding star formation in molecular clouds IV. Column density PDFs from quiescent to massive molecular clouds

Authors: N. Schneider, V. Ossenkopf-Okada, S. Clarke, R. S. Klessen, S. Kabanovic, T. Veltchev, S. Bontemps, S. Dib, T. Csengeri, C. Federrath, J. Di Francesco, F. Motte, Ph. Andre, D. Arzoumanian, J. R. Beattie, L. Bonne, P. Didelon, D. Elia, V. Koenyves, A. Kritsuk, B. Ladjelate, Ph. Myers, S. Pezzuto, J. F. Robitaille, A. Roy , et al. (4 additional authors not shown)

Abstract: We present N-PDFs of 29 Galactic regions obtained from Herschel imaging at high angular resolution, covering diffuse and quiescent clouds, and those showing low-, intermediate-, and high-mass star formation (SF), and characterize the cloud structure using the Delta-variance tool. The N-PDFs are double-log-normal at low column densities, and display one or two power law tails (PLTs) at higher colum… ▽ More We present N-PDFs of 29 Galactic regions obtained from Herschel imaging at high angular resolution, covering diffuse and quiescent clouds, and those showing low-, intermediate-, and high-mass star formation (SF), and characterize the cloud structure using the Delta-variance tool. The N-PDFs are double-log-normal at low column densities, and display one or two power law tails (PLTs) at higher column densities. For diffuse, quiescent, and low-mass SF clouds, we propose that the two log-normals arise from the atomic and molecular phase, respectively. For massive clouds, we suggest that the first log-normal is built up by turbulently mixed H2 and the second one by compressed (via stellar feedback) molecular gas. Nearly all clouds have two PLTs with slopes consistent with self-gravity, where the second one can be flatter or steeper than the first one. A flatter PLT could be caused by stellar feedback or other physical processes that slow down collapse and reduce the flow of mass toward higher densities. The steeper slope could arise if the magnetic field is oriented perpendicular to the LOS column density distribution. The first deviation point (DP), where the N-PDF turns from log-normal into a PLT, shows a clustering around values of a visual extinction of AV (DP1) around 2-5. The second DP, which defines the break between the two PLTs, varies strongly. Using the Delta-variance, we observe that the AV value, where the slope changes between the first and second PLT, increases with the characteristic size scale in the variance spectrum. We conclude that at low column densities, atomic and molecular gas is turbulently mixed, while at high column densities, the gas is fully molecular and dominated by self-gravity. The best fitting model N-PDFs of molecular clouds is thus one with log-normal low column density distributions, followed by one or two PLTs. △ Less

Submitted 29 July, 2022; originally announced July 2022.

Comments: Astronomy and Astrophysics in presss

Journal ref: A&A 666, A165 (2022)

Cosmic-ray-induced H$_2$ line emission: Astrochemical modeling and implications for JWST observations

Authors: Brandt A. L. Gaches, Shmuel Bialy, Thomas G. Bisbas, Marco Padovani, Daniel Seifried, Stefanie Walch

Abstract: Context: It has been proposed that H$_2$ near-infrared lines may be excited by cosmic rays and allow for a determination of the cosmic-ray ionization rate in dense gas. One-dimensional models show that measuring both the H$_2$ gas column density and H$_2$ line intensity enables a constraint on the cosmic-ray ionization rate as well as the spectral slope of low-energy cosmic-ray protons in the inte… ▽ More Context: It has been proposed that H$_2$ near-infrared lines may be excited by cosmic rays and allow for a determination of the cosmic-ray ionization rate in dense gas. One-dimensional models show that measuring both the H$_2$ gas column density and H$_2$ line intensity enables a constraint on the cosmic-ray ionization rate as well as the spectral slope of low-energy cosmic-ray protons in the interstellar medium (ISM). Aims: We aim to investigate the impact of certain assumptions regarding the H$_2$ chemical models and ISM density distributions on the emission of cosmic-ray induced H$_2$ emission lines. This is of particular importance for utilizing observations of these lines with the James Webb Space Telescope to constrain the cosmic-ray ionization rate. Methods: We compare the predicted emission from cosmic-ray induced, ro-vibrationally excited H$_2$ emission lines for different one- and three-dimensional models with varying assumptions on the gas chemistry and density distribution. Results: We find that the model predictions of the H$_2$ line intensities for the (1-0)S(0), (1-0)Q(2), (1-0)O(2) and (1-0)O(4) transitions at 2.22, 2.41, 2.63 and 3.00 $μ$m, respectively, are relatively independent of the astro-chemical model and the gas density distribution when compared against the H$_2$ column density, making them robust tracer of the cosmic-ray ionization rate. Conclusions: We recommend the use of ro-vibrational H$_2$ line emission in combination with estimation of the cloud's H$_2$ column density, to constrain the ionization rate and the spectrum of low energy cosmic-rays. △ Less

Submitted 6 July, 2022; v1 submitted 20 June, 2022; originally announced June 2022.

Comments: Accepted by A&A

Journal ref: A&A 664, A150 (2022)

CO and [CII] line emission of molecular clouds -- the impact of stellar feedback and non-equilibrium chemistry

Authors: S. Ebagezio, D. Seifried, S. Walch, P. C. Nürnberger, T. E. Rathjen, T. Naab

Abstract: We analyse synthetic $^{12}$CO, $^{13}$CO, and [CII] emission maps of simulated molecular clouds of the SILCC-Zoom project, which include an on-the-fly evolution of H$_2$, CO, and C$^+$. We use simulations of hydrodynamical and magnetohydrodynamical clouds, both with and without stellar feedback. We introduce a novel post-processing of the C$^+$ abundance using CLOUDY, to account for further ioniz… ▽ More We analyse synthetic $^{12}$CO, $^{13}$CO, and [CII] emission maps of simulated molecular clouds of the SILCC-Zoom project, which include an on-the-fly evolution of H$_2$, CO, and C$^+$. We use simulations of hydrodynamical and magnetohydrodynamical clouds, both with and without stellar feedback. We introduce a novel post-processing of the C$^+$ abundance using CLOUDY, to account for further ionization states of carbon due to stellar radiation. We report the first self-consistent synthetic emission maps of [CII] in feedback bubbles, largely devoid of emission inside them, as recently found in observations. The C$^+$ mass is only poorly affected by stellar feedback but the [CII] luminosity increases by $50 - 85$ per cent compared to runs without feedback. Furthermore, we investigate the capability of the CO/[CII] line ratio as a tracer of the amount of H$_2$ in the clouds and their evolutionary stage. We obtain, for both $^{12}$CO and $^{13}$CO, no clear trend of the luminosity ratio, $L_\mathrm{CO}/L_\mathrm{[CII]}$. It can therefore \textit{not} be used as a reliable measure of the H$_2$ mass fraction. We note a monotonic relation between $L_\mathrm{CO}/L_\mathrm{[CII]}$ and the H$_2$ fraction when considering the ratio for individual pixels of our synthetic maps, but with large scatter. Moreover, we show that assuming chemical equilibrium results in an overestimation of H$_2$ and CO masses of up to 110 and 30 per cent, respectively, and in an underestimation of H and C$^+$ masses of 65 and 7 per cent, respectively. In consequence, $L_\mathrm{CO}$ would be overestimated by up to 50 per cent, and $L_\mathrm{C[II]}$ be underestimated by up to 35 per cent. Hence, the assumption of chemical equilibrium in molecular cloud simulations introduces intrinsic errors of a factor of up to $\sim2$ in chemical abundances, luminosities and luminosity ratios. △ Less

Submitted 11 July, 2022; v1 submitted 13 June, 2022; originally announced June 2022.

Comments: 18 pages, 18 figures including appendix, modified. Submitted to MNRAS

Tree-based solvers for adaptive mesh refinement code FLASH -- III: a novel scheme for radiation pressure on dust and gas and radiative transfer from diffuse sources

Authors: A. Klepitko, S. Walch, R. Wünsch, D. Seifried, F. Dinnbier, S. Haid

Abstract: Radiation is an important contributor to the energetics of the interstellar medium, yet its transport is difficult to solve numerically. We present a novel approach towards solving radiative transfer of diffuse sources via backwards ray tracing. Here we focus on the radiative transfer of infrared radiation and the radiation pressure on dust. The new module, \textsc{TreeRay/RadPressure}, is an exte… ▽ More Radiation is an important contributor to the energetics of the interstellar medium, yet its transport is difficult to solve numerically. We present a novel approach towards solving radiative transfer of diffuse sources via backwards ray tracing. Here we focus on the radiative transfer of infrared radiation and the radiation pressure on dust. The new module, \textsc{TreeRay/RadPressure}, is an extension to the novel radiative transfer method \textsc{TreeRay} implemented in the grid-based MHD code {\sc Flash}. In \textsc{TreeRay/RadPressure}, every cell and every star particle is a source of infrared radiation. We also describe how gas, dust and radiation are coupled via a chemical network. This allows us to compute the local dust temperature in thermal equilibrium, leading to a significantly improvement over the classical grey approximation. In several tests, we demonstrate that the scheme produces the correct radiative intensities as well as the correct momentum input by radiation pressure. Subsequently, we apply our new scheme to model massive star formation from a collapsing, turbulent core of 150 ${\rm M}_\odot$. We include the effects of both, ionizing and infrared radiation on the dynamics of the core. We find that the newborn massive star prevents fragmentation in its proximity due to radiative heating. Over time, dust and radiation temperature equalize, while the gas temperature can be either warmer due to shock heating or colder due to insufficient dust-gas coupling. Compared to gravity, the effects of radiation pressure are insignificant for the stellar mass on the simulated time scale in this work. △ Less

Submitted 5 February, 2023; v1 submitted 19 April, 2022; originally announced April 2022.

Comments: 25 pages, 19 figures, accepted by MNRAS

SILCC-Zoom: the dynamic balance in molecular cloud substructures

Authors: S. Ganguly, S. Walch, S. D. Clarke, D. Seifried

Abstract: How molecular clouds fragment and create the dense structures which go on to form stars is an open question. We investigate the relative importance of different energy terms (kinetic, thermal, magnetic, and gravity - both self-gravity and tidal forces) for the formation and evolution of molecular clouds and their sub-structures based on the SILCC-Zoom simulations. These simulations follow the self… ▽ More How molecular clouds fragment and create the dense structures which go on to form stars is an open question. We investigate the relative importance of different energy terms (kinetic, thermal, magnetic, and gravity - both self-gravity and tidal forces) for the formation and evolution of molecular clouds and their sub-structures based on the SILCC-Zoom simulations. These simulations follow the self-consistent formation of cold molecular clouds down to scales of 0.1 pc from the diffuse supernova-driven interstellar medium in a stratified galactic disc. We study the time evolution of seven molecular clouds (five with magnetic fields and two without) for 1.5-2 Myr. Using a dendrogram, we identify hierarchical 3D sub-structures inside the clouds with the aim to understand their dynamics and distinguish between the theories of gravo-turbulent fragmentation and global hierarchical collapse. The virial analysis shows that the dense gas is indeed dominated by the interplay of gravity and turbulence, while magnetic fields and thermal pressure are only important for fluffy, atomic structures. Over time, gravitationally bound sub-structures emerge from a marginally bound medium (viral ratio $1 \leq α_{\rm vir}^{\rm vol} <2$) as a result of large-scale supernova-driven inflows rather than global collapse. A detailed tidal analysis shows that the tidal tensor is highly anisotropic. Yet the tidal forces are generally not strong enough to disrupt either large-scale or dense sub-structures but cause their deformation. By comparing tidal and crossing time scales, we find that tidal forces do not seem to be the main driver of turbulence within the molecular clouds. △ Less

Submitted 5 April, 2022; originally announced April 2022.

Comments: 20 pages, 18 figures, submitted to MNRAS

Initial Conditions for Star Formation: A Physical Description of the Filamentary ISM

Authors: A. Hacar, S Clark, F. Heitsch, J. Kainulainen, G. Panopoulou, D. Seifried, R. Smith

Abstract: The interstellar medium contains filamentary structure over a wide range of scales. Understanding the role of this structure, both as a conduit of gas across the scales and a diagnostic tool of local physics, is a major focus of star formation studies. We review recent progress in studying filamentary structure in the ISM, interpreting its properties in terms of physical processes, and exploring f… ▽ More The interstellar medium contains filamentary structure over a wide range of scales. Understanding the role of this structure, both as a conduit of gas across the scales and a diagnostic tool of local physics, is a major focus of star formation studies. We review recent progress in studying filamentary structure in the ISM, interpreting its properties in terms of physical processes, and exploring formation and evolution scenarios. We include structures from galactic-scale filaments to tenth-of-a-parsec scale filaments, comprising both molecular and atomic structures, from both observational and theoretical perspectives. In addition to the literature overview, we assemble a large amount of catalogue data from different surveys and provide the most comprehensive census of filamentary structures to date. Our census consists of 22 803 filamentary structures, facilitating a holistic perspective and new insights. We use our census to conduct a meta-analysis, leading to a description of filament properties over four orders of magnitudes in length and eight in mass. Our analysis emphasises the hierarchical and dynamical nature of filamentary structures. Filaments do not live in isolation, nor they generally resemble static structures close to equilibrium. We propose that accretion during filament formation and evolution sets some of the key scaling properties of filaments. This highlights the role of accretion during filament formation and evolution and also in setting the initial conditions for star formation. Overall, the study of filamentary structures during the past decade has been observationally driven. While great progress has been made on measuring the basic properties of filaments, our understanding of their formation and evolution is clearly lacking. In this context, we identify a number of directions and questions we consider most pressing for the field. △ Less

Submitted 17 March, 2022; originally announced March 2022.

Comments: 36 pages, 10 figures. To appear in Protostars and Planets VII, Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura

HyGAL: Characterizing the Galactic ISM with observations of hydrides and other small molecules -- I. Survey description and a first look toward W3(OH), W3 IRS5 and NGC 7538 IRS1

Authors: A. M. Jacob, D. A. Neufeld, P. Schilke, H. Wiesemeyer, W. Kim, S. Bialy, M. Busch, D. Elia, E. Falgarone, M. Gerin, B. Godard, R. Higgins, P. Hennebelle, N. Indriolo, D. C. Lis, K. M. Menten, A. Sanchez-Monge, V. Ossenkopf-Okada, M. R. Rugel, D. Seifried, P. Sonnentrucker, S. Walch, M. Wolfire, F. Wyrowski, V. Valdivia

Abstract: The HyGAL SOFIA legacy program surveys six hydride molecules -- ArH+, OH+, H2O+, SH, OH, and CH -- and two atomic constituents -- C+ and O -- within the diffuse interstellar medium (ISM) by means of absorption-line spectroscopy toward 25 bright Galactic background continuum sources. This detailed spectroscopic study is designed to exploit the unique value of specific hydrides as tracers and probes… ▽ More The HyGAL SOFIA legacy program surveys six hydride molecules -- ArH+, OH+, H2O+, SH, OH, and CH -- and two atomic constituents -- C+ and O -- within the diffuse interstellar medium (ISM) by means of absorption-line spectroscopy toward 25 bright Galactic background continuum sources. This detailed spectroscopic study is designed to exploit the unique value of specific hydrides as tracers and probes of different phases of the ISM, as demonstrated by recent studies with the Herschel Space Observatory. The observations performed under the HyGAL program will allow us to address several questions related to the lifecycle of molecular material in the ISM and the physical processes that impact its phase transition, such as: (1) What is the distribution function of the H2 fraction in the ISM? (2) How does the ionization rate due to low-energy cosmic-rays vary within the Galaxy? (3) What is the nature of interstellar turbulence, and what mechanisms lead to its dissipation? This overview discusses the observing strategy, synergies with ancillary and archival observations, the data reduction and analysis schemes adopted; and presents the first results obtained toward three of the survey targets, W3(OH), W3IRS5 and NGC7538IRS1. Robust measurements of the column densities of these hydrides -- obtained through widespread observations of absorption lines-- help address the questions raised, and there is a timely synergy between these observations and the development of theoretical models, particularly pertaining to the formation of H2 within the turbulent ISM. The provision of enhanced HyGAL data products will therefore serve as a legacy for future ISM studies. △ Less

Submitted 10 February, 2022; originally announced February 2022.

Comments: 33 pages, 17 figures, 7 tables. Accepted for publication in the Astrophysical Journal

Protostellar Outflows: a window to the past

Authors: P. F. Rohde, S. Walch, D. Seifried, A. P. Whitworth, S. D. Clarke

Abstract: During the early phases of low-mass star formation, episodic accretion causes the ejection of high-velocity outflow bullets, which carry a fossil record of the driving protostar's accretion history. We present 44 SPH simulations of $1\,\mathrm{M}_{\odot}$ cores, covering a wide range of initial conditions, and follow the cores for five free-fall times. Individual protostars are represented by sink… ▽ More During the early phases of low-mass star formation, episodic accretion causes the ejection of high-velocity outflow bullets, which carry a fossil record of the driving protostar's accretion history. We present 44 SPH simulations of $1\,\mathrm{M}_{\odot}$ cores, covering a wide range of initial conditions, and follow the cores for five free-fall times. Individual protostars are represented by sink particles, and the sink particles launch episodic outflows using a subgrid model. The Optics algorithm is used to identify individual episodic bullets within the outflows. The parameters of the overall outflow and the individual bullets are then used to estimate the age and energetics of the outflow, and the accretion events that triggered it; and to evaluate how reliable these estimates are, if observational uncertainties and selection effects (like inclination) are neglected. Of the commonly used methods for estimating outflow ages, it appears that those based on the length and speed of advance of the lobe are the most reliable in the early phases of evolution, and those based on the width of the outflow cavity and the speed of advance are most reliable during the later phases. We describe a new method that is almost as accurate as these methods, and reliable throughout the evolution. In addition we show how the accretion history of the protostar can be accurately reconstructed from the dynamics of the bullets if each lobe contains at least two bullets. The outflows entrain about ten times more mass than originally ejected by the protostar. △ Less

Submitted 9 January, 2022; originally announced January 2022.

Comments: 20 pages, 19 figures, published in MNRAS

Self-absorption in [CII], $^{12}$CO, and HI in RCW120. Building up a geometrical and physical model of the region

Authors: S. Kabanovic, N. Schneider, V. Ossenkopf-Okada, F. Falasca, R. Güsten, J. Stutzki, R. Simon, C. Buchbender, L. Anderson, L. Bonne, C. Guevara, R. Higgins, B. Koribalski, M. Luisi, M. Mertens, Y. Okada, M. Röllig, D. Seifried, M. Tiwari, F. Wyrowski, A. Zavagno, A. G. G. M. Tielens

Abstract: Revealing the 3D dynamics of HII regions and their associated molecular clouds is important for understanding the longstanding problem as to how stellar feedback affects the density structure and kinematics of the interstellar medium. We employed observations of the HII region RCW 120 in [CII], observed within the SOFIA legacy program FEEDBACK, and the $^{12}$CO and $^{13}$CO (3$\to$2) lines, obta… ▽ More Revealing the 3D dynamics of HII regions and their associated molecular clouds is important for understanding the longstanding problem as to how stellar feedback affects the density structure and kinematics of the interstellar medium. We employed observations of the HII region RCW 120 in [CII], observed within the SOFIA legacy program FEEDBACK, and the $^{12}$CO and $^{13}$CO (3$\to$2) lines, obtained with APEX. In addition we used HI data from the Southern Galactic Plane Survey. Two radiative transfer models were used to fit the observed data. A line profile analysis with the 1D non-LTE radiative transfer code SimLine proves that the CO emission cannot stem from a spherically symmetric molecular cloud configuration. With a two-layer multicomponent model, we then quantified the amount of warm background and cold foreground gas. There is a deficit of CO emission along the line-of-sight toward the center of the HII region which indicates that the HII region is associated with a flattened molecular cloud. Self-absorption in the CO line may hide signatures of infalling and expanding molecular gas. The [CII] emission arises from an expanding [CII] bubble and from the PDRs. A significant part of [CII] emission is absorbed in a cool (~60-100 K), low-density (<500 cm$^{-3}$) atomic foreground layer with a thickness of a few parsec. We propose that the RCW 120 HII region formed in a flattened molecular cloud and is now bursting out of its parental cloud. The compressed surrounding molecular layer formed a torus around the spherically expanding HII bubble. This scenario can possibly be generalized for other HII bubbles and would explain the observed "flat" structure of molecular clouds associated with HII bubbles. We suggest that the [CII] absorption observed in many star-forming regions is at least partly caused by low-density, cool, HI-envelopes surrounding the molecular clouds. △ Less

Submitted 21 April, 2022; v1 submitted 21 December, 2021; originally announced December 2021.

Journal ref: A&A 659, A36 (2022)

The width of Herschel filaments varies with distance

Authors: G. V. Panopoulou, S. E. Clark, A. Hacar, F. Heitsch, J. Kainulainen, E. Ntormousi, D. Seifried, R. J. Smith

Abstract: Context: Filamentary structures in nearby molecular clouds have been found to exhibit a characteristic width of 0.1 pc, as observed in dust emission. Understanding the origin of this universal width has become a topic of central importance in the study of molecular cloud structure and the early stages of star formation. Aims: We investigate how the recovered widths of filaments depend on the dista… ▽ More Context: Filamentary structures in nearby molecular clouds have been found to exhibit a characteristic width of 0.1 pc, as observed in dust emission. Understanding the origin of this universal width has become a topic of central importance in the study of molecular cloud structure and the early stages of star formation. Aims: We investigate how the recovered widths of filaments depend on the distance from the observer, by using previously published results from the Herschel Gould Belt Survey. Methods: We obtained updated estimates on the distances to nearby molecular clouds observed with Herschel by using recent results based on 3D dust extinction mapping and Gaia. We examined the widths of filaments from individual clouds separately, as opposed to treating them as a single population. We used these per-cloud filament widths to search for signs of variation amongst the clouds of the previously published study. Results: We find a significant dependence of the mean per-cloud filament width with distance. The distribution of mean filament widths for nearby clouds is incompatible with that of farther away clouds. The mean per-cloud widths scale with distance approximately as 4--5 times the beam size. We examine the effects of resolution by performing a convergence study of a filament profile in the Herschel image of the Taurus Molecular Cloud. We find that resolution can severely affect the shapes of radial profiles over the observed range of distances. Conclusions: We conclude that the data are inconsistent with 0.1 pc being the universal characteristic width of filaments. △ Less

Submitted 8 August, 2022; v1 submitted 15 November, 2021; originally announced November 2021.

Comments: An erratum to this manuscript has been published here: https://www.aanda.org/articles/aa/full_html/2022/07/aa42281e-21/aa42281e-21.html

Journal ref: A&A 657, L13 (2022)

On the accuracy of HI observations in molecular clouds -- More cold HI than thought?

Authors: D. Seifried, H. Beuther, S. Walch, J. Syed, J. D. Soler, P. Girichidis, R. Wünsch

Abstract: We present a study of the cold atomic hydrogen (HI) content of molecular clouds simulated within the SILCC-Zoom project. We produce synthetic observations of HI at 21 cm including HI self-absorption (HISA) and observational effects. We find that HI column densities, $N_\textrm{HI}$, of $\gtrsim$10$^{22}$ cm$^{-2}$ are frequently reached in molecular clouds with HI temperatures as low as $\sim$10 K… ▽ More We present a study of the cold atomic hydrogen (HI) content of molecular clouds simulated within the SILCC-Zoom project. We produce synthetic observations of HI at 21 cm including HI self-absorption (HISA) and observational effects. We find that HI column densities, $N_\textrm{HI}$, of $\gtrsim$10$^{22}$ cm$^{-2}$ are frequently reached in molecular clouds with HI temperatures as low as $\sim$10 K. Hence, HISA observations assuming a fixed HI temperature tend to underestimate the amount of cold HI in molecular clouds by a factor of 3 - 10 and produce an artificial upper limit of $N_\textrm{HI}$ around 10$^{21}$ cm$^{-2}$. We thus argue that the cold HI mass in clouds could be a factor of a few higher than previously estimated. Also $N_\textrm{HI}$-PDFs obtained from HISA observations might be subject to observational biases and should be considered with caution. The underestimation of cold HI in HISA observations is due to both the large HI temperature variations and the effect of noise in regions of high optical depth. We find optical depths of cold HI around 1 - 10 making optical depth corrections essential. We show that the high HI column densities ($\gtrsim$10$^{22}$ cm$^{-2}$) can in parts be attributed to the occurrence of up to 10 individual HI-H$_2$ transitions along the line of sight. This is also reflected in the spectra, necessitating Gaussian decomposition algorithms for their analysis. However, also for a single HI-H$_2$ transition, $N_\textrm{HI}$ frequently exceeds 10$^{21}$ cm$^{-2}$, challenging 1D, semi-analytical models. This is due to non-equilibrium chemistry effects and the fact that HI-H$_2$ transition regions usually do not possess a 1-dimensional geometry. Finally, we show that the HI gas is moderately supersonic with Mach numbers of a few. The corresponding non-thermal velocity dispersion can be determined via HISA observations within a factor of $\sim$2. △ Less

Submitted 3 March, 2022; v1 submitted 22 September, 2021; originally announced September 2021.

Comments: 20 pages including appendix, 24 figures, accepted for publication in MNRAS

Tree-based solvers for adaptive mesh refinement code FLASH -- II: radiation transport module TreeRay

Authors: Richard Wünsch, Stefanie Walch, František Dinnbier, Daniel Seifried, Sebastian Haid, Andre Klepitko, Anthony P. Whitworth, Jan Palouš

Abstract: The treatment of radiative transfer with multiple radiation sources is a critical challenge in simulations of star formation and the interstellar medium. In this paper we present the novel TreeRay method for solving general radiative transfer problems, based on reverse ray tracing combined with tree-based accelerated integration. We implement TreeRay in the adaptive mesh refinement code FLASH, as… ▽ More The treatment of radiative transfer with multiple radiation sources is a critical challenge in simulations of star formation and the interstellar medium. In this paper we present the novel TreeRay method for solving general radiative transfer problems, based on reverse ray tracing combined with tree-based accelerated integration. We implement TreeRay in the adaptive mesh refinement code FLASH, as a module of the tree solver developed by Wünsch et al. However, the method itself is independent of the host code and can be implemented in any grid based or particle based hydrodynamics code. A key advantage of TreeRay is that its computational cost is independent of the number of sources, making it suitable for simulations with many point sources (e.g. massive star clusters) as well as simulations where diffuse emission is important. A very efficient communication and tree-walk strategy enables TreeRay to achieve almost ideal parallel scalings. TreeRay can easily be extended with sub-modules to treat radiative transfer at different wavelengths and to implement related physical processes. Here, we focus on ionising (EUV) radiation and use the On-the-Spot approximation to test the method and its parameters. The ability to set the tree solver time step independently enables the speedy calculation of radiative transfer in a multi-phase interstellar medium, where the hydrodynamic time step is typically limited by the sound speed of the hot gas produced in stellar wind bubbles or supernova remnants. We show that complicated simulations of star clusters with feedback from multiple massive stars become feasible with TreeRay. △ Less

Submitted 20 May, 2021; originally announced May 2021.

Comments: 26 pages, 21 figures, accepted by MNRAS

SILCC VI -- Multi-phase ISM structure, stellar clustering, and outflows with supernovae, stellar winds, ionising radiation and cosmic rays

Authors: Tim-Eric Rathjen, Thorsten Naab, Philipp Girichidis, Stefanie Walch, Richard Wünsch, František Dinnbier, Daniel Seifried, Ralf S. Klessen, Simon C. O. Glover

Abstract: We present simulations of the multi-phase interstellar medium (ISM) at solar neighbourhood conditions including thermal and non-thermal ISM processes, star cluster formation, and feedback from massive stars: stellar winds, hydrogen ionising radiation computed with the novel TreeRay radiative transfer method, supernovae (SN), and the injection of cosmic rays (CR). N-body dynamics is computed with a… ▽ More We present simulations of the multi-phase interstellar medium (ISM) at solar neighbourhood conditions including thermal and non-thermal ISM processes, star cluster formation, and feedback from massive stars: stellar winds, hydrogen ionising radiation computed with the novel TreeRay radiative transfer method, supernovae (SN), and the injection of cosmic rays (CR). N-body dynamics is computed with a 4th-order Hermite integrator. We systematically investigate the impact of stellar feedback on the self-gravitating ISM with magnetic fields, CR advection and diffusion and non-equilibrium chemical evolution. SN-only feedback results in strongly clustered star formation with very high star cluster masses, a bi-modal distribution of the ambient SN densities, and low volume-filling factors (VFF) of warm gas, typically inconsistent with local conditions. Early radiative feedback prevents an initial starburst, reduces star cluster masses and outflow rates. Furthermore, star formation rate surface densities of $Σ_{\dot{M}_\star} = 1.4-5.9 \times 10^{-3}$ $\mathrm{M}_\odot\,\mathrm{yr}^{-1}\,\mathrm{kpc}^{-2}$, VFF$_\mathrm{warm} = 60-80$ per cent as well as thermal, kinetic, magnetic, and cosmic ray energy densities of the model including all feedback mechanisms agree well with observational constraints. On the short, 100 Myr, timescales investigated here, CRs only have a moderate impact on star formation and the multi-phase gas structure and result in cooler outflows, if present. Our models indicate that at low gas surface densities SN-only feedback only captures some characteristics of the star-forming ISM and outflows/inflows relevant for regulating star formation. Instead, star formation is regulated on star cluster scales by radiation and winds from massive stars in clusters, whose peak masses agree with solar neighbourhood estimates. △ Less

Submitted 29 March, 2021; v1 submitted 25 March, 2021; originally announced March 2021.

Comments: 23 pages, 16 figures, accepted by MNRAS, corrected an incorrect reference

Synthetic CO emission and the $X_{\rm CO}$ factor of young molecular clouds: a convergence study

Authors: Elisabeth M. A. Borchert, Stefanie Walch, Daniel Seifried, Seamus D. Clarke, Annika Franeck, Pierre Nürnberger

Abstract: The properties of synthetic CO emission from 3D simulations of forming molecular clouds are studied within the SILCC-Zoom project. Since the time scales of cloud evolution and molecule formation are comparable, the simulations include a live chemical network. Two sets of simulations with an increasing spatial resolution (d$x=3.9$ pc to d$x=0.06$ pc) are used to investigate the convergence of the s… ▽ More The properties of synthetic CO emission from 3D simulations of forming molecular clouds are studied within the SILCC-Zoom project. Since the time scales of cloud evolution and molecule formation are comparable, the simulations include a live chemical network. Two sets of simulations with an increasing spatial resolution (d$x=3.9$ pc to d$x=0.06$ pc) are used to investigate the convergence of the synthetic CO emission, which is computed by post-processing the simulation data with the RADMC-3D radiative transfer code. To determine the excitation conditions, it is necessary to include atomic hydrogen and helium alongside H$_2$, which increases the resulting CO emission by ~7-26 per cent. Combining the brightness temperature of $^{12}$CO and $^{13}$CO, we compare different methods to estimate the excitation temperature, the optical depth of the CO line and hence, the CO column density. An intensity-weighted average excitation temperature results in the most accurate estimate of the total CO mass. When the pixel-based excitation temperature is used to calculate the CO mass, it is over-/underestimated at low/high CO column densities where the assumption that $^{12}$CO is optically thick while $^{13}$CO is optically thin is not valid. Further, in order to obtain a converged total CO luminosity and hence <$X_{\rm CO}$> factor, the 3D simulation must have d$x\lesssim0.1$ pc. The <$X_{\rm CO}$> evolves over time and differs for the two clouds; yet pronounced differences with numerical resolution are found. Since high column density regions with a visual extinction larger than 3~mag are not resolved for d$x\gtrsim 1$~pc, in this case the H$_2$ mass and CO luminosity both differ significantly from the higher resolution results and the local $X_{\rm CO}$ is subject to strong noise. Our calculations suggest that synthetic CO emission maps are only converged for simulations with d$x\lesssim 0.1$ pc. △ Less

Submitted 19 November, 2021; v1 submitted 1 February, 2021; originally announced February 2021.

Comments: 23 pages, 22 figures, accepted for publication in MNRAS

Dissecting the super-critical filaments embedded in the 0.5 pc subsonic region of Barnard 5

Authors: Anika Schmiedeke, Jaime E. Pineda, Paola Caselli, Héctor G. Arce, Gary A Fuller, Alyssa A. Goodman, María José Maureira, Stella S. R. Offner, Dominique Segura-Cox, Daniel Seifried

Abstract: We characterize in detail the two ~0.3 pc long filamentary structures found within the subsonic region of Barnard 5. We use combined GBT and VLA observations of the molecular lines NH$_3$(1,1) and (2,2) at a resolution of 1800 au, as well as JCMT continuum observations at 850 and 450 $μ$m at a resolution of 4400 au and 3000 au, respectively. We find that both filaments are highly super-critical wi… ▽ More We characterize in detail the two ~0.3 pc long filamentary structures found within the subsonic region of Barnard 5. We use combined GBT and VLA observations of the molecular lines NH$_3$(1,1) and (2,2) at a resolution of 1800 au, as well as JCMT continuum observations at 850 and 450 $μ$m at a resolution of 4400 au and 3000 au, respectively. We find that both filaments are highly super-critical with a mean mass per unit length, $M/L$, of ~80 M$_\odot$ pc$^{-1}$, after background subtraction, with local increases reaching values of ~150 M$_\odot$ pc$^{-1}$. This would require a magnetic field strength of ~500 $μ$G to be stable against radial collapse. We extract equidistant cuts perpendicular to the spine of the filament and fit a modified Plummer profile as well as a Gaussian to each of the cuts. The filament widths (deconvolved FWHM) range between 6500-7000 au (~0.03 pc) along the filaments. This equals ~2.0 times the radius of the flat inner region. We find an anti-correlation between the central density and this flattening radius, suggestive of contraction. Further, we also find a strong correlation between the power-law exponent at large radii and the flattening radius. We note that the measurements of these three parameters fall in a plane and derive their empirical relation. Our high-resolution observations provide direct constraints of the distribution of the dense gas within super-critical filaments showing pre- and protostellar activity. △ Less

Submitted 1 January, 2021; originally announced January 2021.

Comments: Accepted for publication in ApJ. 23 pages, 14 figures, 3 tables

Magnetized filamentary gas flows feeding the young embedded cluster in Serpens South

Authors: Thushara G. S. Pillai, Dan P. Clemens, Stefan Reissl, Philip C. Myers, Jens Kauffmann, Enrique Lopez-Rodriguez, Felipe O. Alves, Gabriel P. Franco, Jonathan D. Henshaw, Karl M. Menten, Fumitaka Nakamura, Daniel Seifried, Koji Sugitani, Helmut Wiesemeyer

Abstract: Observations indicate that molecular clouds are strongly magnetized, and that magnetic fields influence the formation of stars. A key observation supporting the conclusion that molecular clouds are significantly magnetized is that the orientation of their internal structure is closely related to that of the magnetic field. At low column densities the structure aligns parallel with the field, where… ▽ More Observations indicate that molecular clouds are strongly magnetized, and that magnetic fields influence the formation of stars. A key observation supporting the conclusion that molecular clouds are significantly magnetized is that the orientation of their internal structure is closely related to that of the magnetic field. At low column densities the structure aligns parallel with the field, whereas at higher column densities, the gas structure is typically oriented perpendicular to magnetic fields, with a transition at visual extinctions $A_V\gtrsim{}3~\rm{}mag$. Here we use far-infrared polarimetric observations from the HAWC+ polarimeter on SOFIA to report the discovery of a further transition in relative orientation, i.e., a return to parallel alignment at $A_V\gtrsim{}21~\rm{}mag$ in parts of the Serpens South cloud. This transition appears to be caused by gas flow and indicates that magnetic supercriticality sets in near $A_V\gtrsim{}21~\rm{}mag$, allowing gravitational collapse and star cluster formation to occur even in the presence of relatively strong magnetic fields. △ Less

Submitted 29 September, 2020; originally announced September 2020.

Comments: 16 pages, 4 figures, Published in Nature Astronomy (August 2020). This is the authors' version before final edits. Link to the NA publication: https://www.nature.com/articles/s41550-020-1172-6

Journal ref: Nature Astronomy, August 2020

The impact of episodic outflow feedback on stellar multiplicity and the star formation efficiency

Authors: P. F. Rohde, S. Walch, S. D. Clarke, D. Seifried, A. P. Whitworth, A. Klepitko

Abstract: The accretion of material onto young protostars is accompanied by the launching of outflows. Observations show that accretion, and therefore also outflows, are episodic. However, the effects of episodic outflow feedback on the core-scale are not well understood. We have performed 88 Smoothed Particle Hydrodynamic simulations of turbulent dense $1 \, \mathrm{M}_{\odot}$ cores, to study the influenc… ▽ More The accretion of material onto young protostars is accompanied by the launching of outflows. Observations show that accretion, and therefore also outflows, are episodic. However, the effects of episodic outflow feedback on the core-scale are not well understood. We have performed 88 Smoothed Particle Hydrodynamic simulations of turbulent dense $1 \, \mathrm{M}_{\odot}$ cores, to study the influence of episodic outflow feedback on the stellar multiplicity and the star formation efficiency (SFE). Protostars are represented by sink particles, which use a sub-grid model to capture stellar evolution, inner-disc evolution, episodic accretion and the launching of outflows. By comparing simulations with and without episodic outflow feedback, we show that simulations with outflow feedback reproduce the binary statistics of young stellar populations, including the relative proportions of singles, binaries, triples, etc. and the high incidence of twin binaries with $q\geq 0.95$; simulations without outflow feedback do not. Entrainment factors (the ratio between total outflowing mass and initially ejected mass) are typically $\sim 7\pm 2$, but can be much higher if the total mass of stars formed in a core is low and/or outflow episodes are infrequent. By decreasing both the mean mass of the stars formed and the number of stars formed, outflow feedback reduces the SFE by about a factor of 2 (as compared with simulations that do not include outflow feedback). △ Less

Submitted 24 September, 2020; originally announced September 2020.

Comments: This is a pre-copyedited, author-produced PDF of an article accepted for publication in MNRAS following peer review. 19 pages, 12 figures

Magnetic fields in star-forming systems (II): examining dust polarization, the Zeeman effect, and the Faraday rotation measure as magnetic field tracers

Authors: Stefan Reissl, Amelia M. Stutz, Ralf S. Klessen, Daniel Seifried, Stefanie Walch

Abstract: The degree to which the formation and evolution of clouds and filaments in the interstellar medium is regulated by magnetic fields remains an open question. Yet the fundamental properties of the fields (strength and 3D morphology) are not readily observable. We investigate the potential for recovering magnetic field information from dust polarization, the Zeeman effect, and the Faraday rotation me… ▽ More The degree to which the formation and evolution of clouds and filaments in the interstellar medium is regulated by magnetic fields remains an open question. Yet the fundamental properties of the fields (strength and 3D morphology) are not readily observable. We investigate the potential for recovering magnetic field information from dust polarization, the Zeeman effect, and the Faraday rotation measure ($RM$) in a SILCC-Zoom magnetohydrodynamic (MHD) filament simulation. The object is analyzed at the onset of star formation, and it is characterized by a line-mass of about M/L $\sim 63\ M_{\odot}\ pc^{-1}$ out to a radius of $1\,$pc and a kinked 3D magnetic field morphology. We generate synthetic observations via POLARIS radiative transfer (RT) post-processing, and compare with an analytical model of helical or kinked field morphology to help interpreting the inferred observational signatures. We show that the tracer signals originate close to the filament spine. We find regions along the filament where the angular-dependency with the line-of-sight (LOS) is the dominant factor and dust polarization may trace the underlying kinked magnetic field morphology. We also find that reversals in the recovered magnetic field direction are not unambiguously associated to any particular morphology. Other physical parameters, such as density or temperature, are relevant and sometimes dominant compared to the magnetic field structure in modulating the observed signal. We demonstrate that the Zeeman effect and the $RM$ recover the line-of-sight magnetic field strength to within a factor 2.1 - 3.4. We conclude that the magnetic field morphology may not be unambiguously determined in low-mass systems by observations of dust polarization, Zeeman effect, or $RM$, whereas the field strengths can be reliably recovered. △ Less

Submitted 9 September, 2020; originally announced September 2020.

Comments: 22 pages, 17 figures, 3 tables

From parallel to perpendicular -- On the orientation of magnetic fields in molecular clouds

Authors: D. Seifried, S. Walch, M. Weis, S. Reissl, J. D. Soler, R. S. Klessen, P. R. Joshi

Abstract: We present synthetic dust polarization maps of simulated molecular clouds (MCs) with the goal to systematically explore the origin of the relative orientation of the magnetic field ($\bf{B}$) with respect to the MC sub-structures identified in density ($n$; 3D) and column density ($N$; 2D). The polarization maps are generated with the radiative transfer code POLARIS, including self-consistently ca… ▽ More We present synthetic dust polarization maps of simulated molecular clouds (MCs) with the goal to systematically explore the origin of the relative orientation of the magnetic field ($\bf{B}$) with respect to the MC sub-structures identified in density ($n$; 3D) and column density ($N$; 2D). The polarization maps are generated with the radiative transfer code POLARIS, including self-consistently calculated efficiencies for radiative torque alignment. The MCs are formed in two sets of 3D MHD simulations: in (i) colliding flows (CF), and (ii) the SILCC-Zoom simulations. In 3D, for the CF simulations with an initial field strength below $\sim$5 $μ$G, $\bf{B}$ is oriented parallel or randomly with respect to the $n$-structures. For CF runs with stronger initial fields and 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_\text{trans}$ $\simeq$ 10$^2$ - 10$^3$ cm$^{-3}$. We suggest that this flip happens if the MC'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. In 2D, we use the Projected Rayleigh Statistics (PRS) to study the orientation of $\bf{B}$. If present, the flip in orientation occurs at $N_\text{trans}$ $\simeq$ 10$^{21 - 21.5}$ cm$^{-2}$, similar to the observed transition value from sub- to supercritical magnetic fields in the ISM. However, projection effects can reduce the power of the PRS method: Depending on the MC or LOS, the projected maps of the SILCC-Zoom simulations do not always show the flip, although expected from the 3D morphology. Such projection effects can explain the variety of recently observed field configurations, in particular within a single MC. Finally, we do not find a correlation between the observed orientation of $\bf{B}$ and the $N$-PDF. △ Less

Submitted 25 August, 2020; v1 submitted 28 February, 2020; originally announced March 2020.

Comments: 20 pages, 12 figures, accepted for publication in MNRAS

Physical Properties of the star-forming clusters in NGC 6334

Authors: Mahya Sadaghiani, Álvaro Sánchez-Monge, Peter Schilke, Hauyu Baobab Liu, Seamus Clarke, Qizhou Zhang, Josep Miquel Girart, Daniel Seifried, Atefeh Aghababaei, Hua-bai Li, Carmen Juárez, Kwok Sun Tang

Abstract: We have studied the high-mass star-forming complex NGC 6334 with ALMA in the continuum emission at a frequency of 87.6 GHz, achieving a spatial resolution of 1300 au. Detecting 142 compact sources distributed over the whole observed area, we then used machine learning algorithms to group the compact cores in different clusters. A total of four main clusters were identified: NGC 6334-E, NGC 6334-I,… ▽ More We have studied the high-mass star-forming complex NGC 6334 with ALMA in the continuum emission at a frequency of 87.6 GHz, achieving a spatial resolution of 1300 au. Detecting 142 compact sources distributed over the whole observed area, we then used machine learning algorithms to group the compact cores in different clusters. A total of four main clusters were identified: NGC 6334-E, NGC 6334-I, NGC 6334-I(N) and NGC 6334-I(NW). The typical separations between cluster members (4000-12000 au) together with the core masses (0.2-100 $M_{\odot}$) are in agreement with turbulent fragmentation at scales of 0.1 pc. We find that the CMFs (core mass functions) show an excess of massive cores compared to the IMF. Typical uncertainties in temperature or unresolved multiplicity may mimic the observed excess of massive cores. Evidence of mass segregation in NGC 6334-I and NGC 6334-I(N), with the most massive cores located closer to the center, was found. Correlations between the physical properties of the clusters and their evolutionary stage were searched for, and found to show a larger separation between cores in the more evolved clusters, favouring the role of gas expulsion and stellar ejection over time. From these results, we suggest that NGC 6334-I(N) was primordially segregated, while NGC 6334-I may have become mass-segregated over time due to dynamical effects. Finally, the lack of massive cores in the most evolved cluster suggests that the gas reservoir may be already exhausted, while the less evolved clusters still have access to a larger gas mass reservoir. In general, the fragmentation in NGC 6334 at large scales (about 1 pc) is governed by turbulent pressure, while at smaller scales (a few hundred au), thermal pressure regulates the fragmentation process. △ Less

Submitted 25 November, 2019; v1 submitted 15 November, 2019; originally announced November 2019.

Comments: Accepted for publication in Astronomy & Astrophysics. 20 pages with 15 figures and 3 tables, plus Appendices with extra figures and tables. (Abstract modified)

Journal ref: A&A 635, A2 (2020)

On the extraction of power-law parts of the probability density functions in star-forming clouds

Authors: Todor V. Veltchev, Philipp Girichidis, Sava Donkov, Nicola Schneider, Orlin Stanchev, Lyubov Marinkova, Daniel Seifried, Ralf S. Klessen

Abstract: We present a new approach to extract the power-law part of a density/column-density probability density function (rho-pdf/N-pdf) in star-forming clouds. It is based on the mathematical method bPLFIT of Virkar & Clauset (2014) and assesses the power-law part of an arbitrary distribution, without any assumptions about the other part of this distribution. The slope and deviation point are derived as… ▽ More We present a new approach to extract the power-law part of a density/column-density probability density function (rho-pdf/N-pdf) in star-forming clouds. It is based on the mathematical method bPLFIT of Virkar & Clauset (2014) and assesses the power-law part of an arbitrary distribution, without any assumptions about the other part of this distribution. The slope and deviation point are derived as averaged values as the number of bins is varied. Neither parameter is sensitive to spikes and other local features of the tail. This adapted bPLFIT method is applied to two different sets of data from numerical simulations of star-forming clouds at scales 0.5 and 500 pc and displays rho-pdf and N-pdf evolution in agreement with a number of numerical and theoretical studies. Applied to Herschel data on the regions Aquila and Rosette, the method extracts pronounced power-law tails, consistent with those seen in simulations of evolved clouds. △ Less

Submitted 1 August, 2019; originally announced August 2019.

Comments: Accepted for publication in MNRAS; 15 pages, 13 figures

X-raying molecular clouds with a short flare: probing statistics of gas density and velocity fields

Authors: I. Khabibullin, E. Churazov, R. Sunyaev, C. Federrath, D. Seifried, S. Walch

Abstract: We take advantage of a set of molecular cloud simulations to demonstrate a possibility to uncover statistical properties of the gas density and velocity fields using reflected emission of a short (with duration much less than the cloud's light-crossing time) X-ray flare. Such situation is relevant for the Central Molecular Zone of our Galaxy where several clouds get illuminated by a $\sim110$ yr-o… ▽ More We take advantage of a set of molecular cloud simulations to demonstrate a possibility to uncover statistical properties of the gas density and velocity fields using reflected emission of a short (with duration much less than the cloud's light-crossing time) X-ray flare. Such situation is relevant for the Central Molecular Zone of our Galaxy where several clouds get illuminated by a $\sim110$ yr-old flare from the supermassive black hole Sgr A*. Due to shortness of the flare ($Δt\lesssim1.6$ yrs), only a thin slice ($Δz\lesssim0.5$ pc) of the molecular gas contributes to the X-ray reflection signal at any given moment, and its surface brightness effectively probes the local gas density. This allows reconstructing the density probability distribution function over a broad range of scales with virtually no influence of attenuation, chemo-dynamical biases and projection effects. Such measurement is key to understanding the structure and star-formation potential of the clouds evolving under extreme conditions in the CMZ. For cloud parameters similar to the currently brightest in X-ray reflection molecular complex Sgr A, the sensitivity level of the best available data is sufficient only for marginal distinction between solenoidal and compressive forcing of turbulence. Future-generation X-ray observatories with large effective area and high spectral resolution will dramatically improve on that by minimising systematic uncertainties due to contaminating signals. Furthermore, measurement of the iron fluorescent line centroid with sub-eV accuracy in combination with the data on molecular line emission will allow direct investigation of the gas velocity field. △ Less

Submitted 4 May, 2020; v1 submitted 27 June, 2019; originally announced June 2019.

Comments: 21 pages, 18 figures, MNRAS, accepted

SILCC-Zoom: H$_2$ and CO-dark gas in molecular clouds -- The impact of feedback and magnetic fields

Authors: D. Seifried, S. Haid, S. Walch, E. M. A. Borchert, T. G. Bisbas

Abstract: We analyse the CO-dark molecular gas content of simulated molecular clouds from the SILCC-Zoom project. The simulations reach a resolution of 0.1 pc and include H$_2$ and CO formation, radiative stellar feedback and magnetic fields. CO-dark gas is found in regions with local visual extinctions $A_\text{V, 3D} \sim$ 0.2 - 1.5, number densities of 10 - 10$^{3}$ cm$^{-3}$ and gas temperatures of few… ▽ More We analyse the CO-dark molecular gas content of simulated molecular clouds from the SILCC-Zoom project. The simulations reach a resolution of 0.1 pc and include H$_2$ and CO formation, radiative stellar feedback and magnetic fields. CO-dark gas is found in regions with local visual extinctions $A_\text{V, 3D} \sim$ 0.2 - 1.5, number densities of 10 - 10$^{3}$ cm$^{-3}$ and gas temperatures of few 10 K - 100 K. CO-bright gas is found at number densities above 300 cm$^{-3}$ and temperatures below 50 K. The CO-dark gas fractions range from 40% to 95% and scale inversely with the amount of well-shielded gas ($A_\text{V, 3D}$ $\gtrsim$ 1.5), which is smaller in magnetised molecular clouds. We show that the density, chemical abundances and $A_\text{V, 3D}$ along a given line-of-sight cannot be properly determined from projected quantities. As an example, pixels with a projected visual extinction of $A_\text{V, 2D} \simeq$ 2.5 - 5 can be both, CO-bright or CO-dark, which can be attributed to the presence or absence of strong density enhancements along the line-of-sight. By producing synthetic CO(1-0) emission maps of the simulations with RADMC-3D, we show that about 15 - 65\% of the H$_2$ is in regions with intensities below the detection limit. Our clouds have $X_\text{CO}$-factors around 1.5 $\times$ 10$^{20}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$ with a spread of up to a factor $\sim$ 4, implying a similar uncertainty in the derived total H$_2$ masses and even worse for individual pixels. Based on our results, we suggest a new approach to determine the H$_2$ mass, which relies on the availability of CO(1-0) emission and $A_\text{V, 2D}$ maps. It reduces the uncertainty of the clouds' overall H$_2$ mass to a factor of $\lesssim$ 1.8 and for individual pixels, i.e. on sub-pc scales, to a factor of $\lesssim$ 3. △ Less

Submitted 17 December, 2019; v1 submitted 3 June, 2019; originally announced June 2019.

Comments: 20 pages, 20 figures, accepted for publication in MNRAS

Astro2020: The Cycling of Matter from the Interstellar Medium to Stars and back

Authors: Robert Simon, Nicola Schneider, Frank Bigiel, Volker Ossenkopf-Okada, Yoko Okada, Doug Johnstone, Peter Schilke, Gordon Stacey, Markus Röllig, Alvaro Sanchez-Monge, Daniel Seifried, Juergen Stutzki, Frank Bertoldi, Christof Buchbender, Michel Fich, Terry Herter, Ronan Higgins, Thomas Nikola

Abstract: Understanding the matter cycle in the interstellar medium of galaxies from the assembly of clouds to star formation and stellar feedback remains an important and exciting field in comtemporary astrophysics. Many open questions regarding cloud and structure formation, the role of turbulence, and the relative importance of the various feedback processes can only be addressed with observations of spe… ▽ More Understanding the matter cycle in the interstellar medium of galaxies from the assembly of clouds to star formation and stellar feedback remains an important and exciting field in comtemporary astrophysics. Many open questions regarding cloud and structure formation, the role of turbulence, and the relative importance of the various feedback processes can only be addressed with observations of spectrally resolved lines. We here stress the importance of two specific sets of lines: the finestructure lines of atomic carbon as a tracer of the dark molecular gas and mid-J CO lines as tracers of the warm, active molecular gas in regions of turbulence dissipation and feedback. The observations must cover a wide range of environments (i.e., physical conditions), which will be achieved by large scale surveys of Galactic molecular clouds, the Galactic Center, the Magellanic clouds, and nearby galaxies. To date, such surveys are completely missing and thus constitute an important science opportunity for the next decade and beyond. For the successful interpretation of the observations, it will be essential to combine them with results from (chemical) modelling and simulations of the interstellar medium. △ Less

Submitted 26 March, 2019; v1 submitted 14 March, 2019; originally announced March 2019.

Comments: 7 pages, 1 figure

On the resolution requirements for modelling molecular gas formation in solar neighbourhood conditions

Authors: P. R. Joshi, S. Walch, D. Seifried, S. C. O. Glover, S. D. Clarke, M. Weis

Abstract: The formation of molecular hydrogen (H$_2$) and carbon monoxide (CO) is sensitive to the volume and column density distribution of the turbulent interstellar medium. In this paper, we study H$_2$ and CO formation in a large set of hydrodynamical simulations of periodic boxes with driven supersonic turbulence, as well as in colliding flows with the \textsc{Flash} code. The simulations include a non… ▽ More The formation of molecular hydrogen (H$_2$) and carbon monoxide (CO) is sensitive to the volume and column density distribution of the turbulent interstellar medium. In this paper, we study H$_2$ and CO formation in a large set of hydrodynamical simulations of periodic boxes with driven supersonic turbulence, as well as in colliding flows with the \textsc{Flash} code. The simulations include a non-equilibrium chemistry network, gas self-gravity, and diffuse radiative transfer. We investigate the spatial resolution required to obtain a converged H$_2$ and CO mass fraction and formation history. From the numerical tests we find that H$_2$ converges at a spatial resolution of $\lesssim0.2$~pc, while the required resolution for CO convergence is $\lesssim 0.04$~pc in gas with solar metallicity which is subject to a solar neighbourhood interstellar radiation field. We derive two critical conditions from our numerical results: the simulation has to at least resolve the densities at which (1) the molecule formation time in each cell in the computational domain is equal to the dissociation time, and (2) the formation time is equal to the the typical cell crossing time. For both H$_2$ and CO, the second criterion is more restrictive. The formulae we derive can be used to check whether molecule formation is converged in any given simulation. △ Less

Submitted 10 December, 2018; originally announced December 2018.

Comments: 22 pages, 21 figures, submitted

Evolution of Hubble wedges in episodic protostellar outflows

Authors: P. F. Rohde, S. Walch, D. Seifried, A. P. Whitworth, S. D. Clarke

Abstract: Young low-mass protostars undergo short phases of high accretion and outburst activity leading to lumpy outflows. Recent observations have shown that the position-velocity and mass-velocity diagrams of such outflows exhibit individual bullet-like features; some of these bullets subscribe to a `Hubble Law' velocity relation, and others are manifest as `Hubble wedges'. In order to explore the origin… ▽ More Young low-mass protostars undergo short phases of high accretion and outburst activity leading to lumpy outflows. Recent observations have shown that the position-velocity and mass-velocity diagrams of such outflows exhibit individual bullet-like features; some of these bullets subscribe to a `Hubble Law' velocity relation, and others are manifest as `Hubble wedges'. In order to explore the origin of these features, we have developed a new episodic outflow model for the SPH code {\sc gandalf}, which mimics the accretion and ejection behaviour of FU Ori type stars. We apply this model to simulations of star formation, invoking two types of initial conditions: spherically symmetric cores in solid-body rotation with $ρ\propto r^{-2}$, and spherically symmetric turbulent cores with density proportional to the density of a Bonnor-Ebert sphere. For a wide range of model parameters, we find that episodic outflows lead to self-regulation of the ejected mass and momentum, and we achieve acceptable results, even with relatively low resolution. Using this model, we find that recently ejected outflow bullets produce a `Hubble wedge' in the position-velocity relation. However, once such a bullet hits the leading shock front, it decelerates and aligns with older bullets to form a `Hubble-law'. Bullets can be identified as bumps in the mass-velocity relation, which can be fit with a power-law, $dM/d\upsilon_{_{\rm RAD}}\propto\upsilon_{_{\rm RAD}}^{-1.5}$. △ Less

Submitted 10 December, 2018; v1 submitted 8 November, 2018; originally announced November 2018.

Comments: Accepted (MNRAS) - 2018 November 30. 19 pages, 10 figures

MSC Class: 85-08

SILCC-ZOOM: The early impact of ionizing radiation on forming molecular clouds

Authors: Sebastian Haid, Stefanie Walch, Daniel Seifried, Richard Wuensch, Frantisek Dinnbier, Thorsten Naab

Abstract: As part of the SILCC-ZOOM project we present our first sub-parsec resolution radiation-hydrodynamic simulations of two molecular clouds self-consistently forming from a turbulent, multi-phase ISM. The clouds have similar initial masses of few 10$^4$ M$_{\odot}$, escape velocities of ~5 km s$^{-1}$, and a similar initial energy budget. We follow the formation of star clusters with a sink based mode… ▽ More As part of the SILCC-ZOOM project we present our first sub-parsec resolution radiation-hydrodynamic simulations of two molecular clouds self-consistently forming from a turbulent, multi-phase ISM. The clouds have similar initial masses of few 10$^4$ M$_{\odot}$, 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. Photo-ionizing 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 effciency is considerably reduced by a factor of ~4 to global cloud values of < 10 % as the mass accretion of sinks that host massive stars is terminated after <1 Myr. Despite the low effciency, 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 substructure and in particular on the amount of gas at high visual extinction. The damage of radiation done to the highly shielded cloud (MC1) 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 effciency of molecular clouds △ Less

Submitted 10 December, 2018; v1 submitted 18 October, 2018; originally announced October 2018.

Synthetic [CII] emission maps of a simulated molecular cloud in formation

Authors: A. Franeck, S. Walch, D. Seifried, S. D. Clarke, V. Ossenkopf-Okada, S. C. O. Glover, R. S. Klessen, P. Girichidis, T. Naab, R. Wünsch, P. C. Clark, E. Pellegrini, T. Peters

Abstract: The C$^{+}$ ion is an important coolant of interstellar gas, and so the [CII] fine structure line is frequently observed in the interstellar medium. However, the physical and chemical properties of the [CII]-emitting gas are still unclear. We carry out non-LTE radiative transfer simulations with RADMC-3D to study the [CII] line emission from a young, turbulent molecular cloud before the onset of s… ▽ More The C$^{+}$ ion is an important coolant of interstellar gas, and so the [CII] fine structure line is frequently observed in the interstellar medium. However, the physical and chemical properties of the [CII]-emitting gas are still unclear. We carry out non-LTE radiative transfer simulations with RADMC-3D to study the [CII] line emission from a young, turbulent molecular cloud before the onset of star formation, using data from the SILCC-Zoom project. The [CII] emission is optically thick over 40% of the observable area with $I_{[\textrm{CII}]} > 0.5$ K km s$^{-1}$. To determine the physical properties of the [CII] emitting gas, we treat the [CII] emission as optically thin. We find that the [CII] emission originates primarily from cold, moderate density gas ($40 \lesssim T \lesssim 65$ K and $50 \lesssim n \lesssim 440$ cm$^{-3}$), composed mainly of atomic hydrogen and with an effective visual extinction between $\sim 0.50$ and $\sim 0.91$. Gas dominated by molecular hydrogen contributes only $\lesssim$20% of the total [CII] line emission. Thus, [CII] is not a good tracer for CO-dark H$_2$ at this early phase in the cloud's lifetime. We also find that the total gas, H and C$^+$ column densities are all correlated with the integrated [CII] line emission, with power law slopes ranging from 0.5 to 0.7. Further, the median ratio between the total column density and the [CII] line emission is $Y_{\rm CII}\approx 1.1 \times 10^{21}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$, and $Y_{\rm CII}$ scales with $I_{[\textrm{CII}]}^{-0.3}$. We expect $Y_{\rm CII}$ to change in environments with a lower or higher radiation field than simulated here. △ Less

Submitted 27 September, 2018; originally announced September 2018.

Comments: 27 pages, 23 figures, accepted for publication in MNRAS

The SILCC project - V. The impact of magnetic fields on the chemistry and the formation of molecular clouds

Authors: Philipp Girichidis, Daniel Seifried, Thorsten Naab, Thomas Peters, Stefanie Walch, Richard Wünsch, Simon C. O. Glover, Ralf S. Klessen

Abstract: Magnetic fields are ubiquitously observed in the interstellar medium (ISM) of present-day star-forming galaxies with dynamically relevant energy densities. Using three-dimensional magneto-hydrodynamic (MHD) simulations of the supernova (SN) driven ISM in the flux-freezing approximation (ideal MHD) we investigate the impact of the magnetic field on the chemical and dynamical evolution of the gas, f… ▽ More Magnetic fields are ubiquitously observed in the interstellar medium (ISM) of present-day star-forming galaxies with dynamically relevant energy densities. Using three-dimensional magneto-hydrodynamic (MHD) simulations of the supernova (SN) driven ISM in the flux-freezing approximation (ideal MHD) we investigate the impact of the magnetic field on the chemical and dynamical evolution of the gas, fragmentation and the formation of molecular clouds. We follow the chemistry with a network of six species (H$^{+}$, H, H$_2$, C$^+$, CO, free electrons) including local shielding effects. We find that magnetic fields thicken the disc by a factor of a few to a scale height of $\sim100\,\mathrm{pc}$, delay the formation of dense (and molecular) gas by $\sim25\,\mathrm{Myr}$ and result in differently shaped gas structures. The magnetised gas fragments into fewer clumps, which are initially at subcritical mass-to-flux ratios, $M/Φ\approx0.3(M/Φ)_\mathrm{crit}$, and accrete gas preferentially parallel to the magnetic field lines until supercritial mass-to-flux ratios of up to order 10 are reached. The accretion rates onto molecular clouds scale with $\dot{M}\propto M^{1.5}$. The median of the inter-cloud velocity dispersion is $\sim2-5\,\mathrm{km\,s}^{-1}$ and lower than the internal velocity dispersion in the clouds ($\sim3-7\,\mathrm{km\,s}^{-1}$). However, individual cloud-cloud collisions occur at speeds of a few $10\,\mathrm{km\,s}^{-1}$. △ Less

Submitted 15 August, 2018; originally announced August 2018.

Comments: 31 pages, 26 figures, MNRAS accepted

The relative impact of photoionizing radiation and stellar winds on different environments

Authors: Sebastian Haid, Stefanie Walch, Daniel Seifried, Richard Wünsch, František Dinnbier, Thorsten Naab

Abstract: Photoionizing radiation and stellar winds from massive stars deposit energy and momentum into the interstellar medium (ISM). They might disperse the local ISM, change its turbulent multi-phase structure, and even regulate star formation. Ionizing radiation dominates the massive stars' energy output, but the relative effect of winds might change with stellar mass and the properties of the ambient I… ▽ More Photoionizing radiation and stellar winds from massive stars deposit energy and momentum into the interstellar medium (ISM). They might disperse the local ISM, change its turbulent multi-phase structure, and even regulate star formation. Ionizing radiation dominates the massive stars' energy output, but the relative effect of winds might change with stellar mass and the properties of the ambient ISM. We present simulations of the interaction of stellar winds and ionizing radiation of 12, 23, and 60 M$_{\odot}$ stars within a cold neutral (CNM, $n_{0}$ = 100 cm$^{-3}$), warm neutral (WNM, $n_{0}$ = 1, 10 cm$^{-3}$) or warm ionized (WIM, $n_{0}$ = 0.1 cm$^{-3}$) medium. The FLASH simulations adopt the novel tree-based radiation transfer algorithm TreeRay. With the On-the-Spot approximation and a temperature-dependent recombination coefficient, it is coupled to a chemical network with radiative heating and cooling. In the homogeneous CNM, the total momentum injection ranges from 1.6$\times$10$^{4}$ to 4$\times$10$^{5}$ M$_{\odot}$ km s$^{-1}$ and is always dominated by the expansion of the ionized H$_{\text{II}}$ region. In the WIM, stellar winds dominate (2$\times$10$^{2}$ to 5$\times$10$^{3}$ M$_{\odot}$ km s$^{-1}$), while the input from radiation is small ($\sim$ 10$^{2}$ M$_{\odot}$ km s$^{-1}$). The WNM ($n_{0}$ = 1 cm$^{-3}$) is a transition regime. Energetically, stellar winds couple more efficiently to the ISM ($\sim$ 0.1 percent of wind luminosity) than radiation ($<$ 0.001 percent of ionizing luminosity). For estimating the impact of massive stars, the strongly mass-dependent ratios of wind to ionizing luminosity and the properties of the ambient medium have to be considered. △ Less

Submitted 26 April, 2018; originally announced April 2018.

SILCC-Zoom: Polarisation and depolarisation in molecular clouds

Authors: D. Seifried, S. Walch, S. Reissl, J. C. Ibáñez-Mejía

Abstract: We present synthetic dust polarisation maps of 3D magneto-hydrodynamical simulations of molecular clouds before the onset of stellar feedback. The clouds are modelled within the SILCC-Zoom project and are embedded in their galactic environment. The radiative transfer is carried out with POLARIS for wavelengths from 70 $μ$m to 3 mm at a resolution of 0.12 pc, and includes self-consistently calculat… ▽ More We present synthetic dust polarisation maps of 3D magneto-hydrodynamical simulations of molecular clouds before the onset of stellar feedback. The clouds are modelled within the SILCC-Zoom project and are embedded in their galactic environment. The radiative transfer is carried out with POLARIS for wavelengths from 70 $μ$m to 3 mm at a resolution of 0.12 pc, and includes self-consistently calculated alignment efficiencies for radiative torque alignment. We explore the reason of the observed depolarisation in the center of molecular clouds: We find that dust grains remain well aligned even at high densities ($n$ $>$ 10$^3$ cm$^{-3}$) and visual extinctions ($A_\text{V}$ $>$ 1). The depolarisation is rather caused by strong variations of the magnetic field direction along the LOS due to turbulent motions. The observed magnetic field structure thus resembles best the mass-weighted, line-of-sight averaged field structure. Furthermore, it differs by only a few 1$^\circ$ for different wavelengths and is little affected by the spatial resolution of the synthetic observations. Noise effects can be reduced by convolving the image. Doing so, for $λ$ $\gtrsim$ 160 $μ$m the observed magnetic field traces reliably the underlying field in regions with intensities $I$ $\gtrsim$ 2 times the noise level and column densities above 1 M$_\text{sun}$ pc$^{-2}$. Here, typical deviations are $\lesssim$ 10$^\circ$. The observed structure is less reliable in regions with low polarisation degrees and possibly in regions with large column density gradients. Finally, we show that a simplified and widely used method without self-consistent dust alignment efficiencies can provide a good representation of the observable polarisation structure with deviations below 5$^\circ$. △ Less

Submitted 18 October, 2018; v1 submitted 26 April, 2018; originally announced April 2018.

Comments: 21 pages including appendix, 20 figures, accepted for publication in MNRAS

Non-Equilibrium Chemistry and Destruction of CO by X-ray Flares

Authors: Jonathan Mackey, Stefanie Walch, Daniel Seifried, Simon C. O. Glover, Richard Wünsch, Felix Aharonian

Abstract: Sources of X-rays such as active galactic nuclei and X-ray binaries are often variable by orders of magnitude in luminosity over timescales of years. During and after these flares 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 magnetohydrodynamical simulations.… ▽ More Sources of X-rays such as active galactic nuclei and X-ray binaries are often variable by orders of magnitude in luminosity over timescales of years. During and after these flares 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 magnetohydrodynamical simulations. A static fractal molecular cloud is irradiated with X-rays of different intensity, and the chemical and thermal evolution of the cloud are studied. For a simulated $10^5$ M$_\odot$ fractal cloud an X-ray flux $<0.01$ erg cm$^{-2}$ s$^{-1}$ allows the cloud to remain molecular, whereas most of the CO and H$_2$ are destroyed for a flux of $>1$ erg cm$^{-2}$ s$^{-1}$. The effects of an X-ray flare, which suddenly increases the X-ray flux by $10^5 \times$ are then studied. A cloud exposed to a bright flare has 99% of its CO destroyed in 10-20 years, whereas it takes $>10^3$ years for 99% of the H$_2$ to be destroyed. CO is primarily destroyed by locally generated far-UV emission from collisions between non-thermal electrons and H$_2$; He$^+$ only becomes an important destruction agent when the CO abundance is already very small. After the flare is over, CO re-forms and approaches its equilibrium abundance after $10^3-10^5$ years. This implies that molecular clouds close to Sgr A$^*$ in the Galactic Centre may still be out of chemical equilibrium, and we predict the existence of clouds near flaring X-ray sources in which CO has been mostly destroyed but H is fully molecular. △ Less

Submitted 25 March, 2019; v1 submitted 27 March, 2018; originally announced March 2018.

Comments: Accepted for publication in MNRAS; this version has some additions following the refereeing process

Is molecular cloud turbulence driven by external supernova explosions?

Authors: D. Seifried, S. Walch, S. Haid, P. Girichidis, T. Naab

Abstract: We present high-resolution ($\sim$ 0.1 pc), hydrodynamical and magnetohydrodynamical simulations to investigate whether the observed level of molecular cloud (MC) turbulence can be generated and maintained by external supernova (SN) explosions. The MCs are formed self-consistently within their large-scale galactic environment following the non-equilibrium formation of H$_2$ and CO including (self-… ▽ More We present high-resolution ($\sim$ 0.1 pc), hydrodynamical and magnetohydrodynamical simulations to investigate whether the observed level of molecular cloud (MC) turbulence can be generated and maintained by external supernova (SN) explosions. The MCs are formed self-consistently within their large-scale galactic environment following the non-equilibrium formation of H$_2$ and CO including (self-) shielding and important heating and cooling processes. The MCs inherit their initial level of turbulence from the diffuse ISM, where turbulence is injected by SN explosions. However, by systematically exploring the effect of individual SNe going off outside the clouds, we show that at later stages the importance of SN driven turbulence is decreased significantly. This holds for different MC masses as well as MCs with and without magnetic fields. The SN impact also decreases rapidly with larger distances. Nearby SNe ($d$ $\sim$ 25 pc) boost the turbulent velocity dispersions of the MC by up to 70 per cent (up to a few km s$^{-1}$). For $d$ $>$ 50 pc, however, their impact decreases fast with increasing $d$ and is almost negligible. For all probed distances the gain in velocity dispersion decays rapidly within a few 100 kyr. This is significantly shorter than the average timescale for a MC to be hit by a nearby SN under solar neighbourhood conditions ($\sim$ 2 Myr). Hence, at these conditions SNe are not able to sustain the observed level of MC turbulence. However, in environments with high gas surface densities and SN rates like the Central Molecular Zone, observed elevated MC dispersions could be triggered by external SNe. △ Less

Submitted 5 March, 2018; v1 submitted 3 February, 2018; originally announced February 2018.

Comments: 10 pages, 11 figures, accepted for publication in ApJ, updated to final version

SILCC-Zoom: The dynamical and chemical evolution of molecular clouds

Authors: D. Seifried, S. Walch, P. Girichidis, T. Naab, R. Wünsch, R. S. Klessen, S. C. O. Glover, T. Peters, P. Clark

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 H$_2$ and CO including (self-) shielding. The two clouds are assem… ▽ More 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 H$_2$ and CO including (self-) shielding. The two clouds are assembled within a few million years with mass growth rates of up to $\sim$ 10$^{-2}$ M$_\mathrm{sun}$ yr$^{-1}$ and final masses of $\sim$ 50 000 M$_\mathrm{sun}$. A spatial resolution of $\lesssim$ 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, H$_2$ 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 $\gtrsim$ 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 \geq$ 300 cm$^{-3}$, H$_2$ is also found for $n \lesssim 30$ cm$^{-3}$ due to turbulent mixing and becomes dominant at column densities around 30 - 50 M$_\mathrm{sun}$ pc$^{-2}$. The CO-to-H$_2$ ratio steadily increases within the first 2 Myr whereas $X_\mathrm{CO}\simeq$ 1 - 4 $\times$ 10$^{20}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$ is approximately constant since the CO(1-0) line quickly becomes optically thick. △ Less

Submitted 6 September, 2017; v1 submitted 21 April, 2017; originally announced April 2017.

Comments: 24 pages, 22 figures, accepted for publication in MNRAS

The origin of dust polarization in molecular outflows

Authors: S. Reissl, D. Seifried, S. Wolf, R. Banerjee, R. S. Klessen

Abstract: Aims: In this paper we present a case study to investigate conditions necessary to detect a characteristic magnetic field substructure embedded in a large-scale field. A helical magnetic field with a surrounding hourglass shaped field is expected from theoretical predictions and self-consistent magnetohydrodynamical (MHD) simulations to be present in the specific case of protostellar outflows. Hen… ▽ More Aims: In this paper we present a case study to investigate conditions necessary to detect a characteristic magnetic field substructure embedded in a large-scale field. A helical magnetic field with a surrounding hourglass shaped field is expected from theoretical predictions and self-consistent magnetohydrodynamical (MHD) simulations to be present in the specific case of protostellar outflows. Hence, such an outflow environment is the perfect for our study. Methodes: We present synthetic polarisation maps in the infrared and millimeter regime of protostellar outflows performed with the newly developed RT and polarisation code POLARIS. The code, as the first, includes a self-consistent description of various alignement mechanism like the imperfect Davis-Greenstein (IDG) and the radiative torque (RAT) alignment. We investigate for which effects the grain size distribution, and applied alignement mechanism have. Results: We find that the IDG mechanism cannot produce any measurable polarization degree (< 1 %) whereas RAT alignment produced polarization degrees of a few 1 %. Furthermore, we developed a method to identify the origin of the polarization. We show that the helical magnetic field in the outflow can only be observed close to the outflow axis and at its tip, whereas in the surrounding regions the hourglass field in the foreground dominates the polarization. Furthermore, the polarization degree in the outflow lobe is lower than in the surroundings in agreement with observations. We also find that the orientation of the polarization vector flips around a few 100 micron due to the transition from dichroic extinction to thermal re-emission. Hence, in order to avoid ambiguities when interpreting polarization data, we suggest to observed in the far-infrared and mm regime. Finally, we show that with ALMA it is possible to observe the polarization emerging from protostellar outflows. △ Less

Submitted 8 March, 2017; originally announced March 2017.

Comments: 15 pages, 9 figures

Journal ref: A&A 603, A71 (2017)

Modelling the chemistry of star-forming filaments - II. Testing filament characteristics with synthetic observations

Authors: D. Seifried, Á. Sánchez-Monge, S. Suri, S. Walch

Abstract: We present synthetic continuum and $^{13}$CO and C$^{18}$O line emission observations of dense and cold filaments. The filaments are dynamically evolved using 3D-MHD simulations that include one of the largest on-the-fly chemical networks used to date, which models the detailed evolution of H$_2$ and CO. We investigate the reliability of observable properties, in particular filament mass and width… ▽ More We present synthetic continuum and $^{13}$CO and C$^{18}$O line emission observations of dense and cold filaments. The filaments are dynamically evolved using 3D-MHD simulations that include one of the largest on-the-fly chemical networks used to date, which models the detailed evolution of H$_2$ and CO. We investigate the reliability of observable properties, in particular filament mass and width, under different simulation conditions like magnetic field orientation and cosmic ray ionisation rate. We find that filament widths of $\sim$0.1 pc can be probed with both line and continuum emission observations with a high accuracy (deviations $\leq$ 20%). However, the width of more narrow filaments can be significantly overestimated by up to a factor of a few. Masses obtained via the dust emission are accurate within a few percent whereas the masses inferred from molecular line emission observations deviate from the actual mass by up to a factor of 10 and show large differences for different $J$ transitions. The inaccurate estimate of filament masses and widths of narrow filaments using molecular line observations can be attributed to (i) the non-isothermal state of the filaments, (ii) optical depth effects, and (iii) the subthermally excited state of CO, while inclination effects and opacity correction only influence the obtained masses and widths by less than 50%. Both, mass and width estimates, can be improved by using two isotopes to correct for the optical depth. Since gas and dust temperature generally differ (by up to 25 K), the filaments appear more gravitationally unstable if the (too low) dust temperature is used for the stability analysis. △ Less

Submitted 1 March, 2017; v1 submitted 10 February, 2017; originally announced February 2017.

Comments: 18 pages, 14 figures, accepted by MNRAS, updated to final version

The turbulent life of dust grains in the supernova-driven, multi-phase interstellar medium

Authors: Thomas Peters, Svitlana Zhukovska, Thorsten Naab, Philipp Girichidis, Stefanie Walch, Simon C. O. Glover, Ralf S. Klessen, Paul C. Clark, Daniel Seifried

Abstract: Dust grains are an important component of the interstellar medium (ISM) of galaxies. We present the first direct measurement of the residence times of interstellar dust in the different ISM phases, and of the transition rates between these phases, in realistic hydrodynamical simulations of the multi-phase ISM. Our simulations include a time-dependent chemical network that follows the abundances of… ▽ More Dust grains are an important component of the interstellar medium (ISM) of galaxies. We present the first direct measurement of the residence times of interstellar dust in the different ISM phases, and of the transition rates between these phases, in realistic hydrodynamical simulations of the multi-phase ISM. Our simulations include a time-dependent chemical network that follows the abundances of H^+, H, H_2, C^+ and CO and take into account self-shielding by gas and dust using a tree-based radiation transfer method. Supernova explosions are injected either at random locations, at density peaks, or as a mixture of the two. For each simulation, we investigate how matter circulates between the ISM phases and find more sizeable transitions than considered in simple mass exchange schemes in the literature. The derived residence times in the ISM phases are characterised by broad distributions, in particular for the molecular, warm and hot medium. The most realistic simulations with random and mixed driving have median residence times in the molecular, cold, warm and hot phase around 17, 7, 44 and 1 Myr, respectively. The transition rates measured in the random driving run are in good agreement with observations of Ti gas-phase depletion in the warm and cold phases in a simple depletion model, although the depletion in the molecular phase is under-predicted. ISM phase definitions based on chemical abundance rather than temperature cuts are physically more meaningful, but lead to significantly different transition rates and residence times because there is no direct correspondence between the two definitions. △ Less

Submitted 20 October, 2016; originally announced October 2016.

Comments: submitted to MNRAS, movies https://hera.ph1.uni-koeln.de/~silcc/

Journal ref: Monthly Notices of the Royal Astronomical Society 467 (2017), 4322-4342

Supernova-blast waves in wind-blown bubbles, turbulent, and power-law ambient media

Authors: Sebastian Haid, Stefanie Walch, Thorsten Naab, Daniel Seifried, Jonathan Mackey, Andrea Gatto

Abstract: Supernova (SN) blast waves inject energy and momentum into the interstellar medium (ISM), control its turbulent multiphase structure and the launching of galactic outflows. Accurate modelling of the blast wave evolution is therefore essential for ISM and galaxy formation simulations. We present an efficient method to compute the input of momentum, thermal energy, and the velocity distribution of t… ▽ More Supernova (SN) blast waves inject energy and momentum into the interstellar medium (ISM), control its turbulent multiphase structure and the launching of galactic outflows. Accurate modelling of the blast wave evolution is therefore essential for ISM and galaxy formation simulations. We present an efficient method to compute the input of momentum, thermal energy, and the velocity distribution of the shock-accelerated gas for ambient media with uniform (and with stellar wind blown bubbles), power-law, and turbulent density distributions. Assuming solar metallicity cooling, the blast wave evolution is followed to the beginning of the momentum conserving snowplough phase. The model recovers previous results for uniform ambient media. The momentum injection in wind-blown bubbles depend on the swept-up mass and the efficiency of cooling, when the blast wave hits the wind shell. For power-law density distributions with $n(r) \sim$ $r^{-2}$ (for $n(r) > n_{_{\rm floor}}$) the amount of momentum injection is solely regulated by the background density $n_{_{\rm floor}}$ and compares to $n_{_{\rm uni}}$ = $n_{_{\rm floor}}$. However, in turbulent ambient media with log-normal density distributions the momentum input can increase by a factor of 2 (compared to the homogeneous case) for high Mach numbers. The average momentum boost can be approximated as $p_{_{\rm turb}}/\mathrm{p_{_{0}}}\ =23.07\, \left(\frac{n_{_{0,\rm turb}}}{1\,{\rm cm}^{-3}}\right)^{-0.12} + 0.82 (\ln(1+b^{2}\mathcal{M}^{2}))^{1.49}\left(\frac{n_{_{0,\rm turb}}}{1\,{\rm cm}^{-3}}\right)^{-1.6}$. The velocity distributions are broad as gas can be accelerated to high velocities in low-density channels. The model values agree with results from recent, computationally expensive, three-dimensional simulations of SN explosions in turbulent media. △ Less

Submitted 15 April, 2016; originally announced April 2016.

Supernova Feedback in Molecular Clouds: Global Evolution and Dynamics

Authors: Bastian Körtgen, Daniel Seifried, Robi Banerjee, Enrique Vázquez-Semadeni, Manuel Zamora-Avilés

Abstract: We use magnetohydrodynamical simulations of converging warm neutral medium flows to analyse the formation and global evolution of magnetised and turbulent molecular clouds subject to supernova feedback from massive stars. We show that supernova feedback alone fails to disrupt entire, gravitationally bound, molecular clouds, but is able to disperse small--sized (~10 pc) regions on timescales of les… ▽ More We use magnetohydrodynamical simulations of converging warm neutral medium flows to analyse the formation and global evolution of magnetised and turbulent molecular clouds subject to supernova feedback from massive stars. We show that supernova feedback alone fails to disrupt entire, gravitationally bound, molecular clouds, but is able to disperse small--sized (~10 pc) regions on timescales of less than 1 Myr. Efficient radiative cooling of the supernova remnant as well as strong compression of the surrounding gas result in non-persistent energy and momentum input from the supernovae. However, if the time between subsequent supernovae is short and they are clustered, large hot bubbles form that disperse larger regions of the parental cloud. On longer timescales, supernova feedback increases the amount of gas with moderate temperatures (T~300-3000 K). Despite its inability to disrupt molecular clouds, supernova feedback leaves a strong imprint on the star formation process. We find an overall reduction of the star formation efficiency by a factor of 2 and of the star formation rate by roughly factors of 2-4. △ Less

Submitted 31 March, 2016; originally announced March 2016.

Comments: 16 pages, 12 figures (2 in appendix), revised version, submitted to MNRAS