FUV and X-ray irradiated protoplanetary disks: a grid of models I. The disk structure

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)Comment: 47 pages, accepted by Astronomy and Astrophysics, a high resolution version of the paper is located at http://www.astro.rug.nl/~meijerink/disk_paperI_xrays.pd

[OI] disk emission in the Taurus star forming region

The structure of protoplanetary disks is thought to be linked to the temperature and chemistry of their dust and gas. Whether the disk is flat or flaring depends on the amount of radiation that it absorbs at a given radius, and on the efficiency with which this is converted into thermal energy. The understanding of these heating and cooling processes is crucial to provide a reliable disk structure for the interpretation of dust continuum emission and gas line fluxes. Especially in the upper layers of the disk, where gas and dust are thermally decoupled, the infrared line emission is strictly related to the gas heating/cooling processes. We aim to study the thermal properties of the disk in the oxygen line emission region, and to investigate the relative importance of X-ray (1-120 Angstrom) and far-UV radiation (FUV, 912-2070 Angstrom) for the heating balance there. We use [OI] 63 micron line fluxes observed in a sample of protoplanetary disks of the Taurus/Auriga star forming region and compare it to the model predictions presented in our previous work. The data were obtained with the PACS instrument on board the Herschel Space Observatory as part of the Herschel Open Time Key Program GASPS (GAS in Protoplanetary diskS), published in Howard et al. (2013). Our theoretical grid of disk models can reproduce the [OI] absolute fluxes and predict a correlation between [OI] and the sum Lx+Lfuv. The data show no correlation between the [OI] line flux and the X-ray luminosity, the FUV luminosity or their sum. The data show that the FUV or X-ray radiation has no notable impact on the region where the [OI] line is formed. This is in contrast with what is predicted from our models. Possible explanations are that the disks in Taurus are less flaring than the hydrostatic models predict, and/or that other disk structure aspects that were left unchanged in our models are important. ..abridged..Comment: 9 pages, accepted for publication in A&

Uncertainties in water chemistry in disks: An application to TW Hya

Context. This paper discusses the sensitivity of water lines to chemical processes and radiative transfer for the protoplanetary disk around TW Hya. The study focuses on the Herschel spectral range in the context of new line detections with the PACS instrument from the Gas in Protoplanetary Systems project (GASPS). Aims. The paper presents an overview of the chemistry in the main water reservoirs in the disk around TW Hya. It discusses the limitations in the interpretation of observed water line fluxes. Methods. ... (abbreviated) Results. We report new line detections of p-H2O (3_22-2_11) at 89.99 micron and CO J=18-17 at 144.78 micron for the disk around TW Hya. Disk modeling shows that the far-IR fine structure lines ([OI], [CII]) and molecular submm lines are very robust to uncertainties in the chemistry, while the water line fluxes can change by factors of a few. The water lines are optically thick, sub-thermally excited and can couple to the background continuum radiation field. The low-excitation water lines are also sensitive to uncertainties in the collision rates, e.g. with neutral hydrogen. The gas temperature plays an important role for the [OI] fine structure line fluxes, the water line fluxes originating from the inner disk as well as the high excitation CO, CH+ and OH lines. Conclusions. Due to their sensitivity on chemical input data and radiative transfer, water lines have to be used cautiously for understanding details of the disk structure. Water lines covering a wide range of excitation energies provide access to the various gas phase water reservoirs (inside and outside the snow line) in protoplanetary disks and thus provide important information on where gas-phase water is potentially located. Experimental and/or theoretical collision rates for H2O with atomic hydrogen are needed to diminish uncertainties from water line radiative transfer.Comment: accepted for publication in A&

Thermochemical modelling of brown dwarf discs

RCH acknowledges funding by the Austrian Science Fund (FWF): project number P24790.The physical properties of brown dwarf discs, in terms of their shapes and sizes, are still largely unexplored by observations. ALMA has by far the best capabilities to observe these discs in sub-mm CO lines and dust continuum, while also spatially resolving some discs. To what extent brown dwarf discs are similar to scaled-down T Tauri discs is currently unknown, and this work is a step towards establishing a relationship through the eventual modelling of future observations. We use observations of the brown dwarf disc ρ Oph 102 to infer a fiducial model around which we build a small grid of brown dwarf disc models, in order to model the CO, HCN, and HCO+ line fluxes and the chemistry which drives their abundances. These are the first brown dwarf models to be published which relate detailed, 2D radiation thermochemical disc models to observational data. We predict that moderately extended ALMA antenna configurations will spatially resolve CO line emission around brown dwarf discs, and that HCN and HCO+ will be detectable in integrated flux, following our conclusion that the flux ratios of these molecules to CO emission are comparable to that of T Tauri discs. These molecules have not yet been observed in sub-mm wavelengths in a brown dwarf disc, yet they are crucial tracers of the warm surface-layer gas and of ionization in the outer parts of the disc. We present the prediction that if the physical and chemical processes in brown dwarf discs are similar to those that occur in T Tauri discs-as our models suggest-then the same diagnostics that are used for T Tauri discs can be used for brown dwarf discs (such as HCN and HCO+ lines that have not yet been observed in the sub-mm), and that these lines should be observable with ALMA. Through future observations, either confirmation (or refutation) of these ideas about brown dwarf disc chemistry will have strong implications for our understanding of disc chemistry, structure, and subsequent planet formation in brown dwarf discs.Publisher PDFPeer reviewe

Mid-IR spectra of pre-main sequence Herbig stars : an explanation for the non-detections of water lines

The research leading to these results has received funding from the European Union Seventh Framework Programme FP7-2011 under grant agreement No. 284405.Context. The mid-IR detection rate of water lines in disks around Herbig stars disks is about 5%, while it is around 50% for disks around T Tauri stars. The reason for this is still unclear. Aims. In this study, we want to find an explanation for the different detection rates between low mass and high mass pre-main-sequence stars in the mid-IR regime. Methods. We ran disk models with stellar parameters adjusted to spectral types B9 through M2, using the radiation thermo-chemical disk modelling code ProDiMo. We explored also a small parameter space around a standard disk model, considering dust-to-gas mass ratio, disk gas mass, mixing coefficient for dust settling, flaring index, dust maximum size, and size power law distribution index. We produced convolved spectra at the resolution of Spitzer, IRS, JWST MIRI, and VLT VISIR spectrographs. We applied random noise derived from typical Spitzer spectra for a direct comparison with observations. Results. The strength of the mid-IR water lines correlates directly with the luminosity of the central star. The models show that it is possible to suppress the water emission; however, current observations are not sensitive enough to detect mid-IR lines in disks for most of the explored parameters. The presence of noise in the spectra, combined with the high continuum flux (noise level is proportional to the continuum flux), is the most likely explanation for the non-detections towards Herbig stars. Conclusions. Mid-IR spectra with resolution higher than 20 000 are needed to investigate water in protoplanetary disks. Intrinsic differences in disk structure, such as inner gaps, gas-to-dust ratio, dust size and distribution, and inner disk scale height, between Herbig and T Tauri star disks are able to explain a lower water detection rate in disks around Herbig stars.Publisher PDFPeer reviewe

Herschel/PACS observations of young sources in Taurus: the far-infrared counterpart of optical jets

Observations of the atomic and molecular line emission associated with jets and outflows emitted by young stellar objects can be used to trace the various evolutionary stages they pass through as they evolve to become main sequence stars. To understand the relevance of atomic and molecular cooling in shocks, and how accretion and ejection efficiency evolves with the source evolutionary state, we will study the far-infrared counterparts of bright optical jets associated with Class I and II sources in Taurus (T Tau, DG Tau A, DG Tau B, FS Tau A+B, and RW Aur). We have analysed Herschel/PACS observations of a number of atomic ([OI]63um, 145um, [CII]158um) and molecular (high-J CO, H2O, OH) lines, collected within the OTKP GASPS. To constrain the origin of the detected lines we have compared the FIR emission maps with the emission from optical-jets and millimetre-outflows, and the line fluxes and ratios with predictions from shock and disk models. All of the targets are associated with extended emission in the atomic lines correlated with the direction of the optical jet/mm-outflow. The atomic lines can be excited in fast dissociative J-shocks. The molecular emission, on the contrary, originates from a compact region, that is spatially and spectrally unresolved. Slow C- or J- shocks with high pre-shock densities reproduce the observed H2O and high-J CO lines; however, the disk and/or UV-heated outflow cavities may contribute to the emission. While the cooling is dominated by CO and H2O lines in Class 0 sources, [OI] becomes an important coolant as the source evolves and the environment is cleared. The cooling and mass loss rates estimated for Class II and I sources are one to four orders of magnitude lower than for Class 0 sources. This provides strong evidence to indicate that the outflow activity decreases as the source evolves.Comment: 18 pages, 9 figures, accepted by A&

Dietary assessment of British police force employees: A description of diet record coding procedures and cross-sectional evaluation of dietary energy intake reporting (The Airwave Health Monitoring Study)

Objectives: Dietary intake is a key aspect of occupational health. To capture the characteristics of dietary behaviour that is affected by occupational environment that may affect disease risk, a collection of prospective multiday dietary records is required. The aims of this paper are to: (1) collect multiday dietary data in the Airwave Health Monitoring Study, (2) describe the di etary coding procedures applied and (3) investigate the plausibility of dietary reporting in this occupational cohort. Design: A dietary coding protocol for this large-scale study was developed to minimise coding error rate. Participants (n 4412) who completed 7-day food records were included for cross-sectional analyses. Energy intake (EI) misreporting was estimated using the Goldberg method. Multivariate logistic regression models were applied to determine participant characteristics associated with EI misreporting. Setting: British police force employees enrolled (2007-2012) into the Airwave Health Monitoring Study. Results: The mean code error rate per food diary was 3.7% (SD 3.2%). The strongest predictors of EI under-reporting were body mass index (BMI) and physical activity. Compared with participants with BMI 30 kg/m2 had increased odds of being classified as under-reporting EI (men OR 5.20 95% CI 3.92 to 6.89; women OR 2.66 95% CI 1.85 to 3.83). Men and women in the highest physical activity category compared with the lowest were also more likely to be classified as underreporting (men OR 3.33 95% CI 2.46 to 4.50; women OR 4.34 95% CI 2.91 to 6.55). Conclusions: A reproducible dietary record coding procedure has been developed to minimise coding error in complex 7-day diet diaries. The prevalence of EI under-reporting is comparable with existing national UK cohorts and, in agreement with previous studies, classification of under-reporting was biased towards specific subgroups of participants

Continuum and line modeling of disks around young stars II. Line diagnostics for GASPS from the DENT grid

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. ... abbreviatedComment: 15 pages, 25 figures, accepted for publication in A&

GASPS observations of Herbig Ae/Be stars with PACS/Herschel. The atomic and molecular content of their protoplanetary discs

We observed a sample of 20 representative Herbig Ae/Be stars and five A-type debris discs with PACS onboard of Herschel. The observations were done in spectroscopic mode, and cover far-IR lines of [OI], [CII], CO, CH+, H2O and OH. We have a [OI]63 micron detection rate of 100% for the Herbig Ae/Be and 0% for the debris discs. [OI]145 micron is only detected in 25%, CO J=18-17 in 45% (and less for higher J transitions) of the Herbig Ae/Be stars and for [CII] 157 micron, we often found spatially variable background contamination. We show the first detection of water in a Herbig Ae disc, HD 163296, which has a settled disc. Hydroxyl is detected as well in this disc. CH+, first seen in HD 100546, is now detected for the second time in a Herbig Ae star, HD 97048. We report fluxes for each line and use the observations as line diagnostics of the gas properties. Furthermore, we look for correlations between the strength of the emission lines and stellar or disc parameters, such as stellar luminosity, UV and X-ray flux, accretion rate, PAH band strength, and flaring. We find that the stellar UV flux is the dominant excitation mechanism of [OI]63 micron, with the highest line fluxes found in those objects with a large amount of flaring and greatest PAH strength. Neither the amount of accretion nor the X-ray luminosity has an influence on the line strength. We find correlations between the line flux of [OI]63 micron and [OI]145 micron, CO J = 18-17 and [OI]6300 \AA, and between the continuum flux at 63 micron and at 1.3 mm, while we find weak correlations between the line flux of [OI]63 micron and the PAH luminosity, the line flux of CO J = 3-2, the continuum flux at 63 micron, the stellar effective temperature and the Brgamma luminosity. (Abbreviated version)Comment: 20 pages, 29 figures, accepted by Astronomy and Astrophysic