Abstract
The European FP7 project DIANA has performed a coherent analysis of a large set of observational data of 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, which probe the dust, polycyclic aromatic hydrocarbons (PAHs) and the gas in these objects. We define and apply a standardized modeling procedure to fit these data by state-of-the-art modeling codes (
ProDiMo
,
MCFOST
,
MCMax
), solving continuum and line radiative transfer (RT), disk chemistry, and the heating and cooling balance for both the gas and the dust. 3D diagnostic RT tools (e.g., FLiTs) are eventually used to predict all available observations from the same disk model, the DIANA-standard model. Our aim is to determine the physical parameters of the disks, such as total gas and dust masses, the dust properties, the disk shape, and the chemical structure in these disks. 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. It allows us to break some of the degeneracies arising from pure Spectral Energy Distribution (SED) modeling. 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. 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. 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. However, some line fluxes may deviate by a larger factor, and sometimes we find puzzling data which the models cannot reproduce. Some of these issues are 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 to the community through an online database (
http://www.univie.ac.at/diana
).
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. Episodic accretion plays an important role in the evolution of young stars. Although it has been under investigation for a long time, the origin of such episodic accretion events is not yet ...understood. Aims. We investigate the dust and gas emission of a sample of young outbursting sources in the infrared to get a better understanding of their properties and circumstellar material, and we use the results in a further work to model the objects. Methods. We used Herschel data, from our PI program of 12 objects and complemented with archival observations to obtain the spectral energy distributions (SEDs) and spectra of our targets. We report here the main characteristics of our sample, focussing on the SED properties and on the gas emission lines detected in the PACS and SPIRE spectra. Results. The SEDs of our sample show the diversity of the outbursting sources, with several targets showing strong emission in the far-infrared from the embedded objects. Most of our targets reside in a complex environment, which we discuss in detail. We detected several atomic and molecular lines, in particular rotational CO emission from several transitions from J = 38−37 to J = 4−3. We constructed rotational diagrams for the CO lines, and derived in three domains of assumed local thermodynamic equilibrium (LTE) temperatures and column densities, ranging mainly between 0−100 K and 400−500 K. We confirm correlation in our sample between intense CO J = 16−15 emission and the column density of the warm domain of CO, N(warm). We notice a strong increase in luminosity of HH 381 IRS and a weaker increase for PP 13 S, which shows the beginning of an outburst.
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.
Context. Consistent modeling of protoplanetary disks requires the simultaneous solution of both continuum and line radiative transfer, heating and cooling balance between dust and gas and, of course, ...chemistry. Such models depend on panchromatic observations that can provide a complete description of the physical and chemical properties and energy balance of protoplanetary systems. Along these lines, we present a homogeneous, panchromatic collection of data on a sample of 85 T Tauri and Herbig Ae objects for which data cover a range from X-rays to centimeter wavelengths. Datasets consist of photometric measurements, spectra, along with results from the data analysis such as line fluxes from atomic and molecular transitions. Additional properties resulting from modeling of the sources such as disk mass and shape parameters, dust size, and polycyclic aromatic hydrocarbon (PAH) properties are also provided for completeness. Aim. The purpose of this data collection is to provide a solid base that can enable consistent modeling of the properties of protoplanetary disks. To this end, we performed an unbiased collection of publicly available data that were combined to homogeneous datasets adopting consistent criteria. Targets were selected based on both their properties and the availability of data. Methods. Data from more than 50 different telescopes and facilities were retrieved and combined in homogeneous datasets directly from public data archives or after being extracted from more than 100 published articles. X-ray data for a subset of 56 sources represent an exception as they were reduced from scratch and are presented here for the first time. Results. Compiled datasets, along with a subset of continuum and emission-line models are stored in a dedicated database and distributed through a publicly accessible online system. All datasets contain metadata descriptors that allow us to track them back to their original resources. The graphical user interface of the online system allows the user to visually inspect individual objects but also compare between datasets and models. It also offers to the user the possibility to download any of the stored data and metadata for further processing.
Accurate measurements of the physical structure of protoplanetary discs are critical inputs for planet formation models. These constraints are traditionally established via complex modelling of ...continuum and line observations. Instead, we present an empirical framework to locate the CO isotopologue emitting surfaces from high spectral and spatial resolution ALMA observations. We apply this framework to the disc surrounding IM Lupi, where we report the first direct, i.e. model independent, measurements of the radial and vertical gradients of temperature and velocity in a protoplanetary disc. The measured disc structure is consistent with an irradiated self-similar disc structure, where the temperature increases and the velocity decreases towards the disc surface. We also directly map the vertical CO snow line, which is located at about one gas scale height at radii between 150 and 300 au, with a CO freeze-out temperature of 21 ± 2 K. In the outer disc (>300 au), where the gas surface density transitions from a power law to an exponential taper, the velocity rotation field becomes significantly sub-Keplerian, in agreement with the expected steeper pressure gradient. The sub-Keplerian velocities should result in a very efficient inward migration of large dust grains, explaining the lack of millimetre continuum emission outside of 300 au. The sub-Keplerian motions may also be the signature of the base of an externally irradiated photo-evaporative wind. In the same outer region, the measured CO temperature above the snow line decreases to ≈15 K because of the reduced gas density, which can result in a lower CO freeze-out temperature, photo-desorption, or deviations from local thermodynamic equilibrium.
Context. In protoplanetary disks micron-size dust grains coagulate to form larger structures with complex shapes and compositions. The coagulation process changes the absorption and scattering ...properties of particles in the disk in significant ways. To properly interpret observations of protoplanetary disks and to place these observations in the context of the first steps of planet formation, it is crucial to understand the optical properties of these complex structures. Aims. We derive the optical properties of dust aggregates using detailed computations of aggregate structures and compare these computationally demanding results with approximate methods that are cheaper to compute in practice. In this way we wish to understand the merits and problems of approximate methods and define the context in which they can or cannot be used to analyze observations of objects where significant grain growth is taking place. Methods. For the detailed computations we used the discrete dipole approximation (DDA), a method able to compute the interaction of light with a complexly shaped, inhomogeneous particle. We compared the results to those obtained using spherical and irregular, homogeneous and inhomogeneous particles. Results. While no approximate method properly reproduces all characteristics of large dust aggregates, the thermal properties of dust can be analyzed using irregularly shaped, porous, inhomogeneous grains. The asymmetry of the scattering phase function is a good indicator of aggregate size, while the degree of polarization is probably determined by the size of the constituent particles. Optical properties derived from aggregates significantly differ from the most frequently used standard (“astronomical silicate” in spherical grains). We outline a computationally fast and relatively accurate method that can be used for a multiwavelength analysis of aggregate dust in protoplanetary disks.
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.
HD 163296 is a young star surrounded by a planet-forming disk that shows clear signatures of dust gaps and rings; likely an indication of ongoing planet formation. We use the radiation thermochemical ...disk code ProDiMo to investigate the impact of dust/gas gaps on the temperature, chemistry and observables. Furthermore, we model high spatial resolution gas and dust observation of HD 163296 (ALMA/DSHARP). Our first results indicate that features in the observed radial intensity profile of the 12CO line are a consequence of the dust gaps and do not require gas depletion. Those preliminary results indicate that self-consistent modelling of the gas (chemistry, heating/cooling) and dust is necessary to accurately infer the degree of gas depletion within dust gaps. Such information is crucial to understand the processes that generate the disk substructure and their relation to planet formation.
PDF
