Context. The majority of the Wolf–Rayet (WR) stars represent the stripped cores of evolved massive stars who lost most of their hydrogen envelope. Wind stripping in single stars is expected to be ...inefficient in producing WR stars in metal-poor environments such as the Small Magellanic Cloud (SMC). While binary interaction can also produce WR stars at low metallicity, it is puzzling that the fraction of WR binaries appears to be about 40%, independent of the metallicity. Aim. We aim to use the recently determined physical properties of the twelve known SMC WR stars to explore their possible formation channels through comparisons with stellar models. Methods. We used the MESA stellar evolution code to construct two grids of stellar models with SMC metallicity. One of these consists of models of rapidly rotating single stars, which evolve in part or completely chemically homogeneously. In a second grid, we analyzed core helium burning stellar models assuming constant hydrogen and helium gradients in their envelopes. Results. We find that chemically homogeneous evolution is not able to account for the majority of the WR stars in the SMC. However, in particular the apparently single WR star SMC AB12, and the double WR system SMC AB5 (HD 5980) appear consistent with this channel. We further find a dichotomy in the envelope hydrogen gradients required to explain the observed temperatures of the SMC WR stars. Shallow gradients are found for the WR stars with O star companions, while much steeper hydrogen gradients are required to understand the group of hot apparently single WR stars. Conclusions. The derived shallow hydrogen gradients in the WR component of the WR+O star binaries are consistent with predictions from binary models where mass transfer occurs early, in agreement with their binary properties. Since the hydrogen profiles in evolutionary models of massive stars become steeper with time after the main sequence, we conclude that most of the hot (Teff > 60 kK ) apparently single WR stars lost their envelope after a phase of strong expansion, e.g., as the result of common envelope evolution with a lower mass companion. The so far undetected companions, either main sequence stars or compact objects, are then expected to still be present. A corresponding search might identify the first immediate double black hole binary progenitor with masses as high as those detected in GW150914.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
Context. The evolution of massive stars is strongly influenced by internal mixing processes such as semiconvection, convective core overshooting, and rotationally induced mixing. None of these ...processes are currently well constrained. Aims. We investigate models for massive stars in the Small Magellanic Cloud (SMC), for which stellar-wind mass loss is less important than for their metal-rich counterparts. We aim to constrain the various mixing efficiencies by comparing model results to observations. Methods. For this purpose, we use the stellar-evolution code MESA to compute more than 60 grids of detailed evolutionary models for stars with initial masses of 9…100 M⊙, assuming different combinations of mixing efficiencies of the various processes in each grid. Our models evolve through core hydrogen and helium burning, such that they can be compared with the massive main sequence and supergiant population of the SMC. Results. We find that for most of the combinations of the mixing efficiencies, models in a wide mass range spend core-helium burning either only as blue supergiants, or only as red supergiants. The latter case corresponds to models that maintain a shallow slope of the hydrogen/helium (H/He) gradient separating the core and the envelope of the models. Only a small part of the mixing parameter space leads to models that produce a significant number of blue and red supergiants, which are both in abundance in the SMC. Some of our grids also predict a cut-off in the number of red supergiants above log L/L⊙ = 5…5.5. Interestingly, these models contain steep H/He gradients, as is required to understand the hot, hydrogen-rich Wolf-Rayet stars in the SMC. We find that unless it is very fast, rotation has a limited effect on the H/He profiles in our models. Conclusions. While we use specific implementations of the considered mixing processes, they comprehensively probe the two first-order structural parameters, the core mass and the H/He gradient in the core-envelope interface. Our results imply that in massive stars, mixing during the main-sequence evolution leads to a moderate increase in the helium core masses, and also that the H/He gradients above the helium cores become very steep. Our model grids can be used to further refine the various mixing efficiencies with the help of future observational surveys of the massive stars in the SMC, and thereby help to considerably reduce the uncertainties in models of massive star evolution.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
Stripped-envelope stars form in binary systems after losing mass through Roche-lobe overflow. They bear astrophysical significance as sources of UV and ionizing radiation in older stellar populations ...and, if sufficiently massive, as stripped supernova progenitors. Binary evolutionary models predict that they are common, but only a handful of subdwarfs with B-type companions are known. The question is whether a large population of such systems has evaded detection as a result of biases, or whether the model predictions are wrong. We reanalyze the well-studied post-interaction binary φ Persei. Recently, new data have improved the orbital solution of the system, which contains an ~1.2M⊙ stripped-envelope star and a rapidly rotating ~9.6M⊙ Be star. We compare with an extensive grid of evolutionary models using a Bayesian approach and constrain the initial masses of the progenitor to 7.2 ± 0.4M⊙ for the stripped star and 3.8 ± 0.4M⊙ for the Be star. The system must have evolved through near-conservative mass transfer. These findings are consistent with earlier studies. The age we obtain, 57 ± 9 Myr, is in excellent agreement with the age of the α Persei cluster. We note that neither star was initially massive enough to produce a core-collapse supernova, but mass exchange pushed the Be star above the mass threshold. We find that the subdwarf is overluminous for its mass by almost an order of magnitude, compared to the expectations for a helium core burning star. We can only reconcile this if the subdwarf resides in a late phase of helium shell burning, which lasts only 2–3% of the total lifetime as a subdwarf. Assuming continuous star formation implies that up to ~50 less evolved, dimmer subdwarfs exist for each system similar to φ Persei, but have evaded detection so far. Our findings can be interpreted as a strong indication that a substantial population of stripped-envelope stars indeed exists, but has so far evaded detection because of observational biases and lack of large-scale systematic searches.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
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.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
Context.
Binary evolution can result in fast-rotating stars through accretion of angular momentum during mass-transfer phases. These fast-rotating stars are predicted to be observable as Be stars. ...Current models struggle to produce a satisfactory description of Be star populations, even though numerous pieces of observational evidence indicate that the accretion process might be the dominant Be formation channel.
Aims.
Given the distinct uncertainties in detailed binary evolution calculations, we investigate a rigorous and model-independent upper limit for the production of Be stars through binary interaction and aim to confront this limit with observations of Be stars in young star clusters.
Methods.
Using extreme assumptions, we calculate the number ratio of post-interaction to pre-interaction binary systems in a coeval population. This ratio describes an upper limit to Be star formation through mass transfer. A detailed comparison is made between our derived upper limit and relevant observations of Be stars, which allows us to probe several aspects of binary star physics.
Results.
We find that in coeval populations, binary interaction can at most account for one-third of all main-sequence stars being Be stars. Near the cluster turn-off region, this limit appears to be realised in the clusters studied. Away from the turn-off, a good fit to the observed Be fraction as a function of mass is obtained by applying simple assumptions about which systems undergo unstable mass-transfer produces.
Conclusions.
We find that assuming distinct physics, binary evolution alone can in principle match the high numbers of Be stars that are observed in open clusters. Whether the required binary physics is realised in nature remains to be investigated.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
Context.
The recent gravitational wave measurements have demonstrated the existence of stellar mass black hole binaries. It is essential for our understanding of massive star evolution to identify ...the contribution of binary evolution to the formation of double black holes.
Aims.
A promising way to progress is investigating the progenitors of double black hole systems and comparing predictions with local massive star samples, such as the population in 30 Doradus in the Large Magellanic Cloud (LMC).
Methods.
With this purpose in mind, we analysed a large grid of detailed binary evolution models at LMC metallicity with initial primary masses between 10 and 40
M
⊙
, and identified the model systems that potentially evolve into a binary consisting of a black hole and a massive main-sequence star. We then derived the observable properties of such systems, as well as peculiarities of the OB star component.
Results.
We find that ∼3% of the LMC late-O and early-B stars in binaries are expected to possess a black hole companion when stars with a final helium core mass above 6.6
M
⊙
are assumed to form black holes. While the vast majority of them may be X-ray quiet, our models suggest that these black holes may be identified in spectroscopic binaries, either by large amplitude radial velocity variations (≳50 km s
−1
) and simultaneous nitrogen surface enrichment, or through a moderate radial velocity (≳10 km s
−1
) and simultaneous rapid rotation of the OB star. The predicted mass ratios are such that main-sequence companions can be excluded in most cases. A comparison to the observed OB+WR binaries in the LMC, Be and X-ray binaries, and known massive black hole binaries supports our conclusion.
Conclusions.
We expect spectroscopic observations to be able to test key assumptions in our models, with important implications for massive star evolution in general and for the formation of double black hole mergers in particular.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
Context.
The majority of massive stars are part of binary systems. In about a quarter of these, the companions are so close that mass transfer occurs while they undergo core hydrogen burning, first ...on the thermal and then on the nuclear timescale. The nuclear timescale mass transfer leads to observational counterparts: the semi-detached so-called massive Algol binaries. These systems may provide urgently needed tests of the physics of mass transfer. However, comprehensive model predictions for these systems are sparse.
Aims.
We use a large grid of detailed evolutionary models of short-period massive binaries and follow-up population synthesis calculations to derive probability distributions of the observable properties of massive Algols and their descendants.
Methods.
Our results are based on ∼10 000 binary model sequences calculated with the stellar evolution code MESA, using a metallicity suitable for the Large Magellanic Cloud (LMC), covering initial donor masses between 10
M
⊙
and 40
M
⊙
and initial orbital periods above 1.4 d. These models include internal differential rotation and magnetic angular momentum transport, non-conservative mass and angular momentum transfer between the binary components, and time-dependent tidal coupling.
Results.
Our models imply ∼30, or ∼3% of the ∼1000, core hydrogen burning O-star binaries in the LMC to be currently in the semi-detached phase. Our donor models are up to 25 times more luminous than single stars of an identical mass and effective temperature, which agrees with the observed Algols. A comparison of our models with the observed orbital periods and mass ratios implies rather conservative mass transfer in some systems, while a very inefficient one in others. This is generally well reproduced by our spin-dependent mass transfer algorithm, except for the lowest considered masses. The observations reflect the slow increase of the surface nitrogen enrichment of the donors during the semi-detached phase all the way to CNO equilibrium. We also investigate the properties of our models after core hydrogen depletion of the donor star, when these models correspond to Wolf-Rayet or helium+OB star binaries.
Conclusions.
A dedicated spectroscopic survey of massive Algol systems may allow to derive the dependence of the efficiency of thermal timescale mass transfer on the binary parameters, as well as the efficiency of semiconvective mixing in the stellar interior. This would be a crucial step towards reliable binary models up to the formation of supernovae and compact objects.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
ABSTRACT
We aim to identify and characterize binary systems containing red supergiant (RSG) stars in the Small Magellanic Cloud (SMC) using a newly available ultraviolet (UV) point source catalogue ...obtained using the Ultraviolet Imaging Telescope (UVIT) on board AstroSat. We select a sample of 560 SMC RSGs based on photometric and spectroscopic observations at optical wavelengths and cross-match this with the far-UV point source catalogue using the UVIT F172M filter, finding 88 matches down to mF172M = 20.3 ABmag, which we interpret as hot companions to the RSGs. Stellar parameters (luminosities, effective temperatures, and masses) for both components in all 88 binary systems are determined and we find mass distributions in the ranges 6.1 < M/M⊙ < 22.3 for RSGs and 3.7 <M/M⊙ < 15.6 for their companions. The most massive RSG binary system in the SMC has a combined mass of 32 ± 4 M⊙, with a mass ratio (q) of 0.92. By simulating observing biases, we find an intrinsic multiplicity fraction of $18.8\, \pm \, 1.5\, {{\ \rm per\ cent}}$ for mass ratios in the range 0.3 < q < 1.0 and orbital periods approximately in the range $3 \lt \log P \rm days \lt 8$. By comparing our results with those of a similar mass on the main sequence, we determine the fraction of single stars to be ∼20 per cent and argue that the orbital period distribution declines rapidly beyond log P ∼ 3.5. We study the mass-ratio distribution of RSG binary systems and find that a uniform distribution best describes the data below 14 M⊙. Above 15 M⊙, we find a lack of high mass ratio systems.
Context.
Massive star evolution at low metallicity is closely connected to many fields in high-redshift astrophysics, but is poorly understood so far. Because of its metallicity of ∼0.2
Z
⊙
, its ...proximity, and because it is currently forming stars, the Small Magellanic Cloud (SMC) is a unique laboratory in which to study metal-poor massive stars.
Aims.
We seek to improve the understanding of this topic using available SMC data and a comparison to stellar evolution predictions.
Methods.
We used a recent catalog of spectral types in combination with
Gaia
magnitudes to calculate temperatures and luminosities of bright SMC stars. By comparing these with literature studies, we tested the validity of our method, and using
Gaia
data, we estimated the completeness of stars in the catalog as a function of luminosity. This allowed us to obtain a nearly complete view of the most luminous stars in the SMC. We also calculated the extinction distribution, the ionizing photon production rate, and the star formation rate.
Results.
Our results imply that the SMS hosts only ∼30 very luminous main-sequence stars (
M
≥ 40
M
⊙
;
L
≳ 3 ⋅ 10
5
L
⊙
), which are far fewer than expected from the number of stars in the luminosity range 3 ⋅ 10
4
<
L
/
L
⊙
< 3 ⋅ 10
5
and from the typically quoted star formation rate in the SMC. Even more striking, we find that for masses above
M
≳ 20
M
⊙
, stars in the first half of their hydrogen-burning phase are almost absent. This mirrors a qualitatively similar peculiarity that is known for the Milky Way and Large Magellanic Cloud. This amounts to a lack of hydrogen-burning counterparts of helium-burning stars, which is more pronounced for higher luminosities. We derived the H I ionizing photon production rate of the current massive star population. It agrees with the H
α
luminosity of the SMC.
Conclusions.
We argue that a declining star formation rate or a steep initial mass function are unlikely to be the sole explanations for the dearth of young bright stars. Instead, many of these stars might be embedded in their birth clouds, although observational evidence for this is weak. We discuss implications for the role that massive stars played in cosmic reionization, and for the top end of the initial mass function.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
Context.
The origin of the observed population of Wolf-Rayet (WR) stars in low-metallicity galaxies, such as the Small Magellanic Cloud (SMC), is not yet understood. Standard, single-star ...evolutionary models predict that WR stars should stem from very massive O-type star progenitors, but these are very rare. On the other hand, binary evolutionary models predict that WR stars could originate from primary stars in close binaries.
Aims.
We conduct an analysis of the massive O star, AzV 14, to spectroscopically determine its fundamental and stellar wind parameters, which are then used to investigate evolutionary paths from the O-type to the WR stage with stellar evolutionary models.
Methods.
Multi-epoch UV and optical spectra of AzV 14 are analyzed using the non-local thermodynamic equilibrium (LTE) stellar atmosphere code PoWR. An optical TESS light curve was extracted and analyzed using the PHOEBE code. The obtained parameters are put into an evolutionary context, using the MESA code.
Results.
AzV 14 is a close binary system with a period of
P
= 3.7058 ± 0.0013 d. The binary consists of two similar main sequence stars with masses of
M
1, 2
≈ 32
M
⊙
. Both stars have weak stellar winds with mass-loss rates of log
Ṁ
/(
M
⊙
yr
−1
) = −7.7 ± 0.2. Binary evolutionary models can explain the empirically derived stellar and orbital parameters, including the position of the AzV 14 components on the Hertzsprung-Russell diagram, revealing its current age of 3.3 Myr. The model predicts that the primary will evolve into a WR star with
T
eff
≈ 100 kK, while the secondary, which will accrete significant amounts of mass during the first mass transfer phase, will become a cooler WR star with
T
eff
≈ 50 kK. Furthermore, WR stars that descend from binary components that have accreted significant amount of mass are predicted to have increased oxygen abundances compared to other WR stars. This model prediction is supported by a spectroscopic analysis of a WR star in the SMC.
Conclusions.
Inspired by the binary evolutionary models, we hypothesize that the populations of WR stars in low-metallicity galaxies may have bimodal temperature distributions. Hotter WR stars might originate from primary stars, while cooler WR stars are the evolutionary descendants of the secondary stars if they accreted a significant amount of mass. These results may have wide-ranging implications for our understanding of massive star feedback and binary evolution channels at low metallicity.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK