Dust Settling Research Articles

Young stars inϵChamaleontis and their disks: disk evolution in sparse associations

(abridge) The nearby young stellar association Epsilon Cha association has an estimated age of 3-5 Myr, making it an ideal laboratory to study the disk dissipation process and provide empirical constraints on the timescale of planet formation. We combine the available literature data with our Spitzer IRS spectroscopy and VLT/VISIR imaging data. The very low mass stars USNO-B120144.7 and 2MASS J12005517 show globally depleted spectral energy distributions pointing at strong dust settling. 2MASS J12014343 may have a disk with a very specific inclination where the central star is effectively screened by the cold outer parts of a flared disk but the 10 micron radiation of the warm inner disk can still reach us. We find the disks in sparse stellar associations are dissipated more slowly than those in denser (cluster) environments. We detect C_{2}H_{2} rovibrational band around 13.7 micron on the IRS spectrum of USNO-B120144.7. We find strong signatures of grain growth and crystallization in all Epsilon Cha members with 10 micron features detected in their IRS spectra. We combine the dust properties derived in the Epsilon Cha sample with those found using identical or similar methods in the MBM 12, Coronet cluster, Eta Cha associations, and in the cores to disks (c2d) legacy program. We find that disks around low-mass young stars show a negative radial gradient in the mass-averaged grain size and mass fraction of crystalline silicates. A positive correlation exists between the mass-averaged grain sizes of amorphous silicates and the accretion rates if the latter is above ~10^{-9} Msun/yr, possibly indicating that those disks are sufficiently turbulent to prevent grains of several microns in size to sink into the disk interior.

Modeling the onset of photosynthesis after the Chicxulub asteroid impact

We do a preliminary modelling of the photosynthetic rates of phytoplankton at the very beginning of the Paleogene, just after the impact of the Chicxulub asteroid, which decisively contributed to the last known mass extinction of the Phanerozoic eon. We assume the worst possible scenario from the photobiological point of view: an already clear atmosphere with no ozone, as the timescale for soot and dust settling (years) is smaller than that of the full ozone regeneration (decades). Even in these conditions we show that most phytoplankton species would have had reasonable potential for photosynthesis in all the three main optical ocean water types. This modelling could help explain why the recovery of phytoplankton was relatively rapid after the huge environmental stress of that asteroid impact. In a more general scope, it also reminds us of the great resilience of the unicellular biosphere against huge environmental perturbations.

DoAr 33: a good candidate for revealing dust growth and settling in protoplanetary disks

Aims. We aim to evaluate the evolutionary stage of the circumstellar disk around DoAr 33, a T Tauri star in the Ophiuchus molecular cloud and a promising target for follow-up observations to find signs of dust evolution in protoplanetary disks.

Dust grain growth and settling in initial gaseous giant protoplanets

Dust grain growth and settling time inside initial gaseous giant protoplanets in the mass range 0.3 to 5 Jovian masses, formed by gravitational instability, have been investigated. We have determined the distribution of thermodynamic and physical variables inside the protoplanets solving the structure equations assuming their gas blobs to be fully convective and with this distribution we have calculated growth and settling time of grains with different initial sizes (10−2 cm ≤ r0 ≤ 1 cm). The results of our calculations are found to be in good agreement with those obtained by different approaches.

Impacts on iron solubility in the mineral dust by processes in the source region and the atmosphere: A review

Mineral dust is a complex entity containing a range of iron minerals including poorly crystalline to crystalline iron oxides to clay minerals. Important progress has been made to characterize iron mineralogical compositions in the dust recently. These include the quantification of the content of crystalline hematite and goethite, which appear to show a regional variation in North Africa as a result of the differences in the degree of chemical weathering. Fractional Fe solubility (dissolved to total iron, FFS) in the atmospheric aerosols has been reported to range from 0.1% to 80%. However, FFS is usually less than 0.5% in the non-atmospherically-processed dust, suggesting that FFS can be enhanced by atmospheric processes. One of the atmospheric processes, gravitational settling of dust, which has been previously hypothesized to cause the abovementioned enhancement of FFS during dust transport has been shown to be insignificant. Cycling of dust particles in the clouds, in which pH is usually higher than 4, and in the aerosol phase, in which pH is usually substantially lower, can significantly affect iron speciation and FFS. Laboratory experiments showed that a significant amount of iron (>0.5%) can only be solubilized in the dust when pH is lower than 4. These laboratory data suggest that acid processing rather than cloud processing might be a prime mechanism to cause an increase in FFS in the dust during transport. Further laboratory studies, field measurements, and modelling are needed to increase the ability of models to quantify the atmospheric processing of iron in the dust.

Transition disks: four candidates for ongoing giant planet formation in Ophiuchus

A large set of Spitzer-selected transitional disks in the Ophiuchus molecular cloud was recently examined by Cieza et al (2010), and 4 of the targets were identified as (giant) planet-forming candidates based on the morphology of their Spectral Energy Distributions (SEDs), the apparent lack of stellar companions, and evidence for accretion. Here we characterize the structures of these disks modeling their optical, infrared and (sub)millimeter SEDs. We use the Monte Carlo radiative transfer package RADMC to construct a parametric model of the dust distribution in a flared disk with an inner cavity and calculate the temperature structure consistent with the density profile, in thermal equilibrium with the irradiating star. For each object, we conducted a Bayesian exploration of the parameter space generating Monte Carlo Markov Chains (MCMC) that allow the identification of the best fit parameters and to constrain their range of statistical confidence. Our calculations point to the presence of evacuated cavities with radii ~2-8 AU, consistent with having been carved by embedded giant planets. We found parameter values consistent with those previously given in the literature, indicating a mild degree of grain growth and dust settling, which deserves to be investigated with further modeling and follow up observations. Resolved images with (sub)millimeter interferometers are required to break some of the degeneracies of the models and better constrain the physical properties of these fascinating disks.

ON THE TRANSITIONAL DISK CLASS: LINKING OBSERVATIONS OF T TAURI STARS AND PHYSICAL DISK MODELS

Two decades ago "transitional disks" described spectral energy distributions (SEDs) of T Tauri stars with small near-IR excesses, but significant mid- and far-IR excesses. Many inferred this indicated dust-free holes in disks, possibly cleared by planets. Recently, this term has been applied disparately to objects whose Spitzer SEDs diverge from the expectations for a typical full disk. Here we use irradiated accretion disk models to fit the SEDs of 15 such disks in NGC 2068 and IC 348. One group has a "dip" in infrared emission while the others' continuum emission decreases steadily at all wavelengths. We find that the former have an inner disk hole or gap at intermediate radii in the disk and we call these objects "transitional" and pre-transitional" disks, respectively. For the latter group, we can fit these SEDs with full disk models and find that millimeter data are necessary to break the degeneracy between dust settling and disk mass. We suggest the term "transitional" only be applied to objects that display evidence for a radical change in the disk's radial structure. Using this definition, we find that transitional and pre-transitional disks tend to have lower mass accretion rates than full disks and that transitional disks have lower accretion rates than pre-transitional disks. These reduced accretion rates onto the star could be linked to forming planets. Future observations of transitional and pre-transitional disks will allow us to better quantify the signatures of planet formation in young disks.

Probing the turbulent mixing strength in protoplanetary disks across the stellar mass range: no significant variations

Dust settling and grain growth are the first steps in the planet-formation process in protoplanetary disks. These disks are observed around stars with different spectral types, and there are indications that the disks around lower mass stars are significantly flatter, which could indicate that they settle and evolve faster, or in a different way. We aim to test this assumption by modeling the median spectral energy distributions (SEDs) of three samples of protoplanetary disks: around Herbig stars, T Tauri stars and brown dwarfs. We focus on the turbulent mixing strength to avoid a strong observational bias from disk and stellar properties that depend on stellar mass. We generated SEDs with the radiative transfer code MCMax, using a hydrostatic disk structure and settling the dust in a self-consistent way with the alpha-prescription to probe the turbulent mixing strength. We are able to fit all three samples with a disk with the same input parameters, scaling the inner edge to the dust evaporation radius and disk mass to millimeter photometry. The Herbig stars require a special treatment for the inner rim regions, while the T-Tauri stars require viscous heating, and the brown dwarfs lack a good estimate of the disk mass because only few millimeter detections exist. We find that the turbulent mixing strength does not vary across the stellar mass range for a fixed grain size distribution and gas-to-dust ratio. Regions with the same temperature have a self-similar vertical structure independent of stellar mass, but regions at the same distance from the central star appear more settled in disks around lower mass stars. We find a relatively low turbulent mixing strength of alpha = 10^(-4) for a standard grain size distribution, but our results are also consistent with alpha = 0.01 for a grain size distribution with fewer small grains or a lower gas-to-dust ratio.

Grain charging in protoplanetary discs

Recent work identified a growth barrier for dust coagulation that originates in the electric repulsion between colliding particles. Depending on its charge state, dust material may have the potential to control key processes towards planet formation such as MHD (magnetohydrodynamic) turbulence and grain growth which are coupled in a two-way process. We quantify the grain charging at different stages of disc evolution and differentiate between two very extreme cases: compact spherical grains and aggregates with fractal dimension D_f = 2. Applying a simple chemical network that accounts for collisional charging of grains, we provide a semi-analytical solution. This allowed us to calculate the equilibrium population of grain charges and the ionisation fraction efficiently. The grain charging was evaluated for different dynamical environments ranging from static to non-stationary disc configurations. The results show that the adsorption/desorption of neutral gas-phase heavy metals, such as magnesium, effects the charging state of grains. The greater the difference between the thermal velocities of the metal and the dominant molecular ion, the greater the change in the mean grain charge. Agglomerates have more negative excess charge on average than compact spherical particles of the same mass. The rise in the mean grain charge is proportional to N**(1/6) in the ion-dust limit. We find that grain charging in a non-stationary disc environment is expected to lead to similar results. The results indicate that the dust growth and settling in regions where the dust growth is limited by the so-called "electro-static barrier" do not prevent the dust material from remaining the dominant charge carrier.

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.

A HERSCHEL SEARCH FOR COLD DUST IN BROWN DWARF DISKS: FIRST RESULTS

We report initial results from a {\it Herschel} program to search for far-infrared emission from cold dust around a statistically significant sample of young brown dwarfs. The first three objects in our survey are all detected at 70\micron, and we report the first detection of a brown dwarf at 160\micron. The flux densities are consistent with the presence of substantial amounts of cold dust in the outer disks around these objects. We modeled the SED's with two different radiative transfer codes. We find that a broad range of model parameters provides a reasonable fit to the SED's, but that the addition of our 70\micron, and especially the 160\micron\ detection enables strong lower limits to be placed on the disk masses since most of the mass is in the outer disk. We find likely disk masses in the range of a few $\times 10^{-6}$ to $10^{-4}$ \msun. Our models provide a good fit to the SED's and do not require dust settling.

A single-dish survey of the HCO+, HCN, and CN emission toward the T Tauri disk population in Taurus

(Abridged) As the stellar X-ray and UV light penetration of a protoplanetary disk depends sensitively on the dust properties, trace molecular species like HCO+, HCN, and CN are expected to show marked differences from photoprocessing effects as the dust content in the disk evolves. We investigate the evolution of the UV irradiation of the molecular gas in a sample of classical T Tauri stars in Taurus that exhibit a wide range in grain growth and dust settling properties. We obtained HCO+ (J=3-2), HCN (J=3-2), and CN (J=2-1) observations of 13 sources with the JCMT. Our sample has 1.3mm fluxes in excess of 75mJy, indicating the presence of significant dust reservoirs; a range of dust settling as traced through their spectral slopes between 6, 13, and 25 microns; and varying degrees of grain growth as extrapolated from the strength of their 10-micron silicate emission features. We compare the emission line strengths with the sources' continuum flux and infrared features, and use detailed modeling based on two different model prescriptions to compare typical disk abundances for HCO+, HCN, and CN with the gas-line observations for our sample. We detected HCO+ (3-2) toward 6 disks, HCN (3-2) from 0 disks, and CN (2-1) toward 4 disks. For the complete sample, there is no correlation between the gas-line strengths or their ratios and either the sources' dust continuum flux or infrared slopes.

DUST PROPERTIES AND DISK STRUCTURE OF EVOLVED PROTOPLANETARY DISKS IN Cep OB2: GRAIN GROWTH, SETTLING, GAS AND DUST MASS, AND INSIDE-OUT EVOLUTION

We present Spitzer/IRS spectra of 31 TTS and IRAM/1.3mm observations for 34 low- and intermediate-mass stars in the Cep OB2 region. Including our previously published data, we analyze 56 TTS and the 3 intermediate-mass stars with silicate features in Tr 37 (~4 Myr) and NGC 7160 (~12 Myr). The silicate emission features are well reproduced with a mixture of amorphous (with olivine, forsterite, and silica stoichiometry) and crystalline grains (forsterite, enstatite). We explore grain size and disk structure using radiative transfer disk models, finding that most objects have suffered substantial evolution (grain growth, settling). About half of the disks show inside-out evolution, with either dust-cleared inner holes or a radially-dependent dust distribution, typically with larger grains and more settling in the innermost disk. The typical strong silicate features require nevertheless the presence of small dust grains, and could be explained by differential settling according to grain size, anomalous dust distributions, and/or optically thin dust populations within disk gaps. M-type stars tend to have weaker silicate emission and steeper SEDs than K-type objects. The inferred low dust masses are in a strong contrast with the relatively high gas accretion rates, suggesting global grain growth and/or an anomalous gas to dust ratio. Transition disks (TD) in the Cep OB2 region display strongly processed grains, suggesting that they are dominated by dust evolution and settling. Finally, the presence of rare but remarkable disks with strong accretion at old ages reveals that some very massive disks may still survive to grain growth, gravitational instabilities, and planet formation.

PHOTOELECTRIC CROSS-SECTIONS OF GAS AND DUST IN PROTOPLANETARY DISKS

We provide simple polynomial fits to the X-ray photoelectric cross-sections (0.03 < E < 10keV) for mixtures of gas and dust found in protoplanetary disks. Using the solar elemental abundances of Asplund et al. (2009) we treat the gas and dust components separately, facilitating the further exploration evolutionary processes such as grain settling and gain growth. We find that blanketing due to advanced grain-growth (a_max > 1 micron) can reduce the X-ray opacity of dust appreciably at E_X ~ 1keV, coincident with the peak of typical T Tauri X-ray spectra. However, the reduction of dust opacity by dust settling, which is known to occur in protoplanetary disks, is probably a more significant effect. The absorption of 1-10keV X-rays is dominated by gas opacity once the dust abundance has been reduced to about 1% of its diffuse interstellar value. The gas disk establishes a floor to the opacity at which point X-ray transport becomes insensitive to further dust evolution. Our choice of fitting function follows that of Morrison & McCammon (1983), providing a degree of backward-compatibility.

THE PROPAGATION OF Lyα IN EVOLVING PROTOPLANETARY DISKS

We study the role resonant scattering plays in the transport of Lyman-alpha photons in accreting protoplanetary disk systems subject to varying degrees of dust settling. While the intrinsic stellar FUV spectrum of accreting T Tauri systems may already be dominated by a strong, broad Lyman-alpha line (~80% of the FUV luminosity), we find that resonant scattering further enhances the Lyman-alpha density in the deep molecular layers of the disk. Lyman-alpha is scattered downwards efficiently by the photodissociated atomic hydrogen layer that exists above the molecular disk. In contrast, FUV-continuum photons pass unimpeded through the photodissociation layer, and (forward-)scatter inefficiently off dust grains. Using detailed, adaptive grid Monte Carlo radiative transfer simulations we show that the resulting Lyman-alpha/FUV-continuum photon density ratio is strongly stratified; FUV-continuum dominated in the photodissociation layer, and Lyman-alpha dominated field in the molecular disk. The enhancement is greatest in the interior of the disk (r~1AU) but is also observed in the outer disk (r~100AU). The majority of the total disk mass is shown to be increasingly Lyman-alpha dominated as dust settles towards the midplane.

Statistical Analysis of Effect of a Negative Ionization System on Dust Settling in a Confinement Swine Building

Statistical Analysis of Effect of a Negative Ionization System on Dust Settling in a Confinement Swine Building

THE SPITZER INFRARED SPECTROGRAPH SURVEY OF T TAURI STARS IN TAURUS

We present 161 Spitzer Infrared Spectrograph (IRS) spectra of T Tauri stars and young brown dwarfs in the Taurus star-forming region. All of the targets were selected based on their infrared excess and are therefore surrounded by protoplanetary disks; they form the complete sample of all available IRS spectra of T Tauri stars with infrared excesses in Taurus. We also present the IRS spectra of seven Class 0/I objects in Taurus to complete the sample of available IRS spectra of protostars in Taurus. We use spectral indices that are not significantly affected by extinction to distinguish between envelope- and disk-dominated objects. Together with data from the literature, we construct spectral energy distributions for all objects in our sample. With spectral indices derived from the IRS spectra we infer disk properties such as dust settling and the presence of inner disk holes and gaps. We find a transitional disk frequency, which is based on objects with unusually large 13–31 μm spectral indices indicative of a wall surrounding an inner disk hole, of about 3%, and a frequency of about 20% for objects with unusually large 10 μm features, which could indicate disk gaps. The shape and strength of the 10 μm silicate emission feature suggests weaker 10 μm emission and more processed dust for very low mass objects and brown dwarfs (spectral types M6–M9). These objects also display weaker infrared excess emission from their disks, but do not appear to have more settled disks than their higher-mass counterparts. We find no difference for the spectral indices and properties of the dust between single and multiple systems.

TRANSITIONAL AND PRE-TRANSITIONAL DISKS: GAP OPENING BY MULTIPLE PLANETS?

We use 2D hydrodynamic simulations of viscous disks to examine whether dynamically-interacting multiple giant planets can explain the large gaps (spanning over one order of magnitude in radius) inferred for the transitional and pre-transitional disks around T Tauri stars. In the absence of inner disk dust depletion, we find that it requires three to four giant planets to open up large enough gaps to be consistent with inferences from spectral energy distributions, because the gap width is limited by the tendency of the planets to be driven together into 2:1 resonances. With very strong tidal torques and/or rapid planetary accretion, fewer planets can also generate a large cavity interior to the locally formed gap(s) by preventing outer disk material from moving in. In these cases, however, the reduction of surface density produces a corresponding reduction in the inner disk accretion rate onto the star; this makes it difficult to explain the observed accretion rates of the pre/transitional disks. We find that even with four planets in disks, additional substantial dust depletion is required to explain observed disk gaps/holes. Substantial dust settling and growth, with consequent significant reductions in optical depths, is inferred for typical T Tauri disks in any case, and an earlier history of dust growth is consistent with the hypothesis that pre/transitional disks are explained by the presence of giant planets. We conclude that the depths and widths of gaps, and disk accretion rates in pre/transitional disks cannot be reproduced by a planet-induced gap opening scenario alone. Significant dust depletion is also required within the gaps/holes. Order of magnitude estimates suggest the mass of small dust particles (<1 micron) relative to the gas must be depleted to 1e-5 -- 1e-2 of the interstellar medium value, implying a very efficient mechanism of small dust removal or dust growth.

Dust settling and rapid planetary migration

Planetary migration is essential to explain the observed mass-period relation for exoplanets. Without some stopping mechanism, the tidal, resonant interaction between planets and their gaseous disc generally causes the planets to migrate inward so efficiently that they plunge into the host star within the gaseous disc lifetime ($\sim $ 1-3 Myrs). We investigate planetary migration by analytically calculating the migration rate and time within self-consistently computed, radiatively heated discs around M stars in which the effects of dust settling are included. We show that dust settling lowers the disc temperature and raises the gas density in the mid-plane. This inescapable evolution of disc structure speeds up type I planetary migration for lower mass bodies by up to a factor of about 2. We also examine the effects of dust settling on the gap-opening mass and type II migration, and find that the gap-opening mass is reduced by a factor of 2 and type II migration becomes slower by a factor of 2. While dust settling can somewhat alleviate the problem of planetary migration for more massive planets, the more rapid migration of low mass planets and planetary cores requires a robust slowing mechanism.

Observing dust settling and coagulation in circumstellar discs

Context. Circumstellar discs are expected to be the nursery of planets. Grain growth within these discs is the first step in the planet formation process in the core-accretion gas-capture scenario.