Observational signatures of outside-in accretion bursts in embedded protostars

We present key results from the Herschel Orion Protostar Survey: spectral energy distributions (SEDs) and model fits of 330 young stellar objects, predominantly protostars, in the Orion molecular clouds. This is the largest sample of protostars studied in a single, nearby star formation complex. With near-infrared photometry from 2MASS, mid- and far-infrared data from Spitzer and Herschel, and submillimeter photometry from APEX, our SEDs cover 1.2–870 μm and sample the peak of the protostellar envelope emission at ∼100 μm. Using mid-IR spectral indices and bolometric temperatures, we classify our sample into 92 Class 0 protostars, 125 Class I protostars, 102 flat-spectrum sources, and 11 Class II pre-main-sequence stars. We implement a simple protostellar model (including a disk in an infalling envelope with outflow cavities) to generate a grid of 30,400 model SEDs and use it to determine the best-fit model parameters for each protostar. We argue that far-IR data are essential for accurate constraints on protostellar envelope properties. We find that most protostars, and in particular the flat-spectrum sources, are well fit. The median envelope density and median inclination angle decrease from Class 0 to Class I to flat-spectrum protostars, despite the broad range in best-fit parameters in each of the three categories. We also discuss degeneracies in our model parameters. Our results confirm that the different protostellar classes generally correspond to an evolutionary sequence with a decreasing envelope infall rate, but the inclination angle also plays a role in the appearance, and thus interpretation, of the SEDs.

We present new Keck images at 0.9 μm and OVRO 1.3 mm continuum images of five Class I protostars in the Taurus star-forming region. We analyze these data in conjunction with broadband spectral energy distributions (SEDs) and 8-13 μm spectra from the literature using a Monte Carlo radiative transfer code. By fitting models for the circumstellar dust distributions simultaneously to the scattered light images, millimeter continuum data, and the SEDs, we attempt to distinguish between flared disks, infalling envelopes with outflow cavities, and combinations of disks and envelopes. For each of these circumstellar density distributions, we generate grids of models for varying geometries, dust masses, and accretion rates and determine the best fits by minimizing the residuals between model and data. Comparison of the residuals for best-fit disk, envelope, and disk+envelope models demonstrates that, in general, models incorporating both massive envelopes and massive embedded disks fit the imaging+SED data best. The implied envelope infall rates for these disk+envelope models are generally consistent with infall rates derived by previous investigators, although they are approximately an order of magnitude larger than inner disk accretion rates inferred from recent spectroscopic measurements. In addition, the disk masses inferred from our models are close to or larger than the limit for gravitationally stable disks, indicating that Class I disks may undergo periodic episodes of enhanced accretion, perhaps as a result of gravitational instabilities. An important caveat to these results is that in some cases, no single model can fit all of the imaging and SED data well, suggesting that further refinements to models of the circumstellar dust distributions around Class I sources are necessary. We discuss several potential improvements to the models, as well as new constraints that will become available with upcoming millimeter and infrared facilities.

We present Spitzer IRAC images that indicate the presence of cavities cut into the dense outer envelope surrounding very young pre-main-sequence stars. These young stellar objects (YSOs) characterized by an outflow represent the earliest stages of star formation. Mid-infrared photons thermally created by the central protostar/disk are scattered by dust particles within the outflow cavity itself into the line of sight. We observed this scattered light from 27 nearby, cavity-resolved YSOs and quantified the shape of the outflow cavities. Using the grid models of Robitaille et al., we matched model spectral energy distributions (SEDs) to the observed SEDs of the 27 cataloged YSOs using photometry from IRAC, MIPS, and IRAS. This allows for the estimation of geometric and physical properties such as inclination angle, cavity density, and accretion rate. By using the relative parameter estimates determined by the models, we are able to deduce an evolutionary picture for outflows. Our work supports the concept that cavities widen with time, beginning as a thin jetlike outflow that widens to reveal the central protostar and disk until the protostellar envelope is completely dispersed by outflow and accretion.

Massive stars form whilst they are still embedded in dense envelopes. As a result, the roles of rotation, mass loss and accretion in massive star formation are not well understood. This study evaluates the source of the Q-band, lambda=19.5 microns, emission of massive young stellar objects (MYSOs). This allows us to determine the relative importance of rotation and outflow activity in shaping the circumstellar environments of MYSOs on 1000 AU scales. We obtained diffraction limited mid-infrared images of a sample of 20 MYSOs using the VLT/VISIR and Subaru/COMICS instruments. For these 8 m class telescopes and the sample selected, the diffraction limit, ~0.6", corresponds to approximately 1000 AU. We compare the images and the spectral energy distributions (SEDs) observed to a 2D, axis-symmetric dust radiative transfer model that reproduces VLTI/MIDI observations of the MYSO W33A. We vary the inclination, mass infall rate, and outflow opening angle to simultaneously recreate the behaviour of the sample of MYSOs in the spatial and spectral domains. The mid-IR emission of 70 percent of the MYSOs is spatially resolved. In the majority of cases, the spatial extent of their emission and their SEDs can be reproduced by the W33A model featuring an in-falling, rotating dusty envelope with outflow cavities. There is independent evidence that most of the sources which are not fit by the model are associated with ultracompact HII regions and are thus more evolved. We find that, in general, the diverse 20 micron morphology of MYSOs can be attributed to warm dust in the walls of outflow cavities seen at different inclinations. This implies that the warm dust in the outflow cavity walls dominates the Q-band emission of MYSOs. In turn, this emphasises that outflows are an ubiquitous feature of massive star formation.

We use Herschel spectrophotometry of BHR71, an embedded Class 0 protostar, to provide new constraints on its physical properties. We detect 645 (non-unique) spectral lines among all spatial pixels. At least 61 different spectral lines originate from the central region. A CO rotational diagram analysis shows four excitation temperature components, 43, 197, 397, and 1057 K. Low-J CO lines trace the outflow while the high-J CO lines are centered on the infrared source. The low-excitation emission lines of trace the large-scale outflow, while the high-excitation emission lines trace a small-scale distribution around the equatorial plane. We model the envelope structure using the dust radiative transfer code, hyperion, incorporating rotational collapse, an outer static envelope, outflow cavity, and disk. The evolution of a rotating collapsing envelope can be constrained by the far-infrared/millimeter spectral energy distribution along with the azimuthally averaged radial intensity profile, and the structure of the outflow cavity plays a critical role at shorter wavelengths. Emission at 20–40 μm requires a cavity with a constant-density inner region and a power-law density outer region. The best-fit model has an envelope mass of 19 inside a radius of 0.315 pc and a central luminosity of 18.8 . The time since collapse began is 24,630–44,000 years, most likely around 36,000 years. The corresponding mass infall rate in the envelope (1.2 × 10−5 ) is comparable to the stellar mass accretion rate, while the mass-loss rate estimated from the CO outflow is 20% of the stellar mass accretion rate. We find no evidence for episodic accretion.

Aims. In this paper we aim to determine the structure on 100 AU scales of the massive young stellar object W33A, using interferometric observations in the mid-infrared. This emission could be caused by a variety of elements, for example, the inner protostellar envelope, outflow cavity walls, or a dusty or gaseous accretion disk. Methods. We used the Unit Telescopes of the VLT Interferometer in conjunction with the MIDI instrument to obtain spectrally dispersed visibilities in the N-band on 4 baselines with an angular resolution between 25 and 60 milli-arcsec (equivalent to 95 and 228 AU at 3.8 kpc). The visibility spectra and spectral energy distribution were compared to 2D-axi-symmetric dust radiative transfer models with a geometry that includes a rotationally flattened envelope and outflow cavities. We assumed an O7.5 ZAMS star as the central source, consistent with the observed bolometric luminosity. The observations were compared to models with and without (dusty and gaseous) accretion disks. Results. The visibilities are between 5% and 15%, and the non-spherically symmetric emitting structure has a typical size of 100 AU. A satisfactory model is constructed to reproduce the visibility spectra for each (u ,v ) point. It fits the N-band flux spectrum, the midinfrared slope, the far-infrared peak, and the (sub)mm regime of the SED. It produces a 350 μm morphology consistent with the observations. Conclusions. The mid-infrared emission of W33A on 100 AU scales is dominated by the irradiated walls of the cavity sculpted by the

The mid- to far-infrared emission of the outbursting FU Orionis objects has been attributed either to a flared outer disk or to an infalling envelope. We revisit this issue using detailed radiative transfer calculations to model the recent, high signal-to-noise data from the IRS instrument on the {Spitzer Space Telescope}. In the case of FU Ori, we find that a physically-plausible flared disk irradiated by the central accretion disk matches the observations. Building on our previous work, our accretion disk model with outer disk irradiation by the inner disk reproduces the spectral energy distribution between ~4000 angstroms to ~40 microns. Our model is consistent with near-infrared interferometry but there are some inconsistencies with mid-infared interferometric results. Including the outer disk allows us to refine our estimate of the outer radius of the outbursting, high mass accretion rate disk in FU Ori as ~ 0.5 AU, which is a crucial parameter in assessing theories of the FU Orionis phenomenon. We are able to place an upper limit on the mass infall rate of any remnant envelope infall rate to ~ 7e-7 Msun/yr assuming a centrifugal radius of 200 AU. The FUor BBW 76 is also well modelled by a 0.6 AU inner disk and a flared outer disk. However, V1515 Cyg requires an envelope with an outflow cavity to adequately reproduce the IRS spectrum. In contrast with the suggestion by Green et al., we do not require a flattened envelope to match the observations; the inferred cavity shape is qualitatively consistent with typical protostellar envelopes. This variety of dusty structures suggests that the FU Orionis phase can be present at either early or late stages of protostellar evolution.

We develop and investigate a procedure that accounts for disk reprocessing of photons that originate in the disk itself. Surface temperatures and simple, black body spectral energy distributions (SEDs) of protostellar disks are calculated. In disks that flare with radius, reprocessing of stellar photons results in temperature profiles considerably shallower than r(sup -3/4). Including the disk as a radiation source (as in the case of actively secreting disks) along with the stellar source further flattens the temperature profile. Disks that flare strongly near the star and then smoothly curve over and become shadowed at some distance (decreasing curvature disks) exhibit nearly power-law temperature profiles which result in power-law infrared SEDs with slopes in agreement with typical observations of young stellar objects. Disk models in which the photospheric thickness is controlled by the local opacity and in which the temperature decreases with radius naturally show this shape. Uniformly flaring models do not match observations as well; progressively stronger reprocessing at larger radii leads to SEDs that flatten toward the infrared or even have a second peak at the wavelength corresponding (through the Wien law) to the temperature of the outer edge of the disk. In FU Orionis outbursting systems, the dominant source of energy is the disk itself. The details of the reprocessing depend sensitively on the assumed disk shape and emitted temperature profile. The thermal instability outburst models of Bell Lin reproduce trends in the observed SEDs of Fuors with T varies as r(sup -3/4) in the inner disk (r approx. less than 0.25au corresponding to lambda approx. less than 10 microns) and T varies as r(sup -1/2) in the outer disk. Surface irradiation during outburst and quiescence is compared in the region of planet formation (1 - 10 au). The contrast between the two phases is diminished by the importance of the reprocessing of photons from the relatively high mass flux outer disk (dot-M = 10(exp -5) solar mass/yr) which is present during both outburst and quiescence.

Context. Studies of the dust distribution, composition, and evolution of protoplanetary disks provide clues for understanding planet formation. However, little is known about the innermost regions of disks where telluric planets are expected to form. Aims. We aim constrain the geometry of the inner disk of the T Tauri star RY Lup by combining spectro-photometric data and interferometric observations in the near-infrared (NIR) collected at the Very Large Telescope Interferometer. We use PIONIER data from the ESO archive and GRAVITY data that were obtained in June 2017 with the four 8m telescopes. Methods. We use a parametric disk model and the 3D radiative transfer code MCFOST to reproduce the spectral energy distribution (SED) and match the interferometric observations. MCFOST produces synthetic SEDs and intensity maps at different wavelengths from which we compute the modeled interferometric visibilities and closure phases through Fourier transform. Results. To match the SED from the blue to the millimetric range, our model requires a stellar luminosity of 2.5 L⊙, higher than any previously determined values. Such a high value is needed to accommodate the circumstellar extinction caused by the highly inclined disk, which has been neglected in previous studies. While using an effective temperature of 4800 K determined through high-resolution spectroscopy, we derive a stellar radius of 2.29 R⊙. These revised fundamental parameters, when combined with the mass estimates available (in the range 1.3–1.5 M⊙), lead to an age of 0.5–2.0 Ma for RY Lup, in better agreement with the age of the Lupus association than previous determinations. Our disk model (that has a transition disk geometry) nicely reproduces the interferometric GRAVITY data and is in good agreement with the PIONIER ones. We derive an inner rim location at 0.12 au from the central star. This model corresponds to an inclination of the inner disk of 50°, which is in mild tension with previous determinations of a more inclined outer disk from SPHERE (70° in NIR) and ALMA (67 ± 5°) images, but consistent with the inclination determination from the ALMA CO spectra (55 ± 5°). Increasing the inclination of the inner disk to 70° leads to a higher line-of-sight extinction and therefore requires a higher stellar luminosity of 4.65 L⊙ to match the observed flux levels. This luminosity would translate to a stellar radius of 3.13 R⊙, leading to an age of 2–3 Ma, and a stellarmass of about 2 M⊙, in disagreement with the observed dynamical mass estimate of 1.3–1.5 M⊙. Critically, this high-inclination inner disk model also fails to reproduce the visibilities observed with GRAVITY. Conclusions. The inner dust disk, as traced by the GRAVITY data, is located at a radius in agreement with the dust sublimation radius. An ambiguity remains regarding the respective orientations of the inner and outer disk, coplanar and mildly misaligned, respectively.As our datasets are not contemporary and the star is strongly variable, a deeper investigation will require a dedicated multi-technique observing campaign.

A new class of pre-main sequence objects has been recently identified as pre-transitional disks. They present near-infrared excess coupled to a flux deficit at about 10 microns and a rising mid-infrared and far-infrared spectrum. These features suggest a disk structure with inner and outer dust components, separated by a dust-depleted region (or gap). We here report on the first interferometric observations of the disk around the Herbig Ae star HD 139614. Its infrared spectrum suggests a flared disk, and presents pre-transitional features,namely a substantial near-infrared excess accompanied by a dip around 6 microns and a rising mid-infrared part. In this framework, we performed a study of the spectral energy distribution (SED) and the mid-infrared VLTI/MIDI interferometric data to constrain thespatial structure of the inner dust disk region and assess its possibly multi-component structure. We based our work on a temperature-gradient disk model that includes dust opacity. While we could not reproduce the SED and interferometric visibilities with a one-component disk, a better agreement was obtained with a two-component disk model composed of an optically thin inner disk extending from 0.22 to 2.3 au, a gap, and an outer temperature-gradient disk starting at 5.6 au. Therefore, our modeling favors an extended and optically thin inner dust component and in principle rules out the possibility that the near-infrared excess originates only from a spatially confined region. Moreover, the outer disk is characterized by a very steep temperature profile and a temperature higher than 300 K at its inner edge. This suggests the existence of a warm component corresponding to a scenario where the inner edge of the outer disk is directly illuminated by the central star. This is an expected consequence of the presence of a gap, thus indicative of a pre-transitional structure.

We study protostellar envelope and outflow evolution using Hubble Space Telescope NICMOS or WFC3 images of 304 protostars in the Orion molecular clouds. These near-IR images resolve structures in the envelopes delineated by the scattered light of the central protostars with 80 au resolution, and they complement the 1.2 μm to 870 μm spectral energy distributions (SEDs) obtained with the Herschel Orion Protostar Survey program. Based on their 1.60 μm morphologies, we classify the protostars into five categories: nondetections, point sources without nebulosity, bipolar cavity sources, unipolar cavity sources, and irregulars. We find point sources without associated nebulosity are the most numerous, and show through monochromatic Monte Carlo radiative transfer modeling that this morphology occurs when protostars are observed at low inclinations or have low envelope densities. We also find that the morphology is correlated with the SED-determined evolutionary class, with Class 0 protostars more likely to be nondetections, Class I protostars to show cavities, and flat-spectrum protostars to be point sources. Using an edge detection algorithm to trace the projected edges of the cavities, we fit power laws to the resulting cavity shapes, thereby measuring the cavity half-opening angles and power-law exponents. We find no evidence for the growth of outflow cavities as protostars evolve through the Class I protostar phase, in contradiction with previous studies of smaller samples. We conclude that the decline of mass infall with time cannot be explained by the progressive clearing of envelopes by growing outflow cavities. Furthermore, the low star formation efficiency inferred for molecular cores cannot be explained by envelope clearing alone.

Multi-Wavelength Properties of Submillimeter-Selected Galaxies

The brightest steady sources of radiation in the universe, active galactic nuclei (AGNs), are powered by gas accretion on to a central supermassive black hole (SMBH). The large sizes and accretion rates implicated in AGN accretion discs are expected to lead to gravitational instability and fragmentation, effectively cutting off mass inflow to the SMBH. Radiative feedback from disc-embedded stars has been invoked to yield marginally stable, steady-state solutions in the outer discs. Here, we examine the consequences of this star formation with a semi-analytical model in which stellar-mass black hole (sBH) remnants in the disc provide an additional source of stabilizing radiative feedback. Assuming star formation seeds the embedded sBH population, we model the time-evolving feedback from both stars and the growing population of accreting sBHs. We find that in the outer disc, the luminosity of the sBHs quickly dominates that of their parent stars. However, because sBHs consume less gas than stars to stabilize the disc, the presence of the sBHs enhances the mass flux to the inner disc. As a result, star formation persists over the lifetime of the AGN, damped in the outer disc, but amplified in a narrow ring in the inner disc. Heating from the embedded sBHs significantly modifies the disc’s temperature profile and hardens its spectral energy distribution, and direct emission from the sBHs adds a new hard X-ray component.

As an initial investigation into the long-term evolution of protostellar disks, we explore the conditions required to explain the large outbursts of disk accretion seen in some young stellar objects. We use one-dimensional time-dependent disk models with a phenomenological treatment of the magnetorotational instability (MRI) and gravitational torques to follow disk evolution over long timescales. Comparison with our previous two-dimensional disk model calculations (Zhu et al. 2009b, Z2009b) indicates that the neglect of radial effects and two-dimensional disk structure in the one-dimensional case makes only modest differences in the results; this allows us to use the simpler models to explore parameter space efficiently. We find that the mass infall rates typically estimated for low-mass protostars generally result in AU-scale disk accretion outbursts, as predicted by our previous analysis (Zhu et al. 2009a,Z2009a). We also confirm quasi-steady accretion behavior for high mass infall rates if the values of $\alpha$-parameter for the magnetorotational instability is small, while at this high accretion rate convection from the thermal instability may lead to some variations. We further constrain the combinations of the $\alpha$-parameter and the MRI critical temperature, which can reproduce observed outburst behavior. Our results suggest that dust sublimation may be connected with full activation of the MRI. This is consistent with the idea that small dust captures ions and electrons to suppress the MRI. In a later paper we will explore both long-term outburst and disk evolution with this model, allowing for infall from protostellar envelopes with differing angular momenta.

We present Spitzer and Chandra observations of the nearby (~260 pc) embedded stellar cluster in the Serpens Cloud Core. We observed, using Spitzer's IRAC and MIPS instruments, in six wavelength bands from 3 to 70 ${\mu}m$, to detect thermal emission from circumstellar disks and protostellar envelopes, and to classify stars using color-color diagrams and spectral energy distributions (SEDs). These data are combined with Chandra observations to examine the effects of circumstellar disks on stellar X-ray properties. Young diskless stars were also identified from their increased X-ray emission. We have identified 138 YSOs in Serpens: 22 class 0/I, 16 flat spectrum, 62 class II, 17 transition disk, and 21 class III stars; 60 of which exhibit X-ray emission. Our primary results are the following: 1.) ten protostars detected previously in the sub-millimeter are detected at lambda 1. This may be the result of grain growth through coagulation and/or the accretion of volatiles in the Serpens cloud core.