Search

We investigate how the multiplicity of binary, triple and quadruple star systems changes as the systems evolve from the zero-age main-sequence to the Hubble time. We find the change in multiplicity fractions over time for each data set, identify the number of changes to the orbital configuration and the dominant underlying physical mechanism responsible for each configuration change. Finally, we identify key properties of the binaries which survive the evolution. We use the stellar evolution population synthesis code Multiple Stellar Evolution (MSE) to follow the evolution of $3 \times 10^4$ of each 1+1 binaries, 2+1 triples, 3+1 quadruples and 2+2 quadruples. The coupled stellar and orbital evolution are computed each iteration. The systems are assumed to be isolated and to have formed in situ. We generate data sets for two different black hole natal kick mean velocity distributions (sigma = 10 km/s and sigma = 50 km/s and with and without the inclusion of stellar fly-bys. Our fiducial model has a mean black hole natal kick velocity if sigma = 10 km/s and includes stellar fly-bys. Each system has at least one star with an initial mass larger than 10 solar masses. All data will be publicly available. We find that at the end of the evolution the large majority of systems are single stars in every data set (> 85%). As the number of objects in the initial system increases, so too does the final non-single system fraction. The single fractions of final systems in our fiducial model are 87.8 $\pm$ 0.2 % for the 2 + 2s, 88.8 $\pm$ 0.3 % for the 3 + 1s, 92.3 $\pm$ 0.2 % for the 2 + 1s and 98.9 $\pm$ 0.3 % for the 1 + 1s., Comment: 20 Pages, 10 Figures, 1 Table Submitted to A&A 17 December 2024

Stars stripped of their hydrogen-rich envelopes through binary interaction are thought to be responsible for both hydrogen-poor supernovae and the hard ionizing radiation observed in low-$Z$ galaxies. A population of these stars was recently observed for the first time, but their prevalence remains unknown. In preparation for such measurements, we estimate the mass distribution of hot, stripped stars using a population synthesis code that interpolates over detailed single and binary stellar evolution tracks. We predict that for a constant star-formation rate of $1 \,M_\odot$/yr and regardless of metallicity, a population contains $\sim$30,000 stripped stars with mass $M_{\rm strip}>1M_\odot$ and $\sim$4,000 stripped stars that are sufficiently massive to explode ($M_{\rm strip}>2.6M_\odot$). Below $M_{\rm strip}=5M_\odot$, the distribution is metallicity-independent and can be described by a power law with the exponent $\alpha \sim -2$. At higher masses and lower metallicity ($Z \lesssim 0.002$), the mass distribution exhibits a drop. This originates from the prediction, frequently seen in evolutionary models, that massive low-metallicity stars do not expand substantially until central helium burning or later and therefore cannot form long-lived stripped stars. With weaker line-driven winds at low metallicity, this suggests that neither binary interaction nor wind mass loss can efficiently strip massive stars at low metallicity. As a result, a "helium-star desert" emerges around $M_{\rm strip} =15\, M_\odot$ at $Z=0.002$, covering an increasingly large mass range with decreasing metallicity. We note that these high-mass stars are those that potentially boost a galaxy's He$^+$-ionizing radiation and that participate in the formation of merging black holes. This "helium-star desert" therefore merits further study., Comment: 16 pages, 9 figures, submitted to A&A

Surveys in the Milky Way and Large Magellanic Cloud revealed that the majority of massive stars will interact with companions during their lives. However, knowledge of the binary properties of massive stars at low metallicity, which approaches the conditions of the Early Universe, remains sparse. We present the Binarity at LOw Metallicity (BLOeM) campaign - an ESO large programme designed to obtain 25 epochs of spectroscopy for 929 massive stars in the SMC - the lowest metallicity conditions in which multiplicity is probed to date (Z = 0.2 Zsun). BLOeM will provide (i) the binary fraction, (ii) the orbital configurations of systems with periods P < 3 yr, (iii) dormant OB+BH binaries, and (iv) a legacy database of physical parameters of massive stars at low metallicity. The stars are observed with the LR02 setup of the giraffe instrument of the Very Large Telescope (3960-4570A, resolving power R=6200; typical signal-to-noise ratio S/N=70-100). This paper utilises the first 9 epochs obtained over a three-month time. We describe the survey and data reduction, perform a spectral classification of the stacked spectra, and construct a Hertzsprung-Russell diagram of the sample via spectral-type and photometric calibrations. The sample covers spectral types from O4 to F5, spanning the effective temperature and luminosity ranges 6.5

On the route towards merging neutron stars and stripped-envelope supernovae, binary population synthesis predicts a large number of post-interaction systems with massive stars that have stripped off their outer layers. Yet, observations of such stars in the intermediate-mass regime below the Wolf-Rayet masses are rare. Using X-Shooting ULLYSES (XShootU) data, we discovered three partially stripped star + Be/Oe binaries in the Magellanic Clouds. We analyzed the UV and optical spectra using the PoWR model atmosphere code by superimposing model spectra corresponding to each component. The estimated current masses of the partially stripped stars fall within the intermediate mass range of 4-8 $M_{\odot}$. These objects are overluminous for their stellar masses, matching core He-burning luminosities. Their Be/Oe secondaries have much higher masses than their stripped primaries (mass ratio > 2). All three partially stripped stars show significant nitrogen enrichment and carbon and oxygen depletion on their surfaces. Additionally, one of our sample stars exhibits significant helium enrichment. Our study provides the first comprehensive determination of the wind parameters of partially stripped stars in the intermediate mass range. The wind mass-loss rates of these stars are found to be on the order of $10^{-7} M_\odot$ yr$^{-1}$, which is over ten times higher than that of OB stars of the same luminosity. Current evolutionary models characterizing this phase typically employ OB or WR mass-loss rates, which underestimate or overestimate stripped stars' mass-loss rates by an order of magnitude. Binary evolution models indicate that the observed primaries had initial masses of 12-17 $M_{\odot}$, making them potential candidates for stripped-envelope supernovae that form neutron stars. If they survive the explosion, these systems may become Be X-ray binaries and later double neutron stars., Comment: Accepted for publication in Astronomy & Astrophysics

The quest for quiet or dormant black holes has been ongoing since several decades. Ellipsoidal variables possibly indicate the existence of a very high-mass invisible companion and are thought to be one of the best ways to find such dormant black holes. This, however, is not a panacea as we show here with one example. We indeed report the discovery of a new semi-detached interacting binary, V1315 Cas, discovered as an ellipsoidal variable. Using data from photometric surveys (ASAS-SN, TESS) and high-resolution spectroscopy, we derived a nearly circular orbit with an orbital period of $P_{\rm{orb}}$=34.54 d. The binary system consists of an evolved F-type star primary that is likely still filling its Roche lobe and a B-type star secondary. Using \textsc{phoebe}2, we derived the following masses and radii: for the primary, $M_p =0.84 \pm 0.03 \, M_\odot$ and $R_p =18.51^{+0.12}_{-0.07} \, R_\odot$; for the secondary, $M_s =7.3 \pm 0.3 \,M_\odot$ and $R_s =4.02^{+2.3}_{-2.0}\,R_\odot$. Modeling the evolution of the system with MESA, we found an age of $\sim$7.7e7 years. The system is at the end of a period of rapid non-conservative mass transfer that reversed its mass ratio, while significantly widening its orbit. The primary shows carbon depletion and nitrogen overabundance, indicative of CNO processed material being exposed due to mass transfer. An infrared excess as well as stationary H$\alpha$ emission suggest the presence of a circumstellar or circumbinary disc. V1315 Cas will likely become a detached stripped star binary., Comment: Accepted in MNRAS

Standard binary evolutionary models predict a significant population of core helium-burning stars that lost their hydrogen-rich envelope after mass transfer via Roche-lobe overflow. However, there is a scarcity of observations of such stripped stars in the intermediate mass regime (~1.5 - 8$ M_{\odot}$), which are thought to be prominent progenitors of SN Ib/c. Especially at low metallicity, a significant fraction of these stars is expected to be only partially stripped, retaining a significant amount of hydrogen on their surfaces. For the first time, we discovered a partially stripped massive star in a binary with a Be-type companion located in the Small Magellanic Cloud (SMC) using a detailed spectroscopic analysis. The stripped-star nature of the primary is revealed by the extreme CNO abundance pattern and very high luminosity-to-mass ratio, which suggest that the primary is likely shell-hydrogen burning. Our target SMCSGS-FS 69 is the most luminous and most massive system among the known stripped star + Be binaries, with Mstripped ~3$ M_{\odot}$ and MBe ~17$ M_{\odot}$. Binary evolutionary tracks suggest an initial mass of Mini $\gtrsim 12 M_{\odot}$ for the stripped star and predict it to be in a transition phase towards a hot compact He star, which will eventually produce a stripped-envelope supernova. Our target marks the first representative of a so-far missing evolutionary stage in the formation pathway of Be X-ray binaries and double neutron star mergers., Comment: Accepted for publication as A&A Letter to the Editor

Observations of individual massive stars, super-luminous supernovae, gamma-ray bursts, and gravitational-wave events involving spectacular black-hole mergers, indicate that the low-metallicity Universe is fundamentally different from our own Galaxy. Many transient phenomena will remain enigmatic until we achieve a firm understanding of the physics and evolution of massive stars at low metallicity (Z). The Hubble Space Telescope has devoted 500 orbits to observe 250 massive stars at low Z in the ultraviolet (UV) with the COS and STIS spectrographs under the ULLYSES program. The complementary ``X-Shooting ULLYSES'' (XShootU) project provides enhanced legacy value with high-quality optical and near-infrared spectra obtained with the wide-wavelength coverage X-shooter spectrograph at ESO's Very Large Telescope. We present an overview of the XShootU project, showing that combining ULLYSES UV and XShootU optical spectra is critical for the uniform determination of stellar parameters such as effective temperature, surface gravity, luminosity, and abundances, as well as wind properties such as mass-loss rates in function of Z. As uncertainties in stellar and wind parameters percolate into many adjacent areas of Astrophysics, the data and modelling of the XShootU project is expected to be a game-changer for our physical understanding of massive stars at low Z. To be able to confidently interpret James Webb Space Telescope (JWST) spectra of the first stellar generations, the individual spectra of low Z stars need to be understood, which is exactly where XShootU can deliver., Comment: Accepted in A&A - 35 Pages, 12 Figures, 4 Tables, 2 Large Tables

The origin of the observed population of Wolf-Rayet (WR) stars in low-metallicity (low-Z) 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. 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. Multi-epoch UV and optical spectra of AzV 14 are analyzed using the non-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. AzV 14 is a close binary system consisting of two similar main sequence stars with masses of 32 Msol. Both stars have weak stellar winds with mass-loss rates of log $\dot{M}$ = -7.7. Binary evolutionary models can explain the empirically derived stellar and orbital parameters. The model predicts that the primary will evolve into a WR star with T = 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 = 50 kK and are predicted to have compared to other WR stars increased oxygen abundances. This model prediction is supported by a spectroscopic analysis of a WR star in the SMC. We hypothesize that the populations of WR stars in low-Z 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., Comment: 21 pages (13 main body + 8 appendix), 16 figures, 9 tables

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 for stars beyond 50 solar masses 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 initialised before the donor star reaches its maximum expansion. The only significant difference was found for the mass distribution of massive binary black hole mergers (total mass > 50 solar masses) 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., Comment: Accepted for publication in MNRAS

The LIGO/Virgo gravitational--wave observatories have detected 50 BH-BH coalescences. This sample is large enough to have allowed several recent studies to draw conclusions about the branching ratios between isolated binaries versus dense stellar clusters as the origin of double BHs. It has also led to the exciting suggestion that the population is highly likely to contain primordial black holes. Here we demonstrate that such conclusions cannot yet be robust, because of the large current uncertainties in several key aspects of binary stellar evolution. These include the development and survival of a common envelope, the mass and angular momentum loss during binary interactions, mixing in stellar interiors, pair-instability mass loss and supernova outbursts. Using standard tools such as the population synthesis codes StarTrack and COMPAS and the detailed stellar evolution code MESA, we examine as a case study the possible future evolution of Melnick 34, the most massive known binary star system. We show that, despite its well-known orbital architecture, various assumptions regarding stellar and binary physics predict a wide variety of outcomes: from a close BH-BH binary (which would lead to a potentially detectable coalescence), through a wide BH-BH binary (which might be seen in microlensing observations), or a Thorne-Zytkow object, to a complete disruption of both objects by pair-instability supernovae. Thus since the future of massive binaries is inherently uncertain, sound predictions about the properties of BH-BH systems are highly challenging at this time. Consequently, drawing conclusions about the formation channels for the LIGO/Virgo BH-BH merger population is premature., Comment: comments welcomed

Common-envelope (CE) evolution in massive binary systems is thought to be one of the most promising channels for the formation of compact binary mergers. In the case of merging binary black holes (BBHs), the essential CE phase takes place at a stage when the first BH is already formed and the companion star expands as a supergiant. We study which BH binaries with supergiant companions will evolve through and potentially survive a CE phase. To this end, we compute envelope binding energies from detailed massive stellar models at different evolutionary stages and metallicities. We make multiple physically extreme choices of assumptions that favor easier CE ejection as well as account for recent advancements in mass transfer stability criteria. We find that even with the most optimistic assumptions, a successful CE ejection in BH (and also NS) binaries is only possible if the donor is a massive convective-envelope giant, a red supergiant (RSG). In other words, pre-CE progenitors of BBH mergers are BH binaries with RSG companions. We find that due to its influence on the radial expansion of massive giants, metallicity has an indirect but a very strong effect on the envelope structure and binding energies of RSGs. Our results suggest that merger rates from population synthesis models could be severely overestimated, especially at low metallicity. Additionally, the lack of observed RSGs with luminosities above log($L/L_{\odot}$) = 5.6-5.8, corresponding to stars with $M > 40 M_{\odot}$, puts into question the viability of the CE channel for the formation of the most massive BBH mergers. Either such RSGs elude detection due to very short lifetimes, or they do not exist and the CE channel can only produce BBH systems with total mass $< 50 M_{\odot}$. We discuss an alternative CE scenario, in which a partial envelope ejection is followed by a phase of possibly long and stable mass transfer., Comment: 20 pages + App., accepted for publication in A&A. For $\lambda_{\rm CE}$ fits, see: https://ftp.science.ru.nl/astro/jklencki/

Gaia16aye was a binary microlensing event discovered in the direction towards the northern Galactic disc and was one of the first microlensing events detected and alerted to by the Gaia space mission. Its light curve exhibited five distinct brightening episodes, reaching up to I=12 mag, and it was covered in great detail with almost 25,000 data points gathered by a network of telescopes. We present the photometric and spectroscopic follow-up covering 500 days of the event evolution. We employed a full Keplerian binary orbit microlensing model combined with the motion of Earth and Gaia around the Sun to reproduce the complex light curve. The photometric data allowed us to solve the microlensing event entirely and to derive the complete and unique set of orbital parameters of the binary lensing system. We also report on the detection of the first-ever microlensing space-parallax between the Earth and Gaia located at L2. The properties of the binary system were derived from microlensing parameters, and we found that the system is composed of two main-sequence stars with masses 0.57$\pm$0.05 $M_\odot$ and 0.36$\pm$0.03 $M_\odot$ at 780 pc, with an orbital period of 2.88 years and an eccentricity of 0.30. We also predict the astrometric microlensing signal for this binary lens as it will be seen by Gaia as well as the radial velocity curve for the binary system. Events such as Gaia16aye indicate the potential for the microlensing method of probing the mass function of dark objects, including black holes, in directions other than that of the Galactic bulge. This case also emphasises the importance of long-term time-domain coordinated observations that can be made with a network of heterogeneous telescopes., Comment: accepted for publication in A&A, 24 pages, 10 figures, tables with the data will be available electronically

The distributions of the initial main-sequence binary parameters are one of the key ingredients in obtaining evolutionary predictions for compact binary (BH-BH / BH-NS / NS-NS) merger rates. Until now, such calculations were done under the assumption that initial binary parameter distributions were independent. Here, we implement empirically derived inter-correlated distributions of initial binary parameters primary mass (M1), mass ratio (q), orbital period (P), and eccentricity (e). Unexpectedly, the introduction of inter-correlated initial binary parameters leads to only a small decrease in the predicted merger rates by a factor of 2 $-$ 3 relative to the previously used non-correlated initial distributions. The formation of compact object mergers in the isolated classical binary evolution favors initial binaries with stars of comparable masses (q = 0.5 $-$ 1) at intermediate orbital periods (log P (days) = 2 $-$ 4). New distributions slightly shift the mass ratios towards smaller values with respect to the previously used flat q distribution, which is the dominant effect decreasing the rates. New orbital periods only negligibly increase the number of progenitors. Additionally, we discuss the uncertainty of merger rate predictions associated with possible variations of the massive-star initial mass function (IMF). We argue that evolutionary calculations should be normalized to a star formation rate (SFR) that is obtained from the observed amount of UV light at wavelength 1500{\AA} (SFR indicator). In this case, contrary to recent reports, the uncertainty of the IMF does not affect the rates by more than a factor of 2. Any change to the IMF slope for massive stars requires a change of SFR in a way that counteracts the impact of IMF variations on the merger rates. In contrast, we suggest that the uncertainty in cosmic SFR at low metallicity can be a significant factor at play., Comment: accepted for publication in A&A

All ten LIGO/Virgo binary black hole (BH-BH) coalescences reported from the O1/O2 runs have near zero effective spins. There are only three potential explanations of this fact. If the BH spin magnitudes are large then (i) either both BH spin vectors must be nearly in the orbital plane or (ii) the spin angular momenta of the BHs must be oppositely directed and similar in magnitude. Or, (iii) the BH spin magnitudes are small. We test the third hypothesis within the framework of the classical isolated binary evolution scenario of the BH-BH merger formation. We test three models of angular momentum transport in massive stars: a mildly efficient transport by meridional currents (as employed in the Geneva code), an efficient transport by the Tayler-Spruit magnetic dynamo (as implemented in the MESA code), and a very-efficient transport (as proposed by Fuller et al.) to calculate natal BH spins. We allow for binary evolution to increase the BH spins through accretion and account for the potential spin-up of stars through tidal interactions. Additionally, we update the calculations of the stellar-origin BH masses, include revisions to the history of star formation and to the chemical evolution across cosmic time. We find that we can match simultaneously the observed BH-BH merger rate density, BH masses, and effective spins. Models with efficient angular momentum transport are favored. The updated stellar-mass weighted gas-phase metallicity evolution now used in our models appears to be a key in better reproducing the LIGO/Virgo merger rate estimate. Mass losses during the pair-instability pulsation supernova phase are likely overestimated if the merger GW170729 hosts a BH more massive than 50 Msun. We also estimate rate of BH-NS mergers from recent LIGO/Virgo observations. Our updated models of BH-BH, BH-NS and NS-NS mergers are now publicly available at www.syntheticuniverse.org., Comment: A&A accepted: revised comparison with LIGO/Virgo (41 pages)