Researcher profile

S. Walch

S. Walch contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2025arXiv

ALMAGAL V. Relations between the core populations and the parent clump physical properties

Context. The fragmentation of massive molecular clumps into smaller, potentially star-forming cores plays a key role in the processes of high-mass star formation. The ALMAGAL project offers high-resolution data to investigate these processes across various evolutionary stages in the Galactic plane. Aims. This study aims at correlating the fragmentation properties of massive clumps, obtained from ALMA observations, with their global physical parameters (e.g., mass, surface density, and temperature) and evolutionary indicators (such as luminosity-to-mass ratio and bolometric temperature) obtained from Herschel observations. It seeks to assess whether the cores evolve in number and mass in tandem with their host clumps, and to determine the possible factors influencing the formation of massive cores (M > 24M_\odot). Methods. We analyzed the masses of 6348 fragments, estimated from 1.4 mm continuum data for 1007 ALMAGAL clumps. Leveraging this unprecedentedly large data set, we evaluated statistical relationships between clump parameters, estimated over about 0.1 pc scales, and fragment properties, corresponding to scales of a few 1000 au, while accounting for potential biases related to distance and observational resolution. Our results were further compared with predictions from numerical simulations. Results. The fragmentation level correlates preferentially with clump surface density, supporting a scenario of density-driven fragmentation, whereas it does not show any clear dependence on total clump mass. Both the mass of the most massive core and the core formation efficiency show a broad range and increase on average by an order of magnitude in the intervals spanned by evolutionary indicators such as clump dust temperature and the luminosity-to-mass ratio. This suggests that core growth continues throughout the clump evolution, favoring clump-fed over core-fed theoretical scenarios.

preprint2022arXiv

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

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

preprint2022arXiv

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

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

preprint2022arXiv

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

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.

preprint2022arXiv

Protostellar Outflows: a window to the past

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

preprint2022arXiv

SILCC-Zoom: the dynamic balance in molecular cloud substructures

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

preprint2020arXiv

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

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

preprint2020arXiv

How fast do young star clusters expel their natal gas?: Estimating the upper limit of the gas expulsion time-scale

Formation of massive stars within embedded star clusters starts a complex interplay between their feedback, inflowing gas and stellar dynamics, which often includes close stellar encounters. Hydrodynamical simulations usually resort to substantial simplifications to model embedded clusters. Here, we address the simplification which approximates the whole star cluster by a single sink particle, which completely neglects the internal stellar dynamics. In order to model the internal stellar dynamics, we implement a Hermite predictor-corrector integration scheme to the hydrodynamic code FLASH. As we illustrate by a suite of tests, this integrator significantly outperforms the current leap-frog scheme, and it is able to follow the dynamics of small compact stellar systems without the necessity to soften the gravitational potential. We find that resolving individual massive stars instead of representing the whole cluster by a single energetic source has a profound influence on the gas component: for clusters of mass less than $\lesssim 3 \times 10^3 M_{\odot}$, it slows gas expulsion by a factor of $\approx 5$ to $\approx 1$ Myr, and it results in substantially more complex gas structures. With increasing cluster mass (up to $\approx 3\times 10^3 M_{\odot}$), the gas expulsion time-scale slightly decreases. However, more massive clusters ($\gtrsim 5\times 10^3 M_{\odot}$) are unable to clear their natal gas with photoionising radiation and stellar winds only if they form with a star formation efficiency (SFE) of $1/3$. This implies that the more massive clusters are either cleared with another feedback mechanism or they form with a SFE higher than $1/3$.

preprint2020arXiv

The hierarchical fragmentation of filaments and the role of sub-filaments

Recent observations have revealed the presence of small fibres or sub-filaments within larger filaments. We present a numerical fragmentation study of fibrous filaments investigating the link between cores and sub-filaments using hydrodynamical simulations performed with the moving-mesh code Arepo. Our study suggests that cores form in two environments: (i) as isolated cores, or small chains of cores, on a single sub-filament, or (ii) as an ensemble of cores located at the junction of sub-filaments. We term these isolated and hub cores respectively. We show that these core populations are statistically different from each other. Hub cores have a greater mean mass than isolated cores, and the mass distribution of hub cores is significantly wider than isolated cores. This fragmentation is reminiscent of parsec-scale hub-filament systems, showing that the combination of turbulence and gravity leads to similar fragmentation signatures on multiple scales, even within filaments. Moreover, the fact that fragmentation proceeds through sub-filaments suggests that there exists no characteristic fragmentation length-scale between cores. This is in opposition to earlier theoretical works studying fibre-less filaments which suggest a strong tendency towards the formation of quasi-periodically spaced cores, but in better agreement with observations. We also show tentative signs that global collapse of filaments preferentially form cores at both filament ends, which are more massive and dense than other cores.

preprint2020arXiv

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

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

preprint2019arXiv

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

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