Search

One objective of the Polstar spectropolarimetry mission is to characterize the degree of nonconservative mass transfer that occurs at various stages of binary evolution, from the initial mass reversal to the late Algol phase. The proposed instrument combines spectroscopic and polarimetric capabilities, where the spectroscopy can resolve Doppler shifts in UV resonance lines with 10 km/s precision, and polarimetry can resolve linear polarization with 1e-3 precision or better. The spectroscopy will identify absorption by mass streams seen in projection against the stellar disk as a function of orbital phase, hot accretion spots, as well as scattering from extended splash structures, circumbinary disks, and other flows in and above/below the orbital plane (e.g. jets) that fail to be transferred conservatively. The polarimetry affects more the light coming from material not seen against the stellar disk, allowing the geometry of the scattering to be tracked, resolving ambiguities left by the spectroscopy and light-curve information. For example, nonconservative mass streams ejected in the polar direction will produce polarization of the opposite sign from conservative transfer accreting in the orbital plane. Also, time domain coverage over a range of phases of the binary orbit are well supported by the Polstar observing strategy. Combining these elements will significantly improve our understanding of the mass transfer process and the amount of mass that can escape from the system, an important channel for changing the final mass, and ultimate supernova, of the large number of massive stars found in binaries at close enough separation to undergo interaction., Comment: Refereed paper in Ap&SS 367:Topical Collection, 23 pages, 10 figures

The most massive stars are thought to lose a significant fraction of their mass in a steady wind during the main-sequence and blue supergiant phases. This in turn sets the stage for their further evolution and eventual supernova, and preconditions the surrounding medium for all following events, with consequences for ISM energization, chemical enrichment, and dust formation. Understanding these processes requires accurate observational constraints on the mass-loss rates of the most luminous stars, which can also be used to test theories of stellar wind driving. In the past, mass-loss rates have been characterized via collisional emission processes such as optical Hα and free-free radio emission, but these so-called “density squared” diagnostics require correction in the presence of widespread clumping. Recent observational and theoretical evidence points to the likelihood of a ubiquitously high level of such clumping in hot-star winds, but quantifying its effects requires a deeper understanding of the complex dynamics of radiatively driven winds and their stochastic instabilities. Furthermore, large-scale structures initiating in surface anisotropies and propagating throughout the wind can also affect wind driving and alter mass-loss diagnostics. Time series spectroscopy of high resonance-line opacity in the UV, capable of high resolution and high signal-to-noise, are required to better understand these complex dynamics, and more accurately determine mass-loss rates. The proposed Polstar mission (Scowen et al. 2022, this volume) provides the necessary resolution at the Sobolev (∼10 km s−1) or sound-speed (∼20 km s−1) scale, for over three dozen bright galactic massive stars with signal-to noise an order of magnitude above that of the celebrated MEGA campaign (Massa et al. 1995) of the International Ultraviolet Explorer (IUE), via continuous observations that track propagating structures through the winds in real time. Supporting geometric constraints are provided by the polarimetric capabilities present in all the datasets of such a mission.

Upsilon Sagittarii is a hydrogen-deficient binary that has been suggested to be in its second stage of mass transfer, after the primary has expanded to become a helium supergiant following core helium exhaustion. A tentative identification of the faint companion in the ultraviolet led to mass estimates of both components that made the helium star in Upsilon Sagittarii a prototypical immediate progenitor of a type Ib/c supernova. However, no consistent model for the complex spectrum has been achieved, casting doubt on this interpretation. In the present study we provide for the first time a composite spectral model that fits the ultraviolet data, and clearly identifies the companion as a rapidly rotating slowly moving ~7Msun B-type star, unlike previously suggested. The stripped helium supergiant is less luminous than previous estimates, and, with an estimated mass smaller than 1Msun, is ruled out as a core-collapse supernova progenitor. We provide a detailed binary evolution scenario that explains the temperature and luminosity of the two components as well as the very low gravity (log g ~ 1.0 [cm/s^2]) and extreme hydrogen deficiency of the primary (atmospheric mass fraction XH ~ 0.001). The best-fitting model is a ~5Msun intermediate-mass primary with an initial orbital period of a few days, and a secondary that appears to have gained a significant amount of mass despite its high rotation. We conclude that Upsilon Sagittarii is a key system for testing binary evolution processes, especially envelope stripping and mass accretion., Accepted to MNRAS

We present the results of a four-month, spectroscopic campaign of the Wolf-Rayet dust-making binary, WR137. We detect only small-amplitude, random variability in the CIII5696 emission line and its integrated quantities (radial velocity, equivalent width, skewness, kurtosis) that can be explained by stochastic clumps in the wind of the WC star. We find no evidence of large-scale, periodic variations often associated with Corotating Interaction Regions that could have explained the observed intrinsic continuum polarization of this star. Our moderately high-resolution and high signal-to-noise average Keck spectrum shows narrow double-peak emission profiles in the Halpha, Hbeta, Hgamma, HeII6678 and HeII5876 lines. These peaks have a stable blue-to-red intensity ratio with a mean of 0.997 and a root-mean-square of 0.004, commensurate with the noise level; no variability is found during the entire observing period. We suggest that these profiles arise in a decretion disk around the O9 companion, which is thus an O9e star. The characteristics of the profiles are compatible with those of other Be/Oe stars. The presence of this disk can explain the constant component of the continuum polarization of this system, for which the angle is perpendicular to the plane of the orbit, implying that the rotation axis of the O9e star is aligned with that of the orbit. It remains to be explained why the disk is so stable within the strong ultraviolet radiation field of the O star. We present a binary evolutionary scenario that is compatible with the current stellar and system parameters., 12 pages, 5 Figures

Context: The star cluster R136 inside the LMC hosts a rich population of massive stars, including the most massive stars known. The strong stellar winds of these very luminous stars impact their evolution and the surrounding environment. We currently lack detailed knowledge of the wind structure that is needed to quantify this impact. Aims: To observationally constrain the stellar and wind properties of the massive stars in R136, in particular the parameters related to wind clumping. Methods: We simultaneously analyse optical and UV spectroscopy of 53 O-type and 3 WNh-stars using the FASTWIND model atmosphere code and a genetic algorithm. The models account for optically thick clumps and effects related to porosity and velocity-porosity, as well as a non-void interclump medium. Results: We obtain stellar parameters, surface abundances, mass-loss rates, terminal velocities and clumping characteristics and compare these to theoretical predictions and evolutionary models. The clumping properties include the density of the interclump medium and the velocity-porosity of the wind. For the first time, these characteristics are systematically measured for a wide range of effective temperatures and luminosities. Conclusions: We confirm a cluster age of 1.0-2.5 Myr and derive an initial stellar mass of $\geq 250 {\rm M}_\odot$ for the most massive star in our sample, R136a1. The winds of our sample stars are highly clumped, with an average clumping factor of $f_{\rm cl}=29\pm15$. We find tentative trends in the wind-structure parameters as a function of mass-loss rate, suggesting that the winds of stars with higher mass-loss rates are less clumped. We compare several theoretical predictions to the observed mass-loss rates and terminal velocities and find that none satisfactorily reproduces both quantities. The prescription of Krti\v{c}ka & Kub\'at (2018) matches best the observed mass-loss rates., Comment: Accepted for publication in A&A; Appendix I will not be included in the published version

The hunt for compact objects is on. Rarely seen massive binaries with a compact object are a crucial phase in the evolution towards compact object mergers. In Gaia data release 3 (DR3), the first Gaia astrometric orbital solutions for binary sources will become available. We investigate how many black holes (BH) with massive main-sequence dwarf companions (OB+BH binaries) are expected to be detected as binaries in DR3 and at the end of the nominal 5-yr mission (DR4). We estimate the fraction of identifiable OB+BH binaries and discuss the distributions of the masses of both components and the orbital periods. We study the impact of different BH-formation scenarios. Using tailored models for the massive star population, which assume a direct collapse and no kick upon BH formation (the fiducial case), we estimate the fraction of OB+BH systems that Gaia will detect as binaries. A distance distribution according to that of the second Alma Luminous Star catalogue (ALSII) is assumed. We investigate how many of the systems detected as binaries are identifiable as OB+BH binaries, using a method based on astrometric data. In the fiducial case we conservatively estimate that 77% of the OB+BH binaries in ALSII will be detected as binaries in DR3, of which 89% are identifiable as OB+BH binaries. This leads to a total of around 190 OB+BH binaries, a 20-fold increase in the known sample of OB+BH binaries, covering an uncharted parameter space of long-period binaries. The size and properties of the identifiable OB+BH population will contain crucial observational constraints to improve our understanding of BH formation. In DR4, the detected fraction will increase to 85%, of which 82% will be identifiable. Hence, an additional ~5 systems could be identified, which are expected to have either very short or long periods. The fractions become smaller for different BH-formation scenarios. (truncated), Comment: 22 pages, 18 figures, 2 tables, accepted for publication in A&A

Classical Wolf-Rayet (WR) stars mark an important stage in the late evolution of massive stars. As hydrogen-poor massive stars, these objects have lost their outer layers, while still losing further mass through strong winds indicated by their prominent emission line spectra. Wolf-Rayet stars have been detected in a variety of different galaxies. Their strong winds are a major ingredient of stellar evolution and population synthesis models. Yet, a coherent theoretical picture of their strong mass-loss is only starting to emerge. In particular, the occurrence of WR stars as a function of metallicity (Z) is still far from being understood. To uncover the nature of the complex and dense winds of Wolf-Rayet stars, we employ a new generation of model atmospheres including a consistent solution of the wind hydrodynamics in an expanding non-LTE situation. With this technique, we can dissect the ingredients driving the wind and predict the resulting mass-loss for hydrogen-depleted massive stars. Our modelling efforts reveal a complex picture with strong, non-linear dependencies on the luminosity-to-mass ratio and Z with a steep, but not totally abrupt onset for WR-type winds in helium stars. With our findings, we provide a theoretical motivation for a population of helium stars at low Z, which cannot be detected via WR-type spectral features. Our study of massive He-star atmosphere models yields the very first mass-loss recipe derived from first principles in this regime. Implementing our first findings in stellar evolution models, we demonstrate how traditional approaches tend to overpredict WR-type mass loss in the young Universe., 6 pages, 3 figures, to appear in the proceedings of IAUS 366 "The Origin of Outflows in Evolved Stars"

Stellar-mass black holes are the final remnants of stars born with more than 15 solar masses. Billions are expected to reside in the Local Group, yet only a few are known, mostly detected through X-rays emitted as they accrete material from a companion star. Here, we report on VFTS 243: a massive X-ray-faint binary in the Large Magellanic Cloud. With an orbital period of 10.4 d, it comprises an O-type star of 25 solar masses and an unseen companion of at least nine solar masses. Our spectral analysis excludes a non-degenerate companion at a 5σ confidence level. The minimum companion mass implies that it is a black hole. No other X-ray-quiet black hole is unambiguously known outside our Galaxy. The (near-)circular orbit and kinematics of VFTS 243 imply that the collapse of the progenitor into a black hole was associated with little or no ejected material or black-hole kick. Identifying such unique binaries substantially impacts the predicted rates of gravitational-wave detections and properties of core-collapse supernovae across the cosmos., This research has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (T.S. and H.S., grant agreement 772225: MULTIPLES). T.S. acknowledges support from the European Union’s Horizon 2020 programme under Marie Skłodowska-Curie grant agreement 101024605. This work is based on observations collected at the ESO under programme IDs 182.D-0222, 090.D-0323 and 092.D-0136. L.M. thanks the ESA and the Belgian Federal Science Policy Office (BELSPO) for their support in the framework of the PRODEX programme. P. Marchant acknowledges support from FWO junior postdoctoral fellowship 12ZY520N. We thank J. Hillier for making the CMFGEN code available. C.H. acknowledges support from the KU Leuven Research Council (grant C16/17/007: MAESTRO). P.A.C. acknowledges support from UK Science and Technology Facilities Council research grant ST/V000853/1. V.H.-B. acknowledges the support of the Natural Sciences and Engineering Research Council of Canada (NSERC) through grant RGPIN-2020-05990. M. Gieles acknowledges support from the Ministry of Science and Innovation through a Europa Excelencia grant (EUR2020-112157). The PoWR code, as well as the associated post-processing and visualization tool WRplot, has been developed under the guidance of W.-R. Hamann with substantial contributions from L. Koesterke, G. Gräfener, A. Sander, T.S. and other co-workers and students. This work has received funding from the ERC under the European Union’s Horizon 2020 research and innovation programme (grant agreement 945806). It is also supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy EXC 2181/1-390900948 (the Heidelberg STRUCTURES Excellence Cluster). S.T. acknowledges support from the Netherlands Research Council NWO (grants VENI 639.041.645, VIDI 203.061). A.H. and D.J.L. acknowledge support by the Spanish MCI through grant PGC-2018-0913741-B-C22 and the Severo Ochoa Programme through CEX2019-000920-S. J.M.A. acknowledges support from the Spanish Government Ministerio de Ciencia, Innovación y Universidades through grant PGC2018-095-049-B-C22.

We present the first visual orbit for the nitrogen-rich Wolf-Rayet binary, WR 133 (WN5o + O9I) based on observations made with the CHARA Array and the MIRC-X combiner. This orbit represents the first visual orbit for a WN star and only the third Wolf-Rayet star with a visual orbit. The orbit has a period of 112.8 d, a moderate eccentricity of 0.36, and a separation of $a$= 0.79 mas on the sky. We combine the visual orbit with an SB2 orbit and Gaia parallax to find that the derived masses of the component stars are $M_{\rm WR}$ = $9.3\pm1.6 M_\odot$ and $M_{\rm O}$ = $22.6\pm 3.2 M_\odot$, with the large errors owing to the nearly face-on geometry of the system combined with errors in the spectroscopic parameters. We also derive an orbital parallax that is identical to the {\it Gaia}-determined distance. We present a preliminary spectral analysis and atmosphere models of the component stars, and find the mass-loss rate in agreement with polarization variability and our orbit. However, the derived masses are low compared to the spectral types and spectral model. Given the close binary nature, we suspect that WR 133 should have formed through binary interactions, and represents an ideal target for testing evolutionary models given its membership in the cluster NGC 6871., Accepted to ApJ Letters

Context. To quantitatively decode the information stored within an observed spectrum, detailed modelling of the physical state and accurate radiative transfer solution schemes are required. The accuracy of the model is then typically evaluated by comparing the calculated and observed spectra. In the analysis of stellar spectra, the numerical model often needs to account for binary companions and 3D structures in the stellar envelopes. The enigmatic binary (or multiple) system LB-1 constitutes a perfect example of such a complex multi-D problem. Thus far, the LB-1 system has been indirectly investigated by 1D stellar-atmosphere codes and by spectral disentangling techniques, yielding differing conclusions about the nature of the system (e.g., a B-star and black-hole binary with an accretion disc around the black hole or a stripped-star and Be-star binary system have been proposed). Aims. To improve our understanding of the LB-1 system, we directly modelled the phase-dependent Hα line profiles of this system. To this end, we developed and present a multi-purpose binary-object spectral-synthesis code in 3D (BOSS-3D). Methods. BOSS-3D calculates synthetic line profiles for a given state of the circumstellar material. The standard pz-geometry commonly used for single stars is extended by defining individual coordinate systems for each involved object and by accounting for the appropriate coordinate transformations. The code is then applied to the LB-1 system, considering two main hypotheses, a binary containing a stripped star and Be star, or a B star and a black hole with a disc. Results. Comparing these two scenarios, neither model can reproduce the detailed phase-dependent shape of the Hα line profiles. A satisfactory match with the observations, however, is obtained by invoking a disc around the primary object in addition to the Be-star disc or the black-hole accretion disc. Conclusions. The developed code can be used to model synthetic line profiles for a wide variety of binary systems, ranging from transit spectra of planetary atmospheres, to post-asymptotic giant branch binaries including circumstellar and circumbinary discs and massive-star binaries with stellar winds and disc systems. For the LB-1 system, our modelling provides strong evidence that each object in the system contains a disc-like structure.

BAT99 126 is a multiple system in the Large Magellanic Cloud containing a Wolf-Rayet (WR) star, which has a reported spectroscopic (orbital) period of 25.5 days and a photometric (orbital) period of 1.55 days, and hence is potentially one of the shortest WR binaries known to date. Such short-period binary systems containing a WR star in low-metallicity environments are prime candidate progenitors of black-hole (BH) mergers. By thoroughly analysing the spectroscopic and photometric data, we aim to establish the true multiplicity of BAT99 126, characterise the orbit(s) of the system, measure the physical properties of its individual components, and determine the overall evolutionary status of the system. Using newly acquired high resolution spectra, we measured radial velocities and performed a spectral analysis of the individual components. We estimated the age of the system and derived an evolutionary scenario for the 1.55-day system. BAT99 126 comprises at least four components. The 1.55-day photometric signal originates in an eclipsing O+O binary of spectral types O4 V and O6.5 V, and masses of 36+/-5 MSun and 15+/-2 MSun, respectively. For the WR star, we derived a spectral type WN2.5-3 and reject the previously reported 25.5-d period. We find evidence in the composite spectra of a fourth component, a B1 V star, which shows radial velocity variation. The configurations of the B-type star, WR star, and possible additional undetected components remain unknown. We estimated the age of the system to be 4.2 Myr. The O+O binary likely went through a phase of conservative mass transfer and is currently a near-contact system. We show that BAT99 126 is a multiple - quadruple or higher-order - system with a total initial mass > 160 MSun. The 1.55-day O+O binary most likely will not evolve towards a BH+BH merger, but instead will merge before the components collapse to BHs., 14+2 pages, 12+1 figures, accepted in A&A

ispartof: Nature vol:580 issue:7805 pages:E11-E15 ispartof: location:England status: published

A high fraction of massive stars are found to be binaries but only a few of them are reported as photometrically variable. By studying the populations of SB2 in the 30 Doradus region, we found a subset of them that have photometry from the OGLE project and that display variations in their light curves related to orbital motions. The goal of this study is to determine the dynamical masses and radii of the 26 binary components to investigate the mass-discrepancy problem and to provide an empirical mass-luminosity relation for the LMC. We use the PHOEBE programme to perform a systematic analysis of the OGLE V and I light curves obtained for 13 binary systems in 30 Dor. We adopt Teff, and orbital parameters derived previously to obtain the inclinations of the systems and the parameters of the individual components. Three systems display eclipses in their light curves, while the others only display ellipsoidal variations. We classify two systems as over-contact, five as semi-detached, and four as detached. The two remaining systems have uncertain configurations due to large uncertainties on their inclinations. The fact that systems display ellipsoidal variations has a significant impact on the inclination errors. From the dynamical masses, luminosities, and radii, we provide LMC-based empirical mass-luminosity and mass-radius relations, and we compare them to other relations given for the Galaxy, the LMC, and the SMC. These relations differ for different mass ranges, but do not seem to depend on the metallicity regimes. We also compare the dynamical, spectroscopic, and evolutionary masses of the stars in our sample. While the dynamical and spectroscopic masses agree with each other, the evolutionary masses are systematically higher, at least for stars in semi-detached systems. This suggests that the mass discrepancy can be partly explained by past or ongoing interactions between the stars., Comment: Accepted to A&A The core of the paper 9 pages, 5 pages of Appendix

HR 6819 was recently proposed to be a triple system consisting of an inner B-type giant + black hole binary with an orbital period of 40d and an outer Be tertiary. This interpretation is mainly based on two inferences: that the emission attributed to the outer Be star is stationary, and that the inner star, which is used as mass calibrator for the black hole, is a B-type giant. We re-investigate the properties of HR 6819 by spectral disentangling and an atmosphere analysis of the disentangled spectra to search for a possibly simpler alternative explanation for HR 6819. Disentangling implies that the Be component is not a static tertiary, but rather a component of the binary in the 40-d orbit. The inferred radial velocity amplitudes imply an extreme mass ratio of M_2/M_1 = 15 +/- 3. We infer spectroscopic masses of 0.4$^{+0.3}_{-0.1}$ Msun and 6$^{+5}_{-3}$ Msun for the primary and secondary, which agree well with the dynamical masses for an inclination of i = 32 deg. This indicates that the primary might be a stripped star rather than a B-type giant. Evolutionary modelling suggests that a possible progenitor system would be a tight (P_i ~ 2d) B+B binary system that experienced conservative mass transfer. While the observed nitrogen enrichment of the primary conforms with the predictions of the evolutionary models, we find no indications for the predicted He enrichment. We suggest that HR 6819 is a binary system consisting of a stripped B-type primary and a rapidly-rotating Be star that formed from a previous mass-transfer event. In the framework of this interpretation, HR 6819 does not contain a black hole. Interferometry can distinguish between these two scenarios by providing an independent measurement of the separation between the visible components., Comment: Submitted to A&A, 13 pages (16 figures and 2 tables); 4 pages supplementary material (4 figures and 4 tables). Comments are welcome

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 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. 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. We obtained Swift X-ray data, HST UV spectra, and additional optical spectra for all our targets. 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. 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. Based on the orbital parameters and the further evolution of the donor stars, the investigated HMXBs will presumably form Thorne-\.Zytkow objects in the future., Comment: 44 pages, 13 figures, 15 tables; accepted for publication in A&A

Classical Wolf-Rayet (WR) stars are massive, hydrogen depleted, post main-sequence stars that exhibit emission-line dominated spectra. For a given metallicity Z, stars exceeding a certain initial mass M_single(Z) can reach the WR phase through intrinsic mass-loss (single-star channel). Stars of lower masses can reach the WR phase via binary mass transfer (binary channel). It is commonly assumed that the binary channel dominates the formation of WR stars in environments with low Z such as the SMC and LMC. However, their reported WR binary fractions of 30-40% are comparable to that of the Galaxy. Here, we explain this apparent contradiction by considering the minimum initial mass M_spec(Z) needed for the stripped product to appear as a WR star. We calibrate M_spec(Z) using the lowest-luminosity WR stars in the Clouds and the Galaxy. A range of M_single(Z) values are explored using various evolution codes. We estimate the additional contribution of the binary channel by considering the interval [M_spec(Z), M_single(Z)], which characterises the initial-mass range in which binaries can form additional WR stars. Results: The WR-phenomenon ceases below luminosities of logL = 4.9, 5.25, and 5.6 [Lsun] in the Galaxy, LMC, and SMC, which translates to He-star masses of 7.5, 11, 17 Msun and initial masses of M_spec = 18, 23, 37 Msun. Stripped stars with lower initial masses in the respective galaxies would tend to not appear as WR stars. M_single lies in the range 20-30, 30-60, and > 40 Msun for the Galaxy, LMC, and SMC. We find that that the additional contribution of the binary channel is a non-trivial function of Z that cannot be conclusively claimed to be monotonically increasing with decreasing Z. Hence, one should not a-priori expect that binary interactions become increasingly important in forming WR stars at low Z, or that the WR binary fraction grows with decreasing Z., 10 pages, 5 figures, accepted to A&A on the 12.01.2020. Comments welcome!

The stellar winds of hot stars have an important impact on both stellar and galactic evolution, yet their structure and internal processes are not fully understood in detail. One of the best nearby laboratories for studying such massive stellar winds is the O4I(n)fp star ζ Pup. After briefly discussing existing X-ray observations from Chandra and XMM, we present a simulation of X-ray emission line profile measurements for the upcoming 840 kilosecond Chandra HETGS observation. This simulation indicates that the increased S/N of this new observation will allow several major steps forward in the understanding of massive stellar winds. By measuring X-ray emission line strengths and profiles, we should be able to differentiate between various stellar wind models and map the entire wind structure in temperature and density. This legacy X-ray spectrum of ζ Pup will be a useful benchmark for future X-ray missions.

Massive binaries (MBs) play a crucial role in the Universe. Knowing the distributions of their orbital parameters (OPs) is important for a wide range of topics, from stellar feedback to binary evolution channels, from the distribution of supernova types to gravitational wave progenitors, yet, no direct measurements exist outside the Milky Way. The Tarantula Massive Binary Monitoring was designed to help fill this gap by obtaining multi-epoch radial velocity monitoring of 102 MBs in the 30 Dor. In this paper, we analyse 32 VLT/FLAMES observations of 93 O- and 7 B-type binaries. We performed a Fourier analysis and obtained orbital solutions for 82 systems: 51 single- and 31 double-lined spectroscopic binaries. Overall, the OPs and binary fraction are remarkably similar across the 30 Dor region and compared to existing Galactic samples (GSs). This indicates that within these domains environmental effects are of second order in shaping the properties of MBs. A small difference is found in the distribution of orbital periods (OrbPs), which is slightly flatter (in log space) in 30 Dor than in the Galaxy, although this may be compatible within error estimates and differences in the fitting methodology. Also, OrbPs in 30 Dor can be as short as 1.1 d; somewhat shorter than seen in GSs. Equal mass binaries q>0.95 in 30 Dor are all found outside NGC 2070 the very young and massive cluster at 30 Dor's core. One outstanding exception however is the fact that earliest spectral types (O2-O7) tend to have shorter OrbPs than latter (O9.2-O9.7). Our results point to a relative universality of the incidence rate of MBs and their OPs in the metallicity range from solar ($Z_{\odot}$) to about $0.5Z_{\odot}$. This provides the first direct constraints on MB properties in massive star-forming galaxies at the Universes peak of star formation at redshifts z~1 to 2, which are estimated to have $Z~0.5Z_{\odot}$., 36 pages, 19 figures, 10 tables, published in A\&A

The massive O4.5 V + O5.5 V binary VFTS 352 in the Tarantula nebula is one of the shortest-period and most massive overcontact binaries known. Recent theoretical studies indicate that some of these systems could ultimately lead to the formation of gravitational waves via black hole binary mergers through the chemically homogeneous evolution pathway. By analyzing ultraviolet-optical phase-resolved spectroscopic data, we aim to constrain atmospheric and wind properties that could be later used to confront theoretical predictions from binary evolution. In particular, surface abundances are powerful diagnostics of the evolutionary status, mass transfer and the internal mixing processes. From a set of 32 VLT/FLAMES visual and 8 HST/COS ultraviolet spectra, we used spectral disentangling to separate the primary and secondary components. Using a genetic algorithm wrapped around the NLTE model atmosphere and spectral synthesis code FASTWIND, we perform an 11-parameter optimization to derive the atmospheric and wind parameters of both components, including the surface abundances of He, C, N, O and Si. We find that both components are hotter than expected compared to single-star evolutionary models indicating that additional mixing processes may be at play. However the derived chemical abundances do not show significant indications of mixing when adopting baseline values typical for the system environment., 24 pages, 4 tables, 12 figures, accepted for publication in ApJ

Massive Wolf-Rayet (WR) stars dominate the radiative and mechanical energy budget of galaxies and probe a critical phase in the evolution of massive stars prior to core-collapse. It is not known whether core He-burning WR stars (classical WR, cWR) form predominantly through wind-stripping (w-WR) or binary stripping (b-WR). With spectroscopy of WR binaries so-far largely avoided due to its complexity, our study focuses on the 44 WR binaries / binary candidates of the Large Magellanic Cloud (LMC, metallicity Z~0.5 Zsun), identified on the basis of radial velocity variations, composite spectra, or high X-ray luminosities. Relying on a diverse spectroscopic database, we aim to derive the physical and orbital parameters of our targets, confronting evolution models of evolved massive stars at sub-solar metallicity, and constraining the impact of binary interaction in forming them. Spectroscopy is performed using the Potsdam Wolf-Rayet (PoWR) code and cross-correlation techniques. Disentanglement is performed using the code Spectangular or the shift-and-add algorithm. Evolutionary status is interpreted using the Binary Population and Spectral Synthesis (BPASS) code, exploring binary interaction and chemically-homogeneous evolution. No obvious dichotomy in the locations of apparently-single and binary WN stars on the Hertzsprung-Russell diagram is apparent. According to commonly used stellar evolution models (BPASS, Geneva), most apparently-single WN stars could not have formed as single stars, implying that they were stripped by an undetected companion. Otherwise, it must follow that pre-WR mass-loss/mixing (e.g., during the red supergiant phase) are strongly underestimated in standard stellar evolution models., Comment: accepted to A&A on 10.05.2019; 69 pages (25 main paper + 44 appendix); Corrigendum: Shenar et al. 2020, A&A, 641, 2: An unfortunate typo in the implementation of the "transformed radius" caused errors of up to ~0.5dex in the derived mass-loss rates. This has now been corrected

Comprehensive spectral analyses of the Galactic Wolf-Rayet stars of the nitrogen sequence (i.e.\ the WN subclass) have been performed in a previous paper. However, the distances of these objects were poorly known. Distances have a direct impact on the "absolute" parameters, such as luminosities and mass-loss rates. The recent Gaia Data Release (DR2) of trigonometric parallaxes includes nearly all WN stars of our Galactic sample. In the present paper, we apply the new distances to the previously analyzed Galactic WN stars and rescale the results accordingly. On this basis, we present a revised catalog of 55 Galactic WN stars with their stellar and wind parameters. The correlations between mass-loss rate and luminosity show a large scatter, for the hydrogen-free WN stars as well as for those with detectable hydrogen. The slopes of the $\log L - \log \dot{M}$ correlations are shallower than found previously. The empirical Hertzsprung-Russell diagram (HRD) still shows the previously established dichotomy between the hydrogen-free early WN subtypes that are located on the hot side of the zero-age main sequence (ZAMS), and the late WN subtypes, which show hydrogen and reside mostly at cooler temperatures than the ZAMS (with few exceptions). However, with the new distances, the distribution of stellar luminosities became more continuous than obtained previously. The hydrogen-showing stars of late WN subtype are still found to be typically more luminous than the hydrogen-free early subtypes, but there is a range of luminosities where both subclasses overlap. The empirical HRD of the Galactic single WN stars is compared with recent evolutionary tracks. Neither these single-star evolutionary models nor binary scenarios can provide a fully satisfactory explanation for the parameters of these objects and their location in the HRD., Comment: 11 pages; accepted for publication in Astronomy & Astrophysics

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.

Context. It is now well established that the majority of massive stars reside in multiple systems. However, the effect of multiplicity is not sufficiently understood, resulting in a plethora of uncertainties about the end stages of massive-star evolution. In order to investigate these uncertainties, it is useful to study massive stars just before their demise. Classical Wolf-Rayet stars represent the final end stages of stars at the upper-mass end. The multiplicity fraction of these stars was reported to be ∼0.4 in the Galaxy but no correction for observational biases has been attempted. Aims. The aim of this study is to conduct a homogeneous radial-velocity survey of a magnitude-limited (V ≤ 12) sample of Galactic Wolf-Rayet stars to derive their bias-corrected multiplicity properties. The present paper focuses on 12 northern Galactic carbon-rich (WC) Wolf-Rayet stars observable with the 1.2 m Mercator telescope on the island of La Palma. Methods. We homogeneously measured relative radial velocities (RVs) for carbon-rich Wolf-Rayet stars using cross-correlation. Variations in the derived RVs were used to flag binary candidates. We investigated probable orbital configurations and provide a first correction of observational biases through Monte-Carlo simulations. Results. Of the 12 northern Galactic WC stars in our sample, seven show peak-to-peak RV variations larger than 10 km s−1, which we adopt as our detection threshold. This results in an observed spectroscopic multiplicity fraction of 0.58 with a binomial error of 0.14. In our campaign, we find a clear lack of short-period (P

Accurate stellar parameters of individual objects in binary systems are essential to constrain the effects of binarity on stellar evolution. These parameters serve as a prerequisite to probing existing and future theoretical evolutionary models. We aim to derive the atmospheric parameters of the 31 SB2s in the TMBM sample. This sample, composed of detached, semi-detached and contact systems with at least one of the components classified as an O star, is an excellent test-bed to study how binarity can impact our knowledge of the evolution of massive stars. 32 epochs of FLAMES/GIRAFFE spectra are analysed using spectral disentangling to construct the individual spectra of 62 components. We apply the CMFGEN atmosphere code to determine their stellar parameters and their He, C and N surface abundances. From these properties, we show that the effects of tides on chemical mixing are limited. Components on longer-period orbits show higher nitrogen enrichment at their surface than those on shorter-period orbits, in contrast to expectations of rotational or tidal mixing, implying that other mechanisms play a role in this process. Components filling their Roche lobe are mass donors. They exhibit higher nitrogen content at their surface and rotate more slowly than their companions. By accreting new material, their companions spin faster and are rejuvenated. Their locations in the N-vsini diagram tend to show that binary products are good candidates to populate the two groups of stars (slowly rotating, nitrogen-enriched and rapidly rotating non-enriched) that cannot be reproduced through single-star population synthesis. This sample is the largest sample of binaries to be studied in such a homogeneous way. The study of these objects gives us strong observational constraints to test theoretical binary evolutionary tracks., Comment: accepted to A&A Core of the paper 15 pages, Appendix: 41 pages

Context. HD 93129A was classified as the earliest O-type star in the Galaxy (O2 If*) and is considered as the prototype of its spectral class. However, interferometry shows that this object is a binary system, while recent observations even suggest a triple configuration. None of the previous spectral analyses of this object accounted for its multiplicity. With new high-resolution UV and optical spectra, we have the possibility to reanalyze this key object, taking its binary nature into account for the first time. Aims. We aim to derive the fundamental parameters and the evolutionary status of HD 93129A, identifying the contributions of both components to the composite spectrum Methods. We analyzed UV and optical observations acquired with the Hubble Space Telescope and ESO’s Very Large Telescope. A multiwavelength analysis of the system was performed using the latest version of the Potsdam Wolf-Rayet model atmosphere code. Results. Despite the similar spectral types of the two components, we are able to find signatures from each of the components in the combined spectrum, which allows us to estimate the parameters of both stars. We derive log(L∕L⊙) = 6.15, Teff = 52 kK, and log Ṁ = −4.7 [M⊙ yr−1] for the primary Aa, and log(L∕L⊙) = 5.58, Teff = 45 kK, and log Ṁ = −5.8 [M⊙yr−1] for the secondary Ab. Conclusions. Even when accounting for the binary nature, the primary of HD 93129A is found to be one of the hottest and most luminous O stars in our Galaxy. Based on the theoretical decomposition of the spectra, we assign spectral types O2 If* and O3 III(f*) to components Aa and Ab, respectively. While we achieve a good fit for a wide spectral range, specific spectral features are not fully reproduced. The data are not sufficient to identify contributions from a hypothetical third component in the system.

Metal-poor massive stars are assumed to be progenitors of certain supernovae, gamma-ray bursts, and compact object mergers that might contribute to the early epochs of the Universe with their strong ionizing radiation. However, this assumption remains mainly theoretical because individual spectroscopic observations of such objects have rarely been carried out below the metallicity of the Small Magellanic Cloud. Aims. Here we explore the predictions of the state-of-the-art theories of stellar evolution combined with those of stellar atmospheres about a certain type of metal-poor (0.02 Z·) hot massive stars, the chemically homogeneously evolving stars that we call Transparent Wind Ultraviolet INtense (TWUIN) stars. Methods. We computed synthetic spectra corresponding to a broad range in masses (20-130 M·) and covering several evolutionary phases from the zero-age main-sequence up to the core helium-burning stage. We investigated the influence of mass loss and wind clumping on spectral appearance and classified the spectra according to the Morgan-Keenan (MK) system. Results. We find that TWUIN stars show almost no emission lines during most of their core hydrogen-burning lifetimes. Most metal lines are completely absent, including nitrogen. During their core helium-burning stage, lines switch to emission, and even some metal lines (oxygen and carbon, but still almost no nitrogen) are detected. Mass loss and clumping play a significant role in line formation in later evolutionary phases, particularly during core helium-burning. Most of our spectra are classified as an early-O type giant or supergiant, and we find Wolf-Rayet stars of type WO in the core helium-burning phase. Conclusions. An extremely hot, early-O type star observed in a low-metallicity galaxy could be the result of chemically homogeneous evolution and might therefore be the progenitor of a long-duration gamma-ray burst or a type Ic supernova. TWUIN stars may play an important role in reionizing the Universe because they are hot without showing prominent emission lines during most of their lifetime. © 2019 ESO., This research was supported by grant 16-01116S (GACR). The Astronomical Institute Ondrejov is supported by the project RVO: 67985815. D. Sz. is grateful for the relevant discussion with G. Grafener and N. Langer, as well as to A. Szabo, as usual. 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 number ST/R000565/1. TS acknowledges funding from the European Research Council (ERC) under the European Union's DLV-772225-MULTIPLES Horizon 2020 research and innovation program.

Context. SMC AB 6 is the shortest-period (P = 6.5 d) Wolf-Rayet (WR) binary in the Small Magellanic Cloud. This binary is therefore a key system in the study of binary interaction and formation of WR stars at low metallicity. The WR component in AB 6 was previously found to be very luminous (log L = 6.3 [L⊙]) compared to its reported orbital mass (≈8 M⊙), placing it significantly above the Eddington limit. Aims. Through spectroscopy and orbital analysis of newly acquired optical data taken with the Ultraviolet and Visual Echelle Spectrograph (UVES), we aim to understand the peculiar results reported for this system and explore its evolutionary history. Methods. We measured radial velocities via cross-correlation and performed a spectral analysis using the Potsdam Wolf-Rayet model atmosphere code. The evolution of the system was analyzed using the Binary Population and Spectral Synthesis evolution code. Results. AB 6 contains at least four stars. The 6.5 d period WR binary comprises the WR primary (WN3:h, star A) and a rather rapidly rotating (veq = 265 km s−1) early O-type companion (O5.5 V, star B). Static N III and N IV emission lines and absorption signatures in He lines suggest the presence of an early-type emission line star (O5.5 I(f), star C). Finally, narrow absorption lines portraying a long-term radial velocity variation show the existence of a fourth star (O7.5 V, star D). Star D appears to form a second 140 d period binary together with a fifth stellar member, which is a B-type dwarf or a black hole. It is not clear that these additional components are bound to the WR binary. We derive a mass ratio of MO∕MWR = 2.2 ± 0.1. The WR star is found to be less luminous than previously thought (log L = 5.9 [L⊙]) and, adopting MO = 41 M⊙ for star B, more massive (MWR = 18 M⊙). Correspondingly, the WR star does not exceed the Eddington limit. We derive the initial masses of Mi,WR = 60 M⊙ and Mi,O = 40 M⊙ and an age of 3.9 Myr for the system. The WR binary likely experienced nonconservative mass transfer in the past supported by the relatively rapid rotation of star B. Conclusions. Our study shows that AB 6 is a multiple – probably quintuple – system. This finding resolves the previously reported puzzle of the WR primary exceeding the Eddington limit and suggests that the WR star exchanged mass with its companion in the past.

Stochastically-triggered photospheric light variations reaching $\sim$$40$ mmag peak-to-valley amplitudes have been detected in the O8Iaf supergiant V973 Scorpii as the outcome of two months of high-precision time-resolved photometric observations with the BRIght Target Explorer (BRITE) nanosatellites. The amplitude spectrum of the time series photometry exhibits a pronounced broad bump in the low-frequency regime ($\lesssim$$0.9$ d$^{-1}$) where several prominent frequencies are detected. A time-frequency analysis of the observations reveals typical mode lifetimes of the order of $5-10$ days. The overall features of the observed brightness amplitude spectrum of V973 Sco match well with those extrapolated from two-dimensional hydrodynamical simulations of convectively-driven internal gravity waves randomly excited from deep in the convective cores of massive stars. An alternative or additional possible source of excitation from a subsurface convection zone needs to be explored in future theoretical investigations., 17 pages, 12 figures, 3 tables; Monthly Notices of the Royal Astronomical Society (MNRAS), in press

Clusters or associations of early-type stars are often associated with a 'superbubble' of hot gas. The formation of such superbubbles is caused by the feedback from massive stars. The complex N206 in the Large Magellanic Cloud exhibits a superbubble and a rich massive star population. We observed these massive stars using the FLAMES multi-object spectrograph at ESO-VLT. Available UV spectra from HST, IUE, and FUSE are also used. The spectral analysis is performed with Potsdam Wolf-Rayet (PoWR) model atmospheres. We present the stellar and wind parameters of the OB stars and the two WR binaries in the N206 complex. Twelve percent of the sample show Oe/Be type emission lines, although most of them appear to rotate far below critical. We found eight runaway stars based on their radial velocity. The wind-momentum luminosity relation of our OB sample is consistent with the expectations. The HRD of the OB stars reveals a large age spread (1-30 Myr), suggesting different episodes of star formation in the complex. The youngest stars are concentrated in the inner part of the complex, while the older OB stars are scattered over outer regions. We derived the present day mass function for the entire N206 complex as well as for the cluster NGC2018. Three very massive Of stars are found to dominate the feedback among 164 OB stars in the sample. The two WR winds alone release about as much mechanical luminosity as the whole OB star sample. The cumulative mechanical feedback from all massive stellar winds is comparable to the combined mechanical energy of the supernova explosions that likely occurred in the complex. Accounting also for the WR wind and supernovae, the mechanical input over the last five Myr is ~$2.3\times10^{52}$ erg, which exceeds the current energy content of the complex by more than a factor of five. The morphology of the complex suggests a leakage of hot gas from the superbubble., Comment: Accepted for publication in A&A

The first directly detected gravitational waves (GW 150914) were emitted by two coalescing black holes (BHs) with masses of ~36Msun and ~29Msun. Several scenarios have been proposed to put this detection into an astrophysical context. The evolution of an isolated massive binary system is among commonly considered models. Various groups have performed detailed binary-evolution calculations that lead to BH merger events. However, the question remains open as to whether binary systems with the predicted properties really exist. The aim of this paper is to help observers to close this gap by providing spectral characteristics of massive binary BH progenitors during a phase where at least one of the companions is still non-degenerate. Stellar evolution models predict fundamental stellar parameters. Using these as input for our stellar atmosphere code (PoWR), we compute a set of models for selected evolutionary stages of massive merging BH progenitors at different metallicities. The synthetic spectra obtained from our atmosphere calculations reveal that progenitors of massive BH merger events start their lives as O2-3V stars that evolve to early-type blue supergiants before they undergo core-collapse during the Wolf-Rayet phase. When the primary has collapsed, the remaining system will appear as a wind-fed high-mass X-ray binary. We provide feedback parameters, broad band magnitudes, and spectral templates that should help to identify such binaries in the future. Comparisons of empirically determined mass-loss rates with those assumed by evolution calculations reveal significant differences. The consideration of the empirical mass-loss rates in evolution calculations will possibly entail a shift of the maximum in the predicted binary-BH merger rate to higher metallicities, that is, more candidates should be expected in our cosmic neighborhood than previously assumed., Comment: 64 pages, 30 figures, accepted for publication in Astronomy & Astrophysics, v2: typos corrected

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 press

Massive stars are the key agents of feedback. Consequently, quantitative analysis of massive stars are required to understand how the feedback of these objects shapes/ creates the large scale structures of the ISM. The giant HII region N206 in the Large Magellanic Cloud contains an OB association that powers a X-ray superbubble, serving as an ideal laboratory in this context. We obtained optical spectra with the muti-object spectrograph FLAMES at the ESO-VLT. When possible, the optical spectroscopy was complemented by UV spectra from the HST, IUE, and FUSE archives. Detailed spectral classifications are presented for our sample Of-type stars. For the quantitative spectroscopic analysis we use the Potsdam Wolf-Rayet (PoWR) model atmosphere code. The physical parameters and nitrogen abundances of our sample stars are determined by fitting synthetic spectra to the observations. The stellar and wind parameters of nine Of-type stars are used to construct wind momentum,luminosity relationship. We find that our sample follows a relation close to the theoretical prediction, assuming clumped winds. The most massive star in the N206 association is an Of supergiant which has a very high mass-loss rate. Two objects in our sample reveal composite spectra, showing that the Of primaries have companions of late O subtype. All stars in our sample have an evolutionary age less than 4 million years, with the O2-type star being the youngest. All these stars show a systematic discrepancy between evolutionary and spectroscopic masses. All stars in our sample are nitrogen enriched. Nitrogen enrichment shows a clear correlation with increasing projected rotational velocities. The mechanical energy input from the Of stars alone is comparable to the energy stored in the N206 superbubble as measured from the observed X-ray and H alpha emission., Comment: Accepted for the pubblication in Astronomy & Astrophysics

The recent discovery of a gravitational wave from the merging of two black holes of about 30 solar masses each challenges our incomplete understanding of massive stars and their evolution. Critical ingredients comprise mass-loss, rotation, magnetic fields, internal mixing, and mass transfer in close binary systems. The imperfect knowledge of these factors implies large uncertainties for models of stellar populations and their feedback. In this contribution we summarize our empirical studies of Wolf-Rayet populations at different metallicities by means of modern non-LTE stellar atmosphere models, and confront these results with the predictions of stellar evolution models. At the metallicity of our Galaxy, stellar winds are probably too strong to leave remnant masses as high as 30 solar masses, but given the still poor agreement between evolutionary tracks and observation even this conclusion is debatable. At the low metallicity of the Small Magellanic Cloud, all WN stars which are (at least now) single are consistent with evolving quasi-homogeneously. O and B-type stars, in contrast, seem to comply with standard evolutionary models without strong internal mixing. Close binaries which avoided early merging could evolve quasi-homogeneously and lead to close compact remnants of relatively high masses that merge within a Hubble time., 5 pages, to be published in the Proceedings of the IAU Symposium No. 329 "The lives and death-throes of massive stars"

We report on the first multi-color precision light curve of the bright Wolf-Rayet binary $\gamma^2$ Velorum, obtained over six months with the nanosatellites in the BRITE- Constellation fleet. In parallel, we obtained 488 high-resolution optical spectra of the system. In this first report on the datasets, we revise the spectroscopic orbit and report on the bulk properties of the colliding winds. We find a dependence of both the light curve and excess emission properties that scales with the inverse of the binary separation. When analyzing the spectroscopic properties in combination with the photometry, we find that the phase dependence is caused only by excess emission in the lines, and not from a changing continuum. We also detect a narrow, high-velocity absorption component from the He I $\lambda$5876 transition, which appears twice in the orbit. We calculate smoothed-particle hydrodynamical simulations of the colliding winds and can accurately associate the absorption from He I to the leading and trailing arms of the wind shock cone passing tangentially through our line of sight. The simulations also explain the general strength and kinematics of the emission excess observed in wind lines such as C III $\lambda$5696 of the system. These results represent the first in a series of investigations into the winds and properties of $\gamma^2$ Velorum through multi-technique and multi-wavelength observational campaigns., Comment: 16 pages, 14 figures, additional measurements to be included in online dataset. Accepted for publication in MNRAS

Wolf-Rayet (WR) stars are the most advanced stage in the evolution of the most massive stars. The strong feedback provided by these objects and their subsequent supernova (SN) explosions are decisive for a variety of astrophysical topics such as the cosmic matter cycle. Consequently, understanding the properties of WR stars and their evolution is indispensable. A crucial but still not well known quantity determining the evolution of WR stars is their mass-loss rate. Since the mass loss is predicted to increase with metallicity, the feedback provided by these objects and their spectral appearance are expected to be a function of the metal content of their host galaxy. This has severe implications for the role of massive stars in general and the exploration of low metallicity environments in particular. Hitherto, the metallicity dependence of WR star winds was not well studied. In this contribution, we review the results from our comprehensive spectral analyses of WR stars in environments of different metallicities, ranging from slightly super-solar to SMC-like metallicities. Based on these studies, we derived empirical relations for the dependence of the WN mass-loss rates on the metallicity and iron abundance, respectively., Comment: 5 pages, 4 figures, to be published in the Proceedings of the IAU Symposium No. 329 "The lives and death-throes of massive stars"

Stars which start their lives with spectral types O and early-B are the progenitors of core-collapse supernovae, long gamma-ray bursts, neutron stars, and black holes. These massive stars are the primary sources of stellar feedback in star-forming galaxies. At low metallicities, the properties of massive stars and their evolution are not yet fully explored. Here we report a spectroscopic study of 320 OB stars in the Small Magellanic Cloud. The data, which we obtained with the ESO Very Large Telescope, were analyzed using state-of-the-art stellar atmosphere models. We find that stellar winds of our sample stars are much weaker than theoretically expected. The stellar rotation rates show a bi-modal distribution. The well-populated upper Hertzsprung-Russell diagram including our sample OB stars from SMC Wing as well as additional evolved stars all over SMC from the literature shows a strict luminosity limit. The comparison with single-star evolutionary tracks suggests a dichotomy in the fate of massive stars in the SMC. Only stars with Minit30M$_{\odot}$ appear to stay always hot and might evolve quasi chemically homogeneously, finally collapsing to relatively massive black holes. However, we find no indication that chemical mixing is correlated with rapid rotation. We report extended star-formation episodes in a quiescent low-density region of the Wing, which is progressing stochastically. We measure the key parameters of stellar feedback and establish the links between the rates of star formation and supernovae. Our study reveals that in metal-poor environments the stellar feedback is dominated by core-collapse supernovae in combination with winds and ionizing radiation supplied by a few of the most massive stars., Accepted for publication in Astronomy & Astrophysics

The study of massive stars in different metallicity environments is a central topic of current stellar research. The spectral analysis of massive stars requires adequate model atmospheres. The computation of such models is difficult and time-consuming. Therefore, spectral analyses are greatly facilitated if they can refer to existing grids of models. Here we provide grids of model atmospheres for OB-type stars at metallicities corresponding to the Small and Large Magellanic Clouds, as well as to solar metallicity. In total, the grids comprise 785 individual models. The models were calculated using the state-of-the-art Potsdam Wolf-Rayet (PoWR) model atmosphere code. The parameter domain of the grids was set up using stellar evolution tracks. For all these models, we provide normalized and flux-calibrated spectra, spectral energy distributions, feedback parameters such as ionizing photons, Zanstra temperatures, and photometric magnitudes. The atmospheric structures (the density and temperature stratification) are available as well. All these data are publicly accessible through the PoWR website., 12 pages, 14 figures, accepted for publication in Astronomy & Astrophysics

We present the first SB2 orbital solution and disentanglement of the massive Wolf-Rayet binary R145 (P = 159d) located in the Large Magellanic Cloud. The primary was claimed to have a stellar mass greater than 300Msun, making it a candidate for the most massive star known. While the primary is a known late type, H-rich Wolf-Rayet star (WN6h), the secondary could not be so far unambiguously detected. Using moderate resolution spectra, we are able to derive accurate radial velocities for both components. By performing simultaneous orbital and polarimetric analyses, we derive the complete set of orbital parameters, including the inclination. The spectra are disentangled and spectroscopically analyzed, and an analysis of the wind-wind collision zone is conducted. The disentangled spectra and our models are consistent with a WN6h type for the primary, and suggest that the secondary is an O3.5 If*/WN7 type star. We derive a high eccentricity of e = 0.78 and minimum masses of M1 sin^3 i ~ M2 sin^3 i ~ 13 +- 2 Msun, with q = M2 / M1 = 1.01 +- 0.07. An analysis of emission excess stemming from a wind-wind collision yields a similar inclination to that obtained from polarimetry (i = 39 +- 6deg). Our analysis thus implies M1 = 53^{+40}_{-20} and M2 = 54^{+40}_{-20} Msun, excluding M1 > 300Msun. A detailed comparison with evolution tracks calculated for single and binary stars, as well as the high eccentricity, suggest that the components of the system underwent quasi-homogeneous evolution and avoided mass-transfer. This scenario would suggest current masses of ~ 80 Msun and initial masses of Mi,1 ~ 105 and Mi,2 ~ 90Msun, consistent with the upper limits of our derived orbital masses, and would imply an age of ~2.2 Myr., Accepted for Publication in A&A, 16 pages, 17 figures and 4 tables

WR 148 (HD 197406) is an extreme runaway system considered to be a potential candidate for a short-period (4.3173 d) rare WR + compact object binary. Provided with new high-resolution, high signal-to-noise spectra from the Keck observatory, we determine the orbital parameters for both the primary WR and the secondary, yielding respective projected orbital velocity amplitudes of 88.1 ± 3.8 km s−1 and 79.2 ± 3.1 km s−1 and implying a mass ratio of 1.1 ± 0.1. We then apply the shift-and-add technique to disentangle the spectra and obtain spectra compatible with a WN7ha and an O4-6 star. Considering an orbital inclination of ∼67°, derived from previous polarimetry observations, the system's total mass would be a mere 2–3M⊙⁠, an unprecedented result for a putative massive binary system. However, a system comprising a 37M⊙ secondary (typical mass of an O5V star) and a 33M⊙ primary (given the mass ratio) would infer an inclination of ∼18°. We therefore reconsider the previous methods of deriving the orbital inclination based on time-dependent polarimetry and photometry. While the polarimetric results are inconclusive requiring better data, the photometric results favour low inclinations. Finally, we compute WR 148’s space velocity and retrace the runaway's trajectory back to the Galactic plane (GP). With an ejection velocity of 198 ± 27 km s−1 and a travel time of 4.7 ± 0.8 Myr to reach its current location, WR 148 was most likely ejected via dynamical interactions in a young cluster.

[Context]: Classical supergiant X-ray binaries (SGXBs) and supergiant fast X-ray transients (SFXTs) are two types of high-mass X-ray binaries (HMXBs) that present similar donors but, at the same time, show very different behavior in the X-rays. The reason for this dichotomy of wind-fed HMXBs is still a matter of debate. Among the several explanations that have been proposed, some of them invoke specific stellar wind properties of the donor stars. Only dedicated empiric analysis of the donors¿ stellar wind can provide the required information to accomplish an adequate test of these theories. However, such analyses are scarce. [Aims]: To close this gap, we perform a comparative analysis of the optical companion in two important systems: IGR J17544-2619 (SFXT) and Vela X-1 (SGXB). We analyze the spectra of each star in detail and derive their stellar and wind properties. As a next step, we compare the wind parameters, giving us an excellent chance of recognizing key differences between donor winds in SFXTs and SGXBs. [Methods]: We use archival infrared, optical and ultraviolet observations, and analyze them with the non-local thermodynamic equilibrium (NLTE) Potsdam Wolf-Rayet model atmosphere code. We derive the physical properties of the stars and their stellar winds, accounting for the influence of X-rays on the stellar winds. [Results]: We find that the stellar parameters derived from the analysis generally agree well with the spectral types of the two donors: O9I (IGR J17544-2619) and B0.5Iae (Vela X-1). The distance to the sources have been revised and also agree well with the estimations already available in the literature. In IGR J17544-2619 we are able to narrow the uncertainty to d = 3.0 ± 0.2 kpc. From the stellar radius of the donor and its X-ray behavior, the eccentricity of IGR J17544-2619 is constrained to e< 0.25. The derived chemical abundances point to certain mixing during the lifetime of the donors. An important difference between the stellar winds of the two stars is their terminal velocities (ν∞ = 1500 km s-1 in IGR J17544-2619 and ν∞ = 700 km s-1 in Vela X-1), which have important consequences on the X-ray luminosity of these sources. [Conclusions]: The donors of IGR J17544-2619 and Vela X-1 have similar spectral types as well as similar parameters that physically characterize them and their spectra. In addition, the orbital parameters of the systems are similar too, with a nearly circular orbit and short orbital period. However, they show moderate differences in their stellar wind velocity and the spin period of their neutron star which has a strong impact on the X-ray luminosity of the sources. This specific combination of wind speed and pulsar spin favors an accretion regime with a persistently high luminosity in Vela X-1, while it favors an inhibiting accretion mechanism in IGR J17544-2619. Our study demonstrates that the relative wind velocity is critical in class determination for the HMXBs hosting a supergiant donor, given that it may shift the accretion mechanism from direct accretion to propeller regimes when combined with other parameters., The work of AG-G was supported by the Spanish MICINN under FPI Fellowship BES-2011-050874 associated to the project AYA2010-15431. 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. This work has been partially supported by the Spanish Ministry of Economy and Competitiveness project numbers ESP2013-48637-C2-2P and ESP2014-53672-C3-3-P, the Generalitat Valenciana project number GV2014/088 and the Vicerectorat d’Investigació, Desenvolupament i Innovació de la Universitat d’Alacant under grant GRE12-35. S.M.N. thanks the support of the Spanish unemployment agency, allowing her to continue her scientific collaborations during the critical situation of the Spanish research system. The authors acknowledge the help of the International Space Science Institute at Bern, Switzerland, and the faculty of the European Space Astronomy Centre. A.S. is supported by the Deutsche Forschungsgemeinschaft (DFG) under grant HA 1455/26. Support for MAST for non-HST data is provided by the NASA office of space science via grant NAG5-7584 and by other grants and contracts.

We report on interferometric observations with the CHARA Array of two classical Wolf-Rayet stars in suspected binary systems, namely WR 137 and WR 138. In both cases, we resolve the component stars to be separated by a few milliarcseconds. The data were collected in the H-band, and provide a measure of the fractional flux for both stars in each system. We find that the WR star is the dominant H-band light source in both systems ($f_{\rm WR, 137} = 0.59\pm0.04$; $f_{\rm WR, 138} = 0.67\pm0.01$), which is confirmed through both comparisons with estimated fundamental parameters for WR stars and O dwarfs, as well as through spectral modeling of each system. Our spectral modeling also provides fundamental parameters for the stars and winds in these systems. The results on WR 138 provide evidence that it is a binary system which may have gone through a previous mass-transfer episode to create the WR star. The separation and position of the stars in the WR 137 system together with previous results from the IOTA interferometer provides evidence that the binary is seen nearly edge-on. The possible edge-on orbit of WR 137 aligns well with the dust production site imaged by the Hubble Space Telescope during a previous periastron passage, showing that the dust production may be concentrated in the orbital plane., 11 pages, 4 tables, 7 figures, accepted to MNRAS

During the summer of 2013, a 4-month spectroscopic campaign took place to observe the variabilities in three Wolf-Rayet stars. The spectroscopic data have been analyzed for WR 134 (WN6b), to better understand its behaviour and long-term periodicity, which we interpret as arising from corotating interaction regions (CIRs) in the wind. By analyzing the variability of the He II $\lambda$5411 emission line, the previously identified period was refined to P = 2.255 $\pm$ 0.008 (s.d.) days. The coherency time of the variability, which we associate with the lifetime of the CIRs in the wind, was deduced to be 40 $\pm$ 6 days, or $\sim$ 18 cycles, by cross-correlating the variability patterns as a function of time. When comparing the phased observational grayscale difference images with theoretical grayscales previously calculated from models including CIRs in an optically thin stellar wind, we find that two CIRs were likely present. A separation in longitude of $\Delta \phi \simeq$ 90$^{\circ}$ was determined between the two CIRs and we suggest that the different maximum velocities that they reach indicate that they emerge from different latitudes. We have also been able to detect observational signatures of the CIRs in other spectral lines (C IV $\lambda\lambda$5802,5812 and He I $\lambda$5876). Furthermore, a DAC was found to be present simultaneously with the CIR signatures detected in the He I $\lambda$5876 emission line which is consistent with the proposed geometry of the large-scale structures in the wind. Small-scale structures also show a presence in the wind, simultaneously with the larger scale structures, showing that they do in fact co-exist., Comment: 13 pages, 13 figures, 4 tables, will appear in the Monthly Notices for the Royal Astronomical Society, http://www.astro.umontreal.ca/~emily/CIR_Lifetime_WR134_full.pdf

Massive WR stars are evolved massive stars characterized by strong mass-loss. Hypothetically, they can form either as single stars or as mass donors in close binaries. About 40% of the known WR stars are confirmed binaries, raising the question as to the impact of binarity on the WR population. By performing a spectral analysis of all multiple WR systems in the SMC, we obtain the full set of stellar parameters for each individual component. Mass-luminosity relations are tested, and the importance of the binary evolution channel is assessed. The spectral analysis is performed with the PoWR model atmosphere code by superimposing model spectra that correspond to each component. Evolutionary channels are constrained using the BPASS evolution tool. Significant Hydrogen mass fractions (0.1 - 0.4) are detected in all WN components. A comparison with mass-luminosity relations and evolutionary tracks implies that the majority of the WR stars in our sample are not chemically homogeneous. The WR component in the binary AB 6 is found to be very luminous (Log L ~ 6.3 [Lsun]) given its orbital mass (~10 Msun), presumably because of observational contamination by a third component. Evolutionary paths derived for our objects suggest that Roche lobe overflow had occurred in most systems, affecting their evolution. However, the implied initial masses are large enough for the primaries to have entered the WR phase, regardless of binary interaction. Together with the results for the putatively single SMC WR stars, our study suggests that the binary evolution channel does not dominate the formation of WR stars at SMC metallicity., Accepted for publication in A&A on 30 March, 2016. 15 + 10 pages, 8 + 8 figures

Eclipsing systems of massive stars allow one to explore the properties of their components in great detail. We perform a multi-wavelength, non-LTE analysis of the three components of the massive multiple system $\delta$ Ori A, focusing on the fundamental stellar properties, stellar winds, and X-ray characteristics of the system. The primary's distance-independent parameters turn out to be characteristic for its spectral type (O9.5 II), but usage of the ${\rm \it Hipparcos}$ parallax yields surprisingly low values for the mass, radius, and luminosity. Consistent values follow only if $\delta$ Ori lies at about twice the ${\rm \it Hipparcos}$ distance, in the vicinity of the $\sigma$-Orionis cluster. The primary and tertiary dominate the spectrum and leave the secondary only marginally detectable. We estimate the V-band magnitude difference between primary and secondary to be $\Delta V \approx 2.\!\!^{\rm m}8$. The inferred parameters suggest the secondary is an early B-type dwarf ($\approx$ B1 V), while the tertiary is an early B-type subgiant ($\approx$ B0 IV). We find evidence for rapid turbulent velocities ($\sim 200$ km ${\rm s}^{-1}$) and wind inhomogeneities, partially optically thick, in the primary's wind. The bulk of the X-ray emission likely emerges from the primary's stellar wind ($\log L_{\text{X}} / L_{\text{Bol}} \approx -6.85$), initiating close to the stellar surface at $R_0 \sim 1.1\,R_*$. Accounting for clumping, the mass-loss rate of the primary is found to be $\log \dot{M} \approx -6.4\,[M_\odot\,{\rm yr}^{-1}]$, which agrees with hydrodynamic predictions, and provides a consistent picture along the X-ray, UV, optical and radio spectral domains., Comment: Accepted for Publication in APJ, 21 pages, 16 figures and 2 tables

Wolf-Rayet (WR) stars have a severe impact on their environments owing to their strong ionizing radiation fields and powerful stellar winds. Since these winds are considered to be driven by radiation pressure, it is theoretically expected that the degree of the wind mass-loss depends on the initial metallicity of WR stars. Following our comprehensive studies of WR stars in the Milky Way, M31, and the LMC, we derive stellar parameters and mass-loss rates for all seven putatively single WN stars known in the SMC. Based on these data, we discuss the impact of a low-metallicity environment on the mass loss and evolution of WR stars. The quantitative analysis of the WN stars is performed with the Potsdam Wolf-Rayet (PoWR) model atmosphere code. The physical properties of our program stars are obtained from fitting synthetic spectra to multi-band observations. In all SMC WN stars, a considerable surface hydrogen abundance is detectable. The majority of these objects have stellar temperatures exceeding 75 kK, while their luminosities range from 10^5.5 to 10^6.1 Lsun. The WN stars in the SMC exhibit on average lower mass-loss rates and weaker winds than their counterparts in the Milky Way, M31, and the LMC. By comparing the mass-loss rates derived for WN stars in different Local Group galaxies, we conclude that a clear dependence of the wind mass-loss on the initial metallicity is evident, supporting the current paradigm that WR winds are driven by radiation. A metallicity effect on the evolution of massive stars is obvious from the HRD positions of the SMC WN stars at high temperatures and high luminosities. Standard evolution tracks are not able to reproduce these parameters and the observed surface hydrogen abundances. Homogeneous evolution might provide a better explanation for their evolutionary past., Comment: 18+12 pages; 22+8 figures; accepted for publication in A&A

We present new grids of Potsdam Wolf-Rayet (PoWR ) model atmospheres for Wolf-Rayet stars of the nitrogen sequence (WN stars). The models have been calculated with the latest version of the PoWR stellar atmosphere code for spherical stellar winds. The WN model atmospheres include the non-LTE solutions of the statistical equations for complex model atoms, as well as the radiative transfer equation in the co-moving frame. Iron-line blanketing is treated with the help of the superlevel approach, while wind inhomogeneities are taken into account via optically thin clumps. Three of our model grids are appropriate for Galactic metallicity. The hydrogen mass fraction of these grids is 50%, 20%, and 0%, thus also covering the hydrogen-rich late-type WR stars that have been discovered in recent years. Three grids are adequate for LMC WN stars and have hydrogen fractions of 40%, 20%, and 0%. Recently, additional grids with SMC metallicity and with 60%, 40%, 20%, and 0% hydrogen have been added. We provide contour plots of the equivalent widths of spectral lines that are usually used for classification and diagnostics.

We report on both high-precision photometry from the Microvariability and Oscillations of Stars (MOST) space telescope and ground-based spectroscopy of the triple system δ Ori A, consisting of a binary O9.5II+early-B (Aa1 and Aa2) with P = 5.7 days, and a more distant tertiary (O9 IV years). This data was collected in concert with X-ray spectroscopy from the Chandra X-ray Observatory. Thanks to continuous coverage for three weeks, the MOST light curve reveals clear eclipses between Aa1 and Aa2 for the first time in non-phased data. From the spectroscopy, we have a well-constrained radial velocity (RV) curve of Aa1. While we are unable to recover RV variations of the secondary star, we are able to constrain several fundamental parameters of this system and determine an approximate mass of the primary using apsidal motion. We also detected second order modulations at 12 separate frequencies with spacings indicative of tidally influenced oscillations. These spacings have never been seen in a massive binary, making this system one of only a handful of such binaries that show evidence for tidally induced pulsations., La lista completa de autores que integran el documento puede consultarse en el archivo., Facultad de Ciencias Astronómicas y Geofísicas, Instituto de Astrofísica de La Plata

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., 13 pages, 10 figures, accepted for publication in A&A

With a deep Chandra/HETGS exposure of WR 6, we have resolved emission lines whose profiles show that the X-rays originate from a uniformly expanding spherical wind of high X-ray-continuum optical depth. The presence of strong helium-like forbidden lines places the source of X-ray emission at tens to hundreds of stellar radii from the photosphere. Variability was present in X-rays and simultaneous optical photometry, but neither were correlated with the known period of the system or with each other. An enhanced abundance of sodium revealed nuclear processed material, a quantity related to the evolutionary state of the star. The characterization of the extent and nature of the hot plasma in WR 6 will help to pave the way to a more fundamental theoretical understanding of the winds and evolution of massive stars., Comment: Accepted by the Astrophysical Journal