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.

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.

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.

Massive stars like company. However, low-mass companions have remained extremely difficult to detect at angular separations ($\rho$) smaller than 1" (~1000-3000 au considering typical distance to nearby massive stars) given the large brightness contrast between the companion and the central star. Constraints on the low-mass end of the companions mass-function for massive stars are however needed, for example to help distinguishing between various scenarios for the formation of massive stars. To obtain statistically significant constraint on the presence of low-mass companions beyond the typical detection limit of current surveys ($\Delta \mathrm{mag} \lesssim 5$ at $\rho \lesssim 1$"), we initiated a survey of O and Wolf-Rayet stars in the Carina region using the SPHERE coronagraphic instrument on the VLT. In this first paper, we aim to introduce the survey, to present the methodology and to demonstrate the capability of SPHERE for massive stars using the multiple system QZ~Car. High-contrast imaging techniques, such as angular- and spectral-differential imaging techniques as well as PSF-fitting, were applied to detect and measure the relative flux of companions in each spectral channel of the instrument. We detected 19 sources around the QZ~Car system with detection limits of 9~mag at $\rho > 200$~mas for IFS and as faint as 13~mag at $\rho$ > 1".8 for IRDIS (corresponding to sub-solar masses for potential companions). All but two are reported here for the first time. Based on this proof of concept, we showed that VLT/SPHERE allows us to reach the sub-solar mass regime of the companion mass function. This paves the way for this type of observation with a large sample of massive stars to provide novel constraints on the multiplicity of massive stars in a region of the parameter space that has remained inaccessible so far.

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.

The multiplicity properties of massive stars are one of the important outstanding issues in stellar evolution. Quantifying the binary statistics of all evolutionary phases is essential to paint a complete picture of how and when massive stars interact with their companions, and to determine the consequences of these interactions. We investigate the multiplicity of an almost complete census of red supergiant stars (RSGs) in NGC 330, a young massive cluster in the SMC. Using a combination of multi-epoch HARPS and MUSE spectroscopy, we estimate radial velocities and assess the kinematic and multiplicity properties of 15 RSGs in NGC 330. Radial velocities are estimated to better than +/-100 m/s for the HARPS data. The line-of-sight velocity dispersion for the cluster is estimated as 3.20 +0.69-0.52 km/s. When virial equilibrium is assumed, the dynamical mass of the cluster is log (M{dyn} /M{sun}) = 5.20+/-0.17, in good agreement with previous upper limits. We detect significant radial velocity variability in our multi-epoch observations and distinguish between variations caused by atmospheric activity and those caused by binarity. The binary fraction of NGC 330 RSGs is estimated by comparisons with simulated observations of systems with a range of input binary fractions. In this way, we account for observational biases and estimate the intrinsic binary fraction for RSGs in NGC 330 as f{RSG} = 0.3+/-0.1 for orbital periods in the range 2.3< log P [days] <4.3, with q>0.1. Using the distribution of the luminosities of the RSG population, we estimate the age of NGC 330 to be 45+/-5 Myr and estimate a red straggler fraction of 50%. We estimate the binary fraction of RSGs in NGC 330 and conclude that it appears to be lower than that of main-sequence massive stars, which is expected because interactions between an RSG and a companion are assumed to effectively strip the RSG envelope.

Liu et al. (2019) recently reported the detection of a 68 [+11/-13] solar mass (Msun) black hole (BH) paired with an 8.2 [+0.9/-1.2] Msun B-type sub-giant star in the 78.9-day spectroscopic binary system LB-1. Such a black hole is over twice as massive as any other known stellar-mass black hole with non-compact companions2 and its mass approaches those that result from BH-BH coalescences that are detected by gravitational wave interferometers. Its presence in a solar-like metallicity environment challenges conventional theories of massive binary evolution, stellar winds and core-collapse supernovae, so that more exotic scenarios seem to be needed to explain the existence and properties of LB-1. Here, we show that the observational diagnostics used to derive the BH mass results from the orbital motion of the B-type star, not that of the BH. As a consequence, no evidence for a massive BH remains in the data, therefore solving the existing tension with formation models of such a massive BH at solar metallicity and with theories of massive star evolution in general.

$\gamma$ Cas stars are a $\sim$1% minority among classical Be stars with hard but only moderately strong continuous thermal X-ray flux and mostly very early-B spectral type. The X-ray flux has been suggested to originate from matter accelerated via magnetic disk-star interaction, by a rapidly rotating neutron star (NS) companion via the propeller effect, or by accretion onto a white dwarf (WD) companion. In view of the growing number of identified $\gamma$ Cas stars and the only imperfect matches between these suggestions and the observations, alternative models should be pursued. Two of the three best-observed $\gamma$ Cas stars, $\gamma$ Cas itself and $\pi$ Aqr, have a low-mass companion with low optical flux; interferometry of BZ Cru is inconclusive. Binary-evolution models are examined for their ability to produce such systems. The OB+He-star stage of post-mass transfer binaries, which is otherwise observationally unaccounted, can potentially reproduce many observed properties of $\gamma$ Cas stars. The interaction of the fast wind of helium stars with the disk and/or with the wind of Be stars may give rise to the production of hard X-rays. While not modelling this process, it is shown that the energy budget is favourable, and that the wind velocities may lead to hard X-rays as observed in $\gamma$ Cas stars. Furthermore, their observed number appears to be consistent with the evolutionary models. Within the Be+He-star binary model, the Be stars in $\gamma$ Cas stars are conventional classical Be stars. They are encompassed by O-star+Wolf-Rayet systems towards higher mass, where no stable Be decretion disks exist, and by Be+sdO systems at lower mass where the sdO winds may be too weak to cause the $\gamma$ Cas phenomenon. In decreasing order of the helium-star mass, the descendants could be Be+black-hole, Be+NS or Be+WD binaries.

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.

To understand the nature of transient obscuring outflows in active galactic nuclei, we observed the Seyfert 1 galaxy NGC 3783 on two occasions in December 2016 triggered by Swift monitoring indicating strong soft X-ray absorption in November. We obtained ultraviolet spectra using COS on HST and optical spectra using FEROS on the MPG/ESO 2.2-m telescope that were simultaneous with X-ray spectra from XMM-Newton and NuSTAR. We find new components of broad, blue-shifted absorption associated with Ly$\alpha$, \ion{N}{v}, \ion{Si}{iv}, and \ion{C}{iv} in our COS spectra. The absorption extends from velocities near zero in the rest-frame of the host galaxy to $-6200$ $\rm km~s^{-1}$. These features appear for the first time in NGC 3783 at the same time as the heavy soft X-ray absorption seen in the XMM-Newton X-ray spectra. The X-ray absorption has a column density of $\sim 10^{23}~\rm cm^{-2}$, and it partially covers the X-ray continuum source. The X-ray absorption becomes more transparent in the second observation, as does the UV absorption. Combining the X-ray column densities with the UV spectral observations yields an ionization parameter for the obscuring gas of log $\xi =1.84^{+0.4}_{-0.2}$ $\rm erg~cm~s^{-1}$. Despite the high intensity of the UV continuum in NGC 3783, F(1470 \AA)=$8 \times 10^{-14}~\rm erg~cm^{-2}~s^{-1}~\AA^{-1}$, the well known narrow UV absorption lines are deeper than in earlier observations in unobscured states, and low ionization species such as \ion{C}{iii} appear, indicating that the narrow-line gas is more distant from the nucleus and is being shadowed by the gas producing the obscuration. Despite the high continuum flux levels in our observations of NGC 3783, moderate velocities in the UV broad line profiles have substantially diminished. We suggest that a collapse of the broad line region has led to the outburst and triggered the obscuring event.

Recent studies show that more than 70% of massive stars do not evolve as effectively single stars, but as members of interacting binary systems. The evolution of these stars is thus strongly altered compared to similar but isolated objects. We investigate the occurrence of parsec-scale mid-infrared nebulae around early-type stars. If they exist over a wide range of stellar properties, one possible overarching explanation is non-conservative mass transfer in binary interactions, or stellar mergers. For ~3850 stars (all OBA stars in the Bright Star Catalogue [BSC], Be stars, BeXRBs, and Be+sdO systems), we visually inspect WISE 22 $\mu$m images. Based on nebular shape and relative position, we distinguish five categories: offset bow shocks structurally aligned with the stellar space velocity, unaligned offset bow shocks, and centered, unresolved, and not classified nebulae. In the BSC, we find that 28%, 13%, and 0.4% of all O, B, and A stars, respectively, possess associated infrared (IR) nebulae. Additionally, 34/234 Be stars, 4/72 BeXRBs, and 3/17 Be+sdO systems are associated with IR nebulae. Aligned or unaligned bow shocks result from high relative velocities between star and interstellar medium (ISM) that are dominated by the star or the ISM, respectively. About 13% of the centered nebulae could be bow shocks seen head- or tail-on. For the rest, the data disfavor explanations as remains of parental disks, supernova remnants of a previous companion, and dust production in stellar winds. The existence of centered nebulae also at high Galactic latitudes strongly limits the global risk of coincidental alignments with condensations in the ISM. Mass loss during binary evolution seems a viable mechanism for the formation of at least some of these nebulae. In total, about 29% of the IR nebulae (2% of all OBA stars in the BSC) may find their explanation in the context of binary evolution.

We report the results of spectrophotometric observations of the massive star MN18 revealed via discovery of a bipolar nebula around it with the Spitzer Space Telescope. Using the optical spectrum obtained with the Southern African Large Telescope, we classify this star as B1 Ia. The evolved status of MN18 is supported by the detection of nitrogen overabundance in the nebula, which implies that it is composed of processed material ejected by the star. We analysed the spectrum of MN18 by using the code CMFGEN, obtaining a stellar effective temperature of \approx 21 kK. The star is highly reddened, E(B-V)\approx 2 mag. Adopting an absolute visual magnitude of M_V=-6.8\pm0.5 (typical of B1 supergiants), MN18 has a luminosity of log L/Lsun \approx 5.42\pm0.30, a mass-loss rate of \approx (2.8-4.5)\times10^{-7} Msun/yr, and resides at a distance of \approx 5.6^{+1.5} _{-1.2} kpc. We discuss the origin of the nebula around MN18 and compare it with similar nebulae produced by other blue supergiants in the Galaxy (Sher 25, HD 168625, [SBW2007] 1) and the Large Magellanic Cloud (Sk-69 202). The nitrogen abundances in these nebulae imply that blue supergiants can produce them from the main sequence stage up to the pre-supernova stage. We also present a K-band spectrum of the candidate luminous blue variable MN56 (encircled by a ring-like nebula) and report the discovery of an OB star at \approx 17 arcsec from MN18. The possible membership of MN18 and the OB star of the star cluster Lynga 3 is discussed.