Source author record

Michael M. Dunham

Michael M. Dunham 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

61works
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

61 published item(s)

preprint2022arXiv

Ice features of low-luminosity protostars in near-infrared spectra of AKARI/IRC

We present near-infrared spectra of three low-luminosity protostars and one background star in the Perseus molecular cloud, acquired using the Infrared Camera (IRC) onboard the \textit{AKARI} space telescope. For the comparison with different star-forming environments, we also present spectra of the massive protostar AFGL 7009S, where the protostellar envelope is heated significantly, and the low-mass protostar RNO 91, which is suspected to be undergoing an episodic burst. We detected ice absorption features of \ch{H2O}, \ch{CO2}, and \ch{CO} at all spectra around the wavelengths of 3.05, 4.27, and 4.67 $μ$m, respectively. At least two low-luminosity protostars, we also detected the \ch{XCN} ice feature at 4.62 $μ$m. The presence of the crystalline \ch{H2O} ice and \ch{XCN} ice components indicates that the low-luminosity protostars experienced a hot phase via accretion bursts during the past mass accretion process. We compared the ice abundances of low-luminosity protostars with those of the embedded low-mass protostars and the dense molecular clouds and cores, suggesting that their ice abundances reflect the strength of prior bursts and the timescale after the last burst.

preprint2021arXiv

Evolution and Kinematics of Protostellar Envelopes in the Perseus Molecular Cloud

We present a comprehensive analysis on the evolution of envelopes surrounding protostellar systems in the Perseus molecular cloud using data from the MASSES survey. We focus our attention to the C$^{18}$O(2--1) spectral line, and we characterize the shape, size, and orientation of 54 envelopes and measure their fluxes, velocity gradients, and line widths. To look for evolutionary trends, we compare these parameters to the bolometric temperature Tbol, a tracer of protostellar age. We find evidence that the angular difference between the elongation angle of the C$^{18}$O envelope and the outflow axis direction generally becomes increasingly perpendicular with increasing Tbol, suggesting the envelope evolution is directly affected by the outflow evolution. We show that this angular difference changes at Tbol = $53 \pm 20$ K, which includes the conventional delineation between Class 0 and I protostars of 70K. We compare the C$^{18}$O envelopes with larger gaseous structures in other molecular clouds and show that the velocity gradient increases with decreasing radius ($|\mathcal{G}| \sim R^{-0.72 \pm 0.06}$). From the velocity gradients we show that the specific angular momentum follows a power law fit $J/M \propto R^{1.83 \pm 0.05}$ for scales from 1pc down to $\sim$500 au, and we cannot rule out a possible flattening out at radii smaller than $\sim$1000 au.

preprint2020arXiv

A Study of Millimeter Variability in FUor Objects

FU Orionis objects (FUors) are rapidly-accreting, pre-main sequence objects that are known to exhibit large outbursts at optical and near-infrared wavelengths, with post-eruption, small-scale photometric variability superimposed on longer-term trends. In contrast, little is known about the variability of FUors at longer wavelengths. To explore this further, we observed six FUor objects using the NOrthern Extended Millimeter Array (NOEMA) and for a subset of three objects we obtained coordinated observations with NOEMA and the Lowell Discovery Telescope (LDT). In combination with previously published NOEMA observations from 2014, our 2017 observations of V1735 Cyg provide the first detection of variability in an FUor object at 2.7 mm. In the absence of significant optical variability, we discount the possibility that the mm flux density changed as a result of irradiation from the central disk. In addition, a change in the dust mass due to infall is highly unlikely. A plausible explanation for the change in 2.7 mm flux density is variability in free-free emission due to changes in the object's jet/wind. Thus, it may be that free-free emission in some FUor objects is significant at $\sim$3 mm and must be considered when deriving disk masses in order to help constrain the mechanism responsible for triggering FUor outbursts.

preprint2020arXiv

Kinematic Analysis of a Protostellar Multiple System: Measuring the Protostar Masses and Assessing Gravitational Instability in the Disks of L1448 IRS3B and L1448 IRS3A

We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations towards a compact (230~au separation) triple protostar system, L1448 IRS3B, at 879~\micron with \contbeam~resolution. Spiral arm structure within the circum-multiple disk is well resolved in dust continuum toward IRS3B, and we detect the known wide (2300~au) companion, IRS3A, also resolving possible spiral substructure. Using dense gas tracers, C17O, H13CO$+$, and H13CN, we resolve the Keplerian rotation for both the circum-triple disk in IRS3B and the disk around IRS3A. Furthermore, we use the molecular line kinematic data and radiative transfer modeling of the molecular line emission to confirm that the disks are in Keplerian rotation with fitted masses of $1.19^{+0.13}_{-0.07}$ for IRS3B-ab, $1.51^{+0.06}_{-0.07}$~Msun for IRS3A, and place an upper limit on the central protostar mass for the tertiary IRS3B-c of 0.2~Msun. We measure the mass of the fragmenting disk of IRS3B to be 0.29~Msun from the dust continuum emission of the circum-multiple disk and estimate the mass of the clump surrounding IRS3B-c to be 0.07~Msun. We also find that the disk around IRS3A has a mass of 0.04~Msun. By analyzing the Toomre~Q parameter, we find the IRS3A circumstellar disk is gravitationally stable (Q$>$5), while the IRS3B disk is consistent with a gravitationally unstable disk (Q$<$1) between the radii 200-500~au. This coincides with the location of the spiral arms and the tertiary companion IRS3B-c, supporting the hypothesis that IRS3B-c was formed in situ via fragmentation of a gravitationally unstable disk.

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 CARMA-NRO Orion Survey: Protostellar Outflows, Energetics, and Filamentary Alignment

We identify 45 protostellar outflows in CO maps of the Orion A giant molecular cloud from the CARMA-NRO Orion survey. Our sample includes 11 newly detected outflows. We measure the mass and energetics of the outflows, including material at low-velocities by correcting for cloud contributions. The total momentum and kinetic energy injection rates of outflows is comparable to the turbulent dissipation rate of the cloud. We also compare the outflow position angles to the orientation of C$^{18}$O filaments. We find that the full sample of outflows is consistent with being randomly oriented with respect to the filaments. A subsample of the most reliable measurements shows a moderately perpendicular outflow-filament alignment which may reflect accretion of mass across filaments and onto the protostellar cores.

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

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'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 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

A Catalog of Low-Mass Star-Forming Cores Observed with SHARC-II at 350 microns

We present a catalog of low-mass dense cores observed with the SHARC-II instrument at 350 microns. Our observations have an effective angular resolution of 10", approximately 2.5 times higher than observations at the same wavelength obtained with the Herschel Space Observatory, albeit with lower sensitivity, especially to extended emission. The catalog includes 81 maps covering a total of 164 detected sources. For each detected source, we tabulate basic source properties including position, peak intensity, flux density in fixed apertures, and radius. We examine the uncertainties in the pointing model applied to all SHARC-II data and conservatively find that the model corrections are good to within ~3", approximately 1/3 of the SHARC-II beam. We examine the differences between two array scan modes and find that the instrument calibration, beam size, and beam shape are similar between the two modes. We also show that the same flux densities are measured when sources are observed in the two different modes, indicating that there are no systematic effects introduced into our catalog by utilizing two different scan patterns during the course of taking observations. We find a detection rate of 95% for protostellar cores but only 45% for starless cores, and demonstrate the existence of a SHARC-II detection bias against all but the most massive and compact starless cores. Finally, we discuss the improvements in protostellar classification enabled by these 350 micron observations.

preprint2016arXiv

An ALMA Search for Substructure, Fragmentation, and Hidden Protostars in Starless Cores in Chamaeleon I

We present an Atacama Large Millimeter/submillimeter Array (ALMA) 106 GHz (Band 3) continuum survey of the complete population of dense cores in the Chamaeleon I molecular cloud. We detect a total of 24 continuum sources in 19 different target fields. All previously known Class 0 and Class I protostars in Chamaeleon I are detected, whereas all of the 56 starless cores in our sample are undetected. We show that the Spitzer+Herschel census of protostars in Chamaeleon I is complete, with the rate at which protostellar cores have been misclassified as starless cores calculated as <1/56, or < 2%. We use synthetic observations to show that starless cores collapsing following the turbulent fragmentation scenario are detectable by our ALMA observations when their central densities exceed ~10^8 cm^-3, with the exact density dependent on the viewing geometry. Bonnor-Ebert spheres, on the other hand, remain undetected to central densities at least as high as 10^10 cm^-3. Our starless core non-detections are used to infer that either the star formation rate is declining in Chamaeleon I and most of the starless cores are not collapsing, matching the findings of previous studies, or that the evolution of starless cores are more accurately described by models that develop less substructure than predicted by the turbulent fragmentation scenario, such as Bonnor-Ebert spheres. We outline future work necessary to distinguish between these two possibilities.

preprint2016arXiv

An extremely high velocity molecular jet surrounded by an ionized cavity in the protostellar source Serpens SMM1

We report ALMA observations of a one-sided, high-velocity ($\sim$80 km s$^{-1}$) CO($J = 2 \rightarrow 1$) jet powered by the intermediate-mass protostellar source Serpens SMM1-a. The highly collimated molecular jet is flanked at the base by a wide-angle cavity; the walls of the cavity can be seen in both 4 cm free-free emission detected by the VLA and 1.3 mm thermal dust emission detected by ALMA. This is the first time that ionization of an outflow cavity has been directly detected via free-free emission in a very young, embedded Class 0 protostellar source that is still powering a molecular jet. The cavity walls are ionized either by UV photons escaping from the accreting protostellar source, or by the precessing molecular jet impacting the walls. These observations suggest that ionized outflow cavities may be common in Class 0 protostellar sources, shedding further light on the radiation, outflow, and jet environments in the youngest, most embedded forming stars.

preprint2016arXiv

Circumstellar Disks of the Most Vigorously Accreting Young Stars

Young stellar objects (YSOs) may not accumulate their mass steadily, as was previously thought, but in a series of violent events manifesting themselves as sharp stellar brightening. These events can be caused by fragmentation due to gravitational instabilities in massive gaseous disks surrounding young stars, followed by migration of dense gaseous clumps onto the star. We report our high angular resolution, coronagraphic near-infrared polarization imaging observations using the High Contrast Instrument for the Subaru Next Generation Adaptive Optics (HiCIAO) of the Subaru 8.2 m Telescope, towards four YSOs which are undergoing luminous accretion outbursts. The obtained infrared images have verified the presence of several hundred AUs scale arms and arcs surrounding these YSOs. In addition, our hydrodynamics simulations and radiative transfer models further demonstrate that these observed structures can indeed be explained by strong gravitational instabilities occurring at the beginning of the disk formation phase. The effect of those tempestuous episodes of disk evolution on star and planet formation remains to be understood.

preprint2016arXiv

Imaging the water-snow line during a protostellar outburst

A snow-line is the region of a protoplanetary disk at which a major volatile, such as water or carbon monoxide, reaches its condensation temperature. Snow-lines play a crucial role in disk evolution by promoting the rapid growth of ice-covered grains. Signatures of the carbon monoxide snow-line (at temperatures of around 20 kelvin) have recently been imaged in the disks surrounding the pre-main-sequence stars TW Hydra and HD163296, at distances of about 30 astronomical units (au) from the star. But the water snow-line of a protoplanetary disk (at temperatures of more than 100 kelvin) has not hitherto been seen, as it generally lies very close to the star (less than 5 au away for solar-type stars). Water-ice is important because it regulates the efficiency of dust and planetesimal coagulation, and the formation of comets, ice giants and the cores of gas giants. Here we report ALMA images at 0.03-arcsec resolution (12 au) of the protoplanetary disk around V883 Ori, a protostar of 1.3 solar masses that is undergoing an outburst in luminosity arising from a temporary increase in the accretion rate. We find an intensity break corresponding to an abrupt change in the optical depth at about 42 au, where the elevated disk temperature approaches the condensation point of water, from which we conclude that the outburst has moved the water snow-line. The spectral behaviour across the snow-line confirms recent model predictions: dust fragmentation and the inhibition of grain growth at higher temperatures results in soaring grain number densities and optical depths. As most planetary systems are expected to experience outbursts caused by accretion during their formation our results imply that highly dynamical water snow-lines must be considered when developing models of disk evolution and planet formation.

preprint2016arXiv

Misalignment of Outflow Axes in the Proto-Multiple Systems in Perseus

We investigate the alignment between outflow axes in nine of the youngest binary/multiple systems in the Perseus Molecular Cloud. These systems have typical member spacing larger than 1000 AU. For outflow identification, we use 12CO(2-1) and 12CO(3-2) data from a large survey with the Sub- millimeter Array: Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES). The distribution of outflow orientations in the binary pairs is consistent with random or preferentially anti-aligned distributions, demonstrating that these outflows are misaligned. This result suggests that these systems are possibly formed in environments where the distribution of angular momentum is complex and disordered, and these systems do not come from the same co-rotating structures or from an initial cloud with aligned vectors of angular momentum.

preprint2016arXiv

Revolutionizing our View of Protostellar Multiplicity and Disks: The VLA Nascent Disk and Multiplicity (VANDAM) Survey of the Perseus Molecular Cloud

There is substantial evidence for disk formation taking place during the early stages of star formation and for most stars being born in multiple systems; however, protostellar multiplicity and disk searches have been hampered by low resolution, sample bias, and variable sensitivity. We have conducted an unbiased, high-sensitivity Karl G. Jansky Very Large Array (VLA) survey toward all known protostars (n = 94) in the Perseus molecular cloud (d~230 pc), with a resolution of ~15 AU (0.06") at 8 mm. We have detected candidate protostellar disks toward 17 sources (with 12 of those in the Class 0 stage) and we have found substructure on < 50AU scales for three Class 0 disk candidates, possibly evidence for disk fragmentation. We have discovered 16 new multiple systems (or new components) in this survey; the new systems have separations < 500 AU and 3 by < 30 AU. We also found a bi-modal distribution of separations, with peaks at ~75 AU and ~3000 AU, suggestive of formation through two distinct mechanisms: disk and turbulent fragmentation. The results from this survey demonstrate the necessity and utility of uniform, unbiased surveys of protostellar systems at millimeter and centimeter wavelengths.

preprint2016arXiv

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 momenta. Although the pairs migrate to closer separations, their spins remain partially misaligned. We produce $^{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 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.

preprint2016arXiv

The VLA Nascent Disk and Multiplicity Survey of Perseus Protostars (VANDAM). II. Multiplicity of Protostars in the Perseus Molecular Cloud

We present a multiplicity study of all known protostars (94) in the Perseus molecular cloud from a Karl G. Jansky Very Large Array (VLA) survey at Ka-band (8 mm and 1 cm) and C-band (4 cm and 6.6 cm). The observed sample has a bolometric luminosity range between 0.1 L$_{\odot}$ and $\sim$33 L$_{\odot}$, with a median of 0.7 L$_{\odot}$. This multiplicity study is based on the Ka-band data, having a best resolution of $\sim$0.065" (15 AU) and separations out to $\sim$43" (10000 AU) can be probed. The overall multiplicity fraction (MF) is found to be of 0.40$\pm$0.06 and the companion star fraction (CSF) is 0.71$\pm$0.06. The MF and CSF of the Class 0 protostars are 0.57$\pm$0.09 and 1.2$\pm$0.2, and the MF and CSF of Class I protostars are both 0.23$\pm$0.08. The distribution of companion separations appears bi-modal, with a peak at $\sim$75 AU and another peak at $\sim$3000 AU. Turbulent fragmentation is likely the dominant mechanism on $>$1000 AU scales and disk fragmentation is likely to be the dominant mechanism on $<$200 AU scales. Toward three Class 0 sources we find companions separated by $<$30 AU. These systems have the smallest separations of currently known Class 0 protostellar binary systems. Moreover, these close systems are embedded within larger (50 AU to 400 AU) structures and may be candidates for ongoing disk fragmentation.

preprint2016arXiv

The VLA Nascent Disk and Multiplicity Survey: First Look at Resolved Candidate Disks around Class 0 and I Protostars in the Perseus Molecular Cloud

We present the first dust emission results toward a sample of seven protostellar disk candidates around Class 0 and I sources in the Perseus molecular cloud from the VLA Nascent Disk and Multiplicity (VANDAM) survey with ~0.05'' or 12 AU resolution. To examine the surface brightness profiles of these sources, we fit the Ka-band 8 mm dust-continuum data in the u,v-plane to a simple, parametrized model based on the Shakura-Sunyaev disk model. The candidate disks are well-fit by a model with a disk-shaped profile and have masses consistent with known Class 0 and I disks. The inner-disk surface densities of the VANDAM candidate disks have shallower density profiles compared to disks around more evolved Class II systems. The best-fit model radii of the seven early-result candidate disks are R_c > 10 AU; at 8 mm, the radii reflect lower limits on the disk size since dust continuum emission is tied to grain size and large grains radially drift inwards. These relatively large disks, if confirmed kinematically, are inconsistent with theoretical models where the disk size is limited by strong magnetic braking to < 10 AU at early times.

preprint2015arXiv

Absence of Significant Cool Disks in Young Stellar Objects Exhibiting Repetitive Optical Outbursts

We report Submillimeter Array (SMA) 1.3 mm high angular resolution observations towards the four EXor type outbursting young stellar objects (YSOs) VY Tau, V1118 Ori, V1143 Ori, and NY Ori. The data mostly show low dust masses $M_{dust}$ in the associated circumstellar disks. Among the sources, NY Ori possesses a relatively massive disk with $M_{dust} \sim 9 \times 10^{-4}$ $M_{\odot}$. V1118 Ori has a marginal detection equivalent to $M_{dust} \sim 6 \times 10^{-5}$ $M_{\odot}$. V1143 Ori has a non-detection also equivalent to $M_{dust} < 6 \times 10^{-5}$ $M_{\odot}$. For the nearest source VY Tau, we get a surprising non-detection which provides a stringent upper limit $M_{dust} < 6 \times 10^{-6}$ $M_{\odot}$. We interpret our findings as suggesting that the gas and dust reservoirs that feed the short duration, repetitive optical outbursts seen in some EXors may be limited to the small scale, innermost region of their circumstellar disks. This hot dust may have escaped our detection limits. Follow-up, more sensitive millimeter observations are needed to improve our understanding of the triggering mechanisms of EXor type outbursts.

preprint2015arXiv

Assessing molecular outflows and turbulence in the protostellar cluster Serpens South

Molecular outflows driven by protostellar cluster members likely impact their surroundings and contribute to turbulence, affecting subsequent star formation. The very young Serpens South cluster consists of a particularly high density and fraction of protostars, yielding a relevant case study for protostellar outflows and their impact on the cluster environment. We combined CO $J=1-0$ observations of this region using the Combined Array for Research in Millimeter-wave Astronomy (CARMA) and the Institut de Radioastronomie Millimétrique (IRAM) 30 m single dish telescope. The combined map allows us to probe CO outflows within the central, most active region at size scales of 0.01 pc to 0.8 pc. We account for effects of line opacity and excitation temperature variations by incorporating $^{12}$CO and $^{13}$CO data for the $J=1-0$ and $J=3-2$ transitions (using Atacama Pathfinder Experiment and Caltech Submillimeter Observatory observations for the higher CO transitions), and we calculate mass, momentum, and energy of the molecular outflows in this region. The outflow mass loss rate, force, and luminosity, compared with diagnostics of turbulence and gravity, suggest that outflows drive a sufficient amount of energy to sustain turbulence, but not enough energy to substantially counter the gravitational potential energy and disrupt the clump. Further, we compare Serpens South with the slightly more evolved cluster NGC 1333, and we propose an empirical scenario for outflow-cluster interaction at different evolutionary stages.

preprint2015arXiv

Episodic molecular outflow in the very young protostellar cluster Serpens South

The loss of mass from protostars, in the form of a jet or outflow, is a necessary counterpart to protostellar mass accretion. Outflow ejection events probably vary in their velocity and/or in the rate of mass loss. Such `episodic' ejection events have been observed during the Class 0 protostellar phase (the early accretion stage), and continue during the subsequent class I phase that marks the first one million years of star formation. Previously observed episodic-ejection sources were relatively isolated; however, the most common sites of star formation are clusters. Outflows link protostars with their environment and provide a viable source of turbulence that is necessary for regulating star formation in clusters, but it is not known how an accretion-driven jet or outflow in a clustered environment manifests itself in its earliest stage. This early stage is important in establishing the initial conditions for momentum and energy transfer to the environment as the protostar and cluster evolve. Here we report that an outflow from a very young class 0 protostar, at the hub of the very active and filamentary Serpens South protostellar cluster, shows unambiguous episodic events. The $^{12}$CO (J=2-1) emission from the protostar reveals 22 distinct features of outflow ejecta, the most recent having the highest velocity. The outflow forms bipolar lobes --- one of the first detectable signs of star formation --- which originate from the peak of 1-mm continuum emission. Emission from the surrounding C$^{18}$O envelope shows kinematics consistent with rotation and an infall of material onto the protostar. The data suggest that episodic accretion-driven outflow begins in the earliest phase of protostellar evolution, and that the outflow remains intact in a very clustered environment, probably providing efficient momentum transfer for driving turbulence.

preprint2015arXiv

High Resolution 8 mm and 1 cm Polarization of IRAS 4A from the VLA Nascent Disk and Multiplicity (VANDAM) Survey

Magnetic fields can regulate disk formation, accretion and jet launching. Until recently, it has been difficult to obtain high resolution observations of the magnetic fields of the youngest protostars in the critical region near the protostar. The VANDAM survey is observing all known protostars in the Perseus Molecular Cloud. Here we present the polarization data of IRAS 4A. We find that with ~ 0.2'' (50 AU) resolution at λ = 8.1 and 10.3 mm, the inferred magnetic field is consistent with a circular morphology, in marked contrast with the hourglass morphology seen on larger scales. This morphology is consistent with frozen-in field lines that were dragged in by rotating material entering the infall region. The field morphology is reminiscent of rotating circumstellar material near the protostar. This is the first polarization detection of a protostar at these wavelengths. We conclude from our observations that the dust emission is optically thin with β ~ 1.3, suggesting that mm/cm-sized grains have grown and survived in the short lifetime of the protostar.

preprint2015arXiv

Infrared and Radio observations of a small group of protostellar objects in the molecular core, L1251-C

We present a multi-wavelength observational study of a low-mass star-forming region, L1251-C, with observational results at wavelengths from the near-infrared to the millimeter. Spitzer Space Telescope observations confirmed that IRAS 22343+7501 is a small group of protostellar objects. The extended emission to east-west direction with its intensity peak at the center of L1251A has been detected at 350 and 850 micron with the CSO and JCMT telescopes, tracing dense envelope materials around L1251A. The single-dish data from the KVN and TRAO telescopes show inconsistencies between the intensity peaks of several molecular line emission and that of the continuum emission, suggesting complex distributions of molecular abundances around L1251A. The SMA interferometer data, however, show intensity peaks of CO 2-1 and 13CO 2-1 located at the position of IRS 1, which is both the brightest source in IRAC image and the weakest source in the 1.3 mm dust continuum map. IRS 1 is the strongest candidate for the driving source of the newly detected compact CO 2-1 outflow. Over the whole region (14' by 14') of L125l-C, 3 Class I and 16 Class II sources have been detected, including three YSOs in L1251A. A comparison with the average projected distance among 19 YSOs in L1251-C and that among 3 YSOs in L1251A suggests L1251-C is an example of low-mass cluster formation, where protostellar objects are forming in a small group.

preprint2015arXiv

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 GHz and 345 GHz from the Sub-millimeter 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" (MASSES). Three dust continuum sources, Source B, Source NW, and Source A, are detected at both frequencies. These sources have corresponding emission peaks in C18O (J=2-1), 13CO (J=2-1), and HCO+ (J=4-3), and have offsets with N2D+ (J=3-2) peaks. High angular resolution data from a complimentary 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 C18O emission are found to be perpendicular to within 20 degrees of the outflow directions as revealed by 12CO (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.

preprint2015arXiv

The Bolocam Galactic Plane Survey. XI. Temperatures and Substructure of Galactic Clumps Based on 350 micron Observations

We present 107 maps of continuum emission at 350 microns from Galactic molecular clumps. Observed sources were mainly selected from the Bolocam Galactic Plane Survey (BGPS) catalog, with 3 additional maps covering star forming regions in the outer Galaxy. The higher resolution of the SHARC-II images (8.5" beam) compared with the 1.1 mm images from BGPS (33" beam) allowed us to identify a large population of smaller substructures within the clumps. A catalog is presented for the 1386 sources extracted from the 350 micron maps. The color temperature distribution of clumps based on the two wavelengths has a median of 13.3 K and mean of 16.3 +- 0.4 K, assuming an opacity law index of 1.7. For the structures with good determination of color temperatures, the mean ratio of gas temperature, determined from NH3 observations, to dust color temperature is 0.88 and the median ratio is 0.76. About half the clumps have more than two substructures and 22 clumps have more than 10. The fraction of the mass in dense substructures seen at 350 microns compared to the mass of their parental clump is ~0.19, and the surface densities of these substructures are, on average, 2.2 times those seen in the clumps identified at 1.1 mm. For a well-characterized sample, 88 structures (31%) exceed a surface density of 0.2 g cm^(-2), and 18 (6%) exceed 1.0 g cm^(-2), thresholds for massive star formation suggested by theorists.

preprint2015arXiv

The Spitzer c2d Survey of Large, Nearby, Interstellar Clouds. XII. The Perseus YSO Population as Observed with IRAC and MIPS

The Spitzer Space Telescope mapped the Perseus molecular cloud complex with IRAC and MIPS as part of the c2d Spitzer Legacy project. This paper combines the observations from both instruments giving an overview of low-mass star formation across Perseus from 3.6 to 70 micron. We provide an updated list of young stellar objects with new classifications and source fluxes from previous works, identifying 369 YSOs in Perseus with the Spitzer dataset. By synthesizing the IRAC and MIPS maps of Perseus and building on the work of previous papers in this series (Jorgensen et al. 2006, Rebull et al. 2007), we present a current census of star formation across the cloud and within smaller regions. 67% of the YSOs are associated with the young clusters NGC 1333 and IC 348. The majority of the star formation activity in Perseus occurs in the regions around the clusters, to the eastern and western ends of the cloud complex. The middle of the cloud is nearly empty of YSOs despite containing regions of high visual extinction. The western half of Perseus contains three-quarters of the total number of embedded YSOs (Class 0+I and Flat SED sources) in the cloud and nearly as many embedded YSOs as Class II and III sources. Class II and III greatly outnumber Class 0+I objects in eastern Perseus and IC 348. These results are consistent with previous age estimates for the clusters. Across the cloud, 56% of YSOs and 91% of the Class 0+I and Flat sources are in areas where Av > 5 mag, indicating a possible extinction threshold for star formation.

preprint2015arXiv

Young Stellar Objects in the Gould Belt

We present the full catalog of Young Stellar Objects (YSOs) identified in the 18 molecular clouds surveyed by the Spitzer Space Telescope "cores to disks" (c2d) and "Gould Belt" (GB) Legacy surveys. Using standard techniques developed by the c2d project, we identify 3239 candidate YSOs in the 18 clouds, 2966 of which survive visual inspection and form our final catalog of YSOs in the Gould Belt. We compile extinction corrected SEDs for all 2966 YSOs and calculate and tabulate the infrared spectral index, bolometric luminosity, and bolometric temperature for each object. We find that 326 (11%), 210 (7%), 1248 (42%), and 1182 (40%) are classified as Class 0+I, Flat-spectrum, Class II, and Class III, respectively, and show that the Class III sample suffers from an overall contamination rate by background AGB stars between 25% and 90%. Adopting standard assumptions, we derive durations of 0.40-0.78 Myr for Class 0+I YSOs and 0.26-0.50 Myr for Flat-spectrum YSOs, where the ranges encompass uncertainties in the adopted assumptions. Including information from (sub)millimeter wavelengths, one-third of the Class 0+I sample is classified as Class 0, leading to durations of 0.13-0.26 Myr (Class 0) and 0.27-0.52 Myr (Class I). We revisit infrared color-color diagrams used in the literature to classify YSOs and propose minor revisions to classification boundaries in these diagrams. Finally, we show that the bolometric temperature is a poor discriminator between Class II and Class III YSOs.

preprint2014arXiv

Episodic Accretion in Young Stars

In the last twenty years, the topic of episodic accretion has gained significant interest in the star formation community. It is now viewed as a common, though still poorly understood, phenomenon in low-mass star formation. The FU Orionis objects (FUors) are long-studied examples of this phenomenon. FUors are believed to undergo accretion outbursts during which the accretion rate rapidly increases from typically $10^{-7}$ to a few $10^{-4}$ $M_\odot$ yr$^{-1}$, and remains elevated over several decades or more. EXors, a loosely defined class of pre-main sequence stars, exhibit shorter and repetitive outbursts, associated with lower accretion rates. The relationship between the two classes, and their connection to the standard pre-main sequence evolutionary sequence, is an open question: do they represent two distinct classes, are they triggered by the same physical mechanism, and do they occur in the same evolutionary phases? Over the past couple of decades, many theoretical and numerical models have been developed to explain the origin of FUor and EXor outbursts. In parallel, such accretion bursts have been detected at an increasing rate, and as observing techniques improve each individual outburst is studied in increasing detail. We summarize key observations of pre-main sequence star outbursts, and review the latest thinking on outburst triggering mechanisms, the propagation of outbursts from star/disk to disk/jet systems, the relation between classical EXors and FUors, and newly discovered outbursting sources -- all of which shed new light on episodic accretion. We finally highlight some of the most promising directions for this field in the near- and long-term.

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

On the reliability of protostellar disc mass measurements and the existence of fragmenting discs

We couple non-magnetic, hydrodynamical simulations of collapsing protostellar cores with radiative transfer evolutionary models to generate synthetic observations. We then use these synthetic observations to investigate the extent to which a simple method for measuring protostellar disc masses used in the literature recovers the intrinsic masses of the discs formed in the simulations. We evaluate the effects of contamination from the surrounding core, partially resolving out the disc, optical depth, fixed assumed dust temperatures, inclination, and the dust opacity law. We show that the combination of these effects can lead to disc mass underestimates by up to factors of 2-3 at millimeter wavelengths and up to an order of magnitude or larger at submillimeter wavelengths. The optically thin portions of protostellar discs are generally cooler in the Class I stage than the Class 0 stage since Class I discs are typically larger and more optically thick, and thus more shielded. The observed disc mass distribution closely resembles the intrinsic distribution if this effect is taken into account, especially at millimeter wavelengths where optical depth effects are minimized. Approximately 50%-70% of protostellar discs observed to date with this method are consistent with the masses of the gravitationally unstable discs formed in the simulations, suggesting that at least some protostellar discs are likely sufficiently massive to fragment. We emphasize key future work needed to confirm these results, including assembling larger, less biased samples, and using molecular line observations to distinguish between rotationally supported, Keplerian discs and magnetically supported pseudodiscs.

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.

preprint2014arXiv

The Spitzer Survey of Interstellar Clouds in the Gould Belt. VI. The Auriga-California Molecular Cloud observed with IRAC and MIPS

We present observations of the Auriga-California Molecular Cloud (AMC) at 3.6, 4.5, 5.8, 8.0, 24, 70 and 160 micron observed with the IRAC and MIPS detectors as part of the Spitzer Gould Belt Legacy Survey. The total mapped areas are 2.5 sq-deg with IRAC and 10.47 sq-deg with MIPS. This giant molecular cloud is one of two in the nearby Gould Belt of star-forming regions, the other being the Orion A Molecular Cloud (OMC). We compare source counts, colors and magnitudes in our observed region to a subset of the SWIRE data that was processed through our pipeline. Using color-magnitude and color-color diagrams, we find evidence for a substantial population of 166 young stellar objects (YSOs) in the cloud, many of which were previously unknown. Most of this population is concentrated around the LkHalpha 101 cluster and the filament extending from it. We present a quantitative description of the degree of clustering and discuss the fraction of YSOs in the region with disks relative to an estimate of the diskless YSO population. Although the AMC is similar in mass, size and distance to the OMC, it is forming about 15 - 20 times fewer stars.

preprint2014arXiv

The VLA Nascent Disk And Multiplicity (VANDAM) Survey of Perseus Protostars. Resolving the Sub-Arcsecond Binary System in NGC 1333 IRAS2A

We are conducting a Jansky VLA Ka-band (8 mm and 1 cm) and C-band (4 cm and 6.4 cm) survey of all known protostars in the Perseus Molecular Cloud, providing resolution down to $\sim$0.06'' and $\sim$0.35" in Ka-band and C-band, respectively. Here we present first results from this survey that enable us to examine the source NGC 1333 IRAS2A in unprecedented detail and resolve it into a proto-binary system separated by 0.621"$\pm$0.006" ($\sim$143 AU) at 8 mm, 1 cm, and 4 cm. These 2 sources (IRAS2A VLA1 and VLA2) are likely driving the two orthogonal outflows known to originate from IRAS2A. The brighter source IRAS2A VLA1 is extended perpendicular to its outflow in the VLA data, with a deconvolved size of 0.055" ($\sim$13 AU), possibly tracing a protostellar disk. The recently reported candidate companions (IRAS2A MM2 and MM3) are not detected in either our VLA data, CARMA 1.3 mm data, or SMA 850 $μ$m data. SMA CO ($J=3\rightarrow2$), CARMA CO ($J=2\rightarrow1$), and lower resolution CARMA CO ($J=1\rightarrow0$) observations are used to examine the outflow origins and the nature of the candidate companions to IRAS2A VLA1. The CO ($J=3\rightarrow2$) and ($J=2\rightarrow1$) data show that IRAS2A MM2 is coincident with a bright CO emission spot in the east-west outflow, and IRAS2A MM3 is within the north-south outflow. In contrast, IRAS2A VLA2 lies at the east-west outflow symmetry point. We propose that IRAS2A VLA2 is the driving source of the East-West outflow and a true companion to IRAS2A VLA1, whereas IRAS2A MM2 and MM3 may not be protostellar.

preprint2013arXiv

An Analysis of the Environments of FU Orionis Objects with Herschel

We present Herschel-HIFI, SPIRE, and PACS 50-670 μm imaging and spectroscopy of six FU Orionis-type objects and candidates (FU Orionis, V1735 Cyg, V1515 Cyg, V1057 Cyg, V1331 Cyg, and HBC 722), ranging in outburst date from 1936-2010, from the "FOOSH" (FU Orionis Objects Surveyed with Herschel) program, as well as ancillary results from Spitzer-IRS and the Caltech Submillimeter Observatory. In their system properties (Lbol, Tbol, line emission), we find that FUors are in a variety of evolutionary states. Additionally, some FUors have features of both Class I and II sources: warm continuum consistent with Class II sources, but rotational line emission typical of Class I, far higher than Class II sources of similar mass/luminosity. Combining several classification techniques, we find an evolutionary sequence consistent with previous mid-IR indicators. We detect [O I] in every source at luminosities consistent with Class 0/I protostars, much greater than in Class II disks. We detect transitions of 13CO (J_up of 5 to 8) around two sources (V1735 Cyg and HBC 722) but attribute them to nearby protostars. Of the remaining sources, three (FU Ori, V1515 Cyg, and V1331 Cyg) exhibit only low-lying CO, but one (V1057 Cyg) shows CO up to J = 23 - 22 and evidence for H2O and OH emission, at strengths typical of protostars rather than T Tauri stars. Rotational temperatures for "cool" CO components range from 20-81 K, for ~ 10^50 total CO molecules. We detect [C I] and [N II] primarily as diffuse emission.

preprint2013arXiv

Embedded Protostars in the Dust, Ice, and Gas In Time (DIGIT) Key Program: Continuum SEDs, and an Inventory of Characteristic Far-Infrared Lines from PACS Spectroscopy

We present 50-210 um spectral scans of 30 Class 0/I protostellar sources, obtained with Herschel-PACS, and 0.5-1000 um SEDs, as part of the Dust, Ice, and Gas in Time (DIGIT) Key Program. Some sources exhibit up to 75 H2O lines ranging in excitation energy from 100-2000 K, 12 transitions of OH, and CO rotational lines ranging from J=14-13 up to J=40-39. [O I] is detected in all but one source in the entire sample; among the sources with detectable [O I] are two Very Low Luminosity Objects (VeLLOs). The mean 63/145 um [O I] flux ratio is 17.2 +/- 9.2. The [O I] 63 um line correlates with Lbol, but not with the time-averaged outflow rate derived from low-J CO maps. [C II] emission is in general not local to the source. The sample Lbol increased by 1.25 (1.06) and Tbol decreased to 0.96 (0.96) of mean (median) values with the inclusion of the Herschel data. Most CO rotational diagrams are characterized by two optically thin components (<N> = (0.70 +/- 1.12) x 10^49 total particles). N_CO correlates strongly with Lbol, but neither Trot nor N_CO(warm)/N_CO(hot) correlates with Lbol, suggesting that the total excited gas is related to the current source luminosity, but that the excitation is primarily determined by the physics of the interaction (e.g., UV- heating/shocks). Rotational temperatures for H2O (<Trot> = 194 +/- 85 K) and OH (<Trot> = 183 +/- 117 K) are generally lower than for CO, and much of the scatter in the observations about the best fit is attributed to differences in excitation conditions and optical depths amongst the detected lines.

preprint2013arXiv

SMA Observations of Class 0 Protostars: A High-Angular Resolution Survey of Protostellar Binary Systems

We present high angular resolution 1.3 mm and 850 um dust continuum data obtained with the Submillimeter Array toward 33 Class 0 protostars in nearby clouds (distance < 500 pc), which represents so far the largest survey toward protostellar binary/multiple systems. The median angular resolution in the survey is 2.5 arcsec, while the median linear resolution is approximately 600 AU. Compact dust continuum emission is observed from all sources in the sample. Twenty-one sources in the sample show signatures of binarity/multiplicity, with separations ranging from 50 to 5000 AU. The numbers of singles, binaries, triples, and quadruples in the sample are 12, 14, 5, and 2, respectively. The derived multiplicity frequency (MF) and companion star fraction (CSF) for Class 0 protostars are 0.64+/-0.08 and 0.91+/-0.05, respectively, with no correction for completeness. The derived MF and CSF in this survey are approximately two times higher than the values found in the binary surveys toward Class I YSOs, and approximately three (for MF) and four (for CSF) times larger than the values found among MS stars, with a similar range of separations. Furthermore, the observed fraction of high order multiple systems to binary systems in Class 0 protostars (0.50+/-0.09) is also larger than the fractions found in Class I YSOs (0.31+/-0.07) and MS stars (< 0.2). These results suggest that binary properties evolve as protostars evolve, as predicted by numerical simulations. The distribution of separations for Class 0 protostellar binary systems shows a general trend in which companion star fraction increases with decreasing companion separation. We find that 67%+/-8% of the protobinary systems have circumstellar mass ratios below 0.5, implying that unequal-mass systems are preferred in the process of binary star formation. We suggest an empirical sequential fragmentation picture for binary star formation.

preprint2013arXiv

The Luminosities of Protostars in the Spitzer c2d and Gould Belt Legacy Clouds

Motivated by the long-standing "luminosity problem" in low-mass star formation whereby protostars are underluminous compared to theoretical expectations, we identify 230 protostars in 18 molecular clouds observed by two Spitzer Space Telescope Legacy surveys of nearby star-forming regions. We compile complete spectral energy distributions, calculate Lbol for each source, and study the protostellar luminosity distribution. This distribution extends over three orders of magnitude, from 0.01 Lsun - 69 Lsun, and has a mean and median of 4.3 Lsun and 1.3 Lsun, respectively. The distributions are very similar for Class 0 and Class I sources except for an excess of low luminosity (Lbol < 0.5 Lsun) Class I sources compared to Class 0. 100 out of the 230 protostars (43%) lack any available data in the far-infrared and submillimeter (70 um < wavelength < 850 um) and have Lbol underestimated by factors of 2.5 on average, and up to factors of 8-10 in extreme cases. Correcting these underestimates for each source individually once additional data becomes available will likely increase both the mean and median of the sample by 35% - 40%. We discuss and compare our results to several recent theoretical studies of protostellar luminosities and show that our new results do not invalidate the conclusions of any of these studies. As these studies demonstrate that there is more than one plausible accretion scenario that can match observations, future attention is clearly needed. The better statistics provided by our increased dataset should aid such future work.

preprint2012arXiv

A Systematic Search for Molecular Outflows Toward Candidate Low-Luminosity Protostars and Very Low Luminosity Objects

We present a systematic single-dish search for molecular outflows toward a sample of 9 candidate low-luminosity protostars and 30 candidate Very Low Luminosity Objects (VeLLOs; L_int < 0.1 L_sun). The sources are identified using data from the Spitzer Space Telescope catalogued by Dunham et al. toward nearby (D < 400 pc) star forming regions. Each object was observed in 12CO and 13CO J = 2-1 simultaneously using the sideband separating ALMA Band-6 prototype receiver on the Heinrich Hertz Telescope at 30 arcsecond resolution. Using 5-point grid maps we identify five new potential outflow candidates and make on-the-fly maps of the regions surrounding sources in the dense cores B59, L1148, L1228, and L1165. Of these new outflow candidates, only the map of B59 shows a candidate blue outflow lobe associated with a source in our survey. We also present larger and more sensitive maps of the previously detected L673-7 and the L1251-A IRS4 outflows and analyze their properties in comparison to other outflows from VeLLOs. The accretion luminosities derived from the outflow properties of the VeLLOs with detected CO outflows are higher than the observed internal luminosity of the protostars, indicating that these sources likely had higher accretion rates in the past. The known L1251-A IRS3 outflow is detected but not remapped. We do not detect clear, unconfused signatures of red and blue molecular wings toward the other 31 sources in the survey indicating that large-scale, distinct outflows are rare toward this sample of candidate protostars. Several potential outflows are confused with kinematic structure in the surrounding core and cloud. Interferometric imaging is needed to disentangle large-scale molecular cloud kinematics from these potentially weak protostellar outflows.

preprint2012arXiv

CO2 Ice toward Low-luminosity, Embedded Protostars: Evidence for Episodic Mass Accretion via Chemical History

We present Spitzer IRS spectroscopy of CO2 ice bending mode spectra at 15.2 micrometer toward 19 young stellar objects with luminosity lower than 1 Lsun (3 with luminosity lower than 0.1 Lsun). Ice on dust grain surfaces can encode the history of heating because pure CO2 ice forms only at elevated temperature, T > 20 K, and thus around protostars of higher luminosity. Current internal luminosities of YSOs with L < 1 Lsun do not provide the conditions needed to produce pure CO2 ice at radii where typical envelopes begin. The presence of detectable amounts of pure CO2 ice would signify a higher past luminosity. Many of the spectra require a contribution from a pure, crystalline CO2 component, traced by the presence of a characteristic band splitting in the 15.2 micrometer bending mode. About half of the sources (9 out of 19) in the low luminosity sample have evidence for pure CO2 ice, and six of these have significant double-peaked features, which are very strong evidence of pure CO2 ice. The presence of the pure CO2 ice component indicates that the dust temperature, and hence luminosity of the central star/accretion disk system, must have been higher in the past. An episodic accretion scenario, in which mixed CO-CO2 ice is converted to pure CO2 ice during each high luminosity phase, explains the presence of pure CO2 ice, the total amount of CO2 ice, and the observed residual C18O gas.

preprint2012arXiv

Discovery of A Binary System in IRAM 04191+1522

We present high angular resolution observations of the Class 0 protostar IRAM04191+1522, using the Submillimeter Array (SMA). The SMA 1.3 mm continuum images reveal within IRAM04191+1522 two distinct sources with an angular separation of 7.8\,$\pm$\,0.2$"$. The two continuum sources are located in the southeast-northwest direction, with total gas masses of about 0.011 M_sun and about 0.005 M_sun, respectively. The southeastern source, associated with an infrared source seen in the Spitzer images, is the well-known Class 0 protostar with a bolometric luminosity of about 0.08 L_sun. The newly-discovered northwestern continuum source is not visible in the Spitzer images at wavelengths from 3.6 to 70 micron, and has an extremely low bolometric luminosity (< 0.03 L_sun). Complementary IRAM N2H+(1-0) data that probe the dense gas in the common envelope suggest that the two sources were formed through the rotational fragmentation of an elongated dense core. Furthermore, comparisons between IRAM04191+1522 and other protostars suggest that most cores with binary systems formed therein have ratios of rotational energy to gravitational energy $β_{\rm rot}$ > 1%. This is consistent with theoretical simulations and indicates that the level of rotational energy in a dense core plays an important role in the fragmentation process.

preprint2012arXiv

Revealing The Millimeter Environment of the New FU Orionis Candidate HBC722 with the Submillimeter Array

We present 230 GHz Submillimeter Array continuum and molecular line observations of the newly discovered FUor candidate HBC722. We report the detection of seven 1.3 mm continuum sources in the vicinity of HBC722, none of which correspond to HBC722 itself. We compile infrared and submillimeter continuum photometry of each source from previous studies and conclude that three are Class 0 embedded protostars, one is a Class I embedded protostar, one is a Class I/II transition object, and two are either starless cores or very young, very low luminosity protostars or first hydrostatic cores. We detect a northwest-southeast outflow, consistent with the previous detection of such an outflow in low-resolution, single-dish observations, and note that its axis may be precessing. We show that this outflow is centered on and driven by one of the nearby Class 0 sources rather than HBC722, and find no conclusive evidence that HBC722 itself is driving an outflow. The non-detection of HBC722 in the 1.3 mm continuum observations suggests an upper limit of 0.02 solar masses for the mass of the circumstellar disk. This limit is consistent with typical T Tauri disks and with a disk that provides sufficient mass to power the burst.

preprint2012arXiv

SMA and Spitzer Observations of Bok Glouble CB17: A Candidate First Hydrostatic Core?

We present high angular resolution SMA and Spitzer observations toward the Bok globule CB17. SMA 1.3mm dust continuum images reveal within CB17 two sources with an angular separation of about 21" (about 5250 AU at a distance of 250 pc). The northwestern continuum source, referred to as CB17 IRS, dominates the infrared emission in the Spitzer images, drives a bipolar outflow extending in the northwest-southeast direction, and is classified as a low luminosity Class0/I transition object (L_bol ~ 0.5 L_sun). The southeastern continuum source, referred to as CB17 MMS, has faint dust continuum emission in the SMA 1.3mm observations (about 6 sigma detection; ~3.8 mJy), but is not detected in the deep Spitzer infrared images at wavelengths from 3.6 to 70 micron. Its bolometric luminosity and temperature, estimated from its spectral energy distribution, are less than 0.04 L_sun and 16 K, respectively. The SMA CO(2-1) observations suggest that CB17 MMS may drive a low-velocity molecular outflow (about 2.5 km/s), extending in the east-west direction. Comparisons with prestellar cores and Class0 protostars suggest that CB17 MMS is more evolved than prestellar cores but less evolved than Class0 protostars. The observed characteristics of CB17 MMS are consistent with the theoretical predictions from radiative/magneto hydrodynamical simulations of a first hydrostatic core, but there is also the possibility that CB17 MMS is an extremely low luminosity protostar deeply embedded in an edge-on circumstellar disk. Further observations are needed to study the properties of CB17 MMS and to address more precisely its evolutionary stage.

preprint2011arXiv

Detection of a Bipolar Molecular Outflow Driven by a Candidate First Hydrostatic Core

We present new 230 GHz Submillimeter Array observations of the candidate first hydrostatic core Per-Bolo 58. We report the detection of a 1.3 mm continuum source and a bipolar molecular outflow, both centered on the position of the candidate first hydrostatic core. The continuum detection has a total flux density of 26.6 +/- 4.0 mJy, from which we calculate a total (gas and dust) mass of 0.11 +/- 0.05 Msun and a mean number density of 2.0 +/- 1.6 X 10^7 cm-3. There is some evidence for the existence of an unresolved component in the continuum detection, but longer-baseline observations are required in order to confirm the presence of this component and determine whether its origin lies in a circumstellar disk or in the dense inner envelope. The bipolar molecular outflow is observed along a nearly due east-west axis. The outflow is slow (characteristic velocity of 2.9 km/s), shows a jet-like morphology (opening semi-angles ~8 degrees for both lobes), and extends to the edges of the primary beam. We calculate the kinematic and dynamic properties of the outflow in the standard manner and compare them to several other protostars and candidate first hydrostatic cores with similarly low luminosities. We discuss the evidence both in support of and against the possibility that Per-Bolo 58 is a first hydrostatic core, and we outline future work needed to further evaluate the evolutionary status of this object.

preprint2011arXiv

Disentangling the Environment of the FU Orionis Candidate HBC 722 with Herschel

We analyze the submillimeter emission surrounding the new FU Orionis-type object, HBC 722. We present the first epoch of observations of the active environs of HBC 722, with imaging and spectroscopy from PACS, SPIRE, and HIFI aboard the Herschel Space Observatory, as well as CO J= 2-1 and 350 um imaging (SHARC-II) with the Caltech Submillimeter Observatory. The primary source of submillimeter continuum emission in the region -- 2MASS 20581767+4353310 -- is located 16$\arcsec$ south-southeast of the optical flaring source while the optical and near-IR emission is dominated by HBC 722. A bipolar outflow extends over HBC 722; the most likely driver is the submillimeter source. We detect warm (100 K) and hot (246 K) CO emission in the surrounding region, evidence of outflow-driven heating in the vicinity. The region around HBC 722 itself shows little evidence of heating driven by the new outbursting source itself.

preprint2011arXiv

Resolving the Luminosity Problem in Low-Mass Star Formation

We determine the observational signatures of protostellar cores by coupling two-dimensional radiative transfer calculations with numerical hydrodynamical simulations that predict accretion rates that both decline with time and feature short-term variability and episodic bursts caused by disk gravitational instability and fragmentation. We calculate the radiative transfer of the collapsing cores throughout the full duration of the collapse, using as inputs the core, disk, and protostellar masses, radii, and mass accretion rates predicted by the hydrodynamical simulations. From the resulting spectral energy distributions, we calculate standard observational signatures (bolometric luminosity, bolometric temperature, ratio of bolometric to submillimeter luminosity) to directly compare to observations. We show that the accretion process predicted by these models reproduces the full spread of observed protostars in both Lbol - Tbol and Lbol - core mass space, including very low luminosity objects, provides a reasonable match to the observed protostellar luminosity distribution, and resolves the long-standing luminosity problem. These models predict an embedded phase duration shorter than recent observationally determined estimates (0.12 Myr vs. 0.44 Myr), and a fraction of total time spent in Stage 0 of 23%, consistent with the range of values determined by observations. On average, the models spend 1.3% of their total time in accretion bursts, during which 5.3% of the final stellar mass accretes, with maximum values being 11.8% and 35.5% for the total time and accreted stellar mass, respectively. Time-averaged models that filter out the accretion variability and bursts do not provide as good of a match to the observed luminosity problem, suggesting that the bursts are required.

preprint2011arXiv

The Spitzer c2d Survey of Nearby Dense Cores XI : Infrared and Submillimeter Observations of CB130

We present new observations of the CB130 region, composed of three separate cores. Using the \textit{Spitzer Space Telescope} we detected a Class 0 and a Class II object in one of these, CB130-1. The observed photometric data from \textit{Spitzer} and ground-based telescopes are used to establish the physical parameters of the Class 0 object. SED fitting with a radiative transfer model shows that the luminosity of the Class 0 object is 0.14 $-$ 0.16 L$_{\odot}$, which is a low luminosity for a protostellar object. In order to constrain the chemical characteristics of the core having the low luminosity object, we compare our molecular line observations to models of lines including abundance variations. We tested both ad hoc step function abundance models and a series of self-consistent chemical evolution models. In the chemical evolution models, we consider a continuous accretion model and an episodic accretion model to explore how variable luminosity affects the chemistry. The step function abundance models can match observed lines reasonably well. The best fitting chemical evolution model requires episodic accretion and the formation of CO$_2$ ice from CO ice during the low luminosity periods. This process removes C from the gas phase, providing a much improved fit to the observed gas-phase molecular lines and the CO$_2$ ice absorption feature. Based on the chemical model result, the low luminosity of CB130-1 is explained better as a quiescent stage between episodic accretion bursts rather than being at the first hydrostatic core stage.

preprint2011arXiv

The Spitzer Survey of Interstellar Clouds in the Gould Belt. III. A Multi-Wavelength View of Corona Australis

We present Spitzer Space Telescope IRAC and MIPS observations of a 0.85 deg^2 field including the Corona Australis (CrA) star-forming region. At a distance of 130 pc, CrA is one of the closest regions known to be actively forming stars, particularly within its embedded association, the Coronet. Using the Spitzer data, we identify 51 young stellar objects (YSOs) in CrA which include sources in the well-studied Coronet cluster as well as distributed throughout the molecular cloud. Twelve of the YSOs discussed are new candidates, one of which is located in the Coronet. Known YSOs retrieved from the literature are also added to the list, and a total of 116 candidate YSOs in CrA are compiled. Based on these YSO candidates, the star formation rate is computed to be 12 M_o Myr^-1, similar to that of the Lupus clouds. A clustering analysis was also performed, finding that the main cluster core, consisting of 68 members, is elongated (having an aspect ratio of 2.36), with a circular radius of 0.59 pc and mean surface density of 150 pc^-2. In addition, we analyze outflows and jets in CrA by means of new CO and H_2 data. We present 1.3 mm interferometric continuum observations made with the Submillimeter Array (SMA) covering R CrA, IRS 5, IRS 7, and IRAS 18595-3712 (IRAS 32). We also present multi-epoch H_2 maps and detect jets and outflows, study their proper motions, and identify exciting sources. The Spitzer and ISAAC/VLT observations of IRAS 32 show a bipolar precessing jet, which drives a CO (2-1) outflow detected in the SMA observations. There is also clear evidence for a parsec-scale precessing outflow, E-W oriented, and originating in the SMA 2 region, likely driven by SMA 2 or IRS 7A.

preprint2011arXiv

The Spitzer Survey of Interstellar Clouds in the Gould Belt. IV. Lupus V and VI Observed with IRAC and MIPS

We present Gould's Belt (GB) Spitzer IRAC and MIPS observations of the Lupus V and VI clouds and discuss them in combination with near-infrared (2MASS) data. Our observations complement those obtained for other Lupus clouds within the frame of the Spitzer "Core to Disk" (c2d) Legacy Survey. We found 43 Young Stellar Object (YSO) candidates in Lupus V and 45 in Lupus VI, including 2 transition disks, using the standard c2d/GB selection method. None of these sources was classified as a pre-main sequence star from previous optical, near-IR and X-ray surveys. A large majority of these YSO candidates appear to be surrounded by thin disks (Class III; ~79% in Lupus V and ~87% in Lupus VI). These Class III abundances differ significantly from those observed for the other Lupus clouds and c2d/GB surveyed star-forming regions, where objects with optically thick disks (Class II) dominate the young population. We investigate various scenarios that can explain this discrepancy. In particular, we show that disk photo-evaporation due to nearby OB stars is not responsible for the high fraction of Class III objects. The gas surface densities measured for Lupus V and VI lies below the star-formation threshold (AV {\eqsim}8.6 mag), while this is not the case for other Lupus clouds. Thus, few Myrs older age for the YSOs in Lupus V and VI with respect to other Lupus clouds is the most likely explanation of the high fraction of Class III objects in these clouds, while a higher characteristic stellar mass might be a contributing factor. Better constraints on the age and binary fraction of

preprint2010arXiv

A Candidate Detection of the First Hydrostatic Core

The first hydrostatic core (FHSC) represents a very early phase in the low-mass star formation process, after collapse of the parent core has begun but before a true protostar has formed. This large (few AU), cool (100 K), pressure supported core of molecular hydrogen is expected from theory, but has yet to be observationally verified. Here we present observations of an excellent candidate for the FHSC phase: Per-Bolo 58, a dense core in Perseus that was previously believed to be starless. The 70 micron flux of 65 mJy, from new deep Spitzer MIPS observations, is consistent with that expected for the FHSC. A low signal-to-noise detection at 24 micron leaves open the possibility that Per-Bolo 58 could be a very low luminosity protostar, however. We utilize radiative transfer models to determine the best-fitting FHSC and protostar models to the spectral energy distribution and 2.9 mm visibilities of Per-Bolo 58. The source is consistent with a FHSC with some source of lower opacity through the envelope allowing 24 micron emission to escape; a small outflow cavity or a cavity in the envelope are both possible. While we are unable to rule out the presence of a protostar, if present it would be one of the lowest luminosity protostellar objects yet observed, with an internal luminosity of approximately 0.01 Lsun.

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

Disk and Envelope Structure in Class 0 Protostars: I. The Resolved Massive Disk in Serpens FIRS 1

We present the first results of a program to characterize the disk and envelope structure of typical Class 0 protostars in nearby low-mass star forming regions. We use Spitzer IRS mid-infrared spectra, high resolution CARMA 230 GHz continuum imaging, and 2-D radiative transfer models to constrain the envelope structure, as well as the size and mass of the circum-protostellar disk in Serpens FIRS 1. The primary envelope parameters (centrifugal radius, outer radius, outflow opening angle, and inclination) are well constrained by the spectral energy distribution (SED), including Spitzer IRAC and MIPS photometry, IRS spectra, and 1.1 mm Bolocam photometry. These together with the excellent uv-coverage (4.5-500 klam) of multiple antenna configurations with CARMA allow for a robust separation of the envelope and a resolved disk. The SED of Serpens FIRS 1 is best fit by an envelope with the density profile of a rotating, collapsing spheroid with an inner (centrifugal) radius of approximately 600 AU, and the millimeter data by a large resolved disk with Mdisk~1.0 Msun and Rdisk~300 AU. These results suggest that large, massive disks can be present early in the main accretion phase. Results for the larger, unbiased sample of Class~0 sources in the Perseus, Serpens, and Ophiuchus molecular clouds are needed to determine if relatively massive disks are typical in the Class 0 stage.

preprint2009arXiv

Evolutionary Signatures in the Formation of Low-Mass Protostars. II. Towards Reconciling Models and Observations

A long-standing problem in low-mass star formation is the "luminosity problem," whereby protostars are underluminous compared to the accretion luminosity expected both from theoretical collapse calculations and arguments based on the minimum accretion rate necessary to form a star within the embedded phase duration. Motivated by this luminosity problem, we present a set of evolutionary models describing the collapse of low-mass, dense cores into protostars, using the Young & Evans (2005) model as our starting point. We calculate the radiative transfer of the collapsing cores throughout the full duration of the collapse in two dimensions. From the resulting spectral energy distributions, we calculate standard observational signatures to directly compare to observations. We incorporate several modifications and additions to the original Young & Evans model in an effort to better match observations with model predictions. We find that scattering, 2-D geometry, mass-loss, and outflow cavities all affect the model predictions, as expected, but none resolve the luminosity problem. A cycle of episodic mass accretion, however, can resolve this problem and bring the model predictions into better agreement with observations. Standard assumptions about the interplay between mass accretion and mass loss in our model give star formation efficiencies consistent with recent observations that compare the core mass function (CMF) and stellar initial mass function (IMF). The combination of outflow cavities and episodic mass accretion reduce the connection between observational Class and physical Stage to the point where neither of the two common observational signatures (bolometric temperature and ratio of bolometric to submillimeter luminosity) can be considered reliable indicators of physical Stage.

preprint2009arXiv

Serpens Cluster B and VV Ser Observed With High Spatial Resolution at 70, 160, and 350um

We report on diffraction-limited observations in the far-infrared and sub- millimeter of the Cluster B region of Serpens (G3-G6 Cluster) and of the Herbig Be star to the south, VV Ser. The observations were made with the Spitzer MIPS instrument in fine-scale mode at 70um, in normal mapping mode at 160um (VV Ser only), and the CSO SHARC-II camera at 350um (Cluster B only). We use these data to define the spectral energy distributions of the tightly grouped members of Cluster B, many of whose SEDs peak in the far-infrared. We compare our results to those of the c2d survey of Serpens and to published models for the far-infrared emission from VV Ser. We find that values of Lbol and Tbol calculated with our new photometry show only modest changes from previous values, and that most source SED classifications remain unchanged.

preprint2009arXiv

Star formation in Perseus - V. Outflows detected by HARP

Molecular outflows provide an alternative method of identifying protostellar cores, complementary to recent mid-infrared studies. Continuing our studies of Perseus, we investigate whether all Spitzer-identified protostars, and particularly those with low luminosities, drive outflows, and if any new protostellar cores (perhaps harbouring low-mass sources) can be identified via their outflows alone. We have used the heterodyne array receiver HARP on JCMT to make deep 12CO 3-2 maps of submm cores in Perseus, extending and deepening our earlier study with RxB and bringing the total number of SCUBA cores studied up to 83. Our survey includes 23/25 of the Dunham et al. (2008) Spitzer low-luminosity objects believed to be embedded protostars, including three VeLLOs. All but one of the cores identified as harbouring embedded YSOs have outflows, confirming outflow detections as a good method for identifying protostars. We detect outflows from 20 Spitzer low-luminosity objects. We do not conclusively detect any outflows from IR-quiet cores, though confusion in clustered regions such as NGC1333 makes it impossible to identify all the individual driving sources. This similarity in detection rates despite the difference in search methods and detection limits suggests either that the sample of protostars in Perseus is now complete, or that the existence of an outflow contributes to the Spitzer detectability, perhaps through the contribution of shocked H2 emission in the IRAC bands. For five of the low-luminosity sources, there is no protostellar envelope detected at 350 microns and the Spitzer emission is entirely due to shocks. Additionally, we detect the outflow from IRAS 03282+3035 at 850 microns with SCUBA due to CO line contamination in the continuum passband.

preprint2009arXiv

The Spitzer c2d Survey of Nearby Dense Cores: Jet and Molecular Outflow Associated with a YSO in core A of L1251

A long infrared jet has been discovered by the Spitzer c2d legacy program in core A of L1251. It is associated with a very embedded Class 0 object with an accretion luminosity of about 0.9 Lsun derived by radiative transfer model fitting to the observed SED. Comparing the observed IRAC colors along the infrared jet with those calculated from a model of an admixture of gas with a power-law temperature distribution indicates that the jet is possibly created by a paraboloidal bow shock propagating into the ambient medium of n(H_2)=10^5 cm^{-3}. In addition, the variation of the power-law index along the jet suggests that the portion of hot gas decreases with distance from the jet engine. The molecular outflow in this region has been mapped for the first time using CO data. From the calculated outflow momentum flux, a very strong lower limit to the average accretion luminosity is 3.6 sin i/cos^3 i Lsun, indicative of a decrease in the accretion rate with time.

preprint2008arXiv

The Spitzer c2d Legacy Results: Star Formation Rates and Efficiencies; Evolution and Lifetimes

(Abridged) The c2d Spitzer Legacy project obtained images and photometry with both IRAC and MIPS instruments for five large, nearby molecular clouds. This paper combines information drawn from studies of individual clouds into a combined and updated statistical analysis of star formation rates and efficiencies, numbers and lifetimes for SED classes, and clustering properties. Current star formation efficiencies range from 3% to 6%. Taken together, the five clouds are producing about 260 solar masses of stars per Myr. The star formation surface density is more than an order of magnitude larger than would be predicted from the Kennicutt relation used in extragalactic studies. Measured against the dense gas probed by the maps of dust continuum emission, the efficiencies are much higher, and the current stock of dense cores would be exhausted in 1.8 Myr on average. The derived lifetime for the Class I phase is 0.44 to 0.54 Myr, considerably longer than some estimates. Similarly, the lifetime for the Class 0 SED class, 0.10 to 0.16 Myr, is longer than early estimates. The great majority (90%) of young stars lie within loose clusters with at least 35 members and a stellar density of 1 solar mass per cubic pc. Accretion at the sound speed from an isothermal sphere over the lifetime derived for the Class I phase could build a star of about 0.25 solar masses, given an efficiency of 0.3. Our data confirm and aggravate the "luminosity problem" for protostars. Our results strongly suggest that accretion is time variable, with prolonged periods of very low accretion. Based on a very simple model and this sample of sources, half the mass of a star would be accreted during only 7% of the Class I lifetime, as represented by the eight most luminous objects.

preprint2006arXiv

SHARC-II Mapping of Spitzer c2d Small Clouds and Cores

We present the results of a submillimeter survey of 53 low-mass dense cores with the Submillimeter High Angular Resolution Camera II (SHARC-II). The survey is a follow-up project to the Spitzer Legacy Program ``From Molecular Cores to Planet-Forming Disks'', with the purpose being to create a complete data set of nearby low-mass dense cores from the infrared to the millimeter. We present maps of 52 cores at 350 microns and three cores at 450 microns, two of which were observed at both wavelengths. Of these 52 cores, 41 were detected by SHARC-II: 32 contained one submillimeter source while 9 contained multiple sources. For each submillimeter source detected, we report various source properties including source position, fluxes in various apertures, size, aspect ratio, and position angle. For the 12 cores that were not detected we present upper limits. The sources detected by SHARC-II have, on average, smaller sizes at the 2sigma contours than those derived from longer-wavelength bolometer observations. We conclude that this is not caused by a failure to integrate long enough to detect the full extent of the core; instead it arises primarily from the fact that the observations presented in this survey are insensitive to smoothly varying extended emission. We find that SHARC-II observations of low-mass cores are much better suited to distinguishing between starless and protostellar cores than observations at longer wavelengths. Very Low Luminosity Objects, a new class of objects being discovered by the Spitzer Space Telescope in cores previously classified as starless, look very similar at 350 microns to other cores with more luminous protostars.

preprint2006arXiv

The Spitzer c2d survey of Large, Nearby, Interstellar Clouds. V. Chamaeleon II Observed with IRAC

We present IRAC (3.6, 4.5, 5.8, and 8.0 micron) observations of the Chamaeleon II molecular cloud. The observed area covers about 1 square degree defined by $A_V >2$. Analysis of the data in the 2005 c2d catalogs reveals a small number of sources (40) with properties similar to those of young stellaror substellar objects (YSOs). The surface density of these YSO candidates is low, and contamination by background galaxies appears to be substantial, especially for sources classified as Class I or flat SED. We discuss this problem in some detail and conclude that very few of the candidate YSOs in early evolutionary stages are actually in the Cha II cloud. Using a refined set of criteria, we define a smaller, but more reliable, set of 24 YSO candidates.

preprint2006arXiv

The Spitzer c2d Survey of Nearby Dense Cores: II: Discovery of a Low Luminosity Object in the "Evolved Starless Core" L1521F

We present Spitzer Space Telescope observations of the "evolved starless core" L1521F which reveal the presence of a very low luminosity object (L < 0.07 Lsun). The object, L1521F-IRS, is directly detected at mid-infrared wavelengths (>5 micron) but only in scattered light at shorter infrared wavelengths, showing a bipolar nebula oriented east-west which is probably tracing an outflow cavity. The nebula strongly suggests that L1521F-IRS is embedded in the L1521F core. Thus L1521F-IRS is similar to the recently discovered L1014-IRS and the previously known IRAM 04191 in its substellar luminosity and dense core environment. However these objects differ significantly in their core density, core chemistry, and outflow properties, and some may be destined to be brown dwarfs rather than stars.