Carbon line formation and spectroscopy in O-type stars

The determination of chemical abundances constitutes a fundamental requirement for obtaining a complete picture of a star. Particularly in massive stars, CNO abundances are of prime interest, due to the nuclear CNO-cycle and various mixing processes which bring these elements to the surface. We aim at enabling a reliable carbon spectroscopy for our unified NLTE atmosphere code FASTWIND. We develop a new carbon model atom including CII/III/IV/V, and discuss problems related to carbon spectroscopy in O-type stars. We describe different tests to examine the reliability of our implementation, and investigate which mechanisms influence the carbon ionization balance. By comparing with high-resolution spectra from six O-type stars, we check in how far observational constraints can be reproduced by our new carbon line synthesis. Carbon lines are even more sensitive to a variation of temperature, gravity, and mass-loss rate, than hydrogen/helium lines. We are able to reproduce most of the observed lines from our stellar sample, and to estimate those specific carbon abundances which bring the lines from different ions into agreement. For hot dwarfs and supergiants earlier than O7, X-rays from wind-embedded shocks can impact the synthesized line strengths, particularly for CIV, potentially affecting the abundance determination. We have demonstrated our capability to derive realistic carbon abundances by means of FASTWIND, using our recently developed model atom. We found that complex effects can have a strong influence on the carbon ionization balance in hot stars. For a further understanding, the UV range needs to be explored as well. By means of detailed nitrogen and oxygen model atoms available to use, we will be able to perform a complete CNO abundance analysis for larger samples of massive stars, and to provide constraints on corresponding evolutionary models and aspects.Comment: 22 pages, 16 figures, 6 table

Atmospheric NLTE-Models for the Spectroscopic Analysis of Blue Stars with Winds. III. X-ray emission from wind-embedded shocks

X-rays/EUV radiation emitted from wind-embedded shocks in hot, massive stars can affect the ionization balance in their outer atmospheres, and can be the mechanism responsible for the production of highly ionized species. To allow for these processes in the context of spectral analysis, we have implemented such emission into our unified, NLTE model atmosphere/spectrum synthesis code FASTWIND. The shock structure and corresponding emission is calculated as a function of user-supplied parameters. We account for a temperature and density stratification inside the post-shock cooling zones, calculated for radiative and adiabatic cooling in the inner and outer wind, respectively. The high-energy absorption of the cool wind is considered by adding important K-shell opacities, and corresponding Auger ionization rates have been included into the NLTE network. We tested and verified our implementation carefully against corresponding results from various alternative model atmosphere codes, and studied the effects from shock emission for important ions from He, C, N, O, Si, and P. Surprisingly, dielectronic recombination turned out to play an essential role for the ionization balance of OIV/OV around Teff = 45,000 K. Finally, we investigated the behavior of the mass absorption coefficient, kappa_nu(r), important in the context of X-ray line formation in massive star winds. In almost all considered cases, direct ionization is of major influence, and Auger ionization significantly affects only NVI and OVI. The approximation of a radially constant kappa_nu is justified for r > 1.2 Rstar and lambda < 18 A, and also for many models at longer wavelengths. To estimate the actual value of this quantity, however, the HeII opacities need to be calculated from detailed NLTE modeling, at least for wavelengths longer than 18 to 20 A, and information on the individual CNO abundances has to be present.Comment: accepted by A&

The Discordance of Mass-Loss Estimates for Galactic O-Type Stars

We have determined accurate values of the product of the mass-loss rate and the ion fraction of P^{4+}, Mdot q(P^{4+}), for a sample of 40 Galactic O-type stars by fitting stellar-wind profiles to observations of the P V resonance doublet obtained with FUSE, ORFEUS/BEFS, and Copernicus. When P^{4+} is the dominant ion in the wind, Mdot q(P^{4+}) approximates the mass-loss rate to within a factor of 2. Theory predicts that P^{4+} is the dominant ion in the winds of O7-O9.7 stars, though an empirical estimator suggests that the range from O4-O7 may be more appropriate. However, we find that the mass-loss rates obtained from P V wind profiles are systematically smaller than those obtained from fits to Halpha emission profiles or radio free-free emission by median factors of about 130 (if P^{4+} is dominant between O7 and O9.7) or about 20 (if P^{4+} is dominant between O4 and O7). These discordant measurements can be reconciled if the winds of O stars in the relevant temperature range are strongly clumped on small spatial scales. We use a simplified two-component model to investigate the volume filling factors of the denser regions. This clumping implies that mass-loss rates determined from "density squared" diagnostics have been systematically over-estimated by factors of 10 or more, at least for a subset of O stars. Reductions in the mass-loss rates of this size have important implications for the evolution of massive stars and quantitative estimates of the feedback that hot-star winds provide to their interstellar environments.Comment: 26 pages, 4 figures; accepted for publication in Ap

Atmospheric NLTE-Models for the Spectroscopic Analysis of Blue Stars with Winds. II. Line-Blanketed Models

We present new or improved methods for calculating NLTE, line-blanketed model atmospheres for hot stars with winds (spectral types A to O), with particular emphasis on a fast performance. These methods have been implemented into a previous, more simple version of the model atmosphere code FASTWIND (Santolaya-Rey et al.1997) and allow to spectroscopically analyze rather large samples of massive stars in a reasonable time-scale, using state-of-the-art physics. We describe our (partly approximate) approach to solve the equations of statistical equilibrium for those elements which are primarily responsible for line-blocking and blanketing, as well as an approximate treatment of the line-blocking itself, which is based on a simple statistical approach using suitable means for line opacities and emissivities. Furthermore, we comment on our implementation of a consistent temperature structure. In the second part, we concentrate on a detailed comparison with results from those two codes which have been used in alternative spectroscopical investigations, namely CMFGEN (Hillier & Miller 1998) and WM-Basic (Pauldrach et al. 2001). All three codes predict almost identical temperature structures and fluxes for lambda > 400 A, whereas at lower wavelengths a number of discrepancies are found. Optical H/He lines as synthesized by FASTWIND are compared with results from CMFGEN, obtaining a remarkable coincidence, except for the HeI singlets in the temperature range between 36,000 to 41,000 K for dwarfs and between 31,000 to 35,000 K for supergiants, where CMFGEN predicts much weaker lines. Consequences due to these discrepancies are discussed.Comment: 30 pages incl. 20 figures, accepted by A&

The SILCC project: III. Regulation of star formation and outflows by stellar winds and supernovae

We study the impact of stellar winds and supernovae on the multi-phase interstellar medium using three-dimensional hydrodynamical simulations carried out with FLASH. The selected galactic disc region has a size of (500 pc)$^2$ x $\pm$ 5 kpc and a gas surface density of 10 M$_{\odot}$/pc$^2$. The simulations include an external stellar potential and gas self-gravity, radiative cooling and diffuse heating, sink particles representing star clusters, stellar winds from these clusters which combine the winds from indi- vidual massive stars by following their evolution tracks, and subsequent supernova explosions. Dust and gas (self-)shielding is followed to compute the chemical state of the gas with a chemical network. We find that stellar winds can regulate star (cluster) formation. Since the winds suppress the accretion of fresh gas soon after the cluster has formed, they lead to clusters which have lower average masses (10$^2$ - 10$^{4.3}$ M$_{\odot}$) and form on shorter timescales (10$^{-3}$ - 10 Myr). In particular we find an anti-correlation of cluster mass and accretion time scale. Without winds the star clusters easily grow to larger masses for ~5 Myr until the first supernova explodes. Overall the most massive stars provide the most wind energy input, while objects beginning their evolution as B-type stars contribute most of the supernova energy input. A significant outflow from the disk (mass loading $\gtrsim$ 1 at 1 kpc) can be launched by thermal gas pressure if more than 50% of the volume near the disc mid-plane can be heated to T > 3x10$^5$ K. Stellar winds alone cannot create a hot volume-filling phase. The models which are in best agreement with observed star formation rates drive either no outflows or weak outflows.Comment: 23 pages; submitted to MNRA

NLTE wind models of hot subdwarf stars

We calculate NLTE models of stellar winds of hot compact stars (central stars of planetary nebulae and subdwarf stars). The studied range of subdwarf parameters is selected to cover a large part of these stars. The models predict the wind hydrodynamical structure and provide mass-loss rates for different abundances. Our models show that CNO elements are important drivers of subdwarf winds, especially for low-luminosity stars. We study the effect of X-rays and instabilities on these winds. Due to the line-driven wind instability, a significant part of the wind could be very hot.Comment: 7 pages, to appear in Astrophysics and Space Science. The final publication will be available at springerlink.com

STIS UV spectroscopy of early B supergiants in M31

We analyze STIS spectra in the 1150-1700 Angstrom wavelength range obtained for six early B supergiants in the neighboring galaxy M31. Because of their likely high (nearly solar) abundance, these stars were originally chosen to be directly comparable to their Galactic counterparts, and represent a much-needed addition to our current sample of B-type supergiants, in our efforts to study the dependence of the Wind Momentum-Luminosity Relationship on spectral type and metallicity. As a first step to determine wind momenta we fit the P-Cygni profiles of the resonance lines of N V, Si IV and C IV with standard methods, and derive terminal velocities for all of the STIS targets. From these lines we also derive ionic stellar wind column densities. Our results are compared with those obtained previously in Galactic supergiants, and confirm earlier claims of `normal' wind line intensities and terminal velocities in M31. For half of the sample we find evidence for an enhanced maximum turbulent velocity when compared to Galactic counterparts.Comment: 15 pages, 9 figures, 2 tables. Accepted for publication in The Astrophysical Journa

A Medium Resolution Near-Infrared Spectral Atlas of O and Early B Stars

We present intermediate resolution (R ~ 8,000 - 12,000) high signal-to-noise H- and K-band spectroscopy of a sample of 37 optically visible stars, ranging in spectral type from O3 to B3 and representing most luminosity classes. Spectra of this quality can be used to constrain the temperature, luminosity and general wind properties of OB stars, when used in conjunction with sophisticated atmospheric model codes. Most important is the need for moderately high resolutions (R > 5000) and very high signal-to-noise (S/N > 150) spectra for a meaningful profile analysis. When using near-infrared spectra for a classification system, moderately high signal-to-noise (S/N ~ 100) is still required, though the resolution can be relaxed to just a thousand or two. In the appendix we provide a set of very high quality near-infrared spectra of Brackett lines in six early-A dwarfs. These can be used to aid in the modeling and removal of such lines when early-A dwarfs are used for telluric spectroscopic standards.Comment: 12 pages, 3 tables, 14 figures. AASTex preprint style. To appear in ApJS, November 2005. All spectra are available by contacting M.M. Hanso

On the sensitivity of HeI singlet lines to the FeIV model atom in O stars

Recent calculations and analyses of O star spectra have revealed discrepancies between theory and observations, and between different theoretical calculations, for the strength of optical HeI singlet transitions.We investigate the source of these discrepancies. Using a non-LTE radiative transfer code we have undertaken detailed test calculations for a range of O star properties. Our principal test model has parameters similar to those of the O9V star, 10 Lac. We show that the discrepancies arise from uncertainties in the radiation field in the HeI resonance transition near 584Angs. The radiation field at 584Angs. is influenced by model assumptions, such as the treatment of line-blanketing and the adopted turbulent velocity, and by the FeIV atomic data. It isshown that two FeIV transitions near 584Angs can have a substantial influence on the strength of the HeI singlet transitions. Because of the difficulty of modeling the HeI singlet lines, particularly in stars with solar metalicity, the HeI triplet lines should be preferred in spectral analyses. These lines are much less sensitive to model assumptions.Comment: 7 pages, 9 figures, accepted for publication in A&

An explanation for the curious mass loss history of massive stars: from OB stars, through Luminous Blue Variables to Wolf-Rayet stars

The stellar winds of massive stars show large changes in mass-loss rates and terminal velocities during their evolution from O-star through the Luminous Blue Variable phase to the Wolf-Rayet phase. The luminosity remains approximately unchanged during these phases. These large changes in wind properties are explained in the context of the radiation driven wind theory, of which we consider four different models. They are due to the evolutionary changes in radius, gravity and surface composition and to the change from optically thin (in continuum) line driven winds to optically thick radiation driven winds.Comment: Accepted for publication in Astronomy and Astrophysics (Letter to the Editor