Search

Surveys in the Milky Way and Large Magellanic Cloud revealed that the majority of massive stars will interact with companions during their lives. However, knowledge of the binary properties of massive stars at low metallicity, which approaches the conditions of the Early Universe, remains sparse. We present the Binarity at LOw Metallicity (BLOeM) campaign - an ESO large programme designed to obtain 25 epochs of spectroscopy for 929 massive stars in the SMC - the lowest metallicity conditions in which multiplicity is probed to date (Z = 0.2 Zsun). BLOeM will provide (i) the binary fraction, (ii) the orbital configurations of systems with periods P < 3 yr, (iii) dormant OB+BH binaries, and (iv) a legacy database of physical parameters of massive stars at low metallicity. The stars are observed with the LR02 setup of the giraffe instrument of the Very Large Telescope (3960-4570A, resolving power R=6200; typical signal-to-noise ratio S/N=70-100). This paper utilises the first 9 epochs obtained over a three-month time. We describe the survey and data reduction, perform a spectral classification of the stacked spectra, and construct a Hertzsprung-Russell diagram of the sample via spectral-type and photometric calibrations. The sample covers spectral types from O4 to F5, spanning the effective temperature and luminosity ranges 6.5

Context. B supergiants (BSGs) represent an important connection between the main sequence and more extreme evolutionary stages of massive stars. Additionally, lying toward the cool end of the hot star regime, determining their wind properties is crucial to constrain the evolution and feedback of massive stars as, for instance, they might manifest the bi-stability jump phenomenon. Aims. We undertake a detailed analysis of a representative sample of 18 Small Magellanic Cloud (SMC) BSGs within the ULLYSES and XShootU datasets. Our UV and optical analysis spans BSGs from B0 to B8 - covering the bi-stability jump region. We aim to evaluate their evolutionary status and verify what their wind properties say about the bi-stability jump in a low-metallicity environment. Methods. We used the CMFGEN to model the spectra and photometry (from UV to infrared) of our sample. We compare our results with different evolutionary models, with previous determinations in the literature of OB stars, and with diverging mass-loss recipes at the bi-stability jump. Additionally, we provide the first BSG models in the SMC including X-rays. Results. (i) Within a single-stellar evolution framework, the evolutionary status of early BSGs seem less clear than that of late BSGs, which agree with H-shell burning models. (ii) UV analysis shows evidence that BSGs contain X-rays in their atmospheres, for which we provide constraints. In general, we find higher X-ray luminosity (close to the standard log(L_X/L) ~ -7) for early BSGs. For cooler BSGs, lower values are preferred, log(L_X/L) ~ -8.5. (iii) The obtained mass-loss rates suggest neither a jump nor a monotonic decrease with temperature. Instead, a rather constant trend is observed, which is at odds with the increase found for Galactic BSGs. (iv) The wind velocity behavior with temperature shows a sharp drop at ~19 kK, similar to what is observed for Galactic BSGs., Comment: 33 pages (23+10)+(22 at Zenodo), 34 figures (21+13)+(21 at Zenodo), 7 tables (3+4)+(1 at Zenodo), accepted for publication

We constrain wind parameters of a sample of 18 O-type stars in the LMC, through analysis with stellar atmosphere and wind models including the effects of optically thick clumping. This allows us to determine the most accurate spectroscopic mass-loss and wind structure properties of massive stars at sub-solar metallicity to date and gain insight into the impact of metallicity on massive stellar winds. Combining high signal to noise (S/N) ratio spectroscopy in the UV and optical gives us access to diagnostics of multiple different physical processes in the stellar wind. We produce synthetic spectra using the stellar atmosphere modelling code FASTWIND, and reproduce the observed spectra using a genetic algorithm based fitting technique. We empirically constrain 15 physical parameters associated with the stellar and wind properties, including temperature, surface gravity, surface abundances, rotation, macroturbulence and wind parameters. We find, on average, mass-loss rates a factor of 4-5 lower than those predicted by Vink et al. 2001, in good agreement with predictions from Bjorklund et al. 2021, and the best agreement with those from Krticka et al. 2018. In the 'weak-wind' regime we find mass-loss rates orders of magnitude below any theoretical predictions. We find a positive correlation of clumping factors (fcl) with effective temperature with an average fcl = 14 +- 8 for the full sample. Above 38 kK an average 46 +- 24% of the wind velocity span is covered by clumps and the interclump density is 10-30% of the mean wind. Below an effective temperature of roughly 38 kK there must be additional light leakage for supergiants. For dwarf stars at low temperatures there is a statistical preference for very low clump velocity spans, however it is unclear if this can be physically motivated as there are no clearly observable wind signatures in UV diagnostics., Comment: 12 pages with 6 figures, plus 16 pages and 18 figures in appendix. Accepted for publication in A&A

Gravitational waves from black-hole merging events have revealed a population of extra-galactic BHs residing in short-period binaries with masses that are higher than expected based on most stellar evolution models - and also higher than known stellar-origin black holes in our Galaxy. It has been proposed that those high-mass BHs are the remnants of massive metal-poor stars. Gaia astrometry is expected to uncover many Galactic wide-binary systems containing dormant BHs, which may not have been detected before. The study of this population will provide new information on the BH-mass distribution in binaries and shed light on their formation mechanisms and progenitors. As part of the validation efforts in preparation for the fourth Gaia data release (DR4), we analysed the preliminary astrometric binary solutions, obtained by the Gaia Non-Single Star pipeline, to verify their significance and to minimise false-detection rates in high-mass-function orbital solutions. The astrometric binary solution of one source, Gaia BH3, implies the presence of a 32.70 \pm 0.82 M\odot BH in a binary system with a period of 11.6 yr. Gaia radial velocities independently validate the astrometric orbit. Broad-band photometric and spectroscopic data show that the visible component is an old, very metal-poor giant of the Galactic halo, at a distance of 590 pc. The BH in the Gaia BH3 system is more massive than any other Galactic stellar-origin BH known thus far. The low metallicity of the star companion supports the scenario that metal-poor massive stars are progenitors of the high-mass BHs detected by gravitational-wave telescopes. The Galactic orbit of the system and its metallicity indicate that it might belong to the Sequoia halo substructure. Alternatively, and more plausibly, it could belong to the ED-2 stream, which likely originated from a globular cluster that had been disrupted by the Milky Way., Comment: 23 pages, accepted fro publication in A&A Letters. New version with small fixes

Massive stars (those larger than 8 solar masses at formation) have radiative envelopes that cannot sustain a dynamo, the mechanism that produces magnetic fields in lower-mass stars. Despite this, approximately 7\% of massive stars have observed magnetic fields, the origin of which is debated. We used multi-epoch interferometric and spectroscopic observations to characterize HD 148937, a binary system of two massive stars. We found that only one star is magnetic and that it appears younger than its companion. The system properties and a surrounding bipolar nebula can be reproduced with a model in which two stars merged (in a previous triple system) to produce the magnetic massive star. Our results provide observational evidence that magnetic fields form in at least some massive stars through stellar mergers., Comment: Full paper with supplementary materials. 59 pages, 18 figures

Using the medium resolution spectrograph X-shooter, spectra of 235 OB and Wolf-Rayet (WR) stars in sub-solar metallicity environments have been secured. [...]This second paper focuses on the optical observations of 232 Magellanic Clouds targets. It describes the uniform reduction of the UVB (300 - 560 nm) and VIS (550 - 1020 nm) XShootU data as well as the preparation of advanced data products [...] . The data reduction of the raw data is based on the ESO CPL X-shooter pipeline. We paid particular attention to the determination of the response curves [...] We implemented slit-loss correction, absolute flux calibration, (semi-)automatic rectification to the continuum, and a correction for telluric lines. The spectra of individual epochs were corrected for the barycentric motion, re-sampled and co-added, and the spectra from the two arms were merged into a single flux calibrated spectrum covering the entire optical range with maximum signal-to-noise ratio. [...] We provide three types of data products: (i) two-dimensional spectra for each UVB and VIS exposure; (ii) one-dimensional UVB and VIS spectra before and after response-correction, as well as after applying various processing, including absolute flux calibration, telluric removal, normalisation and barycentric correction; and (iii) co-added flux-calibrated and rectified spectra over the full optical range, for which all available XShootU exposures were combined. For many of the targets, the final signal-to-noise ratio per resolution element is above 200 in both the UVB and the VIS co-added spectra. The reduced data and advanced scientific data products will be made available to the community upon publication of this paper. [...], Comment: 22 pages, 19 figures; accepted for publication in Astronomy & Astrophysics. Links to online tables and databases will be included upon publication

Gaia's readout window strategy is challenged by very dense fields in the sky. Therefore, in addition to standard Gaia observations, full Sky Mapper (SM) images were recorded for nine selected regions in the sky. A new software pipeline exploits these Service Interface Function (SIF) images of crowded fields (CFs), making use of the availability of the full two-dimensional (2D) information. This new pipeline produced half a million additional Gaia sources in the region of the omega Centauri ($\omega$ Cen) cluster, which are published with this Focused Product Release. We discuss the dedicated SIF CF data reduction pipeline, validate its data products, and introduce their Gaia archive table. Our aim is to improve the completeness of the {\it Gaia} source inventory in a very dense region in the sky, $\omega$ Cen. An adapted version of {\it Gaia}'s Source Detection and Image Parameter Determination software located sources in the 2D SIF CF images. We validated the results by comparing them to the public {\it Gaia} DR3 catalogue and external Hubble Space Telescope data. With this Focused Product Release, 526\,587 new sources have been added to the {\it Gaia} catalogue in $\omega$ Cen. Apart from positions and brightnesses, the additional catalogue contains parallaxes and proper motions, but no meaningful colour information. While SIF CF source parameters generally have a lower precision than nominal {\it Gaia} sources, in the cluster centre they increase the depth of the combined catalogue by three magnitudes and improve the source density by a factor of ten. This first SIF CF data publication already adds great value to the {\it Gaia} catalogue. It demonstrates what to expect for the fourth {\it Gaia} catalogue, which will contain additional sources for all nine SIF CF regions.

Context. Strongly lensed quasars are fundamental sources for cosmology. The Gaia space mission covers the entire sky with the unprecedented resolution of $0.18$" in the optical, making it an ideal instrument to search for gravitational lenses down to the limiting magnitude of 21. Nevertheless, the previous Gaia Data Releases are known to be incomplete for small angular separations such as those expected for most lenses. Aims. We present the Data Processing and Analysis Consortium GravLens pipeline, which was built to analyse all Gaia detections around quasars and to cluster them into sources, thus producing a catalogue of secondary sources around each quasar. We analysed the resulting catalogue to produce scores that indicate source configurations that are compatible with strongly lensed quasars. Methods. GravLens uses the DBSCAN unsupervised clustering algorithm to detect sources around quasars. The resulting catalogue of multiplets is then analysed with several methods to identify potential gravitational lenses. We developed and applied an outlier scoring method, a comparison between the average BP and RP spectra of the components, and we also used an extremely randomised tree algorithm. These methods produce scores to identify the most probable configurations and to establish a list of lens candidates. Results. We analysed the environment of 3 760 032 quasars. A total of 4 760 920 sources, including the quasars, were found within 6" of the quasar positions. This list is given in the Gaia archive. In 87\% of cases, the quasar remains a single source, and in 501 385 cases neighbouring sources were detected. We propose a list of 381 lensed candidates, of which we identified 49 as the most promising. Beyond these candidates, the associate tables in this Focused Product Release allow the entire community to explore the unique Gaia data for strong lensing studies further., Comment: 35 pages, 60 figures, accepted for publication by Astronomy and Astrophysics

The third Gaia Data Release (DR3) provided photometric time series of more than 2 million long-period variable (LPV) candidates. Anticipating the publication of full radial-velocity (RV) in DR4, this Focused Product Release (FPR) provides RV time series for a selection of LPVs with high-quality observations. We describe the production and content of the Gaia catalog of LPV RV time series, and the methods used to compute variability parameters published in the Gaia FPR. Starting from the DR3 LPVs catalog, we applied filters to construct a sample of sources with high-quality RV measurements. We modeled their RV and photometric time series to derive their periods and amplitudes, and further refined the sample by requiring compatibility between the RV period and at least one of the $G$, $G_{\rm BP}$, or $G_{\rm RP}$ photometric periods. The catalog includes RV time series and variability parameters for 9\,614 sources in the magnitude range $6\lesssim G/{\rm mag}\lesssim 14$, including a flagged top-quality subsample of 6\,093 stars whose RV periods are fully compatible with the values derived from the $G$, $G_{\rm BP}$, and $G_{\rm RP}$ photometric time series. The RV time series contain a mean of 24 measurements per source taken unevenly over a duration of about three years. We identify the great most sources (88%) as genuine LPVs, with about half of them showing a pulsation period and the other half displaying a long secondary period. The remaining 12% consists of candidate ellipsoidal binaries. Quality checks against RVs available in the literature show excellent agreement. We provide illustrative examples and cautionary remarks. The publication of RV time series for almost 10\,000 LPVs constitutes, by far, the largest such database available to date in the literature. The availability of simultaneous photometric measurements gives a unique added value to the Gaia catalog (abridged), Comment: 36 pages, 38 figures

The most massive stars known to date are R 136 a1, a2, a3, and c within the central cluster R 136a of the Tarantula nebula in the Large Magellanic Cloud (LMC), with reported masses in excess of 150-200$M_\odot$. However, the mass estimation of these stars relies on the assumption that they are single. We collected three epochs of spectroscopy for R 136 a1, a2, a3, and c with the Space Telescope Imaging Spectrograph (STIS) of the Hubble Space Telescope (HST) in the years 2020-2021 to probe potential radial-velocity (RV) variations. We combine these epochs with an additional HST/STIS observation taken in 2012. We use cross-correlation to quantify the RVs, and establish constraints on possible companions to these stars up to periods of ~10 yr. Objects are classified as binaries when the peak-to-peak RV shifts exceed 50 km/s, and when the RV shift is significant with respect to errors. R 136 a1, a2, and a3 do not satisfy the binary criteria and are thus classified as putatively single, although formal peak-to-peak RV variability on the level 40 km/s is noted for a3. Only R 136 c is classified as binary, in agreement with literature. We can generally rule out massive companions (M2 > ~50 Msun) to R 136 a1, a2, and a3 out to orbital periods of < 1 yr (separations < 5 au) at 95% confidence, or out to tens of years (separations < ~100 au) at 50% confidence. Highly eccentric binaries (e > ~0.9) or twin companions with similar spectra could evade detection down to shorter periods (> ~10 d), though their presence is not supported by the relative X-ray faintness of R 136 a1, a2, and a3. We derive a preliminary orbital solution with a 17.2 d period for the X-ray bright binary R 136 c, though more data are needed to conclusively derive its orbit. Our study supports a lower bound of 150-200 $M_\odot$ on the upper-mass limit at LMC metallicity, Comment: Accepted to A&A

The origin of initial rotation rates of stars, and how a star's surface rotational velocity changes during the evolution, either by internal angular momentum transport or due to interactions with a binary companion, remain open questions in stellar astrophysics. Here, we aim to derive the physical parameters and study the distribution of (projected) rotational velocities of B-type stars in the 35 Myr-old, massive cluster NGC 330 in the Small Magellanic Cloud. NGC 330 is in an age range where the number of post-interaction binaries is predicted to be high near the cluster turnoff (TO). We develop a simultaneous photometric and spectroscopic grid-fitting method adjusting atmosphere models on multi-band Hubble Space Telescope photometry and Multi Unit Spectroscopic Explorer spectroscopy. This allows us to homogeneously constrain the physical parameters of over 250 B and Be stars, brighter than mF814W = 18.8 mag. The rotational velocities of Be stars in NGC 330 are significantly higher than the ones of B stars. The rotational velocities vary as a function of the star's position in the color-magnitude diagram, qualitatively following predictions of binary population synthesis. A comparison to younger clusters shows that stars in NGC 330 rotate more rapidly on average. The rotational velocities of the 35 Myr old population in NGC 330 quantitatively agree with predictions for a stellar population that underwent significant binary interactions: the bulk of the B stars could be single stars or primaries in pre-interaction binaries. The rapidly spinning Be stars could be mass and angular momentum gainers in previous interactions, while those Be stars close to the TO may be spun-up single stars. The slowly rotating, apparently single stars above the TO could be merger products. The different vsini-characteristics of NGC 330 compared to younger populations can be understood in this framework., Comment: 18 pages (incl. appendix), 15 figures, 3 tables, accepted for publication in A&A

Observations of individual massive stars, super-luminous supernovae, gamma-ray bursts, and gravitational-wave events involving spectacular black-hole mergers, indicate that the low-metallicity Universe is fundamentally different from our own Galaxy. Many transient phenomena will remain enigmatic until we achieve a firm understanding of the physics and evolution of massive stars at low metallicity (Z). The Hubble Space Telescope has devoted 500 orbits to observe 250 massive stars at low Z in the ultraviolet (UV) with the COS and STIS spectrographs under the ULLYSES program. The complementary ``X-Shooting ULLYSES'' (XShootU) project provides enhanced legacy value with high-quality optical and near-infrared spectra obtained with the wide-wavelength coverage X-shooter spectrograph at ESO's Very Large Telescope. We present an overview of the XShootU project, showing that combining ULLYSES UV and XShootU optical spectra is critical for the uniform determination of stellar parameters such as effective temperature, surface gravity, luminosity, and abundances, as well as wind properties such as mass-loss rates in function of Z. As uncertainties in stellar and wind parameters percolate into many adjacent areas of Astrophysics, the data and modelling of the XShootU project is expected to be a game-changer for our physical understanding of massive stars at low Z. To be able to confidently interpret James Webb Space Telescope (JWST) spectra of the first stellar generations, the individual spectra of low Z stars need to be understood, which is exactly where XShootU can deliver., Comment: Accepted in A&A - 35 Pages, 12 Figures, 4 Tables, 2 Large Tables

The Gaia-ESO Survey is a public spectroscopic survey that has targeted $\gtrsim10^5$ stars covering all major components of the Milky Way from the end of 2011 to 2018, delivering its public final release in May 2022. Unlike other spectroscopic surveys, Gaia-ESO is the only survey that observed stars across all spectral types with dedicated, specialised analyses: from O ($T_\mathrm{eff} \sim 30,000-52,000$~K) all the way to K-M ($\gtrsim$3,500~K). The physics throughout these stellar regimes varies significantly, which has previously prohibited any detailed comparisons between stars of significantly different type. In the final data release (internal data release 6) of the Gaia-ESO Survey, we provide the final database containing a large number of products such as radial velocities, stellar parameters and elemental abundances, rotational velocity, and also, e.g., activity and accretion indicators in young stars and membership probability in star clusters for more than 114,000 stars. The spectral analysis is coordinated by a number of Working Groups (WGs) within the Survey, which specialise in the various stellar samples. Common targets are analysed across WGs to allow for comparisons (and calibrations) amongst instrumental setups and spectral types. Here we describe the procedures employed to ensure all Survey results are placed on a common scale to arrive at a single set of recommended results for all Survey collaborators to use. We also present some general quality and consistency checks performed over all Survey results., Comment: A&A accepted, minor revision, 36 pages, 38 figures

The winds of massive stars have an impact on stellar evolution and on the surrounding medium. The maximum speed reached by these outflows, the terminal wind speed, is a global wind parameter and an essential input for models of stellar atmospheres and feedback. With the arrival of the ULLYSES programme, a legacy UV spectroscopic survey with HST, we have the opportunity to quantify the wind speeds of massive stars at sub-solar metallicity (in the Large and Small Magellanic Clouds, 0.5Z and 0.2Z) at an unprecedented scale. We empirically quantify the wind speeds of a large sample of OB stars, including supergiants, giants, and dwarfs at sub-solar metallicity. Using these measurements, we investigate trends of terminal wind speed with a number of fundamental stellar parameters, namely effective temperature, metallicity, and surface escape velocity. We empirically determined the terminal wind speed for a sample of 149 OB stars in the Magellanic Clouds either by directly measuring the maximum velocity shift of the absorption component of the Civ 1548-1550 line profile, or by fitting synthetic spectra produced using the Sobolev with exact integration method. Stellar parameters were either collected from the literature, obtained using spectral-type calibrations, or predicted from evolutionary models. We find strong trends of terminal wind speed with effective temperature and surface escape speed when the wind is strong enough to cause a saturated P Cygni profile in Civ 1548-1550. We find evidence for a metallicity dependence on the terminal wind speed proportional to Z^0.22+-0.03 when we compared our results to previous Galactic studies. Our results suggest that effective temperature rather than surface escape speed should be used as a straightforward empirical prediction of terminal wind speed and that the observed metallicity dependence is steeper than suggested by earlier works., Comment: 21 pages, 16 figures, 8 tables. Accepted in A&A

Massive star feedback affects the evolution of galaxies, where the most massive stars may have the largest impact. The majority of massive stars are born as members of close binary systems. Here, we investigate detailed evolutionary models of very massive binaries (30$\dots$90$M_{\odot}$) with Large Magellanic Cloud (LMC) metallicity. We identify four effects defying the conventional knowledge of binary evolution, which are all related to the proximity of the models to the Eddington limit. We find that the majority of systems undergo mass transfer during core hydrogen burning. During the ensuing nuclear timescale evolution, many mass donors remain more massive than their companions (``reverse Algols''), and nuclear timescale mass transfer may be interrupted or absent altogether. Furthermore, due to the elevated luminosity-to-mass ratio, many of the core-hydrogen burning donors may develop Wolf-Rayet-type winds, at luminosities where single stars would not. We identify observational counterparts of very massive reverse Algol binaries in the LMC, and discuss their contribution to the observed hydrogen-rich Wolf-Rayet stars. We argue that an understanding of very massive Algol systems is key to predicting the advanced evolution of very massive binaries, including their ability to evolve into observable gravitational wave sources., Comment: 15 pages, 5 figures, proposed for acceptance in A&A. Animation at https://zenodo.org/record/4068428#.Y-SrOtLMKV4

The Gaia-ESO survey sample of massive OB stars in the Carina Nebula consists of 234 stars. The addition of brighter sources from the Galactic O-Star Spectroscopic Survey and additional sources from the literature allows us to create the most complete census of massive OB stars done so far in the region. It contains a total of 316 stars, being 18 of them in the background and four in the foreground. Of the 294 stellar systems in Car OB1, 74 are of O type, 214 are of non-supergiant B type and 6 are of WR or non-O supergiant (II to Ia) spectral class. We identify 20 spectroscopic binary systems with an O-star primary, of which 6 are reported for the first time, and another 18 with a B-star primary, of which 13 are new detections. The average observed double-lined binary fraction of O-type stars in the surveyed region is 0.35, which represents a lower limit. We find a good correlation between the spectroscopic n-qualifier and the projected rotational velocity of the stars. The fraction of candidate runaways among the stars with and without the n-qualifier is 4.4% and 2.4%, respectively, although non resolved double-lined binaries can be contaminating the fast rotators sample., Comment: Accepted for publication in A&A, 62 pages, 8 figures

Two main scenarios have been proposed for origin of massive runaway stars -- dynamical ejection or release from a binary at the first core collapse -- but their relative contribution remains debated. Using two large spectroscopic campaigns towards massive stars in 30 Doradus, we aim to provide observational constraints on the properties of the O-type runaway population in the most massive active star-forming region in the Local group. We use RV measurements of the O-type star populations in 30 Doradus obtained by the VLT-FLAMES Tarantula Survey and the Tarantula Massive Binary Monitoring to identify single and binary O-type runaways. We discuss their rotational properties and qualitatively compare observations with expectations of ejection scenarios. We identify 23 single and one binary O-type runaway objects, most of them outside the main star-forming regions in 30 Doradus. We find an overabundance of rapid rotators (vsini > 200km/s) among the runaway population, providing an explanation of the overabundance of rapidly rotating stars in the 30 Doradus field. Considerations of the projected rotation rates and runaway line-of-sight (los) velocities reveal a conspicuous absence of rapidly rotating (vsini > 210k/ms), fast moving (v_{los} > 60km/s) runaways, and suggest the presence of two different populations of runaway stars: a population of rapidly-spinning but slowly moving runaways and a population of fast moving but slowly rotating ones. These are detected with a ratio close to 2:1 in our sample. We argue that slowly moving but rapidly spinning runaways result from binary ejections, while rapidly moving but slowly spinning runaways could result from dynamical ejections. Given that detection biases will more strongly impact the slow-moving population, our results suggest that the binary evolution scenario dominates the current massive runaway population in 30 Doradus., Comment: Accepted for publication in A&A Letters; 9 pages, 5 figures

Recent evolutionary computations predict that a few percent of massive OB stars in binary systems should have a dormant BH companion. Despite several reported X-ray quiet OB+BH systems over the last couple of years, finding them with certainty remains challenging. These have great importance as they can be gravitational wave (GW) source progenitors, and are landmark systems in constraining supernova kick physics. This work aims to characterise the hidden companions to the single-lined spectroscopic binaries (SB1s) in the B star population of the young open Galactic cluster NGC 6231 to find candidate systems for harbouring compact object companions. With the orbital solutions for each SB1 previously constrained, we applied Fourier spectral disentangling to multi-epoch optical VLT/FLAMES spectra of each target to extract a potential signature of a faint companion, and to identify newly disentangled double-lined spectroscopic binaries (SB2s). For targets where the disentangling does not reveal any signature of a stellar companion, we performed atmospheric and evolutionary modelling on the primary to obtain constraints on the unseen companion. Seven newly classified SB2 systems with mass ratios down to near 0.1 were identified. From the remaining targets, for which no faint companion could be extracted from the spectra, four are found to have companion masses in the predicted mass ranges of neutron stars (NSes) and BHes. Two of these have companion masses between 1 and 3.5 $M_{\odot}$, making them potential hosts of NSes (or lower mass main sequence stars). The other two are between 2.5 to 8 $M_{\odot}$ and 1.6 and 26 $M_{\odot}$, respectively, and so are identified as candidates for harbouring BH companions. However, unambiguous identification of these systems as X-ray quiet compact object harbouring binaries requires follow up observations., Comment: Submitted to A&A

The Gaia-ESO Public Spectroscopic Survey is an ambitious project designed to obtain astrophysical parameters and elemental abundances for 100,000 stars, including large representative samples of the stellar populations in the Galaxy, and a well-defined sample of 60 (plus 20 archive) open clusters. We provide internally consistent results calibrated on benchmark stars and star clusters, extending across a very wide range of abundances and ages. This provides a legacy data set of intrinsic value, and equally a large wide-ranging dataset that is of value for homogenisation of other and future stellar surveys and Gaia's astrophysical parameters. This article provides an overview of the survey methodology, the scientific aims, and the implementation, including a description of the data processing for the GIRAFFE spectra. A companion paper (arXiv:2206.02901) introduces the survey results. Gaia-ESO aspires to quantify both random and systematic contributions to measurement uncertainties. Thus all available spectroscopic analysis techniques are utilised, each spectrum being analysed by up to several different analysis pipelines, with considerable effort being made to homogenise and calibrate the resulting parameters. We describe here the sequence of activities up to delivery of processed data products to the ESO Science Archive Facility for open use. The Gaia-ESO Survey obtained 202,000 spectra of 115,000 stars using 340 allocated VLT nights between December 2011 and January 2018 from GIRAFFE and UVES. The full consistently reduced final data set of spectra was released through the ESO Science Archive Facility in late 2020, with the full astrophysical parameters sets following in 2022., Comment: 38 pages. A&A in press

The spectrum of the Wolf-Rayet (WR) star WR 63 contains spectral lines of two different O stars that show regular radial velocity (RV) variations with amplitudes of ~160 and ~225 km/s on a ~4.03 d period. The light-curve shows two narrow eclipses that are 0.2 mag deep on the same period as the RV changes. On the other hand, our data show no significant RV variations for the WR spectral lines. Those findings are compatible with WR 63 being a triple system composed of two non-interacting late O stars orbiting a WR star on a period larger than 1000 days. The amplitude of the WR spectral line-profile variability reaches 7-8% of the line intensity and seems related to a 0.04 mag periodic photometric variation. Large wind density structures are a possible origin of this variability, but our data are not sufficient to verify this. Our analysis shows that, should the three stars be bound, they would be coeval with an age of about 5.9+/-1.4 Myrs. The distance to the O stars is estimated to be 3.4+/-0.5 kpc. Their dynamical masses are 14.3+/-0.1 and 10.3+/-0.1 M_sol. Using rotating, single star evolutionary tracks, we estimate their initial masses to be 18+/-2 and 16+/-2 M_sol for the primary and the secondary, respectively. Regular spectral monitoring is required in the future to detect RV variations of the WR star that would prove that it is gravitationally bound to the close O+OB system and to determine its mass., Comment: 11 pages, 9 figures, 2 tables

We present a spectroscopic analysis of the GIRAFFE and UVES data collected by the Gaia-ESO survey for the young open cluster NGC 3293. Archive spectra from the same instruments obtained in the framework of the `VLT-FLAMES survey of massive stars' are also analysed. Atmospheric parameters, non-LTE chemical abundances for six elements, or variability information are reported for a total of about 160 B stars spanning a wide range in terms of spectral types (B1 to B9.5) and rotation rate (up to 350 km/s). We take advantage of the multi-epoch observations to detect several binary systems or intrinsically line-profile variables. A deconvolution algorithm is used to infer the current, true (deprojected) rotational velocity distribution. We find a broad, Gaussian-like distribution peaking around 200-250 km/s. Although some stars populate the high-velocity tail, most stars in the cluster appear to rotate far from critical. We discuss the chemical properties of the cluster, including the low occurrence of abundance peculiarities in the late B stars and the paucity of objects showing CN-cycle burning products at their surface. We argue that the former result can largely be explained by the inhibition of diffusion effects because of fast rotation, while the latter is generally in accord with the predictions of single-star evolutionary models under the assumption of a wide range of initial spin rates at the onset of main-sequence evolution. However, we find some evidence for a less efficient mixing in two quite rapidly rotating stars that are among the most massive objects in our sample. Finally, we obtain a cluster age of ~20 Myrs through a detailed, star-to-star correction of our results for the effect of stellar rotation. This is significantly older than previous estimates from turn-off fitting that fully relied on classical, non-rotating isochrones. [abridged], Comment: 29 pages, 24 figures, accepted for publication in A&A

Aims. To search for these rare objects, we study 32 Galactic O-type stars that were reported as SB1s in the literature. In our sample we include Cyg X-1, which is known to host an accreting stellar-mass BH, and HD 74194, a supergiant fast X-ray transient, in order to validate our methodology. The final goal is to characterise the nature of the unseen companions to determine if they are MS stars, stripped helium stars, triples, or compact objects such as neutron stars or stellar-mass BHs. Methods. After measuring radial velocities and deriving orbital solutions for all the systems in our sample, we performed spectral disentangling to extract putative signatures of faint secondary companions from the composite spectra. We derived stellar parameters for the visible stars and estimated the mass ranges of the secondary stars using the binary mass function. Variability observed in the photometric TESS light curves was also searched for indications of the presence of putative companions, degenerate or not. Results. In 17 of the 32 systems reported as SB1s, we extract secondary signatures, down to mass ratios of ~0.15. For the 17 newly detected double-lined spectroscopic binaries (SB2s), we derive physical properties of the individual components and discuss why they have not been detected as such before. Among the remaining systems, we identify nine systems with possible NS or low-mass MS companions. For Cyg X-1 and HD 130298, we are not able to extract any signatures for the companions, and the minimum masses of their companions are estimated to be about 7Msun. Our simulations show that secondaries with such a mass should be detectable from our dataset, no matter their nature: MS stars, stripped helium stars or even triples. While this is expected for Cyg X-1, confirming our methodology, our simulations also strongly suggest that HD 130298 could be another candidate to host a stellar-mass BH., Comment: 20 pages (31 pages of Appendix), 19 figures, A&A in press. The abstract is shorter than in the paper, and some figures, and tables have been cut in the paper

While the number of stellar-mass black holes detected in X-rays or as gravitational wave sources is steadily increasing, the known population remains orders of magnitude smaller than predicted by stellar evolution theory. A significant fraction of stellar-mass black holes is expected to hide in X-ray-quiet binaries where they are paired with a "normal" star. Although a handful of such quiescent black hole candidates have been proposed, the majority have been challenged by follow-up investigations. A confusion that emerged recently concerns binary systems that appear to contain a normal B-type star with an unseen companion, believed to be a black hole. On closer inspection, some of these seemingly normal B-type stars instead turn out to be stars stripped of most of their mass through an interaction with their binary companion, which in at least two cases is a rapidly rotating star rather than a compact object. These contaminants in the search for quiescent black holes are themselves extremely interesting objects as they represent a rare phase of binary evolution, and should be given special attention when searching for binaries hosting black holes in large spectroscopic studies., Comment: 7 pages, 5 figures

In the last 15 years different ground-based spectroscopic surveys have been started (and completed) with the general aim of delivering stellar parameters and elemental abundances for large samples of Galactic stars, complementing Gaia astrometry. Among those surveys, the Gaia-ESO Public Spectroscopic Survey (GES), the only one performed on a 8m class telescope, was designed to target 100,000 stars using FLAMES on the ESO VLT (both Giraffe and UVES spectrographs), covering all the Milky Way populations, with a special focus on open star clusters. This article provides an overview of the survey implementation (observations, data quality, analysis and its success, data products, and releases), of the open cluster survey, of the science results and potential, and of the survey legacy. A companion article (Gilmore et al.) reviews the overall survey motivation, strategy, Giraffe pipeline data reduction, organisation, and workflow. The GES has determined homogeneous good-quality radial velocities and stellar parameters for a large fraction of its more than 110,000 unique target stars. Elemental abundances were derived for up to 31 elements for targets observed with UVES. Lithium abundances are delivered for about 1/3 of the sample. The analysis and homogenisation strategies have proven to be successful; several science topics have been addressed by the Gaia-ESO consortium and the community, with many highlight results achieved. The final catalogue has been released through the ESO archive at the end of May 2022, including the complete set of advanced data products. In addition to these results, the Gaia-ESO Survey will leave a very important legacy, for several aspects and for many years to come., Comment: Accepted for publication in Astronomy and Astrophysics. 30 pages, 30 figures, 4 tables

Two scenarios have been proposed to match the existing observational constraints of the object HR 6819. The system could consist of a close inner B-type giant plus a black hole (BH) binary with an additional Be companion in a wide orbit. Alternatively, it could be a binary composed of a stripped B star and a Be star in a close orbit. Either scenario makes HR 6819 a cornerstone object as the stellar BH closest to Earth, or as an example of an important transitional, non-equilibrium phase for Be stars with solid evidence for its nature. We aimed to distinguish between the two scenarios for HR 6819. Both models predict two luminous stars but with very different angular separations and orbital motions. Therefore, the presence of bright sources in the 1-100 milliarcsec (mas) regime is a key diagnostic for determining the nature of the HR 6819 system. We obtained new high-angular resolution data with VLT/MUSE and VLTI/GRAVITY of HR 6819. The MUSE data are sensitive to bright companions at large scales, whilst the interferometric GRAVITY data are sensitive down to separations on mas scales and large magnitude differences. The MUSE observations reveal no bright companion at large separations and the GRAVITY observations indicate the presence of a stellar companion at an angular separation of ~1.2 mas that moves on the plane of the sky over a timescale compatible with the known spectroscopic 40-day period. We conclude that HR 6819 is a binary system and that no BH is present in the system. The unique nature of HR 6819, and its proximity to Earth make it an ideal system for quantitatively characterising the immediate outcome of binary interaction and probing how Be stars form., Comment: 16 pages, 11 figures, subject of ESO press release eso2204

The Gaia-ESO Survey (GES) is a large public spectroscopic survey that has collected, over a period of 6 years, spectra of ~ 10^5 stars. This survey provides not only the reduced spectra, but also the stellar parameters and abundances resulting from the analysis of the spectra. The GES dataflow is organised in 19 working groups. Working group 13 (WG13) is responsible for the spectral analysis of the hottest stars (O, B and A type, with a formal cut-off of Teff > 7000 K) that were observed as part of GES. We present the procedures and techniques that have been applied to the reduced spectra, in order to determine the stellar parameters and abundances of these stars. The procedure used is similar to that of other working groups in GES. A number of groups (called `Nodes') each independently analyse the spectra, using their state-of-the-art techniques and codes. Specific for the analysis in WG13 is the large temperature range that is covered (Teff = 7000 - 50,000 K), requiring the use of different analysis codes. Most Nodes can therefore only handle part of the data. Quality checks are applied to the results of these Nodes by comparing them to benchmark stars, and by comparing them one to another. For each star the Node values are then homogenised into a single result: the recommended parameters and abundances. Eight Nodes each analysed (part of) the data. In total 17,693 spectra of 6462 stars were analysed, most of them in 37 open star clusters. The homogenisation led to stellar parameters for 5584 stars. Abundances were determined for a more limited number of stars. Elements studied are He, C, N, O, Ne, Mg, Al, Si and Sc. Abundances for at least one of those elements were determined for 292 stars. The hot-star data analysed here, as well as the Gaia-ESO Survey data in general, will be of considerable use in future studies of stellar evolution and open clusters., Comment: 26 pages, 14 figures, accepted for publication in Astronomy & Astrophysics; language-edited version; two appendices merged

About a quarter of massive binary stars undergo mass transfer while both stars burn hydrogen at their cores, first on the thermal and then on the nuclear timescale. The nuclear timescale mass transfer leads to observable counterparts: the semi-detached so-called massive Algol binaries. However, comprehensive model predictions for these systems are sparse. We study them using a large grid of ~10,000 detailed binary evolution models calculated with the stellar evolution code MESA, covering initial donor masses between 10-40 M$_{\odot}$ and initial orbital periods above 1.4 d, at a metallicity suitable for the Large Magellanic Cloud (LMC). Our models imply ~30, or ~3% of the ~1,000 core hydrogen burning O-star binaries in the LMC to be currently in the semi-detached phase. Our donor models are up to 25-times more luminous than single stars of identical mass and effective temperature, which agrees with the observed Algols. A comparison of our models with the observed orbital periods and mass ratios implies rather conservative mass transfer in some systems, while very inefficient one in others. This is generally well reproduced by our spin-dependent mass transfer algorithm, except for the lowest considered masses. The observations reflect the slow increase of the surface nitrogen enrichment of the donors during the semi-detached phase all the way to CNO equilibrium. We also investigate the properties of our models after core hydrogen depletion of the donor star, when these models correspond to Wolf-Rayet/helium+OB star binaries. A dedicated spectroscopic survey of massive Algol systems may allow to derive the dependence of the efficiency of thermal timescale mass transfer on the binary parameters, as well as the efficiency of semiconvective mixing in the stellar interior. This would be a crucial step towards reliable binary models up to the formation of supernovae and compact objects., Comment: Accepted for publication in A&A (December 6, 2021)

Luminous blue variable stars (LBVs) are of great interest in massive-star evolution as they experience very high mass-loss episodes within short periods of time. HR Car is a famous member of this class in the Galaxy. It has a large circumstellar nebula and has also been confirmed as being in a binary system. One means of gaining information about the evolutionary status and physical nature of LBVs is studying their environments. We investigated the stellar content within ~100 pc of HR Car and also its circumstellar nebula. Very Large Telescope (VLT) Multi Unit Spectroscopic Explorer (MUSE) observations of a 2'x2' region around the star highlight the incompleteness of stellar classification for stars with magnitudes of V > 13 mag. Eight B0 to B9 stars have been identified which may lie in close spatial vicinity to HR Car. For a region with a radius of r =1.2 degree (~100 pc at a distance of 4.8 kpc) around HR Car, existing catalogs list several late O-type and early B-type stars, but only one early O-type star. Given the relatively low stellar and nebular masses in the HR Car system, no early O-type stars and only a few late O-type stars would be expected in association with HR Car. Instead, HR Car's location in a point vector diagram suggests that HR Car is not isolated, but is part of a moving group with a population of B-type stars in a spiral arm, and it has not received a strong kick from a supernova explosion of a companion star or a merger event. Potential binary evolution pathways for the HR Car system cannot be fully explored because of the unknown nature of the companion star. Furthermore, the MUSE observations reveal the presence of a fast outflow and "bullets" that have been ejected at intervals of about 400 years. These features may have been caused by recurrent mass transfer in the system., Comment: 11 pages, 5 figures, to be published in A&A

We investigate the impact of optically thick clumping on stellar wind diagnostics in O supergiants and constrain wind parameters associated with porosity in velocity space. This is the first time the effects of optically thick clumping have been investigated for a sample of massive hot stars, using models including a full optically thick clumping description. We re-analyse spectroscopic observations of a sample of eight O supergiants. Using a genetic algorithm wrapper around the NLTE atmosphere code FASTWIND we obtain simultaneous fits to optical and UV spectra and determine photospheric and wind properties and surface abundances. We provide empirical constraints on a number of wind parameters including the clumping factors, mass-loss rates and terminal wind velocities. Additionally, we establish the first systematic empirical constraints on velocity filling factors and interclump densities. These parameters describe clump distribution in velocity-space and density of the interclump medium in physical-space, respectively. We observe a mass-loss rate reduction of a factor of 3.6 compared to theoretical predictions from Vink et al. (2000), and mass-loss rates within a factor 1.4 of predictions from Bj\"orklund et al. (2021). We confirm that including optically thick clumping allows simultaneous fitting of recombination lines and resonance lines, including the unsaturated UV phosphorus lines (Pv 1118-1128), without reducing the phosphorus abundance. We find that, on average, half of the wind velocity field is covered by dense clumps. We also find that these clumps are 25 times denser than the average wind, and that the interclump medium is 3-10 times less dense than the mean wind. The former result agrees well with theoretical predictions, the latter suggests that lateral filling-in of radially compressed gas might be critical for setting the scale of the rarefied interclump matter., Comment: Accepted for publication in A&A; 35 pages, 20 figures

It is well known that massive O-stars are frequently (if not always) found in binary or higher-order multiple systems, but this fact has been less robustly investigated for the lower mass range of the massive stars, represented by the B-type stars. Obtaining the binary fraction and orbital parameter distributions of B-type stars is crucial to understand the impact of multiplicity on the archetypal progenitor of core-collapse supernovae as well as to properly investigate formation channels for gravitational wave progenitors. This work aims to characterise the multiplicity of the B-star population of the young open cluster NGC 6231 through multi-epoch optical spectroscopy of 80 B-type stars. We analyse 31 FLAMES/GIRAFFE observations of 80 B-type stars, monitoring their radial velocities (RVs) and performing a least-squares spectral analysis (Lomb-Scargle) to search for periodicity in those stars with statistically significant variability in their RVs. We constrain an observed spectroscopic binary fraction of $33\pm5$% for the B-type stars of NGC 6231, with a first order bias-correction giving a true spectroscopic binary fraction of $52\pm8$%. Out of 27 B-type binary candidates, we obtained orbital solutions for 20 systems: 15 single-lined (SB1) and 5 double-lined spectroscopic binaries (SB2s). We present these orbital solutions and the orbital parameter distributions associated with them. Our results indicate that Galactic B-type stars are less frequently found in binary systems than their more massive O-type counterparts, but their orbital properties generally resemble those of B- and O-type stars in both the Galaxy and Large Magellanic Cloud.

Context. Luminous Blue Variables (LBVs) are thought to be in a transitory phase between O stars on the main-sequence and the Wolf-Rayet stage. Recent studies suggest that they might be formed through binary interaction. Only a few are known in binary systems but their multiplicity fraction is uncertain. Aims. This study aims at deriving the binary fraction among the Galactic (confirmed and candidate) LBV population. We combine multi-epoch spectroscopy and long-baseline interferometry. Methods. We use cross-correlation to measure their radial velocities. We identify spectroscopic binaries through significant RV variability (larger than 35 km/s). We investigate the observational biases to establish the intrinsic binary fraction. We use CANDID to detect interferometric companions, derive their parameters and positions. Results. We derive an observed spectroscopic binary fraction of 26 %. Considering period and mass ratio ranges from Porb=1 to 1000 days, and q = 0.1-1.0, and a representative set of orbital parameter distributions, we find a bias-corrected binary fraction of 62%. From interferometry, we detect 14 companions out of 18 objects, providing a binary fraction of 78% at projected separations between 1 and 120 mas. From the derived primary diameters, and the distances of these objects, we measure for the first time the exact radii of Galactic LBVs to be between 100 and 650 Rsun, making unlikely to have short-period systems. Conclusions. This analysis shows that the binary fraction among the Galactic LBV population is large. If they form through single-star evolution, their orbit must be initially large. If they form through binary channel that implies that either massive stars in short binary systems must undergo a phase of fully non-conservative mass transfer to be able to sufficiently widen the orbit or that LBVs form through merging in initially binary or triple systems., Comment: 14 pages, 9 figures, 4 tables; Accepted for publication in A&A. (Abstract shortened to respect the arXiv limit of 1920 characters)

Direct dynamical mass measurements of stars with masses above 30 M${}_\odot$ are rare. This is the result of the low yield of the upper initial mass function and the limited number of such systems in eclipsing binaries. Long-period, double-lined spectroscopic binaries that are also resolved astrometrically offer an alternative for obtaining absolute masses of stellar objects. 9 Sgr is one such long-period, high-mass binary. Unfortunately, a large amount of tension exists between its total dynamical mass inferred from radial velocity measurements and that from astrometric data. We obtained the astrometric orbit from VLTI/PIONIER and VLTI/GRAVITY interferometric measurements. Using archival and new spectroscopy, we performed a grid-based spectral disentangling search to constrain the semi-amplitudes of the radial velocity curves. We computed atmospheric parameters and surface abundances by adjusting \textsc{fastwind} atmosphere models and we compared our results with evolutionary tracks computed with the Bonn Evolutionary Code (BEC). Grid spectral disentangling of 9 Sgr supports the presence of a 53 M${}_\odot$ primary and a 39 M${}_\odot$ secondary. Comparison with BEC evolutionary tracks shows the components of 9 Sgr are most likely coeval with an age of roughly 1 Myr. Our analysis clears up the contradiction between mass and orbital inclination estimates reported in previous studies. We detect the presence of significant CNO-processed material at the surface of the primary, suggesting enhanced internal mixing compared to currently implemented in the BEC models. The present measurements provide a high-quality high-mass anchor to validate stellar evolution models and to test the efficiency of internal mixing processes., Comment: 14 pages, 5 figures, 4 tables (+6 figures and 3 tables in appendices). Accepted for publication in Astronomy and Astrophysics

Observations of massive stars in young open clusters (< ~8 Myr) have shown that a majority of them are in binary systems, most of which will interact during their life. Populations of massive stars older than ~20 Myr allow us to probe the outcome of such interactions after many systems have experienced mass and angular momentum transfer. Using multi-epoch integral-field spectroscopy, we investigate the multiplicity properties of the massive-star population in NGC 330 (~40 Myr) in the Small Magellanic Cloud to search for imprints of stellar evolution on the multiplicity properties. From six epochs of VLT/MUSE observations supported by adaptive optics we extract spectra and measure radial velocities for stars brighter than F814W = 19. We identify single-lined spectroscopic binaries through significant RV variability as well as double-lined spectroscopic binaries, and quantify the observational biases for binary detection. The observed spectroscopic binary fraction is 13.2+/-2.0 %. Considering period and mass ratio ranges from log(P)=0.15-3.5, and q = 0.1-1.0, and a representative set of orbital parameter distributions, we find a bias-corrected close binary fraction of 34 +8 -7 %. This seems to decline for the fainter stars, which indicates either that the close binary fraction drops in the B-type domain, or that the period distribution becomes more heavily weighted towards longer orbital periods. Both fractions vary strongly in different regions of the color-magnitude diagram which probably reveals the imprint of the binary history of different groups of stars. We provide the first homogeneous RV study of a large sample of B-type stars at a low metallicity. The overall bias-corrected close binary fraction of B stars in NGC 330 is lower than the one reported for younger Galactic and LMC clusters. More data are needed to establish whether this result from an age or a metallicty effect., Comment: 19 page (incl. appendix), 10 figures, 4 tables, accepted for publication in A&A

R 144 is the visually brightest WR star in the Large Magellanic Cloud (LMC). R 144 was reported to be a binary, making it potentially the most massive binary thus observed. We perform a comprehensive spectral, photometric, orbital, and polarimetric analysis of R 144. R 144 is an eccentric (e=0.51) 74.2-d binary comprising two relatively evolved (age~2 Myr), H-rich WR stars. The hotter primary (WN5/6h, T=50 kK) and the cooler secondary (WN6/7h,T=45kK) have nearly equal masses. The combination of low rotation and H-depletion observed in the system is well reproduced by contemporary evolution models that include boosted mass-loss at the upper-mass end. The systemic velocity of R 144 and its relative isolation suggest that it was ejected as a runaway from the neighbouring R 136 cluster. The optical light-curve shows a clear orbital modulation that can be well explained as a combination of two processes: excess emission stemming from wind-wind collisions and double wind eclipses. Our light-curve model implies an orbital inclination of i=60.4+-1.5deg, resulting in accurately constrained dynamical masses of 74+-4 and 69+-4 Msun. Assuming that both binary components are core H-burning, these masses are difficult to reconcile with the derived luminosities (logL1,2 = 6.44, 6.39 [Lsun]), which correspond to evolutionary masses of the order of 110 and 100Msun, respectively. Taken at face value, our results imply that both stars have high classical Eddington factors of Gamma_e = 0.78+-0.1. If the stars are on the main sequence, their derived radii (~25Rsun) suggest that they are only slightly inflated, even at this high Eddington factor. Alternatively, the stars could be core-He burning, strongly inflated from the regular size of classical Wolf-Rayet stars (~1Rsun), a scenario that could help resolve the observed mass discrepancy., Comment: Accepted to A&A on the April 06 2021. 17 pages + 5 page appendix. Version 28.10.2021: corrected typo in Eq. (2); numerator and denominator were inversed on RHS (only typo; no impact on results)

Vigorous mass loss in the classical Wolf-Rayet (WR) phase is important for the late evolution and final fate of massive stars. We develop spherically symmetric time-dependent and steady-state hydrodynamical models of the radiation-driven wind outflows and associated mass loss from classical WR stars. The simulations are based on combining the opacities typically used in static stellar structure and evolution models with a simple parametrised form for the enhanced line-opacity expected within a supersonic outflow. Our simulations reveal high mass-loss rates initiated in deep and hot optically thick layers around T\approx 200kK. The resulting velocity structure is non-monotonic and can be separated into three phases: i) an initial acceleration to supersonic speeds ii) stagnation and even deceleration, and iii) an outer region of rapid re-acceleration. The characteristic structures seen in converged steady-state simulations agree well with the outflow properties of our time-dependent models. By directly comparing our dynamic simulations to corresponding hydrostatic models, we demonstrate explicitly that the need to invoke extra energy transport in convectively inefficient regions of stellar structure and evolution models is merely an artefact of enforcing a hydrostatic outer boundary. Moreover, the "dynamically inflated" inner regions of our simulations provide a natural explanation for the often-found mismatch between predicted hydrostatic WR radii and those inferred from spectroscopy. Finally, we contrast our simulations with alternative recent WR wind models based on co-moving frame radiative transfer for computing the radiation force. Since CMF transfer currently cannot handle non-monotonic velocity fields, the characteristic deceleration regions found here are avoided in such simulations by invoking an ad-hoc very high degree of clumping., Comment: 15 pages, 9 figures

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

The intriguing binary LS V +22 25 (LB-1) has drawn much attention following claims of it being a single-lined spectroscopic binary with a 79-day orbit comprising a B-type star and a ~70 Msun black hole. Recent analyses have implied that the visible primary star is a stripped He-rich star. However, the nature of the secondary, which was proposed to be a black hole, a neutron star, or a main sequence star, remains unknown. Based on 26 newly acquired spectroscopic observations, we perform an orbital analysis and spectral disentangling of LB-1 to elucidate the nature of the system. Our analysis reveals that LB-1 contains two components of comparable brightness in the optical. The narrow-lined primary, which we estimate to contribute ~55% in the optical, has spectral properties that suggest that it is a stripped star: it has a small spectroscopic mass (~1 Msun) for a B-type star and it is He- and N-rich. The "hidden" secondary, which contributes about 45% of the optical flux, is a rapidly rotating (vsini ~ 300 km/s) B3 V star with a decretion disk -- a Be star. Hence, LB-1 does not contain a compact object. Instead, it is a rare Be binary system consisting of a stripped star (the former mass donor) and a Be star rotating at near its critical velocity (the former mass accretor). This system is a clear example that binary interactions play a decisive role in the production of rapid stellar rotators and Be stars., Comment: Accepted to A&A. 4 pages (1 figure, 1 table); 3 pages of supporting material (7 figures and 1 table)

The recent gravitational wave measurements have demonstrated the existence of stellar mass black hole binaries. It is essential for our understanding of massive star evolution to identify the contribution of binary evolution to the formation of double black holes. A promising way to progress is investigating the progenitors of double black hole systems and comparing predictions with local massive star samples such as the population in 30 Doradus in the Large Magellanic Cloud (LMC). Methods. To this purpose, we analyse a large grid of detailed binary evolution models at LMC metallicity with initial primary masses between 10 and 40 Msun, and identify which model systems potentially evolve into a binary consisting of a black hole and a massive main sequence star. We then derive the observable properties of such systems, as well as peculiarities of the OB star component. We find that about 3% of the LMC late O and early B stars in binaries are expected to possess a black hole companion, when assuming stars with a final helium core mass above 6.6 M to form black holes. While the vast majority of them may be X-ray quiet, our models suggest that these may be identified in spectroscopic binaries, either by large amplitude radial velocity variations ( > 50 km s ) and simultaneous nitrogen surface enrichment, or through a moderate radial velocity ( > 10 km/s ) and simultaneously rapid rotation of the OB star. The predicted mass ratios are such that main sequence companions could be excluded in most cases. A comparison to the observed OB+WR binaries in the LMC, Be/X-ray binaries, and known massive BH binaries supports our conclusion. We expect spectroscopic observations to be able to test key assumptions in our models, with important implications for massive star evolution in general, and for the formation of double-black hole mergers in particular., Comment: 19 pages, 14 figures, Astronomy and Astrophysics, in press

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

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

A majority of massive stars are part of binary systems, a large fraction of which will inevitably interact during their lives. Binary-interaction products (BiPs), i.e. stars affected by such interaction, are expected to be commonly present in stellar populations. BiPs are thus a crucial ingredient in the understanding of stellar evolution. We aim to identify and characterize a statistically significant sample of BiPs by studying clusters of 10-40 Myr, an age at which binary population models predict the abundance of BiPs to be highest. One example of such a cluster is NGC 330 in the Small Magellanic Cloud. Using MUSE WFM-AO observations of NGC 330, we resolve the dense cluster core for the first time and are able to extract spectra of its entire massive star population. We develop an automated spectral classification scheme based on the equivalent widths of spectral lines in the red part of the spectrum. We characterize the massive star content of the core of NGC 330 which contains more than 200 B stars, 2 O stars, 6 A-type supergiants and 11 red supergiants. We find a lower limit on the Be star fraction of $32 \pm 3\%$ in the whole sample. It increases to at least $46 \pm 10\%$ when only considering stars brighter than V=17 mag. We estimate an age of the cluster core between 35 and 40 Myr and a total cluster mass of $88^{+17}_{-18} \times 10^3 M_{\odot}$. We find that the population in the cluster core is different than the population in the outskirts: while the stellar content in the core appears to be older than the stars in the outskirts, the Be star fraction and the observed binary fraction are significantly higher. Furthermore, we detect several BiP candidates that will be subject of future studies., Comment: 19 pages (incl. appendix), 21 figures, 3 tables, accepted for publication in A&A

The 30 Doradus (30 Dor) nebula in the Large Magellanic Cloud (LMC) is the brightest HII region in the Local Group and a prototype starburst similar to those found in high redshift galaxies. It is thus a stepping stone to understand the complex formation processes of stars in starburst regions across the Universe. Here, we have studied the formation history of massive stars in 30 Dor using masses and ages derived for 452 mainly OB stars from the spectroscopic VLT-FLAMES Tarantula Survey (VFTS). We find that stars of all ages and masses are scattered throughout 30 Dor. This is remarkable because it implies that massive stars either moved large distances or formed independently over the whole field of view in relative isolation. We find that both channels contribute to the 30 Dor massive star population. Massive star formation rapidly accelerated about 8 Myr ago, first forming stars in the field before giving birth to the stellar populations in NGC 2060 and NGC 2070. The R136 star cluster in NGC 2070 formed last and, since then, about 1 Myr ago, star formation seems to be diminished with some continuing in the surroundings of R136. Massive stars within a projected distance of 8 pc of R136 are not coeval but show an age range of up to 6 Myr. Our mass distributions are well populated up to $200\,\mathrm{M}_\odot$. The inferred IMF is shallower than a Salpeter-like IMF and appears to be the same across 30 Dor. By comparing our sample of stars to stellar models in the Hertzsprung-Russell diagram, we find evidence for missing physics in the models above $\log L/\mathrm{L}_\odot=6$ that is likely connected to enhanced wind mass loss for stars approaching the Eddington limit. [abridged], Comment: 20 pages (incl. appendix), 15 figures, 2 tables; accepted for publication in A&A

Many massive stars are part of binary or higher multiplicity systems. The present work focusses on two higher multiplicity systems: HD17505A and HD206267A. Determining the fundamental parameters of the components of the inner binary of these systems is mandatory to quantify the impact of binary or triple interactions on their evolution. We analysed high-resolution optical spectra to determine new orbital solutions of the inner binary systems. After subtracting the spectrum of the tertiary component, a spectral disentangling code was applied to reconstruct the individual spectra of the primary and secondary. We then analysed these spectra with the non-LTE model atmosphere code CMFGEN to establish the stellar parameters and the CNO abundances of these stars. The inner binaries of these systems have eccentric orbits with e ~ 0.13 despite their relatively short orbital periods of 8.6 and 3.7 days for HD17505Aa and HD206267Aa, respectively. Slight modifications of the CNO abundances are found in both components of each system. The components of HD17505Aa are both well inside their Roche lobe, whilst the primary of HD206267Aa nearly fills its Roche lobe around periastron passage. Whilst the rotation of the primary of HD206267Aa is in pseudo-synchronization with the orbital motion, the secondary displays a rotation rate that is higher. The CNO abundances and properties of HD17505Aa can be explained by single star evolutionary models accounting for the effects of rotation, suggesting that this system has not yet experienced binary interaction. The properties of HD206267Aa suggest that some intermittent binary interaction might have taken place during periastron passages, but is apparently not operating anymore., Comment: Accepted for publication in A&A

HD54662 is an O-type binary star belonging to the CMa OB1 association. Due to its long-period orbit, this system is an interesting target to test the adiabatic wind shock model. The goal is to improve our knowledge of the orbital and stellar parameters of HD54662 and to analyze its X-ray emission to test the theoretical scaling of the X-ray emission with orbital separation for adiabatic wind shocks. We applied a spectral disentangling code to optical spectra to determine the radial velocities and the individual spectra of each star. The individual spectra were analyzed using the CMFGEN model atmosphere code. We fitted two X-ray spectra using a Markov Chain Monte Carlo algorithm and compared them to the emission expected from adiabatic shocks. We determine an orbital period of 2103.4days, a low orbital eccentricity of 0.11, and a mass ratio m2/m1=0.84. Combined with the orbital inclination inferred in a previous astrometric study, we obtain surprisingly low masses of 9.7 and 8.2Msun. From the individual spectra, we infer O6.5 spectral types for both stars and a brightness ratio of l1/l2~2. The softness of the X-ray spectra, the very small variation of spectral parameters, and the comparison of the X-ray-to-bolometric luminosity ratio with the canonical value for O-type stars allow us to conclude that X-ray emission from the wind interaction region is quite low. We cannot confirm the runaway status previously attributed to HD54662 and we find no X-ray emission associated with the bow shock detected in the infrared. The lack of hard X-ray emission from the wind-shock region suggests that the mass-loss rates are lower than expected and/or that the pre-shock wind velocities are much lower than the terminal wind velocities. The bow shock associated with HD54662 possibly corresponds to a wind-blown arc created by the interaction of the stellar winds with the ionized gas of CMa OB1. (abridged), Comment: Manuscript has been accepted. A&A, in press

Aims. HD166734 is an eccentric eclipsing binary system composed of two supergiant O-type stars, orbiting with a 34.5-day period. In this rare configuration for such stars, the two objects mainly evolve independently, following single-star evolution so far. This system provides a chance to study the individual parameters of two supergiant massive stars and to derive their real masses. Methods. An intensive monitoring was dedicated to HD166734.We analyzed mid- and high-resolution optical spectra to constrain the orbital parameters of this system. We also studied its light curve for the first time, obtained in the VRI filters. Finally, we disentangled the spectra of the two stars and modeled them with the CMFGEN atmosphere code in order to determine the individual physical parameters. Results. HD166734 is a O7.5If+O9I(f) binary. We confirm its orbital period but we revise the other orbital parameters. In comparison to what we found in the literature, the system is more eccentric and, now, the hottest and the most luminous component is also the most massive one. The light curve exhibits only one eclipse and its analysis indicates an inclination of 63.0{\deg} $\pm$ 2.7{\deg}. The photometric analysis provides us with a good estimation of the luminosities of the stars, and therefore their exact positions in the Hertzsprung-Russell diagram. The evolutionary and the spectroscopic masses show good agreement with the dynamical masses of 39.5 Msun for the primary and 33.5 Msun for the secondary, within the uncertainties. The two components are both enriched in helium and in nitrogen and depleted in carbon. In addition, the primary also shows a depletion in oxygen. Their surface abundances are however not different from those derived from single supergiant stars, yielding, for both components, an evolution similar to that of single stars., Comment: 13 pages, 13 figures, A&A accepted

The role of mass and momentum exchanges in close massive binaries is very important in the subsequent evolution of the components. Such exchanges produce several observational signatures such as asynchronous rotation and altered chemical compositions, that remain after the stars detach again. We investigated these effects for the close O-star binary LSS 3074 (O4 f + O6-7 :(f):), which is a good candidate for a past Roche lobe overflow (RLOF) episode because of its very short orbital period, P = 2.185 days, and the luminosity classes of both components. We determined a new orbital solution for the system. We studied the photometric light curves to determine the inclination of the orbit and Roche lobe filling factors of both stars. Using phase-resolved spectroscopy, we performed the disentangling of the optical spectra of the two stars. We then analysed the reconstructed primary and secondary spectra with the CMFGEN model atmosphere code to determine stellar parameters, such as the effective temperatures and surface gravities, and to constrain the chemical composition of the components. We confirm the apparent low stellar masses and radii reported in previous studies. We also find a strong overabundance in nitrogen and a strong carbon and oxygen depletion in both primary and secondary atmospheres, together with a strong enrichment in helium of the primary star. We propose several possible evolutionary pathways through a RLOF process to explain the current parameters of the system. We confirm that the system is apparently in overcontact configuration and has lost a significant portion of its mass to its surroundings. We suggest that some of the discrepancies between the spectroscopic and photometric properties of LSS 3074 could stem from the impact of a strong radiation pressure of the primary.

The Gaia-ESO Survey (GES) is obtaining high quality spectra of $\sim$10$^{5}$ stars in our Galaxy, providing an homogeneous and unique overview of all the main components of the Milky Way, its formation history and the evolution of young, mature and ancient Galactic populations. Our group is in charge of the early-type massive stars that define the youngest population in the survey. In this contribution, we present the results of the quantitative spectroscopic analysis of O-type stars in the Carina Nebula within the Gaia-ESO Survey. For this aim, we have used FASTWIND and CMFGEN stellar atmosphere codes, providing stellar parameters for the current sample (GES data release iDR4)., Comment: 4 pages, 3 figures, to be published in Highlights on Spanish Astrophysics IX, Proceedings of the XII Scientific Meeting of the Spanish Astronomical Society held on July 18-22, 2016, in Bilbao, Spain. F. Figueras, A. S\'anchez-Lavega, S. P\'erez-Hoyos, A. Alonso, S. Arribas, C. Hern\'andez Monteagudo (eds.)

Aims. The magnetic star HD148937 is the only Galactic Of?p star surrounded by a nebula. The structure of this nebula is particularly complex and is composed, from the center out outwards, of a close bipolar ejecta nebula (NGC6164/5), an ellipsoidal wind-blown shell, and a spherically symmetric Stromgren sphere. The exact formation process of this nebula and its precise relation to the star's evolution remain unknown. Methods. We analyzed infrared Spitzer IRS and far-infrared Herschel/PACS observations of the NGC6164/5 nebula. The Herschel imaging allowed us to constrain the global morphology of the nebula. We also combined the infrared spectra with optical spectra of the central star to constrain its evolutionary status. We used these data to derive the abundances in the ejected material. To relate this information to the evolutionary status of the star, we also determined the fundamental parameters of HD148937 using the CMFGEN atmosphere code. Results. The H$\alpha$ image displays a bipolar or "8"-shaped ionized nebula, whilst the infrared images show dust to be more concentrated around the central object. We determine nebular abundance ratios of N/O = 1.06 close to the star, and N/O = 1.54 in the bright lobe constituting NGC6164. Interestingly, the parts of the nebula located further from HD148937 appear more enriched in stellar material than the part located closer to the star. Evolutionary tracks suggest that these ejecta have occured $\sim$1.2-1.3 and $\sim$0.6 Myrs ago, respectively. In addition, we derive abundances of argon for the nebula compatible with the solar values and we find a depletion of neon and sulfur. The combined analyses of the known kinematics and of the new abundances of the nebula suggest either a helical morphology for the nebula, possibly linked to the magnetic geometry, or the occurrence of a binary merger., Comment: 15 pages, 3 pages of appendix, 12 figures, A&A in press

Massive binary systems are important laboratories in which to probe the properties of massive stars and stellar physics in general. In this context, we analysed optical spectroscopy and photometry of the eccentric short-period early-type binary HD 152218 in the young open cluster NGC 6231. We reconstructed the spectra of the individual stars using a separating code. The individual spectra were then compared with synthetic spectra obtained with the CMFGEN model atmosphere code. We furthermore analysed the light curve of the binary and used it to constrain the orbital inclination and to derive absolute masses of 19.8 +/- 1.5 and 15.0 +/- 1.1 solar masses. Combining radial velocity measurements from over 60 years, we show that the system displays apsidal motion at a rate of (2.04^{+.23}_{-.24}) degree/year. Solving the Clairaut-Radau equation, we used stellar evolution models, obtained with the CLES code, to compute the internal structure constants and to evaluate the theoretically predicted rate of apsidal motion as a function of stellar age and primary mass. In this way, we determine an age of 5.8 +/- 0.6 Myr for HD 152218, which is towards the higher end of, but compatible with, the range of ages of the massive star population of NGC 6231 as determined from isochrone fitting., Comment: Accepted for publication in Astronomy & Astrophysics