Showing 1–50 of 82 results for author: Wünsch, R

Radiation transport methods in star formation simulations

Authors: Richard Wünsch

Abstract: Radiation transport plays a crucial role in star formation models, as certain questions within this field cannot be accurately addressed without taking it into account. Given the high complexity of the interstellar medium from which stars form, numerical simulations are frequently employed to model the star formation process. This study reviews recent methods for incorporating radiation transport… ▽ More Radiation transport plays a crucial role in star formation models, as certain questions within this field cannot be accurately addressed without taking it into account. Given the high complexity of the interstellar medium from which stars form, numerical simulations are frequently employed to model the star formation process. This study reviews recent methods for incorporating radiation transport into star formation simulations, discussing them in terms of the used algorithms, treatment of radiation frequency dependence, the interaction of radiation with the gas, and the parallelization of methods for deployment on supercomputers. Broadly, the algorithms fall into two categories: (i) moment-based methods, encompassing the flux-limited diffusion approximation, M1 closure, and variable Eddington tensor methods, and (ii) methods directly solving the radiation transport equation, including forward and reverse ray tracing, characteristics-based methods, and Monte Carlo techniques. Beyond discussing advantages and disadvantages of these methods, the review also lists recent radiation hydrodynamic codes implemented the described methods. △ Less

Submitted 8 March, 2024; originally announced March 2024.

Comments: 22 pages, accepted by Frontiers in Astronomy and Space Sciences

Journal ref: Front. Astron. Space Sci., 01 March 2024, Sec. Stellar and Solar Physics, Volume 11 - 2024

Ionized regions in the central arcsecond of NGC 1068. YJHK spatially resolved spectroscopy

Authors: P. Vermot, B. Barna, S. Ehlerová, M. R. Morris, J. Palous, R. Wünsch

Abstract: Context. Several bright emission line regions have been observed in the central 100 parsecs of the active galaxy NGC 1068. Aims. We aim to determine the properties and ionization mechanism of three regions of NGC 1068: the nucleus (B) and two clouds located at 0.3" and 0.7" north of it (C and D). Methods. We combined SPHERE (0.95 - 1.65 um) and SINFONI (1.5 - 2.45 um) spectra for the three regions… ▽ More Context. Several bright emission line regions have been observed in the central 100 parsecs of the active galaxy NGC 1068. Aims. We aim to determine the properties and ionization mechanism of three regions of NGC 1068: the nucleus (B) and two clouds located at 0.3" and 0.7" north of it (C and D). Methods. We combined SPHERE (0.95 - 1.65 um) and SINFONI (1.5 - 2.45 um) spectra for the three regions B, C, and D. We compared these spectra to several CLOUDY photoionization models and to the MAPPINGS III Library of Fast Radiative Shock Models. Results. The emission line spectra of the three regions are almost identical to each other and contribute to most of the emission line flux in the nuclear region. The emitting media contain multiple phases, the most luminous of which have temperatures ranging from 104.8 K to 106 K. Central photoionization models can reproduce some features of the spectra, but the fast radiative shock model provides the best fit to the data. Conclusions. The similarity between the three regions indicates that they belong to the same class of objects. Based on our comparisons, we conclude that they are shock regions located where the jet of the active galactic nucleus impacts massive molecular clouds. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: A&A, Forthcoming article, accepted for publication

A 3D model for the stellar populations in the nuclei of NGC 1433,NGC 1566, and NGC 1808

Authors: P. Vermot, J. Palouš, B. Barna, S. Ehlerová, M. R. Morris, R. Wünsch

Abstract: Aims. We aim to characterize the properties of the stellar populations in the central few hundred parsecs of nearby galactic nuclei; specifically their age, mass, and 3D geometry. Methods. We use spatially resolved spectroscopic observations of NGC 1433, NGC 1566, and NGC 1808 obtained with SINFONI to constrain a 3D model composed of a spherically symmetric nuclear star cluster (NSC) and an extend… ▽ More Aims. We aim to characterize the properties of the stellar populations in the central few hundred parsecs of nearby galactic nuclei; specifically their age, mass, and 3D geometry. Methods. We use spatially resolved spectroscopic observations of NGC 1433, NGC 1566, and NGC 1808 obtained with SINFONI to constrain a 3D model composed of a spherically symmetric nuclear star cluster (NSC) and an extended thick stellar disk. We computed UV to mid-infrared single stellar population (UMISSP) spectra to determine the age of the stellar populations and construct synthetic observations for our model. To overcome degeneracies between key parameters, we simultaneously fit the spatially resolved line-of-sight velocity, line-of-sight-velocity-dispersion, low-spectral-resolution NIR continuum, and high-spectral-resolution CO absorption features for each pixel. Results. For the three objects, we derive the age and mass of the young and old stellar populations in the NSC and surrounding disk, as well as their 3D geometry: radius for the NSC; thickness, inclination, and position angle for the disk. These results are consistent with published independent measurements when available. Conclusions. The proposed method allows us to derive a consistent 3D model of the stellar populations in nearby galactic centers solely based on a near-infrared IFU observation. △ Less

Submitted 30 May, 2023; originally announced May 2023.

Journal ref: A&A 674, A135 (2023)

Simulations of pre-supernova feedback in spherical clouds

Authors: Michalis Kourniotis, Richard Wünsch, Sergio Martínez-González, Jan Palouš, Guillermo Tenorio-Tagle, Soňa Ehlerová

Abstract: We present a one-dimensional radiation-hydrodynamic model of a spherically symmetric cloud evolving under the influence of the self-gravity and the feedback from a star cluster forming in its centre. On one hand, the model is simple due to its 1D geometry, on the other hand, the feedback includes the ionising radiation, stellar winds and the radiation pressure acting on gas and dust. The star clus… ▽ More We present a one-dimensional radiation-hydrodynamic model of a spherically symmetric cloud evolving under the influence of the self-gravity and the feedback from a star cluster forming in its centre. On one hand, the model is simple due to its 1D geometry, on the other hand, the feedback includes the ionising radiation, stellar winds and the radiation pressure acting on gas and dust. The star cluster is formed from the gas flowing into the cloud centre and the feedback parameters are determined from stellar evolution models and the cluster star forming history. The model is compared to the semi-analytic code WARPFIELD implementing similar physical processes and exploring the scenario that the young cluster R136 in the Large Magellanic Cloud was formed due to re-collapse of the shell formed by the previous generation star cluster. A good qualitative agreement is found, however, $3 - 4$ times higher stellar mass is needed to disrupt the cloud in our model, because it takes into account (contrary to WARPFIELD) self-gravity of the cloud surrounding the shell. We use the model to explore star formation in clouds with different mass, radius and density profile measuring their star formation efficiency (SFE), i.e. the fraction of the cloud mass converted to stars. We found that SFE is a function of a single parameter, $\mathrm{log(SFE)} \propto -n_{hm}^{-0.46}$, with $n_{hm}$ being the cloud mean particle density within its half-mass radius. Furthermore, we found that the feedback efficiency, i.e. a fraction of the feedback energy retained by gas, has a nearly constant value $\sim 10^{-3}$. △ Less

Submitted 15 March, 2023; originally announced March 2023.

Comments: 15 pages including appendix, 11 figures. Accepted for publication in MNRAS

Tree-based solvers for adaptive mesh refinement code FLASH -- IV: An X-ray radiation scheme to couple discrete and diffuse X-ray emission sources to the thermochemistry of the interstellar medium

Authors: Brandt A. L. Gaches, Stefanie Walch, Richard Wünsch, Jonathan Mackey

Abstract: X-ray radiation, in particular radiation between 0.1 keV and 10 keV, is evident from both point-like sources, such as compact objects and T-Tauri young stellar objects, and extended emission from hot, cooling gas, such as in supernova remnants. The X-ray radiation is absorbed by nearby gas, providing a source of both heating and ionization. While protoplanetary chemistry models now often include X… ▽ More X-ray radiation, in particular radiation between 0.1 keV and 10 keV, is evident from both point-like sources, such as compact objects and T-Tauri young stellar objects, and extended emission from hot, cooling gas, such as in supernova remnants. The X-ray radiation is absorbed by nearby gas, providing a source of both heating and ionization. While protoplanetary chemistry models now often include X-ray emission from the central young stellar object, simulations of star-forming regions have yet to include X-ray emission coupled to the chemo-dynamical evolution of the gas. We present an extension of the {\sc TreeRay} reverse raytrace algorithm implemented in the {\sc Flash} magneto-hydrodynamic code which enables the inclusion of X-ray radiation from 0.1 keV $< E_γがんま <$ 100 keV, dubbed {\rm XrayTheSpot}. {\sc XrayTheSpot} allows for the use of an arbitrary number of bins, minimum and maximum energies, and both temperature-independent and temperature-dependent user-defined cross sections, along with the ability to include both point and extended diffuse emission and is coupled to the thermochemical evolution. We demonstrate the method with several multi-bin benchmarks testing the radiation transfer solution and coupling to the thermochemistry. Finally, we show two example star formation science cases for this module: X-ray emission from protostellar accretion irradiating an accretion disk and simulations of molecular clouds with active chemistry, radiation pressure, protostellar radiation feedback from infrared to X-ray radiation. △ Less

Submitted 29 April, 2023; v1 submitted 30 January, 2023; originally announced January 2023.

Comments: Accepted to MNRAS

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 Arecibo Galaxy Environment Survey XII : Optically dark HI clouds in the Leo I Group

Authors: Rhys Taylor, Joachim Koppen, Pavel Jachym, Robert Minchin, Jan Palous, Jessica Rosenberg, Steven Schneider, Richard Wunsch, Boris Deshev

Abstract: Using data from the Arecibo Galaxy Environment Survey, we report the discovery of five HI clouds in the Leo I group without detected optical counterparts. Three of the clouds are found midway between M96 and M95, one is only 10$^{\prime}$ from the south-east side of the well-known Leo Ring, and the fifth is relatively isolated. HI masses range from 2.6$\times$10$^{6}$ - 9.0$\times$10$^{6}$M… ▽ More Using data from the Arecibo Galaxy Environment Survey, we report the discovery of five HI clouds in the Leo I group without detected optical counterparts. Three of the clouds are found midway between M96 and M95, one is only 10$^{\prime}$ from the south-east side of the well-known Leo Ring, and the fifth is relatively isolated. HI masses range from 2.6$\times$10$^{6}$ - 9.0$\times$10$^{6}$M$_{\odot}$ and velocity widths (W50) from 16 - 42 km/s. Although a tidal origin is the most obvious explanation, this formation mechanism faces several challenges. For the most isolated cloud, the difficulties are its distance from neighbouring galaxies and the lack of any signs of disturbance in the HI discs of those systems. Some of the clouds also appear to follow the baryonic Tully-Fisher relation between mass and velocity width for normal, stable galaxies which is not expected if they are tidal in origin. Three clouds are found between M96 and M95 which have no optical counterparts, but have otherwise similar properties and location to the optically detected galaxy LeG 13. While overall we favour a tidal debris scenario to explain the clouds, we cannot rule out a primordial origin. If the clouds were produced in the same event that gave rise to the Leo Ring, they may provide important constraints on any model attempting to explain that structure △ Less

Submitted 22 September, 2022; originally announced September 2022.

Comments: 16 pages (plus 4 page appendix), 1 table, 7 figures, accepted for publication in AJ

Dust Grain Growth & Dusty Supernovae in Low-Metallicity Molecular Clouds

Authors: Sergio Martínez-González, Richard Wünsch, Guillermo Tenorio-Tagle, Sergiy Silich, Dorottya Szécsi, Jan Palouš

Abstract: We present 3-D hydrodynamical models of the evolution of superbubbles powered by stellar winds and supernovae from young coeval massive star clusters within low metallicity ($Z = 0.02$Z$_{\odot}$), clumpy molecular clouds. We explore the initial stages of the superbubble evolution, including the occurrence of pair-instability and core-collapse supernovae. Our aim is to study the occurrence of dust… ▽ More We present 3-D hydrodynamical models of the evolution of superbubbles powered by stellar winds and supernovae from young coeval massive star clusters within low metallicity ($Z = 0.02$Z$_{\odot}$), clumpy molecular clouds. We explore the initial stages of the superbubble evolution, including the occurrence of pair-instability and core-collapse supernovae. Our aim is to study the occurrence of dust grain growth within orbiting dusty clumps, and in the superbubble's swept-up supershell. We also aim to address the survival of dust grains produced by sequential supernovae. The model accounts for the star cluster gravitational potential and self-gravity of the parent cloud. It also considers radiative cooling (including that induced by dust) and a state-of-the-art population synthesis model for the coeval cluster. As shown before, a superbubble embedded into a clumpy medium becomes highly distorted, expanding mostly due to the hot gas streaming through low density channels. Our results indicate that in the case of massive ($\sim10^7$M$_{\odot}$) molecular clouds, hosting a super star cluster ($\sim5.6\times10^5$M$_{\odot}$), grain growth increments the dust mass at a rate $\sim4.8\times10^{-5}$M$_{\odot}$ yr$^{-1}$ during the first $2.5$Myr of the superbubble's evolution, while the net contribution of pair-instability and core-collapse supernovae to the superbubble's dust budget is $\sim1200$M$_{\odot} (M_{SC}/5.6\times10^{5}$M$_{\odot})$, where $M_{SC}$ is the stellar mass of the starburst. Therefore, dust grain growth and dust injection by supernovae lead to create, without invoking a top-heavy initial mass function, massive amounts of dust within low-metallicity star-forming molecular clouds, in accordance with the large dust mass present in galaxies soon after the onset of cosmic reionization. △ Less

Submitted 13 June, 2022; originally announced June 2022.

Comments: 8 pages, 4 figures. Accepted for publication in The Astrophysical Journal

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

X-ray emission from star cluster winds in starburst galaxies

Authors: Annika Franeck, Richard Wünsch, Sergio Martínez-González, Ivana Orlitová, Peter Boorman, Jiří Svoboda, Dorottya Szécsi, Vanesa Douna

Abstract: Inspired by the excess soft X-ray emission recently detected in Green Pea galaxies, we model the soft X-ray emission (0.5 - 2.0 keV) of hot gas from star cluster winds. By combining individual star clusters, we estimate the soft X-ray emission expected from the typically unresolved diffuse hot gas in starburst galaxies, devoid of competing emission from e.g., AGN or other unresolved point sources.… ▽ More Inspired by the excess soft X-ray emission recently detected in Green Pea galaxies, we model the soft X-ray emission (0.5 - 2.0 keV) of hot gas from star cluster winds. By combining individual star clusters, we estimate the soft X-ray emission expected from the typically unresolved diffuse hot gas in starburst galaxies, devoid of competing emission from e.g., AGN or other unresolved point sources. We use stellar models of sub-solar metallicities (0.02 $Z_{\odot}$ and 0.4 $Z_{\odot}$), and take into account supernova explosions for massive stars. For lower metallicities, we find that stellar winds do not contribute significantly ($\lesssim 3$ % of the mechanical energy) to the observed soft X-ray emission of normal star forming galaxies. For higher metallicities and possibly also for larger proportions of massive star clusters in the simulated starburst galaxies, we reproduce well the observed correlation between star formation rate and X-ray luminosity previously reported in the literature. However, we find that no combination of model assumptions is capable of reproducing the substantial soft X-ray emission observed from Green Pea galaxies, indicating that other emission mechanisms (i.e. unusually large quantities of High-/Low-Mass X-ray Binaries, Ultra-Luminous X-ray sources, a modified initial mass function, Intermediate-Mass Black Holes, or AGN) are more likely to be responsible for the X-ray excess. △ Less

Submitted 28 January, 2022; originally announced January 2022.

Comments: Accepted for publication in ApJ. 18 pages, 9 figures

Flash-light on the Ring: hydrodynamic simulations of expandingsupernova shells near supermassive black holes

Authors: B. Barna, J. Palouš, S. Ehlerová, R. Wünsch, M. R. Morris, Pierre Vermot

Abstract: The way supermassive black holes (SMBH) in galactic centers accumulate their mass is not completely determined. At large scales, it is governed by galactic encounters, mass inflows connected to spirals arms and bars, or due to expanding shells from supernova (SN) explosions in the central parts of galaxies. The investigation of the latter process requires an extensive set of gas dynamical simulati… ▽ More The way supermassive black holes (SMBH) in galactic centers accumulate their mass is not completely determined. At large scales, it is governed by galactic encounters, mass inflows connected to spirals arms and bars, or due to expanding shells from supernova (SN) explosions in the central parts of galaxies. The investigation of the latter process requires an extensive set of gas dynamical simulations to explore the muti-dimensional parameter space needed to frame the phenomenon. The aims of this paper are to extend our investigation of the importance of supernovae for inducing accretion onto a SMBH and carry out a comparison between the fully hydrodynamic code Flash and the much less computationally intensive code Ring, which uses the thin shell approximation. We simulate 3D expanding shells in a gravitational potential similar to that of the Galactic Center with a variety of homogeneous and turbulent environments. In homogeneous media, we find convincing agreement between Flash and Ring in the shapes of shells and their equivalent radii throughout their whole evolution until they become subsonic. In highly inhomogeneous, turbulent media, there is also a good agreement of shapes and sizes of shells, and of the times of their first contact with the central 1 pc sphere, where we assume that they join the accretion flow. The comparison supports the proposition that a SN occurring at a galactocentric distance of 5 pc typically drives 1 - 3 $M_\odot$ into the central 1 pc around the galactic center. △ Less

Submitted 22 December, 2021; originally announced December 2021.

Comments: 14 pages, 12 figures, accepted for MNRAS

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

Impact of the external radiation field on the structure and evolution of supernova remnants

Authors: Mario Romero, Yago Ascasibar, Jan Palouš, Richard Wünsch, Mercedes Mollá

Abstract: We carry out 1D hydrodynamical simulations of the evolution of a spherically symmetric supernova remnant (SNR) subject to an external radiation field (ERF) that influences the cooling and heating rates of the gas. We consider homogeneous media with ambient hydrogen number densities $n_{\rm H,0}$ of $0.1$ and $1$ cm$^{-3}$ permeated by an average radiation field including the cosmic microwave, extr… ▽ More We carry out 1D hydrodynamical simulations of the evolution of a spherically symmetric supernova remnant (SNR) subject to an external radiation field (ERF) that influences the cooling and heating rates of the gas. We consider homogeneous media with ambient hydrogen number densities $n_{\rm H,0}$ of $0.1$ and $1$ cm$^{-3}$ permeated by an average radiation field including the cosmic microwave, extragalactic, and Galactic backgrounds, attenuated by an effective column density $N_{\rm H,eff}$ from $10^{18}$ to $10^{21}$~cm$^{-2}$. Our results may be classified into two broad categories: at low $N_{\rm H,eff}$, the ERF presents little absorption in the ultraviolet (ionising) regime, and all the 'unshielded' cases feature an equilibrium temperature $T_{eq} \sim 7000$~K below which the ambient gas cannot cool further. In this scenario, the SNR develops a nearly isothermal shock profile whose shell becomes thicker over time. At higher $N_{\rm H,eff}$, the ERF is heavily absorbed in the UV range, yielding a roughly constant heating function for temperatures $\lesssim 10^4$ K. These `shielded' cases develop a thin, cold and dense shell throughout their evolution. Energy and momentum injection to the medium do not change significantly between both scenarios, albeit luminosity is higher and more uniformly distributed over the shell for unshielded SNR. △ Less

Submitted 8 June, 2021; originally announced June 2021.

Comments: 12 pages, 5 figures, 2 tables. Accepted in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 505, Issue 4, August 2021, pp.5301-5310

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

Extreme adaptive optics astrometry of R136. Searching for high proper motion stars

Authors: Zeinab Khorrami, M. Langlois, F. Vakili, P. C. Clark, A. S. M. Buckner, M. Gonzalez, P. Crowther, R. Wunsch, J. Palous, A. Boccaletti, S. Lumsden, E. Moraux

Abstract: We compared high-contrast near-infrared images of the core of R136 taken by VLT/SPHERE, in two epochs separated by 3.06 years. For the first time we monitored the dynamics of the detected sources in the core of R136 from a ground-based telescope with adaptive optics. The aim of these observations was to search for High prOper Motion cAndidates (HOMAs) in the central region of R136 (r<6") where it… ▽ More We compared high-contrast near-infrared images of the core of R136 taken by VLT/SPHERE, in two epochs separated by 3.06 years. For the first time we monitored the dynamics of the detected sources in the core of R136 from a ground-based telescope with adaptive optics. The aim of these observations was to search for High prOper Motion cAndidates (HOMAs) in the central region of R136 (r<6") where it has been challenging for other instruments. Two bright sources (K<15mag and V<16mag) are located near R136a1 and R136c (massive WR stars) and have been identified as potential HOMAs. These sources have significantly shifted in the images with respect to the mean shift of all reliable detected sources and their neighbours, and six times their own astrometric errors. We calculate their proper motions to be 1.36\pm0.22 mas/yr (321\pm52 km/s) and 1.15\pm0.11 mas/yr (273\pm26 km/s). We discuss different possible scenarios to explain the magnitude of such extreme proper motions, and argue for the necessity to conduct future observations to conclude on the nature of HOMAs in the core of R136. △ Less

Submitted 26 April, 2021; originally announced April 2021.

Comments: Accepted to be published in A&A/Letter. catalogue of the reliable-consistent sources are available online via CDS

Journal ref: A&A 649, L8 (2021)

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

High contrast and resolution near infrared photometry of the core of R136

Authors: Zeinab Khorrami, Maud Langlois, Paul C. Clark, Farrokh Vakili, Anne S. M. Buckner, Marta Gonzalez, Paul Crowther, Richard Wunsch, Jan Palous, Stuart Lumsden, Estelle Moraux

Abstract: We present the sharpest and deepest near infrared photometric analysis of the core of R136, a newly formed massive star cluster at the centre of the 30 Doradus star forming region in the Large Magellanic Cloud. We used the extreme adaptive optics of the SPHERE focal instrument implemented on the ESO Very Large Telescope and operated in its IRDIS imaging mode, for the second time with longer exposu… ▽ More We present the sharpest and deepest near infrared photometric analysis of the core of R136, a newly formed massive star cluster at the centre of the 30 Doradus star forming region in the Large Magellanic Cloud. We used the extreme adaptive optics of the SPHERE focal instrument implemented on the ESO Very Large Telescope and operated in its IRDIS imaging mode, for the second time with longer exposure time in the H- and K filters. Our aim was to (i) increase the number of resolved sources in the core of R136, and (ii) to compare with the first epoch to classify the properties of the detected common sources between the two epochs. Within the field of view (FOV) of 10.8"x12.1" (2.7pc x3.0pc), we detected 1499 sources in both H and K filters, for which 76% of these sources have visual companions closer than 0.2". The larger number of detected sources, enabled us to better sample the mass function (MF). The MF slopes are estimated at ages of 1, 1.5 and 2 Myr, at different radii, and for different mass ranges. The MF slopes for the mass range of 10-300 solar-mass are about 0.3 dex steeper than the mass range of 3-300 solar-mass, for the whole FOV and different radii. Comparing the JHK colours of 790 sources common in between the two epochs, 67% of detected sources in the outer region (r >3") are not consistent with evolutionary models at 1-2 Myr and with extinctions similar to the average cluster value, suggesting an origin from ongoing star formation within 30 Doradus, unrelated to R136. △ Less

Submitted 11 February, 2021; originally announced February 2021.

Comments: 24 pages, 20 figures, 6 tables. Accepted for publication in MNRAS

Can supernova shells feed supermassive black holes in galactic nuclei?

Authors: Jan Palous, Sona Ehlerova, Richrd Wunsch, Mark R. Morris

Abstract: We simulate shells created by supernovae expanding into the interstellar medium (ISM) of the nuclear region of a galaxy, and analyze how the shell evolution is influenced by the supernova (SN) position relative to the galactic center, by the interstellar matter (ISM) density, and by the combined gravitational pull of the nuclear star cluster (NSC) and supermassive black hole (SMBH).We adopted simp… ▽ More We simulate shells created by supernovae expanding into the interstellar medium (ISM) of the nuclear region of a galaxy, and analyze how the shell evolution is influenced by the supernova (SN) position relative to the galactic center, by the interstellar matter (ISM) density, and by the combined gravitational pull of the nuclear star cluster (NSC) and supermassive black hole (SMBH).We adopted simplified hydrodynamical simulations using the infinitesimally thin layer approximation in 3D (code RING) and determined whether and where the shell expansion may bring new gas into the inner parsec around the SMBH. The simulations show that supernovae occurring within a conical region around the rotational axis of the galaxy can feed the central accretion disk surrounding the SMBH. For ambient densities between 10$^3$ and 10$^5$ cm$^{-3}$, the average mass deposited into the central parsec by individual supernovae varies between 10 to 1000 solar masses depending on the ambient density and the spatial distribution of supernova events. Supernova occurring in the aftermath of a starburst event near a galactic center can supply two to three orders of magnitude more mass into the central parsec, depending on the magnitude of the starburst. The deposited mass typically encounters and joins an accretion disk. The fate of that mass is then divided between the growth of the SMBH and an energetically driven outflow from the disk. △ Less

Submitted 29 October, 2020; originally announced October 2020.

Comments: 12 pages, 8 figures, Astronomy and Astrophysiocs, accepted

Journal ref: A&A 644, A72 (2020)

Bonn Optimized Stellar Tracks (BoOST). Simulated Populations of Massive and Very Massive Stars for Astrophysical Applications

Authors: Dorottya Szécsi, Poojan Agrawal, Richard Wünsch, Norbert Langer

Abstract: Massive and very massive stars can play important roles in stellar populations by ejecting strong stellar winds and exploding in energetic phenomena. It is therefore imperative that their behavior be properly accounted for in synthetic model populations. We present nine grids of stellar evolutionary model sequences, together with finely resolved interpolated sequences and synthetic populations,… ▽ More Massive and very massive stars can play important roles in stellar populations by ejecting strong stellar winds and exploding in energetic phenomena. It is therefore imperative that their behavior be properly accounted for in synthetic model populations. We present nine grids of stellar evolutionary model sequences, together with finely resolved interpolated sequences and synthetic populations, of stars with 9-500 Msun and with metallicities ranging from Galactic metallicity down to 1/250 Zsun. The stellar models were computed with the Bonn evolutionary code with consistent physical ingredients, and covering core hydrogen- and core helium-burning phases. The interpolation and population synthesis were performed with our newly developed routine 'synStars'. Eight of the grids represent slowly rotating massive stars with a normal or classical evolutionary path, while one grid represents fast-rotating, chemically homogeneously evolving models. The grids contain data on stellar wind properties such as estimated wind velocity and kinetic energy of the wind, as well as common stellar parameters such as mass, radius, surface temperature, luminosity, mass-loss rate, and surface abundances of 34 isotopes. We also provide estimates of the helium and carbon-oxygen core mass for calculating the mass of stellar remnants. The Bonn Optimized Stellar Tracks (BoOST) project is published as simple tables that include stellar models, interpolated tracks, and synthetic populations. Covering the broadest mass and metallicity range of any published massive star evolutionary model sets to date, BoOST is ideal for further scientific applications such as star formation studies in both low- and high-redshift galaxies. △ Less

Submitted 28 December, 2021; v1 submitted 17 April, 2020; originally announced April 2020.

Comments: Accepted by A&A

Journal ref: A&A 658, A125 (2022)

Faint and fading tails : the fate of stripped HI gas in Virgo cluster galaxies

Authors: Rhys Taylor, Joachim Köppen, Pavel Jáchym, Robert Minchin, Jan Palouš, Richard Wünsch

Abstract: Although many galaxies in the Virgo cluster are known to have lost significant amounts of HI gas, only about a dozen features are known where the HI extends significantly outside its parent galaxy. Previous numerical simulations have predicted that HI removed by ram pressure stripping should have column densities far in excess of the sensitivity limits of observational surveys. We construct a simp… ▽ More Although many galaxies in the Virgo cluster are known to have lost significant amounts of HI gas, only about a dozen features are known where the HI extends significantly outside its parent galaxy. Previous numerical simulations have predicted that HI removed by ram pressure stripping should have column densities far in excess of the sensitivity limits of observational surveys. We construct a simple model to try and quantify how many streams we might expect to detect. This accounts for the expected random orientation of the streams in position and velocity space as well as the expected stream length and mass of stripped HI. Using archival data from the Arecibo Galaxy Environment Survey, we search for any streams which might previously have been missed in earlier analyses. We report the confident detection of ten streams as well as sixteen other less sure detections. We show that these well-match our analytic predictions for which galaxies should be actively losing gas, however the mass of the streams is typically far below the amount of missing HI in their parent galaxies, implying that a phase change and/or dispersal renders the gas undetectable. By estimating the orbital timescales we estimate that dissolution rates of 1-10 Msolar/yr are able to explain both the presence of a few long, massive streams and the greater number of shorter, less massive features. △ Less

Submitted 10 January, 2020; originally announced January 2020.

Comments: Accepted for publication in AJ. 15 pages, 10 figures, 3 tables (plus 10 page appendix with 18 figures)

Supernovae within Pre-existing Wind-Blown Bubbles: Dust Injection vs. Ambient Dust Destruction

Authors: Sergio Martínez-González, Richard Wünsch, Sergiy Silich, Guillermo Tenorio-Tagle, Jan Palouš, Andrea Ferrara

Abstract: By means of 3-D hydrodynamical simulations, here we evaluate the impact that supernova explosions occurring within wind-driven bubbles have on the survival or destruction of dust grains. We consider both, the dust generated within the ejecta and the dust initially present in the ambient gas and later locked-up in the surrounding wind-driven shell. The collision of the supernova blast wave with the… ▽ More By means of 3-D hydrodynamical simulations, here we evaluate the impact that supernova explosions occurring within wind-driven bubbles have on the survival or destruction of dust grains. We consider both, the dust generated within the ejecta and the dust initially present in the ambient gas and later locked-up in the surrounding wind-driven shell. The collision of the supernova blast wave with the wind-driven shell leads to a transmitted shock that moves into the shell and a reflected shock into the ejecta. The transmitted shock is capable of destroying large amounts of the dust locked in the shell, but only if the mass of the wind-driven shell is small, less than a few tens the ejected mass. Conversely, massive wind-driven shells, with several times the ejected mass, lead upon the interaction to strong radiative cooling, which inhibits the Sedov-Taylor phase and weakens the transmitted shock, making it unable to traverse the wind-driven shell. In such a case, the destruction/disruption of the ambient dust is largely inhibited. On the other hand, the SNRs grow rapidly in the very tenuous region excavated by the stellar winds, and thus a large fraction of the dust generated within the ejecta is not efficiently destroyed by the supernova reverse shock, nor by the reflected shock. Our calculations favor a scenario in which core-collapse supernovae within sufficiently massive wind-driven shells supply more dust to the ISM than what they are able to destroy. △ Less

Submitted 12 November, 2019; originally announced November 2019.

Comments: 9 pages, 5 figures. Accepted for publication in The Astrophysical Journal

On the star formation efficiencies and evolution of multiple stellar generations in Globular Clusters

Authors: Guillermo Tenorio-Tagle, Sergiy Silich, Jan Palous, Casiana Muñoz-Tuñón, Richard Wunsch

Abstract: By adopting empirical estimates of the Helium enhancement (Delta Y) between consecutive stellar generations for a sample of Galactic globular clusters (GGC), we uniquely constraint the star formation efficiency of each stellar generation in these stellar systems. In our approach, the star formation efficiency is the central factor that links stellar generations as it defines both their stellar mas… ▽ More By adopting empirical estimates of the Helium enhancement (Delta Y) between consecutive stellar generations for a sample of Galactic globular clusters (GGC), we uniquely constraint the star formation efficiency of each stellar generation in these stellar systems. In our approach, the star formation efficiency is the central factor that links stellar generations as it defines both their stellar mass and the remaining mass available for further star formation, fixing also the amount of matter required to contaminate the next stellar generation. In this way, the star formation efficiency is here shown to be fully defined by the He enhancement between successive stellar generations in a GC. Our approach has also an impact on the evolution of clusters and thus considers the possible loss of stars through evaporation, tidal interactions and stellar evolution. We focus on the present mass ratio between consecutive stellar generations and the present total mass of Galactic globular clusters. Such considerations suffice to determine the relative proportion of stars of consecutive generations that remain today in globular clusters. The latter is also shown to directly depend on the values of Delta Y and thus the He enhancement between consecutive stellar generations in GGC places major constraints on models of star formation and evolution of GC. △ Less

Submitted 7 June, 2019; v1 submitted 23 May, 2019; originally announced May 2019.

Comments: 19 pages, 4 figures, accepted for publication in The Astrophysical Journal

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

Role of supergiants in the formation of globular clusters

Authors: Dorottya Szécsi, Richard Wünsch

Abstract: Multiple stellar populations are observed in almost all globular-clusters, but the origin of this phenomenon is still debated. We investigate the role cool supergiants may have played. To do this, we combine two investigative methods: state-of-the-art massive stellar evolution and calculations of the hydrodynamic structure of the cluster-gas. This approach allows us to study how star-formation in… ▽ More Multiple stellar populations are observed in almost all globular-clusters, but the origin of this phenomenon is still debated. We investigate the role cool supergiants may have played. To do this, we combine two investigative methods: state-of-the-art massive stellar evolution and calculations of the hydrodynamic structure of the cluster-gas. This approach allows us to study how star-formation in young massive clusters depends on the energy- and mass-input of the first-generation of stars, while predicting the chemical composition of the second-generation. We find that the presence of massive (9-500 M$_{\odot}$) metal-poor supergiants in the young cluster leads to a star-formation episode within the first 4 Myr of the cluster's lifetime, that is, before the first core-collapse supernovae explode or the gas is expelled. The stellar winds accumulate in the cluster center, forming the second-generation there. Its composition is predicted to show variations in O & Na abundances, consistently with observations. The abundance of helium is, similarly to other scenarios involving massive stars, higher than what is referred from observations. Supposing dynamical removal of stars from the outskirts of the cluster, or applying a top-heavy initial-mass-function, we can predict a number ratio of the second-generation as high as 20-80%. The effect of metallicity is shown to be important, as the most luminous supergiants are only predicted at low-metallicity, thus limiting $-$ but not excluding $-$ the extent of a polluted second-generation at high-metallicity. These massive stars becoming black-holes suggests globular-clusters hosting gravitational-wave progenitors. Our scenario predicts a correlation between the mass of the cluster and the extent of the multiple population phenomenon. △ Less

Submitted 27 November, 2018; v1 submitted 5 September, 2018; originally announced September 2018.

Comments: Accepted by ApJ, in Press

Pyroclastic Blowout: Dust Survival in Isolated versus Clustered Supernovae

Authors: Sergio Martínez-González, Richard Wünsch, Jan Palouš, Casiana Muñoz-Tuñón, Sergiy Silich, Guillermo Tenorio-Tagle

Abstract: Following the current debate on the fate of SN-condensed dust grains, here we explore by means of three-dimensional hydrodynamical simulations the interaction of dusty supernova remnants (SNRs) with the shocked winds of neighboring massive stars within young massive stellar clusters (SSCs). As a comparison, we have also explored the evolution of supernova remnants in the diffuse ISM with constant… ▽ More Following the current debate on the fate of SN-condensed dust grains, here we explore by means of three-dimensional hydrodynamical simulations the interaction of dusty supernova remnants (SNRs) with the shocked winds of neighboring massive stars within young massive stellar clusters (SSCs). As a comparison, we have also explored the evolution of supernova remnants in the diffuse ISM with constant density. Since the hydrodynamics of SNRs is intimately related to the properties of their immediate environment, the lifecycle of dust grains in SNRs within SSCs is radically different from that in the diffuse ISM. Moreover, off-centered SNRs evolving in the steep density gradient established due to a star cluster wind experience a blowout phase: shell fragmentation due to protruding Rayleigh-Taylor instabilities and the venting of SN ejecta. Our main finding is that clustered SN explosions will cause a net increase in the amount of dust in the surroundings of young massive stellar clusters. Our analysis considers the multiple dust processing resulting from the passage of the SN reverse shock, including its reflection at the SNR's center, the injection of shocked stellar winds within the respective remnant's volume and the effect of secondary forward shocks produced in sequential SN explosions. In the simulations, we have on-the-fly calculated the rates of thermal sputtering and dust-induced radiative cooling provided an initial distribution of grain sizes and dust content. Fast-moving elongated dusty SN ejecta resemble mushroom clouds violently ascending in a stratified atmosphere after volcanic super-eruptions, where the pyroclasts carried by the clouds are wind-driven and eventually accumulate into the vast surroundings. △ Less

Submitted 20 August, 2018; originally announced August 2018.

Comments: 14 pages, 11 figures. Accepted for publication in the Astrophysical Journal

Journal ref: The Astrophysical Journal, Volume 866, Issue 1, article id. 40, 13 pp. (2018)

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

Simulating the evolution of optically dark HI clouds in the Virgo cluster : will no-one rid me of this turbulent sphere ?

Authors: Rhys Taylor, Richard Wünsch, Jan Palouš

Abstract: Most detected neutral atomic hydrogen (HI) at low redshift is associated with optically bright galaxies. However, a handful of HI clouds are known which appear to be optically dark and have no nearby potential progenitor galaxies, making tidal debris an unlikely explanation. In particular, 6 clouds identified by the Arecibo Galaxy Environment Survey are interesting due to the combination of their… ▽ More Most detected neutral atomic hydrogen (HI) at low redshift is associated with optically bright galaxies. However, a handful of HI clouds are known which appear to be optically dark and have no nearby potential progenitor galaxies, making tidal debris an unlikely explanation. In particular, 6 clouds identified by the Arecibo Galaxy Environment Survey are interesting due to the combination of their small size, isolation, and especially their broad line widths atypical of other such clouds. A recent suggestion is that these clouds exist in pressure equilibrium with the intracluster medium, with the line width arising from turbulent internal motions. Here we explore that possibility by using the FLASH code to perform a series of 3D hydro simulations. Our clouds are modelled using spherical Gaussian density profiles, embedded in a hot, low-density gas representing the intracluster medium. The simulations account for heating and cooling of the gas, and we vary the structure and strength of their internal motions. We create synthetic HI spectra, and find that none of our simulations reproduce the observed cloud parameters for longer than about 100 Myr : the clouds either collapse, disperse, or experience rapid heating which would cause ionisation and render them undetectable to HI surveys. While the turbulent motions required to explain the high line widths generate structures which appear to be inherently unstable, making this an unlikely explanation for the observed clouds, these simulations demonstrate the importance of including the intracluster medium in any model seeking to explain the existence of these objects. △ Less

Submitted 9 May, 2018; originally announced May 2018.

Comments: 15 pages, 11 figures. Accepted for publication in MNRAS

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.

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

Tree-based solvers for adaptive mesh refinement code FLASH - I: gravity and optical depths

Authors: Richard Wünsch, Stefanie Walch, František Dinnbier, Anthony P. Whitworth

Abstract: We describe an OctTree algorithm for the MPI-parallel, adaptive mesh-refinement code {\sc FLASH}, which can be used to calculate the gas self-gravity, and also the angle-averaged local optical depth, for treating ambient diffuse radiation. The algorithm communicates to the different processors only those parts of the tree that are needed to perform the tree walk locally. The advantage of this appr… ▽ More We describe an OctTree algorithm for the MPI-parallel, adaptive mesh-refinement code {\sc FLASH}, which can be used to calculate the gas self-gravity, and also the angle-averaged local optical depth, for treating ambient diffuse radiation. The algorithm communicates to the different processors only those parts of the tree that are needed to perform the tree walk locally. The advantage of this approach is a relatively low memory requirement, important in particular for the optical depth calculation, which needs to process information from many different directions. This feature also enables a general tree-based radiation transport algorithm that will be described in a subsequent paper, and delivers excellent scaling up to at least 1500 cores. Boundary conditions for gravity can be either isolated or periodic, and they can be specified in each direction independently, using a newly developed generalisation of the Ewald method. The gravity calculation can be accelerated with the {\em adaptive block update} technique by partially re-using the solution from the previous time-step. Comparison with the {\sc Flash} internal multi-grid gravity solver shows that tree based methods provide a competitive alternative, particularly for problems with isolated or mixed boundary conditions. We evaluate several multipole acceptance criteria (MACs) and identify a relatively simple APE MAC which provides high accuracy at low computational cost. The optical depth estimates are found to agree very well with those of the {\sc RADMC-3D} radiation transport code, with the tree solver being much faster. Our algorithm is available in the standard release of the {\sc FLASH} code in version 4.0 and later. △ Less

Submitted 2 January, 2018; v1 submitted 21 August, 2017; originally announced August 2017.

Comments: 26 pages, 19 figures, 6 tables. Accepted by MNRAS

Globular Cluster formation in a collapsing supershell

Authors: S. Recchi, R. Wünsch, J. Palous, F. Dinnbier

Abstract: Primordial clouds are supposed to host the so-called population III stars. These stars are very massive and completely metal-free. The final stage of the life of population III stars with masses between 130 and 260 solar masses is a very energetic hypernova explosion. A hypernova drives a shock, behind which a spherically symmetric very dense supershell forms, which might become gravitationally un… ▽ More Primordial clouds are supposed to host the so-called population III stars. These stars are very massive and completely metal-free. The final stage of the life of population III stars with masses between 130 and 260 solar masses is a very energetic hypernova explosion. A hypernova drives a shock, behind which a spherically symmetric very dense supershell forms, which might become gravitationally unstable, fragment, and form stars. In this paper we study under what conditions can an expanding supershell become gravitationally unstable and how the feedback of these supershell stars (SSSs) affects its surroundings. We simulate, by means of a 1-D Eulerian hydrocode, the early evolution of the primordial cloud after the hypernova explosion, the formation of SSSs, and the following evolution, once the SSSs start to release energy and heavy elements into the interstellar medium. Our results indicate that a shell, enriched with nucleosynthetic products from SSSs, propagates inwards, towards the center of the primordial cloud. In a time span of a few Myr, this inward-propagating shell reaches a distance of only a few parsec away from the center of the primordial cloud. Its density is extremely high and its temperature very low, thus the conditions for a new episode of star formation are achieved. We study what fraction of these two distinct populations of stars can remain bound and survive until the present day. We study also under what conditions can this process repeat and form multiple stellar populations. We extensively discuss whether the proposed scenario can help to explain some open questions of the formation mechanism of globular clusters. △ Less

Submitted 10 September, 2017; v1 submitted 16 August, 2017; originally announced August 2017.

Comments: 15 pages, 7 figures, published on Ap&SS. Replaced after editorial corrections

Can Dust Injected by SNe Explain the NIR-MIR Excess in Young Massive Stellar Clusters?

Authors: Sergio Martinez-Gonzalez, Richard Wünsch, Jan Palouš

Abstract: We present a physically-motivated model involving the different processes affecting supernova dust grains as they are incorporated into the thermalized medium within young massive star clusters. The model is used to explain the near- to mid-infrared (NIR-MIR) excess found in such clusters and usually modeled as a blackbody with temperature $\sim (400-1000)$ K. In our approach, dust grains are effi… ▽ More We present a physically-motivated model involving the different processes affecting supernova dust grains as they are incorporated into the thermalized medium within young massive star clusters. The model is used to explain the near- to mid-infrared (NIR-MIR) excess found in such clusters and usually modeled as a blackbody with temperature $\sim (400-1000)$ K. In our approach, dust grains are efficiently produced in the clumpy ejecta of core-collapse supernovae, fragmented into small pieces ($\lesssim 0.05$ $μみゅー$m) as they are incorporated into the hot and dense ISM, heated via frequent collisions with electrons and the absorption of energetic photons. Grains with small sizes can more easily acquire the high temperatures ($\sim 1000$ K) required to produce a NIR-MIR excess with respect to the emission of foreground PAHs and starlight. However, the extreme conditions inside young massive clusters make difficult for these small grains to have a persistent manifestation at NIR-MIR wavelengths as they are destroyed by efficient thermal sputtering. Nevertheless, the chances for a persistent manifestation are increased by taking into account that small grains become increasingly transparent to their impinging ions as their size decreases. For an individual SN event, we find that the NIR-MIR excess last longer if the time required to incorporate all the grains into the thermalized medium is also longer, and in some cases, comparable to the characteristic interval between supernova explosions. Our models, can successfully explain the near-infrared excesses found in the star clusters observed in M33 (Relaño et al. 2016) assuming a low heating efficiency and mass-loading. In this scenario, the presence of the NIR-MIR excess is an indication of efficient dust production in SNe and its subsequent destruction △ Less

Submitted 24 May, 2017; originally announced May 2017.

Comments: 13 pages, 8 figures, Accepted for publication in The Astrophysical Journal

Journal ref: The Astrophysical Journal, Volume 843, Issue 2, article id. 95, 12 pp. (2017)

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

Kinematic clues to the origins of starless HI clouds : dark galaxies or tidal debris ?

Authors: R. Taylor, J. I. Davies, P. Jáchym, O. Keenan, R. F. Minchin, J. Palouš, R. Smith, R. Wünsch

Abstract: Isolated HI clouds with no optical counterparts are often taken as evidence for galaxy-galaxy interactions, though an alternative hypothesis is that these are primordial 'dark galaxies' which have not formed stars. Similarly, certain kinematic features in HI streams are also controversial, sometimes taken as evidence of dark galaxies but also perhaps explicable as the result of harassment. We nume… ▽ More Isolated HI clouds with no optical counterparts are often taken as evidence for galaxy-galaxy interactions, though an alternative hypothesis is that these are primordial 'dark galaxies' which have not formed stars. Similarly, certain kinematic features in HI streams are also controversial, sometimes taken as evidence of dark galaxies but also perhaps explicable as the result of harassment. We numerically model the passage of a galaxy through the gravitational field of cluster. The galaxy consists of SPH particles for the gas and n-bodies for the stars and dark matter, while the cluster includes the gravitational effects of substructure using 400 subhalos (the effects of the intracluster medium are ignored). We find that harassment can indeed produce long HI streams and these streams can include kinematic features resembling dark galaxy candidates such as VIRGOHI21. We also show that apparent clouds with diameter < 20 kpc and velocity widths < 50 km/s are almost invariably produced in these simulations, making tidal debris a highly probable explanation. In contrast, we show that the frequency of isolated clouds of the same size but velocity width > 100 km/s is negligible - making this a very unlikely explanation for the observed clouds in the Virgo cluster with these properties. △ Less

Submitted 19 January, 2017; originally announced January 2017.

Comments: 15 pages, 12 figures (26 pages, 25 figures including appendix). Accepted for publication in MNRAS

Fragmentation of vertically stratified gaseous layers: monolithic or coalescence-driven collapse

Authors: František Dinnbier, Richard Wünsch, Anthony P. Whitworth, Jan Palouš

Abstract: We investigate, using 3D hydrodynamic simulations, the fragmentation of pressure-confined, vertically stratified, self-gravitating gaseous layers. The confining pressure is either thermal pressure acting on both surfaces, or thermal pressure acting on one surface and ram-pressure on the other. In the linear regime of fragmentation, the dispersion relation we obtain agrees well with that derived by… ▽ More We investigate, using 3D hydrodynamic simulations, the fragmentation of pressure-confined, vertically stratified, self-gravitating gaseous layers. The confining pressure is either thermal pressure acting on both surfaces, or thermal pressure acting on one surface and ram-pressure on the other. In the linear regime of fragmentation, the dispersion relation we obtain agrees well with that derived by Elmegreen & Elmegreen (1978), and consequently deviates from the dispersion relations based on the thin shell approximation (Vishniac 1983) or pressure assisted gravitational instability (Wünsch et al. 2010). In the non-linear regime, the relative importance of the confining pressure to the self-gravity is a crucial parameter controlling the qualitative course of fragmentation. When confinement of the layer is dominated by external pressure, self- gravitating condensations are delivered by a two-stage process: first the layer fragments into gravitationally bound but stable clumps, and then these clumps coalesce until they assemble enough mass to collapse. In contrast, when external pressure makes a small contribution to confinement of the layer, the layer fragments monolithically into gravitationally unstable clumps and there is no coalescence. This dichotomy persists whether the external pressure is thermal or ram. We apply these results to fragments forming in a shell swept up by an expanding H II region, and find that, unless the swept up gas is quite hot or the surrounding medium has low density, the fragments have low-mass ( ~< 3 M_Sun ), and therefore they are unlikely to spawn stars that are sufficiently massive to promote sequential self-propagating star formation. △ Less

Submitted 10 January, 2017; originally announced January 2017.

Comments: 21 pages, 12 figures, MNRAS in press

The formation of secondary stellar generations in massive young star clusters from rapidly cooling shocked stellar winds

Authors: Richard Wünsch, Jan Palouš, Guillermo Tenorio-Tagle, Soňa Ehlerová

Abstract: We study a model of rapidly cooling shocked stellar winds in young massive clusters and estimate the circumstances under which secondary star formation, out of the reinserted winds from a first stellar generation (1G), is possible. We have used two implementations of the model: a highly idealized computationally inexpensive spherically symmetric semi-analytic model, and a complex three-dimensional… ▽ More We study a model of rapidly cooling shocked stellar winds in young massive clusters and estimate the circumstances under which secondary star formation, out of the reinserted winds from a first stellar generation (1G), is possible. We have used two implementations of the model: a highly idealized computationally inexpensive spherically symmetric semi-analytic model, and a complex three-dimensional radiation-hydrodynamic simulations, and they are in a good mutual agreement. The results confirm our previous findings that in a cluster with 1G mass $10^7$ M$_\odot$ and half-mass radius $2.38$ pc, the shocked stellar winds become thermally unstable, collapse into dense gaseous structures that partially accumulate inside the cluster, self-shield against ionizing stellar radiation and form the second generation (2G) of stars. We have used the semi-analytic model to explore a subset of the parameter space covering a wide range of the observationally poorly constrained parameters: the heating efficiency, $ηいーた_\mathrm{he}$, and the mass loading, $ηいーた_\mathrm{ml}$. The results show that the fraction of the 1G stellar winds accumulating inside the cluster can be larger than $50$ % if $ηいーた_\mathrm{he} \lesssim 10$ % which is suggested by the observations. Furthermore, for low $ηいーた_\mathrm{he}$, the model provides a self-consistent mechanism predicting 2G stars forming only in the central zones of the cluster. Finally, we have calculated the accumulated warm gas emission in the H30$αあるふぁ$ recombination line, analyzed its velocity profile and estimated its intensity for super star clusters in interacting galaxies NGC4038/9 (Antennae) showing that the warm gas should be detectable with ALMA. △ Less

Submitted 15 December, 2016; originally announced December 2016.

Comments: 18 pages, 10 figures, accepted by ApJ

The impact of magnetic fields on the chemical evolution of the supernova-driven ISM

Authors: Anabele-Linda Pardi, Philipp Girichidis, Thorsten Naab, Stefanie Walch, Thomas Peters, Fabian Heitsch, Simon C. O. Glover, Ralf S. Klessen, Richard Wünsch, Andrea Gatto

Abstract: We present three-dimensional magneto-hydrodynamical simulations of the self-gravitating interstellar medium (ISM) in a periodic (256 pc)$^3$ box with a mean number density of 0.5 cm$^{-3}$. At a fixed supernova rate we investigate the multi-phase ISM structure, H$_{2}$ molecule formation and density-magnetic field scaling for varying initial magnetic field strengths (0, $6\times 10^{-3}$, 0.3, 3… ▽ More We present three-dimensional magneto-hydrodynamical simulations of the self-gravitating interstellar medium (ISM) in a periodic (256 pc)$^3$ box with a mean number density of 0.5 cm$^{-3}$. At a fixed supernova rate we investigate the multi-phase ISM structure, H$_{2}$ molecule formation and density-magnetic field scaling for varying initial magnetic field strengths (0, $6\times 10^{-3}$, 0.3, 3 $μみゅー$G). All magnetic runs saturate at mass weighted field strengths of $\sim$ 1 $-$ 3 $μみゅー$G but the ISM structure is notably different. With increasing initial field strengths (from $6\times 10^{-3}$ to 3 $μみゅー$G) the simulations develop an ISM with a more homogeneous density and temperature structure, with increasing mass (from 5% to 85%) and volume filling fractions (from 4% to 85%) of warm (300 K $<$ T $<$ 8000 K) gas, with decreasing volume filling fractions (VFF) from $\sim$ 35% to $\sim$ 12% of hot gas (T $> 10^5$ K) and with a decreasing H$_{2}$ mass fraction (from 70% to $<$ 1%). Meanwhile the mass fraction of gas in which the magnetic pressure dominates over the thermal pressure increases by a factor of 10, from 0.07 for an initial field of $6\times 10^{-3}$ $μみゅー$G to 0.7 for a 3 $μみゅー$G initial field. In all but the simulations with the highest initial field strength self-gravity promotes the formation of dense gas and H$_{2}$, but does not change any other trends. We conclude that magnetic fields have a significant impact on the multi-phase, chemical and thermal structure of the ISM and discuss potential implications and limitations of the model. △ Less

Submitted 2 November, 2016; originally announced November 2016.

Comments: submitted

The SILCC project --- IV. Impact of dissociating and ionising radiation on the interstellar medium and Halpha emission as a tracer of the star formation rate

Authors: Thomas Peters, Thorsten Naab, Stefanie Walch, Simon C. O. Glover, Philipp Girichidis, Eric Pellegrini, Ralf S. Klessen, Richard Wünsch, Andrea Gatto, Christian Baczynski

Abstract: We present three-dimensional radiation-hydrodynamical simulations of the impact of stellar winds, photoelectric heating, photodissociating and photoionising radiation, and supernovae on the chemical composition and star formation in a stratified disc model. This is followed with a sink-based model for star clusters with populations of individual massive stars. Stellar winds and ionising radiation… ▽ More We present three-dimensional radiation-hydrodynamical simulations of the impact of stellar winds, photoelectric heating, photodissociating and photoionising radiation, and supernovae on the chemical composition and star formation in a stratified disc model. This is followed with a sink-based model for star clusters with populations of individual massive stars. Stellar winds and ionising radiation regulate the star formation rate at a factor of ~10 below the simulation with only supernova feedback due to their immediate impact on the ambient interstellar medium after star formation. Ionising radiation (with winds and supernovae) significantly reduces the ambient densities for most supernova explosions to rho < 10^-25 g cm^-3, compared to 10^-23 g cm^-3 for the model with only winds and supernovae. Radiation from massive stars reduces the amount of molecular hydrogen and increases the neutral hydrogen mass and volume filling fraction. Only this model results in a molecular gas depletion time scale of 2 Gyr and shows the best agreement with observations. In the radiative models, the Halpha emission is dominated by radiative recombination as opposed to collisional excitation (the dominant emission in non-radiative models), which only contributes ~1-10 % to the total Halpha emission. Individual massive stars (M >= 30 M_sun) with short lifetimes are responsible for significant fluctuations in the Halpha luminosities. The corresponding inferred star formation rates can underestimate the true instantaneous star formation rate by factors of ~10. △ 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 466 (2017), 3293-3308

Attack of the Flying Snakes : Formation of Isolated HI Clouds By Fragmentation of Long Streams

Authors: R. Taylor, J. I. Davies, P. Jachym, O. Keenan, R. F. Minchin, J. Palous, R. Smith, R. Wunsch

Abstract: The existence of long (> 100 kpc) HI streams and small (< 20 kpc) free-floating HI clouds is well-known. While the formation of the streams has been investigated extensively, and the isolated clouds are often purported to be interaction debris, little research has been done on the formation of optically dark HI clouds that are not part of a larger stream. One possibility is that such features resu… ▽ More The existence of long (> 100 kpc) HI streams and small (< 20 kpc) free-floating HI clouds is well-known. While the formation of the streams has been investigated extensively, and the isolated clouds are often purported to be interaction debris, little research has been done on the formation of optically dark HI clouds that are not part of a larger stream. One possibility is that such features result from the fragmentation of more extended streams, while another idea is that they are primordial, optically dark galaxies. We test the validity of the fragmentation scenario (via harassment) using numerical simulations. In order to compare our numerical models with observations, we present catalogues of both the known long HI streams (42 objects) and free-floating HI clouds suggested as dark galaxy candidates (51 objects). In particular, we investigate whether it is possible to form compact features with high velocity widths (> 100 km/s), similar to observed clouds which are otherwise intriguing dark galaxy candidates. We find that producing such features is possible but extremely unlikely, occurring no more than 0.2% of the time in our simulations. In contrast, we find that genuine dark galaxies could be extremely stable to harassment and remain detectable even after 5 Gyr in the cluster environment (with the important caveat that our simulations only explore harassment and do not yet include the intracluster medium, heating and cooling, or star formation). We also discuss the possibility that such objects could be the progenitors of recently discovered ultra diffuse galaxies. △ Less

Submitted 17 June, 2016; originally announced June 2016.

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

The SILCC project: III. Regulation of star formation and outflows by stellar winds and supernovae

Authors: A. Gatto, S. Walch, T. Naab, P. Girichidis, R. Wünsch, S. C. O. Glover, R. S. Klessen, P. C. Clark, T. Peters, D. Derigs, C. Baczynski, J. Puls

Abstract: We study the impact of stellar winds and supernovae on the multi-phase interstellar medium using three-dimensional hydrodynamical simulations carried out with FLASH. The selected galactic disc region has a size of (500 pc)$^2$ x $\pm$ 5 kpc and a gas surface density of 10 M$_{\odot}$/pc$^2$. The simulations include an external stellar potential and gas self-gravity, radiative cooling and diffuse h… ▽ More We study the impact of stellar winds and supernovae on the multi-phase interstellar medium using three-dimensional hydrodynamical simulations carried out with FLASH. The selected galactic disc region has a size of (500 pc)$^2$ x $\pm$ 5 kpc and a gas surface density of 10 M$_{\odot}$/pc$^2$. The simulations include an external stellar potential and gas self-gravity, radiative cooling and diffuse heating, sink particles representing star clusters, stellar winds from these clusters which combine the winds from indi- vidual massive stars by following their evolution tracks, and subsequent supernova explosions. Dust and gas (self-)shielding is followed to compute the chemical state of the gas with a chemical network. We find that stellar winds can regulate star (cluster) formation. Since the winds suppress the accretion of fresh gas soon after the cluster has formed, they lead to clusters which have lower average masses (10$^2$ - 10$^{4.3}$ M$_{\odot}$) and form on shorter timescales (10$^{-3}$ - 10 Myr). In particular we find an anti-correlation of cluster mass and accretion time scale. Without winds the star clusters easily grow to larger masses for ~5 Myr until the first supernova explodes. Overall the most massive stars provide the most wind energy input, while objects beginning their evolution as B-type stars contribute most of the supernova energy input. A significant outflow from the disk (mass loading $\gtrsim$ 1 at 1 kpc) can be launched by thermal gas pressure if more than 50% of the volume near the disc mid-plane can be heated to T > 3x10$^5$ K. Stellar winds alone cannot create a hot volume-filling phase. The models which are in best agreement with observed star formation rates drive either no outflows or weak outflows. △ Less

Submitted 16 June, 2016; originally announced June 2016.

Comments: 23 pages; submitted to MNRAS

Bimodal regime in young massive clusters leading to subsequent stellar generations

Authors: Richard Wünsch, Jan Palouš, Guillermo Tenorio-Tagle, Casiana Muñoz-Tuñón, Soňa Ehlerová

Abstract: Massive stars in young massive clusters insert tremendous amounts of mass and energy into their surroundings in the form of stellar winds and supernova ejecta. Mutual shock-shock collisions lead to formation of hot gas, filling the volume of the cluster. The pressure of this gas then drives a powerful cluster wind. However, it has been shown that if the cluster is massive and dense enough, it can… ▽ More Massive stars in young massive clusters insert tremendous amounts of mass and energy into their surroundings in the form of stellar winds and supernova ejecta. Mutual shock-shock collisions lead to formation of hot gas, filling the volume of the cluster. The pressure of this gas then drives a powerful cluster wind. However, it has been shown that if the cluster is massive and dense enough, it can evolve in the so--called bimodal regime, in which the hot gas inside the cluster becomes thermally unstable and forms dense clumps which are trapped inside the cluster by its gravity. We will review works on the bimodal regime and discuss the implications for the formation of subsequent stellar generations. The mass accumulates inside the cluster and as soon as a high enough column density is reached, the interior of the clumps becomes self-shielded against the ionising radiation of stars and the clumps collapse and form new stars. The second stellar generation will be enriched by products of stellar evolution from the first generation, and will be concentrated near the cluster center. △ Less

Submitted 15 January, 2016; originally announced January 2016.

Comments: 8 pages, 4 figures, invited talk at IAU Symposium 316 "Formation, evolution, and survival of massive star clusters", C. Charbonnel & A. Nota, eds

Self-shielding clumps in starburst clusters

Authors: Jan Palouš, Richard Wünsch, Soňa Ehlerová, Guillermo Tenorio-Tagle

Abstract: Young and massive star clusters above a critical mass form thermally unstable clumps reducing locally the temperature and pressure of the hot 10$^{7}$~K cluster wind. The matter reinserted by stars, and mass loaded in interactions with pristine gas and from evaporating circumstellar disks, accumulate on clumps that are ionized with photons produced by massive stars. We discuss if they may become s… ▽ More Young and massive star clusters above a critical mass form thermally unstable clumps reducing locally the temperature and pressure of the hot 10$^{7}$~K cluster wind. The matter reinserted by stars, and mass loaded in interactions with pristine gas and from evaporating circumstellar disks, accumulate on clumps that are ionized with photons produced by massive stars. We discuss if they may become self-shielded when they reach the central part of the cluster, or even before it, during their free fall to the cluster center. Here we explore the importance of heating efficiency of stellar winds. △ Less

Submitted 15 January, 2016; originally announced January 2016.

Comments: 2 pages, 4 figures, IAU Symposium 316: Formation, Evolution, and Survival of Massive Star Clusters

Impact of supernova and cosmic-ray driving on the surface brightness of the galactic halo in soft X-rays

Authors: Thomas Peters, Philipp Girichidis, Andrea Gatto, Thorsten Naab, Stefanie Walch, Richard Wünsch, Simon C. O. Glover, Paul C. Clark, Ralf S. Klessen, Christian Baczynski

Abstract: The halo of the Milky Way contains a hot plasma with a surface brightness in soft X-rays of the order $10^{-12}$erg cm$^{-2}$ s$^{-1}$ deg$^{-2}$. The origin of this gas is unclear, but so far numerical models of galactic star formation have failed to reproduce such a large surface brightness by several orders of magnitude. In this paper, we analyze simulations of the turbulent, magnetized, multi-… ▽ More The halo of the Milky Way contains a hot plasma with a surface brightness in soft X-rays of the order $10^{-12}$erg cm$^{-2}$ s$^{-1}$ deg$^{-2}$. The origin of this gas is unclear, but so far numerical models of galactic star formation have failed to reproduce such a large surface brightness by several orders of magnitude. In this paper, we analyze simulations of the turbulent, magnetized, multi-phase interstellar medium including thermal feedback by supernova explosions as well as cosmic-ray feedback. We include a time-dependent chemical network, self-shielding by gas and dust, and self-gravity. Pure thermal feedback alone is sufficient to produce the observed surface brightness, although it is very sensitive to the supernova rate. Cosmic rays suppress this sensitivity and reduce the surface brightness because they drive cooler outflows. Self-gravity has by far the largest effect because it accumulates the diffuse gas in the disk in dense clumps and filaments, so that supernovae exploding in voids can eject a large amount of hot gas into the halo. This can boost the surface brightness by several orders of magnitude. Although our simulations do not reach a steady state, all simulations produce surface brightness values of the same order of magnitude as the observations, with the exact value depending sensitively on the simulation parameters. We conclude that star formation feedback alone is sufficient to explain the origin of the hot halo gas, but measurements of the surface brightness alone do not provide useful diagnostics for the study of galactic star formation. △ Less

Submitted 22 October, 2015; originally announced October 2015.

Comments: ApJ Letters in press

Journal ref: Astrophysical Journal Letters 813 (2015) L27

Launching cosmic-ray-driven outflows from the magnetized interstellar medium

Authors: Philipp Girichidis, Thorsten Naab, Stefanie Walch, Michal Hanasz, Mordecai-Mark Mac Low, Jeremiah P. Ostriker, Andrea Gatto, Thomas Peters, Richard Wünsch, Simon C. O. Glover, Ralf S. Klessen, Paul C. Clark, Christian Baczynski

Abstract: We present a hydrodynamical simulation of the turbulent, magnetized, supernova (SN)-driven interstellar medium (ISM) in a stratified box that dynamically couples the injection and evolution of cosmic rays (CRs) and a self-consistent evolution of the chemical composition. CRs are treated as a relativistic fluid in the advection-diffusion approximation. The thermodynamic evolution of the gas is comp… ▽ More We present a hydrodynamical simulation of the turbulent, magnetized, supernova (SN)-driven interstellar medium (ISM) in a stratified box that dynamically couples the injection and evolution of cosmic rays (CRs) and a self-consistent evolution of the chemical composition. CRs are treated as a relativistic fluid in the advection-diffusion approximation. The thermodynamic evolution of the gas is computed using a chemical network that follows the abundances of H+, H, H2, CO, C+, and free electrons and includes (self-)shielding of the gas and dust. We find that CRs perceptibly thicken the disk with the heights of 90% (70%) enclosed mass reaching ~1.5 kpc (~0.2 kpc). The simulations indicate that CRs alone can launch and sustain strong outflows of atomic and ionized gas with mass loading factors of order unity, even in solar neighborhood conditions and with a CR energy injection per SN of 10^50 erg, 10% of the fiducial thermal energy of an SN. The CR-driven outflows have moderate launching velocities close to the midplane (~100 km/s) and are denser (ρろー~1e-24 - 1e-26 g/cm^3), smoother, and colder than the (thermal) SN-driven winds. The simulations support the importance of CRs for setting the vertical structure of the disk as well as the driving of winds. △ Less

Submitted 6 January, 2016; v1 submitted 24 September, 2015; originally announced September 2015.

Comments: 7 pages, 5 figures, ApJL, accepted

The SILCC (SImulating the LifeCycle of molecular Clouds) project - II. Dynamical evolution of the supernova-driven ISM and the launching of outflows

Authors: Philipp Girichidis, Stefanie Walch, Thorsten Naab, Andrea Gatto, Richard Wünsch, Simon C. O. Glover, Ralf S. Klessen, Paul C. Clark, Thomas Peters, Dominik Derigs, Christian Baczynski

Abstract: The SILCC project (SImulating the Life-Cycle of molecular Clouds) aims at a more self-consistent understanding of the interstellar medium (ISM) on small scales and its link to galaxy evolution. We present three-dimensional (magneto)hydrodynamic simulations of the ISM in a vertically stratified box including self-gravity, an external potential due to the stellar component of the galactic disc, and… ▽ More The SILCC project (SImulating the Life-Cycle of molecular Clouds) aims at a more self-consistent understanding of the interstellar medium (ISM) on small scales and its link to galaxy evolution. We present three-dimensional (magneto)hydrodynamic simulations of the ISM in a vertically stratified box including self-gravity, an external potential due to the stellar component of the galactic disc, and stellar feedback in the form of an interstellar radiation field and supernovae (SNe). The cooling of the gas is based on a chemical network that follows the abundances of H+, H, H2, C+, and CO and takes shielding into account consistently. We vary the SN feedback by comparing different SN rates, clustering and different positioning, in particular SNe in density peaks and at random positions, which has a major impact on the dynamics. Only for random SN positions the energy is injected in sufficiently low-density environments to reduce energy losses and enhance the effective kinetic coupling of the SNe with the gas. This leads to more realistic velocity dispersions (σしぐま_HI ~ 0.8σしぐま_(300-8000K) ~ 10-20km/s, σしぐま_Hαあるふぁ~ 0.6σしぐま_(8000-3e5K) ~ 20-30km/s), and strong outflows with mass loading factors of up to 10 even for solar neighbourhood conditions. Clustered SNe abet the onset of outflows compared to individual SNe but do not influence the net outflow rate. The outflows do not contain any molecular gas and are mainly composed of atomic hydrogen. The bulk of the outflowing mass is dense (ρろー~ 1e-25-1e-24g/cc) and slow (v ~ 20-40km/s) but there is a high-velocity tail of up to v ~ 500km/s with ρろー~ 1e-28-1e-27g/cc. △ Less

Submitted 11 December, 2015; v1 submitted 26 August, 2015; originally announced August 2015.

Comments: 26 pages, MNRAS, accepted

StarBench: The D-type expansion of an HII region

Authors: T. G. Bisbas, T. J. Haworth, R. J. R. Williams, J. Mackey, P. Tremblin, A. C. Raga, S. J. Arthur, C. Baczynski, J. E. Dale, T. Frostholm, S. Geen, T. Haugboelle, D. Hubber, I. T. Iliev, R. Kuiper, J. Rosdahl, D. Sullivan, S. Walch, R. Wuensch

Abstract: StarBench is a project focused on benchmarking and validating different star-formation and stellar feedback codes. In this first StarBench paper we perform a comparison study of the D-type expansion of an HII region. The aim of this work is to understand the differences observed between the twelve participating numerical codes against the various analytical expressions examining the D-type phase o… ▽ More StarBench is a project focused on benchmarking and validating different star-formation and stellar feedback codes. In this first StarBench paper we perform a comparison study of the D-type expansion of an HII region. The aim of this work is to understand the differences observed between the twelve participating numerical codes against the various analytical expressions examining the D-type phase of HII region expansion. To do this, we propose two well-defined tests which are tackled by 1D and 3D grid- and SPH- based codes. The first test examines the `early phase' D-type scenario during which the mechanical pressure driving the expansion is significantly larger than the thermal pressure of the neutral medium. The second test examines the `late phase' D-type scenario during which the system relaxes to pressure equilibrium with the external medium. Although they are mutually in excellent agreement, all twelve participating codes follow a modified expansion law that deviates significantly from the classical Spitzer solution in both scenarios. We present a semi-empirical formula combining the two different solutions appropriate to both early and late phases that agrees with high-resolution simulations to $\lesssim2\%$. This formula provides a much better benchmark solution for code validation than the Spitzer solution. The present comparison has validated the participating codes and through this project we provide a dataset for calibrating the treatment of ionizing radiation hydrodynamics codes. △ Less

Submitted 20 July, 2015; originally announced July 2015.

Comments: 20 pages, 12 Figures, 4 Tables. Accepted for publication in MNRAS. Comments are welcome. Participation in future StarBench tests is also welcome

The SILCC (SImulating the LifeCycle of molecular Clouds) project: I. Chemical evolution of the supernova-driven ISM

Authors: S. K. Walch, P. Girichidis, T. Naab, A. Gatto, S. C. O. Glover, R. Wünsch, R. S Klessen, P. C. Clark, T. Peters, C. Baczynski

Abstract: The SILCC project (SImulating the Life-Cycle of molecular Clouds) aims at a more self-consistent understanding of the interstellar medium (ISM) on small scales and its link to galaxy evolution. We simulate the evolution of the multi-phase ISM in a 500 pc x 500 pc x 10 kpc region of a galactic disc, with a gas surface density of $Σしぐま_{_{\rm GAS}} = 10 \;{\rm M}_\odot/{\rm pc}^2$. The Flash 4.1 simula… ▽ More The SILCC project (SImulating the Life-Cycle of molecular Clouds) aims at a more self-consistent understanding of the interstellar medium (ISM) on small scales and its link to galaxy evolution. We simulate the evolution of the multi-phase ISM in a 500 pc x 500 pc x 10 kpc region of a galactic disc, with a gas surface density of $Σしぐま_{_{\rm GAS}} = 10 \;{\rm M}_\odot/{\rm pc}^2$. The Flash 4.1 simulations include an external potential, self-gravity, magnetic fields, heating and radiative cooling, time-dependent chemistry of H$_2$ and CO considering (self-) shielding, and supernova (SN) feedback. We explore SN explosions at different (fixed) rates in high-density regions (peak), in random locations (random), in a combination of both (mixed), or clustered in space and time (clustered). Only random or clustered models with self-gravity (which evolve similarly) are in agreement with observations. Molecular hydrogen forms in dense filaments and clumps and contributes 20% - 40% to the total mass, whereas most of the mass (55% - 75%) is in atomic hydrogen. The ionised gas contributes <10%. For high SN rates (0.5 dex above Kennicutt-Schmidt) as well as for peak and mixed driving the formation of H$_2$ is strongly suppressed. Also without self-gravity the H$_2$ fraction is significantly lower ($\sim$ 5%). Most of the volume is filled with hot gas ($\sim$90% within $\pm$2 kpc). Only for random or clustered driving, a vertically expanding warm component of atomic hydrogen indicates a fountain flow. Magnetic fields have little impact on the final disc structure. However, they affect dense gas ($n\gtrsim 10\;{\rm cm}^{-3}$) and delay H$_2$ formation. We highlight that individual chemical species, in particular atomic hydrogen, populate different ISM phases and cannot be accurately accounted for by simple temperature-/density-based phase cut-offs. △ Less

Submitted 8 December, 2014; originally announced December 2014.

Comments: 30 pages, 23 figures, submitted to MNRAS. Comments welcome! For movies of the simulations and download of selected Flash data see the SILCC website: http://www.astro.uni-koeln.de/silcc

Modelling the supernova-driven ISM in different environments

Authors: A. Gatto, S. Walch, M. -M. Mac Low, T. Naab, P. Girichidis, S. C. O. Glover, R. Wünsch, R. S. Klessen, P. C. Clark, C. Baczynski, T. Peters, J. P. Ostriker, J. C. Ibáñez-Mejía, S. Haid

Abstract: We use hydrodynamical simulations in a $(256\;{\rm pc})^3$ periodic box to model the impact of supernova (SN) explosions on the multi-phase interstellar medium (ISM) for initial densities $n = 0.5-30$ cm$^{-3}$ and SN rates $1-720$ Myr$^{-1}$. We include radiative cooling, diffuse heating, and the formation of molecular gas using a chemical network. The SNe explode either at random positions, at d… ▽ More We use hydrodynamical simulations in a $(256\;{\rm pc})^3$ periodic box to model the impact of supernova (SN) explosions on the multi-phase interstellar medium (ISM) for initial densities $n = 0.5-30$ cm$^{-3}$ and SN rates $1-720$ Myr$^{-1}$. We include radiative cooling, diffuse heating, and the formation of molecular gas using a chemical network. The SNe explode either at random positions, at density peaks, or both. We further present a model combining thermal energy for resolved and momentum input for unresolved SNe. Random driving at high SN rates results in hot gas ($T\gtrsim 10^6$ K) filling $> 90$% of the volume. This gas reaches high pressures ($10^4 < P/k_\mathrm{B} < 10^7$ K cm$^{-3}$) due to the combination of SN explosions in the hot, low density medium and confinement in the periodic box. These pressures move the gas from a two-phase equilibrium to the single-phase, cold branch of the cooling curve. The molecular hydrogen dominates the mass ($>50$%), residing in small, dense clumps. Such a model might resemble the dense ISM in high-redshift galaxies. Peak driving results in huge radiative losses, producing a filamentary ISM with virtually no hot gas, and a small molecular hydrogen mass fraction ($\ll 1$%). Varying the ratio of peak to random SNe yields ISM properties in between the two extremes, with a sharp transition for equal contributions. The velocity dispersion in HI remains $\lesssim 10$ km s$^{-1}$ in all cases. For peak driving the velocity dispersion in H$_αあるふぁ$ can be as high as $70$ km s$^{-1}$ due to the contribution from young, embedded SN remnants. △ Less

Submitted 8 March, 2015; v1 submitted 31 October, 2014; originally announced November 2014.

Comments: 19 pages, 12 figures, 2 tables. Accepted for publication in MNRAS. Minor revisions to match published version

Exploring GLIMPSE Bubble N107: Multiwavelength Observations and Simulations

Authors: Vojtech Sidorin, Kevin A. Douglas, Jan Palous, Richard Wunsch, Sona Ehlerova

Abstract: Context. Bubble N107 was discovered in the infrared emission of dust in the Galactic Plane observed by the Spitzer Space Telescope (GLIMPSE survey: l ~ 51.0 deg, b ~ 0.1 deg). The bubble represents an example of shell-like structures found all over the Milky Way Galaxy. Aims. We aim to analyse the atomic and molecular components of N107, as well as its radio continuum emission. With the help of… ▽ More Context. Bubble N107 was discovered in the infrared emission of dust in the Galactic Plane observed by the Spitzer Space Telescope (GLIMPSE survey: l ~ 51.0 deg, b ~ 0.1 deg). The bubble represents an example of shell-like structures found all over the Milky Way Galaxy. Aims. We aim to analyse the atomic and molecular components of N107, as well as its radio continuum emission. With the help of numerical simulations, we aim to estimate the bubble age and other parameters which cannot be derived directly from observations. Methods. From the observations of the HI (I-GALFA) and 13CO (GRS) lines we derive the bubble's kinematical distance and masses of the atomic and molecular components. With the algorithm DENDROFIND, we decompose molecular material into individual clumps. From the continuum observations at 1420 MHz (VGPS) and 327 MHz (WSRT), we derive the radio flux density and the spectral index. With the numerical code ring, we simulate the evolution of stellar-blown bubbles similar to N107. Results. The total HI mass associated with N107 is 5.4E3 Msun. The total mass of the molecular component (a mixture of cold gasses of H2, CO, He and heavier elements) is 1.3E5 Msun, from which 4.0E4 Msun is found along the bubble border. We identified 49 molecular clumps distributed along the bubble border, with the slope of the clump mass function of -1.1. The spectral index of -0.30 of a strong radio source located apparently within the bubble indicates nonthermal emission, hence part of the flux likely originates in a supernova remnant, not yet catalogued. The numerical simulations suggest N107 is likely less than 2.25 Myr old. Since first supernovae explode only after 3 Myr or later, no supernova remnant should be present within the bubble. It may be explained if there is a supernova remnant in the direction towards the bubble, however not associated with it. △ Less

Submitted 23 February, 2014; originally announced February 2014.

Comments: 15 pages, 11 figures