• The point at which a protoplanetary disk becomes a debris disk is difficult to identify. To better understand this, here we study the $\sim$40~AU separation binary T~54 in the Chamaeleon I cloud. We derive a K5 spectral type for T~54~A (which dominates the emission of the system) and an age of $\sim$2~Myr. However, the dust disk properties of T~54 are consistent with those of debris disks seen around older and earlier-type stars. At the same time, T~54 has evidence of gas remaining in the disk as indicated by [Ne II], [Ne III], and [O I] line detections. We model the spectral energy distribution of T~54 and estimate that $\sim$3$\times$10$^{-3}$ Earth-masses of small dust grains ($<$0.25~$\mu$m) are present in an optically thin circumbinary disk along with at least $\sim$3$\times$10$^{-7}$ Earth-masses of larger ($>$10~$\mu$m) grains within a circumprimary disk. Assuming a solar-like mixture, we use Ne line luminosities to place a minimum limit on the gas mass of the disk ($\sim$3$\times$10$^{-4}$ Earth-masses) and derive a gas-to-dust mass ratio of $\sim$0.1. We do not detect substantial accretion, but we do see H$\alpha$ in emission in one epoch, suggestive that there may be intermittent dumping of small amounts of matter onto the star. Considering the low dust mass, the presence of gas, and young age of T~54, we conclude that this system is on the bridge between the protoplanetary and debris disk stages.
  • Transitional discs are a class of circumstellar discs around young stars with extensive clearing of dusty material within their inner regions on 10s of au scales. One of the primary candidates for this kind of clearing is the formation of planet(s) within the disc that then accrete or clear their immediate area as they migrate through the disc. Our sample included eight transitional discs. Using the Keck/NIRC2 instrument we utilised the Sparse Aperture Masking (SAM) interferometry technique to search for asymmetries indicative of ongoing planet formation. We searched for close-in companions using both model fitting and interferometric image reconstruction techniques. Using simulated data, we derived diagnostics that helped us to distinguish between point sources and extended asymmetric disc emission. In addition, we investigated the degeneracy between the contrast and separation that appear for marginally resolved companions. We found FP Tau to contain a previously unseen disc wall, and DM Tau, LkHa 330, and TW Hya to contain an asymmetric signal indicative of point source-like emission. We placed upper limits on the contrast of a companion in RXJ1842.9-3532 and V2246 Oph. We ruled the asymmetry signal in RXJ1615.3-3255 and V2062 Oph to be false positives. In the cases where our data indicated a potential companion we computed estimates for the value of $M_c \dot M_c$ and found values in the range of $10^{-5} - 10^{-3} M^2_J yr^{-1}$. We found significant asymmetries in four targets. Of these, three were consistent with companions. We resolved a previously unseen gap in the disc of FP Tau extending inwards from approximately 10 au.
  • We present key results from the Herschel Orion Protostar Survey (HOPS): 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 sub-millimeter photometry from APEX, our SEDs cover 1.2-870 $\mu$m and sample the peak of the protostellar envelope emission at ~100 $\mu$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 30400 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.
  • To characterize the mechanisms of planet formation it is crucial to investigate the properties and evolution of protoplanetary disks around young stars, where the initial conditions for the growth of planets are set. Our goal is to study grain growth in the disk of the young, intermediate mass star HD163296 where dust processing has already been observed, and to look for evidence of growth by ice condensation across the CO snowline, already identified in this disk with ALMA. Under the hypothesis of optically thin emission we compare images at different wavelengths from ALMA and VLA to measure the opacity spectral index across the disk and thus the maximum grain size. We also use a Bayesian tool based on a two-layer disk model to fit the observations and constrain the dust surface density. The measurements of the opacity spectral index indicate the presence of large grains and pebbles ($\geq$1 cm) in the inner regions of the disk (inside $\sim$50 AU) and smaller grains, consistent with ISM sizes, in the outer disk (beyond 150 AU). Re-analysing ALMA Band 7 Science Verification data we find (radially) unresolved excess continuum emission centered near the location of the CO snowline at $\sim$90 AU. Our analysis suggests a grain size distribution consistent with an enhanced production of large grains at the CO snowline and consequent transport to the inner regions. Our results combined with the excess in infrared scattered light found by Garufi et al. (2014) suggests the presence of a structure at 90~AU involving the whole vertical extent of the disk. This could be evidence for small scale processing of dust at the CO snowline.
  • Theoretical models of grain growth predict dust properties to change as a function of protoplanetary disk radius, mass, age and other physical conditions. We lay down the methodology for a multi-wavelength analysis of (sub-)mm and cm continuum interferometric observations to constrain self-consistently the disk structure and the radial variation of the dust properties. The computational architecture is massively parallel and highly modular. The analysis is based on the simultaneous fit in the uv-plane of observations at several wavelengths with a model for the disk thermal emission and for the dust opacity. The observed flux density at the different wavelengths is fitted by posing constraints on the disk structure and on the radial variation of the grain size distribution. We apply the analysis to observations of three protoplanetary disks (AS 209, FT Tau, DR Tau) for which a combination of spatially resolved observations in the range ~0.88mm to ~10mm is available (from SMA, CARMA, and VLA), finding evidence of a decreasing maximum dust grain size (a_max) with radius. We derive large a_max values up to 1 cm in the inner disk between 15 and 30 AU and smaller grains with a_max~1 mm in the outer disk (R > 80AU). In this paper we develop a multi-wavelength analysis that will allow this missing quantity to be constrained for statistically relevant samples of disks and to investigate possible correlations with disk or stellar parameters.
  • Correlations between the accretion luminosity and emission line luminosities (L_acc and L_line) of pre-main sequence (PMS) stars have been published for many different spectral lines, which are used to estimate accretion rates. Despite the origin of those correlations is unknown, this could be attributed to direct or indirect physical relations between the emission line formation and the accretion mechanism. This work shows that all (near-UV/optical/near-IR) L_acc-L_line correlations are the result of the fact that the accretion luminosity and the stellar luminosity (L_star) are correlated, and are not necessarily related with the physical origin of the line. Synthetic and observational data are used to illustrate how the L_acc-L_line correlations depend on the L_acc-L_star relationship. We conclude that because PMS stars show the L_acc-L_star correlation immediately implies that L_acc also correlates with the luminosity of all emission lines, for which the L_acc-L_line correlations alone do not prove any physical connection with accretion but can only be used with practical purposes to roughly estimate accretion rates. When looking for correlations with possible physical meaning, we suggest that L_acc/L_star and L_line/L_star should be used instead of L_acc and L_line. Finally, the finding that L_acc has a steeper dependence on L_star for T-Tauri stars than for intermediate-mass Herbig Ae/Be stars is also discussed. That is explained from the magnetospheric accretion scenario and the different photospheric properties in the near-UV.
  • We present Herschel Space Observatory PACS spectra of T Tauri stars, in which we detect amorphous and crystalline water ice features. Using irradiated accretion disk models, we determine the disk structure and ice abundance in each of the systems. Combining a model-independent comparison of the ice feature strength and disk size with a detailed analysis of the model ice location, we estimate that the ice emitting region is at disk radii >30AU, consistent with a proto-Kuiper belt. Vertically, the ice emits most below the photodesorption zone, consistent with Herschel observations of cold water vapor. The presence of crystallized water ice at a disk location a) colder than its crystallization temperature and b) where it should have been re-amorphized in ~1 Myr suggests that localized generation is occurring; the most likely cause appears to be micrometeorite impact or planetesimal collisions. Based on simple tests with UV models and different ice distributions, we suggest that the SED shape from 20 to 50 micron may probe the location of the water ice snow line in the disk upper layers. This project represents one of the first extra-solar probes of the spatial structure of the cometary ice reservoir thought to deliver water to terrestrial planets.
  • We present spectra of 13 T Tauri stars in the Taurus-Auriga star-forming region showing emission in Spitzer Space Telescope Infrared Spectrograph (IRS) 5-7.5 micron spectra from water vapor and absorption from other gases in these stars' protoplanetary disks. Seven stars' spectra show an emission feature at 6.6 microns due to the nu_2 = 1-0 bending mode of water vapor, with the shape of the spectrum suggesting water vapor temperatures > 500 K, though some of these spectra also show indications of an absorption band, likely from another molecule. This water vapor emission contrasts with the absorption from warm water vapor seen in the spectrum of the FU Orionis star V1057 Cyg. The other six of the thirteen stars have spectra showing a strong absorption band, peaking in strength at 5.6-5.7 microns, which for some is consistent with gaseous formaldehyde (H2CO) and for others is consistent with gaseous formic acid (HCOOH). There are indications that some of these six stars may also have weak water vapor emission. Modeling of these stars' spectra suggests these gases are present in the inner few AU of their host disks, consistent with recent studies of infrared spectra showing gas in protoplanetary disks.
  • For over a decade, the structure of the inner cavity in the transition disk of TW Hydrae has been a subject of debate. Modeling the disk with data obtained at different wavelengths has led to a variety of proposed disk structures. Rather than being inconsistent, the individual models might point to the different faces of physical processes going on in disks, such as dust growth and planet formation. Our aim is to investigate the structure of the transition disk again and to find to what extent we can reconcile apparent model differences. A large set of high-angular-resolution data was collected from near-infrared to centimeter wavelengths. We investigated the existing disk models and established a new self-consistent radiative-transfer model. A genetic fitting algorithm was used to automatize the parameter fitting. Simple disk models with a vertical inner rim and a radially homogeneous dust composition from small to large grains cannot reproduce the combined data set. Two modifications are applied to this simple disk model: (1) the inner rim is smoothed by exponentially decreasing the surface density in the inner ~3 AU, and (2) the largest grains (>100 um) are concentrated towards the inner disk region. Both properties can be linked to fundamental processes that determine the evolution of protoplanetary disks: the shaping by a possible companion and the different regimes of dust-grain growth, respectively. The full interferometric data set from near-infrared to centimeter wavelengths requires a revision of existing models for the TW Hya disk. We present a new model that incorporates the characteristic structures of previous models but deviates in two key aspects: it does not have a sharp edge at 4 AU, and the surface density of large grains differs from that of smaller grains. This is the first successful radiative-transfer-based model for a full set of interferometric data.
  • The dust sublimation walls of disks around T Tauri stars represent a directly observable cross-section through the disk atmosphere and midplane. Their emission properties can probe the grain size distribution and composition of the innermost regions of the disk, where terrestrial planets form. Here we calculate the inner dust sublimation wall properties for four classical T Tauri stars with a narrow range of spectral types and inclination angles and a wide range of mass accretion rates to determine the extent to which the walls are radially curved. Best-fits to the near- and mid-IR excesses are found for curved, 2-layer walls in which the lower layer contains larger, hotter, amorphous pyroxene grains with Mg/(Mg+Fe)=0.6 and the upper layer contains submicron, cooler, mixed amorphous olivine and forsterite grains. As the mass accretion rates decrease from 10^(-8) to 10^(-10) Msol/yr, the maximum grain size in the lower layer decreases from 3 to 0.5 microns. We attribute this to a decrease in fragmentation and turbulent support for micron-sized grains with decreasing viscous heating. The atmosphere of these disks is depleted of dust with dust-gas mass ratios 1x10^(-4) of the ISM value, while the midplane is enhanced to 8 times the ISM value. For all accretion rates, the wall contributes at least half of the flux in the optically thin 10 micron silicate feature. Finally, we find evidence for an iron gradient in the disk, suggestive of that found in our solar system.
  • Using NASA IRTF SpeX data from 0.8 to 4.5 $\mu$m, we determine self-consistently the stellar properties and excess emission above the photosphere for a sample of classical T Tauri stars (CTTS) in the Taurus molecular cloud with varying degrees of accretion. This process uses a combination of techniques from the recent literature as well as observations of weak-line T Tauri stars (WTTS) to account for the differences in surface gravity and chromospheric activity between the TTS and dwarfs, which are typically used as photospheric templates for CTTS. Our improved veiling and extinction estimates for our targets allow us to extract flux-calibrated spectra of the excess in the near-infrared. We find that we are able to produce an acceptable parametric fit to the near-infrared excesses using a combination of up to three blackbodies. In half of our sample, two blackbodies at temperatures of 8000 K and 1600 K suffice. These temperatures and the corresponding solid angles are consistent with emission from the accretion shock on the stellar surface and the inner dust sublimation rim of the disk, respectively. In contrast, the other half requires three blackbodies at 8000, 1800, and 800 K, to describe the excess. We interpret the combined two cooler blackbodies as the dust sublimation wall with either a contribution from the disk surface beyond the wall or curvature of the wall itself, neither of which should have single-temperature blackbody emission. In these fits, we find no evidence of a contribution from optically thick gas inside the inner dust rim.
  • High-energy radiation from T Tauri stars (TTS) influences the amount and longevity of gas in disks, thereby playing a crucial role in the creation of gas giant planets. Here we probe the high-energy ionizing radiation from TTS using high-resolution mid-infrared (MIR) Spitzer IRS Neon forbidden line detections in a sample of disks from IC 348, NGC 2068, and Chamaeleon. We report three new detections of [Ne III] from CS Cha, SZ Cha, and T 54, doubling the known number of [Ne III] detections from TTS. Using [Ne III]-to-[Ne II] ratios in conjunction with X-ray emission measurements, we probe high-energy radiation from TTS. The majority of previously inferred [Ne III]/[Ne II] ratios based on [Ne III] line upper limits are significantly less than 1, pointing to the dominance of either X-ray radiation or soft Extreme-Ultraviolet (EUV) radiation in producing these lines. Here we report the first observational evidence for hard EUV dominated Ne forbidden line production in a T Tauri disk: [Ne III]/[Ne II]~1 in SZ Cha. Our results provide a unique insight into the EUV emission from TTS, by suggesting that EUV radiation may dominate the creation of Ne forbidden lines, albeit in a minority of cases.
  • We present Herschel Space Observatory PACS spectra of GQ Lup, a protoplanetary disk in the Lupus star-forming region. Through SED fitting from 0.3{\mu}m to 1.3mm, we construct a self-consistent model of this system's temperature and density structures, finding that although it is 3 Myr old, its dust has not settled to the midplane substantially. The disk has a radial gradient in both the silicate dust composition and grain size, with large amorphous grains in the upper layers of the inner disk and an enhancement of submicron, crystalline grains in the outer disk. We detect an excess of emission in the Herschel PACS B2A band near 63{\mu}m and model it with a combination of {\sim}15 to 70{\mu}m crystalline water ice grains with a size distribution consistent with ice recondensation-enhanced grain growth and a mass fraction half of that of our solar system. The combination of crystalline water ice and silicates in the outer disk is suggestive of disk-wide heating events or planetesimal collisions. If confirmed, this would be the first detection of water ice by Herschel.
  • We present an Advanced Camera for Surveys/ Solar Blind Channel far-ultraviolet (FUV) study of \h2 gas in 12 weak T Tauri stars in nearby star-forming regions. The sample consists of sources which have no evidence of inner disk dust. Our new FUV spectra show that in addition to the dust, the gas is depleted from the inner disk. This sample is combined with a larger FUV sample of accretors and non-accretors with ages between 1 and 100 Myr, showing that as early as 1--3 Myr, systems both with and without gas are found. Possible mechanisms for depleting gas quickly include viscous evolution, planet formation and photoevaporation by stellar radiation fields. Since these mechanisms alone cannot account for the lack of gas at 1--3 Myr, it is likely that the initial conditions (e.g. initial disk mass or core angular momentum) contribute to the variety of disks observed at any age. We estimate the angular momentum of a cloud needed for most of the mass to fall very close to the central object and compare this to models of the expected distribution of angular momenta. Up to 20% of cloud cores have low enough angular momenta to form disks with the mass close to the star, which would then accrete quickly; this percentage is similar to the fraction of diskless stars in the youngest star forming regions. With our sample, we characterize the chromospheric contribution to the FUV luminosity and find that $L_{FUV}/L_{bol}$ saturates at $\sim10^{-4.1}$.
  • Two decades ago "transitional disks" described spectral energy distributions (SEDs) of T Tauri stars with small near-IR excesses, but significant mid- and far-IR excesses. Many inferred this indicated dust-free holes in disks, possibly cleared by planets. Recently, this term has been applied disparately to objects whose Spitzer SEDs diverge from the expectations for a typical full disk. Here we use irradiated accretion disk models to fit the SEDs of 15 such disks in NGC 2068 and IC 348. One group has a "dip" in infrared emission while the others' continuum emission decreases steadily at all wavelengths. We find that the former have an inner disk hole or gap at intermediate radii in the disk and we call these objects "transitional" and pre-transitional" disks, respectively. For the latter group, we can fit these SEDs with full disk models and find that millimeter data are necessary to break the degeneracy between dust settling and disk mass. We suggest the term "transitional" only be applied to objects that display evidence for a radical change in the disk's radial structure. Using this definition, we find that transitional and pre-transitional disks tend to have lower mass accretion rates than full disks and that transitional disks have lower accretion rates than pre-transitional disks. These reduced accretion rates onto the star could be linked to forming planets. Future observations of transitional and pre-transitional disks will allow us to better quantify the signatures of planet formation in young disks.
  • This work aims to derive accretion rates for a sample of 38 HAeBe stars. We apply magnetospheric accretion (MA) shock modelling to reproduce the observed Balmer excesses. We look for possible correlations with the strength of the Halpha, [OI]6300, and Brgamma emission lines. The median mass accretion rate is 2 x 10^-7 Msun yr^-1 in our sample. The model fails to reproduce the large Balmer excesses shown by the four hottest stars (T* > 12000 K). We derive Macc propto M*^5 and Lacc propto L*^1.2 for our sample, with scatter. Empirical calibrations relating the accretion and the Halpha, [OI]6300, and Brgamma luminosities are provided. The slopes in our expressions are slightly shallower than those for lower mass stars, but the difference is within the uncertainties, except for the [OI]6300 line. The Halpha 10% width is uncorrelated with Macc, unlike for the lower mass regime. The mean Halpha width shows higher values as the projected rotational velocities of HAe stars increase, which agrees with MA. The accretion rate variations in the sample are typically lower than 0.5 dex on timescales of days to months, Our data suggest that the changes in the Balmer excess are uncorrelated to the simultaneous changes of the line luminosities. The Balmer excesses and Halpha line widths of HAe stars can be interpreted within the context of MA, which is not the case for several HBes. The steep trend relating Macc and M* can be explained from the mass-age distribution characterizing HAeBe stars. The line luminosities used for low-mass objects are also valid to estimate typical accretion rates for the intermediate-mass regime under similar empirical expressions. However, we suggest that several of these calibrations are driven by the stellar luminosity.
  • We present 5 to 36 micron mid-infrared spectra of 82 young stars in the ~2 Myr old Chamaeleon I star-forming region, obtained with the Spitzer Infrared Spectrograph (IRS). We have classified these objects into various evolutionary classes based on their spectral energy distributions and the spectral features seen in the IRS spectra. We have analyzed the mid-IR spectra of Class II objects in Chamaeleon I in detail, in order to study the vertical and radial structure of the protoplanetary disks surrounding these stars. We find evidence for substantial dust settling in most protoplanetary disks in Chamaeleon I. We have identified several disks with altered radial structures in Chamaeleon I, among them transitional disk candidates which have holes or gaps in their disks. Analysis of the silicate emission features in the IRS spectra of Class II objects in Chamaeleon I shows that the dust grains in these disks have undergone significant processing (grain growth and crystallization). However, disks with radial holes/gaps appear to have relatively unprocessed grains. We further find the crystalline dust content in the inner (< 1-2 AU) and the intermediate (< 10 AU) regions of the protoplanetary disks to be tightly correlated. We also investigate the effects of accretion and stellar multiplicity on the disk structure and dust properties. Finally, we compare the observed properties of protoplanetary disks in Cha I with those in slightly younger Taurus and Ophiuchus regions and discuss the effects of disk evolution in the first 1-2 Myr.
  • We present a Spitzer IRS study of variability in 14 T Tauri stars in the Taurus and Chamaeleon star-forming regions. The sample is composed of transitional and pre-transitional objects which contain holes and gaps in their disks. We detect variability between 5-38 microns in all but two of our objects on timescales of 2-3 years. Most of the variability observed can be classified as seesaw behavior, whereby the emission at shorter wavelengths varies inversely with the emission at longer wavelengths. For many of the objects we can reasonably reproduce the observed variability using irradiated disk models, particularly by changing the height of the inner disk wall by ~20%. When the inner wall is taller, the emission at the shorter wavelengths is higher since the inner wall dominates the emission at 2-8 microns. The taller inner wall casts a larger shadow on the outer disk wall, leading to less emission at wavelengths beyond 20 microns where the outer wall dominates. We discuss how the possible presence of planets in these disks could lead to warps which cause changes in the height of the inner wall. We also find that crystalline silicates are common in the outer disks of our objects and that in the four disks in the sample with the most crystalline silicates, variability on timescales of 1 week is present. In addition to explaining the infrared variability described above, planets can create shocks and collisions which can crystallize the dust and lead to short timescale variability.
  • Spectral observations of classical T Tauri stars show a wide range of line profiles, many of which reveal signs of matter inflow and outflow. Halpha is the most commonly observed line profile due to its intensity, and it is highly dependent on the characteristics of the surrounding environment of these stars. Our aim is to analyze how the Halpha line profile is affected by the various parameters of our model which contains both the magnetospheric and disk wind contributions to the Halpha flux. We used a dipolar axisymmetric stellar magnetic field to model the stellar magnetosphere and a modified Blandford & Payne model was used in our disk wind region. A three-level atom with continuum was used to calculate the required Hydrogen level populations. We use the Sobolev approximation and a ray-by-ray method to calculate the integrated line profile. Through an extensive study of the model parameter space, we have investigated the contribution of many of the model parameters on the calculated line profiles. Our results show that the Halpha line is strongly dependent on the densities and temperatures inside the magnetosphere and the disk wind region. The bulk of the flux comes, most of the time, from the magnetospheric component for standard classical T Tauri stars parameters, but the disk wind contribution becomes more important as the mass accretion rate, the temperatures and densities inside the disk wind increase. We have also found that most of the disk wind contribution to the Halpha line is emitted at the innermost region of the disk wind. Models that take into consideration both inflow and outflow of matter are a necessity to fully understand and describe classical T Tauri stars.
  • We have analyzed nearly all images of the Taurus star-forming region at 3.6-24um that were obtained during the cryogenic mission of the Spitzer Space Telescope (46 deg^2) and have measured photometry for all known members of the region that are within these data, corresponding to 348 sources. We have classified the members of Taurus according to whether they show evidence of disks and envelopes (classes I, II, and III). The disk fraction in Taurus is 75% for solar-mass stars and declines to 45% for low-mass stars and brown dwarfs (0.01-0.3 M_sun). This dependence on stellar mass is similar to that measured for Cha I, although the disk fraction in Taurus is slightly higher overall, probably because of its younger age (1 vs. 2-3 Myr). In comparison, the disk fraction for solar-mass stars is much lower (20%) in IC 348 and Sigma Ori, which are denser than Taurus and Cha I and are roughly coeval with the latter. These data indicate that disk lifetimes for solar-mass stars are longer in regions that have lower stellar densities. Through an analysis of multiple epochs of photometry that are available for ~200 Taurus members, we find that stars with disks exhibit significantly greater mid-IR variability than diskless stars. Finally, we have used our data in Taurus to refine the criteria for primordial, evolved, and transitional disks. The number ratio of evolved and transitional disks to primordial disks in Taurus is 15/98 for K5-M5, indicating a timescale of 0.15 x tau(primordial)=0.45 Myr for the clearing of the inner regions of optically thick disks. After applying the same criteria to older clusters (2-10 Myr), we find that the proportions of evolved and transitional disks in those populations are consistent with the measurements in Taurus when their star formation histories are properly taken into account. ERRATUM: In Table 7, we inadvertently omitted the spectral type bins in which class II sources were placed in Table 8 based on their bolometric luminosities (applies only to stars that lack spectroscopic classifications). The bins were K6-M3.5 for FT Tau, DK Tau B, and IRAS 04370+2559, M3.5-M6 for IRAS 04200+2759, IT Tau B, and ITG 1, and M6-M8 for IRAS 04325+2402 C. In addition, the values of K_s-[3.6] in Table 13 and Figure 26 for spectral types of M4-M9 are incorrect. We present corrected versions of Table 13 and Figure 26.
  • In the past few years, several disks with inner holes that are empty of small dust grains have been detected and are known as transitional disks. Recently, Spitzer has identified a new class of "pre-transitional disks" with gaps; these objects have an optically thick inner disk separated from an optically thick outer disk by an optically thin disk gap. A near-infrared spectrum provided the first confirmation of a gap in the pre-transitional disk of LkCa 15 by verifying that the near-infrared excess emission in this object was due to an optically thick inner disk. Here we investigate the difference between the nature of the inner regions of transitional and pre-transitional disks using the same veiling-based technique to extract the near-infrared excess emission above the stellar photosphere. We show that the near-infrared excess emission of the previously identified pre-transitional disks of LkCa 15 and UX Tau A in Taurus as well as the newly identified pre-transitional disk of ROX 44 in Ophiuchus can be fit with an inner disk wall located at the dust destruction radius. We also model the broad-band SEDs of these objects, taking into account the effect of shadowing by the inner disk on the outer disk, considering the finite size of the star. The near-infrared excess continua of these three pre-transitional disks, which can be explained by optically thick inner disks, are significantly different from that of the transitional disks of GM Aur, whose near-infrared excess continuum can be reproduced by emission from sub-micron-sized optically thin dust, and DM Tau, whose near-infrared spectrum is consistent with a disk hole that is relatively free of small dust. The structure of pre-transitional disks may be a sign of young planets forming in these disks and future studies of pre-transitional disks will provide constraints to aid in theoretical modeling of planet formation.
  • Upcoming surveys with the Herschel Space Observatory will yield far-IR photometry of large samples of young stellar objects, which will require careful interpretation. We investigate the color and luminosity diagnostics based on Herschel broad-band filters to identify and discern the properties of low-mass protostars. We compute a grid of 2,016 protostars in various physical congurations, present the expected flux densities and flux density ratios for this grid of protostars, and compare Herschel observations of three protostars to the model results. These provide useful constraints on the range of colors and fluxes of protostar in the Herschel filters. We find that Herschel data alone is likely a useful diagnostic of the envelope properties of young stars
  • We present 70 and 160 micron Herschel science demonstration images of a field in the Orion A molecular cloud that contains the prototypical Herbig-Haro objects HH 1 and 2, obtained with the Photodetector Array Camera and Spectrometer (PACS). These observations demonstrate Herschel's unprecedented ability to study the rich population of protostars in the Orion molecular clouds at the wavelengths where they emit most of their luminosity. The four protostars previously identified by Spitzer 3.6-40 micron imaging and spectroscopy are detected in the 70 micron band, and three are clearly detected at 160 microns. We measure photometry of the protostars in the PACS bands and assemble their spectral energy distributions (SEDs) from 1 to 870 microns with these data, Spitzer spectra and photometry, 2MASS data, and APEX sub-mm data. The SEDs are fit to models generated with radiative transfer codes. From these fits we can constrain the fundamental properties of the protostars. We find luminosities in the range 12-84 L_sun and envelope densities spanning over two orders of magnitude. This implies that the four protostars have a wide range of envelope infall rates and evolutionary states: two have dense, infalling envelopes, while the other two have only residual envelopes. We also show the highly irregular and filamentary structure of the cold dust and gas surrounding the protostars as traced at 160 microns.
  • The NGC 1999 reflection nebula features a dark patch with a size of ~10,000 AU, which has been interpreted as a small, dense foreground globule and possible site of imminent star formation. We present Herschel PACS far-infrared 70 and 160mum maps, which reveal a flux deficit at the location of the globule. We estimate the globule mass needed to produce such an absorption feature to be a few tenths to a few Msun. Inspired by this Herschel observation, we obtained APEX LABOCA and SABOCA submillimeter continuum maps, and Magellan PANIC near-infrared images of the region. We do not detect a submillimer source at the location of the Herschel flux decrement; furthermore our observations place an upper limit on the mass of the globule of ~2.4x10^-2 Msun. Indeed, the submillimeter maps appear to show a flux depression as well. Furthermore, the near-infrared images detect faint background stars that are less affected by extinction inside the dark patch than in its surroundings. We suggest that the dark patch is in fact a hole or cavity in the material producing the NGC 1999 reflection nebula, excavated by protostellar jets from the V 380 Ori multiple system.
  • We have used Spitzer/IRAC to conduct a photometric monitoring program of the IC1396A dark globule in order to study the mid-IR (3.6 - 8 micron) variability of the heavily embedded Young Stellar Objects (YSOs) present in that area. We obtained light curves covering a 14 day timespan with a twice daily cadence for 69 YSOs, and continuous light curves with approximately 12 second cadence over 7 hours for 38 YSOs. Typical accuracies for our relative photometry were 1-2% for the long timespan data and a few mmag, corresponding to less than 0.5%, for the 7 hour continuous "staring-mode" data. More than half of the YSOs showed detectable variability, with amplitudes from ~0.05 mag to ~0.2 mag. About thirty percent of the YSOs showed quasi-sinusoidal light curve shapes with apparent periods from 5-12 days and light curve amplitudes approximately independent of wavelength over the IRAC bandpasses. We have constructed models which simulate the time dependent spectral energy distributions of Class I and I I YSOs in order to attempt to explain these light curves. Based on these models, the apparently periodic light curves are best explained by YSO models where one or two high latitude photospheric spots heat the inner wall of the circumstellar disk, and where we view the disk at fairly large inclination angle. Disk inhomogeneities, such as increasing the height where the accretion funnel flows to the stellar hotspot, enhances the light curve modulations. The other YSOs in our sample show a range of light curve shapes, some of which are probably due to varying accretion rate or disk shadowing events. One star, IC1396A-47, shows a 3.5 hour periodic light curve; this object may be a PMS Delta Scuti star.