Showing 1–50 of 114 results for author: Offner, S S R

FORGE'd in FIRE III: The IMF in Quasar Accretion Disks from STARFORGE

Authors: Philip F. Hopkins, Michael Y. Grudic, Kyle Kremer, Stella S. R. Offner, David Guszejnov, Anna L. Rosen

Abstract: Recently, we demonstrated self-consistent formation of strongly-magnetized quasar accretion disks (QADs) from cosmological radiation-magnetohydrodynamic-thermochemical galaxy-star formation simulations, including the full STARFORGE physics shown previously to produce a reasonable IMF under typical ISM conditions. Here we study star formation and the stellar IMF in QADs, on scales from 100 au to 10… ▽ More Recently, we demonstrated self-consistent formation of strongly-magnetized quasar accretion disks (QADs) from cosmological radiation-magnetohydrodynamic-thermochemical galaxy-star formation simulations, including the full STARFORGE physics shown previously to produce a reasonable IMF under typical ISM conditions. Here we study star formation and the stellar IMF in QADs, on scales from 100 au to 10 pc from the SMBH. We show it is critical to include physics often previously neglected, including magnetic fields, radiation, and (proto)stellar feedback. Closer to the SMBH, star formation is suppressed, but the (rare) stars that do form exhibit top-heavy IMFs. Stars can form only in special locations (e.g. magnetic field switches) in the outer QAD. Protostars accrete their natal cores rapidly but then dynamically decouple from the gas and 'wander,' ceasing accretion on timescales ~100 yr. Their jets control initial core accretion, but the ejecta are 'swept up' into the larger-scale QAD flow without much dynamical effect. The strong tidal environment strongly suppresses common-core multiplicity. The IMF shape depends sensitively on un-resolved dynamics of protostellar disks (PSDs), as the global dynamical times can become incredibly short (< yr) and tidal fields are incredibly strong, so whether PSDs can efficiently transport angular momentum or fragment catastrophically at <10 au scales requires novel PSD simulations to properly address. Most analytic IMF models and analogies with planet formation in PSDs fail qualitatively to explain the simulation IMFs, though we discuss a couple of viable models. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: 39 pages, 25 figures, submitted to the Open Journal of Astrophysics. Comments welcome. Additional images and movies from the simulations available at http://www.tapir.caltech.edu/~phopkins/Site/animations/Movies_zoom.html

Probing the physics of star formation (ProPStar): I. First resolved maps of the electron fraction and cosmic-ray ionization rate in NGC 1333

Authors: Jaime E. Pineda, Olli Sipilä, Dominique M. Segura-Cox, Maria Teresa Valdivia-Mena, Roberto Neri, Michael Kuffmeier, Alexei V. Ivlev, Stella S. R. Offner, Maria Jose Maureira, Paola Caselli, Silvia Spezzano, Nichol Cunningham, Anika Schmiedeke, Mike Chen

Abstract: Electron fraction and cosmic-ray ionization rates (CRIR) in star-forming regions are important quantities in astrochemical modeling and are critical to the degree of coupling between neutrals, ions, and electrons, which regulates the dynamics of the magnetic field. However, these are difficult quantities to estimate. We aim to derive the electron fraction and CRIR maps of an active star-forming re… ▽ More Electron fraction and cosmic-ray ionization rates (CRIR) in star-forming regions are important quantities in astrochemical modeling and are critical to the degree of coupling between neutrals, ions, and electrons, which regulates the dynamics of the magnetic field. However, these are difficult quantities to estimate. We aim to derive the electron fraction and CRIR maps of an active star-forming region. We combined observations of the nearby NGC 1333 star-forming region carried out with the NOEMA interferometer and IRAM 30-m single dish to generate high spatial dynamic range maps of different molecular transitions. We used the DCO$^+$ and H$^{13}$CO$^+$ ratio (in addition to complementary data) to estimate the electron fraction and produce cosmic-ray ionization rate maps. We derived the first large-area electron fraction and CRIR resolved maps in a star-forming region, with typical values of $10^{-6.5}$ and $10^{-16.5}$ s$^{-1}$, respectively. The maps present clear evidence of enhanced values around embedded young stellar objects (YSOs). This provides strong evidence for locally accelerated cosmic rays. We also found a strong enhancement toward the northwest region in the map that might be related either to an interaction with a bubble or to locally generated cosmic rays by YSOs. We used the typical electron fraction and derived a MHD turbulence dissipation scale of 0.054 pc, which could be tested with future observations. We found a higher cosmic-ray ionization rate compared to the canonical value for $N({\rm H_2})=10^{21}-10^{23}$ cm$^{-2}$ of $10^{-17}$ s$^{-1}$ in the region, and it is likely generated by the accreting YSOs. The high value of the electron fraction suggests that new disks will form from gas in the ideal-MHD limit. This indicates that local enhancements of $ζ({\rm H_2})$, due to YSOs, should be taken into account in the analysis of clustered star formation. △ Less

Submitted 25 February, 2024; originally announced February 2024.

Comments: 13 pages, 13 Figures. Accepted to A&A. Abridged abstract. GitHub repo with code available at https://github.com/jpinedaf/NGC1333_NOEMA_Ions

Journal ref: A&A 686, A162 (2024)

Infall of material onto the filaments in Barnard 5

Authors: Spandan Choudhury, Jaime E. Pineda, Paola Caselli, Michael Chun-Yuan Chen, Stella S. R. Offner, Maria Teresa Valdivia-Mena

Abstract: Aims. We aim to study the structure and kinematics of the two filaments inside the subsonic core Barnard 5 in Perseus using high-resolution ($\approx$ 2400 au) NH3 data and a multi-component fit analysis. Methods. We used observations of NH3 (1,1) and (2,2) inversion transitions using the Very Large Array (VLA) and the Green Bank Telescope (GBT). We smoothed the data to a beam of 8'' to reliably… ▽ More Aims. We aim to study the structure and kinematics of the two filaments inside the subsonic core Barnard 5 in Perseus using high-resolution ($\approx$ 2400 au) NH3 data and a multi-component fit analysis. Methods. We used observations of NH3 (1,1) and (2,2) inversion transitions using the Very Large Array (VLA) and the Green Bank Telescope (GBT). We smoothed the data to a beam of 8'' to reliably fit multiple velocity components towards the two filamentary structures identified in B5. Results. Along with the core and cloud components, which dominate the flux in the line of sight, we detected two components towards the two filaments showing signs of infall. We also detected two additional components that can possibly trace new material falling into the subsonic core of B5. Conclusions. Following comparison with previous simulations of filament formation scenarios in planar geometry, we conclude that either the formation of the B5 filaments is likely to be rather cylindrically symmetrical or the filaments are magnetically supported. We also estimate infall rates of $1.6\times10^{-4}\,M_\odot\,yr^{-1}$ and $1.8\times10^{-4}\,M_\odot\,yr^{-1}$ (upper limits) for the material being accreted onto the two filaments. At these rates, the filament masses can change significantly during the core lifetime. We also estimate an upper limit of $3.5\times10^{-5}\,M_\odot\,yr^{-1}$ for the rate of possible infall onto the core itself. Accretion of new material onto cores indicates the need for a significant update to current core evolution models, where cores are assumed to evolve in isolation. △ Less

Submitted 28 December, 2023; originally announced December 2023.

Comments: 15 pages, 13 figures, accepted for publication on Astronomy and Astrophysics. The scripts used for the analysis presented in this paper can be found at https://github.com/SpandanCh/Barnard5_filaments

Co-Evolution of Stars and Gas: Using Analysis of Synthetic Observations to Investigate the Star-Gas Correlation in STARFORGE

Authors: Samuel Millstone, Robert Gutermuth, Stella S. R. Offner, Riwaj Pokhrel, Michael Y. Grudić

Abstract: We explore the relationship between stellar surface density and gas surface density (the star-gas or S-G correlation) in a 20,000 M$_{\odot}$ simulation from the STAR FORmation in Gaseous Environments (STARFORGE) Project. We create synthetic observations based on the Spitzer and Herschel telescopes by modeling active galactic nuclei contamination, smoothing based on angular resolution, cropping th… ▽ More We explore the relationship between stellar surface density and gas surface density (the star-gas or S-G correlation) in a 20,000 M$_{\odot}$ simulation from the STAR FORmation in Gaseous Environments (STARFORGE) Project. We create synthetic observations based on the Spitzer and Herschel telescopes by modeling active galactic nuclei contamination, smoothing based on angular resolution, cropping the field-of-view, and removing close neighbors and low-mass sources. We extract star-gas properties such as the dense gas mass fraction, the Class II:I ratio, and the S-G correlation ($Σ_{\rm YSO}/Σ_{\rm gas}$) from the simulation and compare them to observations of giant molecular clouds, young clusters, and star-forming regions, as well as to analytical models. We find that the simulation reproduces trends in the counts of young stellar objects and the median slope of the S-G correlation. This implies that the S-G correlation is not simply the result of observational biases but is in fact a real effect. However, other statistics, such as the Class II:I ratio and dense gas mass fraction, do not always match observed equivalents in nearby clouds. This motivates further observations covering the full simulation age range and more realistic modeling of cloud formation. △ Less

Submitted 17 October, 2023; originally announced October 2023.

Comments: 25 pages, 16 figures. To be published in The Astrophysical Journal

Stellar Populations in STARFORGE: The Origin and Evolution of Star Clusters and Associations

Authors: Juan P. Farias, Stella S. R. Offner, Michael Y. Grudić, Dávid Guszejnov, Anna L. Rosen

Abstract: Most stars form in highly clustered environments within molecular clouds, but eventually disperse into the distributed stellar field population. Exactly how the stellar distribution evolves from the embedded stage into gas-free associations and (bound) clusters is poorly understood. We investigate the long-term evolution of stars formed in the STARFORGE simulation suite -- a set of radiation-magne… ▽ More Most stars form in highly clustered environments within molecular clouds, but eventually disperse into the distributed stellar field population. Exactly how the stellar distribution evolves from the embedded stage into gas-free associations and (bound) clusters is poorly understood. We investigate the long-term evolution of stars formed in the STARFORGE simulation suite -- a set of radiation-magnetohydrodynamic simulations of star-forming turbulent clouds that include all key stellar feedback processes inherent to star formation. We use Nbody6++GPU to follow the evolution of the young stellar systems after gas removal. We use HDBSCAN to define stellar groups and analyze the stellar kinematics to identify the true bound star clusters. The conditions modeled by the simulations, i.e., global cloud surface densities below 0.15 g cm$^{-2}$,, star formation efficiencies below 15%, and gas expulsion timescales shorter than a free fall time, primarily produce expanding stellar associations and small clusters. The largest star clusters, which have $\sim$1000 bound members, form in the densest and lowest velocity dispersion clouds, representing $\sim$32 and 39% of the stars in the simulations, respectively. The cloud's early dynamical state plays a significant role in setting the classical star formation efficiency versus bound fraction relation. All stellar groups follow a narrow mass-velocity dispersion power law relation at 10 Myr with a power law index of 0.21. This correlation result in a distinct mass-size relationship for bound clusters. We also provide valuable constraints on the gas dispersal timescale during the star formation process and analyze the implications for the formation of bound systems. △ Less

Submitted 20 September, 2023; originally announced September 2023.

Comments: 20 Pages, 10 figures, submitted to MNRAS

Flow of gas detected from beyond the filaments to protostellar scales in Barnard 5

Authors: M. T. Valdivia-Mena, J. E. Pineda, D. M. Segura-Cox, P. Caselli, A. Schmiedeke, S. Choudhury, S. S. R. Offner, R. Neri, A. Goodman, G. A. Fuller

Abstract: The infall of gas from outside natal cores has proven to feed protostars after the main accretion phase (Class 0). This changes our view of star formation to a picture that includes asymmetric accretion (streamers), and a larger role of the environment. However, the connection between streamers and the filaments that prevail in star-forming regions is unknown. We investigate the flow of material t… ▽ More The infall of gas from outside natal cores has proven to feed protostars after the main accretion phase (Class 0). This changes our view of star formation to a picture that includes asymmetric accretion (streamers), and a larger role of the environment. However, the connection between streamers and the filaments that prevail in star-forming regions is unknown. We investigate the flow of material toward the filaments within Barnard 5 (B5) and the infall from the envelope to the protostellar disk of the embedded protostar B5-IRS1. Our goal is to follow the flow of material from the larger, dense core scale, to the protostellar disk scale. We present new HC$_3$N line data from the NOEMA and 30m telescopes covering the coherence zone of B5, together with ALMA H$_2$CO and C$^{18}$O maps toward the protostellar envelope. We fit multiple Gaussian components to the lines so as to decompose their individual physical components. We investigate the HC$_3$N velocity gradients to determine the direction of chemically-fresh gas flow. At envelope scales, we use a clustering algorithm to disentangle the different kinematic components within H$_2$CO emission. At dense core scales, HC$_3$N traces the infall from the B5 region toward the filaments. HC$_3$N velocity gradients are consistent with accretion toward the filament spines plus flow along them. We found a $\sim2800$ au streamer in H$_2$CO emission which is blueshifted with respect to the protostar and deposits gas at outer disk scales. The strongest velocity gradients at large scales curve toward the position of the streamer at small scales, suggesting a connection between both flows. Our analysis suggests that the gas can flow from the dense core to the protostar. This implies that the mass available for a protostar is not limited to its envelope, and can receiving chemically-unprocessed gas after the main accretion phase. △ Less

Submitted 1 August, 2023; v1 submitted 26 July, 2023; originally announced July 2023.

Comments: 25 pages, 27 figures, accepted for publication on Astronomy and Astrophysics. The scripts used for analysis can be seen at https://github.com/tere-valdivia/Barnard_5_infall

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

Alignment of dense molecular core morphology and velocity gradients with ambient magnetic fields

Authors: A. Pandhi, R. K. Friesen, L. Fissel, J. E. Pineda, P. Caselli, M. C-Y. Chen, J. Di Francesco, A. Ginsburg, H. Kirk, P. C. Myers, S. S. R. Offner, A. Punanova, F. Quan, E. Redaelli, E. Rosolowsky, S. Scibelli, Y. M. Seo, Y. Shirley

Abstract: Studies of dense core morphologies and their orientations with respect to gas flows and the local magnetic field have been limited to only a small sample of cores with spectroscopic data. Leveraging the Green Bank Ammonia Survey alongside existing sub-millimeter continuum observations and Planck dust polarization, we produce a cross-matched catalogue of 399 dense cores with estimates of core morph… ▽ More Studies of dense core morphologies and their orientations with respect to gas flows and the local magnetic field have been limited to only a small sample of cores with spectroscopic data. Leveraging the Green Bank Ammonia Survey alongside existing sub-millimeter continuum observations and Planck dust polarization, we produce a cross-matched catalogue of 399 dense cores with estimates of core morphology, size, mass, specific angular momentum, and magnetic field orientation. Of the 399 cores, 329 exhibit 2D $\mathrm{v}_\mathrm{LSR}$ maps that are well fit with a linear gradient, consistent with rotation projected on the sky. We find a best-fit specific angular momentum and core size relationship of $J/M \propto R^{1.82 \pm 0.10}$, suggesting that core velocity gradients originate from a combination of solid body rotation and turbulent motions. Most cores have no preferred orientation between the axis of core elongation, velocity gradient direction, and the ambient magnetic field orientation, favouring a triaxial and weakly magnetized origin. We find, however, strong evidence for a preferred anti-alignment between the core elongation axis and magnetic field for protostellar cores, revealing a change in orientation from starless and prestellar populations that may result from gravitational contraction in a magnetically-regulated (but not dominant) environment. We also find marginal evidence for anti-alignment between the core velocity gradient and magnetic field orientation in the L1228 and L1251 regions of Cepheus, suggesting a preferred orientation with respect to magnetic fields may be more prevalent in regions with locally ordered fields. △ Less

Submitted 24 July, 2023; originally announced July 2023.

Comments: 33 pages, 28 figures, accepted to MNRAS

Does God play dice with star clusters?

Authors: Michael Y. Grudić, Stella S. R. Offner, Dávid Guszejnov, Claude-André Faucher-Giguère, Philip F. Hopkins

Abstract: When a detailed model of a stellar population is unavailable, it is most common to assume that stellar masses are independently and identically distributed according to some distribution: the universal initial mass function (IMF). However, stellar masses resulting from causal, long-ranged physics cannot be truly random and independent, and the IMF may vary with environment. To compare stochastic s… ▽ More When a detailed model of a stellar population is unavailable, it is most common to assume that stellar masses are independently and identically distributed according to some distribution: the universal initial mass function (IMF). However, stellar masses resulting from causal, long-ranged physics cannot be truly random and independent, and the IMF may vary with environment. To compare stochastic sampling with a physical model, we run a suite of 100 STARFORGE radiation magnetohydrodynamics simulations of low-mass star cluster formation in $2000M_\odot$ clouds that form $\sim 200$ stars each on average. The stacked IMF from the simulated clouds has a sharp truncation at $\sim 28 M_\odot$, well below the typically-assumed maximum stellar mass $M_{\rm up} \sim 100-150M_\odot$ and the total cluster mass. The sequence of star formation is not totally random: massive stars tend to start accreting sooner and finish later than the average star. However, final cluster properties such as maximum stellar mass and total luminosity have a similar amount of cloud-to-cloud scatter to random sampling. Therefore stochastic sampling does not generally model the stellar demographics of a star cluster as it is forming, but may describe the end result fairly well, if the correct IMF -- and its environment-dependent upper cutoff -- are known. △ Less

Submitted 7 December, 2023; v1 submitted 30 June, 2023; originally announced July 2023.

Comments: Accepted for publication in Open Journal of Astrophysics

ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): A forming quadruple system with continuum `ribbons' and intricate outflows

Authors: Qiu-yi Luo, Tie Liu, Aaron T. Lee, Stella S. R. Offner, James di Francesco, Doug Johnstone, Mika Juvela, Paul F. Goldsmith, Sheng-Li Qin, Xiaofeng Mai, Xun-chuan Liu, Patricio Sanhueza, Feng-Wei Xu, Ken'ichi Tatematsu, Somnath Dutta, Huei-Ru Vivien Chen, Shanghuo Li, Aiyuan Yang, Sheng-Yuan Liu, Chin-Fei Lee, Naomi Hirano, Chang Won Lee, Dipen Sahu, Hsien Shang, Shih-Ying Hsu , et al. (9 additional authors not shown)

Abstract: One of the most poorly understood aspects of low-mass star formation is how multiple-star systems are formed. Here we present the results of Atacama Large Millimeter/submillimeter Array (ALMA) Band-6 observations towards a forming quadruple protostellar system, G206.93-16.61E2, in the Orion B molecular cloud. ALMA 1.3 mm continuum emission reveals four compact objects, of which two are Class I you… ▽ More One of the most poorly understood aspects of low-mass star formation is how multiple-star systems are formed. Here we present the results of Atacama Large Millimeter/submillimeter Array (ALMA) Band-6 observations towards a forming quadruple protostellar system, G206.93-16.61E2, in the Orion B molecular cloud. ALMA 1.3 mm continuum emission reveals four compact objects, of which two are Class I young stellar objects (YSOs), and the other two are likely in prestellar phase. The 1.3 mm continuum emission also shows three asymmetric ribbon-like structures that are connected to the four objects, with lengths ranging from $\sim$500 au to $\sim$2200 au. By comparing our data with magneto-hydrodynamic (MHD) simulations, we suggest that these ribbons trace accretion flows and also function as gas bridges connecting the member protostars. Additionally, ALMA CO J=2-1 line emission reveals a complicated molecular outflow associated with G206.93-16.61E2 with arc-like structures suggestive of an outflow cavity viewed pole-on. △ Less

Submitted 13 June, 2023; originally announced June 2023.

Comments: The paper was accepted by APJL

A Census of Outflow to Magnetic Field Orientations in Nearby Molecular Clouds

Authors: Duo Xu, Stella S. R. Offner, Robert Gutermuth, Jonathan C. Tan

Abstract: We define a sample of 200 protostellar outflows showing blue and redshifted CO emission in the nearby molecular clouds Ophiuchus, Taurus, Perseus and Orion to investigate the correlation between outflow orientations and local, but relatively large-scale, magnetic field directions traced by Planck 353 GHz dust polarization. At high significance (p~1e-4), we exclude a random distribution of relative… ▽ More We define a sample of 200 protostellar outflows showing blue and redshifted CO emission in the nearby molecular clouds Ophiuchus, Taurus, Perseus and Orion to investigate the correlation between outflow orientations and local, but relatively large-scale, magnetic field directions traced by Planck 353 GHz dust polarization. At high significance (p~1e-4), we exclude a random distribution of relative orientations and find that there is a preference for alignment of projected plane of sky outflow axes with magnetic field directions. The distribution of relative position angles peaks at ~30deg and exhibits a broad dispersion of ~50deg. These results indicate that magnetic fields have dynamical influence in regulating the launching and/or propagation directions of outflows. However, the significant dispersion around perfect alignment orientation implies that there are large measurement uncertainties and/or a high degree of intrinsic variation caused by other physical processes, such as turbulence or strong stellar dynamical interactions. Outflow to magnetic field alignment is expected to lead to a correlation in the directions of nearby outflow pairs, depending on the degree of order of the field. Analyzing this effect we find limited correlation, except on relatively small scales < 0.5 pc. Furthermore, we train a convolutional neural network to infer the inclination angle of outflows with respect to the line of sight and apply it to our outflow sample to estimate their full 3D orientations. We find that the angles between outflow pairs in 3D space also show evidence of small-scale alignment. △ Less

Submitted 7 November, 2022; originally announced November 2022.

Comments: ApJ Accepted

The cosmic ray ionisation and $γ$-ray budgets of star-forming galaxies

Authors: Mark R. Krumholz, Roland M. Crocker, Stella S. R. Offner

Abstract: Cosmic rays in star-forming galaxies are a dominant source of both diffuse $γ$-ray emission and ionisation in gas too deeply shielded for photons to penetrate. Though the cosmic rays responsible for $γ$-rays and ionisation are of different energies, they are produced by the same star formation-driven sources, and thus galaxies' star formation rates, $γ$-ray luminosities, and ionisation rates shoul… ▽ More Cosmic rays in star-forming galaxies are a dominant source of both diffuse $γ$-ray emission and ionisation in gas too deeply shielded for photons to penetrate. Though the cosmic rays responsible for $γ$-rays and ionisation are of different energies, they are produced by the same star formation-driven sources, and thus galaxies' star formation rates, $γ$-ray luminosities, and ionisation rates should all be linked. In this paper we use up-to-date cross-section data to determine this relationship, finding that cosmic rays in a galaxy of star formation rate $\dot{M}_*$ and gas depletion time $t_\mathrm{dep}$ produce a maximum primary ionisation rate $ζ\approx 1\times 10^{-16} (t_\mathrm{dep}/\mbox{Gyr})^{-1}$ s$^{-1}$ and a maximum $γ$-ray luminosity $L_γ\approx 4\times 10^{39} (\dot{M}_*/\mathrm{M}_\odot\mbox{ yr}^{-1})$ erg s$^{-1}$ in the 0.1 - 100 GeV band. These budgets imply either that the ionisation rates measured in Milky Way molecular clouds include a significant contribution from local sources that elevate them above the Galactic mean, or that CR-driven ionisation in the Milky Way is enhanced by sources not linked directly to star formation. Our results also imply that ionisation rates in starburst systems are only moderately enhanced compared to those in the Milky Way. Finally, we point out that measurements of $γ$-ray luminosities can be used to place constraints on galactic ionisation budgets in starburst galaxies that are nearly free of systematic uncertainties on the details of cosmic ray acceleration. △ Less

Submitted 7 February, 2023; v1 submitted 7 November, 2022; originally announced November 2022.

Comments: 19 pages, 8 figures, accepted for publication in MNRAS

SPYGLASS. II. The Multi-Generational and Multi-Origin Star Formation History of Cepheus Far North

Authors: Ronan Kerr, Adam L. Kraus, Simon J. Murphy, Daniel M. Krolikowski, Stella S. R. Offner, Benjamin M. Tofflemire, Aaron C. Rizzuto

Abstract: Young stellar populations provide a record of past star formation, and by establishing their members' dynamics and ages, it is possible to reconstruct the full history of star formation events. Gaia has greatly expanded the number of accessible stellar populations, with one of the most notable recently-discovered associations being Cepheus Far North (CFN), a population containing hundreds of membe… ▽ More Young stellar populations provide a record of past star formation, and by establishing their members' dynamics and ages, it is possible to reconstruct the full history of star formation events. Gaia has greatly expanded the number of accessible stellar populations, with one of the most notable recently-discovered associations being Cepheus Far North (CFN), a population containing hundreds of members spanning over 100 pc. With its proximity (d $\lesssim$ 200 pc), apparent substructure, and relatively small population, CFN represents a manageable population to study in depth, with enough evidence of internal complexity to produce a compelling star formation story. Using Gaia astrometry and photometry combined with additional spectroscopic observations, we identify over 500 candidate CFN members spread across 7 subgroups. Combining ages from isochrones, asteroseismology, dynamics, and lithium depletion, we produce well-constrained ages for all seven subgroups, revealing a largely continuous 10 Myr star formation history in the association. By tracing back the present-day populations to the time of their formation, we identify two spatially and dynamically distinct nodes in which stars form, one associated with $β$ Cephei which shows mostly co-spatial formation, and one associated with EE Draconis with a more dispersed star formation history. This detailed view of star formation demonstrates the complexity of the star formation process, even in the smallest of regions. △ Less

Submitted 15 December, 2022; v1 submitted 26 September, 2022; originally announced September 2022.

Comments: Accepted to ApJ; 34 pages, 15 figures, 6 tables in two-column AASTEX63 format

Effects of the environment on the multiplicity properties of stars in the STARFORGE simulations

Authors: Dávid Guszejnov, Aman N. Raju, Stella S. R. Offner, Michael Y. Grudić, Claude-André Faucher-Giguère, Philip F. Hopkins, Anna L. Rosen

Abstract: Most observed stars are part of a multiple star system, but the formation of such systems and the role of environment and various physical processes is still poorly understood. We present a suite of radiation-magnetohydrodynamic simulations of star-forming molecular clouds from the STARFORGE project that include stellar feedback with varied initial surface density, magnetic fields, level of turbul… ▽ More Most observed stars are part of a multiple star system, but the formation of such systems and the role of environment and various physical processes is still poorly understood. We present a suite of radiation-magnetohydrodynamic simulations of star-forming molecular clouds from the STARFORGE project that include stellar feedback with varied initial surface density, magnetic fields, level of turbulence, metallicity, interstellar radiation field, simulation geometry and turbulent driving. In our fiducial cloud the raw simulation data reproduces the observed multiplicity fractions for Solar-type and higher mass stars, similar to previous works. However, after correcting for observational incompleteness the simulation under-predicts these values. The discrepancy is likely due to the lack of disk fragmentation, as the simulation only resolves multiples that form either through capture or core fragmentation. The raw mass distribution of companions is consistent with randomly drawing from the initial mass function for the companions of $>1\,\mathrm{M_\odot}$ stars, however, accounting for observational incompleteness produces a flatter distribution similar to observations. We show that stellar multiplicity changes as the cloud evolves and anti-correlates with stellar density. This relationship also explains most multiplicity variations between runs, i.e., variations in the initial conditions that increase stellar density (increased surface density, reduced turbulence) decrease multiplicity. While other parameters, such as metallicity, interstellar radiation, and geometry significantly affect the star formation history or the IMF, varying them produces no clear trend in stellar multiplicity properties. △ Less

Submitted 28 September, 2022; v1 submitted 4 August, 2022; originally announced August 2022.

Comments: 20 pages, 21 figures, submitted to MNRAS

Effects of the environment and feedback physics on the initial mass function of stars in the STARFORGE simulations

Authors: Dávid Guszejnov, Michael Y. Grudić, Stella S. R. Offner, Claude-André Faucher-Giguère, Philip F. Hopkins, Anna L. Rosen

Abstract: One of the key mysteries of star formation is the origin of the stellar initial mass function (IMF). The IMF is observed to be nearly universal in the Milky Way and its satellites, and significant variations are only inferred in extreme environments, such as the cores of massive elliptical galaxies. In this work we present simulations from the STARFORGE project that are the first cloud-scale RMHD… ▽ More One of the key mysteries of star formation is the origin of the stellar initial mass function (IMF). The IMF is observed to be nearly universal in the Milky Way and its satellites, and significant variations are only inferred in extreme environments, such as the cores of massive elliptical galaxies. In this work we present simulations from the STARFORGE project that are the first cloud-scale RMHD simulations that follow individual stars and include all relevant physical processes. The simulations include detailed gas thermodynamics, as well as stellar feedback in the form of protostellar jets, stellar radiation, winds and supernovae. In this work we focus on how stellar radiation, winds and supernovae impact star-forming clouds. Radiative feedback plays a major role in quenching star formation and disrupting the cloud, however the IMF peak is predominantly set by protostellar jet physics. We find the effect of stellar winds is minor, and supernovae occur too late}to affect the IMF or quench star formation. We also investigate the effects of initial conditions on the IMF. The IMF is insensitive to the initial turbulence, cloud mass and cloud surface density, even though these parameters significantly shape the star formation history of the cloud, including the final star formation efficiency. The characteristic stellar mass depends weakly on metallicity and the interstellar radiation field. Finally, while turbulent driving and the level of magnetization strongly influences the star formation history, they only influence the high-mass slope of the IMF. △ Less

Submitted 2 August, 2022; v1 submitted 20 May, 2022; originally announced May 2022.

Comments: 24 pages, 20 figures, extra figures at https://github.com/guszejnovdavid/STARFORGE_IMF_paper_extra_plots

The Origin and Evolution of Multiple Star Systems

Authors: Stella S. R. Offner, Maxwell Moe, Kaitlin M. Kratter, Sarah I. Sadavoy, Eric L. N. Jensen, John J. Tobin

Abstract: Observational advances over the last decade have enabled high-resolution, interferometric studies of forming multiple systems, statistical surveys of multiplicity in star-forming regions, and new insights into disk evolution and planetary architectures in these systems. In this review, we compile the results of observational and theoretical studies of stellar multiplicity. We summarize the populat… ▽ More Observational advances over the last decade have enabled high-resolution, interferometric studies of forming multiple systems, statistical surveys of multiplicity in star-forming regions, and new insights into disk evolution and planetary architectures in these systems. In this review, we compile the results of observational and theoretical studies of stellar multiplicity. We summarize the population statistics spanning system evolution from the protostellar phase through the main-sequence phase and evaluate the influence of the local environment. In short, most stars are born in multiple stellar systems, most main sequence stars are members of multiple systems, but most star systems are single. We describe current models for the origin of stellar multiplicity and review the landscape of numerical simulations and assess their consistency with observations. We review the properties of disks and discuss the impact of multiplicity on planet formation and system architectures. Finally, we summarize open questions and discuss the technical requirements for future observational and theoretical progress. △ Less

Submitted 7 December, 2022; v1 submitted 18 March, 2022; originally announced March 2022.

Comments: 36 pages, 9 figures, Protostars and Planets VII, Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura

Cluster assembly and the origin of mass segregation in the STARFORGE simulations

Authors: Dávid Guszejnov, Carleen Markey, Stella S. R. Offner, Michael Y. Grudić, Claude-André Faucher-Giguère, Anna L. Rosen, Philip F. Hopkins

Abstract: Stars form in dense, clustered environments, where feedback from newly formed stars eventually ejects the gas, terminating star formation and leaving behind one or more star clusters. Using the STARFORGE simulations, it is possible to simulate this process in its entirety within a molecular cloud, while explicitly evolving the gas radiation and magnetic fields and following the formation of indivi… ▽ More Stars form in dense, clustered environments, where feedback from newly formed stars eventually ejects the gas, terminating star formation and leaving behind one or more star clusters. Using the STARFORGE simulations, it is possible to simulate this process in its entirety within a molecular cloud, while explicitly evolving the gas radiation and magnetic fields and following the formation of individual, low-mass stars. We find that individual star-formation sites merge to form ever larger structures, while still accreting gas. Thus clusters are assembled through a series of mergers. During the cluster assembly process a small fraction of stars are ejected from their clusters; we find no significant difference between the mass distribution of the ejected stellar population and that of stars inside clusters. The star-formation sites that are the building blocks of clusters start out mass segregated with one or a few massive stars at their center. As they merge the newly formed clusters maintain this feature, causing them to have mass-segregated substructures without themselves being centrally condensed. The merged clusters relax to a centrally condensed mass segregated configuration through dynamical interactions between their members, but this process does not finish before feedback expels the remaining gas from the cluster. In the simulated runs the gas-free clusters then become unbound and break up. We find that turbulent driving and a periodic cloud geometry can significantly reduce clustering and prevent gas expulsion. Meanwhile, the initial surface density and level of turbulence have little qualitative effect on cluster evolution, despite the significantly different star formation histories. △ Less

Submitted 19 May, 2022; v1 submitted 5 January, 2022; originally announced January 2022.

Comments: 18 pages, 15 figures, submitted to MNRAS

The dynamics and outcome of star formation with jets, radiation, winds, and supernovae in concert

Authors: Michael Y. Grudić, Dávid Guszejnov, Stella S. R. Offner, Anna L. Rosen, Aman N. Raju, Claude-André Faucher-Giguère, Philip F. Hopkins

Abstract: We analyze the first giant molecular cloud (GMC) simulation to follow the formation of individual stars and their feedback from jets, radiation, winds, and supernovae, using the STARFORGE framework in the GIZMO code. We evolve the GMC for $\sim 9 \rm Myr$, from initial turbulent collapse to dispersal by feedback. Protostellar jets dominate feedback momentum initially, but radiation and winds cause… ▽ More We analyze the first giant molecular cloud (GMC) simulation to follow the formation of individual stars and their feedback from jets, radiation, winds, and supernovae, using the STARFORGE framework in the GIZMO code. We evolve the GMC for $\sim 9 \rm Myr$, from initial turbulent collapse to dispersal by feedback. Protostellar jets dominate feedback momentum initially, but radiation and winds cause cloud disruption at $\sim 8\%$ star formation efficiency (SFE), and the first supernova at $8.3 \rm Myr$ comes too late to influence star formation significantly. The per-freefall SFE is dynamic, accelerating from 0 to $\sim 18\%$ before dropping quickly to <1%, but the estimate from YSO counts compresses it to a narrower range. The primary cluster forms hierarchically and condenses to a brief ($\sim 1\,\mathrm{Myr}$) compact ($\sim 1 \rm pc$) phase, but does not virialize before the cloud disperses, and the stars end as an unbound expanding association. The initial mass function resembles the Chabrier (2005) form with a high-mass slope $α=-2$ and a maximum mass of $55 M_\odot$. Stellar accretion takes $\sim 400 \rm kyr$ on average, but $\gtrsim 1\rm Myr$ for $>10 M_\odot$ stars, so massive stars finish growing latest. The fraction of stars in multiples increases as a function of primary mass, as observed. Overall, the simulation much more closely resembles reality, compared to variations which neglect different feedback physics entirely. But more detailed comparison with synthetic observations is necessary to constrain the theoretical uncertainties. △ Less

Submitted 24 March, 2022; v1 submitted 3 January, 2022; originally announced January 2022.

Comments: Accepted for publication in MNRAS. 17 pages, 10 figures. Comments welcome. Some data products and plotting code available at https://github.com/mikegrudic/StarforgeFullPhysics. Visit http://www.starforge.space to learn more about the STARFORGE Project. Simulation movie available at http://www.starforge.space/M2e4_fullphysics_flyaround.mp4

The VLA/ALMA Nascent Disk and Multiplicity (VANDAM) Survey of Orion Protostars V. A Characterization of Protostellar Multiplicity

Authors: John J. Tobin, Stella S. R. Offner, Kaitlin M. Kratter, S. Thomas Megeath, Patrick D. Sheehan, Leslie W. Looney, Ana Karla Diaz-Rodriguez, Mayra Osorio, Guillem Anglada, Sarah I. Sadavoy, Elise Furlan, Dominique Segura-Cox, Nicole Karnath, Merel L. R. van 't Hoff, Ewine F. van Dishoeck, Zhi-Yun Li, Rajeeb Sharma, Amelia M. Stutz, Lukasz Tychoniec

Abstract: We characterize protostellar multiplicity in the Orion molecular clouds using ALMA 0.87~mm and VLA 9~mm continuum surveys toward 328 protostars. These observations are sensitive to projected spatial separations as small as $\sim$20~au, and we consider source separations up to 10$^4$~au as potential companions. The overall multiplicity fraction (MF) and companion fraction (CF) for the Orion protost… ▽ More We characterize protostellar multiplicity in the Orion molecular clouds using ALMA 0.87~mm and VLA 9~mm continuum surveys toward 328 protostars. These observations are sensitive to projected spatial separations as small as $\sim$20~au, and we consider source separations up to 10$^4$~au as potential companions. The overall multiplicity fraction (MF) and companion fraction (CF) for the Orion protostars are 0.30$\pm$0.03 and 0.44$\pm$0.03, respectively, considering separations from 20 to 10$^4$~au. The MFs and CFs are corrected for potential contamination by unassociated young stars using a probabilistic scheme based on the surface density of young stars around each protostar. The companion separation distribution as a whole is double peaked and inconsistent with the separation distribution of solar-type field stars, while the separation distribution of Flat Spectrum protostars is consistent solar-type field stars. The multiplicity statistics and companion separation distributions of the Perseus star-forming region are consistent with those of Orion. Based on the observed peaks in the Class 0 separations at $\sim$100~au and $\sim$10$^3$~au, we argue that multiples with separations $<$500~au are likely produced by both disk fragmentation and turbulent fragmentation with migration, and those at $\ga$10$^3$~au result primarily from turbulent fragmentation. We also find that MFs/CFs may rise from Class 0 to Flat Spectrum protostars between 100 and 10$^3$~au in regions of high YSO density. This finding may be evidence for migration of companions from $>$10$^3$~au to $<$10$^3$~au, and that some companions between 10$^3$ and 10$^4$~au must be (or become) unbound. △ Less

Submitted 8 November, 2021; originally announced November 2021.

Comments: 76 pages, 20 Figures, 10 Tables, accepted by the Astrophysical Journal

Less wrong: a more realistic initial condition for simulations of turbulent molecular clouds

Authors: Henry B. Lane, Michael Y. Grudić, Dávid Guszejnov, Stella S. R. Offner, Claude-André Faucher-Giguère, Anna L. Rosen

Abstract: Simulations of isolated giant molecular clouds (GMCs) are an important tool for studying the dynamics of star formation, but their turbulent initial conditions (ICs) are uncertain. Most simulations have either initialized a velocity field with a prescribed power spectrum on a smooth density field (failing to model the full structure of turbulence) or "stirred" turbulence with periodic boundary con… ▽ More Simulations of isolated giant molecular clouds (GMCs) are an important tool for studying the dynamics of star formation, but their turbulent initial conditions (ICs) are uncertain. Most simulations have either initialized a velocity field with a prescribed power spectrum on a smooth density field (failing to model the full structure of turbulence) or "stirred" turbulence with periodic boundary conditions (which may not model real GMC boundary conditions). We develop and test a new GMC simulation setup (called TURBSPHERE) that combines advantages of both approaches: we continuously stir an isolated cloud to model the energy cascade from larger scales, and use a static potential to confine the gas. The resulting cloud and surrounding envelope achieve a quasi-equilibrium state with the desired hallmarks of supersonic ISM turbulence (e.g. density PDF and a $\sim k^{-2}$ velocity power spectrum), whose bulk properties can be tuned as desired. We use the final stirred state as initial conditions for star formation simulations with self-gravity, both with and without continued driving and protostellar jet feedback, respectively. We then disentangle the respective effects of the turbulent cascade, simulation geometry, external driving, and gravity/MHD boundary conditions on the resulting star formation. Without external driving, the new setup obtains results similar to previous simple spherical cloud setups, but external driving can suppress star formation considerably in the new setup. Periodic box simulations with the same dimensions and turbulence parameters form stars significantly slower, highlighting the importance of boundary conditions and the presence or absence of a global collapse mode in the results of star formation calculations. △ Less

Submitted 27 October, 2021; originally announced October 2021.

Comments: Submitted to MNRAS. 11 pages, 8 figures. See http://www.starforge.space/turbsphere_sigma_gas.mp4 for an animated version of Figure 2

Are massive dense clumps truly sub-virial? A new analysis using Gould Belt ammonia data

Authors: Ayushi Singh, Christopher D. Matzner, Rachel K. Friesen, Peter G. Martin, Jaime E. Pineda, Erik W. Rosolowsky, Felipe Alves, Ana Chacón-Tanarro, Hope How-Huan Chen, Michael Chun-Yuan Chen, Spandan Choudhury, James Di Francesco, Jared Keown, Helen Kirk, Anna Punanova, Youngmin Seo, Yancy Shirley, Adam Ginsburg, Stella S. R. Offner, Héctor G. Arce, Paola Caselli, Alyssa A. Goodman, Philip C. Myers, Elena Redaelli

Abstract: Dynamical studies of dense structures within molecular clouds often conclude that the most massive clumps contain too little kinetic energy for virial equilibrium, unless they are magnetized to an unexpected degree. This raises questions about how such a state might arise, and how it might persist long enough to represent the population of massive clumps. In an effort to re-examine the origins of… ▽ More Dynamical studies of dense structures within molecular clouds often conclude that the most massive clumps contain too little kinetic energy for virial equilibrium, unless they are magnetized to an unexpected degree. This raises questions about how such a state might arise, and how it might persist long enough to represent the population of massive clumps. In an effort to re-examine the origins of this conclusion, we use ammonia line data from the Green Bank Ammonia Survey and Planck-calibrated dust emission data from Herschel to estimate the masses and kinetic and gravitational energies for dense clumps in the Gould Belt clouds. We show that several types of systematic error can enhance the appearance of low kinetic-to-gravitational energy ratios: insufficient removal of foreground and background material; ignoring the kinetic energy associated with velocity differences across a resolved cloud; and over-correcting for stratification when evaluating the gravitational energy. Using an analysis designed to avoid these errors, we find that the most massive Gould Belt clumps harbor virial motions, rather than sub-virial ones. As a byproduct, we present a catalog of masses, energies, and virial energy ratios for 85 Gould Belt clumps. △ Less

Submitted 11 August, 2021; originally announced August 2021.

Comments: Submitted to ApJ

TIMES II: Investigating the Relation Between Turbulence and Star-forming Environments in Molecular Clouds

Authors: Hyeong-Sik Yun, Jeong-Eun Lee, Neal J. Evans II, Stella S. R. Offner, Mark H. Heyer, Jungyeon Cho, Brandt A. L. Gaches, Yao-Lun Yang, How-Huan Chen, Yunhee Choi, Yong-Hee Lee, Giseon Baek, Minho Choi, Jongsoo Kim, Hyunwoo Kang, Seokho Lee, Ken'ichi Tatematsu

Abstract: We investigate the effect of star formation on turbulence in the Orion A and Ophiuchus clouds using principal component analysis (PCA). We measure the properties of turbulence by applying PCA on the spectral maps in $^{13}$CO, C$^{18}$O, HCO$^+$ $J=$1$-$0, and CS $J=$2$-$1. First, the scaling relations derived from PCA of the $^{13}$CO maps show that the velocity difference ($δv$) for a given spat… ▽ More We investigate the effect of star formation on turbulence in the Orion A and Ophiuchus clouds using principal component analysis (PCA). We measure the properties of turbulence by applying PCA on the spectral maps in $^{13}$CO, C$^{18}$O, HCO$^+$ $J=$1$-$0, and CS $J=$2$-$1. First, the scaling relations derived from PCA of the $^{13}$CO maps show that the velocity difference ($δv$) for a given spatial scale ($L$) is the highest in the integral shaped filament (ISF) and L1688, where the most active star formation occurs in the two clouds. The $δv$ increases with the number density and total bolometric luminosity of the protostars in the sub-regions. Second, in the ISF and L1688 regions, the $δv$ of C$^{18}$O, HCO$^+$, and CS are generally higher than that of $^{13}$CO, which implies that the dense gas is more turbulent than the diffuse gas in the star-forming regions; stars form in dense gas, and dynamical activities associated with star formation, such as jets and outflows, can provide energy into the surrounding gas to enhance turbulent motions. △ Less

Submitted 28 July, 2021; originally announced July 2021.

Comments: 29 pages, 17 figures, 7 tables, accepted for publication in ApJ

Blowing Bubbles around Intermediate-Mass Stars: Feedback from Main-Sequence Winds is not Enough

Authors: Anna L. Rosen, Stella S. R. Offner, Michael M. Foley, Laura A. Lopez

Abstract: Numerous spherical ``shells" have been observed in young star-forming environments that host low- and intermediate-mass stars. These observations suggest that these shells may be produced by isotropic stellar wind feedback from young main-sequence stars. However, the driving mechanism for these shells remains uncertain because the momentum injected by winds is too low to explain their sizes and dy… ▽ More Numerous spherical ``shells" have been observed in young star-forming environments that host low- and intermediate-mass stars. These observations suggest that these shells may be produced by isotropic stellar wind feedback from young main-sequence stars. However, the driving mechanism for these shells remains uncertain because the momentum injected by winds is too low to explain their sizes and dynamics due to their low mass-loss rates. However, these studies neglect how the wind kinetic energy is transferred to the ISM and instead assume it is instantly lost via radiation, suggesting that these shells are momentum-driven. Intermediate-mass stars have fast ($v_w \gtrsim 1000$ km/s) stellar winds and therefore the energy injected by winds should produce energy-driven adiabatic wind bubbles that are larger than momentum-driven wind bubbles. Here, we explore if energy-driven wind feedback can produce the observed shells by performing a series of 3D magneto-hydrodynamic simulations of wind feedback from intermediate-mass and high-mass stars that are placed in a magnetized, turbulent molecular cloud. We find that, for the high-mass stars modeled, energy-driven wind feedback produces $\sim$pc scale wind bubbles in molecular clouds that agree with the observed shell sizes but winds from intermediate-mass stars can not produce similar shells because of their lower mass-loss rates and velocities. Therefore, such shells must be driven by other feedback processes inherent to low- and intermediate-mass star formation. △ Less

Submitted 26 July, 2021; originally announced July 2021.

Comments: 17 pages, 5 figures. Submitted to the Astrophysical Journal

TIMES I: a Systematic Observation in Multiple Molecular Lines Toward the Orion A and Ophiuchus Clouds

Authors: Hyeong-Sik Yun, Jeong-Eun Lee, Yunhee Choi, Neal J. Evans II, Stella S. R. Offner, Mark H. Heyer, Brandt A. L. Gaches, Yong-Hee Lee, Giseon Baek, Minho Choi, Hyunwoo Kang, Seokho Lee, Ken'ichi Tatematsu, Yao-Lun Yang, How-Huan Chen, Youngung Lee, Jae Hoon Jung, Changhoon Lee, Jungyeon Cho

Abstract: We have used the Taeduk Radio Astronomy Observatory to observe the Orion A and Ophiuchus clouds in the $J=$1$-$0 lines of $^{13}$CO, C$^{18}$O, HCN, HCO$^+$, and N$_2$H$^+$ and the $J=$2$-$1 line of CS. The fully sampled maps with uniform noise levels are used to create moment maps. The variations of the line intensity and velocity dispersion with total column density, derived from dust emission m… ▽ More We have used the Taeduk Radio Astronomy Observatory to observe the Orion A and Ophiuchus clouds in the $J=$1$-$0 lines of $^{13}$CO, C$^{18}$O, HCN, HCO$^+$, and N$_2$H$^+$ and the $J=$2$-$1 line of CS. The fully sampled maps with uniform noise levels are used to create moment maps. The variations of the line intensity and velocity dispersion with total column density, derived from dust emission maps, are presented and compared to previous work. The CS line traces dust column density over more than one order of magnitude, and the N$_2$H$^+$ line best traces the highest column density regime ($\log(N_\mathrm{H_2}$) $>$ 22.8). Line luminosities, integrated over the cloud, are compared to those seen in other galaxies. The HCO$^+$-to-HCN luminosity ratio in the Orion A cloud is similar to that of starburst galaxies, while that in the Ophiuchus cloud is in between those of active galactic nuclei and starburst galaxies. △ Less

Submitted 8 June, 2021; originally announced June 2021.

Comments: 62 pages, 53 figures, 6 tables, accepted for publication in ApJS

Stars with Photometrically Young Gaia Luminosities Around the Solar System (SPYGLASS) I: Mapping Young Stellar Structures and their Star Formation Histories

Authors: Ronan Kerr, Aaron C. Rizzuto, Adam L. Kraus, Stella S. R. Offner

Abstract: Young stellar associations hold a star formation record that can persist for millions of years, revealing the progression of star formation long after the dispersal of the natal cloud. To identify nearby young stellar populations that trace this progression, we have designed a comprehensive framework for the identification of young stars, and use it to identify $\sim$3$\times 10^4$ candidate young… ▽ More Young stellar associations hold a star formation record that can persist for millions of years, revealing the progression of star formation long after the dispersal of the natal cloud. To identify nearby young stellar populations that trace this progression, we have designed a comprehensive framework for the identification of young stars, and use it to identify $\sim$3$\times 10^4$ candidate young stars within a distance of 333 pc using Gaia DR2. Applying the HDBSCAN clustering algorithm to this sample, we identify 27 top-level groups, nearly half of which have little to no presence in previous literature. Ten of these groups have visible substructure, including notable young associations such as Orion, Perseus, Taurus, and Sco-Cen. We provide a complete subclustering analysis on all groups with substructure, using age estimates to reveal each region's star formation history. The patterns we reveal include an apparent star formation origin for Sco-Cen along a semicircular arc, as well as clear evidence for sequential star formation moving away from that arc with a propagation speed of $\sim$4 km s$^{-1}$ ($\sim$4 pc Myr$^{-1}$). We also identify earlier bursts of star formation in Perseus and Taurus that predate current, kinematically identical active star-forming events, suggesting that the mechanisms that collect gas can spark multiple generations of star formation, punctuated by gas dispersal and cloud regrowth. The large spatial scales and long temporal scales on which we observe star formation offer a bridge between the processes within individual molecular clouds and the broad forces guiding star formation at galactic scales. △ Less

Submitted 16 June, 2021; v1 submitted 19 May, 2021; originally announced May 2021.

Comments: Accepted to ApJ; 55 pages, 32 figures, 7 tables in two-column AASTEX63 format

The Single-Cloud Star Formation Relation

Authors: Riwaj Pokhrel, Robert A. Gutermuth, Mark R. Krumholz, Christoph Federrath, Mark Heyer, Shivan Khullar, S. Thomas Megeath, Philip C. Myers, Stella S. R. Offner, Judith L. Pipher, William J. Fischer, Thomas Henning, Joseph L. Hora

Abstract: One of the most important and well-established empirical results in astronomy is the Kennicutt-Schmidt (KS) relation between the density of interstellar gas and the rate at which that gas forms stars. A tight correlation between these quantities has long been measured at galactic scales. More recently, using surveys of YSOs, a KS relationship has been found within molecular clouds relating the sur… ▽ More One of the most important and well-established empirical results in astronomy is the Kennicutt-Schmidt (KS) relation between the density of interstellar gas and the rate at which that gas forms stars. A tight correlation between these quantities has long been measured at galactic scales. More recently, using surveys of YSOs, a KS relationship has been found within molecular clouds relating the surface density of star formation to the surface density of gas; however, the scaling of these laws varies significantly from cloud to cloud. In this Letter, we use a recently developed, high-accuracy catalog of young stellar objects from $\textit{Spitzer}$ combined with high-dynamic-range gas column density maps of twelve nearby ($<$1.5 kpc) molecular clouds from $\textit{Herschel}$ to re-examine the KS relation within individual molecular clouds. We find a tight, linear correlation between clouds' star formation rate per unit area and their gas surface density normalized by the gas free-fall time. The measured intracloud KS relation, which relates star formation rate to the volume density, extends over more than two orders of magnitude within each cloud and is nearly identical in each of the twelve clouds, implying a constant star formation efficiency per free-fall time $ε_{\rm ff}\approx 0.026$. The finding of a universal correlation within individual molecular clouds, including clouds that contain no massive stars or massive stellar feedback, favors models in which star formation is regulated by local processes such as turbulence or stellar feedback such as protostellar outflows, and disfavors models in which star formation is regulated only by galaxy properties or supernova feedback on galactic scales. △ Less

Submitted 9 April, 2021; originally announced April 2021.

Comments: Accepted in ApJL. 16 pages, 4 figures

Transition from Coherent Cores to Surrounding Cloud in L1688

Authors: Spandan Choudhury, Jaime E. Pineda, Paola Caselli, Stella S. R. Offner, Erik Rosolowsky, Rachel K. Friesen, Elena Redaelli, Ana Chacón-Tanarro, Yancy Shirley, Anna Punanova, Helen Kirk

Abstract: Stars form in cold dense cores showing subsonic velocity dispersions. The parental molecular clouds display higher temperatures and supersonic velocity dispersions. The transition from core to cloud has been observed in velocity dispersion, but temperature and abundance variations are unknown. We aim to study the transition from cores to ambient cloud in temperature and velocity dispersion using a… ▽ More Stars form in cold dense cores showing subsonic velocity dispersions. The parental molecular clouds display higher temperatures and supersonic velocity dispersions. The transition from core to cloud has been observed in velocity dispersion, but temperature and abundance variations are unknown. We aim to study the transition from cores to ambient cloud in temperature and velocity dispersion using a single tracer. We use NH3 (1,1) and (2,2) maps in L1688 from the Green Bank Ammonia Survey, smoothed to 1', and determine the physical properties from fits. We identify the coherent cores and study the changes in temperature and velocity dispersion from cores to the surrounding cloud. We obtain a kinetic temperature map tracing the extended cloud, improving from previous maps tracing mostly the cores. The cloud is 4-6 K warmer than the cores, and shows a larger velocity dispersion (diff. = 0.15-0.25 km/s). Comparing to Herschel-based measurements, we find that cores show kinetic temperature $\approx$1.8 K lower than the dust temperature; while the gas temperature is higher than the dust temperature in the cloud. We find an average p-NH3 fractional abundance (with respect to H2) of $(4.2\pm0.2) \times 10^{-9}$ towards the coherent cores, and $(1.4\pm0.1) \times 10^{-9}$ outside the core boundaries. Using stacked spectra, we detect two components, one narrow and one broad, towards cores and their neighbourhoods. We find the turbulence in the narrow component to be correlated to the size of the structure (Pearson-r=0.54). With these unresolved regional measurements, we obtain a turbulence-size relation of $σ_{v,NT}\propto r^{0.5}$, similar to previous findings using multiple tracers. We discover that the subsonic component extends up to 0.15 pc beyond the typical coherent boundaries, unveiling larger extents of the coherent cores and showing gradual transition to coherence over ~0.2 pc. △ Less

Submitted 12 February, 2021; originally announced February 2021.

Comments: 37 pages, 33 figures, 1 table. Accepted for publication in A&A

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

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

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

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

Submitted 1 January, 2021; originally announced January 2021.

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

Application of Convolutional Neural Networks to Identify Protostellar Outflows in CO Emission

Authors: Duo Xu, Stella S. R. Offner, Robert Gutermuth, Colin Van Oort

Abstract: We adopt the deep learning method CASI-3D (Convolutional Approach to Structure Identification-3D) to identify protostellar outflows in molecular line spectra. We conduct magneto-hydrodynamics simulations that model forming stars that launch protostellar outflows and use these to generate synthetic observations. We apply the 3D radiation transfer code RADMC-3D to model 12CO (J=1-0) line emission fr… ▽ More We adopt the deep learning method CASI-3D (Convolutional Approach to Structure Identification-3D) to identify protostellar outflows in molecular line spectra. We conduct magneto-hydrodynamics simulations that model forming stars that launch protostellar outflows and use these to generate synthetic observations. We apply the 3D radiation transfer code RADMC-3D to model 12CO (J=1-0) line emission from the simulated clouds. We train two CASI-3D models: ME1 is trained to predict only the position of outflows, while MF is trained to predict the fraction of the mass coming from outflows in each voxel. The two models successfully identify all 60 previously visually identified outflows in Perseus. Additionally, CASI-3D finds 20 new high-confidence outflows. All of these have coherent high-velocity structures, and 17 of them have nearby young stellar objects, while the remaining three are outside the Spitzer survey coverage. The mass, momentum and energy of individual outflows in Perseus predicted by model MF is comparable to the previous estimations. This similarity is due to a cancelation in errors: previous calculations missed outflow material with velocities comparable to the cloud velocity, however, they compensate for this by over-estimating the amount of mass at higher velocities that has contamination from non-outflow gas. We show outflows likely driven by older sources have more high-velocity gas compared to those driven by younger sources. △ Less

Submitted 23 October, 2020; originally announced October 2020.

Comments: Revised version submitted to ApJ, comments welcome

STARFORGE: Toward a comprehensive numerical model of star cluster formation and feedback

Authors: Michael Y. Grudić, Dávid Guszejnov, Philip F. Hopkins, Stella S. R. Offner, Claude-André Faucher-Giguère

Abstract: We present STARFORGE (STAR FORmation in Gaseous Environments): a new numerical framework for 3D radiation MHD simulations of star formation that simultaneously follow the formation, accretion, evolution, and dynamics of individual stars in massive giant molecular clouds (GMCs) while accounting for stellar feedback, including jets, radiative heating and momentum, stellar winds, and supernovae. We u… ▽ More We present STARFORGE (STAR FORmation in Gaseous Environments): a new numerical framework for 3D radiation MHD simulations of star formation that simultaneously follow the formation, accretion, evolution, and dynamics of individual stars in massive giant molecular clouds (GMCs) while accounting for stellar feedback, including jets, radiative heating and momentum, stellar winds, and supernovae. We use the GIZMO code with the MFM mesh-free Lagrangian MHD method, augmented with new algorithms for gravity, timestepping, sink particle formation and accretion, stellar dynamics, and feedback coupling. We survey a wide range of numerical parameters/prescriptions for sink formation and accretion and find very small variations in star formation history and the IMF (except for intentionally-unphysical variations). Modules for mass-injecting feedback (winds, SNe, and jets) inject new gas elements on-the-fly, eliminating the lack of resolution in diffuse feedback cavities otherwise inherent in Lagrangian methods. The treatment of radiation uses GIZMO's radiative transfer solver to track 5 frequency bands (IR, optical, NUV, FUV, ionizing), coupling direct stellar emission and dust emission with gas heating and radiation pressure terms. We demonstrate accurate solutions for SNe, winds, and radiation in problems with known similarity solutions, and show that our jet module is robust to resolution and numerical details, and agrees well with previous AMR simulations. STARFORGE can scale up to massive ($>10^5 M_\odot $) GMCs on current supercomputers while predicting the stellar ($\gtrsim 0.1 M_\odot$) range of the IMF, permitting simulations of both high- and low-mass cluster formation in a wide range of conditions. △ Less

Submitted 27 February, 2024; v1 submitted 21 October, 2020; originally announced October 2020.

Comments: Fix error in Eq. 44 for wind mass loss rate

STARFORGE: The effects of protostellar outflows on the IMF

Authors: Dávid Guszejnov, Michael Y. Grudić, Philip F. Hopkins, Stella S. R. Offner, Claude-André Faucher-Giguère

Abstract: The initial mass function (IMF) of stars is a key quantity affecting almost every field of astrophysics, yet it remains unclear what physical mechanisms determine it. We present the first runs of the STARFORGE project, using a new numerical framework to follow the formation of individual stars in giant molecular clouds (GMCs) using the GIZMO code. Our suite include runs with increasingly complex p… ▽ More The initial mass function (IMF) of stars is a key quantity affecting almost every field of astrophysics, yet it remains unclear what physical mechanisms determine it. We present the first runs of the STARFORGE project, using a new numerical framework to follow the formation of individual stars in giant molecular clouds (GMCs) using the GIZMO code. Our suite include runs with increasingly complex physics, starting with isothermal ideal magnetohydrodynamic (MHD) and then adding non-isothermal thermodynamics and protostellar outflows. We show that without protostellar outflows the resulting stellar masses are an order of magnitude too high, similar to the result in the base isothermal MHD run. Outflows disrupt the accretion flow around the protostar, allowing gas to fragment and additional stars to form, thereby lowering the mean stellar mass to a value similar to that observed. The effect of jets upon global cloud evolution is most pronounced for lower-mass GMCs and dense clumps, so while jets can disrupt low-mass clouds, they are unable to regulate star formation in massive GMCs, as they would turn an order unity fraction of the mass into stars before unbinding the cloud. Jets are also unable to stop the runaway accretion of massive stars, which could ultimately lead to the formation of stars with masses $\mathrm{>500\,M_\odot}$. Although we find that the mass scale set by jets is insensitive to most cloud parameters (i.e., surface density, virial parameter), it is strongly dependent on the momentum loading of the jets (which is poorly constrained by observations) as well the the temperature of the parent cloud, which predicts slightly larger IMF variations than observed. We conclude that protostellar jets play a vital role in setting the mass scale of stars, but additional physics are necessary to reproduce the observed IMF. △ Less

Submitted 22 January, 2021; v1 submitted 21 October, 2020; originally announced October 2020.

Comments: 18 pages, 13 figures, submitted to MNRAS

Ubiquitous $\rm NH_3$ supersonic component in L1688 coherent cores

Authors: Spandan Choudhury, Jaime E. Pineda, Paola Caselli, Adam Ginsburg, Stella S. R. Offner, Erik Rosolowsky, Rachel K. Friesen, Felipe O. Alves, Ana Chacón-Tanarro, Anna Punanova, Elena Redaelli, Helen Kirk, Philip C. Myers, Peter G. Martin, Yancy Shirley, Michael Chun-Yuan Chen, Alyssa A. Goodman, James Di Francesco

Abstract: Context : Star formation takes place in cold dense cores in molecular clouds. Earlier observations have found that dense cores exhibit subsonic non-thermal velocity dispersions. In contrast, CO observations show that the ambient large-scale cloud is warmer and has supersonic velocity dispersions. Aims : We aim to study the ammonia ($\rm NH_3$) molecular line profiles with exquisite sensitivity tow… ▽ More Context : Star formation takes place in cold dense cores in molecular clouds. Earlier observations have found that dense cores exhibit subsonic non-thermal velocity dispersions. In contrast, CO observations show that the ambient large-scale cloud is warmer and has supersonic velocity dispersions. Aims : We aim to study the ammonia ($\rm NH_3$) molecular line profiles with exquisite sensitivity towards the coherent cores in L1688 in order to study their kinematical properties in unprecedented detail. Methods : We used $\rm NH_3$ (1,1) and (2,2) data from the first data release (DR1) in the Green Bank Ammonia Survey (GAS). We first smoothed the data to a larger beam of 1' to obtain substantially more extended maps of velocity dispersion and kinetic temperature, compared to the DR1 maps. We then identified the coherent cores in the cloud and analysed the averaged line profiles towards the cores. Results : For the first time, we detected a faint (mean $\rm NH_3$(1,1) peak brightness $<$0.25 K in $T_{MB}$), supersonic component towards all the coherent cores in L1688. We fitted two components, one broad and one narrow, and derived the kinetic temperature and velocity dispersion of each component. The broad components towards all cores have supersonic linewidths ($\mathcal{M}_S \ge 1$). This component biases the estimate of the narrow dense core component's velocity dispersion by $\approx$28% and the kinetic temperature by $\approx$10%, on average, as compared to the results from single-component fits. Conclusions : Neglecting this ubiquitous presence of a broad component towards all coherent cores causes the typical single-component fit to overestimate the temperature and velocity dispersion. This affects the derived detailed physical structure and stability of the cores estimated from $\rm NH_3$ observations. △ Less

Submitted 20 July, 2020; v1 submitted 14 July, 2020; originally announced July 2020.

Comments: Accepted for publication in Astronomy & Astrophysics on 06/07/2020. 15 pages, 16 figures, 1 table. Language edits from previous version

Journal ref: A&A 640, L6 (2020)

Turbulence, Coherence and Collapse: Three Phases for Core Evolution

Authors: Stella S. R. Offner, Josh Taylor, Carleen Markey, Hope How-Huan Chen, Jaime E. Pineda, Alyssa A. Goodman, Andreas Burkert, Adam Ginsburg, Spandan Choudhury

Abstract: We study the formation, evolution and collapse of dense cores by tracking structures in a magnetohydrodynamic simulation of a star-forming cloud. We identify cores using the dendrogram algorithm and utilize machine learning techniques, including Neural Gas prototype learning and Fuzzy $c$-means clustering, to analyze the density and velocity dispersion profiles of cores together with six bulk prop… ▽ More We study the formation, evolution and collapse of dense cores by tracking structures in a magnetohydrodynamic simulation of a star-forming cloud. We identify cores using the dendrogram algorithm and utilize machine learning techniques, including Neural Gas prototype learning and Fuzzy $c$-means clustering, to analyze the density and velocity dispersion profiles of cores together with six bulk properties. We produce a 2-d visualization using a Uniform Manifold Approximation and Projection (UMAP), which facilitates the connection between physical properties and three partially-overlapping phases: i) unbound turbulent structures (Phase I), ii) coherent cores that have low turbulence (Phase II), and iii) bound cores, many of which become protostellar (Phase III). Within Phase II we identify a population of long-lived coherent cores that reach a quasi-equilibrium state. Most prestellar cores form in Phase II and become protostellar after evolving into Phase III. Due to the turbulent cloud environment, the initial core properties do not uniquely predict the eventual evolution, i.e., core evolution is stochastic, and cores follow no one evolutionary path. The phase lifetimes are 1.0$\pm$0.1$\times$10$^5$ yr, 1.3$\pm$0.2$\times$10$^5$ yr, and 1.8$\pm$0.3$\times$10$^5$ yr for Phase I, II, and III, respectively. We compare our results to NH$_3$ observations of dense cores. Known coherent cores predominantly map into Phase II, while most turbulent pressure-confined cores map to Phase I or III. We predict that a significant fraction of observed starless cores have unresolved coherent regions and that $\gtrsim 20$% of observed starless cores will not form stars. Measurements of core radial profiles, in addition to the usual bulk properties, will enable more accurate predictions of core evolution. △ Less

Submitted 26 September, 2022; v1 submitted 12 June, 2020; originally announced June 2020.

Comments: 26 pages, 16 figures. Accepted to MNRAS

The Origin of the Stellar Mass Distribution and Multiplicity

Authors: Yueh-Ning Lee, Stella S. R. Offner, Patrick Hennebelle, Philippe André, Hans Zinnecker, Javier Ballesteros-Paredes, Shu-ichiro Inutsuka, J. M. Diederik Kruijssen

Abstract: In this chapter, we review some historical understanding and recent advances on the Initial Mass Function (IMF) and the Core Mass Function (CMF), both in terms of observations and theories. We focus mostly on star formation in clustered environment since this is suggested by observations to be the dominant mode of star formation. The statistical properties and the fragmentation behaviour of turbul… ▽ More In this chapter, we review some historical understanding and recent advances on the Initial Mass Function (IMF) and the Core Mass Function (CMF), both in terms of observations and theories. We focus mostly on star formation in clustered environment since this is suggested by observations to be the dominant mode of star formation. The statistical properties and the fragmentation behaviour of turbulent gas is discussed, and we also discuss the formation of binaries and small multiple systems. △ Less

Submitted 10 June, 2020; originally announced June 2020.

Comments: Review chapter published in Space Science Reviews

Zooming in on Individual Star Formation: Low- and High-mass Stars

Authors: Anna L. Rosen, Stella S. R. Offner, Sarah I. Sadavoy, Asmita Bhandare, Enrique Vázquez-Semadeni, Adam Ginsburg

Abstract: Star formation is a multi-scale, multi-physics problem ranging from the size scale of molecular clouds ($\sim$10s pc) down to the size scales of dense prestellar cores ($\sim$0.1 pc) that are the birth sites of stars. Several physical processes like turbulence, magnetic fields and stellar feedback, such as radiation pressure and outflows, are more or less important for different stellar masses and… ▽ More Star formation is a multi-scale, multi-physics problem ranging from the size scale of molecular clouds ($\sim$10s pc) down to the size scales of dense prestellar cores ($\sim$0.1 pc) that are the birth sites of stars. Several physical processes like turbulence, magnetic fields and stellar feedback, such as radiation pressure and outflows, are more or less important for different stellar masses and size scales. During the last decade a variety of technological and computing advances have transformed our understanding of star formation through the use of multi-wavelength observations, large scale observational surveys, and multi-physics multi-dimensional numerical simulations. Additionally, the use of synthetic observations of simulations have provided a useful tool to interpret observational data and evaluate the importance of various physical processes on different scales in star formation. Here, we review these recent advancements in both high- ($M \gtrsim 8 \, M_{\rm \odot}$) and low-mass star formation. △ Less

Submitted 15 May, 2020; originally announced May 2020.

Comments: 53 pages, 7 figures, To appear in Space Science Reviews (submitted January 31, 2020; accepted May 7, 2020), topical collection on Star formation

Physical Processes in Star Formation

Authors: Philipp Girichidis, Stella S. R. Offner, Alexei G. Kritsuk, Ralf S. Klessen, Patrick Hennebelle, J. M. Diederik Kruijssen, Martin G. H. Krause, Simon C. O. Glover, Marco Padovani

Abstract: Star formation is a complex multi-scale phenomenon that is of significant importance for astrophysics in general. Stars and star formation are key pillars in observational astronomy from local star forming regions in the Milky Way up to high-redshift galaxies. From a theoretical perspective, star formation and feedback processes (radiation, winds, and supernovae) play a pivotal role in advancing o… ▽ More Star formation is a complex multi-scale phenomenon that is of significant importance for astrophysics in general. Stars and star formation are key pillars in observational astronomy from local star forming regions in the Milky Way up to high-redshift galaxies. From a theoretical perspective, star formation and feedback processes (radiation, winds, and supernovae) play a pivotal role in advancing our understanding of the physical processes at work, both individually and of their interactions. In this review we will give an overview of the main processes that are important for the understanding of star formation. We start with an observationally motivated view on star formation from a global perspective and outline the general paradigm of the life-cycle of molecular clouds, in which star formation is the key process to close the cycle. After that we focus on the thermal and chemical aspects in star forming regions, discuss turbulence and magnetic fields as well as gravitational forces. Finally, we review the most important stellar feedback mechanisms. △ Less

Submitted 13 May, 2020; originally announced May 2020.

Comments: 86 pages, 4 figures. To appear in Space Science Reviews, topical collection Star formation

Star-Gas Surface Density Correlations in Twelve Nearby Molecular Clouds I: Data Collection and Star-Sampled Analysis

Authors: Riwaj Pokhrel, Robert A. Gutermuth, Sarah K. Betti, Stella S. R. Offner, Philip C. Myers, S. Thomas Megeath, Alyssa D. Sokol, Babar Ali, Lori Allen, Tom S. Allen, Michael M. Dunham, William J. Fischer, Thomas Henning, Mark Heyer, Joseph L. Hora, Judith L. Pipher, John J. Tobin, Scott J. Wolk

Abstract: We explore the relation between the stellar mass surface density and the mass surface density of molecular hydrogen gas in twelve nearby molecular clouds that are located at $<$1.5 kpc distance. The sample clouds span an order of magnitude range in mass, size, and star formation rates. We use thermal dust emission from $Herschel$ maps to probe the gas surface density and the young stellar objects… ▽ More We explore the relation between the stellar mass surface density and the mass surface density of molecular hydrogen gas in twelve nearby molecular clouds that are located at $<$1.5 kpc distance. The sample clouds span an order of magnitude range in mass, size, and star formation rates. We use thermal dust emission from $Herschel$ maps to probe the gas surface density and the young stellar objects from the most recent $Spitzer$ Extended Solar Neighborhood Archive (SESNA) catalog to probe the stellar surface density. Using a star-sampled nearest neighbor technique to probe the star-gas surface density correlations at the scale of a few parsecs, we find that the stellar mass surface density varies as a power-law of the gas mass surface density, with a power-law index of $\sim$2 in all the clouds. The consistent power-law index implies that star formation efficiency is directly correlated with gas column density, and no gas column density threshold for star formation is observed. We compare the observed correlations with the predictions from an analytical model of thermal fragmentation, and with the synthetic observations of a recent hydrodynamic simulation of a turbulent star-forming molecular cloud. We find that the observed correlations are consistent for some clouds with the thermal fragmentation model and can be reproduced using the hydrodynamic simulations. △ Less

Submitted 11 May, 2020; originally announced May 2020.

Comments: 39 pages, 10 figures, 3 tables. Accepted in ApJ

The Physics of Star Cluster Formation and Evolution

Authors: Martin G. H. Krause, Stella S. R. Offner, Corinne Charbonnel, Mark Gieles, Ralf S. Klessen, Enrique Vazquez-Semadeni, Javier Ballesteros-Paredes, Philipp Girichidis, J. M. Diederik Kruijssen, Jacob L. Ward, Hans Zinnecker

Abstract: Star clusters form in dense, hierarchically collapsing gas clouds. Bulk kinetic energy is transformed to turbulence with stars forming from cores fed by filaments. In the most compact regions, stellar feedback is least effective in removing the gas and stars may form very efficiently. These are also the regions where, in high-mass clusters, ejecta from some kind of high-mass stars are effectively… ▽ More Star clusters form in dense, hierarchically collapsing gas clouds. Bulk kinetic energy is transformed to turbulence with stars forming from cores fed by filaments. In the most compact regions, stellar feedback is least effective in removing the gas and stars may form very efficiently. These are also the regions where, in high-mass clusters, ejecta from some kind of high-mass stars are effectively captured during the formation phase of some of the low mass stars and effectively channeled into the latter to form multiple populations. Star formation epochs in star clusters are generally set by gas flows that determine the abundance of gas in the cluster. We argue that there is likely only one star formation epoch after which clusters remain essentially clear of gas by cluster winds. Collisional dynamics is important in this phase leading to core collapse, expansion and eventual dispersion of every cluster. We review recent developments in the field with a focus on theoretical work. △ Less

Submitted 15 May, 2020; v1 submitted 2 May, 2020; originally announced May 2020.

Comments: 53 pages, 7 figures, preprint of a review to appear in Space Science Reviews Topical collection "Star formation". Changes in this version: updated comments line to to acceptance of the paper by the journal

The highly variable time evolution of star-forming cores identified with dendrograms

Authors: Rachel A. Smullen, Kaitlin M. Kratter, Stella S. R. Offner, Aaron T. Lee, Hope How-Huan Chen

Abstract: We investigate the time evolution of dense cores identified in molecular cloud simulations using dendrograms, which are a common tool to identify hierarchical structure in simulations and observations of star formation. We develop an algorithm to link dendrogram structures through time using the three-dimensional density field from magnetohydrodynamical simulations, thus creating histories for all… ▽ More We investigate the time evolution of dense cores identified in molecular cloud simulations using dendrograms, which are a common tool to identify hierarchical structure in simulations and observations of star formation. We develop an algorithm to link dendrogram structures through time using the three-dimensional density field from magnetohydrodynamical simulations, thus creating histories for all dense cores in the domain. We find that the population-wide distributions of core properties are relatively invariant in time, and quantities like the core mass function match with observations. Despite this consistency, an individual core may undergo large (>40%), stochastic variations due to the redefinition of the dendrogram structure between timesteps. This variation occurs independent of environment and stellar content. We identify a population of short-lived (<200 kyr) overdensities masquerading as dense cores that may comprise ~20% of any time snapshot. Finally, we note the importance of considering the full history of cores when interpreting the origin of the initial mass function; we find that, especially for systems containing multiple stars, the core mass defined by a dendrogram leaf in a snapshot is typically less than the final system stellar mass. This work reinforces that there is no time-stable density contour that defines a star-forming core. The dendrogram itself can induce significant structure variation between timesteps due to small changes in the density field. Thus, one must use caution when comparing dendrograms of regions with different ages or environment properties because differences in dendrogram structure may not come solely from the physical evolution of dense cores. △ Less

Submitted 2 April, 2020; originally announced April 2020.

Comments: 20 pages, 17 figures. Submitted to MNRAS

Relative Alignment between Dense Molecular Cores and Ambient Magnetic Field: The Synergy of Numerical Models and Observations

Authors: Che-Yu Chen, Erica A. Behrens, Jasmin E. Washington, Laura M. Fissel, Rachel K. Friesen, Zhi-Yun Li, Jaime E. Pineda, Adam Ginsburg, Helen Kirk, Samantha Scibelli, Felipe Alves, Elena Redaelli, Paola Caselli, Anna Punanova, James Di Francesco, Erik Rosolowsky, Stella S. R. Offner, Peter G. Martin, Ana Chacón-Tanarro, Hope H. -H. Chen, Michael C. -Y. Chen, Jared Keown, Youngmin Seo, Yancy Shirley, Hector G. Arce , et al. (4 additional authors not shown)

Abstract: The role played by magnetic field during star formation is an important topic in astrophysics. We investigate the correlation between the orientation of star-forming cores (as defined by the core major axes) and ambient magnetic field directions in 1) a 3D MHD simulation, 2) synthetic observations generated from the simulation at different viewing angles, and 3) observations of nearby molecular cl… ▽ More The role played by magnetic field during star formation is an important topic in astrophysics. We investigate the correlation between the orientation of star-forming cores (as defined by the core major axes) and ambient magnetic field directions in 1) a 3D MHD simulation, 2) synthetic observations generated from the simulation at different viewing angles, and 3) observations of nearby molecular clouds. We find that the results on relative alignment between cores and background magnetic field in synthetic observations slightly disagree with those measured in fully 3D simulation data, which is partly because cores identified in projected 2D maps tend to coexist within filamentary structures, while 3D cores are generally more rounded. In addition, we examine the progression of magnetic field from pc- to core-scale in the simulation, which is consistent with the anisotropic core formation model that gas preferably flow along the magnetic field toward dense cores. When comparing the observed cores identified from the GBT Ammonia Survey (GAS) and Planck polarization-inferred magnetic field orientations, we find that the relative core-field alignment has a regional dependence among different clouds. More specifically, we find that dense cores in the Taurus molecular cloud tend to align perpendicular to the background magnetic field, while those in Perseus and Ophiuchus tend to have random (Perseus) or slightly parallel (Ophiuchus) orientations with respect to the field. We argue that this feature of relative core-field orientation could be used to probe the relative significance of the magnetic field within the cloud. △ Less

Submitted 24 March, 2020; originally announced March 2020.

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

Impact of low-energy cosmic rays on star formation

Authors: Marco Padovani, Alexei V. Ivlev, Daniele Galli, Stella S. R. Offner, Nick Indriolo, Donna Rodgers-Lee, Alexandre Marcowith, Philipp Girichidis, Andrei M. Bykov, J. M. Diederik Kruijssen

Abstract: In recent years, exciting developments have taken place in the identification of the role of cosmic rays in star-forming environments. Observations from radio to infrared wavelengths and theoretical modelling have shown that low-energy cosmic rays (<1 TeV) play a fundamental role in shaping the chemical richness of the interstellar medium, determining the dynamical evolution of molecular clouds. I… ▽ More In recent years, exciting developments have taken place in the identification of the role of cosmic rays in star-forming environments. Observations from radio to infrared wavelengths and theoretical modelling have shown that low-energy cosmic rays (<1 TeV) play a fundamental role in shaping the chemical richness of the interstellar medium, determining the dynamical evolution of molecular clouds. In this review we summarise in a coherent picture the main results obtained by observations and by theoretical models of propagation and generation of cosmic rays, from the smallest scales of protostars and circumstellar discs, to young stellar clusters, up to Galactic and extragalactic scales. We also discuss the new fields that will be explored in the near future thanks to new generation instruments, such as: CTA, for the $γ$-ray emission from high-mass protostars; SKA and precursors, for the synchrotron emission at different scales; and ELT/HIRES, JWST, and ARIEL, for the impact of cosmic rays on exoplanetary atmospheres and habitability. △ Less

Submitted 24 February, 2020; originally announced February 2020.

Comments: Submitted to Space Science Reviews topical collection Star formation

Can magnetized turbulence set the mass scale of stars?

Authors: David Guszejnov, Michael Y. Grudić, Philip F. Hopkins, Stella S. R. Offner, Claude-André Faucher-Giguère

Abstract: Understanding the evolution of self-gravitating, isothermal, magnetized gas is crucial for star formation, as these physical processes have been postulated to set the initial mass function (IMF). We present a suite of isothermal magnetohydrodynamic (MHD) simulations using the GIZMO code, that resolve the formation of individual stars in giant molecular clouds (GMCs), spanning a range of Mach numbe… ▽ More Understanding the evolution of self-gravitating, isothermal, magnetized gas is crucial for star formation, as these physical processes have been postulated to set the initial mass function (IMF). We present a suite of isothermal magnetohydrodynamic (MHD) simulations using the GIZMO code, that resolve the formation of individual stars in giant molecular clouds (GMCs), spanning a range of Mach numbers found in observed GMCs. As in past works, the mean and median stellar masses are sensitive to numerical resolution, because they are sensitive to low-mass stars that contribute a vanishing fraction of the overall stellar mass. The {\em mass-weighted} median stellar mass $M_\mathrm{50}$ becomes insensitive to resolution once turbulent fragmentation is well-resolved. Without imposing Larson-like scaling laws, our simulations find $M_\mathrm{50} \propto M_\mathrm{0} \mathcal{M}^{-3} α_\mathrm{turb} \mathrm{SFE}^{1/3}$ for GMC mass $M_\mathrm{0}$, sonic Mach number $\mathcal{M}$, virial parameter $α_\mathrm{turb}$, and star formation efficiency $\mathrm{SFE}=M_\mathrm{\star}/M_\mathrm{0}$. This fit agrees well with previous IMF results from the RAMSES, ORION2, and SphNG codes. Although $M_\mathrm{50}$ has no significant dependence on the magnetic field strength at the cloud scale, MHD is necessary to prevent a fragmentation cascade that results in non-convergent stellar masses. For initial conditions and SFE similar to star-forming GMCs in our Galaxy, we predict $M_\mathrm{50}$ to be $>20 M_{\odot}$, an order of magnitude larger than observed ($\sim 2 M_\odot$), together with an excess of brown dwarfs. Moreover, $M_\mathrm{50}$ is sensitive to initial cloud properties and evolves strongly in time within a given cloud, predicting much larger IMF variations than are observationally allowed. We conclude that physics beyond MHD turbulence and gravity are necessary ingredients for the IMF. △ Less

Submitted 26 October, 2020; v1 submitted 4 February, 2020; originally announced February 2020.

Comments: Note that this version has erratum at the end (does not affect results) 18 pages, 17 figures

Application of Convolutional Neural Networks to Identify Stellar Feedback Bubbles in CO Emission

Authors: Duo Xu, Stella S. R. Offner, Robert Gutermuth, Colin Van Oort

Abstract: We adopt the deep learning method CASI (Convolutional Approach to Shell Identification) and extend it to 3D (CASI-3D) to identify signatures of stellar feedback in molecular line spectra, such as 13CO. We adopt magneto-hydrodynamics simulations that study the impact of stellar winds in a turbulent molecular cloud as an input to generate synthetic observations. We apply the 3D radiation transfer co… ▽ More We adopt the deep learning method CASI (Convolutional Approach to Shell Identification) and extend it to 3D (CASI-3D) to identify signatures of stellar feedback in molecular line spectra, such as 13CO. We adopt magneto-hydrodynamics simulations that study the impact of stellar winds in a turbulent molecular cloud as an input to generate synthetic observations. We apply the 3D radiation transfer code radmc-3d to model 13CO (J=1-0) line emission from the simulated clouds. We train two CASI-3d models: ME1 predicts only the position of feedback, while MF predicts the fraction of the mass coming from feedback in each voxel. We adopt 75% of the synthetic observations as the training set and assess the accuracy of the two models with the remaining data. We demonstrate that model ME1 identifies bubbles in simulated data with 95% accuracy, and model MF predicts the bubble mass within 4% of the true value. We use bubbles previously visually identified in Taurus in 13CO to validate the models and show both perform well on the highest confidence bubbles. We apply our two models on the full 98 deg2 FCRAO 13CO survey of the Taurus cloud. Models ME1 and MF predict feedback gas mass of 2894 M and 302 M, respectively. When including a correction factor for missing energy due to the limited velocity range of the 13CO data cube, model ME1 predicts feedback kinetic energies of 4.0*1e46 ergs and 1.5*1e47 ergs with/without subtracting the cloud velocity gradient. Model MF predicts feedback kinetic energy of 9.6*1e45 ergs and 2.8*1e46 ergs with/without subtracting the cloud velocity gradient. Model ME1 predicts bubble locations and properties consistent with previous visually identified bubbles. However, model MF demonstrates that feedback properties computed based on visual identifications are significantly overestimated due to line of sight confusion and contamination from background and foreground gas. △ Less

Submitted 13 January, 2020; originally announced January 2020.

Comments: ApJ Accepted

The VLA/ALMA Nascent Disk and Multiplicity (VANDAM) Survey of Orion Protostars. A Statistical Characterization of Class 0 and I Protostellar Disks

Authors: John J. Tobin, Patrick Sheehan, S. Thomas Megeath, Ana Karla Diaz-Rodriguez, Stella S. R. Offner, Nadia M. Murillo, Merel van 't Hoff, Ewine F. van Dishoeck, Mayra Osorio, Guillem Anglada, Elise Furlan, Amelia M. Stutz, Nickalas Reynolds, Nicole Karnath, William J. Fischer, Magnus Persson, Leslie W. Looney, Zhi-Yun Li, Ian Stephens, Claire J. Chandler, Erin Cox, Michael M. Dunham, Lukasz Tychoniec, Mihkel Kama, Kaitlin Kratter , et al. (11 additional authors not shown)

Abstract: We have conducted a survey of 328 protostars in the Orion molecular clouds with ALMA at 0.87 mm at a resolution of $\sim$0.1" (40 au), including observations with the VLA at 9 mm toward 148 protostars at a resolution of $\sim$0.08" (32 au). This is the largest multi-wavelength survey of protostars at this resolution by an order of magnitude. We use the dust continuum emission at 0.87 mm and 9 mm t… ▽ More We have conducted a survey of 328 protostars in the Orion molecular clouds with ALMA at 0.87 mm at a resolution of $\sim$0.1" (40 au), including observations with the VLA at 9 mm toward 148 protostars at a resolution of $\sim$0.08" (32 au). This is the largest multi-wavelength survey of protostars at this resolution by an order of magnitude. We use the dust continuum emission at 0.87 mm and 9 mm to measure the dust disk radii and masses toward the Class 0, Class I, and Flat Spectrum protostars, characterizing the evolution of these disk properties in the protostellar phase. The mean dust disk radii for the Class 0, Class I, and Flat Spectrum protostars are 44.9$^{+5.8}_{-3.4}$, 37.0$^{+4.9}_{-3.0}$, and 28.5$^{+3.7}_{-2.3}$ au, respectively, and the mean protostellar dust disk masses are 25.9$^{+7.7}_{-4.0}$, 14.9$^{+3.8}_{-2.2}$, 11.6$^{+3.5}_{-1.9}$ Earth masses, respectively. The decrease in dust disk masses is expected from disk evolution and accretion, but the decrease in disk radii may point to the initial conditions of star formation not leading to the systematic growth of disk radii or that radial drift is keeping the dust disk sizes small. At least 146 protostellar disks (35% out of 379 detected 0.87 mm continuum sources plus 42 non-detections) have disk radii greater than 50 au in our sample. These properties are not found to vary significantly between different regions within Orion. The protostellar dust disk mass distributions are systematically larger than that of Class II disks by a factor of $>$4, providing evidence that the cores of giant planets may need to at least begin their formation during the protostellar phase. △ Less

Submitted 13 January, 2020; originally announced January 2020.

Comments: 107 pages, 21 Figures, 11 Tables, accepted to ApJ. Version with all source figures: https://www.cv.nrao.edu/~jtobin/Orion-disks-fullfigs.pdf Reduced data available from https://dataverse.harvard.edu/dataverse/VANDAMOrion

Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) -- Full Data Release

Authors: Ian W. Stephens, Tyler L. Bourke, Michael M. Dunham, Philip C. Myers, Riwaj Pokhrel, John J. Tobin, Héctor G. Arce, Sarah I. Sadavoy, Eduard I. Vorobyov, Jaime E. Pineda, Stella S. R. Offner, Katherine I. Lee, Lars E. Kristensen, Jes K. Jørgensen, Mark A. Gurwell, Alyssa A. Goodman

Abstract: We present and release the full dataset for the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. This survey used the Submillimeter Array (SMA) to image the 74 known protostars within the Perseus molecular cloud. The SMA was used in two array configurations to capture outflows for scales $>$30$^{\prime\prime}$ ($>$9000 au) and to probe scales down to $\sim$1… ▽ More We present and release the full dataset for the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. This survey used the Submillimeter Array (SMA) to image the 74 known protostars within the Perseus molecular cloud. The SMA was used in two array configurations to capture outflows for scales $>$30$^{\prime\prime}$ ($>$9000 au) and to probe scales down to $\sim$1$^{\prime\prime}$ ($\sim$300 au). The protostars were observed with the 1.3 mm and 850 $μ$m receivers simultaneously to detect continuum at both wavelengths and molecular line emission from CO(2-1), $^{13}$CO(2-1), C$^{18}$O(2-1), N$_2$D$^+$(3-2), CO(3-2), HCO$^+$(4-3), and H$^{13}$CO$^+$(4-3). Some of the observations also used the SMA's recently upgraded correlator, SWARM, whose broader bandwidth allowed for several more spectral lines to be observed (e.g., SO, H$_2$CO, DCO$^+$, DCN, CS, CN). Of the main continuum and spectral tracers observed, 84% of the images and cubes had emission detected. The median C$^{18}$O(2-1) linewidth is $\sim$1.0 km s$^{-1}$, which is slightly higher than those measured with single-dish telescopes at scales of 3000-20000 au. Of the 74 targets, six are suggested to be first hydrostatic core candidates, and we suggest that L1451-mm is the best candidate. We question a previous continuum detection toward L1448 IRS2E. In the SVS13 system, SVS13A certainly appears to be the most evolved source, while SVS13C appears to be hotter and more evolved than SVS13B. The MASSES survey is the largest publicly available interferometric continuum and spectral line protostellar survey to date, and is largely unbiased as it only targets protostars in Perseus. All visibility ($uv$) data and imaged data are publicly available at https://dataverse.harvard.edu/dataverse/full_MASSES/. △ Less

Submitted 19 November, 2019; originally announced November 2019.

Comments: Accepted to ApJS

The Formation and Evolution of Wide-Orbit Stellar Multiples In Magnetized Clouds

Authors: Aaron T. Lee, Stella S. R. Offner, Kaitlin M. Kratter, Rachel A. Smullen, Pak Shing Li

Abstract: Stars rarely form in isolation. Nearly half of the stars in the Milky Way have a companion, and this fraction increases in star-forming regions. However, why some dense cores and filaments form bound pairs while others form single stars remains unclear. We present a set of three-dimensional, gravo-magnetohydrodynamic simulations of turbulent star-forming clouds, aimed at understanding the formatio… ▽ More Stars rarely form in isolation. Nearly half of the stars in the Milky Way have a companion, and this fraction increases in star-forming regions. However, why some dense cores and filaments form bound pairs while others form single stars remains unclear. We present a set of three-dimensional, gravo-magnetohydrodynamic simulations of turbulent star-forming clouds, aimed at understanding the formation and evolution of multiple-star systems formed through large scale (>~$10^3$ AU) turbulent fragmentation. We investigate three global magnetic field strengths, with global mass-to-flux ratios of $μ_φ$=2, 8, and 32. The initial separations of protostars in multiples depends on the global magnetic field strength, with stronger magnetic fields (e.g., $μ_φ$=2) suppressing fragmentation on smaller scales. The overall multiplicity fraction (MF) is between 0.4-0.6 for our strong and intermediate magnetic field strengths, which is in agreement with observations. The weak field case has a lower fraction. The MF is relatively constant throughout the simulations, even though stellar densities increase as collapse continues. While the MF rarely exceeds 60% in all three simulations, over 80% of all protostars are part of a binary system at some point. We additionally find that the distribution of binary spin mis-alignment angles is consistent with a randomized distribution. In all three simulations, several binaries originate with wide separations and dynamically evolve to <~ $10^2$ AU separations. We show that a simple model of mass accretion and dynamical friction with the gas can explain this orbital evolution. △ Less

Submitted 18 November, 2019; originally announced November 2019.

Comments: Accepted to ApJ. 28 pages, 23 figures, comments from the community welcomed

Evolution of giant molecular clouds across cosmic time

Authors: David Guszejnov, Michael Y. Grudić, Stella S. R. Offner, Michael Boylan-Kolchin, Claude-André Faucher-Giguère, Andrew Wetzel, Samantha M. Benincasa, Sarah Loebman

Abstract: Giant molecular clouds (GMCs) are well-studied in the local Universe, however, exactly how their properties vary during galaxy evolution is poorly understood due to challenging resolution requirements, both observational and computational. We present the first time-dependent analysis of giant molecular clouds in a Milky Way-like galaxy and an LMC-like dwarf galaxy of the FIRE-2 (Feedback In Realis… ▽ More Giant molecular clouds (GMCs) are well-studied in the local Universe, however, exactly how their properties vary during galaxy evolution is poorly understood due to challenging resolution requirements, both observational and computational. We present the first time-dependent analysis of giant molecular clouds in a Milky Way-like galaxy and an LMC-like dwarf galaxy of the FIRE-2 (Feedback In Realistic Environments) simulation suite, which have sufficient resolution to predict the bulk properties of GMCs in cosmological galaxy formation self-consistently. We show explicitly that the majority of star formation outside the galactic center occurs within self-gravitating gas structures that have properties consistent with observed bound GMCs. We find that the typical cloud bulk properties such as mass and surface density do not vary more than a factor of 2 in any systematic way after the first Gyr of cosmic evolution within a given galaxy from its progenitor. While the median properties are constant, the tails of the distributions can briefly undergo drastic changes, which can produce very massive and dense self-gravitating gas clouds. Once the galaxy forms, we identify only two systematic trends in bulk properties over cosmic time: a steady increase in metallicity produced by previous stellar populations and a weak decrease in bulk cloud temperatures. With the exception of metallicity we find no significant differences in cloud properties between the Milky Way-like and dwarf galaxies. These results have important implications for cosmological star and star cluster formation and put especially strong constraints on theories relating the stellar initial mass function to cloud properties. △ Less

Submitted 24 February, 2020; v1 submitted 2 October, 2019; originally announced October 2019.

Comments: 15 pages, 11 figures accepted by MNRAS

The VLA/ALMA Nascent Disk and Multiplicity (VANDAM) Survey of Orion Protostars I. Identifying and Characterizing the Protostellar Content of the OMC2-FIR4 and OMC2-FIR3 Regions

Authors: John J. Tobin, S. Thomas Megeath, Merel van 't Hoff, Ana Karla Diaz-Rodriguez, Nickalas Reynolds, Mayra Osorio, Guillem Anglada, Elise Furlan, Nicole Karnath, Stella S. R. Offner, Patrick Sheehan, Sarah I. Sadavoy, Amelia M. Stutz, William J. Fischer, Mihkel Kama, Magnus Persson, James Di Francesco, Leslie W. Looney, Dan M. Watson, Zhi-Yun Li, Ian Stephens, Claire J. Chandler, Erin Cox, Michael M. Dunham, Kaitlin Kratter , et al. (9 additional authors not shown)

Abstract: We present ALMA (0.87~mm) and VLA (9~mm) observations toward OMC2-FIR4 and OMC2-FIR3 within the Orion integral-shaped filament that are thought to be the nearest regions of intermediate mass star formation. We characterize the continuum sources within these regions on $\sim$40~AU (0\farcs1) scales and associated molecular line emission at a factor of $\sim$30 better resolution than previous observ… ▽ More We present ALMA (0.87~mm) and VLA (9~mm) observations toward OMC2-FIR4 and OMC2-FIR3 within the Orion integral-shaped filament that are thought to be the nearest regions of intermediate mass star formation. We characterize the continuum sources within these regions on $\sim$40~AU (0\farcs1) scales and associated molecular line emission at a factor of $\sim$30 better resolution than previous observations at similar wavelengths. We identify six compact continuum sources within OMC2-FIR4, four in OMC2-FIR3, and one additional source just outside OMC2-FIR4. This continuum emission is tracing the inner envelope and/or disk emission on less than 100~AU scales. HOPS-108 is the only protostar in OMC2-FIR4 that exhibits emission from high-excitation transitions of complex organic molecules (e.g., methanol and other lines) coincident with the continuum emission. HOPS-370 in OMC2-FIR3 with L~$\sim$~360~\lsun, also exhibits emission from high-excitation methanol and other lines. The methanol emission toward these two protostars is indicative of temperatures high enough to thermally evaporate methanol from icy dust grains; overall these protostars have characteristics similar to hot corinos. We do not identify a clear outflow from HOPS-108 in \twco, but find evidence of interaction between the outflow/jet from HOPS-370 and the OMC2-FIR4 region. The multitude of observational constraints indicate that HOPS-108 is likely a low to intermediate-mass protostar in its main mass accretion phase and it is the most luminous protostar in OMC2-FIR4. The high resolution data presented here are essential for disentangling the embedded protostars from their surrounding dusty environments and characterizing them. △ Less

Submitted 1 October, 2019; originally announced October 2019.

Comments: 32 Pages, 10 Figures, 5 Tables, accepted to ApJ

KFPA Examinations of Young STellar Object Natal Environments (KEYSTONE): Hierarchical Ammonia Structures in Galactic Giant Molecular Clouds

Authors: Jared Keown, James Di Francesco, Erik Rosolowsky, Ayushi Singh, Charles Figura, Helen Kirk, L. D. Anderson, Michael Chun-Yuan Chen, Davide Elia, Rachel Friesen, Adam Ginsburg, A. Marston, Stefano Pezzuto, Eugenio Schisano, Sylvain Bontemps, Paola Caselli, Hong-Li Liu, Steven Longmore, Frederique Motte, Philip C. Myers, Stella S. R. Offner, Patricio Sanhueza, Nicola Schneider, Ian Stephens, James Urquhart , et al. (1 additional authors not shown)

Abstract: We present initial results from the K-band focal plane array Examinations of Young STellar Object Natal Environments (KEYSTONE) survey, a large project on the 100-m Green Bank Telescope mapping ammonia emission across eleven giant molecular clouds at distances of $0.9-3.0$ kpc (Cygnus X North, Cygnus X South, M16, M17, MonR1, MonR2, NGC2264, NGC7538, Rosette, W3, and W48). This data release includ… ▽ More We present initial results from the K-band focal plane array Examinations of Young STellar Object Natal Environments (KEYSTONE) survey, a large project on the 100-m Green Bank Telescope mapping ammonia emission across eleven giant molecular clouds at distances of $0.9-3.0$ kpc (Cygnus X North, Cygnus X South, M16, M17, MonR1, MonR2, NGC2264, NGC7538, Rosette, W3, and W48). This data release includes the NH$_3$ (1,1) and (2,2) maps for each cloud, which are modeled to produce maps of kinetic temperature, centroid velocity, velocity dispersion, and ammonia column density. Median cloud kinetic temperatures range from $11.4\pm2.2$ K in the coldest cloud (MonR1) to $23.0\pm6.5$ K in the warmest cloud (M17). Using dendrograms on the NH$_3$ (1,1) integrated intensity maps, we identify 856 dense gas clumps across the eleven clouds. Depending on the cloud observed, $40-100\%$ of the clumps are aligned spatially with filaments identified in H$_2$ column density maps derived from SED-fitting of dust continuum emission. A virial analysis reveals that 523 of the 835 clumps ($\sim63\%$) with mass estimates are bound by gravity alone. We find no significant difference between the virial parameter distributions for clumps aligned with the dust-continuum filaments and those unaligned with filaments. In some clouds, however, hubs or ridges of dense gas with unusually high mass and low virial parameters are located within a single filament or at the intersection of multiple filaments. These hubs and ridges tend to host water maser emission, multiple 70$μ$m-detected protostars, and have masses and radii above an empirical threshold for forming massive stars. △ Less

Submitted 29 August, 2019; v1 submitted 27 August, 2019; originally announced August 2019.

Comments: Accepted for publication in ApJ

Impact of Cosmic-Ray Feedback on Accretion and Chemistry in Circumstellar Disks

Authors: S. S. R. Offner, B. A. L. Gaches, J. Holdship

Abstract: We use the gas-grain chemistry code UCLCHEM to explore the impact of cosmic-ray feedback on the chemistry of circumstellar disks. We model the attenuation and energy losses of the cosmic-rays as they propagate outwards from the star and also consider ionization due to stellar radiation and radionuclides. For accretion rates typical of young stars, $\dot M_* \sim 10^{-9}-10^{-6}$ M_\odot yr$^{-1}$,… ▽ More We use the gas-grain chemistry code UCLCHEM to explore the impact of cosmic-ray feedback on the chemistry of circumstellar disks. We model the attenuation and energy losses of the cosmic-rays as they propagate outwards from the star and also consider ionization due to stellar radiation and radionuclides. For accretion rates typical of young stars, $\dot M_* \sim 10^{-9}-10^{-6}$ M_\odot yr$^{-1}$, we show that cosmic rays accelerated by the stellar accretion shock produce a cosmic-ray ionization rate at the disk surface $ζ\gtrsim 10^{-15}$ s$^{-1}$, at least an order of magnitude higher than the ionization rate associated with the Galactic cosmic-ray background. The incident cosmic-ray flux enhances the disk ionization at intermediate to high surface densities ($Σ> 10$ g cm$^{-2}$) particularly within 10 au of the star. We find the dominant ions are C$^+$, S$^+$ and Mg$^+$ in the disk surface layers, while the H$_3^+$ ion dominates at surface densities above 1.0 g cm$^{-2}$. We predict the radii and column densities at which the magneto-rotational instability (MRI) is active in T Tauri disks and show that ionization by cosmic-ray feedback extends the MRI-active region towards the disk mid-plane. However, the MRI is only active at the mid-plane of a minimum mass solar nebula disk if cosmic-rays propagate diffusively ($ζ\propto r^{-1}$) away from the star. The relationship between accretion, which accelerates cosmic rays, the dense accretion columns, which attenuate cosmic rays, and the MRI, which facilitates accretion, create a cosmic-ray feedback loop that mediates accretion and may produce luminosity variability. △ Less

Submitted 10 September, 2019; v1 submitted 21 August, 2019; originally announced August 2019.

Comments: 19 pages, 9 figures, accepted to ApJ

The Astrochemical Impact of Cosmic Rays in Protoclusters II: CI-to-H$_2$ and CO-to-H$_2$ Conversion Factors

Authors: Brandt A. L. Gaches, Stella S. R. Offner, Thomas G. Bisbas

Abstract: We utilize a modified astrochemistry code which includes cosmic ray attenuation in-situ to quantify the impact of different cosmic-ray models on the CO-to-H$_2$ and CI-to-H$_2$ conversion factors, $X_{\rm CO}$ and $X_{\rm CI}$, respectively. We consider the impact of cosmic rays accelerated by accretion shocks, and show that clouds with star formation efficiencies greater than 2\% have… ▽ More We utilize a modified astrochemistry code which includes cosmic ray attenuation in-situ to quantify the impact of different cosmic-ray models on the CO-to-H$_2$ and CI-to-H$_2$ conversion factors, $X_{\rm CO}$ and $X_{\rm CI}$, respectively. We consider the impact of cosmic rays accelerated by accretion shocks, and show that clouds with star formation efficiencies greater than 2\% have $X_{\rm CO} = (2.5 \pm 1)\times10^{20}$ cm$^{-2}$(K km s$^{-1}$)$^{-1}$, consistent with Milky Way observations. We find that changing the cosmic ray ionization rate from external sources from the canonical $ζ\approx 10^{-17}$ to $ζ\approx 10^{-16}$ s$^{-1}$, which better represents observations in diffuse gas, reduces $X_{\rm CO}$ by 0.2 dex for clusters with surface densities below 3 g cm$^{-2}$. We show that embedded sources regulate $X_{\rm CO}$ and decrease its variance across a wide range of surface densities and star formation efficiencies. Our models reproduce the trends of a decreased $X_{\rm CO}$ in extreme cosmic ray environments. $X_{\rm CI}$ has been proposed as an alternative to $X_{\rm CO}$ due to its brightness at high redshifts. The inclusion of internal cosmic ray sources leads to 1.2 dex dispersion in $X_{\rm CI}$ ranging from $2\times10^{20} < X_{\rm CI} < 4\times10^{21}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$. We show that $X_{\rm CI}$ is highly sensitive to the underlying cosmic ray model. △ Less

Submitted 19 August, 2019; originally announced August 2019.

Comments: 9 pages, 3 figures. ApJ Accepted