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

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

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

On the binary nature of massive blue hypergiants: high-resolution X-ray spectroscopy suggests that Cyg OB2 12 is a colliding wind binary

The blue hypergiant Cyg OB2-12 (B3Ia+) is a representative member of the class of very massive stars in a poorly understood evolutionary stage. We obtained its high-resolution X-ray spectrum using Chandra observatory. PoWR model atmospheres were calculated to provide realistic wind opacities and to establish the wind density structure. We find that collisional de-excitation is the dominant mechanism de-populating the metastable upper levels of the forbidden lines of the He-like ions SiXIV and MgXII. Comparison between the model and observations reveals that X-ray emission is produced in a dense plasma, which could reside only at the photosphere or in a colliding wind zone between binary components. The observed X-ray spectra are well fitted by thermal plasma models, with average temperatures in excess of 10 MK. The wind speed in Cyg OB2-12 is not high enough to power such high temperatures, but the collision of two winds in a binary system can be sufficient. We used archival data to investigate the X-ray properties of other blue hypergiants. In general, stars of this class are not detected as X-rays sources. We suggest that our new Chandra observations of Cyg OB2-12 can be best explained if Cyg OB2-12 is a colliding wind binary possessing a late O-type companion. This makes Cyg OB2-12 only the second binary system among the 16 known Galactic hypergiants. This low binary fraction indicates that the blue hypergiants are likely products of massive binary evolution during which they either accreted a significant amount of mass or already merged with their companion.Comment: accepted to Ap

The Galactic WC and WO stars: The impact of revised distances from Gaia DR2 and their role as massive black hole progenitors

Wolf-Rayet stars of the carbon sequence (WC stars) are an important cornerstone in the late evolution of massive stars before their core collapse. As core-helium burning, hydrogen-free objects with huge mass-loss, they are likely the last observable stage before collapse and thus promising progenitor candidates for type Ib/c supernovae. Their strong mass-loss furthermore provides challenges and constraints to the theory of radiatively driven winds. Thus, the determination of the WC star parameters is of major importance for several astrophysical fields. With Gaia DR2, for the first time parallaxes for a large sample of Galactic WC stars are available, removing major uncertainties inherent to earlier studies. In this work, we re-examine the sample from Sander et al. (2012) to derive key properties of the Galactic WC population. All quantities depending on the distance are updated, while the underlying spectral analyses remain untouched. Contrasting earlier assumptions, our study yields that WC stars of the same subtype can significantly vary in absolute magnitude. With Gaia DR2, the picture of the Galactic WC population becomes more complex: We obtain luminosities ranging from log L = 4.9 to 6.0 with one outlier having log L = 4.7. This indicates that the WC stars are likely formed from a broader initial mass range than previously assumed. We obtain mass-loss rates ranging between log Mdot = -5.1 and -4.1, with Mdot propto L^0.68 and a linear scaling of the modified wind momentum with luminosity. We discuss the implications for stellar evolution, including unsolved issues regarding the need of envelope inflation to address the WR radius problem, and the open questions in regard to the connection of WR stars with Gamma-ray bursts. WC and WO stars are progenitors of massive black holes, collapsing either silently or in a supernova that most-likely has to be preceded by a WO stage.Comment: 19 pages, 13 figures, 6 tables; A&A, v2: version in pres

The stellar and wind parameters of six prototypical HMXBs and their evolutionary status

Context. High-mass X-ray binaries (HMXBs) are exceptional astrophysical laboratories that offer a rare glimpse into the physical processes that govern accretion on compact objects, massive-star winds, and stellar evolution. In a subset of the HMXBs, the compact objects accrete matter solely from winds of massive donor stars. These so-called wind-fed HMXBs are divided in persistent (classical) HMXBs and supergiant fast X-ray transients (SFXTs) according to their X-ray properties. While it has been suggested that this dichotomy depends on the characteristics of stellar winds, they have been poorly studied. Aims. With this investigation, we aim to remedy this situation by systematically analyzing donor stars of wind-fed HMXBs that are observable in the UV, concentrating on those with neutron star (NS) companions. Methods. We obtained Swift X-ray data, HST UV spectra, and additional optical spectra for all our targets. The spectral analysis of our program stars was carried out with the Potsdam Wolf–Rayet model atmosphere code. Results. Our multi-wavelength approach allows us to provide stellar and wind parameters for six donor stars (four wind-fed systems and two OBe X-ray binaries). The wind properties are in line with the predictions of the line-driven wind theory. Based on the abundances, three of the donor stars are in an advanced evolutionary stage, while for some of the stars, the abundance pattern indicates that processed material might have been accreted. When passing by the NS in its tight orbit, the donor star wind has not yet reached its terminal velocity but it is still significantly slower; its speed is comparable with the orbital velocity of the NS companion. There are no systematic differences between the two types of wind-fed HMXBs (persistent versus transients) with respect to the donor stars. For the SFXTs in our sample, the orbital eccentricity is decisive for their transient X-ray nature. The dichotomy of wind-fed HMXBs studied in this work is primarily a result of the orbital configuration, while in general it is likely that it reflects a complex interplay between the donor-star parameters, the orbital configuration, and the NS properties. Based on the orbital parameters and the further evolution of the donor stars, the investigated HMXBs will presumably form Thorne–Żytkow objects in the future.The first author of this work (R.H.) is supported by the Deutsche Forschungsgemeinschaft (DFG) under grant HA 1455/28-1. L.M.O. acknowledges support from the Verbundforschung grant 50 OR 1809. J.M.T. acknowledges the research grant ESP2017-85691-P. A.A.C.S. is supported by the Deutsche Forschungsgemeinschaft (DFG) under grant HA 1455/26. F.F., K.S., and A.B. are grateful for support from STScI Grant HST-GO-13703.002-A. T.S. acknowledges support from the European Research Council (ERC) under the European Union’s DLV-772225-MULTIPLES Horizon 2020 research and innovation programme

Phase-resolved spectroscopic analysis of the eclipsing black hole X-ray binary M33 X-7: System properties, accretion, and evolution

M33 X-7 is the only known eclipsing black hole high mass X-ray binary. The system is reported to contain a very massive O supergiant donor and a massive black hole in a short orbit. The high X-ray luminosity and its location in the metal-poor galaxy M33 make it a unique laboratory for studying the winds of metal-poor donor stars with black hole companions and it helps us to understand the potential progenitors of black hole mergers. Using phase-resolved simultaneous HST- and XMM-Newton-observations, we traced the interaction of the stellar wind with the black hole. We observed a strong Hatchett-McCray effect in M33 X-7 for the full range of wind velocities. Our comprehensive spectroscopic investigation of the donor star (X-ray+UV+optical) yields new stellar and wind parameters for the system that differ significantly from previous estimates. In particular, the masses of the components are considerably reduced to ≈38 M⊙ for the O-star donor and ≈11.4 M⊙ for the black hole. The O giant is overfilling its Roche lobe and shows surface He enrichment. The donor shows a densely clumped wind with a mass-loss rate that matches theoretical predictions. An extended ionization zone is even present during the eclipse due to scattered X-ray photons. The X-ray ionization zone extends close to the photosphere of the donor during inferior conjunction. We investigated the wind-driving contributions from different ions and the changes in the ionization structure due to X-ray illumination. Toward the black hole, the wind is strongly quenched due to strong X-ray illumination. For this system, the standard wind-fed accretion scenario alone cannot explain the observed X-ray luminosity, pointing toward an additional mass overflow, which is in line with our acceleration calculations. The X-ray photoionization creates an He II emission region around the system emitting ∼1047 ph s−1. We computed binary evolutionary tracks for the system using MESA. Currently, the system is transitioning toward an unstable mass transfer phase, possibly resulting in a common envelope of the black hole and the O-star donor. Since the mass ratio is q ≳ 3.3 and the period is short, the system is unlikely to survive the common envelope, but will rather merge.VR acknowledges support by the Deutsches Zentrum fsür Luft- und Raumfahrt (DLR) under grant 50 OR 1912. VR and AACS acknowledge support by the Deutsche Forschungsgemeinschaft (DFG – German Research Foundation) in the form of an Emmy Noether Research Group (grant number SA4064/1-1, PI Sander). DP acknowledges financial support by the Deutsches Zentrum für Luft und Raumfahrt (DLR) grant FKZ 50 OR 2005

High-resolution X-ray spectroscopy of supergiant HMXB 4U 1700−37 during the compact object eclipse

We present an analysis of the first observation of the iconic high-mass X-ray binary 4U 1700−37 with the Chandra High-Energy Transmission Gratings (HETGs) during an X-ray eclipse. The goal of the observation was to study the structure/physical conditions in the clumpy stellar wind through high-resolution spectroscopy. We find the following: (a) Emission-line brightness from K-shell transitions, corresponding to near-neutral species, directly correlates with continuum illumination. However, these lines do not greatly diminish during eclipse. This is readily explained if fluorescence K α emission comes from the bulk of the wind. (b) The highly ionized Fe xxv and Fe xxvi Ly α diminish during eclipse. Thus, they must be produced in the vicinity of the compact object where log ξ > 3. (c) To describe the emission-line spectrum, the sum of two self-consistent photoionization models with low ionization (log ξ ∼ −1) and high ionization (log ξ ∼ 2.4) is required. From their emission measures, the clump-to-interclump density ratio can be estimated to be nc/ni ∼ 300. To fit the complex He-like Si xiii profile, the plasma requires a broadening with vbulk ∼ 840 km s−1. Reproducing the observed r ≈ f line fluxes requires the addition of a third collisionally ionized plasma. (d) Emission-line widths appear unresolved at the HETG resolution with the exception of silicon. There is no clear radial segregation between (quasi-)neutral and ionized species, consistent with cold wind clumps interspersed in a hot rarefied interclump medium.This research has been funded under the project ESP2017-85691-P. The research leading to these results has received funding from the European Union’s Horizon 2020 Programme under the Activities for the High-Energy Astrophysiscs Domain (AHEAD) project (grant agreement no. 654215). Victoria Grinberg (VG) was supported through the Margarete von Wrangell fellowship by the ESF and the Ministry of Science, Research and the Arts of Baden-Württemberg. Work at LLNL was performed under the auspices of the U.S. Department of Energy under contract no. DE-AC52-07NA27344 and supported through National Aeronautics and Space Administration (NASA) grants to Lawrence Livermore National Laboratory (LLNL). Lida M. Oskinova acknowledges Deutsches Zentrum für Luft und Raumfahrt (DLR) grant FKZ 50 OR 1508 and partial support by the Russian Government Program of Competitive Growth of Kazan Federal University

On the Apparent Absence of Wolf-Rayet+Neutron Star Systems: The Curious Case of WR124

Among the different types of massive stars in advanced evolutionary stages is the enigmatic WN8h type. There are only a few Wolf-Rayet (WR) stars with this spectral type in our Galaxy. It has long been suggested that WN8h-type stars are the products of binary evolution that may harbor neutron stars (NS). One of the most intriguing WN8h stars is the runaway WR 124 surrounded by its magnificent nebula M1-67. We test the presence of an accreting NS companion in WR 124 using ∼100 ks long observations by the Chandra X-ray observatory. The hard X-ray emission from WR 124 with a luminosity of L ∼ 10 erg s is marginally detected. We use the non-local thermodynamic equilibrium stellar atmosphere code PoWR to estimate the WR wind opacity to the X-rays. The wind of a WN8-type star is effectively opaque for X-rays, hence the low X-ray luminosity of WR 124 does not rule out the presence of an embedded compact object. We suggest that, in general, high-opacity WR winds could prevent X-ray detections of embedded NS, and be an explanation for the apparent lack of WR+NS systems.© 2018. The American Astronomical Society. All rights reserved..The authors are indebted to the anonymous referee for a detailed revision of the manuscript and constructive suggestions. Support for this Letter was provided by the National Aeronautics and Space Administration through Chandra award No. G07-18014X issued by the Chandra X-ray Center, which is operated by the Smithsonian Astrophysical Observatory for and on behalf of the NASA contract NAS8-03060. J.A.T., M.A.G., and H.T. are funded by UNAM DGAPA PAPIIT project IA100318. L.M.O. acknowledges support by the DLR grant 50 OR 1508 and partial support by the Russian Government Program of Competitive Growth of Kazan Federal University. A.A.C.S. is supported by the Deutsche Forschungsgemeinschaft (DFG) under grant HA 1455/26 and would like to thank STFC for funding under grant No. ST/R000565/1. Y.-H.C. acknowledges support from the Ministry of Science and Technology of Taiwan. T.S. acknowledges support from the German Verbundforschung (DLR) grant 50 OR 1612. J.M.T. acknowledges support from ESP2017-85691-P

Evidence of Compton cooling during an X-ray flare supports a neutron star nature of the compact object in 4U1700−37

Based on new Chandra X-ray telescope data, we present empirical evidence of plasma Compton cooling during a flare in the non-pulsating massive X-ray binary 4U1700−37. This behaviour might be explained by quasi-spherical accretion on to a slowly rotating magnetized neutron star (NS). In quiescence, the NS in 4U1700−37 is surrounded by a hot radiatively cooling shell. Its presence is supported by the detection of mHz quasi-periodic oscillations likely produced by its convection cells. The high plasma temperature and the relatively low X-ray luminosity observed during the quiescence, point to a small emitting area ∼1 km, compatible with a hotspot on an NS surface. The sudden transition from a radiative to a significantly more efficient Compton cooling regime triggers an episode of enhanced accretion resulting in a flare. During the flare, the plasma temperature drops quickly. The predicted luminosity for such transitions, ∼3 × 1035 erg s−1, is very close to the luminosity of 4U1700−37 during quiescence. The transition may be caused by the accretion of a clump in the stellar wind of the donor star. Thus, a magnetized NS nature of the compact object is strongly favoured.This research has been supported by the grant ESP2014-53672-C3-3P. AB acknowledges support from STScI award 44A-1096046. JJRR acknowledges support from MECD fellowship PRX17/00114