The YMDB catalog: Young massive detached binaries for the determination of high-precision absolute stellar parameters
Pablo Martín-Ravelo Roberto Gamen Julia I. Arias André-Nicolas Chené Rodolfo H. Barbá






















The YMDB catalog







Departamento de Astronomía, Universidad de La Serena, Av. Cisternas 1200 Norte, La Serena, Chile.

    Gemini Observatory/NSF’s NOIRLab, 670 N. A‘ohoku Place, Hilo, HI 96720, USA 

    pablo.martin@noirlab.edu
    Instituto de Astrofísica de La Plata, CONICET–UNLP, Paseo del Bosque s/n, 1900, La Plata, Argentina.

    Facultad de Ciencias Astronómicas y Geofísicas, Universidad Nacional de La Plata, Argentina.
    





 


    Massive stars play a crucial role in the cosmic dynamics and chemical evolution of galaxies. Despite their significance, our understanding of their evolution and properties remains limited. An accurate determination of stellar parameters, such as the mass and radius, is essential for advancing our knowledge. Detached eclipsing binaries (DEBs) are particularly valuable for these determinations due to the minimal interaction between their stellar components, allowing for precise measurements.
    
    This study aims to introduce the Young Massive Detached Binary (YMDB) catalog, designed to address the gap in the high-precision absolute parameter determination for young massive stars. By focusing on DEBs within the spectral range O9-B1, this catalog seeks to provide a reliable database for future astronomical studies and improve our understanding of massive star evolution.
    
    We conducted a photometric analysis of 87 young massive stars in detached eclipsing systems using TESS light curves (LCs) that were processed through a custom pipeline. This analysis involved determining the amplitude of magnitude variations, orbital periods, times of minima, eccentricities, and the presence of apsidal motion and heartbeat phenomena. A thorough literature review was performed to obtain MK spectral classifications. We performed our own spectral classification of 19 systems to support the sample where a new classification was lacking or inconclusive.
    
    The analysis identified 20 previously unreported binary systems, with 13 newly recognized as variable stars. Among the 87 stars examined, 30 are confirmed as YMDB members, and 25 are candidates pending spectral classification. The exclusion of the remaining 32 stars is attributed to unsuitable spectral types or their nondetached binary nature. Notable findings include the identification of new LC classifications, eccentricities in 13 systems, and heartbeat phenomena in several targets.
    
    The YMDB catalog offers a resource of high-quality LCs and reliable stellar classifications, serving as a valuable tool for the astronomical community.




    

    The YMDB catalog: Young massive detached binaries for the determination of high-precision absolute stellar parametersFull version of Tables 2, 3 and 4 are only available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/.
    
    
    Pablo Martín-Ravelo 
    1,2
    Roberto Gamen 
    3,4
    Julia I. Arias 
    1
    André-Nicolas Chené 
    2
    Rodolfo H. Barbá In Memoriam (1962–2021)
    1
    
    Received: 20 June 2024 / Accepted: 05 July 2024
================================================================================================================================================================================================================================================================================================================================








§ INTRODUCTION


Massive stars, defined by their cataclysmic end as core collapse supernovas, are typically those with an initial mass of eight or more solar masses and they fall into the OB spectral classification. These stellar behemoths play a crucial role in cosmic dynamics and chemical evolution, with their supernova events significantly enriching the interstellar medium with heavy elements. Understanding massive stars is fundamental to comprehending stellar evolution, the evolutionary history of galaxies, and the Universe at large. Due to their immense luminosity, they are key observable objects in distant galaxies, making them essential for astronomical studies with both current and forthcoming space and ground-based large telescopes.



Our current understanding of massive star evolution remains limited, as evidenced by the persistent mass discrepancy between empirical measurements from orbital dynamics and theoretical model predictions. This discrepancy has been a longstanding issue in astrophysics <cit.>. Factors inherent in binary systems, such as stellar rotation, which influence orbital ellipticity, synchronization, tidal effects, limb darkening, and radiative propagation, are crucial in this context. Evolutionary models that incorporate rotation show marked differences compared to nonrotational models. Additionally, research by <cit.> highlights significant variations in how current models handle mass loss in these stars.



The determination of reliable masses for massive stars is only achievable in binary systems, particularly through the study of light curves (LCs) and radial velocity (RV) curves of eclipsing binaries. The likelihood of encountering massive stars within binary systems is fairly high, especially during their main sequence phase, the most extended period of a star's lifecycle. Comprehensive spectroscopic surveys reveal that around 75% of main sequence O-type stars <cit.>, and a similar ratio of early B-type stars <cit.>, are part of gravitationally bound systems.



The endemic multiplicity among massive stars offers unique research opportunities, although the method has its limitations. An inherent consequence of binarity is the interaction among its components, which make their evolution paths deviate from those of solitary stars of similar types. This interaction issue was emphasized by <cit.>, who found that at least 71% of O-type stars in binary systems interact with their companions during their lifetimes, with 29% eventually merging into a single entity. Detecting past interactions poses challenges, as evidenced by discrepancies between empirical mass measurements and predictions from atmospheric or stellar evolution models, which often overlook these interactions. 



Detached eclipsing binaries (DEBs) stand as critical objects in this context. These systems, charaterized by a negligible or minimal interaction between their stellar components, perfectly meet the criteria for accurate stellar analysis. Given that most massive stars are part of multiple systems, it is preferable to focus on those with well-documented multiplicity. The best single stars for study are often the components of binary systems, as proposed by <cit.> and further discussed by <cit.>. DEBs provide a unique opportunity to determine the physical properties of stars with a high accuracy, such as their masses, and radii. Moreover, they can be used to determine distances. In this context, DEBs are key objects. Determinations from DEBs form the basis for calibrating crucial relationships among stellar parameters, such as the mass–luminosity relationship.



Among massive stars, those within the O8–B3 spectral-type range are predominant (just as a consequence of the stellar mass distribution). However, only a limited number have had their absolute parameters accurately determined. 
Our research aims to address this gap by examining a selection of DEBs within this spectral range. We intend to analyze their LCs and RV curves to derive precise stellar parameters. 



In this work, we present a curated database of young massive stars, with spectral types around B0 V, in DEBs. This endeavor is the first step to maintain an up-to-date database of empirical absolute parameters for massive stars. To achieve our objectives, we have carefully selected systems from a careful search in the available literature (see details in Sect. <ref>). After that, we generated the LCs of each target in the database using data from the NASA Transiting Exoplanet Survey Satellite (TESS), and analyzed their variations to identify DEBs. This is explained in Sect. <ref>. TESS LCs were also examined for periodicities, to determine eccentricities, ephemerides, and also for the detection of pulsation and/or heartbeat phenomena. Some targets were spectroscopically observed to obtain their spectral classification, which were confusing in the literature. Observations are described in Sect. <ref>. All of the results that populate the Young Massive Detached Binaries (YMDB) catalog are shown in Sect. <ref>.







§ GENERAL METHODOLOGY FOR THE YMDB CATALOG


Our study's approach integrates a thorough search and analysis methodology to identify and evaluate candidates for detailed investigation. Starting with an extensive review of databases and literature, we refined our list of potential candidates by examining their spectral classifications and the likelihood of misclassification, and by quickly reviewing raw TESS LCs to exclude any undoubted nondetached binaries. This process narrowed our initial pool to 87 candidates.

In preparing the LCs, we initially set up background and target masks using a structured pixel grid to minimize light contamination, which is vital for precise LC extraction. We continuously evaluated data quality and implemented polynomial fitting to correct known and recurrent flux pattern variations in the TESS data, such as rollovers, ensuring a cleaner and more accurate representation of the LCs. 

We utilized Gaussian Fitting (GF) and Box Least Square (BLS) methods to accurately determine orbital periods and other temporal features. Further analysis of the LCs involved meticulous visual inspection to confirm the detached nature of systems and identify any additional features such as pulsations or heartbeats.

To enhance the validity of the spectral classifications, new spectroscopic observations were conducted using the Jorge Sahade telescope at CASLEO. These efforts were focused on acquiring low-resolution spectra, which were processed using established methods, to clarify any existing classification uncertainties.




 §.§ Candidates selection
 

Candidates for the study were identified through extensive searches in various databases and literature. These sources included the Spectroscopic Binary Orbits Ninth Catalog <cit.>, OWN Survey <cit.>, IACOB <cit.>, Eclipsing Variables Catalog <cit.>, and the multiplicity of northern O-type spectroscopic systems project <cit.>, among others. Additionally, potential targets suggested by collaborative efforts within the astronomical community were also considered.

The search criteria encompassed a wide range, including systems with components in the O7–B3 III–V spectral class, to account for potential misclassifications. This cautious approach was used to ensure all potentially relevant systems were considered, particularly those that might align with the desired O9-B1 IV-V classification upon closer scrutiny. 
A thorough review of the bibliography was conducted to assess the reliability of the spectral classifications, the availability of LCs, and the reported detachment status in the literature, narrowing down the original list of 339 targets to 186 candidates.

The process continued with an analysis of TESS data to eliminate any candidates showing clear non-Algol type variability (i.e, Beta Lyræ  and W Ursæ  Majoris) reducing the number to 114 candidates. A custom pipeline was then developed for extracting high quality TESS LCs, allowing for the estimation of orbital periods (P) and minima (T0), as well as the identification of eccentric orbits, heartbeats and apsidal motion throughout the TESS sectors. This was achieved using both GF and BLS methods. The candidates were meticulously scrutinized, focusing on the identification of Algol-type variables. Systems not fitting the Algol type criteria were excluded, resulting in 87 stars qualifying as potential candidates which we present in this work.
Details of this task are given in the following.





 §.§ Construction of light curves
 

In constructing the LC, the initial step involved assessing the quality masks for each TESS sector, related to the star, to determine whether the default rejection criteria were adequate. For the majority of LCs, it was found that the default quality mask met the needs effectively. 




  §.§.§ Background mask


To construct a sky background mask for each sector, a 30×30 pixel area centered on the target system is analyzed. Stars within this area, along with their magnitudes, are identified using Gaia DR3 data. The lightkurve tool <cit.> is then employed to generate a mask with threshold zero, specifically excluding pixels affected by stars brighter than a determined magnitude limit. This upper limit is set at 5 magnitudes fainter than the brightest star within a 3-pixel radius of the target system, effectively minimizing contamination from other stars in the extraction of the desired LC.






  §.§.§ Target mask


Following the construction of the sky background mask, a target mask is developed using a similar methodology. The function that delineates boxes around mapped stars based on their brightness is also applied to isolate the Target Pixel File (TPF) surrounding the star of interest. Adjustments to the size and position of this cut, alongside a specified threshold value, are made using the lightkurve tool to craft the target mask. This stage might require iteration, involving comparisons of the resulting LC with published ones and the LCs of individual pixels, to verify the accuracy of the selected parameters and the suitability of the target mask.





  §.§.§ Data selection


The analysis of the unfolded extracted LC involves identifying and flagging problematic data that could compromise the construction of the LC. A visual inspection of the cadence for each TESS sector is essential, as it can reveal areas potentially affected by gaps in data—either due to TESS CCD readouts or data omitted by the quality mask. Data adjacent to these gaps might exhibit a different background profile, making parts of the LC unreliable if the background extraction fails to account for this variance. Additionally, these segments may experience slight magnitude shifts (Δmag) and may require independent normalization, posing challenges for systems with long periods.





  §.§.§ Determination of ephemeris
 

In determining the ephemeris of our studied systems, we incorporate two principal methods: GF and BLS. GF is primarily utilized for its simplicity and effectiveness in identifying eclipse minima by fitting Gaussian profiles, a method similar to the Phase Dispersion Minimization (PDM) approach.

Phase dispersion minimization is traditionally used to detect stable periodic signals by minimizing scatter across phased data, making it ideal for datasets like those from TESS which may include gaps or nonsinusoidal variations. By adapting this methodology, our GF process not only identifies the minima but also estimates the period from these minima's separations, similarly to how PDM assesses periodicity by examining the variance across different data bins.

After initial minima identification, we fold the LC and apply GF iteratively across all sectors, optimizing the ephemeris precision in a manner analogous to refining PDM's phase coverage by adjusting bin overlaps or employing smoother functions in updated PDM versions.

Box least square is used alongside GF for period verification, with both methods critically reviewed against known periods from literature when available in order to cross-verify with established data the reliability of our adapted PDM-inspired techniques in analyzing TESS's time series data.





  §.§.§ Phase-synchronized polynomial fitting


Photometric observations from the TESS mission, similar to those from its predecessor Kepler, are susceptible to instrumental systematic trends that can obscure or distort the intrinsic stellar signals. These systematics arise primarily from spacecraft jitter and other operational imperfections, which introduce noise and trends across different timescales into the captured LCs. 

To mitigate these effects, we employed a polynomial fitting method, analogous to the Pixel Level Decorrelation (PLD) method utilized in the EVEREST pipeline for K2 data. This method has been demonstrated to effectively remove correlated noise due to spacecraft motion by fitting and subtracting systematics directly from the pixel-level data <cit.>. 

While our approach uses a simpler polynomial fitting of lower orders, it shares the core principle of modeling the observed data to isolate and remove instrumental signatures. The primary distinction of our method, which we designate as "Phase-Synchronized Polynomial Fitting" (PSPF), lies in its utilization of the previously known orbital periods to synchronize (fold) phase points within a TESS sector. This synchronization enables the calculation of fits for data points that share the same phase, thus forming a more complex function that better adapts to the unique signal trends of TESS data. Therefore, PSPF is suited for systems with orbital periods shorter than the duration of a TESS sector. The method cannot be applied if the period exceeds the sector length (PLD should be used in those cases), unless the signal is expected to remain constant in different phases of the LC, such as outside the eclipses in highly detached binaries. It becomes particularly reliable when the orbital period is substantially shorter than the sector length, allowing for enhanced signal integrity post-correction through multiple phase folds. 

The correction process begins with the identification of phase points across the LC that are consistent over multiple orbits, meaning without flux variations other than those related to the eclipsing nature of the system or other synchronized pulsations (whose periods are multiple of the orbital period). This leverages the comparative stability of shorter-period systems which may ignore not-synchronized variations if enough points per phase are provided. Each phase point will provide then a value not only for that phase, but for every other phase point in the LC. 
A weighted mean is calculated for these values in each phase point to construct a robust fit model. This model is then used to adjust the LC, with the median of all polynomial corrections applied to derive the final corrected LC. This procedure was only implemented when a consistent pattern of systematic error was evident across the entire sector, ensuring that the corrections made were both meaningful and substantiated by the data.

For about 20% of the candidates in our study that exhibited clear systematic trends, the PSPF was crucial for assembling a clean LC. A polynomial fit of order 1 was mainly used, but the selection of polynomial order was tailored to each target's specific noise characteristics, with higher orders only used seldom and for more complex systematic patterns.

Our PSPF corrections were validated against already corrected LCs, showing a significant reduction in noise and the preservation of intrinsic astrophysical signals and effectively detrending the TESS LCs, akin to the results reported by the EVEREST pipeline when applying PLD to K2 data <cit.>.






 §.§ Light-curves analysis


All generated LCs underwent meticulous inspection to confirm their detached nature, identify any pulsations or heartbeats, and determine eccentricity and apsidal motion. Visual inspection was used to confirm the detached nature of the system, which was unequivocally determined by the distinct transitions marking the beginning and end of eclipses in the LC. Pulsational features, indicative of either intrinsic variability of the components or the presence of a complex multiple system, were also identified through visual inspection of the LCs. In cases of multiple systems, efforts were made to disentangle the LCs using templates crafted from median data selected, when possible, from segments of the LC minimally affected by eclipses or intrinsic pulsations of the main system. The heartbeat phenomenon was identified as a distinct photometric variation occurring during orbital phases with closer eclipses, indicative of the expected periastron. We compared those variations against all forms of heartbeat showed in <cit.>. Detailed descriptions of these cases are provided in Section <ref> and their LCs are shown in Figure <ref>.

For systems where eccentricity was not immediately obvious, we employed GF to precisely measure the centers of the primary and secondary eclipses, using a threshold of Δϕ=0.0002, to determine whether their phase separation differed from 0.5. It is important to note that an exact half-period interval between eclipses does not necessarily imply a circular orbit; the system could be eccentric with a periastron longitude of either 90^∘ or 270^∘. Therefore, values listed in Tables 2, 3, and 4 may imply eccentricity but never circularity. Variations in the fit to the secondary eclipse across different sectors were analyzed, and systems exhibiting changes beyond a certain threshold (Δϕ=0.0002) were flagged for potential apsidal motion. Additionally, since heartbeat phenomena occur exclusively in eccentric binaries, the detection of heartbeat signals in cases where only one eclipse was observed allowed us to identify the system as eccentric. 





 §.§ Spectroscopic observations and classification
 

To address the inconsistencies in spectral classifications found in the existing literature, we undertook the task of obtaining new spectra.
Low-resolution spectra were obtained using the Jorge Sahade telescope at Complejo Astronómico El Leoncito (CASLEO) in San Juan, Argentina, utilizing the REOSC spectrograph in single mode alongside the new Sophia CCD. This equipment yielded spectra with a resolution of R∼1000, covering the wavelength range of 3760–5845 Å  suitable for spectral classification in the MKK system.
Spectral data processing was conducted in the standard way using iraf[NOIRLab IRAF is distributed by the Community Science and Data Center at NSF NOIRLab, which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the U.S. National Science Foundation.] routines.

The resulting spectra are depicted in Fig. <ref>. They have been classified following the criteria described in <cit.>. To summarize, we employed the ratio between He ii λ4541 and He i λ4387, for late O-types, and Si iii λ4552/Si iv λ4089, for early B-types. Further details about spectral classification of OB stars can be found in the references provided.

The new spectral types are presented in the corresponding Tables <ref>, <ref> and <ref>, column ST1, and are indicated with the label “tw”. No ST2 is given due to the low spectral resolution, which prevents the disentangling of both components in the system.









§ RESULTS
 

Among the 87 systems, this study reveals 20 new eclipsing binaries (Table <ref>), including 13 whose variable nature was previously unknown, five targets which were known as photometric variables or spectroscopic binaries, but not as DEB, and two of them were identified as nonthermal contact (because they present different light minima). Moreover, we introduce new LC classifications for 30 systems and report novel findings such as eccentricity in 13 systems, heartbeat features in 17, new types of variability in 11, and debut TESS LC presentations for the majority of the sample. 

The study's results are presented in three distinct tables: confirmed members of the catalog (Table <ref>), candidates pending spectral classification (Table <ref>), and unqualified candidates (Table <ref>). Confirmed members are systems that meet the essential conditions of displaying Algol-type LCs (detached) and having at least one component with a spectral classification within the O9-B1 V range. Candidates are also detached systems that show potential signs of meeting the spectral type criteria but lack definitive confirmation due to incomplete or inconclusive spectral classification. Unqualified systems are those whose LCs do not appear detached (EB or EW), have primary spectral classifications clearly out of range, or both. The secondary component's classification might either be out of range or impossible to detect due to limitations in current spectroscopic data and analysis. While unqualified detached systems with indeterminate secondary classifications are currently unsuitable for the catalog, future higher-resolution or more sensitive spectra could reveal additional stellar features, prompting their reconsideration as candidates.

For each group, we detail findings from our photometric analysis, including delta magnitude (Δ mag), period (P), time of minima (T_0), observed apsidal motion (Apsidal), additional system variability (MultiP), heartbeat phenomena (HB), and eccentricity (e). Additionally, we provide our spectral classifications for certain stars, either to fill gaps where classifications were absent or to verify existing classifications.

Regarding new spectral types, we confirmed HD 298448, CD-28 5257, HD 309036, V* IK Vel as entries to the YMDB (Table <ref>),
and emphasize that certain eclipsing binary systems were excluded from the YMDB catalog because their spectral types fall outside the range studied here. 
An intriguing example is CM CMa (=Gaia DR3 2928505622380096256), which lacks a documented spectral type in the eight publications indexed in the Simbad database. Our April 2023 night spectrum reveals absorptions consistent with an F5 V type, characterized by similar H and K Calcium lines and the presence of the G-band.
Another case is GN Nor (=Gaia DR3 5884730512723833984), which similarly lacks spectroscopical references in Simbad. Our CASLEO spectrum is classified as A0 V, primarily due to the dominance of Balmer lines.
Similar cases are V* GN Nor and * f Vel, whose spectral classifications resulted out of the range for this catalog.
These spectra are shown in Fig. <ref>.

In certain LCs, we observed additional variability beyond the typical eclipsing patterns. These variations, unless identified as heartbeat phenomena, are denoted with 1 in the tables. Therefore, we detected such oscillations in 37 systems, 11 of which are unprecedented.
Most newly discovered pulsating systems have dedicated paragraphs, except 
HD 309036, 
HD 204827, 
all of them well known DEBs. 

Finally, the visual inspection of TESS LCs allowed us to detect apsidal motion. This phenomena is indicated in the tables (with 1) when the times of the secondary eclipses seem to vary among different sectors.




 §.§ New detached eclipsing binaries discovered
 





We offer a detailed insight into select DEB featured within the YMDB.


BD+66 1674: 
Initially noted as a probable RV variable <cit.>, no additional references were found in the available literature. Although <cit.> classified it as O9.7 IV, an analysis, incorporating high-resolution spectra from the Galactic O-Star Spectroscopic Survey <cit.> and the Library of Libraries of Massive-Star High-Resolution Spectra <cit.>, revealed it to be a B0 V + B0 V system (Maíz Apellániz, priv. comm.). Consequently, it has been included in the YMDB.  Its LC displays additional variations beyond the 18-day eclipsing pattern, resembling those of highly eccentric detached double-eclipsing systems. Moreover, a notable bump between eclipses suggests the presence of a heartbeat-like feature (see Fig. <ref>). Considering the pixel size of TESS data, we conducted a search for other targets within a 40 arcsec radius and identified seven Gaia sources. These sources, at least four G-magnitudes fainter than BD+66 1674, complicate the determination of whether the short-period binary is one of these sources or an unresolved component of BD+66 1674.  





BD +66 1675: 
Classified as O7.5 Vz <cit.>, high-resolution spectra obtained within the context of LiLiMaRlin <cit.> revealed BD +66 1675 to be a triple system comprised of O7.5Vz+O8V +B components. A preliminary analysis of 20 spectra indicated that the RV of the spectral lines belonging to the O7.5 Vz star do not vary within errors. However, the O8V + B pair exhibits RV motion in accordance with the photometric period. Moreover, preliminary orbital elements suggest an eccentric orbit, explaining why only one eclipse is observed. This eclipse occurs when the massive star passes in front of the system, making it a secondary eclipse. Also, we detect a heartbeat-like behavior immediately after the eclipses, but only in the Sector 58 (Fig.  <ref>), and short-period variability identified as pulsations.





HD 278236: 
Classified as O9 V by <cit.>, no publications reporting variability or binarity were found. It is noteworthy that the relative position of the secondary eclipse appears to vary in different sectors, indicating apsidal motion (Fig. <ref>). Moreover, apsidal motion is making the secondary eclipse narrow overtime (from sector 17-59 ∼1100 days and from 59-73 equivalent to ∼1400 days, for a total time span of ∼2500 days).





RAFGL 5223: 
Recognized as a Herbig Ae/Be star candidate, its entire bibliography is centered around this characteristic. However, no indications of variability, either spectroscopic or photometric, were identified. To scrutinize the indicated spectral classification in the Simbad database, we downloaded a X-Shooter spectrum (program ID 084.C-0952(A)) from the ESO portal. By evaluating the ratios between He ii and He i pairs, we determined it to be of O9.2 type, with He ii λ4686/He i λ4713 and Si iv λ4089/He i λ4026 ratios indicative of a V luminosity class. Consequently, it has been included in YMDB as an eccentric short-period binary. We note a potential heartbeat feature, manifested as a bump in the orbital phases where the periastron passage is expected.
We also remark a probable apsidal motion, as separation between eclipses seems to vary among Sectors 7 and 33 (∼700 d).



TYC 8174-540-1: 
Initially identified as an OB star by <cit.> and subsequently classified as O9.5 Vn by <cit.>, TYC 8174-540-1's spectral type was verified through our own CASLEO spectrum (see Fig. <ref>). Despite an absence of reports regarding its binary nature, the TESS LC unmistakably depicts both eclipses, which, in fact, are transit and occultation events. This characteristic increases its significance as a benchmark for determining stellar parameters. Notably, the system exhibits eccentricity and appears to manifest a heartbeat phenomenon before entering the primary eclipse.



HD 102475: 
We found no evidence in the existing literature suggesting that it is a binary or variable star. It is only mentioned in general works on spectral types <cit.>, where it is classified as B0.5 II and B1 III, respectively. Therefore, we include this target as a candidate until its spectral type is confirmed with modern spectroscopic analysis. It also exhibits short-period variability, interpreted as pulsations.



HD 114026: 
While we found no evidence suggesting variability, <cit.> identified it as a probable SB2 system. Despite its spectral classification exhibiting significant dispersion, ranging from OB to B2, with <cit.> designating it as a B0.5 V:n type, we prefer to be cautious and include HD 114026 as a candidate for the YMDB.



CPD -63 3284: 
This star lacks a reliable spectral classification, and unfortunately, we were unable to observe it during our run at CASLEO. Initially identified as an OB star by <cit.>, no additional classification information was found. Nevertheless, its LC unmistakably displays double eclipses, which indeed correspond to occultation and transit events. It also exhibits short-period variability, interpreted as pulsations. We include it as a candidate in the YMDB.



UCAC2 5911156: 
Initially recognized as an early-type star by <cit.> and subsequently classified as B0.5 V <cit.>, we acquired a new spectrum leading to a reclassification of its luminosity class to III. This adjustment is based on the almost similar intensities of the absorption lines He i λ4387 and Si iii λ4552. The giant nature of UCAC2 5911156 naturally explains the pulsations noted in the LC shown in Fig.<ref>. The LC exhibits eclipses with a superimposed set of oscillations forming a beating pattern, akin to those initially discovered in HD187091 <cit.>. Notably, these damped oscillations do not align with the orbital period, as the maximum amplitudes do not occur at the same phase. To the best of our knowledge, this is the first identification of the system as a DEB and pulsating. However, a detailed analysis of this intriguing system is beyond the scope of this work. 





HD 142152: 
Initially classified as a B0 III massive star by <cit.> and subsequently confirmed by <cit.>, HD 142152 lacked published spectra. To address this absence, we observed this target using CASLEO and confirmed its spectral classification (see Fig. <ref>). The identification as a probable binary arises from two discrepant radial velocities reported by <cit.>. The TESS LC distinctly illustrates a highly eccentric detached double-eclipsing behavior. Additionally, pulsation-like variations, occurring with 1/8 of the orbital periodicity, and a discernible heartbeat feature are detected. This renders the target exceptionally interesting, despite its exclusion from YMDB due to its luminosity class.



CD -54 6456: 
Prior to this study, there were no records indicating variability or binarity for this system. Classified as O9.5 V <cit.>, our analysis involved six high-resolution spectra, which revealed no RV variations in the features of the O-type star. Consequently, this suggests that the eclipsing binary likely involves another stellar pair within the TESS pixel or is indistinguishable from CD -54 6456. For this reason, we have included it as a candidate in the YMDB.



HD 144918: 
This target exhibits one of the most widely scattered spectral type assignments in the literature, ranging from O5/7 <cit.> to B0 <cit.>, with none originating from a modern study. Consequently, we decided to acquire a new spectrum. The CASLEO spectrum is classified as O8 V, with this classification based on the observation that He ii λ4542 is slightly fainter than He i λ4471. Concerning its binary nature, it is reported as SB2 by <cit.>, yet no period is provided, and no dedicated work was found in the literature. 



BD +55 2722: 
This system is, indeed, a Trapezium-like configuration with components named A (O8 Vz), B (O9.5 V), and C (O7 V(n)z+B), according to the GOSSS <cit.>. Given that the three sources are indistinguishable in TESS data, we attribute the DEB discovery to the C component, thereby designating BD +55 2722 C as a candidate in the YMDB. This classification stands until the spectral type of the secondary component is more precisely determined.  



HD 277878: 
Despite extensive literature search, no evidence of photometric variability was found for this target. Originally identified as an OB or B0-type star, it was recently reclassified as O7 V((f))z based on LAMOST spectra <cit.>, and indicated as SB1. Consequently, we have excluded it from consideration in the YMDB.



* psi02 Ori: 
It is a widely recognized spectroscopic binary <cit.>, and also known to exhibit ellipsoidal variations <cit.>. In the TESS data, its double-eclipsing nature is clearly evident, alongside the ellipsoidal variations. However, due to its spectral type being identified as B1 III + B2 V <cit.>, we have excluded it from consideration in the YMDB.



LS  VI +00   25: 
It was initially identified as an SB1 system by <cit.>, although no photometric variations have been reported to our knowledge based on available literature. With a spectral type of O9.5 V <cit.>, we have included it in the YMDB.



CD-53 6352: 
This star is recognized as double or multiple according to the Washington Double Star Catalog (WDS J16001-5355AB), but no photometric or RV variability has been identified. Analysis of TESS data has not conclusively determined whether star A or B is the DEB. While component A is classified as O7 III <cit.>, the spectral classification for component B is unavailable. As such, we have included it as a candidate in the YMDB. Additionally, its LC exhibits short-period variations compatible with pulsations.



HD 338961: 
No evidence of variability or binary nature was found for this star in the literature consulted. Despite some dispersion in its spectral classification, we consider the determination B0.5 IIInn from <cit.> to be reliable. Therefore, we have excluded it from inclusion in the YMDB.






 §.§ Contact binary systems discovered


We provide details of two targets initially considered for inclusion in the YMDB. However, subsequent analysis of the TESS data revealed that they are contact binaries, leading to their rejection.


LS V +38 12: 
This star was identified as a binary system, based on spectroscopic observations, by <cit.>, who classified the pair as O7 V((f))+ B0 III-V.



HD 305850: 
It is listed as a pulsating star in Simbad, yet we found no references to its periodicity or LC.  Moreover, its spectral type is actually uncertain, requiring an accurate determination.






 §.§ DEBs showing heartbeat phenomena
  






In the YMDB, in addition to TYC 8174-540-1, RAFGL 5223, HD 142152, BD+66 1675 and the subsystem found within BD+66 1674  which are already detailed in <ref>, several other systems exhibit variations in their LCs consistent with the heartbeat phenomenon. 


V* NY Cep: 
It is a well-known DEB <cit.>, characterized by a single eclipse resulting from its high eccentricity and periastron longitude. Notably, a bump is observed in its LC just before entering the eclipse, which we interpret as a probable heartbeat, although other proximity effects could not be discarded. 



HD 99898: 
This star is a previously identified variable, classified as Algol type <cit.>, and subsequently analyzed as a detached eclipsing binary (DEB) displaying apsidal motion <cit.>. However, no heartbeat phenomenon was reported prior to this study.



CD-28 5257: 
It is a recognized eclipsing binary <cit.>, although the reported periodicity is incorrect. Our analysis reveals it to be an eccentric double-eclipsing system, with eclipses corresponding to both transits and occultations. Furthermore, a marginal bump is observed in the orbital phases where the expected periastron occurs, which we identify as a heartbeat effect, made possible by the precision of the TESS data. This star is marked for exhibiting apsidal motion in Table <ref>, as the orbital phases of the secondary transits appear to vary across available Sectors 7, 8, 34, and 61.



V* V725 Car: 
It is a DEB in a highly eccentric orbit, with an RV curve determined for its primary component <cit.>. The new LC presented here reveals a distinct heartbeat effect between the closer eclipses. An integrated analysis of both datasets, spectroscopic and photometric, will enable the determination of stellar parameters for both components.



V* Y Cyg: 
It is a well-established SB2 eclipsing binary <cit.>. A marginal flux increase is observed at orbital phases corresponding to the periastron passage, which can be interpreted as a heartbeat effect, although it could also be attributed to ellipsoidal variations. This system exhibits apsidal motion, as evidenced by the variation in secondary eclipse times (and possibly their depths) between sectors 15 and 41 (∼700 d).



CD -27 4726: 
It is a very well-known DEB. TESS data show an asymmetric light increasing just before and after the (unique visible) primary eclipse. Its shape reminiscences the one detected on the very massive binary system WR 21a <cit.>. Unluckily, the spectral type of this DEB is not completely reliable thus, we add CD -27 4726 to the candidates list.



HD 306096: 
It is also a very studied DEB, but we note a clear heartbeat feature on the orbital phases where the periastron passage is expected.



V* GN Nor: 
This is another well-known system; however, its heartbeat has not been reported to the best of our knowledge.



*f Vel: 
It is a studied DEB system; however, as far as we are aware, its heartbeat phenomenon has not been documented.






 §.§ DEBs with newly determined periods
  

We compared the periods obtained from our analysis of the TESS data with those previously published. Some of the periods we derived were found to be twice as long as the previously reported ones. This discrepancy is likely due to the improved depth and resolution of the eclipses identified in the TESS data, enabling a clearer distinction between primary and secondary eclipses. While we have already discussed these discrepancies for system CD-28 5257 in Section <ref>, we now address those for which we have not yet provided details.


V* KU Car: 
KU Car, reported as B8 in SIMBAD from <cit.> and initially observed with a period of 5.92 days by <cit.>, is reexamined in the author's Master's dissertation <cit.>. The dissertation utilized ASAS-3 data, which revealed a secondary eclipse previously unnoticed, effectively revising the reported period to 2.96 days. Detailed spectral analysis was performed using data from the Galactic O Star Spectroscopic Survey <cit.>, which led to a revised spectral classification of B0.5 V(n). While O'Connell previously suggested an eccentric orbit based on the LC's shape outside the eclipses, this hypothesis was biased by an incomplete understanding of the LC's true nature. Our observations do not indicate eccentricity from the LC; only a comprehensive RV study could resolve this ambiguity. Ongoing spectroscopic analysis aims to refine the temperature and absolute parameters of KU Car.



HD 99630: 
This star is known as DEB, but its periodicity is badly reported in the literature <cit.>. Our analysis of the TESS LC reveals that its period is nearly double the previously reported value. Additionally, it shows double eclipses and high eccentricity. Our high-resolution spectra, obtained during the OWN Survey campaigns, indicate its spectral type is earlier than the B4-B5 determined by <cit.>. The ratio between He i λ4471 and Mg ii λ4481 lines is greater than three, and C ii λ4267 is identified, thus a B3-type is more suitable. The weakness of the metal lines confirms its dwarf class, B3 V.



CD-59 3165: 
It is recognized as a highly eccentric DEB <cit.>. In the TESS data, in addition to exhibiting double-eclipsing behavior, it also displays other eclipses with a periodicity of 3.18205 d (see Fig. <ref>). Given that this star is identified as a double in the WDS catalog (WDS J10348-6013AB), with both stars separated by only 2.3 arcseconds <cit.>, we can not definitively confirm the origin of these additional eclipses. However, our analysis of two spectra obtained during opposite quadratures of the main DEB system as part of the OWN Survey program reveals distinct spectral features for each component. One component exhibits narrow lines, while the other shows broader lines. The ratio of Si iii λ4552 to Si iv λ4089 is nearly unity in the narrow component, suggesting a spectral type of B0.5. Conversely, this ratio is smaller in the broader component, and considering the absence of He ii λ4542, yield a spectral type of B0-0.2. In terms of luminosity class, the narrow component appears to belong to classes III-I, as He ii λ4686 is markedly fainter than He i λ4713. Conversely, both lines are comparable in the broad component, indicative of a class V classification. In Fig. <ref> we show some spectral regions to illustrate these classifications.
In the future, we plan to conduct a more targeted spectral analysis to identify the spectral signatures of both components of the other DEB system. This analysis will also aim to elucidate the source of the additional eclipses observed.










 §.§ Other intriguing systems


We identified a diverse set of interesting systems, each presenting unique photometric characteristics that merit special mention. Firstly, several systems in our sample displayed Tidally Excited Oscillations (TEOs), which are oscillations within a star or system driven by tidal forces due to a close stellar companion. Notable examples include UCAC2 5911156 (detailed in <ref>, Fig. <ref>) along with V* V4386 Sgr, 2MASS J16542949-4139149, * 23 Ori, and V* V1216 Sco.

Additionally, we have successfully disentangled LCs of systems where the combined photometric data initially obscured individual components. Systems such as CD-59 3165 (detailed in <ref>, Fig. <ref>), BD+66 1674 (<ref>, Fig. <ref>) and eta Ori (detailed below, Fig. <ref>). 

We also observed systems whose eclipse characteristics vary significantly over time, adding another layer of complexity to their study. For instance, HD 278236 initially showed flat eclipses typical of total eclipses but over a span of 1500 days, the eclipse profile gradually transitioned to partial, with a smoothing and narrowing that indicates dynamic changes in the system. Additionally, the shifting of the secondary eclipse suggests apsidal motion, highlighting the system’s evolving orbital dynamics. (Fig. <ref>). HD 93683 exhibited a decrease in eclipse depth over 1500 days, suggesting changes in the system's configuration or surrounding material (Fig. <ref>). And BD+66 1675 (detailed in <ref>, Fig. <ref>) presented a heartbeat-like feature post-primary eclipse, which shifted to prior the eclipse over 1150 days of observation.


HD 93683: 
It is recognized as a SB2 system with a Be-type third component <cit.>. Our analysis of the TESS data unveiled intriguing behavior in this system, characterized by an attenuation of its eclipses (Fig. <ref>). This phenomenon could arise from variations in the brightness of the variable Be-star (either the star itself or its surrounding disk), which may dilute the eclipses. Alternatively, it could result from a rare effect known as Zeipel-Lidov-Kozai cycles, induced by the third component in a noncoplanar orbit, leading to changes in the orbital plane relative to our line of sight. As far as we know, this is the very first case reported in massive systems <cit.>.
This multiple system also exhibits several short-period variations, interpreted as pulsations. <cit.> reported one such variation with a period of 2.4 days. Further analysis is needed to determine additional periodicities in the TESS data.





* eta Ori: 
*Eta Ori is an ideal candidate for our catalog as it meets the criteria of a detached binary with spectral types B0.7 V and B1.5: V. <cit.> analyzed TESS data and fitted the LC to obtain absolute parameters for the system, opting not to use data from sector 32 due to its low variability amplitude. In our study, we present LCs using both sectors (Fig. <ref>). We confirm the periodic variations reported in the literature: the primary eclipse cycle at approximately 7.989 days, indicative of the detached eclipsing binary (EB) system, and a shorter cycle of about 0.432 days, associated with pulsational variability. <cit.> initially suggested that the shorter period might indicate a contact binary with a period of 0.864 days. Southworth proposed the configuration as a detached EB with a period of 7.988 days, where one component exhibits g-mode pulsations, alongside a noneclipsing binary with a period of 0.8641 days, showing strong ellipsoidal variations.










§ THE YOUNG MASSIVE DETACHED BINARY CATALOG


This work presents the comprehensive YMDB catalog, derived from the analysis of TESS LCs and spectroscopic data, with additional support from an extensive review of existing literature. 

Eclipsing binaries offer a unique opportunity to determine stellar parameters with high precision, especially when combining LC information with RV data. Detached binaries are particularly valuable due to the minimal interaction between their stellar components, allowing for accurate determinations of stellar parameters. 

Although systems within the O8–B3 spectral-type range are common, few have had their absolute parameters precisely measured. The YMDB catalog addresses this knowledge gap by providing a curated database of young massive detached binaries, facilitating high-precision stellar parameter determinations. 

Through the analysis of TESS LCs for 87 systems with suspected spectral types in the range O9-B1, this study identified 20 new eclipsing binaries, including 13 previously unknown variable systems and 2 nonthermal contact binaries. Additionally, new LC classifications were reported for 30 systems, and novel features such as eccentricity and heartbeat phenomena were discovered in many targets. 

The YMDB catalog offers a reliable resource of high-quality LCs, serving as a valuable asset for the astronomical community. The primary results of this study are documented in Table <ref>, which lists the 30 confirmed members of the YMDB. These systems feature detached LCs and have at least one component with a spectral classification within the specified range. For the 25 systems that show potential yet require further spectroscopic verification, details can be found in Table <ref>. The 32 systems that do not qualify for the catalog due to nondetached configurations or incompatible spectral types are included in Table <ref>

It is our desire to further calibrate this method and publish a semiautomatic pipeline for public use.





















RG acknowledges support from grant PICT 2019-0344. JIA acknowledges the financial support of DIDULS/ULS, through the project PR2324063. This research made use of Lightkurve, a Python package for Kepler and TESS data analysis (Lightkurve Collaboration, 2018). This work made use of Astropy:[http://www.astropy.org] a community-developed core Python package and an ecosystem of tools and resources for astronomy <cit.>.






aa

64
natexlab#1#1

[Albrecht et al.(2011)Albrecht, Winn, Carter, Snellen, &
  de Mooij]2011ApJ...726...68A
Albrecht, S., Winn, J. N., Carter, J. A., Snellen, I. A. G., & de
  Mooij, E. J. W. 2011, , 726, 68

[Alexander et al.(2016)Alexander, Hanes, Povich, &
  McSwain]2016AJ....152..190A
Alexander, M. J., Hanes, R. J., Povich, M. S., & McSwain, M. V. 2016,
  , 152, 190

[Alfonso-Garzón et al.(2012)Alfonso-Garzón, Domingo,
  Mas-Hesse, & Giménez]2012A   A...548A..79A
Alfonso-Garzón, J., Domingo, A., Mas-Hesse, J. M., & Giménez,
  A. 2012, , 548, A79

[Astropy Collaboration et al.(2022)Astropy Collaboration,
  Price-Whelan, Lim, Earl, Starkman, Bradley, Shupe, Patil,
  Corrales, Brasseur, Nöthe, Donath, Tollerud, Morris,
  Ginsburg, Vaher, Weaver, Tocknell, Jamieson, van Kerkwijk,
  Robitaille, Merry, Bachetti, Günther, Aldcroft,
  Alvarado-Montes, Archibald, Bódi, Bapat, Barentsen,
  Bazán, Biswas, Boquien, Burke, Cara, Cara, Conroy,
  Conseil, Craig, Cross, Cruz, D'Eugenio, Dencheva, Devillepoix,
  Dietrich, Eigenbrot, Erben, Ferreira, Foreman-Mackey, Fox,
  Freij, Garg, Geda, Glattly, Gondhalekar, Gordon, Grant,
  Greenfield, Groener, Guest, Gurovich, Handberg, Hart,
  Hatfield-Dodds, Homeier, Hosseinzadeh, Jenness, Jones, Joseph,
  Kalmbach, Karamehmetoglu, Kałuszyński, Kelley, Kern,
  Kerzendorf, Koch, Kulumani, Lee, Ly, Ma, MacBride, Maljaars,
  Muna, Murphy, Norman, O'Steen, Oman, Pacifici, Pascual,
  Pascual-Granado, Patil, Perren, Pickering, Rastogi, Roulston,
  Ryan, Rykoff, Sabater, Sakurikar, Salgado, Sanghi, Saunders,
  Savchenko, Schwardt, Seifert-Eckert, Shih, Jain, Shukla, Sick,
  Simpson, Singanamalla, Singer, Singhal, Sinha, Sipőcz,
  Spitler, Stansby, Streicher, Šumak, Swinbank, Taranu,
  Tewary, Tremblay, de Val-Borro, Van Kooten, Vasović, Verma,
  de Miranda Cardoso, Williams, Wilson, Winkel, Wood-Vasey, Xue,
  Yoachim, Zhang, Zonca, & Astropy Project
  Contributors]2022ApJ...935..167A
Astropy Collaboration, Price-Whelan, A. M., Lim, P. L., et al. 2022,
  , 935, 167

[Avvakumova et al.(2013)Avvakumova, Malkov, &
  Kniazev]2013AN....334..860A
Avvakumova, E. A., Malkov, O. Y., & Kniazev, A. Y. 2013, Astronomische
  Nachrichten, 334, 860

[Barbá et al.(2017)Barbá, Gamen, Arias, &
  Morrell]2017IAUS..329...89B
Barbá, R. H., Gamen, R., Arias, J. I., & Morrell, N. I. 2017, in
  The Lives and Death-Throes of Massive Stars, ed. J. J. Eldridge, J. C.
  Bray, L. A. S. McClelland, & L. Xiao, Vol. 329, 89–96

[Barbá et al.(2022)Barbá, Gamen, Martín-Ravelo,
  Arias, & Morrell]2022MNRAS.516.1149B
Barbá, R. H., Gamen, R. C., Martín-Ravelo, P., Arias, J. I.,
  & Morrell, N. I. 2022, , 516, 1149

[Bassino et al.(1982)Bassino, Dessaunet, Muzzio, &
  Waldhausen]1982MNRAS.201..885B
Bassino, L. P., Dessaunet, V. H., Muzzio, J. C., & Waldhausen, S.
  1982, , 201, 885

[Bodensteiner et al.(2020)Bodensteiner, Shenar, &
  Sana]2020A   A...641A..42B
Bodensteiner, J., Shenar, T., & Sana, H. 2020, , 641, A42

[Borkovits et al.(2022)Borkovits, Rappaport, Toonen, Moe,
  Mitnyan, & Csányi]2022MNRAS.515.3773B
Borkovits, T., Rappaport, S. A., Toonen, S., et al. 2022, , 515,
  3773

[Cannon & Pickering(1921)]1921AnHar..96....1C
Cannon, A. J. & Pickering, E. C. 1921, Annals of Harvard College
  Observatory, 96, 1

[Chini et al.(2012)Chini, Hoffmeister, Nasseri, Stahl, &
  Zinnecker]2012MNRAS.424.1925C
Chini, R., Hoffmeister, V. H., Nasseri, A., Stahl, O., & Zinnecker,
  H. 2012, , 424, 1925

[Crampton(1971)]1971AJ.....76..260C
Crampton, D. 1971, , 76, 260

[Crampton(1972)]1972MNRAS.158...85C
Crampton, D. 1972, , 158, 85

[Crampton & Fisher(1974)]1974PDAO...14..283C
Crampton, D. & Fisher, W. A. 1974, Publications of the Dominion
  Astrophysical Observatory Victoria, 14, 283

[de Mink et al.(2011)de Mink, Langer, &
  Izzard]2011BSRSL..80..543D
de Mink, S. E., Langer, N., & Izzard, R. G. 2011, Bulletin de la Societe
  Royale des Sciences de Liege, 80, 543

[Ekström(2021)]2021FrASS...8...53E
Ekström, S. 2021, Frontiers in Astronomy and Space Sciences, 8, 53

[Feast et al.(1961)Feast, Stoy, Thackeray, &
  Wesselink]1961MNRAS.122..239F
Feast, M. W., Stoy, R. H., Thackeray, A. D., & Wesselink, A. J. 1961,
  , 122, 239

[Feast & Thackeray(1963)]1963MmRAS..68..173F
Feast, M. W. & Thackeray, A. D. 1963, , 68, 173

[Gamen et al.(2007)Gamen, Barbá, Morrell, Arias,
  Maíz Apellániz, Sota, Walborn, &
  Alfaro]2007BAAA...50..105G
Gamen, R., Barbá, R., Morrell, N., et al. 2007, Spectroscopic
  monitoring of Southern Galactic O and WN stars: State of the Art in 2007

[Gamen et al.(2008)Gamen, Barbá, Morrell, Arias, &
  Maíz Apellániz]2008RMxAC..33...54G
Gamen, R., Barbá, R. H., Morrell, N. I., Arias, J., & Maíz
  Apellániz, J. 2008, in Revista Mexicana de Astronomia y Astrofisica
  Conference Series, Vol. 33, Revista Mexicana de Astronomia y Astrofisica
  Conference Series, 54–54

[Garrison et al.(1983)Garrison, Schild, &
  Hiltner]1983ApJS...52....1G
Garrison, R. F., Schild, R. E., & Hiltner, W. A. 1983, , 52, 1

[Georgelin et al.(1973)Georgelin, Georgelin, &
  Roux]1973AnA....25..337G
Georgelin, Y. M., Georgelin, Y. P., & Roux, S. 1973, , 25, 337

[Herrero et al.(1992)Herrero, Kudritzki, Vilchez, Kunze,
  Butler, & Haser]1992A   A...261..209H
Herrero, A., Kudritzki, R. P., Vilchez, J. M., et al. 1992, , 261,
  209

[Houk(1978)]1978mcts.book.....H
Houk, N. 1978, Michigan catalogue of two-dimensional spectral types for the
  HD stars

[Houk & Cowley(1975)]1975mcts.book.....H
Houk, N. & Cowley, A. P. 1975, University of Michigan Catalogue of
  two-dimensional spectral types for the HD stars. Volume I. Declinations -90_
  to -53_ƒ0.

[Khaliullin et al.(2006)Khaliullin, Antipin, &
  Khaliullina]2006AstL...32..772K
Khaliullin, K. F., Antipin, S. V., & Khaliullina, A. I. 2006, Astronomy
  Letters, 32, 772

[Kim et al.(2018)Kim, Kreiner, Zakrzewski, Ogłoza,
  Kim, & Jeong]2018ApJS..235...41K
Kim, C. H., Kreiner, J. M., Zakrzewski, B., et al. 2018, , 235, 41

[Kiminki & Smith(2018)]2018MNRAS.477.2068K
Kiminki, M. M. & Smith, N. 2018, , 477, 2068

[Lee et al.(1993)Lee, Sung, Koch, Hrivnak, Bradstreet,
  Corcoran, Mitchell, & Blitzstein]1993ASPC...38..239L
Lee, W. B., Sung, E. C., Koch, R. H., et al. 1993, in Astronomical
  Society of the Pacific Conference Series, Vol. 38, New Frontiers in Binary
  Star Research, ed. K.-C. Leung & I.-S. Nha, 239

[Li(2021)]2021ApJS..253...54L
Li, G.-W. 2021, , 253, 54

[Lightkurve Collaboration et al.(2018)Lightkurve Collaboration,
  Cardoso, Hedges, Gully-Santiago, Saunders, Cody, Barclay, Hall,
  Sagear, Turtelboom, Zhang, Tzanidakis, Mighell, Coughlin, Bell,
  Berta-Thompson, Williams, Dotson, & Barentsen]2018ascl.soft12013L
Lightkurve Collaboration, Cardoso, J. V. d. M., Hedges, C., et al.
  2018, Lightkurve: Kepler and TESS time series analysis in Python,
  Astrophysics Source Code Library

[Loden et al.(1976)Loden, Loden, Nordstrom, &
  Sundman]1976AnAS...23..283L
Loden, L. O., Loden, K., Nordstrom, B., & Sundman, A. 1976, , 23,
  283

[Lu(1985)]1985PASP...97..428L
Lu, W. 1985, , 97, 428

[Luger et al.(2016)Luger, Agol, Kruse, Barnes, Becker,
  Foreman-Mackey, & Deming]2016AJ....152..100L
Luger, R., Agol, E., Kruse, E., et al. 2016, , 152, 100

[Luger et al.(2018)Luger, Kruse, Foreman-Mackey, Agol, &
  Saunders]2018AJ....156...99L
Luger, R., Kruse, E., Foreman-Mackey, D., Agol, E., & Saunders, N.
  2018, , 156, 99

[Lynga(1964)]1964MeLuS.141....1L
Lynga, G. 1964, Meddelanden fran Lunds Astronomiska Observatorium Serie II,
  141, 1

[Maíz Apellániz et al.(2016)Maíz Apellániz,
  Sota, Arias, Barbá, Walborn, Simón-Díaz, Negueruela,
  Marco, Leão, Herrero, Gamen, & Alfaro]2016ApJS..224....4M
Maíz Apellániz, J., Sota, A., Arias, J. I., et al. 2016,
  , 224, 4

[Maíz Apellániz et al.(2019a)Maíz
  Apellániz, Trigueros Páez, Jiménez Martínez,
  Barbá, Simón-Díaz, Pellerin, Negueruela, & Rodrigo
  Souza Leão]2019hsax.conf..420M
Maíz Apellániz, J., Trigueros Páez, E., Jiménez
  Martínez, I., et al. 2019a, in Highlights on Spanish
  Astrophysics X, ed. B. Montesinos, A. Asensio Ramos, F. Buitrago,
  R. Schödel, E. Villaver, S. Pérez-Hoyos, &
  I. Ordóñez-Etxeberria, 420–420

[Maíz Apellániz et al.(2019b)Maíz
  Apellániz, Trigueros Páez, Negueruela, Barbá,
  Simón-Díaz, Lorenzo, Sota, Gamen, Fariña, Salas,
  Caballero, Morrell, Pellerin, Alfaro, Herrero, Arias, &
  Marco]2019AnA...626A..20M
Maíz Apellániz, J., Trigueros Páez, E., Negueruela, I.,
  et al. 2019b, , 626, A20

[Martín-Ravelo et al.(2021)Martín-Ravelo, Barbá,
  Morrell, Arias, & Gamen]2021mobs.confE..36M
Martín-Ravelo, P., Barbá, Rodolfo, H., Morrell, N., Arias, J.,
  & Gamen, R. 2021, in MOBSTER-1 virtual conference: Stellar Variability as
  a Probe of Magnetic Fields in Massive Stars, 36

[Martins & Palacios(2013)]2013A   A...560A..16M
Martins, F. & Palacios, A. 2013, , 560, A16

[Mason et al.(2001)Mason, Wycoff, Hartkopf, Douglass, &
  Worley]2001AJ....122.3466M
Mason, B. D., Wycoff, G. L., Hartkopf, W. I., Douglass, G. G., &
  Worley, C. E. 2001, , 122, 3466

[Munari & Tomasella(1999)]1999A   A...343..806M
Munari, U. & Tomasella, L. 1999, , 343, 806

[Muzzio & Orsatti(1977)]1977AJ.....82..474M
Muzzio, J. C. & Orsatti, A. M. 1977, , 82, 474

[O'Connell(1951)]1951PRCO....2...85O
O'Connell, D. J. K. 1951, Publications of the Riverview College Observatory,
  2, 85

[Plaskett(1908)]1908ApJ....28..266P
Plaskett, J. S. 1908, , 28, 266

[Pojmanski(1998)]1998AcA....48...35P
Pojmanski, G. 1998, , 48, 35

[Pourbaix et al.(2004)Pourbaix, Tokovinin, Batten, Fekel,
  Hartkopf, Levato, Morrell, Torres, & Udry]2004A   A...424..727P
Pourbaix, D., Tokovinin, A. A., Batten, A. H., et al. 2004, , 424,
  727

[Pozo Nuñez et al.(2019)Pozo Nuñez, Chini, Barr
  Domínguez, Fein, Hackstein, Pietrzyński, &
  Murphy]2019MNRAS.490.5147P
Pozo Nuñez, F., Chini, R., Barr Domínguez, A., et al. 2019,
  , 490, 5147

[Sana et al.(2012)Sana, de Mink, de Koter, Langer,
  Evans, Gieles, Gosset, Izzard, Le Bouquin, &
  Schneider]2012Sci...337..444S
Sana, H., de Mink, S. E., de Koter, A., et al. 2012, Science, 337, 444

[Sawyer(1887)]1887AJ......7..116S
Sawyer, E. F. 1887, , 7, 116

[Shi et al.(2022)Shi, Qian, & Li]2022ApJS..259...50S
Shi, X.-d., Qian, S.-b., & Li, L.-J. 2022, , 259, 50

[Simón-Díaz et al.(2011)Simón-Díaz, Castro,
  Garcia, Herrero, & Markova]2011BSRSL..80..514S
Simón-Díaz, S., Castro, N., Garcia, M., Herrero, A., &
  Markova, N. 2011, Bulletin de la Societe Royale des Sciences de Liege, 80,
  514

[Sota et al.(2014)Sota, Maíz Apellániz, Morrell,
  Barbá, Walborn, Gamen, Arias, & Alfaro]2014ApJS..211...10S
Sota, A., Maíz Apellániz, J., Morrell, N. I., et al. 2014,
  , 211, 10

[Sota et al.(2011)Sota, Maíz Apellániz, Walborn,
  Alfaro, Barbá, Morrell, Gamen, & Arias]2011ApJS..193...24S
Sota, A., Maíz Apellániz, J., Walborn, N. R., et al. 2011,
  , 193, 24

[Southworth & Bowman(2022)]2022MNRAS.513.3191S
Southworth, J. & Bowman, D. M. 2022, , 513, 3191

[Thompson et al.(2012)Thompson, Everett, Mullally,
  Barclay, Howell, Still, Rowe, Christiansen, Kurtz, Hambleton,
  Twicken, Ibrahim, & Clarke]2012ApJ...753...86T
Thompson, S. E., Everett, M., Mullally, F., et al. 2012, , 753, 86

[Trigueros Páez et al.(2021)Trigueros Páez, Barbá,
  Negueruela, Maíz Apellániz, Simón-Díaz, &
  Holgado]2021A   A...655A...4T
Trigueros Páez, E., Barbá, R. H., Negueruela, I., et al. 2021,
  , 655, A4

[Turner(1980)]1980ApJ...235..146T
Turner, D. G. 1980, , 235, 146

[Vijapurkar & Drilling(1993)]1993ApJS...89..293V
Vijapurkar, J. & Drilling, J. S. 1993, , 89, 293

[Walborn & Fitzpatrick(1990)]1990PASP..102..379W
Walborn, N. R. & Fitzpatrick, E. L. 1990, , 102, 379

[Weidner & Vink(2010)]2010A   A...524A..98W
Weidner, C. & Vink, J. S. 2010, , 524, A98

[Welsh et al.(2011)Welsh, Orosz, Aerts, Brown,
  Brugamyer, Cochran, Gilliland, Guzik, Kurtz, Latham, Marcy,
  Quinn, Zima, Allen, Batalha, Bryson, Buchhave, Caldwell,
  Gautier, Howell, Kinemuchi, Ibrahim, Isaacson, Jenkins, Prsa,
  Still, Street, Wohler, Koch, & Borucki]2011ApJS..197....4W
Welsh, W. F., Orosz, J. A., Aerts, C., et al. 2011, , 197, 4













§ FOLDED LIGHT CURVES FOR ALL 87 SYSTEMS ANALYZED IN THIS STUDY


The appendix presents the LCs for all systems analyzed in this study. Figure <ref> shows the LCs for the confirmed systems, Figure <ref> displays the LCs for the candidate systems, and Figure <ref> includes the LCs for the unqualified systems.















