Search

During their asymptotic giant branch, evolution low-mass stars lose a significant fraction of their mass through an intense wind, enriching the interstellar medium with products of nucleosynthesis. We observed the nearby oxygen-rich asymptotic giant branch star IK Tau using the high-resolution HIFI spectrometer onboard Herschel. We report on the first detection of H_2^{16}O and the rarer isotopologues H_2^{17}O and H_2^{18}O in both the ortho and para states. We deduce a total water content (relative to molecular hydrogen) of 6.6x10^{-5}, and an ortho-to-para ratio of 3:1. These results are consistent with the formation of H_2O in thermodynamical chemical equilibrium at photospheric temperatures, and does not require pulsationally induced non-equilibrium chemistry, vaporization of icy bodies or grain surface reactions. High-excitation lines of 12CO, 13CO, 28SiO, 29SiO, 30SiO, HCN, and SO have also been detected. From the observed line widths, the acceleration region in the inner wind zone can be characterized, and we show that the wind acceleration is slower than hitherto anticipated., Comment: 5 pages, 3 figures, Appendix of 3 pages and 3 figures, Astronomy & Astrophysics, in press

CONTEXT: Water vapour maser emission from evolved oxygen-rich stars remains poorly understood. Additional observations, including polarisation studies and simultaneous observation of different maser transitions may ultimately lead to greater insight. AIMS: We have aimed to elucidate the nature and structure of the VY CMa water vapour masers in part by observationally testing a theoretical prediction of the relative strengths of the 620.701 GHz and the 22.235 GHz maser components of ortho water vapour. METHODS: In its high-resolution mode (HRS) the Herschel Heterodyne Instrument for the Infrared (HIFI) offers a frequency resolution of 0.125 MHz, corresponding to a line-of-sight velocity of 0.06 km/s, which we employed to obtain the strength and linear polarisation of maser spikes in the spectrum of VY CMa at 620.701 GHz. Simultaneous ground based observations of the 22.235 GHz maser with the Max-Planck-Institut f\"ur Radioastronomie 100-meter telescope at Effelsberg, provided a ratio of 620.701 GHz to 22.235 GHz emission. RESULTS:We report the first astronomical detection to date of water vapour maser emission at 620.701 GHz. In VY CMa both the 620.701 and the 22.235 GHz polarisation are weak. At 620.701 GHz the maser peaks are superposed on what appears to be a broad emission component, jointly ejected asymmetrically from the star. We observed the 620.701 GHz emission at two epochs 21 days apart, both to measure the potential direction of linearly polarised maser components and to obtain a measure of the longevity of these components. Although we do not detect significant polarisation levels in the core of the line, they rise up to approximately 6% in its wings.

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., Comment: 9 pages, 7 figures, 2 tables, accepted by A&A

Aims: Little information exists on the circumstellar molecular abundance stratifications of many molecules. The aim is to study the circumstellar chemical abundance pattern of 11 molecules and isotopologs ($^{12}$CO, $^{13}$CO, SiS, $^{28}$SiO, $^{29}$SiO, $^{30}$SiO, HCN, CN, CS, SO, SO$_2$) in the oxygen-rich evolved star IK~Tau. Methods: We have performed an in-depth analysis of a large number of molecular emission lines excited in the circumstellar envelope around IK~Tau. The analysis is done based on a non-local thermodynamic equilibrium (non-LTE) radiative transfer analysis, which calculates the temperature and velocity structure in a self-consistent way. The chemical abundance pattern is coupled to theoretical outer wind model predictions including photodestruction and cosmic ray ionization. Not only the integrated line intensities, but also the line shapes, are used as diagnostic tool to study the envelope structure. Results: The deduced wind acceleration is much slower than predicted from classical theories. SiO and SiS are depleted in the envelope, possibly due to the adsorption onto dust grains. For HCN and CS a clear difference with respect to inner wind non-equilibrium predictions is found, either indicating uncertainties in the inner wind theoretical modeling or the possibility that HCN and CS (or the radical CN) participate in the dust formation. The low signal-to-noise profiles of SO and CN prohibit an accurate abundance determination; the modeling of high-excitation SO$_2$ lines is cumbersome, possibly related to line misidentifications or problems with the collisional rates. The SiO isotopic ratios ($^{29}$SiO/$^{28}$SiO and $^{30}$SiO/$^{28}$SiO) point toward an enhancement in $^{28}$SiO compared to results of classical stellar evolution codes. Predictions for H$_2$O lines in the spectral range of the Herschel/HIFI mission are performed. [abbreviated], Comment: 24 pagees, accepted for publication in Astronomy & Astrophysics

The location of B supergiants in the Hertzsprung-Russell diagram (HRD) represents a long-standing problem in massive star evolution. Here we propose their nature may be revealed utilising their rotational properties, and we highlight a steep drop in massive star rotation rates at an effective temperature of 22000 K. We discuss two potential explanations for it. On the one hand, the feature might be due to the end of the main sequence, which could potentially constrain the core overshooting parameter. On the other hand, the feature might be the result of enhanced mass loss at the predicted location of the bi-stability jump. We term this effect "bi-stability breaking" and discuss its potential consequences for the evolution of massive stars., Comment: Accepted by A&A Letters (4 pages, 5 figures); typos corrected

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], Comment: 19 pages, 11 figures, accepted by ApJ

The Tarantula Survey is an ambitious ESO Large Programme that has obtained multi-epoch spectroscopy of over 1,000 massive stars in the 30 Doradus region of the Large Magellanic Cloud. Here we introduce the scientific motivations of the survey and give an overview of the observational sample. Ultimately, quantitative analysis of every star, paying particular attention to the effects of rotational mixing and binarity, will be used to address fundamental questions in both stellar and cluster evolution., Comment: 6 pages, to appear in the proceedings of IAU Symposium 266 "Star Clusters - Basic Galactic Building Blocks Throughout Time and Space", Rio de Janeiro, Brazil, August 2009

We study the convection zones in the outer envelope of hot massive stars which are caused by opacity peaks associated with iron and helium ionization. We determine the occurrence and properties of these convection zones as function of the stellar parameters. We then confront our results with observations of OB stars. A stellar evolution code is used to compute a grid of massive star models at different metallicities. In these models, the mixing length theory is used to characterize the envelope convection zones. We find the iron convection zone (FeCZ) to be more prominent for lower surface gravity, higher luminosity and higher initial metallicity. It is absent for luminosities below about $10^{3.2}\Lsun$, $10^{3.9}\Lsun$, and $10^{4.2}$\Lsun$ for the Galaxy, LMC and SMC, respectively. We map the strength of the FeCZ on the Hertzsprung-Russell diagram for three metallicities, and compare this with the occurrence of observational phenomena in O stars: microturbulence, non-radial pulsations, wind clumping, and line profile variability. The confirmation of all three trends for the FeCZ as function of stellar parameters by empirical microturbulent velocities argues for a physical connection between sub-photospheric convective motions and small scale stochastic velocities in the photosphere of O- and B-type stars. We further suggest that clumping in the inner parts of the winds of OB stars could be caused by the same mechanism, and that magnetic fields produced in the FeCZ could appear at the surface of OB stars as diagnosed by discrete absorption components in ultraviolet absorption lines., Comment: Accepted for publication in Astronomy and Astrophysics

In this paper we present two efficient implementations of the diffusion approximation to be employed in Monte Carlo computations of radiative transfer in dusty media of massive circumstellar disks. The aim is to improve the accuracy of the computed temperature structure and to decrease the computation time. The accuracy, efficiency and applicability of the methods in various corners of parameter space are investigated. The effects of using these methods on the vertical structure of the circumstellar disk as obtained from hydrostatic equilibrium computations are also addressed. Two methods are presented. First, an energy diffusion approximation is used to improve the accuracy of the temperature structure in highly obscured regions of the disk, where photon counts are low. Second, a modified random walk approximation is employed to decrease the computation time. This modified random walk ensures that the photons that end up in the high-density regions can quickly escape to the lower density regions, while the energy deposited by these photons in the disk is still computed accurately. A new radiative transfer code, MCMax, is presented in which both these diffusion approximations are implemented. These can be used simultaneously to increase both computational speed and decrease statistical noise. We conclude that the diffusion approximations allow for fast and accurate computations of the temperature structure, vertical disk structure and observables of very optically thick circumstellar disks., Comment: Accepted for publication in A&A

We have previously analysed the spectra of 135 early B-type stars in the LMC and found several groups of stars that have chemical compositions that conflict with the theory of rotational mixing. Here we extend this study to Galactic and SMC metallicities with the analysis of ~50 Galactic and ~100 SMC early B-type stars with rotational velocities up to ~300km/s. The surface nitrogen abundances are utilised as a probe of the mixing process. In the SMC, we find a population of slowly rotating nitrogen-rich stars amongst the early B type core-hydrogen burning stars, similar to the LMC. In the Galactic sample we find no significant enrichment amongst the core hydrogen-burning stars, which appears to be in contrast with the expectation from both rotating single-star and close binary evolution models. However, only a small number of the rapidly rotating stars have evolved enough to produce a significant nitrogen enrichment, and these may be analogous to the non-enriched rapid rotators previously found in the LMC sample. Finally, in each metallicity regime, a population of highly enriched supergiants is observed, which cannot be the immediate descendants of core-hydrogen burning stars. Their abundances are, however, compatible with them having gone through a previous red supergiant phase. Together, these observations paint a complex picture of the nitrogen enrichment in massive main sequence and supergiant stellar atmospheres, where age and binarity cause crucial effects. Whether rotational mixing is required to understand our results remains an open question at this time, but could be answered by identifying the true binary fraction in those groups of stars that do not agree with single-star evolutionary models (abridged)., Comment: Accepted paper - 86 pages with tables and figures

We present a parameter study of self-consistent models of protoplanetary disks around Herbig AeBe stars. We use the code developed by Dullemond and Dominik, which solves the 2D radiative transfer problem including an iteration for the vertical hydrostatic structure of the disk. This grid of models will be used for several studies on disk emission and mineralogy in followup papers. In this paper we take a first look on the new models, compare them with previous modeling attempts and focus on the effects of various parameters on the overall structure of the SED that leads to the classification of Herbig AeBe stars into two groups, with a flaring (group I) or self-shadowed (group II) SED. We find that the parameter of overriding importance to the SED is the total mass in grains smaller than 25um, confirming the earlier results by Dullemond and Dominik. All other parameters studied have only minor influences, and will alter the SED type only in borderline cases. We find that there is no natural dichotomy between group I and II. From a modeling point of view, the transition is a continuous function of the small dust mass. We also show that moderate grain growth produces spectra with weak or no 1um feature, both for flaring (Group I) and non-flaring (Group II) sources. The fact that sources with weak features have been found mostly in Group I sources is therefore surprising and must be due to observational biases or evolutionary effects., Comment: 12 pages, 10 figures, Astronomy and Astrophysics, in press

We present results from the first extensive study of convection zones in the envelopes of hot massive stars, which are caused by opacity peaks associated with iron and helium ionization. These convective regions can be located very close to the stellar surface. Recent observations of microturbulence in massive stars from the VLT-Flames survey are in good agreement with our predictions concerning the occurrence and the strength of sub-surface convection in hot stars. We argue further that convection close to the surface may trigger clumping at the base of the stellar wind of massive stars., Comment: to appear in Comm. in Astroseismology - Proceedings of the 38th LIAC/HELAS-ESTA/BAG, 2008

We have used VLT FLAMES data to constrain the uncertain physics of rotational mixing in stellar evolution models. We have simulated a population of single stars and find two groups of observed stars that cannot be explained: (1) a group of fast rotating stars which do not show evidence for rotational mixing and (2) a group of slow rotators with strong N enrichment. Binary effects and fossil magnetic fields may be considered to explain those two groups. We suggest that the element boron could be used to distinguish between rotational mixing and the binary scenario. Our single star population simulations quantify the expected amount of boron in fast and slow rotators and allow a comparison with measured nitrogen and boron abundances in B-stars., Comment: to appear in Comm. in Astroseismology - Contribution to the Proceedings of the 38th LIAC, 2008

In order to deduce properties of dust in astrophysical environments where dust growth through aggregation is important, knowledge of the way aggregated particles interact with radiation, and what information is encoded in the thermal radiation they emit, is needed. The emission characteristics are determined by the size and structure of the aggregate and the composition and shape of the constituents. We thus aim at performing computations of compositionally inhomogeneous aggregates composed of irregularly shaped constituents. In addition we aim at developing an empirical recipe to compute the optical properties of such aggregates in a fast and accurate manner. We performed CDA computations for aggregates of irregularly shaped particles with various compositions. The constituents of the aggregate are assumed to be in the Rayleigh regime, and in addition we assume that the dominant interaction of the aggregate constituents is through dipole-dipole interactions. We computed the spectral structure of the emission efficiency in the 10 micron region for aggregates with 30% amorphous carbon and 70% silicates by volume with various fractions of crystalline and amorphous components. We find that the spectral appearance of the various components of the aggregate are very different and depend on their abundances. Most notably, materials that have a very low abundance appear spectroscopically as if they were in very small grains, while more abundant materials appear, spectroscopically to reside in larger grains. We construct a fast empirical approximate method, based on the idea of an effective medium approximation, to construct the spectra for these aggregates which almost perfectly reproduces the more exact computations. This new method is fast enough to be easily implemented in fitting procedures., Comment: Accepted for publication in A&A

MWC 297 is a young massive nearby B[e] star. The central star has a large projected rotational velocity of 350 km/s. Despite the wealth of published observations, the nature of this object and its dust-rich surroundings is not well understood. With the present paper, we shed light on the geometrical structure of the circumstellar matter which produces the near- to mid-infrared flux excess, and construct an overall image of the source's appearance and evolutionary status. The H-, K- and N-band brightness distribution of MWC 297 is probed with the ESO interferometric spectrographs AMBER and MIDI. We have obtained visibility measurements on 3 AMBER and 12 MIDI baselines, covering a wide range of spatial frequencies. We have reconstructed the brightness distribution in H, K and N with a geometric model consisting of three Gaussian disks with different extent and brightness temperature. This model can account for the entire near- to mid-IR emission of MWC 297. The near- and mid-IR emission, including the silicate emission at 10 micron, emanates from a very compact region (FWHM < 1.5 AU) around the central star. We argue that the circumstellar matter in the MWC 297 system is organized in a disk, seen under moderate (i < 40 deg) inclination. The disk displays no inner emission-free gap at the resolution of our interferometric data. The low inclination of the disk implies that the actual rotational velocity of the star exceeds its critical velocity. We discuss the impact of this result in terms of the formation of high-mass stars, and the main-sequence evolution of classical Be stars., Comment: accepted by Astronomy & Astrophysics (16/04/2008), 14 pages, 8 figures

Near-infrared imaging surveys of high-mass star-forming regions reveal an amazingly complex interplay between star formation and the environment (Churchwell et al. 2006; Alvarez et al. 2004). By means of near-IR spectroscopy the embedded massive young stars can be characterized and placed in the context of their birth site. However, so far spectroscopic surveys have been hopelessly incomplete, hampering any systematic study of these very young massive stars. New integral field instrumentation available at ESO has opened the possibility to take a huge step forward by obtaining a full spectral inventory of the youngest massive stellar populations in star-forming regions currently accessible. Simultaneously, the analysis of the extended emission allows the characterization of the environmental conditions. The Formation and Early Evolution of Massive Stars (FEMS) collaboration aims at setting up a large observing campaign to obtain a full census of the stellar content, ionized material, outflows and PDR's over a sample of regions that covers a large parameter space. Complementary radio, mm and infrared observations will be used for the characterization of the deeply embedded population. For the first eight regions we have obtained 40 hours of SINFONI observations. In this contribution, we present the first results on three regions that illustrate the potential of this strategy., Comment: To appear in ASP Conf. Proceedings of "Massive Star Formation: Observations confront Theory", H. Beuther et al. (eds.), held in Heidelberg, September 2007

The VLT-FLAMES Survey of Massive Stars was an ESO Large Programme to understand rotational mixing and stellar mass-loss in different metallicity environments, in order to better constrain massive star evolution. We gathered high-quality spectra of over 800 stars in the Galaxy and in the Magellanic Clouds. A sample of this size is unprecedented, enabled by the first high-resolution, wide-field, multi-object spectrograph on an 8-m telescope. We developed spectral analysis techniques that, in combination with non-LTE, line-blanketed model atmospheres, were used to quantitatively characterise every star. The large sample, combined with the theoretical developments, has produced exciting new insights into the evolution of the most massive stars., Comment: 7 pages, from the March 2008 edition of The ESO Messenger (vol. 131, p25). A fully formatted version is available at: http://www.eso.org/sci/publications/messenger/

Rotation has become an important element in evolutionary models of massive stars, specifically via the prediction of rotational mixing. Here, we study a sample of stars, including rapid rotators, to constrain such models and use nitrogen enrichments as a probe of the mixing process. Chemical compositions (C, N, O, Mg and Si) have been estimated for 135 early B-type stars in the Large Magellanic Cloud with projected rotational velocities up to ~300km/s using a non-LTE TLUSTY model atmosphere grid. Evolutionary models, including rotational mixing, have been generated attempting to reproduce these observations by adjusting the overshooting and rotational mixing parameters and produce reasonable agreement with 60% of our core hydrogen burning sample. We find (excluding known binaries) a significant population of highly nitrogen enriched intrinsic slow rotators vsini less than 50km/s incompatible with our models ~20% of the sample). Furthermore, while we find fast rotators with enrichments in agreement with the models, the observation of evolved (log g less than 3.7dex) fast rotators that are relatively unenriched (a further ~20% of the sample) challenges the concept of rotational mixing. We also find that 70% of our blue supergiant sample cannot have evolved directly from the hydrogen burning main-sequence. We are left with a picture where invoking binarity and perhaps fossil magnetic fields are required to understand the surface properties of a population of massive main sequence stars., Comment: ApJL. 10 pages, 1 figure. Updated to match accepted version

We study the evolution and fate of solar composition supermassive stars in the mass range 60 - 1000 Msun. Our study is relevant both for very massive objects observed in young stellar complexes as well as supermassive stars that may potentially form through runaway stellar collisions. We predict the outcomes of stellar evolution employing a mass-loss prescription that is consistent with the observed Hertzsprung-Russell Diagram location of the most massive stars. We compute a series of stellar models with an appropriately modified version of the Eggleton evolutionary code. We find that super-massive stars with initial masses up to 1000 Msun end their lives as objects less massive than ~150 Msun. These objects are expected to collapse into black holes (with M < 70 Msun) or explode as pair-instability supernovae. We argue that if ultraluminous X-ray sources (ULXs) contain intermediate-mass black holes, these are unlikely to be the result of runaway stellar collisions in the cores of young clusters., Comment: 26 pages, 15 figures, accepted by A&A. small language corrections

Low and intermediate mass stars lose a significant fraction of their mass through a dust-driven wind during the Asymptotic Giant Branch (AGB) phase. Recent studies show that winds from late-type stars are far from being smooth. Mass-loss variations occur on different time scales, from years to tens of thousands of years. The variations appear to be particularly prominent towards the end of the AGB evolution. The occurrence, amplitude and time scale of these variations are still not well understood. The goal of our study is to gain insight into the structure of the circumstellar envelope (CSE) of WX Psc and map the possible variability of the late-AGB mass-loss phenomenon. We have performed an in-depth analysis of the extreme infrared AGB star WX Psc by modeling (1) the CO J=1-0 through 7-6 rotational line profiles and the full spectral energy distribution (SED) ranging from 0.7 to 1300 micron. We hence are able to trace a geometrically extended region of the CSE. Both mass-loss diagnostics bear evidence of the occurrence of mass-loss modulations during the last ~2000 yr. In particular, WX Psc went through a high mass-loss phase (Mdot~5e-5 Msun/yr) some 800 yr ago. This phase lasted about 600 yr and was followed by a long period of low mass loss (Mdot~5e-8 Msun/yr). The present day mass-loss rate is estimated to be ~6e-6 Msun/yr. The AGB star WX Psc has undergone strong mass-loss rate variability on a time scale of several hundred years during the last few thousand years. These variations are traced in the strength and profile of the CO rotational lines and in the SED. We have consistently simulated the behaviour of both tracers using radiative transfer codes that allow for non-constant mass-loss rates., Comment: 12 pages, accepted for publication in A&A

We present the results of Monte Carlo mass-loss predictions for massive stars covering a wide range of stellar parameters. We critically test our predictions against a range of observed mass-loss rates -- in light of the recent discussions on wind clumping. We also present a model to compute the clumping-induced polarimetric variability of hot stars and we compare this with observations of Luminous Blue Variables, for which polarimetric variability is larger than for O and Wolf-Rayet stars. Luminous Blue Variables comprise an ideal testbed for studies of wind clumping and wind geometry, as well as for wind strength calculations, and we propose they may be direct supernova progenitors., Comment: 3 pages, 3 figures, to appear in the proceedings of workshop 'Clumping in Hot Star Winds', eds. W.-R. Hamann, A. Feldmeier, & L. Oskinova

[Abridged] We present a comprehensive study of the metallicity dependence of the mass-loss rates in stationary stellar winds of hot massive stars. Assuming a power-law dependence of mass loss on metallicity, Mdot \propto Z^{m}, and adopting a theoretical relation between the terminal velocity and metallicity, v_inf \propto Z^{0.13} (Leitherer et al.), we find m = 0.83 +/- 0.16 for non-clumped outflows from an analysis of the wind momentum luminosity relation (WLR) for stars more luminous than 10^{5.2} Lsun. Within the errors, this result agrees with the prediction of m = 0.69 +/- 0.10 from Vink et al. For the high luminosity stars we find the mass loss rates to be greater than the predictions, implying wind clumping factors in their line-forming regions of ~4. For lower luminosity stars, the winds are so weak that their strengths cannot be reliably derived from optical lines, and one must rely on analysis of UV lines. In the low-luminosity domain the Galactic WLR is found to be much steeper than expected from theory, leading to a discrepancy between UV mass-loss rates and the predictions by a factor 100 at luminosities of L ~ 10^{4.75} Lsun, the origin of which is unknown. We emphasize that even if the current mass-loss rates of hot luminous stars are overestimated as a result of wind clumping, the degree of clumping would likely to be independent of metallicity, so the scalings derived in this study are expected to remain correct., Comment: Accepted by A&A, 17 pages

Pair creation supernovae (PCSN) are thought to be produced from very massive low metallicity stars. The spectacularly bright SN 2006gy does show signatures expected from PCSNe. Here, we investigate the metallicity threshold below which PCSN can form and estimate their occurrence rate. We perform stellar evolution calculations for stars of 150$\mso$ and 250$\mso$ of low metallicity (Z$_{\odot}$/5 and Z$_{\odot}$/20), and analyze their mass loss rates. We find that the bifurcation between quasi-chemically homogeneous evolution for fast rotation and conventional evolution for slower rotation, which has been found earlier for massive low metallicity stars, persists in the mass range considered here. Consequently, there are two separate PCSN progenitor types: (I) Fast rotators produce PCSNe from very massive Wolf-Rayet stars, and (II) Slower rotators that generate PCSNe in hydrogen-rich massive yellow hypergiants. We find that hydrogen-rich PCSNe could occur at metallicities as high as Z$_{\odot}$/3, which -- assuming standard IMFs are still valid to estimate their birth rates -- results in a rate of about one PCSN per 1000 supernovae in the local universe, and one PCSN per 100 supernovae at a redshift of $z=5$. PCSNe from WC-type Wolf-Rayet stars are restricted to much lower metallicity., Comment: submitted to A&A

[Abridged] We have studied the optical spectra of 28 O- and early B-type stars in the Large Magellanic Cloud, 22 of which are associated with the young star-forming region N11. Stellar parameters are determined using an automated fitting method, combining the stellar atmosphere code FASTWIND with the genetic-algorithm optimisation routine PIKAIA. Results for stars in the LH9 and LH10 associations of N11 are consistent with a sequential star formation scenario, in which activity in LH9 triggered the formation of LH10. Our sample contains four stars of spectral type O2, of which the hottest is found to be ~49-54 kK (cf. ~45-46 kK for O3 stars). The masses of helium-enriched dwarfs and giants are systematically lower than those implied by non-rotating evolutionary tracks. We interpret this as evidence for efficient rotationally-enhanced mixing, leading to the surfacing of primary helium and to an increase of the stellar luminosity. This result is consistent with findings for SMC stars by Mokiem et al. For bright giants and supergiants no such mass-discrepancy is found, implying that these stars follow tracks of modestly (or non-)rotating objects. Stellar mass-loss properties were found to be intermediate to those found in massive stars in the Galaxy and the SMC, and comparisons with theoretical predictions at LMC metallicity yielded good agreement over the luminosity range of our targets, i.e. 5.0 < log L/L(sun) < 6.1.

Massive stars and supernovae (SNe) have a huge impact on their environment. Despite their importance, a comprehensive knowledge of which massive stars produce which SNe is hitherto lacking. We use a Monte Carlo method to predict the mass-loss rates of massive stars in the Hertzsprung-Russell Diagram (HRD) covering all phases from the OB main sequence, the unstable Luminous Blue Variable (LBV) stage, to the final Wolf-Rayet (WR) phase. Although WR produce their own metals, a strong dependence of the mass-loss rate on the initial iron abundance is found at sub-solar metallicities (1/10 -- 1/100 solar). This may present a viable mechanism to prevent the loss of angular momentum by stellar winds, which could inhibit GRBs occurring at solar metallicities -- providing a significant boost to the collapsar model. Furthermore, we discuss recently reported quasi-sinusoidal modulations in the radio lightcurves of SNe 2001ig and 2003bg. We show that both the sinusoidal behaviour and the recurrence timescale of these modulations are consistent with the predicted mass-loss behaviour of LBVs. We discuss potential ramifications for the ``Conti'' scenario for massive star evolution., Comment: 8 pages, 6 figures, conference talk. To be published in the proceedings of "Mass loss from stars and the evolution of stellar clusters"

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

We report on observations of circumstellar disks around young stars that have been obtained with the MIDI instrument, which is mounted on the VLT Interferometer and operates in the 10 micrometer atmospheric window. The maximum spatial resolution of 5 milli-arcsec corresponds to sub-AU scales at the distance to nearby star formation regions. Thus, we can study the disks on the spatial scales at which important processes occur, such as accretion, dust processing, and planet formation. The main results obtained so far can be summarized as follows: 1. The measured interferometric visibilities are in good qualitative agreement with those predicted by models of circumstellar disks. In particular, a predicted correlation between the strength of the far-infrared excess and the spatial structure of the disk is confirmed by direct measurements; 2. In several objects strong evidence for deviations from circular symmetry is present, indicating that an inclined disk is indeed the dominant component seen in the mid-infrared; 3. The dust properties are not uniform over the disk, but are instead a strong function of distance to the central star. The dust in the innermost disk regions is observed to be more ``processed'' than the dust further out, both in Herbig Ae star disks and in those around T-Tauri stars., Comment: 19 pages, 12 figures, invited talk at SPIE conference "Astronomical telescopes and instrumentation", Orlando, may 2006

We have studied the optical spectra of a sample of 31 O- and early B-type stars in the Small Magellanic Cloud, 21 of which are associated with the young massive cluster NGC 346. Stellar parameters are determined using an automated fitting method. Comparison with predictions of stellar evolution that account for stellar rotation does not result in a unique age, though most stars are best represented by an age of 1-3 Myr. The present day v_sini distribution of the 21 dwarf stars in our sample is consistent with an underlying rotational velocity (v_r) distribution that can be characterised by a mean velocity of about 160-190 km/s and an effective half width of 100-150 km/s. The v_r distribution must include a small percentage of slowly rotating stars. If predictions of the time evolution of the equatorial velocity for massive stars within the environment of the SMC are correct, the young age of the cluster implies that this underlying distribution is representative for the initial rotational velocity distribution. The location in the Hertzsprung-Russell diagram of the stars showing helium enrichment is in qualitative agreement with evolutionary tracks accounting for rotation, but not for those ignoring v_r. The mass loss rates of the SMC objects having luminosities of log L/L_sun > 5.4 are in excellent agreement with predictions. However, for lower luminosity stars the winds are too weak to determine M_dot accurately from the optical spectrum. Two of three spectroscopically classified Vz stars from our sample are located close to the theoretical zero age main sequence, as expected., Comment: 35 pages, accepted for publication in A&A

Context: Mass loss plays a dominant role in the evolution of low mass stars while they are on the Asymptotic Giant Branch (AGB). The gas and dust ejected during this phase are a major source in the mass budget of the interstellar medium. Recent studies have pointed towards the importance of variations in the mass-loss history of such objects. Aims: By modelling the full line profile of low excitation CO lines emitted in the circumstellar envelope, we can study the mass-loss history of AGB stars. Methods: We have developed a non-LTE radiative transfer code, which calculates the velocity structure and gas kinetic temperature of the envelope in a self-consistent way. The resulting structure of the envelope provides the input for the molecular line radiative calculations which are evaluated in the comoving frame. The code allows for the implementation of modulations in the mass-loss rate. This code has been benchmarked against other radiative transfer codes and is shown to perform well and efficiently. Results: We illustrate the effects of varying mass-loss rates in case of a superwind phase. The model is applied to the well-studied case of VY CMa. We show that both the observed integrated line strengths as the spectral structure present in the observed line profiles, unambiguously demonstrate that this source underwent a phase of high mass loss (~ 3.2E-4 Msun/yr) some 1000 yr ago. This phase took place for some 100 yr, and was preceded by a low mass-loss phase (~ 1E-6 Msun/yr) taking some 800 yr. The current mass-loss rate is estimated to be in the order of 8E-5 Msun/yr. Conclusions: In this paper, we demonstrate that both the relative strength of the CO rotational line profiles and the (non)-occurrence of spectral structure in the profile offer strong diagnostics to pinpoint the mass-loss history., Comment: 18 pages, 10 figures, accepted for publication in "Astronomy & Astrophysics"

We model the 10 micron absorption spectra of nonspherical particles composed of amorphous silicate. We consider two classes of particles, compact ones and fractal aggregates composed of homogeneous spheres. For the compact particles we consider Gaussian random spheres with various degrees of non-sphericity. For the fractal aggregates we compute the absorption spectra for various fractal dimensions. The 10 micron spectra are computed for ensembles of these particles in random orientation using the well-known Discrete Dipole Approximation. We compare our results to spectra obtained when using volume equivalent homogeneous spheres and to those computed using a porous sphere approximation. We conclude that, in general, nonspherical particles show a spectral signature that is similar to that of homogeneous spheres with a smaller material volume. This effect is overestimated when approximating the particles by porous spheres with the same volume filling fraction. For aggregates with fractal dimensions typically predicted for cosmic dust, we show that the spectral signature characteristic of very small homogeneous spheres (with a volume equivalent radius r_V<0.5 micron) can be detected even in very large particles. We conclude that particle sizes are underestimated when using homogeneous spheres to model the emission spectra of astronomical sources. In contrast, the particle sizes are severely overestimated when using equivalent porous spheres to fit observations of 10 micron silicate emission., Comment: Accepted for publication in A&A

We have performed a pilot study of mass loss predictions for late-type Wolf-Rayet (WR) stars as a function of metal abundance, over a range between 10^{-5} < (Z/Zsun) < 10. We find that the winds of nitrogen-rich Wolf-Rayet stars are dominated by iron lines, with a dependence of mass loss on Z similar to that of massive OB stars. For more evolved, carbon-rich, WR stars the wind strength is also found to be dependent on the Fe abundance, so that they depend on the chemical environment of the host galaxy, but with a mass loss metallicity dependence that is less steep than for OB stars. Our finding that WR mass loss is Z-dependent is a new one, with important consequences for black hole formation and X-ray population studies in external galaxies. A further finding of our study is that the Z dependence of C-rich WR stars becomes weaker at metallicities below Z/Zsun < 1/10, and mass loss no longer declines once the metal abundance drops below (Z/Zsun) = 10^{-3}. This is the result of an increased importance of radiative driving by intermediate mass elements, such as carbon. In combination with rapid rotation and/or proximity to the Eddington limit -- likely to be relevant for massive Population III stars -- this effect may indicate a role for mass loss in the appearance and evolution of these objects, as well as a potential role for stellar winds in enriching the intergalactic medium of the early Universe., Comment: Accepted by Astronomy & Astrophysics (11 pages)

We present the first automated fitting method for the quantitative spectroscopy of O- and early B-type stars with stellar winds. The method combines the non-LTE stellar atmosphere code FASTWIND from Puls et al. (2005) with the genetic algorithm based optimizing routine PIKAIA from Charbonneau (1995), allowing for a homogeneous analysis of upcoming large samples of early-type stars (e.g. Evans et al. 2005). In this first implementation we use continuum normalized optical hydrogen and helium lines to determine photospheric and wind parameters. We have assigned weights to these lines accounting for line blends with species not taken into account, lacking physics, and/or possible or potential problems in the model atmosphere code. We find the method to be robust, fast, and accurate. Using our method we analysed seven O-type stars in the young cluster Cyg OB2 and five other Galactic stars with high rotational velocities and/or low mass loss rates (including 10lac, zeta Oph, and tau Sco) that have been studied in detail with a previous version of FASTWIND. The fits are found to have a quality that is comparable or even better than produced by the classical ``by eye'' method. We define errorbars on the model parameters based on the maximum variations of these parameters in the models that cluster around the global optimum. Using this concept, for the investigated dataset we are able to recover mass-loss rates down to ~6d-8 Msun/yr to within an error of a factor of two. Comparison of our derived spectroscopic masses with those derived from stellar evolutionary models are in very good agreement, i.e. based on the limited sample that we have studied we do not find indications for a mass discrepancy., Comment: 24 pages, 19 figures, accepted by A&A

We discuss the composition and size distribution of the dust in the coma of comet Hale-Bopp. We do this by fitting simultaneously the infrared emission spectrum measured by the infrared space observatory (ISO) and the measured degree of linear polarization of scattered light at various phase angles and 12 different wavelengths. The effects of particle shape on the modeled optical properties of the dust grains are taken into account. We constrain our fit by forcing the abundances of the major rock forming chemical elements to be solar. The infrared spectrum at long wavelengths reveals that large grains are needed in order to fit the spectral slope. The size and shape distribution we employ allows us to estimate the sizes of the crystalline silicates. The ratios of the strength of various forsterite features show that the crystalline silicate grains in Hale-Bopp must be submicron sized. We exclude the presence of large crystalline silicate grains in the coma. Because of this lack of large crystalline grains combined with the fact that we do need large amorphous grains to fit the emission spectrum at long wavelengths, we need only approximately 4% of crystalline silicates by mass. After correcting for possible hidden crystalline material included in large amorphous grains, our best estimate of the total mass fraction of crystalline material is approximately 7.5%, significantly lower than deduced in previous studies in which the typical derived crystallinity is 20-30%. The implications of this on the possible origin and evolution of the comet are discussed. The crystallinity we observe in Hale-Bopp is consistent with the production of crystalline silicates in the inner solar system by thermal annealing and subsequent radial mixing to the comet forming region., Comment: Accepted for publication in Icarus

We provide a theoretical foundation for the statistical approach for computing the absorption properties of particles in the Rayleigh domain. We present a general method based on the Discrete Dipole Approximation (DDA) to compute the absorption and scattering properties of particles in the Rayleigh domain. The method allows to separate the geometrical aspects of a particle from its material properties. Doing the computation of the optical properties of a particle once, provides them for any set of refractive indices, wavelengths and orientations. This allows for fast computations of e.g. absorption spectra of arbitrarily shaped particles. Other practical applications of the method are in the interpretation of atmospheric and radar measurements as well as computations of the scattering matrix of small particles as a function of the scattering angle. In the statistical approach, the optical properties of irregularly shaped particles are represented by the average properties of an ensemble of particles with simple shapes. We show that the absorption cross section of an ensemble of arbitrarily shaped particles with arbitrary orientations can always be uniquely represented by the average absorption cross section of an ensemble of spheroidal particles with the same composition and fixed orientation. This proves for the first time that the statistical approach is generally viable in the Rayleigh domain., Comment: Accepted for publication in JQSRT

We introduce a new survey of massive stars in the Galaxy and the Magellanic Clouds using the Fibre Large Array Multi-Element Spectrograph (FLAMES) instrument at the Very Large Telescope (VLT). Here we present observations of 269 Galactic stars with the FLAMES-Giraffe Spectrograph (R ~ 25,000), in fields centered on the open clusters NGC 3293, NGC 4755 and NGC 6611. These data are supplemented by a further 50 targets observed with the Fibre-Fed Extended Range Optical Spectrograph (FEROS, R = 48,000). Following a description of our scientific motivations and target selection criteria, the data reduction methods are described; of critical importance the FLAMES reduction pipeline is found to yield spectra that are in excellent agreement with less automated methods. Spectral classifications and radial velocity measurements are presented for each star, with particular attention paid to morphological peculiarities and evidence of binarity. These observations represent a significant increase in the known spectral content of NGC 3293 and NGC 4755, and will serve as standards against which our subsequent FLAMES observations in the Magellanic Clouds will be compared., Comment: 26 pages, 9 figures (reduced size). Accepted for publication in A&A. A copy with full res. figures is available from http://www.ing.iac.es/~cje/flames_mw.ps.gz. Minor changes following correction of proofs

We present spectroscopic observations of a large sample of Herbig Ae stars in the 10 micrometer spectral region. We perform compositional fits of the spectra based on properties of homogeneous as well as inhomogeneous spherical particles, and derive the mineralogy and typical grain sizes of the dust responsible for the 10 $\mu$m emission. Several trends are reported that can constrain theoretical models of dust processing in these systems: i) none of the sources consists of fully pristine dust comparable to that found in the interstellar medium, ii) all sources with a high fraction of crystalline silicates are dominated by large grains, iii) the disks around more massive stars (M >~ 2.5 M_sun, L >~ 60 L_sun) have a higher fraction of crystalline silicates than those around lower mass stars, iv) in the subset of lower mass stars (M <~ 2.5 M_sun) there is no correlation between stellar parameters and the derived crystallinity of the dust. The correlation between the shape and strength of the 10 micron silicate feature reported by van Boekel et al. (2003) is reconfirmed with this larger sample. The evidence presented in this paper is combined with that of other studies to present a likely scenario of dust processing in Herbig Ae systems. We conclude that the present data favour a scenario in which the crystalline silicates are produced in the innermost regions of the disk, close to the star, and transported outward to the regions where they can be detected by means of 10 micron spectroscopy. Additionally, we conclude that the final crystallinity of these disks is reached very soon after active accretion has stopped., Comment: 22 pages, 14 figures, accepted for publication in A&A. Note: this submission was replaced on 26.04.2005: we used incorrect terminology in figure 6 and the discussion of this figure. The vertical axis label of figure 6 has been corrected and now reads "Normalized 11.3/9.8 Flux Ratio", in the discussion of this figure (section 4.2) "continuum subtracted" has been replaced by "normalized"

In this paper we study the combined effects of size and shape of small solid state particles on the absorption, emission and scattering characteristics. We use the statistical approach to calculate these optical properties. In this approach the average optical properties of an ensemble of particles in random orientation are represented by the average optical properties of an ensemble of simple shapes. The validity of this approach is studied in detail for a uniform distribution of hollow spheres where the fractional volume of the central inclusion is varied. We apply the results to two different areas of interest, i) infrared spectroscopy and ii) polarization of scattered light. The effects of particle size and shape on the optical characteristics are discussed. We compare the results using the distribution of hollow spheres with those obtained by using randomly oriented spheroids. Also we compare the results with observations and laboratory measurements. The distribution of hollow spheres is very successful in reproducing laboratory measurements of the scattering angle distribution of the degree of linear polarization for incident unpolarized light of randomly oriented irregular quartz particles. Furthermore, we show that we are able to derive the size distribution of dust grains by fitting the measured degree of linear polarization using computational result for hollow spheres. It is shown that the distribution of hollow spheres is a powerful tool for studying light scattering, absorption and emission by cosmic dust grains and in particular when large numbers of particle parameters need to be considered since the computational demand of the distribution of hollow spheres is small., Comment: 12 pages, A&A, in press

We present observations of the famous OH/IR star OH26.5+0.6 obtained using the Mid-Infrared Interferometric Instrument MIDI at the European Southern Observatory (ESO) Very Large Telescope Interferometer VLTI. The emission of the dusty envelope, spectrally dispersed at a resolution of 30 from 8 to 13.5 micron, appears resolved by a single dish UT telescope. In particular the angular diameter increases strongly within the silicate absorption band. Moreover an acquisition image taken at 8.7 micron exhibits, after deconvolution, a strong asymmetry. The axis ratio is 0.75+/-0.07 with the FWHM of the major and minor axis which are 286mas and 214mas respectively. The measured PA angle, 95 degrees +/-6 degrees is reminiscent of the asymmetry in the OH maser emission detected at 1612MHz by Bowers & Johnston (1990) for this star. In interferometric mode the UT1-UT3 102m baseline was employed to detect the presence of the star. No fringes have been found with a detection threshold estimated to be of the order of 1% of the total flux of the source, i.e. 5-8 Jy. These observations were carried out during the phase of maximum luminosity of the star, when the dust shell is more diluted and therefore the chance to detect the central source maximized. We modeled the dusty environment based on the work of Justannont et al. (1996). In particular, the failure to detect fringes provides strong constraints on the opacities in the inner regions of the dust shell or in the close vicinity of the star., Comment: Accepted in A&A

The core of the nebula surrounding Eta Carinae has been observed with the VLT Adaptive Optics system NACO and with the interferometer VLTI/MIDI to constrain spatially and spectrally the warm dusty environment and the central object. In particular, narrow-band images at 3.74 and 4.05 micron reveal the butterfly shaped dusty environment close to the central star with unprecedented spatial resolution. A void whose radius corresponds to the expected sublimation radius has been discovered around the central source. Fringes have been obtained in the Mid-IR which reveal a correlated flux of about 100Jy situated 0.3" south-east of the photocenter of the nebula at 8.7 micron, which corresponds with the location of the star as seen in other wavelengths. This correlated flux is partly attributed to the central object, and these observations provide an upper limit for the SED of the central source from 2.2 to 13.5 micron. Moreover, we have been able to spectrally disperse the signal from the nebula itself at PA=318 degree, i.e. in the direction of the bipolar nebula 310 degree) within the MIDI field of view of 3". A large amount of corundum (Al2O3) is discovered, peaking at 0.6-1.2" south-east from the star, whereas the dust content of the Weigelt blobs is dominated b silicates. We discuss the mechanisms of dust formation which are closely related to the geometry of this Butterfly nebulae.

We discuss an interesting feature of the distribution of luminous blue variables on the H-R diagram, and we propose a connection with the bistability jump in the winds of early-type supergiants. There appears to be a deficiency of quiescent LBVs on the S Dor instability strip at luminosities between log L/Lsun = 5.6 and 5.8. The upper boundary, is also where the temperature-dependent S Dor instability strip intersects the bistability jump at about 21,000 K. Due to increased opacity, winds of early-type supergiants are slower and denser on the cool side of the bistability jump, and we postulate that this may trigger optically-thick winds that inhibit quiescent LBVs from residing there. We conduct numerical simulations of radiation-driven winds for a range of temperatures, masses, and velocity laws at log L/Lsun=5.7 to see what effect the bistability jump should have. We find that for relatively low stellar masses the increase in wind density at the bistability jump leads to the formation of a modest to strong pseudo photosphere -- enough to make an early B-type star appear as a yellow hypergiant. Thus, the proposed mechanism will be most relevant for LBVs that are post-red supergiants. Yellow hypergiants like IRC+10420 and rho Cas occupy the same luminosity range as the ``missing'' LBVs, and show apparent temperature variations at constant luminosity. If these yellow hypergiants do eventually become Wolf-Rayet stars, we speculate that they may skip the normal LBV phase, at least as far as their apparent positions on the HR diagram are concerned., Comment: 20 pages, 4 figs, accepted by ApJ

We use a grid of 30 line-blanketed unified stellar photosphere and wind models for O-type stars; computed with the code CMFGEN in order to evaluate its potential in the near-infrared spectral domain. The grid includes dwarfs, giants and supergiants. We analyse the equivalent width behaviour of the 20 strongest lines of hydrogen and helium in spectral windows that can be observed using ground-based instrumentation and compare the results with observations. Our main findings are that: i) HeI/HeII line ratios in the J, H and K bands correlate well with the optical ratio employed in spectral classification, and can therefore be used to determine the spectral type; ii) in supergiant stars the transition from the stellar photosphere to the wind follows a shallower density gradient than the standard approach followed in our models, which can be mimicked by adopting a lower gravity in our prescription of the density stratification. iii) the Brackett gamma line poses a number of peculiar problems which partly might be related to wind clumping, and iv) the Brackett alpha line is an excellent mass-loss indicator. For the first and last item we provide quantitative calibrations., Comment: 14 pages, 7 figures, accepted by A&A

We present a systematic study of the effect of metallicity on the stellar spectral energy distribution (SED) of O main sequence (dwarf) stars, focussing on the hydrogen and helium ionizing continua, and on the optical and near-IR lines used for spectral classification. The spectra are based on non-LTE line blanketed atmosphere models with stellar winds calculated using the CMFGEN code of Hillier & Miller (1998). We draw the following conclusions. First, we find that the total number of Lyman photons emitted is almost independent of line blanketing effects and metallicity for a given effective temperature. This is because the flux that is blocked by the forest of metal lines at wavelengths shorter than 600 Angstrom is redistributed mainly within the Lyman continuum. Second, the spectral type, as defined by the ratio of the equivalent widths of HeI 4471 Angstrom and HeII 4542 Angstrom, is shown to depend noticeably on the microturbulent velocity in the atmosphere, on metallicity and, within the luminosity class of dwarfs, on gravity. Third, we confirm the decrease in the effective temperature for a given spectral type due to the inclusion of line blanketing recently found by e.g. Martins et al. (2002). Finally, we find that the SED below ~450 Angstrom is highly dependent on metallicity. This is reflected in the behaviour of nebular fine-structure line ratios such as [NeIII]/[NeII] 15.5/12.8 and [ArIII]/[ArII] 9.0/7.0 micron. This dependence complicates the use of these nebular ratios as diagnostic tools for the effective temperature determination of the ionizing stars in HII regions and for age dating of starburst regions in galaxies., Comment: 16 pages, 11 figures, accepted by A&A

We re-examine the nature of NGC2024-IRS2 in light of the recent discovery of the late O-type star, IRS2b, located 5 arcsec from IRS2. Using L-band spectroscopy, we set a lower limit of Av = 27.0 mag on the visual extinction towards IRS2. Arguments based on the nature of the circumstellar material, favor an Av of 31.5 mag. IRS2 is associated with the UCHII region G206.543-16.347 and the infrared source IRAS 05393-0156. We show that much of the mid-infrared emission towards IRS2, as well as the far infrared emission peaking at ~ 100 micron, do not originate in the direct surroundings of IRS2, but instead from an extended molecular cloud. Using new K-, L- and L'-band spectroscopy and a comprehensive set of infrared and radio continuum measurements from the literature, we apply diagnostics based on the radio slope, the strength of the infrared hydrogen recombination lines, and the presence of CO band-heads to constrain the nature and spatial distribution of the circumstellar material of IRS2. Using simple gaseous and/or dust models of prescribed geometry, we find strong indications that the infrared flux originating in the circumstellar material of IRS2 is dominated by emission from a dense gaseous disk with a radius of about 0.6 AU. At radio wavelengths the flux density distribution is best described by a stellar wind recombining at a radius of about 100 AU. Although NGC2024/IRS2 shares many similarities with BN-like objects, we do not find evidence for the presence of a dust shell surrounding this object. Therefore, IRS2 is likely more evolved., Comment: 15 pages, 10 figures. Accepted for publication in Astronomy & Astrophysics

We present submillimeter observations of rotational transitions of carbon monoxide from J = 2 -> 1 up to 7 -> 6 for a sample of Asymptotic Giant Branch stars and red supergiants. It is the first time that the high transitions J = 6 -> 5 and 7 -> 6 are included in such a study. With line radiative transfer calculations, we aim to determine the mass-loss history of these stars by fitting the CO line intensities. We find that the observed line intensities of the high transitions, including the J = 4 -> 3 transition, are significantly lower than the predicted values. We conclude that the physical structure of the outflow of Asymptotic Giant Branch stars is more complex than previously thought. In order to understand the observed line intensities and profiles, a physical structure with a variable mass-loss rate and/or a gradient in stochastic gas velocity is required. A case study of the AGB star WX Psc is performed. We find that the CO line strengths may be explained by variations in mass-loss on time scales similar to those observed in the separated arc-like structures observed around post-AGB stars. In addition, a gradient in the stochastic velocity may play a role. Until this has been sorted out fully, any mass loss determinations based upon single CO lines will remain suspect., Comment: accepted for publication by A&A

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., Comment: 7 pages, 4 figures. Figure 1 included as jpeg. Accepted for publication in Astronomy & Astrophysics

We present a detailed study of the circumstellar envelope of the post-asymptotic giant branch ``21 micron object'' HD 56126. We build a detailed dust radiative transfer model of the circumstellar envelope in order to derive the dust composition and mass, and the mass-loss history of the star. To model the emission of the dust we use amorphous carbon, hydrogenated amorphous carbon, magnesium sulfide and titanium carbide. We present a detailed parametrisation of the optical properties of hydrogenated amorphous carbon as a function of H/C content. The mid-infrared imaging and spectroscopy is best reproduced by a single dust shell from 1.2 to 2.6 arcsec radius around the central star. This shell originates from a short period during which the mass-loss rate exceeded 10^(-4) M_sun/yr. We find that the strength of the ``21'' micron feature poses a problem for the TiC identification. The low abundance of Ti requires very high absorption cross-sections in the ultraviolet and visible wavelength range to explain the strength of the feature. Other nano-crystalline metal carbides should be considered as well. We find that hydrogenated amorphous carbon in radiative equilibrium with the local radiation field does not reach a high enough temperature to explain the strength of the 3.3-3.4 and 6-9 micron hydrocarbon features relative to the 11-17 micron hydrocarbon features. We propose that the carriers of these hydrocarbon features are not in radiative equilibrium but are transiently heated to high temperature. We find that 2 per cent of the dust mass is required to explain the strength of the ``30'' micron feature, which fits well within the measured atmospheric abundance of Mg and S. This further strengthens the MgS identification of the ``30'' micron feature., Comment: 20 Pages, 10 Figures, accepted for publication in Astronomy and Astrophysics