@article{RamiaramanantsoaMoffatHarmonetal.2018,
  author    = {Ramiaramanantsoa, Tahina and Moffat, Anthony F. J. and Harmon, Robert and Ignace, R. and St-Louis, Nicole and Vanbeveren, Dany and Shenar, Tomer and Pablo, Herbert and Richardson, Noel D. and Howarth, Ian D. and Stevens, Ian R. and Piaulet, Caroline and St-Jean, Lucas and Eversberg, Thomas and Pigulski, Andrzej and Popowicz, Adam and Kuschnig, Rainer and Zoclonska, Elzbieta and Buysschaert, Bram and Handler, Gerald and Weiss, Werner W. and Wade, Gregg A. and Rucinski, Slavek M. and Zwintz, Konstanze and Luckas, Paul and Heathcote, Bernard and Cacella, Paulo and Powles, Jonathan and Locke, Malcolm and Bohlsen, Terry and Chen{\´e}, Andr{\´e}-Nicolas and Miszalski, Brent and Waldron, Wayne L. and Kotze, Marissa M. and Kotze, Enrico J. and B{\"o}hm, Torsten},
  title     = {BRITE-Constellation high-precision time-dependent photometry of the early O-type supergiant zeta Puppis unveils the photospheric drivers of its small- and large-scale wind structures},
  series = {Monthly notices of the Royal Astronomical Society},
  volume    = {473},
  journal   = {Monthly notices of the Royal Astronomical Society},
  number    = {4},
  publisher = {Oxford Univ. Press},
  address   = {Oxford},
  issn      = {0035-8711},
  doi       = {10.1093/mnras/stx2671},
  pages     = {5532 -- 5569},
  year      = {2018},
  abstract  = {From 5.5 months of dual-band optical photometric monitoring at the 1 mmag level, BRITE-Constellation has revealed two simultaneous types of variability in the O4I(n)fp star ζ Puppis: one single periodic non-sinusoidal component superimposed on a stochastic component. The monoperiodic component is the 1.78-d signal previously detected by Coriolis/Solar Mass Ejection Imager, but this time along with a prominent first harmonic. The shape of this signal changes over time, a behaviour that is incompatible with stellar oscillations but consistent with rotational modulation arising from evolving bright surface inhomogeneities. By means of a constrained non-linear light-curve inversion algorithm, we mapped the locations of the bright surface spots and traced their evolution. Our simultaneous ground-based multisite spectroscopic monitoring of the star unveiled cyclical modulation of its He ii λ4686 wind emission line with the 1.78-d rotation period, showing signatures of corotating interaction regions that turn out to be driven by the bright photospheric spots observed by BRITE. Traces of wind clumps are also observed in the He ii λ4686 line and are correlated with the amplitudes of the stochastic component of the light variations probed by BRITE at the photosphere, suggesting that the BRITE observations additionally unveiled the photospheric drivers of wind clumps in ζ Pup and that the clumping phenomenon starts at the very base of the wind. The origins of both the bright surface inhomogeneities and the stochastic light variations remain unknown, but a subsurface convective zone might play an important role in the generation of these two types of photospheric variability.},
  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{RamachandranHamannHainichetal.2018,
  author    = {Ramachandran, Varsha and Hamann, Wolf-Rainer and Hainich, Rainer and Oskinova, Lida and Shenar, Tomer and Sander, Andreas Alexander Christoph and Todt, Helge Tobias and Gallagher, John S.},
  title     = {Stellar population of the superbubble N206 in the LMC II. Parameters of the OB and WR stars, and the total massive star feedback},
  series = {Astronomy and astrophysics : an international weekly journal},
  volume    = {615},
  journal   = {Astronomy and astrophysics : an international weekly journal},
  publisher = {EDP Sciences},
  address   = {Les Ulis},
  issn      = {1432-0746},
  doi       = {10.1051/0004-6361/201832816},
  pages     = {72},
  year      = {2018},
  abstract  = {Context. 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 (LMC) exhibits a superbubble and a rich massive star population. Aims. Our goal is to perform quantitative spectral analyses of all massive stars associated with the N206 superbubble in order to determine their stellar and wind parameters. We compare the superbubble energy budget to the stellar energy input and discuss the star formation history of the region. Results. We present the stellar and wind parameters of the OB stars and the two Wolf-Rayet (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 Hertzsprung-Russell diagram (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. The total ionizing photon flux produced by all massive stars in the N206 complex is Q(0) approximate to 5 x 10(50) s(-1), and the mechanical luminosity of their stellar winds amounts to L-mec = 1.7 x 10(38) erg s(-1). 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 approximate to 2.3 x 10(52) erg. Conclusions. The N206 complex in the LMC has undergone star formation episodes since more than 30 Myr ago. From the spectral analyses of its massive star population, we derive a current star formation rate of 2.2 x 10(-3) M-circle dot yr(-1). From the combined input of mechanical energy from all stellar winds, only a minor fraction is emitted in the form of X-rays. The corresponding input accumulated over a long time also 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.},
  language  = {en}
}
@article{KubatovaSzecsiSanderetal.2019,
  author    = {Kubatova, Brankica and Szecsi, D. and Sander, Andreas Alexander Christoph and Kubat, Jiř{\´i} and Tramper, F. and Krticka, Jiri and Kehrig, C. and Hamann, Wolf-Rainer and Hainich, Rainer and Shenar, Tomer},
  title     = {Low-metallicity massive single stars with rotation},
  series = {Astronomy and astrophysics : an international weekly journal},
  volume    = {623},
  journal   = {Astronomy and astrophysics : an international weekly journal},
  publisher = {EDP Sciences},
  address   = {Les Ulis},
  issn      = {1432-0746},
  doi       = {10.1051/0004-6361/201834360},
  pages     = {32},
  year      = {2019},
  abstract  = {Context. 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(circle dot)) 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-circle dot) 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.},
  language  = {en}
}
@article{HamannGraefenerLiermannetal.2019,
  author    = {Hamann, Wolf-Rainer and Gr{\"a}fener, G. and Liermann, A. and Hainich, Rainer and Sander, Andreas Alexander Christoph and Shenar, Tomer and Ramachandran, Varsha and Todt, Helge Tobias and Oskinova, Lida},
  title     = {The Galactic WN stars revisited},
  series = {Astronomy and astrophysics : an international weekly journal},
  volume    = {625},
  journal   = {Astronomy and astrophysics : an international weekly journal},
  publisher = {EDP Sciences},
  address   = {Les Ulis},
  issn      = {1432-0746},
  doi       = {10.1051/0004-6361/201834850},
  pages     = {11},
  year      = {2019},
  abstract  = {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 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.},
  language  = {en}
}
@article{MunozMoffatHilletal.2017,
  author    = {Munoz, Melissa and Moffat, Anthony F. J. and Hill, Grant M. and Shenar, Tomer and Richardson, Noel D. and Pablo, Herbert and St-Louis, Nicole and Ramiaramanantsoa, Tahina},
  title     = {WR 148: identifying the companion of an extreme runaway massive binary},
  series = {Monthly notices of the Royal Astronomical Society},
  volume    = {467},
  journal   = {Monthly notices of the Royal Astronomical Society},
  publisher = {Oxford Univ. Press},
  address   = {Oxford},
  issn      = {0035-8711},
  doi       = {10.1093/mnras/stw2283},
  pages     = {3105 -- 3121},
  year      = {2017},
  abstract  = {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\&\#8201;km\&\#8201;s\&\#8722;1 and 79.2 ± 3.1\&\#8201;km\&\#8201;s\&\#8722;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 \&\#8764;67°, derived from previous polarimetry observations, the system's total mass would be a mere 2-3M\&\#8857;\&\#8288;, an unprecedented result for a putative massive binary system. However, a system comprising a 37M\&\#8857; secondary (typical mass of an O5V star) and a 33M\&\#8857; primary (given the mass ratio) would infer an inclination of \&\#8764;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\&\#8201;km\&\#8201;s\&\#8722;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.},
  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{OskinovaHuenemoerderHamannetal.2017,
  author    = {Oskinova, Lida and Huenemoerder, D. P. and Hamann, Wolf-Rainer and Shenar, Tomer and Sander, Andreas Alexander Christoph and Ignace, R. and Todt, Helge Tobias and Hainich, Rainer},
  title     = {On the Binary Nature of Massive Blue Hypergiants: High-resolution X-Ray Spectroscopy Suggests That Cyg OB2 12 is a Colliding Wind Binary},
  series = {The astrophysical journal : an international review of spectroscopy and astronomical physics},
  volume    = {845},
  journal   = {The astrophysical journal : an international review of spectroscopy and astronomical physics},
  publisher = {IOP Publ. Ltd.},
  address   = {Bristol},
  issn      = {0004-637X},
  doi       = {10.3847/1538-4357/aa7e79},
  pages     = {11},
  year      = {2017},
  abstract  = {The blue hypergiant Cyg OB2 12 (B3Ia(+)) is a representative member of the class of very massive stars in a poorly understood evolutionary stage. We obtained its high-resolution X-ray spectrum using the Chandra observatory. PoWR model atmospheres were calculated to provide realistic wind opacities and to establish the wind density structure. We find that collisional de-excitation is the dominant mechanism depopulating the metastable upper levels of the forbidden lines of the He-like ions Si XIV and Mg XII. Comparison between the model and observations reveals that X-ray emission is produced in a dense plasma, which could reside only at the photosphere or in a colliding wind zone between binary components. The observed X-ray spectra are well-fitted by thermal plasma models, with average temperatures in excess of 10 MK. The wind speed in Cyg OB2 12 is not high enough to power such high temperatures, but the collision of two winds in a binary system can be sufficient. We used archival data to investigate the X-ray properties of other blue hypergiants. In general, stars of this class are not detected as X-ray sources. We suggest that our new Chandra observations of Cyg OB2 12 can be best explained if Cyg OB2 12 is a colliding wind binary possessing a late O-type companion. This makes Cyg OB2 12 only the second binary system among the 16 known Galactic hypergiants. This low binary fraction indicates that the blue hypergiants are likely products of massive binary evolution during which they either accreted a significant amount of mass or already merged with their companions.},
  language  = {en}
}
@article{GimenezGarciaShenarTorrejonetal.2016,
  author    = {Gimenez-Garcia, Ana and Shenar, Tomer and Torrejon, J. M. and Oskinova, Lida and Martinez-Nunez, S. and Hamann, Wolf-Rainer and Rodes-Roca, J. J. and Gonz{\´a}lez-Galan, A. and Alonso-Santiago, J. and Gonz{\´a}lez-Fern{\´a}ndez, C. and Bernabeu, Guillermo and Sander, Andreas Alexander Christoph},
  title     = {Measuring the stellar wind parameters in IGR J17544-2619 and Vela X-1 constrains the accretion physics in supergiant fast X-ray transient and classical supergiant X-ray binaries},
  series = {Siberian Mathematical Journal},
  volume    = {591},
  journal   = {Siberian Mathematical Journal},
  publisher = {EDP Sciences},
  address   = {Les Ulis},
  issn      = {1432-0746},
  doi       = {10.1051/0004-6361/201527551},
  pages     = {25},
  year      = {2016},
  abstract  = {Aims. To close this gap, we perform a comparative analysis of the optical companion in two important systems: IGR J175442619 (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 (v(infinity) = 1500 km s(-1) in IGR J17544-2619 and v(infinity) = 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.},
  language  = {en}
}
@article{RichardsonShenarRoyLoubieretal.2016,
  author    = {Richardson, Noel D. and Shenar, Tomer and Roy-Loubier, Olivier and Schaefer, Gail and Moffat, Anthony F. J. and St-Louis, Nicole and Gies, Douglas R. and Farrington, Chris and Hill, Grant M. and Williams, Peredur M. and Gordon, Kathryn and Pablo, Herbert and Ramiaramanantsoa, Tahina},
  title     = {The CHARA Array resolves the long-period Wolf-Rayet binaries WR 137 and WR 138},
  series = {Monthly notices of the Royal Astronomical Society},
  volume    = {461},
  journal   = {Monthly notices of the Royal Astronomical Society},
  publisher = {Oxford Univ. Press},
  address   = {Oxford},
  issn      = {0035-8711},
  doi       = {10.1093/mnras/stw1585},
  pages     = {4115 -- 4124},
  year      = {2016},
  abstract  = {We report on interferometric observations with the CHARAArray of two classical Wolf-Rayet (WR) 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 (fWR, 137 = 0.59 +/- 0.04; fWR, 138 = 0.67 +/- 0.01), which is confirmed through both comparisons with estimated fundamental parameters for WR stars and O dwarfs, as well as through spectral modelling of each system. Our spectral modelling 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 edgeon. 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.},
  language  = {en}
}