Search

Westerlund 1 (Wd 1) is one of the most relevant star clusters in the Milky Way to study massive star formation, although it is still poorly known. Here, we used photometric and spectroscopic data to model the eclipsing binary W36, showing that its spectral type is O6.5 III + O9.5 IV, hotter and more luminous than thought before. Its distance $d_{\rm W36}$ $=$ 4.03$\pm$0.25 kpc agrees, within the errors, with three recent Gaia-EDR3-based distances reported in Paper I, Beasor & Davies, and by Negueruela's group. However, they follow different approaches to fix the zero-points for red sources such as those in Wd 1 and to select the best approach, we used an accurate modelling of W36. The weighted mean distance of our parallax (Paper I) and binary distances results in $d_{\rm wd1}$ = 4.05 $\pm$0.20 kpc, with an unprecedented accuracy of 5%. We adopted isochrones based on the Geneva code with supersolar abundances to infer the age of W36B as 6.4 $\pm$ 0.7 Myr. This object seems to be part of the prolific star formation burst represented by OB giants and supergiants that occurred at 7.1 $\pm$ 0.5 Myr ago, which coincides with the recently published PMS isochrone with age 7.2 Myr. Other BA-type luminous evolved stars and Yellow Hypergiants spread in the age range of 8--11 Myr. The four Red Supergiants discussed in paper I represent the oldest population of the cluster with an age of 10.7 $\pm$ 1 Myr. The multiple episodes of star formation in Wd 1 are reminiscent of that reported for the R136/30 Dor LMC cluster., 15 pages, 6 figures, 2 appendix (with one table)

Context: The star cluster R136 inside the LMC hosts a rich population of massive stars, including the most massive stars known. The strong stellar winds of these very luminous stars impact their evolution and the surrounding environment. We currently lack detailed knowledge of the wind structure that is needed to quantify this impact. Aims: To observationally constrain the stellar and wind properties of the massive stars in R136, in particular the parameters related to wind clumping. Methods: We simultaneously analyse optical and UV spectroscopy of 53 O-type and 3 WNh-stars using the FASTWIND model atmosphere code and a genetic algorithm. The models account for optically thick clumps and effects related to porosity and velocity-porosity, as well as a non-void interclump medium. Results: We obtain stellar parameters, surface abundances, mass-loss rates, terminal velocities and clumping characteristics and compare these to theoretical predictions and evolutionary models. The clumping properties include the density of the interclump medium and the velocity-porosity of the wind. For the first time, these characteristics are systematically measured for a wide range of effective temperatures and luminosities. Conclusions: We confirm a cluster age of 1.0-2.5 Myr and derive an initial stellar mass of $\geq 250 {\rm M}_\odot$ for the most massive star in our sample, R136a1. The winds of our sample stars are highly clumped, with an average clumping factor of $f_{\rm cl}=29\pm15$. We find tentative trends in the wind-structure parameters as a function of mass-loss rate, suggesting that the winds of stars with higher mass-loss rates are less clumped. We compare several theoretical predictions to the observed mass-loss rates and terminal velocities and find that none satisfactorily reproduces both quantities. The prescription of Krti\v{c}ka & Kub\'at (2018) matches best the observed mass-loss rates., Comment: Accepted for publication in A&A; Appendix I will not be included in the published version

In this paper we present a catalog of 4584 eclipsing binaries observed during the first two years (26 sectors) of the TESS survey. We discuss selection criteria for eclipsing binary candidates, detection of hither-to unknown eclipsing systems, determination of the ephemerides, the validation and triage process, and the derivation of heuristic estimates for the ephemerides. Instead of keeping to the widely used discrete classes, we propose a binary star morphology classification based on a dimensionality reduction algorithm. Finally, we present statistical properties of the sample, we qualitatively estimate completeness, and discuss the results. The work presented here is organized and performed within the TESS Eclipsing Binary Working Group, an open group of professional and citizen scientists; we conclude by describing ongoing work and future goals for the group. The catalog is available from http://tessEBs.villanova.edu and from MAST., 32 pages, 21 figures, accepted to ApJ Supplement Series; comments welcome

Stellar-mass black holes are the final remnants of stars born with more than 15 solar masses. Billions are expected to reside in the Local Group, yet only a few are known, mostly detected through X-rays emitted as they accrete material from a companion star. Here, we report on VFTS 243: a massive X-ray-faint binary in the Large Magellanic Cloud. With an orbital period of 10.4 d, it comprises an O-type star of 25 solar masses and an unseen companion of at least nine solar masses. Our spectral analysis excludes a non-degenerate companion at a 5σ confidence level. The minimum companion mass implies that it is a black hole. No other X-ray-quiet black hole is unambiguously known outside our Galaxy. The (near-)circular orbit and kinematics of VFTS 243 imply that the collapse of the progenitor into a black hole was associated with little or no ejected material or black-hole kick. Identifying such unique binaries substantially impacts the predicted rates of gravitational-wave detections and properties of core-collapse supernovae across the cosmos., This research has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (T.S. and H.S., grant agreement 772225: MULTIPLES). T.S. acknowledges support from the European Union’s Horizon 2020 programme under Marie Skłodowska-Curie grant agreement 101024605. This work is based on observations collected at the ESO under programme IDs 182.D-0222, 090.D-0323 and 092.D-0136. L.M. thanks the ESA and the Belgian Federal Science Policy Office (BELSPO) for their support in the framework of the PRODEX programme. P. Marchant acknowledges support from FWO junior postdoctoral fellowship 12ZY520N. We thank J. Hillier for making the CMFGEN code available. C.H. acknowledges support from the KU Leuven Research Council (grant C16/17/007: MAESTRO). P.A.C. acknowledges support from UK Science and Technology Facilities Council research grant ST/V000853/1. V.H.-B. acknowledges the support of the Natural Sciences and Engineering Research Council of Canada (NSERC) through grant RGPIN-2020-05990. M. Gieles acknowledges support from the Ministry of Science and Innovation through a Europa Excelencia grant (EUR2020-112157). The PoWR code, as well as the associated post-processing and visualization tool WRplot, has been developed under the guidance of W.-R. Hamann with substantial contributions from L. Koesterke, G. Gräfener, A. Sander, T.S. and other co-workers and students. This work has received funding from the ERC under the European Union’s Horizon 2020 research and innovation programme (grant agreement 945806). It is also supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy EXC 2181/1-390900948 (the Heidelberg STRUCTURES Excellence Cluster). S.T. acknowledges support from the Netherlands Research Council NWO (grants VENI 639.041.645, VIDI 203.061). A.H. and D.J.L. acknowledge support by the Spanish MCI through grant PGC-2018-0913741-B-C22 and the Severo Ochoa Programme through CEX2019-000920-S. J.M.A. acknowledges support from the Spanish Government Ministerio de Ciencia, Innovación y Universidades through grant PGC2018-095-049-B-C22.

Given that mergers are often invoked to explain many exotic phenomena in massive star evolution, understanding the evolutionary phase directly preceding a merger, the overcontact phase, is of crucial importance. Despite its importance, large uncertainties exist in our understanding of the evolution of massive overcontact binaries. We aim to provide robust observational constraints on the future dynamical evolution of massive overcontact systems by measuring the rate at which the periods change for a sample of six such objects. Furthermore, we aim to investigate whether the periods of unequal mass systems show higher rates of change than their equal mass counterparts as theoretical models predict. Using archival photometric data from various ground- and space-based missions covering up to ~40 years, we measure the periods of each system over several smaller time spans. We then fit a linear regression through the measured periods to determine the rate at which the period is changing over the entire data set. We find that all of the stars in our sample have very small period changes and that there does not seem to be a correlation with the mass ratio. This implies that the orbital periods for these systems are stable on the nuclear timescale, and that the unequal mass systems may not equalize as expected. When comparing our results with population synthesis distributions, we find large discrepancies between the expected mass ratios and period stabilities. We find that these discrepancies can be mitigated to a degree by removing systems with shorter initial periods, suggesting that the observed sample of overcontact systems may originate from binary systems with longer initial orbital periods., Comment: 11 pages, 6 (+5) Figures, accepted for publication in A&A

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

Barium (Ba) dwarfs and CH subgiants are the less evolved analogues of Ba and CH giants. They are F-to G-Type main-sequence stars polluted with heavy elements by their binary companions when the companion was on the asymptotic giant branch (AGB). This companion is now a white dwarf that in most cases cannot be directly detected. We present a large systematic study of 60 objects classified as Ba dwarfs or CH subgiants. Combining radial-velocity measurements from HERMES and SALT high-resolution spectra with radial-velocity data from CORAVEL and CORALIE, we determine the orbital parameters of 27 systems. We also derive their masses by comparing their location in the Hertzsprung-Russell diagram with evolutionary models. We confirm that Ba dwarfs and CH subgiants are not at different evolutionary stages, and that they have similar metallicities, despite their different names. Additionally, Ba giants appear significantly more massive than their main-sequence analogues. This is likely due to observational biases against the detection of hotter main-sequence post-mass-Transfer objects. Combining our spectroscopic orbits with the HIPPARCOS astrometric data, we derive the orbital inclination and the mass of the WD companion for four systems. Since this cannot be done for all systems in our sample yet (but should be possible with upcoming Gaia data releases), we also analyse the mass-function distribution of our binaries. We can model this distribution with very narrow mass distributions for the two components and random orbital orientations on the sky. Finally, based on BINSTAR evolutionary models, we suggest that the orbital evolution of low-mass Ba systems can be affected by a second phase of interactions along the red giant branch of the Ba star, which impact the eccentricities and periods of the giants., SCOPUS: ar.j, info:eu-repo/semantics/published

O-stars are known to experience a wide range of variability mechanisms originating at both their surface and their near-core regions. Characterization and understanding of this variability and its potential causes are integral for evolutionary calculations. We use a new extensive high-resolution spectroscopic data set to characterize the variability observed in both the spectroscopic and space-based photometric observations of the O+B eclipsing binary HD~165246. We present an updated atmospheric and binary solution for the primary component, involving a high level of microturbulence ($13_{-1.3}^{+1.0}\,$km\,s$^{-1}$) and a mass of $M_1=23.7_{-1.4}^{+1.1}$~M$_{\odot}$, placing it in a sparsely explored region of the Hertzsprung-Russell diagram. Furthermore, we deduce a rotational frequency of $0.690\pm 0.003\,$d$^{-1}$ from the combined photometric and line-profile variability, implying that the primary rotates at 40\% of its critical Keplerian rotation rate. We discuss the potential explanations for the overall variability observed in this massive binary, and discuss its evolutionary context., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society, 14 pages, 10 figures, 6 tables, 1 appendix (4 additional 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. The accuracy of the model is then typically evaluated by comparing the calculated and observed spectra. 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. Thus far, the LB-1 system has been indirectly investigated by 1D stellar-atmosphere codes and by spectral disentangling techniques, yielding differing conclusions about the nature of the system (e.g., a B-star and black-hole binary with an accretion disc around the black hole or a stripped-star and Be-star binary system have been proposed). Aims. To improve our understanding of the LB-1 system, we directly modelled the phase-dependent Hα line profiles of this system. To this end, we developed and present 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α 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.

As potential progenitors of several exotic phenomena including gravitational wave sources, magnetic stars, and Be stars, close massive binary systems probe a crucial area of the parameter space in massive star evolution. Despite the importance of these systems, large uncertainties regarding the nature and efficiency of the internal mixing mechanisms still exist. In this work, we aim to provide robust observational constraints on the internal mixing processes by spectroscopically analyzing a sample of three massive overcontact binaries at different metallicities. Using optical phase-resolved spectroscopic data, we perform an atmosphere analysis using more traditional 1D techniques and using state-of-the-art 3D techniques. We compare and contrast the assumptions and results of each technique and investigate how the assumptions affect the final derived atmospheric parameters. We find that in all three cases, both components of system are highly overluminous indicating either efficient internal mixing of helium or previous non-conservative mass transfer. However, we do not find strong evidence of helium or CNO surface abundance changes usually associated with mixing. Additionally, we find that in unequal mass systems, the measured effective temperature and luminosity of the less massive component places it very close to the more massive component on the Hertzsprung-Russell diagram. These results were obtained independently using both of the techniques mentioned above, which suggests that these measurements are robust. The observed discrepancies between the temperature and the surface abundance measurements when compared to theoretical expectations indicate that unaccounted for additional physical mechanisms may be at play., Comment: Accepted for publication in A&A; 28 pages, 20 figures

Eclipsing binary systems with components that pulsate in gravity modes ($g$ modes) allow for simultaneous and independent constraints of the chemical mixing profiles of stars. We aim to assemble a sample of $g$-mode pulsators in detached eclipsing binaries with the purpose of finding good candidates for future evolutionary and asteroseismic modelling. In addition, we present a case study of the eclipsing binary KIC9850387, identified as our most promising candidate. We selected all of the detached eclipsing binaries in the Kepler eclipsing binary catalogue and performed a visual inspection to determine the presence and density of $g$ modes, and the presence of $g$-mode period-spacing patterns in their frequency spectra. We then characterised our sample based on their $g$-mode pulsational parameters and binary and atmospheric parameters. A spectroscopic follow-up of KIC9850387 was then performed. We then performed spectral disentangling followed by atmospheric modelling and abundance analysis for the primary star. We utilised an iterative approach to simultaneously optimise the pulsational and eclipse models, and subsequently performed an analysis of the pressure- ($p$-) and $g$-mode pulsational frequencies. We compiled a sample of 93 Kepler eclipsing binary stars with $g$-mode pulsating components and identified clear $g$-mode period-spacing patterns in the frequency spectra of seven of these systems. We found that the $g$-mode pulsational parameters and the binary and atmospheric parameters of our sample are weakly correlated at best. We find that the eclipsing binary KIC9850387 is a double-lined spectroscopic binary in a near-circular orbit with an intermediate-mass primary, and a solar-like secondary. We find $\ell=1$ and $\ell=2$ period-spacing patterns in the frequency spectrum of KIC9850387 spanning more than ten radial orders each., 27 pages, 14 figures, accepted for publication in A&A

PHOEBE 2 is a Python package for modeling the observables of eclipsing star systems, but until now has focused entirely on the forward-model -- that is, generating a synthetic model given fixed values of a large number of parameters describing the system and the observations. The inverse problem, obtaining orbital and stellar parameters given observational data, is more complicated and computationally expensive as it requires generating a large set of forward-models to determine which set of parameters and uncertainties best represent the available observational data. The process of determining the best solution and also of obtaining reliable and robust uncertainties on those parameters often requires the use of multiple algorithms, including both optimizers and samplers. Furthermore, the forward-model of PHOEBE has been designed to be as physically robust as possible, but is computationally expensive compared to other codes. It is useful, therefore, to use whichever code is most efficient given the reasonable assumptions for a specific system, but learning the intricacies of multiple codes presents a barrier to doing this in practice. Here we present the 2.3 release of PHOEBE (publicly available from http://phoebe-project.org) which introduces a general framework for defining and handling distributions on parameters, and utilizing multiple different estimation, optimization, and sampling algorithms. The presented framework supports multiple forward-models, including the robust model built into PHOEBE itself., accepted for publication in ApJS

Massive stars play an important role in the mechanical and chemical evolution of galaxies. Understanding the internal processes of these stars is vital to our understanding of their evolution and eventual end products. Deformations from spherical geometry are common for massive stars; however, the tools that are currently available for the study of these systems are almost exclusively one-dimensional. We present a new spectroscopic analysis tool tailored for massive stars that deviate from spherical symmetry. This code (entitled SPAMMS) is a spectroscopic patch model that takes the three-dimensional surface geometry of the system into account to produce spectral profiles at given phases and orientations. In using the Wilson-Devinney-like code PHOEBE in combination with the non local thermodynamic equilibrium (NLTE) radiative transfer code FASTWIND, we created a three-dimensional mesh that represents the surface geometry of our system and we assigned FASTWIND emergent intensity line profiles to each mesh triangle, which take the local parameters such as temperature, surface gravity, and radius into account. These line profiles were then integrated across the visible surface, where their flux contribution and radial velocity are taken into account, thus returning a final line profile for the visible surface of the system at a given phase. We demonstrate that SPAMMS can accurately reproduce the morphology of observed spectral line profiles for overcontact systems. Additionally, we show how line profiles of rapidly-rotating single stars differ when taking rotational distortion into account, and the effects that these can have on the determined parameters. Finally, we demonstrate the code's ability to reproduce the Rossiter-Mclaughlin and Struve-Sahade effects., 10 pages, 10 figures, accepted for publication in A&A

HR 6819 was recently proposed to be a triple system consisting of an inner B-type giant + black hole binary with an orbital period of 40d and an outer Be tertiary. This interpretation is mainly based on two inferences: that the emission attributed to the outer Be star is stationary, and that the inner star, which is used as mass calibrator for the black hole, is a B-type giant. We re-investigate the properties of HR 6819 by spectral disentangling and an atmosphere analysis of the disentangled spectra to search for a possibly simpler alternative explanation for HR 6819. Disentangling implies that the Be component is not a static tertiary, but rather a component of the binary in the 40-d orbit. The inferred radial velocity amplitudes imply an extreme mass ratio of M_2/M_1 = 15 +/- 3. We infer spectroscopic masses of 0.4$^{+0.3}_{-0.1}$ Msun and 6$^{+5}_{-3}$ Msun for the primary and secondary, which agree well with the dynamical masses for an inclination of i = 32 deg. This indicates that the primary might be a stripped star rather than a B-type giant. Evolutionary modelling suggests that a possible progenitor system would be a tight (P_i ~ 2d) B+B binary system that experienced conservative mass transfer. While the observed nitrogen enrichment of the primary conforms with the predictions of the evolutionary models, we find no indications for the predicted He enrichment. We suggest that HR 6819 is a binary system consisting of a stripped B-type primary and a rapidly-rotating Be star that formed from a previous mass-transfer event. In the framework of this interpretation, HR 6819 does not contain a black hole. Interferometry can distinguish between these two scenarios by providing an independent measurement of the separation between the visible components., Comment: Submitted to A&A, 13 pages (16 figures and 2 tables); 4 pages supplementary material (4 figures and 4 tables). Comments are welcome

The massive O4.5 V + O5.5 V binary VFTS 352 in the Tarantula nebula is one of the shortest-period and most massive overcontact binaries known. Recent theoretical studies indicate that some of these systems could ultimately lead to the formation of gravitational waves via black hole binary mergers through the chemically homogeneous evolution pathway. By analyzing ultraviolet-optical phase-resolved spectroscopic data, we aim to constrain atmospheric and wind properties that could be later used to confront theoretical predictions from binary evolution. In particular, surface abundances are powerful diagnostics of the evolutionary status, mass transfer and the internal mixing processes. From a set of 32 VLT/FLAMES visual and 8 HST/COS ultraviolet spectra, we used spectral disentangling to separate the primary and secondary components. Using a genetic algorithm wrapped around the NLTE model atmosphere and spectral synthesis code FASTWIND, we perform an 11-parameter optimization to derive the atmospheric and wind parameters of both components, including the surface abundances of He, C, N, O and Si. We find that both components are hotter than expected compared to single-star evolutionary models indicating that additional mixing processes may be at play. However the derived chemical abundances do not show significant indications of mixing when adopting baseline values typical for the system environment., 24 pages, 4 tables, 12 figures, accepted for publication in ApJ

The Kepler Mission has provided unprecedented, nearly continuous photometric data of $\sim$200,000 objects in the $\sim$105 deg$^{2}$ field of view from the beginning of science operations in May of 2009 until the loss of the second reaction wheel in May of 2013. The Kepler Eclipsing Binary Catalog contains information including but not limited to ephemerides, stellar parameters and analytical approximation fits for every known eclipsing binary system in the Kepler Field of View. Using Target Pixel level data collected from Kepler in conjunction with the Kepler Eclipsing Binary Catalog, we identify false positives among eclipsing binaries, i.e. targets that are not eclipsing binaries themselves, but are instead contaminated by eclipsing binary sources nearby on the sky and show eclipsing binary signatures in their light curves. We present methods for identifying these false positives and for extracting new light curves for the true source of the observed binary signal. For each source, we extract three separate light curves for each quarter of available data by optimizing the signal-to-noise ratio, the relative percent eclipse depth and the flux eclipse depth. We present 289 new eclipsing binaries in the Kepler Field of View that were not targets for observation, and these have been added to the Catalog. An online version of this Catalog with downloadable content and visualization tools is maintained at http://keplerEBs.villanova.edu., Comment: 21 pages, 9 figures

The primary Kepler Mission provided nearly continuous monitoring of ~200,000 objects with unprecedented photometric precision. We present the final catalog of eclipsing binary systems within the 105 square degree Kepler field of view. This release incorporates the full extent of the data from the primary mission (Q0-Q17 Data Release). As a result, new systems have been added, additional false positives have been removed, ephemerides and principal parameters have been recomputed, classifications have been revised to rely on analytical models, and eclipse timing variations have been computed for each system. We identify several classes of systems including those that exhibit tertiary eclipse events, systems that show clear evidence of additional bodies, heartbeat systems, systems with changing eclipse depths, and systems exhibiting only one eclipse event over the duration of the mission. We have updated the period and galactic latitude distribution diagrams and included a catalog completeness evaluation. The total number of identified eclipsing and ellipsoidal binary systems in the Kepler field of view has increased to 2878, 1.3% of all observed Kepler targets. An online version of this catalog with downloadable content and visualization tools is maintained at http://keplerEBs.villanova.edu., Comment: 52 pages, 14 figures, aastex

Accurate stellar parameters of individual objects in binary systems are essential to constrain the effects of binarity on stellar evolution. These parameters serve as a prerequisite to probing existing and future theoretical evolutionary models. We aim to derive the atmospheric parameters of the 31 SB2s in the TMBM sample. This sample, composed of detached, semi-detached and contact systems with at least one of the components classified as an O star, is an excellent test-bed to study how binarity can impact our knowledge of the evolution of massive stars. 32 epochs of FLAMES/GIRAFFE spectra are analysed using spectral disentangling to construct the individual spectra of 62 components. We apply the CMFGEN atmosphere code to determine their stellar parameters and their He, C and N surface abundances. From these properties, we show that the effects of tides on chemical mixing are limited. Components on longer-period orbits show higher nitrogen enrichment at their surface than those on shorter-period orbits, in contrast to expectations of rotational or tidal mixing, implying that other mechanisms play a role in this process. Components filling their Roche lobe are mass donors. They exhibit higher nitrogen content at their surface and rotate more slowly than their companions. By accreting new material, their companions spin faster and are rejuvenated. Their locations in the N-vsini diagram tend to show that binary products are good candidates to populate the two groups of stars (slowly rotating, nitrogen-enriched and rapidly rotating non-enriched) that cannot be reproduced through single-star population synthesis. This sample is the largest sample of binaries to be studied in such a homogeneous way. The study of these objects gives us strong observational constraints to test theoretical binary evolutionary tracks., Comment: accepted to A&A Core of the paper 15 pages, Appendix: 41 pages

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