Source author record

Amelia M. Stutz

Amelia M. Stutz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

30works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

30 published item(s)

preprint2022arXiv

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions -- XI. From inflow to infall in hub-filament systems

We investigate the presence of hub-filament systems in a large sample of 146 active proto-clusters, using H$^{13}$CO$^{+}$ J=1-0 molecular line data obtained from the ATOMS survey. We find that filaments are ubiquitous in proto-clusters, and hub-filament systems are very common from dense core scales ($\sim$0.1 pc) to clump/cloud scales ($\sim$1-10 pc). The proportion of proto-clusters containing hub-filament systems decreases with increasing dust temperature ($T_d$) and luminosity-to-mass ratios ($L/M$) of clumps, indicating that stellar feedback from H{\sc ii} regions gradually destroys the hub-filament systems as proto-clusters evolve. Clear velocity gradients are seen along the longest filaments with a mean velocity gradient of 8.71 km s$^{-1}$pc$^{-1}$ and a median velocity gradient of 5.54 km s$^{-1}$pc$^{-1}$. We find that velocity gradients are small for filament lengths larger than $\sim$1~pc, probably hinting at the existence of inertial inflows, although we cannot determine whether the latter are driven by large-scale turbulence or large-scale gravitational contraction. In contrast, velocity gradients below $\sim$1~pc dramatically increase as filament lengths decrease, indicating that the gravity of the hubs or cores starts to dominate gas infall at small scales. We suggest that self-similar hub-filament systems and filamentary accretion at all scales may play a key role in high-mass star formation.

preprint2022arXiv

Pre-main Sequence Brackett Emitters in the APOGEE DR17 Catalog: Line Strengths and Physical Properties of Accretion Columns

Very young (t $\lesssim$ 10 Myrs) stars possess strong magnetic fields that channel ionized gas from the interiors of their circumstellar discs to the surface of the star. Upon impacting the stellar surface, the shocked gas recombines and emits hydrogen spectral lines. To characterize the density and temperature of the gas within these accretion streams, we measure equivalent widths of Brackett (Br) 11-20 emission lines detected in 1101 APOGEE spectra of 326 likely pre-main sequence accretors. For sources with multiple observations, we measure median epoch-to-epoch line strength variations of 10% in Br11 and 20% in Br20. We also fit the measured line ratios to predictions of radiative transfer models by Kwan & Fischer. We find characteristic best-fit electron densities of $n_e$ = 10$^{11} - 10^{12}$ cm$^{-3}$, and excitation temperatures that are inversely correlated with electron density (from T$\sim$5000 K for $n_e \sim 10^{12}$ cm$^{-3}$, to T$\sim$12500 K at $n_e \sim 10^{11}$ cm$^{-3}$). These physical parameters are in good agreement with predictions from modelling of accretion streams that account for the hydrodynamics and radiative transfer within the accretion stream. We also present a supplementary catalog of line measurements from 9733 spectra of 4255 Brackett emission line sources in the APOGEE DR17 dataset.

preprint2022arXiv

Stellar Properties for a Comprehensive Collection of Star Forming Regions in the SDSS APOGEE-2 Survey

The Sloan Digital Sky Survey IV (SDSS-IV) APOGEE-2 primary science goal was to observe red giant stars throughout the Galaxy to study its dynamics, morphology, and chemical evolution. The APOGEE instrument, a high-resolution 300 fiber H-band (1.55-1.71 micron) spectrograph, is also ideal to study other stellar populations in the Galaxy, among which are a number of star forming regions and young open clusters. We present the results of the determination of six stellar properties ($T_{eff}$, $\log{g}$, [Fe/H], $L/L_\odot$, $M/M_\odot$, and ages) for a sample that is composed of 3360 young stars, of sub-solar to super-solar types, in sixteen Galactic star formation and young open cluster regions. Those sources were selected by using a clustering method that removes most of the field contamination. Samples were also refined by removing targets affected by various systematic effects of the parameter determination. The final samples are presented in a comprehensive catalog that includes all six estimated parameters. This overview study also includes parameter spatial distribution maps for all regions and Hertzprung-Russell ($L/L_\odot$ vs. $T_{eff}$) diagrams. This study serves as a guide for detailed studies on individual regions, and paves the way for the future studies on the global properties of stars in the pre-main sequence phase of stellar evolution using more robust samples.

preprint2022arXiv

The Rate, Amplitude and Duration of Outbursts from Class 0 Protostars in Orion

At least half of a protostar's mass is accreted in the Class 0 phase, when the central protostar is deeply embedded in a dense, infalling envelope. We present the first systematic search for outbursts from Class 0 protostars in the Orion clouds. Using photometry from Spitzer/IRAC spanning 2004 to 2017, we detect three outbursts from Class 0 protostars with $\ge 2$ mag changes at 3.6 or 4.5 $μ$m. This is comparable to the magnitude change of a known protostellar FU Ori outburst. Two are newly detected bursts from the protostars HOPS 12 and 124. The number of detections implies that Class 0 protostars burst every 438 yr, with a 95% confidence interval of 161 to 1884 yr. Combining Spitzer and WISE/NEOWISE data spanning 2004-2019, we show that the bursts persist for more than nine years with significant variability during each burst. Finally, we use $19-100$ $μ$m photometry from SOFIA, Spitzer and Herschel to measure the amplitudes of the bursts. Based on the burst interval, a duration of 15 yr, and the range of observed amplitudes, 3-100% of the mass accretion during the Class 0 phase occurs during bursts. In total, we show that bursts from Class 0 protostars are as frequent, or even more frequent, than those from more evolved protostars. This is consistent with bursts being driven by instabilities in disks triggered by rapid mass infall. Furthermore, we find that bursts may be a significant, if not dominant, mode of mass accretion during the Class 0 phase.

preprint2021arXiv

Final Targeting Strategy for the SDSS-IV APOGEE-2N Survey

APOGEE-2 is a dual-hemisphere, near-infrared (NIR), spectroscopic survey with the goal of producing a chemo-dynamical mapping of the Milky Way Galaxy. The targeting for APOGEE-2 is complex and has evolved with time. In this paper, we present the updates and additions to the initial targeting strategy for APOGEE-2N presented in Zasowski et al. (2017). These modifications come in two implementation modes: (i) "Ancillary Science Programs" competitively awarded to SDSS-IV PIs through proposal calls in 2015 and 2017 for the pursuit of new scientific avenues outside the main survey, and (ii) an effective 1.5-year expansion of the survey, known as the Bright Time Extension, made possible through accrued efficiency gains over the first years of the APOGEE-2N project. For the 23 distinct ancillary programs, we provide descriptions of the scientific aims, target selection, and how to identify these targets within the APOGEE-2 sample. The Bright Time Extension permitted changes to the main survey strategy, the inclusion of new programs in response to scientific discoveries or to exploit major new datasets not available at the outset of the survey design, and expansions of existing programs to enhance their scientific success and reach. After describing the motivations, implementation, and assessment of these programs, we also leave a summary of lessons learned from nearly a decade of APOGEE-1 and APOGEE-2 survey operations. A companion paper, Santana et al. (submitted), provides a complementary presentation of targeting modifications relevant to APOGEE-2 operations in the Southern Hemisphere.

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

Extended HNCO, SiO, and HC$_{3}$N emission in 43 southern star-forming regions

We have selected 43 southern massive star-forming regions to study the spatial distribution of HNCO 4$_{04}$-3$_{03}$, SiO 2-1 and HC$_{3}$N 10-9 line emission and to investigate their spatial association with the dust emission. The morphology of HNCO 4$_{04}$-3$_{03}$ and HC$_{3}$N 10-9 agrees well with the dust emission. HC$_{3}$N 10-9 tends to originate from more compact regions than HNCO 4$_{04}$-3$_{03}$ and SiO 2-1. We divided our sources into three groups: those in the Central Molecular Zone (CMZ), those associated with bubbles (Bubble), and the remaining sources, which are termed 'normal star forming regions' (NMSFR). These three groups, subdivided into three different categories with respect to line widths, integrated intensities, and column densities, hint at the presence of different physical and chemical processes. We find that the dust temperature $T_{\rm d}$, and the abundance ratios of $N_{\rm HNCO}/N_{\rm SiO}$ and $N_{\rm HNCO}/N_{\rm HC3N}$ show a decreasing trend towards the central dense regions of CMZ sources, while $N_{\rm HC3N}/N_{\rm SiO}$ moves into the opposite direction. Moreover, a better agreement is found between $T_{\rm d}$ and $N_{\rm HC3N}/N_{\rm SiO}$ in Bubble and NMSFR category sources. Both outflow and inflow activities have been found in eight of the sixteen bubble and NMSFR sources. The low outflow detection rate indicates that in these sources the SiO 2-1 line wing emission is either below our sensitivity limit or that the bulk of the SiO emission may be produced by the expansion of an H{\sc\,ii} region or supernova remnant, which has pushed molecular gas away forming a shock and yielding SiO.

preprint2020arXiv

Filament rotation in the California L1482 cloud

We analyze the gas mass distribution, the gas kinematics, and the young stellar object (YSO) content of the California Molecular Cloud (CMC) L1482 filament. We derive a Gaia DR2 YSO distance of 511$^{+17}_{-16}$ pc. We derive scale-free power-laws for the mean gas line-mass (M/L) profiles; we calculate the gravitational potential and field profiles consistent with these. We present IRAM 30 m C$^{18}$O (1-0) (and other tracers) position-velocity (PV) diagrams that exhibit complex velocity twisting and turning structures. We find a rotational profile in C$^{18}$O perpendicular to the southern filament ridgeline. The profile is regular, confined ($r\lesssim0.4$ pc), anti-symmetric, and to first order linear with a break at $r\sim0.25$ pc. The timescales of the inner (outer) gradients are $\sim$0.7 (6.0) Myr. We show that the centripetal force, compared to gravity, increases toward the break; when the ratio of forces approaches unity, the profile turns over, just before filament breakup is achieved. The timescales and relative roles of gravity to rotation indicate that the structure is stable, long lived ($\sim$ a few times 6 Myr), and undergoing outside-in evolution. Moreover, this filament has practically no star formation, a perpendicular Planck plane-of-the-sky (POS) magnetic field morphology, and POS "zig-zag" morphology, which together with the rotation profile lead to the suggestion that the 3D shape is a corkscrew filament with a helical magnetic field. These results, combined with results in Orion and G035.39-00.33, suggest evolution toward higher densities as rotating filaments shed angular momentum. Thus, magnetic fields may be an essential feature of high-mass (M $\sim10^5$ M$_{\odot}$) cloud filament evolution toward cluster formation.

preprint2020arXiv

The Herschel Orion Protostar Survey: Far-Infrared Photometry and Colors of Protostars and Their Variations across Orion A and B

The degree to which the properties of protostars are affected by environment remains an open question. To investigate this, we look at the Orion A and B molecular clouds, home to most of the protostars within 500 pc. At ~400 pc, Orion is close enough to distinguish individual protostars across a range of environments in terms of both the stellar and gas projected densities. As part of the Herschel Orion Protostar Survey (HOPS), we used the Photodetector Array Camera and Spectrometer (PACS) to map 108 partially overlapping square fields with edge lengths of 5 arcmin or 8 arcmin and measure the 70 micron and 160 micron flux densities of 338 protostars within them. In this paper we examine how these flux densities and their ratio depend on evolutionary state and environment within the Orion complex. We show that Class 0 protostars occupy a region of the 70 micron flux density versus 160 micron to 70 micron flux density ratio diagram that is distinct from their more evolved counterparts. We then present evidence that the Integral-Shaped Filament (ISF) and Orion B contain protostars with more massive envelopes than those in the more sparsely populated LDN 1641 region. This can be interpreted as evidence for increasing star formation rates in the ISF and Orion B or as a tendency for more massive envelopes to be inherited from denser birth environments. We also provide technical details about the map-making and photometric procedures used in the HOPS program.

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" (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.

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.

preprint2016arXiv

Slingshot Mechanism in Orion: Kinematic Evidence For Ejection of Protostars by Filaments

By comparing 3 constituents of Orion A (gas, protostars, and pre-main-sequence stars), both morhologically and kinematically, we derive the following. The gas surface density near the integral-shaped filament (ISF) is well represented by a power law, Sigma(b)=72 Msun/pc^2(b/pc)^{-5/8} for our entire range, 0.05<b/pc<8.5, of distance from the filament ridge. Essentially all protostars lie on the ISF or other filament ridges, while almost all pre-main-sequence stars do not. Combined with the fact that protostars move <1 kms relative to the filaments while stars move several times faster, this implies that protostellar accretion is terminated by a slingshot ejection from the filaments. The ISF is the 3rd in a series of star bursts that are progressively moving south, with separations of a few Myr in time and 3 pc in space. This, combined with the filament's observed undulations (spatial and velocity), suggests that repeated propagation of transverse waves thru the filament is progressively digesting the material that formerly connected Orion A and B into stars in discrete episodes. We construct an axially symmetric gas density profile rho(r)=16 Msun/pc^3(r/pc)^{-13/8}. The model implies that the observed magnetic fields are supercritical on scales of the observed undulations, suggesting that the filament's transverse waves are magnetically induced. Because the magnetic fields are subcritical on scales of the filament on larger scales, the system as a whole is relatively stable and long lived. Protostellar ejection occurs because the gas accelerates away from the protostars, not the other way around. The model also implies that the ISF is kinematically young, which is consistent with other lines of evidence. The southern filament has a broken power law, which matches the ISF profile for 2.5<b/pc<8.5, but is shallower closer in. It is also kinematically older than the ISF.

preprint2015arXiv

Evolution of Mass Outflow in Protostars

We have surveyed 84 Class 0, Class I, and flat-spectrum protostars in mid-infrared [Si II], [Fe II] and [S I] line emission, and 11 of these in far-infrared [O I] emission. We use the results to derive their mass outflow rates. Thereby we observe a strong correlation of mass outflow rates with bolometric luminosity, and with the inferred mass accretion rates of the central objects, which continues through the Class 0 range the trend observed in Class II young stellar objects. Along this trend from large to small mass-flow rates, the different classes of young stellar objects lie in the sequence Class 0 -- Class I/flat-spectrum -- Class II, indicating that the trend is an evolutionary sequence in which mass outflow and accretion rates decrease together with increasing age, while maintaining rough proportionality. The survey results include two which are key tests of magnetocentrifugal outflow-acceleration mechanisms: the distribution of the outflow/accretion branching ratio b, and limits on the distribution of outflow speeds. Neither rule out any of the three leading outflow-acceleration, angular-momentum-ejection mechanisms, but they provide some evidence that disk winds and accretion-powered stellar winds (APSWs) operate in many protostars. An upper edge observed in the branching-ratio distribution is consistent with the upper bound of b = 0.6 found in models of APSWs, and a large fraction (0.31) of the sample have branching ratio sufficiently small that only disk winds, launched on scales as large as several AU, have been demonstrated to account for them.

preprint2015arXiv

HOPS 383: An Outbursting Class 0 Protostar in Orion

We report the dramatic mid-infrared brightening between 2004 and 2006 of HOPS 383, a deeply embedded protostar adjacent to NGC 1977 in Orion. By 2008, the source became a factor of 35 brighter at 24 microns with a brightness increase also apparent at 4.5 microns. The outburst is also detected in the submillimeter by comparing APEX/SABOCA to SCUBA data, and a scattered-light nebula appeared in NEWFIRM K_s imaging. The post-outburst spectral energy distribution indicates a Class 0 source with a dense envelope and a luminosity between 6 and 14 L_sun. Post-outburst time-series mid- and far-infrared photometry shows no long-term fading and variability at the 18% level between 2009 and 2012. HOPS 383 is the first outbursting Class 0 object discovered, pointing to the importance of episodic accretion at early stages in the star formation process. Its dramatic rise and lack of fading over a six-year period hint that it may be similar to FU Ori outbursts, although the luminosity appears to be significantly smaller than the canonical luminosities of such objects.

preprint2015arXiv

Measurement of [HDCO]/[H2CO] Ratios in the Envelopes of Extremely Cold Protostars in Orion

We present observations of HDCO and H2CO emission toward a sample of 15 Class 0 protostars in the Orion A and B clouds. Of these, eleven protostars are Herschel-identified PACS Bright Red Sources (PBRS) and four are previously identified protostars. Our observations revealed the chemical properties of the PBRS envelope for the first time. The column densities of HDCO and H2CO are derived from single dish observations at an angular resolution of ~20 arcsec (~8400 AU). The degree of deuteration in H2CO ([HDCO]/[H2CO]) was estimated to range from 0.03 to 0.31. The deuterium fractionation of most PBRS (70%) is similar to that of the non-PBRS sources. Three PBRS (30%) exhibit high deuterium fractionation, larger than 0.15. The large variation of the deuterium fractionation of H2CO in the whole PBRS sample may reflect the diversity in the initial conditions of star forming cores. There is no clear correlation between the [HDCO]/[H2CO] ratio and the evolutionary sequence of protostars.

preprint2014arXiv

Characterizing the Youngest Herschel-detected Protostars I. Envelope Structure Revealed by CARMA Dust Continuum Observations

We present CARMA 2.9 mm dust continuum emission observations of a sample of 14 Herschel-detected Class 0 protostars in the Orion A and B molecular clouds, drawn from the PACS Bright Red Sources (PBRS) sample (Stutz et al.). These objects are characterized by very red 24 \micron\ to 70 \micron\ colors and prominent submillimeter emission, suggesting that they are very young Class 0 protostars embedded in dense envelopes. We detect all of the PBRS in 2.9 mm continuum emission and emission from 4 protostars and 1 starless core in the fields toward the PBRS; we also report 1 new PBRS source. The ratio of 2.9 mm luminosity to bolometric luminosity is higher by a factor of $\sim$5 on average, compared to other well-studied protostars in the Perseus and Ophiuchus clouds. The 2.9 mm visibility amplitudes for 6 of the 14 PBRS are very flat as a function of uv-distance, with more than 50\% of the source emission arising from radii $<$ 1500 AU. These flat visibility amplitudes are most consistent with spherically symmetric envelope density profiles with $ρ$~$\propto$~R$^{-2.5}$. Alternatively, there could be a massive unresolved structure like a disk or a high-density inner envelope departing from a smooth power-law. The large amount of mass on scales $<$ 1500 AU (implying high average central densities) leads us to suggest that that the PBRS with flat visibility amplitude profiles are the youngest PBRS and may be undergoing a brief phase of high mass infall/accretion and are possibly among the youngest Class 0 protostars. The PBRS with more rapidly declining visibility amplitudes still have large envelope masses, but could be slightly more evolved.

preprint2014arXiv

Line Profiles of Cores within Clusters. III. What is the most reliable tracer of core collapse in dense clusters?

Recent observational and theoretical investigations have emphasised the importance of filamentary networks within molecular clouds as sites of star formation. Since such environments are more complex than those of isolated cores, it is essential to understand how the observed line profiles from collapsing cores with non-spherical geometry are affected by filaments. In this study, we investigate line profile asymmetries by performing radiative transfer calculations on hydrodynamic models of three collapsing cores that are embedded in filaments. We compare the results to those that are expected for isolated cores. We model the five lowest rotational transition line (J = 1-0, 2-1, 3-2, 4-3, and 5-4) of both optically thick (HCN, HCO$^+$) as well as optically thin (N$_2$H$^+$, H$^{13}$CO$^+$) molecules using constant abundance laws. We find that less than 50% of simulated (1-0) transition lines show blue infall asymmetries due to obscuration by the surrounding filament. However, the fraction of collapsing cores that have a blue asymmetric emission line profile rises to 90% when observed in the (4-3) transition. Since the densest gas towards the collapsing core can excite higher rotational states, upper level transitions are more likely to produce blue asymmetric emission profiles. We conclude that even in irregular, embedded cores one can trace infalling gas motions with blue asymmetric line profiles of optically thick lines by observing higher transitions. The best tracer of collapse motions of our sample is the (4-3) transition of HCN, but the (3-2) and (5-4) transitions of both HCN and HCO$^+$ are also good tracers.

preprint2014arXiv

Molecular Outflows Driven by Low-Mass Protostars. I. Correcting for Underestimates When Measuring Outflow Masses and Dynamical Properties

We present a survey of 28 molecular outflows driven by low-mass protostars, all of which are sufficiently isolated spatially and/or kinematically to fully separate into individual outflows. Using a combination of new and archival data from several single-dish telescopes, 17 outflows are mapped in CO (2-1) and 17 are mapped in CO (3-2), with 6 mapped in both transitions. For each outflow, we calculate and tabulate the mass, momentum, kinetic energy, mechanical luminosity, and force assuming optically thin emission in LTE at an excitation temperature of 50 K. We show that all of the calculated properties are underestimated when calculated under these assumptions. Taken together, the effects of opacity, outflow emission at low velocities confused with ambient cloud emission, and emission below the sensitivities of the observations increase outflow masses and dynamical properties by an order of magnitude, on average, and factors of 50-90 in the most extreme cases. Different (and non-uniform) excitation temperatures, inclination effects, and dissociation of molecular gas will all work to further increase outflow properties. Molecular outflows are thus almost certainly more massive and energetic than commonly reported. Additionally, outflow properties are lower, on average, by almost an order of magnitude when calculated from the CO (3-2) maps compared to the CO (2-1) maps, even after accounting for different opacities, map sensitivities, and possible excitation temperature variations. It has recently been argued in the literature that the CO (3-2) line is subthermally excited in outflows, and our results support this finding.

preprint2014arXiv

The Evolution of Protostars: Insights from Ten Years of Infrared Surveys with Spitzer and Herschel

Stars form from the gravitational collapse of dense molecular cloud cores. In the protostellar phase, mass accretes from the core onto a protostar, likely through an accretion disk, and it is during this phase that the initial masses of stars and the initial conditions for planet formation are set. Over the past decade, new observational capabilities provided by the Spitzer Space Telescope and Herschel Space Observatory have enabled wide-field surveys of entire star-forming clouds with unprecedented sensitivity, resolution, and infrared wavelength coverage. We review resulting advances in the field, focusing both on the observations themselves and the constraints they place on theoretical models of star formation and protostellar evolution. We also emphasize open questions and outline new directions needed to further advance the field.

preprint2013arXiv

A Herschel and APEX Census of the Reddest Sources in Orion: Searching for the Youngest Protostars

We perform a census of the reddest, and potentially youngest, protostars in the Orion molecular clouds using data obtained with the PACS instrument onboard the Herschel Space Observatory and the LABOCA and SABOCA instruments on APEX as part of the Herschel Orion Protostar Survey (HOPS). A total of 55 new protostar candidates are detected at 70 um and 160 um that are either too faint (m24 > 7 mag) to be reliably classified as protostars or undetected in the Spitzer/MIPS 24 um band. We find that the 11 reddest protostar candidates with log (lambda F_lambda 70) / (lambda F_lambda 24) > 1.65 are free of contamination and can thus be reliably explained as protostars. The remaining 44 sources have less extreme 70/24 colors, fainter 70 um fluxes, and higher levels of contamination. Taking the previously known sample of Spitzer protostars and the new sample together, we find 18 sources that have log (lambda F_lambda 70) / (lambda F_lambda 24) > 1.65; we name these sources "PACS Bright Red sources", or PBRs. Our analysis reveals that the PBRs sample is composed of Class 0 like sources characterized by very red SEDs (T_bol < 45 K) and large values of sub-millimeter fluxes (L_smm/L_bol > 0.6%). Modified black-body fits to the SEDs provide lower limits to the envelope masses of 0.2 M_sun to 2 M_sun and luminosities of 0.7 L_sun to 10 L_sun. Based on these properties, and a comparison of the SEDs with radiative transfer models of protostars, we conclude that the PBRs are most likely extreme Class 0 objects distinguished by higher than typical envelope densities and hence, high mass infall rates.

preprint2013arXiv

G048.66-0.29: Physical State of an Isolated Site of Massive Star Formation

We present continuum observations of the infrared dark cloud (IRDC) G48.66-0.22 (G48) obtained with Herschel, Spitzer, and APEX, in addition to several molecular line observations. The Herschel maps are used to derive temperature and column density maps of G48 using a model based on a modified blackbody. We find that G48 has a relatively simple structure and is relatively isolated; thus this IRDC provides an excellent target to study the collapse and fragmentation of a filamentary structure in the absence of complicating factors such as strong external feedback. The derived temperature structure of G48 is clearly non-isothermal from cloud to core scale. The column density peaks are spatially coincident with the lowest temperatures (~ 17.5 K) in G48. A total cloud mass of ~390Msun is derived from the column density maps. By comparing the luminosity-to-mass ratio of 13 point sources detected in the Herschel/PACS bands to evolutionary models, we find that two cores are likely to evolve into high-mass stars (M>8 Msun). The derived mean projected separation of point sources is smaller than in other IRDCs but in good agreement with theoretical predications for cylindrical collapse. We detect several molecular species such as CO, HCO+, HCN, HNC and N2H+. CO is depleted by a factor of ~3.5 compared to the expected interstellar abundance, from which we conclude that CO freezes out in the central region. Furthermore, the molecular clumps, associated with the sub-millimeter peaks in G48, appear to be gravitationally unbound or just pressure confined. The analysis of critical line masses in G48 show that the entire filament is collapsing, overcoming any internal support.

preprint2013arXiv

HOPS 136: An Edge-On Orion Protostar Near the End of Envelope Infall

Edge-on protostars are valuable for understanding the disk and envelope properties of embedded young stellar objects, since the disk, envelope, and envelope cavities are all distinctly visible in resolved images and well constrained in modeling. Comparing 2MASS, WISE, Spitzer, Herschel, APEX, and IRAM photometry and limits from 1.25 to 1200 microns, Spitzer spectroscopy from 5 to 40 microns, and high-resolution Hubble imaging at 1.60 and 2.05 microns to radiative transfer modeling, we determine envelope and disk properties for the Class I protostar HOPS 136, an edge-on source in Orion's Lynds 1641 region. The source has a bolometric luminosity of 0.8 L_sun, a bolometric temperature of 170 K, and a ratio of submillimeter to bolometric luminosity of 0.8%. Via modeling, we find a total luminosity of 4.7 L_sun (larger than the observed luminosity due to extinction by the disk), an envelope mass of 0.06 M_sun, and a disk radius and mass of 450 AU and 0.002 M_sun. The stellar mass is highly uncertain but is estimated to fall between 0.4 and 0.5 M_sun. To reproduce the flux and wavelength of the near-infrared scattered-light peak in the spectral energy distribution, we require 5.4 x 10^-5 M_sun of gas and dust in each cavity. The disk has a large radius and a mass typical of more evolved T Tauri stars in spite of the significant remaining envelope. HOPS 136 appears to be a key link between the protostellar and optically revealed stages of star formation.

preprint2012arXiv

Dust SEDs in the era of Herschel and Planck: a Hierarchical Bayesian fitting technique

We present a hierarchical Bayesian method for fitting infrared spectral energy distributions (SEDs) of dust emission to observed fluxes. Under the standard assumption of optically thin single temperature (T) sources the dust SED as represented by a power--law modified black body is subject to a strong degeneracy between T and the spectral index beta. The traditional non-hierarchical approaches, typically based on chi-square minimization, are severely limited by this degeneracy, as it produces an artificial anti-correlation between T and beta even with modest levels of observational noise. The hierarchical Bayesian method rigorously and self-consistently treats measurement uncertainties, including calibration and noise, resulting in more precise SED fits. As a result, the Bayesian fits do not produce any spurious anti-correlations between the SED parameters due to measurement uncertainty. We demonstrate that the Bayesian method is substantially more accurate than the chi-square fit in recovering the SED parameters, as well as the correlations between them. As an illustration, we apply our method to Herschel and sub millimeter ground-based observations of the star-forming Bok globule CB244. This source is a small, nearby molecular cloud containing a single low-mass protostar and a starless core. We find that T and beta are weakly positively correlated -- in contradiction with the chi-square fits, which indicate a T-beta anti-correlation from the same data-set. Additionally, in comparison to the chi-square fits the Bayesian SED parameter estimates exhibit a reduced range in values.

preprint2012arXiv

KELT-1b: A Strongly Irradiated, Highly Inflated, Short Period, 27 Jupiter-mass Companion Transiting a mid-F Star

We present the discovery of KELT-1b, the first transiting low-mass companion from the wide-field Kilodegree Extremely Little Telescope-North (KELT-North) survey. The V=10.7 primary is a mildly evolved, solar-metallicity, mid-F star. The companion is a low-mass brown dwarf or super-massive planet with mass of 27.23+/-0.50 MJ and radius of 1.110+0.037-0.024 RJ, on a very short period (P=1.21750007) circular orbit. KELT-1b receives a large amount of stellar insolation, with an equilibrium temperature assuming zero albedo and perfect redistribution of 2422 K. Upper limits on the secondary eclipse depth indicate that either the companion must have a non-zero albedo, or it must experience some energy redistribution. Comparison with standard evolutionary models for brown dwarfs suggests that the radius of KELT-1b is significantly inflated. Adaptive optics imaging reveals a candidate stellar companion to KELT-1, which is consistent with an M dwarf if bound. The projected spin-orbit alignment angle is consistent with zero stellar obliquity, and the vsini of the primary is consistent with tidal synchronization. Given the extreme parameters of the KELT-1 system, we expect it to provide an important testbed for theories of the emplacement and evolution of short-period companions, and theories of tidal dissipation and irradiated brown dwarf atmospheres.

preprint2012arXiv

KELT-2Ab: A Hot Jupiter Transiting the Bright (V=8.77) Primary Star of a Binary System

We report the discovery of KELT-2Ab, a hot Jupiter transiting the bright (V=8.77) primary star of the HD 42176 binary system. The host is a slightly evolved late F-star likely in the very short-lived "blue-hook" stage of evolution, with $\teff=6148\pm48{\rm K}$, $\log{g}=4.030_{-0.026}^{+0.015}$ and $\feh=0.034\pm0.78$. The inferred stellar mass is $M_*=1.314_{-0.060}^{+0.063}$\msun\ and the star has a relatively large radius of $R_*=1.836_{-0.046}^{+0.066}$\rsun. The planet is a typical hot Jupiter with period $4.11379\pm0.00001$ days and a mass of $M_P=1.524\pm0.088$\mj\ and radius of $R_P=1.290_{-0.050}^{+0.064}$\rj. This is mildly inflated as compared to models of irradiated giant planets at the $\sim$4 Gyr age of the system. KELT-2A is the third brightest star with a transiting planet identified by ground-based transit surveys, and the ninth brightest star overall with a transiting planet. KELT-2Ab's mass and radius are unique among the subset of planets with $V<9$ host stars, and therefore increases the diversity of bright benchmark systems. We also measure the relative motion of KELT-2A and -2B over a baseline of 38 years, robustly demonstrating for the first time that the stars are bound. This allows us to infer that KELT-2B is an early K-dwarf. We hypothesize that through the eccentric Kozai mechanism KELT-2B may have emplaced KELT-2Ab in its current orbit. This scenario is potentially testable with Rossiter-McLaughlin measurements, which should have an amplitude of $\sim$44 m s$^{-1}$.

preprint2012arXiv

Line Profiles of Cores within Clusters: I. The Anatomy of a Filament

Observations are revealing the ubiquity of filamentary structures in molecular clouds. As cores are often embedded in filaments, it is important to understand how line profiles from such systems differ from those of isolated cores. We perform radiative transfer calculations on a hydrodynamic simulation of a molecular cloud in order to model line emission from collapsing cores embedded in filaments. We model two optically thick lines, CS(2-1) and HCN(1-0), and one optically thin line, N2H+(1-0), from three embedded cores. In the hydrodynamic simulation, gas self-gravity, turbulence, and bulk flows create filamentary regions within which cores form. Though the filaments have large dispersions, the N2H+(1-0) lines indicate subsonic velocities within the cores. We find that the observed optically thick line profiles of CS(2-1) and HCN(1-0) vary drastically with viewing angle. In over 50% of viewing angles, there is no sign of a blue asymmetry, an idealised signature of infall motions in an isolated spherical collapsing core. Profiles which primarily trace the cores, with little contribution from the surrounding filament, are characterised by a systematically higher HCN(1-0) peak intensity. The N2H+(1-0) lines do not follow this trend. We demonstrate that red asymmetric profiles are also feasible in the optically thick lines, due to emission from the filament or one-sided accretion flows onto the core. We conclude that embedded cores may frequently undergo collapse without showing a blue asymmetric profile, and that observational surveys including filamentary regions may underestimate the number of collapsing cores if based solely on profile shapes of optically thick lines.

preprint2012arXiv

Multiwavelength Observations of V2775 Ori, an Outbursting Protostar in L 1641: Exploring the Edge of the FU Orionis Regime

Individual outbursting young stars are important laboratories for studying the physics of episodic accretion and the extent to which this phenomenon can explain the luminosity distribution of protostars. We present new and archival data for V2775 Ori (HOPS 223), a protostar in the L 1641 region of the Orion molecular clouds that was discovered by Caratti o Garatti et al. (2011) to have recently undergone an order-of-magnitude increase in luminosity. Our near-infrared spectra of the source have strong blueshifted He I 10830 absorption, strong H2O and CO absorption, and no H I emission, all typical of FU Orionis sources. With data from IRTF, 2MASS, HST, Spitzer, WISE, Herschel, and APEX that span from 1 to 70 microns pre-outburst and from 1 to 870 microns post-outburst, we estimate that the outburst began between 2005 April and 2007 March. We also model the pre- and post-outburst spectral energy distributions of the source, finding it to be in the late stages of accreting its envelope with a disk-to-star accretion rate that increased from about 2x10^-6 M_sun/yr to about 10^-5 M_sun/yr during the outburst. The post-outburst luminosity at the epoch of the FU Orionis-like near-IR spectra is 28 L_sun, making V2775 Ori the least luminous documented FU Orionis outburster with a protostellar envelope. The existence of low-luminosity outbursts supports the notion that a range of episiodic accretion phenomena can partially explain the observed spread in protostellar luminosities.

preprint2010arXiv

Observational Constraints on Submillimeter Dust Opacity

Infrared extinction maps and submillimeter dust continuum maps are powerful probes of the density structure in the envelope of star-forming cores. We make a direct comparison between infrared and submillimeter dust continuum observations of the low-mass Class 0 core, B335, to constrain the ratio of submillimeter to infrared opacity (\kaprat) and the submillimeter opacity power-law index ($κ\propto λ^{-β}$). Using the average value of theoretical dust opacity models at 2.2 \micron, we constrain the dust opacity at 850 and 450 \micron . Using new dust continuum models based upon the broken power-law density structure derived from interferometric observations of B335 and the infall model derived from molecular line observations of B335, we find that the opacity ratios are $\frac{κ_{850}}{κ_{2.2}} = (3.21 - 4.80)^{+0.44}_{-0.30} \times 10^{-4}$ and $\frac{κ_{450}}{κ_{2.2}} = (12.8 - 24.8)^{+2.4}_{-1.3} \times 10^{-4}$ with a submillimeter opacity power-law index of $β_{smm} = (2.18 - 2.58)^{+0.30}_{-0.30}$. The range of quoted values are determined from the uncertainty in the physical model for B335. For an average 2.2 \micron\ opacity of $3800 \pm 700$ cm$^2$g$^{-1}$, we find a dust opacity at 850 and 450 \micron\ of $κ_{850} = (1.18 - 1.77)^{+0.36}_{-0.24}$ and $κ_{450} = (4.72 - 9.13)^{+1.9}_{-0.98}$ cm$^2$g$^{-1}$ of dust. These opacities are from $(0.65 - 0.97) κ^{\rm{OH}5}_{850}$ of the widely used theoretical opacities of Ossenkopf and Henning for coagulated ice grains with thin mantles at 850\micron.

preprint2010arXiv

The Spitzer c2d Survey of Nearby Dense Cores. IX. Discovery of a Very Low Luminosity Object Driving a Molecular Outflow in the Dense Core L673-7

We present new infrared, submillimeter, and millimeter observations of the dense core L673-7 and report the discovery of a low-luminosity, embedded Class 0 protostar driving a molecular outflow. L673-7 is seen in absorption against the mid-infrared background in 5.8, 8, and 24 micron Spitzer images, allowing for a derivation of the column density profile and total enclosed mass of L673-7, independent of dust temperature assumptions. Estimates of the core mass from these absorption profiles range from 0.2-4.5 solar masses. Millimeter continuum emission indicates a mass of about 2 solar masses, both from a direct calculation assuming isothermal dust and from dust radiative transfer models constrained by the millimeter observations. We use dust radiative transfer models to constrain the internal luminosity of L673-7, defined to be the luminosity of the central source and excluding the luminosity from external heating, to be 0.01-0.045 solar luminosities, with 0.04 solar luminosities the most likely value. L673-7 is thus classified as a very low luminosity object (VeLLO), and is among the lowest luminosity VeLLOs yet studied. We calculate the kinematic and dynamic properties of the molecular outflow in the standard manner, and we show that the expected accretion luminosity based on these outflow properties is greater than or equal to 0.36 solar luminosities. The discrepancy between this expected accretion luminosity and the internal luminosity derived from dust radiative transfer models indicates that the current accretion rate is much lower than the average rate over the lifetime of the outflow. Although the protostar embedded within L673-7 is consistent with currently being substellar, it is unlikely to remain as such given the substantial mass reservoir remaining in the core.

preprint2009arXiv

Spitzer and HHT observations of starless cores: masses and environments

We present Spitzer observations of a sample of 12 starless cores selected to have prominent 24 micron shadows. The Spitzer images show 8 and 24 micron shadows and in some cases 70 micron shadows; these spatially resolved absorption features trace the densest regions of the cores. We have carried out a 12CO (2-1) and 13CO (2-1) mapping survey of these cores with the Heinrich Hertz Telescope (HHT). We use the shadow features to derive optical depth maps. We derive molecular masses for the cores and the surrounding environment; we find that the 24 micron shadow masses are always greater than or equal to the molecular masses derived in the same region, a discrepancy likely caused by CO freeze--out onto dust grains. We combine this sample with two additional cores that we studied previously to bring the total sample to 14 cores. Using a simple Jeans mass criterion we find that ~ 2/3 of the cores selected to have prominent 24 micron shadows are collapsing or near collapse, a result that is supported by millimeter line observations. Of this subset at least half have indications of 70 micron shadows. All cores observed to produce absorption features at 70 micron are close to collapse. We conclude that 24 micron shadows, and even more so the 70 micron ones, are useful markers of cloud cores that are approaching collapse.