@article{AlmeidaSanaTayloretal.2017,
  author    = {Almeida, Leonardo A. and Sana, H. and Taylor, W. and Barb{\´a}, Rodolfo and Bonanos, Alceste Z. and Crowther, Paul and Damineli, Augusto and de Koter, A. and de Mink, Selma E. and Evans, C. J. and Gieles, Mark and Grin, Nathan J. and H{\´e}nault-Brunet, V. and Langer, Norbert and Lennon, D. and Lockwood, Sean and Ma{\´i}z Apell{\´a}niz, Jes{\´u}s and Moffat, A. F. J. and Neijssel, C. and Norman, C. and Ram{\´i}rez-Agudelo, O. H. and Richardson, N. D. and Schootemeijer, Abel and Shenar, Tomer and Soszyński, Igor and Tramper, Frank and Vink, J. S.},
  title     = {The tarantula massive binary monitoring},
  series = {Astronomy and astrophysics : an international weekly journal},
  volume    = {598},
  journal   = {Astronomy and astrophysics : an international weekly journal},
  publisher = {EDP Sciences},
  address   = {Les Ulis},
  issn      = {1432-0746},
  doi       = {10.1051/0004-6361/201629844},
  pages     = {36},
  year      = {2017},
  abstract  = {Context: Massive binaries play a crucial role in the Universe. Knowing the distributions of their orbital parameters is important for a wide range of topics from stellar feedback to binary evolution channels and from the distribution of supernova types to gravitational wave progenitors, yet no direct measurements exist outside the Milky Way. Aims: The Tarantula Massive Binary Monitoring project was designed to help fill this gap by obtaining multi-epoch radial velocity (RV) monitoring of 102 massive binaries in the 30 Doradus region. Methods: In this paper we analyze 32 FLAMES/GIRAFFE observations of 93 O- and 7 B-type binaries. We performed a Fourier analysis and obtained orbital solutions for 82 systems: 51 single-lined (SB1) and 31 double-lined (SB2) spectroscopic binaries. Results: Overall, the binary fraction and orbital properties across the 30 Doradus region are found to be similar to existing Galactic samples. This indicates that within these domains environmental effects are of second order in shaping the properties of massive binary systems. A small difference is found in the distribution of orbital periods, which is slightly flatter (in log space) in 30 Doradus than in the Galaxy, although this may be compatible within error estimates and differences in the fitting methodology. Also, orbital periods in 30 Doradus can be as short as 1.1 d, somewhat shorter than seen in Galactic samples. Equal mass binaries (q> 0.95) in 30 Doradus are all found outside NGC 2070, the central association that surrounds R136a, the very young and massive cluster at 30 Doradus's core. Most of the differences, albeit small, are compatible with expectations from binary evolution. One outstanding exception, however, is the fact that earlier spectral types (O2-O7) tend to have shorter orbital periods than later spectral types (O9.2-O9.7). Conclusions: Our results point to a relative universality of the incidence rate of massive binaries and their orbital properties in the metallicity range from solar (Z⊙) to about half solar. This provides the first direct constraints on massive binary properties in massive star-forming galaxies at the Universe's peak of star formation at redshifts z ~ 1 to 2 which are estimated to have Z ~ 0.5 Z⊙.},
  language  = {en}
}
@article{CorcoranNicholsPabloetal.2015,
  author    = {Corcoran, Michael F. and Nichols, Joy S. and Pablo, Herbert and Shenar, Tomer and Pollock, Andy M. T. and Waldron, Wayne L. and Moffat, Anthony F. J. and Richardson, Noel D. and Russell, Christopher M. P. and Hamaguchi, Kenji and Huenemoerder, David P. and Oskinova, Lida and Hamann, Wolf-Rainer and Naze, Yael and Ignace, Richard and Evans, Nancy Remage and Lomax, Jamie R. and Hoffman, Jennifer L. and Gayley, Kenneth and Owocki, Stanley P. and Leutenegger, Maurice and Gull, Theodore R. and Hole, Karen Tabetha and Lauer, Jennifer and Iping, Rosina C.},
  title     = {A coordinated X-Ray and optical campaign of the nearest massive eclipsing binary, delta ORIONIS Aa. I. Overview of thr X-Ray spectrum},
  series = {The astrophysical journal : an international review of spectroscopy and astronomical physics},
  volume    = {809},
  journal   = {The astrophysical journal : an international review of spectroscopy and astronomical physics},
  number    = {2},
  publisher = {IOP Publ. Ltd.},
  address   = {Bristol},
  issn      = {0004-637X},
  doi       = {10.1088/0004-637X/809/2/132},
  pages     = {15},
  year      = {2015},
  abstract  = {We present an overview of four deep phase-constrained Chandra HETGS X-ray observations of delta Ori A. Delta Ori A is actually a triple system that includes the nearest massive eclipsing spectroscopic binary, delta Ori Aa, the only such object that can be observed with little phase-smearing with the Chandra gratings. Since the fainter star, delta Ori Aa2, has a much lower X-ray luminosity than the brighter primary (delta Ori Aa1), delta Ori Aa provides a unique system with which to test the spatial distribution of the X-ray emitting gas around delta Ori Aa1 via occultation by the photosphere of, and wind cavity around, the X-ray dark secondary. Here we discuss the X-ray spectrum and X-ray line profiles for the combined observation, having an exposure time of nearly 500 ks and covering nearly the entire binary orbit. The companion papers discuss the X-ray variability seen in the Chandra spectra, present new space-based photometry and ground-based radial velocities obtained simultaneously with the X-ray data to better constrain the system parameters, and model the effects of X-rays on the optical and UV spectra. We find that the X-ray emission is dominated by embedded wind shock emission from star Aa1, with little contribution from the tertiary star Ab or the shocked gas produced by the collision of the wind of Aa1 against the surface of Aa2. We find a similar temperature distribution to previous X-ray spectrum analyses. We also show that the line half-widths are about 0.3-0.5 times the terminal velocity of the wind of star Aa1. We find a strong anti-correlation between line widths and the line excitation energy, which suggests that longer-wavelength, lower-temperature lines form farther out in the wind. Our analysis also indicates that the ratio of the intensities of the strong and weak lines of Fe XVII and Ne X are inconsistent with model predictions, which may be an effect of resonance scattering.},
  language  = {en}
}
@article{HainichOskinovaShenaretal.2018,
  author    = {Hainich, Rainer and Oskinova, Lida and Shenar, Tomer and Marchant Campos, Pablo and Eldridge, J. J. and Sander, Andreas Alexander Christoph and Hamann, Wolf-Rainer and Langer, Norbert and Todt, Helge Tobias},
  title     = {Observational properties of massive black hole binary progenitors},
  series = {Astronomy and astrophysics : an international weekly journal},
  volume    = {609},
  journal   = {Astronomy and astrophysics : an international weekly journal},
  publisher = {EDP Sciences},
  address   = {Les Ulis},
  issn      = {1432-0746},
  doi       = {10.1051/0004-6361/201731449},
  pages     = {62},
  year      = {2018},
  abstract  = {Context: The first directly detected gravitational waves (GW 150914) were emitted by two coalescing black holes (BHs) with masses of ≈ 36 M⊙ and ≈ 29 M⊙. 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. Aims: 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. Methods: Stellar evolution models predict fundamental stellar parameters. Using these as input for our stellar atmosphere code (Potsdam Wolf-Rayet), we compute a set of models for selected evolutionary stages of massive merging BH progenitors at different metallicities. Results: 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. Based on our atmosphere models, we provide feedback parameters, broad band magnitudes, and spectral templates that should help to identify such binaries in the future. Conclusions: While the predicted parameter space for massive BH binary progenitors is partly realized in nature, none of the known massive binaries match our synthetic spectra of massive BH binary progenitors exactly. 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.},
  language  = {en}
}
@article{HainichPasemannTodtetal.2015,
  author    = {Hainich, Rainer and Pasemann, Diana and Todt, Helge Tobias and Shenar, Tomer and Sander, Andreas Alexander Christoph and Hamann, Wolf-Rainer},
  title     = {Wolf-Rayet stars in the Small Magellanic Cloud I. Analysis of the single WN stars},
  series = {Astronomy and astrophysics : an international weekly journal},
  volume    = {581},
  journal   = {Astronomy and astrophysics : an international weekly journal},
  publisher = {EDP Sciences},
  address   = {Les Ulis},
  issn      = {1432-0746},
  doi       = {10.1051/0004-6361/201526241},
  pages     = {30},
  year      = {2015},
  abstract  = {Context. 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. Aims. 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. Methods. 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. Results. 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) L-circle dot. 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. Conclusions. 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.},
  language  = {en}
}
@article{HainichRamachandranShenaretal.2019,
  author    = {Hainich, Rainer and Ramachandran, Varsha and Shenar, Tomer and Sander, Andreas Alexander Christoph and Todt, Helge Tobias and Gruner, David and Oskinova, Lida and Hamann, Wolf-Rainer},
  title     = {PoWR grids of non-LTE model atmospheres for OB-type stars of various metallicities},
  series = {Astronomy and astrophysics : an international weekly journal},
  volume    = {621},
  journal   = {Astronomy and astrophysics : an international weekly journal},
  publisher = {EDP Sciences},
  address   = {Les Ulis},
  issn      = {1432-0746},
  doi       = {10.1051/0004-6361/201833787},
  pages     = {12},
  year      = {2019},
  abstract  = {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.},
  language  = {en}
}
@article{PabloRichardsonMoffatetal.2015,
  author    = {Pablo, Herbert and Richardson, Noel D. and Moffat, Anthony F. J. and Corcoran, Michael and Shenar, Tomer and Benvenuto, Omar and Fuller, Jim and Naze, Yael and Hoffman, Jennifer L. and Miroshnichenko, Anatoly and Apellaniz, Jesus Maiz and Evans, Nancy and Eversberg, Thomas and Gayley, Ken and Gull, Ted and Hamaguchi, Kenji and Hamann, Wolf-Rainer and Henrichs, Huib and Hole, Tabetha and Ignace, Richard and Iping, Rosina and Lauer, Jennifer and Leutenegger, Maurice and Lomax, Jamie and Nichols, Joy and Oskinova, Lida and Owocki, Stan and Pollock, Andy and Russell, Christopher M. P. and Waldron, Wayne and Buil, Christian and Garrel, Thierry and Graham, Keith and Heathcote, Bernard and Lemoult, Thierry and Li, Dong and Mauclaire, Benjamin and Potter, Mike and Ribeiro, Jose and Matthews, Jaymie and Cameron, Chris and Guenther, David and Kuschnig, Rainer and Rowe, Jason and Rucinski, Slavek and Sasselov, Dimitar and Weiss, Werner},
  title     = {A coordinated X-Ray and optical campaign of the nearest massive eclipsing binary, delta ORIONIS Aa. III. Analysis of optical photometric (most) and spectroscopic (ground based) variations},
  series = {The astrophysical journal : an international review of spectroscopy and astronomical physics},
  volume    = {809},
  journal   = {The astrophysical journal : an international review of spectroscopy and astronomical physics},
  number    = {2},
  publisher = {IOP Publ. Ltd.},
  address   = {Bristol},
  issn      = {0004-637X},
  doi       = {10.1088/0004-637X/809/2/134},
  pages     = {11},
  year      = {2015},
  abstract  = {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 delta 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 P > 400 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.},
  language  = {en}
}
@article{SanderFuerstKretschmaretal.2018,
  author    = {Sander, Andreas Alexander Christoph and F{\"u}rst, F. and Kretschmar, P. and Oskinova, Lida and Todt, Helge Tobias and Hainich, Rainer and Shenar, Tomer and Hamann, Wolf-Rainer},
  title     = {Coupling hydrodynamics with comoving frame radiative transfer},
  series = {Astronomy and astrophysics : an international weekly journal},
  volume    = {610},
  journal   = {Astronomy and astrophysics : an international weekly journal},
  publisher = {EDP Sciences},
  address   = {Les Ulis},
  issn      = {1432-0746},
  doi       = {10.1051/0004-6361/201731575},
  pages     = {19},
  year      = {2018},
  abstract  = {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 used the recently updated version of the Potsdam Wolf-Rayet 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 u(infinity) approximate to 600 km s(-1). On the other hand, the wind velocity in the inner region where the NS is located is only approximate to 100 km s(-1), 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.},
  language  = {en}
}
@article{SanderHamannTodtetal.2017,
  author    = {Sander, Andreas Alexander Christoph and Hamann, Wolf-Rainer and Todt, Helge Tobias and Hainich, Rainer and Shenar, Tomer},
  title     = {Coupling hydrodynamics with comoving frame radiative transfer I. A unified approach for OB and WR stars},
  series = {Astronomy and astrophysics : an international weekly journal},
  volume    = {603},
  journal   = {Astronomy and astrophysics : an international weekly journal},
  publisher = {EDP Sciences},
  address   = {Les Ulis},
  issn      = {1432-0746},
  doi       = {10.1051/0004-6361/201730642},
  pages     = {14},
  year      = {2017},
  abstract  = {Context. For more than two decades, stellar atmosphere codes have been used to derive the stellar and wind parameters of massive stars. Although they have become a powerful tool and sufficiently reproduce the observed spectral appearance, they can hardly be used for more than measuring parameters. One major obstacle is their inconsistency between the calculated radiation field and the wind stratification due to the usage of prescribed mass-loss rates and wind-velocity fields. Aims. We present the concepts for a new generation of hydrodynamically consistent non-local thermodynamical equilibrium (nonLTE) stellar atmosphere models that allow for detailed studies of radiation-driven stellar winds. As a first demonstration, this new kind of model is applied to a massive O star. Methods. Based on earlier works, the PoWR code has been extended with the option to consistently solve the hydrodynamic equation together with the statistical equations and the radiative transfer in order to obtain a hydrodynamically consistent atmosphere stratification. In these models, the whole velocity field is iteratively updated together with an adjustment of the mass-loss rate. Results. The concepts for obtaining hydrodynamically consistent models using a comoving-frame radiative transfer are outlined. To provide a useful benchmark, we present a demonstration model, which was motivated to describe the well-studied O4 supergiant zeta Pup. The obtained stellar and wind parameters are within the current range of literature values. Conclusions. For the first time, the PoWR code has been used to obtain a hydrodynamically consistent model for a massive O star. This has been achieved by a profound revision of earlier concepts used for Wolf-Rayet stars. The velocity field is shaped by various elements contributing to the radiative acceleration, especially in the outer wind. The results further indicate that for more dense winds deviations from a standard beta-law occur.},
  language  = {en}
}
@article{SanderShenarHainichetal.2015,
  author    = {Sander, Andreas Alexander Christoph and Shenar, Tomer and Hainich, Rainer and Gimenez-Garcia, Ana and Todt, Helge Tobias and Hamann, Wolf-Rainer},
  title     = {On the consistent treatment of the quasi-hydrostatic layers in hot star atmospheres},
  series = {Astronomy and astrophysics : an international weekly journal},
  volume    = {577},
  journal   = {Astronomy and astrophysics : an international weekly journal},
  publisher = {EDP Sciences},
  address   = {Les Ulis},
  issn      = {0004-6361},
  doi       = {10.1051/0004-6361/201425356},
  pages     = {13},
  year      = {2015},
  abstract  = {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.},
  language  = {en}
}
@article{ShenarOskinovaHamannetal.2015,
  author    = {Shenar, Tomer and Oskinova, Lida and Hamann, Wolf-Rainer and Corcoran, Michael F. and Moffat, Anthony F. J. and Pablo, Herbert and Richardson, Noel D. and Waldron, Wayne L. and Huenemoerder, David P. and Maiz Apellaniz, Jesus and Nichols, Joy S. and Todt, Helge Tobias and Naze, Yael and Hoffman, Jennifer L. and Pollock, Andy M. T. and Negueruela, Ignacio},
  title     = {A coordinated X-Ray and optical campaign of the nearest massive eclipsing binary, delta ORIONIS Aa. IV. A multiwavelength, non-lte spectroscopic analysis},
  series = {The astrophysical journal : an international review of spectroscopy and astronomical physics},
  volume    = {809},
  journal   = {The astrophysical journal : an international review of spectroscopy and astronomical physics},
  number    = {2},
  publisher = {IOP Publ. Ltd.},
  address   = {Bristol},
  issn      = {0004-637X},
  doi       = {10.1088/0004-637X/809/2/135},
  pages     = {20},
  year      = {2015},
  abstract  = {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 Hipparcos parallax yields surprisingly low values for the mass, radius, and luminosity. Consistent values follow only if delta Ori lies at about twice the 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 approximate to 2.(m)8. The inferred parameters suggest that the secondary is an early B-type dwarf (approximate to B1 V), while the tertiary is an early B-type subgiant (approximate to B0 IV). We find evidence for rapid turbulent velocities (similar to 200 km 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 (logL(X)/L-Bol approximate to -6.85), initiating close to the stellar surface at R-0 similar to 1.1 R-*. Accounting for clumping, the mass-loss rate of the primary is found to be log (M) over dot approximate to -6.4 (M-circle dot yr(-1))., which agrees with hydrodynamic predictions, and provides a consistent picture along the X-ray, UV, optical, and radio spectral domains.},
  language  = {en}
}