No Arabic abstract
Aims. We want to understand the chemistry and physics of disks on the basis of a large unbiased and statistically relevant grid of disk models. One of the main goals is to explore the diagnostic power of various gas emission lines and line ratios for deriving main disk parameters such as the gas mass. Methods. We explore the results of the DENT grid (Disk Evolution with Neat Theory) that consists of 300 000 disk models with 11 free parameters. Through a statistical analysis, we search for correlations and trends in an effort to find tools for disk diagnostic. Results. All calculated quantities like species masses, temperatures, continuum and line fluxes differ by several orders of magnitude across the entire parameter space. The broad distribution of these quantities as a function of input parameters shows the limitation of using a prototype T Tauri or Herbig Ae/Be disk model. The statistical analysis of the DENT grid shows that CO gas is rarely the dominant carbon reservoir in disks. Models with large inner radii (10 times the dust condensation radius) and/or shallow surface density gradients lack massive gas phase water reservoirs. Also, 60% of the disks have gas temperatures averaged over the oxygen mass in the range between 15 and 70 K; the average gas temperatures for CO and O differ by less than a factor two. Studying the observational diagnostics, the [CII] 158 mum fine structure line flux is very sensitive to the stellar UV flux and presence of a UV excess and it traces the outer disk radius (Rout). In the submm, the CO low J rotational lines also trace Rout. Low [OI] 63/145 line ratios (< a few) can be explained with cool atomic O gas in the uppermost surface layers leading to self-absorption in the 63 mum line; this occurs mostly for massive non-flaring, settled disk models without UV excess. ... abbreviated
Methods. We use the recently developed disk code ProDiMo to calculate the physico-chemical structure of protoplanetary disks and apply the Monte-Carlo line radiative transfer code RATRAN to predict observable line profiles and fluxes. We consider a series of Herbig Ae type disk models ranging from 10^-6 M_Sun to 2.2 10^-2 M_Sun (between 0.5 and 700 AU) to discuss the dependency of the line fluxes and ratios on disk mass for otherwise fixed disk parameters. Results. We find the [CII] 157.7 mum line to originate in LTE from the surface layers of the disk, where Tg > Td . The total emission is dominated by surface area and hence depends strongly on disk outer radius. The [OI] lines can be very bright (> 10^-16 W/m^2) and form in slightly deeper and closer regions under non-LTE conditions. The high-excitation [OI] 145.5 mum line, which has a larger critical density, decreases more rapidly with disk mass than the 63.2 mum line. Therefore, the [OI] 63.2 mum/145.5 mum ratio is a promising disk mass indicator, especially as it is independent of disk outer radius for Rout > 200 AU. CO is abundant only in deeper layers A_V >~ 0.05. For too low disk masses (M_disk <~10^-4 M_Sun) the dust starts to become transparent, and CO is almost completely photo-dissociated. For masses larger than that the lines are an excellent independent tracer of disk outer radius and can break the outer radius degeneracy in the [OI] 63.2 mum/[CII]157.7 mum line ratio. Conclusions. The far-IR fine-structure lines of [CII] and [OI] observable with Herschel provide a promising tool to measure the disk gas mass, although they are mainly generated in the atomic surface layers. In spatially unresolved observations, none of these lines carry much information about the inner, possibly hot regions < 30 AU.
We present self-consistent models of gas in optically-thick dusty disks and calculate its thermal, density and chemical structure. The models focus on an accurate treatment of the upper layers where line emission originates, and at radii $gtrsim 0.7$ AU. We present results of disks around $sim 1{rm M}_{odot}$ stars where we have varied dust properties, X-ray luminosities and UV luminosities. We separately treat gas and dust thermal balance, and calculate line luminosities at infrared and sub-millimeter wavelengths from all transitions originating in the predominantly neutral gas that lies below the ionized surface of the disk. We find that the [ArII] 7$mu$m, [NeII] 12.8$mu$m, [FeI] 24$mu$m, [SI] 25$mu$m, [FeII] 26$mu$m, [SiII] 35 $mu$m, [OI] 63$mu$m and pure rotational lines of H$_2$, H$_2$O and CO can be quite strong and are good indicators of the presence and distribution of gas in disks. We apply our models to the disk around the nearby young star, TW Hya, and find good agreement between our model calculations and observations. We also predict strong emission lines from the TW Hya disk that are likely to be detected by future facilities. A comparison of CO observations with our models suggests that the gas disk around TW Hya may be truncated to $sim 120 $ AU, compared to its dust disk of 174 AU. We speculate that photoevaporation due to the strong stellar FUV field from TW Hya is responsible for the gas disk truncation.
We calculate the emission of protoplanetary disks threaded by a poloidal magnetic field and irradiated by the central star. The radial structure of these disks was studied by Shu and collaborators and the vertical structure was studied by Lizano and collaborators. We consider disks around low mass protostars, T Tauri stars, and FU Ori stars with different mass-to-flux ratios $lambda_{rm sys}$. We calculate the spectral energy distribution and the antenna temperature profiles at 1 mm and 7 mm convolved with the ALMA and VLA beams. We find that disks with weaker magnetization (high values of $lambda_{rm sys}$) emit more than disks with stronger magnetization (low values of $lambda_{rm sys}$). This happens because the former are denser, hotter and have larger aspect ratios, receiving more irradiation from the central star. The level of magnetization also affects the optical depth at millimeter wavelengths, being larger for disks with high $lambda_{rm sys}$. In general, disks around low mass protostars and T Tauri stars are optically thin at 7 mm while disks around FU Ori are optically thick. A qualitative comparison of the emission of these magnetized disks, including heating by an external envelope, with the observed millimeter antenna temperature profiles of HL Tau indicates that large cm grains are required to increase the optical depth and reproduce the observed 7 mm emission at large radii.
We discuss accretion and outflow properties of three very low-mass young stellar objects based on broad-band mid-resolution X-Shooter/VLT spectra. Our targets (FU TauA, 2M 1207-39, and Par-Lup3-4) have spectral types between M5 and M8, ages between 1 Myr and ~10 Myr, and are known to be accreting from previous studies. The final objective of our project is the determination of mass outflow to accretion rate for objects near or within the substellar regime as a probe for the T Tauri phase of brown dwarfs and the investigation of variability in the accretion and outflow processes.
Extreme ultraviolet (EUV, 13.6 eV $< h u lta 100$ eV) and X-rays in the 0.1-2 keV band can heat the surfaces of disks around young, low mass stars to thousands of degrees and ionize species with ionization potentials greater than 13.6 eV. Shocks generated by protostellar winds can also heat and ionize the same species close to the star/disk system. These processes produce diagnostic lines (e.g., [NeII] 12.8 $mu$m and [OI] 6300 AA) that we model as functions of key parameters such as EUV luminosity and spectral shape, X-ray luminosity and spectral shape, and wind mass loss rate and shock speed. Comparing our models with observations, we conclude that either internal shocks in the winds or X-rays incident on the disk surfaces often produce the observed [NeII] line, although there are cases where EUV may dominate. Shocks created by the oblique interaction of winds with disks are unlikely [NeII] sources because these shocks are too weak to ionize Ne. Even if [NeII] is mainly produced by X-rays or internal wind shocks, the neon observations typically place upper limits of $lta 10^{42}$ s$^{-1}$ on the EUV photon luminosity of these young low mass stars. The observed [OI] 6300 AA line has both a low velocity component (LVC) and a high velocity component. The latter likely arises in internal wind shocks. For the former we find that X-rays likely produce more [OI] luminosity than either the EUV layer, the transition layer between the EUV and X-ray layer, or the shear layer where the protostellar wind shocks and entrains disk material in a radial flow across the surface of the disk. Our soft X-ray models produce [OI] LVCs with luminosities up to $10^{-4}$ L$_odot$, but may not be able to explain the most luminous LVCs.