X-ray Diagnostics of Grain Depletion in Matter Accreting onto T Tauri Stars

Recent analysis of high resolution Chandra X-ray spectra has shown that the Ne/O abundance ratio is remarkably constant in stellar coronae. Based on this result, we point out the utility of the Ne/O ratio as a discriminant for accretion-related X-rays from T Tauri stars, and for probing the measure of grain-depletion of the accreting material in the inner disk. We apply the Ne/O diagnostic to the classical T Tauri stars BP Tau and TW Hya--the two stars found to date whose X-ray emission appears to originate, at least in part, from accretion activity. We show that TW Hya appears to be accreting material which is significantly depleted in O relative to Ne. In constrast, BP Tau has an Ne/O abundance ratio consistent with that observed for post-T Tauri stars. We interpret this result in terms of the different ages and evolutionary states of the circumstellar disks of these stars. In the young BP Tau disk (age 0.6 Myr) dust is still present near the disk corotation radius and can be ionized and accreted, re-releasing elements depleted onto grains. In the more evolved TW Hya disk (age 10 Myr), evidence points to ongoing coagulation of grains into much larger bodies, and possibly planets, that can resist the drag of inward-migrating gas, and accreting gas is consequently depleted of grain-forming elements.Comment: 13 pages, 1 Figure, ApJ Letters, in pres

Emission Line Variability of the Accreting Young Brown Dwarf 2MASSW J1207334-393254: From Hours to Years

We have obtained a series of high-resolution optical spectra for the brown dwarf 2MASSW J1207334-393254 (2M1207) using the ESO Very Large Telescope with the UVES spectrograph during two consecutive observing nights (time resolution of ~12 min) and the Magellan Clay telescope with the MIKE spectrograph. Combined with previously published results, these data allow us to investigate changes in the emission line spectrum of 2M1207 on timescales of hours to years. Most of the emission line profiles of 2M1207 are broad, in particular that of Halpha, indicating that the dominant fraction of the emission must be attributed to disk accretion rather than to magnetic activity. From the Halpha 10% width we deduce a relatively stable accretion rate between 10^(-10.1...-9.8) Msun/yr for two nights of consecutive observations. Therefore, either the accretion stream is nearly homogeneous over (sub-)stellar longitude or the system is seen face-on. Small but significant variations are evident throughout our near-continuous observation, and they reach a maximum after ~8 h, roughly the timescale on which maximum variability is expected across the rotation cycle. Together with past measurements, we confirm that the accretion rate of 2M1207 varies by more than one order of magnitude on timescales of months to years. Such variable mass accretion yields a plausible explanation for the observed spread in the accretion rate vs. mass diagram. The magnetic field required to drive the funnel flow is on the order of a few hundred G. Despite the obvious presence of a magnetic field, no radio nor X-ray emission has been reported for 2M1207. Possibly strong accretion suppresses magnetic activity in brown dwarfs, similar to the findings for higher mass T Tauri stars.Comment: accepted for publication in Ap

X-rays from accretion shocks in T Tauri stars: The case of BP Tau

We present an XMM-Newton observation of the classical T Tauri star BP Tau. In the XMM-Newton RGS spectrum the O {\sc vii} triplet is clearly detected with a very weak forbidden line indicating high plasma densities and/or a high UV flux environment. At the same time concurrent UV data point to a small hot spot filling factor suggesting an accretion funnel shock as the site of the X-ray and UV emission. Together with the X-ray data on TW Hya these new observations suggest such funnels to be a general feature in classical T Tauri stars.Comment: 4 pages, 4 figures, accepted by A&

On the origin of ionising photons emitted by T Tauri stars

We address the issue of the production of Lyman continuum photons by T Tauri stars, in an attempt to provide constraints on theoretical models of disc photoionisation. By treating the accretion shock as a hotspot on the stellar surface we show that Lyman continuum photons are produced at a rate approximately three orders of magnitude lower than that produced by a corresponding black body, and that a strong Lyman continuum is only emitted for high mass accretion rates. When our models are extended to include a column of material accreting on to the hotspot we find that the accretion column is extremely optically thick to Lyman continuum photons. Further, we find that radiative recombination of hydrogen atoms within the column is not an efficient means of producing photons with energies greater than 13.6eV, and find that an accretion column of any conceivable height suppresses the emission of Lyman continuum photons to a level below or comparable to that expected from the stellar photosphere. The photospheric Lyman continuum is itself much too weak to affect disc evolution significantly, and we find that the Lyman continuum emitted by an accretion shock is similarly unable to influence disc evolution significantly. This result has important consequences for models which use photoionisation as a mechanism to drive the dispersal of circumstellar discs, essentially proving that an additional source of Lyman continuum photons must exist if disc photoionisation is to be significant.Comment: 6 pages, 4 figures. Accepted for publication in MNRA

Profiles of Strong Permitted Lines in Classical T Tauri Stars

We present a spectral analysis of 30 T Tauri stars observed with the Hamilton echelle spectrograph over more than a decade. One goal is to test magnetospheric accretion model predictions. Observational evidence previously published supporting the model, such as emission line asymmetry and a high frequency of redshifted absorption components, are considered. We also discuss the relation between different line forming regions and search for good accretion rate indicators. In this work we confirm several important points of the models, such as the correlation between accretion and outflow, broad emission components that are mostly central or slightly blueshifted and only the occasional presence of redshifted absorption. We also show, however, that the broad emission components supposedly formed in the magnetospheric accretion flow only partially support the models. Unlike the predictions, they are sometimes redshifted, and are mostly found to be symmetric. The published theoretical profiles do not have a strong resemblance to our observed ones. We emphasize the need for accretion models to include a strong turbulent component before their profiles will match the observations. The effects of rotation, and the outflow components, will also be needed to complete the picture.Comment: 25 pages including 9 figures, 3 tables, accepted for publication in the Astronomical Journa

Characterizing the IYJ Excess Continuum Emission in T Tauri Stars

We present the first characterization of the excess continuum emission of accreting T Tauri stars between optical and near-infrared wavelengths. With nearly simultaneous spectra from 0.48 to 2.4 microns acquired with HIRES and NIRSPEC on Keck and SpeX on the IRTF, we find significant excess continuum emission throughout this region, including the I, Y, and J bands, which are usually thought to diagnose primarily photospheric emission. The IYJ excess correlates with the excess in the V band, attributed to accretion shocks in the photosphere, and the excess in the K band, attributed to dust in the inner disk near the dust sublimation radius, but it is too large to be an extension of the excess from these sources. The spectrum of the excess emission is broad and featureless, suggestive of blackbody radiation with a temperature between 2200 and 5000 K. The luminosity of the IYJ excess is comparable to the accretion luminosity inferred from modeling the blue and ultraviolet excess emission and may require reassessment of disk accretion rates. The source of the IYJ excess is unclear. In stars of low accretion rate, the size of the emitting region is consistent with cooler material surrounding small hot accretion spots in the photosphere. However, for stars with high accretion rates, the projected area is comparable to or exceeds that of the stellar surface. We suggest that at least some of the IYJ excess emission arises in the dust-free gas inside the dust sublimation radius in the disk.Comment: Accepted to ApJ, 31 pages, 21 figure

The 2008 extreme outburst of the young eruptive variable star Ex Lupi

In early 2008, the young low-mass star EX Lupi, the prototype of the EXor class of eruptive variables, optically brightened by over five magnitudes for a period of 7 months. The previous time a change of such amplitude had been observed in EX Lup was over 50 years ago. In this Letter we present new optical and near-IR high resolution spectroscopy of EX~Lup during the 2008 outburst. We investigate the physical characteristics of the outburst both soon after it began and some four months later, and consider the energetics and kinematics observed. Emission line strengths, widths, and profiles changed significantly between the two observations. Also, modeling of the 2.2935 um CO overtone bandhead emission suggests that an inner gap in the circumstellar gas disk around the star may be present and it is from the inner edge of the gas disk that the CO overtone emission probably arises. We derive a mass accretion luminosity and rate during the extreme outburst of ~2+-0.5~Lsun and ~2+-0.5x10^-7 Msun yr^-1, respectively, which suggests that this outburst was indeed one of the strongest witnessed in EX Lup, yet not as strong as those observed in FU Orionis stars.Comment: 5 pages, 4 figure

The spectral energy distribution of self-gravitating protostellar disks

The long wavelength emission of protostellar objects is commonly attributed to a disk of gas and dust around the central protostar. In the first stages of disk accretion or in the case of high mass protostars, the disk mass is likely to be sufficiently large, so that the disk self-gravity may have an impact on the dynamics and the emission properties of the disk. In this paper we describe the spectral energy distribution (SED) produced by a simple, non-flaring, self-gravitating accretion disk model. Self-gravity is included in the calculation of the rotation curve of the disk and in the energy balance equation, as a term of effective heating related to Jeans instability. In order to quantify in detail the requirements on the mass of the disk and on the accretion rate posed on the models by realistic situations, we compare the SEDs produced by these models with the observed SEDs of a small sample of well-studied protostellar objects. We find that relatively modest disks - even lighter than the central star - can lead to an interesting fit to the infrared SED of the FU Orionis objects considered, while in the case of T Tauri stars the required parameters fall outside the range suggested as acceptable by the general theoretical and observational scenario. On the basis of the present results, we may conclude that the contribution of a self-gravitating disk is plausible in several cases (in particular, for FU Orionis objects) and that, in the standard irradiation dominated disk scenario, it would help softening the requirements encountered by Keplerian accretion models.Comment: 26 pages, 7 figures, accepted by A&

HST measures of Mass Accretion Rates in the Orion Nebula Cluster

The present observational understanding of the evolution of the mass accretion rates (Macc) in pre-main sequence stars is limited by the lack of accurate measurements of Macc over homogeneous and large statistical samples of young stars. Such observational effort is needed to properly constrain the theory of star formation and disk evolution. Based on HST/WFPC2 observations, we present a study of Macc for a sample of \sim 700 sources in the Orion Nebula Cluster, ranging from the Hydrogen-burning limit to M\ast \sim 2M\odot. We derive Macc from both the U-band excess and the H{\alpha} luminosity (LH{\alpha}), after determining empirically both the shape of the typical accretion spectrum across the Balmer jump and the relation between the accretion luminosity (Lacc) and LH{\alpha}, that is Lacc/L\odot = (1.31\pm0.03)\cdotLH{\alpha}/L\odot + (2.63\pm 0.13). Given our large statistical sample, we are able to accurately investigate relations between Macc and the parameters of the central star such as mass and age. We clearly find Macc to increase with stellar mass, and decrease over evolutionary time, but we also find strong evidence that the decay of Macc with stellar age occurs over longer timescales for more massive PMS stars. Our best fit relation between these parameters is given by: log(Macc/M\odot\cdotyr)=(-5.12 \pm 0.86) -(0.46 \pm 0.13) \cdot log(t/yr) -(5.75 \pm 1.47)\cdot log(M\ast/M\odot) + (1.17 \pm 0.23)\cdot log(t/yr) \cdot log(M\ast/M\odot). These results also suggest that the similarity solution model could be revised for sources with M\ast > 0.5M\odot. Finally, we do not find a clear trend indicating environmental effects on the accretion properties of the sources.Comment: 17 pages, 15 figures, accepted for publication in Ap

The global structure of thin, stratified "alpha"-discs and the reliability of the one layer approximation

We report the results of a systematic comparison between the vertically averaged model and the vertically explicit model of steady state, Keplerian, optically thick "alpha"-discs. The simulations have concerned discs currently found in three different systems: dwarf novae, young stellar objects and active galactic nuclei. In each case, we have explored four decades of accretion rates and almost the whole disc area (except the narrow region where the vertically averaged model has degenerate solutions). We find that the one layer approach gives a remarkably good estimate of the main physical quantities in the disc, and specially the temperature at the equatorial plane which is accurate to within 30% for cases considered. The major deviations (by a factor < 4) are observed on the disc half-thickness. The sensitivity of the results to the "alpha"-parameter value has been tested for 0.001 < alpha < 0.1 and appears to be weak. This study suggests that the ``precision'' of the vertically averaged model which is easy to implement should be sufficient in practice for many astrophysical applications.Comment: 4 pages, PostScript. Accepted in Astronomy & Astrophysic