Showing 1–50 of 67 results for author: Peters, T

Gravitational collapse at low to moderate Mach numbers: The relationship between star formation efficiency and the fraction of mass in the massive object

Authors: Jorge Saavedra-Bastidas, Dominik R. G. Schleicher, Ralf S. Klessen, Sunmyon Chon, Kazuyuki Omukai, Thomas Peters, Lewis R. Prole, Bastián Reinoso, Rafeel Riaz, Paulo Solar

Abstract: The formation of massive objects via gravitational collapse is relevant both for explaining the origin of the first supermassive black holes and in the context of massive star formation. Here, we analyze simulations of the formation of massive objects pursued by different groups and in various environments, concerning the formation of supermassive black holes, primordial stars, as well as present-… ▽ More The formation of massive objects via gravitational collapse is relevant both for explaining the origin of the first supermassive black holes and in the context of massive star formation. Here, we analyze simulations of the formation of massive objects pursued by different groups and in various environments, concerning the formation of supermassive black holes, primordial stars, as well as present-day massive stars. We focus particularly on the regime of small virial parameters, i.e., low ratios of the initial kinetic to gravitational energy, low to moderate Mach numbers, and the phase before feedback is very efficient. We compare the outcomes of collapse under different conditions using dimensionless parameters, particularly the star formation efficiency ε_*, the fraction f_* of mass in the most massive object relative to the total stellar mass, and the fraction f_{\rm tot} of mass of the most massive object as a function of the total mass. We find that in all simulations analyzed here, f_{\rm tot} increases as a function of ε_*, although the steepness of the increase depends on the environment. The relation between f_* and ε_* is found to be more complex and also strongly depends on the number of protostars present at the beginning of the simulations. We show that a collision parameter, estimated as the ratio of the system size divided by the typical collision length, allows us to approximately characterize whether collisions will be important. We analyze the statistical correlation between the dimensionless quantities using the Spearman coefficient and confirm via a machine learning analysis that good predictions of f_* can be obtained from ε_* together with a rough estimate of the collision parameter. This suggests that a good estimate of the mass of the most massive object can be obtained once the maximum efficiency for a given environment is known. △ Less

Submitted 2 August, 2024; originally announced August 2024.

Comments: 13 pages, 8 figures, 2 tables. Submitted to A&A

Density distributions, magnetic field structures and fragmentation in high-mass star formation

Authors: H. Beuther, C. Gieser, J. D. Soler, Q. Zhang, R. Rao, D. Semenov, Th. Henning, R. Pudritz, T. Peters, P. Klaassen, M. T. Beltran, A. Palau, T. Moeller, K. G. Johnston, H. Zinnecker, J. Urquhart, R. Kuiper, A. Ahmadi, A. Sanchez-Monge, S. Feng, S. Leurini, S. E. Ragan

Abstract: Methods: Observing the large pc-scale Stokes I mm dust continuum emission with the IRAM 30m telescope and the intermediate-scale (<0.1pc) polarized submm dust emission with the Submillimeter Array toward a sample of 20 high-mass star-forming regions allows us to quantify the dependence of the fragmentation behaviour of these regions depending on the density and magnetic field structures. Results… ▽ More Methods: Observing the large pc-scale Stokes I mm dust continuum emission with the IRAM 30m telescope and the intermediate-scale (<0.1pc) polarized submm dust emission with the Submillimeter Array toward a sample of 20 high-mass star-forming regions allows us to quantify the dependence of the fragmentation behaviour of these regions depending on the density and magnetic field structures. Results: We infer density distributions n~r^{-p} of the regions with typical power-law slopes p around ~1.5. There is no obvious correlation between the power-law slopes of the density structures on larger clump scales (~1pc) and the number of fragments on smaller core scales (<0.1pc). Comparing the large-scale single-dish density profiles to those derived earlier from interferometric observations at smaller spatial scales, we find that the smaller-scale power-law slopes are steeper, typically around ~2.0. The flattening toward larger scales is consistent with the star-forming regions being embedded in larger cloud structures that do not decrease in density away from a particular core. Regarding the magnetic field, for several regions it appears aligned with filamentary structures leading toward the densest central cores. Furthermore, we find different polarization structures with some regions exhibiting central polarization holes whereas other regions show polarized emission also toward the central peak positions. Nevertheless, the polarized intensities are inversely related to the Stokes I intensities. We estimate magnetic field strengths between ~0.2 and ~4.5mG, and we find no clear correlation between magnetic field strength and the fragmentation level of the regions. Comparison of the turbulent to magnetic energies shows that they are of roughly equal importance in this sample. The mass-to-flux ratios range between ~2 and ~7, consistent with collapsing star-forming regions. △ Less

Submitted 20 November, 2023; originally announced November 2023.

Comments: Accepted for Astronomy & Astrophysics, 14 pages, 14 figures plus appendices, also download option at https://www2.mpia-hd.mpg.de/homes/beuther/papers.html

Feasibility of Passive Sounding of Uranian Moons using Uranian Kilometric Radiation

Authors: Andrew Romero-Wolf, Gregor Steinbruegge, Julie Castillo-Rogez, Corey J. Cochrane, Tom A. Nordheim, Karl L. Mitchell, Natalie S. Wolfenbarger, Dustin M. Schroeder, Sean T. Peters

Abstract: We present a feasibility study for passive sounding of Uranian icy moons using Uranian Kilometric Radio (UKR) emissions in the 100 - 900 kHz band. We provide a summary description of the observation geometry, the UKR characteristics, and estimate the sensitivity for an instrument analogous to the Cassini Radio Plasma Wave Science (RPWS) but with a modified receiver digitizer and signal processing… ▽ More We present a feasibility study for passive sounding of Uranian icy moons using Uranian Kilometric Radio (UKR) emissions in the 100 - 900 kHz band. We provide a summary description of the observation geometry, the UKR characteristics, and estimate the sensitivity for an instrument analogous to the Cassini Radio Plasma Wave Science (RPWS) but with a modified receiver digitizer and signal processing chain. We show that the concept has the potential to directly and unambiguously detect cold oceans within Uranian satellites and provide strong constraints on the interior structure in the presence of warm or no oceans. As part of a geophysical payload, the concept could therefore have a key role in the detection of oceans within the Uranian satellites. The main limitation of the concept is coherence losses attributed to the extended source size of the UKR and dependence on the illumination geometry. These factors represent constraints on the tour design of a future Uranus mission in terms of flyby altitudes and encounter timing. △ Less

Submitted 5 May, 2023; originally announced May 2023.

Kinematics and stability of high-mass protostellar disk candidates at sub-arcsecond resolution -- Insights from the IRAM NOEMA large program CORE

Authors: Aida Ahmadi, H. Beuther, F. Bosco, C. Gieser, S. Suri, J. C. Mottram, R. Kuiper, Th. Henning, Á. Sánchez-Monge, H. Linz, R. E. Pudritz, D. Semenov, J. M. Winters, T. Möller, M. T. Beltrán, T. Csengeri, R. Galván-Madrid, K. G. Johnston, E. Keto, P. D. Klaassen, S. Leurini, S. N. Longmore, S. L. Lumsden, L. T. Maud, L. Moscadelli , et al. (6 additional authors not shown)

Abstract: The fragmentation mode of high-mass molecular clumps and the accretion processes that form the most massive stars ($M\gtrsim 8M_\odot$) are still not well understood. To this end, we have undertaken a large observational program (CORE) making use of interferometric observations from the Northern Extended Millimetre Array (NOEMA) for a sample of 20 luminous ($L>10^4L_\odot$) protostellar objects in… ▽ More The fragmentation mode of high-mass molecular clumps and the accretion processes that form the most massive stars ($M\gtrsim 8M_\odot$) are still not well understood. To this end, we have undertaken a large observational program (CORE) making use of interferometric observations from the Northern Extended Millimetre Array (NOEMA) for a sample of 20 luminous ($L>10^4L_\odot$) protostellar objects in the 1.37 mm wavelength regime in both continuum and line emission, reaching $\sim$0.4" resolution (800 au at 2 kpc). Using the dense gas tracer CH$_3$CN, we find velocity gradients across 13 cores perpendicular to the directions of bipolar molecular outflows, making them excellent disk candidates. Specific angular momentum ($j$) radial profiles are on average $\sim10^{-3}$ km /s pc and follow $j \propto r^{1.7}$, consistent with a poorly resolved rotating and infalling envelope/disk model. Fitting the velocity profiles with a Keplerian model, we find protostellar masses in the range of $\sim 10-25$ $M_\odot$. Modelling the level population of CH$_3$CN lines, we present temperature maps and find median gas temperatures in the range $70-210$ K. We create Toomre $Q$ maps to study the stability of the disks and find almost all (11 of 13) disk candidates to be prone to fragmentation due to gravitational instabilities at the scales probed by our observations. In particular, disks with masses greater than $\sim10-20\%$ of the mass of their host (proto)stars are Toomre unstable, and more luminous protostellar objects tend to have disks that are more massive and hence more prone to fragmentation. Our finings show that most disks around high-mass protostars are prone to disk fragmentation early in their formation due to their high disk to stellar mass ratio. This impacts the accretion evolution of high-mass protostars which will have significant implications for the formation of the most massive stars. △ Less

Submitted 3 May, 2023; v1 submitted 28 April, 2023; originally announced May 2023.

Comments: 27 pages, 12 figures, 6 appendices - accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 677, A171 (2023)

Clustered Formation of Massive Stars within an Ionized Rotating Disk

Authors: Roberto Galván-Madrid, Qizhou Zhang, Andrés Izquierdo, Charles J. Law, Thomas Peters, Eric Keto, Hauyu Baobab Liu, Paul T. P. Ho, Adam Ginsburg, Carlos Carrasco-González

Abstract: We present ALMA observations with a 800 au resolution and radiative-transfer modelling of the inner part ($r\approx6000$ au) of the ionized accretion flow around a compact star cluster in formation at the center of the luminous ultra-compact (UC) HII region G10.6-0.4. We modeled the flow with an ionized Keplerian disk with and without radial motions in its outer part, or with an external Ulrich en… ▽ More We present ALMA observations with a 800 au resolution and radiative-transfer modelling of the inner part ($r\approx6000$ au) of the ionized accretion flow around a compact star cluster in formation at the center of the luminous ultra-compact (UC) HII region G10.6-0.4. We modeled the flow with an ionized Keplerian disk with and without radial motions in its outer part, or with an external Ulrich envelope. The MCMC fits to the data give total stellar masses $M_\star$ from 120 to $200~M_\odot$, with much smaller ionized-gas masses $M_\mathrm{ion-gas} = 0.2$ to $0.25~M_\odot$. The stellar mass is distributed within the gravitational radius $R_g\approx 1000$ to 1500 au, where the ionized gas is bound. The viewing inclination angle from the face-on orientation is $i = 49$ to $56~°$. Radial motions at radii $r > R_g$ converge to $v_{r,0} \approx 8.7$ km/s, or about the speed of sound of ionized gas, indicating that this gas is marginally unbound at most. From additional constraints on the ionizing-photon rate and far-IR luminosity of the region, we conclude that the stellar cluster consists of a few massive stars with $M_\mathrm{star} = 32$ to $60~M_\odot$, or one star in this range of masses accompanied by a population of lower-mass stars. Any active accretion of ionized gas onto the massive (proto)stars is residual. The inferred cluster density is very large, comparable to that reported at similar scales in the Galactic Center. Stellar interactions are likely to occur within the next Myr. △ Less

Submitted 11 December, 2022; originally announced December 2022.

Comments: Accepted to ApJ Letters, December 7 2022

The sharp ALMA view of infall and outflow in the massive protocluster G31.41+0.31

Authors: M. T. Beltrán, V. M. Rivilla, R. Cesaroni, D. Galli, L. Moscadelli, A. Ahmadi, H. Beuther, S. Etoka, C. Goddi, P. D. Klaassen, R. Kuiper, M. S. N. Kumar, A. Lorenzani, T. Peters, Á. Sánchez-Monge, P. Schilke, F. van der Tak, S. Vig

Abstract: Context. To better understand the formation of high-mass stars, it is fundamental to investigate how matter accretes onto young massive stars, how it is ejected, and how all this differs from the low-mass case. The massive protocluster G31.41+0.31 is the ideal target to study all these processes because observations at millimeter and centimeter wavelengths have resolved the emission of the Main co… ▽ More Context. To better understand the formation of high-mass stars, it is fundamental to investigate how matter accretes onto young massive stars, how it is ejected, and how all this differs from the low-mass case. The massive protocluster G31.41+0.31 is the ideal target to study all these processes because observations at millimeter and centimeter wavelengths have resolved the emission of the Main core into at least four massive dust continuum sources, named A, B, C, and D, within 1" or 0.018 pc, and have identified signatures of infall and several outflows associated with the core. Aims. We study the interplay between infall and outflow in G31.41+0.31 by investigating at a spatial resolution of a few 100 au their properties and their possible impact on the core. Methods. We carried out molecular line observations of typical high-density tracers, such as CH3CN or H2CO, and shock and outflow tracers, such as SiO, with ALMA at 1.4 mm that achieved an angular resolution of 0.09" (340 au). Results. The observations have revealed inverse P-Cygni profiles in CH3CN and H2CO toward the four sources embedded in the Main core, suggesting that all of them are undergoing collapse. The infall rates, estimated from the red-shifted absorption are on the order of 1E-2 Msun/yr. The individual infall rates imply that the accretion timescale of the Main core is an order of magnitude smaller than its rotation timescale. This confirms that rotating toroids such as the G31 Main core are non-equilibrium, transient collapsing structures that need to be constantly replenished with fresh material from a large-scale reservoir. For sources B, C, and D, the infall could be accelerating inside the sources, while for source A, the presence of a second emission component complicates the interpretation. The SiO observations have revealed the presence of at least six outflows in the G31.41+0.31 star-forming region, ... △ Less

Submitted 25 January, 2022; originally announced January 2022.

Comments: 17 pages, 9 figures, 2 tables. Accepted by A&A

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

Clustered star formation at early evolutionary stages. Physical and chemical analysis of the young star-forming regions ISOSS J22478+6357 and ISOSS J23053+5953

Authors: C. Gieser, H. Beuther, D. Semenov, S. Suri, J. D. Soler, H. Linz, J. Syed, Th. Henning, S. Feng, T. Möller, A. Palau, J. M. Winters, M. T. Beltrán, R. Kuiper, L. Moscadelli, P. Klaassen, J. S. Urquhart, T. Peters, S. N. Longmore, Á. Sánchez-Monge, R. Galván-Madrid, R. E. Pudritz, K. G. Johnston

Abstract: We aim to characterize the physical and chemical properties of fragmented cores during the earliest evolutionary stages in the very young star-forming regions ISOSS J22478+6357 and ISOSS J23053+5953. NOEMA 1.3 mm data are used in combination with archival mid- and far-infrared observations to construct and fit the SEDs of individual fragmented cores. The radial density profiles are inferred from t… ▽ More We aim to characterize the physical and chemical properties of fragmented cores during the earliest evolutionary stages in the very young star-forming regions ISOSS J22478+6357 and ISOSS J23053+5953. NOEMA 1.3 mm data are used in combination with archival mid- and far-infrared observations to construct and fit the SEDs of individual fragmented cores. The radial density profiles are inferred from the 1.3 mm continuum visibility profiles and the radial temperature profiles are estimated from H2CO rotation temperature maps. Molecular column densities are derived with the line fitting tool XCLASS. The physical and chemical properties are combined by applying the physical-chemical model MUSCLE in order to constrain the chemical timescales of a few line-rich cores. The morphology and spatial correlations of the molecular emission are analyzed using the HOG method. The mid-infrared data show that both regions contain a cluster of young stellar objects. Bipolar molecular outflows are observed in the CO 2-1 transition toward the strong mm cores indicating protostellar activity. We find strong molecular emission of SO, SiO, H2CO, and CH3OH in locations which are not associated with the mm cores. These shocked knots can be either associated with the bipolar outflows or, in the case of ISOSS J23053+5953, with a colliding flow that creates a large shocked region between the mm cores. The mean chemical timescale of the cores is lower (20 000 yr) compared to that of the sources of the more evolved CORE sample (60 000 yr). With the HOG method, we find that the spatial emission of species tracing the extended emission and of shock-tracing molecules are well correlated within transitions of these groups. △ Less

Submitted 14 October, 2021; v1 submitted 5 October, 2021; originally announced October 2021.

Comments: 34 pages, 17 figures, accepted for publication in A&A

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

Disk fragmentation in high-mass star formation. High-resolution observations towards AFGL 2591-VLA 3

Authors: S. Suri, H. Beuther, C. Gieser, A. Ahmadi, Á. Sánchez-Monge, J. M. Winters, H. Linz, Th. Henning, M. T. Beltrán, F. Bosco, R. Cesaroni, T. Csengeri, S. Feng, M. G. Hoare, K. G. Johnston, P. Klaasen, R. Kuiper, S. Leurini, S. Longmore, S. Lumsden, L. Maud, L. Moscadelli, T. Möller, A. Palau, T. Peters , et al. (7 additional authors not shown)

Abstract: Increasing evidence suggests that, similar to their low-mass counterparts, high-mass stars form through a disk-mediated accretion process. At the same time, formation of high-mass stars still necessitates high accretion rates, and hence, high gas densities, which in turn can cause disks to become unstable against gravitational fragmentation. We study the kinematics and fragmentation of the disk ar… ▽ More Increasing evidence suggests that, similar to their low-mass counterparts, high-mass stars form through a disk-mediated accretion process. At the same time, formation of high-mass stars still necessitates high accretion rates, and hence, high gas densities, which in turn can cause disks to become unstable against gravitational fragmentation. We study the kinematics and fragmentation of the disk around the high-mass star forming region AFGL 2591-VLA 3 which was hypothesized to be fragmenting based on the observations that show multiple outflow directions. We use a new set of high-resolution (0.19 arcsec) IRAM/NOEMA observations at 843 micron towards VLA 3 which allow us to resolve its disk, characterize the fragmentation, and study its kinematics. In addition to the 843 micron continuum emission, our spectral setup targets warm dense gas and outflow tracers such as HCN, HC$_3$N and SO$_2$, as well as vibrationally excited HCN lines. The high resolution continuum and line emission maps reveal multiple fragments with subsolar masses within the inner 1000 AU of VLA 3. Furthermore, the velocity field of the inner disk observed at 843 micron shows a similar behavior to that of the larger scale velocity field studied in the CORE project at 1.37 mm. We present the first observational evidence for disk fragmentation towards AFGL 2591-VLA 3, a source that was thought to be a single high-mass core. While the fragments themselves are low-mass, the rotation of the disk is dominated by the protostar with a mass of 10.3$\pm 1.8~M_{\odot}$. These data also show that NOEMA Band 4 can obtain the highest currently achievable spatial resolution at (sub-)mm wavelengths in observations of strong northern sources. △ Less

Submitted 10 September, 2021; originally announced September 2021.

Comments: 16 pages, 11 figures, Accepted for publication in A&A

Journal ref: A&A 655, A84 (2021)

Fragmentation and kinematics in high-mass star formation: CORE-extension targeting two very young high-mass star-forming regions

Authors: H. Beuther, C. Gieser, S. Suri, H. Linz, P. Klaassen, D. Semenov, J. M. Winters, Th. Henning, J. D. Soler, J. S. Urquhart, J. Syed, S . Feng, T. Moeller, M. T. Beltran, A. Sanchez-Monge, S. N. Longmore, T. Peters, J. Ballesteros-Paredes, P. Schilke, L. Moscadelli, A. Palau, R. Cesaroni, S. Lumsden, R. Pudritz, F. Wyrowski , et al. (2 additional authors not shown)

Abstract: Context: The formation of high-mass star-forming regions from their parental gas cloud and the subsequent fragmentation processes lie at the heart of star formation research. Aims: We aim to study the dynamical and fragmentation properties at very early evolutionary stages of high-mass star formation. Methods: Employing the NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30m telescope, we… ▽ More Context: The formation of high-mass star-forming regions from their parental gas cloud and the subsequent fragmentation processes lie at the heart of star formation research. Aims: We aim to study the dynamical and fragmentation properties at very early evolutionary stages of high-mass star formation. Methods: Employing the NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30m telescope, we observed two young high-mass star-forming regions, ISOSS22478 and ISOSS23053, in the 1.3mm continuum and spectral line emission at a high angular resolution (~0.8''). Results: We resolved 29 cores that are mostly located along filament-like structures. Depending on the temperature assumption, these cores follow a mass-size relation of approximately M~r^2.0, corresponding to constant mean column densities. However, with different temperature assumptions, a steeper mass-size relation up to M~r^3.0, which would be more likely to correspond to constant mean volume densities, cannot be ruled out. The correlation of the core masses with their nearest neighbor separations is consistent with thermal Jeans fragmentation. We found hardly any core separations at the spatial resolution limit, indicating that the data resolve the large-scale fragmentation well. Although the kinematics of the two regions appear very different at first sight - multiple velocity components along filaments in ISOSS22478 versus a steep velocity gradient of more than 50km/s/pc in ISOSS23053 - the findings can be explained within the framework of a dynamical cloud collapse scenario. Conclusions: While our data are consistent with a dynamical cloud collapse scenario and subsequent thermal Jeans fragmentation, the importance of additional environmental properties, such as the magnetization of the gas or external shocks triggering converging gas flows, is nonetheless not as well constrained and would require future investigation. △ Less

Submitted 6 April, 2021; originally announced April 2021.

Comments: 16 pages, 16 figures, accepted for Astronomy and Astrophysics, a higher-resolution version can also be found at https://www.mpia.de/homes/beuther/papers.html

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

Fragmentation in the massive G31.41+0.31 protocluster

Authors: M. T. Beltrán, V. M. Rivilla, R. Cesaroni, L. T. Maud, D. Galli, L. Moscadelli, A. Lorenzani, A. Ahmadi, H. Beuther, T. Csengeri, S. Etoka, C. Goddi, P. D. Klaassen, R. Kuiper, M. S. N. Kumar, T. Peters, Á. Sánchez-Monge, P. Schilke, F. van der Tak, S. Vig, H. Zinnecker

Abstract: Context. ALMA observations at 1.4 mm and 0.2'' (750au) angular resolution of the Main core in the high-mass star forming region G31.41+0.31 have revealed a puzzling scenario: on the one hand, the continuum emission looks very homogeneous and the core appears to undergo solid-body rotation, suggesting a monolithic core stabilized by the magnetic field; on the other hand, rotation and infall speed u… ▽ More Context. ALMA observations at 1.4 mm and 0.2'' (750au) angular resolution of the Main core in the high-mass star forming region G31.41+0.31 have revealed a puzzling scenario: on the one hand, the continuum emission looks very homogeneous and the core appears to undergo solid-body rotation, suggesting a monolithic core stabilized by the magnetic field; on the other hand, rotation and infall speed up toward the core center, where two massive embedded free-free continuum sources have been detected, pointing to an unstable core having undergone fragmentation. Aims. To establish whether the Main core is indeed monolithic or its homogeneous appearance is due to a combination of large dust opacity and low angular resolution, we carried out millimeter observations at higher angular resolution and different wavelengths. Methods. We carried out ALMA observations at 1.4 mm and 3.5 mm that achieved angular resolutions of 0.1''(375 au) and 0.075'' (280 au), respectively. VLA observations at 7 mm and 1.3 cm at even higher angular resolution, 0.05'' (190 au) and 0.07'' (260 au), respectively, were also carried out to better study the nature of the free-free continuum sources detected in the core. Results. The millimeter continuum emission of the Main core has been clearly resolved into at least four sources, A, B, C, and D, within 1'', indicating that the core is not monolithic. The deconvolved radii of the dust emission of the sources, estimated at 3.5 mm, are 400-500au, their masses range from 15 to 26 Msun, and their number densities are several 1E9 cm-3. Sources A and B, located closer to the center of the core and separated by 750 au, are clearly associated with two free-free continuum sources, likely thermal radio jets, and are the brightest in the core. The spectral energy distribution of these two sources and their masses and sizes are similar and suggest a common origin. △ Less

Submitted 8 March, 2021; originally announced March 2021.

Comments: 18 pages,9 figures, 5 tables. Accepted by A&A

Journal ref: A&A 648, A100 (2021)

The physical and chemical structure of high-mass star-forming regions. Unraveling chemical complexity with the NOEMA large program "CORE"

Authors: C. Gieser, H. Beuther, D. Semenov, A. Ahmadi, S. Suri, T. Möller, M. T. Beltran, P. Klaassen, Q. Zhang, J. S. Urquhart, Th. Henning, S. Feng, R. Galván-Madrid, V. de Souza Magalhães, L. Moscadelli, S. Longmore, S. Leurini, R. Kuiper, T. Peters, K. M. Menten, T. Csengeri, G. Fuller, F. Wyrowski, S. Lumsden, Á. Sánchez-Monge , et al. (8 additional authors not shown)

Abstract: We use sub-arcsecond resolution ($\sim$0.4$''$) observations with NOEMA at 1.37 mm to study the dust emission and molecular gas of 18 high-mass star-forming regions. We combine the derived physical and chemical properties of individual cores in these regions to estimate their ages. The temperature structure of these regions are determined by fitting H2CO and CH3CN line emission. The density profil… ▽ More We use sub-arcsecond resolution ($\sim$0.4$''$) observations with NOEMA at 1.37 mm to study the dust emission and molecular gas of 18 high-mass star-forming regions. We combine the derived physical and chemical properties of individual cores in these regions to estimate their ages. The temperature structure of these regions are determined by fitting H2CO and CH3CN line emission. The density profiles are inferred from the 1.37 mm continuum visibilities. The column densities of 11 different species are determined by fitting the emission lines with XCLASS. Within the 18 observed regions, we identify 22 individual cores with associated 1.37 mm continuum emission and with a radially decreasing temperature profile. We find an average temperature power-law index of q = 0.4$\pm$0.1 and an average density power-law index of p = 2.0$\pm$0.2 on scales on the order of several 1 000 au. Comparing these results with values of p derived in the literature suggest that the density profiles remain unchanged from clump to core scales. The column densities relative to N(C18O) between pairs of dense gas tracers show tight correlations. We apply the physical-chemical model MUSCLE to the derived column densities of each core and find a mean chemical age of $\sim$60 000 yrs and an age spread of 20 000-100 000 yrs. With this paper we release all data products of the CORE project available at https://www.mpia.de/core. The CORE sample reveals well constrained density and temperature power-law distributions. Furthermore, we characterize a large variety in molecular richness that can be explained by an age spread confirmed by our physical-chemical modeling. The hot molecular cores show the most emission lines, but we also find evolved cores at an evolutionary stage, in which most molecules are destroyed and thus the spectra appear line-poor again. △ Less

Submitted 23 February, 2021; originally announced February 2021.

Comments: 29 pages, 14 figures, accepted for publication in A&A

Journal ref: A&A 648, A66 (2021)

Multi-scale view of star formation in IRAS 21078+5211: From clump fragmentation to disk wind

Authors: L. Moscadelli, H. Beuther, A. Ahmadi, C. Gieser, F. Massi, R. Cesaroni, Á. Sánchez-Monge, F. Bacciotti, M. T. Beltrán, T. Csengeri, R. Galván-Madrid, Th. Henning, P. D. Klaassen, R. Kuiper, S. Leurini, S. N. Longmore, L. T. Maud, T. Möller, A. Palau, T. Peters, R. E. Pudritz, A. Sanna, D. Semenov, J. S. Urquhart, J. M. Winters , et al. (1 additional authors not shown)

Abstract: In the massive star-forming region IRAS 21078+5211, a highly fragmented cluster (0.1~pc in size) of molecular cores is observed, located at the density peak of an elongated (1~pc in size) molecular cloud. A small (1~km/s per 0.1~pc) LSR velocity (Vlsr) gradient is detected across the axis of the molecular cloud. Assuming we are observing a mass flow from the harboring cloud to the cluster, we deri… ▽ More In the massive star-forming region IRAS 21078+5211, a highly fragmented cluster (0.1~pc in size) of molecular cores is observed, located at the density peak of an elongated (1~pc in size) molecular cloud. A small (1~km/s per 0.1~pc) LSR velocity (Vlsr) gradient is detected across the axis of the molecular cloud. Assuming we are observing a mass flow from the harboring cloud to the cluster, we derive a mass infall rate of about 10^{-4}~M_{sun}~yr^{-1}. The most massive cores (labeled 1, 2, and 3) are found at the center of the cluster, and these are the only ones that present a signature of protostellar activity in terms of emission from high-excitation molecular lines or a molecular outflow. We reveal an extended (size about 0.1~pc), bipolar collimated molecular outflow emerging from core 1. We believe this is powered by a (previously discovered) compact (size <= 1000~au) radio jet, ejected by a YSO embedded in core 1 (named YSO-1), since the molecular outflow and the radio jet are almost parallel and have a comparable momentum rate. By means of high-excitation lines, we find a large (14~km/s over 500~au) Vlsr gradient at the position of YSO-1, oriented approximately perpendicular to the radio jet. Assuming this is an edge-on, rotating disk and fitting a Keplerian rotation pattern, we determine the YSO-1 mass to be 5.6+/-2.0~M_{sun}. The water masers (previously observed with VLBI) emerge within 100-300~au from YSO-1 and are unique tracers of the jet kinematics. Their three-dimensional (3D) velocity pattern reveals that the gas flows along, and rotates about, the jet axis. We show that the 3D maser velocities are fully consistent with the magneto-centrifugal disk-wind models predicting a cylindrical rotating jet. Under this hypothesis, we determine the jet radius to be about 16~au and the corresponding launching radius and terminal velocity to be about 2.2~au and 200~km/s, respectively. △ Less

Submitted 9 February, 2021; originally announced February 2021.

Comments: 25 pages, 16 figures

Journal ref: A&A 647, A114 (2021)

Time-Variable Radio Recombination Line Emission in W49A

Authors: C. G. De Pree, D. J. Wilner, L. E. Kristensen, R. Galván-Madrid, W. M. Goss, R. S. Klessen, M. -M. Mac Low, T. Peters, A. Robinson, S. Sloman, M. Rao

Abstract: We present new Jansky Very Large Array (VLA) images of the central region of the W49A star-forming region at 3.6~cm and at 7~mm at resolutions of 0\farcs15 (1650 au) and 0\farcs04 (440 au), respectively. The 3.6~cm data reveal new morphological detail in the ultracompact \ion{H}{2} region population, as well as several previously unknown and unresolved sources. In particular, source A shows elonga… ▽ More We present new Jansky Very Large Array (VLA) images of the central region of the W49A star-forming region at 3.6~cm and at 7~mm at resolutions of 0\farcs15 (1650 au) and 0\farcs04 (440 au), respectively. The 3.6~cm data reveal new morphological detail in the ultracompact \ion{H}{2} region population, as well as several previously unknown and unresolved sources. In particular, source A shows elongated, edge-brightened, bipolar lobes, indicative of a collimated outflow, and source E is resolved into three spherical components. We also present VLA observations of radio recombination lines at 3.6~cm and 7~mm, and IRAM Northern Extended Millimeter Array (NOEMA) observations at 1.2~mm. Three of the smallest ultracompact \ion{H}{2} regions (sources A, B2 and G2) all show broad kinematic linewidths, with $Δ$V$_{FWHM}\gtrsim$40~km~s$^{-1}$. A multi-line analysis indicates that broad linewidths remain after correcting for pressure broadening effects, suggesting the presence of supersonic flows. Substantial changes in linewidth over the 21 year time baseline at both 3.6 cm and 7 mm are found for source G2. At 3.6 cm, the linewidth of G2 changed from 31.7$\pm$1.8 km s$^{-1}$ to 55.6$\pm$2.7 km s$^{-1}$, an increase of $+$23.9$\pm$3.4 km s$^{-1}$. The G2 source was previously reported to have shown a 3.6~cm continuum flux density decrease of 40\% between 1994 and 2015. This source sits near the center of a very young bipolar outflow whose variability may have produced these changes. △ Less

Submitted 9 October, 2020; originally announced October 2020.

Comments: 18 pages, 7 figures, Accepted for publication in The Astronomical Journal

Multi-wavelength modelling of the circumstellar environment of the massive proto-star AFGL 2591 VLA 3

Authors: F. A. Olguin, M. G. Hoare, K. G. Johnston, F. Motte, H. -R. V. Chen, H. Beuther, J. C. Mottram, A. Ahmadi, C. Gieser, D. Semenov, T. Peters, A. Palau, P. D. Klaassen, R. Kuiper, Á. Sánchez-Monge, Th. Henning

Abstract: We have studied the dust density, temperature and velocity distributions of the archetypal massive young stellar object (MYSO) AFGL 2591. Given its high luminosity ($L=2 \times 10^5$ L$_\odot$) and distance ($d=3.3$ kpc), AFGL 2591 has one of the highest $\sqrt{L}/d$ ratio, giving better resolved dust emission than any other MYSO. As such, this paper provides a template on how to use resolved mult… ▽ More We have studied the dust density, temperature and velocity distributions of the archetypal massive young stellar object (MYSO) AFGL 2591. Given its high luminosity ($L=2 \times 10^5$ L$_\odot$) and distance ($d=3.3$ kpc), AFGL 2591 has one of the highest $\sqrt{L}/d$ ratio, giving better resolved dust emission than any other MYSO. As such, this paper provides a template on how to use resolved multi-wavelength data and radiative transfer to obtain a well-constrained 2-D axi-symmetric analytic rotating infall model. We show for the first time that the resolved dust continuum emission from Herschel 70 $μ$m observations is extended along the outflow direction, whose origin is explained in part from warm dust in the outflow cavity walls. However, the model can only explain the kinematic features from CH$_3$CN observations with unrealistically low stellar masses ($<15$ M$_\odot$), indicating that additional physical processes may be playing a role in slowing down the envelope rotation. As part of our 3-step continuum and line fitting, we have identified model parameters that can be further constrained by specific observations. High-resolution mm visibilities were fitted to obtain the disc mass (6 M$_\odot$) and radius (2200 au). A combination of SED and near-IR observations were used to estimate the luminosity and envelope mass together with the outflow cavity inclination and opening angles. △ Less

Submitted 12 August, 2020; originally announced August 2020.

Comments: 34 pages, 22 figures, 15 tables, 4 appendices. Accepted for publication in MNRAS

Chemical complexity in high-mass star formation: An observational and modeling case study of the AFGL 2591 VLA 3 hot core

Authors: C. Gieser, D. Semenov, H. Beuther, A. Ahmadi, J. C. Mottram, Th. Henning, M. Beltran, L. T. Maud, F. Bosco, S. Leurini, T. Peters, P. Klaassen, R. Kuiper, S. Feng, J. S. Urquhart, L. Moscadelli, T. Csengeri, S. Lumsden, J. M. Winters, S. Suri, Q. Zhang, R. Pudritz, A. Palau, K. M. Menten, R. Galvan-Madrid , et al. (8 additional authors not shown)

Abstract: We present a detailed observational and modeling study of the hot core VLA 3 in the high-mass star-forming region AFGL 2591, which is a target region of the NOrthern Extended Millimeter Array (NOEMA) large program CORE. Using NOEMA observations at 1.37 mm with an angular resolution of ~0."42 (1 400 au at 3.33 kpc), we derived the physical and chemical structure of the source. We modeled the observ… ▽ More We present a detailed observational and modeling study of the hot core VLA 3 in the high-mass star-forming region AFGL 2591, which is a target region of the NOrthern Extended Millimeter Array (NOEMA) large program CORE. Using NOEMA observations at 1.37 mm with an angular resolution of ~0."42 (1 400 au at 3.33 kpc), we derived the physical and chemical structure of the source. We modeled the observed molecular abundances with the chemical evolution code MUSCLE (MUlti Stage ChemicaL codE). Results. With the kinetic temperature tracers CH3CN and H2CO we observe a temperature distribution with a power-law index of q = 0.41+-0.08. Using the visibilities of the continuum emission we derive a density structure with a power-law index of p = 1.7+-0.1. The hot core spectra reveal high molecular abundances and a rich diversity in complex molecules. The majority of the molecules have an asymmetric spatial distribution around the forming protostar(s), which indicates a complex physical structure on scales < 1 400 au. Using MUSCLE, we are able to explain the observed molecular abundance of 10 out of 14 modeled species at an estimated hot core chemical age of ~21 100 years. In contrast to the observational analysis, our chemical modeling predicts a lower density power-law index of p < 1.4. Reasons for this discrepancy are discussed. Conclusions. Combining high spatial resolution observations with detailed chemical modeling allows us to derive a concise picture of the physical and chemical structure of the famous AFGL 2591 hot core. The next steps are to conduct a similar analysis for the whole CORE sample, and then use this analysis to constrain the chemical diversity in high-mass star formation to a much greater depth. △ Less

Submitted 11 October, 2019; originally announced October 2019.

Comments: 22 pages, 13 figures, accepted for publication in A&A

Journal ref: A&A 631, A142 (2019)

Fragmentation, rotation and outflows in the high-mass star-forming region IRAS 23033+5951. A case study of the IRAM NOEMA large program CORE

Authors: F. Bosco, H. Beuther, A. Ahmadi, J. C. Mottram, R. Kuiper, H. Linz, L. Maud, J. M. Winters, T. Henning, S. Feng, T. Peters, D. Semenov, P. D. Klaassen, P. Schilke, J. S. Urquhart, M. T. Beltrán, S. L. Lumsden, S. Leurini, L. Moscadelli, R. Cesaroni, Á. Sánchez-Monge, A. Palau, R. Pudritz, F. Wyrowski, S. Longmore , et al. (1 additional authors not shown)

Abstract: The formation process of high-mass stars (>8M$_\odot$) is poorly constrained, particularly, the effects of clump fragmentation creating multiple systems and the mechanism of mass accretion onto the cores. We study the fragmentation of dense gas clumps, and trace the circumstellar rotation and outflows by analyzing observations of the high-mass (~500M$_\odot$) star-forming region IRAS 23033+5951. U… ▽ More The formation process of high-mass stars (>8M$_\odot$) is poorly constrained, particularly, the effects of clump fragmentation creating multiple systems and the mechanism of mass accretion onto the cores. We study the fragmentation of dense gas clumps, and trace the circumstellar rotation and outflows by analyzing observations of the high-mass (~500M$_\odot$) star-forming region IRAS 23033+5951. Using the Northern Extended Millimeter Array (NOEMA) in three configurations and the IRAM 30-m single-dish telescope at 220GHz, we probe the gas and dust emission at an angular resolution of ~0.45arcsec, corresponding to 1900au. In the mm continuum emission, we identify a protostellar cluster with at least four mm-sources, where three of them show a significantly higher peak intensity well above a signal-to-noise ratio of 100. Hierarchical fragmentation from large to small spatial scales is discussed. Two fragments are embedded in rotating structures and drive molecular outflows, traced by $^{13}$CO (2-1) emission. The velocity profiles across two of the cores are similar to Keplerian but are missing the highest velocity components close to the center of rotation, which is a common phenomena from observations like these, and other rotation scenarios are not excluded entirely. Position-velocity diagrams suggest protostellar masses of ~6 and 19M$_\sun$. Rotational temperatures from fitting CH$_3$CN ($12_K-11_K$) spectra are used for estimating the gas temperature and by that the disk stability against gravitational fragmentation, utilizing Toomre's $Q$ parameter. [We] identify only one candidate disk to be unstable against gravitational instability caused by axisymmetric perturbations. The dominant sources cover different evolutionary stages within the same maternal gas clump. The appearance of rotation and outflows of the cores are similar to those found in low-mass star-forming regions. △ Less

Submitted 9 July, 2019; originally announced July 2019.

Comments: 23 pages, 16 figures, 5 tables, accepted for publication in A&A

Journal ref: A&A 629, A10 (2019)

Substructures in the Keplerian disc around the O-type (proto)star G17.64+0.16

Authors: L. T. Maud, R. Cesaroni, M. S. N. Kumar, V. M. Rivilla, A. Ginsburg, P. D. Klaassen, D. Harsono, A. Sanchez-Monge, A. Ahmadi, V. Allen, M. T. Beltran, H. Beuther, R. Galvan-Madrid, C. Goddi, M. G. Hoare, M. R. Hogerheijde, K. G. Johnston, R. Kuiper, L. Moscadelli, T. Peters, L. Testi, F. F. S. van der Tak, W. J. de Wit

Abstract: We present the highest angular resolution (20x15mas - 44x33au) Atacama Large Millimeter/sub-millimeter Array (ALMA) observations currently possible of the proto-O-star G17.64+0.16 in Band 6. The Cycle 5 observations with baselines out to 16km probes scales <50au and reveal the rotating disc around G17.64+0.16, a massive forming O-type star. The disc has a ring-like enhancement in the dust emission… ▽ More We present the highest angular resolution (20x15mas - 44x33au) Atacama Large Millimeter/sub-millimeter Array (ALMA) observations currently possible of the proto-O-star G17.64+0.16 in Band 6. The Cycle 5 observations with baselines out to 16km probes scales <50au and reveal the rotating disc around G17.64+0.16, a massive forming O-type star. The disc has a ring-like enhancement in the dust emission, especially visible as arc structures to the north and south. The Keplerian kinematics are most prominently seen in the vibrationally excited water line, H2O (Eu=3461.9K). The mass of the central source found by modelling the Keplerian rotation is consistent with 45+/-10Mo. The H30alpha (231.9GHz) radio-recombination line and the SiO (5-4) molecular line were detected at up to the 10 sigma$ level. The estimated disc mass is 0.6-2.6Mo under the optically thin assumption. Analysis of the Toomre Q parameter, in the optically thin regime, indicates that the disc stability is highly dependent on temperature. The disc currently appears stable for temperatures >150K, this does not preclude that the substructures formed earlier through disc fragmentation. △ Less

Submitted 15 June, 2019; originally announced June 2019.

Comments: 9 Total, inc references, appendix, long author list, 4 figures (inc appendix)

Journal ref: A&A 627, L6 (2019)

IRAS23385+6053: An embedded massive cluster in the making

Authors: R. Cesaroni, H. Beuther, A. Ahmadi, M. T. Beltran, T. Csengeri, R. Galvan-Madrid, C. Gieser, T. Henning, K. G. Johnston, P. D. Klaassen, R. Kuiper, S. Leurini, H. Linz, S. Longmore, S. L. Lumsden, L. T. Maud, L. Moscadelli, J. C. Mottram, A. Palau, T. Peters, R. E. Pudritz, A. Sanchez-Monge, P. Schilke, D. Semenov, S. Suri , et al. (4 additional authors not shown)

Abstract: This study is part of the project ``CORE'', an IRAM/NOEMA large program consisting of observations of the millimeter continuum and molecular line emission towards 20 selected high-mass star forming regions. We focus on IRAS23385+6053, which is believed to be the least evolved source of the CORE sample. The observations were performed at ~1.4 mm and employed three configurations of NOEMA and additi… ▽ More This study is part of the project ``CORE'', an IRAM/NOEMA large program consisting of observations of the millimeter continuum and molecular line emission towards 20 selected high-mass star forming regions. We focus on IRAS23385+6053, which is believed to be the least evolved source of the CORE sample. The observations were performed at ~1.4 mm and employed three configurations of NOEMA and additional single-dish maps, merged with the interferometric data to recover the extended emission. Our correlator setup covered a number of lines from well-known hot core tracers and a few outflow tracers. The angular (~0.45"$-$0.9") and spectral (0.5 km/s) resolutions were sufficient to resolve the clump in IRAS23385+6053 and investigate the existence of large-scale motions due to rotation, infall, or expansion. We find that the clump splits into six distinct cores when observed at sub-arcsecond resolution. These are identified through their 1.4 mm continuum and molecular line emission. We produce maps of the velocity, line width, and rotational temperature from the methanol and methyl cyanide lines, which allow us to investigate the cores and reveal a velocity and temperature gradient in the most massive core. We also find evidence of a bipolar outflow, possibly powered by a low-mass star. We present the tentative detection of a circumstellar self-gravitating disk lying in the most massive core and powering a large-scale outflow previously known in the literature. In our scenario, the star powering the flow is responsible for most of the luminosity of IRAS23385+6053 (~$3000~L_\odot$). The other cores, albeit with masses below the corresponding virial masses, appear to be accreting material from their molecular surroundings and are possibly collapsing or on the verge of collapse. We conclude that we are observing a sample of star-forming cores that is bound to turn into a cluster of massive stars. △ Less

Submitted 27 May, 2019; originally announced May 2019.

Journal ref: A&A 627, A68 (2019)

Young massive star cluster formation in the Galactic Centre is driven by global gravitational collapse of high-mass molecular clouds

Authors: A. T. Barnes, S. N. Longmore, A. Avison, Y. Contreras, A. Ginsburg, J. D. Henshaw, J. M. Rathborne, D. L. Walker, J. Alves, J. Bally, C. Battersby, M. T. Beltrán, H. Beuther, G. Garay, L. Gomez, J. Jackson, J. Kainulainen, J. M. D. Kruijssen, X. Lu, E. A. C. Mills, J. Ott, T. Peters

Abstract: Young massive clusters (YMCs) are the most compact, high-mass stellar systems still forming at the present day. The precursor clouds to such systems are, however, rare due to their large initial gas mass reservoirs and rapid dispersal timescales due to stellar feedback. Nonetheless, unlike their high-z counterparts, these precursors are resolvable down to the sites of individually forming stars, a… ▽ More Young massive clusters (YMCs) are the most compact, high-mass stellar systems still forming at the present day. The precursor clouds to such systems are, however, rare due to their large initial gas mass reservoirs and rapid dispersal timescales due to stellar feedback. Nonetheless, unlike their high-z counterparts, these precursors are resolvable down to the sites of individually forming stars, and hence represent the ideal environments in which to test the current theories of star and cluster formation. Using high angular resolution (1$^{\prime\prime}$ / 0.05pc) and sensitivity ALMA observations of two YMC progenitor clouds in the Galactic Centre, we have identified a suite of molecular line transitions -- e.g. c-C$_{3}$H$_{2} $($7-6$) -- that are believed to be optically thin, and reliably trace the gas structure in the highest density gas on star-forming core scales. We conduct a virial analysis of the identified core and proto-cluster regions, and show that half of the cores (5/10) and both proto-clusters are unstable to gravitational collapse. This is the first kinematic evidence of global gravitational collapse in YMC precursor clouds at such an early evolutionary stage. The implications are that if these clouds are to form YMCs, then they likely do so via the "conveyor-belt" mode, whereby stars continually form within dispersed dense gas cores as the cloud undergoes global gravitational collapse. The concurrent contraction of both the cluster-scale gas and embedded (proto)stars ultimately leads to the high (proto)stellar density in YMCs. △ Less

Submitted 14 March, 2019; originally announced March 2019.

Comments: Accepted for publication in MNRAS. 22 pages (+22 pages of appendix), 10 (+21 in appendix) figures, 5 (+1 in appendix) tables

Chasing discs around O-type (proto)stars - ALMA evidence for an SiO disc and disc wind from G17.64+0.16

Authors: L. T. Maud, R. Cesaroni, M. S. N. Kumar, F. F. S. van der Tak, V. Allen, M. G. Hoare, P. D. Klaassen, D. Harsono, M. R. Hogerheijde, Á. Sánchez-Monge, P. Schilke, A. Ahmadi, M. T. Beltrán, H. Beuther, T. Csengeri, S. Etoka, G. Fuller, R. Galván-Madrid, C. Goddi, Th. Henning, K. G. Johnston, R. Kuiper, S. Lumsden, L. Moscadelli, J. C. Mottram , et al. (6 additional authors not shown)

Abstract: We present high angular resolution 0.2 arcsec continuum and molecular emission line Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of G17.64+0.16 in Band 6 (220GHz) taken as part of a campaign in search of circumstellar discs around (proto)-O-stars. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peak… ▽ More We present high angular resolution 0.2 arcsec continuum and molecular emission line Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of G17.64+0.16 in Band 6 (220GHz) taken as part of a campaign in search of circumstellar discs around (proto)-O-stars. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. We detect SiO (5-4) emission that is marginally resolved and elongated in a direction perpendicular to the large-scale outflow seen in the 13CO (2-1) line using the main ALMA array in conjunction with the Atacama Compact Array (ACA). Morphologically, the SiO appears to represent a disc-like structure. Using parametric models we show that the position-velocity profile of the SiO is consistent with the Keplerian rotation of a disc around an object between 10-30Mo in mass, only if there is also radial expansion from a separate structure. The radial motion component can be interpreted as a disc wind from the disc surface. Models with a central stellar object mass between 20 and 30Mo are the most consistent with the stellar luminosity (100000 Lo) and indicative of an O-type star. The H30a millimetre recombination line (231.9GHz) is also detected, but spatially unresolved, and is indicative of a very compact, hot, ionised region co-spatial with the dust continuum core. Accounting for all observables, we suggest that G17.64 is consistent with a O-type young stellar object in the final stages of protostellar assembly, driving a wind, but that has not yet developed into a compact HII region. The existance and detection of the disc in G17.64 is likely related to its isolated and possibly more evolved nature, traits which may underpin discs in similar sources. △ Less

Submitted 9 October, 2018; originally announced October 2018.

Comments: 23 pages, 13 fig (inc 3 Appendix figures)

Journal ref: A&A 620, A31 (2018)

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

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

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

Submitted 27 September, 2018; originally announced September 2018.

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

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

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

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

Submitted 15 August, 2018; originally announced August 2018.

Comments: 31 pages, 26 figures, MNRAS accepted

Core fragmentation and Toomre stability analysis of W3(H2O): A case study of the IRAM NOEMA large program CORE

Authors: A. Ahmadi, H. Beuther, J. C. Mottram, F. Bosco, H. Linz, Th. Henning, J. M. Winters, R. Kuiper, R. Pudritz, Á. Sánchez-Monge, E. Keto, M. Beltran, S. Bontemps, R. Cesaroni, T. Csengeri, S. Feng, R. Galvan-Madrid, K. G. Johnston, P. Klaassen, S. Leurini, S. N. Longmore, S. Lumsden, L. T. Maud, K. M. Menten, L. Moscadelli , et al. (8 additional authors not shown)

Abstract: The fragmentation mode of high-mass molecular clumps and the properties of the central rotating structures surrounding the most luminous objects have yet to be comprehensively characterised. Using the IRAM NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30-m telescope, the CORE survey has obtained high-resolution observations of 20 well-known highly luminous star-forming regions in the 1.3… ▽ More The fragmentation mode of high-mass molecular clumps and the properties of the central rotating structures surrounding the most luminous objects have yet to be comprehensively characterised. Using the IRAM NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30-m telescope, the CORE survey has obtained high-resolution observations of 20 well-known highly luminous star-forming regions in the 1.37 mm wavelength regime in both line and dust continuum emission. We present the spectral line setup of the CORE survey and a case study for W3(H2O). At ~0.35" (700 AU at 2 kpc) resolution, the W3(H2O) clump fragments into two cores (West and East), separated by ~2300 AU. Velocity shifts of a few km/s are observed in the dense-gas tracer, CH3CN, across both cores, consistent with rotation and perpendicular to the directions of two bipolar outflows, one emanating from each core. The kinematics of the rotating structure about W3(H2O) W shows signs of differential rotation of material, possibly in a disk-like object. The observed rotational signature around W3(H2O) E may be due to a disk-like object, an unresolved binary (or multiple) system, or a combination of both. We fit the emission of CH3CN (12-11) K = 4-6 and derive a gas temperature map with a median temperature of ~165 K across W3(H2O). We create a Toomre Q map to study the stability of the rotating structures against gravitational instability. The rotating structures appear to be Toomre unstable close to their outer boundaries, with a possibility of further fragmentation in the differentially-rotating core W3(H2O) W. Rapid cooling in the Toomre-unstable regions supports the fragmentation scenario. Combining millimeter dust continuum and spectral line data toward the famous high-mass star-forming region W3(H2O), we identify core fragmentation on large scales, and indications for possible disk fragmentation on smaller spatial scales. △ Less

Submitted 1 August, 2018; originally announced August 2018.

Comments: 23 pages, 26 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 618, A46 (2018)

Flux Density Variations at 3.6 cm in the Massive Star-Forming Region W49A

Authors: C. G. De Pree, R. Galván-Madrid, W. M. Goss, R. S. Klessen, M. -M. Mac Low, T. Peters, D. Wilner, J. Bates, T. Melo, B. Presler-Marshall, R. Webb-Forgus

Abstract: A number of ultracompact H II regions in Galactic star forming environments have been observed to vary significantly in radio flux density on timescales of 10-20 years. Theory predicted that such variations should occur when the accretion flow that feeds a young massive star becomes unstable and clumpy. We have targeted the massive star-forming region W49A with the Karl G. Jansky Very Large Array… ▽ More A number of ultracompact H II regions in Galactic star forming environments have been observed to vary significantly in radio flux density on timescales of 10-20 years. Theory predicted that such variations should occur when the accretion flow that feeds a young massive star becomes unstable and clumpy. We have targeted the massive star-forming region W49A with the Karl G. Jansky Very Large Array (VLA) for observations at 3.6 cm with the B-configuration at 0.8'' resolution, to compare to nearly identical observations taken almost 21 years earlier (February 2015 and August 1994). Most of the sources in the crowded field of ultracompact and hypercompact H II regions exhibit no significant changes over this time period. However, one source, W49A/G2, decreased by 20% in peak intensity (from 71+/-4 mJy/beam to 57+/-3 mJy/beam), and 40% in integrated flux (from 0.109+/-0.011 Jy to 0.067+/-0.007 Jy), where we cite 5 sigma errors in peak intensity, and 10% errors in integrated flux. We present the radio images of the W49A region at the two epochs, the difference image that indicates the location of the flux density decrease, and discuss explanations for the flux density decrease near the position of W49A/G2. △ Less

Submitted 27 July, 2018; originally announced July 2018.

Comments: 9 pages, 2 figures, accepted for publication, The Astrophysical Journal Letters, July 26, 2018

Fragmentation and disk formation during high-mass star formation: The IRAM NOEMA (Northern Extended Millimeter Array) large program CORE

Authors: H. Beuther, J. C. Mottram, A. Ahmadi, F. Bosco, H. Linz, Th. Henning, P. Klaassen, J. M. Winters, L. T. Maud, R. Kuiper, D. Semenov, C. Gieser, T. Peters, J. S. Urquhart, R. Pudritz, S. E. Ragan, S. Feng, E. Keto, S. Leurini, R. Cesaroni, M. Beltran, A. Palau, A. Sanchez-Monge, R. Galvan-Madrid, Q. Zhang , et al. (8 additional authors not shown)

Abstract: Aims: We aim to understand the fragmentation as well as the disk formation, outflow generation and chemical processes during high-mass star formation on spatial scales of individual cores. Methods: Using the IRAM Northern Extended Millimeter Array (NOEMA) in combination with the 30m telescope, we have observed in the IRAM large program CORE the 1.37mm continuum and spectral line emission at high… ▽ More Aims: We aim to understand the fragmentation as well as the disk formation, outflow generation and chemical processes during high-mass star formation on spatial scales of individual cores. Methods: Using the IRAM Northern Extended Millimeter Array (NOEMA) in combination with the 30m telescope, we have observed in the IRAM large program CORE the 1.37mm continuum and spectral line emission at high angular resolution (~0.4'') for a sample of 20 well-known high-mass star-forming regions with distances below 5.5kpc and luminosities larger than 10^4Lsun. Results: We present the overall survey scope, the selected sample, the observational setup and the main goals of CORE. Scientifically, we concentrate on the mm continuum emission on scales on the order of 1000AU. We detect strong mm continuum emission from all regions, mostly due to the emission from cold dust. The fragmentation properties of the sample are diverse. We see extremes where some regions are dominated by a single high-mass core whereas others fragment into as many as 20 cores. A minimum-spanning-tree analysis finds fragmentation at scales on the order of the thermal Jeans length or smaller suggesting that turbulent fragmentation is less important than thermal gravitational fragmentation. The diversity of highly fragmented versus singular regions can be explained by varying initial density structures and/or different initial magnetic field strengths. Conclusions: The smallest observed separations between cores are found around the angular resolution limit which indicates that further fragmentation likely takes place on even smaller spatial scales. The CORE project with its numerous spectral line detections will address a diverse set of important physical and chemical questions in the field of high-mass star formation. △ Less

Submitted 3 May, 2018; originally announced May 2018.

Comments: 29 pages, 17 figures, Astronomy & Astrophysics in press, for a higher-resolution version of the paper see http://www.mpia.de/homes/beuther/papers.html

Journal ref: A&A 617, A100 (2018)

Ionised gas kinematics in bipolar H II regions

Authors: Hannah S. Dalgleish, Steven N. Longmore, Thomas Peters, Jonathan D. Henshaw, Joshua L. Veitch-Michaelis, James S. Urquhart

Abstract: Stellar feedback plays a fundamental role in shaping the evolution of galaxies. Here we explore the use of ionised gas kinematics in young, bipolar H II regions as a probe of early feedback in these star-forming environments. We have undertaken a multiwavelength study of a young, bipolar H II region in the Galactic disc, G$316.81-0.06$, which lies at the centre of a massive ($\sim10^3$ M… ▽ More Stellar feedback plays a fundamental role in shaping the evolution of galaxies. Here we explore the use of ionised gas kinematics in young, bipolar H II regions as a probe of early feedback in these star-forming environments. We have undertaken a multiwavelength study of a young, bipolar H II region in the Galactic disc, G$316.81-0.06$, which lies at the centre of a massive ($\sim10^3$ M$_{\odot}$) infrared-dark cloud filament. It is still accreting molecular gas as well as driving a $\sim 0.2$ pc ionised gas outflow perpendicular to the filament. Intriguingly, we observe a large velocity gradient ($47.81 \pm 3.21$ km s$^{-1}$ pc$^{-1}$) across the ionised gas in a direction perpendicular to the outflow. This kinematic signature of the ionised gas shows a reasonable correspondence with the simulations of young H II regions. Based on a qualitative comparison between our observations and these simulations, we put forward a possible explanation for the velocity gradients observed in G$316.81-0.06$. If the velocity gradient perpendicular to the outflow is caused by rotation of the ionised gas, then we infer that this rotation is a direct result of the initial net angular momentum in the natal molecular cloud. If this explanation is correct, this kinematic signature should be common in other young (bipolar) H II regions. We suggest that further quantitative analysis of the ionised gas kinematics of young H II regions, combined with additional simulations, should improve our understanding of feedback at these early stages. △ Less

Submitted 1 May, 2018; originally announced May 2018.

Comments: 14 pages and 9 figures. Accepted for publication in MNRAS

Accelerating infall and rotational spin-up in the hot molecular core G31.41+0.31

Authors: M. T. Beltrán, R. Cesaroni, V. M. Rivilla, Á. Sánchez-Monge, L. Moscadelli, A. Ahmadi, V. Allen, H. Beuther, S. Etoka, D. Galli, R. Galván-Madrid, C. Goddi, K. G. Johnston, A. Kölligan, R. Kuiper, M. S. N. Kumar, L. T. Maud, J. C. Mottram, T. Peters, P. Schilke, L. Testi, F. van der Tak, C. M. Walmsley

Abstract: As part of our effort to search for circumstellar disks around high-mass stellar objects, we observed the well-known core G31.41+0.31 with ALMA at 1.4 mm with an angular resolution of~0.22" (~1700 au). The dust continuum emission has been resolved into two cores namely Main and NE. The Main core, which has the stronger emission and is the more chemically rich, has a diameter of ~5300 au, and is as… ▽ More As part of our effort to search for circumstellar disks around high-mass stellar objects, we observed the well-known core G31.41+0.31 with ALMA at 1.4 mm with an angular resolution of~0.22" (~1700 au). The dust continuum emission has been resolved into two cores namely Main and NE. The Main core, which has the stronger emission and is the more chemically rich, has a diameter of ~5300 au, and is associated with two free-free continuum sources. The Main core looks featureless and homogeneous in dust continuum emission and does not present any hint of fragmentation. Each transition of CH3CN and CH3OCHO, both ground and vibrationally excited, as well as those of CH3CN isotopologues, shows a clear velocity gradient along the NE-SW direction, with velocity linearly increasing with distance from the center, consistent with solid-body rotation. However, when comparing the velocity field of transitions with different upper level energies, the rotation velocity increases with increasing energy of the transition, which suggests that the rotation speeds up towards the center. Spectral lines towards the dust continuum peak show an inverse P-Cygni profile that supports the existence of infall in the core. The infall velocity increases with the energy of the transition suggesting that the infall is accelerating towards the center of the core, consistent with gravitational collapse. Despite the monolithic appearance of the Main core, the presence of red-shifted absorption, the existence of two embedded free-free sources at the center, and the rotational spin-up are consistent with an unstable core undergoing fragmentation with infall and differential rotation due to conservation of angular momentum. Therefore, the most likely explanation for the monolithic morphology is that the large opacity of the dust emission prevents the detection of any inhomogeneity in the core. △ Less

Submitted 21 June, 2018; v1 submitted 14 March, 2018; originally announced March 2018.

Comments: 28 pages, 16 figures, and 7 tables. Accepted for publication in A&A

Radiative transfer calculations of the diffuse ionised gas in disc galaxies with cosmic ray feedback

Authors: Bert Vandenbroucke, Kenneth Wood, Philipp Girichidis, Alex Hill, Thomas Peters

Abstract: The large vertical scale heights of the diffuse ionised gas (DIG) in disc galaxies are challenging to model, as hydrodynamical models including only thermal feedback seem to be unable to support gas at these heights. In this paper, we use a three dimensional Monte Carlo radiation transfer code to post-process disc simulations of the Simulating the Life-Cycle of Molecular Clouds (SILCC) project tha… ▽ More The large vertical scale heights of the diffuse ionised gas (DIG) in disc galaxies are challenging to model, as hydrodynamical models including only thermal feedback seem to be unable to support gas at these heights. In this paper, we use a three dimensional Monte Carlo radiation transfer code to post-process disc simulations of the Simulating the Life-Cycle of Molecular Clouds (SILCC) project that include feedback by cosmic rays. We show that the more extended discs in simulations including cosmic ray feedback naturally lead to larger scale heights for the DIG which are more in line with observed scale heights. We also show that including a fiducial cosmic ray heating term in our model can help to increase the temperature as a function of disc scale height, but fails to reproduce observed DIG nitrogen and sulphur forbidden line intensities. We show that, to reproduce these line emissions, we require a heating mechanism that affects gas over a larger density range than is achieved by cosmic ray heating, which can be achieved by fine tuning the total luminosity of ionising sources to get an appropriate ionising spectrum as a function of scale height. This result sheds a new light on the relation between forbidden line emissions and temperature profiles for realistic DIG gas distributions. △ Less

Submitted 21 February, 2018; originally announced February 2018.

Comments: 14 pages, 11 figures, accepted for publication in Montly Notices of the Royal Astronomical Society

The evolution of young HII regions

Authors: P. D. Klaassen, K. G. Johnston, J. S. Urquhart, J. C. Mottram, T. Peters, R. Kuiper, H. Beuther, F. F. S. van der Tak, C. Goddi

Abstract: High-mass stars form in much richer environments than those associated with isolated low-mass stars, and once they reach a certain mass, produce ionised (HII) regions. The formation of these pockets of ionised gas are unique to the formation of high-mass stars (M $>8$ M$_\odot$), and present an excellent opportunity to study the final stages of accretion, which could include accretion through the… ▽ More High-mass stars form in much richer environments than those associated with isolated low-mass stars, and once they reach a certain mass, produce ionised (HII) regions. The formation of these pockets of ionised gas are unique to the formation of high-mass stars (M $>8$ M$_\odot$), and present an excellent opportunity to study the final stages of accretion, which could include accretion through the HII region itself. This study of the dynamics of the gas on both sides of these ionisation boundaries in very young HII regions aims to quantify the relationship between the HII regions and their immediate environments.We present high-resolution ($\sim$ 0.5$"$) ALMA observations of nine HII regions selected from the Red MSX Source (RMS) survey with compact radio emission and bolometric luminosities greater than 10$^4$ L$_\odot$. We focus on the initial presentation of the data, including initial results from the radio recombination line H29$α$, some complementary molecules, and the 256 GHz continuum emission. Of the six (out of nine) regions with H29$α$ detections, two appear to have cometary morphologies with velocity gradients across them, and two appear more spherical with velocity gradients suggestive of infalling ionised gas. The remaining two were either observed at low resolution or had signals that were too weak to draw robust conclusions. We also present a description of the interactions between the ionised and molecular gas (as traced by CS (J=5-4)), often (but not always) finding theHII region had cleared its immediate vicinity of molecules. Of our sample of nine, the observations of the two clusters expected to have the youngest HII regions (from previous radio observations) are suggestive of having infalling motions in the H29$α$ emission, which could be indicative of late stage accretion onto the stars despite the presence of an HII region. △ Less

Submitted 13 December, 2017; originally announced December 2017.

Comments: Accepted for publication in A&A. 25 pages, 15 figures, 9 tables, with Table A.2 made available electronically via the CDS

Journal ref: A&A 611, A99 (2018)

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

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

Abstract: We present 3D "zoom-in" simulations of the formation of two molecular clouds out of the galactic interstellar medium. We model the clouds - identified from the SILCC simulations - with a resolution of up to 0.06 pc using adaptive mesh refinement in combination with a chemical network to follow heating, cooling, and the formation of H$_2$ and CO including (self-) shielding. The two clouds are assem… ▽ More We present 3D "zoom-in" simulations of the formation of two molecular clouds out of the galactic interstellar medium. We model the clouds - identified from the SILCC simulations - with a resolution of up to 0.06 pc using adaptive mesh refinement in combination with a chemical network to follow heating, cooling, and the formation of H$_2$ and CO including (self-) shielding. The two clouds are assembled within a few million years with mass growth rates of up to $\sim$ 10$^{-2}$ M$_\mathrm{sun}$ yr$^{-1}$ and final masses of $\sim$ 50 000 M$_\mathrm{sun}$. A spatial resolution of $\lesssim$ 0.1 pc is required for convergence with respect to the mass, velocity dispersion, and chemical abundances of the clouds, although these properties also depend on the cloud definition such as based on density thresholds, H$_2$ or CO mass fraction. To avoid grid artefacts, the progressive increase of resolution has to occur within the free-fall time of the densest structures (1 - 1.5 Myr) and $\gtrsim$ 200 time steps should be spent on each refinement level before the resolution is progressively increased further. This avoids the formation of spurious, large-scale, rotating clumps from unresolved turbulent flows. While CO is a good tracer for the evolution of dense gas with number densities $n \geq$ 300 cm$^{-3}$, H$_2$ is also found for $n \lesssim 30$ cm$^{-3}$ due to turbulent mixing and becomes dominant at column densities around 30 - 50 M$_\mathrm{sun}$ pc$^{-2}$. The CO-to-H$_2$ ratio steadily increases within the first 2 Myr whereas $X_\mathrm{CO}\simeq$ 1 - 4 $\times$ 10$^{20}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$ is approximately constant since the CO(1-0) line quickly becomes optically thick. △ Less

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

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

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

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

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

Submitted 2 November, 2016; originally announced November 2016.

Comments: submitted

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

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

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

Submitted 20 October, 2016; originally announced October 2016.

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

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

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

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

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

Submitted 20 October, 2016; originally announced October 2016.

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

Journal ref: Monthly Notices of the Royal Astronomical Society 466 (2017), 3293-3308

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

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

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

Submitted 16 June, 2016; originally announced June 2016.

Comments: 23 pages; submitted to MNRAS

Simulating the Formation of Massive Protostars: I. Radiative Feedback and Accretion Disks

Authors: Mikhail Klassen, Ralph Pudritz, Rolf Kuiper, Thomas Peters, Robi Banerjee

Abstract: We present radiation hydrodynamic simulations of collapsing protostellar cores with initial masses of 30, 100, and 200 M$_{\odot}$. We follow their gravitational collapse and the formation of a massive protostar and protostellar accretion disk. We employ a new hybrid radiative feedback method blending raytracing techniques with flux-limited diffusion for a more accurate treatment of the temperatur… ▽ More We present radiation hydrodynamic simulations of collapsing protostellar cores with initial masses of 30, 100, and 200 M$_{\odot}$. We follow their gravitational collapse and the formation of a massive protostar and protostellar accretion disk. We employ a new hybrid radiative feedback method blending raytracing techniques with flux-limited diffusion for a more accurate treatment of the temperature and radiative force. In each case, the disk that forms becomes Toomre-unstable and develops spiral arms. This occurs between 0.35 and 0.55 freefall times and is accompanied by an increase in the accretion rate by a factor of 2-10. Although the disk becomes unstable, no other stars are formed. In the case of our 100 and 200 M$_{\odot}$ simulation, the star becomes highly super-Eddington and begins to drive bipolar outflow cavities that expand outwards. These radiatively-driven bubbles appear stable, and appear to be channeling gas back onto the protostellar accretion disk. Accretion proceeds strongly through the disk. After 81.4 kyr of evolution, our 30 M$_{\odot}$ simulation shows a star with a mass of 5.48 M$_{\odot}$ and a disk of mass 3.3 M$_{\odot}$, while our 100 M$_{\odot}$ simulation forms a 28.8 M$_{\odot}$ mass star with a 15.8 M$_{\odot}$ disk over the course of 41.6 kyr, and our 200 M$_{\odot}$ simulation forms a 43.7 M$_{\odot}$ star with an 18 M$_{\odot}$ disk in 21.9 kyr. In the absence of magnetic fields or other forms of feedback, the masses of the stars in our simulation do not appear limited by their own luminosities. △ Less

Submitted 7 April, 2016; v1 submitted 23 March, 2016; originally announced March 2016.

Comments: 24 pages, 14 figures. Accepted to The Astrophysical Journal

Impact of initial models and variable accretion rates on the pre-main-sequence evolution of massive and intermediate-mass stars and the early evolution of HII regions

Authors: Lionel Haemmerlé, Thomas Peters

Abstract: Massive star formation requires the accretion of gas at high rate while the star is already bright. Its actual luminosity depends sensitively on the stellar structure. We compute pre-main-sequence tracks for massive and intermediate-mass stars with variable accretion rates and study the evolution of stellar radius, effective temperature and ionizing luminosity, starting at $2\,M_\odot$ with convec… ▽ More Massive star formation requires the accretion of gas at high rate while the star is already bright. Its actual luminosity depends sensitively on the stellar structure. We compute pre-main-sequence tracks for massive and intermediate-mass stars with variable accretion rates and study the evolution of stellar radius, effective temperature and ionizing luminosity, starting at $2\,M_\odot$ with convective or radiative structures. The radiative case shows a much stronger swelling of the protostar for high accretion rates than the convective case. For radiative structures, the star is very sensitive to the accretion rate and reacts quickly to accretion bursts, leading to considerable changes in photospheric properties on timescales as short as 100 - 1000 yr. The evolution for convective structures is much less influenced by the instantaneous accretion rate, and produces a monotonically increasing ionizing flux that can be many orders of magnitude smaller than in the radiative case. For massive stars, it results in a delay of the HII region expansion by up to 10,000 yr. In the radiative case, the HII region can potentially be engulfed by the star during the swelling, which never happens in the convective case. We conclude that the early stellar structure has a large impact on the radiative feedback during the pre-main-sequence evolution of massive protostars and introduces an important uncertainty that should be taken into account. Because of their lower effective temperatures, our convective models may hint at a solution to an observed discrepancy between the luminosity distribution functions of massive young stellar objects and compact HII regions. △ Less

Submitted 1 March, 2016; originally announced March 2016.

Direct detection of precursors of gas giants formed by gravitational instability with the Atacama Large Millimetre/sub-millimetre Array

Authors: Lucio Mayer, Thomas Peters, Jaime E. Pineda, James Wadsley

Abstract: Phases of gravitational instability are expected in the early phases of disk evolution, when the disk mass is still a substantial fraction of the mass of the star. Disk fragmentation into sub-stellar objects could occur in the cold exterior part of the disk. Direct detection of massive gaseous clumps on their way to collapse into gas giant planets would offer an unprecedented test of the disk inst… ▽ More Phases of gravitational instability are expected in the early phases of disk evolution, when the disk mass is still a substantial fraction of the mass of the star. Disk fragmentation into sub-stellar objects could occur in the cold exterior part of the disk. Direct detection of massive gaseous clumps on their way to collapse into gas giant planets would offer an unprecedented test of the disk instability model. Here we use state-of-the-art 3D radiation-hydro simulations of disks undergoing fragmentation into massive gas giants, post-processed with the RADMC-3D ray-tracing code to produce dust continuum emission maps. These are then fed into the Common Astronomy Software Applications (CASA) ALMA simulator. The synthetic maps show that both overdense spiral arms and actual clumps at different stages of collapse can be detected with the Atacama Large Millimetre/sub-millimetre Array (ALMA) in the full configuration at the distance of the Ophiuchus star forming region (125 pc). The detection of clumps is particularly effective at shorter wavelengths (690 GHz) combining two resolutions with multi-scale clean. Furthermore, we show that a flux-based estimate of the mass of a protoplanetary clump can be from comparable to a factor of 3 higher than the gravitationally bound clump mass. The estimated mass depends on the assumed opacity, and on the gas temperature, which should be set using the input of radiation-hydro simulations. We conclude that ALMA has the capability to detect "smoking gun" systems that are a signpost of the disk instability model for gas giant planet formation. △ Less

Submitted 15 February, 2016; originally announced February 2016.

Evidence of Short Timescale Flux Density Variations of UC HII regions in Sgr B2 Main and North

Authors: C. G. De Pree, T. Peters, M. -M. Mac Low, D. J. Wilner, W. M. Goss, R. Galván-Madrid, E. R. Keto, R. S. Klessen, A. Monsrud

Abstract: We have recently published observations of significant flux density variations at 1.3 cm in HII regions in the star forming regions Sgr B2 Main and North (De Pree et al. 2014). To further study these variations, we have made new 7 mm continuum and recombination line observations of Sgr B2 at the highest possible angular resolution of the Karl G. Jansky Very Large Array (VLA). We have observed Sgr… ▽ More We have recently published observations of significant flux density variations at 1.3 cm in HII regions in the star forming regions Sgr B2 Main and North (De Pree et al. 2014). To further study these variations, we have made new 7 mm continuum and recombination line observations of Sgr B2 at the highest possible angular resolution of the Karl G. Jansky Very Large Array (VLA). We have observed Sgr B2 Main and North at 42.9 GHz and at 45.4 GHz in the BnA configuration (Main) and the A configuration (North). We compare these new data to archival VLA 7 mm continuum data of Sgr B2 Main observed in 2003 and Sgr B2 North observed in 2001. We find that one of the 41 known ultracompact and hypercompact HII regions in Sgr B2 (K2-North) has decreased $\sim$27% in flux density from 142$\pm$14 mJy to 103$\pm$10 mJy (2.3$σ$) between 2001 and 2012. A second source, F3c-Main has increased $\sim$30% in flux density from 82$\pm$8 mJy to 107 $\pm$11 mJy (1.8$σ$) between 2003 and 2012. F3c-Main was previously observed to increase in flux density at 1.3 cm over a longer time period between 1989 and 2012 (De Pree et al. 2014). An observation of decreasing flux density, such as that observed in K2-North, is particularly significant since such a change is not predicted by the classical hypothesis of steady expansion of HII regions during massive star accretion. Our new observations at 7 mm, along with others in the literature, suggest that the formation of massive stars occurs through time-variable and violent accretion. △ Less

Submitted 16 November, 2015; originally announced November 2015.

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

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

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

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

Submitted 22 October, 2015; originally announced October 2015.

Comments: ApJ Letters in press

Journal ref: Astrophysical Journal Letters 813 (2015) L27

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

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

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

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

Comments: 7 pages, 5 figures, ApJL, accepted

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

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

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

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

Comments: 26 pages, MNRAS, accepted

Mach number study of supersonic turbulence: The properties of the density field

Authors: Lukas Konstandin, Wolfram Schmidt, Philipp Girichidis, Thomas Peters, Rahul Shetty, Ralf S. Klessen

Abstract: We model driven, compressible, isothermal, turbulence with Mach numbers ranging from the subsonic ($\mathcal{M} \approx 0.65$) to the highly supersonic regime ($\mathcal{M}\approx 16 $). The forcing scheme consists both solenoidal (transverse) and compressive (longitudinal) modes in equal parts. We find a relation $σ_{s}^2 = \mathrm{b}\log{(1+\mathrm{b}^2\mathcal{M}^2)}$ between the Mach number an… ▽ More We model driven, compressible, isothermal, turbulence with Mach numbers ranging from the subsonic ($\mathcal{M} \approx 0.65$) to the highly supersonic regime ($\mathcal{M}\approx 16 $). The forcing scheme consists both solenoidal (transverse) and compressive (longitudinal) modes in equal parts. We find a relation $σ_{s}^2 = \mathrm{b}\log{(1+\mathrm{b}^2\mathcal{M}^2)}$ between the Mach number and the standard deviation of the logarithmic density with $\mathrm{b} = 0.457 \pm 0.007$. The density spectra follow $\mathcal{D}(k,\,\mathcal{M}) \propto k^{ζ(\mathcal{M})}$ with scaling exponents depending on the Mach number. We find $ζ(\mathcal{M}) = α\mathcal{M}^β$ with a coefficient $α$ that varies slightly with resolution, whereas $β$ changes systematically. We extrapolate to the limit of infinite resolution and find $α= -1.91 \pm 0.01,\, β=-0.30\pm 0.03$. The dependence of the scaling exponent on the Mach number implies a fractal dimension $D=2+0.96 \mathcal{M}^{-0.30}$. We determine how the scaling parameters depend on the wavenumber and find that the density spectra are slightly curved. This curvature gets more pronounced with increasing Mach number. We propose a physically motivated fitting formula $\mathcal{D}(k) = \mathcal{D}_0 k^{ζk^η}$ by using simple scaling arguments. The fit reproduces the spectral behaviour down to scales $k\approx 80$. The density spectrum follows a single power-law $η= -0.005 \pm 0.01$ in the low Mach number regime and the strongest curvature $η= -0.04 \pm 0.02$ for the highest Mach number. These values of $η$ represent a lower limit, as the curvature increases with resolution. △ Less

Submitted 11 June, 2015; originally announced June 2015.

Comments: 12 pages, 11 figures, MNRAS, submitted

Radiation Hydrodynamics using Characteristics on Adaptive Decomposed Domains for Massively Parallel Star Formation Simulations

Authors: Lars Buntemeyer, Robi Banerjee, Thomas Peters, Mikhail Klassen, Ralph E. Pudritz

Abstract: We present an algorithm for solving the radiative transfer problem on massively parallel computers using adaptive mesh refinement and domain decomposition. The solver is based on the method of characteristics which requires an adaptive raytracer that integrates the equation of radiative transfer. The radiation field is split into local and global components which are handled separately to overcome… ▽ More We present an algorithm for solving the radiative transfer problem on massively parallel computers using adaptive mesh refinement and domain decomposition. The solver is based on the method of characteristics which requires an adaptive raytracer that integrates the equation of radiative transfer. The radiation field is split into local and global components which are handled separately to overcome the non-locality problem. The solver is implemented in the framework of the magneto-hydrodynamics code FLASH and is coupled by an operator splitting step. The goal is the study of radiation in the context of star formation simulations with a focus on early disc formation and evolution. This requires a proper treatment of radiation physics that covers both the optically thin as well as the optically thick regimes and the transition region in particular. We successfully show the accuracy and feasibility of our method in a series of standard radiative transfer problems and two 3D collapse simulations resembling the early stages of protostar and disc formation. △ Less

Submitted 30 April, 2015; v1 submitted 19 January, 2015; originally announced January 2015.

Comments: 33 pages, 15 figures, prepared for submission to New Astronomy. Comments welcome

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

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

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

Submitted 8 December, 2014; originally announced December 2014.

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

Modelling the supernova-driven ISM in different environments

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

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

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

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

The physics of volume rendering

Authors: Thomas Peters

Abstract: Radiation transfer is an important topic in several physical disciplines, probably most prominently in astrophysics. Computer scientists use radiation transfer, among other things, for the visualisation of complex data sets with direct volume rendering. In this note, I point out the connection between physical radiation transfer and volume rendering, and I describe an implementation of direct volu… ▽ More Radiation transfer is an important topic in several physical disciplines, probably most prominently in astrophysics. Computer scientists use radiation transfer, among other things, for the visualisation of complex data sets with direct volume rendering. In this note, I point out the connection between physical radiation transfer and volume rendering, and I describe an implementation of direct volume rendering in the astrophysical radiation transfer code RADMC-3D. I show examples for the use of this module on analytical models and simulation data. △ Less

Submitted 22 October, 2014; originally announced October 2014.

Journal ref: Eur. J. Phys. 35 (2014) 065028

A general hybrid radiation transport scheme for star formation simulations on an adaptive grid

Authors: Mikhail Klassen, Rolf Kuiper, Ralph E. Pudritz, Thomas Peters, Robi Banerjee, Lars Buntemeyer

Abstract: Radiation feedback plays a crucial role in the process of star formation. In order to simulate the thermodynamic evolution of disks, filaments, and the molecular gas surrounding clusters of young stars, we require an efficient and accurate method for solving the radiation transfer problem. We describe the implementation of a hybrid radiation transport scheme in the adaptive grid-based FLASH genera… ▽ More Radiation feedback plays a crucial role in the process of star formation. In order to simulate the thermodynamic evolution of disks, filaments, and the molecular gas surrounding clusters of young stars, we require an efficient and accurate method for solving the radiation transfer problem. We describe the implementation of a hybrid radiation transport scheme in the adaptive grid-based FLASH general magnetohydrodynamics code. The hybrid scheme splits the radiative transport problem into a raytracing step and a diffusion step. The raytracer captures the first absorption event, as stars irradiate their environments, while the evolution of the diffuse component of the radiation field is handled by a flux-limited diffusion (FLD) solver. We demonstrate the accuracy of our method through a variety of benchmark tests including the irradiation of a static disk, subcritical and supercritical radiative shocks, and thermal energy equilibration. We also demonstrate the capability of our method for casting shadows and calculating gas and dust temperatures in the presence of multiple stellar sources. Our method enables radiation-hydrodynamic studies of young stellar objects, protostellar disks, and clustered star formation in magnetized, filamentary environments. △ Less

Submitted 15 October, 2014; originally announced October 2014.

Comments: 16 pages, 15 figures, accepted to ApJ

Fractal basins of escape and the formation of spiral arms in a galactic potential with a bar

Authors: Andreas Ernst, Thomas Peters

Abstract: We investigate the dynamics in the close vicinity of and within the critical area in a 2D effective galactic potential with a bar of Zotos. We have calculated Poincaré surfaces of section and the basins of escape. In both the Poincaré surfaces of section and the basins of escape we find numerical evidence for the existence of a separatrix which hinders orbits from escaping out of the bar region. W… ▽ More We investigate the dynamics in the close vicinity of and within the critical area in a 2D effective galactic potential with a bar of Zotos. We have calculated Poincaré surfaces of section and the basins of escape. In both the Poincaré surfaces of section and the basins of escape we find numerical evidence for the existence of a separatrix which hinders orbits from escaping out of the bar region. We present numerical evidence for the similarity between spiral arms of barred spiral galaxies and tidal tails of star clusters. △ Less

Submitted 2 July, 2014; originally announced July 2014.

Comments: 12 pages, 10 figures, accepted by MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc. 443, 2579-2589 (2014)

Low-metallicity star formation: Relative impact of metals and magnetic fields

Authors: Thomas Peters, Dominik R. G. Schleicher, Rowan J. Smith, Wolfram Schmidt, Ralf S. Klessen

Abstract: Low-metallicity star formation poses a central problem of cosmology, as it determines the characteristic mass scale and distribution for the first and second generations of stars forming in our Universe. Here, we present a comprehensive investigation assessing the relative impact of metals and magnetic fields, which may both be present during low-metallicity star formation. We show that the presen… ▽ More Low-metallicity star formation poses a central problem of cosmology, as it determines the characteristic mass scale and distribution for the first and second generations of stars forming in our Universe. Here, we present a comprehensive investigation assessing the relative impact of metals and magnetic fields, which may both be present during low-metallicity star formation. We show that the presence of magnetic fields generated via the small-scale dynamo stabilises the protostellar disc and provides some degree of support against fragmentation. In the absence of magnetic fields, the fragmentation timescale in our model decreases by a factor of ~10 at the transition from Z=0 to Z>0, with subsequently only a weak dependence on metallicity. Similarly, the accretion timescale of the cluster is set by the large-scale dynamics rather than the local thermodynamics. In the presence of magnetic fields, the primordial disc can become completely stable, therefore forming only one central fragment. At Z>0, the number of fragments is somewhat reduced in the presence of magnetic fields, though the shape of the mass spectrum is not strongly affected in the limits of the statistical uncertainties. The fragmentation timescale, however, increases by roughly a factor of 3 in the presence of magnetic fields. Indeed, our results indicate comparable fragmentation timescales in primordial runs without magnetic fields and Z>0 runs with magnetic fields. △ Less

Submitted 2 June, 2014; originally announced June 2014.

Comments: MNRAS in press

Journal ref: Monthly Notices of the Royal Astronomical Society 442 (2014), 3112-3126

Collective outflow from a small multiple stellar system

Authors: Thomas Peters, Pamela D. Klaassen, Mordecai-Mark Mac Low, Martin Schrön, Christoph Federrath, Michael D. Smith, Ralf S. Klessen

Abstract: The formation of high-mass stars is usually accompanied by powerful protostellar outflows. Such high-mass outflows are not simply scaled-up versions of their lower-mass counterparts, since observations suggest that the collimation degree degrades with stellar mass. Theoretically, the origins of massive outflows remain open to question because radiative feedback and fragmentation of the accretion f… ▽ More The formation of high-mass stars is usually accompanied by powerful protostellar outflows. Such high-mass outflows are not simply scaled-up versions of their lower-mass counterparts, since observations suggest that the collimation degree degrades with stellar mass. Theoretically, the origins of massive outflows remain open to question because radiative feedback and fragmentation of the accretion flow around the most massive stars, with M > 15 M_Sun, may impede the driving of magnetic disk winds. We here present a three-dimensional simulation of the early stages of core fragmentation and massive star formation that includes a subgrid-scale model for protostellar outflows. We find that stars that form in a common accretion flow tend to have aligned outflow axes, so that the individual jets of multiple stars can combine to form a collective outflow. We compare our simulation to observations with synthetic H_2 and CO observations and find that the morphology and kinematics of such a collective outflow resembles some observed massive outflows, such as Cepheus A and DR 21. We finally compare physical quantities derived from simulated observations of our models to the actual values in the models to examine the reliability of standard methods for deriving physical quantities, demonstrating that those methods indeed recover the actual values to within a factor of 2-3. △ Less

Submitted 22 April, 2014; originally announced April 2014.

Comments: ApJ in press, movies: http://www.ita.uni-heidelberg.de/~tpeters/collective.shtml

Journal ref: The Astrophysical Journal 788 (2014), 14