Researcher profile

Stella S. R. Offner

Stella S. R. Offner contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

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

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

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

Building this graph slice

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

Published work

26 published item(s)

preprint2023arXiv

Infall of material onto the filaments in Barnard 5

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.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2021arXiv

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

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.

preprint2021arXiv

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

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 &#34;stirred&#34; 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.

preprint2021arXiv

STARFORGE: The effects of protostellar outflows on the IMF

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.

preprint2021arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

Efficient methanol production on the dark side of a prestellar core

We present ALMA maps of the starless molecular cloud core Ophiuchus/H-MM1 in the lines of deuterated ammonia (ortho-NH2D), methanol (CH3OH), and sulphur monoxide (SO). The dense core is seen in NH2D emission, whereas the CH3OH and SO distributions form a halo surrounding the core. Because methanol is formed on grain surfaces, its emission highlights regions where desorption from grains is particularly efficient. Methanol and sulphur monoxide are most abundant in a narrow zone that follows the eastern side of the core. This side is sheltered from the stronger external radiation field coming from the west. We show that photodissociation on the illuminated side can give rise to an asymmetric methanol distribution, but that the stark contrast observed in H-MM1 is hard to explain without assuming enhanced desorption on the shaded side. The region of the brightest emission has a wavy structure that rolls up at one end. This is the signature of Kelvin-Helmholtz instability occurring in sheared flows. We suggest that in this zone, methanol and sulphur are released as a result of grain-grain collisions induced by shear vorticity.

preprint2020arXiv

Evolution of giant molecular clouds across cosmic time

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.

preprint2020arXiv

Impact of low-energy cosmic rays on star formation

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.

preprint2020arXiv

Physical Processes in Star Formation

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

The Origin of the Stellar Mass Distribution and Multiplicity

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.

preprint2020arXiv

The Physics of Star Cluster Formation and Evolution

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.

preprint2020arXiv

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

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&#34; (40 au), including observations with the VLA at 9 mm toward 148 protostars at a resolution of $\sim$0.08&#34; (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.

preprint2020arXiv

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

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&#39; 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&#39;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.

preprint2020arXiv

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

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.

preprint2019arXiv

Droplets II: Internal Velocity Structures and Potential Rotational Motions in Pressure-dominated Coherent Structures

We present an analysis of the internal velocity structures of the newly identified sub-0.1 pc coherent structures, droplets, in L1688 and B18. By fitting 2D linear velocity fields to the observed maps of velocity centroids, we determine the magnitudes of linear velocity gradients and examine the potential rotational motions that could lead to the observed velocity gradients. The results show that the droplets follow the same power-law relation between the velocity gradient and size found for larger-scale dense cores. Assuming that rotational motion giving rise to the observed velocity gradient in each core is a solid-body rotation of a rotating body with a uniform density, we derive the &#34;net rotational motions&#34; of the droplets. We find a ratio between rotational and gravitational energies, $β$, of $\sim 0.046$ for the droplets, and when including both droplets and larger-scale dense cores, we find $β\sim 0.039$. We then examine the alignment between the velocity gradient and the major axis of each droplet, using methods adapted from the histogram of relative orientations (HRO) introduced by Soler et al. (2013). We find no definitive correlation between the directions of velocity gradients and the elongations of the cores. Lastly, we discuss physical processes other than rotation that may give rise to the observed velocity field.

preprint2019arXiv

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

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&#39;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/.

preprint2019arXiv

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

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.

preprint2019arXiv

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

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.