Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userFUV and X-ray irradiated protoplanetary disks: a grid of models II - Gas diagnostic line emission
298
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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)
rate research
Read More
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)
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.
High resolution spectra of typical wind diagnostics ([OI] 6300 r{A} and other forbidden emission lines) can often be decomposed into multiple components: high-velocity components with blueshifts up to several 100 km/s are usually attributed to fast jets, while narrow (NLVC) and broad (BLVC) low-velocity components are believed to trace slower disc winds. Under the assumption that the line-broadening is dominated by Keplerian rotation, several studies have found that the BLVCs should trace gas launched between 0.05 and 0.5 au and correlations between the properties of BLVCs and NLVCs have been interpreted as evidence for the emission tracing an extended MHD wind and not a photoevaporative wind. We calculated synthetic line profiles obtained from detailed photoionisation calculations of an X-ray photoevaporation model and a simple MHD wind model and analyzed the emission regions of different diagnostic lines and the resulting spectral profiles. The photoevaporation model reproduces most of the observed NLVCs but not the BLVCs or HVCs. The MHD model is able to reproduce all components but produces Keplerian double peaks at average inclinations that are rarely observed. The combination of MHD and photoevaporative winds could solve this problem. Our results suggest that the Gaussian decomposition does not allow for a clear distinction of flux from different wind regions and that the line broadening is often dominated by the velocity gradient in the outflow rather than by Keplerian rotation. We show that observed correlations between BLVC and NLVC do not necessarily imply a common origin in an extended MHD wind.
We study the PAH emission from protoplanetary disks. First, we discuss the dependence of the PAH band ratios on the hardness of the absorbed photons and the temperature of the stars. We show that the photon energy together with a varying degree of the PAH hydrogenation accounts for most of the observed PAH band ratios without the need to change the ionization degree of the molecules. We present an accurate treatment of stochastic heated grains in a vectorized three dimensional Monte Carlo dust radiative transfer code. The program is verified against results using ray tracing techniques. Disk models are presented for T Tauri and Herbig Ae stars. Particular attention is given to the photo-dissociation of the molecules. We consider beside PAH destruction also the survival of the molecules by vertical mixing within the disk. By applying typical X-ray luminosities the model accounts for the low PAH detection probability observed in T Tauri and the high PAH detection statistics found in Herbig Ae disks. Spherical halos above the disks are considered. We show that halos reduce the observed PAH band-to-continuum ratios when observed at high inclination. Finally, mid-IR images of disks around Herbig Ae disks are presented. We show that they are easier to resolve when PAH emission dominate.
ALMA observations of protoplanetary disks confirm earlier indications that there is a clear difference between the dust and gas radial extents. The origin of this difference is still debated, with both radial drift of the dust and optical depth effects suggested in the literature. In this work, the feedback of realistic dust particle distributions onto the gas chemistry and molecular emissivity is investigated, with a particular focus on CO isotopologues. The radial dust grain size distribution is determined using dust evolution models that include growth, fragmentation and radial drift. A new version of the code DALI is used to take into account how dust surface area and density influence the disk thermal structure, molecular abundances and excitation. The difference of dust and gas radial sizes is largely due to differences in the optical depth of CO lines and millimeter continuum, without the need to invoke radial drift. The effect of radial drift is primarily visible in the sharp outer edge of the continuum intensity profile. The gas outer radius probed by $^{12}$CO emission can easily differ by a factor of $sim 2$ between the models for a turbulent $alpha$ ranging between typical values. Grain growth and settling concur in thermally decoupling the gas and dust components, due to the low collision rate with large grains. As a result, the gas can be much colder than the dust at intermediate heights, reducing the CO excitation and emission, especially for low turbulence values. Also, due to disk mid-plane shadowing, a second CO thermal desorption (rather than photodesorption) front can occur in the warmer outer mid-plane disk. The models are compared to ALMA observations of HD 163296 as a test case. In order to reproduce the observed CO snowline of the system, a binding energy for CO typical of ice mixtures needs to be used rather than the lower pure CO value.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
