X-Shooter study of accretion in ρ\rho-Ophiucus: very low-mass stars and brown dwarfs

We present new VLT/X-Shooter optical and NIR spectra of a sample of 17 candidate young low-mass stars and BDs in the rho-Ophiucus cluster. We derived SpT and Av for all the targets, and then we determined their physical parameters. All the objects but one have M*<0.6 Msun, and 8 have mass below or close to the hydrogen-burning limit. Using the intensity of various emission lines present in their spectra, we determined the Lacc and Macc for all the objects. When compared with previous works targeting the same sample, we find that, in general, these objects are not as strongly accreting as previously reported, and we suggest that the reason is our more accurate estimate of the photospheric parameters. We also compare our findings with recent works in other slightly older star-forming regions to investigate possible differences in the accretion properties, but we find that the accretion properties for our targets have the same dependence on the stellar and substellar parameters as in the other regions. This leads us to conclude that we do not find evidence for a different dependence of Macc with M* when comparing low-mass stars and BDs. Moreover, we find a similar small (1 dex) scatter in the Macc-M* relation as in some of our recent works in other star-forming regions, and no significant differences in Macc due to different ages or properties of the regions. The latter result suffers, however, from low statistics and sample selection biases in the current studies. The small scatter in the Macc-M* correlation confirms that Macc in the literature based on uncertain photospheric parameters and single accretion indicators, such as the Ha width, can lead to a scatter that is unphysically large. Our studies show that only broadband spectroscopic surveys coupled with a detailed analysis of the photospheric and accretion properties allows us to properly study the evolution of disk accretion rates.Comment: accepted for publication in Astronomy & Astrophysics. Abstract shortened to fit arXiv constraint

The gas temperature in the surface layers of protoplanetary disks

Models for the structure of protoplanetary disks have so far been based on the assumption that the gas and the dust temperature are equal. The gas temperature, an essential ingredient in the equations of hydrostatic equilibrium of the disk, is then determined from a continuum radiative transfer calculation, in which the continuum opacity is provided by the dust. It has been long debated whether this assumption still holds in the surface layers of the disk, where the dust infrared emission features are produced. In this paper we compute the temperature of the gas in the surface layers of the disk in a self-consistent manner. The gas temperature is determined from a heating-cooling balance equation in which processes such as photoelectric heating, dissociative heating, dust-gas thermal heat exchange and line cooling are included. The abundances of the dominant cooling species such as CO, C, C+ and O are determined from a chemical network based on the atomic species H, He, C, O, S, Mg, Si, Fe (Kamp & Bertoldi 2000). The underlying disk models to our calculations are the models of Dullemond, van Zadelhoff & Natta (2002). We find that in general the dust and gas temperature are equal to withing 10% for A_V >~ 0.1, which is above the location of the `super-heated surface layer' in which the dust emission features are produced (e.g. Chiang & Goldreich 1997). High above the disk surface the gas temperature exceeds the dust temperature and can can become -- in the presence of polycyclic aromatic hydrocarbons -- as high as 600 K at a radius of 100 AU. This is a region where CO has fully dissociated, but a significant fraction of hydrogen is still in molecular form. The densities are still high enough for non-negligible H_2 emission to be produced.....(see paper for full abstract)Comment: 28 pages, 8 figures, accepted for publication in Ap

The Onset of Planet Formation in Brown Dwarf Disks

The onset of planet formation in protoplanetary disks is marked by the growth and crystallization of sub-micron-sized dust grains accompanied by dust settling toward the disk mid-plane. Here we present infrared spectra of disks around brown dwarfs and brown dwarf candidates. We show that all three processes occur in such cool disks in a way similar or identical to that in disks around low- and intermediate-mass stars. These results indicate that the onset of planet formation extends to disks around brown dwarfs, suggesting that planet formation is a robust process occurring in most young circumstellar disks.Comment: Published in Science 2005, vol 310, 834; 3 pages in final format, 4 figures + 8 pages Supporting Online Material. For final typeset, see http://www.sciencemag.org/cgi/content/abstract/310/5749/834?eto

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

An extensive VLT/X-Shooter library of photospheric templates of pre-main sequence stars

Studies of the formation and evolution of young stars and their disks rely on the knowledge of the stellar parameters of the young stars. The derivation of these parameters is commonly based on comparison with photospheric template spectra. Furthermore, chromospheric emission in young active stars impacts the measurement of mass accretion rates, a key quantity to study disk evolution. Here we derive stellar properties of low-mass pre-main sequence stars without disks, which represent ideal photospheric templates for studies of young stars. We also use these spectra to constrain the impact of chromospheric emission on the measurements of mass accretion rates. The spectra in reduced, flux-calibrated, and corrected for telluric absorption form are made available to the community. We derive the spectral type for our targets by analyzing the photospheric molecular features present in their VLT/X-Shooter spectra by means of spectral indices and comparison of the relative strength of photospheric absorption features. We also measure effective temperature, gravity, projected rotational velocity, and radial velocity from our spectra by fitting them with synthetic spectra with the ROTFIT tool. The targets have negligible extinction and spectral type from G5 to M8. We perform synthetic photometry on the spectra to derive the typical colors of young stars in different filters. We measure the luminosity of the emission lines present in the spectra and estimate the noise due to chromospheric emission in the measurements of accretion luminosity in accreting stars. We provide a calibration of the photospheric colors of young PMS stars as a function of their spectral type in a set of standard broad-band optical and near-infrared filters. For stars with masses of ~ 1.5Msun and ages of ~1-5 Myr, the chromospheric noise converts to a limit of measurable mass accretion rates of ~ 3x10^-10 Msun/yr.Comment: Accepted for publication on Astronomy & Astrophysics. The spectra of the photospheric templates will be uploaded to Vizier, but are already available on request. Abstract shortened for arxiv constraints. Language edited versio

On the gas content of transitional disks: a VLT/X-Shooter study of accretion and winds

Transitional disks (TDs) are thought to be a late evolutionary stage of protoplanetary disks with dust depleted inner regions. The mechanism responsible for this depletion is still under debate. To constrain the models it is mandatory to have a good understanding of the properties of the gas content of the inner disk. Using X-Shooter broad band -UV to NIR- medium resolution spectroscopy we derive the stellar, accretion, and wind properties of a sample of 22 TDs. The analysis of these properties allows us to put strong constraints on the gas content in a region very close to the star (<0.2 AU) which is not accessible with any other observational technique. We fit the spectra with a self-consistent procedure to derive simultaneously SpT,Av,and mass accretion rates (Macc) of the targets. From forbidden emission lines we derive the wind properties of the targets. Comparing our findings to values for cTTs, we find that Macc and wind properties of 80% of the TDs in our sample, which is strongly biased towards strongly accreting objects, are comparable to those of cTTs. Thus, there are (at least) some TDs with Macc compatible with those of cTTs, irrespective of the size of the dust inner hole.Only in 2 cases Macc are much lower, while the wind properties are similar. We do not see any strong trend of Macc with the size of the dust depleted cavity, nor with the presence of a dusty optically thick disk close to the star. In the TDs in our sample there is a gas rich inner disk with density similar to that of cTTs disks. At least for some TDs, the process responsible of the inner disk clearing should allow for a transfer of gas from the outer disk to the inner region. This should proceed at a rate that does not depend on the physical mechanism producing the gap seen in the dust emission and results in a gas density in the inner disk similar to that of unperturbed disks around stars of similar mass.Comment: Accepted on Astronomy & Astrophysics. Abstract shortened to fit arXiv constraint

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&

Trapping dust particles in the outer regions of protoplanetary disks

Aims. We attempt to explain grain growth to mm sized particles and their retention in the outer regions of protoplanetary disks, as observed at sub-mm and mm wavelengths, by investigating whether strong inhomogeneities in the gas density profiles can decelerate excessive radial drift and help the dust particles to grow. Methods. We use coagulation/fragmentation and disk-structure models, to simulate the evolution of dust in a bumpy surface density profile, which we mimic with a sinusoidal disturbance. For different values of the amplitude and length scale of the bumps, we investigate the ability of this model to produce and retain large particles on million-year timescales. In addition, we compare the pressure inhomogeneities considered in this work with the pressure profiles that come from magnetorotational instability. Using the Common Astronomy Software Applications ALMA simulator, we study whether there are observational signatures of these pressure inhomogeneities that can be seen with ALMA. Results. We present the conditions required to trap dust particles and the corresponding calculations predicting the spectral slope in the mm-wavelength range, to compare with current observations. Finally, we present simulated images using different antenna configurations of ALMA at different frequencies, to show that the ring structures will be detectable at the distances of either the Taurus Auriga or Ophiucus star-forming regions

Testing the theory of grain growth and fragmentation by millimeter observations of protoplanetary disks

Context. Observations at sub-millimeter and mm wavelengths will in the near future be able to resolve the radial dependence of the mm spectral slope in circumstellar disks with a resolution of around a few AU at the distance of the closest star-forming regions. Aims. We aim to constrain physical models of grain growth and fragmentation by a large sample of (sub-)mm observations of disks around pre-main sequence stars in the Taurus-Auriga and Ophiuchus star-forming regions. Methods. State-of-the-art coagulation/fragmentation and disk-structure codes are coupled to produce steady-state grain size distributions and to predict the spectral slopes at (sub-)mm wavelengths. Results. This work presents the first calculations predicting the mm spectral slope based on a physical model of grain growth. Our models can quite naturally reproduce the observed mm-slopes, but a simultaneous match to the observed range of flux levels can only be reached by a reduction of the dust mass by a factor of a few up to about 30 while keeping the gas mass of the disk the same. This dust reduction can either be due to radial drift at a reduced rate or during an earlier evolutionary time (otherwise the predicted fluxes would become too low) or due to efficient conversion of dust into larger, unseen bodies.Comment: Accepted for publication in A&A Letters. 5 pages, 3 figure

Passive irradiated circumstellar disks with an inner hole

A model for irradiated dust disks around Herbig Ae stars is proposed. The model is based on the flaring disk model of Chiang & Goldreich (1997), but with the central regions of the disk removed. The inner rim of the disk is puffed up and is much hotter than the rest of the disk, because it is directly exposed to the stellar flux. If located at the dust evaporation radius, its reemitted flux produces a conspicuous bump in the SED which peaks at 2-3 micron. We propose that this emission is the explanation for the near-infrared bump observed in the SEDs of Herbig Ae stars. We study for which stellar parameters this bump would be observable, and find that it is the case for Herbig Ae stellar parameters but not for T-Tauri stars, confirming what is found from the observations. We also study the effects of the shadow cast by the inner rim over the rest of the flaring disk. The shadowed region can be quite large, and under some circumstances the entire disk may lie in the shadow. This shadowed region will be much cooler than an unshadowed flaring disk, since its only heating sources are radial radiative diffusion and possible indirect sources of irradiation. Under certain special circumstances the shadowing effect can suppress, or even completely eliminate, the 10 micron emission feature from the spectrum, which might explain the anomalous SEDs of some isolated Herbig Ae stars in the sample of Meeus et al. (2001). At much larger radii the disk emerges from the shadow, and continues as a flaring disk towards the outer edge. The complete model, including structure of the inner edge, shadowed region and the flared outer part, is described in detail in this paper, and we show examples of the general behavior of the model for varying parameters.Comment: Accepted for publication in the Astrophysical Journa