Spectroscopic diagnostic for the mineralogy of large dust grains

We examine the thermal infrared spectra of large dust grains of different chemical composition and mineralogy. Strong resonances in the optical properties result in detectable spectral structure even when the grain is much larger than the wavelength at which it radiates. We apply this to the thermal infrared spectra of compact amorphous and crystalline silicates. The weak resonances of amorphous silicates at 9.7 and 18 micron virtually disappear for grains larger than about 10 micron. In contrast, the strong resonances of crystalline silicates produce emission dips in the infrared spectra of large grains; these emission dips are shifted in wavelength compared to the emission peaks commonly seen in small crystalline silicate grains. We discuss the effect of a fluffy or compact grain structure on the infrared emission spectra of large grains, and apply our theory to the dust shell surrounding Vega.Comment: Submitted to A&A Letter

FU Orionis - The MIDI/VLTI Perspective

We present the first mid-infrared interferometric measurements of FU Orionis. We clearly resolve structures that are best explained with an optically thick accretion disk. A simple accretion disk model fits the observed SED and visibilities reasonably well and does not require the presence of any additional structure such as a dusty envelope. The inclination and also the position angle of the disk can be constrained from the multibaseline interferometric observations. Our disk model is in general agreement with most published near-infrared interferometric measurements. From the shape and strength of the 8-13 micrometer spectrum the dust composition of the accretion disk is derived for the first time. We conclude that most dust particles are amorphous and already much larger than those typically observed in the ISM. Although the high accretion rate of the system provides both, high temperatures out to large radii and an effective transport mechanism to distribute crystalline grains, we do not see any evidence for crystalline silicates neither in the total spectrum nor in the correlated flux spectra from the inner disk regions. Possible reasons for this non-detection are mentioned. All results are discussed in context with other high-spatial resolution observations of FU Ori and other FU Ori objects. We also address the question whether FU Ori is in a younger evolutionary stage than a classical TTauri star.Comment: 41 pages (aastex style), 11 figures, 8 tables, accepted by Ap

Ab initio phonon dispersion curves and interatomic force constants of barium titanate

The phonon dispersion curves of cubic BaTiO_3 have been computed within a first-principles approach and the results compared to the experimental data. The curves obtained are very similar to those reported for KNbO_3 by Yu and Krakauer [Phys. Rev. Lett. 74, 4067 (1995)]. They reveal that correlated atomic displacements along chains are at the origin of the ferroelectric instability. A simplified model illustrates that spontaneous collective displacements will occur when a dozen of aligned atoms are coupled. The longitudinal interatomic force constant between nearest neighbour Ti and O atoms is relatively weak in comparison to that between Ti atoms in adjacent cells. The small coupling between Ti and O displacements seems however necessary to reproduce a ferroelectric instability.Comment: 12 pages, 4 figure

Dust crystallinity in protoplanetary disks: the effect of diffusion/viscosity ratio

The process of turbulent radial mixing in protoplanetary disks has strong relevance to the analysis of the spatial distribution of crystalline dust species in disks around young stars and to studies of the composition of meteorites and comets in our own solar system. A debate has gone on in the recent literature on the ratio of the effective viscosity coefficient $\nu$ (responsible for accretion) to the turbulent diffusion coefficient $D$ (responsible for mixing). Numerical magneto-hydrodynamic simulations have yielded values between $\nu/D\simeq 10$ (Carballido, Stone & Pringle, 2005) and $\nu/D\simeq 0.85$ (Johansen & Klahr, 2005}). Here we present two analytic arguments for the ratio $\nu/D=1/3$ which are based on elegant, though strongly simplified assumptions. We argue that whichever of these numbers comes closest to reality may be determined {\em observationally} by using spatially resolved mid-infrared measurements of protoplanetary disks around Herbig stars. If meridional flows are present in the disk, then we expect less abundance of crystalline dust in the surface layers, a prediction which can likewise be observationally tested with mid-infrared interferometers.Comment: 9 pages, 5 figures, accepted for publication in A&

C2D Spitzer-IRS spectra of disks around T Tauri stars V. Spectral decomposition

(Abridged) Dust particles evolve in size and lattice structure in protoplanetary disks, due to coagulation, fragmentation and crystallization, and are radially and vertically mixed in disks. This paper aims at determining the mineralogical composition and size distribution of the dust grains in disks around 58 T Tauri stars observed with Spitzer/IRS. We present a spectral decomposition model that reproduces the IRS spectra over the full spectral range. The model assumes two dust populations: a warm component responsible for the 10\mu m emission arising from the disk inner regions and a colder component responsible for the 20-30\mu m emission, arising from more distant regions. We show evidence for a significant size distribution flattening compared to the typical MRN distribution, providing an explanation for the usual boxy 10\mu m feature profile generally observed. We reexamine the crystallinity paradox, observationally identified by Olofsson et al. (2009), and we find a simultaneous enrichment of the crystallinity in both the warm and cold regions, while grain sizes in both components are uncorrelated. Our modeling results do not show evidence for any correlations between the crystallinity and either the star spectral type, or the X-ray luminosity (for a subset of the sample). The size distribution flattening may suggests that grain coagulation is a slightly more effective process than fragmentation in disk atmospheres, and that this imbalance may last over most of the T Tauri phase. This result may also point toward small grain depletion via strong stellar winds or radiation pressure in the upper layers of disk. The non negligible cold crystallinity fractions suggests efficient radial mixing processes in order to distribute crystalline grains at large distances from the central object, along with possible nebular shocks in outer regions of disks that can thermally anneal amorphous grains

Determining the forsterite abundance of the dust around Asymptotic Giant Branch stars

Aims. We present a diagnostic tool to determine the abundance of the crystalline silicate forsterite in AGB stars surrounded by a thick shell of silicate dust. Using six infrared spectra of high mass-loss oxygen rich AGB stars we obtain the forsterite abundance of their dust shells. Methods. We use a monte carlo radiative transfer code to calculate infrared spectra of dust enshrouded AGB stars. We vary the dust composition, mass-loss rate and outer radius. We focus on the strength of the 11.3 and the 33.6 \mu m forsterite bands, that probe the most recent (11.3 \mu m) and older (33.6 \mu m) mass-loss history of the star. Simple diagnostic diagrams are derived, allowing direct comparison to observed band strengths. Results. Our analysis shows that the 11.3 \mu m forsterite band is a robust indicator for the forsterite abundance of the current mass-loss period for AGB stars with an optically thick dust shell. The 33.6 \mu m band of forsterite is sensitive to changes in the density and the geometry of the emitting dust shell, and so a less robust indicator. Applying our method to six high mass-loss rate AGB stars shows that AGB stars can have forsterite abundances of 12% by mass and higher, which is more than the previously found maximum abundance of 5%.Comment: Accepted for publication in A&

The shape and composition of interstellar silicate grains

We investigate the composition and shape distribution of silicate dust grains in the interstellar medium. The effect of the amount of magnesium in the silicate lattice is studied. We fit the spectral shape of the interstellar 10 mu extinction feature as observed towards the galactic center. We use very irregularly shaped coated and non-coated porous Gaussian Random Field particles as well as a statistical approach to model shape effects. For the dust materials we use amorphous and crystalline silicates with various composition and SiC. The results of our analysis of the 10 mu feature are used to compute the shape of the 20 mu silicate feature and to compare this with observations. By using realistic particle shapes we are, for the first time, able to derive the magnesium fraction in interstellar silicates. We find that the interstellar silicates are highly magnesium rich (Mg/(Fe+Mg)>0.9) and that the stoichiometry lies between pyroxene and olivine type silicates. This composition is not consistent with that of the glassy material found in GEMS in interplanetary dust particles indicating that these are, in general, not unprocessed remnants from the interstellar medium. Also, we find a significant fraction of SiC (~3%). We discuss the implications of our results for the formation and evolutionary history of cometary and circumstellar dust. We argue that the fact that crystalline silicates in cometary and circumstellar grains are almost purely magnesium silicates is a natural consequence of our findings that the amorphous silicates from which they were formed were already magnesium rich.Comment: Accepted for publication in A&

EX Lupi in quiescence

EX Lup is the prototype of EXors, a subclass of low-mass pre-main sequence stars whose episodic eruptions are attributed to temporarily increased accretion. In quiescence the optical and near-infrared properties of EX Lup cannot be distinguished from those of normal T Tau stars. Here we investigate whether it is the circumstellar disk structure which makes EX Lup an atypical Class II object. During outburst the disk might undergo structural changes. Our characterization of the quiescent disk is intended to serve as a reference to study the physical changes related to one of EX Lupi's strongest known eruptions in 2008 Jan-Sep. We searched the literature for photometric and spectroscopic observations including ground-based, IRAS, ISO and Spitzer data. After constructing the optical-infrared spectral energy distribution (SED), we compared it with the typical SEDs of other young stellar objects and modeled it using the Monte Carlo radiative transfer code RADMC. A mineralogical decomposition of the 10 micron silicate emission feature and also the description of the optical and near-infrared spectra were performed. The SED is in general similar to that of a typical T Tauri star, though above 7 micron EX Lup emits higher flux. The quiescent phase data suggest low level variability in the optical-mid-infrared domain. Integrating the optical and infrared fluxes we derived a bolometric luminosity of 0.7 L_Sun. The 10 micron silicate profile could be fitted by a mixture consisting of amorphous silicates, no crystalline silicates were found. A modestly flaring disk model with a total mass of 0.025 M_Sun and an outer radius of 150 AU was able to reproduce the observed SED. The derived inner radius of 0.2 AU is larger than the sublimation radius, and this inner gap sets EX Lup aside from typical T Tauri stars.Comment: Accepted for publication in Astronomy and Astrophysics, 10 pages, 6 figure

Dust amorphization in protoplanetary disks

High-energy irradiation of the circumstellar material might impact the structure and the composition of a protoplanetary disk and hence the process of planet formation. In this paper, we present a study on the possible influence of the stellar irradiation, indicated by X-ray emission, on the crystalline structure of the circumstellar dust. The dust crystallinity is measured for 42 class II T Tauri stars in the Taurus star-forming region using a decomposition fit of the 10 micron silicate feature, measured with the Spitzer IRS instrument. Since the sample includes objects with disks of various evolutionary stages, we further confine the target selection, using the age of the objects as a selection parameter. We correlate the X-ray luminosity and the X-ray hardness of the central object with the crystalline mass fraction of the circumstellar dust and find a significant anti-correlation for 20 objects within an age range of approx. 1 to 4.5 Myr. We postulate that X-rays represent the stellar activity and consequently the energetic ions of the stellar winds which interact with the circumstellar disk. We show that the fluxes around 1 AU and ion energies of the present solar wind are sufficient to amorphize the upper layer of dust grains very efficiently, leading to an observable reduction of the crystalline mass fraction of the circumstellar, sub-micron sized dust. This effect could also erase other relations between crystallinity and disk/star parameters such as age or spectral type.Comment: accepted for publication by A&

Resolving the compact dusty discs around binary post-AGB stars using N-band interferometry

We present the first mid-IR long baseline interferometric observations of the circumstellar matter around binary post-AGB stars. Two objects, SX Cen and HD 52961, were observed using the VLTI/MIDI instrument during Science Demonstration Time. Both objects are known binaries for which a stable circumbinary disc is proposed to explain the SED characteristics. This is corroborated by our N-band spectrum showing a crystallinity fraction of more than 50 % for both objects, pointing to a stable environment where dust processing can occur. Surprisingly, the dust surrounding SX Cen is not resolved in the interferometric observations providing an upper limit of 11 mas (or 18 AU at the distance of this object) on the diameter of the dust emission. This confirms the very compact nature of its circumstellar environment. The dust emission around HD 52961 originates from a very small but resolved region, estimated to be ~ 35 mas at 8 micron and ~ 55 mas at 13 micron. These results confirm the disc interpretation of the SED of both stars. In HD 52961, the dust is not homogeneous in its chemical composition: the crystallinity is clearly concentrated in the hotter inner region. Whether this is a result of the formation process of the disc, or due to annealing during the long storage time in the disc is not clear.Comment: 12 pages, 10 figures, accepted for publication in A &