Source author record

L. B. F. M. Waters

L. B. F. M. Waters appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.SR astro-ph.GA astro-ph astro-ph.EP astro-ph.IM

Catalog footprint

What is connected

79works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2025arXiv

Multi-band infrared imaging reveals dusty spiral arcs around the binary B[e] star 3 Puppis

3 Puppis is the brightest known B[e] star. Recent work classifies this A-type object as a supergiant, yet the impact of its binarity on the circumstellar environment (CE) remains hard to characterize. To resolve its dusty region at 5-10 mas, we obtained mid-IR interferometric observations with VLTI/MATISSE over 3-12 μm. Because the (u,v) coverage supports imaging, we introduce a statistical interferometric-imaging workflow based on MiRA to generate averaged images: this systematic approach enables the selection of an optimal set of reconstructions, improving the robustness and fidelity of the recovered features. We also use SPARCO, an independent tool well suited to bright central objects embedded in fainter extended emission. Images from both tools in the L, M, and N bands agree and reveal an asymmetric, elongated feature ~17 mas (~10 au at 631 pc) southeast of the star with ~20% density contrast. A second northwest asymmetry and a skewed inner rim are detected. Simple geometric modelling, guided by the MATISSE images, constrains the morphology, location, and flux of the CE and its asymmetries. The images are consistent with earlier VLTI measurements but expose a more complex CE with large-scale clumps in the southeast and northwest parts of the disc. Hydrodynamic modelling indicates that tidal spiral-wake perturbations from the central binary, dynamically excited at Lindblad resonances in the circumbinary disc, best explain the radial extent and curvature of the elongated structures seen in all bands.

preprint2022arXiv

Forming planets around stars with non-solar elemental composition

Stars in the solar neighbourhood have refractory element ratios slightly different from the Sun. It is unclear how much the condensation of solids and thus the composition of planets forming around these stars is affected. We aim to understand the impact of changing the ratios of refractory elements Mg, Si, and Fe within the range observed in solar type stars within 150~pc on the composition of planets forming around them. We use the GGchem code to simulate the condensation of solids in protoplanetary disks with a Minimum Mass Solar Nebula around main sequence G-type stars in the Solar neighbourhood. We extract the stellar elemental composition from the Hypatia database. We find that a lower Mg/Si ratio shifts the condensation sequence from forsterite (Mg$_2$SiO$_4$) and SiO to enstatite (MgSiO$_3$) and quartz (SiO$_2$); a lower Fe/S ratio leads to the formation of FeS and FeS$_2$ and little or no Fe-bearing silicates. Ratios of refractory elements translate directly from the gas phase to the condensed phase for $T\,<\,1000$~K. However, ratios with respect to volatile elements (e.g.\ oxygen and sulphur) in the condensates -- the building blocks of planets -- differ from the original stellar composition. Our study shows that the composition of planets crucially depends on the abundances of the stellar system under investigation. Our results can have important implications for planet interiors, which depend strongly on the degree of oxidation and the sulphur abundance.

preprint2022arXiv

MATISSE, the VLTI mid-infrared imaging spectro-interferometer

Context:Optical interferometry is at a key development stage. ESO's VLTI has established a stable, robust infrastructure for long-baseline interferometry for general astronomical observers. The present second-generation instruments offer a wide wavelength coverage and improved performance. Their sensitivity and measurement accuracy lead to data and images of high reliability. Aims:We have developed MATISSE, the Multi AperTure mid-Infrared SpectroScopic Experiment, to access high resolution imaging in a wide spectral domain and explore topics such: stellar activity and mass loss; planet formation and evolution in the gas and dust disks around young stars; accretion processes around super massive black holes in AGN. Methods:The instrument is a spectro-interferometric imager covering three atmospheric bands (L,M,N) from 2.8 to 13.0 mu, combining four optical beams from the VLTI's telscopes. Its concept, related observing procedure, data reduction and calibration approach are the product of 30 years of instrumental research. The instrument utilizes a multi-axial beam combination that delivers spectrally dispersed fringes. The signal provides the following quantities at several spectral resolutions: photometric flux, coherent fluxes, visibilities, closure phases, wavelength differential visibilities and phases, and aperture-synthesis imaging. Results:We provide an overview of the physical principle of the instrument and its functionalities, the characteristics of the delivered signal, a description of the observing modes and of their performance limits. An ensemble of data and reconstructed images are illustrating the first acquired key observations. Conclusion:The instrument has been in operation at Cerro Paranal, ESO, Chile since 2018, and has been open for science use by the international community since April 2019. The first scientific results are being published now.

preprint2022arXiv

The disk of FU Orionis viewed with MATISSE/VLTI: first interferometric observations in $L$ and $M$ bands

The disk of FU Orionis is marginally resolved with MATISSE, suggesting that the region emitting in the thermal infrared is rather compact. An upper limit of $\sim1.3\pm0.1$ mas (in $L$) can be given for the diameter of the disk region probed in the $L$ band, corresponding to 0.5 au at the adopted Gaia EDR3 distance. This represents the hot, gaseous region of the accretion disk. The $N$-band data indicate that the dusty passive disk is silicate-rich. Only the innermost region of said dusty disk is found to emit strongly in the $N$ band, and it is resolved at an angular size of $\sim5$ mas, which translates to a diameter of about 2 au. The observations therefore place stringent constraints for the outer radius of the inner accretion disk. Dust radiative transfer simulations with RADMC-3D provide adequate fits to the spectral energy distribution from the optical to the submillimeter and to the interferometric observables when opting for an accretion rate $\dot{M}\sim 2\times 10^{-5}\, M_\odot$ yr$^{-1}$ and assuming $M_*=0.6\, M_\odot$. Most importantly, the hot inner accretion disk's outer radius can be fixed at 0.3 au. The outer radius of the dusty disk is placed at 100 au, based on constraints from scattered-light images in the literature. The dust mass contained in the disk is $2.4\times10^{-4}\, M_\odot$, and for a typical gas-to-dust ratio of 100, the total mass in the disk is approximately 0.02 $M_\odot$. We did not find any evidence for a nearby companion in the current interferometric data, and we tentatively explored the case of disk misalignment. For the latter, our modeling results suggest that the disk orientation is similar to that found in previous imaging studies by ALMA. Should there be an asymmetry in the very compact, inner accretion disk, this might be resolved at even smaller spatial scales ($\leq1$ mas).

preprint2022arXiv

The Infrared Database of Extragalactic Observables from Spitzer. -- II. The Database & The Diagnostic Power of Crystalline Silicate Features in Galaxy Spectra

We present the Infrared Database of Extragalactic Observables from Spitzer (IDEOS), a homogeneous, publicly available, database of 77 fitted mid-infrared observables in the 5.4-36um range, comprising measurements for 3335 galaxies observed in the low-resolution staring mode of the Infrared Spectrometer onboard the Spitzer Space Telescope. Among the included observables are PAH fluxes and their equivalent widths, the strength of the 9.8um silicate feature, emission line fluxes, solid-state features, rest frame continuum fluxes, synthetic photometry, and a mid-infrared spectral classification. The IDEOS spectra were selected from the Cornell Atlas of Spitzer-IRS Sources. To our surprise we have detected at a >95% confidence level crystalline silicates in the spectra of 786 IDEOS galaxies. The detections range from single band detections to detections of all fitted crystalline bands (16, 19, 23, 28 and 33um). We find the strength of the crystalline silicate bands to correlate with the amorphous silicate strength, and the change from an emission to an absorption feature to occur at higher obscuration as the wavelength of the crystalline silicate band is longer. These observed characteristics are consistent with an origin for the amorphous and crystalline silicate features in a centrally heated dust geometry, either an edge-on disk or a cocoon. We find the 23 and 33um crystalline silicate bands to be well-suited to classify the obscuration level of galactic nuclei, even in the presence of strong circumnuclear star formation. Based on our detection statistics, we conclude that crystalline silicates are a common component of the interstellar medium of galactic nuclei.

preprint2021arXiv

ATOMIUM: ALMA tracing the origins of molecules in dust forming oxygen rich M-type stars: Motivation, sample, calibration, and initial results

This overview paper presents ATOMIUM, a Large Programme in Cycle 6 with the Atacama Large Millimeter-submillimeter Array (ALMA). The goal of ATOMIUM is to understand the dynamics and the gas phase and dust formation chemistry in the winds of evolved asymptotic giant branch (AGB) and red supergiant (RSG) stars. A more general aim is to identify chemical processes applicable to other astrophysical environments. 17 oxygen-rich AGB and RSG stars spanning a range in (circum)stellar parameters and evolutionary phases were observed in a homogeneous observing strategy allowing for an unambiguous comparison. Data were obtained between 213.83 and 269.71 GHz at high (0.025-0.050 arcsec), medium (0.13-0.24 arcsec), and low (about 1 arcsec) angular resolution. The sensitivity per 1.3 km/s channel was 1.5-5 mJy/beam. 13 molecules were designated as primary molecules in the survey: CO, SiO, AlO, AlOH, TiO, TiO2, HCN, SO, SO2, SiS, CS, H2O, and NaCl. The scientific motivation, survey design, sample properties, data reduction, and an overview of the data products are described; and we highlight one scientific result - the wind kinematics of the ATOMIUM sources. The ATOMIUM sources often have a slow wind acceleration, and a fraction of the gas reaches a velocity which can be up to a factor of two times larger than previously reported terminal velocities assuming isotropic expansion, and the wind kinematic profiles establish that the radial velocity described by the momentum equation for a spherical wind structure cannot capture the complexity of the velocity field. In 15 sources, some molecular transitions other than 12CO v=0 J=2-1 reach a higher outflow velocity, with a spatial emission zone that is often greater than 30 stellar radii, but much less than the extent of CO. Binary interaction with a (sub)stellar companion might (partly) explain the non-monotonic behaviour of the projected velocity field.

preprint2021arXiv

VLTI-MATISSE L- and N-band aperture-synthesis imaging of the unclassified B[e] star FS Canis Majoris

Context: FS Canis Majoris (FS CMa, HD 45677) is an unclassified B[e] star surrounded by an inclined dust disk. The evolutionary stage of FS CMa is still debated. Perpendicular to the circumstellar disk, a bipolar outflow was detected. Infrared aperture-synthesis imaging provides us with a unique opportunity to study the disk structure. Aims: Our aim is to study the intensity distribution of the disk of FS CMa in the mid-infrared L and N bands. Methods: We performed aperture-synthesis imaging of FS CMa with the MATISSE instrument (Multi AperTure mid-Infrared SpectroScopic Experiment) in the low spectral resolution mode to obtain images in the L and N bands. We computed radiative transfer models that reproduce the L- and N-band intensity distributions of the resolved disks. Results: We present L- and N-band aperture-synthesis images of FS CMa reconstructed in the wavelength bands of 3.4-3.8 and 8.6-9.0 micrometer. In the L-band image, the inner rim region of an inclined circumstellar disk and the central object can be seen with a spatial resolution of 2.7 milliarcsec (mas). An inner disk cavity with an angular diameter of 6x12mas is resolved. The L-band disk consists of a bright northwestern (NW) disk region and a much fainter southeastern (SE) region. The images suggest that we are looking at the bright inner wall of the NW disk rim, which is on the far side of the disk. In the N band, only the bright NW disk region is seen. In addition to deriving the inclination and the inner disk radius, fitting the reconstructed brightness distributions via radiative transfer modeling allows one to constrain the innermost disk structure, in particular the shape of the inner disk rim.

preprint2020arXiv

Inner dusty envelope of the AGB stars W~Hydrae, SW\,Virginis, and R~Crateris using SPHERE/ZIMPOL

Stars with initial masses between $\sim0.8$ and 8~$M_\odot$ present copious mass loss during the asymptotic giant branch (AGB) at the end of their lives. The accepted mass-loss mechanism requires radiation pressure acting on dust grains that form in the extended AGB stellar atmospheres. The details of this process are not yet well understood, however. Using the extreme-adaptive-optics imager and polarimeter SPHERE/ZIMPOL, we observed light polarised by grains around W\,Hya, SW\,Vir, and R\,Crt, which have mass-loss rates between 10$^{-7}$ and 10$^{-6}~M_\odot~{\rm yr^{-1}}$. We find the distribution of dust to be asymmetric around the three targets. A biconical morphology is seen for R Crt, with a position angle that is very similar to those inferred from interferometric observations of maser emission and of mid-infrared continuum emission. The cause of the biconical outflow cannot be directly inferred from the ZIMPOL data. The dust grains polarise light more efficiently at 0.65~$μ$m for R\,Crt and SW\,Vir and at 0.82~$μ$m for W\,Hya. This indicates that at the time of the observations, the grains around SW\,Vir and R\,Crt had sizes $< 0.1~μ$m, while those around W\,Hya were larger, with sizes $\gtrsim 0.1~μ$m. The asymmetric distribution of dust around R\,Crt makes the interpretation more uncertain for this star, however. We find that polarised light is produced already from within the visible photosphere of W~Hya, which we reproduce using models with an inner dust shell that is optically thick to scattering. The radial profile of the polarised light observed around W\,Hya reveal a steep dust density profile. We find the wind-acceleration region of W\,Hya to extend to at least $\sim 7~R_\star$, in agreement with theoretical predictions of acceleration up to $\sim 12~R_\star$.

preprint2020arXiv

The asymmetric inner disk of the Herbig Ae star HD 163296 in the eyes of VLTI/MATISSE: evidence for a vortex?

Context. The inner few au region of planet-forming disks is a complex environment. High angular resolution observations have a key role in understanding the disk structure and the dynamical processes at work. Aims. In this study we aim to characterize the mid-infrared brightness distribution of the inner disk of the young intermediate-mass star HD 163296, from VLTI/MATISSE observations. Methods. We use geometric models to fit the data. Our models include a smoothed ring, a flat disk with inner cavity, and a 2D Gaussian. The models can account for disk inclination and for azimuthal asymmetries as well. We also perform numerical hydro-dynamical simulations of the inner edge of the disk. Results. Our modeling reveals a significant brightness asymmetry in the L-band disk emission. The brightness maximum of the asymmetry is located at the NW part of the disk image, nearly at the position angle of the semimajor axis. The surface brightness ratio in the azimuthal variation is $3.5 \pm 0.2$. Comparing our result on the location of the asymmetry with other interferometric measurements, we confirm that the morphology of the $r<0.3$ au disk region is time-variable. We propose that this asymmetric structure, located in or near the inner rim of the dusty disk, orbits the star. For the physical origin of the asymmetry, we tested a hypothesis where a vortex is created by Rossby wave instability, and we find that a unique large scale vortex may be compatible with our data. The half-light radius of the L-band emitting region is $0.33\pm 0.01$ au, the inclination is ${52^\circ}^{+5^\circ}_{-7^\circ}$, and the position angle is $143^\circ \pm 3^\circ$. Our models predict that a non-negligible fraction of the L-band disk emission originates inside the dust sublimation radius for $μ$m-sized grains. Refractory grains or large ($\gtrsim 10\ μ$m-sized) grains could be the origin for this emission.

preprint2016arXiv

A proposed new diagnostic for Herbig disc geometry: FWHM versus J of CO ro-vibrational lines

Aims. CO ro-vibrational lines observed from Herbig group II discs are seen to be often broad while the same lines observed from group I discs are often narrow. This difference is not well understood. In this paper we explore the underlying cause for this difference and provide a pathway for a better understanding of the geometry and structure of the inner discs around Herbig Ae/Be stars. Methods. High-spectral-resolution infrared spectra of CO ro-vibrational emission from six Herbig Ae/Be candidate stars were taken with CRIRES (CRyogenic high-resolution InfraRed Echelle Spectrograph) at the VLT (Very Large Telescope). From these spectra, we produce individual and co-added CO ro-vibrational line profiles. We investigate line profile shape differences, and we explore the FWHM (Full Width Half Maximum) variations with J quantum number in the context of disc geometry. Furthermore, we put our new sources into the context of earlier observed sources to study a large sample. For comparison, we also investigate the FWHM variations with J of modelled CO ro-vibrational lines from two typical disc geometries produced with the thermo-chemical disc modelling code ProDiMo. Results. - abbreviated - Comparing dust and gas inner disc geometries (inferred by SED and CO ro-vibrational emission) for the expanded sample of observed Herbig discs, we find no clear correspondence between the SED groups of the sources and their inner CO radius. Conclusions. The FWHM versus J is a potential new gas diagnostic for the inner disc with e.g. a constant FWHM versus J indicating the presence of a large gas hole or gap. Both models and observations indicate the potential of this new diagnostic. - abbreviated -

preprint2016arXiv

Shadows cast on the transition disk of HD 135344B. Multiwavelength VLT/SPHERE polarimetric differential imaging

The protoplanetary disk around the F-type star HD 135344B (SAO 206462) is in a transition stage and shows many intriguing structures both in scattered light and thermal (sub-)millimeter emission which are possibly related to planet formation processes. We study the morphology and surface brightness of the disk in scattered light to gain insight into the innermost disk regions, the formation of protoplanets, planet-disk interactions traced in the surface and midplane layers, and the dust grain properties of the disk surface. We have carried out high-contrast polarimetric differential imaging (PDI) observations with VLT/SPHERE and obtained polarized scattered light images with ZIMPOL in R- and I-band and with IRDIS in Y- and J-band. The scattered light images reveal with unprecedented angular resolution and sensitivity the spiral arms as well as the 25 au cavity of the disk. Multiple shadow features are discovered on the outer disk with one shadow only being present during the second observation epoch. A positive surface brightness gradient is observed in the stellar irradiation corrected images in southwest direction possibly due to an azimuthally asymmetric perturbation of the temperature and/or surface density by the passing spiral arms. The disk integrated polarized flux, normalized to the stellar flux, shows a positive trend towards longer wavelengths which we attribute to large aggregate dust grains in the disk surface. Part of the the non-azimuthal polarization signal in the Uphi image of the J-band observation could be the result of multiple scattering in the disk. The detected shadow features and their possible variability have the potential to provide insight into the structure of and processes occurring in the innermost disk regions.

preprint2016arXiv

Study of the inner dust envelope and stellar photosphere of the AGB star R Doradus using SPHERE/ZIMPOL

We use high-angular-resolution images obtained with SPHERE/ZIMPOL to study the photosphere, the warm molecular layer, and the inner wind of the close-by oxygen-rich AGB star R Doradus. We present observations in filters V, cntH$α$, and cnt820 and investigate the surface brightness distribution of the star and of the polarised light produced in the inner envelope. Thanks to second-epoch observations in cntH$α$, we are able to see variability on the stellar photosphere. We find that in the first epoch the surface brightness of R Dor is asymmetric in V and cntH$α$, the filters where molecular opacity is stronger, while in cnt820 the surface brightness is closer to being axisymmetric. The second-epoch observations in cntH$α$ show that the morphology of R Dor changes completely in a timespan of 48 days to a more axisymmetric and compact configuration. The polarised intensity is asymmetric in all epochs and varies by between a factor of 2.3 and 3.7 with azimuth for the different images. We fit the radial profile of the polarised intensity using a spherically symmetric model and a parametric description of the dust density profile, $ρ(r)=ρ_\circ r^{-n}$. On average, we find exponents of $- 4.5 \pm 0.5$ that correspond to a much steeper density profile than that of a wind expanding at constant velocity. The dust densities we derive imply an upper limit for the dust-to-gas ratio of $\sim 2\times10^{-4}$ at 5.0 $R_\star$. Given the uncertainties in observations and models, this value is consistent with the minimum values required by wind-driving models for the onset of a wind, of $\sim 3.3\times10^{-4}$. However, if the steep density profile we find extends to larger distances from the star, the dust-to-gas ratio will quickly become too small for the wind of R Dor to be driven by the grains that produce the scattered light.

preprint2016arXiv

The abundance and thermal history of water ice in the disk surrounding HD142527 from the DIGIT Herschel Key Program

The presence or absence of ice in protoplanetary disks is of great importance for the formation of planets. By enhancing the solid surface density and increasing the sticking efficiency, ice catalyzes the rapid formation of planetesimals and decreases the time scale for giant planet core accretion. Aims: In this paper we analyse the composition of the outer disk around the Herbig star HD~142527. We focus on the composition of the water ice, but also analyse the abundances of previously proposed minerals. Methods: We present new Herschel far infrared spectra and a re-reduction of archival data from the Infrared Space Observatory (ISO). We model the disk using full 3D radiative transfer to obtain the disk structure. Also, we use an optically thin analysis of the outer disk spectrum to obtain firm constraints on the composition of the dust component. Results: The water ice in the disk around HD~142527 contains a large reservoir of crystalline water ice. We determine the local abundance of water ice in the outer disk (i.e. beyond 130\,AU). The re-reduced ISO spectrum differs significantly from that previously published, but matches the new Herschel spectrum at their common wavelength range. In particular, we do not detect any significant contribution from carbonates or hydrous silicates, in contrast to earlier claims. Conclusions: The amount of water ice detected in the outer disk requires $\sim80\,$\% of the oxygen atoms. This is comparable to the water ice abundance in the outer solar system, in comets and in dense interstellar clouds. The water ice is highly crystalline while the temperatures where we detect it are too low to crystallize the water on relevant time scales. We discuss the implications of this finding.

preprint2016arXiv

Variability in the CO ro-vibrational lines from HD163296

We present for the first time a direct comparison of multi-epoch (2001-2002 and 2012) CO ro-vibrational emission lines from HD163296. We find that both the line shapes and the FWHM (Full Width Half Maximum) differ between these two epochs. The FWHM of the median observed line profiles are 10-25 km/s larger in the earlier epoch, and confirmed double peaks are only present in high J lines from 2001-2002. The line wings of individual transitions are similar in the two epochs making an additional central component in the later epoch a likely explanation for the single peaks and the lower FWHM. Variations in NIR brightness have been reported and could be linked to the observed variations. Additionally, we use the thermo chemical disc code ProDiMo to compare for the first time the line shapes, peak separations, FWHM, and line fluxes, to those observed. The ProDiMo model reproduces the peak separations, and low and mid J line fluxes well. The FWHM however, are over predicted and high J line fluxes are under predicted. We propose that a variable non-Keplerian component of the CO ro-vibrational emission, such as a disc wind or an episodic accretion funnel, is causing the difference between the two data sets collected at different epochs, and between model and observations. Additional CO ro-vibrational line detections (with CRIRES/VLT or NIRSPEC/Keck) or [Ne II] line observations with VISIR/VLT could help to clarify the cause of the variability.

preprint2015arXiv

Consistent dust and gas models for protoplanetary disks. I. Disk shape, dust settling, opacities, and PAHs

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. We propose new standard dust opacities for disk models, we present a simplified treatment of PAHs sufficient to reproduce the PAH emission features, and we suggest using a simple treatment of dust settling. We roughly adjust parameters to obtain a model that predicts typical Class II T Tauri star continuum and line observations. We systematically study the impact of each model parameter (disk mass, disk extension and shape, dust settling, dust size and opacity, gas/dust ratio, etc.) on all continuum and line observables, in particular on the SED, mm-slope, continuum visibilities, and emission lines including [OI] 63um, high-J CO lines, (sub-)mm CO isotopologue lines, and CO fundamental ro-vibrational lines. We find that evolved dust properties (large grains) often needed to fit the SED, have important consequences for disk chemistry and heating/cooling balance, leading to stronger emission lines in general. Strong dust settling and missing disk flaring have similar effects on continuum observations, but opposite effects on far-IR gas emission lines. PAH molecules can shield the gas from stellar UV radiation because of their strong absorption and negligible scattering opacities. The observable millimetre-slope of the SED can become significantly more gentle in the case of cold disk midplanes, which we find regularly in our T Tauri models. We propose to use line observations of robust chemical tracers of the gas, such as O, CO, and H2, as additional constraints to determine some key properties of the disks, such as disk shape and mass, opacities, and the dust/gas ratio, by simultaneously fitting continuum and line observations.

preprint2015arXiv

Dusty wind of W Hya. Multi-wavelength modelling of the present-day and recent mass-loss

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\,$μ$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.

preprint2015arXiv

Herschel observations of extreme OH/IR stars - the isotopic ratios of oxygen as a sign-post for the stellar mass

Aim: The late stages of stellar evolution are mainly governed by the mass of the stars. Low- and intermediate-mass stars lose copious amounts of mass during the asymptotic giant branch (AGB) which obscure the central star making it difficult to study the stellar spectra and determine the stellar mass. In this study, we present observational data that can be used to determine lower limits to the stellar mass. Method: Spectra of nine heavily reddened AGB stars taken by the Herschel Space Observatory display numerous molecular emission lines. The strongest emission lines are due to H2O. We search for the presence of isotopologues of H2O in these objects. Result: We detected the 16O and 17O isotopologues of water in these stars, but lines due to H2^{18}O are absent. The lack of 18O is predicted by a scenario where the star has undergone hot-bottom burning which preferentially destroys 18O relative to 16O and 17O. From stellar evolution calculations, this process is thought to occur when the stellar mass is above 5 Msun for solar metallicity. Hence, observations of different isotopologues of H2O can be used to help determine the lower limit to the initial stellar mass. Conclusion: From our observations, we deduce that these extreme OH/IR stars are intermediate-mass stars with masses of >= 5 Msun. Their high mass-loss rates of ~ 1.0e-4 Msun/yr may affect the enrichment of the interstellar medium and the overall chemical evolution of our Galaxy.

preprint2015arXiv

Micron-sized forsterite grains in the pre-planetary nebula of IRAS 17150-3224 - Searching for clues on the mysterious evolution of massive AGB stars

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.

preprint2015arXiv

The Mid-Infrared Instrument for the James Webb Space Telescope, VI: The Medium Resolution Spectrometer

We describe the design and performance of the Medium Resolution Spectrometer (MRS) for the JWST-MIRI instrument. The MRS incorporates four coaxial spectral channels in a compact opto-mechanical layout that generates spectral images over fields of view up to 7.7 X 7.7 arcseconds in extent and at spectral resolving powers ranging from 1,300 to 3,700. Each channel includes an all-reflective integral field unit (IFU): an 'image slicer' that reformats the input field for presentation to a grating spectrometer. Two 1024 X 1024 focal plane arrays record the output spectral images with an instantaneous spectral coverage of approximately one third of the full wavelength range of each channel. The full 5 to 28.5 micron spectrum is then obtained by making three exposures using gratings and pass-band-determining filters that are selected using just two three-position mechanisms. The expected on-orbit optical performance is presented, based on testing of the MIRI Flight Model and including spectral and spatial coverage and resolution. The point spread function of the reconstructed images is shown to be diffraction limited and the optical transmission is shown to be consistent with the design expectations.

preprint2015arXiv

The structure of disks around intermediate-mass young stars from mid-infrared interferometry. Evidence for a population of group II disks with gaps

The disks around Herbig Ae/Be stars are commonly divided into group I and group II based on their far-infrared spectral energy distribution, and the common interpretation for that is flared and flat disks. Recent observations suggest that many flaring disks have gaps, whereas flat disks are thought to be gapless. The different groups of objects can be expected to have different structural signatures in high-angular-resolution data. Over the past 10 years, the MIDI instrument on the Very Large Telescope Interferometer has collected observations of several tens of protoplanetary disks. We model the large set of observations with simple geometric models. A population of radiative-transfer models is synthesized for interpreting the mid-infrared signatures. Objects with similar luminosities show very different disk sizes in the mid-infrared. Restricting to the young objects of intermediate mass, we confirm that most group I disks are in agreement with being transitional. We find that several group II objects have mid-infrared sizes and colors overlapping with sources classified as group I, transition disks. This suggests that these sources have gaps, which has been demonstrated for a subset of them. This may point to an intermediate population between gapless and transition disks. Flat disks with gaps are most likely descendants of flat disks without gaps. Gaps, potentially related to the formation of massive bodies, may therefore even develop in disks in a far stage of grain growth and settling. The evolutionary implications of this new population could be twofold. Either gapped flat disks form a separate population of evolved disks, or some of them may further evolve into flaring disks with large gaps. The latter transformation may be governed by the interaction with a massive planet, carving a large gap and dynamically exciting the grain population in the disk.

preprint2014arXiv

Dust composition and mass-loss return from the luminous blue variable R71 in the LMC

We present an analysis of mid-and far-infrared (IR) spectrum and spectral energy distribution (SED) of the LBV R71 in the LMC.This work aims to understand the overall contribution of high-mass LBVs to the total dust-mass budget of the interstellar medium (ISM) of the LMC and compare this with the contribution from low-mass asymptotic giant branch (AGB) stars. As a case study, we analyze the SED of R71. We compiled all the available photometric and spectroscopic observational fluxes from various telescopes for a wide wavelength range (0.36 -- 250\,$μ$m). We determined the dust composition from the spectroscopic data, and derived the ejected dust mass, dust mass-loss rate, and other dust shell properties by modeling the SED of R71. We noted nine spectral features in the dust shell of R71 by analyzing Spitzer spectroscopic data. Among these, we identified three new crystalline silicate features. We computed our model spectrum by using 3D radiative transfer code MCMax. Our model calculation shows that dust is dominated by amorphous silicates, with some crystalline silicates, metallic iron, and a very tiny amount of polycyclic aromatic hydrocarbon (PAH) molecules. The presence of both silicates and PAHs indicates that the dust has a mixed chemistry. We derived a dust mass of 0.01 M$_\odot$, from which we arrive at a total ejected mass of $\approx$ 5 M$_\odot$. This implies a time-averaged dust mass-loss rate of 2.5$\times$10$^{-6}$ M$_\odot$\,yr$^{-1}$ with an explosion about 4000 years ago. We assume that the other five confirmed dusty LBVs in the LMC loose mass at a similar rate, and estimate the total contribution to the mass budget of the LMC to be $\approx$ 10$^{-5}$ M$_\odot$\,yr$^{-1}$, which is comparable to the contribution by all the AGB stars in the LMC. Based on our analysis on R71, we speculate that LBVs as a class may be an important dust source in the ISM of the LMC.

preprint2014arXiv

Herschel/PACS observations of the 69 $μm$ band of crystalline olivine around evolved stars

We present 48 Herschel/PACS spectra of evolved stars in the wavelength range of 67-72 $μ$m. This wavelength range covers the 69 $μ$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 $μ$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 $μ$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 $μ$m band of crystalline olivine ($\text{Mg}_{(2-2x)}\text{Fe}_{(2x)}\text{SiO}_{4}$). The 69 $μ$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 $μ$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 $μ$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...

preprint2014arXiv

Location and sizes of forsterite grains in protoplanetary disks: interpretation from the Herschel DIGIT programme

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.

preprint2014arXiv

Polycyclic aromatic hydrocarbon ionization as a tracer of gas flows through protoplanetary disk gaps

Planet-forming disks of gas and dust around young stars contain polycyclic aromatic hydrocarbons (PAHs). We aim to characterize how the charge state of PAHs can be used as a probe of flows of gas through protoplanetary gaps. In this context, our goal is to understand the PAH spectra of four transitional disks. In addition, we want to explain the observed correlation between PAH ionization (traced by the 6.2/11.3 feature ratio) and the disk mass (traced by the 1.3 mm luminosity). We implement a model to calculate the charge state of PAHs in the radiative transfer code MCMax. The emission spectra and ionization balance are calculated. A benchmark modeling grid is presented that shows how PAH ionization and luminosity behave as a function of star and disk properties. The PAH ionization is most sensitive to ultraviolet (UV) radiation and the electron density. In optically thick disks, where the UV field is low and the electron density is high, PAHs are predominantly neutral. Ionized PAHs trace low-density optically thin disk regions where the UV field is high and the electron density is low. Such regions are characteristic of gas flows through the gaps of transitional disks. We demonstrate that fitting the PAH spectra of four transitional disks requires a contribution of ionized PAHs in gas flows through the gap. The PAH spectra of transitional disks can be understood as superpositions of neutral and ionized PAHs. For HD97048, neutral PAHs in the optically thick disk dominate the spectrum. In the cases of HD169142, HD135344B and Oph IRS 48, small amounts of ionized PAHs located in the gas flows through the gap are strong contributors to the total PAH luminosity. The observed trend between the disk mass and PAH ionization may imply that lower-mass disks have larger gaps. Ionized PAHs in gas flows through these gaps contribute strongly to their spectra.

preprint2014arXiv

Surprising detection of an equatorial dust lane on the AGB star IRC+10216

Understanding the formation of planetary nebulae remains elusive because in the preceding asymtotic giant branch (AGB) phase these stars are heavily enshrouded in an optically thick dusty envelope. To further understand the morphology of the circumstellar environments of AGB stars we observe the closest carbon-rich AGB star IRC+10216 in scattered light. When imaged in scattered light at optical wavelengths, IRC+10216 surprisingly shows a narrow equatorial density enhancement, in contrast to the large-scale spherical rings that have been imaged much further out. We use radiative transfer models to interpret this structure in terms of two models: firstly, an equatorial density enhancement, commonly observed in the more evolved post-AGB stars, and secondly, in terms of a dust rings model, where a local enhancement of mass-loss creates a spiral ring as the star rotates. We conclude that both models can be used to reproduce the dark lane in the scattered light images, which is caused by an equatorially density enhancement formed by dense dust rather than a bipolar outflow as previously thought. We are unable to place constraints on the formation of the equatorial density enhancement by a binary system.

preprint2014arXiv

The structure of disks around Herbig Ae/Be stars as traced by CO ro-vibrational emission

We study the emission and absorption of CO ro-vibrational lines in the spectra of intermediate mass pre-main-sequence stars with the aim to determine both the spatial distribution of the CO gas and its physical properties. We also aim to correlate CO emission properties with disk geometry. Using high-resolution spectra containing fundamental and first overtone CO ro-vibrational emission, observed with CRIRES on the VLT, we probe the physical properties of the circumstellar gas by studying its kinematics and excitation conditions. We detect and spectrally resolve CO fundamental ro-vibrational emission in 12 of the 13 stars observed, and in two cases in absorption. Keeping in mind that we studied a limited sample, we find that the physical properties and spatial distribution of the CO gas correlate with disk geometry. Flaring disks show highly excited CO fundamental emission up to v$_u$ = 5, while self-shadowed disks show CO emission that is not as highly excited. Rotational temperatures range between ~250-2000 K. The $^{13}$CO rotational temperatures are lower than those of $^{12}$CO. The vibrational temperatures in self-shadowed disks are similar to or slightly below the rotational temperatures, suggesting that thermal excitation or IR pumping is important in these lines. In flaring disks the vibrational temperatures reach as high as 6000 K, suggesting fluorescent pumping. Using a simple kinematic model we show that the CO inner radius of the emitting region is $\approx$10 au for flaring disks and $\leq$ 1 au for self-shadowed disks. (abstract truncated due to 1920 character limit for arXiv)

preprint2014arXiv

The wind of W Hya as seen by Herschel - I. The CO envelope

Asymptotic giant branch (AGB) stars lose their envelopes by means of a stellar wind whose driving mechanism is not understood well. Characterizing the composition and thermal and dynamical structure of the outflow provides constraints that are essential for understanding AGB evolution, including the rate of mass loss and isotopic ratios. We modeled the CO emission from the wind of the low mass-loss rate oxygen-rich AGB star W Hya using data obtained by the HIFI, PACS, and SPIRE instruments onboard the Herschel Space Observatory and ground-based telescopes. $^{12}$CO and $^{13}$CO lines are used to constrain the intrinsic $^{12}$C/$^{13}$C ratio from resolved HIFI lines. The acceleration of the outflow up to about 5.5 km/s is quite slow and can be represented by a beta-type velocity law with index 5. Beyond this point, acceleration up the terminal velocity of 7 km/s is faster. Using the J=10-9, 9-8, and 6-5 transitions, we find an intrinsic $^{12}$C/$^{13}$C ratio of $18\pm10$ for W Hya, where the error bar is mostly due to uncertainties in the $^{12}$CO abundance and the stellar flux around 4.6 $μ$m. To match the low-excitation CO lines, these molecules need to be photo-dissociated at about 500 stellar radii. The radial dust emission intensity profile measured by PACS images at 70 $μ$m shows substantially stronger emission than our model predicts beyond 20 arcsec. The initial slow acceleration of the wind implies inefficient wind driving in the lower part of the envelope. The final injection of momentum in the wind might be the result of an increase in the opacity thanks to the late condensation of dust species. The derived intrinsic isotopologue ratio for W Hya is consistent with values set by the first dredge-up and suggestive of an initial mass of 2 M$_\odot$ or more. However, the uncertainty in the main-sequence mass derived based on this isotopologic ratio is large.

preprint2014arXiv

The wind of W Hya as seen by Herschel. II. The molecular envelope of W Hya

The evolution of low- and intermediate-mass stars on the asymptotic giant branch (AGB) is mainly controlled by the rate at which these stars lose mass in a stellar wind. Understanding the driving mechanism and strength of the stellar winds of AGB stars and the processes enriching their surfaces with products of nucleosynthesis are paramount to constraining AGB evolution and predicting the chemical evolution of galaxies. In a previous paper we have constrained the structure of the outflowing envelope of W Hya using spectral lines of the $^{12}$CO molecule. Here we broaden this study by modelling an extensive set of H$_{2}$O and $^{28}$SiO lines observed by the three instruments on board Herschel using a state-of-the-art molecular excitation and radiative transfer code. The oxygen isotopic ratios and the $^{28}$SiO abundance profile can be connected to the initial stellar mass and to crucial aspects of dust formation at the base of the stellar wind, respectively. The modelling of H$_{2}$O and $^{28}$SiO confirms the properties of the envelope model of W Hya derived from $^{12}$CO lines. We find an H$_2$O ortho-to-para ratio of 2.5\,$^{+2.5}_{-1.0}$, consistent with what is expected for an AGB wind. The O$^{16}$/O$^{17}$ ratio indicates that W Hya has an initial mass of about 1.5 M$_\odot$. Although the ortho- and para-H$_{2}$O lines observed by HIFI appear to trace gas of slightly different physical properties, a turbulence velocity of $0.7\pm0.1$ km s$^{-1}$ fits the HIFI lines of both spin isomers and those of $^{28}$SiO well. The ortho- and para-H$_2^{16}$O and $^{28}$SiO abundances relative to H$_{2}$ are $(6^{+3}_{-2}) \times 10^{-4}$, $(3^{+2}_{-1}) \times 10^{-4}$, and $(3.3\pm 0.8)\times 10^{-5}$, respectively. Assuming a solar silicon-to-carbon ratio, the $^{28}$SiO line emission model is consistent with about one-third of the silicon atoms being locked up in dust particles.

preprint2013arXiv

CO ro-vibrational lines in HD100546: A search for disc asymmetries and the role of fluorescence

We have studied the emission of CO ro-vibrational lines in the disc around the Herbig Be star HD100546 with the final goal of using these lines as a diagnostic to understand inner disc structure in the context of planet formation. High-resolution IR spectra of CO ro-vibrational emission at eight different position angles were taken with CRIRES at the VLT. From these spectra flux tables, CO ro-vibrational line profiles, and population diagrams were produced. We have investigated variations in the line profile shapes and line strengths as a function of slit position angle. We used the thermochemical disc modelling code ProDiMo based on the chemistry, radiation field, and temperature structure of a previously published model for HD100546. Comparing observations and the model, we investigated the possibility of disc asymmetries, the excitation mechanism (UV fluorescence), the geometry, and physical conditions of the inner disc. The observed CO ro-vibrational lines are largely emitted from the inner rim of the outer disc at 10-13 AU. The line shapes are similar for all v levels and line fluxes from all vibrational levels vary only within one order of magnitude. All line profile asymmetries and variations can be explained with a symmetric disc model to which a slit correction and pointing offset is applied. Because the angular size of the CO emitting region (10-13 AU) and the slit width are comparable the line profiles are very sensitive to the placing of the slit. The model reproduces the line shapes and the fluxes of the v=1-0 lines as well as the spatial extent of the CO ro-vibrational emission. It does not reproduce the observed band ratios of 0.5-0.2 with higher vibrational bands. We find that lower gas volume densities at the surface of the inner rim of the outer disc can make the fluorescence pumping more effcient and reproduce the observed band ratios.

preprint2013arXiv

HIFISTARS Herschel/HIFI observations of VY Canis Majoris. Molecular-line inventory of the envelope around the largest known star

To gain insight into the physical conditions and kinematics of the warm (100-1000 K) gas around the red hyper-giant VY CMa, we performed sensitive high spectral resolution observations of molecular lines in the sub-mm/FIR using the HIFI instrument of the Herschel Space Observatory. We observed CO, H2O, and other molecular species, sampling excitation energies from a few tens to a few thousand K. These observations are part of the Herschel Guaranteed Time Key Program HIFISTARS. We detected the J=6-5, J=10-9, and J=16-15 lines of 12CO and 13CO at about 100, 300, and 750K above the ground state (and the 13CO J=9-8 line). These lines are crucial for improving the modelling of the internal layers of the envelope around VY CMa. We also detected 27 lines of H2O and its isotopomers, and 96 lines of species such as NH3, SiO, SO, SO2 HCN, OH and others, some of them originating from vibrationally excited levels. Three lines were not unambiguously assigned. Our observations confirm that VY CMa's envelope must consist of two or more detached components. The molecular excitation in the outer layers is significantly lower than in the inner ones, resulting in strong self-absorbed profiles in molecular lines that are optically thick in this outer envelope, for instance, low-lying lines of H2O. Except for the most abundant species, CO and H2O, most of the molecular emission detected at these sub-mm/FIR wavelengths arise from the central parts of the envelope. The spectrum of VY CMa is very prominent in vibrationally excited lines, which are caused by the strong IR pumping present in the central regions. Compared with envelopes of other massive evolved stars, VY CMa's emission is particularly strong in these vibrationally excited lines, as well as in the emission from less abundant species such as H13CN, SO, and NH3.

preprint2013arXiv

Identifying gaps in flaring Herbig Ae/Be disks using spatially resolved mid-infrared imaging. Are all group I disks transitional?

*Context The evolution of young massive protoplanetary disks toward planetary systems is expected to include the formation of gaps and the depletion of dust and gas. *Aims A special group of flaring disks around Herbig Ae/Be stars do not show prominent silicate emission features. We focus our attention on four key Herbig Ae/Be stars to understand the structural properties responsible for the absence of silicate feature emission. *Methods We investigate Q- and N-band images taken with Subaru/COMICS, Gemini South/T-ReCS and VLT/VISIR. Our radiative transfer modeling solutions require a separation of inner- and outer- disks by a large gap. From this we characterize the radial density structure of dust and PAHs in the disk. *Results The inner edge of the outer disk has a high surface brightness and a typical temperature between ~100-150 K and therefore dominates the emission in the Q-band. We derive radii of the inner edge of the outer disk of 34, 23, 30 and 63 AU for HD97048, HD169142, HD135344B and Oph IRS 48 respectively. For HD97048 this is the first detection of a disk gap. The continuum emission in the N-band is not due to emission in the wings of PAHs. This continuum emission can be due to VSGs or to thermal emission from the inner disk. We find that PAH emission is not always dominated by PAHs on the surface of the outer disk. *Conclusions. The absence of silicate emission features is due to the presence of large gaps in the critical temperature regime. Many, if not all Herbig disks with Spectral Energy Distribution (SED) classification `group I' are disks with large gaps and can be characterized as (pre-) transitional. An evolutionary path from the observed group I to the observed group II sources seems no longer likely. Instead, both might derive from a common ancestor.

preprint2013arXiv

OH/IR stars and their superwinds as observed by the Herschel Space Observatory

Aim : In order to study the history of mass loss in extreme OH/IR stars, we observed a number of these objects using CO as a tracer of the density and temperature structure of their circumstellar envelopes. Method : Combining CO observations from the Herschel Space Observatory with those from the ground, we trace mass loss rates as a function of radius in five extreme OH/IR stars. Using radiative transfer modelling, we modelled the dusty envelope as well as the CO emission. The high-rotational transitions of CO indicate that they originate in a dense superwind region close to the star while the lower transitions tend to come from a more tenuous outer wind which is a result of the mass loss since the early AGB phase. Result : The models of the circumstellar envelopes around these stars suggest that they have entered a superwind phase in the past 200 - 500 years. The low 18O/17O (~ 0.1 compared to the solar abundance ratio of ~ 5) and 12C/13C (3-30 cf. the solar value of 89) ratios derived from our study support the idea that these objects have undergone hot-bottom burning and hence that they are massive M >= 5 solar-mass AGB stars.

preprint2013arXiv

RCW36: characterizing the outcome of massive star formation

Massive stars play a dominant role in the process of clustered star formation, with their feedback into the molecular cloud through ionizing radiation, stellar winds and outflows. The formation process of massive stars is poorly constrained because of their scarcity, the short formation timescale and obscuration. By obtaining a census of the newly formed stellar population, the star formation history of the young cluster and the role of the massive stars within it can be unraveled. We aim to reconstruct the formation history of the young stellar population of the massive star-forming region RCW 36. We study several dozens of individual objects, both photometrically and spectroscopically, look for signs of multiple generations of young stars and investigate the role of the massive stars in this process. We obtain a census of the physical parameters and evolutionary status of the young stellar population. Using a combination of near-infrared photometry and spectroscopy we estimate ages and masses of individual objects. We identify the population of embedded young stellar objects (YSO) by their infrared colors and emission line spectra. RCW 36 harbors a stellar population of massive and intermediate-mass stars located around the center of the cluster. Class 0/I and II sources are found throughout the cluster. The central population has a median age of 1.1 +/- 0.6 Myr. Of the stars which could be classified, the most massive ones are situated in the center of the cluster. The central cluster is surrounded by filamentary cloud structures; within these, some embedded and accreting YSOs are found. Our age determination is consistent with the filamentary structures having been shaped by the ionizing radiation and stellar winds of the central massive stars. The formation of a new generation of stars is ongoing, as demonstrated by the presence of embedded protostellar clumps, and two exposed jets.

preprint2013arXiv

The 69 micron forsterite band in spectra of protoplanetary disks - Results from the Herschel DIGIT programme

Context: We have analysed Herschel-PACS spectra of 32 circumstellar disks around Herbig Ae/Be and T-Tauri stars obtained within the Herschel key programme DIGIT. In this paper we focus on the 69mu emission band of the crystalline silicate forsterite. Aims: This work provides an overview of the 69mu forsterite bands in the DIGIT sample. We aim to derive the temperature and composition of the forsterite grains. With this information, constraints can be placed on the spatial distribution of the forsterite in the disk and its formation history. Methods: Position and shape of the 69mu band are used to derive the temperature and composition of the dust by comparison to laboratory spectra of that band. We combine our data with existing Spitzer IRS spectra to compare the presence and strength of the 69mu band to the forsterite bands at shorter wavelengths. Results: A total of 32 sources have been observed, 8 of them show a 69mu emission band that can be attributed to forsterite. With the exception of the T-Tauri star AS205, all of the detections are for disks associated with Herbig Ae/Be stars. Most of the forsterite grains that give rise to the 69mu bands are warm (~100-200 K) and iron-poor (less than ~2% iron). Only AB-Aur requires approximately 3-4% of iron. Conclusions: Our findings support the hypothesis that the forsterite grains form through an equilibrium condensation process at high temperatures. The connection between the strength of the 69 and 33mu bands shows that at least part of the emission in these bands originates from the same dust grains. Further, any model that explains the PACS and the Spitzer IRS observations must take the effects of a wavelength dependent optical depth into account. We find indications of a correlation of the detection rate of the 69mu band with the spectral type of the host stars. However, our sample is too small to obtain a definitive result.

preprint2013arXiv

The color dependent morphology of the post-AGB star HD161796

Context. Many protoplanetary nebulae show strong asymmetries in their surrounding shell, pointing to asymmetries during the mass loss phase. Questions concerning the origin and the onset of deviations from spherical symmetry are important for our understanding of the evolution of these objects. Here we focus on the circumstellar shell of the post-AGB star HD 161796. Aims. We aim at detecting signatures of an aspherical outflow, as well as to derive the properties of it. Methods. We use the imaging polarimeter ExPo (the extreme polarimeter), a visitor instrument at the William Herschel Telescope, to accurately image the dust shell surrounding HD 161796 in various wavelength filters. Imaging polarimetry allows us to separate the faint, polarized, light from circumstellar material from the bright, unpolarized, light from the central star. Results. The shell around HD 161796 is highly aspherical. A clear signature of an equatorial density enhancement can be seen. This structure is optically thick at short wavelengths and changes its appearance to optically thin at longer wavelengths. In the classification of the two different appearances of planetary nebulae from HST images it changes from being classified as DUPLEX at short wavelengths to SOLE at longer wavelengths. This strengthens the interpretation that these two appearances are manifestations of the same physical structure. Furthermore, we find that the central star is hotter than often assumed and the relatively high observed reddening is due to circumstellar rather than interstellar extinction.

preprint2013arXiv

The problematically short superwind of OH/IR stars - Probing the outflow with the 69 μm spectral band of forsterite

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 μm. The 69 μ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 μ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 μ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.

preprint2013arXiv

The VLTI/MIDI survey of massive young stellar objects - Sounding the inner regions around intermediate- and high-mass young stars using mid-infrared interferometry

We aim to characterize the distribution and composition of circumstellar material around young massive stars, and to investigate exactly which physical structures in these objects are probed by long-baseline mid-infrared interferometric observations. We used the two-telescope interferometric instrument MIDI of the Very Large Telescope Interferometer of the European Southern Observatory to observe a sample of 24 intermediate- and high-mass young stellar objects in the N band (8-13 micron). We had successful fringe detections for 20 objects, and present spectrally-resolved correlated fluxes and visibility levels for projected baselines of up to 128 m. We fit the visibilities with geometric models to derive the sizes of the emitting regions, as well as the orientation and elongation of the circumstellar material. Fourteen objects in the sample show the 10 micron silicate feature in absorption in the total and correlated flux spectra. For 13 of these objects, we were able to fit the correlated flux spectra with a simple absorption model, allowing us to constrain the composition and absorptive properties of the circumstellar material. Nearly all of the massive young stellar objects observed show significant deviations from spherical symmetry at mid-infrared wavelengths. In general, the mid-infrared emission can trace both disks and outflows, and in many cases it may be difficult to disentangle these components on the basis of interferometric data alone, because of the sparse spatial frequency coverage normally provided by current long-baseline interferometers. For the majority of the objects in this sample, the absorption occurs on spatial scales larger than those probed by MIDI. Finally, the physical extent of the mid-infrared emission around these sources is correlated with the total luminosity, albeit with significant scatter.

preprint2012arXiv

Comet-like mineralogy of olivine crystals in an extrasolar proto-Kuiper belt

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.

preprint2012arXiv

Direct imaging of a massive dust cloud around R Coronae Borealis

We present recent polarimetric images of the highly variable star R CrB using ExPo and archival WFPC2 images from the HST. We observed R CrB during its current dramatic minimum where it decreased more than 9 mag due to the formation of an obscuring dust cloud. Since the dust cloud is only in the line-of-sight, it mimics a coronograph allowing the imaging of the star's circumstellar environment. Our polarimetric observations surprisingly show another scattering dust cloud at approximately 1.3" or 2000 AU from the star. We find that to obtain a decrease in the stellar light of 9 mag and with 30% of the light being reemitted at infrared wavelengths (from R CrB's SED) the grains in R CrB's circumstellar environment must have a very low albedo of approximately 0.07%. We show that the properties of the dust clouds formed around R CrB are best fitted using a combination of two distinct populations of grains size. The first are the extremely small 5 nm grains, formed in the low density continuous wind, and the second population of large grains (~0.14 μm) which are found in the ejected dust clouds. The observed scattering cloud, not only contains such large grains, but is exceptionally massive compared to the average cloud.

preprint2012arXiv

Herschel/HIFI observations of red supergiants and yellow hypergiants: I. Molecular inventory

Red supergiant stars (RSGs) and yellow hypergiant stars (YHGs) are believed to be the high-mass counterparts of stars in the AGB and early post-AGB phases. We study the mass-loss in the post main-sequence evolution of massive stars, through the properties of their envelopes in the intermediate and warm gas layers. These are the regions where the acceleration of the gas takes place and the most recent mass-loss episodes can be seen. We used the HIFI instrument on-board the Herschel Space Observatory to observe sub-mm and FIR transitions of CO, water, and their isotopologues in a sample of two RSGs (NML Cyg and Betelgeuse) and two YHGs (IRC+10420 and AFGL 2343) stars. We present an inventory of the detected lines and analyse the information revealed by their spectral profiles. On the basis of the results presented in an earlier study, we model the CO and 13CO emission in IRC+10420 and compare it to a set of lines ranging from the mm, to the FIR. Red supergiants have stronger high-excitation lines than the YHGs, indicating that they harbour dense and hot inner shells contributing to these transitions. Consequently, these high-J lines in RSGs originate from acceleration layers that have not yet reached the circumstellar terminal velocity and have narrower profiles than their flat-topped lower-J counterparts. The YHGs tend to lack this inner component, in line with the picture of detached, hollow envelopes derived from studies at longer wavelengths. NH3 is only detected in two sources (NML Cyg, IRC+10420), which are also observed to be the strongest water-line emitters of the studied sample. In contrast, OH is detected in all sources and does not seem to correlate with the water line intensities. We show that the IRC+10420 model derived solely from mm low-J CO transitions is capable of reproducing the high-J transitions when the temperature in the inner shell is simply lowered by about 30%.

preprint2012arXiv

Mid-IR imaging of the transitional disk of HD169142: Measuring the size of the gap

The disk around the Herbig Ae star HD\,169142 was imaged and resolved at 18.8 and 24.5\,$μ$m using Subaru/COMICS. We interpret the observations using a 2D radiative transfer model and find evidence for the presence of a large gap. The MIR images trace dust that emits at the onset of the strong rise in the spectral energy distribution (SED) at 20\,$μ$m, therefore are very sensitive to the location and characteristics of the inner wall of the outer disk and its dust. We determine the location of the wall to be 23$^{+3}_{-5}$\,AU from the star. An extra component of hot dust must exist close to the star. We find that a hydrostatic optically thick inner disk does not produce enough flux in the NIR and an optically thin geometrically thick component is our solution to fit the SED. Considering the recent findings of gaps and holes in a number of Herbig Ae/Be group I disks, we suggest that such disk structures may be common in group I sources. Classification as group I should be considered a support for classification as a transitional disk, though improved imaging surveys are needed to support this speculation.

preprint2012arXiv

Observational evidence for composite grains in an AGB outflow: MgS in the extreme carbon star LL Peg

The broad 30 μm feature in carbon stars is commonly attributed to MgS dust particles. However, reproducing the 30 μ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 μ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 λ< 10 μ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 μ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 μm feature when it is part of a composite grain population that has optical properties representative of an ensemble of particle shapes.

preprint2012arXiv

On the physical structure of IRC+10216. Ground-based and Herschel observations of CO and CCH

The C-rich AGB star IRC+10216 undergoes strong mass loss, and quasi-periodic density enhancements in the circumstellar matter have been reported. CO is ubiquitous in the CSE, while CCH emission comes from a spatially confined shell. With the IRAM 30m telescope and Herschel/HIFI, we recently detected unexpectedly strong emission from the CCH N=4-3, 6-5, 7-6, 8-7, and 9-8 transitions, challenging the available chemical and physical models. We aim to constrain the physical properties of IRC+10216's CSE, including the effect of episodic mass loss on the observed emission. In particular, we aim to determine the excitation region and conditions of CCH and to reconcile these with interferometric maps of the N=1-0 transition. Via radiative-transfer modelling, we provide a physical description of the CSE, constrained by the SED and a sample of 20 high-resolution and 29 low-resolution CO lines. We further present detailed radiative-transfer analysis of CCH. Assuming a distance of 150pc, the SED is modelled with a stellar luminosity of 11300Lsun and a dust-mass-loss rate of 4.0\times10^{-8}Msun/yr. Based on the analysis of 20 high resolution CO observations, an average gas-mass-loss rate for the last 1000yrs of 1.5\times10^{-5}Msun/yr is derived. This gives a gas-to-dust-mass ratio of 375, typical for an AGB star. The gas kinetic temperature throughout the CSE is described by 3 powerlaws: it goes as r^{-0.58} for r<9R*, as r^{-0.40} for 9<=r<=65R*, and as r^{-1.20} for r>65R*. This model successfully describes all 49 CO lines. We show the effect of wind-density enhancements on the CCH-abundance profile, and the good agreement of the model with the CCH N=1-0 transition and with the lines observed with the 30m telescope and HIFI. We report on the importance of radiative pumping to the vibrationally excited levels of CCH and the significant effect this has on the excitation of all levels of the CCH-molecule.

preprint2011arXiv

First firm spectral classification of an early-B pre-main-sequence star: B275 in M17

The optical to near-infrared (300 - 2500 nm) spectrum of the candidate massive Young Stellar Object (YSO) B275, embedded in the star forming region M17, has been observed with X-shooter on the ESO Very Large Telescope. The spectrum includes both photospheric absorption lines and emission features (H and Ca triplet emission lines, 1st and 2nd overtone CO bandhead emission), as well as an infrared excess indicating the presence of a (flaring) circumstellar disk. The strongest emission lines are double-peaked with a peak separation ranging between 70 and 105 km s^-1 and provide information on the physical structure of the disk. The underlying photospheric spectrum is classified as B6-B7, significantly cooler than a previous estimate based on modeling of the spectral energy distribution. This discrepancy is solved by allowing for a larger stellar radius (i.e. a bloated star), and thus positioning the star above the main sequence. This constitutes the first firm spectral classification of an early-B pre-main-sequence (PMS) star. We discuss the position of B275 in the Hertzsprung-Russell diagram in terms of PMS evolution. Although the position is consistent with PMS tracks of heavily accreting protostars (Macc > 10^-5 Msun yr^-1), the fact that the photosphere of the object is detectable suggests that the current mass-accretion rate is not very high.

preprint2011arXiv

HIFI detection of HF in the carbon star envelope IRC +10216

We report the detection of emission in the J=1-0 rotational transition of hydrogen fluoride (HF), together with observations of the J=1-0 to J=3-2 rotational lines of H35Cl and H37Cl, towards the envelope of the carbon star IRC +10216. High-sensitivity, high-spectral resolution observations have been carried out with the HIFI instrument on board Herschel, allowing us to resolve the line profiles and providing insights into the spatial distribution of the emission. Our interpretation of the observations, with the use of radiative transfer calculations, indicates that both HF and HCl are formed in the inner regions of the envelope close to the AGB star. Thermochemical equilibrium calculations predict HF and HCl to be the major reservoirs of fluorine and chlorine in the atmospheres of AGB stars. The abundances relative to H2 derived for HF and HCl, 8e-9 and 1e-7 respectively, are substantially lower than those predicted by thermochemical equilibrium, indicating that F and Cl are likely affected by significant depletion onto dust grains, although some chlorine may be in the form of atomic Cl. The H35Cl/H37Cl abundance ratio is 3.3 +/- 0.3. The low abundance derived for HF in IRC +10216 makes it likely that the fluorine abundance is not enhanced over the solar value by nucleosynthesis in the AGB star, although this conclusion may not be robust because the HF abundance we derive is a lower limit to the elemental abundance of F. These observations suggest that both HF and HCl should be detectable through low J rotational transitions in other evolved stars.

preprint2011arXiv

Low abundance, strong features: Window-dressing crystalline forsterite in the disk wall of HD 100546

Forsterite is one of the crystalline dust species that is often observed in protoplanetary disks and solar system comets. Being absent in the interstellar medium, it must be produced during the disk lifetime. It can therefore serve as a tracer of dust processing and disk evolution, which can lead to a better understanding of the physical processes occurring in the disk, and possibly planet formation. However, the connection of these processes with the overall disk crystallinity remains unclear. We aim to characterize the forsterite abundance and spatial distribution in the disk of the Herbig Be star HD 100546, to investigate if a connection exists with the large disk gap. We use a 2D radiative transfer code, MCMax, to model the circumstellar dust around HD 100546. We use VISIR Q-band imaging to probe the outer disk geometry and mid-infrared features to model the spatial distribution of forsterite. The temperature-dependent shape of the 69 micron feature observed with Herschel PACS is used as a critical tool to constrain this distribution. We find a crystalline mass fraction of 40 - 60 %, located close to the disk wall between 13 and 20 AU, and possibly farther out at the disk surface. The forsterite is in thermal contact with the other dust species. We put an upper limit on the iron content of forsterite of 0.3 %. Optical depth effects play a key role in explaining the observed forsterite features, hiding warm forsterite from view at short wavelengths. The disk wall acts as a showcase: it displays a localized high abundance of forsterite, which gives rise to a high observed crystallinity, while the overall mass fraction of forsterite is a factor of ten lower.

preprint2011arXiv

Sequential star formation in IRAS 06084-0611 (GGD 12-15): From intermediate-mass to high-mass stars

Context. The formation and early evolution of high- and intermediate-mass stars towards the main sequence involves the interplay of stars in a clustered and highly complex environment. To obtain a full census of this interaction, the Formation and Early evolution of Massive Stars (FEMS) collaboration studies a well-selected sample of 10 high-mass star-forming regions. Aims. In this study we examine the stellar content of the high-mass star-forming region centered on IRAS 06084-0611 in the Monoceros R2 cloud. Methods. Using the near-infrared H- and K-band spectra from the VLT/SINFONI instrument on the ESO Very Large Telescope (VLT)and photometric near-infrared NTT/SOFI, 2MASS and Spitzer/IRAC data, we were able to determine the spectral types for the most luminous stars in the cluster. Results. Two very young and reddened massive stars have been detected by SINFONI: a massive Young Stellar Object (YSO) con- sistent with an early-B spectral type and a Herbig Be star. Furthermore, stars of spectral type G and K are detected while still in the Pre-Main Sequence (PMS) phase. We derive additional properties such as temperatures, extinctions, radii and masses. We present a Hertzsprung-Russell diagram and find most objects having intermediate masses between \sim1.5-2.5 M\odot. For these stars we derive a median cluster age of \sim4 Myr. Conclusions. Using Spitzer/IRAC data we confirm earlier studies that the younger class 0/I objects are centrally located while the class II objects are spread out over a larger area, with rough scale size radii of \sim0.5 pc and \sim1.25 pc respectively. Moreover, the presence of a massive YSO, an ultracompact H ii region and highly reddened objects in the center of the cluster suggest a much younger age of < 1 Myr. A possible scenario for this observation would be sequential star formation along the line of sight; from a cluster of intermediate-mass to high-mass stars.

preprint2011arXiv

The complex circumstellar environment of HD142527

The recent findings of gas giant planets around young A-type stars suggest that disks surrounding Herbig Ae/Be stars will develop planetary systems. An interesting case is HD142527, for which previous observations revealed a complex circumstellar environment and an unusually high ratio of infrared to stellar luminosity. Its properties differ considerably from other Herbig Ae/Be stars. This suggests that the disk surrounding HD142527 is in an uncommon evolutionary stage. We aim for a better understanding of the geometry and evolutionary status of the circumstellar material around the Herbig Ae/Be star HD142527. We map the composition and spatial distribution of the dust around HD142527. We analyze SEST and ATCA millimeter data, VISIR N and Q-band imaging and spectroscopy. We gather additional relevant data from the literature. We use the radiative transfer code MCMax to construct a model of the geometry and density structure of the circumstellar matter, which fits all of the observables satisfactorily. We find that the disk of HD142527 has three geometrically distinct components separated by a disk gap running from 30 to 130 AU. There is a geometrically flat inner disk running from 0.3 AU up to 30 AU; an optically thin halo-like component of dust in the inner disk regions; and a massive self-shadowed outer disk running from 130 AU up to 200 AU. We derived a total dust mass in small grains of 10^-3 M0 and a vertical height of the inner wall of the outer disk of h = 60 AU. ....

preprint2011arXiv

XMM-Newton observation of 4U 1820-30: Broad band spectrum and the contribution of the cold interstellar medium

We present the analysis of the bright X-ray binary 4U 1820-30, based mainly on XMM-Newton-RGS data, but using complementary data from XMM-Epic, Integral, and Chandra-HETG, to investigate different aspects of the source. The broad band continuum is well fitted by a classical combination of black body and Comptonized emission. The continuum shape and the high flux of the source (L/L_Edd\sim0.16) are consistent with a "high state" of the source. We do not find significant evidence of iron emission at energies >=6.4 keV. The soft X-ray spectrum contain a number of absorption features. Here we focus on the cold-mildly ionized gas. The neutral gas column density is N_H\sim1.63x10^21 cm^-2. The detailed study of the oxygen and iron edge reveals that those elements are depleted, defined here as the ratio between dust and the total ISM cold phase, by a factor 0.20\pm0.02 and 0.87\pm0.14, respectively. Using the available dust models, the best fit points to a major contribution of Mg-rich silicates, with metallic iron inclusion. Although we find that a large fraction of Fe is in dust form, the fit shows that Fe-rich silicates are disfavored. The measured Mg:Fe ratio is 2.0\pm0.3. Interestingly, this modeling may point to a well studied dust constituent (GEMS), sometimes proposed as a silicate constituent in our Galaxy. Oxygen and iron are found to be slightly over- and under-abundant, respectively (1.23 and 0.85 times the solar value) along this line of sight. We also report the detection of two absorption lines, tentatively identified as part of an outflow of mildly ionized gas (ξ\sim-0.5) at a velocity of \sim1200 km/s.

preprint2010arXiv

A high resolution line survey of IRC+10216 with Herschel. First results: Detection of warm silicon dicarbide SiC2

We present the first results of a high-spectral-resolution survey of the carbon-rich evolved star IRC+10216 that was carried out with the HIFI spectrometer onboard Herschel. This survey covers all HIFI bands, with a spectral range from 488 to 1901GHz. In this letter we focus on the band-1b spectrum, in a spectral range 554.5-636.5GHz, where we identified 130 spectral features with intensities above 0.03 K and a signal-to-noise ratio >5. Detected lines arise from HCN, SiO, SiS, CS, CO, metal-bearing species and, surprisingly, silicon dicarbide (SiC2). We identified 55 SiC2 transitions involving energy levels between 300 and 900 K. By analysing these rotational lines, we conclude that SiC2 is produced in the inner dust formation zone, with an abundance of ~2x10^-7 relative to molecular hydrogen. These SiC2 lines have been observed for the first time in space and have been used to derive an SiC2 rotational temperature of ~204 K and a source-averaged column density of ~6.4x10^15 cm^-2. Furthermore, the high quality of the HIFI data set was used to improve the spectroscopic rotational constants of SiC2.

preprint2010arXiv

An Inner Gaseous Disk around the Herbig Be Star MWC 147

We present high-spectral-resolution, optical spectra of the Herbig Be star MWC 147, in which we spectrally resolve several emission lines, including the [O I] lines at 6300 and 6363°. Their highly symmetric, double-peaked line profiles indicate that the emission originates in a rotating circumstellar disk. We deconvolve the Doppler-broadened [O I] emission lines to obtain a measure of emission as a function of distance from the central star. The resulting radial surface brightness profiles are in agreement with a disk structure consisting of a flat, inner, gaseous disk and a flared, outer, dust disk. The transition between these components at 2 to 3 AU corresponds to the estimated dust sublimation radius. The width of the double-peaked Mg II line at 4481°suggests that the inner disk extends to at least 0.10 AU, close to the corotation radius.

preprint2010arXiv

Determining the forsterite abundance of the dust around AGB stars

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

preprint2010arXiv

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 μm forsterite bands, that probe the most recent (11.3 μm) and older (33.6 μ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 μ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 μ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%.

preprint2010arXiv

Dust evolution in protoplanetary disks around Herbig Ae/Be stars - The Spitzer view

In this paper we present mid-infrared spectra of a comprehensive set of Herbig Ae/Be stars observed with the Spitzer Space Telescope. The signal-to-noise ratio of these spectra is very high, ranging between about a hundred and several hundreds. During the analysis of these data we tested the validity of standard protoplanetary dust models and studied grain growth and crystal formation. On the basis of the analyzed spectra, the major constituents of protoplanetary dust around Herbig Ae/Be stars are amorphous silicates with olivine and pyroxene stoichiometry, crystalline forsterite and enstatite and silica. No other solid state features, indicating other abundant dust species, are present in the Spitzer spectra. Deviations of the synthetic spectra from the observations are most likely related to grain shape effects and uncertainties in the iron content of the dust grains. Our analysis revealed that larger grains are more abundant in the disk atmosphere of flatter disks than in that of flared disks, indicating that grain growth and sedimentation decrease the disk flaring. We did not find, however, correlations between the value of crystallinity and any of the investigated system parameters. Our analysis shows that enstatite is more concentrated toward the warm inner disk than forsterite, in contrast to predictions of equilibrium condensation models. None of the three crystal formation mechanisms proposed so far can alone explain all our findings. It is very likely that all three play at least some role in the formation of crystalline silicates.

preprint2010arXiv

Dust, Ice and Gas in Time (DIGIT) Herschel program first results: A full PACS-SED scan of the gas line emission in protostar DK Cha

DK Cha is an intermediate-mass star in transition from an embedded configuration to a star plus disk stage. We aim to study the composition and energetics of the circumstellar material during this pivotal stage. Using the Range Scan mode of PACS on the Herschel Space Observatory, we obtained a spectrum of DK Cha from 55 to 210 micron as part of the DIGIT Key Program. Almost 50 molecular and atomic lines were detected, many more than the 7 lines detected in ISO-LWS. Nearly the entire ladder of CO from J=14-13 to 38-37 (E_u/k = 4080 K), water from levels as excited as E_u/k = 843 K, and OH lines up to E_u/k = 290 K were detected. The continuum emission in our PACS SED scan matches the flux expected from a model consisting of a star, a surrounding disk of 0.03 Solar mass, and an envelope of a similar mass, supporting the suggestion that the object is emerging from its main accretion stage. Molecular, atomic, and ionic emission lines in the far-infrared reveal the outflow's influence on the envelope. The inferred hot gas can be photon-heated, but some emission could be due to C-shocks in the walls of the outflow cavity.

preprint2010arXiv

First results of the Herschel Key Program 'Dust, Ice and Gas in Time': Dust and Gas Spectroscopy of HD 100546

We present far-infrared spectroscopic observations, taken with the Photodetector Array Camera and Spectrometer (PACS) on the Herschel Space Observatory, of the protoplanetary disk around the pre-main-sequence star HD 100546. These observations are the first within the DIGIT Herschel key program, which aims to follow the evolution of dust, ice, and gas from young stellar objects still embedded in their parental molecular cloud core, through the final pre-main-sequence phases when the circumstellar disks are dissipated. Our aim is to improve the constraints on temperature and chemical composition of the crystalline olivines in the disk of HD 100546 and to give an inventory of the gas lines present in its far-infrared spectrum. The 69 μ\m feature is analyzed in terms of position and shape to derive the dust temperature and composition. Furthermore, we detected 32 emission lines from five gaseous species and measured their line fluxes. The 69 μ\m emission comes either from dust grains with ~70 K at radii larger than 50 AU, as suggested by blackbody fitting, or it arises from ~200 K dust at ~13 AU, close to the midplane, as supported by radiative transfer models. We also conclude that the forsterite crystals have few defects and contain at most a few percent iron by mass. Forbidden line emission from [CII] at 157 μ\m and [OI] at 63 and 145 μ\m, most likely due to photodissociation by stellar photons, is detected. Furthermore, five H2O and several OH lines are detected. We also found high-J rotational transition lines of CO, with rotational temperatures of ~300 K for the transitions up to J=22-21 and T~800 K for higher transitions.

preprint2010arXiv

HD 95881: A gas rich to gas poor transition disk?

Context. Based on the far infrared excess the Herbig class of stars is divided into a group with flaring circumstellar disks (group I) and a group with flat circumstellar disks (group II). Dust sedimentation is generally proposed as an evolution mechanism to transform flaring disks into flat disks. Theory predicts that during this process the disks preserve their gas content, however observations of group II Herbig Ae stars demonstrate a lack of gas. Aims. We map the spatial distribution of the gas and dust around the group II Herbig Ae star HD 95881. Methods. We analyze optical photometry, Q-band imaging, infrared spectroscopy, and K and N-band interferometric spectroscopy. We use a Monte Carlo radiative transfer code to create a model for the density and temperature structure which quite accurately reproduces all the observables. Results. We derive a consistent picture in which the disk consists of a thick puffed up inner rim and an outer region which has a flaring gas surface and is relatively void of 'visible' dust grains. Conclusions. HD 95881 is in a transition phase from a gas rich flaring disk to a gas poor self-shadowed disk.

preprint2010arXiv

PACS spectroscopy of OH/IR stars

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.

preprint2010arXiv

Sequential Star Formation in RCW 34: A Spectroscopic Census of the Stellar Content of High-mass Star-forming Regions

We present VLT/SINFONI integral field spectroscopy of RCW 34 along with Spitzer/IRAC photometry of the surroundings. RCW 34 consists of three different regions. A large bubble has been detected on the IRAC images in which a cluster of intermediate- and low-mass class II objects is found. At the northern edge of this bubble, an HII region is located, ionized by 3 OB stars. Intermediate mass stars (2 - 3 Msun) are detected of G- and K- spectral type. These stars are still in the pre-main sequence (PMS) phase. North of the HII region, a photon-dominated region is present, marking the edge of a dense molecular cloud traced by H2 emission. Several class 0/I objects are associated with this cloud, indicating that star formation is still taking place. The distance to RCW 34 is revised to 2.5 +- 0.2 kpc and an age estimate of 2 - 1 Myrs is derived from the properties of the PMS stars inside the HII region. The most likely scenario for the formation of the three regions is that star formation propagates from South to North. First the bubble is formed, produced by intermediate- and low-mass stars only, after that, the HII region is formed from a dense core at the edge of the molecular cloud, resulting in the expansion as a champagne flow. More recently, star formation occurred in the rest of the molecular cloud. Two different formation scenarios are possible: (a) The bubble with the cluster of low- and intermediate mass stars triggered the formation of the O star at the edge of the molecular cloud which in turn induces the current star-formation in the molecular cloud. (b) An external triggering is responsible for the star-formation propagating from South to North. [abridged]

preprint2010arXiv

SPITZER-IRS spectral fitting of discs around binary post-AGB stars - Corrigendum

Recently, we have discovered an error in our Monte-Carlo spectral fitting routine, more specifically where the errors on the fluxes were rescaled to get a reduced chi2 of 1. The rescaled errors were too big, resulting in too wide a range of good fits in our 100 step Monte-Carlo routine. This problem affects Figs. 7-9 and Tables A.1, A.2 in Gielen et al. (2008), Table 3 in Gielen et al. (2009a), and Table 4 in Gielen et al. (2009b). We corrected for this error and present the new values and errors in the tables below. The new values and errors nearly all fall within the old error range. Our best chi2 values and overall former scientific results are not affected. With these new errors some possible new trends in the dust parameters might be observed. These will be discussed in an upcoming paper where we extend the sample presented in Gielen et al. (2008) with newly obtained SPITZER-IRS data.

preprint2010arXiv

Spitzer's view on aromatic and aliphatic hydrocarbon emission in Herbig Ae stars

The chemistry of astronomical hydrocarbons, responsible for the well-known infrared emission features detected in a wide variety of targets, remains enigmatic. Here we focus on the group of young intermediate-mass Herbig Ae stars. We have analyzed the aliphatic and polycyclic aromatic hydrocarbon (PAH) emission features in the infrared spectra of a sample of 53 Herbig Ae stars, obtained with the Infrared Spectrograph aboard the Spitzer Space Telescope. We confirm that the PAH-to-stellar luminosity ratio is higher in targets with a flared dust disk. However, a few sources with a flattened dust disk still show relatively strong PAH emission. Since PAH molecules trace the gas disk, this indicates that gas disks may still be flared, while the dust disk has settled due to grain growth. There are indications that the strength of the 11.3 um feature also depends on dust disk structure, with flattened disks being less bright in this feature. We confirm that the CC bond features at 6.2 and 7.8 um shift to redder wavelengths with decreasing stellar effective temperature. Moreover, we show that this redshift is accompanied by a relative increase of aliphatic CH emission and a decrease of the aromatic 8.6 um CH feature strength. Cool stars in our sample are surrounded by hydrocarbons with a high aliphatic/aromatic CH ratio and a low aromatic CH/CC ratio, and vice versa for the hot stars. We conclude that, while the overall hydrocarbon emission strength depends on the dust disk's geometry, the relative differences seen in the IR emission features in disks around Herbig Ae stars are mainly due to chemical differences of the hydrocarbon molecules induced by the stellar UV field. Strong UV flux reduces the aliphatic component and emphasizes the spectral signature of the aromatic molecules in the IR spectra.

preprint2010arXiv

The 9.7 and 18 um silicate absorption profiles towards diffuse and molecular cloud lines-of-sight

Studying the composition of dust in the interstellar medium (ISM) is crucial in understanding the cycle of dust in our galaxy. The mid-infrared spectral signature of amorphous silicates, the most abundant dust species in the ISM, is studied in different lines-of-sight through the Galactic plane, thus probing different conditions in the ISM. We have analysed 10 spectra from the Spitzer archive, of which 6 lines-of-sight probe diffuse interstellar medium material and 4 probe molecular cloud material. The 9.7 um silicate absorption features in 7 of these spectra were studied in terms of their shape and strength. In addition, the shape of the 18 um silicate absorption features in 4 of the diffuse sightline spectra were analysed. The 9.7 um silicate absorption bands in the diffuse sightlines show a strikingly similar band shape. This is also the case for all but one of the 18 um silicate absorption bands observed in diffuse lines-of-sight. The 9.7 um bands in the 4 molecular sightlines show small variations in shape. These modest variations in the band shape are inconsistent with the interpretation of the large variations in τ_9.7/E(J-K) between diffuse and molecular sightlines in terms of silicate grain growth. Instead, we suggest that the large changes in τ_9.7 / E(J-K) must be due to changes in E(J-K).

preprint2010arXiv

The lunar phases of dust grains orbiting Fomalhaut

Optical images of the nearby star Fomalhaut show a ring of dust orbiting the central star. This dust is in many respects expected to be similar to the zodiacal dust in the solar system. The ring displays a clear brightness asymmetry, attributed to asymmetric scattering of the central starlight by the circumstellar dust grains. Recent measurements show that the bright side of the Fomalhaut ring is oriented away from us. This implies that the grains in this system scatter most of the light in the backward direction, in sharp contrast to the forward-scattering nature of the grains in the solar system. In this letter, we show that grains considerably larger than those dominating the solar system zodiacal dust cloud provide a natural explanation for the apparent backward scattering behavior. In fact, we see the phases of the dust grains in the same way as we can observe the phases of the Moon and other large solar system bodies. We outline how the theory of the scattering behavior of planetesimals can be used to explain the Fomalhaut dust properties. This indicates that the Fomalhaut dust ring is dominated by very large grains. The material orbiting Fomalhaut, which is at the transition between dust and planetesimals, can, with respect to their optical behavior, best be described as micro-asteroids.

preprint2010arXiv

The widespread occurence of water vapor in the circumstellar envelopes of carbon-rich AGB stars: first results from a survey with Herschel/HIFI

We report the preliminary results of a survey for water vapor in a sample of eight C stars with large mid-IR continuum fluxes: V384 Per, CIT 6, V Hya, Y CVn, IRAS 15194-5115, V Cyg, S Cep, and IRC+40540. This survey, performed using the HIFI instrument on board the Herschel Space Observatory, entailed observations of the lowest transitions of both ortho- and para-water: the 556.936 GHz 1(10)-1(01) and 1113.343 GHz 1(11)-0(00) transitions, respectively. Water vapor was unequivocally detected in all eight of the target stars. Prior to this survey, IRC+10216 was the only carbon-rich AGB star from which thermal water emissions had been discovered, in that case with the use of the Submillimeter Wave Astronomy Satellite (SWAS). Our results indicate that IRC+10216 is not unusual, except insofar as its proximity to Earth leads to a large line flux that was detectable with SWAS. The water spectral line widths are typically similar to those of CO rotational lines, arguing against the vaporization of a Kuiper belt analog (Ford & Neufeld 2001) being the general explanation for water vapor in carbon-rich AGB stars. There is no apparent correlation between the ratio of the integrated water line fluxes to the 6.3 micron continuum flux - a ratio which measures the water outflow rate - and the total mass-loss rate for the stars in our sample.

preprint2009arXiv

A binary engine fuelling HD87643' s complex circumstellar environment, using AMBER/VLTI

Context. The star HD 87643, exhibiting the "B[e] phenomenon", has one of the most extreme infrared excesses for this object class. It harbours a large amount of both hot and cold dust, and is surrounded by an extended reflection nebula. Aims. One of our major goals was to investigate the presence of a companion in HD87643. In addition, the presence of close dusty material was tested through a combination of multi-wavelength high spatial 5Aresolution observations. Methods. We observed HD 87643 with high spatial resolution techniques, using the near-IR AMBER/VLTI interferometer with baselines ranging from 60 m to 130 m and the mid-IR MIDI/VLTI interferometer with baselines ranging from 25 m to 65 m. These observations are complemented by NACO/VLT adaptive-optics-corrected images in the K and L-bands, ESO-2.2m optical Wide-Field Imager large-scale images in the B, V and R-bands, Results. We report the direct detection of a companion to HD 87643 by means of image synthesis using the AMBER/VLTI instrument. The presence of the companion is confirmed by the MIDI and NACO data, although with a lower confidence. The companion is separated by ~ 34 mas with a roughly north-south orientation. The period must be large (several tens of years) and hence the orbital parameters are not determined yet. Binarity with high eccentricity might be the key to interpreting the extreme characteristics of this system, namely a dusty circumstellar envelope around the primary, a compact dust nebulosity around the binary system and a complex extended nebula witnessing past violent ejections.

preprint2009arXiv

A dam around the Water Fountain Nebula? The dust shell of IRAS16342-3814 spatially resolved with VISIR/VLT

Either by collimating a fast stellar wind or by driving a jet via accretion in the central system, dusty torii or stable disks may be crucial ingredients for the shaping of PNe. We study the dust distribution in the very young Proto-Planetary Nebule (PPN) IRAS16342-3814, also known as the Water Fountain Nebula, which is known to show strong bipolar characteristics in the shape of two reflection lobes, and high-velocity collimated molecular outlfows. We use the new Mid-IR (MIR) instrument VISIR on the Very Large Telescope (VLT) both in imaging and spectroscopy mode at wavelengths from 8 to 13 micron. We present the first spatially resolved MIR observations of a dusty evolved star obtained with VISIR and find that the improved spatial resolution contradicts previous claims of an elliptical brightness distribution at the heart of IRAS16342: we find the waist region to be dark even in the MIR. We show that the filling angle of the obscuring dust lane, which is made mostly of amorphous silicates, is very large, possibly even close to a spherically symmetric superwind as seen in OH/IR stars. We conclude that, in contrast to the multitude of recent dusty-disk detections in Post-AGB stars and PNe, IRAS16342 does not show this extreme equatorial density enhancement, at least not on the scale of the dusty environment which lends the object its IR appearance. Rather, it appears that the observed precessing jets are shaping the bipolar nature in the remains of a spherically symmetric AGB superwind.

preprint2009arXiv

Analysis of the infrared spectra of the peculiar post-AGB stars EPLyr and HD52961

Aim: We aim to study in detail the peculiar mineralogy and structure of the circumstellar environment of two binary post-AGB stars, EPLyr and HD52961. Both stars were selected from a larger sample of evolved disc sources observed with Spitzer and show unique solid-state and gas features in their infrared spectra. Moreover, they show a very small infrared excess in comparison with the other sample stars. Methods: The different dust and gas species are identified on the basis of high-resolution Spitzer-IRS spectra. We fit the full spectrum to constrain grain sizes and temperature distributions in the discs. This, combined with our broad-band spectral energy distribution and interferometric measurements, allows us to study the physical structure of the disc, using a self-consistent 2D radiative-transfer disc model. Results: We find that both stars have strong emission features due to CO_2 gas, dominated by ^{12}C^{16}O_2, but with clear ^{13}C^{16}O_2 and even ^{16}O^{12}C^{18}O isotopic signatures. Crystalline silicates are apparent in both sources but proved very hard to model. EP Lyr also shows evidence of mixed chemistry, with emission features of the rare class-C PAHs. Whether these PAHs reside in the oxygen-rich disc or in a carbon-rich outflow is still unclear. With the strongly processed silicates, the mixed chemistry and the low ^{12}C/^{13}C ratio, EP Lyr resembles some silicate J-type stars, although the depleted photosphere makes nucleosynthetic signatures difficult to probe. We find that the disc environment of both sources is, to a first approximation, well modelled with a passive disc, but additional physics such as grain settling, radial dust distributions, and an outflow component must be included to explain the details of the observed spectral energy distributions in both stars.

preprint2009arXiv

Evidence for CO depletion in the inner regions of gas-rich protoplanetary disks

We investigate the physical properties and spatial distribution of Carbon Monoxide (CO) gas in the disks around the Herbig Ae/Be stars HD 97048 and HD 100546. Using high-spectral-resolution 4.588-4.715 $μ$m spectra containing fundamental CO emission taken with CRIRES on the VLT, we probe the circumstellar gas and model the kinematics of the emission lines. By using spectro-astrometry on the spatially resolved targets, we constrain the physical size of the emitting regions in the disks. We resolve, spectrally and spatially, the emission of the $^{13}$CO v(1-0) vibrational band and the $^{12}$CO $v=1-0, v=2-1, v=3-2$ and $v=4-3$ vibrational bands in both targets, as well as the $^{12}$CO $v=5-4$ band in HD 100546. Modeling of the CO emission with a homogeneous disk in Keplerian motion, yields a best fit with an inner and outer radius of the CO emitting region of 11 and $\geq$ 100 AU for HD 97048. HD 100546 is not fit well with our model, but we derive a lower limit on the inner radius of 8 AU. The fact that gaseous [OI] emission was previously detected in both targets at significantly smaller radii suggests that CO may be effectively destroyed at small radii in the surface layers of these disks

preprint2009arXiv

Spitzer/IRAC view of Sh 2-284: Searching for evidence of triggered star formation in an isolated region in the outer Milky Way

Using Spitzer/IRAC observations of a region to be observed by the CoRoT satellite, we have unraveled a new complex star-forming region at low metallicity in the outer Galaxy. We perform a study of S284 in order to outline the chain of events in this star-forming region. We used four-band Spitzer/IRAC photometry as well as Halpha imaging obtained with INT/WFC. Combining these data with the optical photometry obtained in the frame of CoRoTs preparation and the 2MASS catalog we analysed the properties and distribution of young stellar objects (YSOs) associated with point-like sources. We also studied the SEDs of regions of extended emission, complementing our dataset with IRAS and MSX data. We find that S284 is unique in several ways: it is very isolated at the end of a spiral arm and both the diffuse dust and ionized emission are remarkably symmetric. We have partially resolved the central clusters of the three bubbles present in this region. Despite the different scales present in its multiple-bubble morphology, our study points to a very narrow spread of ages among the powering high-mass clusters. In contrast, the particular sawtooth structure of the extended emission at the rim of each ionized bubble harbours either small lower-mass clusters with a younger stellar population or individual young reddened protostars. In particular, triggered star formation is considered to be at work in these regions.

preprint2005arXiv

ESO-VLT and Spitzer spectroscopy of IRAS05328-6827: a massive young stellar object in the Large Magellanic Cloud

We present the first thermal-infrared spectra of an extra-galactic Young Stellar Object (YSO), IRAS05328-6827 in the HII region LHA 120-N148 in the Large Magellanic Cloud. The observed and modelled spectral energy distribution reveals a massive YSO, M=20 Msun, which is heavily-embedded and probably still accreting. The reduced dust content as a consequence of the lower metallicity of the LMC allows a unique view into this object, and together with a high C/O ratio may be responsible for the observed low abundance of water ice and relatively high abundances of methanol and CO2 ices.

preprint2003arXiv

Identification of the ionizing source of NGC 2024

We propose the late-O, early-B star IRS2b as the ionizing source of the Flame Nebula (NGC 2024). It has been clear that such a hot, massive star must be present in this heavily obscured region, and now it has been identified. New near-infrared photometry shows that IRS2b is the most luminous and hottest star in the young star cluster embedded in the center of NGC 2024. The near-infrared observations (5' x 5') cover ~90 % of the HII region detected in radio continuum radiation, making the probability very low that the ionizing star is not present in the field. A K-band spectrum of IRS2b obtained with ISAAC on the Very Large Telescope indicates that the spectral type of IRS2b is in the range O8V - B2V. Additional arguments indicate that its spectral type is likely closer to O8 than to B2. The corresponding amount of ionizing radiation is consistent with published radio continuum and recombination line observations.

preprint2001arXiv

The composition of circumstellar gas and dust in 51 Oph

We analyze ISO archive data of the nearby bright emission-line star 51 Oph, previously classified as a proto-planetary system similar to beta Pic. The infrared spectrum reveals the presence of gas-phase emission bands of hot (approx 850 K) CO, CO_2, H_2O and NO. In addition to this, partially crystalline silicate dust is present. The solid-state bands and the energy distribution are indicative of dust that has formed recently, rather than of debris dust. The presence of hot molecular gas and the composition of the circumstellar dust are highly unusual for a young star and are reminiscent of what is found around evolved (AGB) stars, although we exclude the possibility of 51 Oph belonging to this group. We suggest several explanations for the nature of 51 Oph, including a recent episode of mass loss from a Be star, and the recent destruction of a planet-sized body around a young star.

preprint2000arXiv

ISO spectroscopy of circumstellar dust in 14 Herbig Ae/Be systems: towards an understanding of dust processing

We present Infrared Space Observatory (ISO) spectra of fourteen isolated Herbig Ae/Be (HAEBE) stars, to study the characteristics of their circumstellar dust. These spectra show large star-to-star differences, in the emission features of both carbon-rich and oxygen-rich dust grains. The IR spectra were combined with photometric data ranging from the UV through the optical into the sub-mm region. We defined two key groups, based upon the spectral shape of the infrared region. The derived results can be summarized as follows: (1) the continuum of the IR to sub-mm region of all stars can be reconstructed by the sum of a power-law and a cool component, which can be represented by a black body. Possible locations for these components are an optically thick, geometrically thin disc (power-law component) and an optically thin flared region (black body); (2) all stars have a substantial amount of cold dust around them, independent of the amount of mid-IR excess they show; (3) also the near-IR excess is unrelated to the mid-IR excess, indicating different composition/location of the emitting material; (4) remarkably, some sources lack the silicate bands; (5) apart from amorphous silicates, we find evidence for crystalline silicates in several stars, some of which are new detections; (6) PAH bands are present in at least 50% of our sample, and their appearance is slightly different from PAHs in the ISM; (7) PAH bands are, with one exception, not present in sources which only show a power-law continuum in the IR; their presence is unrelated to the presence of the silicate bands; (8) the dust in HAEBE stars shows strong evidence for coagulation; this dust processing is unrelated to any of the central star properties (such as age, spectral type and activity).

preprint1999arXiv

The central stars of the planetary nebulae NGC 7027 and NGC 6543

Infrared spectra of NGC 7027 and NGC 6543 ranging from 2.4 to 45 micron were obtained with the Short Wavelength Spectrometer on board the Infrared Space Observatory. A first analysis of these spectra, with the aid of photo-ionization models, is presented. We report the first detection of the Ar VI 4.53 micron and Ne VI 7.65 micron lines in the spectrum of NGC 7027. When compared with older observations it is clear that the Ar VI line and possibly also other lines have increased in strength since 1981. We argue that a likely explanation for this variability is a change in the spectral energy distribution of the central star, possibly an increase in effective temperature. However, this result needs to be confirmed by further observations. We also report a non-detection of the O IV 25.9 micron line and the first detection of the Na III 7.32 micron line in the spectrum of NGC 6543. The non-detection is not expected based on a blackbody approximation for the spectrum of the central star. The ionization threshold for O(3+) is just beyond the He II limit, and the absence of this line shows that the stellar flux drops at least by a factor 350 at the He II limit. Modeling the O IV line may prove to be a valuable test for atmosphere models.

preprint1999arXiv

The spectral evolution of post-AGB stars

A parameter study of the spectral evolution of a typical post-AGB star, with particular emphasis on the evolution of the infrared colours, is presented. The models are based on the latest evolutionary tracks for hydrogen burning post-AGB stars. For such tracks the evolutionary rate is very dependent on the assumed mass loss rate as a function of time. We investigate this effect by modifying the mass loss prescription. The newly calculated evolutionary rates and density distributions are used to model the spectral evolution of a post-AGB star with the photo-ionization code Cloudy, including dust in the radiative transfer. Different assumptions for the dust properties and dust formation are considered. It is shown that by varying these parameters in a reasonable way, entirely different paths are followed in the IRAS colour-colour diagram. First of all, the effects of the evolution of the central star on the expanding dust shell can not be neglected. Also the dust properties and the definition of the end of the AGB phase have an important effect. The model tracks show that objects occupying the same location in the IRAS colour-colour diagram can have a different evolutionary past, and therefore the position in the IRAS colour-colour diagram alone can not a priori give a unique determination of the evolutionary status of an object. An alternative colour-colour diagram, the K-[12] vs. [12]-[25] diagram, is presented. The tracks in this diagram seem less affected by particulars of the grain emission. This diagram may be a valuable additional tool for studying post-AGB evolution.

preprint1997arXiv

Obscured Asymptotic Giant Branch stars in the Magellanic Clouds III. New IRAS counterparts

We have searched for near-infrared stellar counterparts of IRAS point sources in the Large Magellanic Cloud (LMC), in J and K-bands. This resulted in the detection of 21 counterparts, of which 19 are new discoveries. Using colour-magnitude and colour-colour diagrams, we identify 13 Asymptotic Giant Branch (AGB) stars with thick circumstellar dust envelopes, 7 possible early post-AGB stars or stars recovering from a thermal pulse, and 1 red supergiant or foreground star. For 10 of the IRAS targets we do not succeed in detecting and/or identifying a near-infrared counterpart. We serendipitously detect 14 other red sources, of which 2 are known Long Period Variables, and a few galaxies. The near-infrared and optical colours of the galaxies may indicate considerable interstellar extinction through the LMC, as much as A_V about 2-4 mag. The relative number of AGB carbon stars over oxygen stars is shown to decrease as the luminosity increases. Yet amongst the faintest mass-losing AGB stars oxygen-rich stars still exist, which puts constraints on current convection theories that predict the occurrence of third dredge-up and Hot Bottom Burning. We investigate the nature of some LMC stars that have infrared properties very similar to suspected Galactic post-AGB stars.

preprint1996arXiv

Circumstellar C2, CN, and CH+ in the optical spectra of post-AGB stars

We present optical high-resolution spectra of a sample of sixteen post-AGB stars and IRC +10216. Of the post-AGB stars, ten show C2 Phillips and Swan and CN Red System absorption, one CH+ emission, one CH+ absorption, and four without any molecules. We find typically Trot=43-399, 155-202, and 18-50 K, log N = 14.90-15.57, 14.35, and 15.03-16.47 cm-2 for C2, CH+, and CN respectively, and 0.6<N(CN)/N(C2)<11.2. We did not detect isotopic lines, which places a lower limit on the isotope ratio of 12C/13C>20. The presence of C2 and CN absorption is correlated with cold dust (Tdust<300K) and the presence of CH+ with hot dust (Tdust>300K). All objects with the unidentified 21mum emission feature exhibit C2 and CN absorption, but not all objects with C2 and CN detections exhibit a 21mum feature. The derived expansion velocity, ranging from 5 to 44 km/s, is the same as that derived from CO millimeter line emission. This unambiguously proves that these lines are of circumstellar origin and are formed in the AGB ejecta (circumstellar shell expelled during the preceding AGB phase). Furthermore there seems to be a relation between the C2 molecular column density and the expansion velocity, which is attributed to the fact that a higher carbon abundance of the dust leads to a more efficient acceleration of the AGB wind. Using simple assumptions for the location of the molecular lines and molecular abundances, mass-loss rates have been derived from the molecular absorption lines and are comparable to those obtained from CO emission lines and the infrared excess.

preprint1995arXiv

A study on the nature of the peculiar supergiant HD101584

We present a study of low- and high-resolution ultraviolet, high-resolution optical CAT/CES spectra and ultraviolet, optical and infrared photometry of the peculiar supergiant HD101584. From the photometry we learn that the ultraviolet and optical energy distribution cannot be fitted in a consistent way and we need a model in which the UV and optical energy distribution are formed by different gas. The Geneva photometry is best fitted to a B9II Kurucz model, Teff=12000+-1000K and log g=3.0 +-1.0, with an extinction of E(B-V)=0.49+-0.05. The observed spectral features in the spectrum of HD101584 are classified in eight different categories based on the velocity, shape of profile and the identification. The high-excitation HeI(chi=20.87eV), NII(chi=18.40eV), CII (chi=14.39eV) and NI (chi=10.29eV) optical absorption lines are formed in the photosphere of a late B-star (e.g. B8-9I-II). These absorption lines show radial velocity variations which are attributed to binary motion, with the secondary being a white dwarf or a low-mass main sequence object. The low-excitation P-Cygni lines in the optical and UV are formed in the wind. The number density of absorption lines in the UV is so large that the wind spectrum acts as an iron curtain in front of the B-star. The terminal velocity of the wind of $v_{\infty}=100\pm30$ km~s$^{-1}~$ is consistent with the star being a low-mass post-AGB star and the low effective gravity is attributed to the presence of a nearby, unseen, secondary. We estimate a mass-loss rate of M$\approx 10^{-8}$~M$_{\odot}$~yr$^{-1}$. Narrow absorption and emission lines are observed which are formed in a circumsystem disk with a typical radius of $10^{2} R_{\ast}$.

preprint1995arXiv

Obscured AGB stars in the Magellanic Clouds II. Near-infrared and mid-infrared counterparts

We have carried out an infrared search for obscured AGB stars in the Magellanic Clouds. The survey uncovered a number of obscured AGB stars as well as some supergiants with infrared excess. We present photometry of the sources and discuss the colour diagrams and bolometric luminosities. Most of the AGB stars are luminous, often close to the classical limit of $M_{\rm bol}=-7.1$. To determine whether the stars are oxygen-rich or carbon-rich, we have acquired narrow-band mid-infrared photometry with the ESO TIMMI camera for several sources. All but one are found to show the silicate feature and therefore to have oxygen-rich dust: the colours of the remaining source are consistent with either an oxygen-rich or a carbon-rich nature. A method to distinguish carbon and oxygen stars based on H$-$K versus K$-$[12] colours is presented. We discuss several methods of calculating the mass-loss rate: for the AGB stars the mass-loss rates vary between approximately 5 times 10**-4 and 5 times 10**-6 solar masses per year, depending on assumed dust-to-gas mass ratio. We present a new way to calculate mass-loss rates from the OH-maser emission. We find no evidence for a correlation of the mass-loss rates with luminosity in these obscured stars. Neither do the mass-loss rates for the LMC and SMC stars differ in any clear systematic way from each other. Expansion velocities appear to be slightly lower in the LMC than in the Galaxy. Period determinations are discussed for two sources: the periods are comparable to those of the longer-period galactic OH/IR stars. All of the luminous stars for which periods are available, have significantly higher luminosities than predicted from the period--luminosity relations.

L. B. F. M. Waters

What is connected

Connect this record

See the researcher in context

Building this map preview

79 published item(s)

Multi-band infrared imaging reveals dusty spiral arcs around the binary B[e] star 3 Puppis

Forming planets around stars with non-solar elemental composition

MATISSE, the VLTI mid-infrared imaging spectro-interferometer

The disk of FU Orionis viewed with MATISSE/VLTI: first interferometric observations in $L$ and $M$ bands

The Infrared Database of Extragalactic Observables from Spitzer. -- II. The Database & The Diagnostic Power of Crystalline Silicate Features in Galaxy Spectra

ATOMIUM: ALMA tracing the origins of molecules in dust forming oxygen rich M-type stars: Motivation, sample, calibration, and initial results

VLTI-MATISSE L- and N-band aperture-synthesis imaging of the unclassified B[e] star FS Canis Majoris

Inner dusty envelope of the AGB stars W~Hydrae, SW\,Virginis, and R~Crateris using SPHERE/ZIMPOL

The asymmetric inner disk of the Herbig Ae star HD 163296 in the eyes of VLTI/MATISSE: evidence for a vortex?

A proposed new diagnostic for Herbig disc geometry: FWHM versus J of CO ro-vibrational lines

Shadows cast on the transition disk of HD 135344B. Multiwavelength VLT/SPHERE polarimetric differential imaging

Study of the inner dust envelope and stellar photosphere of the AGB star R Doradus using SPHERE/ZIMPOL

The abundance and thermal history of water ice in the disk surrounding HD142527 from the DIGIT Herschel Key Program

Variability in the CO ro-vibrational lines from HD163296

Consistent dust and gas models for protoplanetary disks. I. Disk shape, dust settling, opacities, and PAHs

Dusty wind of W Hya. Multi-wavelength modelling of the present-day and recent mass-loss

Herschel observations of extreme OH/IR stars - the isotopic ratios of oxygen as a sign-post for the stellar mass

Micron-sized forsterite grains in the pre-planetary nebula of IRAS 17150-3224 - Searching for clues on the mysterious evolution of massive AGB stars

The Mid-Infrared Instrument for the James Webb Space Telescope, VI: The Medium Resolution Spectrometer

The structure of disks around intermediate-mass young stars from mid-infrared interferometry. Evidence for a population of group II disks with gaps

Dust composition and mass-loss return from the luminous blue variable R71 in the LMC

Herschel/PACS observations of the 69 $μm$ band of crystalline olivine around evolved stars

Location and sizes of forsterite grains in protoplanetary disks: interpretation from the Herschel DIGIT programme

Polycyclic aromatic hydrocarbon ionization as a tracer of gas flows through protoplanetary disk gaps

Surprising detection of an equatorial dust lane on the AGB star IRC+10216

The structure of disks around Herbig Ae/Be stars as traced by CO ro-vibrational emission

The wind of W Hya as seen by Herschel - I. The CO envelope

The wind of W Hya as seen by Herschel. II. The molecular envelope of W Hya

CO ro-vibrational lines in HD100546: A search for disc asymmetries and the role of fluorescence

HIFISTARS Herschel/HIFI observations of VY Canis Majoris. Molecular-line inventory of the envelope around the largest known star

Identifying gaps in flaring Herbig Ae/Be disks using spatially resolved mid-infrared imaging. Are all group I disks transitional?

OH/IR stars and their superwinds as observed by the Herschel Space Observatory

RCW36: characterizing the outcome of massive star formation

The 69 micron forsterite band in spectra of protoplanetary disks - Results from the Herschel DIGIT programme

The color dependent morphology of the post-AGB star HD161796

The problematically short superwind of OH/IR stars - Probing the outflow with the 69 μm spectral band of forsterite

The VLTI/MIDI survey of massive young stellar objects - Sounding the inner regions around intermediate- and high-mass young stars using mid-infrared interferometry

Comet-like mineralogy of olivine crystals in an extrasolar proto-Kuiper belt

Direct imaging of a massive dust cloud around R Coronae Borealis

Herschel/HIFI observations of red supergiants and yellow hypergiants: I. Molecular inventory

Mid-IR imaging of the transitional disk of HD169142: Measuring the size of the gap

Observational evidence for composite grains in an AGB outflow: MgS in the extreme carbon star LL Peg

On the physical structure of IRC+10216. Ground-based and Herschel observations of CO and CCH

First firm spectral classification of an early-B pre-main-sequence star: B275 in M17

HIFI detection of HF in the carbon star envelope IRC +10216

Low abundance, strong features: Window-dressing crystalline forsterite in the disk wall of HD 100546

Sequential star formation in IRAS 06084-0611 (GGD 12-15): From intermediate-mass to high-mass stars

The complex circumstellar environment of HD142527

XMM-Newton observation of 4U 1820-30: Broad band spectrum and the contribution of the cold interstellar medium

A high resolution line survey of IRC+10216 with Herschel. First results: Detection of warm silicon dicarbide SiC2

An Inner Gaseous Disk around the Herbig Be Star MWC 147

Determining the forsterite abundance of the dust around AGB stars

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

Dust evolution in protoplanetary disks around Herbig Ae/Be stars - The Spitzer view

Dust, Ice and Gas in Time (DIGIT) Herschel program first results: A full PACS-SED scan of the gas line emission in protostar DK Cha

First results of the Herschel Key Program 'Dust, Ice and Gas in Time': Dust and Gas Spectroscopy of HD 100546

HD 95881: A gas rich to gas poor transition disk?

PACS spectroscopy of OH/IR stars

Sequential Star Formation in RCW 34: A Spectroscopic Census of the Stellar Content of High-mass Star-forming Regions

SPITZER-IRS spectral fitting of discs around binary post-AGB stars - Corrigendum

Spitzer's view on aromatic and aliphatic hydrocarbon emission in Herbig Ae stars

The 9.7 and 18 um silicate absorption profiles towards diffuse and molecular cloud lines-of-sight

The lunar phases of dust grains orbiting Fomalhaut

The widespread occurence of water vapor in the circumstellar envelopes of carbon-rich AGB stars: first results from a survey with Herschel/HIFI

A binary engine fuelling HD87643' s complex circumstellar environment, using AMBER/VLTI

A dam around the Water Fountain Nebula? The dust shell of IRAS16342-3814 spatially resolved with VISIR/VLT

Analysis of the infrared spectra of the peculiar post-AGB stars EPLyr and HD52961

Evidence for CO depletion in the inner regions of gas-rich protoplanetary disks

Spitzer/IRAC view of Sh 2-284: Searching for evidence of triggered star formation in an isolated region in the outer Milky Way

ESO-VLT and Spitzer spectroscopy of IRAS05328-6827: a massive young stellar object in the Large Magellanic Cloud

Identification of the ionizing source of NGC 2024

The composition of circumstellar gas and dust in 51 Oph

ISO spectroscopy of circumstellar dust in 14 Herbig Ae/Be systems: towards an understanding of dust processing

The central stars of the planetary nebulae NGC 7027 and NGC 6543