The effect of scattering on the structure and SED of protoplanetary disks

In this paper we investigate how the inclusion of scattering of the stellar radiation into a passive flaring disk model affects its structure and spectral energy distribution, and whether neglecting it could significantly decrease the model reliability. In order to address these questions we construct a detailed 1+1D vertical structure model in which the scattering properties of the dust can be varied. Models are presented with and without dust scattering, and for different albedos and phase functions. It is found that scattering has the effect of reducing the disk temperature at all heights, so that the disk "shrinks", i.e., the the density at all intermediate heights decreases. However, this effect in most cases is more than compensated by the increase of the total extinction (absorption + scattering) cross section, so that the surface scale height increases, and images in scattered light will see a slightly thicker disk. The integrated infrared emission decreases as the albedo increases, because an increasing part of the flux captured by the disk is reflected away instead of absorbed and reprocessed. The reduction of the infrared thermal emission of the disk is stronger at short wavelengths (near infrared) and practically negligible at millimeter wavelengths. For relatively low albedo (alb <~ 0.5), or for strongly forward-peaked scattering (g roughly >0.8), the infrared flux reduction is relatively small.Comment: Accepted for publication in Astronomy & Astrophysic

On the Evolution and Survival of Protoplanets Embedded in a Protoplanetary Disk

We model the evolution of a Jupiter-mass protoplanet formed by the disk instability mechanism at various radial distances accounting for the presence of the disk. Using three different disk models, it is found that a newly-formed Jupiter-mass protoplanet at radial distance of $\lesssim$ 5-10 AU cannot undergo a dynamical collapse and evolve further to become a gravitational bound planet. We therefore conclude that {\it giant planets, if formed by the gravitational instability mechanism, must form and remain at large radial distances during the first $\sim$ 10$^5-10^6$ years of their evolution}. The minimum radial distances in which protoplanets of 1 Saturn-mass, 3 and 5 Jupiter-mass protoplanets can evolve using a disk model with $\dot{M}=10^{-6} M_{Sun}/yr$ and $\alpha=10^{-2}$ are found to be 12, 9, and 7 AU, respectively. The effect of gas accretion on the planetary evolution of a Jupiter-mass protoplanet is also investigated. It is shown that gas accretion can shorten the pre-collapse timescale substantially. Our study suggests that the timescale of the pre-collapse stage does not only depend on the planetary mass, but is greatly affected by the presence of the disk and efficient gas accretion.Comment: 26 pages, 2 tables, 10 figures. Accepted for publication in Ap

A Spatially Resolved Inner Hole in the Disk around GM Aurigae

We present 0.3 arcsec resolution observations of the disk around GM Aurigae with the Submillimeter Array (SMA) at a wavelength of 860 um and with the Plateau de Bure Interferometer at a wavelength of 1.3 mm. These observations probe the distribution of disk material on spatial scales commensurate with the size of the inner hole predicted by models of the spectral energy distribution. The data clearly indicate a sharp decrease in millimeter optical depth at the disk center, consistent with a deficit of material at distances less than ~20 AU from the star. We refine the accretion disk model of Calvet et al. (2005) based on the unresolved spectral energy distribution (SED) and demonstrate that it reproduces well the spatially resolved millimeter continuum data at both available wavelengths. We also present complementary SMA observations of CO J=3-2 and J=2-1 emission from the disk at 2" resolution. The observed CO morphology is consistent with the continuum model prediction, with two significant deviations: (1) the emission displays a larger CO J=3-2/J=2-1 line ratio than predicted, which may indicate additional heating of gas in the upper disk layers; and (2) the position angle of the kinematic rotation pattern differs by 11 +/- 2 degrees from that measured at smaller scales from the dust continuum, which may indicate the presence of a warp. We note that photoevaporation, grain growth, and binarity are unlikely mechanisms for inducing the observed sharp decrease in opacity or surface density at the disk center. The inner hole plausibly results from the dynamical influence of a planet on the disk material. Warping induced by a planet could also potentially explain the difference in position angle between the continuum and CO data sets.Comment: 12 pages, 6 figures, accepted for publication in Ap

The Effects of UV Continuum and Lyman alpha Radiation on the Chemical Equilibrium of T Tauri Disks

We show in this Letter that the spectral details of the FUV radiation fields have a large impact on the chemistry of protoplanetary disks surrounding T Tauri stars. We show that the strength of a realistic stellar FUV field is significantly lower than typically assumed in chemical calculations and that the radiation field is dominated by strong line emission, most notably Lyman alpha radiation. The effects of the strong Lyman alpha emission on the chemical equilibrium in protoplanetary disks has previously been unrecognized. We discuss the impact of this radiation on molecular observations in the context of a radiative transfer model that includes both direct attenuation and scattering. In particular, Lyman alpha radiation will directly dissociate water vapor and may contribute to the observed enhancements of CN/HCN in disks.Comment: 14 pages, 4 figures, accepted by ApJ Letter

Vertical structure models of T Tauri and Herbig Ae/Be disks

In this paper we present detailed models of the vertical structure (temperature and density) of passive irradiated circumstellar disks around T Tauri and Herbig Ae/Be stars. In contrast to earlier work, we use full frequency- and angle-dependent radiative transfer instead of the usual moment equations. We find that this improvement of the radiative transfer has strong influence on the resulting vertical structure of the disk, with differences in temperature as large as 70 %. However, the spectral energy distribution (SED) is only mildly affected by this change. In fact, the SED compares reasonably well with that of improved versions of the Chiang & Goldreich (CG) model. This shows that the latter is a reasonable model for the SED, in spite of its simplicity. It also shows that from the SED alone, little can be learned about the vertical structure of a passive circumstellar disk. The molecular line emission from these disks is more sensitive to the vertical temperature and density structure, and we show as an example how the intensity and profiles of various CO lines depend on the adopted disk model. The models presented in this paper can also serve as the basis of theoretical studies of e.g. dust coagulation and settling in disks.Comment: 12 pages, 15 figures, accepted for publication in A&

Empirical Constraints on Turbulence in Protoplanetary Accretion Disks

We present arcsecond-scale Submillimeter Array observations of the CO(3-2) line emission from the disks around the young stars HD 163296 and TW Hya at a spectral resolution of 44 m/s. These observations probe below the ~100 m/s turbulent linewidth inferred from lower-resolution observations, and allow us to place constraints on the turbulent linewidth in the disk atmospheres. We reproduce the observed CO(3-2) emission using two physical models of disk structure: (1) a power-law temperature distribution with a tapered density distribution following a simple functional form for an evolving accretion disk, and (2) the radiative transfer models developed by D'Alessio et al. that can reproduce the dust emission probed by the spectral energy distribution. Both types of models yield a low upper limit on the turbulent linewidth (Doppler b-parameter) in the TW Hya system (<40 m/s), and a tentative (3-sigma) detection of a ~300 m/s turbulent linewidth in the upper layers of the HD 163296 disk. These correspond to roughly <10% and 40% of the sound speed at size scales commensurate with the resolution of the data. The derived linewidths imply a turbulent viscosity coefficient, alpha, of order 0.01 and provide observational support for theoretical predictions of subsonic turbulence in protoplanetary accretion disks.Comment: 18 pages, 9 figures, accepted for publication in Ap

Millimeter imaging of HD 163296: probing the disk structure and kinematics

We present new multi-wavelength millimeter interferometric observations of the Herbig Ae star HD 163296 obtained with the IRAM/PBI, SMA and VLA arrays both in continuum and in the 12CO, 13CO and C18O emission lines. Gas and dust properties have been obtained comparing the observations with self-consistent disk models for the dust and CO emission. The circumstellar disk is resolved both in the continuum and in CO. We find strong evidence that the circumstellar material is in Keplerian rotation around a central star of 2.6 Msun. The disk inclination with respect to the line of sight is 46+-4 deg with a position angle of 128+-4 deg. The slope of the dust opacity measured between 0.87 and 7 mm (beta=1) confirms the presence of mm/cm-size grains in the disk midplane. The dust continuum emission is asymmetric and confined inside a radius of 200 AU while the CO emission extends up to 540 AU. The comparison between dust and CO temperature indicates that CO is present only in the disk interior. Finally, we obtain an increasing depletion of CO isotopomers from 12CO to 13CO and C18O. We argue that these results support the idea that the disk of HD 163296 is strongly evolved. In particular, we suggest that there is a strong depletion of dust relative to gas outside 200 AU; this may be due to the inward migration of large bodies that form in the outer disk or to clearing of a large gap in the dust distribution by a low mass companion.Comment: Accepted for publication on A&A, 16 page

A Spitzer IRS Survey of NGC 1333: Insights into disk evolution from a very young cluster

We report on the {\lambda} = 5-36{\mu}m Spitzer Infrared Spectrograph spectra of 79 young stellar objects in the very young nearby cluster NGC 1333. NGC 1333's youth enables the study of early protoplanetary disk properties, such as the degree of settling as well as the formation of gaps and clearings. We construct spectral energy distributions (SEDs) using our IRS data as well as published photometry and classify our sample into SED classes. Using "extinction-free" spectral indices, we determine whether the disk, envelope, or photosphere dominates the spectrum. We analyze the dereddened spectra of objects which show disk dominated emission using spectral indices and properties of silicate features in order to study the vertical and radial structure of protoplanetary disks in NGC 1333. At least nine objects in our sample of NGC 1333 show signs of large (several AU) radial gaps or clearings in their inner disk. Disks with radial gaps in NGC 1333 show more-nearly pristine silicate dust than their radially continuous counterparts. We compare properties of disks in NGC 1333 to those in three other well studied regions, Taurus-Auriga, Ophiuchus and Chamaeleon I, and find no difference in their degree of sedimentation and dust processing.Comment: 67 pages, 20 figures, accepted to The Astrophysical Journal Supplement Serie

Spitzer observations of the Orion OB1 association: disk census in the low mass stars

We present new Spitzer Space Telescope observations of two fields in the Orion OB1 association. We report here IRAC/MIPS observations for 115 confirmed members and 41 photometric candidates of the ~10 Myr 25 Orionis aggregate in the OB1a subassociation, and 106 confirmed members and 65 photometric candidates of the 5 Myr region located in the OB1b subassociation. The 25 Orionis aggregate shows a disk frequency of 6% while the field in the OB1b subassociation shows a disk frequency of 13%. Combining IRAC, MIPS and 2MASS photometry we place stars bearing disks in several classes: stars with optically thick disks (class II systems), stars with an inner transitional disks (transitional disk candidates) and stars with "evolved disks"; the last exhibit smaller IRAC/MIPS excesses than class II systems. In all, we identify 1 transitional disk candidate in the 25 Orionis aggregate and 3 in the OB1b field; this represents ~10% of the disk bearing stars, indicating that the transitional disk phase can be relatively fast. We find that the frequency of disks is a function of the stellar mass, suggesting a maximum around stars with spectral type M0. Comparing the infrared excess in the IRAC bands among several stellar groups we find that inner disk emission decays with stellar age, showing a correlation with the respective disk frequencies. The disk emission at the IRAC and MIPS bands in several stellar groups indicates that disk dissipation takes place faster in the inner region of the disks. Comparison with models of irradiated accretion disks, computed with several degrees of settling, suggests that the decrease in the overall accretion rate observed in young stellar groups is not sufficient to explain the weak disk emission observed in the IRAC bands for disk bearing stars with ages 5 Myr or older.Comment: Accepted in the Astrophysical Journa