Showing 1–50 of 78 results for author: Marchant, P

X-Shooting ULLYSES: Massive Stars at Low Metallicity

Authors: Jorick S. Vink, Paul Crowther, Alex Fullerton, Miriam Garcia, Fabrice Martins, Nidia Morrell, Lida Oskinova, Nicole St. Louis, Asif ud-Doula, Andreas Sander, Hugues Sana, Jean-Claude Bouret, Brankica Kubatova, Pablo Marchant, Lucimara P. Martins, Aida Wofford, Jacco van Loon, O. Grace Telford, Ylva Götberg, Dominic Bowman, Christi Erba, Venu Kalari, The XShootU Collaboration

Abstract: The Hubble Space Telescope has devoted 500 orbits to observing 250 massive stars with low metallicity in the ultraviolet (UV) range within the framework of the ULLYSES program. The X-Shooting ULLYSES (XShootU) project enhances the legacy value of this UV dataset by providing high-quality optical and near-infrared spectra, which are acquired using the wide-wavelength-coverage X-shooter spectrograph… ▽ More The Hubble Space Telescope has devoted 500 orbits to observing 250 massive stars with low metallicity in the ultraviolet (UV) range within the framework of the ULLYSES program. The X-Shooting ULLYSES (XShootU) project enhances the legacy value of this UV dataset by providing high-quality optical and near-infrared spectra, which are acquired using the wide-wavelength-coverage X-shooter spectrograph at ESO's Very Large Telescope. XShootU emphasises the importance of combining UV with optical spectra for the consistent determination of key stellar parameters such as effective temperature, surface gravity, luminosity, abundances, and wind characteristics including mass-loss rates as a function of metallicity. Since uncertainties in these parameters have implications across various branches of astrophysics, the data and modelling generated by the XShootU project are poised to significantly advance our understanding of massive stars at low metallicity. This is particularly crucial for confidently interpreting JWST data of the earliest stellar generations, making XShootU a unique resource for comprehending individual spectra of low-metallicity stars. △ Less

Submitted 30 April, 2024; originally announced May 2024.

Comments: 6 pages, 6 figures. ESO Large Programme Overview

Journal ref: ESO Messenger, 2024

A magnetic massive star has experienced a stellar merger

Authors: A. J. Frost, H. Sana, L. Mahy, G. Wade, J. Barron, J. -B. Le Bouquin, A. Mérand, F. R. N. Schneider, T. Shenar, R. H. Barbá, D. M. Bowman, M. Fabry, A. Farhang, P. Marchant, N. I. Morrell, J. V. Smoker

Abstract: Massive stars (those larger than 8 solar masses at formation) have radiative envelopes that cannot sustain a dynamo, the mechanism that produces magnetic fields in lower-mass stars. Despite this, approximately 7\% of massive stars have observed magnetic fields, the origin of which is debated. We used multi-epoch interferometric and spectroscopic observations to characterize HD 148937, a binary sys… ▽ More Massive stars (those larger than 8 solar masses at formation) have radiative envelopes that cannot sustain a dynamo, the mechanism that produces magnetic fields in lower-mass stars. Despite this, approximately 7\% of massive stars have observed magnetic fields, the origin of which is debated. We used multi-epoch interferometric and spectroscopic observations to characterize HD 148937, a binary system of two massive stars. We found that only one star is magnetic and that it appears younger than its companion. The system properties and a surrounding bipolar nebula can be reproduced with a model in which two stars merged (in a previous triple system) to produce the magnetic massive star. Our results provide observational evidence that magnetic fields form in at least some massive stars through stellar mergers. △ Less

Submitted 15 April, 2024; originally announced April 2024.

Comments: Full paper with supplementary materials. 59 pages, 18 figures

Journal ref: Science, 11 Apr 2024, Vol 384, Issue 6692, pp. 214-217

Analytic approximations for massive close post-mass transfer binary systems

Authors: Christoph Schürmann, Norbert Langer, Joana A. Kramer, Pablo Marchant, Chen Wang, Koushik Sen

Abstract: Massive binary evolution models are needed to predict massive star populations in star forming galaxies, the supernova diversity, and the number and properties of gravitational wave sources. Such models are often computed using so called rapid binary evolution codes, which approximate the evolution of the binary components based on detailed single star models. However, about one third of the inter… ▽ More Massive binary evolution models are needed to predict massive star populations in star forming galaxies, the supernova diversity, and the number and properties of gravitational wave sources. Such models are often computed using so called rapid binary evolution codes, which approximate the evolution of the binary components based on detailed single star models. However, about one third of the interacting massive binary stars undergo mass transfer during core hydrogen burning (Case A mass transfer), whose outcome is difficult to derive from single star models. Here, we use a large grid of detailed binary evolution models for primaries in the initial mass range 10 to 40 Solar masses of LMC and SMC composition, to derive analytic fits for the key quantities needed in rapid binary evolution codes, i.e., the duration of core hydrogen burning, and the resulting donor star mass. Systems with shorter orbital periods produce up to 50% lighter stripped donors and have a up to 30% larger lifetime than wider systems. We find that both quantities depend strongly on the initial binary orbital period, but that the initial mass ratio and the mass transfer efficiency of the binary have little impact on the outcome. Our results are easily parameterisable and can be used to capture the effects of Case A mass transfer more accurately in rapid binary evolution codes. △ Less

Submitted 12 April, 2024; originally announced April 2024.

Comments: 10 pages

Investigating the Chemically Homogeneous Evolution Channel and its Role in the Formation of the Enigmatic Binary Black Hole Progenitor Candidate HD 5980

Authors: K. Sharpe, L. A. C. van Son, S. E. de Mink, R. Farmer, P. Marchant, G. Koenigsberger

Abstract: Chemically homogeneous evolution (CHE) is a promising channel for forming massive binary black holes. The enigmatic, massive Wolf-Rayet (WR) binary HD 5980 A&B has been proposed to have formed through this channel. We investigate this claim by comparing its observed parameters with CHE models. Using MESA, we simulate grids of close massive binaries then use a Bayesian approach to compare them with… ▽ More Chemically homogeneous evolution (CHE) is a promising channel for forming massive binary black holes. The enigmatic, massive Wolf-Rayet (WR) binary HD 5980 A&B has been proposed to have formed through this channel. We investigate this claim by comparing its observed parameters with CHE models. Using MESA, we simulate grids of close massive binaries then use a Bayesian approach to compare them with the stars' observed orbital period, masses, luminosities, and hydrogen surface abundances. The most probable models, given the observational data, have initial periods ~3 days, widening to the present-day ~20 day orbit as a result of mass loss -- correspondingly, they have very high initial stellar masses ($\gtrsim$150 M$_\odot$). We explore variations in stellar wind-mass loss and internal mixing efficiency, and find that models assuming enhanced mass-loss are greatly favored to explain HD 5980, while enhanced mixing is only slightly favoured over our fiducial assumptions. Our most probable models slightly underpredict the hydrogen surface abundances. Regardless of its prior history, this system is a likely binary black hole progenitor. We model its further evolution under our fiducial and enhanced wind assumptions, finding that both stars produce black holes with masses ~19-37 M$_\odot$. The projected final orbit is too wide to merge within a Hubble time through gravitational waves alone. However, the system is thought to be part of a 2+2 hierarchical multiple. We speculate that secular effects with the (possible) third and fourth companions may drive the system to promptly become a gravitational-wave source. △ Less

Submitted 19 February, 2024; originally announced February 2024.

Comments: 23 pages, 9 figures, 2 tables, 2 appendices. Accepted for publication in ApJ

Modelling time-dependent convective penetration in 1D stellar evolution

Authors: Cole Johnston, Mathias Michielsen, Evan H. Anders, Mathieu Renzo, Matteo Cantiello, P. Marchant, Jared A. Goldberg, Richard H. D. Townsend, Gautham Sabhahit, Adam S. Jermyn

Abstract: 1D stellar evolution calculations produce uncertain predictions for quantities like the age, core mass, core compactness, and nucleo-synthetic yields; a key source of uncertainty is the modeling of interfaces between regions that are convectively stable and those that are not. Theoretical and numerical work has demonstrated that there should be numerous processes adjacent to the convective boundar… ▽ More 1D stellar evolution calculations produce uncertain predictions for quantities like the age, core mass, core compactness, and nucleo-synthetic yields; a key source of uncertainty is the modeling of interfaces between regions that are convectively stable and those that are not. Theoretical and numerical work has demonstrated that there should be numerous processes adjacent to the convective boundary that induce chemical and angular momentum transport, as well as modify the thermal structure of the star. One such process is called convective penetration, wherein vigorous convection extends beyond the nominal convective boundary and alters both the composition and thermal structure. In this work, we incorporate the process of convective penetration in stellar evolution calculations using the stellar evolution software instrument mesa. We implement convective penetration according to the description presented by Anders et al. (2022a) to calculate a grid of models from the pre main sequence to He core depletion. The extent of the convective penetration zone is self-consistently calculated at each time step without introducing new free parameters. We find both a substantial penetration zone in all models with a convective core and observable differences to global stellar properties such as the luminosity and radius. We preset how the predicted radial extent of the penetration zone scales with the total stellar mass, age and the metallicity of the star. We discuss our results in the context of existing numerical and observational studies. △ Less

Submitted 13 December, 2023; originally announced December 2023.

Comments: Accepted for publication in The Astrophysical Journal, 16 pages, 8 figures, accompanying repository (https://github.com/mathren/dbconvpen)

An upper limit on the spins of merging binary black holes formed through binary evolution

Authors: Pablo Marchant, Philipp Podsiadlowski, Ilya Mandel

Abstract: As gravitational wave detectors improve, observations of black hole (BH) mergers will provide the joint distribution of their masses and spins. This will be a critical benchmark to validate formation scenarios. Merging binary BHs formed through isolated binary evolution require both components to be stripped of their hydrogen envelopes before core-collapse. The rotation rates of such stripped star… ▽ More As gravitational wave detectors improve, observations of black hole (BH) mergers will provide the joint distribution of their masses and spins. This will be a critical benchmark to validate formation scenarios. Merging binary BHs formed through isolated binary evolution require both components to be stripped of their hydrogen envelopes before core-collapse. The rotation rates of such stripped stars are constrained by their surface critical rotation, restricting their angular momentum content at core-collapse. We use stripped star models at low metallicities ($Z_\odot/10$, $Z_\odot/50$ and $Z_\odot/250$) to determine the spins of BHs produced by critically rotating stellar progenitors. To study how such progenitors can arise, we consider their formation through chemically homogeneous evolution (CHE). We use a semianalytical model to study the final spins of CHE binaries, and compare our results against available detailed population synthesis models. We find that above BH masses of $\simeq 25M_\odot$, the dimensionless spin of critically rotating stripped stars ($a = Jc/(GM^2$)) is below unity. This results in an exclusion region at high chirp masses and effective spins that cannot be populated by binary evolution. CHE can produce binaries where both BHs hit this limit, producing a pile-up at the boundary of the excluded region. Highly spinning BHs arise from very low-metallicity CHE systems with short delay times, which merge at higher redshifts. On the other hand, the contribution of CHE to merging binary BHs in the third observing run of the LVK collaboration is expected to be dominated by systems with low spins ($χ_\mathrm{eff}<0.5$) which merge near redshift zero. Owing to its higher projected sensitivity and runtime, the fourth observing run of the LVK collaboration can potentially place constraints on the high spin population and the existence of a limit set by critical rotation. △ Less

Submitted 23 November, 2023; originally announced November 2023.

Comments: Submitted to A&A. Abstract abridged

The Evolution of Massive Binary Stars

Authors: Pablo Marchant, Julia Bodensteiner

Abstract: Massive stars play a major role in the evolution of their host galaxies, and serve as important probes of the distant Universe. It has been established that the majority of massive stars reside in close binaries and will interact with their companion stars during their lifetime. Such interactions drastically alter their life cycles and complicate our understanding of their evolution, but are also… ▽ More Massive stars play a major role in the evolution of their host galaxies, and serve as important probes of the distant Universe. It has been established that the majority of massive stars reside in close binaries and will interact with their companion stars during their lifetime. Such interactions drastically alter their life cycles and complicate our understanding of their evolution, but are also responsible for the production of interesting and exotic interaction products. - Extensive observation campaigns with well-understood detection sensitivities have allowed to convert the observed properties into intrinsic characteristics, facilitating a direct comparison to theory. - Studies of large samples of massive stars in our Galaxy and the Magellanic Clouds have unveiled new types of interaction products, providing critical constraints on the mass transfer phase and the formation of compact objects. - The direct detection of gravitational waves has revolutionized the study of stellar mass compact objects, providing a new window to study massive star evolution. Their formation processes are, however, still unclear. The known sample of compact object mergers will grow by orders of magnitude in the coming decade, turning into the best understood astrophysical population. △ Less

Submitted 3 November, 2023; originally announced November 2023.

Comments: This is the authors' submitted version. To be published with significant updates on the Annual Reviews of Astronomy and Astrophysics. Please refer to the published version once available, which will become available under this url <https://www.annualreviews.org/journal/astro>

Forming merging double compact objects with stable mass transfer

Authors: Annachiara Picco, Pablo Marchant, Hugues Sana, Gijs Nelemans

Abstract: Merging double compact objects (CO) represent the inferred sources of every detected gravitational wave (GW) signal, thus modeling their progenitors is important to constrain stellar evolution theory. Stable mass transfer (MT) between a donor star and a black hole is one of the proposed tightening mechanisms to form binary black holes merging within the age of the universe. We aim to assess the po… ▽ More Merging double compact objects (CO) represent the inferred sources of every detected gravitational wave (GW) signal, thus modeling their progenitors is important to constrain stellar evolution theory. Stable mass transfer (MT) between a donor star and a black hole is one of the proposed tightening mechanisms to form binary black holes merging within the age of the universe. We aim to assess the potential of stable non conservative mass transfer to produce the pairings of COs: black holes (BHs), neutron stars (NSs) and white dwarfs (WDs). We study the conditions required for mass transfer between a star and a CO to be stable and to lead to merging binary COs. We use published results for the response of the stellar radii to rapid mass loss; coupled with analytical models for orbital evolution, we determine the boundary for unstable MT and the post interaction properties of binaries undergoing stable MT. We investigate the impact of angular momentum loss prescription in the hardening by accounting for isotropic re emission from the accretor vicinity and mass outflow from the Lagrangian point L2. Stable MT in systems with a CO + Roche lobe filling star, in the limit of isotropic re-emission, is shown to be able to form any pair of merging COs apart from WD + BH. Considering mass outflow from L2, the resulting parameter space for GW progenitors is shifted towards smaller initial mass ratios, ruling out the formation of NS + NS while allowing the production of merging WD + BH pairs. We compare our results with single-degenerate binaries and find that conditions for stable MT to operate are present in nature. We show that stable MT in the isotropic re-emission limit can produce merging binary BHs with mass ratios consistent with the majority of inferred sources of the third Gravitational Wave Transient Catalogue. Angular momentum loss from L2 lifts the achievable final mass ratio. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: 18 pages, 17 figures

Journal ref: A&A 681, A31 (2024)

A massive helium star with a sufficiently strong magnetic field to form a magnetar

Authors: T. Shenar, G. A. Wade, P. Marchant, S. Bagnulo, J. Bodensteiner, D. M. Bowman, A. Gilkis, N. Langer, A. -N. Chene, L. Oskinova, T. Van Reeth, H. Sana, N. St-Louis, A. Soares de Oliveira, H. Todt, S. Toonen

Abstract: Magnetars are highly magnetized neutron stars; their formation mechanism is unknown. Hot helium-rich stars with spectra dominated by emission lines are known as Wolf-Rayet stars. We observe the binary system HD 45166 using spectropolarimetry, finding that it contains a Wolf-Rayet star with a mass of 2 solar masses and a magnetic field of 43 kilogauss. Stellar evolution calculations indicate that t… ▽ More Magnetars are highly magnetized neutron stars; their formation mechanism is unknown. Hot helium-rich stars with spectra dominated by emission lines are known as Wolf-Rayet stars. We observe the binary system HD 45166 using spectropolarimetry, finding that it contains a Wolf-Rayet star with a mass of 2 solar masses and a magnetic field of 43 kilogauss. Stellar evolution calculations indicate that this component will explode as a type Ib or IIb supernova, and the strong magnetic field favors a magnetar remnant. We propose that the magnatized Wolf-Rayet star formed by the merger of two lower mass helium stars. △ Less

Submitted 16 August, 2023; originally announced August 2023.

Comments: Published in Science on the 18 August 2023. Radial velocities, spectra, and software available in: https://zenodo.org/record/8042656 ESO press release: www.eso.org/public/news/eso2313

The B-type Binaries Characterisation Programme II. VFTS 291: A stripped star from a recent mass transfer phase

Authors: J. I. Villaseñor, D. J. Lennon, A. Picco, T. Shenar, P. Marchant, N. Langer, P. L. Dufton, F. Nardini, C. J. Evans, J. Bodensteiner, S. E. de Mink, Y. Götberg, I. Soszyński, W. D. Taylor, H. Sana

Abstract: Recent studies of massive binaries with putative black hole companions have uncovered a phase of binary evolution that has not been observed before, featuring a bloated stripped star that very recently ceased transferring mass to a main-sequence companion. In this study, we focus on the candidate system VFTS 291, a binary with an orbital period of 108 d and a high semi-amplitude velocity (… ▽ More Recent studies of massive binaries with putative black hole companions have uncovered a phase of binary evolution that has not been observed before, featuring a bloated stripped star that very recently ceased transferring mass to a main-sequence companion. In this study, we focus on the candidate system VFTS 291, a binary with an orbital period of 108 d and a high semi-amplitude velocity ($K_1=93.7\pm0.2$ km s$^{-1}$). Through our analysis of the disentangled spectra of the two components, together with dynamical and evolutionary arguments, we identify a narrow-lined star of ~1.5-2.5 $M_\odot$ dominating the spectrum, and an early B-type main-sequence companion of $13.2\pm1.5$ $M_\odot$. The low mass of the narrow-lined star, and the high mass ratio, suggest that VFTS 291 is a post-mass-transfer system, with the narrow-lined star being bloated and stripped of its hydrogen-rich envelope, sharing many similarities with other recently discovered stripped stars. Our finding is supported by our detailed binary evolution models, which indicate that the system can be well explained by an initial configuration consisting of an 8.1 $M_\odot$ primary with an 8 $M_\odot$ companion in a 7 d orbital period. While some open questions remain, particularly concerning the surface helium enrichment of the stripped star and the rotational velocity of the companion, we expect that high-resolution spectroscopy may help reconcile our estimates with theory. Our study highlights the importance of multi-epoch spectroscopic surveys to identify and characterize binary interaction products, and provides important insights into the evolution of massive binary stars. △ Less

Submitted 11 September, 2023; v1 submitted 15 July, 2023; originally announced July 2023.

Comments: 25 pages, 17 figures. Published by MNRAS, Advance Access

X-Shooting ULLYSES: massive stars at low metallicity. I. Project Description

Authors: Jorick S. Vink, A. Mehner, P. A. Crowther, A. Fullerton, M. Garcia, F. Martins, N. Morrell, L. M. Oskinova, N. St-Louis, A. ud-Doula, A. A. C. Sander, H. Sana, J. -C. Bouret, B. Kubatova, P. Marchant, L. P. Martins, A. Wofford, J. Th. van Loon, O. Grace Telford, Y. Gotberg, D. M. Bowman, C. Erba, V. M. Kalari, M. Abdul-Masih, T. Alkousa , et al. (56 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 1 June, 2023; v1 submitted 10 May, 2023; originally announced May 2023.

Comments: Accepted in A&A - 35 Pages, 12 Figures, 4 Tables, 2 Large Tables

Journal ref: A&A 675, A154 (2023)

ALMA detection of CO rotational line emission in red supergiant stars of the massive young star cluster RSGC1 -- Determination of a new mass-loss rate prescription for red supergiants

Authors: Leen Decin, Anita M. S. Richards, Pablo Marchant, Hugues Sana

Abstract: [Abridged] Aim: We aim to derive a new mass-loss rate prescription for RSGs that is not afflicted with some uncertainties inherent in preceding studies. Methods: We have observed CO rotational line emission towards a sample of RSGs in the open cluster RSGC1 that all are of similar initial mass. The ALMA CO(2-1) line detections allow to retrieve the gas mass-loss rates (Mdot_CO). In contrast to mas… ▽ More [Abridged] Aim: We aim to derive a new mass-loss rate prescription for RSGs that is not afflicted with some uncertainties inherent in preceding studies. Methods: We have observed CO rotational line emission towards a sample of RSGs in the open cluster RSGC1 that all are of similar initial mass. The ALMA CO(2-1) line detections allow to retrieve the gas mass-loss rates (Mdot_CO). In contrast to mass-loss rates derived from the analysis of dust spectral features (Mdot_SED), the data allow a direct determination of the wind velocity and no uncertain dust-to-gas correction factor is needed. Results: Five RSGs in RSGC1 have been detected in CO(2-1). The retrieved Mdot_CO values are systematically lower than Mdot_SED. Although only five RSGs in RSGC1 have been detected, the data allow to propose a new mass-loss rate relation for M-type red supergiants that is dependent on luminosity and initial mass. The new mass-loss rate relation is based on the new Mdot_CO values for the RSGs in RSGC1 and on prior Mdot_SED values for RSGs in 4 clusters, including RSGC1. The new Mdot-prescription yields a good prediction for the mass-loss rate of some well-known Galactic RSGs that are observed in multiple CO rotational lines, including alpha Ori, mu Cep and VX Sgr. However, there are indications that a stronger, potentially eruptive, mass-loss process - different than captured by our new mass-loss rate prescription - is occurring during some fraction of the RSG lifetime. Implementing a lower mass-loss rate in evolution codes for massive stars has important consequences for the nature of their end-state. A reduction of the RSG mass-loss rate implies that quiescent RSG mass loss is not enough to strip a single star's hydrogen-rich envelope. Upon core-collapse such single stars would explode as RSG. △ Less

Submitted 13 November, 2023; v1 submitted 16 March, 2023; originally announced March 2023.

Comments: 26 pages, submitted to Astronomy & Astrophysics, minor revisions already implemented

Journal ref: A&A 681, A17 (2024)

Modeling contact binaries, II. The effect of energy transfer

Authors: Matthias Fabry, Pablo Marchant, Norbert Langer, Hugues Sana

Abstract: Context. It is common for massive stars to engage in binary interaction. In close binaries, the components can enter a contact phase, where both stars overflow their respective Roche lobes simultaneously. While there exist observational constraints on the stellar properties of such systems, the most detailed stellar evolution models that feature a contact phase are not fully reconcilable with thos… ▽ More Context. It is common for massive stars to engage in binary interaction. In close binaries, the components can enter a contact phase, where both stars overflow their respective Roche lobes simultaneously. While there exist observational constraints on the stellar properties of such systems, the most detailed stellar evolution models that feature a contact phase are not fully reconcilable with those measurements. Aims. We aim to consistently model contact phases of binary stars in a 1D stellar evolution code. To this end, we develop the methodology to account for energy transfer in the common contact layers. Methods. We implement an approximative model for energy transfer between the components of a contact binary based on the von Zeipel theorem in the stellar evolution code MESA. We compare structure and evolution models with and without this transfer and analyze the implications for the observable properties of the contact phase. Results. Implementing energy transfer helps eliminating baroclinicity in the common envelope between the components of a contact binary, which, if present, would drive strong thermal flows. We find that accounting for energy transfer in massive contact binaries significantly alters the mass ratio evolution and can extend the lifetime of an unequal mass ratio contact system. △ Less

Submitted 16 March, 2023; originally announced March 2023.

Comments: 11 pages, 5 figures. Accepted to A&A

Journal ref: A&A 672, A175 (2023)

A spectroscopic multiplicity survey of Galactic Wolf-Rayet stars. III. The northern late-type nitrogen-rich sample

Authors: K. Dsilva, T. Shenar, H. Sana, P. Marchant

Abstract: Massive stars are powerful cosmic engines. In the phases immediately preceding core-collapse, massive stars in the Galaxy with $M_i \gtrsim 20$ $M_{\odot}$ may appear as classical Wolf-Rayet (WR) stars. As the final contribution of a homogeneous RV survey, this work constrains the multiplicity properties of northern Galactic late-type nitrogen-rich Wolf-Rayet (WNL) stars. We compare their intrinsi… ▽ More Massive stars are powerful cosmic engines. In the phases immediately preceding core-collapse, massive stars in the Galaxy with $M_i \gtrsim 20$ $M_{\odot}$ may appear as classical Wolf-Rayet (WR) stars. As the final contribution of a homogeneous RV survey, this work constrains the multiplicity properties of northern Galactic late-type nitrogen-rich Wolf-Rayet (WNL) stars. We compare their intrinsic binary fraction and orbital period distribution to the carbon-rich (WC) and early-type nitrogen-rich (WNE) populations from previous works. We obtained high-resolution spectra of the complete magnitude-limited sample of 11 Galactic WNL stars with the Mercator telescope on the island of La Palma. We used cross-correlation to measure relative RVs and flagged binary candidates based on the peak-to-peak RV dispersion. By using Monte Carlo sampling and a Bayesian framework, we computed the three-dimensional likelihood and one-dimensional posteriors for the upper period cut-off, power-law index, and intrinsic binary fraction. Adopting a threshold $C$ of 50 km s$^{-1}$, our Bayesian analysis produced an intrinsic fraction of $0.42\substack{+0.15 \\ -0.17}$ for the parent WNL population alongside distributions for the power-law index and the orbital periods. The observed period distribution of Galactic WN and WC binaries from the literature is in agreement with what is found. The period distribution of Galactic WN binaries peaks at $P{\sim}1$-$10$d and that of the WC population at $P{\sim}5000\,$d. This shift cannot be reconciled by orbital evolution due to mass loss or mass transfer. At long periods, the evolutionary sequence O($\xrightarrow{}$LBV)$\xrightarrow{}$WN$\xrightarrow{}$WC seems feasible. The high frequency of short-period WN binaries compared to WC binaries suggests that they either tend to merge or that the WN components in these binaries rarely evolve into WC stars in the Galaxy. △ Less

Submitted 13 December, 2022; originally announced December 2022.

Comments: Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 674, A88 (2023)

The initial spin distribution of B-type stars revealed by the split main sequences of young star clusters

Authors: Chen Wang, Ben Hastings, Abel Schootemeijer, Norbert Langer, Selma E. de Mink, Julia Bodensteiner, Antonino Milone, Stephen Justham, Pablo Marchant

Abstract: Spectroscopic observations of stars in young open clusters have revealed evidence for a dichotomous distribution of stellar rotational velocities, with 10-30% of stars rotating slowly and the remaining 70-90% rotating fairly rapidly. At the same time, high-precision multiband photometry of young star clusters shows a split main sequence band, which is again interpreted as due to a spin dichotomy.… ▽ More Spectroscopic observations of stars in young open clusters have revealed evidence for a dichotomous distribution of stellar rotational velocities, with 10-30% of stars rotating slowly and the remaining 70-90% rotating fairly rapidly. At the same time, high-precision multiband photometry of young star clusters shows a split main sequence band, which is again interpreted as due to a spin dichotomy. Recent papers suggest that extreme rotation is required to retrieve the photometric split. Our new grids of MESA models and the prevalent SYCLIST models show, however, that initial slow (0-35% of the linear Keplerian rotation velocities) and intermediate (50-65% of the Keplerian rotation velocities) rotation are adequate to explain the photometric split. These values are consistent with the recent spectroscopic measurements of cluster and field stars, and are likely to reflect the birth spin distributions of upper main-sequence stars. A fraction of the initially faster-rotating stars may be able to reach near-critical rotation at the end of their main-sequence evolution and produce Be stars in the turn-off region of young star clusters. However, we find that the presence of Be stars up to two magnitudes below the cluster turnoff advocates for a crucial role of binary interaction in creating Be stars. We argue that surface chemical composition measurements may help distinguish these two Be star formation channels. While only the most rapidly rotating, and therefore nitrogen-enriched, single stars can evolve into Be stars, slow pre-mass-transfer rotation and inefficient accretion allows for mild or no enrichment even in critically rotating accretion-induced Be stars. Our results shed new light on the origin of the spin distribution of young and evolved B-type main sequence stars. △ Less

Submitted 7 December, 2022; v1 submitted 28 November, 2022; originally announced November 2022.

Comments: 18 pages, 14 figures. Accepted for publication in A&A

Journal ref: A&A 670, A43 (2023)

Detection of single-degenerate massive binaries with Gaia: The impact of blue supergiants, triples, mass precision, and high-precision parallax requirements

Authors: S. Janssens, T. Shenar, H. Sana, P. Marchant

Abstract: Hundreds of black holes with massive main-sequence companions (OB+BHs) might be identified from Gaia astrometry with the Astrometric Mass-Ratio Function (AMRF). We investigate the impact of blue supergiant companions (BSG) instead of dwarfs and the presence of additional companions in the system that are unresolved by Gaia on the astrometric identification of OB+BHs. We also explore how accurate t… ▽ More Hundreds of black holes with massive main-sequence companions (OB+BHs) might be identified from Gaia astrometry with the Astrometric Mass-Ratio Function (AMRF). We investigate the impact of blue supergiant companions (BSG) instead of dwarfs and the presence of additional companions in the system that are unresolved by Gaia on the astrometric identification of OB+BHs. We also explore how accurate the primary mass needs to be constrained. Moreover, we assess how the high-precision publishing constraints of astrometric binary orbits in the latest Gaia data release DR3 impact the detection of OB+BHs. We establish a BSG mass-magnitude relation and compute BSG AMRF curves. From a mock population of non- and single-degenerate massive binaries and triples with OB or BSG primaries, we asses the false-positive identification fraction and the effect of the BSG AMRF curves. We compare the number of stars with astrometric DR3 orbits in the second Alma Luminous Star catalogue (ALSII) with new predictions on the OB+BH detection using the conservative DR3 publishing criteria. BSG primaries and triples are not expected to impact the false-positive identification fractions significantly. However, if the evolutionary stage of the primary star is unknown, the usage of the BSG curves is still recommended. This also reduces the OB+BH identification fraction significantly. The primary mass does not need to be known to benefit from the high identification fraction of OB+BHs, while keeping the fraction of false-positives low. We find no OB+BH candidates in the ALS II among the DR3 orbits. This null-detection cannot be attributed to the underlying BH-formation scenario, but rather to the stringent DR3 selection criteria. To infer the BH-formation scenario with Gaia, we propose that the constraint on the relative parallax precision in DR4 should be 95% less conservative than the DR3 criterion.(Abridged) △ Less

Submitted 10 November, 2022; originally announced November 2022.

Comments: 15 pages (+ 5 appendix), 10 figures, 3 tables, submitted to A&A

Journal ref: A&A 670, A79 (2023)

The Nature of Unseen Companions in Massive Single-Line Spectroscopic Binaries

Authors: Hugues Sana, Michael Abdul-Masih, Gareth Banyard, Julia Bodensteiner, Dominic M. Bowman, Karan Dsilva, C. Eldridge, Matthias Fabry, Abigail J. Frost, Calum Hawcroft, Soetkin Janssens, Laurent Mahy, Pablo Marchant, Norbert Langer, Timothy Van Reeth, Koushik Sen, Tomer Shenar

Abstract: Massive stars are predominantly found in binaries and higher order multiples. While the period and eccentricity distributions of OB stars are now well established across different metallicity regimes, the determination of mass-ratios has been mostly limited to double-lined spectroscopic binaries. As a consequence, the mass-ratio distribution remains subject to significant uncertainties. Open quest… ▽ More Massive stars are predominantly found in binaries and higher order multiples. While the period and eccentricity distributions of OB stars are now well established across different metallicity regimes, the determination of mass-ratios has been mostly limited to double-lined spectroscopic binaries. As a consequence, the mass-ratio distribution remains subject to significant uncertainties. Open questions include the shape and extent of the companion mass-function towards its low-mass end and the nature of undetected companions in single-lined spectroscopic binaries. In this contribution, we present the results of a large and systematic analysis of a sample of over 80 single-lined O-type spectroscopic binaries (SB1s) in the Milky Way and in the Large Magellanic Cloud (LMC). We report on the developed methodology, the constraints obtained on the nature of SB1 companions, the distribution of O star mass-ratios at LMC metallicity and the occurrence of quiescent OB+black hole binaries. △ Less

Submitted 13 September, 2022; originally announced September 2022.

Comments: To appear in the proceedings of IAUS361: Massive stars near and far; 6 pages, 1 figure

A synthetic population of Wolf-Rayet stars in the LMC based on detailed single and binary star evolution models

Authors: D. Pauli, N. Langer, D. R. Aguilera-Dena, C. Wang, P. Marchant

Abstract: Without doubt, mass transfer in close binary systems contributes to the populations of Wolf-Rayet (WR) stars in the Milky Way and the Magellanic Clouds. However, the binary formation channel is so far not well explored. We want to remedy this by exploring large grids of detailed binary and single star evolution models computed with the publicly available MESA code, for a metallicity appropriate fo… ▽ More Without doubt, mass transfer in close binary systems contributes to the populations of Wolf-Rayet (WR) stars in the Milky Way and the Magellanic Clouds. However, the binary formation channel is so far not well explored. We want to remedy this by exploring large grids of detailed binary and single star evolution models computed with the publicly available MESA code, for a metallicity appropriate for the Large Magellanic Cloud (LMC). The binary models are calculated through Roche-lobe overflow and mass transfer, until the initially more massive star exhausts helium in its core. We distinguish models of WR and helium stars based on the estimated stellar wind optical depth. We use these models to build a synthetic WR population, assuming constant star formation. Our models can reproduce the WR population of the LMC to significant detail, including the number and luminosity functions of the main WR subtypes. We find that for binary fractions of 100% (50%), all LMC WR stars below $10^6\,L_{\odot}$ ($10^{5.7}\,L_{\odot}$) are stripped binary mass donors. We also identify several insightful mismatches. With a single star fraction of 50\%, our models produce too many yellow supergiants, calling either for a larger initial binary fraction, or for enhanced mass-loss near the Humphreys-Davidson limit. Our models predict more long-period WR binaries than observed, arguably due to an observational bias towards short periods. Our models also underpredict the shortest-period WR binaries, which may have implications for understanding the progenitors of double black hole mergers. The fraction of binary produced WR stars may be larger than often assumed, and outline the risk to mis-calibrate stellar physics when only single star models are used to reproduce the observed WR stars. △ Less

Submitted 22 August, 2022; originally announced August 2022.

Comments: 15 pages + 13 pages appendix, 14 figures, 2 tables

Journal ref: A&A 667, A58 (2022)

Modules for Experiments in Stellar Astrophysics (MESA): Time-Dependent Convection, Energy Conservation, Automatic Differentiation, and Infrastructure

Authors: Adam S. Jermyn, Evan B. Bauer, Josiah Schwab, R. Farmer, Warrick H. Ball, Earl P. Bellinger, Aaron Dotter, Meridith Joyce, Pablo Marchant, Joey S. G. Mombarg, William M. Wolf, Tin Long Sunny Wong, Giulia C. Cinquegrana, Eoin Farrell, R. Smolec, Anne Thoul, Matteo Cantiello, Falk Herwig, Odette Toloza, Lars Bildsten, Richard H. D. Townsend, F. X. Timmes

Abstract: We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). The new auto_diff module implements automatic differentiation in MESA, an enabling capability that alleviates the need for hard-coded analytic expressions or finite difference approximations. We significantly enhance the treatment of the growth and decay of convection in MES… ▽ More We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). The new auto_diff module implements automatic differentiation in MESA, an enabling capability that alleviates the need for hard-coded analytic expressions or finite difference approximations. We significantly enhance the treatment of the growth and decay of convection in MESA with a new model for time-dependent convection, which is particularly important during late-stage nuclear burning in massive stars and electron degenerate ignition events. We strengthen MESA's implementation of the equation of state, and we quantify continued improvements to energy accounting and solver accuracy through a discussion of different energy equation features and enhancements. To improve the modeling of stars in MESA we describe key updates to the treatment of stellar atmospheres, molecular opacities, Compton opacities, conductive opacities, element diffusion coefficients, and nuclear reaction rates. We introduce treatments of starspots, an important consideration for low-mass stars, and modifications for superadiabatic convection in radiation-dominated regions. We describe new approaches for increasing the efficiency of calculating monochromatic opacities and radiative levitation, and for increasing the efficiency of evolving the late stages of massive stars with a new operator split nuclear burning mode. We close by discussing major updates to MESA's software infrastructure that enhance source code development and community engagement. △ Less

Submitted 30 December, 2022; v1 submitted 7 August, 2022; originally announced August 2022.

Comments: 50 pages, 29 figures; Accepted to ApJS

Constraining the overcontact phase in massive binary evolution -- II. Period stability of known O+O overcontact systems

Authors: Michael Abdul-Masih, Ana Escorza, Athira Menon, Laurent Mahy, Pablo Marchant

Abstract: Given that mergers are often invoked to explain many exotic phenomena in massive star evolution, understanding the evolutionary phase directly preceding a merger, the overcontact phase, is of crucial importance. Despite its importance, large uncertainties exist in our understanding of the evolution of massive overcontact binaries. We aim to provide robust observational constraints on the future dy… ▽ More Given that mergers are often invoked to explain many exotic phenomena in massive star evolution, understanding the evolutionary phase directly preceding a merger, the overcontact phase, is of crucial importance. Despite its importance, large uncertainties exist in our understanding of the evolution of massive overcontact binaries. We aim to provide robust observational constraints on the future dynamical evolution of massive overcontact systems by measuring the rate at which the periods change for a sample of six such objects. Furthermore, we aim to investigate whether the periods of unequal mass systems show higher rates of change than their equal mass counterparts as theoretical models predict. Using archival photometric data from various ground- and space-based missions covering up to ~40 years, we measure the periods of each system over several smaller time spans. We then fit a linear regression through the measured periods to determine the rate at which the period is changing over the entire data set. We find that all of the stars in our sample have very small period changes and that there does not seem to be a correlation with the mass ratio. This implies that the orbital periods for these systems are stable on the nuclear timescale, and that the unequal mass systems may not equalize as expected. When comparing our results with population synthesis distributions, we find large discrepancies between the expected mass ratios and period stabilities. We find that these discrepancies can be mitigated to a degree by removing systems with shorter initial periods, suggesting that the observed sample of overcontact systems may originate from binary systems with longer initial orbital periods. △ Less

Submitted 2 August, 2022; originally announced August 2022.

Comments: 11 pages, 6 (+5) Figures, accepted for publication in A&A

Journal ref: A&A 666, A18 (2022)

Identifying quiescent compact objects in massive Galactic single-lined spectroscopic binaries

Authors: L. Mahy, H. Sana, T. Shenar, K. Sen, N. Langer, P. Marchant, M. Abdul-Masih, G. Banyard, J. Bodensteiner, D. M. Bowman, K. Dsilva, M. Fabry, C. Hawcroft, S. Janssens, T. Van Reeth, C. Eldridge

Abstract: Aims. To search for these rare objects, we study 32 Galactic O-type stars that were reported as SB1s in the literature. In our sample we include Cyg X-1, which is known to host an accreting stellar-mass BH, and HD 74194, a supergiant fast X-ray transient, in order to validate our methodology. The final goal is to characterise the nature of the unseen companions to determine if they are MS stars, s… ▽ More Aims. To search for these rare objects, we study 32 Galactic O-type stars that were reported as SB1s in the literature. In our sample we include Cyg X-1, which is known to host an accreting stellar-mass BH, and HD 74194, a supergiant fast X-ray transient, in order to validate our methodology. The final goal is to characterise the nature of the unseen companions to determine if they are MS stars, stripped helium stars, triples, or compact objects such as neutron stars or stellar-mass BHs. Methods. After measuring radial velocities and deriving orbital solutions for all the systems in our sample, we performed spectral disentangling to extract putative signatures of faint secondary companions from the composite spectra. We derived stellar parameters for the visible stars and estimated the mass ranges of the secondary stars using the binary mass function. Variability observed in the photometric TESS light curves was also searched for indications of the presence of putative companions, degenerate or not. Results. In 17 of the 32 systems reported as SB1s, we extract secondary signatures, down to mass ratios of ~0.15. For the 17 newly detected double-lined spectroscopic binaries (SB2s), we derive physical properties of the individual components and discuss why they have not been detected as such before. Among the remaining systems, we identify nine systems with possible NS or low-mass MS companions. For Cyg X-1 and HD 130298, we are not able to extract any signatures for the companions, and the minimum masses of their companions are estimated to be about 7Msun. Our simulations show that secondaries with such a mass should be detectable from our dataset, no matter their nature: MS stars, stripped helium stars or even triples. While this is expected for Cyg X-1, confirming our methodology, our simulations also strongly suggest that HD 130298 could be another candidate to host a stellar-mass BH. △ Less

Submitted 15 July, 2022; originally announced July 2022.

Comments: 20 pages (31 pages of Appendix), 19 figures, A&A in press. The abstract is shorter than in the paper, and some figures, and tables have been cut in the paper

An X-ray quiet black hole born with a negligible kick in a massive binary within the Large Magellanic Cloud

Authors: Tomer Shenar, Hugues Sana, Laurent Mahy, Kareem El-Badry, Pablo Marchant, Norbert Langer, Calum Hawcroft, Matthias Fabry, Koushik Sen, Leonardo A. Almeida, Michael Abdul-Masih, Julia Bodensteiner, Paul A. Crowther, Mark Gieles, Mariusz Gromadzki, Vincent Henault-Brunet, Artemio Herrero, Alex de Koter, Patryk Iwanek, Szymon Kozłowski, Daniel J. Lennon, Jesus Maız Apellaniz, Przemysław Mroz, Anthony F. J. Moffat, Annachiara Picco , et al. (13 additional authors not shown)

Abstract: Stellar-mass black holes are the final remnants of stars born with more than 15 solar masses. Billions are expected to reside in the Local Group, yet only few are known, mostly detected through X-rays emitted as they accrete material from a companion star. Here, we report on VFTS 243: a massive X-ray faint binary in the Large Magellanic Cloud. With an orbital period of 10.4-d, it comprises an O-ty… ▽ More Stellar-mass black holes are the final remnants of stars born with more than 15 solar masses. Billions are expected to reside in the Local Group, yet only few are known, mostly detected through X-rays emitted as they accrete material from a companion star. Here, we report on VFTS 243: a massive X-ray faint binary in the Large Magellanic Cloud. With an orbital period of 10.4-d, it comprises an O-type star of 25 solar masses and an unseen companion of at least nine solar masses. Our spectral analysis excludes a non-degenerate companion at a 5-sigma confidence level. The minimum companion mass implies that it is a black hole. No other X-ray quiet black hole is unambiguously known outside our Galaxy. The (near-)circular orbit and kinematics of VFTS 243 imply that the collapse of the progenitor into a black hole was associated with little or no ejected material or black-hole kick. Identifying such unique binaries substantially impacts the predicted rates of gravitational-wave detections and properties of core-collapse supernovae across the Cosmos. △ Less

Submitted 15 July, 2022; originally announced July 2022.

Comments: Accepted to Nature Astronomy, 64 pages, 15 figures, 2 tables; ESO press release: https://www.eso.org/public/news/eso2210/; Nat Asr paper URL: https://www.nature.com/articles/s41550-022-01730-y

The Tarantula Massive Binary Monitoring VI: Characterisation of hidden companions in 51 single-lined O-type binaries, a flat mass-ratio distribution, and black-hole binary candidates

Authors: Tomer Shenar, Hugues Sana, Laurent Mahy, Jesus Maiz Apellaniz, Paul A. Crowther, Mariusz Gromadzki, Artemio Herrero, Norbert Langer, Pablo Marchant, Fabian R. N. Schneider, Koushik Sen, Igor Soszynski, S. Toonen

Abstract: We aim to hunt for massive binaries hosting a black hole companion (OB+BH) and establish the natal mass-ratio distribution of massive stars at the subsolar metallicity environment of the Large Magellanic Cloud (LMC). We use the shift-and-add grid disentangling technique to characterize the hidden companions in 51 SB1 O-type and evolved B-type binaries in the LMC monitored in the framework of the T… ▽ More We aim to hunt for massive binaries hosting a black hole companion (OB+BH) and establish the natal mass-ratio distribution of massive stars at the subsolar metallicity environment of the Large Magellanic Cloud (LMC). We use the shift-and-add grid disentangling technique to characterize the hidden companions in 51 SB1 O-type and evolved B-type binaries in the LMC monitored in the framework of the Tarantula Massive Binary Monitoring (TMBM). Out of the 51 SB1 systems, 43 (84%) are found to have non-degenerate stellar companions, of which 28 are confident detections, and 15 are less certain (SB1: or SB2:). Of these 43 targets, one is found to be a triple (VFTS 64), and two are found to be quadruples (VFTS 120, 702). The remaining eight targets (16%) retain an SB1 classification. Aside from the unambiguous case of VFTS 243, analysed in detailed in a separate paper, we identify two additional OB+BH candidates: VFTS 514 and VFTS 779. Additional black holes may be present in the sample but at a lower probability. Our study firmly establishes a virtually flat natal mass-ratio distribution for O-type stars at LMC metallicity, covering the entire mass-ratio range (0.05 < q < 1) and periods in the range 0 < log P < 3 [d]. The nature of the OB+BH candidates should be verified through future monitoring, but the frequency of OB+BH candidates is generally in line with recent predictions at LMC metallicity. △ Less

Submitted 5 October, 2022; v1 submitted 15 July, 2022; originally announced July 2022.

Comments: 41 pages (14 main article + 27 appendix), accepted to A&A

Journal ref: A&A 665, A148 (2022)

Detecting Stripped Stars While Searching for Quiescent Black Holes

Authors: J. Bodensteiner, M. Heida, M. Abdul-Masih, D. Baade, G. Banyard, D. M. Bowman, M. Fabry, A. Frost, L. Mahy, P. Marchant, A. Mérand, M. Reggiani, Th. Rivinius, H. Sana, F. Selman, T. Shenar

Abstract: While the number of stellar-mass black holes detected in X-rays or as gravitational wave sources is steadily increasing, the known population remains orders of magnitude smaller than predicted by stellar evolution theory. A significant fraction of stellar-mass black holes is expected to hide in X-ray-quiet binaries where they are paired with a "normal" star. Although a handful of such quiescent bl… ▽ More While the number of stellar-mass black holes detected in X-rays or as gravitational wave sources is steadily increasing, the known population remains orders of magnitude smaller than predicted by stellar evolution theory. A significant fraction of stellar-mass black holes is expected to hide in X-ray-quiet binaries where they are paired with a "normal" star. Although a handful of such quiescent black hole candidates have been proposed, the majority have been challenged by follow-up investigations. A confusion that emerged recently concerns binary systems that appear to contain a normal B-type star with an unseen companion, believed to be a black hole. On closer inspection, some of these seemingly normal B-type stars instead turn out to be stars stripped of most of their mass through an interaction with their binary companion, which in at least two cases is a rapidly rotating star rather than a compact object. These contaminants in the search for quiescent black holes are themselves extremely interesting objects as they represent a rare phase of binary evolution, and should be given special attention when searching for binaries hosting black holes in large spectroscopic studies. △ Less

Submitted 1 July, 2022; originally announced July 2022.

Comments: 7 pages, 5 figures

Journal ref: The Messenger, 2022, vol. 186, p. 3-9

A spectroscopic multiplicity survey of Galactic Wolf-Rayet stars: II. The northern WNE sequence

Authors: K. Dsilva, T. Shenar, H. Sana, P. Marchant

Abstract: Most massive stars reside in multiple systems that will interact over the course of their lifetime. Classical Wolf-Rayet (WR) stars represent the final end stages of stellar evolution at the upper-mass end. As part of a homogeneous, magnitude-limited ($V\leq12$) spectroscopic survey of northern Galactic WR stars, this paper aims to establish the observed and intrinsic multiplicity properties of th… ▽ More Most massive stars reside in multiple systems that will interact over the course of their lifetime. Classical Wolf-Rayet (WR) stars represent the final end stages of stellar evolution at the upper-mass end. As part of a homogeneous, magnitude-limited ($V\leq12$) spectroscopic survey of northern Galactic WR stars, this paper aims to establish the observed and intrinsic multiplicity properties of the early-type nitrogen-rich WR population (WNE). We obtained high-resolution spectroscopic time series of the complete magnitude-limited sample of 16 WNE stars observable with the 1.2 m Mercator telescope at La Palma, typically providing a time base of about two to eight years. We measured relative radial velocities (RVs) using cross-correlation and used RV variations to flag binary candidates. Adopting a peak-to-peak RV variability threshold of 50 km/s as a criterion found an observed multiplicity fraction of 0.44$\pm$0.12. Using an updated Monte Carlo method with a Bayesian framework, we calculated the three-dimensional likelihood for the intrinsic binary fraction, the maximum period, and the power-law index for the period distribution for the WNE population. We also re-derived multiplicity parameters for the Galactic WC population. We found an intrinsic multiplicity fraction of $0.56\substack{+0.20 \\ -0.15}$ for the parent WNE population. For the Galactic WC population, we re-derive an intrinsic multiplicity fraction of $0.96\substack{+0.04 \\ -0.22}$. The derived multiplicity parameters for the WNE population are quite similar to those derived for main-sequence O binaries but differ from those of the WC population. The significant shift in the WC period distribution towards longer periods is too large to be explained via expansion of the orbit due to stellar winds, and we discuss possible implications of our results. △ Less

Submitted 26 April, 2022; originally announced April 2022.

Comments: 17 pages, 12 figures

Journal ref: A&A 664, A93 (2022)

Active Learning for Computationally Efficient Distribution of Binary Evolution Simulations

Authors: Kyle Akira Rocha, Jeff J. Andrews, Christopher P. L. Berry, Zoheyr Doctor, Aggelos K. Katsaggelos, Juan Gabriel Serra Pérez, Pablo Marchant, Vicky Kalogera, Scott Coughlin, Simone S. Bavera, Aaron Dotter, Tassos Fragos, Konstantinos Kovlakas, Devina Misra, Zepei Xing, Emmanouil Zapartas

Abstract: Binary stars undergo a variety of interactions and evolutionary phases, critical for predicting and explaining observed properties. Binary population synthesis with full stellar-structure and evolution simulations are computationally expensive requiring a large number of mass-transfer sequences. The recently developed binary population synthesis code POSYDON incorporates grids of MESA binary star… ▽ More Binary stars undergo a variety of interactions and evolutionary phases, critical for predicting and explaining observed properties. Binary population synthesis with full stellar-structure and evolution simulations are computationally expensive requiring a large number of mass-transfer sequences. The recently developed binary population synthesis code POSYDON incorporates grids of MESA binary star simulations which are then interpolated to model large-scale populations of massive binaries. The traditional method of computing a high-density rectilinear grid of simulations is not scalable for higher-dimension grids, accounting for a range of metallicities, rotation, and eccentricity. We present a new active learning algorithm, psy-cris, which uses machine learning in the data-gathering process to adaptively and iteratively select targeted simulations to run, resulting in a custom, high-performance training set. We test psy-cris on a toy problem and find the resulting training sets require fewer simulations for accurate classification and regression than either regular or randomly sampled grids. We further apply psy-cris to the target problem of building a dynamic grid of MESA simulations, and we demonstrate that, even without fine tuning, a simulation set of only $\sim 1/4$ the size of a rectilinear grid is sufficient to achieve the same classification accuracy. We anticipate further gains when algorithmic parameters are optimized for the targeted application. We find that optimizing for classification only may lead to performance losses in regression, and vice versa. Lowering the computational cost of producing grids will enable future versions of POSYDON to cover more input parameters while preserving interpolation accuracies. △ Less

Submitted 16 September, 2022; v1 submitted 30 March, 2022; originally announced March 2022.

Comments: 21 pages, 10 figures, ApJ in press

Journal ref: Astrophysical Journal; 938(1):64(15); 2022

Modeling overcontact binaries, I. The effect of tidal deformation

Authors: M. Fabry, P. Marchant, H. Sana

Abstract: In the realm of massive stars, strong binary interaction is commonplace. One extreme case are overcontact systems, which is expected to be part of the evolution of all stars evolving towards a merger, and is hypothesized to play a role in the formation of binary black holes. However, important simplifications are made to model the evolution of overcontact binaries, namely that tidal or rotational… ▽ More In the realm of massive stars, strong binary interaction is commonplace. One extreme case are overcontact systems, which is expected to be part of the evolution of all stars evolving towards a merger, and is hypothesized to play a role in the formation of binary black holes. However, important simplifications are made to model the evolution of overcontact binaries, namely that tidal or rotational deformation is frequently ignored. Yet, both observation and theory show that overcontact stars are heavily tidally deformed, leaving a potentially important effect on the outer layers unaccounted for in models. In this work we develop the methodology to represented tidally deformed stars. Using numerical methods, we compute the structure correction factors to the 1D spherical stellar structure equations due to the binary Roche potential, and compare them to existing results and the structure corrections of single rotating stars. We implement the new structure correction factors into the stellar evolution code MESA and explore several case studies. We compare the differences between our simulations when no rotation is included, when we treat rotation using single star corrections (i.e. only accounting for centrifugal deformation) or when we use tidal deformation. We find that ignoring rotation in deformed detached eclipsing binaries can produce a radius discrepancy of up to 5%. The difference between tidal and single star centrifugal distortion models is more benign at 1%, showing that single rotating star models are a suitable approximation of tidally deformed stars in a binary system. In overcontact configurations, we find a similar 5% variation in surface properties as a result of tidal distortion with respect to non-rotating models, showing that it is inappropriate to model binary stars that fill their Roche lobe significantly, as non-rotating. △ Less

Submitted 17 February, 2022; originally announced February 2022.

Comments: 12 pages, 9 figures, accepted to appear in Astronomy & Astrophysics

Journal ref: A&A 661, A123 (2022)

Stellar mergers as the origin of the blue main-sequence band in young star clusters

Authors: Chen Wang, Norbert Langer, Abel Schootemeijer, Antonino Milone, Ben Hastings, Xiao-Tian Xu, Julia Bodensteiner, Hugues Sana, Norberto Castro, D. J. Lennon, Pablo Marchant, A. de Koter, Selma E. de Mink

Abstract: Recent high-quality Hubble Space Telescope (HST) photometry shows that the main sequences (MS) stars of young star clusters form two discrete components in the color-magnitude diagram (CMD). Based on their distribution in the CMD, we show that stars of the blue MS component can be understood as slow rotators originating from stellar mergers. We derive the masses of the blue MS stars, and find that… ▽ More Recent high-quality Hubble Space Telescope (HST) photometry shows that the main sequences (MS) stars of young star clusters form two discrete components in the color-magnitude diagram (CMD). Based on their distribution in the CMD, we show that stars of the blue MS component can be understood as slow rotators originating from stellar mergers. We derive the masses of the blue MS stars, and find that they follow a nearly flat mass function, which supports their unusual formation path. Our results imply that the cluster stars gain their mass in two different ways, by disk accretion leading to rapid rotation, contributing to the red MS, or by binary merger leading to slow rotation and populating the blue MS. We also derive the approximate merger time of the individual stars of the blue MS component, and find a strong early peak in the merger rate, with a lower level merger activity prevailing for tens of Myr. This supports recent binary formation models, and explains new velocity dispersion measurements for members of young star clusters. Our findings shed new light on the origin of the bi-modal mass, spin, and magnetic field distributions of main-sequence stars. △ Less

Submitted 11 February, 2022; originally announced February 2022.

Comments: Published in Nature Astronomy online. Link to the paper:https://rdcu.be/cGLpH

Uncovering astrometric black hole binaries with massive main-sequence companions with Gaia

Authors: S. Janssens, T. Shenar, H. Sana, S. Faigler, N. Langer, P. Marchant, T. Mazeh, C. Schürmann, S. Shahaf

Abstract: The hunt for compact objects is on. Rarely seen massive binaries with a compact object are a crucial phase in the evolution towards compact object mergers. In Gaia data release 3 (DR3), the first Gaia astrometric orbital solutions for binary sources will become available. We investigate how many black holes (BH) with massive main-sequence dwarf companions (OB+BH binaries) are expected to be detect… ▽ More The hunt for compact objects is on. Rarely seen massive binaries with a compact object are a crucial phase in the evolution towards compact object mergers. In Gaia data release 3 (DR3), the first Gaia astrometric orbital solutions for binary sources will become available. We investigate how many black holes (BH) with massive main-sequence dwarf companions (OB+BH binaries) are expected to be detected as binaries in DR3 and at the end of the nominal 5-yr mission (DR4). We estimate the fraction of identifiable OB+BH binaries and discuss the distributions of the masses of both components and the orbital periods. We study the impact of different BH-formation scenarios. Using tailored models for the massive star population, which assume a direct collapse and no kick upon BH formation (the fiducial case), we estimate the fraction of OB+BH systems that Gaia will detect as binaries. A distance distribution according to that of the second Alma Luminous Star catalogue (ALSII) is assumed. We investigate how many of the systems detected as binaries are identifiable as OB+BH binaries, using a method based on astrometric data. In the fiducial case we conservatively estimate that 77% of the OB+BH binaries in ALSII will be detected as binaries in DR3, of which 89% are identifiable as OB+BH binaries. This leads to a total of around 190 OB+BH binaries, a 20-fold increase in the known sample of OB+BH binaries, covering an uncharted parameter space of long-period binaries. The size and properties of the identifiable OB+BH population will contain crucial observational constraints to improve our understanding of BH formation. In DR4, the detected fraction will increase to 85%, of which 82% will be identifiable. Hence, an additional ~5 systems could be identified, which are expected to have either very short or long periods. The fractions become smaller for different BH-formation scenarios. (truncated) △ Less

Submitted 11 November, 2021; originally announced November 2021.

Comments: 22 pages, 18 figures, 2 tables, accepted for publication in A&A

Journal ref: A&A 658, A129 (2022)

Detailed models of interacting short-period massive binary stars

Authors: K. Sen, N. Langer, P. Marchant, A. Menon, S. E. de Mink, A. Schootemeijer, C. Schürmann, L. Mahy, B. Hastings, K. Nathaniel, H. Sana, C. Wang, X. T. Xu

Abstract: About a quarter of massive binary stars undergo mass transfer while both stars burn hydrogen at their cores, first on the thermal and then on the nuclear timescale. The nuclear timescale mass transfer leads to observable counterparts: the semi-detached so-called massive Algol binaries. However, comprehensive model predictions for these systems are sparse. We study them using a large grid of ~10,00… ▽ More About a quarter of massive binary stars undergo mass transfer while both stars burn hydrogen at their cores, first on the thermal and then on the nuclear timescale. The nuclear timescale mass transfer leads to observable counterparts: the semi-detached so-called massive Algol binaries. However, comprehensive model predictions for these systems are sparse. We study them using a large grid of ~10,000 detailed binary evolution models calculated with the stellar evolution code MESA, covering initial donor masses between 10-40 M$_{\odot}$ and initial orbital periods above 1.4 d, at a metallicity suitable for the Large Magellanic Cloud (LMC). Our models imply ~30, or ~3% of the ~1,000 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 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 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/helium+OB star binaries. 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. △ Less

Submitted 9 December, 2021; v1 submitted 5 November, 2021; originally announced November 2021.

Comments: Accepted for publication in A&A (December 6, 2021)

Journal ref: A&A 659, A98 (2022)

Binary Black Hole Formation with Detailed Modeling: Stable Mass Transfer Leads to Lower Merger Rates

Authors: Monica Gallegos-Garcia, Christopher P L Berry, Pablo Marchant, Vicky Kalogera

Abstract: Rapid binary population synthesis codes are often used to investigate the evolution of compact-object binaries. They typically rely on analytical fits of single-star evolutionary tracks and parameterized models for interactive phases of evolution (e.g., mass-transfer on thermal timescale, determination of dynamical instability, and common envelope) that are crucial to predict the fate of binaries.… ▽ More Rapid binary population synthesis codes are often used to investigate the evolution of compact-object binaries. They typically rely on analytical fits of single-star evolutionary tracks and parameterized models for interactive phases of evolution (e.g., mass-transfer on thermal timescale, determination of dynamical instability, and common envelope) that are crucial to predict the fate of binaries. These processes can be more carefully implemented in stellar structure and evolution codes such as MESA. To assess the impact of such improvements, we compare binary black hole mergers as predicted in models with the rapid binary population synthesis code COSMIC to models ran with MESA simulations through mass transfer and common-envelope treatment. We find that results significantly differ in terms of formation paths, the orbital periods and mass ratios of merging binary black holes, and consequently merger rates. While common-envelope evolution is the dominant formation channel in COSMIC, stable mass transfer dominates in our MESA models. Depending upon the black hole donor mass, and mass-transfer and common-envelope physics, at sub-solar metallicity COSMIC overproduces the number of binary black hole mergers by factors of 2--35 with a significant fraction of them having merger times orders of magnitude shorter than the binary black holes formed when using detailed MESA models. Therefore we find that some binary black hole merger rate predictions from rapid population syntheses of isolated binaries may be overestimated by factors of ~5--500. We conclude that the interpretation of gravitational-wave observations requires the use of detailed treatment of these interactive binary phases. △ Less

Submitted 5 November, 2021; v1 submitted 12 July, 2021; originally announced July 2021.

Comments: 17 pages (13 main text), 5 figures, ApJ, in press

Probing the progenitors of spinning binary black-hole mergers with long gamma-ray bursts

Authors: Simone S. Bavera, Tassos Fragos, Emmanouil Zapartas, Enrico Ramirez-Ruiz, Pablo Marchant, Luke Z. Kelley, Michael Zevin, Jeff J. Andrews, Scott Coughlin, Aaron Dotter, Konstantinos Kovlakas, Devina Misra, Juan G. Serra-Perez, Ying Qin, Kyle A. Rocha, Jaime Román-Garza, Nam H. Tran, Zepei Xing

Abstract: Long-duration gamma-ray bursts are thought to be associated with the core-collapse of massive, rapidly spinning stars and the formation of black holes. However, efficient angular momentum transport in stellar interiors, currently supported by asteroseismic and gravitational-wave constraints, leads to predominantly slowly-spinning stellar cores. Here, we report on binary stellar evolution and popul… ▽ More Long-duration gamma-ray bursts are thought to be associated with the core-collapse of massive, rapidly spinning stars and the formation of black holes. However, efficient angular momentum transport in stellar interiors, currently supported by asteroseismic and gravitational-wave constraints, leads to predominantly slowly-spinning stellar cores. Here, we report on binary stellar evolution and population synthesis calculations, showing that tidal interactions in close binaries not only can explain the observed sub-population of spinning, merging binary black holes but also lead to long gamma-ray bursts at the time of black-hole formation. Given our model calibration against the distribution of isotropic-equivalent energies of luminous long gamma-ray bursts, we find that ~10% of the GWTC-2 reported binary black holes had a luminous long gamma-ray burst associated with their formation, with GW190517 and GW190719 having a probability of ~85% and ~60%, respectively, being among them. Moreover, given an assumption about their average beaming fraction, our model predicts the rate density of long gamma-ray bursts, as a function of redshift, originating from this channel. For a constant beaming fraction $f_\mathrm{B}\sim 0.05$ our model predicts a rate density comparable to the observed one, throughout the redshift range, while, at redshift $z \in [0,2.5]$, a tentative comparison with the metallicity distribution of observed LGRB host galaxies implies that between 20% to 85% of the observed long gamma-ray bursts may originate from progenitors of merging binary black holes. The proposed link between a potentially significant fraction of observed, luminous long gamma-ray bursts and the progenitors of spinning binary black-hole mergers allows us to probe the latter well outside the horizon of current-generation gravitational wave observatories, and out to cosmological distances. △ Less

Submitted 3 December, 2021; v1 submitted 30 June, 2021; originally announced June 2021.

Comments: Accepted for publication in A&A Letters, 12 pages, 6 figures, 1 table

Resolving the dynamical mass tension of the massive binary 9 Sagittarii

Authors: M. Fabry, C. Hawcroft, A. J. Frost, L. Mahy, P. Marchant, J-B. Le Bouquin, H. Sana

Abstract: Direct dynamical mass measurements of stars with masses above 30 M${}_\odot$ are rare. This is the result of the low yield of the upper initial mass function and the limited number of such systems in eclipsing binaries. Long-period, double-lined spectroscopic binaries that are also resolved astrometrically offer an alternative for obtaining absolute masses of stellar objects. 9 Sgr is one such lon… ▽ More Direct dynamical mass measurements of stars with masses above 30 M${}_\odot$ are rare. This is the result of the low yield of the upper initial mass function and the limited number of such systems in eclipsing binaries. Long-period, double-lined spectroscopic binaries that are also resolved astrometrically offer an alternative for obtaining absolute masses of stellar objects. 9 Sgr is one such long-period, high-mass binary. Unfortunately, a large amount of tension exists between its total dynamical mass inferred from radial velocity measurements and that from astrometric data. We obtained the astrometric orbit from VLTI/PIONIER and VLTI/GRAVITY interferometric measurements. Using archival and new spectroscopy, we performed a grid-based spectral disentangling search to constrain the semi-amplitudes of the radial velocity curves. We computed atmospheric parameters and surface abundances by adjusting \textsc{fastwind} atmosphere models and we compared our results with evolutionary tracks computed with the Bonn Evolutionary Code (BEC). Grid spectral disentangling of 9 Sgr supports the presence of a 53 M${}_\odot$ primary and a 39 M${}_\odot$ secondary. Comparison with BEC evolutionary tracks shows the components of 9 Sgr are most likely coeval with an age of roughly 1 Myr. Our analysis clears up the contradiction between mass and orbital inclination estimates reported in previous studies. We detect the presence of significant CNO-processed material at the surface of the primary, suggesting enhanced internal mixing compared to currently implemented in the BEC models. The present measurements provide a high-quality high-mass anchor to validate stellar evolution models and to test the efficiency of internal mixing processes. △ Less

Submitted 20 May, 2021; originally announced May 2021.

Comments: 14 pages, 5 figures, 4 tables (+6 figures and 3 tables in appendices). Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 651, A119 (2021)

Constraining the Overcontact Phase in Massive Binary Evolution I. Mixing in V382 Cyg, VFTS 352, and OGLE SMC-SC10 108086

Authors: Michael Abdul-Masih, Hugues Sana, Calum Hawcroft, Leonardo A. Almeida, Sarah A. Brands, Selma E. de Mink, Stephen Justham, Norbert Langer, Laurent Mahy, Pablo Marchant, Athira Menon, Joachim Puls, Jon Sundqvist

Abstract: As potential progenitors of several exotic phenomena including gravitational wave sources, magnetic stars, and Be stars, close massive binary systems probe a crucial area of the parameter space in massive star evolution. Despite the importance of these systems, large uncertainties regarding the nature and efficiency of the internal mixing mechanisms still exist. In this work, we aim to provide rob… ▽ More As potential progenitors of several exotic phenomena including gravitational wave sources, magnetic stars, and Be stars, close massive binary systems probe a crucial area of the parameter space in massive star evolution. Despite the importance of these systems, large uncertainties regarding the nature and efficiency of the internal mixing mechanisms still exist. In this work, we aim to provide robust observational constraints on the internal mixing processes by spectroscopically analyzing a sample of three massive overcontact binaries at different metallicities. Using optical phase-resolved spectroscopic data, we perform an atmosphere analysis using more traditional 1D techniques and using state-of-the-art 3D techniques. We compare and contrast the assumptions and results of each technique and investigate how the assumptions affect the final derived atmospheric parameters. We find that in all three cases, both components of system are highly overluminous indicating either efficient internal mixing of helium or previous non-conservative mass transfer. However, we do not find strong evidence of helium or CNO surface abundance changes usually associated with mixing. Additionally, we find that in unequal mass systems, the measured effective temperature and luminosity of the less massive component places it very close to the more massive component on the Hertzsprung-Russell diagram. These results were obtained independently using both of the techniques mentioned above, which suggests that these measurements are robust. The observed discrepancies between the temperature and the surface abundance measurements when compared to theoretical expectations indicate that unaccounted for additional physical mechanisms may be at play. △ Less

Submitted 16 April, 2021; v1 submitted 15 April, 2021; originally announced April 2021.

Comments: Accepted for publication in A&A; 28 pages, 20 figures

The Tarantula Massive Binary Monitoring V. R 144: a wind-eclipsing binary with a total mass > 140 Msun

Authors: T. Shenar, H. Sana, P. Marchant, B. Pablo, N. Richardson, A. F. J. Moffat, T. Van Reeth, R. H. Barba, D. M. Bowman, P. Broos, P. A. Crowther, J. S. Clark, A. de Koter, S. E. de Mink, K. Dsilva, G. Graefener, I. D. Howarth, N. Langer, L. Mahy, J. Maiz Apellaniz, A. M. T. Pollock, F. R. N. Schneider, L. Townsley, J. S. Vink

Abstract: R 144 is the visually brightest WR star in the Large Magellanic Cloud (LMC). R 144 was reported to be a binary, making it potentially the most massive binary thus observed. We perform a comprehensive spectral, photometric, orbital, and polarimetric analysis of R 144. R 144 is an eccentric (e=0.51) 74.2-d binary comprising two relatively evolved (age~2 Myr), H-rich WR stars. The hotter primary (W… ▽ More R 144 is the visually brightest WR star in the Large Magellanic Cloud (LMC). R 144 was reported to be a binary, making it potentially the most massive binary thus observed. We perform a comprehensive spectral, photometric, orbital, and polarimetric analysis of R 144. R 144 is an eccentric (e=0.51) 74.2-d binary comprising two relatively evolved (age~2 Myr), H-rich WR stars. The hotter primary (WN5/6h, T=50 kK) and the cooler secondary (WN6/7h,T=45kK) have nearly equal masses. The combination of low rotation and H-depletion observed in the system is well reproduced by contemporary evolution models that include boosted mass-loss at the upper-mass end. The systemic velocity of R 144 and its relative isolation suggest that it was ejected as a runaway from the neighbouring R 136 cluster. The optical light-curve shows a clear orbital modulation that can be well explained as a combination of two processes: excess emission stemming from wind-wind collisions and double wind eclipses. Our light-curve model implies an orbital inclination of i=60.4+-1.5deg, resulting in accurately constrained dynamical masses of 74+-4 and 69+-4 Msun. Assuming that both binary components are core H-burning, these masses are difficult to reconcile with the derived luminosities (logL1,2 = 6.44, 6.39 [Lsun]), which correspond to evolutionary masses of the order of 110 and 100Msun, respectively. Taken at face value, our results imply that both stars have high classical Eddington factors of Gamma_e = 0.78+-0.1. If the stars are on the main sequence, their derived radii (~25Rsun) suggest that they are only slightly inflated, even at this high Eddington factor. Alternatively, the stars could be core-He burning, strongly inflated from the regular size of classical Wolf-Rayet stars (~1Rsun), a scenario that could help resolve the observed mass discrepancy. △ Less

Submitted 28 October, 2021; v1 submitted 7 April, 2021; originally announced April 2021.

Comments: Accepted to A&A on the April 06 2021. 17 pages + 5 page appendix. Version 28.10.2021: corrected typo in Eq. (2); numerator and denominator were inversed on RHS (only typo; no impact on results)

Journal ref: A&A 650, A147 (2021)

The role of mass transfer and common envelope evolution in the formation of merging binary black holes

Authors: Pablo Marchant, Kaliroë M. W. Pappas, Monica Gallegos-Garcia, Christopher P. L. Berry, Ronald E. Taam, Vicky Kalogera, Philipp Podsiadlowski

Abstract: As the number of observed merging binary black holes (BHs) grows, accurate models are required to disentangle multiple formation channels. In models with isolated binaries, important uncertainties remain regarding the stability of mass transfer (MT) and common-envelope (CE) evolution. To study some of these uncertainties, we have computed simulations using MESA of a $30M_\odot$, low metallicity (… ▽ More As the number of observed merging binary black holes (BHs) grows, accurate models are required to disentangle multiple formation channels. In models with isolated binaries, important uncertainties remain regarding the stability of mass transfer (MT) and common-envelope (CE) evolution. To study some of these uncertainties, we have computed simulations using MESA of a $30M_\odot$, low metallicity ($Z_\odot/10$) star with a BH companion. We developed a prescription to compute MT rates including possible outflows from outer Lagrangian points, and a method to self-consistently determine the core-envelope boundary in the case of CE evolution. We find that binaries survive a CE only if unstable MT happens after the formation of a deep convective envelope, resulting in a narrow range (0.2 dex) in period for envelope ejection. All cases where interaction is initiated with a radiative envelope have large binding energies ($\sim 10^{50}$ erg), and merge during CE even under the assumption that all the internal and recombination energy of the envelope, as well as the energy from an inspiral, is used for ejection. This is independent of core helium ignition for the donor, a condition under which various rapid-population synthesis calculations assume a successful ejection is possible. Moreover, we find that the critical mass ratio for instability is such that for periods between $\sim 1-1000$ days merging binary BHs can be formed via stable MT. A large fraction of these systems overflow their L$_2$ equipotential, in which case we find stable MT produces merging binary BHs even under extreme assumptions of mass and angular momentum outflows. Our conclusions are limited to the study of one donor star, but suggest that population synthesis calculations overestimate the formation rate of merging binary BHs produced by CE evolution, and that stable MT could dominate the rate from isolated binaries. △ Less

Submitted 16 March, 2021; originally announced March 2021.

Comments: 22 pages (including 4 pages of appendices). Abstract abridged for submission to the arXiv. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 650, A107 (2021)

Dynamically inflated wind models of classical Wolf-Rayet stars

Authors: L. G. Poniatowski, J. O. Sundqvist, N. D. Kee, S. P. Owocki, P. Marchant, L. Decin, A. de Koter, L. Mahy, H. Sana

Abstract: Vigorous mass loss in the classical Wolf-Rayet (WR) phase is important for the late evolution and final fate of massive stars. We develop spherically symmetric time-dependent and steady-state hydrodynamical models of the radiation-driven wind outflows and associated mass loss from classical WR stars. The simulations are based on combining the opacities typically used in static stellar structure an… ▽ More Vigorous mass loss in the classical Wolf-Rayet (WR) phase is important for the late evolution and final fate of massive stars. We develop spherically symmetric time-dependent and steady-state hydrodynamical models of the radiation-driven wind outflows and associated mass loss from classical WR stars. The simulations are based on combining the opacities typically used in static stellar structure and evolution models with a simple parametrised form for the enhanced line-opacity expected within a supersonic outflow. Our simulations reveal high mass-loss rates initiated in deep and hot optically thick layers around T\approx 200kK. The resulting velocity structure is non-monotonic and can be separated into three phases: i) an initial acceleration to supersonic speeds ii) stagnation and even deceleration, and iii) an outer region of rapid re-acceleration. The characteristic structures seen in converged steady-state simulations agree well with the outflow properties of our time-dependent models. By directly comparing our dynamic simulations to corresponding hydrostatic models, we demonstrate explicitly that the need to invoke extra energy transport in convectively inefficient regions of stellar structure and evolution models is merely an artefact of enforcing a hydrostatic outer boundary. Moreover, the "dynamically inflated" inner regions of our simulations provide a natural explanation for the often-found mismatch between predicted hydrostatic WR radii and those inferred from spectroscopy. Finally, we contrast our simulations with alternative recent WR wind models based on co-moving frame radiative transfer for computing the radiation force. Since CMF transfer currently cannot handle non-monotonic velocity fields, the characteristic deceleration regions found here are avoided in such simulations by invoking an ad-hoc very high degree of clumping. △ Less

Submitted 10 December, 2020; originally announced December 2020.

Comments: 15 pages, 9 figures

Journal ref: A&A 647, A151 (2021)

Detailed evolutionary models of massive contact binaries: I. Model grids and synthetic populations for the Magellanic Clouds

Authors: Athira Menon, Norbert Langer, Selma E. de Mink, Stephen Justham, Koushik Sen, Dorottya Szécsi, Alex de Koter, Michael Abdul-Masih, Hugues Sana, Laurent Mahy, Pablo Marchant

Abstract: The majority of close massive binary stars with initial periods of a few days experience a contact phase, in which both stars overflow their Roche lobes simultaneously. We perform the first dedicated study of the evolution of massive contact binaries and provide a comprehensive prediction of their observed properties. We compute 2790 detailed binary models for the Large and Small Magellanic Clouds… ▽ More The majority of close massive binary stars with initial periods of a few days experience a contact phase, in which both stars overflow their Roche lobes simultaneously. We perform the first dedicated study of the evolution of massive contact binaries and provide a comprehensive prediction of their observed properties. We compute 2790 detailed binary models for the Large and Small Magellanic Clouds each, assuming a conservative mass transfer. The initial parameter space for both grids span total masses from 20 to 80$\,\textrm{M}_{\odot}$, orbital periods of 0.6 to 2 days and mass ratios of 0.6 to 1.0. We find that models that remain in contact over nuclear timescales evolve towards equal masses, echoing the mass ratios of their observed counterparts. Ultimately, the fate of our nuclear-timescale models is to merge on the main sequence. Our predicted period-mass ratio distributions of O-type contact binaries are similar for both galaxies, and we expect 10 such systems together in both Magellanic Clouds. While we can largely reproduce the observed distribution, we over-estimate the population of equal-mass contact binaries. This situation is somewhat remedied if we also account for binaries that are approaching contact. Our theoretical distributions work particularly well for contact binaries with periods $<$2 days and total masses $\lessapprox45\,\textrm{M}_{\odot}$. We expect stellar winds, non-conservative mass transfer and envelope inflation to have played a role in the formation of the more massive and longer-period contact binaries. △ Less

Submitted 4 August, 2021; v1 submitted 26 November, 2020; originally announced November 2020.

Comments: 23 pages, 11 figures, accepted for publication in MNRAS

One Channel to Rule Them All? Constraining the Origins of Binary Black Holes using Multiple Formation Pathways

Authors: Michael Zevin, Simone S. Bavera, Christopher P. L. Berry, Vicky Kalogera, Tassos Fragos, Pablo Marchant, Carl L. Rodriguez, Fabio Antonini, Daniel E. Holz, Chris Pankow

Abstract: The second LIGO-Virgo catalog of gravitational wave transients has more than quadrupled the observational sample of binary black holes. We analyze this catalog using a suite of five state-of-the-art binary black hole population models covering a range of isolated and dynamical formation channels and infer branching fractions between channels as well as constraints on uncertain physical processes t… ▽ More The second LIGO-Virgo catalog of gravitational wave transients has more than quadrupled the observational sample of binary black holes. We analyze this catalog using a suite of five state-of-the-art binary black hole population models covering a range of isolated and dynamical formation channels and infer branching fractions between channels as well as constraints on uncertain physical processes that impact the observational properties of mergers. Given our set of formation models, we find significant differences between the branching fractions of the underlying and detectable populations, and that the diversity of detections suggests that multiple formation channels are at play. A mixture of channels is strongly preferred over any single channel dominating the detected population: an individual channel does not contribute to more than $\simeq 70\%$ of the observational sample of binary black holes. We calculate the preference between the natal spin assumptions and common envelope efficiencies in our models, favoring natal spins of isolated black holes of $\lesssim 0.1$, and marginally preferring common envelope efficiencies of $\gtrsim 2.0$ while strongly disfavoring highly inefficient common envelopes. We show that it is essential to consider multiple channels when interpreting gravitational wave catalogs, as inference on branching fractions and physical prescriptions becomes biased when contributing formation scenarios are not considered or incorrect physical prescriptions are assumed. Although our quantitative results can be affected by uncertain assumptions in model predictions, our methodology is capable of including models with updated theoretical considerations and additional formation channels. △ Less

Submitted 6 April, 2021; v1 submitted 19 November, 2020; originally announced November 2020.

Comments: 27 pages (14 pages main text + 13 pages appendices/references), 8 figures, 1 table, published in ApJ

Journal ref: The Astrophysical Journal 910, 152 (2021)

BAT99 126: A multiple Wolf-Rayet system in the Large Magellanic Cloud with a massive near-contact binary

Authors: Soetkin Janssens, Tomer Shenar, Laurent Mahy, Pablo Marchant, Hugues Sana, Julia Bodensteiner

Abstract: BAT99 126 is a multiple system in the Large Magellanic Cloud containing a Wolf-Rayet (WR) star, which has a reported spectroscopic (orbital) period of 25.5 days and a photometric (orbital) period of 1.55 days, and hence is potentially one of the shortest WR binaries known to date. Such short-period binary systems containing a WR star in low-metallicity environments are prime candidate progenitors… ▽ More BAT99 126 is a multiple system in the Large Magellanic Cloud containing a Wolf-Rayet (WR) star, which has a reported spectroscopic (orbital) period of 25.5 days and a photometric (orbital) period of 1.55 days, and hence is potentially one of the shortest WR binaries known to date. Such short-period binary systems containing a WR star in low-metallicity environments are prime candidate progenitors of black-hole (BH) mergers. By thoroughly analysing the spectroscopic and photometric data, we aim to establish the true multiplicity of BAT99 126, characterise the orbit(s) of the system, measure the physical properties of its individual components, and determine the overall evolutionary status of the system. Using newly acquired high resolution spectra, we measured radial velocities and performed a spectral analysis of the individual components. We estimated the age of the system and derived an evolutionary scenario for the 1.55-day system. BAT99 126 comprises at least four components. The 1.55-day photometric signal originates in an eclipsing O+O binary of spectral types O4 V and O6.5 V, and masses of 36+/-5 MSun and 15+/-2 MSun, respectively. For the WR star, we derived a spectral type WN2.5-3 and reject the previously reported 25.5-d period. We find evidence in the composite spectra of a fourth component, a B1 V star, which shows radial velocity variation. The configurations of the B-type star, WR star, and possible additional undetected components remain unknown. We estimated the age of the system to be 4.2 Myr. The O+O binary likely went through a phase of conservative mass transfer and is currently a near-contact system. We show that BAT99 126 is a multiple - quadruple or higher-order - system with a total initial mass > 160 MSun. The 1.55-day O+O binary most likely will not evolve towards a BH+BH merger, but instead will merge before the components collapse to BHs. △ Less

Submitted 18 November, 2020; originally announced November 2020.

Comments: 14+2 pages, 12+1 figures, accepted in A&A

Journal ref: A&A 646, A33 (2021)

The impact of mass-transfer physics on the observable properties of field binary black hole populations

Authors: Simone S. Bavera, Tassos Fragos, Michael Zevin, Christopher P. L. Berry, Pablo Marchant, Jeff J. Andrews, Scott Coughlin, Aaron Dotter, Konstantinos Kovlakas, Devina Misra, Juan G. Serra-Perez, Ying Qin, Kyle A. Rocha, Jaime Román-Garza, Nam H. Tran, Emmanouil Zapartas

Abstract: We study the impact of mass-transfer physics on the observable properties of binary black hole populations formed through isolated binary evolution. We investigate the impact of mass-accretion efficiency onto compact objects and common-envelope efficiency on the observed distributions of $χ_{eff}$, $M_{chirp}$ and $q$. We find that low common envelope efficiency translates to tighter orbits post c… ▽ More We study the impact of mass-transfer physics on the observable properties of binary black hole populations formed through isolated binary evolution. We investigate the impact of mass-accretion efficiency onto compact objects and common-envelope efficiency on the observed distributions of $χ_{eff}$, $M_{chirp}$ and $q$. We find that low common envelope efficiency translates to tighter orbits post common envelope and therefore more tidally spun up second-born black holes. However, these systems have short merger timescales and are only marginally detectable by current gravitational-waves detectors as they form and merge at high redshifts ($z\sim 2$), outside current detector horizons. Assuming Eddington-limited accretion efficiency and that the first-born black hole is formed with a negligible spin, we find that all non-zero $χ_{eff}$ systems in the detectable population can come only from the common envelope channel as the stable mass-transfer channel cannot shrink the orbits enough for efficient tidal spin-up to take place. We find the local rate density ($z\simeq 0.01$) for the common envelope channel is in the range $\sim 17-113~Gpc^{-3}yr^{-1}$ considering a range of $α_{CE} \in [0.2,5.0]$ while for the stable mass transfer channel the rate density is $\sim 25~Gpc^{-3}yr^{-1}$. The latter drops by two orders of magnitude if the mass accretion onto the black hole is not Eddington limited because conservative mass transfer does not shrink the orbit as efficiently as non-conservative mass transfer does. Finally, using GWTC-2 events, we constrain the lower bound of branching fraction from other formation channels in the detected population to be $\sim 0.2$. Assuming all remaining events to be formed through either stable mass transfer or common envelope channels, we find moderate to strong evidence in favour of models with inefficient common envelopes. △ Less

Submitted 15 February, 2021; v1 submitted 30 October, 2020; originally announced October 2020.

Comments: 26 pages, 13 figures, accepted for publication in A&A

Journal ref: A&A 647, A153 (2021)

Chemically Homogeneous Evolution: A rapid population synthesis approach

Authors: Jeff Riley, Ilya Mandel, Pablo Marchant, Ellen Butler, Kaila Nathaniel, Coenraad Neijssel, Spencer Shortt, Alejandro Vigna-Gomez

Abstract: We explore chemically homogeneous evolution (CHE) as a formation channel for massive merging binary black holes (BBHs). We develop methods to include CHE in a rapid binary population synthesis code, Compact Object Mergers: Population Astrophysics and Statistics (COMPAS), which combines realistic models of binary evolution with cosmological models of the star-formation history of the Universe. For… ▽ More We explore chemically homogeneous evolution (CHE) as a formation channel for massive merging binary black holes (BBHs). We develop methods to include CHE in a rapid binary population synthesis code, Compact Object Mergers: Population Astrophysics and Statistics (COMPAS), which combines realistic models of binary evolution with cosmological models of the star-formation history of the Universe. For the first time, we simultaneously explore conventional isolated binary star evolution under the same set of assumptions. This approach allows us to constrain population properties and make simultaneous predictions about the gravitational-wave detection rates of BBH mergers for the CHE and conventional formation channels. The overall mass distribution of detectable BBHs is consistent with existing gravitational-wave observations. We find that the CHE channel may yield up to ~70% of all gravitational-wave detections of BBH mergers coming from isolated binary evolution. △ Less

Submitted 2 May, 2021; v1 submitted 29 September, 2020; originally announced October 2020.

Comments: 15 pages, 16 figures

Progenitors of Type IIb Supernovae: II. Observable Properties

Authors: Niharika Sravan, Pablo Marchant, Vassiliki Kalogera, Dan Milisavljevic, Raffaella Margutti

Abstract: Type IIb supernovae (SNe IIb) present a unique opportunity for investigating the evolutionary channels and mechanisms governing the evolution of stripped-envelope SN progenitors due to a variety of observational constraints available. Comparison of these constraints with the full distribution of theoretical properties not only help ascertain the prevalence of observed properties in nature, but can… ▽ More Type IIb supernovae (SNe IIb) present a unique opportunity for investigating the evolutionary channels and mechanisms governing the evolution of stripped-envelope SN progenitors due to a variety of observational constraints available. Comparison of these constraints with the full distribution of theoretical properties not only help ascertain the prevalence of observed properties in nature, but can also reveal currently unobserved populations. In this follow-up paper, we use the large grid of models presented in Sravan et al. 2019 to derive distributions of single and binary SNe IIb progenitor properties and compare them to constraints from three independent observational probes: multi-band SN light-curves, direct progenitor detections, and X-ray/radio observations. Consistent with previous work, we find that while current observations exclude single stars as SN IIb progenitors, SN IIb progenitors in binaries can account for them. We also find that the distributions indicate the existence of an unobserved dominant population of binary SNe IIb at low metallicity that arise due to mass transfer initiated on the Hertzsprung Gap. In particular, our models indicate the existence of a group of highly stripped (envelope mass ~0.1-0.2 M_sun) progenitors that are compact (<50 R_sun) and blue (T_eff <~ 10^5K) with ~10^4.5-10^5.5 L_sun and low density circumstellar mediums. As discussed in Sravan et al. 2019, this group is necessary to account for SN IIb fractions and likely exist regardless of metallicity. The detection of the unobserved populations indicated by our models would support weak stellar winds and inefficient mass transfer in SN IIb progenitors. △ Less

Submitted 14 September, 2020; originally announced September 2020.

Comments: Resubmitted to the Astrophysical Journal after incorporating suggestions from the referee

Luminous supernovae associated with ultra-long gamma-ray bursts from hydrogen-free progenitors extended by pulsational pair-instability

Authors: Takashi J. Moriya, Pablo Marchant, Sergei I. Blinnikov

Abstract: We show that the luminous supernovae (SNe) associated with ultra-long gamma-ray bursts (GRBs) can be related to the slow cooling from the explosions of hydrogen-free progenitors extended by pulsational pair-instability. In the accompanying paper (Marchant & Moriya 2020), we have shown that some rapidly-rotating hydrogen-free GRB progenitors that experience pulsational pair-instability can keep an… ▽ More We show that the luminous supernovae (SNe) associated with ultra-long gamma-ray bursts (GRBs) can be related to the slow cooling from the explosions of hydrogen-free progenitors extended by pulsational pair-instability. In the accompanying paper (Marchant & Moriya 2020), we have shown that some rapidly-rotating hydrogen-free GRB progenitors that experience pulsational pair-instability can keep an extended structure caused by pulsational pair-instability until the core collapse. Such progenitors have large radii exceeding 10 Rsun and they sometimes reach beyond 1000 Rsun at the time of the core collapse. They are, therefore, promising progenitors of ultra-long GRBs. We here perform the light-curve modeling of the explosions of one extended hydrogen-free progenitor with a radius of 1962 Rsun. Thanks to the large progenitor radius, the ejecta experience slow cooling after the shock breakout and they become rapidly evolving (<~ 10 days) luminous (>~ 1e43 erg/s) SNe in optical even without the energy input from the 56Ni nuclear decay when the explosion energy is more than 1e52 erg. The 56Ni decay energy input can affect the light curves after the optical light-curve peak and make the light-curve decay slow when the 56Ni mass is around 1 Msun. They also have fast photospheric velocity above 10,000 km/s and hot photospheric temperature above 10,000 K at around the peak luminosity. We find that the rapid rise and luminous peak found in the optical light curve of SN 2011kl, which is associated with the ultra-long GRB 111209A, can be explained as the cooling phase of the extended progenitor. The ultra-long GRB progenitors proposed in Marchant & Moriya (2020) can explain both the ultra-long GRB duration and the accompanying SN properties. When the GRB jet is off-axis or choked, the luminous SNe could be observed as fast blue optical transients without accompanying GRBs. (abridged) △ Less

Submitted 4 September, 2020; v1 submitted 13 July, 2020; originally announced July 2020.

Comments: 5 pages, 5 figures, accepted by Astronomy & Astrophysics Letters

Journal ref: Astronomy & Astrophysics, Volume 641, id.L10, 5 pp. (2020)

The impact of stellar rotation on the black hole mass-gap from pair-instability supernovae

Authors: Pablo Marchant, Takashi Moriya

Abstract: Models of pair-instability supernovae (PISNe) predict a gap in black hole (BH) masses between $\sim 45M_\odot-120M_\odot$, which is referred to as the upper BH mass-gap. With the advent of gravitational-wave astrophysics it has become possible to test this prediction, and there is an important associated effort to understand what theoretical uncertainties modify the boundaries of this gap. In this… ▽ More Models of pair-instability supernovae (PISNe) predict a gap in black hole (BH) masses between $\sim 45M_\odot-120M_\odot$, which is referred to as the upper BH mass-gap. With the advent of gravitational-wave astrophysics it has become possible to test this prediction, and there is an important associated effort to understand what theoretical uncertainties modify the boundaries of this gap. In this work we study the impact of rotation on the hydrodynamics of PISNe, which leave no compact remnant, as well as the evolution of pulsational-PISNe (PPISNe), which undergo thermonuclear eruptions before forming a compact object. We perform simulations of non-rotating and rapidly-rotating stripped helium stars in a metal poor environment $(Z_\odot/50)$ in order to resolve the lower edge of the upper mass-gap. We find that the outcome of our simulations is dependent on the efficiency of angular momentum transport, with models that include efficient coupling through the Spruit-Tayler dynamo shifting the lower edge of the mass-gap upwards by $\sim 4\%$, while simulations that do not include this effect shift it upwards by $\sim 15\%$. From this, we expect the lower edge of the upper mass-gap to be dependent on BH spin, which can be tested as the number of observed BH mergers increases. Moreover, we show that stars undergoing PPISNe have extended envelopes ($R\sim 10-1000~R_\odot$) at iron-core collapse, making them promising progenitors for ultra-long gamma-ray bursts. △ Less

Submitted 13 July, 2020; originally announced July 2020.

Comments: 9 pages, including 4 pages of appendices

Journal ref: A&A 640, L18 (2020)

A spectroscopic multiplicity survey of Galactic Wolf-Rayet stars. I. The northern WC sequence

Authors: Karan Dsilva, Tomer Shenar, Hugues Sana, Pablo Marchant

Abstract: 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… ▽ More 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 ${\sim}0.4$ in the Galaxy but no correction for observational biases has been attempted. The aim of this study is to conduct a homogeneous radial-velocity survey of a magnitude-limited ($V$ $\leq 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.2m Mercator telescope on the island of La Palma. 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. 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~$<~\sim$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. △ Less

Submitted 24 June, 2020; originally announced June 2020.

Journal ref: A&A 641, A26 (2020)

Is HR 6819 a triple system containing a black hole? -- An alternative explanation

Authors: J. Bodensteiner, T. Shenar, L. Mahy, M. Fabry, P. Marchant, M. Abdul-Masih, G. Banyard, D. M. Bowman, K. Dsilva, A. J. Frost, C. Hawcroft, M. Reggiani, H. Sana

Abstract: HR 6819 was recently proposed to be a triple system consisting of an inner B-type giant + black hole binary with an orbital period of 40d and an outer Be tertiary. This interpretation is mainly based on two inferences: that the emission attributed to the outer Be star is stationary, and that the inner star, which is used as mass calibrator for the black hole, is a B-type giant. We re-investigate t… ▽ More HR 6819 was recently proposed to be a triple system consisting of an inner B-type giant + black hole binary with an orbital period of 40d and an outer Be tertiary. This interpretation is mainly based on two inferences: that the emission attributed to the outer Be star is stationary, and that the inner star, which is used as mass calibrator for the black hole, is a B-type giant. We re-investigate the properties of HR 6819 by spectral disentangling and an atmosphere analysis of the disentangled spectra to search for a possibly simpler alternative explanation for HR 6819. Disentangling implies that the Be component is not a static tertiary, but rather a component of the binary in the 40-d orbit. The inferred radial velocity amplitudes imply an extreme mass ratio of M_2/M_1 = 15 +/- 3. We infer spectroscopic masses of 0.4$^{+0.3}_{-0.1}$ Msun and 6$^{+5}_{-3}$ Msun for the primary and secondary, which agree well with the dynamical masses for an inclination of i = 32 deg. This indicates that the primary might be a stripped star rather than a B-type giant. Evolutionary modelling suggests that a possible progenitor system would be a tight (P_i ~ 2d) B+B binary system that experienced conservative mass transfer. While the observed nitrogen enrichment of the primary conforms with the predictions of the evolutionary models, we find no indications for the predicted He enrichment. We suggest that HR 6819 is a binary system consisting of a stripped B-type primary and a rapidly-rotating Be star that formed from a previous mass-transfer event. In the framework of this interpretation, HR 6819 does not contain a black hole. Interferometry can distinguish between these two scenarios by providing an independent measurement of the separation between the visible components. △ Less

Submitted 18 June, 2020; originally announced June 2020.

Comments: Submitted to A&A, 13 pages (16 figures and 2 tables); 4 pages supplementary material (4 figures and 4 tables). Comments are welcome

Journal ref: A&A 641, A43 (2020)

The 'hidden' companion in LB-1 unveiled by spectral disentangling

Authors: T. Shenar, J. Bodensteiner, M. Abdul-Masih, M. Fabry, L. Mahy, P. Marchant, G. Banyard, D. M. Bowman, K. Dsilva, C. Hawcroft, M. Reggiani, H. Sana

Abstract: 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 Msun black hole. Recent analyses have implied that the visible primary star is a stripped He-rich star. However, the nature of the secondary, which was proposed to be a black hole, a neutron star, or a main sequenc… ▽ More 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 Msun black hole. Recent analyses have implied that the visible primary star is a stripped He-rich star. However, the nature of the secondary, which was proposed to be a black hole, a neutron star, or a main sequence star, remains unknown. Based on 26 newly acquired spectroscopic observations, we perform an orbital analysis and spectral disentangling of LB-1 to elucidate the nature of the system. Our analysis reveals 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 Msun) for a B-type star and it is He- and N-rich. The "hidden" secondary, which contributes about 45% of the optical flux, is a rapidly rotating (vsini ~ 300 km/s) B3 V star with a decretion disk -- a Be star. Hence, 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. △ Less

Submitted 11 July, 2020; v1 submitted 27 April, 2020; originally announced April 2020.

Comments: Accepted to A&A. 4 pages (1 figure, 1 table); 3 pages of supporting material (7 figures and 1 table)

Journal ref: A&A 639, L6 (2020)

Cosmic Rates of Black Hole Mergers and Pair-Instability Supernovae from Chemically Homogeneous Binary Evolution

Authors: Lise du Buisson, Pablo Marchant, Philipp Podsiadlowski, Chiaki Kobayashi, Filipe B. Abdalla, Philip Taylor, Ilya Mandel, Selma E. de Mink, Takashi J. Moriya, Norbert Langer

Abstract: During the first three observing runs of the Advanced gravitational-wave detector network, the LIGO/Virgo collaboration detected several black hole binary (BHBH) mergers. As the population of detected BHBH mergers grows, it will become possible to constrain different channels for their formation. Here we consider the chemically homogeneous evolution (CHE) channel in close binaries, by performing p… ▽ More During the first three observing runs of the Advanced gravitational-wave detector network, the LIGO/Virgo collaboration detected several black hole binary (BHBH) mergers. As the population of detected BHBH mergers grows, it will become possible to constrain different channels for their formation. Here we consider the chemically homogeneous evolution (CHE) channel in close binaries, by performing population synthesis simulations that combine realistic binary models with detailed cosmological calculations of the chemical and star-formation history of the Universe. This allows us to constrain population properties, as well as cosmological and aLIGO/aVirgo detection rates of BHBH mergers formed through this pathway. We predict a BHBH merger rate at redshift zero of $5.8 \, \textrm{Gpc}^{-3} \textrm{yr}^{-1}$ through the CHE channel, to be compared with aLIGO/aVirgo's measured rate of ${53.2}_{-28.2}^{+55.8} \, \text{Gpc}^{-3} \text{yr}^{-1}$, and find that eventual merger systems have BH masses in the range $17 - 43 \, \textrm{M}_{\odot}$ below the pair-instability supernova (PISN) gap, and $>124 \, \textrm{M}_{\odot}$ above the PISN gap. We investigate effects of momentum kicks during black hole formation, and calculate cosmological and magnitude limited PISN rates. We also study the effects of high-redshift deviations in the star formation rate. We find that momentum kicks tend to increase delay times of BHBH systems, and our magnitude limited PISN rate estimates indicate that current deep surveys should be able to detect such events. Lastly, we find that our cosmological merger rate estimates change by at most $\sim 8\%$ for mild deviations of the star formation rate in the early Universe, and by up to $\sim 40\%$ for extreme deviations. △ Less

Submitted 8 November, 2020; v1 submitted 26 February, 2020; originally announced February 2020.

Comments: 20 pages, 18 figures, 4 tables. References added, text clarifications added, results and figures unchanged

Sensitivity of the lower-edge of the pair instability black hole mass gap to the treatment of time dependent convection

Authors: M. Renzo, R. J. Farmer, S. Justham, S. E. de Mink, Y. Götberg, P. Marchant

Abstract: 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 timescales shorter than the convective turn-over timescale. This invalidates the steady-state assumption on w… ▽ More 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 timescales shorter than the convective turn-over timescale. 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 \texttt{MESA}: (i) using the default implementation of \cite{paxton:18} and (ii) solving an additional equation accounting for the timescale 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 $\sim$\,$5\,M_\odot$. The differences are much smaller for the progenitors which determine the maximum mass of BHs below the gap. This prediction is robust at $M_{\rm BH, max}\simeq 48\,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. △ Less

Submitted 19 February, 2020; originally announced February 2020.

Comments: 7 pages + 1 appendix, accepted in MNRAS, online results at https://zenodo.org/record/3406320