Search

The cosmic origin of the elements, the fundamental chemical building blocks of the universe, is still uncertain. Binary interactions play a key role in the evolution of many massive stars, yet their impact on chemical yields is poorly understood. Using the MESA stellar evolution code, we predict the chemical yields ejected in wind mass loss and the supernovae of single and binary-stripped stars. We do this with a large 162-isotope nuclear network at solar metallicity. We find that binary-stripped stars are more effective producers of the elements than single stars, due to their increased mass loss and an increased chance to eject their envelopes during a supernova. This increased production by binaries varies across the periodic table, with F and K being more significantly produced by binary-stripped stars than single stars. We find that the 12C/13C could be used as an indicator of the conservativeness of mass transfer, as 13C is preferentially ejected during mass transfer while 12C is preferentially ejected during wind mass loss. We identify a number of gamma-ray-emitting radioactive isotopes that may be used to help constrain progenitor and explosion models of core-collapse supernovae with next-generation gamma-ray detectors. For single stars we find that 44V and 52Mn are strong probes of the explosion model, while for binary-stripped stars it is 48Cr. Our findings highlight that binary-stripped stars are not equivalent to two single stars and that detailed stellar modeling is needed to predict their final nucleosynthetic yields. [ABSTRACT FROM AUTHOR]

Gravitational-wave detections are now probing the black hole (BH) mass distribution, including the predicted pair-instability mass gap. These data require robust quantitative predictions, which are challenging to obtain. The most massive BH progenitors experience episodic mass ejections on time-scales shorter than the convective turnover time-scale. This invalidates the steady-state assumption on which the classic mixing length theory relies. We compare the final BH masses computed with two different versions of the stellar evolutionary code |$\tt{MESA}$|⁠ : (i) using the default implementation of Paxton et al. (2018) and (ii) solving an additional equation accounting for the time-scale for convective deceleration. In the second grid, where stronger convection develops during the pulses and carries part of the energy, we find weaker pulses. This leads to lower amounts of mass being ejected and thus higher final BH masses of up to ∼ |$5\, \mathrm{M}_\odot$|⁠. The differences are much smaller for the progenitors that determine the maximum mass of BHs below the gap. This prediction is robust at |$M_{\rm BH, max}\simeq 48\, \mathrm{M}_\odot$|⁠ , at least within the idealized context of this study. This is an encouraging indication that current models are robust enough for comparison with the present-day gravitational-wave detections. However, the large differences between individual models emphasize the importance of improving the treatment of convection in stellar models, especially in the light of the data anticipated from the third generation of gravitational-wave detectors. [ABSTRACT FROM AUTHOR]

We present a detailed study of stellar rotation in the massive 1.5 Gyr old cluster NGC 1846 in the Large Magellanic Cloud. Similar to other clusters at this age, NGC 1846 shows an extended main-sequence turn-off (eMSTO), and previous photometric studies have suggested it could be bimodal. In this study, we use MUSE integral-field spectroscopy to measure the projected rotational velocities (v sin  i) of around |$1400$| stars across the eMSTO and along the upper main sequence of NGC 1846. We measure v sin  i values up to |$\sim 250\, {\rm km\, s^{-1}}$| and find a clear relation between the v sin  i of a star and its location across the eMSTO. Closer inspection of the distribution of rotation rates reveals evidence for a bimodal distribution, with the fast rotators centred around |$v\sin i=140\, {\rm km\, s^{-1}}$| and the slow rotators centred around |$v\sin i=60\, {\rm km\, s^{-1}}$|⁠. We further observe a lack of fast rotating stars along the photometric binary sequence of NGC 1846, confirming results from the field that suggest that tidal interactions in binary systems can spin-down stars. However, we do not detect a significant difference in the binary fractions of the fast and slowly rotating sub-populations. Finally, we report on the serendipitous discovery of a planetary nebula associated with NGC 1846. [ABSTRACT FROM AUTHOR]

Stars stripped of their envelopes from interaction with a binary companion emit a significant fraction of their radiation as ionizing photons. They are potentially important stellar sources of ionizing radiation, however, they are still often neglected in spectral synthesis simulations or simulations of stellar feedback. In anticipating the large datasets of galaxy spectra from the upcoming James Webb Space Telescope, we modeled the radiative contribution from stripped stars by using detailed evolutionary and spectral models. We estimated their impact on the integrated spectra and specifically on the emission rates of H I-, He I-, and He II-ionizing photons from stellar populations. We find that stripped stars have the largest impact on the ionizing spectrum of a population in which star formation halted several Myr ago. In such stellar populations, stripped stars dominate the emission of ionizing photons, mimicking a younger stellar population in which massive stars are still present. Our models also suggest that stripped stars have harder ionizing spectra than massive stars. The additional ionizing radiation, with which stripped stars contribute affects observable properties that are related to the emission of ionizing photons from stellar populations. In co-eval stellar populations, the ionizing radiation from stripped stars increases the ionization parameter and the production efficiency of hydrogen ionizing photons. They also cause high values for these parameters for about ten times longer than what is predicted for massive stars. The effect on properties related to non-ionizing wavelengths is less pronounced, such as on the ultraviolet continuum slope or stellar contribution to emission lines. However, the hard ionizing radiation from stripped stars likely introduces a characteristic ionization structure of the nebula, which leads to the emission of highly ionized elements such as O2+ and C3+. We, therefore, expect that the presence of stripped stars affects the location in the BPT diagram and the diagnostic ratio of O III to O II nebular emission lines. Our models are publicly available through CDS database and on the STARBURST99 website. [ABSTRACT FROM AUTHOR]

How very massive stars form is still an open question in astrophysics. VFTS682 is among the most massive stars known, with an inferred initial mass of ≳ |$150\,\text{M}_\odot$|⁠. It is located in 30 Doradus at a projected distance of 29 pc from the central cluster R136. Its apparent isolation led to two hypotheses: either it formed in relative isolation or it was ejected dynamically from the cluster. We investigate the kinematics of VFTS682 as obtained by Gaia and Hubble Space Telescope astrometry. We derive a projected velocity relative to the cluster of |$38\pm 17{\, \mathrm{km\ s^{-1}}}$| (1σ confidence interval). Although the error bars are substantial, two independent measures suggest that VFTS682 is a runaway ejected from the central cluster. This hypothesis is further supported by a variety of circumstantial clues. The central cluster is known to harbour other stars more massive than |$150\, \text{M}_\odot$| of similar spectral type and recent astrometric studies on VFTS16 and VFTS72 provide direct evidence that the cluster can eject some of its most massive members, in agreement with theoretical predictions. If future data confirm the runaway nature, this would make VFTS682 the most massive runaway star known to date. [ABSTRACT FROM AUTHOR]

Most massive stars, the progenitors of core-collapse supernovae, are in close binary systems and may interact with their companion through mass transfer or merging. We undertake a population synthesis study to compute the delay-time distribution of core-collapse supernovae, that is, the supernova rate versus time following a starburst, taking into account binary interactions. We test the systematic robustness of our results by running various simulations to account for the uncertainties in our standard assumptions. We find that a significant fraction, 15+9-8%, of core-collapse supernovae are "late", that is, they occur 50-200 Myr after birth, when all massive single stars have already exploded. These late events originate predominantly from binary systems with at least one, or, in most cases, with both stars initially being of intermediate mass (4-8 M⊙). The main evolutionary channels that contribute often involve either the merging of the initially more massive primary star with its companion or the engulfment of the remaining core of the primary by the expanding secondary that has accreted mass at an earlier evolutionary stage. Also, the total number of core-collapse supernovae increases by 14+15-14% because of binarity for the same initial stellar mass. The high rate implies that we should have already observed such late core-collapse supernovae, but have not recognized them as such. We argue that φ Persei is a likely progenitor and that eccentric neutron star - white dwarf systems are likely descendants. Late events can help explain the discrepancy in the delay-time distributions derived from supernova remnants in the Magellanic Clouds and extragalactic type Ia events, lowering the contribution of prompt Ia events. We discuss ways to test these predictions and speculate on the implications for supernova feedback in simulations of galaxy evolution. [ABSTRACT FROM AUTHOR]

Context: The VLT-FLAMES Tarantula Survey has observed hundreds of O-type stars in the 30 Doradus region of the Large Magellanic Cloud (LMC). Aims: We study the properties of a statistically significant sample of O-type dwarfs in the same star-forming region and test the latest atmospheric and evolutionary models of the early main-sequence phase of massive stars. Methods: We performed quantitative spectroscopic analysis of 105 apparently single O-type dwarfs. To determine stellar and wind parameters, we used the iacob-gbat package, an automatic procedure based on a large grid of atmospheric models that are calculated with the fastwind code. This package was developed for the analysis of optical spectra of O-type stars. In addition to classical techniques, we applied the Bayesian bonnsai tool to estimate evolutionary masses. Results: We provide a new calibration of effective temperature vs. spectral type for O-type dwarfs in the LMC, based on our homogeneous analysis of the largest sample of such objects to date and including all spectral subtypes. Good agreement with previous results is found, although the sampling at the earliest subtypes could be improved. Rotation rates and helium abundances are studied in an evolutionary context. We find that most of the rapid rotators (v sin i > 300 km s-1) in our sample have masses below ~25 M⊙ and intermediate rotation-corrected gravities (3.9 < log gc < 4.1). Such rapid rotators are scarce at higher gravities (i.e. younger ages) and absent at lower gravities (larger ages). This is not expected from theoretical evolutionary models, and does not appear to be due to a selection bias in our sample. We compare the estimated evolutionary and spectroscopic masses, finding a trend that the former is higher for masses below ~20 M⊙. This can be explained as a consequence of limiting our sample to the O-type stars, and we see no compelling evidence for a systematic mass discrepancy. For most of the stars in the sample we were unable to estimate the wind-strength parameter (hence mass-loss rates) reliably, particularly for objects with lower luminosity (log L/L⊙ ≲ 5.1). Only with ultraviolet spectroscopy will we be able to undertake a detailed investigation of the wind properties of these dwarfs. [ABSTRACT FROM AUTHOR]

Context. Theoretically, rotation-induced chemical mixing in massive stars has far reaching evolutionary consequences, affecting the sequence of morphological phases, lifetimes, nucleosynthesis, and supernova characteristics. Aims. Using a sample of 72 presumably single O-type giants to supergiants observed in the context of the VLT-FLAMES Tarantula Survey (VFTS), we aim to investigate rotational mixing in evolved core-hydrogen burning stars initially more massive than 15M๏ by analysing their surface nitrogen abundances. Methods. Using stellar and wind properties derived in a previous VFTS study we computed synthetic spectra for a set of up to 21 Nii-v lines in the optical spectral range, using the non-LTE atmosphere code FASTWIND.We constrained the nitrogen abundance by fitting the equivalent widths of relatively strong lines that are sensitive to changes in the abundance of this element. Given the quality of the data, we constrained the nitrogen abundance in 38 cases; for 34 stars only upper limits could be derived, which includes almost all stars rotating at 3e sin i > 200 km s-1. Results. We analysed the nitrogen abundance as a function of projected rotation rate 3e sin i and confronted it with predictions of rotational mixing. We found a group of N-enhanced slowly-spinning stars that is not in accordance with predictions of rotational mixing in single stars. Among O-type stars with (rotation-corrected) gravities less than log gc = 3.75 this group constitutes 30-40 percent of the population. We found a correlation between nitrogen and helium abundance which is consistent with expectations, suggesting that, whatever the mechanism that brings N to the surface, it displays CNO-processed material. For the rapidly-spinning O-type stars we can only provide upper limits on the nitrogen abundance, which are not in violation with theoretical expectations. Hence, the data cannot be used to test the physics of rotation induced mixing in the regime of high spin rates. Conclusions. While the surface abundances of 60-70 percent of presumed single O-type giants to supergiants behave in conformity with expectations, at least 30-40 percent of our sample can not be understood in the current framework of rotational mixing for single stars. Even though we have excluded stars showing radial velocity variations, of our sample may have remained contaminated by postinteraction binary products. Hence, it is plausible that effects of binary interaction need to be considered to understand their surface properties. Alternatively, or in conjunction, the effects of magnetic fields or alternative mass-loss recipes may need to be invoked. [ABSTRACT FROM AUTHOR]

We test the predictions of spectral synthesis models based on seven different massive-star prescriptions against Legacy ExtraGalactic UV Survey (LEGUS) observations of eight young massive clusters in two local galaxies, NGC 1566 and NGC 5253, chosen because predictions of all seven models are available at the published galactic metallicities. The high angular resolution, extensive cluster inventory, and full near-ultraviolet to near-infrared photometric coverage make the LEGUS data set excellent for this study. We account for both stellar and nebular emission in the models and try two different prescriptions for attenuation by dust. From Bayesian fits of model libraries to the observations, we find remarkably low dispersion in the median E(B - V) (~0.03 mag), stellar masses (~104M⊙), and ages (~1 Myr) derived for individual clusters using different models, although maximum discrepancies in these quantities can reach 0.09 mag and factors of 2.8 and 2.5, respectively. This is for ranges in median properties of 0.05-0.54 mag, 1.8-10 ? 104M⊙, and 1.6-40 Myr spanned by the clusters in our sample. In terms of best fit, the observations are slightly better reproduced by models with interacting binaries and least well reproduced by models with single rotating stars. Our study provides a first quantitative estimate of the accuracies and uncertainties of the most recent spectral synthesis models of young stellar populations, demonstrates the good progress of models in fitting high-quality observations, and highlights the needs for a larger cluster sample and more extensive tests of the model parameter space. [ABSTRACT FROM AUTHOR]

Context. Oxygen sequence Wolf-Rayet (WO) stars are a very rare stage in the evolution of massive stars. Their spectra show strong emission lines of helium-burning products, in particular highly ionized carbon and oxygen. The properties of WO stars can be used to provide unique constraints on the (post-)helium burning evolution of massive stars, and their remaining lifetimes and the expected properties of their supernovae. Aims. We aim to homogeneously analyze the currently known presumed-single WO stars to obtain the key stellar and outflow properties and to constrain their evolutionary state. Methods. We use the line-blanketed non-local thermal equilibrium atmosphere code CMFGEN to model the X-Shooter spectra of the WO stars and to deduce the atmospheric parameters. We calculate dedicated evolutionary models to determine the evolutionary state of the stars. Results. The WO stars have extremely high temperatures that range from 150 kK to 210 kK, and very low surface helium mass fractions that range from 44% down to 14%. Their properties can be reproduced by evolutionary models with helium zero-age main sequence masses of MHe,ini = 15-25 M0 that exhibit a fairly strong (a few times 10-5 M0 yr-1 ), homogeneous (fc > 0.3) stellar wind. Conclusions. WO stars represent the final evolutionary stage of stars with estimated initial masses of Mini = 40-60 M0. They are post core-helium burning and predicted to explode as type Ic supernovae within a few thousand years. [ABSTRACT FROM AUTHOR]

Context. Model atmosphere analyses have been previously undertaken for both Galactic and extragalactic B-type supergiants. By contrast, little attention has been given to a comparison of the properties of single supergiants and those that are members of multiple systems. Aims. Atmospheric parameters and nitrogen abundances have been estimated for all the B-type supergiants identified in the VLT-FLAMES Tarantula survey. These include both single targets and binary candidates. The results have been analysed to investigate the role of binarity in the evolutionary history of supergiants. Methods. tlusty non-local thermodynamic equilibrium (LTE) model atmosphere calculations have been used to determine atmospheric parameters and nitrogen abundances for 34 single and 18 binary supergiants. Effective temperatures were deduced using the silicon balance technique, complemented by the helium ionisation in the hotter spectra. Surface gravities were estimated using Balmer line profiles and microturbulent velocities deduced using the silicon spectrum. Nitrogen abundances or upper limits were estimated from the N?ii spectrum. The effects of a flux contribution from an unseen secondary were considered for the binary sample. Results. We present the first systematic study of the incidence of binarity for a sample of B-type supergiants across the theoretical terminal age main sequence (TAMS). To account for the distribution of effective temperatures of the B-type supergiants it may be necessary to extend the TAMS to lower temperatures. This is also consistent with the derived distribution of mass discrepancies, projected rotational velocities and nitrogen abundances, provided that stars cooler than this temperature are post-red supergiant objects. For all the supergiants in the Tarantula and in a previous FLAMES survey, the majority have small projected rotational velocities. The distribution peaks at about 50 km?s-1 with 65% in the range 30 km?s-1 ≤ ve sin i ≤ 60 km?s-1. About ten per cent have larger vesini (≥100 km?s-1), but surprisingly these show little or no nitrogen enhancement. All the cooler supergiants have low projected rotational velocities of ≤70 km?s-1 and high nitrogen abundance estimates, implying that either bi-stability braking or evolution on a blue loop may be important. Additionally, there is a lack of cooler binaries, possibly reflecting the small sample sizes. Single-star evolutionary models, which include rotation, can account for all of the nitrogen enhancement in both the single and binary samples. The detailed distribution of nitrogen abundances in the single and binary samples may be different, possibly reflecting differences in their evolutionary history. Conclusions. The first comparative study of single and binary B-type supergiants has revealed that the main sequence may be significantly wider than previously assumed, extending to Teff = 20?000 K. Some marginal differences in single and binary atmospheric parameters and abundances have been identified, possibly implying non-standard evolution for some of the sample. This sample as a whole has implications for several aspects of our understanding of the evolutionary status of blue supergiants. [ABSTRACT FROM AUTHOR]

Context. With growing evidence for the existence of very massive stars at subsolar metallicity, there is an increased need for corresponding stellar evolution models. Aims. We present a dense model grid with a tailored input chemical composition appropriate for the Large Magellanic Cloud (LMC). Methods. We use a one-dimensional hydrodynamic stellar evolution code, which accounts for rotation, transport of angular momentum by magnetic fields, and stellar wind mass loss to compute our detailed models. We calculate stellar evolution models with initial masses from 70 to 500 M⊙ and with initial surface rotational velocities from 0 to 550 kms-1, covering the core-hydrogen burning phase of evolution. Results. We find our rapid rotators to be strongly influenced by rotationally induced mixing of helium, with quasi-chemically homogeneous evolution occurring for the fastest rotating models. Above 160 M⊙, homogeneous evolution is also established through mass loss, producing pure helium stars at core hydrogen exhaustion independent of the initial rotation rate. Surface nitrogen enrichment is also found for slower rotators, even for stars that lose only a small fraction of their initial mass. For models above ~150 M⊙ at zero age, and for models in the whole considered mass range later on, we find a considerable envelope inflation due to the proximity of these models to their Eddington limit. This leads to a maximum ZAMS surface temperature of ~56000 K, at ~180 M⊙, and to an evolution of stars in the mass range 50 M⊙ . . . 100 M⊙a to the regime of luminous blue variables in the Hertzsprung-Russell diagram with high internal Eddington factors. Inflation also leads to decreasing surface temperatures during the chemically homogeneous evolution of stars above ~180 M⊙ . Conclusions. The cool surface temperatures due to the envelope inflation in our models lead to an enhanced mass loss, which prevents stars at LMC metallicity from evolving into pair-instability supernovae. The corresponding spin-down will also prevent very massive LMC stars to produce long-duration gamma-ray bursts, which might, however, originate from lower masses. [ABSTRACT FROM AUTHOR]

The VLT-FLAMES Tarantula Survey is an ESO Large Program from which we have obtained multi-epoch optical spectroscopy of over 800 massive stars in the 30 Doradus region of the Large Magellanic Cloud. This unprecedented dataset is being used to address outstanding questions in how massive stars evolve from the early main sequence to their deaths as core collapse supernovae. Here we focus on the rotation properties of the population of presumably single O stars and use binary population synthesis predictions to show that the rapid rotators among this population likely are post-interaction binaries. The same type of population synthesis can be used to study the mass function of massive young clusters. We argue - on the basis of predictions for the Arches and Quintuplet clusters - that a sizable fraction of the very massive WNh stars in 30 Doradus may also have such a binary interaction history. We single out the WNh star discovered in the VFTS, VFTS682, and discuss its properties. [ABSTRACT FROM AUTHOR]

Context. The Tarantula Nebula in the Large Magellanic Cloud is our closest view of a starburst region and is the ideal environment to investigate important questions regarding the formation, evolution and final fate of the most massive stars. Aims. We analyze the multiplicity properties of the massive O-type star population observed through multi-epoch spectroscopy in the framework of the VLT-FLAMES Tarantula Survey. With 360 O-type stars, this is the largest homogeneous sample of massive stars analyzed to date. Methods. We use multi-epoch spectroscopy and variability analysis to identify spectroscopic binaries. We also use a Monte-Carlo method to correct for observational biases. By modeling simultaneously the observed binary fraction, the distributions of the amplitudes of the radial velocity variations and the distribution of the time scales of these variations, we constrain the intrinsic current binary fraction and period and mass-ratio distributions. Results. We observe a spectroscopic binary fraction of 0:35±0:03, which corresponds to the fraction of objects displaying statistically significant radial velocity variations with an amplitude of at least 20 km s-1.We compute the intrinsic binary fraction to be 0:51±0:04. We adopt power-laws to describe the intrinsic period and mass-ratio distributions: f (log10 P=d) ~ (log10 P=d)π (with log10 P=d in the range 0.15-3.5) and f (q) ~ q" with 0:1 ⩽ q = M2=M1 ⩽ 1:0. The power-law indexes that best reproduce the observed quantities are π = -0:45 ± 0:30 and ⩽ = -1:0 ± 0:4. The period distribution that we obtain thus favours shorter period systems compared to an Öpik law (π = 0). The mass ratio distribution is slightly skewed towards low mass ratio systems but remains incompatible with a random sampling of a classical mass function (⩽ = -2:35). The binary fraction seems mostly uniform across the field of view and independent of the spectral types and luminosity classes. The binary fraction in the outer region of the field of view (r > 7:80, i.e. π117 pc) and among the O9.7 I/II objects are however significantly lower than expected from statistical fluctuations. The observed and intrinsic binary fractions are also lower for the faintest objects in our sample (Ks > 15:5 mag), which results from observational effects and the fact that our O star sample is not magnitude-limited but is defined by a spectral-type cutoff. We also conclude that magnitude-limited investigations are biased towards larger binary fractions. Conclusions. Using the multiplicity properties of the O stars in the Tarantula region and simple evolutionary considerations, we estimate that over 50% of the current O star population will exchange mass with its companion within a binary system. This shows that binary interaction is greatly affecting the evolution and fate of massive stars, and must be taken into account to correctly interpret unresolved populations of massive stars. [ABSTRACT FROM AUTHOR]

We report the detection of a strong, organized magnetic field in the secondary component of the massive O8III/I+O7.5V/III double-lined spectroscopic binary system HD 47129 (Plaskett's star) in the context of the Magnetism in Massive Stars survey. Eight independent Stokes V observations were acquired using the Echelle SpectroPolarimetric Device for the Observations of Stars (ESPaDOnS) spectropolarimeter at the Canada–France–Hawaii Telescope and the Narval spectropolarimeter at the Télescope Bernard Lyot. Using least-squares deconvolution we obtain definite detections of signal in Stokes V in three observations. No significant signal is detected in the diagnostic null (N) spectra. The Zeeman signatures are broad and track the radial velocity of the secondary component; we therefore conclude that the rapidly rotating secondary component is the magnetized star. Correcting the polarized spectra for the line and continuum of the (sharp-lined) primary, we measured the longitudinal magnetic field from each observation. The longitudinal field of the secondary is variable and exhibits extreme values of −810 ± 150 and +680 ± 190 G, implying a minimum surface dipole polar strength of 2850 ± 500 G. In contrast, we derive an upper limit (3σ) to the primary's surface magnetic field of 230 G. The combination of a strong magnetic field and rapid rotation leads us to conclude that the secondary hosts a centrifugal magnetosphere fed through a magnetically confined wind. We revisit the properties of the optical line profiles and X-ray emission – previously interpreted as a consequence of colliding stellar winds – in this context. We conclude that HD 47129 represents a heretofore unique stellar system – a close, massive binary with a rapidly rotating, magnetized component – that will be a rich target for further study. [ABSTRACT FROM PUBLISHER]

The presence of a nearby companion alters the evolution of massive stars in binary systems, leading to phenomena such as stellar mergers, x-ray binaries, and gamma-ray bursts. Unambiguous constraints on the fraction of massive stars affected by binary interaction were lacking. We simultaneously measured all relevant binary characteristics in a sample of Galactic massive 0 stars and quantified the frequency and nature of binary interactions. More than 70% of all massive stars will exchange mass with a companion, leading to a binary merger in one-third of the cases. These numbers greatly exceed previous estimates and imply that binary interaction dominates the evolution of massive stars, with implications for populations of massive stars and their supernovae. [ABSTRACT FROM AUTHOR]

The hydrogen-rich outer layers of massive stars can be removed by interactions with a binary companion. Theoretical models predict that this stripping produces a population of hot helium stars of ~2 to 8 solar masses (M☉), however, only one such system has been identified thus far. We used ultraviolet photometry to identify potential stripped helium stars then investigated 25 of them using optical spectroscopy. We identified stars with high temperatures (~60,000 to 100,000 kelvin), high surface gravities, and hydrogen-depleted surfaces; 16 stars also showed binary motion. These properties match expectations for stars with initial masses of 8 to 25 M☉ that were stripped by binary interaction. Their masses fall in the gap between subdwarf helium stars and Wolf-Rayet stars. We propose that these stars could be progenitors of stripped-envelope supernovae. [ABSTRACT FROM AUTHOR]