@article{KrtičkaFeldmeier2021, author = {Krtička, Jiř{\´i} and Feldmeier, Achim}, title = {Stochastic light variations in hot stars from wind instability}, series = {Astronomy and astrophysics : an international weekly journal / European Southern Observatory (ESO)}, volume = {648}, journal = {Astronomy and astrophysics : an international weekly journal / European Southern Observatory (ESO)}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {1432-0746}, doi = {10.1051/0004-6361/202040148}, pages = {9}, year = {2021}, abstract = {Context Line-driven wind instability is expected to cause small-scale wind inhomogeneities, X-ray emission, and wind line profile variability. The instability can already develop around the sonic point if it is initiated close to the photosphere due to stochastic turbulent motions. In such cases, it may leave its imprint on the light curve as a result of wind blanketing. Aims We study the photometric signatures of the line-driven wind instability. Methods We used line-driven wind instability simulations to determine the wind variability close to the star. We applied two types of boundary perturbations: a sinusoidal one that enables us to study in detail the development of the instability and a stochastic one given by a Langevin process that provides a more realistic boundary perturbation. We estimated the photometric variability from the resulting mass-flux variations. The variability was simulated assuming that the wind consists of a large number of independent conical wind sectors. We compared the simulated light curves with TESS light curves of OB stars that show stochastic variability. Results We find two typical signatures of line-driven wind instability in photometric data: a knee in the power spectrum of magnitude fluctuations, which appears due to engulfment of small-scale structure by larger structures, and a negative skewness of the distribution of fluctuations, which is the result of spatial dominance of rarefied regions. These features endure even when combining the light curves from independent wind sectors. Conclusions The stochastic photometric variability of OB stars bears certain signatures of the line-driven wind instability. The distribution function of observed photometric data shows negative skewness and the power spectra of a fraction of light curves exhibit a knee. This can be explained as a result of the line-driven wind instability triggered by stochastic base perturbations.}, language = {en} } @article{SanderVinkHamann2019, author = {Sander, Andreas Alexander Christoph and Vink, Jorick S. and Hamann, Wolf-Rainer}, title = {Driving classical Wolf-Rayet winds}, series = {Monthly notices of the Royal Astronomical Society}, volume = {491}, journal = {Monthly notices of the Royal Astronomical Society}, number = {3}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/stz3064}, pages = {4406 -- 4425}, year = {2019}, abstract = {Classical Wolf-Rayet (cWR) stars are at a crucial evolutionary stage for constraining the fates of massive stars. The feedback of these hot, hydrogen-depleted stars dominates their surrounding by tremendous injections of ionizing radiation and kinetic energy. The strength of a Wolf-Rayet (WR) wind decides the eventual mass of its remnant, likely a massive black hole. However, despite their major influence and importance for gravitational wave detection statistics, WR winds are particularly poorly understood. In this paper, we introduce the first set of hydrodynamically consistent stellar atmosphere models for cWR stars of both the carbon (C) and the nitrogen (N) sequence, i.e. WC and WN stars, as a function of stellar luminosity-to-mass ratio (or Eddington Gamma) and metallicity. We demonstrate the inapplicability of the CAK wind theory for cWR stars and confirm earlier findings that their winds are launched at the (hot) iron (Fe) opacity peak. For log Z/Z(circle dot) > -2, Fe is also the main accelerator throughout the wind. Contrasting previous claims of a sharp lower mass-loss limit forWR stars, we obtain a smooth transition to optically thin winds. Furthermore, we find a strong dependence of the mass-loss rates on Eddington Gamma, both at solar and subsolar metallicity. Increases inWCcarbon and oxygen abundances turn out to slightly reduce the predicted mass-loss rates. Calculations at subsolar metallicities indicate that below the metallicity of the Small Magellanic Cloud, WR mass-loss rates decrease much faster than previously assumed, potentially allowing for high black hole masses even in the local Universe.}, language = {en} } @article{IrrgangGeierKreuzeretal.2020, author = {Irrgang, Andreas and Geier, Stephan and Kreuzer, Simon and Pelisoli, Ingrid Domingos and Heber, Ulrich}, title = {A stripped helium star in the potential black hole binary LB-1}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {633}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {0004-6361}, doi = {10.1051/0004-6361/201937343}, pages = {9}, year = {2020}, abstract = {Context The recently claimed discovery of a massive (M-BH = 68(-13)(+11) M-circle dot) black hole in the Galactic solar neighborhood has led to controversial discussions because it severely challenges our current view of stellar evolution. Aims A crucial aspect for the determination of the mass of the unseen black hole is the precise nature of its visible companion, the B-type star LSV +22 25. Because stars of different mass can exhibit B-type spectra during the course of their evolution, it is essential to obtain a comprehensive picture of the star to unravel its nature and, thus, its mass. Methods To this end, we study the spectral energy distribution of LSV +22 25 and perform a quantitative spectroscopic analysis that includes the determination of chemical abundances for He, C, N, O, Ne, Mg, Al, Si, S, Ar, and Fe. Results Our analysis clearly shows that LSV +22 25 is not an ordinary main sequence B-type star. The derived abundance pattern exhibits heavy imprints of the CNO bi-cycle of hydrogen burning, that is, He and N are strongly enriched at the expense of C and O. Moreover, the elements Mg, Al, Si, S, Ar, and Fe are systematically underabundant when compared to normal main-sequence B-type stars. We suggest that LSV +22 25 is a stripped helium star and discuss two possible formation scenarios. Combining our photometric and spectroscopic results with the Gaia parallax, we infer a stellar mass of 1.1 +/- 0.5 M-circle dot. Based on the binary system's mass function, this yields a minimum mass of 2-3 M-circle dot for the compact companion, which implies that it may not necessarily be a black hole but a massive neutron- or main sequence star. Conclusions The star LSV +22 25 has become famous for possibly having a very massive black hole companion. However, a closer look reveals that the star itself is a very intriguing object. Further investigations are necessary for complete characterization of this object.}, language = {en} } @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{SchoellerHubrigFossatietal.2017, author = {Sch{\"o}ller, Markus and Hubrig, Swetlana and Fossati, L. and Carroll, Thorsten Anthony and Briquet, Maryline and Oskinova, Lida and J{\"a}rvinen, S. and Ilyin, Ilya and Castro, N. and Morel, T. and Langer, N. and Przybilla, N. and Nieva, M. -F. and Kholtygin, A. F. and Sana, H. and Herrero, A. and Barba, R. H. and de Koter, A.}, title = {B fields in OB stars (BOB)}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {599}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, organization = {BOB Collaboration}, issn = {1432-0746}, doi = {10.1051/0004-6361/201628905}, pages = {11}, year = {2017}, abstract = {Aims. The B fields in OB stars (BOB) Collaboration is based on an ESO Large Programme to study the occurrence rate, properties, and ultimately the origin of magnetic fields in massive stars. Methods. In the framework of this program, we carried out low-resolution spectropolarimetric observations of a large sample of massive stars using FORS2 installed at the ESO VLT 8m telescope. Results. We determined the magnetic field values with two completely independent reduction and analysis pipelines. Our in-depth study of the magnetic field measurements shows that differences between our two pipelines are usually well within 3 sigma errors. From the 32 observations of 28 OB stars, we were able to monitor the magnetic fields in CPD -57 degrees 3509 and HD164492C, confirm the magnetic field in HD54879, and detect a magnetic field in CPD -62 degrees 2124. We obtain a magnetic field detection rate of 6 +/- 3\% for the full sample of 69 OB stars observed with FORS 2 within the BOB program. For the preselected objects with a nu sin i below 60 km s(-1), we obtain a magnetic field detection rate of 5 +/- 5\%. We also discuss X-ray properties and multiplicity of the objects in our FORS2 sample with respect to the magnetic field detections.}, language = {en} } @article{RamachandranHainichHamannetal.2017, author = {Ramachandran, Varsha and Hainich, Rainer and Hamann, Wolf-Rainer and Oskinova, Lida and Shenar, T. and Sander, Andreas Alexander Christoph and Todt, Helge Tobias and Gallagher, John S.}, title = {Stellar population of the superbubble N206 in the LMC I. Analysis of the Of-type stars}, 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/201731093}, pages = {26}, year = {2017}, abstract = {Context. Massive stars severely influence their environment by their strong ionizing radiation and by the momentum and kinetic energy input provided by their stellar winds and supernovae. Quantitative analyses of massive stars are required to understand how their feedback creates and shapes large scale structures of the interstellar medium. The giant H II region N206 in the Large Magellanic Cloud contains an OB association that powers a superbubble filled with hot X-ray emitting gas, serving as an ideal laboratory in this context. Aims. We aim to estimate stellar and wind parameters of all OB stars in N206 by means of quantitative spectroscopic analyses. In this first paper, we focus on the nine Of-type stars located in this region. We determine their ionizing flux and wind mechanical energy. The analysis of nitrogen abundances in our sample probes rotational mixing. Methods. We obtained optical spectra with the multi-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 used the Potsdam Wolf-Rayet model atmosphere code. We determined the physical parameters and nitrogen abundances of our sample stars by fitting synthetic spectra to the observations. Results. The stellar and wind parameters of nine Of-type stars, which are largely derived from spectral analysis 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 that 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 of less than 4 million yr, 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. Conclusions. 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.}, 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{PillitteriFossatiRodriguezetal.2017, author = {Pillitteri, Ignazio and Fossati, Luca and Rodriguez, N. Castro and Oskinova, Lida and Wolk, Scott J.}, title = {Detection of magnetic field in the B2 star rho Ophiuchi A with ESO FORS2}, 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/201732078}, pages = {4}, year = {2017}, abstract = {Circumstantial evidence suggests that magnetism and enhanced X-ray emission are likely correlated in early B-type stars: similar fractions of them (similar to 10\%) are strong and hard X-ray sources and possess strong magnetic fields. It is also known that some B-type stars have spots on their surface. Yet up to now no X-ray activity associated with spots on early-type stars was detected. In this Letter we report the detection of a magnetic field on the B2V star rho Oph A. Previously, we assessed that the X-ray activity of this star is associated with a surface spot, herewith we establish its magnetic origin. We analyze spectra of rho Oph A obtained with the FORS2 spectrograph at ESO Very Large Telescope (VLT) at two epochs, and detect a longitudinal component of the magnetic field of the order of similar to 500 G in one of the datasets. The detection of the magnetic field only at one epoch can be explained by stellar rotation which is also invoked to explain observed periodic X-ray activity. From archival HARPS ESO VLT high resolution spectra we derived the fundamental stellar parameters of rho Oph A and further constrained its age. We conclude that rho Oph A provides strong evidence for the presence of active X-ray emitting regions on young magnetized early type stars.}, 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{AbdallaAbramowskiAharonianetal.2018, author = {Abdalla, Hassan E. and Abramowski, A. and Aharonian, Felix A. and Benkhali, F. Ait and Akhperjanian, A. G. and Andersson, T. and Anguener, E. O. and Arakawa, M. and Arrieta, M. and Aubert, P. and Backes, M. and Balzer, A. and Barnard, M. and Becherini, Y. and Tjus, J. Becker and Berge, D. and Bernhard, S. and Bernloehr, K. and Blackwell, R. and Bottcher, M. and Boisson, C. and Bolmont, J. and Bordas, Pol and Bregeon, J. and Brun, F. and Brun, P. and Bryan, M. and Buechele, M. and Bulik, T. and Capasso, M. and Carr, J. and Casanova, Sabrina and Cerruti, M. and Chakraborty, N. and Chalme-Calvet, R. and Chaves, R. C. G. and Chen, A. and Chevalier, J. and Chretien, M. and Coffaro, M. and Colafrancesco, S. and Cologna, G. and Condon, B. and Conrad, J. and Cui, Y. and Davids, I. D. and Decock, J. and Degrange, B. and Deil, C. and Devin, J. and deWilt, P. and Dirson, L. and Djannati-Atai, A. and Domainko, W. and Donath, A. and Dutson, K. and Dyks, J. and Edwards, T. and Egberts, Kathrin and Eger, P. and Ernenwein, J. -P. and Eschbach, S. and Farnier, C. and Fegan, S. and Fernandes, M. V. and Fiasson, A. and Fontaine, G. and Foerster, A. and Funk, S. and Fuessling, M. and Gabici, S. and Gajdus, M. and Gallant, Y. A. and Garrigoux, T. and Giavitto, G. and Giebels, B. and Glicenstein, J. F. and Gottschall, D. and Goyal, A. and Grondin, M. -H. and Hahn, J. and Haupt, M. and Hawkes, J. and Heinzelmann, G. and Henri, G. and Hermann, G. and Hervet, O. and Hinton, J. A. and Hofmann, W. and Hoischen, Clemens and Holler, M. and Horns, D. and Ivascenko, A. and Iwasaki, H. and Jacholkowska, A. and Jamrozy, M. and Janiak, M. and Jankowsky, D. and Jankowsky, F. and Jingo, M. and Jogler, T. and Jouvin, L. and Jung-Richardt, I. and Kastendieck, M. A. and Katarzynski, K. and Katsuragawa, M. and Katz, U. and Kerszberg, D. and Khangulyan, D. and Khelifi, B. and Kieffer, M. and King, J. and Klepser, S. and Klochkov, D. and Kluzniak, W. and Kolitzus, D. and Komin, Nu. and Kosack, K. and Krakau, S. and Kraus, M. and Kruger, P. P. and Laffon, H. and Lamanna, G. and Lau, J. and Lees, J. -P. and Lefaucheur, J. and Lefranc, V. and Lemiere, A. and Lemoine-Goumard, M. and Lenain, J. -P. and Leser, Eva and Lohse, T. and Lorentz, M. and Liu, R. and Lopez-Coto, R. and Lypova, I. and Marandon, V. and Marcowith, A. and Mariaud, C. and Marx, R. and Maurin, G. and Maxted, N. and Mayer, M. and Meintjes, P. J. and Meyer, M. and Mitchell, A. M. W. and Moderski, R. and Mohamed, M. and Mohrmann, L. and Mora, K. and Moulin, E. and Murach, T. and Nakashima, S. and de Naurois, M. and Niederwanger, F. and Niemiec, J. and Oakes, L. and Odaka, H. and Oettl, S. and Ohm, S. and Ostrowski, M. and Oya, I. and Padovani, M. and Panter, M. and Parsons, R. D. and Pekeur, N. W. and Pelletier, G. and Perennes, C. and Petrucci, P. -O. and Peyaud, B. and Piel, Q. and Pita, S. and Poon, H. and Prokhorov, D. and Prokoph, H. and Puehlhofer, G. and Punch, M. and Quirrenbach, A. and Raab, S. and Reimer, A. and Reimer, O. and Renaud, M. and de los Reyes, R. and Richter, S. and Rieger, F. and Romoli, C. and Rowell, G. and Rudak, B. and Rulten, C. B. and Sahakian, V. and Saito, S. and Salek, D. and Sanchez, D. A. and Santangelo, A. and Sasaki, M. and Schlickeiser, R. and Schussler, F. and Schulz, A. and Schwanke, U. and Schwemmer, S. and Seglar-Arroyo, M. and Settimo, M. and Seyffert, A. S. and Shafi, N. and Shilon, I. and Simoni, R. and Sol, H. and Spanier, F. and Spengler, G. and Spies, F. and Stawarz, L. and Steenkamp, R. and Stegmann, Christian Michael and Stycz, K. and Sushch, I. and Takahashi, T. and Tavernet, J. -P. and Tavernier, T. and Taylor, A. M. and Terrier, R. and Tibaldo, L. and Tiziani, D. and Tluczykont, M. and Trichard, C. and Tsuji, N. and Tuffs, R. and Uchiyama, Y. and van der Walt, D. J. and van Eldik, C. and van Rensburg, C. and van Soelen, B. and Vasileiadis, G. and Veh, J. and Venter, C. and Viana, A. and Vincent, P. and Vink, J. and Voisin, F. and Voelk, H. J. and Vuillaume, T. and Wadiasingh, Z. and Wagner, S. J. and Wagner, P. and Wagner, R. M. and White, R. and Wierzcholska, A. and Willmann, P. and Woernlein, A. and Wouters, D. and Yang, R. and Zabalza, V. and Zaborov, D. and Zacharias, M. and Zanin, R. and Zdziarski, A. A. and Zech, Alraune and Zefi, F. and Ziegler, A. and Zywucka, N.}, title = {Systematic search for very-high-energy gamma-ray emission from bow shocks of runaway stars}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {612}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, organization = {HESS Collaboration}, issn = {1432-0746}, doi = {10.1051/0004-6361/201630151}, pages = {6}, year = {2018}, abstract = {Context. Runaway stars form bow shocks by ploughing through the interstellar medium at supersonic speeds and are promising sources of non-thermal emission of photons. One of these objects has been found to emit non-thermal radiation in the radio band. This triggered the development of theoretical models predicting non-thermal photons from radio up to very-high-energy (VHE, E >= 0.1 TeV) gamma rays. Subsequently, one bow shock was also detected in X-ray observations. However, the data did not allow discrimination between a hot thermal and a non-thermal origin. Further observations of different candidates at X-ray energies showed no evidence for emission at the position of the bow shocks either. A systematic search in the Fermi-LAT energy regime resulted in flux upper limits for 27 candidates listed in the E-BOSS catalogue. Aims. Here we perform the first systematic search for VHE gamma-ray emission from bow shocks of runaway stars. Methods. Using all available archival H.E.S.S. data we search for very-high-energy gamma-ray emission at the positions of bow shock candidates listed in the second E-BOSS catalogue release. Out of the 73 bow shock candidates in this catalogue, 32 have been observed with H.E.S.S. Results. None of the observed 32 bow shock candidates in this population study show significant emission in the H.E.S.S. energy range. Therefore, flux upper limits are calculated in five energy bins and the fraction of the kinetic wind power that is converted into VHE gamma rays is constrained. Conclusions. Emission from stellar bow shocks is not detected in the energy range between 0.14 and 18 TeV. The resulting upper limits constrain the level of VHE gamma-ray emission from these objects down to 0.1-1\% of the kinetic wind energy.}, language = {en} }