Search

Wolf-Rayet stars (WRs) are thought to be the final evolutionary stages of massive stars and immediate progenitors of stellar-mass black holes. Their multiplicity forms an important anchor point in single and binary population models for predicting gravitational-wave progenitors. Recent spectroscopic campaigns have suggested incompatible multiplicity fractions and period distributions for N- and C-rich Galactic WRs (WNs and WCs) at short as well as long orbital periods, in contradiction with evolutionary model predictions. In this work, we employed infrared interferometry using the $K$-band instrument GRAVITY at the VLTI to investigate the multiplicity of WRs at long periods and explore the nature of their companions. We present a survey of 39 Galactic WRs, including 11 WN, 15 WC and 13 H-rich WN (WNh) stars. We detected wide companions with GRAVITY for only four stars: WR 48, WR 89, WR 93 and WR 115. Combining with spectroscopic studies, we arrived at multiplicity fractions of $f^{\rm WN}_{\rm obs} = 0.55\pm0.15$, $f^{\rm WC}_{\rm obs} = 0.40\pm0.13$ and $f^{\rm WNh}_{\rm obs} = 0.23\pm0.12$. In addition, we also found other features in the GRAVITY dataset such as (i) a diffuse extended component in over half the WR sample; (ii) five known spectroscopic binaries resolved in differential phase data and (iii) spatially resolved winds in four stars: WR 16, WR 31a, WR 78 and WR 110. Our survey reveals a lack of intermediate (few 100s d) and long- (few years to decades) period WR systems. The 200-d peak in the period distributions of WR+OB and BH+OB binaries predicted by Case B mass-transfer binary evolution models is not seen in our data. The rich companionship of their O-type progenitors in this separation range suggest that the WR progenitor stars expand and interact with their companions, most likely through unstable mass-transfer, resulting in either a short-period system or a merger., Comment: Submitted to A&A. Comments welcome

Significant progress has been made toward understanding the formation of massive ($M > 8~$M$_{\odot}$) binaries in close orbits. For example, the detection of a very low velocity dispersion among the massive stars in the young region M17 and the measurement of a positive trend of velocity dispersion with age in Galactic clusters. The velocity dispersion observed in M17 could be explained either by the lack of binaries among the stars in this region or by larger binary separations than typically observed, but with a binary fraction similar to other young Galactic clusters. The latter implies that over time, the binary components migrate toward each other. We aim to determine the origin of the strikingly low velocity dispersion by determining the observed and intrinsic binary fraction of massive stars in M17 through multi-epoch spectroscopy. We performed a multi-epoch spectroscopic survey consisting of three epochs separated by days and months. We determine the radial velocity of each star at each epoch by fitting the stellar absorption profiles. We determine an observed binary fraction of 27% and an intrinsic binary fraction of 87%, consistent with that of other Galactic clusters. We conclude that the low velocity dispersion is due to a large separation among the young massive binaries in M17. Our result is in agreement with a migration scenario in which massive stars are born in binaries or higher order systems at large separation and harden within the first million years of evolution. Such an inward migration may either be driven by interaction with a remnant accretion disk, with other young stellar objects present in the system or by dynamical interactions within the cluster. Our results imply that possibly both dynamical interactions and binary evolution are key processes in the formation of gravitational wave sources., Comment: Accepted for publication in A&A. 18 pages, 8 figures. Appendices A and C will be available in electronic form at the CDS

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

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

In order to predict the black hole mass distributions at high redshift, we need to understand whether very massive single stars ($M>40$ M$_\odot$) at low metallicity $Z$ lose their hydrogen-rich envelopes, like their metal-rich counterparts, or whether a binary companion is required to achieve this. To test this, we undertake a deep spectroscopic search for binary companions of the seven apparently single Wolf-Rayet (WR) stars in the Small Magellanic Cloud (SMC; $Z \simeq 1/5 Z_\odot$). For each of them, we acquired six high-quality VLT-UVES spectra spread over 1.5 years. By using the narrow N V lines in these spectra, we monitor radial velocity (RV) variations to search for binary motion. We find low RV variations between 6 and 23 km/s for the seven WR stars, with a median standard deviation of $5$ km/s. Our Monte Carlo simulations imply probabilities below ~5% for any of our target WR stars to have a binary companion more massive than ~5 M$_\odot$ at orbital periods of less than a year. We estimate that the probability that all our target WR stars have companions with orbital periods shorter than 10 yr is below ~10$^{-5}$, and argue that the observed modest RV variations may originate from intrinsic atmosphere or wind variability. Our findings imply that metal-poor massive stars born with $M \gtrsim 40$ M$_\odot$ can lose most of their hydrogen-rich envelopes via stellar winds or eruptive mass loss, which strongly constrains their initial mass - black hole mass relation. We also identify two of our seven target stars (SMC AB1 and SMC AB11) as runaway stars with a peculiar radial velocity of ~80 km/s. Moreover, with all five previously detected WR binaries in the SMC exhibiting orbital periods of below 20 d, a puzzling absence of intermediate-to-long-period WR binaries has emerged, with strong implications for the outcome of massive binary interaction at low metallicity., Comment: Accepted for publication in Astronomy & Astrophysics

The spin evolution of stars in close binary systems can be strongly affected by tides. We investigate the rotational synchronisation of the stellar components for 69 SB1 systems and 14 SB2 B-type systems in the 30 Doradus region of the Large Magellanic Cloud using observations from the VFTS and BBC surveys. Their orbital periods range from a few to a few hundred days, while estimated primary masses for these systems are in the range 5-20 Msun with mass ratio ranges of approximately q=0.03-0.5 and q=0.6-1.0 for the SB1 and SB2 systems, respectively. Projected rotational velocities of the stellar components have been compared with their synchronous velocities derived from the orbital periods. We find that effectively all systems with orbital period of more than 10 days must be asynchronous, whilst all the systems with periods of less than 3 days are likely synchronised. In terms of the stellar fractional radius (r), our results imply that all systems with r<0.1 are asynchronous, with those having r>0.2 probably being synchronised. For the apparently synchronised systems our results are more consistent with synchronisation at the mean orbital angular velocity rather than with that at periastron., Comment: 13 pages, 5 figures. submitted to A&A

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

Given the high incidence of binaries among mature field massive stars, it is clear that multiplicity is an inevitable outcome of high-mass star formation. Understanding how massive multiples form requires the study of the birth environments of massive stars, covering the innermost to outermost regions. We aim to detect and characterise low-mass companions around massive young stellar objects (MYSOs) during and shortly after their formation phase. To investigate large spatial scales, we carried out an $L'$-band high-contrast direct imaging survey seeking low-mass companions (down to $L_{\text{bol}}\approx 10 L_{\odot}$, or late A-type) around thirteen previously identified MYSOs using the VLT/NACO instrument. From those images, we looked for the presence of companions on a wide orbit, covering scales from 300 to 56,000 au. Detection limits were determined for all targets and we tested the gravitational binding to the central object based on chance projection probabilities. We have discovered a total of thirty-nine potential companions around eight MYSOs, the large majority of which have never been reported to date. We derived a multiplicity frequency (MF) of $62\pm13$% and a companion fraction (CF) of $3.0\pm0.5$. The derived MF and CF are compared to other studies for similar separation ranges. The comparisons are effective for a fixed evolutionary stage spanning a wide range of masses and vice versa. We find an increased MF and CF compared to the previous studies targeting MYSOs, showing that the statement in which multiplicity scales with primary mass also extends to younger evolutionary stages. The separations at which the companions are found and their location with relation to the primary star allow us to discuss the implications for the massive star formation theories., Comment: 24 pages, 2 Tables, 5 Figures and 4 appendices (containing figures)

Over the past decades, libraries of stellar spectra have been used in a large variety of science cases, including as sources of reference spectra for a given object or a given spectral type. Despite the existence of large libraries and the increasing number of projects of large-scale spectral surveys, there is to date only one very high-resolution spectral library offering spectra from a few hundred objects from the southern hemisphere (UVES-POP) . We aim to extend the sample, offering a finer coverage of effective temperatures and surface gravity with a uniform collection of spectra obtained in the northern hemisphere. Between 2010 and 2020, we acquired several thousand echelle spectra of bright stars with the Mercator-HERMES spectrograph located in the Roque de Los Muchachos Observatory in La Palma, whose pipeline offers high-quality data reduction products. We have also developed methods to correct for the instrumental response in order to approach the true shape of the spectral continuum. Additionally, we have devised a normalisation process to provide a homogeneous normalisation of the full spectral range for most of the objects. We present a new spectral library consisting of 3256 spectra covering 2043 stars. It combines high signal-to-noise and high spectral resolution over the entire range of effective temperatures and luminosity classes. The spectra are presented in four versions: raw, corrected from the instrumental response, with and without correction from the atmospheric molecular absorption, and normalised (including the telluric correction)., Comment: 17 pages, 18 figures Preview and access to the library: https://www.royer.se/melchiors.html

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

Context. MWC 656 was reported as the first known Be star with a black-hole (BH) companion in a 60 d period. The mass of the proposed BH companion is estimated to be between 4 - 7 MSun. This estimate is based on radial velocity (RV) measurements derived from the Fe ii 4583 emission line of the Be star disc and from the He ii 4686 emission line, assumed to be formed in a disc around the putative BH. Aims. Using new high-resolution spectroscopic data, we investigate whether MWC 656 truly contains a BH. Methods. We used the cross-correlation method to calculate the RVs of both the Be star and the He ii 4686 emission line and we derive a new orbital solution. We also performed disentangling to look for the spectral signature of a companion. Results. We derive an orbital period of 59.028 +- 0.011 d and a mass ratio q = M_Heii/M_Be = 0.12 +- 0.03, much lower than the previously reported q = 0.41 +- 0.07. Adopting a mass of the Be star of M_Be = 7.8 +- 2.0MSun, the companion has a mass of 0.94 +- 0.34MSun. For the upper limit of M_Be = 16MSun and q = 0.15, the companion has a mass 2.4MSun. Performing disentangling on mock spectra shows that the spectral signature of a non-degenerate stellar companion with such a low mass cannot be retrieved using our data. Conclusions. Our measurements do not support the presence of a BH companion in MWC 656. The derived upper limit on the mass of the companion rather indicates that it is a neutron star, a white dwarf, or a hot helium star. Far-UV data will help to reject or confirm a hot helium-star companion., Comment: 6 pages, 2 + 6 figures, 3 + 1 tables, accepted for publication in A&A

Magnetars are highly magnetized neutron stars; their formation mechanism is unknown. Hot helium-rich stars with spectra dominated by emission lines are known as Wolf-Rayet stars. We observe the binary system HD 45166 using spectropolarimetry, finding that it contains a Wolf-Rayet star with a mass of 2 solar masses and a magnetic field of 43 kilogauss. Stellar evolution calculations indicate that this component will explode as a type Ib or IIb supernova, and the strong magnetic field favors a magnetar remnant. We propose that the magnatized Wolf-Rayet star formed by the merger of two lower mass helium stars., Comment: Published in Science on the 18 August 2023. Radial velocities, spectra, and software available in: https://zenodo.org/record/8042656 ESO press release: www.eso.org/public/news/eso2313

Recent studies of massive binaries with putative black hole companions have uncovered a phase of binary evolution that has not been observed before, featuring a bloated stripped star that very recently ceased transferring mass to a main-sequence companion. In this study, we focus on the candidate system VFTS 291, a binary with an orbital period of 108 d and a high semi-amplitude velocity ($K_1=93.7\pm0.2$ km s$^{-1}$). Through our analysis of the disentangled spectra of the two components, together with dynamical and evolutionary arguments, we identify a narrow-lined star of ~1.5-2.5 $M_\odot$ dominating the spectrum, and an early B-type main-sequence companion of $13.2\pm1.5$ $M_\odot$. The low mass of the narrow-lined star, and the high mass ratio, suggest that VFTS 291 is a post-mass-transfer system, with the narrow-lined star being bloated and stripped of its hydrogen-rich envelope, sharing many similarities with other recently discovered stripped stars. Our finding is supported by our detailed binary evolution models, which indicate that the system can be well explained by an initial configuration consisting of an 8.1 $M_\odot$ primary with an 8 $M_\odot$ companion in a 7 d orbital period. While some open questions remain, particularly concerning the surface helium enrichment of the stripped star and the rotational velocity of the companion, we expect that high-resolution spectroscopy may help reconcile our estimates with theory. Our study highlights the importance of multi-epoch spectroscopic surveys to identify and characterize binary interaction products, and provides important insights into the evolution of massive binary stars., Comment: 25 pages, 17 figures. Published by MNRAS, Advance Access

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 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

The 'Great Dimming' of the prototypical red supergiant Betelgeuse, which occurred between December 2019 and April 2020, gives us unprecedented insight into the processes occurring on the stellar surface and in the inner wind of this type of star. In particular it may bring further understanding of their dust nucleation and mass loss processes. Here, we present and analyse VLTI/MATISSE observations in the N-band (8 - 13 $\mu$m) taken near the brightness minimum in order to assess the status of the dusty circumstellar environment. We explore the compatibility of a dust clump obscuring the star with our mid-infrared interferometric observations using continuum 3D radiative transfer modelling, and probe the effect of adding multiple clumps close to the star on the observables. We also test the viability of a large cool spot on the stellar surface without dust present in the ambient medium. Using the visibility data, we derive a uniform disk diameter of 59.02 $\pm$ 0.64 mas in the spectral range 8 to 8.75 $\mu$m. We find that both the dust clump and the cool spot models are compatible with the data. Further to this, we note that the extinction and emission of our localised dust clump in the line of sight of the star, directly compensate each other making the clump undetectable in the spectral energy distribution and visibilities. The lack of infrared brightening during the 'Great Dimming' therefore does not exclude extinction due to a dust clump as one of the possible mechanisms. The visibilities can be reproduced by a spherical wind with dust condensing at 13 stellar radii and a dust mass-loss rate of (2.1 - 4.9) $\times$ 10$^{-10}$ $\mathit{M}_{\odot} {\rm yr}^{-1}$, however, in order to reproduce the complexity of the observed closure phases, additional surface features or dust clumps would be needed., Comment: 13 Pages, accepted for publication 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

Massive stars are powerful cosmic engines. In the phases immediately preceding core-collapse, massive stars in the Galaxy with $M_i \gtrsim 20$ $M_{\odot}$ may appear as classical Wolf-Rayet (WR) stars. As the final contribution of a homogeneous RV survey, this work constrains the multiplicity properties of northern Galactic late-type nitrogen-rich Wolf-Rayet (WNL) stars. We compare their intrinsic binary fraction and orbital period distribution to the carbon-rich (WC) and early-type nitrogen-rich (WNE) populations from previous works. We obtained high-resolution spectra of the complete magnitude-limited sample of 11 Galactic WNL stars with the Mercator telescope on the island of La Palma. We used cross-correlation to measure relative RVs and flagged binary candidates based on the peak-to-peak RV dispersion. By using Monte Carlo sampling and a Bayesian framework, we computed the three-dimensional likelihood and one-dimensional posteriors for the upper period cut-off, power-law index, and intrinsic binary fraction. Adopting a threshold $C$ of 50 km s$^{-1}$, our Bayesian analysis produced an intrinsic fraction of $0.42\substack{+0.15 \\ -0.17}$ for the parent WNL population alongside distributions for the power-law index and the orbital periods. The observed period distribution of Galactic WN and WC binaries from the literature is in agreement with what is found. The period distribution of Galactic WN binaries peaks at $P{\sim}1$-$10$d and that of the WC population at $P{\sim}5000\,$d. This shift cannot be reconciled by orbital evolution due to mass loss or mass transfer. At long periods, the evolutionary sequence O($\xrightarrow{}$LBV)$\xrightarrow{}$WN$\xrightarrow{}$WC seems feasible. The high frequency of short-period WN binaries compared to WC binaries suggests that they either tend to merge or that the WN components in these binaries rarely evolve into WC stars in the Galaxy., Comment: Accepted for publication in Astronomy and Astrophysics

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

Hundreds of black holes with massive main-sequence companions (OB+BHs) might be identified from Gaia astrometry with the Astrometric Mass-Ratio Function (AMRF). We investigate the impact of blue supergiant companions (BSG) instead of dwarfs and the presence of additional companions in the system that are unresolved by Gaia on the astrometric identification of OB+BHs. We also explore how accurate the primary mass needs to be constrained. Moreover, we assess how the high-precision publishing constraints of astrometric binary orbits in the latest Gaia data release DR3 impact the detection of OB+BHs. We establish a BSG mass-magnitude relation and compute BSG AMRF curves. From a mock population of non- and single-degenerate massive binaries and triples with OB or BSG primaries, we asses the false-positive identification fraction and the effect of the BSG AMRF curves. We compare the number of stars with astrometric DR3 orbits in the second Alma Luminous Star catalogue (ALSII) with new predictions on the OB+BH detection using the conservative DR3 publishing criteria. BSG primaries and triples are not expected to impact the false-positive identification fractions significantly. However, if the evolutionary stage of the primary star is unknown, the usage of the BSG curves is still recommended. This also reduces the OB+BH identification fraction significantly. The primary mass does not need to be known to benefit from the high identification fraction of OB+BHs, while keeping the fraction of false-positives low. We find no OB+BH candidates in the ALS II among the DR3 orbits. This null-detection cannot be attributed to the underlying BH-formation scenario, but rather to the stringent DR3 selection criteria. To infer the BH-formation scenario with Gaia, we propose that the constraint on the relative parallax precision in DR4 should be 95% less conservative than the DR3 criterion.(Abridged), Comment: 15 pages (+ 5 appendix), 10 figures, 3 tables, submitted to A&A

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 ESA/KU Leuven CubeSpec mission is specifically designed to provide low-cost space-based high-resolution optical spectroscopy. Here we highlight the science requirements and capabilities of CubeSpec. The primary science goal is to perform pulsation mode identification from spectroscopic line profile variability and empower asteroseismology of massive stars., Comment: Authors' submitted version of poster proceedings paper for IAU Symposium 361: Massive Stars Near and Far, held in Ballyconnell, Ireland, 9-13 May 2022

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

Most massive stars reside in multiple systems that will interact over the course of their lifetime. Classical Wolf-Rayet (WR) stars represent the final end stages of stellar evolution at the upper-mass end. As part of a homogeneous, magnitude-limited ($V\leq12$) spectroscopic survey of northern Galactic WR stars, this paper aims to establish the observed and intrinsic multiplicity properties of the early-type nitrogen-rich WR population (WNE). We obtained high-resolution spectroscopic time series of the complete magnitude-limited sample of 16 WNE stars observable with the 1.2 m Mercator telescope at La Palma, typically providing a time base of about two to eight years. We measured relative radial velocities (RVs) using cross-correlation and used RV variations to flag binary candidates. Adopting a peak-to-peak RV variability threshold of 50 km/s as a criterion found an observed multiplicity fraction of 0.44$\pm$0.12. Using an updated Monte Carlo method with a Bayesian framework, we calculated the three-dimensional likelihood for the intrinsic binary fraction, the maximum period, and the power-law index for the period distribution for the WNE population. We also re-derived multiplicity parameters for the Galactic WC population. We found an intrinsic multiplicity fraction of $0.56\substack{+0.20 \\ -0.15}$ for the parent WNE population. For the Galactic WC population, we re-derive an intrinsic multiplicity fraction of $0.96\substack{+0.04 \\ -0.22}$. The derived multiplicity parameters for the WNE population are quite similar to those derived for main-sequence O binaries but differ from those of the WC population. The significant shift in the WC period distribution towards longer periods is too large to be explained via expansion of the orbit due to stellar winds, and we discuss possible implications of our results., Comment: 17 pages, 12 figures

In accelerating and supersonic media, the interaction of photons with spectral lines can be of ultimate importance. However, fully accounting for such line forces currently can only be done by specialised codes in 1-D steady-state flows. More general cases and higher dimensions require alternative approaches. We presented a comprehensive and fast method for computing the radiation line-force using tables of spectral line-strength distribution parameters, which can be applied in arbitrary (multi-D, time-dependent) simulations, including those accounting for the line-deshadowing instability, to compute the appropriate opacities. We assumed local thermodynamic equilibrium to compute a flux-weighted line opacity from $>4$ million spectral lines. We derived the spectral line strength and tabulated the corresponding line-distribution parameters for a range of input densities $\rho\in[10^{-20},10^{-10}]gcm^{-3}$ and temperatures $T\in[10^4,10^{4.7}]K$. We found that the variation of the line distribution parameters plays an essential role in setting the wind dynamics in our models. In our benchmark study, we also found a good overall agreement between the O-star mass-loss rates of our models and those derived from steady-state studies using more detailed radiative transfer. Our models reinforce that self-consistent variation of the line-distribution parameters is important for the dynamics of line-driven flows. Within a well-calibrated O-star regime, our results support the proposed methodology. In practice, utilising the provided tables, yielded a factor $>100$ speed-up in computational time compared to specialised 1-D model-atmosphere codes of line-driven winds, which constitutes an important step towards efficient multi-D simulations. We conclude that our method and tables are ready to be exploited in various radiation-hydrodynamic simulations where the line force is important.

Spectroscopic multiplicity surveys of O stars in young clusters and OB associations have revealed that a large portion ($\sim$ 70%) of these massive stars (M$_{i}$ $\gt$ 15 $M_{\odot}$) belong to close and short-period binaries (physical separation d $\lt$few au). Despite the recent and significant progress, the formation mechanisms leading to such close massive multiple systems remain to be elucidated. As a result, young massive close binaries (or higher-order multiple systems) are unique laboratories to figure out the pairing mechanism of high-mass stars. We present the first VLTI/GRAVITY observations of six young O stars in the M17 star-forming region ($\lesssim$ 1 Myr) and two additional foreground stars. From the interferometric model fitting of visibility amplitudes and closure phases, we search for companions and measure their positions and flux ratios. Combining the resulting magnitude difference with atmosphere models and evolutionary tracks, we further constrain the masses of the individual components. All of the six high-mass stars are in multiple systems, leading to a multiplicity fraction (MF) of 100%, yielding a 68% confidence interval of 94-100%. We detect a total number of 9 companions with separations up to 120 au. Including previously identified spectroscopic companions, the companion fraction of the young O stars in our sample reaches 2.3$\pm$0.6. The derived masses span a wide range from 2.5 to 50 $M_{\odot}$, with a great tendency towards high-mass companions. While based on a modest sample, our results clearly indicate that the origin of the high degree of multiplicity is rooted in their star formation mechanism. No clear evidence for one of the competing concepts of massive star formation (core accretion or competitive accretion) could be found. However, our results are compatible with migration as a scenario for the formation of close massive binaries., Comment: 20 pages, 14 figures, 4 tables, 3 appendices

The evolutionary status of Be-type stars remains unclear, with both single-star and binary pathways having been proposed. Here, VFTS spectroscopy of 73 Be-type stars, in the spectral-type range, B0--B3, is analysed to estimate projected rotational velocities, radial velocities and stellar parameters. They are found to be rotating faster than the corresponding VFTS B-type sample but simulations imply that their projected rotational velocities are inconsistent with them all rotating at near critical velocities. The de-convolution of the projected rotational velocities estimates leads to a mean rotational velocity estimate of 320-350 km/s, approximately 100 km/s larger than that for the corresponding B-type sample. There is a dearth of targets with rotational velocities less than 0.4 of the critical velocity, with a broad distribution reaching up to critical rotation. Our best estimate for the mean or median of the rotational velocitiy is 0.68 of the critical velocity. Rapidly-rotating B-type stars are more numerous than their Be-type counterparts, whilst the observed frequency of Be-type stars identified as binary systems is significantly lower than that for normal B-type stars, consistent with their respective radial-velocity dispersions. The semi-amplitudes for the Be-type binaries are also smaller. Similar results are found for a SMC Be-type sample centred on NGC346 with no significant differences being found between the two samples. These results are compared with the predictions of single and binary stellar evolutionary models for Be-type stars. Assuming that a single mechanism dominated the production of classical Be-type stars, our comparison would favour a binary evolutionary history., Comment: Accepted for publication in MNRAS, Appendicesand Tables 2 and 4 published as supplementary material

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

In the realm of massive stars, strong binary interaction is commonplace. One extreme case are overcontact systems, which is expected to be part of the evolution of all stars evolving towards a merger, and is hypothesized to play a role in the formation of binary black holes. However, important simplifications are made to model the evolution of overcontact binaries, namely that tidal or rotational deformation is frequently ignored. Yet, both observation and theory show that overcontact stars are heavily tidally deformed, leaving a potentially important effect on the outer layers unaccounted for in models. In this work we develop the methodology to represented tidally deformed stars. Using numerical methods, we compute the structure correction factors to the 1D spherical stellar structure equations due to the binary Roche potential, and compare them to existing results and the structure corrections of single rotating stars. We implement the new structure correction factors into the stellar evolution code MESA and explore several case studies. We compare the differences between our simulations when no rotation is included, when we treat rotation using single star corrections (i.e. only accounting for centrifugal deformation) or when we use tidal deformation. We find that ignoring rotation in deformed detached eclipsing binaries can produce a radius discrepancy of up to 5%. The difference between tidal and single star centrifugal distortion models is more benign at 1%, showing that single rotating star models are a suitable approximation of tidally deformed stars in a binary system. In overcontact configurations, we find a similar 5% variation in surface properties as a result of tidal distortion with respect to non-rotating models, showing that it is inappropriate to model binary stars that fill their Roche lobe significantly, as non-rotating., Comment: 12 pages, 9 figures, accepted to appear in Astronomy & Astrophysics

Almost all massive stars are part of a binary system. Given the wide range of separations at which these companions are found, several observational techniques have been adopted to characterize them, but contrasts greater than 4 in the H-band have never been reached between 0".1 and 1". We used VLT/SPHERE to observe simultaneously with the IRDIS and IFS sub-systems 18 O-type stars within 6 kpc and ages between 1-5 Myrs to probe the existence of stellar companions in the angular separation range from 0".15 to 6" down to very low mass ratios. The IFS YJH- band observations have allowed us to probe the presence of sub-solar companions in a 1".7x1".7 field-of-view down to magnitude limits of deltaH=10 at 0".4. In the wider 12"x12" IRDIS field-of-view, we reached contrasts of deltaK=12 at 1", enabling us to look for even fainter companions. This paper presents five newly discovered intermediate (<1") separation companions, three of which are smaller than 0.2M_sun. If confirmed by future analyses of proper motions, these new detections represent the lowest-mass companions ever found around O-type stars. Assuming that all sources detected within 1" are physically bound, the observed fraction of companions for O-type stars between 0".15 and 0".9 is 0.39+/-0.15, whereas it increases to 1.6+/-0.3 in the separation range from 0".9 to 6". These findings clearly support the notion that massive stars form almost exclusively in multiple systems, and that larger AO-assisted coronagraphic surveys are crucial in placing constraints on the multiplicity properties of massive star companions in regions of the parameter space that have previously gone unexplored, and demonstrate that the companion mass function is populated down to the lowest stellar masses., Comment: Accepted for publication in Astronomy & Astrophysics. 14 pages, 12 figures

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

Context. There is currently a niche for providing high-cadence, high resolution, time-series optical spectroscopy from space, which can be filled by using a low-cost cubesat mission. The Belgian-led ESA CubeSpec mission is specifically designed to provide space-based, low-cost spectroscopy with specific capabilities that can be optimised for a particular science need. Approved as an ESA in-orbit demonstrator, the CubeSpec satellite's primary science objective will focus on obtaining high-cadence, high resolution optical spectroscopic data to facilitate asteroseismology of pulsating massive stars. Aims. In this first paper, we aim to search for pulsating massive stars suitable for the CubeSpec mission, specifically $\beta$ Cep stars, which typically require time series spectroscopy to identify the geometry of their pulsation modes. Methods. Based on the science requirements needed to enable asteroseismology of massive stars with the capabilities of CubeSpec's spectrograph, we combine a literature study for pulsation with the analysis of recent high-cadence time series TESS photometry to classify the variability for stars brighter than V < 4 mag and between O9 and B3 in spectral type. Results. Among the 90 stars that meet our magnitude and spectral type requirements, we identify 23 promising $\beta$ Cep stars with high-amplitude (non-)radial pulsation modes with frequencies below 7 d$^{-1}$. Using further constraints on projected rotational velocities, pulsation amplitudes and number of pulsation modes, we devise a prioritised target list for the CubeSpec mission according to its science requirements and the potential of the targets for asteroseismology. The full target catalogue further provide a modern TESS-based review of line profile and photometric variability properties among bright O9-B3 stars., Comment: 14 pages, accepted for publication in A&A

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

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)

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.

We present results from the B-type Binaries Characterisation (BBC) programme, a multi-epoch spectroscopic study of 88 early B-type binary candidates in the 30 Doradus region of the Large Magellanic Cloud (LMC). From radial-velocity analysis of 29 observational epochs we confirm the binary status of 64 of our targets, comprising 50 SB1 and 14 SB2 B-type binaries. A further 20 systems (classified as SB1*) show clear signs of periodicity but with more tentative periods. Orbital solutions are presented for these 84 systems, providing the largest homogeneous sample to date of the binary properties of early B-type stars. Our derived orbital-period distribution is generally similar to those for samples of more massive (O-type) binaries in both the LMC and the Galaxy. This similarity with the properties of the more massive O-type binaries is important as early B-type stars are expected to account for the majority of core-collapse supernovae. Differences in the period distributions of the different samples start to increase above 4 d, and are also present between the earliest (B0-0.7) and later-type (B1-2.5) systems within the BBC sample, although further study is required to understand if this is an observational bias or a real physical effect. We have examined the semi-amplitude velocities and orbital periods of our sample to identify potential candidates that could hide compact companions. Comparing with probability distributions of finding black hole companions to OB-type stars from a recent theoretical study, we have found 16 binaries in the higher probability region that warrant further study., Comment: 42 pages (main body: 28 pages, supplementary material: 14 pages), 20 figures. Accepted for publication in MNRAS

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)

Context. Red supergiants are observed to undergo vigorous mass-loss. However, to date, no theoretical model has succeeded in explaining the origins of these objects' winds. This strongly limits our understanding of red supergiant evolution and Type II-P and II-L supernova progenitor properties. Aims. We examine the role that vigorous atmospheric turbulence may play in initiating and determining the mass-loss rates of red supergiant stars. Methods. We analytically and numerically solve the equations of conservation of mass and momentum, which we later couple to an atmospheric temperature structure, to obtain theoretically motivated mass-loss rates. We then compare these to state-of-the-art empirical mass-loss rate scaling formulae as well as observationally inferred mass-loss rates of red supergiants. Results. We find that the pressure due to the characteristic turbulent velocities inferred for red supergiants is sufficient to explain the mass-loss rates of these objects in the absence of the normally employed opacity from circumstellar dust. Motivated by this initial success, we provide a first theoretical and fully analytic mass-loss rate prescription for red supergiants. We conclude by highlighting some intriguing possible implications of these rates for future studies of stellar evolution, especially in light of the lack of a direct dependence on metallicity., Comment: 14 pages, 9 figures, 2 tables

The processes by which red supergiants lose mass are not fully understood thus-far and their mass-loss rates lack theoretical constraints. The ambient surroundings of the nearby M0.5 Iab star Antares offers an ideal environment to obtain detailed empirical information on the outflow properties at its onset, and hence indirectly, on the mode(s) of mass loss. We present and analyse optical VLT/SPHERE/ZIMPOL polarimetric imaging with angular resolution down to 23 milli-arcsec, sufficient to spatially resolve both the stellar disk and its direct surroundings. We detect a conspicuous feature in polarised intensity that we identify as a clump containing dust, which we characterise through 3D radiative transfer modelling. The clump is positioned behind the plane of the sky, therefore has been released from the backside of the star, and its inner edge is only 0.3 stellar radii above the surface. The current dust mass in the clump is $1.3^{+0.2}_{-1.0} \times 10^{-8}$ M$_{\odot}$, though its proximity to the star implies that dust nucleation is probably still ongoing. The ejection of clumps of gas and dust makes a non-negligible contribution to the total mass lost from the star which could possibly be linked to localised surface activity such as convective motions or non-radial pulsations., Comment: 15 pages, accepted for publication in MNRAS

The majority of massive stars ($>8$ $\rm{M_{\odot}}$) in OB associations are found in close binary systems. Nonetheless, the formation mechanism of these close massive binaries is not understood yet. Using literature data, we measured the radial-velocity dispersion ($\sigma_\mathrm{RV}$) as a proxy for the close binary fraction in ten OB associations in the Galaxy and the Large Magellanic Cloud, spanning an age range from 1 to 6 Myrs. We find a positive trend of this dispersion with the cluster's age, which is consistent with binary hardening. Assuming a universal binary fraction of $f_\mathrm{bin}$ = 0.7, we converted the $\sigma_\mathrm{RV}$ behavior to an evolution of the minimum orbital period $P_\mathrm{cutoff}$ from $\sim$9.5 years at 1 Myr to $\sim$1.4 days for the oldest clusters in our sample at $\sim$6 Myr. Our results suggest that binaries are formed at larger separations, and they harden in around 1 to 2 Myrs to produce the period distribution observed in few million year-old OB binaries. Such an inward migration may either be driven by an interaction with a remnant accretion disk or with other young stellar objects present in the system. Our findings constitute the first empirical evidence in favor of migration as a scenario for the formation of massive close binaries., Comment: 7 pages. 5 figures. Accepted for publication in Astronomy & Astrophysics

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

Binary interactions dominate the evolution of massive stars, but their role is less clear for low- and intermediate-mass stars. The evolution of a spherical wind from an asymptotic giant branch (AGB) star into a nonspherical planetary nebula (PN) could be due to binary interactions. We observed a sample of AGB stars with the Atacama Large Millimeter/submillimeter Array (ALMA) and found that their winds exhibit distinct nonspherical geometries with morphological similarities to planetary nebulae (PNe). We infer that the same physics shapes both AGB winds and PNe; additionally, the morphology and AGB mass-loss rate are correlated. These characteristics can be explained by binary interaction. We propose an evolutionary scenario for AGB morphologies that is consistent with observed phenomena in AGB stars and PNe., Comment: 19 pages main journal, 97 pages Supplementary Information

About 20% of all B-type stars are classical Be stars. The Be phenomenon is strongly correlated with rapid rotation, the origin of which remains unclear. It may be rooted in single- or binary-star evolution. In the framework of the binary channel, the initially more massive star transfers mass and angular momentum to the original secondary, which becomes a Be star. The system then evolves into a Be binary with a post-main-sequence companion, which may later be disrupted in a supernova event. Hence, if the binary channel dominates the formation of Be stars, one may expect a strong lack of close Be binaries with main sequence (MS) companions. Through an extensive, star-by-star review of the literature of a magnitude-limited sample of Galactic early-type Be stars, we investigate whether Be binaries with MS companions are known to exist. Our sample is constructed from the BeSS database and cross-matched with all available literature on the individual stars. Out of an initial list of 505 Be stars, we compile a final sample of 287 Galactic Be stars earlier than B1.5 with V<=12 mag. Out of those, 13 objects were reported as Be binaries with known post-MS companions and 11 as binaries with unknown, uncertain or debated companions. We find no confirmed reports of Be binaries with MS companions. For the remaining 263 targets, no significant reports of multiplicity exist in the literature, implying that they are either Be binaries with faint companions, or truly single. The clear lack of reported MS companions to Be stars, which stands in contrast to the high number of detected B+B MS binaries, strongly supports the hypothesis that early-type Be stars are binary interaction products that spun up after mass and angular momentum transfer from a companion star. Taken at face value, our results may suggest that a large majority of the early-type Be stars have formed through binary mass-transfer., Comment: 15 pages (incl. appendix), 6 figures, 3 tables, accepted for publication in A&A