Engineering topics

Offner, Stella S. R.

Publications and source records attributed to Offner, Stella S. R..

18 records

DOE OSTI2017

Impact of Protostellar Outflows on Turbulence and Star Formation Efficiency in Magnetized Dense Cores

The star-forming efficiency of dense gas is thought to be set within cores by outflow and radiative feedback. We use magnetohydrodynamic simulations to investigate the relation between protostellar outflow evolution, turbulence, and star formation efficiency. We model the collapse and evolution of isolated dense cores for ≳0.5 Myr including the effects of turbulence, radiation transfer, and both radiation and outflow feedback from forming protostars. We show that outflows drive and maintain turbulence in the core environment even with strong initial fields. The star formation efficiency decreases with increasing field strength, and the final efficiencies are 15%–40%. The Stage 0 lifetime, during which the protostellar mass is lower than that of the dense envelope, increases proportionally with the initial magnetic field strength and ranges from ∼0.1 to 0.4 Myr. The average accretion rate is well represented by a tapered turbulent core model, which is a function of the final protostellar mass and is independent of the magnetic field strength. By tagging material launched in the outflow, we demonstrate that the outflow entrains about three times the actual launched gas mass, a ratio that remains roughly constant in time regardless of the initial magnetic field strength. However, turbulent driving increases for stronger fields since momentum is more efficiently imparted to non-outflow material. The protostellar outflow momentum is highest during the first 0.1 Myr and declines thereafter by a factor of ≳10 as the accretion rate diminishes.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

The Green Bank Ammonia Survey: Dense Cores under Pressure in Orion A

We use data on gas temperature and velocity dispersion from the Green Bank Ammonia Survey and core masses and sizes from the James Clerk Maxwell Telescope Gould Belt Survey to estimate the virial states of dense cores within the Orion A molecular cloud. Surprisingly, we find that almost none of the dense cores are sufficiently massive to be bound when considering only the balance between self-gravity and the thermal and non-thermal motions present in the dense gas. Including the additional pressure binding imposed by the weight of the ambient molecular cloud material and additional smaller pressure terms, however, suggests that most of the dense cores are pressure-confined.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Alignment between Protostellar Outflows and Filamentary Structure

We present new Submillimeter Array (SMA) observations of CO(2–1) outflows toward young, embedded protostars in the Perseus molecular cloud as part of the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. For 57 Perseus protostars, we characterize the orientation of the outflow angles and compare them with the orientation of the local filaments as derived from Herschel observations. We find that the relative angles between outflows and filaments are inconsistent with purely parallel or purely perpendicular distributions. Instead, the observed distribution of outflow-filament angles are more consistent with either randomly aligned angles or a mix of projected parallel and perpendicular angles. A mix of parallel and perpendicular angles requires perpendicular alignment to be more common by a factor of ∼3. Our results show that the observed distributions probably hold regardless of the protostar’s multiplicity, age, or the host core’s opacity. These observations indicate that the angular momentum axis of a protostar may be independent of the large-scale structure. We discuss the significance of independent protostellar rotation axes in the general picture of filament-based star formation.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

The Green Bank Ammonia Survey: First Results of NH{sub 3} Mapping of the Gould Belt

We present an overview of the first data release (DR1) and first-look science from the Green Bank Ammonia Survey (GAS). GAS is a Large Program at the Green Bank Telescope to map all Gould Belt star-forming regions with A{sub V}≳7 mag visible from the northern hemisphere in emission from NH{sub 3} and other key molecular tracers. This first release includes the data for four regions in the Gould Belt clouds: B18 in Taurus, NGC 1333 in Perseus, L1688 in Ophiuchus, and Orion A North in Orion. We compare the NH{sub 3} emission to dust continuum emission from Herschel and find that the two tracers correspond closely. We find that NH{sub 3} is present in over 60% of the lines of sight with A{sub V}≳7 mag in three of the four DR1 regions, in agreement with expectations from previous observations. The sole exception is B18, where NH{sub 3} is detected toward ∼40% of the lines of sight with A{sub V}≳7 mag. Moreover, we find that the NH{sub 3} emission is generally extended beyond the typical 0.1 pc length scales of dense cores. We produce maps of the gas kinematics, temperature, and NH{sub 3} column densities through forward modeling of the hyperfine structure of the NH{sub 3} (1, 1) and (2, 2) lines. We show that the NH{sub 3} velocity dispersion, σ{sub v}, and gas kinetic temperature, T {sub K}, vary systematically between the regions included in this release, with an increase in both the mean value and the spread of σ{sub v} and T {sub K} with increasing star formation activity. The data presented in this paper are publicly available (https://dataverse.harvard.edu/dataverse/GAS-DR1).

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

AN EXPLORATION OF THE STATISTICAL SIGNATURES OF STELLAR FEEDBACK

All molecular clouds are observed to be turbulent, but the origin, means of sustenance, and evolution of the turbulence remain debated. One possibility is that stellar feedback injects enough energy into the cloud to drive observed motions on parsec scales. Recent numerical studies of molecular clouds have found that feedback from stars, such as protostellar outflows and winds, injects energy and impacts turbulence. We expand upon these studies by analyzing magnetohydrodynamic simulations of molecular clouds, including stellar winds, with a range of stellar mass-loss rates and magnetic field strengths. We generate synthetic {sup 12}CO(1–0) maps assuming that the simulations are at the distance of the nearby Perseus molecular cloud. By comparing the outputs from different initial conditions and evolutionary times, we identify differences in the synthetic observations and characterize these using common astrostatistics. We quantify the different statistical responses using a variety of metrics proposed in the literature. We find that multiple astrostatistics, including the principal component analysis, the spectral correlation function, and the velocity coordinate spectrum (VCS), are sensitive to changes in stellar mass-loss rates and/or time evolution. A few statistics, including the Cramer statistic and VCS, are sensitive to the magnetic field strength. These findings demonstrate that stellar feedback influences molecular cloud turbulence and can be identified and quantified observationally using such statistics.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

ALMA Cycle 1 observations of the HH 46/47 molecular outflow: structure, entrainment, and core impact

We present Atacama Large Millimeter/sub-millimeter Array Cycle 1 observations of the HH 46/47 molecular outflow using combined 12 m array and Atacama Compact Array observations. The improved angular resolution and sensitivity of our multi-line maps reveal structures that help us study the entrainment process in much more detail and allow us to obtain more precise estimates of outflow properties than in previous observations. We use {sup 13}CO (1–0) and C{sup 18}O (1–0) emission to correct for the {sup 12}CO (1–0) optical depth to accurately estimate the outflow mass, momentum, and kinetic energy. This correction increases the estimates of the mass, momentum, and kinetic energy by factors of about 9, 5, and 2, respectively, with respect to estimates assuming optically thin emission. The new {sup 13}CO and C{sup 18}O data also allow us to trace denser and slower outflow material than that traced by the {sup 12}CO maps, and they reveal an outflow cavity wall at very low velocities (as low as 0.2 km s{sup −1} with respect to the core’s central velocity). Adding the slower material traced only by {sup 13}CO and C{sup 18}O, there is another factor of three increase in the mass estimate and 50% increase in the momentum estimate. The estimated outflow properties indicate that the outflow is capable of dispersing the parent core within the typical lifetime of the embedded phase of a low-mass protostar and that it is responsible for a core-to-star efficiency of 1/4 to 1/3. We find that the outflow cavity wall is composed of multiple shells associated with a series of jet bow-shock events. Within about 3000 au of the protostar the {sup 13}CO and C{sup 18}O emission trace a circumstellar envelope with both rotation and infall motions, which we compare with a simple analytic model. The CS (2–1) emission reveals tentative evidence of a slowly moving rotating outflow, which we suggest is entrained not only poloidally but also toroidally by a disk wind that is launched from relatively large radii from the source.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

THE TURBULENT ORIGIN OF OUTFLOW AND SPIN MISALIGNMENT IN MULTIPLE STAR SYSTEMS

The protostellar outflows of wide-separation forming binaries frequently appear misaligned. We use magneto-hydrodynamic simulations to investigate the alignment of protostellar spin and molecular outflows for forming binary pairs. We show that the protostellar pairs, which form from turbulent fragmentation within a single parent core, have randomly oriented angular momentum. Although the pairs migrate to closer separations, their spins remain partially misaligned. We produce {sup 12}CO(2-1) synthetic observations of the simulations and characterize the outflow orientation in the emission maps. The CO-identified outflows exhibit a similar random distribution and are also statistically consistent with the observed distribution of molecular outflows. We conclude that the observed misalignment provides a clear signature of binary formation via turbulent fragmentation. The persistence of misaligned outflows and stellar spins following dynamical evolution may provide a signature of binary origins for more evolved multiple star systems.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

Mass assembly of stellar systems and their evolution with the SMA (MASSES). multiplicity and the physical environment in L1448N

We present continuum and molecular line observations at 230 and 345 GHz from the Submillimeter Array (SMA) toward three protostars in the Perseus L1448N region. The data are from the large project “Mass Assembly of Stellar Systems and their Evolution with the SMA.” Three dust continuum sources, Source B, Source NW, and Source A, are detected at both frequencies. These sources have corresponding emission peaks in C{sup 18}O (J=2→1), {sup 13}CO (J=2→1), and HCO{sup +} (J=4→3), and have offsets with N{sub 2}D{sup +} (J=3→2) peaks. High angular resolution data from a complementary continuum survey with the Karl G. Jansky Very Large Array show that Source B is associated with three 8 mm continuum objects, Source NW with two, and Source A remains single. These results suggest that multiplicity in L1448N exists at different spatial scales from a few thousand AU to <100 AU. Velocity gradients in each source obtained from two-dimensional fits to the SMA C{sup 18}O emission are found to be perpendicular to within 20° of the outflow directions as revealed by {sup 12}CO (J=2→1). We have observed that Sources B and NW with multiplicity have higher densities than Source A without multiplicity. This suggests that thermal Jeans fragmentation can be relevant in the fragmentation process. However, we have not observed a difference in the ratio between rotational and gravitational energy between sources with and without multiplicity. We also have not observed a trend between non-thermal velocity dispersions and the level of fragmentation. Our study has provided the first direct and comprehensive comparison between multiplicity and core properties in low-mass protostars, although based on small number statistics.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

Impact of winds from intermediate-mass stars on molecular cloud structure and turbulence

Observations of nearby molecular clouds detect “shells,” which are likely caused by winds from young main sequence stars. However, the progenitors of these observed features are not well characterized and the mass-loss rates inferred from the gas kinematics are several orders of magnitude greater than those predicted by atomic line-driven stellar wind models. We use magnetohydrodynamic simulations to model winds launching within turbulent molecular clouds and explore the impact of wind properties on cloud morphology and turbulence. We find that winds do not produce clear features in turbulent statistics such as the Fourier spectra of density and momentum but do impact the Fourier velocity spectrum. The density and velocity distribution functions, especially as probed by CO spectral lines, strongly indicate the presence and influence of winds. We show that stellar mass-loss rates for individual stars must be m-dot {sub w}≳10{sup −7} M{sub ⊙} yr{sup −1}, similar to those estimated from observations, to reproduce shell properties. Consequently, we conclude that B and A-type main sequence stars have mass-loss rates several orders of magnitude larger that those predicted by models or that young stars are more variable than expected due to magnetic activity or accretion.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

Predictions for observing protostellar outflows with ALMA

Protostellar outflows provide a means to probe the accretion process of forming stars and their ability to inject energy into their surroundings. However, conclusions based on outflow observations depend upon the degree of accuracy with which their properties can be estimated. We examine the quality of Atacama Large Millimeter/submillimeter Array (ALMA) observations of protostellar outflows by producing synthetic {sup 12}CO(1-0) and {sup 13}CO(1-0) observations of numerical simulations. We use various ALMA configurations, observational parameters, and outflow inclinations to assess how accurately different assumptions and setups can recover underlying properties. We find that more compact arrays and longer observing times can improve the mass and momentum recovery by a factor of two. During the first ∼0.3 Myr of evolution, {sup 12}CO(1-0) is optically thick, even for velocities |v|⩾1 km s{sup −1}, and outflow mass is severely underestimated without an optical depth correction. Likewise, {sup 13}CO(1-0) is optically thick during the first ≃0.1 Myr. However, underestimation due to shorter observing time, missing flux, and optical depth are partially offset by the assumption of local thermodynamic equilibrium and higher excitation temperatures. Overall, we expect that full ALMA {sup 13}CO(1-0) observations of protostellar sources within 500 pc with observing times ≳1 hr and assumed excitation temperatures of T<20 K will reliably measure mass and line of sight momentum to within 20%.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

ASTROCHEMICAL CORRELATIONS IN MOLECULAR CLOUDS

We investigate the spectral correlations between different species used to observe molecular clouds. We use hydrodynamic simulations and a full chemical network to study the abundances of over 150 species in typical Milky Way molecular clouds. We perform synthetic observations in order to produce emission maps of a subset of these tracers. We study the effects of different lines of sight and spatial resolution on the emission distribution and perform a robust quantitative comparison of the species to each other. We use the Spectral Correlation Function (SCF), which quantifies the root mean squared difference between spectra separated by some length scale, to characterize the structure of the simulated cloud in position-position-velocity (PPV) space. We predict the observed SCF for a broad range of observational tracers, and thus identify homologous species. In particular, we show that the pairs C and CO, C{sup +} and CN, and NH{sub 3} and H{sub 2}CS have very similar SCFs. We measure the SCF slope variation as a function of beam size for all species and demonstrate that the beam size has a distinct effect on different species emission. However, for beams of up to 10'', placing the cloud at 1 kpc, the change is not large enough to move the SCF slopes into different regions of parameter space. The results from this study provide observational guidance for choosing the best tracer to probe various cloud length scales.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Investigations of protostellar outflow launching and gas entrainment: Hydrodynamic simulations and molecular emission

We investigate protostellar outflow evolution, gas entrainment, and star formation efficiency using radiation-hydrodynamic simulations of isolated, turbulent low-mass cores. We adopt an X-wind launching model, in which the outflow rate is coupled to the instantaneous protostellar accretion rate and evolution. We vary the outflow collimation angle from θ = 0.01-0.1 and find that even well-collimated outflows effectively sweep up and entrain significant core mass. The Stage 0 lifetime ranges from 0.14-0.19 Myr, which is similar to the observed Class 0 lifetime. The star formation efficiency of the cores spans 0.41-0.51. In all cases, the outflows drive strong turbulence in the surrounding material. Although the initial core turbulence is purely solenoidal by construction, the simulations converge to approximate equipartition between solenoidal and compressive motions due to a combination of outflow driving and collapse. When compared to simulation of a cluster of protostars, which is not gravitationally centrally condensed, we find that the outflows drive motions that are mainly solenoidal. The final turbulent velocity dispersion is about twice the initial value of the cores, indicating that an individual outflow is easily able to replenish turbulent motions on sub-parsec scales. We post-process the simulations to produce synthetic molecular line emission maps of {sup 12}CO, {sup 13}CO, and C{sup 18}O and evaluate how well these tracers reproduce the underlying mass and velocity structure.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Synthetic observations of the evolution of starless cores in a molecular cloud simulation: Comparisons with JCMT data and predictions for alma

Interpreting the nature of starless cores has been a prominent goal in star formation for many years. In order to characterize the evolutionary stages of these objects, we perform synthetic observations of a numerical simulation of a turbulent molecular cloud. We find that nearly all cores that we detect are associated with filaments and eventually form protostars. We conclude that observed starless cores that appear Jeans unstable are only marginally larger than their respective Jeans masses (within a factor of three). We note single dish observations such as those performed with the James Clerk Maxwell Telescope appear to miss significant core structure on small scales because of beam averaging. Finally, we predict that interferometric observations with Atacama Large Millimeter Array Cycle 1 will resolve the important small-scale structure, which has so far been missed by millimeter wavelength observations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

MODELING THE ATOMIC-TO-MOLECULAR TRANSITION AND CHEMICAL DISTRIBUTIONS OF TURBULENT STAR-FORMING CLOUDS

We use 3D-PDR, a three-dimensional astrochemistry code for modeling photodissociation regions (PDRs), to post-process hydrodynamic simulations of turbulent, star-forming clouds. We focus on the transition from atomic to molecular gas, with specific attention to the formation and distribution of H, C{sup +}, C, H{sub 2}, and CO. First, we demonstrate that the details of the cloud chemistry and our conclusions are insensitive to the simulation spatial resolution, to the resolution at the cloud edge, and to the ray angular resolution. We then investigate the effect of geometry and simulation parameters on chemical abundances and find weak dependence on cloud morphology as dictated by gravity and turbulent Mach number. For a uniform external radiation field, we find similar distributions to those derived using a one-dimensional PDR code. However, we demonstrate that a three-dimensional treatment is necessary for a spatially varying external field, and we caution against using one-dimensional treatments for non-symmetric problems. We compare our results with the work of Glover et al., who self-consistently followed the time evolution of molecule formation in hydrodynamic simulations using a reduced chemical network. In general, we find good agreement with this in situ approach for C and CO abundances. However, the temperature and H{sub 2} abundances are discrepant in the boundary regions (A{sub v} {<=} 5), which is due to the different number of rays used by the two approaches.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

OBSERVING SIMULATED PROTOSTARS WITH OUTFLOWS: HOW ACCURATE ARE PROTOSTELLAR PROPERTIES INFERRED FROM SEDs?

The properties of unresolved protostars and their local environment are frequently inferred from spectral energy distributions (SEDs) using radiative transfer modeling. In this paper, we use synthetic observations of realistic star formation simulations to evaluate the accuracy of properties inferred from fitting model SEDs to observations. We use ORION, an adaptive mesh refinement (AMR) three-dimensional gravito-radiation-hydrodynamics code, to simulate low-mass star formation in a turbulent molecular cloud including the effects of protostellar outflows. To obtain the dust temperature distribution and SEDs of the forming protostars, we post-process the simulations using HYPERION, a state-of-the-art Monte Carlo radiative transfer code. We find that the ORION and HYPERION dust temperatures typically agree within a factor of two. We compare synthetic SEDs of embedded protostars for a range of evolutionary times, simulation resolutions, aperture sizes, and viewing angles. We demonstrate that complex, asymmetric gas morphology leads to a variety of classifications for individual objects as a function of viewing angle. We derive best-fit source parameters for each SED through comparison with a pre-computed grid of radiative transfer models. While the SED models correctly identify the evolutionary stage of the synthetic sources as embedded protostars, we show that the disk and stellar parameters can be very discrepant from the simulated values, which is expected since the disk and central source are obscured by the protostellar envelope. Parameters such as the stellar accretion rate, stellar mass, and disk mass show better agreement, but can still deviate significantly, and the agreement may in some cases be artificially good due to the limited range of parameters in the set of model SEDs. Lack of correlation between the model and simulation properties in many individual instances cautions against overinterpreting properties inferred from SEDs for unresolved protostellar sources.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2011

RADIATION-HYDRODYNAMIC SIMULATIONS OF PROTOSTELLAR OUTFLOWS: SYNTHETIC OBSERVATIONS AND DATA COMPARISONS

We present results from three-dimensional, self-gravitating, radiation-hydrodynamic simulations of low-mass protostellar outflows. We construct synthetic observations in {sup 12}CO in order to compare with observed outflows and evaluate the effects of beam resolution and outflow orientation on inferred outflow properties. To facilitate the comparison, we develop a quantitative prescription for measuring outflow opening angles. Using this prescription, we demonstrate that, in both simulations and synthetic observations, outflow opening angles broaden with time similarly to observed outflows. However, the interaction between the outflowing gas and the turbulent core envelope produces significant asymmetry between the redshifted and blueshifted outflow lobes. We find that applying a velocity cutoff may result in outflow masses that are underestimated by a factor five or more, and masses derived from optically thick CO emission further underpredict the mass of the high-velocity gas by a factor of 5-10. Derived excitation temperatures indicate that outflowing gas is hotter than the ambient gas with temperature rising over time, which is in agreement with the simulation gas temperatures. However, excitation temperatures are otherwise not well correlated with the actual gas temperature.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2011

ON THE RELIABILITY OF STELLAR AGES AND AGE SPREADS INFERRED FROM PRE-MAIN-SEQUENCE EVOLUTIONARY MODELS

We revisit the problem of low-mass pre-main-sequence stellar evolution and its observational consequences for where stars fall on the Hertzsprung-Russell diagram (HRD). In contrast to most previous work, our models follow stars as they grow from small masses via accretion, and we perform a systematic study of how the stars' HRD evolution is influenced by their initial radius, by the radiative properties of the accretion flow, and by the accretion history, using both simple idealized accretion histories and histories taken from numerical simulations of star cluster formation. We compare our numerical results to both non-accreting isochrones and to the positions of observed stars in the HRD, with a goal of determining whether both the absolute ages and the age dispersions inferred from non-accreting isochrones are reliable. We show that non-accreting isochrones can sometimes overestimate stellar ages for more massive stars (those with effective temperatures above {approx}3500 K), thereby explaining why non-accreting isochrones often suggest a systematic age difference between more and less massive stars in the same cluster. However, we also find the only way to produce a similar overestimate for the ages of cooler stars is if these stars grow from {approx}0.01 M{sub sun} seed protostars that are an order of magnitude smaller than predicted by current theoretical models, and if the size of the seed protostar correlates systematically with the final stellar mass at the end of accretion. We therefore conclude that, unless both of these conditions are met, inferred ages and age spreads for cool stars are reliable, at least to the extent that the observed bolometric luminosities and temperatures are accurate. Finally, we note that the time dependence of the mass accretion rate has remarkably little effect on low-mass stars' evolution on the HRD, and that such time dependence may be neglected for all stars except those with effective temperatures above {approx}4000 K.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2011

THE PROTOSTELLAR LUMINOSITY FUNCTION

The protostellar luminosity function (PLF) is the present-day luminosity function of the protostars in a region of star formation. It is determined using the protostellar mass function in combination with a stellar evolutionary model that provides the luminosity as a function of instantaneous and final stellar mass. In 2010, McKee and Offner considered three main accretion models: the isothermal sphere (IS) model, the turbulent core (TC) model, and an approximation of the competitive accretion (CA) model. We also consider the effect of an accretion rate that tapers off linearly in time and an accelerating star formation rate. For each model, we characterize the luminosity distribution using the mean, median, maximum, ratio of the median to the mean, standard deviation of the logarithm of the luminosity, and the fraction of very low luminosity objects. We compare the models with bolometric luminosities observed in local star-forming regions and find that models with an approximately constant accretion time, such as the TC and CA models, appear to agree better with observation than those with a constant accretion rate, such as the IS model. We show that observations of the mean protostellar luminosity in these nearby regions of low-mass star formation suggest a mean star formation time of 0.3 {+-} 0.1 Myr. Such a timescale, together with some accretion that occurs non-radiatively and some that occurs in high-accretion, episodic bursts, resolves the classical 'luminosity problem' in low-mass star formation, in which observed protostellar luminosities are significantly less than predicted. An accelerating star formation rate is one possible way of reconciling the observed star formation time and mean luminosity. Future observations will place tighter constraints on the observed luminosities, star formation time, and episodic accretion, enabling better discrimination between star formation models and clarifying the influence of variable accretion on the PLF.

79 ASTRONOMY AND ASTROPHYSICS↗