No Arabic abstract
Context. Planets are thought to eventually form from the mostly gaseous (~99% of the mass) disks around young stars. The density structure and chemical composition of protoplanetary disks are affected by the incident radiation field at optical, FUV, and X-ray wavelengths, as well as by the dust properties. Aims. The effect of FUV and X-rays on the disk structure and the gas chemical composition are investigated. This work forms the basis of a second paper, which discusses the impact on diagnostic lines of, e.g., C+, O, H2O, and Ne+ observed with facilities such as Spitzer and Herschel. Methods. A grid of 240 models is computed in which the X-ray and FUV luminosity, minimum grain size, dust size distribution, and surface density distribution are varied in a systematic way. The hydrostatic structure and the thermo-chemical structure are calculated using ProDiMo. Results. The abundance structure of neutral oxygen is stable to changes in the X-ray and FUV luminosity, and the emission lines will thus be useful tracers of the disk mass and temperature. The C+ abundance distribution is sensitive to both X-rays and FUV. The radial column density profile shows two peaks, one at the inner rim and a second one at a radius r=5-10 AU. Ne+ and other heavy elements have a very strong response to X-rays, and the column density in the inner disk increases by two orders of magnitude from the lowest (LX = 1e29 erg/s) to the highest considered X-ray flux (LX = 1e32 erg/s). FUV confines the Ne+ ionized region to areas closer to the star at low X-ray luminosities (LX = 1e29 erg/s). H2O abundances are enhanced by X-rays due to higher temperatures in the inner disk and higher ionization fractions in the outer disk. The line fluxes and profiles are affected by the effects on these species, thus providing diagnostic value in the study of FUV and X-ray irradiated disks around T Tauri stars. (abridged)
Most of the mass in protoplanetary disks is in the form of gas. The study of the gas and its diagnostics is of fundamental importance in order to achieve a detailed description of the thermal and chemical structure of the disk. The radiation from the central star (from optical to X-ray wavelengths) and viscous accretion are the main source of energy and dominates the disk physics and chemistry in its early stages. This is the environment in which the first phases of planet formation will proceed. We investigate how stellar and disk parameters impact the fine-structure cooling lines [NeII], [ArII], [OI], [CII] and H2O rotational lines in the disk. These lines are potentially powerful diagnostics of the disk structure and their modelling permits a thorough interpretation of the observations carried out with instrumental facilities such as Spitzer and Herschel. Following Aresu et al. (2011), we computed a grid of 240 disk models, in which the X-ray luminosity, UV-excess luminosity, minimum dust grain size, dust size distribution power law and surface density distribution power law, are systematically varied. We solve self-consistently for the disk vertical hydrostatic structure in every model and apply detailed line radiative transfer to calculate line fluxes and profiles for a series of well known mid- and far-infrared cooling lines. The [OI] 63 micron line flux increases with increasing FUV luminosity when Lx < 1e30 erg/s, and with increasing X-ray luminosity when LX > 1e30 erg/s. [CII] 157 micron is mainly driven by FUV luminosity via C+ production, X-rays affect the line flux to a lesser extent. [NeII] 12.8 micron correlates with X-rays; the line profile emitted from the disk atmosphere shows a double-peaked component, caused by emission in the static disk atmosphere, next to a high velocity double-peaked component, caused by emission in the very inner rim. (abridged)
There is growing theoretical and observational evidence that protoplanetary disc evolution may be significantly affected by the canonical levels of far ultraviolet (FUV) radiation found in a star forming environment, leading to substantial stripping of material from the disc outer edge even in the absence of nearby massive stars. In this paper we perform the first full radiation hydrodynamic simulations of the flow from the outer rim of protoplanetary discs externally irradiated by such intermediate strength FUV fields, including direct modelling of the photon dominated region (PDR) which is required to accurately compute the thermal properties. We find excellent agreement between our models and the semi-analytic models of Facchini et al. (2016) for the profile of the flow itself, as well as the mass loss rate and location of their critical radius. This both validates their results (which differed significantly from prior semi-analytic estimates) and our new numerical method, the latter of which can now be applied to elements of the problem that the semi--analytic approaches are incapable of modelling. We also obtain the composition of the flow, but given the simple geometry of our models we can only hint at some diagnostics for future observations of externally irradiated discs at this stage. We also discuss the potential for these models as benchmarks for future photochemical-dynamical codes.
We present an open access grid of 3930 calculations of externally evaporating protoplanetary discs. This spans a range of disc sizes (1-400AU), disc masses, UV field strengths (10-10$^4$G$_0$) and stellar masses (0.05-1.9M$_odot$). The grid is publicly available for download, and offers a means of cheaply including external photoevaporation in disc evolutionary calculations. It can also be queried using an online tool for quick estimates of instantaneous mass loss rates (e.g for convenient evaluation of real observed systems). The `FRIED grid itself illustrates that for discs around stars $leq0.3$M$_odot$ external photoevaporation is effective down to small radii ($<50$AU) down to UV fields at least as weak as 10G$_0$. At the other end of the scale, in a $10^4$G$_0$ environment photoevaporation is effective down to 1AU even for stellar masses at least as high as 1.9M$_odot$. We also illustrate in which regimes CO survives in the photoevaporative outflow for significant mass loss rates; marking a system a good candidate to detect external photoevaporation in weak-intermediate UV environments through sub-Keplerian rotation. Finally we make illustrative mass loss rate estimates for discs in Taurus based on the Guilloteau et al. (2011) star-disc parameters, finding that around half are expected to have both significant mass loss and retain CO in the photoevaporative outflow.
We present new models for the X-ray photoevaporation of circumstellar discs which suggest that the resulting mass loss (occurring mainly over the radial range 10-40 AU) may be the dominant dispersal mechanism for gas around low mass pre-main sequence stars, contrary to the conclusions of previous workers. Our models combine use of the MOCASSIN Monte Carlo radiative transfer code and a self-consistent solution of the hydrostatic structure of the irradiated disc. We estimate the resulting photoevaporation rates assuming sonic outflow at the surface where the gas temperature equals the local escape temperature and derive mass loss rates of ~10^{-9} M_sun/yr, typically a factor 2-10 times lower than the corresponding rates in our previous work (Ercolano et al., 2008) where we did not adjust the density structure of the irradiated disc. The somewhat lower rates, and the fact that mass loss is concentrated towards slightly smaller radii, result from the puffing up of the heated disc at a few AU which partially screens the disc at tens of AU. (.....) We highlight the fact that X-ray photoevaporation has two generic advantages for disc dispersal compared with photoevaporation by extreme ultraviolet (EUV) photons that are only modestly beyond the Lyman limit: the demonstrably large X-ray fluxes of young stars even after they have lost their discs and the fact that X-rays are effective at penetrating much larger columns of material close to the star (abridged).
We study atomic line diagnostics of the inner regions of protoplanetary disks with our model of X-ray irradiated disk atmospheres which was previously used to predict observable levels of the NeII and NeIII fine-structure transitions at 12.81 and 15.55mum. We extend the X-ray ionization theory to sulfur and calculate the fraction of sulfur in S, S+, S2+ and sulfur molecules. For the DAlessio generic T Tauri star disk, we find that the SI fine-structure line at 25.55mum is below the detection level of the Spitzer Infrared Spectrometer (IRS), in large part due to X-ray ionization of atomic S at the top of the atmosphere and to its incorporation into molecules close to the mid-plane. We predict that observable fluxes of the SII 6718/6732AA forbidden transitions are produced in the upper atmosphere at somewhat shallower depths and smaller radii than the neon fine-structure lines. This and other forbidden line transitions, such as the OI 6300/6363AA and the CI 9826/9852AA lines, serve as complementary diagnostics of X-ray irradiated disk atmospheres. We have also analyzed the potential role of the low-excitation fine-structure lines of CI, CII, and OI, which should be observable by SOFIA and Herschel.