Showing 1–50 of 69 results for author: Thi, W -

HCN as a probe of the inner disk in a candidate proto-brown dwarf

Authors: B. Riaz, W. -F. Thi, M. N. Machida

Abstract: The detection of Keplerian rotation is rare among Class 0 protostellar systems. We have investigated the high-density tracer HCN as a probe of the inner disk in a Class 0 proto-brown dwarf candidate. Our ALMA high angular resolution observations show the peak in the HCN (3-2) line emission arises from a compact component near the proto-brown dwarf with a small bar-like structure and a deconvolved… ▽ More The detection of Keplerian rotation is rare among Class 0 protostellar systems. We have investigated the high-density tracer HCN as a probe of the inner disk in a Class 0 proto-brown dwarf candidate. Our ALMA high angular resolution observations show the peak in the HCN (3-2) line emission arises from a compact component near the proto-brown dwarf with a small bar-like structure and a deconvolved size of $\sim$50 au. Radiative transfer modelling indicates that this HCN feature is tracing the innermost, dense regions in the proto-brown dwarf where a small Keplerian disk is expected to be present. The limited velocity resolution of the observations, however, makes it difficult to confirm the rotational kinematics of this feature. A brightening in the HCN emission towards the core center suggests that HCN can survive in the gas phase in the inner, dense regions of the proto-brown dwarf. In contrast, modelling of the HCO$^{+}$ (3-2) line emission indicates that it originates from the outer pseudo-disk/envelope region and is centrally depleted. HCN line emission can reveal the small-scale structures and can be an efficient observational tool to study the inner disk properties in such faint compact objects where spatially resolving the disk is nearly impossible. △ Less

Submitted 18 May, 2024; originally announced May 2024.

Comments: Accepted in MNRAS Letters

2D disc modelling of the JWST line spectrum of EX Lupi

Authors: P. Woitke, W. -F. Thi, A. M. Arabhavi, I. Kamp, A. Kospal, P. Abraham

Abstract: We introduce a number of new theoretical approaches and improvements to the thermo-chemical disc modelling code ProDiMo to better predict and analyse the JWST line spectra of protoplanetary discs. We develop a new line escape probability method for disc geometries, a new scheme for dust settling, and discuss how to apply UV molecular shielding factors to photorates in 2D disc geometry. We show tha… ▽ More We introduce a number of new theoretical approaches and improvements to the thermo-chemical disc modelling code ProDiMo to better predict and analyse the JWST line spectra of protoplanetary discs. We develop a new line escape probability method for disc geometries, a new scheme for dust settling, and discuss how to apply UV molecular shielding factors to photorates in 2D disc geometry. We show that these assumptions are crucial for the determination of the gas heating/cooling rates and discuss how they affect the predicted molecular concentrations and line emissions. We apply our revised 2D models to the protoplanetary disc around the T Tauri star EX Lupi in quiescent state. We calculate infrared line emission spectra between 5 and 20 mic by CO, H2O, OH, CO2, HCN, C2H2 and H2, including lines of atoms and ions, using our full 2D predictions of molecular abundances, dust opacities, gas and dust temperatures. We develop a disc model with a slowly increasing surface density structure around the inner rim that can simultaneously fit the spectral energy distribution, the overall shape of the JWST spectrum of EX Lupi, and the main observed molecular characteristics in terms of column densities, emitting areas and molecular emission temperatures, which all result from one consistent disc model. The spatial structure of the line emitting regions of the different molecules is discussed. High abundances of HCN and C2H2 are caused in the model by stellar X-ray irradiation of the gas around the inner rim. △ Less

Submitted 30 November, 2023; originally announced November 2023.

Comments: accepted by A&A, 25 pages, 15 figures, 7 tables

First observations of warm and cold methanol in Class 0/I proto-brown dwarfs

Authors: B. Riaz, W. -F. Thi, M. N. Machida

Abstract: We present results from the first molecular line survey to search for the fundamental complex organic molecule, methanol (CH$_{3}$OH), in 14 Class 0/I proto-brown dwarfs (proto-BDs). IRAM 30-m observations over the frequency range of 92-116 GHz and 213-280 GHz have revealed emission in 14 CH$_{3}$OH transition lines, at upper state energy level, E$_{upper}\sim$7-49 K, and critical densities,… ▽ More We present results from the first molecular line survey to search for the fundamental complex organic molecule, methanol (CH$_{3}$OH), in 14 Class 0/I proto-brown dwarfs (proto-BDs). IRAM 30-m observations over the frequency range of 92-116 GHz and 213-280 GHz have revealed emission in 14 CH$_{3}$OH transition lines, at upper state energy level, E$_{upper}\sim$7-49 K, and critical densities, $n_{crit}$ of 10$^{5}$ to 10$^{9}$ cm$^{-3}$. The most commonly detected lines are at E$_{upper} <$ 20 K, while 11 proto-BDs also show emission in the higher excitation lines at E$_{upper}\sim$21-49 K and $n_{crit}\sim$10$^{5}$ to 10$^{8}$ cm$^{-3}$. In comparison with the brown dwarf formation models, the high excitation lines likely probe the warm ($\sim$25-50 K) corino region at $\sim$10-50 au in the proto-BDs, while the low-excitation lines trace the cold ($<$ 20 K) gas at $\sim$50-150 au. The column density for the cold component is an order of magnitude higher than the warm component. The CH$_{3}$OH ortho-to-para ratios range between $\sim$0.3-2.3. The volume-averaged CH$_{3}$OH column densities show a rise with decreasing bolometric luminosity among the proto-BDs, with the median column density higher by a factor of $\sim$3 compared to low-mass protostars. Emission in high-excitation (E$_{upper}>$ 25 K) CH$_{3}$OH lines together with the model predictions suggest that a warm corino is present in $\sim$78\% of the proto-BDs in our sample. The remaining show evidence of only the cold component, possibly due to the absence of a strong, high-velocity jet that can stir up the warm gas around it. △ Less

Submitted 8 May, 2023; originally announced May 2023.

Comments: Accepted in MNRAS

Simulation of CH$_3$OH ice UV photolysis under laboratory conditions

Authors: W. R. M. Rocha, P. Woitke, S. Pilling, W. -F. Thi, J. K. Jørgensen, L. E. Kristensen, G. Perotti, I. Kamp

Abstract: Methanol is the most complex molecule securely identified in interstellar ices and is a key chemical species for understanding chemical complexity in astrophysical environments. Important aspects of the methanol ice photochemistry are still unclear such as the branching ratios and photo-dissociation cross-sections at different temperatures and irradiation fluxes. This work aims at a quantitative a… ▽ More Methanol is the most complex molecule securely identified in interstellar ices and is a key chemical species for understanding chemical complexity in astrophysical environments. Important aspects of the methanol ice photochemistry are still unclear such as the branching ratios and photo-dissociation cross-sections at different temperatures and irradiation fluxes. This work aims at a quantitative agreement between laboratory experiments and astrochemical modelling of the CH3OH ice UV photolysis. This work also allows us to better understand which processes govern the methanol ice photochemistry present in laboratory experiments. We use ProDiMo to simulate the conditions of laboratory measurements. The simulations start with simple chemistry consisting only of methanol ice and helium to mimic the residual gas in the experimental chamber. A surface chemical network enlarged by photo-dissociation reactions is used to study the chemical reactions within the ice. Additionally, different surface chemistry parameters (surface competition, tunnelling, thermal diffusion and reactive desorption) are adopted to check those that reproduce the experimental results. The chemical models with ProDiMo can reproduce the methanol ice destruction via UV photodissociation at temperatures of 20, 30, 50 and 70 K as observed in the experiments. We note that the results are sensitive to different branching ratios after photolysis and to the mechanisms of reactive desorption. In the simulations of a molecular cloud at 20 K, we observed an increase in the methanol gas abundance of one order of magnitude, with a similar decrease in the solid-phase abundance. Comprehensive astrochemical models provide new insights into laboratory experiments as the quantitative understanding of the processes that govern the reactions within the ice. Ultimately, those insights can help to better interpret astronomical observations. △ Less

Submitted 8 February, 2023; originally announced February 2023.

Comments: Accepted for publication in A&A on 08-Feb-2023

Journal ref: A&A 673, A70 (2023)

Interpreting molecular hydrogen and atomic oxygen line emission of T Tauri disks with photoevaporative disk-wind models

Authors: Ch. Rab, M. Weber, T. Grassi, B. Ercolano, G. Picogna, P. Caselli, W. -F. Thi, I. Kamp, P. Woitke

Abstract: Winds in protoplanetary disks play an important role in their evolution and dispersal. However, what physical process is driving the winds is still unclear (i.e. magnetically vs thermally driven), and can only be understood by directly confronting theoretical models with observational data. We use hydrodynamic photoevaporative disk-wind models and post-process them with a thermo-chemical model to… ▽ More Winds in protoplanetary disks play an important role in their evolution and dispersal. However, what physical process is driving the winds is still unclear (i.e. magnetically vs thermally driven), and can only be understood by directly confronting theoretical models with observational data. We use hydrodynamic photoevaporative disk-wind models and post-process them with a thermo-chemical model to produce synthetic observables for the o-H$_2$ at 2.12 micron and [OI] at 0.63 micron spectral lines and directly compare the results to a sample of observations. Our photoevaporative disk-wind model is consistent with the observed signatures of the blueshifted narrow low-velocity component (NLVC), which is usually associated with slow disk winds, for both tracers. Only for one out of seven targets that show blueshifted NLVCs does the photoevaporative model fail to explain the observed line kinematics. Our results also indicate that interpreting spectral line profiles by simple methods, such as the thin-disk approximation, to determine the line emitting region can yield misleading conclusions. The photoevaporative disk-wind models are largely consistent with the studied observational data set, but it is not possible to clearly discriminate between different wind-driving mechanisms. Further improvements to the models, such as consistent modelling of the dynamics and chemistry and detailed modelling of individual targets would be beneficial. Furthermore, a direct comparison of magnetically driven disk-wind models to the observational data set is necessary in order to determine whether or not spatially unresolved observations of multiple wind tracers are sufficient to discriminate between theoretical models. △ Less

Submitted 18 November, 2022; v1 submitted 27 October, 2022; originally announced October 2022.

Comments: 14 pages, 13 figures, accepted for publication in A&A

Journal ref: A&A 668, A154 (2022)

Mixing and diffusion in protoplanetary disc chemistry

Authors: P. Woitke, A. M. Arabhavi, I. Kamp, W. -F. Thi

Abstract: We develop a simple iterative scheme to include vertical turbulent mixing and diffusion in ProDiMo thermo-chemical models for protoplanetary discs. The models are carefully checked for convergence toward the time-independent solution of the reaction-diffusion equations, as e.g. used in exoplanet atmosphere models. A series of five T Tauri disc models is presented where we vary the mixing parameter… ▽ More We develop a simple iterative scheme to include vertical turbulent mixing and diffusion in ProDiMo thermo-chemical models for protoplanetary discs. The models are carefully checked for convergence toward the time-independent solution of the reaction-diffusion equations, as e.g. used in exoplanet atmosphere models. A series of five T Tauri disc models is presented where we vary the mixing parameter α mix from 0 to 0.01 and take into account (a) the radiative transfer feedback of the opacities of icy grains that are mixed upward and (b) the feedback of the changing molecular abundances on the gas temperature structure caused by exothermic reactions and increased line heating/cooling. We see considerable changes of the molecular and ice concentrations in the disc. The most abundant species (H2, CH4, CO, the neutral atoms in higher layers, and the ices in the midplane) are transported both up and down, and at the locations where these abundant chemicals finally decompose, for example by photo processes, the release of reaction products has important consequences for all other molecules. This generally creates a more active chemistry, with a richer mixture of ionised, atomic, molecular and ice species and new chemical pathways that are not relevant in the unmixed case. We discuss the impact on three spectral observations caused by mixing and find that (i) icy grains can reach the observable disc surface where they cause ice absorption and emission features at IR to far-IR wavelengths, (ii) mixing increases the concentrations of certain neutral molecules observable by mid-IR spectroscopy, in particular OH, HCN and C2H2, and (iii) mixing can change the optical appearance of CO in ALMA line images and channel maps, where strong mixing would cause the CO molecules to populate the distant midplane. △ Less

Submitted 25 September, 2022; originally announced September 2022.

Comments: 19 pages, 13 Figures, accepted by A&A

Journal ref: A&A 668, A164 (2022)

A survey of HDCO and D$_{2}$CO towards Class 0/I proto-brown dwarfs

Authors: B. Riaz, W. -F. Thi

Abstract: Deuterium fractionation can constrain the physical and chemical conditions at the early stage of brown dwarf formation. We present IRAM 30m observations over a wide frequency range of 213-279 GHz of singly and doubly deuterated species of formaldehyde (HDCO and D$_{2}$CO) towards Class 0/I proto-brown dwarfs (proto-BDs). Multiple low-excitation HDCO and D$_{2}$CO transition lines with upper energy… ▽ More Deuterium fractionation can constrain the physical and chemical conditions at the early stage of brown dwarf formation. We present IRAM 30m observations over a wide frequency range of 213-279 GHz of singly and doubly deuterated species of formaldehyde (HDCO and D$_{2}$CO) towards Class 0/I proto-brown dwarfs (proto-BDs). Multiple low-excitation HDCO and D$_{2}$CO transition lines with upper energy level $\leq$ 40 K are detected. The D$_{2}$CO/HDCO, HDCO/H$_{2}$CO, D$_{2}$CO/H$_{2}$CO abundance ratios range between 0.01 and 2.5 for the proto-BDs, similar to the range seen in low-mass protostars. The highest ratios of D$_{2}$CO/HDCO $\sim$1.3-2.5 are measured for two Stage 0 proto-BDs. These objects could possess a warm corino, similar to the few hot corino cases reported among Class 0 protostars. The mean D$_{2}$CO/HDCO, D$_{2}$CO/H$_{2}$CO, and HDCO/H$_{2}$CO ratios for the proto-BDs are comparatively higher than the range predicted by the current gas-grain chemical models, indicating that HDCO and D$_{2}$CO are formed via grain surface reactions in the dense and cold interiors of the proto-BDs at an early formation stage. △ Less

Submitted 6 June, 2022; originally announced June 2022.

Comments: Accepted in MNRAS

Deuterium chemistry and D/H ratios in Class 0/I proto-brown dwarfs

Authors: B. Riaz, W. -F. Thi

Abstract: We have conducted the first extensive observational survey of several deuterated species in 16 Class 0/I proto-brown dwarfs (proto-BDs) and 4 Class Flat/Class II brown dwarfs. Observations were obtained with the IRAM 30m telescope in the DCO$^{+}$ (3-2), DCN (3-2), DNC (3-2), and N$_{2}$D$^{+}$ (3-2) lines. The DCO$^{+}$/H$^{13}$CO$^{+}$, DCN/H$^{13}$CN, and DNC/HN$^{13}$C ratios are comparatively… ▽ More We have conducted the first extensive observational survey of several deuterated species in 16 Class 0/I proto-brown dwarfs (proto-BDs) and 4 Class Flat/Class II brown dwarfs. Observations were obtained with the IRAM 30m telescope in the DCO$^{+}$ (3-2), DCN (3-2), DNC (3-2), and N$_{2}$D$^{+}$ (3-2) lines. The DCO$^{+}$/H$^{13}$CO$^{+}$, DCN/H$^{13}$CN, and DNC/HN$^{13}$C ratios are comparatively higher and show a narrower range than the DCO$^{+}$/HCO$^{+}$, DCN/HCN, and DNC/HNC ratios, respectively. The mean D/H ratios for the proto-BDs derived from these molecules are in the range of $\sim$0.02-3. Both low-temperature gas-phase ion-molecule deuteron transfer and grain surface reactions are required to explain the enhanced deuterium fractionation. The very dense and cold ($n_{H_{2}}$ $\geq$10$^{6}$ cm$^{-3}$, T $\leq$10 K) interior of the proto-BDs provide the suitable conditions for efficient deuterium fractionation in these cores. There is no correlation between the D/H ratios and the CO depletion factor, with the exception of the DCN/HCN ratios that show a strong anti-correlation possibly due to the difference in the peak emitting regions of the DCN and HCN molecules. Over a wide range in the bolometric luminosities spanning $\sim$0.002--40 L$_{\odot}$, we find a trend of higher DCO$^{+}$/HCO$^{+}$ (r = -0.7) and DCN/HCN (r = -0.6) ratios, nearly constant DNC/HNC (r = -0.4) and DNC/HN$^{13}$C (r = -0.1) ratios, lower N$_{2}$D$^{+}$/N$_{2}$H$^{+}$ ratios (r = 0.6) in the proto-BDs compared to protostars. Only one Class II brown dwarf shows emission in the DCO$^{+}$ (3-2) line. No correlation is seen between the D/H ratios and the evolutionary stage. △ Less

Submitted 18 January, 2022; originally announced January 2022.

Comments: Accepted in MNRAS

First CH$_{3}$D detection in Class 0/I proto-brown dwarfs: constraints on CH$_{4}$ abundances

Authors: B. Riaz, W. -F. Thi

Abstract: We report the first detection in the ${\it J}_{\it K}$ = 1$_{0}$ - 0$_{0}$ rotational transition line of CH$_{3}$D towards three Class 0/I proto-brown dwarfs (proto-BDs) from IRAM 30 m observations. Assuming a rotational temperature of 25 K, the CH$_{3}$D abundances (relative to H$_{2}$) are in the range of (2.3--14.5) $\times$ 10$^{-7}$. The CH$_{4}$ abundances derived from the CH$_{3}$D abundanc… ▽ More We report the first detection in the ${\it J}_{\it K}$ = 1$_{0}$ - 0$_{0}$ rotational transition line of CH$_{3}$D towards three Class 0/I proto-brown dwarfs (proto-BDs) from IRAM 30 m observations. Assuming a rotational temperature of 25 K, the CH$_{3}$D abundances (relative to H$_{2}$) are in the range of (2.3--14.5) $\times$ 10$^{-7}$. The CH$_{4}$ abundances derived from the CH$_{3}$D abundances and assuming the DCO$^+$/HCO$^+$ ratios are in the range of (0.05--4.8) $\times $10$^{-5}$. The gas-phase formation of CH$_{3}$D via CH$_2$D$^+$ is enhanced at high densities of 10$^{8}$--10$^{10}$ cm$^{-3}$ and our observations are likely probing the innermost dense and warm regions in proto-BDs. Thermal and/or non-thermal desorption can return the CH$_{3}$D and CH$_{4}$ molecules formed at an early stage on grain surfaces to the gas-phase. The gas phase abundances are indicative of warm carbon-chain chemistry in proto-BDs where carbon-chain molecules are synthesized in a lukewarm ($\sim$20-30 K) region close to the central source. △ Less

Submitted 18 January, 2022; originally announced January 2022.

Comments: Accepted in MNRAS Letters

The GRAVITY young stellar object survey III. The dusty disk of RY Lup

Authors: GRAVITY Collaboration, Y. -I. Bouarour, K. Perraut, F. Ménard, W. Brandner, A. Caratti o Garatti, P. Caselli, E. van Dishoeck, C. Dougados, R. Garcia-Lopez, R. Grellmann, T. Henning, L. Klarmann, L. Labadie, A. Natta, J. Sanchez-Bermudez, W. -F. Thi, P. T. de Zeeuw, A. Amorim, M. Bauböck, M. Benisty, J. -P. Berger, Y. Clenet, V. Coudé du Foresto, G. Duvert , et al. (33 additional authors not shown)

Abstract: We use PIONIER data from the ESO archive and GRAVITY data that were obtained in June 2017 with the four 8m telescopes. We use a parametric disk model and the 3D radiative transfer code MCFOST to reproduce the Spectral Energy Distribution and match the interferometric observations. To match the SED , our model requires a stellar luminosity of 2.5 Lsun, higher than any previously determined values.… ▽ More We use PIONIER data from the ESO archive and GRAVITY data that were obtained in June 2017 with the four 8m telescopes. We use a parametric disk model and the 3D radiative transfer code MCFOST to reproduce the Spectral Energy Distribution and match the interferometric observations. To match the SED , our model requires a stellar luminosity of 2.5 Lsun, higher than any previously determined values. Such a high value is needed to accommodate the circumstellar extinction caused by the highly inclined disk, which has been neglected in previous studies. While using an effective temperature of 4800 K determined through high-resolution spectroscopy, we derive a stellar radius of 2.29 Rsun. These revised fundamental parameters, when combined with the mass estimates available , lead to an age of 0.5-2.0 Ma for RY Lup, in better agreement with the age of the Lupus association than previous determinations. Our disk model nicely reproduces the interferometric GRAVITY data and is in good agreement with the PIONIER ones. We derive an inner rim location at 0.12~au from the central star. This model corresponds to an inclination of the inner disk of 50deg, which is in mild tension with previous determinations of a more inclined outer disk from SPHERE (70 deg in NIR) and ALMA(67 $\pm$5 deg) images, but consistent with the inclination determination from the ALMA CO spectra (55$\pm$5deg). Increasing the inclination of the inner disk to 70 deg leads to a higher line-of-sight extinction and therefore requires a higher stellar luminosity of 4.65 Lsun to match the observed flux levels. This luminosity would translate to a stellar radius of 3.13~Rsun, leading to an age of 2-3~Ma, and a stellar mass of about 2 Msun, in disagreement with the observed dynamical mass estimate of 1.3-1.5 Msun. Critically, this high-inclination inner disk model also fails to reproduce the visibilities observed with GRAVITY. △ Less

Submitted 19 August, 2020; originally announced August 2020.

Comments: Accepted for publication in A&A; 11 pages, 5 figures, 3 tables

Journal ref: A&A 642, A162 (2020)

Interpreting high spatial resolution line observations of planet-forming disks with gaps and rings -- The case of HD 163296

Authors: Ch. Rab, I. Kamp, C. Dominik, C. Ginski, G. A. Muro-Arena, W. -F. Thi, L. B. F. M. Waters, P. Woitke

Abstract: Spatially resolved continuum observations of planet-forming disks show prominent ring and gap structures in their dust distribution. However, the picture from gas observations is much less clear and constraints on the radial gas density structure (i.e. gas gaps) remain rare and uncertain. We want to investigate the importance of thermo-chemical processes for the interpretation of high-spatial-reso… ▽ More Spatially resolved continuum observations of planet-forming disks show prominent ring and gap structures in their dust distribution. However, the picture from gas observations is much less clear and constraints on the radial gas density structure (i.e. gas gaps) remain rare and uncertain. We want to investigate the importance of thermo-chemical processes for the interpretation of high-spatial-resolution gas observations of planet-forming disks and their impact on derived gas properties. We apply the radiation thermo-chemical disk code ProDiMo (PROtoplanetary DIsk MOdel) to model self-consistently the dust and gas disk of HD 163296, using the DSHARP gas and dust observations. With this model we investigate the impact of dust gaps and gas gaps, considering chemistry and heating/cooling processes, on the observables and the derived gas properties. We find distinct peaks in the radial line intensity profiles of the CO line data of HD 163296 at the location of the dust gaps. Our model indicates that those peaks are not only a consequence of a gas temperature increase within the gaps but are mainly caused by the absorption of line emission from the back side of the disk by the dust rings. For two of the three prominent dust gaps in HD 163296, we find that thermo-chemical effects are negligible for deriving density gradients via measurements of the rotation velocity. However, for the gap with the highest dust depletion, the temperature gradient can be dominant and needs to be considered to derive accurate gas density profiles. Self-consistent gas and dust thermo-chemical modelling in combination with high-quality observations of multiple molecules are necessary to accurately derive gas gap depths and shapes. This is crucial to determine the origin of gaps and rings in planet-forming disks and to improve the mass estimates of forming planets if they are the cause of the gap. △ Less

Submitted 3 September, 2020; v1 submitted 13 August, 2020; originally announced August 2020.

Comments: 12 pages, 9 figures, accepted for publication in A&A

Journal ref: A&A 642, A165 (2020)

A family portrait of disk inner rims around Herbig Ae/Be stars: Hunting for warps, rings, self shadowing, and misalignments in the inner astronomical units

Authors: J. Kluska, J. -P. Berger, F. Malbet, B. Lazareff, M. Benisty, J. -B. Le Bouquin, O. Absil, F. Baron, A. Delboulbé, G. Duvert, A. Isella, L. Jocou, A. Juhasz, S. Kraus, R. Lachaume, F. Ménard, R. Millan-Gabet, J. D. Monnier, T. Moulin, K. Perraut, C. Pinte, S. Rochat, F. Soulez, M. Tallon, W. -F. Thi , et al. (3 additional authors not shown)

Abstract: The innermost astronomical unit in protoplanetary disks is a key region for stellar and planet formation, as exoplanet searches have shown a large occurrence of close-in planets that are located within the first au around their host star. We aim to reveal the morphology of the disk inner rim using near-infrared interferometric observations with milli-arcsecond resolution provided by infrared inter… ▽ More The innermost astronomical unit in protoplanetary disks is a key region for stellar and planet formation, as exoplanet searches have shown a large occurrence of close-in planets that are located within the first au around their host star. We aim to reveal the morphology of the disk inner rim using near-infrared interferometric observations with milli-arcsecond resolution provided by infrared interferometry. We provide reconstructed images of 15 objects selected from the Herbig AeBe survey carried out with PIONIER at the VLTI, using SPARCO. We find that 40% of the systems are centrosymmetric at the angular resolution of the observations. For the rest of the objects, we find evidence for asymmetric emission due to moderate-to-strong inclination of a disk-like structure for 30% of the objects and noncentrosymmetric morphology due to a nonaxisymmetric and possibly variable environment (30%). Among the systems with a disk-like structure, 20% show a resolved dust-free cavity. The image reconstruction process is a powerful tool to reveal complex disk inner rim morphologies. At the angular resolution reached by near-infrared interferometric observations, most of the images are compatible with a centrally peaked emission (no cavity). For the most resolved targets, image reconstruction reveals morphologies that cannot be reproduced by generic parametric models. Moreover, the nonaxisymmetric disks show that the spatial resolution probed by optical interferometers makes the observations of the near-infrared emission sensitive to temporal evolution with a time-scale down to a few weeks. The evidence of nonaxisymmetric emission that cannot be explained by simple inclination and radiative transfer effects requires alternative explanations, such as a warping of the inner disks. Interferometric observations can, therefore, be used to follow the evolution of the asymmetry of those disks at a sub-au scale. △ Less

Submitted 3 April, 2020; originally announced April 2020.

Comments: Accepted in A&A, 24pages, 22 figures

Journal ref: A&A 636, A116 (2020)

A model exploration of NIR ro-vibrational CO emission as a tracer of inner cavities in protoplanetary disks

Authors: S. Antonellini, A. Banzatti, I. Kamp, W. -F. Thi, P. Woitke

Abstract: Near-IR observations of protoplanetary disks provide information about the properties of the inner disk. High resolution spectra of abundant molecules such as CO can be used to determine the disk structure in the warm inner parts. The $v2/v1$ ro-vibrational ratio of $v_{1-0}$ and $v_{2-1}$ transitions has been recently observed to follow distinct trends with the CO emitting radius, in a sample of… ▽ More Near-IR observations of protoplanetary disks provide information about the properties of the inner disk. High resolution spectra of abundant molecules such as CO can be used to determine the disk structure in the warm inner parts. The $v2/v1$ ro-vibrational ratio of $v_{1-0}$ and $v_{2-1}$ transitions has been recently observed to follow distinct trends with the CO emitting radius, in a sample of TTauri and Herbig disks; these trends have been empirically interpreted as due to inner disk depletion from gas and dust. In this work we use existing thermo-chemical disk models to explore the interpretation of these observed trends in ro-vibrational CO emission. We use the radiation thermo-chemical code ProDiMo, exploring a set of previously published models with different disk properties and varying one parameter at a time: the inner radius, the dust-to-gas mass ratio, the gas mass. In addition, we use models where we change the surface density power law index, and employ a larger set of CO ro-vibrational levels, including also fluorescence from the first electronic state. We investigate these models for both TTauri and Herbig star disks. Finally, we include a set of DIANA models for individual TTauri and Herbig disks which were constructed to reproduce a large set of multi-wavelength observations. △ Less

Submitted 26 March, 2020; originally announced March 2020.

Journal ref: A&A 637, A29 (2020)

The infrared line-emitting regions of T Tauri protoplanetary disks

Authors: A. J. Greenwood, I. Kamp, L. B. F. M. Waters, P. Woitke, W. -F. Thi

Abstract: Mid-infrared molecular line emission detected with the Spitzer Space Telescope is often interpreted using slab models. However, we need to understand the mid-infrared line emission in 2D disk models, such that we gain information about from where the lines are being emitted and under which conditions, such that we gain information about number densities, temperatures, and optical depths in both th… ▽ More Mid-infrared molecular line emission detected with the Spitzer Space Telescope is often interpreted using slab models. However, we need to understand the mid-infrared line emission in 2D disk models, such that we gain information about from where the lines are being emitted and under which conditions, such that we gain information about number densities, temperatures, and optical depths in both the radial and vertical directions. In this paper, we introduce a series of 2D thermochemical models of a prototypical T Tauri protoplanetary disk, in order to examine how sensitive the line-emitting regions are to changes in the UV and X-ray fluxes, the disk flaring angle, dust settling, and the dust-to-gas ratio. These all affect the heating of the inner disk, and thus can affect the mid-infrared spectral lines. Using the ProDiMo and FLiTs codes, we produce a series of 2D thermochemical disk models. We find that there is often a significant difference between the gas and dust temperatures in the line emitting regions, and we illustrate that the size of the line emitting regions is relatively robust against changes in the stellar and disk parameters (namely, the UV and X-ray fluxes, the flaring angle, and dust settling). These results demonstrate the potential for localized variations in the line-emitting region to greatly affect the resulting spectra and line fluxes, and the necessity of allowing for such variations in our models. △ Less

Submitted 11 October, 2019; originally announced October 2019.

Comments: 39 pages, 39 figures, 4 tables

The effects of dust evolution on disks in the mid-IR

Authors: A. J. Greenwood, I. Kamp, L. B. F. M. Waters, P. Woitke, W. -F. Thi

Abstract: In this paper, we couple together the dust evolution code two-pop-py with the thermochemical disk modelling code ProDiMo. We create a series of thermochemical disk models that simulate the evolution of dust over time from 0.018 Myr to 10 Myr, including the radial drift, growth, and settling of dust grains. We examine the effects of this dust evolution on the mid-infrared gas emission, focussing on… ▽ More In this paper, we couple together the dust evolution code two-pop-py with the thermochemical disk modelling code ProDiMo. We create a series of thermochemical disk models that simulate the evolution of dust over time from 0.018 Myr to 10 Myr, including the radial drift, growth, and settling of dust grains. We examine the effects of this dust evolution on the mid-infrared gas emission, focussing on the mid-infrared spectral lines of C2H2, CO2, HCN, NH3, OH, and H2O that are readily observable with Spitzer and the upcoming E-ELT and JWST. The addition of dust evolution acts to increase line fluxes by reducing the population of small dust grains. We find that the spectral lines of all species except C2H2 respond strongly to dust evolution, with line fluxes increasing by more than an order of magnitude across the model series as the density of small dust grains decreases over time. The C2H2 line fluxes are extremely low due to a lack of abundance in the infrared line-emitting regions, despite C2H2 being commonly detected with Spitzer, suggesting that warm chemistry in the inner disk may need further investigation. Finally, we find that the CO2 flux densities increase more rapidly than the other species as the dust disk evolves. This suggests that the flux ratios of CO2 to other species may be lower in disks with less-evolved dust populations. △ Less

Submitted 1 April, 2019; v1 submitted 29 March, 2019; originally announced March 2019.

Comments: 13 pages, 9 figures, accepted in A&A

Journal ref: A&A 626, A6 (2019)

4MOST Survey Strategy Plan

Authors: G. Guiglion, C. Battistini, C. P. M. Bell, T. Bensby, T. Boller, C. Chiappini, J. Comparat, N. Christlieb, R. Church, M. -R. L. Cioni, L. Davies, T. Dwelly, R. S. de Jong, S. Feltzing, A. Gueguen, L. Howes, M. Irwin, I. Kushniruk, M. I Lam, J. Liske, R. McMahon, A. Merloni, P. Norberg, A. S. G. Robotham, O. Schnurr , et al. (8 additional authors not shown)

Abstract: The current status of and motivation for the 4MOST Survey Strategy, as developed by the Consortium science team, are presented here. Key elements of the strategy are described, such as sky coverage, number of visits and total exposure times in different parts of the sky, and how to deal with different observing conditions. The task of organising the strategy is not simple, with many different surv… ▽ More The current status of and motivation for the 4MOST Survey Strategy, as developed by the Consortium science team, are presented here. Key elements of the strategy are described, such as sky coverage, number of visits and total exposure times in different parts of the sky, and how to deal with different observing conditions. The task of organising the strategy is not simple, with many different surveys that have vastly different target brightnesses and densities, sample completeness levels, and signal-to-noise requirements. We introduce here a number of concepts that we will use to ensure all surveys are optimised. Astronomers who are planning to submit a Participating Survey proposal are strongly encouraged to read this article and any relevant 4MOST Survey articles in this issue of The Messenger such that they can optimally complement and benefit from the planned surveys of the 4MOST Consortium. △ Less

Submitted 1 April, 2019; v1 submitted 6 March, 2019; originally announced March 2019.

Comments: Part of the 4MOST issue of The Messenger, published in preparation of 4MOST Community Workshop, see http://www.eso.org/sci/meetings/2019/4MOST.html

Journal ref: The Messenger, 2019, 175, 17

4MOST Scientific Operations

Authors: C. J. Walcher, M. Banerji, C. Battistini, C. P. M. Bell, O. Bellido-Tirado, T. Bensby, J. M. Bestenlehner, T. Boller, J. Brynnel, A. Casey, C. Chiappini, N. Christlieb, R. Church, M. -R. L. Cioni, S. Croom, J. Comparat, L. J. M. Davies, R. S. de Jong, T. Dwelly, H. Enke, S. Feltzing, D. Feuillet, M. Fouesneau, D. Ford, S. Frey , et al. (43 additional authors not shown)

Abstract: The 4MOST instrument is a multi-object spectrograph that will address Galactic and extragalactic science cases simultaneously by observing targets from a large number of different surveys within each science exposure. This parallel mode of operation and the survey nature of 4MOST require some distinct 4MOST-specific operational features within the overall operations model of ESO. The main feature… ▽ More The 4MOST instrument is a multi-object spectrograph that will address Galactic and extragalactic science cases simultaneously by observing targets from a large number of different surveys within each science exposure. This parallel mode of operation and the survey nature of 4MOST require some distinct 4MOST-specific operational features within the overall operations model of ESO. The main feature is that the 4MOST Consortium will deliver, not only the instrument, but also contractual services to the user community, which is why 4MOST is also described as a facility. This white paper concentrates on information particularly useful to answering the forthcoming Call for Letters of Intent. △ Less

Submitted 1 April, 2019; v1 submitted 6 March, 2019; originally announced March 2019.

Comments: Part of the 4MOST issue of The Messenger, published in preparation of 4MOST Community Workshop, see http://www.eso.org/sci/meetings/2019/4MOST.html

Journal ref: The Messenger, 2019, 175, 12

4MOST: Project overview and information for the First Call for Proposals

Authors: R. S. de Jong, O. Agertz, A. Agudo Berbel, J. Aird, D. A. Alexander, A. Amarsi, F. Anders, R. Andrae, B. Ansarinejad, W. Ansorge, P. Antilogus, H. Anwand-Heerwart, A. Arentsen, A. Arnadottir, M. Asplund, M. Auger, N. Azais, D. Baade, G. Baker, S. Baker, E. Balbinot, I. K. Baldry, M. Banerji, S. Barden, P. Barklem , et al. (313 additional authors not shown)

Abstract: We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolut… ▽ More We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = λ/Δλ\sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations. △ Less

Submitted 1 April, 2019; v1 submitted 6 March, 2019; originally announced March 2019.

Comments: Part of the 4MOST issue of The Messenger, published in preparation of 4MOST Community Workshop, see http://www.eso.org/sci/meetings/2019/4MOST.html

Journal ref: The Messenger, 2019, 175, 3

Observing the gas component of circumplanetary disks around wide-orbit planet-mass companions in the (sub)mm regime

Authors: Ch. Rab, I. Kamp, C. Ginski, N. Oberg, G. A. Muro-Arena, C. Dominik, L. B. F. M. Waters, W. -F. Thi, P. Woitke

Abstract: Several detections of wide-orbit planet-mass/sub-stellar companions around young solar-like stars were reported in the last decade. The origin of those possible planets is still unclear but accretion tracers and VLT/SPHERE observations indicate that they are surrounded by circumplanetary material or even a circumplanetary disk. We want to investigate if the gas component of disks around wide-orbit… ▽ More Several detections of wide-orbit planet-mass/sub-stellar companions around young solar-like stars were reported in the last decade. The origin of those possible planets is still unclear but accretion tracers and VLT/SPHERE observations indicate that they are surrounded by circumplanetary material or even a circumplanetary disk. We want to investigate if the gas component of disks around wide-orbit companions is detectable with current and future (sub)mm telescopes and what constraints such gas observations can provide on the nature of the circumplanetary material and on the mass of the companion. We applied the radiation thermo-chemical disk code ProDiMo to model the dust and gas component of passive circumplanetary disks and produced realistic synthetic observables. We considered different companion properties, disk parameters and radiative environments and compared the resulting synthetic observables to telescope sensitivities and to existing dust observations. The main criterion for a successful detection is the size of the circumplanetary disk. At a distance of about 150 pc, a circumplanetary disk with an outer radius of about 10 au is detectable with ALMA in about 6 hours in optically thick CO lines. Other aspects such as the companion's luminosity, disk inclination and background radiation fields are also relevant and should be considered to optimize the observing strategy for detection experiments. For most of the known wide-orbit planet-mass companions, their maximum theoretical disk size of one third of the Hill radius would be sufficient to allow detection of CO lines. It is therefore feasible to detect their gas disks and constrain the mass of the companion through the kinematic signature. Even in the case of non-detections such observations will provide stringent constraints on disk size and gas mass, information crucial for formation theories. △ Less

Submitted 28 March, 2019; v1 submitted 11 February, 2019; originally announced February 2019.

Comments: 13 pages, 11 figures, accepted for publication in A&A

Journal ref: A&A 624, A16 (2019)

Constraining the stellar energetic particle flux in young solar-like stars

Authors: Ch. Rab, M. Padovani, M. Güdel, I. Kamp, W. -F. Thi, P. Woitke

Abstract: Anomalies in the abundance measurements of short lived radionuclides in meteorites indicate that the protosolar nebulae was irradiated by a large number of energetic particles ($E\gtrsim10\,$MeV), often called solar cosmic rays. The particle flux of the contemporary Sun cannot explain these anomalies, but, similar to \mbox{T Tauri} stars, the young Sun was more active and probably produced enough… ▽ More Anomalies in the abundance measurements of short lived radionuclides in meteorites indicate that the protosolar nebulae was irradiated by a large number of energetic particles ($E\gtrsim10\,$MeV), often called solar cosmic rays. The particle flux of the contemporary Sun cannot explain these anomalies, but, similar to \mbox{T Tauri} stars, the young Sun was more active and probably produced enough high energy particles. However, the stellar particle (SP) flux of young stars is essentially unknown. We model the impact of high-energy ionization sources on the chemistry of the circumstellar environment (disks and envelopes). The model includes X-ray radiative transfer and makes use of particle transport models to calculate the individual molecular hydrogen ionization rates. We study the impact on the chemistry via the ionization tracers HCO$^+$ and N$_2$H$^+$. We argue that spatially resolved observations of those molecules combined with detailed models allow for disentangling the contribution of the individual high-energy ionization sources and to put constraints on the SP flux in young stars. △ Less

Submitted 3 February, 2019; originally announced February 2019.

Comments: 2 pages, 1 figure, accepted for IAUS345 "Origins: From the Protosun to the First Steps of Life" proceedings

Consistent dust and gas models for protoplanetary disks III. Models for selected objects from the FP7 DIANA project

Authors: P. Woitke, I. Kamp, S. Antonellini, F. Anthonioz, C. Baldovin-Saveedra, A. Carmona, O. Dionatos, C. Dominik, J. Greaves, M. Güdel, J. D. Ilee, A. Liebhardt, F. Menard, M. Min, C. Pinte, C. Rab, L. Rigon, W. -F. Thi, N. Thureau, L. B. F. M. Waters

Abstract: The European FP7 project DIANA has performed a coherent analysis of a large set of observations from protoplanetary disks by means of thermo-chemical disk models. The collected data include extinction-corrected stellar UV and X-ray input spectra (as seen by the disk), photometric fluxes, low and high resolution spectra, interferometric data, emission line fluxes, line velocity profiles and line ma… ▽ More The European FP7 project DIANA has performed a coherent analysis of a large set of observations from protoplanetary disks by means of thermo-chemical disk models. The collected data include extinction-corrected stellar UV and X-ray input spectra (as seen by the disk), photometric fluxes, low and high resolution spectra, interferometric data, emission line fluxes, line velocity profiles and line maps. We define and apply a standardized modelling procedure to simultaneously fit all these data by state-of-the-art modelling codes (ProDiMo, MCFOST, MCMax) which solve the continuum and line radiative transfer, disk chemistry, and the heating & cooling balance for both the gas and the dust. We allow for up to two radial disk zones to obtain our best-fitting models that have about 20 free parameters. This approach is novel and unique in its completeness and level of consistency. In this paper, we present the results from pure SED fitting for 27 objects and from the all inclusive DIANA-standard models for 14 objects. We fit most infrared to millimeter emission line fluxes within a factor better than 3, simultaneously with SED, PAH features and radial brightness profiles extracted from images at various wavelengths. Our analysis shows a number of Herbig Ae and T Tauri stars with very cold and massive outer disks which are situated at least partly in the shadow of a tall and gas-rich inner disk. The disk masses derived are often in excess to previously published values, since these disks are partially optically thick even at millimeter wavelength and so cold that they emit less than in the Rayleigh-Jeans limit. Some line observations cannot be reproduced by the models, probably caused by foreground cloud absorption or object variability. Our data collection, the fitted physical disk parameters as well as the full model output are available at an online database (http://www.univie.ac.at/diana). △ Less

Submitted 6 December, 2018; originally announced December 2018.

Comments: 70 pages, 35 figures, 22 tables, accepted by PASP

Chemical tracers in proto-brown dwarfs: CN, HCN, and HNC observations

Authors: B. Riaz, W. -F. Thi, P. Caselli

Abstract: We present results from a study of nitrogen chemistry in Class 0/I proto-brown dwarfs (proto-BDs). We have used the IRAM 30 m telescope to observe the CN (2-1), HCN (3-2), and HNC (3-2) lines in 7 proto-BDs. All proto-BDs show a large CN/HCN abundance ratio of >20, and a HNC/HCN abundance ratio close to or larger than unity. The enhanced CN/HCN ratios can be explained by high UV flux originating f… ▽ More We present results from a study of nitrogen chemistry in Class 0/I proto-brown dwarfs (proto-BDs). We have used the IRAM 30 m telescope to observe the CN (2-1), HCN (3-2), and HNC (3-2) lines in 7 proto-BDs. All proto-BDs show a large CN/HCN abundance ratio of >20, and a HNC/HCN abundance ratio close to or larger than unity. The enhanced CN/HCN ratios can be explained by high UV flux originating from an active accretion zone in the proto-BDs. The larger than unity HNC/HCN ratio for the proto-BDs is likely caused by a combination of low temperature and high density. Both CN and HNC show a flat distribution with CO, indicating that these species can survive in regions where CO is depleted. We have investigated the correlations in the molecular abundances of these species for the proto-BDs with Class 0/I protostars. We find tentative trends of CN (HCN) abundances being about an order of magnitude higher (lower) in the proto-BDs compared to protostars. HNC for the proto-BDs shows a nearly constant abundance unlike the large spread of ~2 orders of magnitude seen for the protostars. Also notable is a rise in the HNC/HCN abundance ratio for the lowest luminosity objects, suggesting that this ratio is higher under low-temperature environments. None of the relatively evolved Class Flat/Class II brown dwarfs in our sample show emission in HNC. The HNC molecule can be considered as an efficient tracer to search and identify early stage sub-stellar mass objects. △ Less

Submitted 26 September, 2018; originally announced September 2018.

Comments: Accepted in MNRAS

Modelling mid-infrared molecular emission lines from T Tauri stars

Authors: P. Woitke, M. Min, W. -F. Thi, C. Roberts, A. Carmona, I. Kamp, F. Menard, C. Pinte

Abstract: We introduce a new modelling framework called FLiTs to simulate infrared line emission spectra from protoplanetary discs. This paper focuses on the mid-IR spectral region between 9.7 um to 40 um for T Tauri stars. The generated spectra contain several tens of thousands of molecular emission lines of H2O, OH, CO, CO2, HCN, C2H2, H2 and a few other molecules, as well as the forbidden atomic emission… ▽ More We introduce a new modelling framework called FLiTs to simulate infrared line emission spectra from protoplanetary discs. This paper focuses on the mid-IR spectral region between 9.7 um to 40 um for T Tauri stars. The generated spectra contain several tens of thousands of molecular emission lines of H2O, OH, CO, CO2, HCN, C2H2, H2 and a few other molecules, as well as the forbidden atomic emission lines of SI, SII, SIII, SiII, FeII, NeII, NeIII, ArII and ArIII. In contrast to previously published works, we do not treat the abundances of the molecules nor the temperature in the disc as free parameters, but use the complex results of detailed 2D ProDiMo disc models concerning gas and dust temperature structure, and molecular concentrations. FLiTs computes the line emission spectra by ray tracing in an efficient, fast and reliable way. The results are broadly consistent with R=600 Spitzer/IRS observational data of T Tauri stars concerning line strengths, colour, and line ratios. In order to achieve that agreement, however, we need to assume either a high gas/dust mass ratio of order 1000, or the presence of illuminated disc walls at distances of a few au. The molecules in these walls cannot be photo-dissociated easily by UV because of the large densities in the walls which favour their re-formation. Most observable molecular emission lines are found to be optically thick, rendering a standard analysis with column densities difficult. We find that the difference between gas and dust temperatures in the disc surface is important for the line formation. We briefly discuss the effects of C/O ratio and choice of chemical rate network on these results. Our analysis offers new ways to infer the chemical and temperature structure of T Tauri discs from future JWST/MIRI observations, and to possibly detect secondary illuminated disc walls based on their specific mid-IR molecular signature. △ Less

Submitted 25 July, 2018; v1 submitted 16 July, 2018; originally announced July 2018.

Comments: accepted by A&A, 13 figures, 4 tables

Journal ref: A&A 618, A57 (2018)

X-ray radiative transfer in protoplanetary disks - The role of dust and X-ray background fields

Authors: Ch. Rab, M. Güdel, P. Woitke, I. Kamp, W. -F. Thi, M. Min, G. Aresu, R. Meijerink

Abstract: The X-ray luminosities of T Tauri stars are about two to four orders of magnitude higher than the luminosity of the contemporary Sun. As these stars are born in clusters, their disks are not only irradiated by their parent star but also by an X-ray background field produced by the cluster members. We aim to quantify the impact of X-ray background fields produced by young embedded clusters on the c… ▽ More The X-ray luminosities of T Tauri stars are about two to four orders of magnitude higher than the luminosity of the contemporary Sun. As these stars are born in clusters, their disks are not only irradiated by their parent star but also by an X-ray background field produced by the cluster members. We aim to quantify the impact of X-ray background fields produced by young embedded clusters on the chemical structure of disks. Further, we want to investigate the importance of the dust for X-ray radiative transfer in disks. We present a new X-ray radiative transfer module for the radiation thermo-chemical disk code ProDiMo, which includes X-ray scattering and absorption by both the gas and dust component. For the X-ray radiative transfer, we consider irradiation by the star and by X-ray background fields. To study the impact of X-rays on the chemical structure of disks we use the well-established disk ionization tracers N2H+ and HCO+. For evolved dust populations, X-ray opacities are mostly dominated by the gas; only for photon energies $E\gtrsim5-10\,$keV, dust opacities become relevant. Consequently, the local disk X-ray radiation field is only affected in dense regions close to the disk midplane. X-ray background fields can dominate the local X-ray disk ionization rate for disk radii $r\gtrsim20\,$au. However, the N2H+ and HCO+ column densities are only significantly affected in case of low cosmic-ray ionization rates, or if the background flux is at least a factor of ten higher than the flux level expected for clusters typical for the solar vicinity. Observable signatures of X-ray background fields in low-mass star-formation regions, like Taurus, are only expected for cluster members experiencing a strong X-ray background field. For the majority of the cluster members, the X-ray background field has only little impact on the disk chemical structure. △ Less

Submitted 20 November, 2017; originally announced November 2017.

Comments: 16 pages, 12 figures, accepted for publication in A&A

Journal ref: A&A 609, A91 (2018)

The chemistry of episodic accretion in embedded objects. 2D radiation thermo-chemical models of the post-burst phase

Authors: Ch. Rab, V. Elbakyan, E. Vorobyov, M. Güdel, O. Dionatos, M. Audard, I. Kamp, W. -F. Thi, P. Woitke, A. Postel

Abstract: Episodic accretion is an important process in the evolution of young stars and their environment. The observed strong luminosity bursts of young stellar objects likely have a long lasting impact on the chemical evolution of the disk and envelope structure. We want to investigate observational signatures of the chemical evolution in the post-burst phase for embedded sources. With such signatures it… ▽ More Episodic accretion is an important process in the evolution of young stars and their environment. The observed strong luminosity bursts of young stellar objects likely have a long lasting impact on the chemical evolution of the disk and envelope structure. We want to investigate observational signatures of the chemical evolution in the post-burst phase for embedded sources. With such signatures it is possible to identify targets that experienced a recent luminosity burst. We present a new model for episodic accretion chemistry based on the 2D, radiation thermo-chemical disk code ProDiMo. We have extended ProDiMo with a proper treatment for envelope structures. For a representative Class I model, we calculated the chemical abundances in the post-burst phase and produced synthetic observables like intensity maps and radial profiles. During a burst many chemical species, like CO, sublimate from the dust surfaces. As the burst ends they freeze out again (post-burst phase). This freeze-out happens from inside-out due to the radial density gradient in the disk and envelope structure. This inside-out freeze-out produces clear observational signatures in spectral line emission, like rings and distinct features in the slope of radial intensity profiles. We fitted synthetic C18O J=2-1 observations with single and two component fits and find that post-burst images are much better matched by the latter. Comparing the quality of such fits allows identification of post-burst targets in a model-independent way. Our models confirm that it is possible to identify post-burst objects from spatially resolved CO observations. However, to derive proper statistics, like frequencies of bursts, from observations it is important to consider aspects like the inclination and structure of the target and also dust properties as those have a significant impact on the freeze-out timescale. △ Less

Submitted 1 June, 2017; v1 submitted 10 May, 2017; originally announced May 2017.

Comments: 17 pages, 6 figures, 2 tables. Accepted for publication in A&A

Journal ref: A&A 604, A15 (2017)

Stellar energetic particle ionization in protoplanetary disks around T Tauri stars

Authors: Ch. Rab, M. Güdel, M. Padovani, I. Kamp, W. -F. Thi, P. Woitke, G. Aresu

Abstract: Anomalies in the abundance measurements of short lived radionuclides in meteorites indicate that the protosolar nebulae was irradiated by a high amount of energetic particles (E$\gtrsim$10 MeV). The particle flux of the contemporary Sun cannot explain these anomalies. However, similar to T Tauri stars the young Sun was more active and probably produced enough high energy particles to explain those… ▽ More Anomalies in the abundance measurements of short lived radionuclides in meteorites indicate that the protosolar nebulae was irradiated by a high amount of energetic particles (E$\gtrsim$10 MeV). The particle flux of the contemporary Sun cannot explain these anomalies. However, similar to T Tauri stars the young Sun was more active and probably produced enough high energy particles to explain those anomalies. We want to study the interaction of stellar energetic particles with the gas component of the disk and identify possible observational tracers of this interaction. We use a 2D radiation thermo-chemical protoplanetary disk code to model a disk representative for T Tauri stars. We use a particle energy distribution derived from solar flare observations and an enhanced stellar particle flux proposed for T Tauri stars. For this particle spectrum we calculate the stellar particle ionization rate throughout the disk with an accurate particle transport model. We study the impact of stellar particles for models with varying X-ray and cosmic-ray ionization rates. We find that stellar particle ionization has a significant impact on the abundances of the common disk ionization tracers HCO$^+$ and N$_2$H$^+$, especially in models with low cosmic-ray ionization rates. In contrast to cosmic rays and X-rays, stellar particles cannot reach the midplane of the disk. Therefore molecular ions residing in the disk surface layers are more affected by stellar particle ionization than molecular ions tracing the cold layers/midplane of the disk. Spatially resolved observations of molecular ions tracing different vertical layers of the disk allow to disentangle the contribution of stellar particle ionization from other competing ionization sources. Modeling such observations with a model like the one presented here allows to constrain the stellar particle flux in disks around T Tauri stars. △ Less

Submitted 24 February, 2017; originally announced February 2017.

Comments: 15 pages, 11 figures, 5 tables. Accepted for publication in A&A

Journal ref: A&A 603, A96 (2017)

Thermochemical modelling of brown dwarf discs

Authors: A. J. Greenwood, I. Kamp, L. B. F. M. Waters, P. Woitke, W. -F. Thi, Ch. Rab, G. Aresu, M. Spaans

Abstract: The physical properties of brown dwarf discs, in terms of their shapes and sizes, are still largely unexplored by observations. 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… ▽ More The physical properties of brown dwarf discs, in terms of their shapes and sizes, are still largely unexplored by observations. 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. △ Less

Submitted 20 February, 2017; v1 submitted 15 February, 2017; originally announced February 2017.

Comments: Accepted for publication in A&A. 11 pages, 10 figures, 3 tables

Journal ref: A&A 601, A44 (2017)

Structure of Herbig AeBe disks at the milliarcsecond scale A statistical survey in the H band using PIONIER-VLTI

Authors: B. Lazareff, J. -P. Berger, J. Kluska, J. -B. Le Bouquin, M. Benisty, F. Malbet, C. Koen, C. Pinte, W. -F. Thi, O. Absil, F. Baron, A. Delboulbé, G. Duvert, A. Isella, L. Jocou, A. Juhasz, S. Kraus, R. Lachaume, F. Ménard, R. Millan-Gabet, J. D. Monnier, T. Moulin, K. Perraut, S. Rochat, F. Soulez , et al. (4 additional authors not shown)

Abstract: Context. It is now generally accepted that the near-infrared excess of Herbig AeBe stars originates in the dust of a circumstellar disk. Aims. The aims of this article are to infer the radial and vertical structure of these disks at scales of order one au, and the properties of the dust grains. Methods. The program objects (51 in total) were observed with the H-band (1.6micron) PIONIER/VLTI interf… ▽ More Context. It is now generally accepted that the near-infrared excess of Herbig AeBe stars originates in the dust of a circumstellar disk. Aims. The aims of this article are to infer the radial and vertical structure of these disks at scales of order one au, and the properties of the dust grains. Methods. The program objects (51 in total) were observed with the H-band (1.6micron) PIONIER/VLTI interferometer. The largest baselines allowed us to resolve (at least partially) structures of a few tenths of an au at typical distances of a few hundred parsecs. Dedicated UBVRIJHK photometric measurements were also obtained. Spectral and 2D geometrical parameters are extracted via fits of a few simple models: ellipsoids and broadened rings with azimuthal modulation. Model bias is mitigated by parallel fits of physical disk models. Sample statistics were evaluated against similar statistics for the physical disk models to infer properties of the sample objects as a group. Results. We find that dust at the inner rim of the disk has a sublimation temperature Tsub~1800K. A ring morphology is confirmed for approximately half the resolved objects; these rings are wide delta_r>=0.5. A wide ring favors a rim that, on the star-facing side, looks more like a knife edge than a doughnut. The data are also compatible with a the combination of a narrow ring and an inner disk of unspecified nature inside the dust sublimation radius. The disk inner part has a thickness z/r~0.2, flaring to z/r~0.5 in the outer part. We confirm the known luminosity-radius relation; a simple physical model is consistent with both the mean luminosity-radius relation and the ring relative width; however, a significant spread around the mean relation is present. In some of the objects we find a halo component, fully resolved at the shortest interferometer spacing, that is related to the HAeBe class. △ Less

Submitted 25 November, 2016; originally announced November 2016.

Comments: 41 pages. 23 figures. Four appendices

Journal ref: A&A 599, A85 (2017)

The mid-IR water and silicate relation in protoplanetary disks

Authors: S. Antonellini, J. Bremer, I. Kamp, P. Riviere-Marichalar, F. Lahuis, W. -F. Thi, P. Woitke, R. Meijerink, G. Aresu, M. Spaans

Abstract: Mid-IR water lines from protoplanetary disks around T Tauri stars have a detection rate of 50\%. Models have identified multiple physical properties of disks such as dust-to-gas mass ratio, dust size power law distribution, disk gas mass, disk inner radius, and disk scale height as potential explanation for the current detection rate. We search for a connection between mid-IR water line fluxes and… ▽ More Mid-IR water lines from protoplanetary disks around T Tauri stars have a detection rate of 50\%. Models have identified multiple physical properties of disks such as dust-to-gas mass ratio, dust size power law distribution, disk gas mass, disk inner radius, and disk scale height as potential explanation for the current detection rate. We search for a connection between mid-IR water line fluxes and the strength of the 10~$μ$m silicate feature. We analyse observed water line fluxes from three blends and compute the 10~$μ$m silicate feature strength from Spitzer spectra. We use a series of published models, exploring disk dust and gas properties, and the effects of different stars. The models also show that the increasing stellar luminosity enhance simultaneously the strength of this dust feature and the water lines fluxes. No correlation is found between the observed mid-IR water lines and the 10~$μ$m silicate. Our sample shows the same difference in the peak strength between amorphous and crystalline silicates that was noted in earlier studies, but our models do not support this intrinsic difference in silicate peak strength. Individual properties of our model series are not able to reproduce the most extreme observations, suggesting that more complex dust properties are required. A parametrized settling prescription is able to boost the peak strength by a factor 2 for the standard model. Water line fluxes are unrelated to the composition of the dust. The pronounced regular trends seen in the model results are washed out in the data due to the larger diversity in stellar and disk properties compared to our model series. The disks with with weaker mid-IR water line fluxes are depleted in gas or enhanced in dust in the inner 10~au. In the case of gas depleted disks, settling produces very strong 10~$μ$m silicate features, with strong peak strength. △ Less

Submitted 31 October, 2016; originally announced October 2016.

Journal ref: A&A 597, A72 (2017)

HCO+ Detection of Dust-Depleted Gas in the Inner Hole of the LkCa 15 Pre-Transitional Disk

Authors: E. Drabek-Maunder, S. Mohanty, J. Greaves, I. Kamp, R. Meijerink, M. Spaans, W. -F. Thi, P. Woitke

Abstract: LkCa 15 is an extensively studied star in the Taurus region known for its pre-transitional disk with a large inner cavity in dust continuum and normal gas accretion rate. The most popular hypothesis to explain the LkCa 15 data invokes one or more planets to carve out the inner cavity, while gas continues to flow across the gap from the outer disk onto the central star. We present spatially unresol… ▽ More LkCa 15 is an extensively studied star in the Taurus region known for its pre-transitional disk with a large inner cavity in dust continuum and normal gas accretion rate. The most popular hypothesis to explain the LkCa 15 data invokes one or more planets to carve out the inner cavity, while gas continues to flow across the gap from the outer disk onto the central star. We present spatially unresolved HCO+ J=4-3 observations of the LkCa 15 disk from the JCMT and model the data with the ProDiMo code. We find that: (1) HCO+ line-wings are clearly detected, certifying the presence of gas in the cavity within <50 AU of the star. (2) Reproducing the observed line-wing flux requires both a significant suppression of cavity dust (by a factor >10^4 compared to the ISM) and a substantial increase in the gas scale-height within the cavity (H_0/R_0 ~ 0.6). An ISM dust-to-gas ratio (d:g=10^-2) yields too little line-wing flux regardless of the scale-height or cavity gas geometry, while a smaller scale-height also under predicts the flux even with a reduced d:g. (3) The cavity gas mass is consistent with the surface density profile of the outer disk extended inwards to the sublimation radius (corresponding to mass M_d ~ 0.03 M_sun), and masses lower by a factor >10 appear to be ruled out. △ Less

Submitted 19 September, 2016; originally announced September 2016.

Comments: ApJ accepted, 46 pages (including the appendix), 18 figures

A proposed new diagnostic for Herbig disc geometry: FWHM versus J of CO ro-vibrational lines

Authors: R. P. Hein Bertelsen, I. Kamp, G. van der Plas, M. E. van den Ancker, L. B. F. M. Waters, W. -F. Thi, P. Woitke

Abstract: Aims. CO ro-vibrational lines observed from Herbig group II discs are seen to be often broad while the same lines observed from group I discs are often narrow. This difference is not well understood. In this paper we explore the underlying cause for this difference and provide a pathway for a better understanding of the geometry and structure of the inner discs around Herbig Ae/Be stars. Methods.… ▽ More Aims. CO ro-vibrational lines observed from Herbig group II discs are seen to be often broad while the same lines observed from group I discs are often narrow. This difference is not well understood. In this paper we explore the underlying cause for this difference and provide a pathway for a better understanding of the geometry and structure of the inner discs around Herbig Ae/Be stars. Methods. High-spectral-resolution infrared spectra of CO ro-vibrational emission from six Herbig Ae/Be candidate stars were taken with CRIRES (CRyogenic high-resolution InfraRed Echelle Spectrograph) at the VLT (Very Large Telescope). From these spectra, we produce individual and co-added CO ro-vibrational line profiles. We investigate line profile shape differences, and we explore the FWHM (Full Width Half Maximum) variations with J quantum number in the context of disc geometry. Furthermore, we put our new sources into the context of earlier observed sources to study a large sample. For comparison, we also investigate the FWHM variations with J of modelled CO ro-vibrational lines from two typical disc geometries produced with the thermo-chemical disc modelling code ProDiMo. Results. - abbreviated - Comparing dust and gas inner disc geometries (inferred by SED and CO ro-vibrational emission) for the expanded sample of observed Herbig discs, we find no clear correspondence between the SED groups of the sources and their inner CO radius. Conclusions. The FWHM versus J is a potential new gas diagnostic for the inner disc with e.g. a constant FWHM versus J indicating the presence of a large gas hole or gap. Both models and observations indicate the potential of this new diagnostic. - abbreviated - △ Less

Submitted 11 March, 2016; originally announced March 2016.

Comments: accepted for publication in A&A

Journal ref: A&A 590, A98 (2016)

Variability in the CO ro-vibrational lines from HD163296

Authors: Rosina P. Hein Bertelsen, I. Kamp, G. van der Plas, M. E. van den Ancker, L. B. F. M. Waters, W. -F. Thi, P. Woitke

Abstract: We present for the first time a direct comparison of multi-epoch (2001-2002 and 2012) CO ro-vibrational emission lines from HD163296. We find that both the line shapes and the FWHM (Full Width Half Maximum) differ between these two epochs. The FWHM of the median observed line profiles are 10-25 km/s larger in the earlier epoch, and confirmed double peaks are only present in high J lines from 2001-… ▽ More We present for the first time a direct comparison of multi-epoch (2001-2002 and 2012) CO ro-vibrational emission lines from HD163296. We find that both the line shapes and the FWHM (Full Width Half Maximum) differ between these two epochs. The FWHM of the median observed line profiles are 10-25 km/s larger in the earlier epoch, and confirmed double peaks are only present in high J lines from 2001-2002. The line wings of individual transitions are similar in the two epochs making an additional central component in the later epoch a likely explanation for the single peaks and the lower FWHM. Variations in NIR brightness have been reported and could be linked to the observed variations. Additionally, we use the thermo chemical disc code ProDiMo to compare for the first time the line shapes, peak separations, FWHM, and line fluxes, to those observed. The ProDiMo model reproduces the peak separations, and low and mid J line fluxes well. The FWHM however, are over predicted and high J line fluxes are under predicted. We propose that a variable non-Keplerian component of the CO ro-vibrational emission, such as a disc wind or an episodic accretion funnel, is causing the difference between the two data sets collected at different epochs, and between model and observations. Additional CO ro-vibrational line detections (with CRIRES/VLT or NIRSPEC/Keck) or [Ne II] line observations with VISIR/VLT could help to clarify the cause of the variability. △ Less

Submitted 16 January, 2016; originally announced January 2016.

Comments: 13 pages, 10 figures, 3 tables. Accepted to MNRAS

Consistent dust and gas models for protoplanetary disks. I. Disk shape, dust settling, opacities, and PAHs

Authors: P. Woitke, M. Min, C. Pinte, W. -F. Thi, I. Kamp, C. Rab, F. Anthonioz, S. Antonellini, C. Baldovin-Saavedra, A. Carmona, C. Dominik, O. Dionatos, J. Greaves, M. Güdel, J. D. Ilee, A. Liebhart, F. Ménard, L. Rigon, L. B. F. M. Waters, G. Aresu, R. Meijerink, M. Spaans

Abstract: We propose a set of standard assumptions for the modelling of Class II and III protoplanetary disks, which includes detailed continuum radiative transfer, thermo-chemical modelling of gas and ice, and line radiative transfer from optical to cm wavelengths. We propose new standard dust opacities for disk models, we present a simplified treatment of PAHs sufficient to reproduce the PAH emission feat… ▽ More We propose a set of standard assumptions for the modelling of Class II and III protoplanetary disks, which includes detailed continuum radiative transfer, thermo-chemical modelling of gas and ice, and line radiative transfer from optical to cm wavelengths. We propose new standard dust opacities for disk models, we present a simplified treatment of PAHs sufficient to reproduce the PAH emission features, and we suggest using a simple treatment of dust settling. We roughly adjust parameters to obtain a model that predicts typical Class II T Tauri star continuum and line observations. We systematically study the impact of each model parameter (disk mass, disk extension and shape, dust settling, dust size and opacity, gas/dust ratio, etc.) on all continuum and line observables, in particular on the SED, mm-slope, continuum visibilities, and emission lines including [OI] 63um, high-J CO lines, (sub-)mm CO isotopologue lines, and CO fundamental ro-vibrational lines. We find that evolved dust properties (large grains) often needed to fit the SED, have important consequences for disk chemistry and heating/cooling balance, leading to stronger emission lines in general. Strong dust settling and missing disk flaring have similar effects on continuum observations, but opposite effects on far-IR gas emission lines. PAH molecules can shield the gas from stellar UV radiation because of their strong absorption and negligible scattering opacities. The observable millimetre-slope of the SED can become significantly more gentle in the case of cold disk midplanes, which we find regularly in our T Tauri models. We propose to use line observations of robust chemical tracers of the gas, such as O, CO, and H2, as additional constraints to determine some key properties of the disks, such as disk shape and mass, opacities, and the dust/gas ratio, by simultaneously fitting continuum and line observations. △ Less

Submitted 11 November, 2015; originally announced November 2015.

Comments: accepted by A&A, 36 pages, 30 figures, 5 tables

Journal ref: A&A 586, A103 (2016)

Mid-IR spectra of Pre-Main Sequence Herbig stars: an explanation for the non-detections of water lines

Authors: S. Antonellini, I. Kamp, F. Lahuis, P. Woitke, W. -F. Thi, R. Meijerink, G. Aresu, M. Spaans, M. Güdel, A. Liebhart

Abstract: 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 TTauri stars. The reason for this is still unclear. In this study, we want to find an explanation for the different detection rates between low mass and high mass pre-main-sequence stars (PMSs) in the mid-IR regime. We run disk models with stellar parameters adjusted t… ▽ More 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 TTauri stars. The reason for this is still unclear. In this study, we want to find an explanation for the different detection rates between low mass and high mass pre-main-sequence stars (PMSs) in the mid-IR regime. We run disk models with stellar parameters adjusted to spectral types B9 through M2, using the radiation thermo-chemical disk modeling code ProDiMo. We produce convolved spectra at the resolution of Spitzer IRS, JWST MIRI and VLT VISIR spectrographs. We apply random noise derived from typical Spitzer spectra for a direct comparison with observations. The strength of the mid-IR water lines correlates directly with the luminosity of the central star. We explored 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. 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. Mid-IR spectra with resolution higher than 20000 are needed to investigate water in protoplanetary disks. Intrinsic differences in disk structure, e.g. inner gaps, gas-to-dust ratio, dust size and distribution, and inner disk scale height, between Herbig and TTauri star disks are able to explain a lower water detection rate in disks around Herbig stars. △ Less

Submitted 27 October, 2015; originally announced October 2015.

Journal ref: A&A 585, A61 (2016)

Understanding the water emission in the mid- and far-IR from protoplanetary disks around T~Tauri stars

Authors: S. Antonellini, I. Kamp, P. Riviere-Marichalar, R. Meijerink, P. Woitke, W. -F. Thi, M. Spaans, G. Aresu, E. K. H. Lee

Abstract: We investigate which properties of protoplanetary disks around T Tauri stars affect the physics and chemistry in the regions where mid- and far-IR water lines originate and their respective line fluxes. We search for diagnostics for future observations. With the code ProDiMo, we build a series of models exploring a large parameter space, computing rotational and rovibrational transitions of water… ▽ More We investigate which properties of protoplanetary disks around T Tauri stars affect the physics and chemistry in the regions where mid- and far-IR water lines originate and their respective line fluxes. We search for diagnostics for future observations. With the code ProDiMo, we build a series of models exploring a large parameter space, computing rotational and rovibrational transitions of water in nonlocal thermodynamic equilibrium (non-LTE). We select a sample of transitions in the mid- IR regime and the fundamental ortho and para water transitions in the far-IR. We investigate the chemistry and the local physical conditions in the line emitting regions. We calculate Spitzer spectra for each model and compare far-IR and mid-IR lines. In addition, we use mid-IR colors to tie the water line predictions to the dust continuum. Parameters affecting the water line fluxes in disks by more than a factor of three are : the disk gas mass, the dust-to-gas mass ratio, the dust maximum grain size, ISM(InterStellarMedium) UV radiation field, the mixing parameter of Dubrulle settling, the disk flaring parameter, and the dust size distribution. The first four parameters affect the mid-IR lines much more than the far-IR lines. A key driver behind water spectroscopy is the dust opacity, which sets the location of the water line emitting region. We identify three types of parameters. Parameters, such as dust-to-gas ratio, ISM radiation field, and dust size distribution, affect the mid-IR lines more, while the far-IR transitions are more affected by the flaring index. The gas mass greatly affects lines in both regimes. Higher spectral resolution and line sensitivities, like from the James Webb Space Telescope, are needed to detect a statistically relevant sample of individual water lines to distinguish further between these types of parameters. △ Less

Submitted 6 October, 2015; originally announced October 2015.

Journal ref: A&A 582, A105 (2015)

Grain size limits derived from 3.6 μm and 4.5 μm coreshine

Authors: J. Steinacker, M. Andersen, W. -F. Thi, R. Paladini, M. Juvela, A. Bacmann, V. -M. Pelkonen, L. Pagani, C. Lefèvre, Th. Henning, A. Noriega-Crespo

Abstract: Recently discovered scattered light from molecular cloud cores in the wavelength range 3-5 μm (called "coreshine") seems to indicate the presence of grains with sizes above 0.5 μm. We aim to analyze 3.6 and 4.5 μm coreshine from molecular cloud cores to probe the largest grains in the size distribution. We analyzed dedicated deep Cycle 9 Spitzer IRAC observations in the 3.6 and 4.5 μm bands for a… ▽ More Recently discovered scattered light from molecular cloud cores in the wavelength range 3-5 μm (called "coreshine") seems to indicate the presence of grains with sizes above 0.5 μm. We aim to analyze 3.6 and 4.5 μm coreshine from molecular cloud cores to probe the largest grains in the size distribution. We analyzed dedicated deep Cycle 9 Spitzer IRAC observations in the 3.6 and 4.5 μm bands for a sample of 10 low-mass cores. We used a new modeling approach based on a combination of ratios of the two background- and foreground-subtracted surface brightnesses and observed limits of the optical depth. The dust grains were modeled as ice-coated silicate and carbonaceous spheres. We discuss the impact of local radiation fields with a spectral slope differing from what is seen in the DIRBE allsky maps. For the cores L260, ecc806, L1262, L1517A, L1512, and L1544, the model reproduces the data with maximum grain sizes around 0.9, 0.5, 0.65, 1.5, 0.6, and > 1.5 μm, respectively. The maximum coreshine intensities of L1506C, L1439, and L1498 in the individual bands require smaller maximum grain sizes than derived from the observed distribution of band ratios. Additional isotropic local radiation fields with a spectral shape differing from the DIRBE map shape do not remove this discrepancy. In the case of Rho Oph 9, we were unable to reliably disentangle the coreshine emission from background variations and the strong local PAH emission. Considering surface brightness ratios in the 3.6 and 4.5 μm bands across a molecular cloud core is an effective method of disentangling the complex interplay of structure and opacities when used in combination with observed limits of the optical depth. △ Less

Submitted 17 August, 2015; originally announced August 2015.

Comments: 23 pages, 18 figures, accepted for publication in A&A

Signatures of warm carbon monoxide in protoplanetary discs observed with Herschel SPIRE

Authors: M. H. D. van der Wiel, D. A. Naylor, I. Kamp, F. Ménard, W. -F. Thi, P. Woitke, G. Olofsson, K. M. Pontoppidan, J. Di Francesco, A. M. Glauser, J. S. Greaves, R. J. Ivison

Abstract: Molecular gas constitutes the dominant mass component of protoplanetary discs. To date, these sources have not been studied comprehensively at the longest far-infrared and shortest submillimetre wavelengths. This paper presents Herschel SPIRE FTS spectroscopic observations toward 18 protoplanetary discs, covering the entire 450-1540 GHz (666-195 $μ$m) range at R~400-1300. The spectra reveal clear… ▽ More Molecular gas constitutes the dominant mass component of protoplanetary discs. To date, these sources have not been studied comprehensively at the longest far-infrared and shortest submillimetre wavelengths. This paper presents Herschel SPIRE FTS spectroscopic observations toward 18 protoplanetary discs, covering the entire 450-1540 GHz (666-195 $μ$m) range at R~400-1300. The spectra reveal clear detections of the dust continuum and, in six targets, a significant amount of spectral line emission primarily attributable to $^{12}$CO rotational lines. Other targets exhibit little to no detectable spectral lines. Low signal-to-noise detections also include signatures from $^{13}$CO, [CI] and HCN. For completeness, we present upper limits of non-detected lines in all targets, including low-energy transitions of H2O and CH$^+$ molecules. The ten $^{12}$CO lines that fall within the SPIRE FTS bands trace energy levels of ~50-500 K. Combined with lower and higher energy lines from the literature, we compare the CO rotational line energy distribution with detailed physical-chemical models, for sources where these are available and published. Our 13CO line detections in the disc around Herbig Be star HD 100546 exceed, by factors of ~10-30, the values predicted by a model that matches a wealth of other observational constraints, including the SPIRE $^{12}$CO ladder. To explain the observed $^{12}$CO/$^{13}$CO ratio, it may be necessary to consider the combined effects of optical depth and isotope selective (photo)chemical processes. Considering the full sample of 18 objects, we find that the strongest line emission is observed in discs around Herbig Ae/Be stars, although not all show line emission. In addition, two of the six T Tauri objects exhibit detectable $^{12}$CO lines in the SPIRE range. △ Less

Submitted 22 August, 2014; originally announced August 2014.

Comments: Accepted for publication in MNRAS. 17 pages, 6 figures

Herschel evidence for disk flattening or gas depletion in transitional disks

Authors: J. T. Keane, I. Pascucci, C. Espaillat, P. Woitke, S. Andrews, I. Kamp, W. -F. Thi, G. Meeus, W. R. F. Dent

Abstract: Transitional disks are protoplanetary disks characterized by reduced near- and mid-infrared emission with respect to full disks. This characteristic spectral energy distribution indicates the presence of an optically thin inner cavity within the dust disk believed to mark the disappearance of the primordial massive disk. We present new Herschel Space Observatory PACS spectra of [OI] 63 micron for… ▽ More Transitional disks are protoplanetary disks characterized by reduced near- and mid-infrared emission with respect to full disks. This characteristic spectral energy distribution indicates the presence of an optically thin inner cavity within the dust disk believed to mark the disappearance of the primordial massive disk. We present new Herschel Space Observatory PACS spectra of [OI] 63 micron for 21 transitional disks. Our survey complements the larger Herschel GASPS program "Gas in Protoplanetary Systems" (Dent et al. 2013) by quadrupling the number of transitional disks observed with PACS at this wavelength. [OI] 63 micron traces material in the outer regions of the disk, beyond the inner cavity of most transitional disks. We find that transitional disks have [OI] 63 micron line luminosities two times fainter than their full disk counterparts. We self consistently determine various stellar properties (e.g. bolometric luminosity, FUV excess, etc.) and disk properties (e.g. disk dust mass, etc.) that could influence the [OI] 63 micron line luminosity and we find no correlations that can explain the lower [OI] 63 micron line luminosities in transitional disks. Using a grid of thermo-chemical protoplanetary disk models, we conclude that either transitional disks are less flared than full disks or they possess lower gas-to-dust ratios due to a depletion of gas mass. This result suggests that transitional disks are more evolved than their full disk counterparts, possibly even at large radii. △ Less

Submitted 2 April, 2014; originally announced April 2014.

Comments: Accepted for publication in ApJ; 52 pages, 16 figures, 8 tables

Disk evolution, element abundances and cloud properties of young gas giant planets

Authors: Ch. Helling, P. Woitke, P. B. Rimmer, I. Kamp, W. -F. Thi, R. Meijerink

Abstract: We discuss the chemical pre-conditions for planet formation, in terms of gas and ice abundances in a protoplanetary disk, as function of time and position, and the resulting chemical composition and cloud properties in the atmosphere when young gas giant planets form. Large deviations between the abundances of the host star and its gas giants seem likely to occur if the planet formation follows th… ▽ More We discuss the chemical pre-conditions for planet formation, in terms of gas and ice abundances in a protoplanetary disk, as function of time and position, and the resulting chemical composition and cloud properties in the atmosphere when young gas giant planets form. Large deviations between the abundances of the host star and its gas giants seem likely to occur if the planet formation follows the core-accretion scenario. These deviations stem from the separate evolution of gas and dust in the disk, where the dust forms the planet cores, followed by the final run-away accretion of the left-over gas. ProDiMo protoplanetary disk models are used to predict the chemical evolution of gas and ice in the midplane. We find that cosmic rays play a crucial role in slowly un-blocking the CO, where the liberated oxygen forms water, which then freezes out quickly. Therefore, the C/O ratio in the gas phase is found to gradually increase with time, in a region bracketed by the water and CO ice-lines. In this regions, C/O is found to approach unity after about 5 Myrs, scaling with the cosmic ray ionisation rate. We then explore how the atmospheric chemistry and cloud properties in young gas giants are affected when the non-solar C/O ratios predicted by the disk models are assumed. The DRIFT cloud formation model is applied to study the formation of atmospheric clouds under the influence of varying primordial element abundances and its feedback onto the local gas. We demonstrate that element depletion by cloud formation plays a crucial role in converting an oxygen-rich atmosphere gas into carbon-rich gas when non-solar, primordial element abundances are considered as suggested by disk model. △ Less

Submitted 18 March, 2014; originally announced March 2014.

Comments: 32 pages, accepted for publication in the OA journal 'Life', Special issue 'Planet Formation and the Rise of Life'. Life (ISSN 2075-1729)

Journal ref: Life 2014, volume 4, issue 2, pages 142-173

[OI] disk emission in the Taurus star forming region

Authors: G. Aresu, I. Kamp, R. Meijerink, M. Spaans, S. Vicente, L. Podio, P. Woitke, F. Menard, W. -F. Thi, M. Gudel, A. Liebhart

Abstract: 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… ▽ More 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.. △ Less

Submitted 11 February, 2014; originally announced February 2014.

Comments: 9 pages, accepted for publication in A&A

Detecting scattered light from low-mass molecular cores at 3.6 $μ$m - Impact of global effects on the observation of coreshine

Authors: J. Steinacker, M. Andersen, W. -F. Thi, A. Bacmann

Abstract: Recently discovered scattered light at 3-5 $μ$m from low-mass cores (so-called "coreshine") reveals the presence of grains around 1 $μ$m, which is larger than the grains found in the low-density interstellar medium. But only about half of the 100+ cores investigated so far show the effect. This prompts further studies on the origin of this detection rate. From the 3D continuum radiative transfer e… ▽ More Recently discovered scattered light at 3-5 $μ$m from low-mass cores (so-called "coreshine") reveals the presence of grains around 1 $μ$m, which is larger than the grains found in the low-density interstellar medium. But only about half of the 100+ cores investigated so far show the effect. This prompts further studies on the origin of this detection rate. From the 3D continuum radiative transfer equation, we derive the expected scattered light intensity from a core placed in an arbitrary direction seen from Earth. We use the approximation of single scattering, consider extinction up to 2nd-order Taylor approximation, and neglect spatial gradients in the dust size distribution. The impact of the directional characteristics of the scattering on the detection of scattered light from cores is calculated for a given grain size distribution, and local effects like additional radiation field components are discussed. The surface brightness profiles of a core with a 1D density profile are calculated for various Galactic locations, and the results are compared to the approximate detection limits. We find that for optically thin radiation and a constant size distribution, a simple limit for detecting scattered light from a low-mass core can be derived that holds for grains with sizes smaller than 0.5 $μ$m. The extinction by the core prohibits detection in bright parts of the Galactic plane, especially near the Galactic center. For scattered light received from low-mass cores with grain sizes beyond 0.5 $μ$m, the directional characteristics of the scattering favors the detection of scattered light above and below the Galactic center, and to some extent near the Galactic anti-center. We identify the local incident radiation field as the major unknown causing deviations from this simple scheme. △ Less

Submitted 7 February, 2014; originally announced February 2014.

Comments: 10 pages, 10 figures, accepted by Astronomy & Astrophysics

Probing the gaseous disk of T Tau N with CN 5-4 lines

Authors: L. Podio, I. Kamp, C. Codella, B. Nisini, G. Aresu, S. Brittain, S. Cabrit, C. Dougados, C. Grady, R. Meijerink, G. Sandell, M. Spaans, W. -F. Thi, G. J. White, P. Woitke

Abstract: We present spectrally resolved Herschel/HIFI observations of the young multiple system T Tau in atomic and molecular lines. While CO, H2O, [C II], and SO lines trace the envelope and the outflowing gas up to velocities of 33 km/s with respect to systemic, the CN 5-4 hyperfine structure lines at 566.7, 566.9 GHz show a narrow double-peaked profile centered at systemic velocity, consistent with an o… ▽ More We present spectrally resolved Herschel/HIFI observations of the young multiple system T Tau in atomic and molecular lines. While CO, H2O, [C II], and SO lines trace the envelope and the outflowing gas up to velocities of 33 km/s with respect to systemic, the CN 5-4 hyperfine structure lines at 566.7, 566.9 GHz show a narrow double-peaked profile centered at systemic velocity, consistent with an origin in the outer region of the compact disk of T Tau N. Disk modeling of the T Tau N disk with the thermo-chemical code ProDiMo produces CN line fluxes and profiles consistent with the observed ones and constrain the size of the gaseous disk (R_out = 110 (+10, -20) AU) and its inclination (i = 25 \pm 5 degree). The model indicates that the CN lines originate in a disk upper layer at 40-110 AU from the star, which is irradiated by the stellar UV field and heated up to temperatures of 50-700 K. With respect to previously observed CN 2-1 millimeter lines, the CN 5-4 lines appear to be less affected by envelope emission, due to their larger critical density and excitation temperature. Hence, high-J CN lines are a unique confusion-free tracer of embedded disks, such as the disk of T Tau N. △ Less

Submitted 30 January, 2014; originally announced January 2014.

Comments: 6 pages, 3 figures, accepted for publication on ApJ Letters

CO ro-vibrational lines in HD100546: A search for disc asymmetries and the role of fluorescence

Authors: R. P. Hein Bertelsen, I. Kamp, M. Goto, G. van der Plas, W. -F. Thi, L. B. F. M. Waters, M. E. van den Ancker, P. Woitke

Abstract: We have studied the emission of CO ro-vibrational lines in the disc around the Herbig Be star HD100546 with the final goal of using these lines as a diagnostic to understand inner disc structure in the context of planet formation. High-resolution IR spectra of CO ro-vibrational emission at eight different position angles were taken with CRIRES at the VLT. From these spectra flux tables, CO ro-vibr… ▽ More We have studied the emission of CO ro-vibrational lines in the disc around the Herbig Be star HD100546 with the final goal of using these lines as a diagnostic to understand inner disc structure in the context of planet formation. High-resolution IR spectra of CO ro-vibrational emission at eight different position angles were taken with CRIRES at the VLT. From these spectra flux tables, CO ro-vibrational line profiles, and population diagrams were produced. We have investigated variations in the line profile shapes and line strengths as a function of slit position angle. We used the thermochemical disc modelling code ProDiMo based on the chemistry, radiation field, and temperature structure of a previously published model for HD100546. Comparing observations and the model, we investigated the possibility of disc asymmetries, the excitation mechanism (UV fluorescence), the geometry, and physical conditions of the inner disc. The observed CO ro-vibrational lines are largely emitted from the inner rim of the outer disc at 10-13 AU. The line shapes are similar for all v levels and line fluxes from all vibrational levels vary only within one order of magnitude. All line profile asymmetries and variations can be explained with a symmetric disc model to which a slit correction and pointing offset is applied. Because the angular size of the CO emitting region (10-13 AU) and the slit width are comparable the line profiles are very sensitive to the placing of the slit. The model reproduces the line shapes and the fluxes of the v=1-0 lines as well as the spatial extent of the CO ro-vibrational emission. It does not reproduce the observed band ratios of 0.5-0.2 with higher vibrational bands. We find that lower gas volume densities at the surface of the inner rim of the outer disc can make the fluorescence pumping more effcient and reproduce the observed band ratios. △ Less

Submitted 21 November, 2013; originally announced November 2013.

Comments: 20 pages, 21 figures

Journal ref: A&A, Volume 561,A102,2014

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

Authors: I. Kamp, W. -F. Thi, G. Meeus, P. Woitke, C. Pinte, R. Meijerink, M. Spaans, I. Pascucci, G. Aresu, W. R. F. Dent

Abstract: 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… ▽ More 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. △ Less

Submitted 8 August, 2013; originally announced August 2013.

Comments: accepted for publication in A&A

Herschel Observations of Gas and Dust in the Unusual 49 Ceti Debris Disk

Authors: A. Roberge, I. Kamp, B. Montesinos, W. R. F. Dent, G. Meeus, J. K. Donaldson, J. Olofsson, A. Moor, J. -C. Augereau, C. Howard, C. Eiroa, W. -F. Thi, D. R. Ardila, G. Sandell, P. Woitke

Abstract: We present far-IR/sub-mm imaging and spectroscopy of 49 Ceti, an unusual circumstellar disk around a nearby young A1V star. The system is famous for showing the dust properties of a debris disk, but the gas properties of a low-mass protoplanetary disk. The data were acquired with the Herschel Space Observatory PACS and SPIRE instruments, largely as part of the "Gas in Protoplanetary Systems" (GASP… ▽ More We present far-IR/sub-mm imaging and spectroscopy of 49 Ceti, an unusual circumstellar disk around a nearby young A1V star. The system is famous for showing the dust properties of a debris disk, but the gas properties of a low-mass protoplanetary disk. The data were acquired with the Herschel Space Observatory PACS and SPIRE instruments, largely as part of the "Gas in Protoplanetary Systems" (GASPS) Open Time Key Programme. Disk dust emission is detected in images at 70, 160, 250, 350, and 500 μm; 49 Cet is significantly extended in the 70 μm image, spatially resolving the outer dust disk for the first time. Spectra covering small wavelength ranges centered on eight atomic and molecular emission lines were obtained, including [OI] 63 μm and [CII] 158 μm. The CII line was detected at the 5σ level - the first detection of atomic emission from the disk. No other emission lines were seen, despite the fact that the OI line is the brightest one observed in Herschel protoplanetary disk spectra (Meeus et al. 2012; Dent et al. 2013). We present an estimate of the amount of circumstellar atomic gas implied by the CII emission. The new far-IR/sub-mm data fills in a large gap in the previous spectral energy distribution (SED) of 49 Cet. A simple model of the new SED confirms the two-component structure of the disk: warm inner dust and cold outer dust that produces most of the observed excess. Finally, we discuss preliminary thermochemical modeling of the 49 Cet gas/dust disk and our attempts to match several observational results simultaneously. Although we are not yet successful in doing so, our investigations shed light on the evolutionary status of the 49 Cet gas, which might not be primordial gas but rather secondary gas coming from comets. △ Less

Submitted 13 May, 2013; originally announced May 2013.

Comments: 11 pages, 4 figures, accepted for publication in The Astrophysical Journal

Journal ref: ApJ 771 (2013) 69

Water vapor in the protoplanetary disk of DG Tau

Authors: L. Podio, I. Kamp, C. Codella, S. Cabrit, B. Nisini, C. Dougados, G. Sandell, J. P. Williams, L. Testi, W. -F. Thi, P. Woitke, R. Meijerink, M. Spaans, G. Aresu, F. Menard, C. Pinte

Abstract: Water is key in the evolution of protoplanetary disks and the formation of comets and icy/water planets. While high excitation water lines originating in the hot inner disk have been detected in several T Tauri stars (TTSs), water vapor from the outer disk, where most of water ice reservoir is stored, was only reported in the closeby TTS TW Hya. We present spectrally resolved Herschel/HIFI observa… ▽ More Water is key in the evolution of protoplanetary disks and the formation of comets and icy/water planets. While high excitation water lines originating in the hot inner disk have been detected in several T Tauri stars (TTSs), water vapor from the outer disk, where most of water ice reservoir is stored, was only reported in the closeby TTS TW Hya. We present spectrally resolved Herschel/HIFI observations of the young TTS DG Tau in the ortho- and para- water ground-state transitions at 557, 1113 GHz. The lines show a narrow double-peaked profile, consistent with an origin in the outer disk, and are ~19-26 times brighter than in TW Hya. In contrast, CO and [C II] lines are dominated by emission from the envelope/outflow, which makes H2O lines a unique tracer of the disk of DG Tau. Disk modeling with the thermo-chemical code ProDiMo indicates that the strong UV field, due to the young age and strong accretion of DG Tau, irradiates a disk upper layer at 10-90 AU from the star, heating it up to temperatures of 600 K and producing the observed bright water lines. The models suggest a disk mass of 0.015-0.1 Msun, consistent with the estimated minimum mass of the solar nebula before planet formation, and a water reservoir of ~1e2-1e3 Earth oceans in vapour, and ~100 times larger in the form of ice. Hence, this detection supports the scenario of ocean delivery on terrestrial planets by impact of icy bodies forming in the outer disk. △ Less

Submitted 6 February, 2013; originally announced February 2013.

Comments: 6 pages, 2 figures, 3 tables, letter accepted for publication on ApJ

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

Authors: R. Meijerink, G. Aresu, I. Kamp, M. Spaans, W. -F. Thi, P. Woitke

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… ▽ 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) △ Less

Submitted 24 August, 2012; originally announced August 2012.

Comments: 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.pdf

An icy Kuiper-Belt around the young solar-type star HD 181327

Authors: J. Lebreton, J. -C. Augereau, W. -F. Thi, A. Roberge, J. Donaldson, G. Schneider, S. T. Maddison, F. Ménard, P. Riviere-Marichalar, G. S. Mathews, I. Kamp, C. Pinte, W. R. F. Dent, D. Barrado, G. Duchêne, J. -F. Gonzalez, C. A. Grady, G. Meeus, E. Pantin, J. P. Williams, P. Woitke

Abstract: HD 181327 is a young F5/F6V star belonging to the Beta Pictoris moving group (12 Myr). It harbors an optically thin belt of circumstellar material at 90 AU. We aim to study the dust properties in the belt in details, and to constrain the gas-to-dust ratio. We obtained far-IR observations with the Herschel/PACS instrument, and 3.2 mm observations with the ATCA array. The geometry of the belt is con… ▽ More HD 181327 is a young F5/F6V star belonging to the Beta Pictoris moving group (12 Myr). It harbors an optically thin belt of circumstellar material at 90 AU. We aim to study the dust properties in the belt in details, and to constrain the gas-to-dust ratio. We obtained far-IR observations with the Herschel/PACS instrument, and 3.2 mm observations with the ATCA array. The geometry of the belt is constrained with newly reduced HST/NICMOS images that break the degeneracy between the disk geometry and the dust properties. We use the radiative transfer code GRaTer to compute a large grid of models, and we identify the grain models that best reproduce the Spectral Energy Distribution through a Bayesian analysis. We attempt to detect the [OI] and [CII] lines with PACS spectroscopy, providing observables to our photochemical code ProDiMo. The HST observations confirm that the dust is confined in a narrow belt. The continuum is detected in the far-IR with PACS and the disk is resolved with both PACS and ATCA. A medium integration of the gas spectral lines only provides upper limits on the line fluxes. We show that the HD 181327 dust disk consists of micron-sized grains of porous amorphous silicates and carbonaceous material surrounded by an important layer of ice, for a total dust mass of 0.05 M\oplus (up to 1 mm). We discuss evidences that the grains are fluffy aggregates. The upper limits on the gas lines do not provide unambiguous constraints: only if the PAH abundance is high, the gas mass must be lower than 17 M\oplus. Despite the weak constraints on the gas disk, the age of HD 181327 and the properties of the dust disk suggest that it has passed the stage of gaseous planets formation. The dust reveals a population of icy planetesimals, similar to the primitive Kuiper Belt, that may be a source for the future delivery of water and volatiles onto forming terrestrial planets. △ Less

Submitted 14 December, 2011; originally announced December 2011.

Comments: 17 pages, 12 figures. Accepted for publication in Astronomy & Astrophysics (2011 december 7)

Gas modelling in the disc of HD 163296

Authors: I. Tilling, P. Woitke, G. Meeus, A. Mora, B. Montesinos, P. Riviere-Marichalar, C. Eiroa, W. -F. Thi, A. Isella, A. Roberge, C. Martin-Zaidi, I. Kamp, C. Pinte, G. Sandell, W. D. Vacca, F. Ménard, I. Mendigutía, G. Duchêne, W. R. F. Dent, G. Aresu, R. Meijerink, M. Spaans

Abstract: We present detailed model fits to observations of the disc around the Herbig Ae star HD 163296. This well-studied object has an age of ~ 4 Myr, with evidence of a circumstellar disc extending out to ~ 540AU. We use the radiation thermo-chemical disc code ProDiMo to model the gas and dust in the circumstellar disc of HD 163296, and attempt to determine the disc properties by fitting to observationa… ▽ More We present detailed model fits to observations of the disc around the Herbig Ae star HD 163296. This well-studied object has an age of ~ 4 Myr, with evidence of a circumstellar disc extending out to ~ 540AU. We use the radiation thermo-chemical disc code ProDiMo to model the gas and dust in the circumstellar disc of HD 163296, and attempt to determine the disc properties by fitting to observational line and continuum data. These include new Herschel/PACS observations obtained as part of the open-time key program GASPS (Gas in Protoplanetary Systems), consisting of a detection of the [OI]63mic line and upper limits for several other far infrared lines. We complement this with continuum data and ground-based observations of the 12CO 3-2, 2-1 and 13CO J=1-0 line transitions, as well as the H2 S(1) transition. We explore the effects of stellar ultraviolet variability and dust settling on the line emission, and on the derived disc properties. Our fitting efforts lead to derived gas/dust ratios in the range 9-100, depending on the assumptions made. We note that the line fluxes are sensitive in general to the degree of dust settling in the disc, with an increase in line flux for settled models. This is most pronounced in lines which are formed in the warm gas in the inner disc, but the low excitation molecular lines are also affected. This has serious implications for attempts to derive the disc gas mass from line observations. We derive fractional PAH abundances between 0.007 and 0.04 relative to ISM levels. Using a stellar and UV excess input spectrum based on a detailed analysis of observations, we find that the all observations are consistent with the previously assumed disc geometry. △ Less

Submitted 10 November, 2011; originally announced November 2011.

Comments: 20 pages, 15 figures. Accepted for publication in A&A

Detection of CH+ emission from the disc around HD100546

Authors: W. -F. Thi, F. Ménard, G. Meeus, C. Martin-Zaïdi, P. Woitke, E. Tatulli, M. Benisty, I. Kamp, I. Pascucci, C. Pinte, C. A. Grady, S. Brittain, G. J. White, C. D. Howard, G. Sandell, C. Eiroa

Abstract: Despite its importance in the thermal-balance of the gas and in the determination of primeval planetary atmospheres, the chemistry in protoplanetary discs remains poorly constrained with only a handful of detected species. We observed the emission from disc around the Herbig Be star HD 100546 with the PACS instrument in the spectroscopic mode on board the Herschel Space Telescope as part of the Ga… ▽ More Despite its importance in the thermal-balance of the gas and in the determination of primeval planetary atmospheres, the chemistry in protoplanetary discs remains poorly constrained with only a handful of detected species. We observed the emission from disc around the Herbig Be star HD 100546 with the PACS instrument in the spectroscopic mode on board the Herschel Space Telescope as part of the Gas in Protoplanetary Systems (GASPS) programme and used archival data from the DIGIT programme to search for the rotational emission of CH+. We detected in both datasets an emission line centred at 72.16 micron that most likely corresponds to the J=5-4 rotational emission of CH+. The J=3-2 and 6-5 transitions are also detected albeit with lower confidence. Other CH+ rotational lines in the PACS observations are blended with water lines. A rotational diagram analysis shows that the CH+ gas is warm at 323 (+2320/-151) K with a mass of 3e-14-5e-12 M_Sun. We modelled the CH+ chemistry with the chemo-physical code ProDiMo using a disc density structure and grain parameters that match continuum observations and near- and mid-infrared interferometric data. The model suggests that CH+ is most abundant at the location of the disc rim at 10-13 AU from the star where the gas is warm, consistent with previous observations of hot CO gas emission. △ Less

Submitted 13 April, 2011; v1 submitted 12 April, 2011; originally announced April 2011.

Comments: 5 pages, 5 figures, A & A letter, accepted, version after language edition