With the number of confirmed rocky exoplanets increasing steadily, their characterisation and the search for exoplanetary biospheres is becoming an increasingly urgent issue in astrobiology. We aim to investigate the possibility of characterising an exoplanet (in terms of habitability, geology, presence of life etc.) by studying material ejected from the surface during an impact event. For given parameters characterising the impact event, we estimate the escaping mass and assess its subsequent collisional evolution in a circumstellar orbit, assuming a Sun-like host star. We calculate the fractional luminosity of the dust as a function of time after the impact event and study its detectability with current and future instrumentation. We consider the possibility to constrain the dust composition, giving information on the geology or the presence of a biosphere. As examples, we investigate whether calcite, silica or ejected microorganisms could be detected. For a 20 km diameter impactor, we find that the dust mass escaping the exoplanet is roughly comparable to the zodiacal dust. The collisional evolution is best modelled by considering two independent dust populations, a spalled population consisting of non-melted ejecta evolving on timescales of millions of years, and dust recondensed from melt or vapour evolving on much shorter timescales. While the presence of dust can potentially be inferred with current telescopes, studying its composition requires advanced instrumentation not yet available. The direct detection of biological matter turns out to be extremely challenging. Despite considerable difficulties (small dust masses, noise such as exozodiacal dust etc.), studying dusty material ejected from an exoplanetary surface might become an interesting complement to atmospheric studies in the future.

The debris disk around $\beta$~Pictoris is known to contain gas. Previous ALMA observations revealed a CO belt at $\sim$85 au with a distinct clump, interpreted as a location of enhanced gas production. Photodissociation converts CO into C and O within $\sim$50 years. We resolve CI emission at 492 GHz using ALMA and study its spatial distribution. CI shows the same clump as seen for CO. This is surprising, as C is expected to quickly spread in azimuth. We derive a low C mass (between $5\times10^{-4}$ and $3.1\times10^{-3}$ M$_\oplus$), indicating that gas production started only recently (within $\sim$5000 years). No evidence is seen for an atomic accretion disk inwards of the CO belt, perhaps because the gas did not yet have time to spread radially. The fact that C and CO share the same asymmetry argues against a previously proposed scenario where the clump is due to an outward migrating planet trapping planetesimals in an resonance; nor can the observations be explained by an eccentric planetesimal belt secularly forced by a planet. Instead, we suggest that the dust and gas disks should be eccentric. Such a configuration, we further speculate, might be produced by a recent tidal disruption event. Assuming that the disrupted body has had a CO mass fraction of 10%, its total mass would be $\gtrsim$3 $M_\mathrm{Moon}$.

Meteorites contain minerals from Solar System asteroids with different properties (like size, presence of water, core formation). We provide new mid-IR transmission spectra of powdered meteorites to obtain templates of how mid-IR spectra of asteroidal debris would look like. This is essential for interpreting mid-IR spectra of past and future space observatories, like the James Webb Space Telescope. We show that the transmission spectra of wet and dry chondrites, carbonaceous and ordinary chondrites and achondrite and chondrite meteorites are distinctly different in a way one can distinguish in astronomical mid-IR spectra. The two observables that spectroscopically separate the different meteorites groups (and thus the different types of parent bodies) are the pyroxene-olivine feature strength ratio and the peak shift of the olivine spectral features due to an increase in the iron concentration of the olivine.

We present far-infrared and submillimeter maps from the Herschel Space Observatory and the James Clerk Maxwell Telescope of the debris disk host star AU Microscopii. Disk emission is detected at 70, 160, 250, 350, 450, 500 and 850 micron. The disk is resolved at 70, 160 and 450 micron. In addition to the planetesimal belt, we detect thermal emission from AU Mic's halo for the first time. In contrast to the scattered light images, no asymmetries are evident in the disk. The fractional luminosity of the disk is $3.9 \times 10^{-4}$ and its mm-grain dust mass is 0.01 MEarth (+/- 20%). We create a simple spatial model that reconciles the disk SED as a blackbody of 53 +/- 2 K (a composite of 39 and 50 K components) and the presence of small (non-blackbody) grains which populate the extended halo. The best fit model is consistent with the "birth ring" model explored in earlier works, i.e., an edge-on dust belt extending from 8.8-40 AU, but with an additional halo component with an $r^{-1.5}$ surface density profile extending to the limits of sensitivity (140 AU). We confirm that AU Mic does not exert enough radiation force to blow out grains. For stellar mass loss rates of 10-100x solar, compact (zero porosity) grains can only be removed if they are very small, consistently with previous work, if the porosity is 0.9, then grains approaching 0.1 micron can be removed via corpuscular forces (i.e., the stellar wind).

We aim to constrain the structure of the circumstellar material around the post-AGB binary and RV Tauri pulsator AC Her. We want to constrain the spatial distribution of the amorphous as well as of the crystalline dust. We present very high-quality mid-IR interferometric data that were obtained with MIDI/VLTI. We analyse the MIDI data and the full SED, using the MCMax radiative transfer code, to find a good structure model of AC Her's circumbinary disk. We include a grain size distribution and midplane settling of dust self-consistently. The spatial distribution of crystalline forsterite in the disk is investigated with the mid-IR features, the 69~$\mu$m band and the 11.3~$\mu$m signatures in the interferometric data. All the data are well fitted. The inclination and position angle of the disk are well determined at i=50+-8 and PA=305+-10. We firmly establish that the inner disk radius is about an order of magnitude larger than the dust sublimation radius. Significant grain growth has occurred, with mm-sized grains being settled to the midplane of the disk. A large dust mass is needed to fit the sub-mm fluxes. By assuming {\alpha}=0.01, a good fit is obtained with a small grain size power law index of 3.25, combined with a small gas/dust ratio <10. The resulting gas mass is compatible with recent estimates employing direct gas diagnostics. The spatial distribution of the forsterite is different from the amorphous dust, as more warm forsterite is needed in the surface layers of the inner disk. The disk in AC Her is very evolved, with its small gas/dust ratio and large inner hole. Mid-IR interferometry offers unique constraints, complementary to mid-IR features, for studying the mineralogy in disks. A better uv coverage is needed to constrain in detail the distribution of the crystalline forsterite in AC Her, but we find strong similarities with the protoplanetary disk HD100546.

Context: AGB stars lose a large percentage of their mass in a dust-driven wind. This creates a circumstellar envelope, which can be studied through thermal dust emission and molecular emission lines. In the case of high mass-loss rates, this study is complicated by the high optical depths and the intricate coupling between gas and dust radiative transfer characteristics. An important aspect of the physics of gas-dust interactions is the strong influence of dust on the excitation of several molecules, including H2O. Aims: The dust and gas content of the envelope surrounding the high mass-loss rate OH/IR star OH 127.8+0.0, as traced by Herschel observations, is studied, with a focus on the H2O content and the dust-to-gas ratio. We report detecting a large number of H2O vapor emission lines up to J = 9 in the Herschel data, for which we present the measured line strengths. Methods: The treatments of both gas and dust species are combined using two numerical radiative transfer codes. The method is illustrated for both low and high mass-loss-rate sources. Specifically, we discuss different ways of assessing the dust-to-gas ratio: 1) from the dust thermal emission spectrum and the CO molecular gas line strengths; 2) from the momentum transfer from dust to gas and the measured gas terminal velocity; and 3) from the determination of the required amount of dust to reproduce H2O lines for a given H2O vapor abundance. These three diagnostics probe different zones of the outflow, for the first time allowing an investigation of a possible radial dependence of the dust-to-gas ratio. Results: ... Continued in article.

We aim at determining the spatial distribution of the gas and dust in star-forming regions and address their relative abundances in quantitative terms. We also examine the dust opacity exponent beta for spatial and/or temporal variations. Using mapping observations of the very dense rho Oph A core, we examined standard 1D and non-standard 3D methods to analyse data of far-infrared and submillimeter (submm) continuum radiation. The resulting dust surface density distribution can be compared to that of the gas. The latter was derived from the analysis of accompanying molecular line emission, observed with Herschel from space and with APEX from the ground. As a gas tracer we used N2H+, which is believed to be much less sensitive to freeze-out than CO and its isotopologues. Radiative transfer modelling of the N2H+(J=3-2) and (J=6-5) lines with their hyperfine structure explicitly taken into account provides solutions for the spatial distribution of the column density N(H2), hence the surface density distribution of the gas. The gas-to-dust mass ratio is varying across the map, with very low values in the central regions around the core SM 1. The global average, =88, is not far from the canonical value of 100, however. In rho Oph A, the exponent beta of the power-law description for the dust opacity exhibits a clear dependence on time, with high values of 2 for the envelope-dominated emission in starless Class -1 sources to low values close to 0 for the disk-dominated emission in Class III objects. beta assumes intermediate values for evolutionary classes in between. Since beta is primarily controlled by grain size, grain growth mostly occurs in circumstellar disks. The spatial segregation of gas and dust, seen in projection toward the core centre, probably implies that, like C18O, also N2H+ is frozen onto the grains.

Low- and intermediate-mass stars go through a period of intense mass-loss at the end of their lives in a phase known as the asymptotic giant branch (AGB). During the AGB a significant fraction of their initial mass is expelled in a stellar wind. This process controls the final stages of their evolution and contributes to the chemical evolution of galaxies. However, the wind-driving mechanism of AGB stars is not yet well understood, especially so for oxygen-rich sources. Characterizing both the present-day mass-loss and wind structure and the evolution of the mass-loss rate of such stars is paramount to advancing our understanding of this processes. We modelled the dust envelope of W Hya using an advanced radiative transfer code. The dust model was analysed in the light of a previously calculated gas-phase wind model and compared to measurements available in the literature, such as infrared spectra, infrared images, and optical scattered light fractions. We find that the dust spectrum of W Hya can partly be explained by a gravitationally bound dust shell that probably is responsible for most of the amorphous Al$_2$O$_3$ emission. The composition of the large ($\sim$\,0.3\,$\mu$m) grains needed to explain the scattered light cannot be constrained, but probably is dominated by silicates. Silicate emission in the thermal infrared was found to originate from beyond 40 AU from the star and we find that they need to have substantial near-infrared opacities to be visible at such large distances. The increase in near-infrared opacity of the dust at these distances roughly coincides with a sudden increase in expansion velocity as deduced from the gas-phase CO lines. Finally, the recent mass loss of W Hya is confirmed to be highly variable and we identify a strong peak in the mass-loss rate that occurred about 3500 years ago and lasted for a few hundred years.

We study the grain properties and location of the forsterite crystals in the circumstellar environment of the pre-planetary nebula (PPN) IRAS 17150-3224 in order to learn more about the as yet poorly understood evolutionary phase prior to the PPN. We use the best-fit model for IRAS 17150-3224 of Meixner et al. (2002) and add forsterite to this model. We investigate different spatial distributions and grain sizes of the forsterite crystals in the circumstellar environment. We compare the spectral bands of forsterite in the mid-infrared and at 69 micrometre in radiative transport models to those in ISO-SWS and Herschel/PACS observations. We can reproduce the non-detection of the mid-infrared bands and the detection of the 69 micrometre feature with models where the forsterite is distributed in the whole outflow, in the superwind region, or in the AGB-wind region emitted previous to the superwind, but we cannot discriminate between these three models. To reproduce the observed spectral bands with these three models, the forsterite crystals need to be dominated by a grain size population of 2 micrometre up to 6 micrometre. We hypothesise that the large forsterite crystals were formed after the superwind phase of IRAS 17150-3224, where the star developed an as yet unknown hyperwind with an extremely high mass-loss rate (10^-3 Msol/yr). The high densities of such a hyperwind could be responsible for the efficient grain growth of both amorphous and crystalline dust in the outflow. Several mechanisms are discussed that might explain the lower-limit of 2 micrometre found for the forsterite grains, but none are satisfactory. Among the mechanisms explored is a possible selection effect due to radiation pressure based on photon scattering on micron-sized grains.

Binaries with circumbinary disks are commonly found among optically bright post-AGB stars. Although clearly linked to binary interaction processes, the formation, evolution and fate of these disks are still badly understood. Due to their compactness, interferometric techniques are required to resolve them. Here, we discuss our high-quality multiwavelength interferometric data of two prototypical yet very different post-AGB binaries, AC and 89 Herculis, as well as the modeling thereof with radiative transfer models. A detailed account of the data and models of both objects is published in three separate papers elsewhere; here we focus on comparing the modeling results for the two objects. In particular we discuss the successes and limitations of the models which were developed for protoplanetary disks around young stars. We conclude that multiwavelength high-angular-resolution observations and radiative transfer disk models are indispensible to understand these complex interacting objects and their place in the grand scheme of the (binary) evolution of low and intermediate mass stars.

$\tau$ Ceti is a nearby, mature G-type star very similar to our Sun, with a massive Kuiper Belt analogue (Greaves et al. 2004) and possible multiplanet system (Tuomi et al. 2013) that has been compared to our Solar System. We present Herschel Space Observatory images of the debris disk, finding the disk is resolved at 70 and 160 microns, and marginally resolved at 250 microns. The Herschel images and infrared photometry from the literature are best modelled using a wide dust annulus with an inner edge between 1-10 AU and an outer edge at ~55 AU, inclined from face-on by 35$\pm$10 degrees, and with no significant azimuthal structure. We model the proposed tightly-packed planetary system of five super-Earths and find that the innermost dynamically stable disk orbits are consistent with the inner edge found by the observations. The photometric modelling, however, cannot rule out a disk inner edge as close to the star as 1 AU, though larger distances produce a better fit to the data. Dynamical modelling shows that the 5 planet system is stable with the addition of a Neptune or smaller mass planet on an orbit outside 5 AU, where the Tuomi et al. analysis would not have detected a planet of this mass.

The spectra of protoplanetary disks contain mid- and far- infrared emission features produced by forsterite dust grains. The spectral features contain information about the forsterite temperature, chemical composition and grain size. We aim to characterize how the 23 and 69 micron features can be used to constrain the physical locations of forsterite in disks. We check for consistency between two independent forsterite temperature measurements: the 23/69 feature strength ratio and the shape of the 69 micron band. We performed radiative transfer modeling to study the effect of disk properties to the forsterite spectral features. Temperature-dependent forsterite opacities were considered in self-consistent models to compute forsterite emission from protoplanetary disks. Modelling grids are presented to study the effects of grain size, disk gaps, radial mixing and optical depth to the forsterite features. Independent temperature estimates derived from the 23/69 feature strength ratio and the 69 micron band shape are most inconsistent for HD141569 and Oph IRS 48. A case study of the disk of HD141569 shows two solutions to fit the forsterite spectrum. A model with T ~ 40 K, iron-rich (~0-1 % Fe) and 1 micron forsterite grains, and a model with warmer (T ~ 100 K), iron-free, and larger (10 micron) grains. We find that for disks with low upper limits of the 69 micron feature (most notably in flat, self-shadowed disks), the forsterite must be hot, and thus close to the star. We find no correlation between disk gaps and the presence or absence of forsterite features. We argue that the 69 micron feature of the evolved transitional disks HD141569 and Oph IRS 48 is most likely a tracer of larger (i.e. ~10 micron) forsterite grains.

We present 48 Herschel/PACS spectra of evolved stars in the wavelength range of 67-72 $\mu$m. This wavelength range covers the 69 $\mu$m band of crystalline olivine ($\text{Mg}_{2-2x}\text{Fe}_{(2x)}\text{SiO}_{4}$). The width and wavelength position of this band are sensitive to the temperature and composition of the crystalline olivine. Our sample covers a wide range of objects: from high mass-loss rate AGB stars (OH/IR stars, $\dot M \ge 10^{-5}$ M$_\odot$/yr), through post-AGB stars with and without circumbinary disks, to planetary nebulae and even a few massive evolved stars. The goal of this study is to exploit the spectral properties of the 69 $\mu$m band to determine the composition and temperature of the crystalline olivine. Since the objects cover a range of evolutionary phases, we study the physical and chemical properties in this range of physical environments. We fit the 69 $\mu$m band and use its width and position to probe the composition and temperature of the crystalline olivine. For 27 sources in the sample, we detected the 69 $\mu$m band of crystalline olivine ($\text{Mg}_{(2-2x)}\text{Fe}_{(2x)}\text{SiO}_{4}$). The 69 $\mu$m band shows that all the sources produce pure forsterite grains containing no iron in their lattice structure. The temperature of the crystalline olivine as indicated by the 69 $\mu$m band, shows that on average the temperature of the crystalline olivine is highest in the group of OH/IR stars and the post-AGB stars with confirmed Keplerian disks. The temperature is lower for the other post-AGB stars and lowest for the planetary nebulae. A couple of the detected 69 $\mu$m bands are broader than those of pure magnesium-rich crystalline olivine, which we show can be due to a temperature gradient in the circumstellar environment of these stars. continued...

The unparalleled photometric data obtained by NASA's Kepler space telescope led to an improved understanding of red giant stars and binary stars. Seismology allows us to constrain the properties of red giants. In addition to eclipsing binaries, eccentric non-eclipsing binaries, exhibiting ellipsoidal modulations, have been detected with Kepler. We aim to study the properties of eccentric binary systems containing a red giant star and derive the parameters of the primary giant component. We apply asteroseismic techniques to determine masses and radii of the primary component of each system. For a selected target, light and radial velocity curve modelling techniques are applied to extract the parameters of the system. The effects of stellar on the binary system are studied. The paper presents the asteroseismic analysis of 18 pulsating red giants in eccentric binary systems, for which masses and radii were constrained. The orbital periods of these systems range from 20 to 440days. From radial velocity measurements we find eccentricities between e=0.2 to 0.76. As a case study we present a detailed analysis of KIC5006817. From seismology we constrain the rotational period of the envelope to be at least 165 d, roughly twice the orbital period. The stellar core rotates 13 times faster than the surface. From the spectrum and radial velocities we expect that the Doppler beaming signal should have a maximum amplitude of 300ppm in the light curve. Through binary modelling, we determine the mass of the secondary component to be 0.29$\pm$0.03\,$M_\odot$. For KIC5006817 we exclude pseudo-synchronous rotation of the red giant with the orbit. The comparison of the results from seismology and modelling of the light curve shows a possible alignment of the rotational and orbital axis at the 2$\sigma$ level. Red giant eccentric systems could be progenitors of cataclysmic variables and hot subdwarf B stars.

Spectra of OH/IR stars show prominent spectral bands of crystalline olivine (Mg$_{(2-2x)}$Fe$_{(2x)}$SiO$_{4}$). To learn more about the timescale of the outflows of OH/IR stars, we study the spectral band of crystalline olivine at 69 {\mu}m. The 69 {\mu}m band is of interest because its width and peak wavelength position are sensitive to the grain temperature and to the exact composition of the crystalline olivine. With Herschel/PACS, we observed the 69 {\mu}m band in the outflow of 14 OH/IR stars. By comparing the crystalline olivine features of our sample with those of model spectra, we determined the size of the outflow and its crystalline olivine abundance. The temperature indicated by the observed 69 {\mu}m bands can only be reproduced by models with a geometrically compact superwind ($R_{\rm{SW}}\lesssim$ 2500 AU = 1400 R$_{*}$).This means that the superwind started less than 1200 years ago (assuming an outflow velocity of 10 km/s). The small amount of mass lost in one superwind and the high progenitor mass of the OH/IR stars introduce a mass loss and thus evolutionary problem for these objects, which has not yet been understood.

The space-missions MOST, CoRoT, and Kepler deliver a huge amount of high-quality photometric data suitable to study numerous pulsating stars. Our ultimate goal is a detection and analysis of an extended sample of Gamma Dor-type pulsating stars with the aim to search for observational evidence of non-uniform period spacings and rotational splittings of gravity modes in main-sequence stars typically twice as massive as the Sun. We applied an automated supervised photometric classification method to select a sample of 69 Gamma Doradus candidate stars. We used an advanced method to extract the Kepler light curves from the pixel data information using custom masks. For 36 of the stars, we obtained high-resolution spectroscopy with the HERMES spectrograph installed at the Mercator telescope. We find that all stars for which spectroscopic estimates of Teff and logg are available fall into the region of the HR diagram where the Gamma Dor and Delta Sct instability strips overlap. The stars cluster in a 700 K window in effective temperature, logg measurements suggest luminosity class IV-V. From the Kepler photometry, we identify 45 Gamma Dor-type pulsators, 14 Gamma Dor/Delta Sct hybrids, and 10 stars which are classified as "possibly Gamma Dor/Delta Sct hybrid pulsators". The results of photometric and spectroscopic classifications according to the type of variability are in perfect agreement. We find a clear correlation between the spectroscopically derived vsini and the frequencies of independent pulsation modes and show that it has nothing to do with rotational modulation of the stars but is related to their stellar pulsations. Our sample and frequency determinations offer a good starting point for seismic modelling of slow to moderately rotating Gamma Dor stars.

Some planetary systems harbour debris disks containing planetesimals such as asteroids and comets. Collisions between such bodies produce small dust particles, the spectral features of which reveal their composition and, hence, that of their parent bodies. A measurement of the composition of olivine crystals has been done for the protoplanetary disk HD100546 and for olivine crystals in the warm inner parts of planetary systems. The latter compares well with the iron-rich olivine in asteroids (x<0.29). In the cold outskirts of the Beta Pictoris system, an analogue to the young Solar System, olivine crystals were detected but their composition remained undetermined, leaving unknown how the composition of the bulk of Solar System cometary olivine grains compares with that of extrasolar comets. Here we report the detection of the 69-micrometre-wavelength band of olivine crystals in the spectrum of Beta Pictoris. Because the disk is optically thin, we can associate the crystals with an extrasolar proto-Kuiper belt a distance of 15-45 astronomical units from the star (one astronomical unit is the Sun-Earth distance), determine their magnesium-rich composition (x=0.01+/-0.001) and show that they make up 3.6+/-1.0 per cent of the total dust mass. These values are strikingly similar to those for the dust emitted by the most primitive comets in the Solar System, even though Beta Pictoris is more massive and more luminous and has a different planetary system architecture.

The broad 30 \mu m feature in carbon stars is commonly attributed to MgS dust particles. However, reproducing the 30 \mu m feature with homogeneous MgS grains would require much more sulfur relative to the solar abundance. Direct gas-phase condensation of MgS occurs at a low efficiency. Precipitation of MgS on SiC precursor grains provides a more efficient formation mechanism, such that the assumption of homogeneous MgS grains may not be correct. Using a Monte Carlo-based radiative transfer code, we aim to model the 30 \mu m feature of the extreme carbon star LL Peg with MgS dust particles. We find that for LL Peg this modeling is insensitive to the unknown MgS optical properties at \lambda < 10 \mu m. When MgS is allowed to be in thermal contact with amorphous carbon and SiC, the amount of MgS required to reproduce the strength of 30 \mu m feature agrees with the solar abundance of sulfur, thereby resolving the reported MgS mass problem. We conclude that MgS is a valid candidate to be the carrier of the 30 \mu m feature when it is part of a composite grain population that has optical properties representative of an ensemble of particle shapes.

We report the discovery of low-amplitude gravity-mode oscillations in the massive binary star V380 Cyg, from 180 d of Kepler custom-aperture space photometry and 5 months of high-resolution high signal-to-noise spectroscopy. The new data are of unprecedented quality and allowed to improve the orbital and fundamental parameters for this binary. The orbital solution was subtracted from the photometric data and led to the detection of periodic intrinsic variability with frequencies of which some are multiples of the orbital frequency and others are not. Spectral disentangling allowed the detection of line-profile variability in the primary. With our discovery of intrinsic variability interpreted as gravity mode oscillations, V380 Cyg becomes an important laboratory for future seismic tuning of the near-core physics in massive B-type stars.

Solid state spectroscopy continues to be an important source of information on the mineralogical composition and physical properties of dust grains both in space and on planetary surfaces. With only a few exceptions, artificially produced or natural terrestrial analog materials, rather than 'real' cosmic dust grains, are the subject of solid state astrophysics. The Jena laboratory has provided a large number of data sets characterizing the UV, optical and infrared properties of such cosmic dust analogs. The present paper highlights recent developments and results achieved in this context, focussing on 'non-standard conditions' such as very low temperatures, very high temperatures and very long wavelengths.

Observations of high-excitation molecular emission lines can greatly increase our understanding of AGB winds, as they trace the innermost regions of the circumstellar envelope. The PACS spectrometer on-board the Herschel Space Telescope, provides for the first time the spectral resolution and sensitivity necessary to trace these lines. We report on the first modelling efforts of a PACS spectral scan for the OH/IR star V669 Cas. Central to our methodology is the consistent treatment of both dust and gas by using a line radiative transfer and a continuum radiative transfer code conjointly. Water emission lines are found to be extremely sensitive to the dust-to-gas ratio, emphasizing the need of consistent modelling for dust and gas.

We present a diagnostic tool to determine the forsterite abundance of the dust ejected by AGB stars. Our method is based on a comparison between the observed strength of spectral bands of forsterite and model calculations. We show that the 11.3 {\mu}m forsterite band is a robust indicator of 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. We apply this method to six high mass-loss rate AGB stars, showing that AGB stars can have forsterite abundances of 12% by mass and higher, which is more than the previously found maximum abundance of 5%.

In this article we present the detection of the 69 {\mu}m band of the crystalline olivine forsterite within the MESS key program of Herschel. We determine the temperature of the forsterite grains by fitting the 69 {\mu}m band.

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%.