Context.
About 20% of all B-type stars are classical Be stars – stars whose spectra imply the presence of a circumstellar decretion disk. The disk phenomenon is strongly correlated with rapid ...rotation, the origin of which remains unclear. It may be rooted in single- or binary-star evolution. In the framework of the binary channel, the initially more massive star transfers mass and angular momentum to the original secondary, which becomes a Be star. The system then evolves into a Be binary with a post-main-sequence companion, which, depending on the companion mass, may later be disrupted in a supernova event. Hence, if the binary channel dominates the formation of Be stars, one may expect a strong lack of close Be binaries with main sequence (MS) companions.
Aims.
We want to test the prediction of the binary channel. Through an extensive, star-by-star review of the literature of a magnitude-limited sample of Galactic early-type Be stars, we investigate whether Be binaries with MS companions are known to exist.
Methods.
Our sample is constructed from the BeSS database and cross-matched with all available literature on the individual stars. Archival and amateur spectra are used to verify the existing literature when conflicting reports are found.
Results.
Out of an initial list of 505 Be stars, we compile a final sample of 287 Galactic Be stars earlier than B1.5 with
V
≤ 12 mag. Out of those, 13 objects were reported as Be binaries with known post-MS companions (i.e., compact objects or helium stars) and 11 as binaries with unknown, uncertain or debated companions. We find no confirmed reports of Be binaries with MS companions. For the remaining 263 targets, no significant reports of multiplicity exist in the literature, implying that they are either Be binaries with faint companions, or truly single.
Conclusions.
The clear lack of reported MS companions to Be stars, which stands in contrast to the high number of detected B+B MS binaries, strongly supports the hypothesis that early-type Be stars are binary interaction products that spun up after mass and angular momentum transfer from a companion star. Taken at face value, our results may suggest that a large majority of the early-type Be stars have formed through binary mass-transfer.
Context.
The intriguing binary
LS V +22 25
(LB-1) has drawn much attention following claims of it being a single-lined spectroscopic binary with a 79-day orbit comprising a B-type star and a ≈70
M
⊙
...black hole – the most massive stellar black hole reported to date. Subsequent studies demonstrated a lack of evidence for a companion of such great mass. Recent analyses have implied that the primary star is a stripped He-rich star with peculiar sub-solar abundances of heavy elements, such as Mg and Fe. However, the nature of the secondary, which was proposed to be a black hole, a neutron star, or a main sequence star, remains unknown.
Aims.
Based on 26 newly acquired spectroscopic observations secured with the HERMES and FEROS spectrographs covering the orbit of the system, we perform an orbital analysis and spectral disentangling of LB-1 to elucidate the nature of the system.
Methods.
To derive the radial velocity semi-amplitude
K
2
of the secondary and extract the spectra of the two components, we used two independent disentangling methods: the shift-and-add technique and Fourier disentangling with FDBinary. We used atmosphere models to constrain the surface properties and abundances.
Results.
Our disentangling and spectral analysis shows that LB-1 contains two components of comparable brightness in the optical. The narrow-lined primary, which we estimate to contribute ≈55% in the optical, has spectral properties that suggest that it is a stripped star: it has a small spectroscopic mass (≈1
M
⊙
) for a B-type star and it is He- and N-rich. Unlike previous reports, the abundances of heavy elements are found to be solar. The “hidden” secondary, which contributes about 45% of the optical flux, is a rapidly rotating (
v
sin
i
≈ 300 km s
−1
) B3 V star with a decretion disk – a Be star. As a result of its rapid rotation and dilution, the photospheric absorption lines of the secondary are not readily apparent in the individual observations. We measure a semi-amplitude for this star of
K
2
= 11.2 ± 1.0 km s
−1
and adopting a mass of
M
2
= 7 ± 2
M
⊙
typical for B3 V stars, we derive an orbital mass for the stripped primary of
M
1
= 1.5 ± 0.4
M
⊙
. The orbital inclination of 39 ± 4° implies a near-critical rotation for the Be secondary (
v
eq
≈ 470 km s
−1
).
Conclusions.
LB-1 does not contain a compact object. Instead, it is a rare Be binary system consisting of a stripped star (the former mass donor) and a Be star rotating at near its critical velocity (the former mass accretor). This system is a clear example that binary interactions play a decisive role in the production of rapid stellar rotators and Be stars.
Context. It is now well established that the majority of massive stars reside in multiple systems. However, the effect of multiplicity is not sufficiently understood, resulting in a plethora of ...uncertainties about the end stages of massive-star evolution. In order to investigate these uncertainties, it is useful to study massive stars just before their demise. Classical Wolf-Rayet stars represent the final end stages of stars at the upper-mass end. The multiplicity fraction of these stars was reported to be ∼0.4 in the Galaxy but no correction for observational biases has been attempted. Aims. The aim of this study is to conduct a homogeneous radial-velocity survey of a magnitude-limited (V ≤ 12) sample of Galactic Wolf-Rayet stars to derive their bias-corrected multiplicity properties. The present paper focuses on 12 northern Galactic carbon-rich (WC) Wolf-Rayet stars observable with the 1.2 m Mercator telescope on the island of La Palma. Methods. We homogeneously measured relative radial velocities (RVs) for carbon-rich Wolf-Rayet stars using cross-correlation. Variations in the derived RVs were used to flag binary candidates. We investigated probable orbital configurations and provide a first correction of observational biases through Monte-Carlo simulations. Results. Of the 12 northern Galactic WC stars in our sample, seven show peak-to-peak RV variations larger than 10 km s−1, which we adopt as our detection threshold. This results in an observed spectroscopic multiplicity fraction of 0.58 with a binomial error of 0.14. In our campaign, we find a clear lack of short-period (P < ∼100 d), indicating that a large number of Galactic WC binaries likely reside in long-period systems. Finally, our simulations show that at the 10% significance level, the intrinsic multiplicity fraction of the Galactic WC population is at least 0.72.
Context.
Classical Wolf-Rayet (WR) stars are massive, hydrogen-depleted, post main-sequence stars that exhibit emission-line dominated spectra. For a given metallicity
Z
, stars exceeding a certain ...initial mass
M
single
WR
(Z) can reach the WR phase through intrinsic mass-loss or eruptions (single-star channel). In principle, stars of lower masses can reach the WR phase via stripping through binary interactions (binary channel). Because winds become weaker at low
Z
, it is commonly assumed that the binary channel dominates the formation of WR stars in environments with low metallicity such as the Small and Large Magellanic Clouds (SMC, LMC). However, the reported WR binary fractions of 30−40% in the SMC (
Z
= 0.002) and LMC (
Z
= 0.006) are comparable to that of the Galaxy (
Z
= 0.014), and no evidence for the dominance of the binary channel at low
Z
could be identified observationally. Here, we explain this apparent contradiction by considering the minimum initial mass
M
spec
WR
(Z) needed for the stripped product to appear as a WR star.
Aims.
By constraining
M
spec
WR
(Z) and
M
single
WR
(Z), we estimate the importance of binaries in forming WR stars as a function of
Z
.
Methods.
We calibrated
M
spec
WR
using the lowest-luminosity WR stars in the Magellanic Clouds and the Galaxy. A range of
M
single
WR
values were explored using various evolution codes. We estimated the additional contribution of the binary channel by considering the interval
M
spec
WR
(Z),
M
single
WR
(Z), which characterizes the initial-mass range in which the binary channel can form additional WR stars.
Results.
The WR-phenomenon ceases below luminosities of log
L
≈ 4.9, 5.25, and 5.6
L
⊙
in the Galaxy, the LMC, and the SMC, respectively, which translates to minimum He-star masses of 7.5, 11, 17
M
⊙
and minimum initial masses of
M
spec
WR
= 18, 23, 37
M
⊙
. Stripped stars with lower initial masses in the respective galaxies would tend not to appear as WR stars. The minimum mass necessary for self-stripping,
M
single
WR
(Z), is strongly model-dependent, but it lies in the range 20−30, 30−60, and ≳40
M
⊙
for the Galaxy, LMC, and SMC, respectively. We find that that the additional contribution of the binary channel is a non-trivial and model-dependent function of
Z
that cannot be conclusively claimed to be monotonically increasing with decreasing
Z
.
Conclusions.
The WR spectral appearance arises from the presence of strong winds. Therefore, both
M
spec
WR
and
M
single
WR
increase with decreasing metallicity. Considering this, we show that one should not a-priori expect that binary interactions become increasingly important in forming WR stars at low
Z
, or that the WR binary fraction grows with decreasing
Z
.
The study of massive stars in different metallicity environments is a central topic of current stellar research. The spectral analysis of massive stars requires adequate model atmospheres. The ...computation of such models is difficult and time-consuming. Therefore, spectral analyses are greatly facilitated if they can refer to existing grids of models. Here we provide grids of model atmospheres for OB-type stars at metallicities corresponding to the Small and Large Magellanic Clouds, as well as to solar metallicity. In total, the grids comprise 785 individual models. The models were calculated using the state-of-the-art Potsdam Wolf-Rayet (PoWR) model atmosphere code. The parameter domain of the grids was set up using stellar evolution tracks. For all these models, we provide normalized and flux-calibrated spectra, spectral energy distributions, feedback parameters such as ionizing photons, Zanstra temperatures, and photometric magnitudes. The atmospheric structures (the density and temperature stratification) are available as well. All these data are publicly accessible through the PoWR website.
Context. Spectroscopic analysis remains the most common method to derive masses of massive stars, the most fundamental stellar parameter. While binary orbits and stellar pulsations can provide much ...sharper constraints on the stellar mass, these methods are only rarely applicable to massive stars. Unfortunately, spectroscopic masses of massive stars heavily depend on the detailed physics of model atmospheres. Aims. We demonstrate the impact of a consistent treatment of the radiative pressure on inferred gravities and spectroscopic masses of massive stars. Specifically, we investigate the contribution of line and continuum transitions to the photospheric radiative pressure. We further explore the effect of model parameters, e.g., abundances, on the deduced spectroscopic mass. Lastly, we compare our results with the plane-parallel TLUSTY code, commonly used for the analysis of massive stars with photospheric spectra. Methods. We calculate a small set of O-star models with the Potsdam Wolf-Rayet (PoWR) code using different approaches for the quasi-hydrostatic part. These models allow us to quantify the effect of accounting for the radiative pressure consistently. We further use PoWR models to show how the Doppler widths of line profiles and abundances of elements such as iron affect the radiative pressure, and, as a consequence, the derived spectroscopic masses. Results. Our study implies that errors on the order of a factor of two in the inferred spectroscopic mass are to be expected when neglecting the contribution of line and continuum transitions to the radiative acceleration in the photosphere. Usage of implausible microturbulent velocities, or the neglect of important opacity sources such as Fe, may result in errors of approximately 50% in the spectroscopic mass. A comparison with TLUSTY model atmospheres reveals a very good agreement with PoWR at the limit of low mass-loss rates.
Context.
In the era of gravitational wave astrophysics and with the precise astrometry of billions of stellar sources, the hunt for compact objects is more alive than ever. Rarely seen massive ...binaries with a compact object are a crucial phase in the evolution towards compact object mergers. With the upcoming third
Gaia
data release (DR3), the first
Gaia
astrometric orbital solutions for binary sources will become available, potentially revealing many such binaries.
Aims.
We investigate how many black holes (BHs) with massive main-sequence dwarf companions (OB+BH binaries) are expected to be detected as binaries in
Gaia
DR3 and at the end of the nominal 5-year mission. We estimate how many of those are identifiable as OB+BH binaries and discuss the distributions of the masses of both components as well as of their orbital periods. We also explore how different BH-formation scenarios affect these distributions.
Methods.
We apply observational constraints to tailored models for the massive star population, which assume a direct collapse and no kick upon BH formation, to estimate the fraction of OB+BH systems that will be detected as binaries by
Gaia
, and consider these the fiducial results. These OB+BH systems follow a distance distribution according to that of the second Alma Luminous Star catalogue (ALS II). We use a method based on astrometric data to identify binaries with a compact object and investigate how many of the systems detected as binaries are identifiable as OB+BH binaries. Different scenarios for BH natal kicks and supernova mechanisms are explored and compared to the fiducial results.
Results.
In the fiducial case we conservatively estimate that 77% of the OB+BH binaries in the ALS II will be detected as binaries in DR3, of which 89% will be unambiguously identifiable as OB+BH binaries. By the end of the nominal 5-year mission, the detected fraction will increase to 85%, of which 82% will be identifiable. The 99% confidence intervals on these fractions are of the order of a few percent. These fractions become smaller for different BH-formation scenarios.
Conclusions.
Assuming direct collapse and no natal kick, we expect to find around 190 OB+BH binaries in
Gaia
DR3 among the sources in the ALS II, which increases the known sample of OB+BH binaries by more than a factor of 20 and covers an uncharted parameter space of long-period binaries (10 ≲
P
≲ 1000 d). Our results further show that the size and properties of the OB+BH population that is identifiable using
Gaia
DR3 will contain crucial observational constraints that will help us improve our understanding of BH formation. An additional ∼5 OB+BH binaries could be identified at the end of the nominal 5-year mission, which are expected to have either very short (
P
≲ 10 d) or long periods (
P
≳ 1000 d).
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.
Context. Massive Wolf-Rayet (WR) stars are evolved massive stars (M sub(i)> or = 20 M sub(middot in circle)) characterized by strong mass-loss. Hypothetically, they can form either as single stars or ...as mass donors in close binaries. About 40% of all known WR stars are confirmed binaries, raising the question as to the impact of binarity on the WR population. Studying WR binaries is crucial in this context, and furthermore enable one to reliably derive the elusive masses of their components, making them indispensable for the study of massive stars. Aims. By performing a spectral analysis of all multiple WR systems in the Small Magellanic Cloud (SMC), we obtain the full set of stellar parameters for each individual component. Mass-luminosity relations are tested, and the importance of the binary evolution channel is assessed. Methods. The spectral analysis is performed with the Potsdam Wolf-Rayet (PoWR) model atmosphere code by superimposing model spectra that correspond to each component. Evolutionary channels are constrained using the Binary Population and Spectral Synthesis (BPASS) evolution tool. Results. Significant hydrogen mass fractions (0.1 <X sub(H)< 0.4) are detected in all WN components. A comparison with mass-luminosity relations and evolutionary tracks implies that the majority of the WR stars in our sample are not chemically homogeneous. The WR component in the binary AB6 is found to be very luminous (logLapproximate 6.3 L sub(middot in circle)) given its orbital mass (approximate10 M sub(middot in circle)), presumably because of observational contamination by a third component. Evolutionary paths derived for our objects suggest that Roche lobe overflow had occurred in most systems, affecting their evolution. However, the implied initial masses (> or =60 M sub(middot in circle)) are large enough for the primaries to have entered the WR phase, regardless of binary interaction. Conclusions. Together with the results for the putatively single SMC WR stars, our study suggests that the binary evolution channel does not dominate the formation of WR stars at SMC metallicity.
ABSTRACT
Massive stars are the progenitors of black holes and neutron stars, the mergers of which can be detected with gravitational waves (GW). The expansion of massive stars is one of the key ...factors affecting their evolution in close binary systems, but it remains subject to large uncertainties in stellar astrophysics. For population studies and predictions of GW sources, the stellar expansion is often simulated with the analytic formulae from Hurley et al. (2000). These formulae need to be extrapolated and are often considered outdated. In this work, we present five different prescriptions developed from 1D stellar models to constrain the maximum expansion of massive stars. We adopt these prescriptions to investigate how stellar expansion affects mass transfer interactions and in turn the formation of GW sources. We show that limiting radial expansion with updated 1D stellar models, when compared to the use of Hurley et al. (2000) radial expansion formulae, does not significantly affect GW source properties (rates and masses). This is because most mass transfer events leading to GW sources are initialized in our models before the donor star reaches its maximum expansion. The only significant difference was found for the mass distribution of massive binary black hole mergers (Mtot > 50 M⊙) formed from stars that may evolve beyond the Humphreys–Davidson limit, whose radial expansion is the most uncertain. We conclude that understanding the expansion of massive stars and the origin of the Humphrey–Davidson limit is a key factor for the study of GW sources.