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 1 au, and the properties of the dust grains. Methods. The program objects (51 in total) were observed with the H-band (1.6 mu m) 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 T sub(sub)approximate 1800 K. A ring morphology is confirmed for approximately half the resolved objects; these rings are wide delta r/r> or = 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 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/rapproximate 0.2, flaring to z/rapproximate 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.
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. The first paper of this series focuses on the assumptions about the shape of the disk, the dust opacities, dust settling, and polycyclic aromatic hydrocarbons (PAHs). In particular, we propose new standard dust opacities for disk models, we present a simplified treatment of PAHs in radiative equilibrium which is sufficient to reproduce the PAH emission features, and we suggest using a simple yet physically justified treatment of dust settling. We propose to use line observations of robust chemical tracers of the gas, such as O, CO, and H2, as additional constraints to determine a number of key properties of the disks, such as disk shape and mass, opacities, and the dust/gas ratio, by simultaneously fitting continuum and line observations.
Context.
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.
Aims.
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 the derived gas properties.
Methods.
We applied the radiation thermo-chemical disk code P
RO
D
I
M
O
(PROtoplanetary DIsk MOdel) to model the dust and gas disk of HD 163296 self-consistently, using the DSHARP (Disk Substructure at High Angular Resolution) gas and dust observations. With this model we investigated the impact of dust gaps and gas gaps on the observables and the derived gas properties, considering chemistry, and heating and cooling processes.
Results.
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.
Conclusions.
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.
We introduce a new modelling framework including the Fast Line Tracer (FLITS) to simulate infrared line emission spectra from protoplanetary discs. This paper focusses on the mid-IR spectral region ...between 9.7 and 40 μm 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 S I, S II, S III, Si II, Fe II, Ne II, Ne III, Ar II, and Ar III. 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, for example, due to disc–planet interactions. These walls are irradiated and heated by the star which causes the molecules to emit strongly in the mid-IR. The molecules in the walls cannot be photodissociated easily by UV because of the large densities in the walls favouring their re-formation. Most observable molecular emission lines are found to be optically thick. An abundance analysis is hence not straightforward, and the results of simple slab models concerning molecular column densities can be misleading. We find that the difference between gas and dust temperatures in the disc surface is important for the line formation. The mid-IR emission features of different molecules probe the gas temperature at different depths in the disc, along the following sequence: OH (highest)–CO–H2O and CO2–HCN–C2H2 (deepest), just where these molecules start to become abundant. 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 James Webb Space Telescope (JWST)/MIRI observations, and to possibly detect secondary illuminated disc walls based on their specific mid-IR molecular signature.
Context. 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 ≳ 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. Aims. We aim to study the interaction of stellar energetic particles with the gas component of the disk (i.e. ionization of molecular hydrogen) and identify possible observational tracers of this interaction. Methods. We used a 2D radiation thermo-chemical protoplanetary disk code to model a disk representative for T Tauri stars. We used 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 calculated the stellar particle ionization rate throughout the disk with an accurate particle transport model. We studied the impact of stellar particles for models with varying X-ray and cosmic-ray ionization rates. Results. We find that stellar particle ionization has a significant impact on the abundances of the common disk ionization tracers HCO+ and N2H+, especially in models with low cosmic-ray ionization rates (e.g. 10-19 s-1 for molecular hydrogen). 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 and midplane of the disk. Conclusions. 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. Modelling such observations with a model like the one presented here allows to constrain the stellar particle flux in disks around T Tauri stars.
Aims: We aim to define a small and large chemical network which can be used for the quantitative simultaneous analysis of molecular emission from the near-IR to the submm. We also aim to revise ...reactions of excited molecular hydrogen, which are not included in UMIST, to provide a homogeneous database for future applications. Methods: We have used the thermo-chemical disk modeling code ProDiMo and a standard T Tauri disk model to evaluate the impact of various chemical networks, reaction rate databases and sets of adsorption energies on a large sample of chemical species and emerging line fluxes from the near-IR to the submm wavelength range. Results: We find large differences in the masses and radial distribution of ice reservoirs when considering freeze-out on bare or polar ice coated grains. Most strongly the ammonia ice mass and the location of the snow line (water) change. As a consequence molecules associated to the ice lines such as N2H+ change their emitting region; none of the line fluxes in the sample considered here changes by more than 25% except CO isotopologues, CN and N2H+ lines. The three-body reaction N+H2+M plays a key role in the formation of water in the outer disk. Besides that, differences between the UMIST 2006 and 2012 database change line fluxes in the sample considered here by less than a factor of two (a subset of low excitation CO and fine structure lines stays even within 25%); exceptions are OH, CN, HCN, HCO+ and N2H+ lines. However, different networks such as OSU and KIDA 2011 lead to pronounced differences in the chemistry inside 100 au and thus affect emission lines from high excitation CO, OH and CN lines. H2 is easily excited at the disk surface and state-to-state reactions enhance the abundance of CH+ and to a lesser extent HCO+. For sub-mm lines of HCN, N2H+ and HCO+, a more complex larger network is recommended. Conclusions: More work is required to consolidate data on key reactions leading to the formation of water, molecular ions such as HCO+ and N2H+ as well as the nitrogen chemistry. This affects many of the key lines used in the interpretation of disk observations. Differential analysis of various disk models using the same chemical input data will be more robust than the interpretation of absolute fluxes.
Aims. Brown dwarfs are covered by dust cloud layers which cause inhomogeneous surface features and move below the observable tau = 1 level during the object's evolution. The cloud layers have a ...strong influence on the structure and spectral appearance of brown dwarfs and extra-solar planets, e.g. by providing high local opacities and by removing condensable elements from the atmosphere causing a sub-solar metalicity in the atmosphere. We aim at understanding the formation of cloud layers in quasi-static substellar atmospheres that consist of dirty grains composed of numerous small islands of different solid condensates. Methods. The time-dependent description is a kinetic model describing nucleation, growth and evaporation. It is extended to treat gravitational settling and is applied to the static-stationary case of substellar model atmospheres. From the solution of the dust moments, we determine the grain size distribution function approximately which, together with the calculated material volume fractions, provides the basis for applying effective medium theory and Mie theory to calculate the opacities of the composite dust grains. Results. The cloud particles in brown dwarfs and hot giant-gas planets are found to be small in the high atmospheric layers (a approx 0.01 mum), and are composed of a rich mixture of all considered condensates, in particular MgSiO sub(3) s, Mg sub(2) SiO sub(4) s and SiO sub(2) s. As the particles settle downward, they increase in size and reach several 100 mum in the deepest layers. The more volatile parts of the grains evaporate and the particles stepwise purify to form composite particles of high-temperature condensates in the deeper layers, mainly made of Fes and Al sub(2) O sub(3) s. The gas phase abundances of the elements involved in the dust formation process vary by orders of magnitudes throughout the atmosphere. The grain size distribution is found to be relatively broad in the upper atmospheric layers but strongly peaked in the deeper layers. This reflects the cessation of the nucleation process at intermediate heights. The spectral appearance of the cloud layers in the mid IR (7-20 mum) is close to a grey body with only weak broad features of a few percent, mainly caused by MgSiO sub(3) s, and Mg sub(2) SiO sub(4) s. These features are, nevertheless, a fingerprint of the dust in the higher atmospheric layers that can be probed by observations. Conclusions. Our models predict that the gas phase depletion is much weaker than phase-equilibrium calculations in the high atmospheric layers. Because of the low densities, the dust formation process is incomplete there, which results in considerable amounts of left-over elements that might produce stronger and broader neutral metallic lines.
We study the hydrostatic density structure of the inner disc rim around Herbig Ae stars using the thermo-chemical hydrostatic code prodimo. We compare the spectral energy distributions (SEDs) and ...images from our hydrostatic disc models to that from prescribed density structure discs. The 2D continuum radiative transfer in prodimo includes isotropic scattering. The dust temperature is set by the condition of radiative equilibrium. In the thermal-decoupled case, the gas temperature is governed by the balance between various heating and cooling processes. The gas and dust interact thermally via photoelectrons, radiatively, and via gas accommodation on grain surfaces. As a result, the gas is much hotter than in the thermo-coupled case, where the gas and dust temperatures are equal, reaching a few thousands K in the upper disc layers and making the inner rim higher. A physically motivated density drop at the inner radius ('soft-edge') results in rounded inner rims, which appear ring-like in near-infrared images. The combination of lower gravity pull and hot gas beyond ∼1 au results in a disc atmosphere that reaches a height over radius ratio z/r of 0.2, while this ratio is 0.1 only in the thermo-coupled case. This puffed-up disc atmosphere intercepts larger amount of stellar radiation, which translates into enhanced continuum emission in the 3-30 μm wavelength region from hotter grains at ∼500 K. We also consider the effect of disc mass and grain size distribution on the SEDs self-consistently feeding those quantities back into the gas temperature, chemistry and hydrostatic equilibrium computation.
PDF
