Aims.It is commonly assumed that the massive winds of AGB stars are dust-driven and pulsation-enhanced. However, detailed frequency-dependent dynamical models that can explain the observed magnitudes ...of mass loss rates and outflow velocities have been published so far only for C-stars. This letter reports on first results of similar models for oxygen-rich AGB stars. The aim is to provide a better understanding of the wind driving mechanism, the dust condensation sequence, and the role of pulsations. Methods.New dynamical models for dust-driven winds of oxygen-rich AGB stars are presented which include frequency-dependent Monte Carlo radiative transfer by means of a sparse opacity distribution technique and a time-dependent treatment of the nucleation, growth and evaporation of inhomogeneous dust grains composed of a mixture of Mg2SiO4, SiO2, Al2O3, TiO2, and solid Fe. Results.The frequency-dependent treatment of radiative transfer reveals that the gas is cold close to the star ($700{-}900\,$K at $1.5{-}2\,R_\star$) which facilitates the nucleation process. The dust temperatures are strongly material-dependent, with differences as large as 1000 K for different pure materials, which has an important influence on the dust formation sequence. Two dust layers are formed in the dynamical models: almost pure glassy Al2O3 close to the star ($r \ga 1.5\,R_\star$) and the more opaque Fe-poor Mg-Fe-silicates further out. Solid Fe and Fe-rich silicates are found to be the only condensates that can efficiently absorb the stellar light in the near IR. Consequently, they play a key role in the wind driving mechanism and act as a thermostat. Only small amounts of Fe can be incorporated into the grains, because otherwise the grains become too hot. Thus, the models reveal almost no mass loss, and no dust shells. Conclusions.The observed dust sequence Al2O3 $\to$ Fe-poor Mg-Fe-silicates for oxygen-rich AGB stars having low $\to$ high mass loss rates is in agreement with the presented model and can be understood as follows: Al2O3 is present in the extended atmosphere of the star below the wind acceleration region, also without mass loss. The Mg-Fe-silicates form further out and, therefore, their amount depends on the mass loss rate. The driving mechanism of oxygen-rich AGB stars is still an unsolved puzzle.
New axisymmetric (2D) models for dust-driven winds of C-stars are presented which include hydrodynamics with radiation pressure on dust, equilibrium chemistry and time-dependent dust formation with ...coupled grey Monte Carlo radiative transfer. Considering the simplest case without stellar pulsation (hydrostatic inner boundary condition) these models reveal a more complex picture of the dust formation and wind acceleration than earlier published spherically symmetric (1D) models. The so-called exterior κ-mechanism causes radial oscillations with short phases of active dust formation between longer phases without appreciable dust formation, just like in the 1D models. However, in 2D geometry, the oscillations can be out-of-phase at different places above the stellar atmosphere which results in the formation of dust arcs or smaller caps that only occupy a certain fraction of the total solid angle. These dust structures are accelerated outward by radiation pressure, expanding radially and tangentially like mushroom clouds, while dust-poor matter is falling back towards the star at other places. A highly dynamical and turbulent dust formation zone is created in this way, which again leads to inhomogeneous dust production. Further away from the star, flow instabilities (e.g. Rayleigh-Taylor) have time to fragment the outward moving arcs and shells to produce numerous small-scale cloud-like sub-structures.
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. WASP-18b is an ultra-hot Jupiter with a temperature difference of up to 2500 K between day and night. Such giant planets begin to emerge as a planetary laboratory for understanding cloud ...formation and gas chemistry in well-tested parameter regimes in order to better understand planetary mass loss and for linking observed element ratios to planet formation and evolution. Aims. We aim to understand where clouds form, their interaction with the gas-phase chemistry through depletion and enrichment, the ionisation of the atmospheric gas, and the possible emergence of an ionosphere on ultra-hot Jupiters. Methods. We used 1D profiles from a 3D atmosphere simulation for WASP-18b as input for kinetic cloud formation and gas-phase chemical equilibrium calculations. We solved our kinetic cloud formation model for these 1D profiles, which sample the atmosphere of WASP-18b at 16 different locations along the equator and in the mid-latitudes. We derived the gas-phase composition consistently. Results. The dayside of WASP-18b emerges as completely cloud-free as a result of the very high atmospheric temperatures. In contrast, the nightside is covered in geometrically extended and chemically heterogeneous clouds with dispersed particle size distributions. The atmospheric C/O ratio increases to >0.7 and the enrichment of the atmospheric gas with cloud particles is ρd/ρgas > 10−3. The clouds that form at the limbs appear located farther inside the atmosphere, and they are the least extended. Not all day- to nightside terminator regions form clouds. The gas phase is dominated by H2, CO, SiO, H2O, H2S, CH4, and SiS. In addition, the dayside has a substantial degree of ionisation that is due to ions such as Na+, K+, Ca+, and Fe+. Al+ and Ti+ are the most abundant of their element classes. We find that WASP-18b, as one example for ultra-hot Jupiters, develops an ionosphere on the dayside.
Dynamic mineral clouds on HD 189733b Lee, G; Dobbs-Dixon, I; Helling, Ch ...
Astronomy and astrophysics (Berlin),
10/2016, Volume:
594
Journal Article
Peer reviewed
Context. Observations of exoplanet atmospheres have revealed the presence of cloud particles in their atmospheres. 3D modelling of cloud formation in atmospheres of extrasolar planets coupled to the ...atmospheric dynamics has long been a challenge. Aims. We investigate the thermo-hydrodynamic properties of cloud formation processes in the atmospheres of hot Jupiter exoplanets. Methods. We simulate the dynamic atmosphere of HD 189733b with a 3D model that couples 3D radiative-hydrodynamics with a kinetic, microphysical mineral cloud formation module designed for RHD/GCM exoplanet atmosphere simulations. Our simulation includes the feedback effects of cloud advection and settling, gas phase element advection and depletion/replenishment and the radiative effects of cloud opacity. We model the cloud particles as a mix of mineral materials which change in size and composition as they travel through atmospheric thermo-chemical environments. All local cloud properties such as number density, grain size and material composition are time-dependently calculated. Gas phase element depletion as a result of cloud formation is included in the model. In situ effective medium theory and Mie theory is applied to calculate the wavelength dependent opacity of the cloud component. Results. We present a 3D cloud structure of a chemically complex, gaseous atmosphere of the hot Jupiter HD 189733b. Mean cloud particle sizes are typically sub-micron (0.01-0.5 mu m) at pressures less than 1 bar with hotter equatorial regions containing the smallest grains. Denser cloud structures occur near terminator regions and deeper (~1 bar) atmospheric layers. Silicate materials such as MgSiO sub(3)s are found to be abundant at mid-high latitudes, while TiO sub(2)s and SiO sub(2)s dominate the equatorial regions. Elements involved in the cloud formation can be depleted by several orders of magnitude. Conclusions. The interplay between radiative-hydrodynamics and cloud kinetics leads to an inhomogeneous, wavelength dependent opacity cloud structure with properties differing in longitude, latitude and depth. This suggests that transit spectroscopy would sample a variety of cloud particles properties (sizes, composition, densities).
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
).
Context.
The analysis of spectral energy distributions (SEDs) of protoplanetary disks to determine their physical properties is known to be highly degenerate. Hence, a full Bayesian analysis is ...required to obtain parameter uncertainties and degeneracies. The main challenge here is computational speed, as one proper full radiative transfer model requires at least a couple of CPU minutes to compute.
Aims.
We performed a full Bayesian analysis for 30 well-known protoplanetary disks to determine their physical disk properties, including uncertainties and degeneracies. To circumvent the computational cost problem, we created neural networks (NNs) to emulate the SED generation process.
Methods.
We created two sets of MCFOST Monte Carlo radiative transfer disk models to train and test two NNs that predict SEDs for continuous and discontinuous disks, with 18 and 26 free model parameters, respectively. A Bayesian analysis was then performed on 30 protoplanetary disks with SED data collected by the FP7-Space DIANA project to determine the posterior distributions of all parameters. We ran this analysis twice, (i) with old distances and additional parameter constraints as used in a previous study, to compare results, and (ii) with updated distances and free choice of parameters to obtain homogeneous and unbiased model parameters. We evaluated the uncertainties in the determination of physical disk parameters from SED analysis, and detected and quantified the strongest degeneracies.
Results.
The NNs are able to predict SEDs within ~1 ms with uncertainties of about 5% compared to the true SEDs obtained by the radiative transfer code. We find parameter values and uncertainties that are significantly different from previous values obtained by
χ
2
fitting. Comparing the global evidence for continuous and discontinuous disks, we find that 26 out of 30 objects are better described by disks that have two distinct radial zones. The analysed sample shows a significant trend for massive disks to have small scale heights, which is consistent with lower midplane temperatures in massive disks. We find that the frequently used analytic relationship between disk dust mass and millimetre-flux systematically underestimates the dust mass for high-mass disks (dust mass ≥10
−4
M
⊙
). We determine how well the dust mass can be determined with our method for different numbers of flux measurements. As a byproduct, we created an interactive graphical tool that instantly returns the SED predicted by our NNs for any parameter combination.
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.
CAI formation in the early Solar System Woitke, P.; Drążkowska, J.; Lammer, H. ...
Astronomy and astrophysics (Berlin),
07/2024, Volume:
687
Journal Article
Peer reviewed
Calcium-aluminium-rich inclusions (CAIs) are the oldest dated solid materials in the Solar System, and are found as light-coloured crystalline ingredients in carbonaceous chondrite meteorites. Their ...formation time is commonly associated with age zero of the Solar System. Nevertheless, the physical and chemical processes that once led to the formation of these submillimetre- to centimetre-sized mineral particles in the early solar nebula are still a matter of debate. In this paper, we propose a pathway to form such inclusions during the earliest phases of disc evolution. We combine 1D viscous disc evolutionary models with 2D radiative transfer, equilibrium condensation, and new dust opacity calculations. We show that the viscous heating associated with the high accretion rates in the earliest evolutionary phases causes the midplane inside of about 0.5 au to heat up to limiting temperatures of about 1500–1700 K, but no further. These high temperatures force all refractory material components of the inherited interstellar dust grains to sublimate – except for a few Al-Ca-Ti oxides, such as Al 2 O 3 , Ca 2 Al 2 SiO 7 , and CaTiO 3 . This is a recurring and very stable result in all our simulations, because these minerals form a natural thermostat. Once the Mg-Fe silicates are gone, the dust becomes more transparent and the heat is more efficiently transported to the disc surface, which prevents further warming. This thermostat mechanism keeps these minerals above their annealing temperature for hundreds of thousands of years, allowing them to form large pure crystalline particles. These particles are dragged out by the viscously spreading disc, and once they reach a distance of about 0.5 au, the silicates recondense on the surface of the Ca-Al-rich particles, adding an amorphous silicate matrix. We estimate that this mechanism of CAI production works during the first 50 000 yr of disc evolution. These particles then continue to move outward and populate the entire disc up to radii of about 50 au, before the accretion rate eventually subsides, the disc cools, and the particles start to drift inwards.