Search

Magnetic fields can drastically change predictions of evolutionary models of massive stars via mass-loss quenching, magnetic braking, and efficient angular momentum transport, which we aim to quantify in this work. We use the MESA software instrument to compute an e xtensiv e main-sequence grid of stellar structure and evolution models, as well as isochrones, accounting for the effects attributed to a surface fossil magnetic field. The grid is densely populated in initial mass (3–60 M �), surface equatorial magnetic field strength (0–50 kG), and metallicity (representative of the Solar neighbourhood and the Magellanic Clouds). We use two magnetic braking and two chemical mixing schemes and compare the model predictions for slowly rotating, nitrogen- enriched (‘Group 2’) stars with observations in the Large Magellanic Cloud. We quantify a range of initial field strengths that allow for producing Group 2 stars and find that typical values (up to a few kG) lead to solutions. Between the subgrids, we find notable departures in surface abundances and evolutionary paths. In our magnetic models, chemical mixing is al w ays less efficient compared to non-magnetic models due to the rapid spin-down. We identify that quasi-chemically homogeneous main sequence evolution by efficient mixing could be prevented by fossil magnetic fields. We recommend comparing this grid of evolutionary models with spectropolarimetric and spectroscopic observations with the goals of (i) revisiting the derived stellar parameters of known magnetic stars, and (ii) observationally constraining the uncertain magnetic braking and chemical mixing schemes.

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

We present the 2-200 um Infrared Space Observatory (ISO) spectrum and 3.8-20 um ISAAC and TIMMI2 images of the extreme OH/IR star IRAS 16342-3814. Amorphous silicate absorption features are seen, together with crystalline silicate absorption features up to almost 45 um. No other OH/IR star is known to have crystalline silicate features in absorption up to these wavelengths. This suggests that IRAS 16342-3814 must have, or recently had, an extremely high mass-loss rate. Preliminary radiative transfer calculations suggest that the mass-loss rate may be as large as 10^{-3} Msun/yr. The 3.8 um ISAAC image shows a bipolar reflection nebula with a dark equatorial waist or torus, similar to that seen in optical Hubble Space Telescope (HST) images. The position angle of the nebula decreases significantly with increasing wavelength, suggesting that the dominant source of emission changes from scattering to thermal emission. Still, even up to 20 um the nebula is oriented approximately along the major axis of the nebula seen in the HST and ISAAC images, suggesting that the torus must be very cold, in agreement with the very red ISO spectrum. The 20 um image shows a roughly spherically symmetric extended halo, approximately 6'' in diameter, which is probably due to a previous phase of mass-loss on the AGB, suggesting a transition from a (more) spherically symmetric to a (more) axial symmetric form of mass-loss at the end of the AGB. We estimate the maximum dust particle sizes in the torus and in the reflection nebula to be 1.3 and 0.09 um respectively. The size of the particles in the torus is large compared to typical ISM values, but in agreement with high mass-loss rate objects like AFGL 4106 and HD161796. We discuss the possible reason for the difference in particle size between the torus and the reflection nebula., Comment: Accepted for publication by A&A

We present a model of the growth of water ice on silicate grains in the circumstellar envelopes of Asymptotic Giant Branch (AGB) stars and Red Super Giants. We consider the growth of ice by gas grain collisions, the thermal evaporation of ice from a grain, and sputtering. Our model contains several improvements compared to earlier models, including a detailed treatment of the effects of sputtering, a detailed calculation of the radiation pressure on the grain, and the treatment of subsonic drift velocities. In terms of drift velocity between the grains and gas in the envelope, we find that the ice formation process can be divided into three regimes: (i) a sputtering dominated regime where ice growth is heavily suppressed, (ii) an intermediate regime with moderately efficient condensation and (iii) a thermally dominated regime where ice formation is highly efficient. Sputtering is the critical factor which determines if ice formation can occur at all. We find that in Red Supergiants, ice formation is suppressed, while the winds of OH/IR stars allow for efficient condensation and can convert significant fractions of the available water vapor (tens of percent) into ice mantles on grains. Population II AGB stars hardly form ice due to their low dust to gas ratios. We also modify an analytical equation describing condensation and depletion (Jura and Morris, 1985) in order to give reasonable results for high and low drift velocities. Initially, ice will condense in crystalline form, but continuing condensation at low temperatures, and damage caused by interstellar UV photons favor the production of amorphous ice as well. We predict that a significant fraction of the ice formed will be amorphous., Comment: Accepted for publication by A&A

We have investigated the spatial distribution, and the properties and chemical composition of the dust orbiting HD 100546. This system is remarkably different from other isolated Herbig Ae/Be stars in both the strength of the mid-IR excess and the composition of the circumstellar dust. We speculate that the formation and spatial distribution of the crystalline dust observed in the HD 100546 system may be linked to the formation of a proto-Jupiter in the disk around HD 100546. Such a proto-Jupiter could gravitationally stir the disk leading to a collisional cascade of asteroidal sized objects producing small crystalline grains, or it could cause shocks by tidal interaction with the disk which might produce crystalline dust grains through flash heating. As shown by Malfait et al. (1998), the infrared spectrum of HD 100546 is very similar to that of C/1995 O1 Hale-Bopp (Crovisier et al. 1997). Using an identical methodology, we have therefore also studied this solar system comet. Both objects have an almost identical grain composition, but with the important difference that the individual dust species in Hale-Bopp are in thermal contact with each other, while this is not the case in HD 100546. This suggests that if similar processes leading to the dust composition as seen in HD 100546 also occurred in our own solar system, that Hale-Bopp formed after the formation of one or more proto-gas giants., Comment: 17 pages, 10 figures, Accepted for publication in Astronomy and Astrophysics

We present radiation-driven wind models for Luminous Blue Variables (LBVs) and predict their mass-loss rates. A comparison between our predictions and the observations of AG Car shows that the variable mass loss behaviour of LBVs is due the recombination/ionisation of Fe IV/III and Fe III/II. We also derive a present-day mass of 35 Msun for AG Car., Comment: 2 pages; in: A Massive Star Odyssey, from Main Sequence to Supernova

We show that the stellar masses implied by our predictions of the wind properties of massive stars are in agreement with masses derived from evolution theory and from direct measurements using spectroscopic binaries, contrary to previous attempts to derive masses from wind theory., Comment: 2 pages; in: A Massive Star Odyssey, from Main Sequence to Supernova

We present radiation-driven wind models for Luminous Blue Variables (LBVs) and predict their mass-loss rates. We study the effects of lower masses and modified abundances in comparison to the normal OB supergiants, and we find that the main difference in mass loss is due to the lower masses of LBVs. In addition, we find that the increase in helium abundance changes the mass-loss properties by small amounts (up to about 0.2 dex in log Mdot), while CNO processing is relatively unimportant for the mass-loss rate. A comparison between our mass loss predictions and the observations is performed for four relatively well-studied LBVs. The comparison shows that (i) the winds of LBVs are driven by radiation pressure on spectral lines, (ii) the variable mass loss behaviour of LBVs during their S Doradus-type variation cycles is explained by changes in the line driving efficiency, notably due to the recombination/ionisation of Fe IV/III and Fe III/II, and finally, (iii) the winds of LBVs can be used to derive their masses, as exemplified by the case of AG Car, for which we derive a present-day mass of 35 Msun., Comment: 12 pages; A&A accepted

We analyse the optical and IR spectra, as well as the spectral energy distribution (UV to mm) of the candidate Herbig Ae star HD100453. This star is particular, as it shows an energy distribution similar to that of other isolated Herbig Ae/Be stars (HAEBEs), but unlike most of them, it does not have a silicate emission feature at 10 um, as is shown in Meeus (2001). We confirm the HAEBE nature of HD100453 through an analysis of its optical spectrum and derived location in the H-R diagram. The IR spectrum of HD100453 is modelled by an optically thin radiative transfer code, from which we derive constraints on the composition, grain-size and temperature distribution of the circumstellar dust. We show that it is both possible to explain the lack of the silicate feature as (1) a grain-size effect - lack of small silicate grains, and (2) a temperature effect - lack of small, hot silicates, as proposed by Dullemond (2001), and discuss both possibilities., Comment: 9 pages, 7 figures; accepted by A&A

We have analysed the infrared hydrogen recombination lines of a sample of well studied hot massive stars observed with the Infrared Space Observatory. Our sample contains stars from several classes of objects, whose circumstellar environment is believed to be dominated by an ionized stellar wind (the Luminous Blue Variables) or by a dense disk-like geometry (Be stars and B[e] stars). We show that hydrogen infrared recombination lines can be used as a diagnostic tool to constrain the geometry of the ionized circumstellar material. The line strengths are sensitive to the density of the emitting gas. High densities result in optically thick lines for which line strengths are only dependent on the emitting surface. Low density gas produces optically thin lines which may be characterized by Menzel case B recombination. The ISO observations show that stellar winds are dominated by optically thin \hi\ recombination lines, while disks are dominated by optically thick lines. Disks and winds are well separated in a diagnostic diagram using the \hu/\bra and the \hu/\pfg line flux ratios. This diagnostic tool is useful to constrain the nature of hot star environments in case they are highly obscured, for instance while they are still embedded in their natal molecular cloud., Comment: 4 pages, 2 figures

We present an atlas of spectra of O- and B-type stars, obtained with the Short Wavelength Spectrometer (SWS) during the Post-Helium program of the Infrared Space Observatory (ISO). This program is aimed at extending the Morgan & Keenan classification scheme into the near-infrared. Later type stars will be discussed in a seperate publication. The observations consist of 57 SWS Post-Helium spectra from 2.4 to 4.1 microns, supplemented with 10 spectra acquired during the nominal mission with a similar observational setting. For B-type stars, this sample provides ample spectral converage in terms of subtype and luminosity class. For O-type stars,the ISO sample is coarse and therefore is complemented with 8 UKIRT L'-band observations. In terms of the presence of diagnostic lines, the L'-band is likely the most promising of the near-infrared atmospheric windows for the study of the physical properties of B stars. Specifically, this wavelength interval contains the Brackett alpha, Pfund gamma, and other Pfund lines which are probes of spectral type, luminosity class and mass loss. Here, we present simple empirical methods based on the lines present in the 2.4 to 4.1 microns interval that allow the determination of: the spectral type of B dwarfs and giants to within two subtypes; the luminosity class of B stars to within two classes; the mass-loss rate of O stars and B supergiants to within 0.25 dex., Comment: 19 pages, 11 Postscript figures, accepted by A&A

We present a fit to the spectral energy distribution of OH 127.8+0.0, a typical asymptotic giant branch star with an optically thick circumstellar dust shell. The fit to the dust spectrum is achieved using non-spherical grains consisting of metallic iron, amorphous and crystalline silicates and water ice. Previous similar attempts have not resulted in a satisfactory fit to the observed spectral energy distributions, mainly because of an apparent lack of opacity in the 3--8 micron region of the spectrum. Non-spherical metallic iron grains provide an identification for the missing source of opacity in the near-infrared. Using the derived dust composition, we have calculated spectra for a range of mass-loss rates in order to perform a consistency check by comparison with other evolved stars. The L-[12 micron] colours of these models correctly predict the mass-loss rate of a sample of AGB stars, strengthening our conclusion that the metallic iron grains dominate the near-infrared flux. We discuss a formation mechanism for non-spherical metallic iron grains., Comment: 10 pages, 6 figures, accepted for publication by A&A

We report the discovery of a double ring structure in the waist of the nebula surrounding eta Carinae. The rings are detected in the mid-IR dust continuum at wavelengths of 7.9, 11.9, 12.9 and 20 micronmeter. The dust in the rings has a temperature of about 300 K. The orientation of the rings is inclined with respect to the axis of the homunculus by either 37 or 58 degrees. The central star is not in the projected centre of the structure defined by the two rings. This geometry is reminiscent of that seen in SN1987A and some planetary nebulae. We discuss several possible origins for this remarkable geometry and its orientation., Comment: Accepted for publication in A&A Letters