Source author record

W. -F. Thi

W. -F. Thi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

36works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

36 published item(s)

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2019arXiv

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

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.

preprint2016arXiv

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

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 -

preprint2016arXiv

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

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.

preprint2016arXiv

The mid-IR water and silicate relation in protoplanetary disks

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.

preprint2016arXiv

Variability in the CO ro-vibrational lines from HD163296

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.

preprint2015arXiv

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

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.

preprint2015arXiv

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

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.

preprint2014arXiv

[OI] disk emission in the Taurus star forming region

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

preprint2014arXiv

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

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.

preprint2014arXiv

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

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.

preprint2014arXiv

Herschel evidence for disk flattening or gas depletion in transitional disks

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.

preprint2014arXiv

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

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.

preprint2014arXiv

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

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.

preprint2013arXiv

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

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.

preprint2013arXiv

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

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.

preprint2013arXiv

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

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

preprint2013arXiv

Water vapor in the protoplanetary disk of DG Tau

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.

preprint2012arXiv

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

Context. Planets are thought to eventually form from the mostly gaseous (~99% of the mass) disks around young stars. The density structure and chemical composition of protoplanetary disks are affected by the incident radiation field at optical, FUV, and X-ray wavelengths, as well as by the dust properties. Aims. The effect of FUV and X-rays on the disk structure and the gas chemical composition are investigated. This work forms the basis of a second paper, which discusses the impact on diagnostic lines of, e.g., C+, O, H2O, and Ne+ observed with facilities such as Spitzer and Herschel. Methods. A grid of 240 models is computed in which the X-ray and FUV luminosity, minimum grain size, dust size distribution, and surface density distribution are varied in a systematic way. The hydrostatic structure and the thermo-chemical structure are calculated using ProDiMo. Results. The abundance structure of neutral oxygen is stable to changes in the X-ray and FUV luminosity, and the emission lines will thus be useful tracers of the disk mass and temperature. The C+ abundance distribution is sensitive to both X-rays and FUV. The radial column density profile shows two peaks, one at the inner rim and a second one at a radius r=5-10 AU. Ne+ and other heavy elements have a very strong response to X-rays, and the column density in the inner disk increases by two orders of magnitude from the lowest (LX = 1e29 erg/s) to the highest considered X-ray flux (LX = 1e32 erg/s). FUV confines the Ne+ ionized region to areas closer to the star at low X-ray luminosities (LX = 1e29 erg/s). H2O abundances are enhanced by X-rays due to higher temperatures in the inner disk and higher ionization fractions in the outer disk. The line fluxes and profiles are affected by the effects on these species, thus providing diagnostic value in the study of FUV and X-ray irradiated disks around T Tauri stars. (abridged)

preprint2011arXiv

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

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.

preprint2011arXiv

Constraining the structure of the planet-forming region in the disk of the Herbig Be star HD 100546

Studying the physical conditions in circumstellar disks is a crucial step toward understanding planet formation. Of particular interest is the case of HD 100546, a Herbig Be star that presents a gap within the first 13 AU of its protoplanetary disk, that may originate in the dynamical interactions of a forming planet. We gathered a large amount of new interferometric data using the AMBER/VLTI instrument in the H- and K-bands to spatially resolve the warm inner disk and constrain its structure. Then, combining these measurements with photometric observations, we analyze the circumstellar environment of HD 100546 in the light of a passive disk model based on 3D Monte-Carlo radiative transfer. Finally, we use hydrodynamical simulations of gap formation by planets to predict the radial surface density profile of the disk and test the hypothesis of ongoing planet formation. The SED and the NIR interferometric data are adequately reproduced by our model. We show that the H- and K-band emissions are coming mostly from the inner edge of the internal dust disk, located near 0.24 AU from the star, i.e., at the dust sublimation radius in our model. We directly measure an inclination of $33^{\circ} \pm 11^{\circ}$ and a position angle of $140^{\circ} \pm 16^{\circ}$ for the inner disk. This is similar to the values found for the outer disk ($i \simeq 42^{\circ}$, $PA \simeq 145^{\circ}$), suggesting that both disks may be coplanar. We finally show that 1 to 8 Jupiter mass planets located at $\sim 8$ AU from the star would have enough time to create the gap and the required surface density jump of three orders of magnitude between the inner and outer disk. However, no information on the amount of matter left in the gap is available, which precludes us from setting precise limits on the planet mass, for now.

preprint2011arXiv

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

Aims. We want to understand the chemistry and physics of disks on the basis of a large unbiased and statistically relevant grid of disk models. One of the main goals is to explore the diagnostic power of various gas emission lines and line ratios for deriving main disk parameters such as the gas mass. Methods. We explore the results of the DENT grid (Disk Evolution with Neat Theory) that consists of 300 000 disk models with 11 free parameters. Through a statistical analysis, we search for correlations and trends in an effort to find tools for disk diagnostic. Results. All calculated quantities like species masses, temperatures, continuum and line fluxes differ by several orders of magnitude across the entire parameter space. The broad distribution of these quantities as a function of input parameters shows the limitation of using a prototype T Tauri or Herbig Ae/Be disk model. The statistical analysis of the DENT grid shows that CO gas is rarely the dominant carbon reservoir in disks. Models with large inner radii (10 times the dust condensation radius) and/or shallow surface density gradients lack massive gas phase water reservoirs. Also, 60% of the disks have gas temperatures averaged over the oxygen mass in the range between 15 and 70 K; the average gas temperatures for CO and O differ by less than a factor two. Studying the observational diagnostics, the [CII] 158 \mum fine structure line flux is very sensitive to the stellar UV flux and presence of a UV excess and it traces the outer disk radius (Rout). In the submm, the CO low J rotational lines also trace Rout. Low [OI] 63/145 line ratios (< a few) can be explained with cool atomic O gas in the uppermost surface layers leading to self-absorption in the 63 \mum line; this occurs mostly for massive non-flaring, settled disk models without UV excess. ... abbreviated

preprint2011arXiv

Detection of CH+ emission from the disc around HD100546

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.

preprint2011arXiv

Gas modelling in the disc of HD 163296

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.

preprint2011arXiv

The unusual protoplanetary disk around the T Tauri star ET Cha

We present new continuum and line observations, along with modelling, of the faint (6-8) Myr old T Tauri star ET Cha belonging to the eta Chamaeleontis cluster. We have acquired HERSCHEL/PACS photometric fluxes at 70 mic and 160 mic, as well as a detection of the [OI] 63 mic fine-structure line in emission, and derived upper limits for some other far-IR OI, CII, CO and o-H2O lines. The HERSCHEL data is complemented by new ANDICAM B-K photometry, new HST/COS and HST/STIS UV-observations, a non-detection of CO J=3-2 with APEX, re-analysis of a UCLES high-resolution optical spectrum showing forbidden emission lines like [OI] 6300A, [SII] 6731A and 6716A, and [NII] 6583A, and a compilation of existing broad-band photometric data. We used the thermo-chemical disk code ProDiMo and the Monte-Carlo radiative transfer code MCFOST to model the protoplanetary disk around ET Cha. Based on these models we can determine the disk dust mass Mdust = (2.E-8 - 5.E-8) Msun, whereas the total disk gas mass is found to be only little constrained, Mgas = (5.E-5 - 3.E-3) Msun. In the models, the disk extends from 0.022 AU (just outside of the co-rotation radius) to only about 10 AU. Larger disks are found to be inconsistent with the CO J=3-2 non-detection. The low velocity component of the [OI] 6300A emission line is consistent with being emitted from the inner disk. The model can also reproduce the line flux of H2 v=1-0 S(1) at 2.122 mic. An additional high-velocity component of the [OI] 6300A emission line, however, points to the existence of an additional jet/outflow of low velocity (40 - 65) km/s with mass loss rate ~1.E-9 Msun/yr. In relation to our low estimations of the disk mass, such a mass loss rate suggests a disk lifetime of only ~(0.05 - 3) Myr, substantially shorter than the cluster age. The evolutionary state of this unusual protoplanetary disk is discussed.

preprint2011arXiv

Water depletion in the disk atmosphere of Herbig AeBe stars

We present high resolution (R = 100,000) L-band spectroscopy of 11 Herbig AeBe stars with circumstellar disks. The observations were obtained with the VLT/CRIRES to detect hot water and hydroxyl radical emission lines previously detected in disks around T Tauri stars. OH emission lines are detected towards 4 disks. The OH P4.5 (1+,1-) doublet is spectrally resolved as well as the velocity profile of each component of the doublet. Its characteristic double-peak profile demonstrates that the gas is in Keplerian rotation and points to an emitting region extending out to ~ 15-30 AU. The OH, emission correlates with disk geometry as it is mostly detected towards flaring disks. None of the Herbig stars analyzed here show evidence of hot water vapor at a sensitivity similar to that of the OH lines. The non-detection of hot water vapor emission indicates that the atmosphere of disks around Herbig AeBe stars are depleted of water molecules. Assuming LTE and optically thin emission we derive a lower limit to the OH/H2O column density ratio > 1 - 25 in contrast to T Tauri disks for which the column density ratio is 0.3 -- 0.4.

preprint2010arXiv

Evidence for episodic warm outflowing CO gas from the intermediate mass young stellar object IRAS 08470-4321

We present a R=10,000 M-band spectrum of LLN19 (IRAS 08470-4321), a heavily embedded intermediate-mass young stellar object located in the Vela Molecular Cloud, obtained with VLT-ISAAC. The data were fitted by a 2-slab cold-hot model and a wind model. The spectrum exhibits deep broad ro-vibrational absorption lines of 12CO v=1<-0 and 13CO v=1<-0. A weak CO ice feature at 4.67 micron is also detected. Differences in velocity indicate that the warm gas is distinct from the cold millimeter emitting gas, which may be associated with the absorption by cooler gas (45K). The outflowing warm gas at 300-400K and with a mass-loss rate varying between 0.48E-7 and 4.2E-7 MSun /yr can explain most of the absorption. Several absorption lines were spectrally resolved in subsequent spectra obtained with the VLT-CRIRES instrument. Multiple absorption substructures in the high-resolution (R=100,000) spectra indicate that the mass-loss is episodic with at least two major events that occurred recently (<28 years). The discrete mass-loss events together with the large turbulent width of the gas (dv=10-12 km/s) are consistent with the predictions of the Jet-Bow shock outflow and the wide-angle wind model. The CO gas/solid column density ratio of 20-100 in the line-of-sight confirms that the circumstellar environment of LLN~19 is warm. We also derive a 12C/13C ratio of 67 +/- 3, consistent with previous measurements in local molecular clouds but not with the higher ratios found in the envelope of other young stellar objects.

preprint2010arXiv

The Herschel view of Gas in Protoplanetary Systems (GASPS). First comparisons with a large grid of models

The Herschel GASPS Key Program is a survey of the gas phase of protoplanetary discs, targeting 240 objects which cover a large range of ages, spectral types, and disc properties. To interpret this large quantity of data and initiate self-consistent analyses of the gas and dust properties of protoplanetary discs, we have combined the capabilities of the radiative transfer code MCFOST with the gas thermal balance and chemistry code ProDiMo to compute a grid of 300 000 disc models (DENT). We present a comparison of the first Herschel/GASPS line and continuum data with the predictions from the DENT grid of models. Our objective is to test some of the main trends already identified in the DENT grid, as well as to define better empirical diagnostics to estimate the total gas mass of protoplanetary discs. Photospheric UV radiation appears to be the dominant gas-heating mechanism for Herbig stars, whereas UV excess and/or X-rays emission dominates for T Tauri stars. The DENT grid reveals the complexity in the analysis of far-IR lines and the difficulty to invert these observations into physical quantities. The combination of Herschel line observations with continuum data and/or with rotational lines in the (sub-)millimetre regime, in particular CO lines, is required for a detailed characterisation of the physical and chemical properties of circumstellar discs.

preprint2010arXiv

X-ray impact on the protoplanetary disks around T Tauri stars

Context: T Tauri stars have X-ray luminosities ranging from L_X = 10^28-10^32 erg/s. These luminosities are similar to UV luminosities (L_UV 10^30-10^31 erg/s) and therefore X-rays are expected to affect the physics and chemistry of the upper layers of their surrounding protoplanetary disks. Aim: The effects and importance of X-rays on the chemical and hydrostatic structure of protoplanetary disks are investigated, species tracing X-ray irradiation (for L_X >= 10^29 erg/s) are identified and predictions for [OI], [CII] and [NII] fine structure line fluxes are provided. Methods: We have implemented X-ray physics and chemistry into the chemo-physical disk code ProDiMo. We include Coulomb heating and H2 ionization as heating processes and primary and secondary ionization due to X-rays in the chemistry. Results: X-rays heat up the gas causing it to expand in the optically thin surface layers. Neutral molecular species are not much affected in their abundance and spatial distribution, but charged species such as N+, OH+, H2O+ and H3O+ show enhanced abundances in the disk surface. Conclusions: Coulomb heating by X-rays changes the vertical structure of the disk, yielding temperatures of ~ 8000 K out to distances of 50 AU. The chemical structure is altered by the high electron abundance in the gas in the disk surface, causing an efficient ion-molecule chemistry. The products of this, OH+, H2O+ and H3O+, are of great interest for observations of low-mass young stellar objects with the Herschel Space Observatory. [OI] (at 63 and 145 mic) and [CII] (at 158 mic) fine structure emission are only affected for L_X > 10^30 erg/s.

preprint2009arXiv

Radiation thermo-chemical models of protoplanetary disks II. Line diagnostics

Methods. We use the recently developed disk code ProDiMo to calculate the physico-chemical structure of protoplanetary disks and apply the Monte-Carlo line radiative transfer code RATRAN to predict observable line profiles and fluxes. We consider a series of Herbig Ae type disk models ranging from 10^-6 M_Sun to 2.2 10^-2 M_Sun (between 0.5 and 700 AU) to discuss the dependency of the line fluxes and ratios on disk mass for otherwise fixed disk parameters. Results. We find the [CII] 157.7 mum line to originate in LTE from the surface layers of the disk, where Tg > Td . The total emission is dominated by surface area and hence depends strongly on disk outer radius. The [OI] lines can be very bright (> 10^-16 W/m^2) and form in slightly deeper and closer regions under non-LTE conditions. The high-excitation [OI] 145.5 mum line, which has a larger critical density, decreases more rapidly with disk mass than the 63.2 mum line. Therefore, the [OI] 63.2 mum/145.5 mum ratio is a promising disk mass indicator, especially as it is independent of disk outer radius for Rout > 200 AU. CO is abundant only in deeper layers A_V >~ 0.05. For too low disk masses (M_disk <~10^-4 M_Sun) the dust starts to become transparent, and CO is almost completely photo-dissociated. For masses larger than that the lines are an excellent independent tracer of disk outer radius and can break the outer radius degeneracy in the [OI] 63.2 mum/[CII]157.7 mum line ratio. Conclusions. The far-IR fine-structure lines of [CII] and [OI] observable with Herschel provide a promising tool to measure the disk gas mass, although they are mainly generated in the atomic surface layers. In spatially unresolved observations, none of these lines carry much information about the inner, possibly hot regions < 30 AU.

preprint2007arXiv

A PDR-Code Comparison Study

We present a comparison between independent computer codes, modeling the physics and chemistry of interstellar photon dominated regions (PDRs). Our goal was to understand the mutual differences in the PDR codes and their effects on the physical and chemical structure of the model clouds, and to converge the output of different codes to a common solution. A number of benchmark models have been created, covering low and high gas densities and far ultraviolet intensities. The benchmark models were computed in two ways: one set assuming constant temperatures, thus testing the consistency of the chemical network and photo-processes, and a second set determining the temperature selfconsistently. We investigated the impact of PDR geometry and agreed on the comparison of results from spherical and plane-parallel PDR models. We identified a number of key processes governing the chemical network which have been treated differently in the various codes, and defined a proper common treatment. We established a comprehensive set of reference models for ongoing and future PDR model bench-marking and were able to increase the agreement in model predictions for all benchmark models significantly.

preprint2006arXiv

Detectability of dirty dust grains in brown dwarf atmospheres

Dust clouds influence the atmospheric structure of brown dwarfs, and they affect the heat transfer and change the gas-phase chemistry. However, the physics of their formation and evolution is not well understood. In this letter, we predict dust signatures and propose a potential observational test of the physics of dust formation in brown dwarf atmosphere based on the spectral features of the different solid components predicted by dust formation theory. A momentum method for the formation of dirty dust grains (nucleation, growth, evaporation, drift) is used in application to a static brown dwarf atmosphere structure to compute the dust grain properties, in particular the heterogeneous grain composition and the grain size. Effective medium and Mie theory are used to compute the extinction of these spherical grains. Dust formation results in grains whose composition differs from that of grains formed at equilibrium. Our kinetic model predicts that solid amorphous SiO2[s] (silica) is one of the most abundant solid component followed by amorphous MgSiO4[s] and MgSiO3[s], while SiO2[s] is absent in equilibrium models because it is a metastable solid. Solid amorphous SiO2[s] possesses a strong broad absorption feature centered at 8.7mum, while amorphous Mg2SiO4[s]/MgSiO3[s] absorb at 9.7mum beside other absorption features at longer wavelength. Those features at lambda < 15mum are detectable in absorption if grains are small (radius < 0.2mum) in the upper atmosphere as suggested by our model. We suggest that the detection of a feature at 8.7mum in deep infrared spectra could provide evidence for non-equilibrium dust formation that yields grains composed of metastable solids in brown dwarf atmospheres. This feature will shift towards 10mum and broaden if silicates (e.g. fosterite) are much more abundant.