The Inner Rim of YSO Disks: Effects of dust grain evolution

Dust-grain growth and settling are the first steps towards planet formation. An understanding of dust physics is therefore integral to a complete theory of the planet formation process. In this paper, we explore the possibility of using the dust evaporation front in YSO disks (`the inner rim') as a probe of the dust physics operating in circumstellar disks. The geometry of the rim depends sensitively on the composition and spatial distribution of dust. Using radiative transfer and hydrostatic equilibrium calculations we demonstrate that dust growth and settling can curve the evaporation front dramatically (from a cylindrical radius of about 0.5 AU in the disk mid-plane to 1.2 AU in the disk upper layers for an A0 star). We compute synthetic images and interferometric visibilities for our representative rim models and show that the current generation of near-IR long-baseline interferometers (VLTI, CHARA) can strongly constrain the dust properties of circumstellar disks, shedding light on the relatively poorly understood processes of grain growth, settling and turbulent mixing.Comment: 26 pages, 9 figures. Accepted for publication in Ap

Strong Near-Infrared Emission Interior to the Dust-Sublimation Radius of Young Stellar Objects MWC275 and AB Aur

Using the longest optical-interferometeric baselines currently available, we have detected strong near-infrared (NIR) emission from inside the dust-destruction radius of Herbig Ae stars MWC275 and AB Aur. Our sub-milli-arcsecond resolution observations unambiguously place the emission between the dust-destruction radius and the magnetospheric co-rotation radius. We argue that this new component corresponds to hot gas inside the dust-sublimation radius, confirming recent claims based on spectrally-resolved interferometry and dust evaporation front modeling.Comment: 12 pages, 4 figures, Accepted for publication in ApJ

Mid-infrared size survey of Young Stellar Objects: Description of Keck segment-tilting experiment and basic results

The mid-infrared properties of pre-planetary disks are sensitive to the temperature and flaring profiles of disks for the regions where planet formation is expected to occur. In order to constrain theories of planet formation, we have carried out a mid-infrared (wavelength 10.7 microns) size survey of young stellar objects using the segmented Keck telescope in a novel configuration. We introduced a customized pattern of tilts to individual mirror segments to allow efficient sparse-aperture interferometry, allowing full aperture synthesis imaging with higher calibration precision than traditional imaging. In contrast to previous surveys on smaller telescopes and with poorer calibration precision, we find most objects in our sample are partially resolved. Here we present the main observational results of our survey of 5 embedded massive protostars, 25 Herbig Ae/Be stars, 3 T Tauri stars, 1 FU Ori system, and 5 emission-line objects of uncertain classification. The observed mid-infrared sizes do not obey the size-luminosity relation found at near-infrared wavelengths and a companion paper will provide further modelling analysis of this sample. In addition, we report imaging results for a few of the most resolved objects, including complex emission around embedded massive protostars, the photoevaporating circumbinary disk around MWC 361A, and the subarcsecond binaries T Tau, FU Ori and MWC 1080.Comment: Accepted by Astrophysical Journal. 38 pages. 9 figure

High Resolution K-band Spectroscopy of MWC 480 and V1331 Cyg

We present high resolution (R=25,000-35,000) K-band spectroscopy of two young stars, MWC 480 and V1331 Cyg. Earlier spectrally dispersed (R=230) interferometric observations of MWC 480 indicated the presence of an excess continuum emission interior to the dust sublimation radius, with a spectral shape that was interpreted as evidence for hot water emission from the inner disk of MWC 480. Our spectrum of V1331 Cyg reveals strong emission from CO and hot water vapor, likely arising in a circumstellar disk. In comparison, our spectrum of MWC 480 appears mostly featureless. We discuss possible ways in which strong water emission from MWC 480 might go undetected in our data. If strong water emission is in fact absent from the inner disk, as our data suggest, the continuum excess interior to the dust sublimation radius that is detected in the interferometric data must have another origin. We discuss possible physical origins for the continuum excess.Comment: 29 pages, 5 figures, to appear in Ap

Benchmark problems for continuum radiative transfer. High optical depths, anisotropic scattering, and polarisation

Solving the continuum radiative transfer equation in high opacity media requires sophisticated numerical tools. In order to test the reliability of such tools, we present a benchmark of radiative transfer codes in a 2D disc configuration. We test the accuracy of seven independently developed radiative transfer codes by comparing the temperature structures, spectral energy distributions, scattered light images, and linear polarisation maps that each model predicts for a variety of disc opacities and viewing angles. The test cases have been chosen to be numerically challenging, with midplane optical depths up 10^6, a sharp density transition at the inner edge and complex scattering matrices. We also review recent progress in the implementation of the Monte Carlo method that allow an efficient solution to these kinds of problems and discuss the advantages and limitations of Monte Carlo codes compared to those of discrete ordinate codes. For each of the test cases, the predicted results from the radiative transfer codes are within good agreement. The results indicate that these codes can be confidently used to interpret present and future observations of protoplanetary discs.Comment: 15 pages, 10 figures, accepted for publication in A&

On the interplay between flaring and shadowing in disks around Herbig Ae/Be stars

Based on the SED, Herbig stars have been categorized into two observational groups, reflecting their overall disk structure: group I members have disks with a higher degree of flaring than their group II counterparts. We investigate the 5-35 um Spitzer IRS spectra of a sample of 13 group I sources and 20 group II sources. We focus on the continuum emission to study the underlying disk geometry. We have determined the [30/13.5] and [13.5/7] continuum flux ratios. The 7-um flux excess with respect to the stellar photosphere is measured, as a marker for the strength of the near-IR emission produced by the inner disk. We have compared our data to self-consistent passive-disk model spectra, for which the same quantities were derived. We confirm the literature result that the difference in continuum emission between group I and II sources can largely be explained by a different amount of small dust grains. However, we report a strong correlation between the [30/13.5] and [13.5/7] flux ratios for Meeus group II sources. Moreover, the [30/13.5] flux ratio decreases with increasing 7-um excess for all targets in the sample. To explain these correlations with the models, we need to introduce an artificial scaling factor for the inner disk height. In roughly 50% of the Herbig Ae/Be stars in our sample, the inner disk must be inflated by a factor 2 to 3 beyond what hydrostatic calculations predict. The total disk mass in small dust grains determines the degree of flaring. We conclude, however, that for any given disk mass in small dust grains, the shadowing of the outer (tens of AU) disk is determined by the scale height of the inner disk (1 AU). The inner disk partially obscures the outer disk, reducing the disk surface temperature. Here, for the first time, we prove these effects observationally.Comment: 4 pages, 3 figures, accepted by A&

A low optical depth region in the inner disk of the HerbigAe star HR5999

Circumstellar disks surrounding young stars are known to be the birthplaces of planets, and the innermost astronomical unit is of particular interest. We present new long-baseline spectro-interferometric observations of the HerbigAe star, HR5999, obtained in the H and K bands with the AMBER instrument at the VLTI, and aim to produce near-infrared images at the sub-AU spatial scale. We spatially resolve the circumstellar material and reconstruct images using the MiRA algorithm. In addition, we interpret the interferometric observations using models that assume that the near-infrared excess is dominated by the emission of a circumstellar disk. We compare the images reconstructed from the VLTI measurements to images obtained using simulated model data. The K-band image reveals three main elements: a ring-like feature located at ~0.65 AU, a low surface brightness region inside, and a central spot. At the maximum angular resolution of our observations (1.3 mas), the ring is resolved while the central spot is only marginally resolved, preventing us from revealing the exact morphology of the circumstellar environment. We suggest that the ring traces silicate condensation, i.e., an opacity change, in a circumstellar disk around HR 5999. We build a model that includes a ring at the silicate sublimation radius and an inner disk of low surface brightness responsible for a large amount of the near-infrared continuum emission. The model successfully fits the SED, visibilities, and closure phases, and provides evidence of a low surface brightness region inside the silicate sublimation radius. This study provides additional evidence that in HerbigAe stars, there is material in a low surface brightness region, probably a low optical depth region, located inside the silicate sublimation radius and of unknown nature.Comment: 11 pages, 10 figure

Revealing the sub-AU asymmetries of the inner dust rim in the disk around the Herbig Ae star R CrA

Models predict that in the innermost AU of the disk around Herbig Ae/Be star, the dust disk forms a "puffed-up" inner rim, which should result in a strongly asymmetric brightness distribution for disks seen under intermediate inclination. Using the VLTI/AMBER long-baseline interferometer, we obtained 24 near-infrared (H- and K-band) spectro-interferometric observations on the Herbig Ae star R CrA. In the derived visibility function, we detect the signatures of an extended (25 mas) and a compact component (5.8 mas), with the compact component contributing about 2/3 of the total flux. The brightness distribution is highly asymmetric, as indicated by the strong closure phases (up to 40 deg) and the detected position angle dependence of the visibilities and closure phases. To interpret these asymmetries, we employ geometric as well as physical models, including a binary model, a skewed ring model, and a puffed-up inner rim model with a vertical or curved rim shape. Our curved puffed-up rim model can reasonably well reproduce the interferometric observables and the SED simultaneously and suggests a luminosity of 29 L_sun and the presence of relatively large (> 1.2 micron) Silicate dust grains. Perpendicular to the disk, two bow shock-like structures appear in the associated reflection nebula NGC 6729, suggesting that the resolved sub-AU size disk is the driving engine of a large-scale outflow. Detecting, for the first time, strong non-localized asymmetries in the inner regions of a Herbig Ae disk, our study supports the existence of a puffed-up inner rim in YSO disks.Comment: 17 pages; 15 figures; Accepted by A&

The Hot Inner Disk of FU Ori

We have constructed a detailed radiative transfer disk model which reproduces the main features of the spectrum of the outbursting young stellar object FU Orionis from ~ 4000 angstrom, to ~ 8 micron. Using an estimated visual extinction Av~1.5, a steady disk model with a central star mass ~0.3 Msun and a mass accretion rate ~ 2e-4 Msun/yr, we can reproduce the spectral energy distribution of FU Ori quite well. With the mid-infrared spectrum obtained by the Infrared Spectrograph (IRS) on board the Spitzer Space Telescope, we estimate that the outer radius of the hot, rapidly accreting inner disk is ~ 1 AU using disk models truncated at this outer radius. Inclusion of radiation from a cooler irradiated outer disk might reduce the outer limit of the hot inner disk to ~ 0.5 AU. In either case, the radius is inconsistent with a pure thermal instability model for the outburst. Our radiative transfer model implies that the central disk temperature Tc > 1000 K out to ~ 0.5 - 1 AU, suggesting that the magnetorotational instability (MRI) can be supported out to that distance. Assuming that the ~ 100 yr decay timescale in brightness of FU Ori represents the viscous timescale of the hot inner disk, we estimate the viscosity parameter (alpha) to be ~ 0.2 - 0.02 in the outburst state, consistent with numerical simulations of MRI in disks. The radial extent of the high mass accretion region is inconsistent with the model of Bell & Lin, but may be consistent with theories incorporating both gravitational instability and MRI.Comment: 32 pages, 10 figures, to appear in the Astrophysical Journa

First visual orbit for the prototypical colliding-wind binary WR 140

Wolf-Rayet stars represent one of the final stages of massive stellar evolution. Relatively little is known about this short-lived phase and we currently lack reliable mass, distance, and binarity determinations for a representative sample. Here we report the first visual orbit for WR 140(=HD193793), a WC7+O5 binary system known for its periodic dust production episodes triggered by intense colliding winds near periastron passage. The IOTA and CHARA interferometers resolved the pair of stars in each year from 2003--2009, covering most of the highly-eccentric, 7.9 year orbit. Combining our results with the recent improved double-line spectroscopic orbit of Fahed et al. (2011), we find the WR 140 system is located at a distance of 1.67 +/- 0.03 kpc, composed of a WR star with M_WR = 14.9 +/- 0.5 Msun and an O star with M_O = 35.9 +/- 1.3 Msun. Our precision orbit yields key parameters with uncertainties times 6 smaller than previous work and paves the way for detailed modeling of the system. Our newly measured flux ratios at the near-infrared H and Ks bands allow an SED decomposition and analysis of the component evolutionary states.Comment: Complete OIFITS dataset included via Data Conservancy Projec