Open Access
Issue
A&A
Volume 667, November 2022
Article Number A100
Number of page(s) 18
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/202244730
Published online 11 November 2022

© A. Schootemeijer et al. 2022

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe-to-Open model. Subscribe to A&A to support open access publication.

1. Introduction

The early Universe contained a small amount of heavy elements, and the same is true for dwarf galaxies (Mateo 1998). As such, nearby dwarf galaxies are uniquely suitable for studying the early generations of massive stars, the stars that could have reionized the Universe (Dayal & Ferrara 2018) and whose feedback controlled early star formation. As a result of the low abundance of heavy elements (referred to as metallicity, Z), massive stars (stars that are at least eight times more massive than the Sun) have weaker radiation-driven stellar winds (Hainich et al. 2015; Vink & Sander 2021). Weak stellar winds lead to low angular momentum loss, and allow massive low-metallicity stars to maintain higher rotation velocities (Mokiem et al. 2007; Ramachandran et al. 2019; Brott et al. 2011) as well as higher masses prior to core collapse, favoring the formation of massive black holes. Interestingly, superluminous supernovae (Lunnan et al. 2014), long-duration gamma-ray bursts (Savaglio et al. 2009), and ultra-luminous X-ray emission from compact sources (Kaaret et al. 2017) primarily occur in low-metallicity dwarf galaxies; furthermore, all have been linked to rapidly rotating massive stars (Aguilera-Dena et al. 2018; Marchant et al. 2017).

Rapid rotation is known to affect the structure, appearance, and evolution of massive stars in various ways. First, the closer stars are to their critical rotation velocity, the more they expand in the equatorial direction. Then, gravity darkening reduces their effective temperatures near the equator (von Zeipel 1924; Espinosa Lara & Rieutord 2011), which makes rotating stars look redder. Second, as stars rotate at or close to their critical rotation velocity, they can acquire a decretion disk (Rivinius et al. 2013). This disk emits more light in the red than in the blue, providing a second mechanism that makes an extremely rapidly rotating source intrinsically redder (Martayan et al. 2010). Because of the emission lines that form in the disk, these stars are classified as OBe stars (Struve 1931). Finally, rotation is thought to induce internal mixing. If efficient enough, this could lead to chemically homogeneous evolution, where stars remain compact and evolve to ever-higher temperatures and luminosities (Maeder 1987; Langer 1992). This type of evolution would open an extra window to gravitational wave events caused by black hole mergers (de Mink & Mandel 2016; Marchant et al. 2016). Also, massive stars experiencing chemically homogeneous evolution could produce high-energy photons (Szécsi et al. 2015), providing a potential explanation for nebular emission of, for example, doubly ionized helium that is observed in high-redshift, metal-poor galaxies (Erb et al. 2010).

In the Milky Way and its satellite dwarf galaxies, the Large Magellanic Cloud (LMC; 0.5 Z Trundle et al. 2007) and the Small Magellanic Cloud (SMC; 0.2 Z Korn et al. 2000), rotation velocities of large samples of massive stars have been inferred spectroscopically from line-broadening (e.g., Mokiem et al. 2006; Ramachandran et al. 2019). These studies found that the average rotation velocity of massive stars modestly increases with decreasing metallicity. This trend might continue towards lower metallicity; but whether it does is uncertain, because of the currently scarce statistics in the more metal-poor star-forming galaxies, which results from the increased cost of spectroscopic observations of these distant systems (Telford et al. 2021; Garcia et al. 2021).

An alternative method to line broadening measurements for gauging stellar rotation, which is cheaper and does not suffer from inclination effects, involves the measurement of the fraction of OBe stars. There, we distinguish between two different OBe star fractions: the OBe star fraction measured in subsamples of stars, and the galaxy-wide OBe star fraction. A galaxy-wide OBe star fraction has so far never been measured. In previous work, one approach to measure OBe star fractions was to use photometry to detect excess infrared (IR) emission from decretion disks of OBe stars in a subsample of stars with known spectral types in the Magellanic Clouds (Bonanos et al. 2009, 2010). An alternative approach is to use narrow-band Hα photometry, which traces Hα emission from the disk of the OBe star in order to measure the OBe star fraction in star clusters of the Magellanic Clouds (Keller et al. 1999; Wisniewski & Bjorkman 2006; Iqbal & Keller 2013) and the Milky Way (McSwain & Gies 2005), as well as in parts of the Andromeda galaxy (Peters et al. 2020). Measured cluster OBe star fractions tend to be higher at low metallicity, but within a galaxy the obtained results depend strongly on the cluster that is picked and which of its stars are included in the measurement (see below). We note that the studied clusters are often older than the maximum lifetime of massive stars – which are the focus of this work–, such that they mainly probe longer-lived low-mass stars. Also, measurements in clusters that are young enough to still contain massive stars might be biased because they can have lost OBe stars produced via the binary channel through supernova kicks (as found by Dallas et al. 2022). Therefore, only a measurement of the galaxy-wide OBe star fraction can provide an unbiased approach. The goals of this work are twofold: (i) we aim to measure this galaxy-wide fraction of bright OBe stars in the Magellanic Clouds, and (ii) to extend the measurement of the galaxy-wide OBe star fraction to subSMC metallicity.

2. Methods and data sets

In this work, we make use of archival broad-band photometry. We use short-wavelength (near-ultraviolet(NUV)/blue optical) data to identify hot stars – which OBe stars are – and longer-wavelength IR data to establish the presence of a disk. In this section, we explain how we construct our data sets, which we analyze in Sect. 3. In Appendix A we provide various tests, where we discuss in detail how complete our data sets are, and how well we can use them to identify OBe stars.

Magellanic Clouds. For both the SMC and the LMC, with CDS-xmatch1 (Boch et al. 2012; Pineau et al. 2020) we cross-correlated sources from UBVI catalogs (Zaritsky et al. 2002, 2004) to Spitzer-SAGE (Meixner et al. 2006) sources within 1″. The resulting data files were then cross-correlated with Gaia EDR3 (Gaia Collaboration 2021) sources, with again a 1″ cross-correlation radius. It can happen that a source from a UBVI catalog is cross-matched to more than one source that would appear in our color–magnitude diagrams (CMDs; i.e., brighter than MGbp = −2, and a listed magnitude for each of the following filters: U, B, Gbp, Grp, J, and [3.6]). In those cases, we take the closest source. As a result, we discard 110 sources in the SMC and 353 in the LMC. These cases are relatively few in number, because the SMC and LMC CMDs contain of the order of tens of thousands of sources. Therefore, we do not expect cross-correlation to the wrong sources to affect our results.

The cross-correlated catalogs are cleaned of foreground sources. For the SMC, the method is similar to the one from Schootemeijer et al. (2021): sources with a Gaiaparallax_over_parallax_error value larger than five are removed, as are sources that do not fit the proper motion criteria. More specifically, for the SMC we remove sources more than 2σ away from its bulk proper motion (μ), for which we adopt the values from Yang et al. (2019): μRA = 0.695 ± 0.240 mas yr−1 (Right Ascension (RA)) and μDec = −1.206 ± 0.140 mas yr−1 (Declination (Dec)). For the LMC, we take the same approach but for its bulk motion we adopt μRA = 1.750 ± 0.375 mas yr−1 and μDec = 0.250 ± 0.375 mas yr−1, based on Fig. 16 from Gaia Collaboration (2018).

Sextans A. For Sextans A, no sufficiently deep Gaia and IR data are available, and ground-based UBVI data lack the required spatial resolution. Therefore, we use data obtained with the Hubble Space Telescope (HST). There is F336W and F439W photometry available2 (Bianchi et al. 2012, hereafter referred to as B12). These filters are centered on similar wavelengths to the U and B filters. B12 also registered photometry in the F555 and F814W filters, which have central wavelengths similar to the GaiaGbp and Grp filters. However, in the F555 and F814W filters additional, deeper observations exist3 (Holtzman et al. 2006, hereafter H06). Both B12 and H06 contain two HST fields in Sextans A, which cover about two-thirds of the galaxy. To separate MS, OBe, and HeB stars in Sextans A, we cross-correlate B12 with H06 using TOPCAT (Taylor 2005), and check the result with the Aladin sky atlas (Bonnarel et al. 2000; Boch & Fernique 2014). We note that with a distance modulus DM = 25.63 (Tammann et al. 2011), corresponding to a distance of 1.3 Mpc, Sextans A is about 20 times further away than the Magellanic Clouds, which makes the cross-correlation a more difficult task. Given a systematic coordinate offset in B12 and H06 data, we shift the B12 data before cross-matching. For data in frame 4 of B12, this shift is 1.321″ in RA and 0.265″ in Dec. For their frame 5, the shift is 0.752″ in RA and 1.087″ in Dec.

Holmberg I and Holmberg II. Our strategy to identify OBe stars in Holmberg I and Holmberg II is the same as for Sextans A, except that we employ data from the the Legacy ExtraGalactic UV Survey4 (LEGUS; Sabbi et al. 2018), which are also obtained with HST. In LEGUS, the F336W, F438W, F555W, and F814W filters are for simplicity referred to as the U, B, V, and I filters, respectively. In this work, we use the same nomenclature for LEGUS data. LEGUS UBVI data cover most of Holmberg I and about one-third of Holmberg II.

3. Results

3.1. Small Magellanic Cloud

We show the CMD of the sources in our SMC data set in the left panel of Fig. 1. On the right, the red supergiant branch can be seen, and on the left there are two bands of stars with blue colors, which we focus on in this work. These two blue bands lie at Gbp − Grp ≈ −0.25 and at Gbp − Grp ≈ 0.

thumbnail Fig. 1.

Left: color magnitude diagram of stars in the SMC, constructed with Gaia filters. Each black dot represents a source. The brightest sources are at the top and the bluest sources are on the left. We adopt a value of 18.91 for the distance modulus (Hilditch et al. 2005). Right: Stacked histograms of the Gbp − Grp color in the absolute magnitude ranges that are written in the top right of each panel. Main sequence stars are shown in magenta, OBe stars in light orange, and He-burning stars in green. The black line shows the number distribution of all the sources combined.

To investigate the nature of the sources in the two blue components, we plot the IR J − [3.6] color against the U − B color in the left panel of Fig. 2 for all sources in our SMC data set brighter than an absolute magnitude of MGbp = −2. Three distinct groups form, which can be interpreted with the aid of the location of stars with known spectral types in the right panel of Fig. 2. About 90% of the sources in the bottom left group (i.e., with the bluest colors) that have known spectral types are O-type or early B-type stars. Such spectral types are typical for massive main sequence stars without decretion disks (hereafter referred to as MS stars, for simplicity). Similarly, about 90% of the sources in the group with a blue U − B color but a red color in the IR (photometric OBe stars) on the right side of Fig. 2 are known to be Oe and early Be stars (see also histograms). Their location in the color–color diagram can be understood from their decretion disks emitting mainly at longer wavelengths, giving OBe sources a much redder color in the IR (Bonanos et al. 2009, 2010), while their near-UV/blue optical U − B color remains largely unaffected. Finally, at the top of the color–color diagrams there is a group of sources with U − B > −0.5. Most of these sources have spectral types of B6 or later (Fig. 2, right), implying that they have effective temperatures of Teff ≲ 15 kK (Schootemeijer et al. 2021). Such effective temperatures explain the relatively red U − B colors (compared to O and early B stars) and imply that the sources in this group are He-burning (HeB) objects, as the MS for massive SMC stars extends only up to ∼20 kK even for extreme assumptions (Schootemeijer et al. 2019, 2021).

thumbnail Fig. 2.

Left: color–color diagram of sources in our SMC data set that are brighter than MGbp = −2. Histograms are shown in the right panel and the top panel, where the top histogram includes only sources bluer than U − B = −0.5. Right: same as the left side, but showing sources from various spectroscopic studies (Martayan et al. 2007; Hunter et al. 2008; Dufton et al. 2019; Bonanos et al. 2010) that are colored for their spectral type.

The spectroscopic confirmation of the nature of the sources in the U − B vs. J − [3.6] diagram implies that it is a very effective tool for identifying the three types of blue massive stars in the SMC: MS (O/early B) stars, HeB (late B+) stars, and OBe stars.

On the right side of the CMD (Fig. 1), we show histograms for MS, HeB, and OBe stars, as identified by their colors in Fig. 2. The bluest band at Gbp − Grp ≈ −0.25 is almost exclusively made up of photometric MS stars. The HeB sources occupy a broad range of colors, while OBe stars concentrate at Gbp − Grp ≈ 0 and give rise to a peak in the number distribution there. This color difference of a quarter of a magnitude between MS and OBe stars was also found in previous studies where filters were centered at similar wavelengths (Martayan et al. 2010; Milone et al. 2018; Bodensteiner et al. 2020a). For sources brighter than MGbp = −5, the number of OBe stars becomes very small, and the parallel feature starts to fade away. The total number of OBe stars shown in the CMD is large: slightly more than 2500 (Table 1).

Table 1.

Number of OBe stars (NOBe) and OBe star fractions (fOBe) in the LMC, SMC, HoI, HoII, and SexA.

3.2. The Large Magellanic Cloud, Holmberg I and II, and Sextans A

We extend our analysis to four more dwarf galaxies, starting with the LMC. As well as the SMC, the LMC is a satellite galaxy of the Milky Way, and is slightly closer. However, with half the Solar metallicity (Trundle et al. 2007) it is more enriched in heavy elements than the SMC. We can apply the same method because the LMC is covered by the data sets that we used for the SMC. We further describe this in Appendix A, where we also discuss the CMD for the LMC (Fig. A.1) and its color–color diagram (Fig. A.2). We find as many as approximately 4500 photometric OBe stars in the LMC that are brighter than MGbp = −2 (Table 1).

Furthermore, we explore the dwarf galaxies Sextans A, Holmberg I (Ho I), and Holmberg II (Ho II), which all have a lower metallicity than the SMC and are much more distant. Again, a more detailed description can be found in Appendix A. The CMDs of these galaxies are shown in Figs. A.3, A.8, and A.10, and the color–color diagrams are shown in Figs. A.4, A.9, and A.11. Because we use different filters to identify OBe star disks in these dwarf galaxies, in Appendix A.1 we provide multiple extra tests to demonstrate that with broadband photometry from HST, we can also distinguish O/early B stars, late B+ stars, and OBe stars from one another. These tests involve stars with known spectral types, and narrow-band Hα photometry. Among these three dwarf galaxies, the most metal-poor one that we explore is Sextans A, with a metallicity Z ≈ 1/10 Z (Kaufer et al. 2004; Garcia et al. 2017). It is a distance of about 1.3 Mpc (Tammann et al. 2011). The final two dwarf galaxies, Ho I and Ho II, are at distances of 3 − 4 Mpc (Sabbi et al. 2018). As such, these objects are beyond the Local Group, and belong to the M81 Group instead (Karachentsev et al. 2002). According to their oxygen abundances, Ho I and Ho II have metallicities of ∼1/6 Z and ∼1/7 Z, respectively (Croxall et al. 2009; Bergemann et al. 2021). In each of these three dwarf galaxies, which are smaller and less massive than the Magellanic Clouds, we identify of the order of 100 OBe stars brighter than Mabs = −3: 62 in Sextans A, 83 in Ho I, and 240 in Ho II (Table 1).

3.3. Magnitude-, mass-, and metallicity-dependent OBe star fractions

To calculate the OBe star fraction in different absolute magnitude bins, we divide the number of OBe stars, NOBe, by the total number of H-burning stars, NOBe + NMS (selected by the color cuts in Figs. 2, A.2, A.4, A.9, and A.11). When constructing Fig. 3, we assumed that the disk adds no significant flux to the Gbp, V, or F555W filter (which is what we find in the Appendix A.4). MS sources in the SMC that are dimmer than Mabs ≈ −3 are not bright enough to be observed in the Spitzer [3.6] filter. For the LMC and Sextans A, the completeness fractions of OBe stars and nonOBe MS stars also start to significantly deviate at a similar absolute magnitude (Appendix A.2). Therefore, we cannot accurately calculate the OBe star fraction at Mabs ≳ −3. At the bright end (Mabs < −5.5), we notice that the fraction of Wolf-Rayet and B[e] stars posing as OBe stars becomes significant, in particular in the LMC (Appendix A.3). Therefore, we only display the absolute magnitude range −3 > Mabs > −5.5 in Fig. 3. The values we find for NOBe and fOBe are listed in Table 1. At the top of Fig. 3, we write in each absolute magnitude interval the average evolutionary mass that has been determined previously (Schootemeijer et al. 2021) for SMC stars of type B5 and earlier (see Appendix B). The obtained masses indicate that the sources shown in Fig. 3 are typically massive stars with masses up to 25 M. Given the spread in evolutionary masses and potential biases, these masses are meant to serve as an indication.

thumbnail Fig. 3.

Galaxy-wide OBe star fraction in the Magellanic Clouds, Sextans A, and Holmberg I and II, measured in absolute magnitude ranges −3 > Mabs > −3.5, −3.5 > Mabs > −4, and so on. The tick symbols at the top of the plot indicate the average evolutionary mass that has been inferred (Schootemeijer et al. 2021) for stars of spectral type B5 and earlier in each absolute magnitude bin in the SMC.

We find that the OBe star fraction is higher in the SMC than in the LMC. In the LMC, fOBe ≈ 0.25 at 10 M, but then it decreases quickly for M ≳ 15 M. In the SMC, fOBe ≈ 0.35 around 10 M, and it decreases to ∼0.20 around 20 M. The relatively large difference in OBe star fraction between the SMC and LMC towards MGbp = −5 could be the result of stellar winds (which are stronger at higher metallicity) spinning down the bright stars more strongly in the LMC. However, for disk angular momentum loss of OBe stars (Okazaki 2001), the metallicity dependence is not as well understood. With further decreasing metallicity in Ho I, Ho II, and Sextans A, we find OBe star fractions that are comparable to what we find in the SMC. This becomes more apparent in Fig. 4, where we plot our galaxy-wide OBe star fractions against metallicity. In all four galaxies in the regime 1/10 ≤ Z/Z ≤ 1/5, about 30% of the stars in the considered magnitude range are OBe stars. This implies that extremely rapid rotation is commonplace in metal-poor environments.

thumbnail Fig. 4.

Metallicity (Z) trends of the OBe star fraction. Top panel: in the low-Z regime, we show the galaxy-wide OBe star fraction calculated in this work in the absolute magnitude interval −3 > Mabs > −5 with large diamonds. Bottom panel: literature studies of stellar subsamples within different galaxies, which targeted cluster and/or field stars: K99 (Keller et al. 1999), M99 (Maeder et al. 1999), M05 (McSwain & Gies 2005), W06 (Wisniewski & Bjorkman 2006), I13 (Iqbal & Keller 2013), Bonanos (Bonanos et al. 2009, 2010), P20 (Peters et al. 2020), and ZB97 (Zorec & Briot 1997). We only show clusters with 20 or more OB+OBe stars, and we do not show multiple measurements for any cluster. In K99, M99, M05, W06, I13, and P20, narrow-band photometry was used to identify OBe stars. For a part of the sample in M99, OBe stars were identified with spectroscopy. Bonanos used IR excess of stars with known spectral types to identify OBe stars. ZB97 used a sample of stars with known spectral types. For reference, we show the results of this work again with a dotted line. The metallicities in Holmberg I and II (Ho I and II) are estimated using oxygen abundances (Croxall et al. 2009). The two highest-metallicity galaxies are our own Milky Way (MW) and Andromeda (M 31). For ZMW, we adopt the Solar metallicity, and we adopt ZM31 = 1.5 Z from P20.

In the SMC and LMC, we can discuss our results in the context of literature values reported for stellar subsamples. However, we caution that a direct comparison is inappropriate because of different sample selection, and time-evolution of cluster OBe star fractions. Indeed, as a result of this Fig. 4 shows strong variation from cluster to cluster within individual galaxies. Also, investigated clusters are often older than 40 Myr, which is roughly the lifetime of an 8 M star. Therefore, they contain mainly B and Be stars below 8 M rather than massive stars, which are the topic of this work. The LMC and SMC OBe star fractions that we find are above the average literature values for cluster studies. They are also higher than results for subsamples of stars with known spectral types (green circles), possibly because spectral type catalogs can biased towards the brighter stars, for which we find lower OBe star fractions (Fig. 3).

3.4. Implications for OBe star formation

Whether OBe stars are mainly produced by single star evolution (van Bever & Vanbeveren 1997) or through binary evolution (Schootemeijer et al. 2018; Klement et al. 2019; Bodensteiner et al. 2020b) where mass transfer from one star to the other leads to spin-up and the production of an OBe star is currently a topic of debate. Here we address this issue, and briefly discuss the high value of fOBe ≈ 0.3 in the SMC from a theoretical perspective. Previous work, which used the rotational velocity distribution of Dufton et al. (2013) as the initial distribution, has shown that the single star channel struggles to produce the (high) observed fraction of OBe stars in SMC star cluster NGC 330 (Hastings et al. 2020). Here we redo the analysis of Hastings et al. (2020), but we assume that the theoretical population undergoes constant star formation (CSF) instead of it being coeval; for details, see Appendix C.1. Figure 5 shows that for CSF, the single star population can match the observed OBe star fraction as long as v/vcrit ≳ 0.7 (vcrit being the breakup velocity) is sufficient for stars to manifest themselves as OBe stars. However, if stars have to rotate at v/vcrit ≳ 0.8 − 0.9 to display the OBe phenomenon (Townsend et al. 2004; Frémat et al. 2005), the single star channel cannot explain the high observed OBe star fractions in the SMC in our simulations.

thumbnail Fig. 5.

Comparison of the observed absolute-magnitude-dependent SMC OBe star fraction with theory. At the top of each panel we show indicative masses, as in Fig. 3. Top: OBe star fraction in a theoretical single star population, for which we explore three different minimum v/vcrit values that stars require to manifest themselves as OBe stars, with v/vcrit being the fraction of the critical rotation velocity. Bottom: OBe star fraction in theoretical binary populations. With blue and orange lines, we show theoretical predictions derived from ComBinE simulations, for different combinations of the mass transfer efficiency β and binary fraction fbin. In the theoretical populations, single stars are assumed to not contribute to producing OBe stars. For each theoretical binary population, we assume that all stars with v/vcrit > 0.7 are OBe stars. In the Appendix, we provide a similar figure, but where v/vcrit > 0.8 and v/vcrit > 0.9 are set as the required rotation velocities for OBe stars in the binary populations (Fig. C.1).

On the other hand, the binary channel also struggles to explain the observed OBe star fraction. In the SMC (and Ho I, Ho II, and Sextans A), the OBe star fraction is so high that it touches the stringent upper limit of fOBe ≲ 1/3 at and below the turnoff in coeval star clusters (Hastings et al. 2021). However, for a galaxy-wide population undergoing CSF, efficient mass transfer can help to increase the OBe star fraction shown in Fig. 3 because it can turn lower mass stars – which are favored by the initial mass function (IMF) – into brighter OBe stars. To test this, we simulate three SMC binary populations with the rapid population synthesis code ComBinE (Kruckow et al. 2018), varying the mass transfer efficiency (Fig. 5). We elaborate on these simulations in Appendix C.2. In the theoretical population, for inefficient to moderately efficient mass transfer (β ≤ 0.5), the binary channel cannot explain an OBe star fraction as high as 0.3, even for 90% of sources being interacting binaries. However, for conservative mass transfer, it does become possible to reproduce fOBe ≈ 0.3. As is the case for the single star simulations, our binary simulations do a good job at explaining the behavior of the OBe star fraction as a function of absolute magnitude, which strengthens our results. In principle, a scenario in which stars in binaries spin up to become OBe stars before mass exchange takes place could help to explain the high observed OBe star fractions. However, recent work found no clear evidence for any MS companion in a sample of over 250 Galactic Be stars (Bodensteiner et al. 2020b), suggesting that OBe star formation before interaction in binaries rarely takes place.

4. Discussion and conclusions

We find that in color–color diagrams of metal-poor dwarf galaxies, MS stars, helium-burning stars, and OBe stars form well-separated groups of sources. We took advantage of this to identify OBe stars based on broad-band photometry, and were able to apply this method to far-away galaxies. We analyzed the Milky Way satellite dwarf galaxies LMC (50 kpc; 1/2 Z) and SMC (60 kpc; 1/5 Z), and three distant galaxies: Sextans A (1.3 Mpc away, at the border of the Local Group; 1/10 Z), and Holmberg I (3.8 Mpc; 1/6 Z) and Holmberg II (3.2 Mpc; 1/7 Z) in the M 81 group. Using archival photometric data, we were able to find of the order of 100 and up to a few thousand OBe stars in each galaxy. This allowed us to calculate galaxy-wide OBe star fractions.

We find that the OBe star fraction increases with decreasing metallicity between the LMC and the SMC, and stays high – around 30% – in Sextans A, Holmberg I, and Holmberg II. In other words, extremely rapidly rotating massive stars are common in low-metallicity dwarf galaxies. These high OBe star fractions shed light on their formation channel. By comparison with SMC population synthesis models, we find that the single-star channel can only reproduce the high observed OBe star fractions if significantly subcritical rotation suffices to display the OBe phenomenon. For the binary channel to produce a sufficient number of bright OBe stars, the binary fraction has to be close to unity and – at the same time – mass transfer has to be efficient.

Our work provides galaxy-wide OBe star fraction measurements in the Magellanic Clouds, and we extended OBe star fraction measurements to subSMC metallicity and beyond the Local Group. Within the now-accessible volume, our method can be applied to even more metal-poor galaxies to further test the behavior of the OBe star fraction at low Z. There are a number of such galaxies that are currently star forming (Garcia et al. 2021), but would require near-UV observations in addition to the existing optical photometry (e.g., ANGST – Dalcanton et al. 2009). Examples of promising targets are Sextans B (1/10 Z–1/15 Z: Kniazev et al. 2005), Leo A (∼1/20 Z: van Zee et al. 2006), and SagDIG (∼1/20 Z: Saviane et al. 2002). Fortunately, such efforts are currently being undertaken (e.g., Sabbi et al. 2020). In the galaxies up to the distance of ∼4 Mpc reached here, OBe stars in the range −3 > Mabs > −5 did not hit the HST detection limit. Therefore, even more distant galaxies could be studied, as long as the data quality remains high and color–color diagrams similar to the ones verified in this work are used to identify OBe stars. Exploiting such galaxies with sub-Sextans A metallicities would allow us to get even closer to measuring the incidence of extremely rapid rotation among massive stars in environments that resemble the primordial Universe.


2

Which we accessed through the Hubble Legacy Archive: https://hla.stsci.edu/hlaview.html

Acknowledgments

We thank the anonymous referee for their careful reading of the manuscript and for useful comments. M.G. gratefully acknowledges support by grants PID2019-105552RB-C41 and MDM-2017-0737 Unidad de Excelencia "María de Maeztu"-Centro de Astrobiología (CSIC-INTA), funded by the Spanish Ministry of Science and Innovation/State Agency of Research MCIN/AEI/10.13039/501100011033.

References

  1. Aguilera-Dena, D. R., Langer, N., Moriya, T. J., & Schootemeijer, A. 2018, ApJ, 858, 115 [NASA ADS] [CrossRef] [Google Scholar]
  2. Arenou, F., Luri, X., Babusiaux, C., et al. 2018, A&A, 616, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Bergemann, M., Hoppe, R., Semenova, E., et al. 2021, MNRAS, 508, 2236 [NASA ADS] [CrossRef] [Google Scholar]
  4. Bianchi, L., Efremova, B., Hodge, P., Massey, P., & Olsen, K. A. G. 2012, AJ, 143, 74 [NASA ADS] [CrossRef] [Google Scholar]
  5. Boch, T., & Fernique, P. 2014, in Astronomical Data Analysis Software and Systems XXIII, eds. N. Manset, & P. Forshay, ASP Conf. Ser., 485, 277 [NASA ADS] [Google Scholar]
  6. Boch, T., Pineau, F., & Derriere, S. 2012, in Astronomical Data Analysis Software and Systems XXI, eds. P. Ballester, D. Egret, & N. P. F. Lorente, ASP Conf. Ser., 461, 291 [NASA ADS] [Google Scholar]
  7. Bodensteiner, J., Sana, H., Mahy, L., et al. 2020a, A&A, 634, A51 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Bodensteiner, J., Shenar, T., & Sana, H. 2020b, A&A, 641, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Bonanos, A. Z., Massa, D. L., Sewilo, M., et al. 2009, AJ, 138, 1003 [Google Scholar]
  10. Bonanos, A. Z., Lennon, D. J., Köhlinger, F., et al. 2010, AJ, 140, 416 [NASA ADS] [CrossRef] [Google Scholar]
  11. Bonnarel, F., Fernique, P., Bienaymé, O., et al. 2000, A&AS, 143, 33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Brott, I., de Mink, S. E., Cantiello, M., et al. 2011, A&A, 530, A115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Choi, J., Dotter, A., Conroy, C., et al. 2016, ApJ, 823, 102 [Google Scholar]
  14. Corsaro, E., Lee, Y.-N., García, R. A., et al. 2017, Nat. Astron., 1, 0064 [NASA ADS] [CrossRef] [Google Scholar]
  15. Croxall, K. V., van Zee, L., Lee, H., et al. 2009, ApJ, 705, 723 [NASA ADS] [CrossRef] [Google Scholar]
  16. Dalcanton, J. J., Williams, B. F., Seth, A. C., et al. 2009, ApJS, 183, 67 [Google Scholar]
  17. Dallas, M. M., Oey, M. S., & Castro, N. 2022, ApJ, 936, 112 [NASA ADS] [CrossRef] [Google Scholar]
  18. Dayal, P., & Ferrara, A. 2018, Phys. Rep., 780, 1 [Google Scholar]
  19. de Mink, S. E., & Mandel, I. 2016, MNRAS, 460, 3545 [NASA ADS] [CrossRef] [Google Scholar]
  20. Dohm-Palmer, R. C., Skillman, E. D., Mateo, M., et al. 2002, AJ, 123, 813 [NASA ADS] [CrossRef] [Google Scholar]
  21. Dotter, A. 2016, ApJS, 222, 8 [Google Scholar]
  22. Dufton, P. L., Langer, N., Dunstall, P. R., et al. 2013, A&A, 550, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Dufton, P. L., Evans, C. J., Hunter, I., Lennon, D. J., & Schneider, F. R. N. 2019, A&A, 626, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Ekström, S., Meynet, G., Maeder, A., & Barblan, F. 2008, A&A, 478, 467 [Google Scholar]
  25. Erb, D. K., Pettini, M., Shapley, A. E., et al. 2010, ApJ, 719, 1168 [Google Scholar]
  26. Espinosa Lara, F., & Rieutord, M. 2011, A&A, 533, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Evans, C. J., Howarth, I. D., Irwin, M. J., Burnley, A. W., & Harries, T. J. 2004, MNRAS, 353, 601 [NASA ADS] [CrossRef] [Google Scholar]
  28. Evans, C. J., Lennon, D. J., Smartt, S. J., & Trundle, C. 2006, A&A, 456, 623 [CrossRef] [EDP Sciences] [Google Scholar]
  29. Evans, C. J., Castro, N., Gonzalez, O. A., et al. 2019, A&A, 622, A129 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Frémat, Y., Zorec, J., Hubert, A. M., & Floquet, M. 2005, A&A, 440, 305 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Gaia Collaboration (Helmi, A., et al.) 2018, A&A, 616, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Gaia Collaboration (Brown, A. G. A., et al.) 2021, A&A, 649, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Garcia, M., Evans, C. J., Bestenlehner, J. M., et al. 2021, Exp. Astron., 51, 887 [NASA ADS] [CrossRef] [Google Scholar]
  34. Garcia, M., Herrero, A., Najarro, F., et al. 2017, in The Lives and Death-Throes of Massive Stars, eds. J. J. Eldridge, J. C. Bray, L. A. S. McClelland, & L. Xiao, 329, 313 [Google Scholar]
  35. Gehan, C., Mosser, B., Michel, E., & Cunha, M. S. 2021, A&A, 645, A124 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Gordon, K. D., Clayton, G. C., Misselt, K. A., Landolt, A. U., & Wolff, M. J. 2003, ApJ, 594, 279 [NASA ADS] [CrossRef] [Google Scholar]
  37. Hainich, R., Pasemann, D., Todt, H., et al. 2015, A&A, 581, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Hastings, B., Wang, C., & Langer, N. 2020, A&A, 633, A165 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Hastings, B., Langer, N., Wang, C., Schootemeijer, A., & Milone, A. P. 2021, A&A, 653, A144 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Hilditch, R. W., Howarth, I. D., & Harries, T. J. 2005, MNRAS, 357, 304 [NASA ADS] [CrossRef] [Google Scholar]
  41. Holtzman, J. A., Afonso, C., & Dolphin, A. 2006, ApJS, 166, 534 [NASA ADS] [CrossRef] [Google Scholar]
  42. Hummel, W., Szeifert, T., Gässler, W., et al. 1999, A&A, 352, L31 [NASA ADS] [Google Scholar]
  43. Hunter, I., Lennon, D. J., Dufton, P. L., et al. 2008, A&A, 479, 541 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Iqbal, S., & Keller, S. C. 2013, MNRAS, 435, 3103 [NASA ADS] [CrossRef] [Google Scholar]
  45. Jackson, R. J., & Jeffries, R. D. 2010, MNRAS, 402, 1380 [NASA ADS] [CrossRef] [Google Scholar]
  46. Kaaret, P., Feng, H., & Roberts, T. P. 2017, ARA&A, 55, 303 [Google Scholar]
  47. Karachentsev, I. D., Dolphin, A. E., Geisler, D., et al. 2002, A&A, 383, 125 [CrossRef] [EDP Sciences] [Google Scholar]
  48. Kaufer, A., Venn, K. A., Tolstoy, E., Pinte, C., & Kudritzki, R.-P. 2004, AJ, 127, 2723 [NASA ADS] [CrossRef] [Google Scholar]
  49. Keller, S. C., Wood, P. R., & Bessell, M. S. 1999, A&AS, 134, 489 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Klement, R., Carciofi, A. C., Rivinius, T., et al. 2019, ApJ, 885, 147 [Google Scholar]
  51. Kniazev, A. Y., Grebel, E. K., Pustilnik, S. A., Pramskij, A. G., & Zucker, D. B. 2005, AJ, 130, 1558 [NASA ADS] [CrossRef] [Google Scholar]
  52. Korn, A. J., Becker, S. R., Gummersbach, C. A., & Wolf, B. 2000, A&A, 353, 655 [NASA ADS] [Google Scholar]
  53. Kruckow, M. U., Tauris, T. M., Langer, N., Kramer, M., & Izzard, R. G. 2018, MNRAS, 481, 1908 [CrossRef] [Google Scholar]
  54. Langer, N. 1992, A&A, 265, L17 [NASA ADS] [Google Scholar]
  55. Lunnan, R., Chornock, R., Berger, E., et al. 2014, ApJ, 787, 138 [NASA ADS] [CrossRef] [Google Scholar]
  56. Macri, L. M., Stanek, K. Z., Bersier, D., Greenhill, L. J., & Reid, M. J. 2006, ApJ, 652, 1133 [NASA ADS] [CrossRef] [Google Scholar]
  57. Maeder, A. 1987, A&A, 178, 159 [NASA ADS] [Google Scholar]
  58. Maeder, A., Grebel, E. K., & Mermilliod, J.-C. 1999, A&A, 346, 459 [NASA ADS] [Google Scholar]
  59. Marchant, P., Langer, N., Podsiadlowski, P., Tauris, T. M., & Moriya, T. J. 2016, A&A, 588, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  60. Marchant, P., Langer, N., Podsiadlowski, P., et al. 2017, A&A, 604, A55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Martayan, C., Frémat, Y., Hubert, A. M., et al. 2006, A&A, 452, 273 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  62. Martayan, C., Frémat, Y., Hubert, A. M., et al. 2007, A&A, 462, 683 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Martayan, C., Baade, D., & Fabregat, J. 2010, A&A, 509, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Mateo, M. L. 1998, ARA&A, 36, 435 [NASA ADS] [CrossRef] [Google Scholar]
  65. McSwain, M. V., & Gies, D. R. 2005, ApJS, 161, 118 [Google Scholar]
  66. Meixner, M., Gordon, K. D., Indebetouw, R., et al. 2006, AJ, 132, 2268 [NASA ADS] [CrossRef] [Google Scholar]
  67. Milone, A. P., Marino, A. F., Di Criscienzo, M., et al. 2018, MNRAS, 477, 2640 [Google Scholar]
  68. Mokiem, M. R., de Koter, A., Evans, C. J., et al. 2006, A&A, 456, 1131 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  69. Mokiem, M. R., de Koter, A., Vink, J. S., et al. 2007, A&A, 473, 603 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Mosser, B., Gehan, C., Belkacem, K., et al. 2018, A&A, 618, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  71. Neugent, K. F., Massey, P., & Morrell, N. 2018, ApJ, 863, 181 [NASA ADS] [CrossRef] [Google Scholar]
  72. Okazaki, A. T. 2001, PASJ, 53, 119 [Google Scholar]
  73. Peters, M., Wisniewski, J. P., Williams, B. F., et al. 2020, AJ, 159, 119 [NASA ADS] [CrossRef] [Google Scholar]
  74. Pineau, F. X., Boch, T., Derrière, S., & Schaaff, A. 2020, in Astronomical Data Analysis Software and Systems XXVII, eds. P. Ballester, J. Ibsen, M. Solar, & K. Shortridge, ASP Conf. Ser., 522, 125 [NASA ADS] [Google Scholar]
  75. Ramachandran, V., Hamann, W. R., Hainich, R., et al. 2018, A&A, 615, A40 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  76. Ramachandran, V., Hamann, W. R., Oskinova, L. M., et al. 2019, A&A, 625, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  77. Rivinius, T., Carciofi, A. C., & Martayan, C. 2013, A&A Rev., 21, 69 [NASA ADS] [CrossRef] [Google Scholar]
  78. Sabbi, E., Calzetti, D., Ubeda, L., et al. 2018, ApJS, 235, 23 [Google Scholar]
  79. Sabbi, E., Adamo, A., Bianchi, L. C., et al. 2020, GULP: Galaxy UV Legacy Project, HST Proposal, Cycle 28, ID. #16316 [Google Scholar]
  80. Sana, H., de Mink, S. E., de Koter, A., et al. 2012, Science, 337, 444 [Google Scholar]
  81. Savaglio, S., Glazebrook, K., & Le Borgne, D. 2009, ApJ, 691, 182 [Google Scholar]
  82. Saviane, I., Rizzi, L., Held, E. V., Bresolin, F., & Momany, Y. 2002, A&A, 390, 59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  83. Schneider, F. R. N., Ohlmann, S. T., Podsiadlowski, P., et al. 2019, Nature, 574, 211 [Google Scholar]
  84. Schootemeijer, A., Götberg, Y., Mink, S. E. D., Gies, D., & Zapartas, E. 2018, A&A, 615, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  85. Schootemeijer, A., Langer, N., Grin, N. J., & Wang, C. 2019, A&A, 625, A132 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  86. Schootemeijer, A., Langer, N., Lennon, D., et al. 2021, A&A, 646, A106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  87. Shenar, T., Gilkis, A., Vink, J. S., Sana, H., & Sander, A. A. C. 2020, A&A, 634, A79 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  88. Struve, O. 1931, ApJ, 73, 94 [Google Scholar]
  89. Szécsi, D., Langer, N., Yoon, S.-C., et al. 2015, A&A, 581, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  90. Tammann, G. A., Reindl, B., & Sandage, A. 2011, A&A, 531, A134 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  91. Taylor, M. B. 2005, in Astronomical Data Analysis Software and Systems XIV, eds. P. Shopbell, M. Britton, & R. Ebert, ASP Conf. Ser., 347, 29 [Google Scholar]
  92. Telford, O. G., Chisholm, J., McQuinn, K. B. W., & Berg, D. A. 2021, ApJ, 922, 191 [NASA ADS] [CrossRef] [Google Scholar]
  93. Townsend, R. H. D., Owocki, S. P., & Howarth, I. D. 2004, MNRAS, 350, 189 [Google Scholar]
  94. Trundle, C., Dufton, P. L., Hunter, I., et al. 2007, A&A, 471, 625 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  95. van Bever, J., & Vanbeveren, D. 1997, A&A, 322, 116 [NASA ADS] [Google Scholar]
  96. van Zee, L., Skillman, E. D., & Haynes, M. P. 2006, ApJ, 637, 269 [Google Scholar]
  97. Vink, J. S., & Sander, A. A. C. 2021, MNRAS, 504, 2051 [NASA ADS] [CrossRef] [Google Scholar]
  98. von Zeipel, H. 1924, MNRAS, 84, 665 [NASA ADS] [CrossRef] [Google Scholar]
  99. Wang, C., Langer, N., Schootemeijer, A., et al. 2020, ApJ, 888, L12 [NASA ADS] [CrossRef] [Google Scholar]
  100. Wisniewski, J. P., & Bjorkman, K. S. 2006, ApJ, 652, 458 [NASA ADS] [CrossRef] [Google Scholar]
  101. Yang, M., Bonanos, A. Z., Jiang, B.-W., et al. 2019, A&A, 629, A91 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  102. Zaritsky, D., Harris, J., Thompson, I. B., Grebel, E. K., & Massey, P. 2002, AJ, 123, 855 [Google Scholar]
  103. Zaritsky, D., Harris, J., Thompson, I. B., & Grebel, E. K. 2004, AJ, 128, 1606 [Google Scholar]
  104. Zorec, J., & Briot, D. 1997, A&A, 318, 443 [NASA ADS] [Google Scholar]

Appendix A: Tests

A.1. Identification of OBe stars

Magellanic Clouds To divide the sources in our CMDs of the SMC (Fig. 1) and the LMC (Fig. A.1) into HeB, MS, and OBe stars, we use color cuts in color–color diagrams (as described in the main text for the SMC). Here we provide a few more details and show how we repeated the method for the LMC.

thumbnail Fig. A.1.

Same as Fig. 1, but for the LMC. The adopted value for the distance modulus is 18.41 (Macri et al. 2006).

We base the location of the color cuts in the color–color diagrams (Figs. 2, A.2) on the distribution of stars with known spectral types. The studies that we use are from a sample of spectroscopic SMC studies (Martayan et al. 2007; Hunter et al. 2008; Dufton et al. 2019; Bonanos et al. 2010, hereafter referred to as SSSS) and a sample of spectroscopic LMC studies (Evans et al. 2006; Martayan et al. 2006; Ramachandran et al. 2018; Bonanos et al. 2009, hereafter referred to as SSLS). To be able investigate the colors and magnitudes of the sources presented in these studies, we cross-correlated them with the Gaia EDR3, UBVI, and Spitzer-SAGE catalogs described above, again adopting a 1" cross-correlation radius.

thumbnail Fig. A.2.

Same as Fig. 2, but for the LMC. The sources with known spectral types are from Martayan et al. (2006), Evans et al. (2006), Ramachandran et al. (2018). The border between MS and OBe stars lies at the same J − [3.6] color as in the SMC, while the border between He-burning stars and MS/OBe stars lies at U − B = −0.3 (for the SMC we adopt U − B = −0.5).

From the Bonanos et al. catalogs, hereafter referred to as B09 (Bonanos et al. 2009) for the LMC catalog and B10 (Bonanos et al. 2010) for the SMC catalog, in Figs.2 and A.2 we show only the stars with spectral type B6 and later (i.e., where we do not make a distinction between stars with and without emission features in their spectrum). We chose this approach because the B09 and B10 catalogs are based on observations that do not provide as accurate a means to distinguish O/early B from OBe stars as the rest of the SSLS and SSSS data, as we show below . The stars that are classified as OBe stars in B09 and B10 are almost always stars for which Hα spectroscopy exists, which is the main diagnostic for the presence of a decretion disk. However, most of the sources in B09 and B10 do not have Hα coverage. To illustrate: about 4000 out of 5000 sources in the B10 catalog are from the 2dF survey (Evans et al. 2004), and from these 4000 sources about 1000 have been observed in Hα. A possible OBe nature would most likely not be recognized for the other 3000 sources.

For the location of the color cut between early and late B stars in the SMC, we choose U − B = −0.5. We do so because above this value most stars with known spectral types are late B stars, while below this value the majority are early B stars (Fig. 2). We notice a systematic shift in U − B color of about 0.2 magnitudes between the SMC and the LMC (see U − B histograms on the left half of Figs. 2 and A.2). Therefore, we choose U − B = −0.3 as the color cut between early and late B-type stars in the LMC.

For the color cut between MS and OBe stars in the SMC, we use J − [3.6]=0.25 because at larger values, the number of OBe stars starts to surpass the number of nonOBe stars (right panel in Fig. 2). We notice that unlike for the U − B color, there is no systematic shift in J − [3.6] color distributions between the SMC and the LMC. This could be explained by the small amount of extinction that takes place at these long wavelengths (Gordon et al. 2003) and the lower sensitivity of the J − [3.6] color (compared to U − B) to temperatures of hot stars. Therefore, we adopt J − [3.6]=0.25 as the color cut for OBe stars, as we do for the SMC. Also, we note that the resemblance of OBe star J − [3.6] colors in the LMC and SMC could be an indication that the disk properties do not dramatically change between these two environments.

Sextans A The CMD constructed with the Sextans A data set is shown in Fig. A.3 and the color–color diagram in in Fig. A.4. In this diagram, a distinct group of stars resides to the right of the MS stars, as was the case in the SMC and LMC color–color diagrams (Figs. 2 and A.2). As the disk causes a smaller excess for the F555W − F814W color than for the J − [3.6] color, the OBe stars and MS stars are relatively close to each other in Fig. A.4. To test our assumption that the three groups in Fig. A.4 correspond to MS, OBe, and HeB stars, we provide three different tests.

thumbnail Fig. A.3.

Same as Fig. 1, but for HST data (Holtzman et al. 2006; Bianchi et al. 2012) for Sextans A. For this galaxy, we identify MS, OBe, and HeB stars using the color–color diagram Fig. A.4. We adopted a distance modulus of DM = 25.63 (Tammann et al. 2011).

As a first test, in Fig. A.5 we make a similar plot for SSSS data in the SMC (the colors shown on both axes are not exactly the same as in Fig. A.4, but the mean wavelengths of the relevant filters are very similar5. This plot shows that, indeed, in the SMC the MS, OBe, and HeB stars are well separated in this diagram, although they are not as segregated as in Fig. 2. We elect to adopt slightly conservative color cuts for OBe stars in Fig. A.4. In Fig. A.5, we see that in the SMC the best location for the color cut lies where the number distribution in the top panels reaches about one-quarter of its maximum value Nmax. In Fig. A.4, we adopt N/Nmax = 1/7.

thumbnail Fig. A.4.

Same as the left side of Fig. A.5, except that the colors are constructed using HST data (Holtzman et al. 2006; Bianchi et al. 2012) for Sextans A. Only sources brighter than an absolute magnitude of MF555W = −3 are shown, for an adopted distance modulus of DM = 25.63 (Tammann et al. 2011).

thumbnail Fig. A.5.

Same as Fig. 2, but on the x-axis is the Gbp − Grp color instead of J − [3.6], and only showing sources brighter than MGbp = −3. Also, as we employ a diagonal color cut (gray line) on the x-axis of the top panel we show the distance to the color cut.

As a second test of the effectiveness of the color cuts in Fig. A.4, we investigate a subset of sources in Sextans A that have been observed with the Hα narrow-band filter F656N on board HST. Similarly to what we did previously, we cross-correlate F656N data from Dohm-Palmer et al. (2002) with our B12 x H06 data product in a 0.35" radius. If we cross-match two sources that have a difference larger than 0.4 magnitudes in F555W between the values from Dohm-Palmer et al. (2002) and H06, we discard the cross-match. We present a color–color diagram that includes a color constructed with F656N filter data (Fig. A.6). In the left panel, we draw a box for sources that have a blue U − B color, and a red F555W − F656N color that indicates that they have Hα in emission. We expect OBe stars to reside in this box. The right panel of Fig. A.6 shows that almost all photometric OBe stars appear to have Hα in emission, as would be expected.

thumbnail Fig. A.6.

Color–color diagrams of Sextans A that are used to test how well we can identify OBe stars there. The catalogs that are used to build photometric colors are provided in parentheses are B12 (Bianchi et al. 2012), D02 (Dohm-Palmer et al. 2002), and H06 (Holtzman et al. 2006). Left: Color–color diagram where the color on the x-axis is constructed using the magnitude in the Hα narrow-band filter F656N. Sources that are in the box on the bottom right have colors that indicate that they are hot stars with Hα emission (i.e., OBe stars). Right: Same axes and color cuts as in Fig. A.4, but we show sources that are in the box for hot sources with Hα emission (left panel) in orange here. The presence of Hα excess in most photometric OBe sources supports their inferred OBe nature.

One might wonder if the use of an optical/near-IR color for Sextans A (Fig. A.4) instead of the J − [3.6] color (e.g., Fig. A.4, for the SMC) to identify OBe stars might systematically affect the derived OBe star fraction. As a third test, we investigate this in the SMC, where both optical/near-IR and J − [3.6] data are available. Figure A.5 is similar to Fig. 2, except that the color on the x-axis is based on optical/near-IR filters. We test the impact that the choice of color–color diagram for identifying OBe stars (Fig. 2 versus Fig. A.5) has on the OBe star fraction that we derive. As we show in Fig. A.7, we find that the two methods result in highly similar OBe fractions for sources dimmer than MGbp = −5. In the brighter-source regime at MGbp < −5, the OBe star fraction is higher according to the method that employs the Gbp − Grp color instead of J − [3.6]. This could be the result of, for example, the very brightest stars being affected by extinction more than average (see Fig. 6 of Schootemeijer et al. 2021), causing them to be unjustly identified as OBe stars. Therefore, we caution that the OBe star fraction in Sextans A in the −5 > Mabs > −5.5 bin (which excluded in Fig. 4) might in reality be lower than what is shown in Fig. 3.

thumbnail Fig. A.7.

OBe star fraction in the SMC as a function of absolute magnitude. The red line shows the result obtained using a color–color cut with the J − [3.6] color on the x-axis (shown in Fig. 2). As a test, with a cyan line we also show the result obtained when the Gbp − Grp color is used on the x-axis instead (shown in Fig. A.5). In the grey shaded area, we do not present an OBe star fraction.

Holmberg I and Holmberg II For Holmberg I, we show the CMD in Fig. A.8 and the color–color diagram in Fig. A.9. For Holmberg II, we show the CMD in Fig. A.10 and the color–color diagram in Fig. A.11. In both Fig. A.9 and A.11, we adopt a magnitude cut of MV = −3. Therefore, we list no number or fraction of OBe stars for these dimmer sources in Table 1. As for Sextans A, for Ho I and Ho II the color cuts are again placed slightly conservatively, at the point where N/Nmax reaches down to the value of 1/7 in the number distributions shown above the color–color diagrams.

thumbnail Fig. A.8.

Same as Fig. A.3, but for HST data (Sabbi et al. 2018) of Holmberg I. The adopted distance modulus is 27.92 (Sabbi et al. 2018).

thumbnail Fig. A.9.

Same as the left panel of Fig. A.5, except that the colors are constructed using HST data (Sabbi et al. 2018) for Holmberg I.

thumbnail Fig. A.10.

Same as Fig. A.3, but for HST data (Sabbi et al. 2018) for Holmberg II. The adopted distance modulus is 27.55 (Sabbi et al. 2018).

thumbnail Fig. A.11.

Same as the left side of Fig. A.5, except that the colors are constructed using HST data (Sabbi et al. 2018) for Holmberg II.

A.2. Completeness

Magellanic Clouds We estimate the completeness of the cross-matched data products by counting the number of sources they contain in different magnitude intervals, and dividing that number by the number of sources in the Gaia EDR3 catalog (after cleaning the Gaia EDR3 catalog of foreground stars in the same fashion). The Gaia EDR3 catalog should be essentially complete (Arenou et al. 2018) in our magnitude range of interest (12 ≲ mGbp ≲ 17). We only consider a part of the sky that is fully covered by the cross-matched data product (e.g., the UB data do no cover the Wing of the SMC).

With Gaia photometry alone, we cannot distinguish MS, OBe, and HeB stars. Therefore, we estimate the completeness in the "MS band" and the "OBe/HeB band". In the SMC, we define the MS band as the sources with Gbp − Grp < −0.15 and for the OBe/HeB band we adopt −0.15 < Gbp − Grp < 0.2. For the LMC, we place the cuts at Gbp − Grp < −0.1 for the MS band and −0.1 < Gbp − Grp < 0.25 for the OBe/HeB band. The result is shown in Fig. A.12. We find that the completeness is about 80%−90% in the magnitude range −3 > MGbp > −5.5 (i.e., the range in which we plot the OBe star fraction in Fig. 3). For the dimmest sources, the completeness drops dramatically, especially for sources in the MS bands. We interpret this as the consequence of the hot MS sources becoming too dim to be observed in the [3.6] band. This drop is deeper for the SMC because its distance modulus is half a magnitude larger (Macri et al. 2006; Hilditch et al. 2005).

thumbnail Fig. A.12.

Completeness fraction estimate for MS sources and sources in the band that contains OBe stars and He-burning sources (OBe/HeB-band; see main text) shown for both the SMC and the LMC.

Sextans A Similar to the SMC and LMC photometric data sets, we estimate the completeness of the cross-matched Sextans A HST data. Here, we adopt F555W − F814W < −0.05 for the MS band and −0.05 < F555W − F814W < 0.2 for the OBe/HeB band. Now, we divide the number of sources in the cross-correlated H06 x B12 catalog by the number of sources with a F555W as well as F814W magnitude in H06. As such, this completeness estimate assumes that the H06 F555W data are complete, which might not be the case. However, we see no clear reason for why the H06 F555W data (that go two magnitudes deeper than the the dimmest sources we consider at MF555W = −2, the limit in our CMD) would be significantly more (or less) complete for MS stars than for OBe stars, which would bias the values we derive for the OBe star fraction. Figure A.13 shows that the cross-correlated data set is about 80% complete compared to the H06 catalog. In the absolute magnitude range −2 > MF555W > −3 (where we do not present an OBe star fraction), fewer sources are retained after the cross-correlation.

thumbnail Fig. A.13.

Completeness fraction estimate for MS sources and sources in the band that contains OBe and He-burning stars (OBe/HeB-band; see main text) shown for Sextans A.

Holmberg I and Holmberg II For Holmberg I and Holmberg II, we adopt a similar approach as described above for Sextans A. To estimate the completeness, here we divide the number of sources that have UBVI data (i.e., in all four filters) by the number of sources that have at least VI data, in different MV intervals. We adopt V − I < −0.1 for the MS band and −0.1 < V − I < 0.25 for the OBe/HeB band in both Holmberg I an II. The resulting magnitude-dependent completeness fractions are shown in Fig. A.14. We find that in the range −3 > MV − 5.5, where we present OBe star fractions, the estimated completeness is about 70%. Most importantly, we obtain very similar values for the completeness in the MS band and the OBe/HeB band, implying that completeness issues do not bias our derived values for the OBe star fraction.

thumbnail Fig. A.14.

Completeness fraction estimate for MS sources and sources in the band that contains OBe and He-burning stars (OBe/HeB-band; see main text) shown for Holmberg II (left) and Holmberg I (right).

A.3. Pollution by Wolf-Rayet and B[e] stars

We use the B09 and B10 catalogs to investigate potential contamination of our OBe stars by Wolf-Rayet (WR) and B[e] stars in the Magellanic Clouds. The LMC B09 catalog contains 89 WR stars in the magnitude range −2 > MGbp > −7. Out of these, 69 (i.e., the vast majority) have colors that would lead to them being labeled OBe stars (see Fig. A.2). The B[e] stars in the magnitude range −2 > MGbp > −7 are fewer in number (eight), and seven of these fall within the OBe color cut (typically they are far redder in the IR than OBe stars, with J − [3.6]> 2). In the top right panel of Fig. A.15, we show the logarithm of the number of photometric OBe stars and the polluting WR and B[e] stars in the LMC. In the bottom right panel, we display the fraction of photometric OBe sources in the LMC that are known to have a WR or B[e] spectral type. This shows that for sources brighter than MGbp = −5, where the photometric OBe stars become few in number, polluting WR and B[e] stars start to make up a significant fraction. At least 20% of the photometric OBe stars in the LMC brighter than MGbp = −5 are known to not be OBe stars. However, for dimmer sources, pollution by WR and B[e] stars plays no significant role. We also note that the current census of WR stars in the LMC and SMC should be close to complete (Neugent et al. 2018). Therefore, we do not expect the WR pollution of photometric OBe stars to be significantly stronger than what is shown in the bottom panels of Fig. A.15.

thumbnail Fig. A.15.

Contamination of photometric OBe stars (‘phot OBe’; as identified by the color cuts in Figs. 2 and A.2). Top panels: Number distributions of absolute magnitude MGbp for various types of sources. The green line shows photometric OBe stars. The orange and blue lines show photometric OBe stars that are known to instead be B[e] and Wolf-Rayet (WR) stars, respectively, and that are in the SMC and LMC spectral type catalogs of Bonanos et al. (2009, 2010). Bottom panels: Fraction of all photometric OBe stars that are known to be WR stars or B[e] stars instead, in different absolute magnitude intervals.

In the SMC (left panels of Fig. A.15), there are relatively fewer WR stars than in the LMC. One of the reasons for this is the weaker winds, which mean that stars need to be more luminous to manifest themselves as WR stars (Shenar et al. 2020). There are also three B[e] stars in the B10 catalog that have OBe-like colors in our magnitude range of interest. For the SMC, it is only in the −6 < MGbp < −5.5 bin that there is significant pollution by WR and B[e] stars (20%).

From the above, we conclude that WR star sources more often than not manifest themselves as photometric OBe stars, but that they are far fewer in number than genuine OBe stars. The exception is that for sources brighter than MGbp = −5.5, the WR plus B[e] star pollution can be at least 20% in both the SMC and the LMC. For this reason, we refrain from presenting an OBe star fraction for sources brighter than Mabs = −5.5.

A.4. Assessing uncertainties in calculating the OBe star fraction

To calculate the OBe star fraction in various absolute magnitude intervals, we divide the number of OBe stars (NOBe) by the total number of H-burning stars, which is the sum of the number of OBe stars plus the number of nonOBe MS stars (NMS): fOBe = NOBe/(NOBe + NMS). We then calculate the statistical error as: σ f OBe = f OBe · 1 f OBe N OBe + N MS $ \sigma_{f_{\mathrm{OBe}}} = \sqrt{ f_{\mathrm{OBe}} \cdot \frac{ 1-f_{\mathrm{OBe}} }{N_{\mathrm{OBe}} + N_{\mathrm{MS}}}} $. This is a lower limit on the true error, as other uncertainties (e.g., the exact location of the MS-OBe border in our color–color diagrams) are not quantifiable, at least not in all our target dwarf galaxies.

In theory, the disk of an OBe star could contribute to the Gbp, V, or F555W flux, biasing our comparison. We do not correct for a disk contribution; below, we justify this choice. In Fig. A.16, we show the distribution of various colors for early B and Be sources in the B10 catalog. This figure shows that these B and Be stars appear very similar in terms of the colors that are constructed with the bluer filters (U − B, U − Gbp, B − Gbp). We interpret this as being the result of an at most small disk contribution in the U, B, and Gbp filters. This makes a correction for the disk contribution to Gbp unnecessary. It is only in the Gbp − Grp and IR colors that OBe stars start to look relatively red. We note that to measure the OBe star fractions we chose relatively blue filters (i.e., Gbp, V, and F555W instead of G/Grp, I, and F814W) on purpose, because there we expected the disk contribution to play a less significant role.

thumbnail Fig. A.16.

Histograms of the different colors of sources that have a B spectral type (Bonanos et al. 2010) (including OBe). Only early B stars (B0 to B3) with an absolute Gbp magnitude of −3 > MGbp > −5 are included. The dash-dotted lines show the median color. The mean wavelength λmean of each filter (top right panel) is taken from http://svo2.cab.inta-csic.es/theory/fps/.

An uncertainty for the OBe fraction could be an apparent lack of young stars that are observable. In the SMC, it was found that massive stars in the first half of their lifetime seem to be making up far less than half of the population (Schootemeijer et al. 2021), possibly because they are hidden by their birth cloud, or perhaps as a result of a declining star-formation rate. Stars at low Z in the second half of their lifetime are more often predicted to be OBe stars —as a result of core contraction in the single star channel (Ekström et al. 2008; Hastings et al. 2020), and binary interaction typically not taking place during the early MS phase (Wang et al. 2020). Therefore, a dearth of stars in the first half of their MS lifetime could cause us to overestimate the OBe star fraction in the SMC, at least compared to a scenario where star formation is constant and stars are observable throughout their lifetime. However, we do not expect the star formation histories of the other dwarf galaxies to follow that of the SMC, indicating that the high OBe star fraction of ∼0.3 is not an artifact caused by star formation history.

Appendix B: Magnitude-mass relations for blue stars

To estimate stellar masses for OBe and nonOBe MS stars from their absolute magnitudes, we look up the evolutionary masses that have been determined for these sources in the SMC by Schootemeijer et al. (2021), to which we refer to below as S21. From S21, we consider only stars that have spectral types B5 and earlier, and we include both nonOBe and OBe stars. In the left panel of Fig. B.1 we plot their evolutionary mass against absolute magnitude, and on the right we show the evolutionary mass distribution in different absolute magnitude intervals. We note that there might be other biases resulting from selection effects in the B10 catalog. However, the B10 catalog is expected to be rather complete —about 80% at the bright end, and roughly 50% for dimmer sources, according to S21. Another possible source of error in the magnitude–mass relations comes from extrapolating them from the SMC to other galaxies and other filters (F555W and V instead of Gbp). We note that a difference of ΔV = 0.5 mag has been inferred between equal-mass massive stars in the Milky Way and at Z/50 (Evans et al. 2019), implying offsets of the order of a tenth of a magnitude in the metallicity range of 0.1 ≤ Z/Z < 0.5 considered in this work. Because of this, in combination with the sizeable scatter of evolutionary masses in each MGbp interval, we caution that the interpretation of our Mabs –evolutionary mass relations for blue stars should not be that each of the MS and OBe sources has exactly the mass given by the relation; instead the masses shown at the top of Fig. 3 are intended to serve as an indication.

thumbnail Fig. B.1.

Left: Diagram showing the evolutionary mass Mevo of stars with spectral types B5 and earlier plotted against absolute magnitude MGbp, for both stars with (OBe) and without (MS) emission features in their spectra. The spectral type catalog is from Bonanos et al. (2010) and the evolutionary masses were determined by Schootemeijer et al. (2021). Right: Histograms that show the evolutionary mass distributions of the sources that are displayed in the left panel, in different MGbp intervals. The average evolutionary masses of sources in these magnitude intervals (i.e., OBe and MS) are written in each panel.

Appendix C: Simulations: comparison with theory in the SMC

C.1. Single stars

We first theoretically predict the OBe star fraction in the SMC under the assumption that all OBe stars are formed through the single-star channel, that is, without the aid of a binary companion. For this, we repeat the star cluster study of Hastings et al. (2020), but now assuming constant star formation instead of a coeval scenario and extending the mass range of the models to 60 M. This upper mass limit sufficiently high because more massive stellar models are brighter than MGbp = −5 their entire lifetime. This analysis adopts an initial rotation distribution based on observations (Dufton et al. 2013). We note that Dufton et al. (2013) assume random spin axis orientation. Alignment of spin axes (as found by Corsaro et al. 2017, in old open clusters) could affect the distribution of Dufton et al. (2013), although other work does support random spin axis orientation (Hummel et al. 1999; Jackson & Jeffries 2010; Mosser et al. 2018; Gehan et al. 2021).

To obtain the absolute Gbp magnitude, we use the effective temperature and luminosity of a model star, and bolometric corrections from MIST6 (Dotter 2016; Choi et al. 2016). The result is shown in Fig. 5.

C.2. Binary stars

We also theoretically predict the OBe star fraction under the assumption that only binary systems produce OBe stars. We simulate a population of binary systems with the rapid binary evolution code ComBinE (Kruckow et al. 2018), and briefly highlight the most important physics assumptions here. We again adopt a constant star formation rate. The initial rotation velocities are chosen such that the rotation period of both stars is equal to the orbital period, as is the case for tidally synchronized systems. We adopt a Salpeter IMF and an initial mass range of 5 < Mini/M < 60. The initial mass ratio distribution is flat and ranges from qini = 0.1 to qini = 1. In the orbital period distribution, the number of systems scales with the logarithm of the orbital period to the power π = −0.55 (Sana et al. 2012). The initial orbital period range is 0.15 < log(Porb,  ini/d)< 3.5. For the extreme initial mass ratios of q < 0.3, we assume that the stars merge as Roche lobe overflow (RLOF) commences. We also assume that RLOF in systems where the donor has a convective envelope (that goes deeper than 10% of the stellar mass) leads to a merger, and that contact systems merge. This merger product does remain in the stellar population. For the merger product we assume that 10% of the combined stellar mass is lost in the merger event. The starting point of the merger product is a single star model of the appropriate mass (i.e., 90% of the combined progenitor star masses) that has burned the same amount of hydrogen as both progenitors stars have burned pre-merger. In our simulations the merger does not produce an OBe star because it is expected that a high amount of angular momentum is lost post-merger (Schneider et al. 2019). All systems that do not meet one of the merger criteria described above are assumed to go through stable RLOF, resulting in the accretor becoming an OBe star. We assume that the material that is lost from the system during RLOF carries the same specific orbital angular momentum as the accretor. Tidally induced changes of the stellar spins are accounted for. We consider these assumptions to be optimistic, that is, favoring a high OBe star fraction from the binary channel.

When we adopt a binary fraction lower than fbin = 1, we assume that the single stars never become OBe stars, and simply calculate the OBe star fraction as fOBe = fbin ⋅ fOBe (fbin = 1), the latter term being the OBe star fraction calculated for a binary fraction of unity. For binaries, we show the magnitude-dependent theoretical OBe star fractions at the bottom half of Fig. 5. These are calculated under the assumption that stars rotating at v/vcrit ≳ 0.7 manifest themselves as OBe stars. Figure C.1 shows the OBe star fractions that we obtain when v/vcrit ≳ 0.8 and 0.9 are set as minimum rotation rates. A higher minimum rotation rate for OBe stars slightly reduces the number of OBe stars in the theoretical binary population, but not nearly as much as for the single stars.

thumbnail Fig. C.1.

Same as the lower panel of Fig. 5, but here v/vcrit > 0.8 and v/vcrit > 0.9 are set as the required rotation velocity stars to be OBe stars in the binary populations.

All Tables

Table 1.

Number of OBe stars (NOBe) and OBe star fractions (fOBe) in the LMC, SMC, HoI, HoII, and SexA.

All Figures

thumbnail Fig. 1.

Left: color magnitude diagram of stars in the SMC, constructed with Gaia filters. Each black dot represents a source. The brightest sources are at the top and the bluest sources are on the left. We adopt a value of 18.91 for the distance modulus (Hilditch et al. 2005). Right: Stacked histograms of the Gbp − Grp color in the absolute magnitude ranges that are written in the top right of each panel. Main sequence stars are shown in magenta, OBe stars in light orange, and He-burning stars in green. The black line shows the number distribution of all the sources combined.

In the text
thumbnail Fig. 2.

Left: color–color diagram of sources in our SMC data set that are brighter than MGbp = −2. Histograms are shown in the right panel and the top panel, where the top histogram includes only sources bluer than U − B = −0.5. Right: same as the left side, but showing sources from various spectroscopic studies (Martayan et al. 2007; Hunter et al. 2008; Dufton et al. 2019; Bonanos et al. 2010) that are colored for their spectral type.

In the text
thumbnail Fig. 3.

Galaxy-wide OBe star fraction in the Magellanic Clouds, Sextans A, and Holmberg I and II, measured in absolute magnitude ranges −3 > Mabs > −3.5, −3.5 > Mabs > −4, and so on. The tick symbols at the top of the plot indicate the average evolutionary mass that has been inferred (Schootemeijer et al. 2021) for stars of spectral type B5 and earlier in each absolute magnitude bin in the SMC.

In the text
thumbnail Fig. 4.

Metallicity (Z) trends of the OBe star fraction. Top panel: in the low-Z regime, we show the galaxy-wide OBe star fraction calculated in this work in the absolute magnitude interval −3 > Mabs > −5 with large diamonds. Bottom panel: literature studies of stellar subsamples within different galaxies, which targeted cluster and/or field stars: K99 (Keller et al. 1999), M99 (Maeder et al. 1999), M05 (McSwain & Gies 2005), W06 (Wisniewski & Bjorkman 2006), I13 (Iqbal & Keller 2013), Bonanos (Bonanos et al. 2009, 2010), P20 (Peters et al. 2020), and ZB97 (Zorec & Briot 1997). We only show clusters with 20 or more OB+OBe stars, and we do not show multiple measurements for any cluster. In K99, M99, M05, W06, I13, and P20, narrow-band photometry was used to identify OBe stars. For a part of the sample in M99, OBe stars were identified with spectroscopy. Bonanos used IR excess of stars with known spectral types to identify OBe stars. ZB97 used a sample of stars with known spectral types. For reference, we show the results of this work again with a dotted line. The metallicities in Holmberg I and II (Ho I and II) are estimated using oxygen abundances (Croxall et al. 2009). The two highest-metallicity galaxies are our own Milky Way (MW) and Andromeda (M 31). For ZMW, we adopt the Solar metallicity, and we adopt ZM31 = 1.5 Z from P20.

In the text
thumbnail Fig. 5.

Comparison of the observed absolute-magnitude-dependent SMC OBe star fraction with theory. At the top of each panel we show indicative masses, as in Fig. 3. Top: OBe star fraction in a theoretical single star population, for which we explore three different minimum v/vcrit values that stars require to manifest themselves as OBe stars, with v/vcrit being the fraction of the critical rotation velocity. Bottom: OBe star fraction in theoretical binary populations. With blue and orange lines, we show theoretical predictions derived from ComBinE simulations, for different combinations of the mass transfer efficiency β and binary fraction fbin. In the theoretical populations, single stars are assumed to not contribute to producing OBe stars. For each theoretical binary population, we assume that all stars with v/vcrit > 0.7 are OBe stars. In the Appendix, we provide a similar figure, but where v/vcrit > 0.8 and v/vcrit > 0.9 are set as the required rotation velocities for OBe stars in the binary populations (Fig. C.1).

In the text
thumbnail Fig. A.1.

Same as Fig. 1, but for the LMC. The adopted value for the distance modulus is 18.41 (Macri et al. 2006).

In the text
thumbnail Fig. A.2.

Same as Fig. 2, but for the LMC. The sources with known spectral types are from Martayan et al. (2006), Evans et al. (2006), Ramachandran et al. (2018). The border between MS and OBe stars lies at the same J − [3.6] color as in the SMC, while the border between He-burning stars and MS/OBe stars lies at U − B = −0.3 (for the SMC we adopt U − B = −0.5).

In the text
thumbnail Fig. A.3.

Same as Fig. 1, but for HST data (Holtzman et al. 2006; Bianchi et al. 2012) for Sextans A. For this galaxy, we identify MS, OBe, and HeB stars using the color–color diagram Fig. A.4. We adopted a distance modulus of DM = 25.63 (Tammann et al. 2011).

In the text
thumbnail Fig. A.4.

Same as the left side of Fig. A.5, except that the colors are constructed using HST data (Holtzman et al. 2006; Bianchi et al. 2012) for Sextans A. Only sources brighter than an absolute magnitude of MF555W = −3 are shown, for an adopted distance modulus of DM = 25.63 (Tammann et al. 2011).

In the text
thumbnail Fig. A.5.

Same as Fig. 2, but on the x-axis is the Gbp − Grp color instead of J − [3.6], and only showing sources brighter than MGbp = −3. Also, as we employ a diagonal color cut (gray line) on the x-axis of the top panel we show the distance to the color cut.

In the text
thumbnail Fig. A.6.

Color–color diagrams of Sextans A that are used to test how well we can identify OBe stars there. The catalogs that are used to build photometric colors are provided in parentheses are B12 (Bianchi et al. 2012), D02 (Dohm-Palmer et al. 2002), and H06 (Holtzman et al. 2006). Left: Color–color diagram where the color on the x-axis is constructed using the magnitude in the Hα narrow-band filter F656N. Sources that are in the box on the bottom right have colors that indicate that they are hot stars with Hα emission (i.e., OBe stars). Right: Same axes and color cuts as in Fig. A.4, but we show sources that are in the box for hot sources with Hα emission (left panel) in orange here. The presence of Hα excess in most photometric OBe sources supports their inferred OBe nature.

In the text
thumbnail Fig. A.7.

OBe star fraction in the SMC as a function of absolute magnitude. The red line shows the result obtained using a color–color cut with the J − [3.6] color on the x-axis (shown in Fig. 2). As a test, with a cyan line we also show the result obtained when the Gbp − Grp color is used on the x-axis instead (shown in Fig. A.5). In the grey shaded area, we do not present an OBe star fraction.

In the text
thumbnail Fig. A.8.

Same as Fig. A.3, but for HST data (Sabbi et al. 2018) of Holmberg I. The adopted distance modulus is 27.92 (Sabbi et al. 2018).

In the text
thumbnail Fig. A.9.

Same as the left panel of Fig. A.5, except that the colors are constructed using HST data (Sabbi et al. 2018) for Holmberg I.

In the text
thumbnail Fig. A.10.

Same as Fig. A.3, but for HST data (Sabbi et al. 2018) for Holmberg II. The adopted distance modulus is 27.55 (Sabbi et al. 2018).

In the text
thumbnail Fig. A.11.

Same as the left side of Fig. A.5, except that the colors are constructed using HST data (Sabbi et al. 2018) for Holmberg II.

In the text
thumbnail Fig. A.12.

Completeness fraction estimate for MS sources and sources in the band that contains OBe stars and He-burning sources (OBe/HeB-band; see main text) shown for both the SMC and the LMC.

In the text
thumbnail Fig. A.13.

Completeness fraction estimate for MS sources and sources in the band that contains OBe and He-burning stars (OBe/HeB-band; see main text) shown for Sextans A.

In the text
thumbnail Fig. A.14.

Completeness fraction estimate for MS sources and sources in the band that contains OBe and He-burning stars (OBe/HeB-band; see main text) shown for Holmberg II (left) and Holmberg I (right).

In the text
thumbnail Fig. A.15.

Contamination of photometric OBe stars (‘phot OBe’; as identified by the color cuts in Figs. 2 and A.2). Top panels: Number distributions of absolute magnitude MGbp for various types of sources. The green line shows photometric OBe stars. The orange and blue lines show photometric OBe stars that are known to instead be B[e] and Wolf-Rayet (WR) stars, respectively, and that are in the SMC and LMC spectral type catalogs of Bonanos et al. (2009, 2010). Bottom panels: Fraction of all photometric OBe stars that are known to be WR stars or B[e] stars instead, in different absolute magnitude intervals.

In the text
thumbnail Fig. A.16.

Histograms of the different colors of sources that have a B spectral type (Bonanos et al. 2010) (including OBe). Only early B stars (B0 to B3) with an absolute Gbp magnitude of −3 > MGbp > −5 are included. The dash-dotted lines show the median color. The mean wavelength λmean of each filter (top right panel) is taken from http://svo2.cab.inta-csic.es/theory/fps/.

In the text
thumbnail Fig. B.1.

Left: Diagram showing the evolutionary mass Mevo of stars with spectral types B5 and earlier plotted against absolute magnitude MGbp, for both stars with (OBe) and without (MS) emission features in their spectra. The spectral type catalog is from Bonanos et al. (2010) and the evolutionary masses were determined by Schootemeijer et al. (2021). Right: Histograms that show the evolutionary mass distributions of the sources that are displayed in the left panel, in different MGbp intervals. The average evolutionary masses of sources in these magnitude intervals (i.e., OBe and MS) are written in each panel.

In the text
thumbnail Fig. C.1.

Same as the lower panel of Fig. 5, but here v/vcrit > 0.8 and v/vcrit > 0.9 are set as the required rotation velocity stars to be OBe stars in the binary populations.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.