The Wolf-Rayet stars in M31: I. Analysis of the late-type WN stars

Context: Comprehensive studies of Wolf-Rayet stars were performed in the past for the Galactic and the LMC population. The results revealed significant differences, but also unexpected similarities between the WR populations of these different galaxies. Analyzing the WR stars in M31 will extend our understanding of these objects in different galactic environments. Aims: The present study aims at the late-type WN stars in M31. The stellar and wind parameters will tell about the formation of WR stars in other galaxies with different metallicity and star formation histories. The obtained parameters will provide constraints to the evolution of massive stars in the environment of M31. Methods: We used the latest version of the Potsdam Wolf-Rayet model atmosphere code to analyze the stars via fitting optical spectra and photometric data. To account for the relatively low temperatures of the late WN10 and WN11 subtypes, our WN models have been extended into this temperature regime. Results: Stellar and atmospheric parameters are derived for all known late-type WN stars in M31 with available spectra. All of these stars still have hydrogen in their outer envelopes, some of them up to 50% by mass. The stars are located on the cool side of the zero age main sequence in the Hertzsprung-Russell diagram, while their luminosities range from $10^5$ to $10^6$ Lsun. It is remarkable that no star exceeds $10^6$ Lsun. Conclusions: If formed via single-star evolution, the late-type WN stars in M31 stem from an initial mass range between 20 and 60 Msun. From the very late-type WN9-11 stars, only one star is located in the S Doradus instability strip. We do not find any late-type WN stars with the high luminosities known in the Milky Way.Comment: 11+11 pages, 13+18 figures, A&A, in pres

The impact of rotation on the line profiles of Wolf-Rayet stars

Massive Wolf-Rayet stars are recognized today to be in a very common, but short, evolutionary phase of massive stars. While our understanding of Wolf-Rayet stars has increased dramatically over the past decades, it remains unclear whether rapid rotators are among them. There are various indications that rapidly rotating Wolf-Rayet stars should exist. Unfortunately, due to their expanding atmospheres, rotational velocities of Wolf-Rayet stars are very difficult to measure. However, recently observed spectra of several Wolf-Rayet stars reveal peculiarly broad and round emission lines. Could these spectra imply rapid rotation? In this work, we model the effects of rotation on the atmospheres of Wolf-Rayet stars. We further investigate whether the peculiar spectra of five Wolf-Rayet stars may be explained with the help of stellar rotation, infer appropriate rotation parameters, and discuss the implications of our results. We make use of the Potsdam Wolf-Rayet (PoWR) non-LTE model atmosphere code. Since the observed spectra of Wolf-Rayet stars are mainly formed in their expanding atmospheres, rotation must be accounted for with a 3D integration scheme of the formal integral. For this purpose, we assume a rotational velocity field consisting of an inner co-rotating domain and an outer domain, where the angular momentum is conserved. We find that rotation can reproduce the unique spectra analyzed here. However, the inferred rotational velocities at the stellar surface are large (~200 km/s), and the inferred co-rotation radii (~10 stellar radii) suggest the existence of very strong photospheric magnetic fields (~20 kG)

Moving inhomogeneous envelopes of stars

Massive stars are extremely luminous and drive strong winds, blowing a large part of their matter into the galactic environment before they finally explode as a supernova. Quantitative knowledge of massive star feedback is required to understand our Universe as we see it. Traditionally, massive stars have been studied under the assumption that their winds are homogeneous and stationary, largely relying on the Sobolev approximation. However, observations with the newest instruments, together with progress in model calculations, ultimately dictate a cardinal change of this paradigm: stellar winds are highly inhomogeneous. Hence, we are now advancing to a new stage in our understanding of stellar winds. Using the foundations laid by V.V. Sobolev and his school, we now update and further develop the stellar spectral analysis techniques. New sophisticated 3-D models of radiation transfer in inhomogeneous expanding media elucidate the physics of stellar winds and improve classical empiric mass-loss rate diagnostics. Applications of these new techniques to multiwavelength observations of massive stars yield consistent and robust stellar wind parameters.Comment: slightly corrected version of the review for the special issue "V.V. Sobolev and his Legacy", Journal of Quantitative Spectroscopy and Radiative Transfe

Coupling hydrodynamics with comoving frame radiative transfer: II. Stellar wind stratification in the high-mass X-ray binary Vela X-1

CONTEXT: Vela X-1, a prototypical high mass X-ray binary (HMXB), hosts a neutron star (NS) in a close orbit around an early-B supergiant donor star. Accretion of the donor star's wind onto the NS powers its strong X-ray luminosity. To understand the physics of HMXBs, detailed knowledge about the donor star winds is required. AIMS: To gain a realistic picture of the donor star in Vela X-1, we constructed a hydrodynamically consistent atmosphere model describing the wind stratification while properly reproducing the observed donor spectrum. To investigate how X-ray illumination affects the stellar wind, we calculated additional models for different X-ray luminosity regimes. METHODS: We use the recently updated version of the PoWR code to consistently solve the hydrodynamic equation together with the statistical equations and the radiative transfer. RESULTS: The wind flow in Vela X-1 is driven by ions from various elements with Fe III and S III leading in the outer wind. The model-predicted mass-loss rate is in line with earlier empirical studies. The mass-loss rate is almost unaffected by the presence of the accreting NS in the wind. The terminal wind velocity is confirmed at $v_\infty \approx 600$ km/s. On the other hand, the wind velocity in the inner region where the NS is located is only $\approx 100$ km/s, which is not expected on the basis of a standard $\beta$-velocity law. In models with an enhanced level of X-rays, the velocity field in the outer wind can be altered. If the X-ray flux is too high, the acceleration breaks down because the ionization increases. CONCLUSIONS: Accounting for radiation hydrodynamics, our Vela X-1 donor atmosphere model reveals a low wind speed at the NS location, and it provides quantitative information on wind driving in this important HMXB.Comment: 19 pages, 10 figures, accepted for publication in Astronomy & Astrophysic

Engineering failure analysis and design optimisation with HiP-HOPS

The scale and complexity of computer-based safety critical systems, like those used in the transport and manufacturing industries, pose significant challenges for failure analysis. Over the last decade, research has focused on automating this task. In one approach, predictive models of system failure are constructed from the topology of the system and local component failure models using a process of composition. An alternative approach employs model-checking of state automata to study the effects of failure and verify system safety properties. In this paper, we discuss these two approaches to failure analysis. We then focus on Hierarchically Performed Hazard Origin & Propagation Studies (HiP-HOPS) - one of the more advanced compositional approaches - and discuss its capabilities for automatic synthesis of fault trees, combinatorial Failure Modes and Effects Analyses, and reliability versus cost optimisation of systems via application of automatic model transformations. We summarise these contributions and demonstrate the application of HiP-HOPS on a simplified fuel oil system for a ship engine. In light of this example, we discuss strengths and limitations of the method in relation to other state-of-the-art techniques. In particular, because HiP-HOPS is deductive in nature, relating system failures back to their causes, it is less prone to combinatorial explosion and can more readily be iterated. For this reason, it enables exhaustive assessment of combinations of failures and design optimisation using computationally expensive meta-heuristics. (C) 2010 Elsevier Ltd. All rights reserved

On the consistent treatment of the quasi-hydrostatic layers in hot star atmospheres

Context. Spectroscopic analysis remains the most common method to derive masses of massive stars, the most fundamental stellar parameter. While binary orbits and stellar pulsations can provide much sharper constraints on the stellar mass, these methods are only rarely applicable to massive stars. Unfortunately, spectroscopic masses of massive stars heavily depend on the detailed physics of model atmospheres. Aims. We demonstrate the impact of a consistent treatment of the radiative pressure on inferred gravities and spectroscopic masses of massive stars. Specifically, we investigate the contribution of line and continuum transitions to the photospheric radiative pressure. We further explore the effect of model parameters, e.g., abundances, on the deduced spectroscopic mass. Lastly, we compare our results with the plane-parallel TLUSTY code, commonly used for the analysis of massive stars with photospheric spectra. Methods. We calculate a small set of O-star models with the Potsdam Wolf-Rayet (PoWR) code using different approaches for the quasi-hydrostatic part. These models allow us to quantify the effect of accounting for the radiative pressure consistently. We further use PoWR models to show how the Doppler widths of line profiles and abundances of elements such as iron affect the radiative pressure, and, as a consequence, the derived spectroscopic masses. Results. Our study implies that errors on the order of a factor of two in the inferred spectroscopic mass are to be expected when neglecting the contribution of line and continuum transitions to the radiative acceleration in the photosphere. Usage of implausible microturbulent velocities, or the neglect of important opacity sources such as Fe, may result in errors of approximately 50% in the spectroscopic mass. A comparison with TLUSTY model atmospheres reveals a very good agreement with PoWR at the limit of low mass-loss rates.The first author of this work (A.S.) is supported by the Deutsche Forschungsgemeinschaft (DFG) under grant HA 1455/22. T.S. is grateful for financial support from the Leibniz Graduate School for Quantitative Spectroscopy in Astrophysics, a joint project of the Leibniz Institute for Astrophysics Potsdam (AIP) and the Institute of Physics and Astronomy of the University of Potsdam. A.S. would like to thank the Aspen Center for Physics and the NSF Grant #1066293 for hospitality during the invention and writing of this paper

A new type of X-ray pulsar

X-ray emission from stars much more massive than the Sun was discovered only 35 years ago. Such stars drive fast stellar winds where shocks can develop, and it is commonly assumed that the X-rays emerge from the shock-heated plasma. Many massive stars additionally pulsate. However, hitherto it was neither theoretically predicted nor observed that these pulsations would affect their X-ray emission. Here we report the discovery of pulsating X-rays from the massive B-type star Xi1 Canis Majoris. This star is a variable of beta Cephei type and has a strong magnetic field. Our observations with the XMM-Newton telescope reveal X-ray pulsations with the same period as the fundamental stellar pulsation. This discovery challenges our understanding of stellar winds from massive stars, their X-ray emission, and their magnetism.Comment: manuscript draft. The revised paper is published in Nature Communication