Michael M. Dunham

Identifying the Low-Luminosity Population of Embedded Protostars in the c2d Observations of Clouds and Cores

We present the results of a search for all embedded protostars with internal luminosities ≤1.0 L☉ in the full sample of nearby, low-mass star-forming regions surveyed by the Spitzer Space Telescope Legacy Project From Molecular Cores to Planet Forming Disks (c2d). The internal luminosity of a source, Lint, is the luminosity of the central source and excludes luminosity arising from external heating. On average, the Spitzer c2d data are sensitive to embedded protostars with -->Lint ≥ 4 × 10−3(d/140 pc)2 L☉, a factor of 25 better than the sensitivity of the Infrared Astronomical Satellite (IRAS) to such objects. We present a set of selection criteria used to identify candidates from the Spitzer data and examine complementary data to decide whether each candidate is truly an embedded protostar. We find a tight correlation between the 70 μm flux and internal luminosity of a protostar, an empirical result based on both observations and detailed two-dimensional radiative transfer models of protostars. We identify 50 embedded protostars with -->Lint ≤ 1.0 L☉; 15 have -->Lint ≤ 0.1 L☉. The intrinsic distribution of source luminosities increases to lower luminosities. While we find sources down to the above sensitivity limit, indicating that the distribution may extend to luminosities lower than probed by these observations, we are able to rule out a continued rise in the distribution below -->Lint = 0.1 L☉. Between 75% and 85% of cores classified as starless prior to being observed by Spitzer remain starless to our luminosity sensitivity; the remaining 15%-25% harbor low-luminosity, embedded protostars. We compile complete spectral energy distributions for all 50 objects and calculate standard evolutionary signatures (Lbol, Tbol, and Lbol/Lsmm) and argue that these objects are inconsistent with the simplest picture of star formation, wherein mass accretes from the core onto the protostar at a constant rate.

Evolutionary Signatures in the Formation of Low-Mass Protostars. II. Toward 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. We use as our starting point the evolutionary model following the inside-out collapse of a singular isothermal sphere as presented by Young & Evans. We calculate the radiative transfer of the collapsing core throughout the full duration of the collapse in two dimensions. From the resulting spectral energy distributions, we calculate standard observational signatures (L bol, T bol, L bol/L smm) 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 include (1) the opacity from scattering in the radiative transfer, (2) a circumstellar disk directly in the two-dimensional radiative transfer, (3) a two-dimensional envelope structure, taking into account the effects of rotation, (4) mass-loss and the opening of outflow cavities, and (5) a simple treatment of episodic mass accretion. We find that scattering, two-dimensional geometry, mass-loss, and outflow cavities all affect the model predictions, as expected, but none resolve the luminosity problem. On the other hand, we find that a cycle of episodic mass accretion similar to that predicted by recent theoretical work 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 and stellar initial mass function. Finally, the combination of outflow cavities and episodic mass accretion reduces the connection between observational class and physical stage to the point where neither of the two commonly used observational signatures (T bol and L bol/L smm) can be considered reliable indicators of physical stage.

The Spitzer c2d Survey of Nearby Dense Cores. I. First Direct Detection of the Embedded Source in IRAM 04191+1522

We report the first detections of the Class 0 protostellar source IRAM 04191+1522 at wavelengths shortward of 60 μm with the Spitzer Space Telescope. We see extended emission in the Spitzer images that suggests the presence of an outflow cavity in the circumstellar envelope. We combine the Spitzer observations with existing data to form a complete data set ranging from 3.6 to 1300 μm and use these data to construct radiative transfer models of the source. We conclude that the internal luminosity of IRAM 04191+1522, defined to be the sum of the luminosity from the internal sources (a star and a disk), is L_(int) = 0.08 ± 0.04 L_⊙, placing it among the lowest luminosity protostars known. Although it was discovered before the launch of the Spitzer Space Telescope, IRAM 04191+1522 falls within a new class of very low luminosity objects being discovered by Spitzer. Unlike the two other well-studied objects in this class, which are associated either with weak, compact outflows or no outflows at all, IRAM 04191+1522 has a well-defined molecular outflow with properties consistent with those expected based on relations derived from higher luminosity (L_(int) ≥ 1 L_⊙) protostars. We discuss the difficulties in understanding IRAM 04191+1522 in the context of the standard model of star formation and suggest a possible explanation for the very low luminosity of this source.

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 μm 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 μm 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 an FHSC with some source of lower opacity through the envelope allowing 24 μm emission to escape; a small outflow cavity and 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 ~0.01 L ☉.

First results of the Herschel key program "Dust, Ice and Gas In Time" (DIGIT): Dust and gas spectroscopy of HD 100546

Context. We present far-infrared spectroscopic observations, taken with the Photodetector Array Camera and Spectrometer (PACS) on the Herschel Space Observatory, of the protoplanetary disk around the pre-main-sequence star HD100546. These observations are the first within the DIGIT Herschel key program, which aims to follow the evolution of dust, ice, and gas from young stellar objects still embedded in their parental molecular cloud core, through the final pre-main-sequence phases when the circumstellar disks are dissipated. Aims. Our aim is to improve the constraints on temperature and chemical composition of the crystalline olivines in the disk of HD100546 and to give an inventory of the gas lines present in its far-infrared spectrum. Methods. The 69 μm feature is analyzed in terms of position and shape to derive the dust temperature and composition. Furthermore, we detected 32 emission lines from five gaseous species and measured their line fluxes. Results. The 69 μm emission comes either from dust grains with ~70 K at radii larger than 50 AU, as suggested by blackbody fitting, or it arises from ~200K dust at ~13 AU, close to the midplane, as supported by radiative transfer models. We also conclude that the forsterite crystals have few defects and contain at most a few percent iron by mass. Forbidden line emission from [C_(II)] at 157 μm and [O_I] at 63 and 145 μm, most likely due to photodissociation by stellar photons, is detected. Furthermore, five H_2O and several OH lines are detected. We also found high-J rotational transition lines of CO, with rotational temperatures of ~300K for the transitions up to J = 22−21 and T ~ 800 K for higher transitions.

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 Infrared Spectrograph (IRS) mid-infrared spectra, high-resolution Combined Array for Research in Millimeter-wave Astronomy (CARMA) 230 GHz continuum imaging, and two-dimensional radiative transfer models to constrain the envelope structure, as well as the size and mass of the circumprotostellar 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 kλ) 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 M ☉ 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.

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 Mand a mean number density of 2.0 ± 1.6 × 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 −1 ), shows a jet-like morphology (opening semi-angles ∼8 ◦ 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.

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 12CO (2-1) and 17 are mapped in 12CO (3-2), with 6 mapped in both transitions. For each outflow, we calculate and tabulate the mass (M flow), momentum (P flow), kinetic energy (E flow), mechanical luminosity (L flow), and force (F flow) assuming optically thin emission in LTE at an excitation temperature, T ex, 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 12CO (3-2) maps compared to the 12CO (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 12CO (3-2) line is subthermally excited in outflows, and our results support this finding.

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 (refs 3, 10), 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 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.

A Triple Protostar System Formed via Fragmentation of a Gravitationally Unstable Disk

Binary and multiple star systems are a frequent outcome of the star formation process and as a result almost half of all stars with masses similar to that of the Sun have at least one companion star. Theoretical studies indicate that there are two main pathways that can operate concurrently to form binary/multiple star systems: large-scale fragmentation of turbulent gas cores and filaments or smaller-scale fragmentation of a massive protostellar disk due to gravitational instability. Observational evidence for turbulent fragmentation on scales of more than 1,000 astronomical units has recently emerged. Previous evidence for disk fragmentation was limited to inferences based on the separations of more-evolved pre-main sequence and protostellar multiple systems. The triple protostar system L1448 IRS3B is an ideal system with which to search for evidence of disk fragmentation as it is in an early phase of the star formation process, it is likely to be less than 150,000 years old and all of the protostars in the system are separated by less than 200 astronomical units. Here we report observations of dust and molecular gas emission that reveal a disk with a spiral structure surrounding the three protostars. Two protostars near the centre of the disk are separated by 61 astronomical units and a tertiary protostar is coincident with a spiral arm in the outer disk at a separation of 183 astronomical units. The inferred mass of the central pair of protostellar objects is approximately one solar mass, while the disk surrounding the three protostars has a total mass of around 0.30 solar masses. The tertiary protostar itself has a minimum mass of about 0.085 solar masses. We demonstrate that the disk around L1448 IRS3B appears susceptible to disk fragmentation at radii between 150 and 320 astronomical units, overlapping with the location of the tertiary protostar. This is consistent with models for a protostellar disk that has recently undergone gravitational instability, spawning one or two companion stars.