Source author record

R. S. Klessen

R. S. Klessen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

45works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

45 published item(s)

preprint2026arXiv

PHANGS-JWST: the largest extragalactic molecular cloud catalog traced by polycyclic aromatic hydrocarbon emission

High-resolution JWST images of nearby spiral galaxies reveal polycyclic aromatic hydrocarbon (PAH) emission structures that trace molecular gas, including CO-dark regions. We identify ISM cloud structures in PHANGS-JWST 7.7 $μ$m PAH maps for 66 galaxies, smoothed to 30 pc and at native resolution, extracting 108,466 and 146,040 clouds, respectively. Molecular properties were inferred using a linear conversion from PAH to CO. Given the tendency for clouds in galaxy centers to overlap in velocity space, we opted to flag these and omit them from the analysis in this work. The remaining clouds correspond to giant molecular clouds, such as those detected in CO(2-1) emission by ALMA, or lower surface density clouds that either fall below the ALMA detection limits of existing maps or genuinely have no molecular counterpart. Cross-matching with ALMA CO maps at 90 pc in 27 galaxies shows that 41 % of PAH clouds have CO associations. The converted molecular properties vary little across environments, but the most massive clouds are preferentially found in spiral arms. Fitting lognormal mass distributions down to $2\times10^{3} M_{\odot}$ shows that spiral arms host the highest-mass clouds, consistent with enhanced formation in arm gravitational potentials. Cloud molecular surface densities decline by a factor of $\sim 1.5-2$ toward $2 - 3 R_{e}$. However, the trend largely varies in individual galaxies, with flat, decreasing, and even no trend as a function of galactocentric radius. Factors like large-scale processes and morphologies might influence the observed trends. We publish two catalogs online, one at the common resolution of 30 pc and another at the native resolution. We expect them to have broad utility for future PAH clouds, molecular clouds, and star formation studies.

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.

preprint2025arXiv

Anisotropy in the carbon monoxide (CO) line emission across the Milky Way's disk

We present a study of the $^{12}$CO(1-0) line emission anisotropy across the Milky Way's disk to examine the effect of stellar feedback and Galactic dynamics on the distribution of the dense interstellar medium. The Hessian matrix method is used to characterize the CO line emission distribution and identify the preferential orientation across line-of-sight velocity channels in the Dame et al. 2001 composite Galactic plane survey, which covers the Galactic latitude range $|b|<5^{\circ}$. The structures sampled with this tracer are predominantly parallel to the Galactic plane toward the inner Galaxy, in clear contrast with the predominantly perpendicular orientation of the structures traced by neutral atomic hydrogen (HI) emission toward the same regions. The analysis of the Galactic plane portions sampled at higher angular resolution with other surveys reveals that the alignment with the Galactic plane is also prevalent at smaller scales. We find no preferential orientation in the CO emission toward the outer Galaxy, in contrast with the preferential alignment with the Galactic plane displayed by HI in that portion of the Milky Way. We interpret these results as the combined effect of the decrease in mid-plane pressure with increasing Galactocentric radius and SN feedback lifting diffuse gas more efficiently than dense gas off the Galactic plane.

preprint2022arXiv

Linking stellar populations to HII regions across nearby galaxies: I. Constraining pre-supernova feedback from young clusters in NGC1672

One of the fundamental factors regulating the evolution of galaxies is stellar feedback. However, we still do not have strong observational constraints on the relative importance of the different feedback mechanisms (e.g. radiation, ionised gas pressure, stellar winds) in driving HII region evolution and molecular cloud disruption. In this letter, we constrain the relative importance of the various feedback mechanisms from young massive star populations by resolving HII regions across the disk of the nearby star-forming galaxy NGC 1672. We combine measurements of ionised gas nebular lines obtained by PHANGS-MUSE, with high-resolution imaging from the HST in both the narrow-band Hα and broad-band filters. We identify a sample of 40 isolated, compact HII regions in the HST Hα image, for which we measure the sizes that were previously unresolved in seeing-limited ground-based observations. Additionally, we identify the ionisation source(s) for each HII region from catalogues produced as part of the PHANGS-HST survey. We find that the HII regions investigated are mildly dominated by thermal or wind pressure, yet their elevation above the radiation pressure is within the expected uncertainty range. We see that radiation pressure provides a substantially higher contribution to the total pressure than previously found in the literature over similar size scales. In general, we find higher pressures within more compact HII regions, which is driven by the inherent size scaling relations of each pressure term, albeit with significant scatter introduced by the variation in the stellar population properties (e.g. luminosity, mass, age, metallicity). For nearby galaxies, here we provide a promising approach that could yield the statistics required to map out how the importance of different stellar feedback mechanisms evolve over the lifetime of an HII region.

preprint2022arXiv

Origin of supermassive black holes in massive metal-poor protoclusters

While large numbers of supermassive black holes have been detected at z>6, their origin is still essentially unclear. Numerical simulations have shown that the conditions for the classical direct collapse scenario are very restrictive and fragmentation is very difficult to be avoided. We thus consider here a more general case of a dense massive protostar cluster at low metallicity (<~ 10^{-3} Z_solar) embedded in gas. We estimate the mass of the central massive object, formed via collisions and gas accretion, considering the extreme cases of a logarithmically flat and a Salpeter-type initial mass function. Objects with masses of at least 10^4 solar could be formed for inefficient radiative feedback, whereas ~10^3 solar mass objects could be formed when the accretion time is limited via feedback. These masses will vary depending on the environment and could be considerably larger, particularly due to the continuous infall of gas into the cloud. As a result, one may form intermediate mass black holes of ~ 10^4 solar masses or more. Upcoming observations with the James Webb Space Telescope (JWST) and other observatories may help to detect such massive black holes and their environment, thereby shedding additional light on such a formation channel.

preprint2021arXiv

The Close AGN Reference Survey (CARS): No obvious signature of AGN feedback on star formation, but subtle trends

[Abridged] Active Galactic Nuclei (AGN) are thought to be responsible for the suppression of star formation in massive ~10$^{10}$ M$_\odot$ galaxies. While this process is a key feature in numerical simulations, it is not yet unambiguously confirmed in observational studies. Characterization of the star formation rate (SFR) in AGN host galaxies is challenging as AGN light contaminates most SFR tracers. We aim to obtain and compare SFR estimates from different tracers for AGN host galaxies in the Close AGN Reference Survey (CARS) to provide new observational insights. We construct integrated panchromatic spectral energy distributions (SED) to measure the FIR luminosity as a tracer for the recent (< 100 Myr) SFR. In addition, we use integral-field unit observation of the CARS targets to employ the H$α$ luminosity decontaminated by AGN excitation as a proxy for the current (< 5 Myr) SFR. We find that significant differences in specific SFR of the AGN host galaxies as compared with the larger galaxy population disappear once cold gas mass, in addition to stellar mass, is used to predict the SFR. We identify individual galaxies with a significant difference in their SFR which can be related to a recent enhancement or decline in their SFR history that might be related to various processes including interactions, gas consumption, outflows and AGN feedback. AGN can occur in various stages of galaxy evolution which makes it difficult to relate the SFR solely to the impact of the AGN. We do not find any strong evidence for global positive or negative AGN feedback in the CARS sample. However, there is tentative evidence that 1) the relative orientation of the AGN engine with respect to the host galaxies might alter the efficiency of AGN feedback and 2) the recent SFH is an additional tool to identify rapid changes in galaxy growth driven by the AGN or other processes.

preprint2020arXiv

Cloud formation in the atomic and molecular phase: HI self absorption (HISA) towards a Giant Molecular Filament

Molecular clouds form from the atomic phase of the interstellar medium. However, characterizing the transition between the atomic and the molecular interstellar medium (ISM) is a difficult observational task. Here we address cloud formation processes by combining HSIA with molecular line data. One scenario proposed by numerical simulations is that the column density probability density functions (N-PDF) evolves from a log-normal shape at early times to a power-law-like shape at later times. In this paper, we study the cold atomic component of the giant molecular filament GMF38a (d=3.4 kpc, length$\sim230$ pc). We identify an extended HISA feature, which is partly correlated with the 13CO emission. The peak velocities of the HISA and 13CO observations agree well on the eastern side of the filament, whereas a velocity offset of approximately 4 km s$^{-1}$ is found on the western side. The sonic Mach number we derive from the linewidth measurements shows that a large fraction of the HISA, which is ascribed to the cold neutral medium (CNM), is at subsonic and transonic velocities. The column density of the CNM is on the order of 10$^{20}$ to 10$^{21}$ cm$^{-2}$. The column density of molecular hydrogen is an order of magnitude higher. The N-PDFs from HISA (CNM), HI emission (WNM+CNM), and 13CO (molecular component) are well described by log-normal functions, which is in agreement with turbulent motions being the main driver of cloud dynamics. The N-PDF of the molecular component also shows a power law in the high column-density region, indicating self-gravity. We suggest that we are witnessing two different evolutionary stages within the filament. The eastern subregion seems to be forming a molecular cloud out of the atomic gas, whereas the western subregion already shows high column density peaks, active star formation and evidence of related feedback processes.

preprint2020arXiv

Formation of SMBH seeds in Pop III star clusters through collisions : the importance of mass loss

Runaway collisions in dense clusters may lead to the formation of supermassive black hole (SMBH) seeds, and this process can be further enhanced by accretion, as recent models of SMBH seed formation in Population III star clusters have shown. This may explain the presence of supermassive black holes already at high redshift, $z>6$. However, in this context, mass loss during collisions was not considered and could play an important role for the formation of the SMBH seed. Here, we study the effect of mass loss, due to collisions of protostars, in the formation and evolution of a massive object in a dense primordial cluster. We consider both constant mass loss fractions as well as analytic models based on the stellar structure of the collision components. Our calculations indicate that mass loss can significantly affect the final mass of the possible SMBH seed. Considering a constant mass loss of 5% for every collision, we can lose between 60-80% of the total mass that is obtained if mass loss were not considered. Using instead analytical prescriptions for mass loss, the mass of the final object is reduced by 15-40%, depending on the accretion model for the cluster we study. Altogether, we obtain masses of the order of $10^4M_{\odot}$, which are still massive enough to be SMBH seeds.

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'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

Less than the sum of its parts: the dust-corrected H$α$ luminosity of star-forming galaxies explored at different spatial resolutions with MaNGA and MUSE

The H$α$ and H$β$ emission line luminosities measured in a single integrated spectrum are affected in non-trivial ways by point-to-point variations in dust attenuation in a galaxy. This work investigates the impact of this variation when estimating global H$α$ luminosities corrected for the presence of dust by a global Balmer decrement. Analytical arguments show that the dust-corrected H$α$ luminosity is always underestimated when using the global H$α$/H$β$ flux ratio to correct for dust attenuation. We measure this effect on 156 face-on star-forming galaxies from the Mapping Nearby Galaxies at APO (MaNGA) survey. At 1-2 kpc spatial resolution, the effect is small but systematic, with the integrated dust-corrected H$α$ luminosity underestimated by $2$-$4$ per cent (and typically not more than by $10$ per cent), and depends on the specific star formation rate of the galaxy. Given the spatial resolution of MaNGA, these are lower limits for the effect. From Multi Unit Spectroscopic Explorer (MUSE) observations of NGC 628 with a resolution of 36 pc we find the discrepancy between the globally and the point-by-point dust-corrected H$α$ luminosity to be $14 \pm 1$ per cent, which may still underestimate the true effect. We use toy models and simulations to show that the true difference depends strongly on the spatial variance of the H$α$/H$β$ flux ratio, and on the slope of the relation between H$α$ luminosity and dust attenuation within a galaxy. Larger samples of higher spatial resolution observations are required to quantify the dependence of this effect as a function of galaxy properties.

preprint2020arXiv

Statistical mass function of prestellar cores from the density distribution of their natal clouds

The mass function of clumps observed in molecular clouds raises interesting theoretical issues, especially in its relation to the stellar initial mass function. We propose a statistical model of the mass function of prestellar cores (CMF), formed in self-gravitating isothermal clouds at a given stage of their evolution. The latter is characterized by the mass-density probability distribution function ($ρ$-PDF), which is a power-law with slope $q$. The variety of MCs is divided in ensembles according to the PDF slope and each ensemble is represented by a single spherical cloud. The cores are considered as elements of self-similar structure typical for fractal clouds and are modeled by spherical objects populating each cloud shell. Our model assumes relations between size, mass and density of the statistical cores. Out of them a core mass-density relationship $ρ\propto m^x$ is derived where $x=1/(1+q)$. We found that $q$ determines the existence or non-existence of a threshold density for core collapse. The derived general CMF is a power law of slope $-1$ while the CMF of gravitationally unstable cores has a slope $(-1 + x/2)$, comparable with the slopes of the high-mass part of the stellar initial mass function and of observational CMFs.

preprint2020arXiv

The effects of a background potential in star cluster evolution: a delay in the relaxation time-scale and runaway collision processes

Runaway stellar collisions in dense star clusters are invoked to explain the presence of very massive stars or blue stragglers in the center of those systems. This process has also been explored for the first star clusters in the Universe and shown to yield stars that may collapse at some points into an intermediate mass black hole. Although the early evolution of star clusters requires the explicit modeling of the gas out of which the stars form, these calculations would be extremely time-consuming and often the effects of the gas can be accurately treated by including a background potential to account for the extra gravitational force. We apply this approximation to model the early evolution of the first dense star clusters formed in the Universe by performing $N$-body simulations, our goal is to understand how the additional gravitational force affects the growth of a very massive star through stellar mergers in the central parts of the star cluster. Our results show that the background potential increases the velocities of the stars, causing an overall delay in the evolution of the clusters and in the runaway growth of a massive star at the center. The population of binary stars is lower due to the increased kinetic energy of the stars, initially reducing the number of stellar collisions, and we show that relaxation processes are also affected. Despite these effects, the external potential enhances the mass of the merger product by a factor $\sim$2 if the collisions are maintained for long times.

preprint2020arXiv

Ubiquitous velocity fluctuations throughout the molecular interstellar medium

The density structure of the interstellar medium (ISM) determines where stars form and release energy, momentum, and heavy elements, driving galaxy evolution. Density variations are seeded and amplified by gas motion, but the exact nature of this motion is unknown across spatial scale and galactic environment. Although dense star-forming gas likely emerges from a combination of instabilities, convergent flows, and turbulence, establishing the precise origin is challenging because it requires quantifying gas motion over many orders of magnitude in spatial scale. Here we measure the motion of molecular gas in the Milky Way and in nearby galaxy NGC 4321, assembling observations that span an unprecedented spatial dynamic range ($10^{-1}{-}10^3$ pc). We detect ubiquitous velocity fluctuations across all spatial scales and galactic environments. Statistical analysis of these fluctuations indicates how star-forming gas is assembled. We discover oscillatory gas flows with wavelengths ranging from $0.3{-}400$ pc. These flows are coupled to regularly-spaced density enhancements that likely form via gravitational instabilities. We also identify stochastic and scale-free velocity and density fluctuations, consistent with the structure generated in turbulent flows. Our results demonstrate that ISM structure cannot be considered in isolation. Instead, its formation and evolution is controlled by nested, interdependent flows of matter covering many orders of magnitude in spatial scale.

preprint2019arXiv

Feedback in W49A diagnosed with radio recombination lines and models

We present images of radio recombination lines (RRLs) at wavelengths around 17 cm from the star-forming region W49A to determine the kinematics of ionized gas in the THOR survey (The HI/OH/Recombination line survey of the inner Milky Way) at an angular resolution of 16.8"x13.8". The distribution of ionized gas appears to be affected by feedback processes from the star clusters in W49A. The velocity structure of the RRLs shows a complex behavior with respect to the molecular gas. We find a shell-like distribution of ionized gas as traced by RRL emission surrounding the central cluster of OB stars in W49A. We describe the evolution of the shell with the recent feedback model code WARPFIELD that includes the important physical processes and has previously been applied to the 30 Doradus region in the Large Magellanic Cloud. The cloud structure and dynamics of W49A are in agreement with a feedback-driven shell that is re-collapsing. The shell may have triggered star formation in other parts of W49A. We suggest that W49A is a potential candidate for star formation regulated by feedback-driven and re-collapsing shells.

preprint2019arXiv

Strong excess Faraday rotation on the Inside of the Sagittarius spiral arm

We present first results for Faraday rotation of compact polarized sources (1 to 2 GHz continuum) in The HI/OH/Recombination line (THOR) survey of the inner Galaxy. In the Galactic longitude range 39 degr < l < 52 degr, we find rotation measures in the range -310 rad/m2 < RM < +4219 rad/m2, with the highest values concentrated within a degree of l = 48 degrees at the Sagittarius arm tangent. Most of the high RMs arise in diffuse plasma, along lines of sight that do not intersect HII regions. For l > 49 degr, RM drops off rapidly, while at l < 47 degr, the mean RM is higher with a larger standard deviation than at l > 49 degr. We attribute the RM structure to the compressed diffuse Warm Ionized Medium in the spiral arm, upstream of the major star formation regions. The Sagittarius arm acts as a significant Faraday screen inside the Galaxy. This has implications for models of the Galactic magnetic field and the expected amount of Faraday rotation of Fast Radio Bursts from their host galaxies. We emphasize the importance of sensitivity to high Faraday depth in future polarization surveys.

preprint2019arXiv

The HI/OH/Recombination line survey of the inner Milky Way (THOR): data release 2 and HI overview

With the $Karl~G.~Jansky$ Very Large Array (VLA) in C-configuration, we observed a large portion of the first Galactic quadrant ($l=14.0-67.4^\circ $ and $\lvert b \rvert \leq 1.25^\circ $) achieving an angular resolution of $\leq 40^{\prime\prime}$. At $L$ Band, the WIDAR correlator at the VLA was set to cover the 21~cm HI line, four OH transitions, a series of H$nα$ radio recombination lines (RRLs; $n=151$ to 186), and eight 128~MHz wide continuum spectral windows (SPWs) simultaneously. The HI emission shows clear filamentary substructures at negative velocities with low velocity crowding. The emission at positive velocities is more smeared-out likely due to higher spatial and velocity crowding of structures at the positive velocities. Comparing to the spiral arm model of the Milky Way, the atomic gas follows the Sagittarius and Perseus Arm well but with significant material in the inter-arm regions. With the C-configuration-only HI+continuum data, we produced a HI optical depth map of the THOR areal coverage from 228 absorption spectra with the nearest-neighbor method. With this $τ$ map, we corrected the HI emission for optical depth and the derived column density is 38% higher than the column density with optically thin assumption. The total HI mass with optical depth correction in the survey region is 4.7$\times10^8~M_\odot$, 31% more than the mass derived assuming the emission is optically thin. If we apply this 31% correction to the whole Milky Way, the total atomic gas mass would be 9.4-10.5$\times 10^9~M_\odot$. Comparing the HI with existing CO data, we find a significant increase in the atomic-to-molecular gas ratio from the spiral arms to the inter-arm regions.

preprint2019arXiv

WARPFIELD Population Synthesis: The physics of (extra-)Galactic star formation and feedback driven cloud structure and emission from sub-to-kpc scales

We present a novel method to model galactic scale star formation and the resulting emission from star clusters and the multi-phase interstellar medium. We combine global parameters, such as SFR and CMF, with {\sc warpfield} which provides a description of the feedback-driven evolution of individual star-forming regions. Our approach includes stellar evolution, stellar winds, radiation pressure, supernovae, all of which couple to the dynamical evolution of the parental cloud in a highly non-linear fashion. The heating of diffuse gas and dust is calculated self-consistently with the age, mass and density dependent escape fractions of the local star-forming regions. From this we construct the interstellar radiation field at any point in the galaxy, and we employ the multi-frequency Monte Carlo radiative transfer code {\sc polaris} to produce synthetic emission maps for the one-to-one comparison with observational data. We demonstrate the capabilities of our approach by applying the method to a Milky Way like galaxy built-up in a high-resolution cosmological MHD simulation. We give three examples. First, we compute the multi-scale distribution of electron $n_{e^-}$ and $T_{e^{-}}$ and synthesize the MW all-sky H$α$ emission. We use a multipole expansion method to show that the resulting maps are consistent with observations. Second, we predict the expected \SIII 9530~Å emission. This line is a key target of several planned large survey programs. It suffers less extinction than other diagnostic lines and provides information about star formation in very dense environments that are otherwise observationally not readily accessible. Third, we explore the effects of differential extinction as seen by an extra-galactic observer, and discuss the consequences for the correct interpretation of \Ha emission as a star-formation rate tracer at different viewing angles.(abridged)

preprint2016arXiv

Continuum sources from the THOR survey between 1 and 2 GHz

We carried out a large program with the Karl G. Jansky Very Large Array (VLA): "THOR: The HI, OH, Recombination line survey of the Milky Way". We observed a significant portion of the Galactic plane in the first quadrant of the Milky Way in the 21cm HI line, 4 OH transitions, 19 radio recombination lines, and continuum from 1 to 2 GHz. In this paper we present a catalog of the continuum sources in the first half of the survey (l=14.0-37.9deg and l=47.1-51.2deg, |b|<1.1deg) at a spatial resolution of 10-25", with a spatially varying noise level of ~0.3-1 mJy/beam. The catalog contains ~4400 sources. Around 1200 of these are spatially resolved, and ~1000 are possible artifacts, given their low signal-to-noise ratios. Since the spatial distribution of the unresolved objects is evenly distributed and not confined to the Galactic plane, most of them are extragalactic. Thanks to the broad bandwidth of the observations from 1 to 2 GHz, we are able to determine a reliable spectral index for ~1800 sources. The spectral index distribution reveals a double-peaked profile with maxima at spectral indices of alpha = -1 and alpha = 0 , corresponding to steep declining and flat spectra, respectively. This allows us to distinguish between thermal and non-thermal emission, which can be used to determine the nature of each source. We examine the spectral index of ~300 known HII regions, for which we find thermal emission with spectral indices around alpha = 0. In contrast, supernova remnants (SNR) show non-thermal emission with alpha = -0.5 and extragalactic objects generally have a steeper spectral index of alpha = -1. Using the spectral index information of the THOR survey, we investigate potential SNR candidates. We classify the radiation of four SNR candidates as non-thermal, and for the first time, we provide strong evidence for the SNR origin of these candidates.

preprint2016arXiv

Understanding star formation in molecular clouds III. Probability distribution functions of molecular lines in Cygnus X

Column density (N) PDFs serve as a powerful tool to characterize the physical processes that influence the structure of molecular clouds. Star-forming clouds can best be characterized by lognormal PDFs for the lower N range and a power-law tail for higher N, commonly attributed to turbulence and self-gravity and/or pressure, respectively. We report here on PDFs obtained from observations of 12CO, 13CO, C18O, CS, and N2H+ in the Cygnus X North region and compare to a PDF derived from dust observations with the Herschel satellite. The PDF of 12CO is lognormal for Av~1-30, but is cut for higher Av due to optical depth effects. The PDFs of C18O and 13CO are mostly lognormal up for Av~1-15, followed by excess up to Av~40. Above that value, all CO PDFs drop, most likely due to depletion. The high density tracers CS and N2H+ exhibit only a power law distribution between Av~15 and 400, respectively. The PDF from dust is lognormal for Av~3-15 and has a power-law tail up to Av~500. Absolute values for the molecular line column densities are, however, rather uncertain due to abundance and excitation temperature variations. Taken the dust PDF face value, we 'calibrate' the molecular line PDF of CS to the one of the dust and determined an abundance [CS]/[H2] of 10^-9. The slopes of the power-law tails of the CS, N2H+, and dust PDFs are consistent with free-fall collapse of filaments and clumps. A quasi static configuration of filaments and clumps can possibly also account for the observed N-PDFs, as long as they have a sufficiently condensed density structure and external ram pressure by gas accretion is provided. The somehow flatter slopes of N2H+ and CS can reflect an abundance change and/or subthermal excitation at low column densities.

preprint2015arXiv

Detection of two power-law tails in the probability distribution functions of massive GMCs

We report the novel detection of complex high-column density tails in the probability distribution functions (PDFs) for three high-mass star-forming regions (CepOB3, MonR2, NGC6334), obtained from dust emission observed with Herschel. The low column density range can be fit with a lognormal distribution. A first power-law tail starts above an extinction (Av) of ~6-14. It has a slope of alpha=1.3-2 for the rho~r^-alpha profile for an equivalent density distribution (spherical or cylindrical geometry), and is thus consistent with free-fall gravitational collapse. Above Av~40, 60, and 140, we detect an excess that can be fitted by a flatter power law tail with alpha>2. It correlates with the central regions of the cloud (ridges/hubs) of size ~1 pc and densities above 10^4 cm^-3. This excess may be caused by physical processes that slow down collapse and reduce the flow of mass towards higher densities. Possible are: 1. rotation, which introduces an angular momentum barrier, 2. increasing optical depth and weaker cooling, 3. magnetic fields, 4. geometrical effects, and 5. protostellar feedback. The excess/second power-law tail is closely linked to high-mass star-formation though it does not imply a universal column density threshold for the formation of (high-mass) stars.

preprint2015arXiv

Fervent: Chemistry-coupled, ionising and non-ionising radiative feedback in magnetohydrodynamical simulations

We introduce a radiative transfer code module for the magnetohydrodynamical adaptive mesh refinement code FLASH 4. It is coupled to an efficient chemical network which explicitly tracks the three hydrogen species H, H_2, H+ as well as C+ and CO. The module is geared towards modeling all relevant thermal feedback processes of massive stars, and is able to follow the non-equilibrium time-dependent thermal and chemical state of the present-day interstellar medium as well as that of dense molecular clouds. We describe in detail the implementation of all relevant thermal stellar feedback mechanisms, i.e. photoelectric, photoionization and H_2 dissociation heating as well as pumping of molecular hydrogen by UV photons. All included radiative feedback processes are extensively tested. We also compare our module to dedicated photon-dominated region (PDR) codes and find good agreement in our modeled hydrogen species once our radiative transfer solution reaches equilibrium. In addition, we show that the implemented radiative feedback physics is insensitive to the spatial resolution of the code and show under which conditions it is possible to obtain well-converged evolution in time. Finally, we briefly explore the robustness of our scheme for treating combined ionizing and non-ionizing radiation.

preprint2015arXiv

Modelling the supernova-driven ISM in different environments

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.

preprint2015arXiv

Star-formation rates from young-star counts and the structure of the ISM across the NGC346/N66 complex in the SMC

The rate at which interstellar gas is converted into stars, and its dependence on environment, is one of the pillars on which our understanding of the visible Universe is build. We present a comparison of the surface density of young stars (Sigma_*) and dust surface density (Sigma_d) across NGC346 (N66) in 115 independent pixels of 6x6 pc^2. We find a correlation between Sigma_* and Sigma_d with a considerable scatter. A power law fit to the data yields a steep relation with an exponent of 2.6+-0.2. We convert Sigma_d to gas surface density (Sigma_g) and Sigma_* to star formation rate (SFR) surface densities (Sigma_SFR), using simple assumptions for the gas-to-dust mass ratio and the duration of star formation. The derived total SFR (4+-1 10^-3 M_sun/yr) is consistent with SFR estimated from the Ha emission integrated over the Ha nebula. On small scales the Sigma_SFR derived using Ha systematically underestimates the count-based Sigma_SFR, by up to a factor of 10. This is due to ionizing photons escaping the area, where the stars are counted. We find that individual 36 pc^2 pixels fall systematically above integrated disc-galaxies in the Schmidt-Kennicutt diagram by on average a factor of ~7. The NGC346 average SFR over a larger area (90 pc radius) lies closer to the relation but remains high by a factor of ~3. The fraction of the total mass (gas plus young stars) locked in young stars is systematically high (~10 per cent) within the central 15 pc and systematically lower outside (2 per cent), which we interpret as variations in star formation efficiency. The inner 15 pc is dominated by young stars belonging to a centrally condensed cluster, while the outer parts are dominated by a dispersed population. Therefore, the observed trend could reflect a change of star formation efficiency between clustered and non-clustered star-formation.

preprint2015arXiv

THOR - The HI, OH, Recombination Line Survey of the Milky Way - The pilot study: HI observations of the giant molecular cloud W43

To study the atomic, molecular and ionized emission of Giant Molecular Clouds (GMCs), we have initiated a Large Program with the VLA: 'THOR - The HI, OH, Recombination Line survey of the Milky Way'. We map the 21cm HI line, 4 OH lines, 19 H_alpha recombination lines and the continuum from 1 to 2 GHz of a significant fraction of the Milky Way (l=15-67deg, |b|<1deg) at ~20" resolution. In this paper, we focus on the HI emission from the W43 star-formation complex. Classically, the HI 21cm line is treated as optically thin with column densities calculated under this assumption. This might give reasonable results for regions of low-mass star-formation, however, it is not sufficient to describe GMCs. We analyzed strong continuum sources to measure the optical depth, and thus correct the HI 21cm emission for optical depth effects and weak diffuse continuum emission. Hence, we are able to measure the HI mass of W43 more accurately and our analysis reveals a lower limit of M=6.6x10^6 M_sun, which is a factor of 2.4 larger than the mass estimated with the assumption of optically thin emission. The HI column densities are as high as N(HI)~150 M_sun/pc^2 ~ 1.9x10^22 cm^-2, which is an order of magnitude higher than for low mass star formation regions. This result challenges theoretical models that predict a threshold for the HI column density of ~10 M_sun/pc^2, at which the formation of molecular hydrogen should set in. By assuming an elliptical layered structure for W43, we estimate the particle density profiles. While at the cloud edge atomic and molecular hydrogen are well mixed, the center of the cloud is strongly dominated by molecular hydrogen. We do not identify a sharp transition between hydrogen in atomic and molecular form. Our results are an important characterization of the atomic to molecular hydrogen transition in an extreme environment and challenges current theoretical models.

preprint2015arXiv

TOPoS: II. On the bimodality of carbon abundance in CEMP stars. Implications on the early chemical evolution of galaxies

In the course of the TOPoS (Turn Off Primordial Stars) survey, aimed at discovering the lowest metallicity stars, we have found several carbon-enhanced metal-poor (CEMP) stars. We here present our analysis of six CEMP stars. Calcium and carbon are the only elements that can be measured in all six stars. The range is -5.0<=[Ca/H]< -2.1 and 7.12<=A(C)<=8.65. For star SDSS J1742+2531 we were able to detect three FeI lines from which we deduced [Fe/H]=-4.80, from four CaII lines we derived [Ca/H]=-4.56, and from synthesis of the G-band we derived A(C)=7.26. For SDSS J1035+0641 we were not able to detect any iron lines, yet we could place a robust (3sigma) upper limit of [Fe/H]< -5.0 and measure the Ca abundance, with [Ca/H]=-5.0, and carbon, A(C)=6.90. No lithium is detected in the spectrum of SDSS J1742+2531 or SDSS J1035+0641, which implies a robust upper limit of A(Li)<1.8 for both stars. Our measured carbon abundances confirm the bimodal distribution of carbon in CEMP stars, identifying a high-carbon band and a low-carbon band. We propose an interpretation of this bimodality according to which the stars on the high-carbon band are the result of mass transfer from an AGB companion, while the stars on the low-carbon band are genuine fossil records of a gas cloud that has also been enriched by a faint supernova (SN) providing carbon and the lighter elements. (Abridged)

preprint2015arXiv

Turbulence and its effect on protostellar disk formation

We analyse simulations of turbulent, magnetised molecular cloud cores focussing on the formation of Class 0 stage protostellar discs and the physical conditions in their surroundings. We show that for a wide range of initial conditions Keplerian discs are formed in the Class 0 stage already. Furthermore, we show that the accretion of mass and angular momentum in the surroundings of protostellar discs occurs in a highly anisotropic manner, by means of a few narrow accretion channels. The magnetic field structure in the vicinity of the discs is highly disordered, revealing field reversals up to distances of 1000 AU. These findings demonstrate that as soon as even mild turbulent motions are included, the classical disc formation scenario of a coherently rotating environment and a well-ordered magnetic field breaks down.

preprint2015arXiv

Understanding star formation in molecular clouds II. Signatures of gravitational collapse of IRDCs

We analyse column density and temperature maps derived from Herschel dust continuum observations of a sample of massive infrared dark clouds (G11.11-0.12, G18.82-0.28, G28.37+0.07, G28.53-0.25). We disentangle the velocity structure of the clouds using 13CO 1-0 and 12CO 3-2 data, showing that these IRDCs are the densest regions in massive giant molecular clouds and not isolated features. The probability distribution function (PDF) of column densities for all clouds have a power-law distribution over all (high) column densities, regardless of the evolutionary stage of the cloud: G11.11-0.12, G18.82-0.28, and G28.37+0.07 contain (proto)-stars, while G28.53-0.25 shows no signs of star formation. This is in contrast to the purely log-normal PDFs reported for near/mid-IR extinction maps. We only find a log-normal distribution for lower column densities, if we perform PDFs of the column density maps of the whole GMC in which the IRDCs are embedded. By comparing the PDF slope and the radial column density profile, we attribute the power law to the effect of large-scale gravitational collapse and to local free-fall collapse of pre- and protostellar cores. Independent from the PDF analysis, we find infall signatures in the spectral profiles of 12CO for G28.37+0.07 and G11.11-0.12, supporting the scenario of gravitational collapse. IRDCs are the densest regions within GMCs, which may be the progenitors of massive stars or clusters. At least some of the IRDCs are probably the same features as ridges (high column density regions with N>1e23 cm-2 over small areas), which were defined for nearby IR-bright GMCs. Because IRDCs are only confined to the densest (gravity dominated) cloud regions, the PDF constructed from this kind of a clipped image does not represent the (turbulence dominated) low column density regime of the cloud.

preprint2014arXiv

Accretion and magnetic field morphology around Class 0 stage protostellar discs

We analyse simulations of turbulent, magnetised molecular cloud cores focussing on the formation of Class 0 stage protostellar discs and the physical conditions in their surroundings. We show that for a wide range of initial conditions Keplerian discs are formed in the Class 0 stage already. In particular, we show that even subsonic turbulent motions reduce the magnetic braking efficiency sufficiently in order to allow rotationally supported discs to form. We therefore suggest that already during the Class 0 stage the fraction of Keplerian discs is significantly higher than 50%, consistent with recent observational trends but significantly higher than predictions based on simulations with misaligned magnetic fields, demonstrating the importance of turbulent motions for the formation of Keplerian discs. We show that the accretion of mass and angular momentum in the surroundings of protostellar discs occurs in a highly anisotropic manner, by means of a few narrow accretion channels. The magnetic field structure in the vicinity of the discs is highly disordered, revealing field reversals up to distances of 1000 AU. These findings demonstrate that as soon as even mild turbulent motions are included, the classical disc formation scenario of a coherently rotating environment and a well-ordered magnetic field breaks down. Hence, it is highly questionable to assess the magnetic braking efficiency based on non-turbulent collapse simulation. We strongly suggest that, in addition to the global magnetic field properties, the small-scale accretion flow and detailed magnetic field structure have to be considered in order to assess the likelihood of Keplerian discs to be present.

preprint2013arXiv

Flickering of 1.3 cm Sources in Sgr B2: Towards a Solution to the Ultracompact HII Region Lifetime Problem

Accretion flows onto massive stars must transfer mass so quickly that they are themselves gravitationally unstable, forming dense clumps and filaments. These density perturbations interact with young massive stars, emitting ionizing radiation, alternately exposing and confining their HII regions. As a result, the HII regions are predicted to flicker in flux density over periods of decades to centuries rather than increasing monotonically in size as predicted by simple Spitzer solutions. We have recently observed the Sgr B2 region at 1.3 cm with the VLA in its three hybrid configurations (DnC, CnB and BnA) at a resolution of 0.25''. These observations were made to compare in detail with matched continuum observations from 1989. At 0.25'' resolution, Sgr B2 contains 41 UC HII regions, 6 of which are hypercompact. The new observations of Sgr B2 allow comparison of relative peak flux densites for the HII regions in Sgr B2 over a 23 year time baseline (1989-2012) in one of the most source-rich massive star forming regions in the Milky Way. The new 1.3 cm continuum images indicate that four of the 41 UC HII regions exhibit significant changes in their peak flux density, with one source (K3) dropping in peak flux density, and the other 3 sources (F10.303, F1 and F3) increasing in peak flux density. The results are consistent with statistical predictions from simulations of high mass star formation, suggesting that they offer a solution to the lifetime problem for ultracompact HII regions.

preprint2013arXiv

Rotation and Internal Structure of Population III Protostars

We analyze the cosmological simulations performed in the recent work of Greif et al. (2012), which followed the early growth and merger history of Pop III stars while resolving scales as small as 0.05 R_sol. This is the first set of cosmological simulations to self-consistently resolve the rotation and internal structure of Pop III protostars. We find that Pop III stars form under significant rotational support which is maintained for the duration of the simulations. The protostellar surfaces spin from ~50% to nearly 100% of Keplerian rotational velocity. These rotation rates persist after experiencing multiple stellar merger events. In the brief time period simulated (~ 10 yr), the protostars show little indication of convective instability, and their properties furthermore show little correlation with the properties of their host minihaloes. If Pop III protostars within this range of environments generally form with high degrees of rotational support, and if this rotational support is maintained for a sufficient amount of time, this has a number of crucial implications for Pop III evolution and nucleosynthesis, as well as the possibility for Pop III pair-instability supernovae, and the question of whether the first stars produced gamma-ray bursts.

preprint2013arXiv

Turbulence-induced disc formation in strongly magnetised cloud cores

We present collapse simulations of strongly magnetised, turbulent molecular cloud cores with masses ranging from 2.6 to 1000 M_sun in order to study the influence of the initial conditions on the turbulence-induced disc formation mechanism proposed recently by Seifried et al. 2012. We find that Keplerian discs are formed in all cases independently of the core mass, the strength of turbulence, or the presence of global rotation. The discs appear within a few kyr after the formation of the protostar, are 50 - 150 AU in size, and have masses between 0.05 and a few 0.1 M_sun. During the formation of the discs the mass-to-flux ratio stays well below the critical value of 10 for Keplerian disc formation. Hence, flux-loss alone cannot explain the formation of Keplerian discs. The formation of rotationally supported discs at such early phases is rather due to the disordered magnetic field structure and due to turbulent motions in the surroundings of the discs, two effects lowering the classical magnetic braking efficiency. Binary systems occurring in the discs are mainly formed via the disc capturing mechanism rather than via disc fragmentation, which is largely suppressed by the presence of magnetic fields.

preprint2012arXiv

Disc formation in turbulent cloud cores: Circumventing the magnetic braking catastrophe

We present collapse simulations of strongly magnetised, 100 M_sun, turbulent cloud cores. Around the protostars formed during the collapse Keplerian discs with typical sizes of up to 100 AU build up in contrast to previous simulations neglecting turbulence. Analysing the condensations in which the discs form, we show that the magnetic flux loss is not sufficient to explain the build-up of Keplerian discs. The average magnetic field is strongly inclined to the disc which might reduce the magnetic braking efficiency. However, the main reason for the reduced magnetic braking efficiency is the highly disordered magnetic field in the surroundings of the discs. Furthermore, due to the lack of a coherently rotating structure in the turbulent environment of the disc no toroidal magnetic field necessary for angular momentum extraction can build up. Simultaneously the angular momentum inflow remains high due to local shear flows created by the turbulent motions. We suggest that the "magnetic braking catastrophe" is an artefact of the idealised non-turbulent initial conditions and that turbulence provides a natural mechanism to circumvent this problem.

preprint2012arXiv

Disc formation in turbulent massive cores: Circumventing the magnetic braking catastrophe

We present collapse simulations of 100 M_{\sun}, turbulent cloud cores threaded by a strong magnetic field. During the initial collapse phase filaments are generated which fragment quickly and form several protostars. Around these protostars Keplerian discs with typical sizes of up to 100 AU build up in contrast to previous simulations neglecting turbulence. We examine three mechanisms potentially responsible for lowering the magnetic braking efficiency and therefore allowing for the formation of Keplerian discs. Analysing the condensations in which the discs form, we show that the build-up of Keplerian discs is neither caused by magnetic flux loss due to turbulent reconnection nor by the misalignment of the magnetic field and the angular momentum. It is rather a consequence of the turbulent surroundings of the disc which exhibit no coherent rotation structure while strong local shear flows carry large amounts of angular momentum. We suggest that the "magnetic braking catastrophe", i.e. the formation of sub-Keplerian discs only, is an artefact of the idealised non-turbulent initial conditions and that turbulence provides a natural mechanism to circumvent this problem.

preprint2012arXiv

Magnetic fields during the early stages of massive star formation - II. A generalised outflow criterion

Numerical simulations of outflows formed during the collapse of 100 M_sun cloud cores are presented. We derive a generalised criterion from MHD wind theory to analyse the launching mechanism of these outflows. The criterion is successfully applied to the whole outflow structure and cases with sub-Keplerian disc rotation. It allows us to decide whether an outflow is driven centrifugally or by the toroidal magnetic pressure. We show that quantities such as the magnetic field line inclination or the ratio of the toroidal to poloidal magnetic field alone are insufficient to determine the driving mechanism of outflows. By performing 12 runs we are able to study the influence of the initial conditions on the properties of outflows around massive protostars in detail. Our simulations reveal a strong effect of the magnetic field strength on the outflow morphology. In runs with weak fields or high rotational energies, well collimated, fast jets are observed whereas for strong fields poorly collimated, low-velocity outflows are found. We show that the occurrence of a fast jet is coupled to the existence of a Keplerian protostellar disc. Despite the very different morphologies all outflows are launched from the discs by centrifugal acceleration with the toroidal magnetic field increasingly contributing to the gas acceleration further away from the discs. The poor collimation of the outflows in runs with strong fields is a consequence of the weak hoop stresses. This in turn is caused by the slow build-up of a toroidal field due to sub-Keplerian disc rotation. The mass ejection/accretion ratios scatter around a mean of 0.3 in accordance with observational and analytical results. We suggest an evolutionary scenario for the earliest stage of massive star formation in which initially poorly collimated outflows develop which progressively get better collimated due to the generation of fast jets.

preprint2012arXiv

The Density Variance--Mach Number Relation in Supersonic Turbulence: I. Isothermal, magnetised gas

It is widely accepted that supersonic, magnetised turbulence plays a fundamental role for star formation in molecular clouds. It produces the initial dense gas seeds out of which new stars can form. However, the exact relation between gas compression, turbulent Mach number, and magnetic field strength is still poorly understood. Here, we introduce and test an analytical prediction for the relation between the density variance and the root-mean-square Mach number in supersonic, isothermal, magnetised turbulent flows. We approximate the density and velocity structure of the interstellar medium as a superposition of shock waves. We obtain the density contrast considering the momentum equation for a single magnetised shock and extrapolate this result to the entire cloud. Depending on the field geometry, we then make three different assumptions based on observational and theoretical constraints: B independent of density, B proportional to the root square of the density and B proportional to the density. We test the analytically derived density variance--Mach number relation with numerical simulations, and find that for B proportional to the root square of the density, the variance in the logarithmic density contrast, $σ_{\ln ρ/ρ_0}^2=\ln[1+b^2\mathscr{M}^2β_0/(β_0+1)]$, fits very well to simulated data with turbulent forcing parameter b=0.4, when the gas is super-Alfvénic. However, this result breaks down when the turbulence becomes trans-Alfvénic or sub-Alfvénic, because in this regime the turbulence becomes highly anisotropic. Our density variance--Mach number relations simplify to the purely hydrodynamic relation as the ratio of thermal to magnetic pressure $β_0$ approaches infinite.

preprint2011arXiv

Collapse, outflows and fragmentation of massive, turbulent and magnetized prestellar barotropic cores

Stars and more particularly massive stars, have a drastic impact on galaxy evolution. Yet the conditions in which they form and collapse are still not fully understood. In particular, the influence of the magnetic field on the collapse of massive clumps is relatively unexplored, it is thus of great relevance in the context of the formation of massive stars to investigate its impact. We perform high resolution, MHD simulations of the collapse of hundred solar masses, turbulent and magnetized clouds, using the adaptive mesh refinement code RAMSES. We compute various quantities such as mass distribution, magnetic field and angular momentum within the collapsing core and study the episodic outflows and the fragmentation that occurs during the collapse. The magnetic field has a drastic impact on the cloud evolution. We find that magnetic braking is able to substantially reduce the angular momentum in the inner part of the collapsing cloud. Fast and episodic outflows are being launched with typical velocities of the order of 3-5 km s$^{-1}$ although the highest velocities can be as high as 30-40 km s$^{-1}$. The fragmentation in several objects, is reduced in substantially magnetized clouds with respect to hydrodynamical ones by a factor of the order of 1.5-2. We conclude that magnetic fields have a significant impact on the evolution of massive clumps. In combination with radiation, magnetic fields largely determine the outcome of massive core collapse. We stress that numerical convergence of MHD collapse is a challenging issue. In particular, numerical diffusion appears to be important at high density therefore possibly leading to an over-estimation of the number of fragments.

preprint2011arXiv

Magnetic fields during the early stages of massive star formation - I. Accretion and disk evolution

We present simulations of collapsing 100 M_\sun mass cores in the context of massive star formation. The effect of variable initial rotational and magnetic energies on the formation of massive stars is studied in detail. We focus on accretion rates and on the question under which conditions massive Keplerian disks can form in the very early evolutionary stage of massive protostars. For this purpose, we perform 12 simulations with different initial conditions extending over a wide range in parameter space. The equations of magnetohydrodynamics (MHD) are solved under the assumption of ideal MHD. We find that the formation of Keplerian disks in the very early stages is suppressed for a mass-to-flux ratio normalised to the critical value μbelow 10, in agreement with a series of low-mass star formation simulations. This is caused by very efficient magnetic braking resulting in a nearly instantaneous removal of angular momentum from the disk. For weak magnetic fields, corresponding to μ> 10, large-scale, centrifugally supported disks build up with radii exceeding 100 AU. A stability analysis reveals that the disks are supported against gravitationally induced perturbations by the magnetic field and tend to form single stars rather than multiple objects. We find protostellar accretion rates of the order of a few 10^-4 M_\sun yr^-1 which, considering the large range covered by the initial conditions, vary only by a factor of ~ 3 between the different simulations. We attribute this fact to two competing effects of magnetic fields. On the one hand, magnetic braking enhances accretion by removing angular momentum from the disk thus lowering the centrifugal support against gravity. On the other hand, the combined effect of magnetic pressure and magnetic tension counteracts gravity by exerting an outward directed force on the gas in the disk thus reducing the accretion onto the protostars.

preprint2011arXiv

The role of stellar collisions for the formation of massive stars

We use direct N-body simulations of gas embedded star clusters to study the importance of stellar collisions for the formation and mass accretion history of high-mass stars. Our clusters start in virial equilibrium as a mix of gas and proto-stars. Proto-stars then accrete matter using different mass accretion rates and the amount of gas is reduced in the same way as the mass of stars increases. During the simulations we check for stellar collisions and we investigate the role of these collisions for the build-up of high-mass stars and the formation of runaway stars. We find that a significant number of collisions only occur in clusters with initial half-mass radii r_h < 0.1 pc. After emerging from their parental gas clouds, such clusters end up too compact compared to observed young, massive open clusters. In addition, collisions lead mainly to the formation of a single runaway star instead of the formation of many high mass stars with a broad mass spectrum. We therefore conclude that massive stars form mainly by gas accretion, with stellar collisions only playing a minor role if any at all. Collisions of stars in the pre-main sequence phase might however contribute to the formation of the most massive stars in the densest star clusters and possibly to the formation of intermediate-mass black holes with masses up to a few 100 Msun.

preprint2010arXiv

Fragmentation in turbulent primordial gas

We report results from numerical simulations of star formation in the early universe that focus on the role of subsonic turbulence, and investigate whether it can induce fragmentation of the gas. We find that dense primordial gas is highly susceptible to fragmentation, even for rms turbulent velocity dispersions as low as 20% of the initial sound speed. The resulting fragments cover over two orders of magnitude in mass, ranging from 0.1 to 40 solar masses. However, our results suggest that the details of the fragmentation depend on the local properties of the turbulent velocity field and hence we expect considerable variations in the resulting stellar mass spectrum in different halos.

preprint2009arXiv

Modelling CO formation in the turbulent interstellar medium

We present results from high-resolution three-dimensional simulations of turbulent interstellar gas that self-consistently follow its coupled thermal, chemical and dynamical evolution, with a particular focus on the formation and destruction of H2 and CO. We quantify the formation timescales for H2 and CO in physical conditions corresponding to those found in nearby giant molecular clouds, and show that both species form rapidly, with chemical timescales that are comparable to the dynamical timescale of the gas. We also investigate the spatial distributions of H2 and CO, and how they relate to the underlying gas distribution. We show that H2 is a good tracer of the gas distribution, but that the relationship between CO abundance and gas density is more complex. The CO abundance is not well-correlated with either the gas number density n or the visual extinction A_V: both have a large influence on the CO abundance, but the inhomogeneous nature of the density field produced by the turbulence means that n and A_V are only poorly correlated. There is a large scatter in A_V, and hence CO abundance, for gas with any particular density, and similarly a large scatter in density and CO abundance for gas with any particular visual extinction. This will have important consequences for the interpretation of the CO emission observed from real molecular clouds. Finally, we also examine the temperature structure of the simulated gas. We show that the molecular gas is not isothermal. Most of it has a temperature in the range of 10--20 K, but there is also a significant fraction of warmer gas, located in low-extinction regions where photoelectric heating remains effective.

preprint2009arXiv

Solenoidal versus compressive turbulence forcing

We analyze the statistics and star formation rate obtained in high-resolution numerical experiments of forced supersonic turbulence, and compare with observations. We concentrate on a systematic comparison of solenoidal (divergence-free) and compressive (curl-free) forcing, which are two limiting cases of turbulence driving. Our results show that for the same RMS Mach number, compressive forcing produces a three times larger standard deviation of the density probability distribution. When self-gravity is included in the models, the star formation rate is more than one order of magnitude higher for compressive forcing than for solenoidal forcing.

preprint2009arXiv

Star Formation at Very Low Metallicity. V. The Greater Importance of Initial Conditions Compared to Metallicity Thresholds

The formation of the first stars out of metal-free gas appears to result in stars at least an order of magnitude more massive than in the present-day case. We here consider what controls the transition from a primordial to a modern initial mass function. It has been proposed that this occurs when effective metal line cooling occurs at a metallicity threshold of Z/Z_sun > 10^{-3.5}. We study the influence of low levels of metal enrichment on the cooling and collapse of initially ionized gas in small protogalactic halos using three-dimensional, smoothed particle hydrodynamics simulations with particle splitting. Our initial conditions represent protogalaxies forming within a previously ionized H ii region that has not yet had time to cool and recombine. These differ considerably from those used in simulations predicting a metallicity threshold, where the gas was initially cold and only partially ionized. In the centrally condensed potential that we study here, a wide variety of initial conditions for the gas yield a monolithic central collapse. Our models show no fragmentation during collapse to number densities as high as 10^5 cm^{-3}, for metallicities reaching as high as 10^{-1} Z_sun in one rotating case, far above the threshold suggested by previous work. Rotation allows for the formation of gravitationally stable gas disks over large fractions of the local Hubble time. Turbulence slows the growth of the central density slightly, but both spherically symmetric and turbulent initial conditions collapse and form a single sink particle. We therefore argue that fragmentation at moderate density depends on the initial conditions for star formation more than on the metal abundances present.

preprint2008arXiv

The Clustering Behavior of Pre-Main Sequence Stars in NGC 346 in the Small Magellanic Cloud

We present evidence that the star-forming region NGC 346/N66 in the Small Magellanic Cloud is the product of hierarchical star formation, probably from more than one star formation event. We investigate the spatial distribution and clustering behavior of the pre-main sequence (PMS) stellar population in the region, using data obtained with Hubble Space Telescope's Advanced Camera for Surveys. By applying the nearest neighbor and minimum spanning tree methods on the rich sample of PMS stars previously discovered in the region we identify ten individual PMS clusters in the area and quantify their structures. The clusters show a wide range of morphologies from hierarchical multi-peak configurations to centrally condensed clusters. However, only about 40 per cent of the PMS stars belong to the identified clusters. The central association NGC 346 is identified as the largest stellar concentration, which cannot be resolved into subclusters. Several PMS clusters are aligned along filaments of higher stellar density pointing away from the central part of the region. The PMS density peaks in the association coincide with the peaks of [OIII] and 8 micron emission. While more massive stars seem to be concentrated in the central association when considering the entire area, we find no evidence for mass segregation within the system itself.

preprint2004arXiv

Stellar Mass Spectra from Non-Isothermal Gravoturbulent Fragmentation

(abridged version) Identifying the processes that determine the initial mass function of stars (IMF) is a fundamental problem in star formation theory. One of the major uncertainties is the exact chemical state of the star forming gas and its influence on the dynamical evolution. Most simulations of star forming clusters use an isothermal equation of state (EOS). We address these issues and study the effect of a piecewise polytropic EOS on the formation of stellar clusters in turbulent, self-gravitating molecular clouds. We increase the polytropic exponent gamma from 0.7 to 1.1 at some chosen density n_c, which we vary from from 4.3x10^4 cm^-3 to 4.3x10^7 cm^-3. The change of thermodynamic state at n_c selects a characteristic mass scale for fragmentation M_ch, which we relate to the peak of the observed IMF. We find a relation M_ch ~ n_c^-0.5, supporting the idea that the distribution of stellar masses largely depends on the thermodynamic state of the star-forming gas.

preprint2000arXiv

Gravitational Collapse in Turbulent Molecular Clouds. II. Magnetohydrodynamical Turbulence

Hydrodynamic supersonic turbulence can only prevent local gravitational collapse if the turbulence is driven on scales smaller than the local Jeans lengths in the densest regions, a very severe requirement (Paper I). Magnetic fields have been suggested to support molecular clouds either magnetostatically or via magnetohydrodynamic (MHD) waves. Whereas the first mechanism would form sheet-like clouds, the second mechanism not only could exert a pressure onto the gas counteracting the gravitational forces, but could lead to a transfer of turbulent kinetic energy down to smaller spatial scales via MHD wave interactions. This turbulent magnetic cascade might provide sufficient energy at small scales to halt local collapse. We test this hypothesis with MHD simulations at resolutions up to 256^3 zones, done with ZEUS-3D. We first derive a resolution criterion for self-gravitating, magnetized gas: in order to prevent collapse of magnetostatically supported regions due to numerical diffusion, the minimum Jeans length must be resolved by four zones. Resolution of MHD waves increases this requirement to roughly six zones. We then find that magnetic fields cannot prevent local collapse unless they provide magnetostatic support. Weaker magnetic fields do somewhat delay collapse and cause it to occur more uniformly across the supported region in comparison to the hydrodynamical case. However, they still cannot prevent local collapse for much longer than a global free-fall time.