The Tarantula Massive Binary Monitoring VI: Characterisation of hidden companions in 51 single-lined O-type binaries, a flat mass-ratio distribution, and black-hole binary candidates

We aim to hunt for massive binaries hosting a black hole companion (OB+BH) and establish the natal mass-ratio distribution of massive stars at the subsolar metallicity environment of the Large Magellanic Cloud (LMC). We use the shift-and-add grid disentangling technique to characterize the hidden companions in 51 SB1 O-type and evolved B-type binaries in the LMC monitored in the framework of the Tarantula Massive Binary Monitoring (TMBM). Out of the 51 SB1 systems, 43 (84%) are found to have non-degenerate stellar companions, of which 28 are confident detections, and 15 are less certain (SB1: or SB2:). Of these 43 targets, one is found to be a triple (VFTS 64), and two are found to be quadruples (VFTS 120, 702). The remaining eight targets (16%) retain an SB1 classification. Aside from the unambiguous case of VFTS 243, analysed in detailed in a separate paper, we identify two additional OB+BH candidates: VFTS 514 and VFTS 779. Additional black holes may be present in the sample but at a lower probability. Our study firmly establishes a virtually flat natal mass-ratio distribution for O-type stars at LMC metallicity, covering the entire mass-ratio range (0.05<q<1) and periods in the range 0<log P<3 [d]. The nature of the OB+BH candidates should be verified through future monitoring, but the frequency of OB+BH candidates is generally in line with recent predictions at LMC metallicity.


Introduction
The number of LIGO/VIRGO gravitational-wave (GW) detections originating in coalescing neutron stars (NSs) and black holes (BHs) is rapidly growing. The number of merger events is approaching the 100 mark (Abbott et al. 2019;The LIGO Scientific Collaboration et al. 2021), and it is expected to grow by Article number, page 1 of 42 arXiv:2207.07674v2 [astro-ph.SR] 5 Oct 2022 orders of magnitude in the coming years. There is an overwhelming international effort aimed at constraining the pathways and rates of these merger events (see Mandel & Broekgaarden 2022, for a recent review). This endeavour is intricately tied with our understanding of the progenitors of BHs and NSs -massive stars -which are born with masses higher than 8 M .
A central problem in the field of massive stars concerns the impact and incidence of stellar multiplicity (e.g. Paczyński 1967;Vanbeveren et al. 1998;Langer 2012;Sana et al. 2012). A fundamental input for binary population syntheses in this context is the initial mass ratio distribution of the binary components in massive binaries. In principle, direct measurements of the mass ratios of spectroscopic binaries can only be obtained from double-lined spectroscopic binaries (SB2). However, SB2 binaries are, by definition, systems in which both stars can be observed, and they are hence biased towards mass ratios close to unity. Moreover, SB2 binaries are most easily spotted when the spectral lines are not blended. For this to happen, the binary components need to exhibit appreciable radial-velocity (RV) amplitudes, which again introduces a bias towards short periods and high-mass companions. Short-period (P 10 d) binaries specifically entail the danger of having undergone an interaction in the past, as testified by short-period Algol-like systems (e.g. Mahy et al. 2020b,a;Janssens et al. 2021;Sen et al. 2022). In contrast, single-lined spectroscopic binaries (SB1) are dominated by binaries with more extreme mass ratios and longer periods, where a past binary interaction can be excluded in the case of two non-degenerate components. They, therefore, provide a vital extension to SB2 binaries for a derivation of the natal mass ratio distributions.
Furthermore, SB1 binaries contain an additional population: that of post-interaction binaries (de Mink et al. 2014). Of special interest are OB+BH binaries. Such binaries represent an intermediate evolutionary phase that appears in the majority of evolutionary sequences leading to the formation of BH+BH mergers (Belczynski et al. 2002;Giacobbo & Mapelli 2018;Kruckow et al. 2018;Marchant et al. 2019). As such, OB+BH binaries provide crucial constraints on binary evolution, supernova physics, and the presence of BH kicks (e.g. Fryer et al. 2012;Vanbeveren et al. 2020;Sen et al. 2021). Langer et al. (2020) predicted that about 2 to 3% of massive stars should host a BH companion, implying that hundreds of such binaries are lurking in Galactic and extragalactic OB-star catalogues. And yet, such binaries are sparsely known, and the ones observed are primarily high-mass X-ray binaries (HMXBs, Corral-Santana et al. 2016). However, the X-ray irradiation, mass transfer, and tidal interaction in X-ray binaries makes their analysis and interpretation complex (e.g. Giménez-García et al. 2016, Hainich et al. 2020. The strong tidal circularisation in these binaries also blurs information stored in their orbit regarding the ejecta mass and kick of the BH progenitor. Identifying X-ray quiescent O+BH binaries is thus highly desirable. Indeed, the HMXB population is expected to represent just a small fraction of a much larger population of X-ray quiescent binaries hosting a BH (Langer et al. 2020;Sen et al. 2021).
A few reports of quiescent OB+BH systems have recently emerged (e.g. MWC 656: Casares et al. 2014;LB-1: Liu et al. 2019; HR 6819: Rivinius et al. 2020;NGC 1850BH1: Saracino et al. 2022NGC2004#115: ?;HD 96670: Gomez & Grindlay 2021). However, most of these discoveries were soon thereafter challenged (e.g. Abdul-Masih et al. 2020;Shenar et al. 2020;Bodensteiner et al. 2020;Frost et al. 2022), and the presence of a BH in MWC 656 and HD 96670 still requires validation. In the upcoming year, Gaia is expected to uncover hundreds of OB+BH binaries in our Galaxy relying on accurate astrometry (Breivik et al. 2017;Mashian & Loeb 2017;Janssens et al. 2022). However, to uncover these populations at subsolar metallicities in neighbouring galaxies such as the Small and Large Magellanic Clouds (SMC, LMC), spectroscopy of SB1 binaries is the method of choice. Shenar et al. (2022) identified VFTS 243 in the LMC as an O+BH binary as part of the study put forth here. Recently, a Galactic O+BH system was also reported (HD 130298, Mahy et al. 2022, accepted to A&A).
Owing to its well constrained distance (d = 49.97 kpc, Pietrzyński et al. 2013), modest reddening (Maíz Apellániz et al. 2014), ongoing star formation, and sub-solar metallicity content (Z ≈ 0.5 Z ), the Tarantula nebula (30 Doradus) in the LMC provides an ideal laboratory for investigating massive binaries. It hosts a rich population of massive OB-type stars (Walborn & Blades 1997;Walborn et al. 2014), including the most massive stars and binaries reported to date (de Koter et al. 1997;Crowther et al. 2010;Hainich et al. 2014;Shenar et al. 2017Shenar et al. , 2021Bestenlehner et al. 2022;Brands et al. 2022). The VLT-FLAMES Tarantula Survey (VFTS, Evans et al. 2011) provided a homogeneous, multi-epoch spectroscopic dataset aimed to establish the stellar and binary properties of the massive-star content in this region. Among the ≈ 800 targets included in the VFTS, 360 were classified as O-type (Evans et al. 2011;Walborn et al. 2014). Of those, 116 (32%) were found to exhibit large radial velocity (RV) variations (∆RV > 20 km s −1 ), making them excellent binary candidates . The Tarantula Massive Binary Monitoring (TMBM, PI: Sana, programme IDs: 090.D-0323, 092.D-0136) campaign acquired 32 additional epochs for the majority of the O-type binary candidates, as well as a handful of evolved B-type stars allocated to spare fibres. In total, 93 binaries were studied by Almeida et al. (2017), resulting in orbital solutions for 31 SB2 and 51 SB1 systems. The SB2s were recently analysed by Mahy et al. (2020a,b).
This paper aims to characterise the hidden companions in a sample of O-type and evolved B-type SB1 binaries in the TMBM sample utilising multi-epoch spectroscopy (Sect. 2). Our goal is twofold: to probe the lower-end of the mass ratio distribution, and to hunt for OB+BH candidates. This is achieved through the implementation of spectral disentangling, described in Sect. 3 and Appendix A. The results are presented in Sect. 4 and discussed in Sect. 5. We provide our final conclusions regarding the mass-ratio distribution and OB+BH population in this sample in Sect. 6.

Sample and observations
Our sample comprises the 51 SB1 binaries identified by Almeida et al. (2017) in the TMBM sample. In Table 1, we provide the list of the 51 targets, along with their previous spectral types (SpT), new SpTs derived from the disentangled spectra (Sect. 3), and the orbital periods. We note that Almeida et al. (2017) defined SB1 binaries as systems for which RVs could not be retrieved for a secondary source, which could either mean that the secondary is too faint to be detected, or that severe line blending prevents a trustworthy measurement of the RVs. Hence, for a few targets in the sample, the original SpTs by Almeida et al. (2017) already suggest the presence of a non-degenerate secondary.
Our analysis relies on 32 epochs of optical spectra acquired with the Fibre Large Array Multi Element Spectrograph (FLAMES/GIRAFFE) mounted on UT2 of the Very Large Telescope (VLT). The data and their reduction are described in detail in Almeida et al. (2017); here, we only repeat the essen- Notes. Provided are VFTS identifiers, previous spectral types (adopted from Almeida et al. 2017), new spectral types from our study, and orbital periods (from Almeida et al. 2017, possibly updated in our study). We also provide the number of epochs used, mean S/N of spectra, mass function, estimates for the mass of the primary from spectral-type calibration and evolution tracks, and the mass range of the companion M 2 . (a) Almeida et al. (2017), unless revised here , (b) Walborn et al. (2014) , (c) Derived using disentangled spectra. colon (:) signifies uncertain classification, see Sota et al. (2011) and Walborn et al. (2014) for meaning of remaining qualifiers , (d) Flux ratio of the primary in the visual (see Sect. 3.3). Typical errors are 10%. , (e) Unless stated otherwise, using spectral-type calibrations from Doran et al. (2013) and available photometry (see text for details). Representative errors are 2000 K (1 spectral type ordinate) and 0.3 dex (including uncertainties on SpT, luminosity class, light ratios, and reddening) on T eff and log L, respectively. , (f) Mean evolutionary mass of all apparently-single stars with similar spectral types in Schneider et al. (2018). , (g) Evolutionary mass obtained with the BONNSAI tool (Brott et al. 2011;Köhler et al. 2015;Schneider et al. 2014) using T eff and log L as input. , (h) M min,2 is evaluated from mass function using M 1,SpT . Errors represent 1σ (68% confidence interval). , (i) X-ray detection with the Chandra Visionary programme T-ReX . Targets with no entry were excluded from the X-ray analysis, either because they have log L < 5.0 [erg s −1 ], or due to complex diffuse background contamination. , (j)  naries, we also use five additional epochs from the original VFTS campaign, described in Evans et al. (2011). Table 1 provides the total number of epochs and average signal-to-noise ratios (S/N), which are typically of the order of 50-100 per wavelength element.
In Table 1, we also provide the mass functions f for each systems. Their calculation is primarily based on the results obtained by Almeida et al. (2017), though we use updated orbital parameters obtained in our study when applicable. The estimated flux ratios of the primaries in the visual f 1 / f tot (V) are also given (see Sect. 3.3). Table 1 also lists the effective temperatures T eff,1 and bolometric luminosities log L 1 of the primaries. Generally, the effective temperature is obtained using SpT calibrations by Doran et al. (2013). For log L 1 , we use bolometric corrections from Doran et al. (2013), V-band magnitudes and B − V colours from Evans et al. (2011), and reddening obtained from the difference between the observed colours and SpT-dependent colour (typically B − V 0 ≈ −0.31 ± 0.01, Conti et al. 2008), assuming a relative extinction of R V = 3.5 . We also account for the flux ratio of the primary. We implicitly assume here that the observed colours are not affected by a potential secondary source, which is valid given the early spectral types of the primaries and their companions. For the few targets that were analysed as single stars in previous studies (Bestenlehner et al. 2014;McEvoy et al. 2015;Ramírez-Agudelo et al. 2017), we adopted previously published T eff and log L unless the secondary is found to contribute more than 20% to the visual flux.
We compute two independent estimates for the mass of the primaries: once based on SpT calibrations (M 1,SpT ) and once based on evolution models (M 1,ev ). To estimate M 1,SpT , we calculate the mean of the evolutionary masses of all apparently-single stars in the VFTS sample (Schneider et al. 2018) sharing a similar spectral types (i.e. differing by not more than one subtype and luminosity class ordinate). For the evolutionary masses, we use T eff and log L estimates as input in the BONNSAI tool 1 (Schneider et al. 2014) using tracks from Brott et al. (2011) and Köhler et al. (2015). As priors, we used the Salpter initial mass function and the initial rotational velocity distribution from Ramírez-Agudelo et al. (2013). While some discrepancies between M 1,SpT and (M 1,ev ) are observed, they all agree within 1σ, and no obvious trends are present. This is not entirely surprising, since we are comparing evolutionary masses computed by two different means, which suffer similar systematics (unlike comparing spectroscopic or dynamical masses, e.g. Herrero et al. 1992) . As M 1,ev suffers from more assumptions than M 1,ev , we adopt M 1 = M 1,SpT for the analysis, but using M 1,ev would not impact our conclusions. From the binary mass function and the estimates on the masses of the primaries, constraints on the mass of the secondary M 2 can be derived, which are also provided in Table 1. The derivation of these constraints is described in Sect. 4. We note that low-mass bloated stripped stars could be contaminating the sample (Irrgang et al. 2020;Shenar et al. 2020), and to rule this out, spectroscopic analyses will be needed. We discuss this in Sect. 5.2.
Finally, Table 1 provides constraints on intrinsic X-ray luminosities obtained with the Chandra Visionary programme T-ReX Crowther et al. (2022). Evidently, only four targets (VFTS 333,404,532,603) have finite detections, all of which smaller than 10 33 erg s −1 and lying below or close to the canonical value L x /L ≈ −7 (Pallavicini et al. 1981;Sana et al. 2006). This stands in contrast to the only HMXB candidate in the 30 Dor region, 1 https://www.astro.uni-bonn.de/stars/bonnsai/ VFTS 399 (log L X = 34.7 [erg s −1 ], Clark et al. 2015). The few X-ray detections are discussed in Appendix B.
In addition to the spectroscopy, we also make use of I-band photometry obtained between Oct 2001 and Mar 2020 with the Optical Gravitational Lensing Experiment data reductions III and IV (OGLE, Udalski 2003;Udalski et al. 2008Udalski et al. , 2015.

Spectral disentangling
To investigate whether non-degenerate secondaries are present in our sample, we utilise the method of spectral disentangling (Bagnuolo & Gies 1991;Simon & Sturm 1994;Hadrava 1995;Ilijic et al. 2004;Sablowski et al. 2019). By exploiting the information stored in the spectral-line variability of a spectral time series of a binary, disentangling allows one to separate the component spectra in a binary and derive the orbital parameters. There are different algorithms of spectral disentangling, operating in either wavelength space or Fourier space. Here, we utilise the iterative shift-and-add scheme (Marchenko et al. 1998;González & Levato 2006;Mahy et al. 2010;Shenar et al. 2017). This technique generally yields comparable results to Fourier disentangling (e.g. Shenar et al. 2020), and, in some cases, appears to be more robust (e.g. Bodensteiner et al. 2020). The method of spectral disentangling is not limited to a specific types of stars: it has been implemented on O-type stars (e.g. Martins et al. 2012;Mahy et al. 2020b), B and Be stars (Hensberge et al. 2000;Bisikalo et al. 2000;Saad et al. 2005), lower-mass stars (e.g. Maceroni et al. 2014), and even Wolf-Rayet stars (e.g. Shenar et al. 2017Shenar et al. , 2018Shenar et al. , 2019Shenar et al. , 2021. The underlying assumption is that the variability of the spectral lines is dominated by Doppler shifts of one or more sources (see also Appendix A.) The spectral separation procedure assumes knowledge of the RVs of both components, RV i and RV 2,i , where i ∈ 1, ..., N runs over all available epochs of observation. Equivalently, the orbital parameters can be assumed, from which the RVs may be computed. The procedure is iterative. Starting with an approximation for the secondary's spectrum at the j th iteration, B j , the approximation for the primary's spectrum at the j th iteration, A j is computed by shifting-and-subtracting B j from all observations according to the available RVs, and co-adding the residual spectra in the frame of reference of the primary. The same procedure is then used to compute B j+1 from A j . The first iteration assumes a flat spectrum for the (fainter) secondary, unless otherwise stated, such that the first approximation for the primary is simply the co-added spectrum. The convergence rate of the separation algorithm and hence the number of iterations depends on the RV amplitudes and line profiles, and may range between a few to hundreds, depending on the line profiles, RV amplitudes, and desired level of accuracy.
In many applications, however, the RVs of one or both components are not constrained. In this case, one can perform the disentangling across the K 1 , K 2 axes (grid disentangling), where K 1 , K 2 are the RV semi-amplitudes of both components. In this context, specific spectral lines (or sets of lines) are disentangled. By computing χ 2 (K 1 , K 2 ) from the residual spectra: one can minimise χ 2 to retrieve the K 1 , K 2 values that best reproduce the data, and separate the composite spectra with these Article number, page 4 of 42 T. Shenar et al.: Characterisation of hidden companions in SB1 binaries in the Tarantula values. In Eq. (1), N λ is the number of wavelength bins in the selected range, N epochs is the number of epochs (cf. Table 1). A i,k , B i,k are the disentangled spectra obtained for K 1 , K 2 , evaluated at λ k , and shifted to the appropriate RVs of the i th epoch. O i is the observed i th spectrum, σ i is the S/N in the continuum region in the vicinity of the spectral line, and ν = N λ (N epochs − 2) is the number of degrees of freedom, since each pixel in each of the two disentangled spectra is considered a free variable. Naturally, one may extend this to other orbital parameters, such as the eccentricity (e) or the orbital period (P), albeit this becomes computationally expensive for the shift-and-add algorithm. However, since these parameters are typically well constrained in SB1 binaries, they can be kept fixed in our study, unless clear deviations are seen that require a revision of the parameters.
It is well known that the shift-and-add technique often results in apparent "emission wings" in the disentangled spectra, which are not of astrophysical origin, but are intrinsic to the method (Quintero et al. 2020). Removing these features typically requires very large number of iterations, making the disentangling procedure inefficient. A workaround is to enforce the disentangled spectra to lie below the continuum for absorptionline spectra, unless indications exist that the spectrum contains emission (Shenar et al. 2019;Quintero et al. 2020). This assumption is adopted throughout the paper, unless otherwise stated.
The method is thoroughly tested by creating mock data of different synthetic binaries, discussed in detail in Appendix A. Generally, the method works well in terms of reproducing the input spectra and RV semi-amplitudes whenever the semiamplitudes are larger than half the resolution element (K 1 , K 2 20 km s −1 ). Fortunately, this holds for the majority of targets. Depending on the S/N, we can retrieve companion spectra down to flux contributions of f 2 / f tot ≈ 1 − 3%, noting that the limit strongly depends on the data quality and the rotational velocity of the companion. In principle, pushing to fainter companions could have been possible using the strong Balmer lines. However, the strong nebular contamination in the TMBM data often make Balmer lines unusable for this purpose (see Sect. 3.4).
Finally, we highlight a few facts that appear to hold generally: (i) it is easier to detect the presence of a companion than to determine K 2 . Therefore, in some cases, we can confidently conclude the presence of a non-degenerate companion even if we cannot determine K 2 . (ii) Due to the large number of degrees of freedom, errors on K 1 , K 2 from disentangling are larger than those obtained from orbital fitting. (iii) The fainter the secondary, the less of an impact K 2 has on the disentangled spectra. In contrast, even small deviations in K 1 can result in spurious features in the disentangled spectrum of the secondary. This typically results in a secondary spectrum that mimics that of the primary, which can lead to an erroneous detection of a non-degenerate companion. We critically assess these facts throughout our analysis, and discuss this for the individual targets in Appendix B.

Identifying non-degenerate companions
In the absence of eclipses, the disentangled spectra can only be retrieved up to a scaling factor that depends on the light contribution of the two components to the total flux, l 1 = f 1 /( f 1 + f 2 ) and l 2 = 1 − l 1 = f 2 /( f 1 + f 2 ). The fainter the companion is, the larger the amplification factor 1/l 2 is. Fortunately, the spectral typing does not depend on these scaling factors.
It is important to realise that the method of disentangling, by nature, always returns a disentangled spectrum for a secondary star. It is then up to the user to decide whether the spectrum is "flat", or whether it is of stellar origin. Whether or not the companion is classified as "non-degenerate" is left as a subjective decision, based on a visual inspection of its spectrum after applying the scaling, as well as comparison to spectral models. We generally avoid relying on Balmer lines due to the strong nebular contamination (Sect. 3.4), unless stated otherwise in Appendix B. We give much weight to the presence of He i and He ii lines (depending on the level of nebular contamination). In most cases, when a non-degenerate companion is reported, the stellar spectrum of the secondary is evident, especially upon magnification through re-scaling.
In quite a few cases, it is not clear whether the observed features are of stellar origin, or whether they are caused by nebular contamination (Sect. 3.4) or sources of intrinsic variability (Appendix A). In these cases, we mark the systems as "uncertain SB2" or "uncertain SB1" (SB2:, SB1:), depending on the "degree of belief" we attribute to the spectral features of the secondary, such as their line profiles and widths, their appearance with respect to those of the primary, and the overall level of nebular contamination. Our reasoning for each target is individually discussed in Appendix B. While automated criteria could be developed to identify and classify the binaries, the presence of nebular lines makes it difficult to establish suitable objective criteria. Given the limited sample size, we prefer to keep the classification subjective, since developing an automated algorithm in this context would anyhow boil down to a subjective choice of criteria. The disentangled spectra are presented in Appendix C and are available online 2 for the independent inspection of the reader.

Spectral classification and light ratios
For the spectral classification, we use the Marxist Ghost Buster (MGB) code (Maíz Apellániz et al. 2012 on the disentangled spectra, which allows the user to interactively compare the observed spectra with the spectral type-luminosity class 2-D grid of spectral standards. We use the latest spectral grid derived from the Galactic O-Star Spectroscopic Survey (GOSSS Maíz Apellániz et al. 2011) that includes stars as late as A0 (Maíz Apellániz et al. in prep.). As the GOSSS grid has a spectral resolution 2500, we degrade our spectra to that value.
Spectral classification of hot stars in the blue-violet is usually done in the 3900-5100 Å range. However, here we rely primarily on the disentangled spectra, which only extends down to 4560 Å. The primary luminosity classification criteria and various qualifiers for most O stars involve He ii λ4686 and neighbouring lines, which are only present in the non-disentangled spectra. The 'z' qualifier (He ii λ4686 in absorption and stronger than He i λ4471 and He ii λ4542) cannot be established as that depends sensitively on the ratio of these lines, and we therefore omit it from the classification scheme Maíz Apellániz et al. 2016). As we are also missing the 4630-4660 Å region in the disentangled spectra, the f and c qualifiers (involving the presence of N and C emission lines in this region) cannot be determined (Walborn et al. 2010;Sota et al. 2011). In this case, we adopted these qualifiers from the original classifications of Walborn et al. (2014), implicitly assuming that the features seen in the non-disentangled spectra belong to the primary. We therefore note that the spectral types provided here reflect those that best match our data (disentangled and non-disentangled combined), and may be revised in the future, when disentangling of the full 3900-5100 Å range might become available.
To obtain estimates for the light ratios and the correct scaling of the disentangled spectra, we rely on calibrations between spectral types and equivalent widths of He ii λ4200, 4542, and He i λ4026, 4144, 4388, and 4471, as long as the equivalent widths are significant and not strongly contaminated by nebular lines. We use equivalent-width formulae provided by Sana et al. (in prep.), which are reminiscent of existing quantitative classification schemes by Conti & Alschuler (1971) and Conti & Leep (1974), but using modern data. If two non-degenerate stars are retrieved from disentangling, we compute the light ratio as a weighted mean between the light ratios obtained from the primary's scaling and secondary's scaling. If the spectral subtype of the secondary cannot be determined, we adopt a light ratio of 5% for putative non-degenerate companions and 3% otherwise.

Impact of nebular lines
Only a limited number of spectral lines are available in the data with sufficient S/N to be useful for the disentangling: the Balmer lines Hδ and Hγ, the He i λ4026, 4144, 4388, and 4471 lines, the He ii λ4200, 4542 lines, and rarely metal lines such as Mg ii λ4481 and Si iii λ4553, due to the limited S/N. Of those, only the He ii lines are free of nebular contamination, but as we are mostly hunting for faint B-type companions, He ii lines are of limited use. While sky spectra were subtracted to reduce the nebular contamination, strong residuals remain, likely due to the variability of the emission across the Tarantula. This problem worsens by the fact that the nebular residuals are not constant, but vary from epoch to epoch in strength and shape in a non-trivial manner, presumably due to the varying weather and seeing conditions between the different epochs. At the moderate resolution of our data (∆λ ≈ 0.7 Å ), the nebular lines cover a significant part of the stellar profiles, and cannot be removed by simple interpolation without significantly biasing the data. To partly account for the nebular contamination, we extended the shift-and-add algorithm to three components, where component C represents the nebular lines (see also Abdul-Masih et al. 2019). Component C is enforced to be static, lie above the continuum, and achieve non-vanishing values only in regions overlapping with the nebular emission. To account for the varying strengths of the nebular lines, we scale the component C with the equivalent widths of the emission cores of Hγ.
While our method for removing the residual nebular-line contamination works reasonably well (see examples in Sect. 3.5 and Appendix B), it does not fully remove the impact of nebular lines. This is because the nebular lines vary in shape and strength, and are at times saturated or over-subtracted. This is especially important when the nebular contamination is very strong (typically in Balmer lines and often in the He i λ4471 line), or when the companion is faint, such that the nebular contamination dominates over the secondary's features. Generally, not accounting for nebular contamination leads to spurious features that are especially apparent in the amplified spectra of the faint secondaries. These features can appear confusingly similar to actual stellar features, though they often exhibit strong asymmetries. As we are not able to fully remove the impact of nebular lines, for some systems, our results are ambiguous in terms of the presence of a non-degenerate companion (see Sect. 3.2).

Examples: VFTS 350 (O+O) and VFTS 779 (B+BH?)
We illustrate the method by using VFTS 350 and 779 as examples. VFTS 350 was classified O8 V by Walborn et al. (2014)  reveals a well constrained minimum of χ 2 (K 1 , K 2 ). We perform a similar analysis for all strong helium lines, avoiding He i λ4471 and the Balmer lines due to strong nebular contamination. We obtain the final K 1 , K 2 values by computing the weighted mean of all measurements, amounting to K 1 = 70.3 ± 3.7 km s −1 , and K 2 = 91 ± 20 km s −1 . We note that the K 1 value derived here exceeds the value of 60.2 km s −1 derived by Almeida et al. (2017). This is a typical bias seen when spectral lines are blended by two components (Bodensteiner et al. 2021;Banyard et al. 2022). Figure 3 shows a comparison between the observations and the disentangled spectra for the derived K 1 , K 2 values at two epochs corresponding to the RV extremes. Upon inspection of the line variability in these two epochs, the presence of a nondegenerate companion becomes readily apparent. Evidently, the secondary exhibits broader lines, making it more difficult to spot. The full disentangled spectra, corrected for line dilution, are shown in Fig. 4. Relying on Sota et al. (2011), we classify the secondary as O9.5 V, and estimated its flux contribution at 25%.
We now repeat this exercise for VFTS 779, which was classified B1 II-Ib by Walborn et al. (2014) and revised to B1 II-III here. Almeida et al. (2017) Fig. 6. As Fig. 2, but for VFTS 779. K 1 is measured at K 1 = 35.5 ± 4.0 km s −1 (consistent with the final measurement of K 1 = 30.8 km s −1 ), but K 2 is not constrained.
with P = 59.9 d and a near-circular orbit. A dynamical spectrum is shown in Fig. 5. Like for VFTS 350, a companion is not readily seen in the dynamical spectrum. In Fig. 6, we show the χ 2 (K 1 , K 2 ) map obtained for the He i λ4338 line. While K 1 can be well constrained, K 2 is fully unconstrained. Indeed, when inspecting the spectral variability of this line at RV extremes (Fig. 7), a companion cannot be identified. The spectral appearance of the disentangled spectra is virtually independent of the value of K 2 . As the spectral separation still requires a K 2 value, we adopt the arbitrary value of K 2 = 3 × K 1 = 93 km s −1 (i.e. assuming the companion is three times less massive). In Fig. 8, the full disentangled spectra are shown. To enhance the spectral features of a putative secondary, we adopt an extreme light ratio of l 2 = f 2 / f tot = 0.03 in the visual. Evidently, the disentangled spectrum is virtually flat, with the exception of features originating in diffuse interstellar bands (DIBs) and nebular contamination. VFTS 779 therefore does not have an identified companion. With a minimum secondary mass of M 2 > 3.9±0.6 M (Table 1), this makes VFTS 779 a prime OB+BH candidate. However, given the that the system can accommodate a secondary as light as ≈ 3 M , we cannot rule out that it is a faint non-degenerate companion (main sequence star or helium star), which would contribute less than 1% to the flux. Indeed, VFTS 243 is the only candidate for which non-degenerate companions can be ruled out

Light curve analysis
The OGLE light curves are used to search for photometric variability in the systems, with focus on eclipses or ellipsoidal variations. We analysed available light curves using the Period04 (Lenz & Breger 2005)  Schwarzenberg-Czerny 1997), and the light curves are phased using these periods and inspected. If a notable signal is observed, we present the light curve and discuss it in Appendix B. While not fundamental for our study, the light curves help uncover several eclipsing binaries and high-order multiples. An in-depth analysis of the light curves is beyond the scope of our study.

Results
We perform grid disentangling (1D or 2D) for all targets in the sample. A detailed account of our analysis for each system is given in Appendix B, and the scaled disentangled spectra are presented in Appendix C. Table 2 provides the final orbital parameters, mass ratios q, estimated secondary masses M 2 , and classifications. Figure 9 summarises our findings visually by showing, for each system, the 68% confidence interval on the mass of the secondary, M 2 . The derivation of these mass ranges is described in Sect. 4.1.
Using spectral disentangling, we can recover non-degenerate companions for the majority of our sample. In total, among the original 51 SB1 binaries, non-degenerate companions were found in 43 systems (84%), of which 28 are considered certain and 15 less certain. For the majority of these new SB2 systems, K 2 RV semi-amplitudes could be derived, albeit often with substantial associated uncertainties. The hidden companions are typically early B-type stars, with some exceptions. As Fig. 9 illustrates, almost all companions with minimum masses above ≈ 5 M could be recovered through disentangling. In the remaining eight SB1 binaries (16%), non-degenerate companions could not be uncovered, including the confirmed O+BH binary VFTS 243 ) and two additional OB+BH candidates, VFTS 514 and VFTS 779. We discuss this sample in Sect. 5.2.

Derivation of mass ratios and mass ranges
For systems in which K 2 could be derived, the computation of q and M 2 along with their 68% confidence intervals follows directly from the inferred mass ratio (K 1 /K 2 ) and the primary mass, for which we adopt M 1 = M 1,SpT with the associated errors. For systems for which K 2 could not be constrained (including all SB1 systems), the computation is more involved. First, we asses the possible range of orbital inclinations for these systems. We assume a maximum inclination of 90 • for all targets. The minimum inclination i crit is estimated by equating the measured equatorial rotation velocity sin i/ sin i crit with the critical rotation velocity eq,crit of the luminous components . For simplicity, we adopt a typical value of the critical equatorial velocity as 500 km s −1 (e.g. Townsend et al. 2004), such that: While we could estimate the critical rotation individually, the errors in spectral-type calibration and disentangling procedure easily dominate over the variance of crit over the relevant spectral range. To estimate sin i, we measure the full-widths half-maxima (FWHM) of the He i λ4388 and He ii λ4542 lines of all spectra through Gaussian fitting, and use sin i-FWHM calibrations provided by Ramírez-Agudelo et al. (2015) to obtain sin i. We only accept measurements for which the peak of the Gaussian exceeds the S/N by 6σ (i.e. the half-maximum In total, 28 out of the 51 SB1 systems in our sample are unambiguously found to host non-degenerate companions, though 31 symbols are plotted due to the presence of higher-order multiples. For the tertiary of VFTS 64, the outer period and inner binary mass are denoted. For the quadruples (VFTS 120 and 702), only the binary pairs are considered, since the outer period is not known, and a primary mass of M 1 = 20 ± 10 M is assumed for VFTS 702B for plotting purposes. The symbols mark the medians of the distributions, and the error bars denote the 68% confidence intervals (see Sect. 4.1). Lower Panel: as upper panel, but for the SB1 and uncertain SB1 and SB2 systems (SB1:, SB2:), which could in principle host a BH. Shaded region depicts the mass regime M 2 > 3 M in which BHs are expected to reside. Along with the confirmed O+BH binary VFTS 243, two other OB+BH candidates are highlighted (VFTS 514 and 779).
is at least 3σ below the continuum). We provide the computed sin i and i crit values in Appendix D. We note that the method neglects other sources of star-dependent broadening mechanisms (pressure broadening, micro-and macroturbulence), such that the sin i values should only be considered approximate. The method generally results in poor lower bounds on the inclination of the order of 10 • , and is therefore hardly constraining. We note that this method may not be valid for binaries containing degenerate secondaries due to kicks, which could cause the stellar spin to mis-align with the orbital spin. However, since the method only impacts the upper bound on the mass, it does not have important consequences in the context of identifying OB+BH candidates.
With the mass functions, estimated primary masses (M 1 = M 1,SpT ), and inclination ranges, a probability density for q can be computed for systems without a K 2 measurement. To do so, we perform a Monte Carlo simulation. We ignore errors on the mass function, which are negligible in this context. We sample a Gaussian distribution for the mass of the primary using the values in Table 1. We sample the inclination in the range i ∈ [i crit , 90 • ] (see Eq. 2), following the distribution P(i) ∝ sin i, which re-sults from the projection of randomly oriented orbital orientations. For given i and M 1 , the mass of the secondary M 2 and the mass ratio q can then be computed. For each SB1 system (or SB2 for which K 2 is not known), we thus obtain probability densities P(M 2 ) and P(q). Using these probability densities, we calculate the mediansM 2 andq and the 68% confidence intervals. These values are given in Table 2, and shown in Fig. 9.

High-order multiples
By combining spectroscopic and photometric information, three systems are found to be higher-order multiples. Hence, Table 2 provides 56 entries in total instead of 51. VFTS 64 is a triple system consisting of an outer O7.5 II tertiary orbiting an eclipsing O6+O6 contact binary (P out = 903 d, P in = 1.3 d). VFTS 120 is found to be a quadruple system comprising an O9.5 IV + B0 V binary with a spectroscopic period of 15.65 d and an eclipsing B+B binary with an photometric period of P = 17.41 d. Finally, VFTS 702 is found to be a quadruple system owing to the presence of two distinct eclipse periods (1.98 d and 2.93 d), but the quality of the spectra does not enable us to determine the spectroscopic properties of the system unambiguously. Since the outer periods of both quadruples could not be determined, Fig. 9 includes 54 entries in total.
Given that high-order multiples were only found owing to eclipses, it is likely that additional ones are present in the sample. The FLAMES and OGLE angular resolutions (0.6" and 0.26", respectively) at LMC distance amount to separations of the order of tens of thousands of au, and orbital periods extending up to log P 8 [d]. In this range, the high-order multiplicity fraction is expected to be larger than 50% (Moe & Di Stefano 2017). However, spectroscopy and photometry are not efficient techniques for establishing the multiplicity of stars in the range log P 3 − 4 [d], where instead interferometry or direct imaging are needed, both of which are challenging in the LMC. We therefore cannot rule out that higher-order multiples contaminate our results. However, bright companions that do not trace the orbital period tend to result in K 2 amplitudes of the order of 0 km s −1 (see, e.g. VFTS 64, Appendix B), and such rare cases are discussed in Appendix B. Relatively faint companions would be interesting to find, but have no implication on our analysis. Hence, it is unlikely that higher-order multiplicity impacts our results.

Mass ratio distribution
In this section, we combine all available constraints on mass ratios coming from the homogeneous VFTS and TMBM surveys to reassess the power index κ of the mass ratio distribution, modelled in the form f (q) ∝ q κ (with 0 < q ≤ 1). As discussed in Sect. 1, while SB2 systems allow for direct measurement of q = K 1 /K 2 , they are biased towards short periods, high component masses, and mass ratios close to unity. Disentangling is helpful in this sense, since it extends the sensitivity to wider periods and more extreme mass ratios. Hence, the results become less bias-dominated when disentangling is used. However, even with disentangling, the bias is not fully removed, since a substantial fraction of the SB1 systems for which no companions were retrieved are likely to host faint non-degenerate companions with even more extreme mass ratios. Finally, given that the binary detection probability becomes low for extreme mass ratios, one can expect a population of binaries with two non-degenerate com-A&A proofs: manuscript no. main Table 2. Orbital parameters and companion mass ranges for the 51 O-type SB1 binaries considered in this work.   Notes. Provided are VFTS identifiers, new spectral types from our study, orbital parameters, mass ratios, probability of the system to host a black hole assuming the companion is a compact object, estimated secondary mass, and individual comments. See text for details. (a) Derived in our study . (b) Derived by Almeida et al. (2017), unless otherwise stated . (c) Values and errors represent the median and 68% confidence intervals of the respective probability density function (see text for details) . (d) The probability the companion is a black hole assuming it is a compact object, computed as the probability that M 2 > 3 M (see Sect. 5.2). Probabilities for uncertain SB1 systems (SB1:), rounded to nearest integer, are provided in parentheses, and are omitted for uncertain SB2 systems (SB2:). . (e) Adopted based on symmetry of light curve ; (f) fixed panions that were classified as presumably single stars in the original VFTS survey ). These undetected binaries need to be accounted for if we are to obtain an unbiased estimate for κ.
In Fig. 10, we show a histogram of the mass ratios measured for SB2 systems in the TMBM sample (Almeida et al. 2017) and newly uncovered SB2 systems in this study. We omit VFTS 603 due to its uncertain nature (see Appendix B), and we omit VFTS 243 due to it being identified as an O+BH system  Full population P > 10 d P > 100 d 10 < P < 100 d 0 < P < 10 d Fig. 10. Results for the bias-corrected mass-ratio distribution of the Otype binaries in the Tarantula region. Upper panel: histograms of the mass ratio distributions of the known SB2, new SB2, SB1, and hidden binaries (see legend). Lower panel: corresponding probability densities of the exponent κ of the mass ratio distribution ( f (q) ∝ q κ ). The analysis is perform for various period cuts (see legend). . We further omit the five SB2 systems classified as post-interaction binaries by Mahy et al. (2020b), as we are interested in the natal mass ratio distribution. In addition, we show the median q values for the SB1 systems, derived as describe in Sect. 4.1. The final population we need to consider comprises non-degenerate binaries that were classified as single stars. To visualise this population in the histograms in Fig. 10, we create a mock population of such binaries by using the binary detection probability function provided by Sana et al. (2013) as a function of q, f det (q). We digitise the function, extend it down to q = 0, where per definition f det (0) = 0, and fit a polynomial through the data points (see Appendix E). For every q bin, the number of undetected binaries is given by 1/ f det (q) × M(q), where M(q) is the number of detected binaries in the respective bin.
The distributions match the expected behaviour of the different samples. The known SB2 systems are preferentially found in the range 0.6 < q < 1, while the new SB2 systems, obtained from disentangling of the SB1 systems in this study, probe more uniformly the range 0.2 < q < 1. The "true SB1" binaries (which also include the potential SB2 systems for which K 2 could not be constrained) cluster around more extreme mass ratios of the order 0 < q < 0.5, which is reasonable since a companion was not detected in them. Finally, the undetected binaries have the highest relative contribution for the lowest q values, as could be anticipated.
Assuming a power-law in the form f (q) ∝ q κ , we now estimate κ in a Bayesian approach by computing the posterior of κ given the set of measurements {q i } as We assume a flat prior on κ, such that P(κ) amounts to a normalisation constant, as does the probability for observing the data P({q i }). We account for the errors on q i by multiplying the mass-ratio distribution q κ with the probability distributions of each q i , P(q i ). The population of undetected binaries is accounted for via multiplication with the binary detection probability, f det (q). The posterior P(κ | {q i }) is thus given by: where i runs over the populations of known SB2, new SB2, and SB1 systems. For SB2 binaries, P(q i ) is modelled as a Gaussian using the values given in Table 2. For SB1 binaries, P(q i ) is directly obtained from our Monte Carlo simulations described above.
The lower panel of Fig. 10 shows the results obtained for P(κ|q i ). The result obtained for the entire sample is κ all = 0.2 ± 0.2, consistent with a uniform distribution (κ = 0). However, we also perform the analysis for various period cuts between 1 and 1 000 d to investigate a possible dependence of κ on the period interval, as advocated by Moe & Di Stefano (2017). While most period cuts yield consistent results in the range −0.2 < κ < 0.2, the group of short-period binaries with P < 10 d exhibits a distinctively different distribution, with κ = 1.2 +0.6 −0.5 , implying that very short period binaries tend to have mass ratios closer to unity. A possible explanation or this result could be that many such short-period binaries at more extreme mass ratios (q ini 0.5) have already interacted and merged (Soberman et al. 1997;Langer et al. 2020), and are therefore not included in our sample. However, it is also possible that the star-formation process at very short orbital separations tends to form equal-mass binaries. This fundamental question should be explored in detail in the future.
The unbiased κ values derived here are consistent with a uniform mass ratio distribution (κ = 0). This is also evident when comparing the observed cumulative distribution of the sample with that of power laws for different κ values (Fig. 11). A flat mass-ratio distribution is in contrast with the value of κ = −1.0 ± 0.4 provided by Sana et al. (2013) using the same sample. However, Sana et al. (2013) only had access to RV variability information. A glance at their figure 6 reveals that RV diagnostics are poorly sensitive to κ. Hence, the discrepancy could be easily masked by statistical and systematic biases. On the other hand, a flat mass ratio distribution is more in line with the findings of Sana et al. (2012) for Galactic O-type stars (κ = −0.1 ± 0.6), and also better matches results for Galactic O-types stars by Kobulnicky & Fryer (2007). . Cumulative distribution of the mass ratios measured for the SB2 and SB1 population, including bias-correction due to the finite binary detection probability (see text for details), defined in the bins 0.05, 0.15, ..., 0.95 (each representing an interval of ∆q = 0.1). The distribution is compared to cumulative functions obtained for κ = −1, 0, 0.2, and 1, defined in the interval 0.0 ≤ q ≤ 0.95, except for the κ = −1 case, which is defined for 0.05 ≤ q ≤ 0.95 to avoid a singularity. Evidently, the observed distribution closely follows a flat mass-ratio κ = 0 (uniform).
It is possible that a few of the SB1 binaries considered here host neutron stars or black holes, and should have been omitted from the sample of SB1 binaries. In Sect. 5.2, we estimate that ≈ 4 − 5 OB+BH are expected to be present in our sample. Such systems have typical mass ratios of the order of 0.5-1 (Langer et al. 2020). Having omitted the confirmed O+BH binary VFTS 243 , and given our sample size, the impact of a few O+BH binaries is likely negligible. While OB+NS binaries would strongly skew the distribution towards extreme mass ratios, such binaries are unlikely to be part of our sample. The RV amplitudes of OB+NS binaries are typically of the order of 10-20 km s −1 or less. Hence, such binaries are very unlikely to have passed the threshold of ∆RV > 20 km s −1 imposed by the VFTS survey to be identified as spectroscopic binaries.

The OB+BH population
Among our sample, only one system could be robustly shown to be an O+BH binary: VFTS 243 . In principle, all systems for which a non-degenerate companion could not be found are potential OB+BH candidates. However, not at the same likelihood. The likelihood that a binary in the sample hosts a BH can be written as P(BH) = P(BH|co) · P(co), where P(co) is the probability that the companion is a compact object, and P(BH|co) is the probability that the companion is a BH assuming it is a compact object. The latter is relatively easy to compute. The mass of the most massive neutron star known slightly exceeds 2 M (Fonseca et al. 2021), and neutron stars are not expected to exist above 3 M given constraints from their equation of state (Oppenheimer & Volkoff 1939;Akmal et al. 1998;Lattimer & Prakash 2004;Oertel et al. 2017) . While the bulk of known BHs in the Galactic neighbourhood exceed 5 M (Özel et al. 2010;Farr et al. 2011;Corral-Santana et al. 2016), the existence of lower-mass BHs is still debated (Kreidberg et al. 2012;Wyrzykowski & Mandel 2020;Lam et al. 2022). Hence, to reasonable approximation, P(BH|co) ≈ P3 P(M 2 > 3 M ). In Sect. 4.1, we outlined our Monte Carlo method for establishing probability distributions for M 2 and q. We can use the probability densities P(M 2 ) to estimate P3.
Less trivial is the computation of P(co), which is the complementary probability of the companion being a non-degenerate companion (assuming it is a single star). In principle, one can perform simulations to narrow down the possible types of companions (e.g. faint main sequence stars, helium stars; see Shenar et al. 2022 andMahy et al. 2022). One would then need to combine this with assumptions regarding the likelihood of these companions, which involves the initial mass function and binary evolution. Given the non-trivial uncertainties and complications involved, we simply consider P3 as a proxy for P(BH). These probabilities are given in Table 2. For uncertain SB2 systems (SB2:), we provide the values in parentheses, since P(co) is systematically lower in this case.
Aside from the confirmed O+BH system VFTS 243, two SB1 systems approach or exceed a probability of 90%: VFTS 514 and 779. We mark these systems as prime OB+BH candidates, noting that non-degenerate companions could only be fully excluded in VFTS 243. It is also possible that BHs hide among the other five SB1 binaries, but at a lower probability. The 15 uncertain SB2 or SB1 binaries (SB1:, SB2:) could also host BHs -especially those with high P3 probabilities. More data will be needed to clarify this question (see below).
We now estimate the number of OB+BH binaries expected to lurk among the 51 SB1 systems considered here. According to Langer et al. (2020), about 2% of an unbiased population of massive OB-type binaries are expected to host BH companions 3 . There are 360 O-type stars 4 included in the VFTS survey, and hence we expect roughly seven O+BH binaries among them. With typical RV semi-amplitudes in excess of 50 km s −1 , these O+BH binaries are all expected to have been flagged as binaries in the VFTS survey. The TMBM follow-up targeted 69% of the flagged binaries (Almeida et al. 2017), hence there should be 4-5 O+BH binaries in the TMBM sample. Naturally, none of those would have been classified as SB2 by Almeida et al. (2017), and so we can conclude that these 4-5 O+BH binaries would be present in our SB1 sample (8-10%). With three promising OB+BH candidates (one of which is unambiguous) and 20 less likely candidates, our results seem qualitatively in line with predictions. However, we refrain from over-interpreting this result until the OB+BH nature of the remaining candidates can be verified.
It is interesting to note that the orbit of the confirmed O+BH binary VFTS 243 is near-circular, implying negligible kick or ejecta for the BH progenitor . This is in contrast to the recent report of a BH in the eccentric (e ≈ 0.5) binary HD 130298, which otherwise shares a similar orbital period and spectral type to VFTS 243 . The orbit of the OB+BH candidate VFTS 779 is near-circular, while VFTS 514 is eccentric. If all of these systems are truly O+BH binaries, this may give evidence that some BHs are born with strong kicks, while others are not, in line with recent predictions (e.g. Fryer et al. 2012).
Future high-resolution spectroscopy should significantly improve the detection thresholds, since it should allow for a better removal of the nebular lines and thus make the Balmer lines usable. The Balmer lines are key diagnostics to reject the presence of companions contributing as little as ≈ 0.5% to the flux. Moreover, high resolution is vital for binaries with RV amplitudes comparable or smaller to the resolution element of FLAMES (K 1 40 km s −1 ). Additionally, UV spectroscopy should help reject the presence of non-degenerate companions through dilution of P-Cygni lines (e.g. Georgiev et al. 2011;Shenar et al. 2016).

Conclusions
In this study, we characterised the nature of the hidden companions in 51 SB1 O-type and evolved B-type binaries identified in the framework of the Tarantula Massive Binary Monitoring (TMBM) campaign (Almeida et al. 2017). We implemented the shift-and-add disentangling algorithm to separate the component spectra and establish the RV amplitudes of the two components, whenever possible. We also investigated OGLE lightcurves in search for eclipses or other orbital modulations. Below, we summarise our main results: 1. Out of the 51 SB1 binaries, 43 are found to host nondegenerate companions: 28 are considered certain, and 15 are considered less certain. We also find one triple (VFTS 64) and two quadruple (VFTS 120 and VFTS 702) systems, and eight eclipsing binaries. The remaining eight targets retain their SB1 classifications. 2. Combining constraints from our sample with previously known SB2 binaries (Almeida et al. 2017), and accounting for detection bias, we derive κ = 0.2 ± 0.2 for the exponent of the mass-ratio distribution f (q) ∝ q κ in the range 0.05 < q < 1 and 1 < P < 1000 d. This drops to κ = −0.2 ± 0.2 if binaries with periods shorter than 10 d (potential post-interaction products) are removed, and increases to κ = 1.2±0.5 for binaries with P < 10 d, which we propose is a result of binary interaction among the tight binaries. Our results are therefore consistent with a flat natal mass-ratio distribution for O-type stars at LMC metallicity. 3. Through a probabilistic approach, we identify three O+BH candidates: VFTS 243, 514, 779, of which VFTS 243 is confirmed as an O+BH binary ). The other five SB1 binaries (VFTS 225,631,645,743,769) and 15 uncertain SB1/SB2 binaries (VFTS 73,171,184,256,332,409,479,603,619,657,736,750,812,827,829) could also be O+BH binaries, but at a lower probability.
We strongly encourage further investigations of the various important sub-samples uncovered in this study: the OB+BH candidates, the remaining SB1 binaries (including uncertain SB2), the eclipsing binaries, and the higher-order multiples. The VFTS and TMBM samples offer unparalleled laboratories to homogeneously study binary interaction and evolution at subsolar metallicity.
Step by step, the nature of each and every star in this population is illuminated, providing indispensable constraints on evolution models of massive single, binary, and multiple stellar systems at subsolar metallicity.

Appendix A: The applicability regime of the shift-and-add method
In the context of our study, disentangling comprises two main steps: the retrieval of the RV amplitudes (K 1 , K 2 ), and the retrieval of the disentangled spectra. The robustness of the disentangling procedure depends on many factors, including the individual S/N, Doppler-space coverage and number of observations, line profiles, light ratios, and RV amplitudes. The underlying assumption of spectral disentangling, as implemented here, is that each observed spectrum is the sum of stellar spectra (two in the case of a binary) that are Doppler-shifted with the respective RVs of the stars, contaminated only by Gaussian noise and static nebular lines that can vary in strength from epoch to epoch (Sect. 3.4). Hence, intrinsic non-Doppler variability can impact the results, especially if it is (quasi-)periodic. In the case of O-type stars, the most important source of contamination is likely stellar pulsations, which can result in apparent peak-to-peak RV variability of up to ≈ 20 km s −1 (Ritchie et al. 2009;Aerts et al. 2009;Sana et al. 2013). Since all binaries in our sample have K 1 > 10 km s −1 per definition, we expect intrinsic variability to play a secondary role for most targets. Moreover, we note that the procedure of disentangling tends to smear-out features that do not follow the orbital period. We therefore do not expect pulsations to impact the qualitative or quantitative results significantly. An exception is when the dominant pulsation period is mistaken for the orbital period, or when they are otherwise comparable. A notable example is VFTS 603, which has a derived RV semi-amplitude of K 1 = 11.4 km s −1 and a period of P = 1.76 d. The low-amplitude motion, combined with the short period, raises doubts that pulsations cause the underlying variability. Disentangling does suggest the presence of two nondegenerate stars, but it is possible that the second source is a spurious result of non-Doppler variability (see Appendix B).
Fully mapping the impact of all parameters and line-profile variability on the output of spectral disentangling is beyond the scope of our study, and would require a dedicated study with prespecified models of line variability. Critical for our study is to establish threshold light ratios below which the spectrum of a non-degenerate secondary would be discerned from the noise. To estimate it, and to test the applicability regime of our method, we run a set of simulations of synthetic noisy composite spectra that mimic the resolving power, S/N, and timing of the TMBM programme. Specifically, we consider a hypothetical binary with P = 100 d, a circular orbit, and a mass ratio of 0.5 (i.e. K 2 = 2 K 1 ). We use synthetic TLUSTY spectra (Lanz & Hubeny 2003Hubeny & Lanz 1995) for the primary and secondary. For the primary, we adopt T 1 = 35 kK, log g = 4.0 [cgs], and sin i = 200 km s −1 . For the secondary, we assume T 2 = 26 kK, log g = 4.0 [cgs], and sin i = 100 km s −1 . We mimic the data quality of VFTS 231, which has varying S/N for the individual spectra between ≈ 10 and 40, with an average of ≈ 30. VFTS 231 is among the 10% worst-quality targets in terms of S/N, and is chosen to ensure that our estimates are conservative.
We constructed 12 mock binaries in total, spanning the values K 1 = 200, 50, and 20 km s −1 (and hence K 2 = 400, 100, and 40 km s −1 ), and secondary light contributions of l 2 = f 2 / f tot (V) = 0.5, 0.1, 0.03, and 0 (the latter corresponding to a black hole secondary). We then apply the 2D shift-and-add disentangling technique to derive K 1 , K 2 and disentangle the spectra using those values. The results are shown in Fig. A.1. Evidently, the retrieval procedure works well across the parameter space, although discrepancies become apparent for lower RV ampli-tudes and fainter secondaries. Figure. A.1 illustrates that, at the low data quality of the simulations, companions contributing as little as ≈ 5% can be retrieved. This improves to 1 − 3% for higher S/N data available for some of our targets.

Appendix B: Comments on individual targets
This appendix provides a detailed commentary on the individual targets and the nature of the unidentified secondaries. We also discuss the light curves of the targets (when available). Unless stated otherwise, cited periods and eccentricities are taken from Almeida et al. (2017).
VFTS 64, O8 II:(f) + (O7 V + O7 V), was reported to be a long-period, eccentric binary with P = 903 d and e = 0.53. Among our sample, it has the largest mass function, with f (m) = 10.7 ± 0.4 M . Spectral disentangling unambiguously shows the presence of a companion. However, disentangling of various He lines results in unrealistically small K 2 values, and the disentangled spectra of the secondary appear to centre around a different systemic velocity than that derived by Almeida et al. (2017). This suggests that this system may be a multiple.
Available OGLE photometry confirms the triple nature of the system. Figure B.1 shows the phased I-band light-curve of VFTS 64. A period of P = 1.327181 ± 0.000003 is found from the light-curve analysis. This period is clearly different than the one derived from RVs (903 d). The light curve suggests that the secondary is in fact an eclipsing contact system consisting of two nearly-identical stars, such that the system is a triple with an outer period of P out = 903 d and an inner period of P in = 1.33 d.
We extended the 2D disentangling technique to account for three components. Denoting with the subscript '2' the contact binary as a whole, we now need to constrain the RV amplitudes K 3 and K 4 of the components of the contact binary. We fix the orbital parameters derived by Almeida et al. (2017) for the primary, including its RV amplitude of K 1 = 57.2 km s −1 . The period and time of primary eclipse P in , T 0,in are fixed to those found from the light curve analysis (Fig. B.1). Given the properties of the light curve, we adopt e in = 0 and ω in = 90 • , such that T 0,in corresponds to zero RV shift. Assuming that the inclinations of the inner and outer orbits are aligned given the relatively low separation between the binary and the tertiary, it can be shown from Newtonian relations that the RV amplitude of the contact binary itself, K 2 , is obtained via We note that by requiring K 2 > 0 km s −1 , Eq. (B.1) imposes the condition K 3 + K 4 420 km s −1 . We then measured K 3 and K 4 by disentangling the strong He ii and He i λ4471 lines (e.g. Fig. B.2), whose weighted mean yields K 3 = 282 ± 36 km s −1 , K 4 = 323 ± 33 km s −1 , yielding in turn K 2 = 41 ± 13. This unique system could be a prototypical progenitor of the triple WR system BAT99 126, which was also found to host a contact binary (Janssens et al. 2021).
VFTS 73, O9.5 IV + B: has a reported period of P = 151 d and an eccentricity of e = 0.20. The spectral variability does not readily suggest the presence of a non-degenerate companion (Fig. B.3    companion and an assumed a 5% light contribution, but we cannot rule out a BH. Analysis of the OGLE light curve did not reveal significant periods. VFTS 86, O9.5 III + O8 IV:n has a reported period of P = 183 d and an eccentricity of e = 0.51. Inspection of the line profiles is suggestive of the presence of two non-degenerate stars in the system, both of which appear to contribute a significant amount to the visual flux (Fig. B.4). For this reason, we resort to 2D disentangling. We used the strong He i lines (which suffer very little nebular contamination) and the two He ii lines. The lines yield broadly consistent results, albeit with large errors. The smallest errors are obtained from the He i λ4388 line, yielding K 1 = 43 ± 8 km s −1 and K 2 = 48 ± 10 km s −1 . With a derived light contribution of 44%, the secondary appears to be an O8 star, that is, of earlier spectral type than the primary. Analysis of the OGLE light curve did not reveal significant frequencies.  VFTS 93, O8.5 V + B0.2: V has a reported period of P = 250 d and an eccentricity of e = 0.20. The primary moves with a very low RV amplitude of K 1 = 11 km s −1 , making it difficult to properly disentangle the spectra. However, careful inspection of the line-profile variability of strong He i lines at extremes suggests the presence of a non-degenerate companion (Fig. B.5). We re-derived K 1 by measuring the RVs of the He ii lines, which are not contaminated by the secondary, and obtain K 1 = 12.8 ± 1.1 km s −1 . Fixing this value, the strong He i lines all result in minima of the order of K 2 = 10−30 km s −1 . A weighted mean of all measurements yields K 2 = 20 ± 9 km s −1 . The spectrum is dominated by He i absorption, and best corresponds to a B0.2 spectral type, with an estimated light contribution of 21%. Analysis of the OGLE light curve did not reveal significant frequencies.     Fig. B.8. Comparison of He i λ4388 spectra of VFTS 171 at RV extremes, along with the disentangled spectra and their sum, as derived for K 1 = 12 km s −1 and K 2 = 3 × K 1 = 24 km s −1 . The disentangled spectra are not scaled by the light ratio in this plot.
VFTS 120, (O9.2 IV + B) + (O9.5 V + O9.5 V) was reported to be have a period of P = 15.7 d and an eccentricity of e = 0.28. Inspection of the He i lines reveals line-profile variations that are suggestive of a second non-degenerate star in the spectrum (Fig. B.6). For this reason, we implemented the 2D disentangling algorithm. However, the solution appeared to be strongly dependent on the lines used. Moreover, the He i λ4388 and He i λ4471 lines, which are the strongest and appear to be hardly contaminated by nebular emission, result in K 2 amplitudes of the order of 0 km s −1 . As in the case of VFTS 64, this could imply the presence of additional stellar components in the system.
Interestingly, Fourier analysis of the OGLE light curve shows a peak at 8.7 d and its harmonics. Folding the light curve with twice this value, that is, P = 17.41 d, reveals a clear signature of eclipses in the system (Fig. B.7). It can be fully ruled out that this period corresponds to the spectroscopic period. Given the similarity of the periods, an eclipsing tertiary can be ruled out. The only conceivable way to explain this is that VFTS 120 is a quadruple system. For this reason, we extended the shiftand-add algorithm to account for four components. To reduce the parameter space, we adopt for the eclipsing binary the photometric period, an eccentricity of 0 (as implied by the ≈ 0.5 phase difference in eclipses), and the time of periastron (defined here as time of primary eclipse). The four RV amplitudes K 1 , K 2 , K 3 , K 4 are then scanned through a χ 2 minimisation. We find K 1 = 92±13 km s −1 , K 2 = 93±17 km s −1 , K 3 = 20±7km s −1 , and K 4 = 21 ± 7 km s −1 . However, the results should be taken with caution, given the assumptions, limited resolution, and data quality available here. The outer period connecting the two binaries cannot be established here.
VFTS 171, O8.5 III:(f)) + B1.5: V, was reported to be a long-period binary with P = 677 d and e = 0.56. The lowamplitude motion of the primary (K 1 = 12 km s −1 ) prevents us from obtaining unambiguous results for this system, and inspection of the line-profile variability is not readily suggestive of a non-degenerate secondary (Fig. B.8). Disentangling the spectra results in faint signatures in Balmer lines, He i lines, and He ii lines. However, we cannot derive K 2 from the relatively flat χ 2 map. We tentatively classify the companion as B1.5: V and estimate its light contribution at 8%, but given the low RV amplitude, we consider the SB2 nature uncertain. The OGLE light curve shows a variability pattern when phased at the orbital period (Fig. B.9), though it is not clear that this pattern is of really astrophysical origin.  rived by Almeida et al. (2017) shows substantial scatter. The lack of line-profile variability at RV extremes does not readily suggest the presence of a non-degenerate companion (Fig. B.10). Given the scatter in the orbital solution, we attempted 2D disentangling across the K 1 , K 2 axes, and probing all He i and He ii lines, as well as Hδ, which is weakly contaminated by nebular lines. The results of the different lines are consistent, but the errors are very large, with final mean averages of K 1 = 29 ± 11km s −1 and K 2 = 61 ± 29 km s −1 . Disentangling the spectrum for these RV amplitudes results in similar spectra for both components, which sheds doubt on the validity of the results, as this is often the case when the adopted K 1 is significantly different than the true one. When disentangling for K 1 = 12.1 km s −1 , the features in the spectrum of the secondary weaken, though they do not fully disappear. Our results are therefore uncertain. We tentatively classify the companion as OB:, and refrain from using our measurements for any further analysis. With a minimum mass of M 2 > 1.7 ± 0.4 M , the nature of the secondary is virtually unconstrained. Higher-resolution spectra will be necessary to disentangle this system. Analysis of the OGLE light curve did not reveal significant frequencies.
VFTS 191, O9.2 V + O9.7: V: has a reported period of P = 359 d and an eccentricity of e = 0.22, although the orbital solution (especially the eccentricity) are uncertain. The spectral variability does not readily suggest the presence of a non-degenerate companion (Fig. B.11). However, grid disentangling consistently points towards the presence of a relatively bright companion entangled in the spectrum. We are not able to constrain a K 2 value using disentangling of the strong He i lines. However, adopting plausible K 2 values of the order of  2 − 3 K 1 ≈ 50 − 70 km s −1 results in a secondary spectrum that appears to belong to a fainter, rapidly rotating O-type star, contributing roughly 13% to the total flux. Analysis of the OGLE light curve did not reveal significant frequencies.
VFTS 201, O9.7 V + B1.5: V has a reported period of P = 15.3 d and eccentricity of e = 0.46. The variability seen in the He i lines (Fig. B.12) is suggestive of a second non-degenerate component in the system. Due to the lack of strong He ii lines and the apparent contribution of the secondary, we resort to 2D disentangling. All He i lines are in broad agreement with each other, but show substantial errors on the RV amplitudes. Disentangling all strong He i simultaneously, we find K 1 = 79 ± 9 km s −1 and K 2 = 105 ± 26 km s −1 . The disentangled spectrum of the secondary, which has an estimated light contribution of roughly 30%, shows strong He i lines and no He ii lines, with only faint evidence for the Mg ii λλ4481.1, 4481.3 doublet, and it is classified as B1.5: V. Analysis of the OGLE light curve did not reveal significant frequencies.
VFTS 225, B0.7 III has a reported period of P = 8.2 d and a virtually circular orbit. Disentangling of the strong He i lines (Fig. B.13) does not reveal a significant minimum for K 2 , and the disentangled spectrum of the secondary is virtually featureless in He i lines, and shows very little signatures in Balmer lines that appear to originate from inaccurate subtractions of the weak nebular lines. The system is thus a bona fide SB1 system. However, the minimum mass imposed by the orbital parameters and the estimated mass of the primary yield M 2 > 1.8 ± 0.3 M , such that the nature of the secondary remains unconstrained.
The OGLE III data show substantial scatter and a longterm magnitude increase that is not supported by the OGLE IV data. We therefore remove them from the photometric analysis. Fourier analysis of the OGLE IV data reveals a peak at the orbital period. The phased light curve is shown in Fig. B.14. VFTS 231, O9.7 V + B1.5: V has a reported period of P = 7.9 d and eccentricity of e = 0.41. Almeida et al. (2017) noticed the presence of a secondary star in the system, but could not derive its RVs. Indeed, comparing line profiles belonging to strong He i lines such as He i λ4026 (Fig. B.15) clearly shows the presence of a non-degenerate companion. We implemented 2D disentangling and find consistent results for all He i lines. A mean average of the measurements obtained for the strong He i lines (with exception of the He i λ4471 line, which suffers strong nebular contamination) yields K 1 = 99 ± 12 and K 2 = 130 ± 25 km s −1 . From the strengths of the spectral lines, we estimate a light contribution of 33% for the secondary. We classify it as B1.5:. Despite the short period of the system, we cannot find significant periodicities in the OGLE light curve. VFTS 243, O7 V(n)((f)) + BH was shown to be a bona fide O+BH binary with a P = 10.4 d period and a near-circular orbit in Shenar et al. (2022), to which we refer for a detailed analysis of the spectra and light curve of this unique binary.
VFTS 256, O7.5 V: +OB: has a period of P = 246 d and an eccentricity of e = 0.63. Inspection of the line-profile variability at RV extremes does not readily reveal the presence of a nondegenerate companion (Fig. B.16). The low-amplitude motion (K 1 = 19.2 km s −1 ) and nebular contamination make disentangling of the system challenging. Given the potential contamination of the secondary, we performed 2D disentangling, but find a consistent K 1 value to that derived by Almeida et al. (2017), and therefore adopt their value. We cannot find a minimum in the χ 2 (K 2 ) map, and the disentangled spectra provide marginal evidence for the presence of a companion at faint He i lines features and Balmer features, which may however be impacted by nebular contamination. Comparably weak He ii features are seen, thought these features may be spurious due to an inaccurate K 1 value. We tentatively classify this as an SB2:. With a minimum mass of M 2 > 4.1 ± 0.5 M , the secondary could be a faint nondegenerate star or a BH. Analysis of the OGLE light curve does not reveal any periodicity at the orbital period, but does show periodicity at P = 341 d. This periodicity may well be related to the yearly cycle, and hence cannot be confirmed to represent an intrinsic period of the system.
VFTS 277, O9 V + B1.5: V has a period of P = 240 d and a high eccentricity of e = 0.93. The mass function and estimated mass of the primary imply a companion of at least ≈ 3 M . The spectral variability of He i lines weakly suggests the presence of a non-degenerate companion (Fig. B.17 sively. Disentangling reveals significant He i absorption for the secondary but no He ii absorption. However, K 2 cannot be constrained with the current data quality. Based on the appearance of the Hδ region and the Mg ii λλ4481.1, 4481.3 doublet, we tentatively classify the companion as B1.5. However, as the spectral appearance of the secondary depends on the adopted K 2 value (here, we use K 1 = 2 × K 2 = 128 km s −1 ), the classification is not certain. Moreover, the presence of a faint, blue-shifted He ii λ4200 absorption in the spectrum of the secondary implies that this system may be a higher-order multiple. Analysis of the OGLE light curve did not reveal significant frequencies.
VFTS 314, O9.7 V(n) + B has a short period of P = 2.5 d and an eccentricity of e = 0.17. Almeida et al. (2017) noted the presence of a secondary in the individual spectra, but could not measure its RVs. Indeed, disentangling of the He i lines (Fig. B.18) results in a relatively well defined minimum around 200 km s −1 ; the most robust result is found from disentangling of the He i λ4026 line, which yields K 2 = 232 ± 16 km s −1 . The corresponding disentangled spectrum shows a spectrum with clear He i absorption and no He ii absorption, but a precise classification of the secondary is difficult due to the lack of metal lines. We estimate that the companion contributes 5% or less to the visual flux.
As could be anticipated from the short period, VFTS 314 is an eclipsing system. The light curve of VFTS 314 is shown in Fig. B.19. From analysis of the light curve, the period can be refined to P = 2.550786 ± 0.000005. Interestingly, the primary and secondary minima do not yield an identical period, suggesting that the system experiences apsidal motion.  VFTS 318, O9.5 V + O9.2 V is a nearly circular binary with a reported period of P = 14 d. Careful inspection of the line variability at RV extremes of He i lines such as He i λ4026 or He i λ4388 (Fig. B.20) suggests the presence of a second nondegenerate star in the binary. For this reason, we implemented 2D disentangling on various He i lines. A weighted mean of the RVs measured from the He i lines yields K 1 = 53±12 km s −1 and K 2 = 93±24 km s −1 . The He ii lines do not yield well-constrained measurements. The value of K 1 is more than twice the amplitude derived by Almeida et al. (2017), a consequence of accounting for the secondary. The two stars a very similar spectral type, and the secondary is estimated to contribute 33% to the visual flux. It is therefore somewhat surprising that the secondary is found to have roughly half the mass of the primary, suggesting that the primary may be an evolved star. Analysis of the OGLE light curve did not yield significant periods.
VFTS 329, O9.5 V(n) + B1: V: has a reported period of P = 7.0 d and eccentricity of e = 0.44. The spectral-line variability of He i lines suggests the presence of a fainter non-degenerate secondary in the spectrum, as is shown for the He i λ4338 line in Fig. B.21. Given the non-negligible contribution, we perform a 2D disentangling for the system. A mean average of the of strong He i, aside from the contaminated He i λ4471 line, yields K 1 = 93 ± 8, K 2 = 131 ± 24 km s −1 . The disentangled spectrum of the secondary, which has an estimated light contribution of 22%, matches a B1 V spectral type. Despite the relatively short period, analysis of the OGLE light curve did not yield significant periods.
VFTS 332, O9 III + O9.2 V is an eccentric long-period binary with P = 1025 d and e = 0.81. The coverage of the spectra in Doppler space is very poor, with only one epoch showing significant RV shifts compared to the others. Therefore, results  from disentangling of this system should be taken with caution.
The results point to the presence of a second, non-degenerate object in the system, although its presence is not readily seen from the spectral variability (Fig. B.22). Disentangling of all He lines results in a weighted mean of K 1 = 79 ± 10 km s −1 and K 2 = 46 ± 13 km s −1 . As in the case of VFTS 184, the similarity of the disentangled spectra suggests that the secondary spectrum may be contaminated by the primary due to a mismatch of the K 1 value. However, regardless of the K 1 value used for disentangling, a stellar spectrum is always obtained for the secondary, suggesting that it is likely real. Taking our results at face value, it appears that the hidden secondary is the more massive star, with a spectral class of roughly O9. Despite of this, the disentangled spectra imply that it contributes only 24% to the total light, suggesting that the primary is perhaps more evolved. On the other hand, the spectral appearance of the secondary also corresponds to an evolved object. These results combined are hard to reconcile with each other, suggesting that the disentangled spectra suffer from crosscontamination between the components. High resolution spectroscopy, obtained during periastron passage, will be required to disentangle this object unambiguously. Analysis of the OGLE light curve did not yield significant periods.
VFTS 333, O9 II((f)) + O6.5 V: is an eccentric long-period binary with P = 980 d and e = 0.75,. Unlike VFTS 332, the spectra offer a good coverage in Doppler space. A non-degenerate companion can be seen at the RV extreme of strong He i lines (Fig. B.23). Disentangling of the strong He ii and He i lines yields values between 20-40 km s −1 for K 2 , and suggest that the secondary is relatively bright. We therefore resort to 2D disentangling. In all cases, K 1 is found to be of the order of 40 km s −1 , with a weighted mean of K 1 = 43.0±2.5 and K 2 = 34±14 km s −1 .  All measurements suggests that the object previously identified as the primary is in fact the less massive companion. We estimate it to contribute 42% to the flux based on the strengths of the spectral lines. Analysis of the OGLE light curve did not reveal significant frequencies.
VFTS 350, O8.5 V + O9.5 V was thoroughly discussed in Sect. 3.5. Analysis of the OGLE light curve does not reveal any significant frequencies.
VFTS 386, O9 V(n) + B1 V: has a reported period of P = 20 d, an eccentricity of e = 0.25, and a low-amplitude motion of K 1 = 14.3 km s −1 for the primary. However, the orbital solution established by Almeida et al. (2017) exhibits substantial scatter, with a reduced χ 2 of 9.5. While the presence of a non-degenerate companion is not readily clear from the line profile variability, comparison of the variability in He i and He ii lines reveals that this binary hosts two non-degenerate companions (Fig B.24). Given the low amplitude motion and the lack of significant variability in He i lines, we do not proceed with 2D disentangling. Instead, we re-derive the orbital solution of the system by focusing only on the He ii lines. The revised orbital solution yields a reduced χ 2 of 2.4 (Fig. B.25). The orbital parameters are then P = 20.45 ± 0.02 d, T 0 = 54816.5 ± 2.8 [JD -2400000], K 1 = 30.8 ± 2.1 km s −1 , and e = 0.15 ± 0.05. We then proceed with 1D disentangling of the He i lines, fixing the orbital parameters derived above. The best constraints are obtained from the He i λ4388 (Fig. B.26) and He i λ4471 lines, which yield K 2 = 64 ± 12 km s −1 . The spectral type of the secondary corresponds to a B1 star, and we estimate its light contribution to be 16%. Analysis of the OGLE light curve did not yield significant periodicities. VFTS 390, O5.5 V:((fc)) + O9.7: V: has a reported period of 22 d and an eccentricity of e = 0.49. Disentangling of the spectra proved to be very challenging due to very strong nebular-line contamination and the fact that the secondary appears to possess very weak lines (Fig. B.27). While disentangling of various He i lines supports K 2 values in the range 100 − 250 km s −1 , the errors are very large. A weighted mean of the RVs obtained for He i λ4026 and He i λ4388 yields K 2 = 137 ± 34 km s −1 , but the nebular emission and very low S/N may entail systemic errors. The secondary's spectrum matches a O9.7 V star, and has an estimated light contribution of 15%. Interestingly, the system shows a single eclipse (Fig. B.28  degenerate companion in the system, and allows us to slightly refine the ephemeris. From the orbital configuration, the observed eclipse occurs when the hotter, brighter primary eclipses the cooler, fainter secondary. Due to the orbital configuration, a second eclipse is not seen. VFTS 404, O3.5: V:((fc)) + O5 V: has a reported period of P = 146 d and an eccentricity of e = 0.72. Almeida et al. (2017) noted the presence of a non-degenerate companion, but could not derive its RVs. Indeed, the secondary is clearly seen in both He i and He ii lines (Fig. B.29). Given the significant contribution of the secondary to these lines and the fact that the components are constantly blended, it is very likely that the RVs measured by Almeida et al. (2017) for the primary were strongly impacted by the secondary. Therefore, we utilise 2D disentangling for this system. For this, we use the He ii lines, which are free of nebular contamination and are stronger than the He i lines. A weighted mean yields K 1 = 98 ± 11 km s −1 and K 2 = 106 ± 11 km s −1 . We note that the K 1 value derive here is roughly three times larger than derived by Almeida et al. (2017), which is a result of the secondary's anti-phase motion that leads to an apparent lower amplitude of motion. The companion is estimated to contribute 42% to the visual flux, and is classified O5 V. Analysis of the OGLE light curve did not reveal significant frequencies. We note that VFTS 404 is one of the four targets for which X-rays are detected, and exhibits the highest X-ray luminosity in our sample (log L X = 32.84 [erg s −1 ]). This luminosity lies slightly above the canonical value of L X /L ≈ −7, perhaps due to the presence of colliding winds in the system. The disentangled spectrum shows faint signatures in He i lines (Fig. B.30), though nebular contamination questions the validity of these results. If these signatures are real, then the light contribution of the secondary is ≈ 0.05. We tentatively classify the secondary as B:, and the binary as SB2:. With a minimum mass of M 2 > 7.6 ± 1.1 M , the secondary could be an early B-type star, a helium star, or a BH. Analysis of the OGLE light curve did not reveal significant frequencies.
VFTS 429, O7 V: + B1: V: has a reported period of P = 30 d and an eccentricity of e = 0.56. This system is clearly SB2, as is readily seen from the spectral variability of strong He i lines. However, the spectra are strongly contaminated by nebular lines, making the derivation of K 2 challenging. The least contaminated line clearly showing both stars isolated is the He i λ4144 line (Fig. B.31). A weighted mean of the strong He i lines (with exception of He i λ4471 due to very strong nebular contamination) yields K 2 = 143 ± 27 km s −1 . Judging by its spectral appearance, we classify the secondary as B1: V:, and estimate its light contribution in the visual at 10%. Similar to VFTS 409, phasing the OGLE light curve of VFTS 429 with its spectroscopic ephemeris reveals a single eclipse (Fig. B.32). We used it to slightly refine the orbital period. Unlike VFTS 409 however, the eclipse in VFTS 429 occurs when the cooler B0 V secondary eclipses the hotter O7.5 V primary.
VFTS 440, O6: V:(f) + O8 V, is a long-period binary with P = 1019 d and e = 0.28. The low RV amplitude of the pri- mary (K 1 = 11.6 km s −1 ) and the strong nebular contamination make the disentangling of this binary challenging. We first implemented the 2D disentangling of the He ii lines to ensure that the low RV amplitude is not impacted by a secondary contributing significantly to the flux. We find K 1 that agree well with the value derived by Almeida et al. (2017), and we therefore fix it and continue with 1D disentangling across the K 2 axis. While not readily apparent, careful inspection of the He ii λ4542 line at RV extremes suggests the presence of a faint secondary (Fig. B.33). Disentangling of the He ii lines results in poor constraints on K 2 , since all values in the range [0,100] km s −1 result in acceptable χ 2 values (within 1σ). However, disentangling suggests that the He i λ4009 line is present only in the secondary. We therefore use this line for disentangling as well. A weighted mean of the poorly constrained He ii measurements and the He i λ4009 line yields K 2 = 29 ± 18 km s −1 . The RV amplitude of the secondary is thus poorly constrained. We classify the companion as O8 V and estimate it contributes 16% to the total flux. Analysis of the OGLE light curve did not reveal significant frequencies.
VFTS 441, O9.2 V + B0.5 V has a reported period and eccentricity of P = 6.9 d and e = 0.22. Inspection of the He i lines unambiguously reveals this system as an SB2 binary, albeit the strong nebular contamination of the He i lines makes the disentangling challenging. Because of the presence of a nonnegligible secondary, we implement the 2D disentangling technique. The line least contaminated by nebular lines is He i λ4144, shown in Fig. B.34 at RV extremes. A weighted mean of all He i lines but He i λ4471 (due to dominating nebular contamination) yields K 1 = 73 ± 5 and K 2 = 108 ± 8 km s −1 . The disentangled spectrum of the secondary, estimated to contribute 21% to the visual flux, shows no evidence for He ii absorption. Its spectral appearance matches a B0.5 spectral type. Despite of the short or- bital period, analysis of the OGLE light curve of VFTS 441 did not reveal significant periods. VFTS 475, O9.7 V + B0 V has a reported period of P = 4.1 d and an eccentricity of e = 0.57. The spectral variability of He i lines implies the presence of a non-degenerate secondary in the system (Fig. B.35), and we therefore implement 2D disentangling. Plausible constraints are obtained with the He i λ4388 and He i λ4475 lines, albeit with substantial errors due to the low S/N of the data. A weighted mean of the measurements yields K 1 = 135±33, K 2 = 169±58 km s −1 . The He i λ4026 line results in low K 2 amplitudes, which do not agree with the other lines. The reason for this discrepancy is not clear, and could imply that the system is a higher-order multiple. The spectral appearance of the disentangled spectrum off the secondary matches a B0 V spectral type, and we estimate its light contribution at 35%. The OGLE light curve, phased with the orbital period, shows a clear periodic signature, potentially originating in ellipsoidal variations and irradiation effects (Fig. B.36).
VFTS 479, O4.5 V((fc)) + B: has a reported period of P = 14.7 d and an eccentricity of e = 0.31. The spectra are very noisy and contaminated by nebular emission, such that disentangling of the system is challenging (Fig. B.37). Disentangling of He i lines results in K 2 values of the order of 200 km s −1 with very large errors. A weighted mean of the measurements obtained for the He i λ4026, He i λ4144, and He i λ4388 lines yields K 2 = 238 ± 96 km s −1 , but given the large error, we refrain from providing this value. The disentangled spectrum of the secondary exhibits He i features, but no He ii features. We therefore tentatively classify this system as SB2:, and the secondary as B:, although we cannot exclude that the nebular lines may bias our interpretation. We estimate the light contribution of the secondary at ≈ 5%. With a minimum mass of M 2 > 12.9 +4.7 −3.8 M , the secondary is likely a B0 star, though we cannot fully rule out a BH due to the nebular contamination. The OGLE light curve, phased with the orbital period, is shown in Fig. B.38, showing a very faint periodic signature.
VFTS 481, O8.5 V + O9.7: V: has a reported period of P = 142 d and a high eccentricity of e = 0.93. The presence of a non-degenerate companion is barely seen in the wings of strong He i lines such as He i λ4388 (Fig. B.39). Disentangling supports the presence of a non-degenerate companion. The K 2 value is difficult to constrain given the faintness of the companion. A weighted mean of the strong He i lines yields K 2 = 303 ± 50 km s −1 . The corresponding disentangled spec-  trum approx. matches an O9.7 V star, with an estimated light contribution of 10%. The OGLE light curve does not reveal any significant periods. VFTS 514, O9.7 V, has a reported period of P = 185 d and eccentricity of e = 0.41. No clear minimum for K 2 is obtained, regardless of the disentangled line. The spectrum of the secondary appears to be featureless (Fig. B.40), with exception of the Balmer lines, which are strongly contaminated by nebular lines and are therefore not considered. We can rule out a companion contributing more than ≈ 5% to the visual light, corresponding roughly to B3 V. With a minimum mass of M 2 > 5.3 ± 1.8 M , we cannot rule out that the secondary is a faint non-degenerate star, but given the high BH probability (Table 2), we consider VFTS 514 an O+BH candidate. The OGLE light curve did not reveal any significant frequencies.  Fig. B.43. Comparison of He ii λ4026 spectra of VFTS 603 at RV extremes, along with the disentangled spectra and their sum, as derived for the best-fitting K 1 , K 2 values of 45, 49 km s −1 , respectively. The disentangled spectra are not scaled by the light ratio in this plot.
VFTS 532, O3.5: V:((f*)) + B III, has a reported period of P = 5.8 d and an eccentricity of e = 0.46. Given the very early spectral class of the primary, its orbital solution and RV amplitude of K 1 = 34.4 km s −1 can be well constrained using the He ii lines, which are not impacted by nebular contamination. Similarly to VFTS 386, the He ii lines exhibit clear Doppler motion, while the He i lines appear to be relatively static, which indicates that the secondary contributes mainly to He i lines (Fig. B.41). Because of the low S/N, spectral disentangling of various He i lines does not yield a clear minimum for K 2 , and the results are generally line-dependent, with the K 2 values lying in the range 40 -100 km s −1 . A weighted mean of the measurements of the He i λ4026, He i λ4388, and He i λ4471 lines yields K 2 = 81±32 km s −1 . The disentangled spectrum of the secondary matches best a B-type star, but the exact subtype is difficult to retrieve.
Analysis of the OGLE light curve yields a period of P = 5.796223±0.000002 d, and reveals the unmistakable presence of two eclipses in the system (Fig. B.42) . This further supports the fact that the system hosts two non-degenerate stars. The eclipse ratio implies that the secondary is much cooler than the primary, which agrees with the estimates spectral type. However, a B2 V star would contribute a mere 1% to the flux, and would exhibit RV amplitudes that are 4 − 5 times larger than derived here. Therefore, the secondary is likely a giant star, and we classify it as such. We note that VFTS 532 is one of the four objects in our sample for which X-rays were detected (Table 1).
VFTS 603, O4 III:(fc) +OB:, was reported as a short period binary of P = 1.76 with an eccentricity of e = 0.11. The primary has a very low reported RV amplitude of 11 km s −1 , and the presence of a non-degenerate secondary is not readily evident (Fig. B.43). However, the orbital solution obtained by Almeida et al. (2017) has a fairly high reduced χ 2 value of 6.9. Indeed, 1D disentangling of the strong He ii lines results in a clear residual for the companion that appears distorted. We therefore resorted to 2D disentangling. We focus on the He ii lines, including the He ii λ4026 line, as the remaining lines are too weak and are dominated by nebular emission. However, we do not get consistent results from all lines. The K 1 amplitude is found consistently to be in the range 30 − 50 km s −1 , but K 2 ranges from vanishingly small values for He ii λ4542 to ≈ 50 km s −1 for He i λ4026. The latter is the only line that gives plausible results for the secondary, with K 1 = 45 ± 8 and K 2 = 49 ± 7 km s −1 . These amplitudes are peculiarly low for such a short period, resulting in very small minimum masses of M 1 sin 3 i = 0.08 M and M 2 sin 3 i = 0.07 M , and implying a very small inclination of i ≈ 6 • . While the inclination is not likely (the probability of i ≤ 6.3 • is 0.6%), there is an appreciable probability for this occurring ones in our sample of 51 stars (≈ 25%). Given the short period and low RV amplitudes, it is possible that intrinsic pulsations led to a spurious classification of this object as a binary, and specifically an SB2. Hence, we do not readily adopt our solution in Tables 1 2. Higher-resolution data will be needed to verify the nature of this binary system.
Interestingly, despite of the short period, the OGLE lightcurve does not reveal any clear periodicities beyond integer numbers that are not of astrophysical origin. Moreover, this target is one of the four targets of our sample for which X-rays were detected, albeit at rather low luminosity (Table 1).
VFTS 613, O9 V + O7.5 V, was reported to have a period of P = 69 d and an eccentricity of e = 0.35. However, the RV measurements by Almeida et al. (2017) shows substantial scatter with respect to their derived orbital solution. Inspection of the spectral line variability suggests that the system is in fact SB2, consisting of two similar stars (Fig. B.44). We therefore performed 2D grid disentangling of various lines. The only lines to yield sensible constraints are the He i λ4388 and He ii λ4542 lines, and a weighted mean of the measured RV amplitudes yields K 1 = 66 ± 15km s −1 and K 2 = 96 ± 39km s −1 . The spectral appearance of the secondary suggests that, while less massive, it is hotter than the primary, with an estimated spectral type of O7 and visual light contribution of 49%. This may imply that the primary is an evolved giant or supergiant, though the spectra quality does not enable us to establish this empirically. There exists no OGLE data for this object.
VFTS 619, O8: V, was reported to have P = 14.5 d and e = 0.08. The spectra are strongly contaminated with nebular lines, which makes the disentangling challenging. We cannot find any clear minima in the χ 2 maps of K 2 . The disentangled spectra of the secondary plausible K 2 values yields faint signatures for the secondary, but they are likely the result of the nebular lines (e.g. Fig. B.45). We therefore tentatively classify this system as SB1.
With a minimum mass of 3.8 ± 0.4 M , the companion may be a B-type star, a helium star, or a BH. The OGLE light curve did not reveal any significant frequencies.
VFTS 631, O9.7 V, was reported to have a virtually circular orbit and an orbital period of P = 5.4 d. The spectral variability is not readily suggestive of the presence of a non-degenerate companion. Disentangling of the He i λ4388 line (Fig. B.46) suggests a χ 2 minimum at K 2 = 285 ± 33 km s −1 , but it is barely significant and therefore disregarded. Disentangling of the other He i lines does not result in a clearly defined minimum. Since the nebular contamination is relatively low, our conservative fiducial simulations suggest that we could probe down companions with brightness ratios down to ≈ 3%. Walborn et al. (2014) estimated M V = −3.92 mag for the system (dominated by the primary), such that our threshold would correspond to a companion brightness of M V ≈ −0.1 mag, which matches roughly B6-8 V, or ≈ 3 − 5 M . With a minimum mass of 3.1 ± 0.4 M , the companion could therefore be a late type B or an early type A star, a helium star, or a BH. A Fourier analysis shows a peak at the orbital period, and the phased light curve shows a periodic pattern at a low amplitude when phased at the orbital period (Fig. B.47).
VFTS 645, O9.5 V, has a reported period of 12.5 d and eccentricity of e = 0.23. The disentangling of various strong He i lines such as He i λ4026 (Fig. B.48)  ular contamination. We therefore classify this system as SB1. Walborn et al. (2014) estimated a visual absolute magnitude of M V = −3.86 mag for the system, which is dominated by the primary. Adopting our conservative threshold estimate of 3% for the light contribution for an unseen companion, we can rule out dwarfs with spectral types as late as B6-8 V. With a minimum mass of M 2 2.6 ± 0.5 M , the companion could be an A or Btype dwarf, a helium star, a NS or a BH. The OGLE light curve does not show any significant periods. VFTS 657, O7 II:((f)) + OB:, has a reported period of P = 63 d and eccentricity of e = 0.48. The strong nebular line contamination makes it difficult to determine whether a non-degenerate companion is causing the observed variability (Fig. B.49), and also hampers our ability to disentangle the system robustly. Nonetheless, the He i λ4388 and He i λ4144 lines yield a consistent minimum for K 2 , and the weighted mean of the measurements yields K 2 = 169 ± 38 km s −1 . The spectrum of the companion best matches an early B-type star, but its line coincide with nebular lines, making the SB2 nature of the system questionable. With a minimum mass of 9.2 ± 2.7 M , the companion is either a non-degenerate early-type star or a BH. Analysis of the OGLE light curve did not reveal significant frequencies.
VFTS 702, O8 V(n) + OB + (OB + OB), was reported to be a virtually circular short-period binary with P = 1.98 d. The spectra suffer from low S/N and very strong nebular contamination, hampering our ability to disentangle them robustly. However, Fourier analysis of the OGLE light curve reveals the presence of two eclipsing binaries in the system. The period of the first one closely matches the spectroscopic period derived by Almeida et al. (2017), refined to P = 1.981546 ± 0.000030 d (Fig. B.50). The second period is P = 2.932818 ± 0.000025 d. In     B.50, we also show the light curve phased with this period, after subtracting the 1.98 d signal. Given the short periods, it is clear that the system must be a doubly-eclipsing quadruple system. From the comparable strengths of the eclipses, it is likely that all components are OB-type stars. Examination of available HST images does not reveal any nearby source within 0.1" (roughly 5 000 au), suggesting that the two binaries are gravitationally bound. Unfortunately, the data are not sufficient for disentangling the system, let alone with four components. VFTS 733, O7.5 V + B1 II, was classified as an O9.7p binary with a period of P = 5.9 d and a virtually-circular orbit.  Inspection of the line profiles readily reveals that this system is SB2. In fact, the He ii lines belong almost solely to hotter primary, while the He i and a multitude of metal lines belong to the cooler secondary. Therefore, deriving the nature and orbital solution of this system does not require spectral disentangling, as the RVs can be well measured from isolated lines. Here, we use the He ii λ4542 and the Si iii λ4553 for the hotter primary and cooler secondary, respectively (Fig. B.51). The phased RV measurements and orbital solution are shown in Fig. B.52. We find K 1 = 39.4 ± 1.9 km s −1 and K 2 = 102.2 ± 0.9 km s −1 for the primary and secondary, respectively. The less massive secondary appears to contribute roughly half of the flux in the visual, suggesting that it is evolved (luminosity class III-I). This is supported by the spectral appearance of the secondary. The low minimum masses of M 1 sin 3 = 1.26 M and M 2 sin 3 i = 0.49 M imply a low inclination of i ≈ 20 • . Eclipses are therefore not to be expected. The OGLE light curve shows a periodic pattern when phased at the orbital period (Fig. B.53).
VFTS 736, O9.5 V + B:, has a reported period of P = 68.8 d and eccentricity of e = 0.09. Disentangling of various He i lines does not yield a consistent solution for K 2 , and the spectral variability is not readily suggestive of a non-degenerate companion (Fig. B.54) Fig. B.56. Comparison of He i λ4388 spectra of VFTS 750 at RV extremes, along with the disentangled spectra and their sum, as derived for K 1 = 29 km s −1 and K 2 = 3 × K 1 = 87 km s −1 . The disentangled spectra are not scaled by the light ratio in this plot.
lines, but these may be the result of nebular contamination. As M 2 > 3.9 ± 0.5 M , the companion could either be an B-or Atype star, a helium star, or a black hole. Analysis of the OGLE light curve did not reveal significant frequencies.
VFTS 743, O9.5 V((n)), has a reported period of P = 14.9 d and a virtually circular orbit. The spectral variability and disentangled spectrum of the secondary for various K 2 values do not support the presence of a non-degenerate companion in the system (Fig. B.55). Some faint He i features, as well as Balmer features, are seen in the disentangled spectrum of the secondary, but their appearance suggests that they originate primarily in nebular contamination, and we therefore classify the system as SB1. Given the minimum mass of M 2 > 2.1 ± 0.3 M , the nature of the secondary is virtually unconstrained. Analysis of the OGLE light curve did not reveal significant frequencies.  VFTS 750, O9.5 V + B:, is reported to be an eccentric long period binary with P = 417 d and e = 0.78 . The spectral variability is not readily suggestive of the presence of a nondegenerate companion (Fig. B.56). When disentangling for plausible K 2 values in the range 2-5 K 1 , faint features of He i can be seen. Since the nebular contamination is modest in this star, it is likely that these features are real. However, due to the lowamplitude motion of the primary, the relatively poor Doppler coverage, and the faintness of the potential companion, we cannot retrieve an unambiguous value for K 2 . We therefore tentatively classify the companion as B: and the system as SB2:. With a minimum mass of M 2 > 5.6 ± 1.2 M , the companion could be a non-degenerate object or a BH. Analysis of the OGLE light curve did not reveal significant frequencies.
VFTS 769, O9.7 V, has a short period of P = 2.4 d and a virtually circular orbit. The spectral variability of strong He i lines such as He i λ4388 (Fig. B.57) is not readily suggestive of the presence of a non-degenerate companion. The disentangled spectra of the secondary for various K 2 values are suggestive of Balmer absorption and hints of He i lines, but we cannot exclude that they are the result of nebular contamination. With a minimum mass of M 2 > 1.9 ± 0.2 M , the nature of the secondary is unconstrained. Analysis of the OGLE light curve reveals a significant peak at the orbital period of P = 2.365644 d and its harmonics. The phased light curve is shown in Fig. B.58.
VFTS 779, B1 III, has a reported period of P = 59.9 d and a virtually circular orbit. The spectral variability of the He i does not readily suggest the presence of a non-degenerate secondary (e.g. Fig. B.59) . Disentangling of various He i lines results in nearly-flat χ 2 (K 2 ) maps, though an insignificant tendency to-  wards values in the range K 2 ≈ 2 − 4 × K 2 is noted. The disentangled spectrum of the secondary for such K 2 values shows very faint He i absorption in some of the lines, though it is not clear whether this is due to low levels of nebular contamination. Given the relatively high quality of the data in this case, we can exclude companions contributing more than 2 − 3% to the total light, which would correspond roughly to a B3. We classify this system as SB1. As M 2 > 3.9 ± 0.6 M , the companion is either a mid-to late-type B dwarf, a He star, or a BH. Given the BH probability (Table 2), we consider this system a good B+BH candidate. Analysis of the OGLE light curve did not reveal significant frequencies.
VFTS 802, O7 V: + O8 Vn, has a reported period of P = 182 d and an eccentricity of e = 0.60. The presence of a non-degenerate companion is readily seen in a multitude of He lines as a shallow, broad feature (Fig. B.60). Given the strong contribution of the secondary, we perform 2D Disentangling of the strong He i and He ii lines. A mean average yields: K 1 = 70.4 ± 1.9 km s −1 , K 2 = 76 ± 14 km s −1 . The RV amplitude of the secondary is more difficult to constrain due to the substantial rotational broadening of its spectrum. The companion is estimated to contribute 36% of the visual flux, and is classified as O8. Analysis of the OGLE light curve did not reveal significant frequencies.
VFTS 810, O9.7 V + B1 V, has a reported period of P = 15.7 d with an eccentricity of e = 0.68. The presence of a nondegenerate companion was already recognised by (Almeida et al. 2017), who however could not measure its RVs. The lines of the secondary become deblended in certain epochs, making its presence readily evident (Fig. B.61). 2D Disentangling of the strong He i lines yield a weighted mean of K 1 = 108 ± 9 km s −1 and K 2 = 162 ± 27 km s −1 . The companion is estimated to contribute Wavelength [Å] Comparison of He i λ4471 spectra of VFTS 827 at RV extremes, along with the disentangled spectra and their sum, as derived for K 1 = 25 km s −1 and K 2 = 2 × K 1 = 50 km s −1 . The disentangled spectra are not scaled by the light ratio in this plot.
29% of the visual flux, and is classified as B1. Analysis of the OGLE light curve did not reveal significant frequencies.
VFTS 812, O4 V((fc)), has a reported period of P = 17.3 d and eccentricity of e = 0.62. No sign for a companion is seen in He ii lines (e.g. Fig. B.62), and He i lines are contaminated by nebular lines, hampering our ability to use them robustly. No clear value for K 2 is retrieved from disentangling. The disentangled spectra of the secondary obtained for plausible K 2 values show Balmer lines and some hints of He i lines, but not sufficiently clear to be able to classify this system SB2 unambiguously. With M 2 ≥ 5.1 ± 0.7 M , the companion may be a latetype OB star or a BH. Analysis of the OGLE light curve did not reveal significant frequencies.
VFTS 827, B1.5 III, has a reported period and eccentricity of P = 43 d and e = 0.24. The spectral line variability is not suggestive of the presence of a second non-degenerate companion (e.g. Fig. B.63) Fig. B.64. Comparison of He i λ4471 spectra of VFTS 829 at RV extremes, along with the disentangled spectra and their sum, as derived for K 1 = 13 km s −1 and K 2 = 4 × K 1 = 52 km s −1 . The disentangled spectra are not scaled by the light ratio in this plot.

4390 4395
Wavelength favours K 2 values in the range K 2 70 km s −1 . Disentangling for various K 2 values does not result a clear stellar spectrum for the secondary. However, strong Balmer lines are seen in the disentangled spectrum of the secondary. While the nebular contamination in this case is very low, the narrow widths of the Balmer lines in the disentangled spectrum of the secondary raise the suspicion that they originate in nebular contamination. We therefore classify the system as SB1:. With a minimum mass of M 2 > 2.5 ± 0.5 M , the companion might also a late B-type star or an A-type star, and we cannot fully rule out that it is a NS or a BH. No OGLE light curve is available for this object.
VFTS 829, B1.5 III, has a reported period and eccentricity of P = 203 d and e = 0.27. The spectral line variability is not suggestive of the presence of a second non-degenerate companion (e.g. Fig. B.64) , and disentangling of various He i results in relatively flat χ 2 (K 2 ) maps. The disentangled spectrum of the secondary for K 2 values in the range 2 − 10 × K 1 is fully free of features with the exception of weak Balmer absorption that could be the result of nebular-line contamination. Given the low levels of nebular-line contamination and the relatively high S/N ratios of the spectra, we classify this system tentatively as SB1:. However, the low minimum mass of M 2 > 2.0 ± 0.5 M leaves the nature of the secondary virtually unconstrained.
VFTS 887, O9.7: V: + O9.5: V, has a reported period of P = 2.7 d and a virtually circular orbit. The spectral line variability (Fig. B.65) clearly reveals the presence of two components in the spectrum. We therefore utilise 2D disentangling on various He i and He ii lines, which yield consistent results. A weighted mean yields K 1 = 105 ± 7 km s −1 and K 2 = 98 ± 12 km s −1 , implying that the object identified as the primary by Almeida et al. (2017) is the slightly less massive component. The disentangled spectra correspond to O9.7 V and O9.5 V spectral types, respectively,  and are estimated to contribute 45% to the visual flux. Despite the short period, the OGLE light curve does not show a clear period pattern (Fig. B.66).