-
Formation of Wind-Fed Black Hole High-mass X-ray Binaries: The Role of Roche-lobe-Overflow Post Black-Hole Formation
Authors:
Zepei Xing,
Tassos Fragos,
Emmanouil Zapartas,
Tom M. Kwan,
Lixin Dai,
Ilya Mandel,
Matthias U. Kruckow,
Max Briel,
Jeff J. Andrews,
Simone S. Bavera,
Seth Gossage,
Konstantinos Kovlakas,
Kyle A. Rocha,
Meng Sun,
Philipp M. Srivastava
Abstract:
The three dynamically confirmed wind-fed black hole high-mass X-ray binaries (BH-HMXBs) are suggested to all contain a highly spinning black hole (BH). However, based on the theories of efficient angular momentum transport inside the stars, we expect that the first-born BHs in binary systems should have low spins, which is consistent with gravitational-wave observations. As a result, the origin of…
▽ More
The three dynamically confirmed wind-fed black hole high-mass X-ray binaries (BH-HMXBs) are suggested to all contain a highly spinning black hole (BH). However, based on the theories of efficient angular momentum transport inside the stars, we expect that the first-born BHs in binary systems should have low spins, which is consistent with gravitational-wave observations. As a result, the origin of the high BH spins measured in wind-fed BH-HMXBs remains a mystery. In this paper, we conduct a binary population synthesis study on wind-fed BH-HMXBs at solar metallicity with the use of the newly developed code POSYDON, considering three scenarios for BH accretion: Eddington-limited, moderately super-Eddington, and fully conservative accretion. Taking into account the conditions for accretion-disk formation, we find that regardless of the accretion model, these systems are more likely to have already experienced a phase of Roche-lobe overflow after the BH formation. To account for the extreme BH spins, highly conservative accretion onto BHs is required, when assuming the accreted material carries the specific angular momentum at the innermost stable orbit. Besides, in our simulations we found that the systems with donor stars within the mass range of $10-20\,M_{\odot}$ are prevalent, posing a challenge in explaining simultaneously all observed properties of the BH-HMXB in our Galaxy, Cygnus X-1, and potentially hinting that the accretion efficiency onto non-degenerate stars, before the formation of the BH, is also more conservative than assumed in our simulations.
△ Less
Submitted 28 June, 2024;
originally announced July 2024.
-
To Be or not to Be: the role of rotation in modeling Galactic Be X-ray Binaries
Authors:
Kyle Akira Rocha,
Vicky Kalogera,
Zoheyr Doctor,
Jeff J. Andrews,
Meng Sun,
Seth Gossage,
Simone S. Bavera,
Tassos Fragos,
Konstantinos Kovlakas,
Matthias U. Kruckow,
Devina Misra,
Philipp M. Srivastava,
Zepei Xing,
Emmanouil Zapartas
Abstract:
Be X-ray binaries (Be-XRBs) are crucial in understanding high-mass X-ray binaries, featuring a rapidly rotating Be star and a neutron star companion in an eccentric orbit, intermittently accreting material from the Be star's decretion disk. Originating from binary stellar evolution, Be-XRBs are of significant interest to binary population synthesis (BPS) studies, encapsulating the physics of super…
▽ More
Be X-ray binaries (Be-XRBs) are crucial in understanding high-mass X-ray binaries, featuring a rapidly rotating Be star and a neutron star companion in an eccentric orbit, intermittently accreting material from the Be star's decretion disk. Originating from binary stellar evolution, Be-XRBs are of significant interest to binary population synthesis (BPS) studies, encapsulating the physics of supernovae, common envelope, and mass transfer (MT). Using the POSYDON BPS code, employing pre-computed grids of detailed binary stellar evolution models, we investigate the Galactic Be-XRB population. POSYDON incorporates stellar rotation self-consistently during MT phases, enabling a detailed examination of the rotational distribution of Be stars. Our fiducial BPS and Be-XRB model align well with the orbital properties of Galactic Be-XRBs, emphasizing the role of rotational constraints. Our modeling reveals a bimodal rotational distribution of Be-XRB-like systems, in excellent agreement with literature values. All Be-XRBs undergo an MT phase before the first compact object forms, with over half experiencing a second MT phase from a stripped helium companion (Case BB). Computing rotationally-limited MT efficiencies and applying them to our population, we find that the majority of Be-XRBs have undergone highly non-conservative MT (beta ~ 0.15). Our study underscores the importance of detailed angular momentum modeling during MT in interpreting Be-XRB populations, emphasizing this population as a key probe for the stability and efficiency of MT in interacting binaries.
△ Less
Submitted 11 March, 2024;
originally announced March 2024.
-
Exploring the nature of ultra-luminous X-ray sources across stellar population ages using detailed binary evolution calculations
Authors:
Devina Misra,
Konstantinos Kovlakas,
Tassos Fragos,
Jeff J. Andrews,
Simone S. Bavera,
Emmanouil Zapartas,
Zepei Xing,
Aaron Dotter,
Kyle Akira Rocha,
Philipp M. Srivastava,
Meng Sun
Abstract:
Ultra-luminous X-ray sources (ULXs) are sources observed to exceed the Eddington limit of a stellar-mass black hole (BH). A fraction of ULX sources show X-ray pulses which are evidence for accreting neutron stars (NSs). Theoretical studies have suggested that NSs dominate the compact objects of intrinsic ULXs, even though the majority of observed sample is non-pulsating, implying that X-ray pulses…
▽ More
Ultra-luminous X-ray sources (ULXs) are sources observed to exceed the Eddington limit of a stellar-mass black hole (BH). A fraction of ULX sources show X-ray pulses which are evidence for accreting neutron stars (NSs). Theoretical studies have suggested that NSs dominate the compact objects of intrinsic ULXs, even though the majority of observed sample is non-pulsating, implying that X-ray pulses from many NS ULXs are unobservable. We use POSYDON to generate and study X-ray binary populations spanning starburst ages 5 to 1000Myr. Following theoretical predictions for the alignment of the NS spin axis with the accretion disc, we estimate the required accreted mass in ULXs so that the alignment suppresses observable X-ray pulses. While the properties of ULXs are sensitive to model assumptions, there are certain trends that the populations follow. Young and old stellar populations are dominated by BH and NS accretors, respectively. The donors go from massive H-rich main-sequence (MS) stars in young populations (<100Myr) to low-mass post-MS H-rich stars in older populations (>100Myr), with stripped He-rich giant stars dominating the populations at around 100Myr. In addition, we find that NS ULXs exhibit stronger geometrical beaming than BH ULXs, leading to an under-representation of NS accretors in observed populations. Coupled with our finding that X-ray pulses are suppressed in at least 60% of the NS ULXs, we suggest that the observed fraction of ULXs with detectable X-ray pulses is very small, in agreement with observations. This study investigates the effects of age on ULXs as well as the effects of different model assumptions on ULX demographics. We show that geometrical beaming and the mass-accretion phase are critical aspects of understanding ULX observations. Our results suggest that even though most ULXs have accreting NSs, those with observable X-ray pulses would be very few.
△ Less
Submitted 20 December, 2023; v1 submitted 27 September, 2023;
originally announced September 2023.
-
From ZAMS to Merger: Detailed Binary Evolution Models of Coalescing Neutron Star-Black Hole Systems at Solar Metallicity
Authors:
Zepei Xing,
Simone S. Bavera,
Tassos Fragos,
Matthias U. Kruckow,
Jaim Román-Garza,
Jeff J. Andrews,
Aaron Dotter,
Konstantinos Kovlakas,
Devina Misra,
Philipp M. Srivastava,
Kyle A. Rocha,
Meng Sun,
Emmanouil Zapartas
Abstract:
Neutron star $-$ black hole (NSBH) merger events bring us new opportunities to constrain theories of stellar and binary evolution, and understand the nature of compact objects. In this work, we investigate the formation of merging NSBH binaries at solar metallicity by performing a binary population synthesis study of merging NSBH binaries with the newly developed code POSYDON. The latter incorpora…
▽ More
Neutron star $-$ black hole (NSBH) merger events bring us new opportunities to constrain theories of stellar and binary evolution, and understand the nature of compact objects. In this work, we investigate the formation of merging NSBH binaries at solar metallicity by performing a binary population synthesis study of merging NSBH binaries with the newly developed code POSYDON. The latter incorporates extensive grids of detailed single and binary evolution models, covering the entire evolution of a double compact object progenitor. We explore the evolution of NSBHs originating from different formation channels, which in some cases differ from earlier studies performed with rapid binary population synthesis codes. Then, we present the population properties of merging NSBH systems and their progenitors such as component masses, orbital features, and BH spins, and investigate the model uncertainties in our treatment of common envelope (CE) evolution and core-collapse process. We find that at solar metallicity, under the default model assumptions, most of the merging NSBHs have BH masses in a range of $3-11\,M{_\odot}$ and chirp masses within $1.5-4\,M{_\odot}$. Independently of our model variations, the BH always forms first with dimensionless spin parameter $\lesssim 0.2$, which is correlated to the initial binary orbital period. Some BHs can subsequently spin up moderately ($χ_{\rm BH} \lesssim 0.4$) due to mass transfer, which we assume to be Eddington limited. Binaries that experienced CE evolution rarely demonstrate large tilt angles. Conversely, approximately $40\%$ of the binaries that undergo only stable mass transfer without CE evolution contain an anti-aligned BH. Finally, accounting for uncertainties in both the population modeling and the NS equation of state, we find that $0-18.6\%$ of NSBH mergers may be accompanied by an electromagnetic counterpart.
△ Less
Submitted 18 September, 2023;
originally announced September 2023.
-
The formation of $30\,M_\odot$ merging black holes at solar metallicity
Authors:
Simone S. Bavera,
Tassos Fragos,
Emmanouil Zapartas,
Jeff J. Andrews,
Vicky Kalogera,
Christopher P. L. Berry,
Matthias Kruckow,
Aaron Dotter,
Konstantinos Kovlakas,
Devina Misra,
Kyle A. Rocha,
Philipp M. Srivastava,
Meng Sun,
Zepei Xing
Abstract:
The maximum mass of black holes formed in isolated binaries is determined by stellar winds and the interactions between the binary components. We consider for the first time fully self-consistent detailed stellar structure and binary evolution calculations in population-synthesis models and a new, qualitatively different picture emerges for the formation of black-hole binaries, compared to studies…
▽ More
The maximum mass of black holes formed in isolated binaries is determined by stellar winds and the interactions between the binary components. We consider for the first time fully self-consistent detailed stellar structure and binary evolution calculations in population-synthesis models and a new, qualitatively different picture emerges for the formation of black-hole binaries, compared to studies employing rapid population synthesis models. We find merging binary black holes can form with a non-negligible rate ($\sim 4\times10^{-7}\,M_\odot^{-1}$) at solar metallicity. Their progenitor stars with initial masses $\gtrsim 50\,M_\odot$ do not expand to supergiant radii, mostly avoiding significant dust-driven or luminous blue variable winds. Overall, the progenitor stars lose less mass in stellar winds, resulting in black holes as massive as $\sim 30\,M_\odot$, and, approximately half of them avoid a mass-transfer episode before forming the first-born black hole. Finally, binaries with initial periods of a few days, some of which may undergo episodes of Roche-lobe overflow mass transfer, result in mildly spinning first-born black holes, $χ_\mathrm{BH1} \lesssim 0.2$, assuming efficient angular-momentum transport.
△ Less
Submitted 21 December, 2022;
originally announced December 2022.
-
Searching for candidates of coalescing binary black holes formed through chemically homogeneous evolution in GWTC-3
Authors:
Ying Qin,
Yuan-Zhu Wang,
Simone S. Bavera,
Shichao Wu,
Georges Meynet,
Yi-Ying Wang,
Rui-Chong Hu,
Jin-Ping Zhu,
Dong-Hong Wu,
Xin-Wen Shu,
Fang-Kun Peng,
Han-Feng Song,
Da-Ming Wei
Abstract:
The LIGO, Virgo, and KAGRA (LVK) collaboration has announced 90 coalescing binary black holes (BBHs) with $p_{\rm astro} > 50\%$ to date, however, the origin of their formation channels is still an open scientific question. Given various properties of BBHs (BH component masses and individual spins) inferred using the default priors by the LVK, independent groups have been trying to explain the for…
▽ More
The LIGO, Virgo, and KAGRA (LVK) collaboration has announced 90 coalescing binary black holes (BBHs) with $p_{\rm astro} > 50\%$ to date, however, the origin of their formation channels is still an open scientific question. Given various properties of BBHs (BH component masses and individual spins) inferred using the default priors by the LVK, independent groups have been trying to explain the formation of the BBHs with different formation channels. Of all formation scenarios, the chemically homogeneous evolution (CHE) channel has stood out with distinguishing features, namely, nearly-equal component masses and preferentially high individual spins aligned with the orbital angular momentum. We perform Bayesian inference on the BBH events officially reported in GWTC-3 with astrophysically-predicted priors representing different formation channels of the isolated binary evolution (CEE: common-envelope evolution channel; CHE; SMT: stable mass transfer). Given assumed models, we report strong evidence for GW190517\_055101 being most likely to have formed through the CHE channel. Assuming the BBH events in the subsample are all formed through one of the isolated binary evolution channels, we obtain the lower limits on the local merger rate density of these channels at $11.45 ~\mathrm{Gpc^{-3}~yr^{-1}}$ (CEE), $0.18 ~\mathrm{Gpc^{-3}~yr^{-1}}$ (CHE), and $0.63 ~\mathrm{Gpc^{-3}~yr^{-1}}$ (SMT) at $90\%$ credible level.
△ Less
Submitted 23 December, 2022; v1 submitted 10 November, 2022;
originally announced November 2022.
-
A Black Hole Kicked At Birth: MAXI J1305-704
Authors:
Chase Kimball,
Sam Imperato,
Vicky Kalogera,
Kyle A. Rocha,
Zoheyr Doctor,
Jeff J. Andrews,
Aaron Dotter,
Emmanouil Zapartas,
Simone S. Bavera,
Konstantinos Kovlakas,
Tassos Fragos,
Phillip M. Srivastava,
Devina Misra,
Meng Sun,
Zepei Xing
Abstract:
When a compact object is formed in a binary, any mass lost during core collapse will impart a kick on the binary's center of mass. Asymmetries in this mass loss or neutrino emission would impart an additional natal kick on the remnant black hole or neutron star, whether it was formed in a binary or in isolation. While it is well established that neutron stars receive natal kicks upon formation, it…
▽ More
When a compact object is formed in a binary, any mass lost during core collapse will impart a kick on the binary's center of mass. Asymmetries in this mass loss or neutrino emission would impart an additional natal kick on the remnant black hole or neutron star, whether it was formed in a binary or in isolation. While it is well established that neutron stars receive natal kicks upon formation, it is unclear whether black holes do as well. Here, we consider the low-mass X-ray binary MAXI J1305-704, which has been reported to have a space velocity $\gtrsim$ 200 km/s. In addition to integrating its trajectory to infer its velocity upon formation of its black hole, we account for recent estimates of its period, black hole mass, mass ratio, and donor effective temperature from photometric and spectroscopic observations. We find that if MAXI J1305-704 formed via isolated binary evolution in the thick Galactic disk, then the supernova that formed its black hole imparted a natal kick of at least 70 km/s while ejecting less than $\simeq 1$ M$_\odot$ with 95% confidence assuming uninformative priors on mass loss and natal kick velocity.
△ Less
Submitted 19 July, 2023; v1 submitted 3 November, 2022;
originally announced November 2022.
-
Investigating the Lower Mass Gap with Low Mass X-ray Binary Population Synthesis
Authors:
Jared C. Siegel,
Ilia Kiato,
Vicky Kalogera,
Christopher P. L. Berry,
Thomas J. Maccarone,
Katelyn Breivik,
Jeff J. Andrews,
Simone S. Bavera,
Aaron Dotter,
Tassos Fragos,
Konstantinos Kovlakas,
Devina Misra,
Kyle A. Rocha,
Philipp M. Srivastava,
Meng Sun,
Zepei Xing,
Emmanouil Zapartas
Abstract:
Mass measurements from low-mass black hole X-ray binaries (LMXBs) and radio pulsars have been used to identify a gap between the most massive neutron stars (NSs) and the least massive black holes (BHs). BH mass measurements in LMXBs are typically only possible for transient systems: outburst periods enable detection via all-sky X-ray monitors, while quiescent periods enable radial-velocity measure…
▽ More
Mass measurements from low-mass black hole X-ray binaries (LMXBs) and radio pulsars have been used to identify a gap between the most massive neutron stars (NSs) and the least massive black holes (BHs). BH mass measurements in LMXBs are typically only possible for transient systems: outburst periods enable detection via all-sky X-ray monitors, while quiescent periods enable radial-velocity measurements of the low-mass donor. We quantitatively study selection biases due to the requirement of transient behavior for BH mass measurements. Using rapid population synthesis simulations (COSMIC), detailed binary stellar-evolution models (MESA), and the disk instability model of transient behavior, we demonstrate that transient-LMXB selection effects introduce observational biases, and can suppress mass-gap BHs in the observed sample. However, we find a population of transient LMXBs with mass-gap BHs form through accretion-induced collapse of a NS during the LMXB phase, which is inconsistent with observations. These results are robust against variations of binary evolution prescriptions. The significance of this accretion-induced collapse population depends upon the maximum NS birth mass $M_\mathrm{ NS, birth-max}$. To reflect the observed dearth of low-mass BHs, COSMIC and MESA models favor $M_\mathrm{ NS, birth-max} \lesssim2M_{\odot}$. In the absence of further observational biases against LMXBs with mass-gap BHs, our results indicate the need for additional physics connected to the modeling of LMXB formation and evolution.
△ Less
Submitted 25 July, 2023; v1 submitted 14 September, 2022;
originally announced September 2022.
-
X-ray luminosity function of high-mass X-ray binaries: Studying the signatures of different physical processes using detailed binary evolution calculations
Authors:
Devina Misra,
Konstantinos Kovlakas,
Tassos Fragos,
Margaret Lazzarini,
Simone S. Bavera,
Bret D. Lehmer,
Andreas Zezas,
Emmanouil Zapartas,
Zepei Xing,
Jeff J. Andrews,
Aaron Dotter,
Kyle A. Rocha,
Philipp M. Srivastava,
Meng Sun
Abstract:
The ever-expanding observational sample of X-ray binaries (XRBs) makes them excellent laboratories for constraining binary evolution theory. Such constraints can be obtained by studying the effects of various physical assumptions on synthetic X-ray luminosity functions (XLFs) and comparing to observed XLFs. In this work, we focus on high-mass XRBs (HMXBs) and study the effects on the XLF of variou…
▽ More
The ever-expanding observational sample of X-ray binaries (XRBs) makes them excellent laboratories for constraining binary evolution theory. Such constraints can be obtained by studying the effects of various physical assumptions on synthetic X-ray luminosity functions (XLFs) and comparing to observed XLFs. In this work, we focus on high-mass XRBs (HMXBs) and study the effects on the XLF of various, poorly-constrained assumptions regarding physical processes such as the common-envelope phase, the core-collapse, and wind-fed accretion. We use the new binary population synthesis code POSYDON, which employs extensive pre-computed grids of detailed stellar structure and binary evolution models, to simulate the evolution of binaries. We generate 96 synthetic XRB populations corresponding to different combinations of model assumptions. The generated HMXB XLFs are feature-rich, deviating from the commonly assumed single-power law. We find a break in our synthetic XLF at luminosity $\sim 10^{38}$ erg s$^{-1}$, similar to observed XLFs. However, we find also a general overabundance of XRBs (up to a factor of $\sim$10 for certain model parameter combinations) driven primarily by XRBs with black hole accretors. Assumptions about the transient behavior of Be-XRBs, asymmetric supernova kicks, and common-envelope physics can significantly affect the shape and normalization of our synthetic XLFs. We find that less well-studied assumptions regarding the circularization of the orbit at the onset of Roche-lobe overflow and criteria for the formation of an X-ray emitting accretion disk around wind-accreting black holes can also impact our synthetic XLFs. Our study reveals the importance of large-scale parameter studies, highlighting the power of XRBs in constraining binary evolution theory.
△ Less
Submitted 14 March, 2023; v1 submitted 12 September, 2022;
originally announced September 2022.
-
The $χ_\mathrm{eff}-z$ correlation of field binary black hole mergers and how 3G gravitational-wave detectors can constrain it
Authors:
Simone S. Bavera,
Maya Fishbach,
Michael Zevin,
Emmanouil Zapartas,
Tassos Fragos
Abstract:
Understanding the origin of merging binary black holes is currently one of the most pressing quests in astrophysics. We show that if isolated binary evolution dominates the formation mechanism of merging binary black holes, one should expect a correlation between the effective spin parameter, $χ_\mathrm{eff}$, and the redshift of the merger, $z$, of binary black holes. This correlation comes from…
▽ More
Understanding the origin of merging binary black holes is currently one of the most pressing quests in astrophysics. We show that if isolated binary evolution dominates the formation mechanism of merging binary black holes, one should expect a correlation between the effective spin parameter, $χ_\mathrm{eff}$, and the redshift of the merger, $z$, of binary black holes. This correlation comes from tidal spin-up systems preferentially forming and merging at higher redshifts due to the combination of weaker orbital expansion from low metallicity stars given their reduced wind mass loss rate, delayed expansion and have smaller maximal radii during the supergiant phase compared to stars at higher metallicity. As a result, these tightly bound systems merge with short inspiral times. Given our fiducial model of isolated binary evolution, we show that the origin of a $χ_\mathrm{eff}-z$ correlation in the detectable LIGO--Virgo binary black hole population is different from the intrinsic population, which will become accessible only in the future by third-generation gravitational-wave detectors such as Einstein Telescope and Cosmic Explorer. Finally, we compare our model predictions with population predictions based on the current catalog of binary black hole mergers and find that current data favor a positive correlation of $χ_\mathrm{eff}-z$ as predicted by our model of isolated binary evolution.
△ Less
Submitted 27 June, 2022; v1 submitted 6 April, 2022;
originally announced April 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.
-
Astrophysics with the Laser Interferometer Space Antenna
Authors:
Pau Amaro Seoane,
Jeff Andrews,
Manuel Arca Sedda,
Abbas Askar,
Quentin Baghi,
Razvan Balasov,
Imre Bartos,
Simone S. Bavera,
Jillian Bellovary,
Christopher P. L. Berry,
Emanuele Berti,
Stefano Bianchi,
Laura Blecha,
Stephane Blondin,
Tamara Bogdanović,
Samuel Boissier,
Matteo Bonetti,
Silvia Bonoli,
Elisa Bortolas,
Katelyn Breivik,
Pedro R. Capelo,
Laurentiu Caramete,
Federico Cattorini,
Maria Charisi,
Sylvain Chaty
, et al. (134 additional authors not shown)
Abstract:
The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery…
▽ More
The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe.
△ Less
Submitted 25 May, 2023; v1 submitted 11 March, 2022;
originally announced March 2022.
-
Suspicious Siblings: The Distribution of Mass and Spin Across Component Black Holes in Isolated Binary Evolution
Authors:
Michael Zevin,
Simone S. Bavera
Abstract:
The LIGO and Virgo gravitational-wave detectors have uncovered binary black hole systems with definitively nonzero spins, as well as systems with significant spin residing in the more massive black hole of the pair. We investigate the ability of isolated binary evolution in forming such highly spinning, asymmetric-mass systems through both accretion onto the first-born black hole and tidal spin-up…
▽ More
The LIGO and Virgo gravitational-wave detectors have uncovered binary black hole systems with definitively nonzero spins, as well as systems with significant spin residing in the more massive black hole of the pair. We investigate the ability of isolated binary evolution in forming such highly spinning, asymmetric-mass systems through both accretion onto the first-born black hole and tidal spin-up of the second-born black hole using a rapid population synthesis approach with detailed considerations of spin-up through tidal interactions. Even with the most optimistic assumptions regarding the efficiency at which an accreting star receives material from a donor, we find that it is difficult to form systems with significant mass asymmetry and moderate or high spins in the primary black hole component. Assuming efficient angular momentum transport within massive stars and Eddington-limited accretion onto black holes, we find that $< 1.5\%$ of systems in the underlying binary black hole population have a primary black hole spin greater than 0.2 and a mass asymmetry of greater than 2:1 in our most optimistic models, with most models finding that this criteria is only met in $\sim 0.01\%$ of systems. The production of systems with significant mass asymmetries and spin in the primary black hole component is thus an unlikely byproduct of isolated evolution unless highly super-Eddington accretion is invoked or angular momentum transport in massive stars is less efficient than typically assumed.
△ Less
Submitted 8 July, 2022; v1 submitted 4 March, 2022;
originally announced March 2022.
-
POSYDON: A General-Purpose Population Synthesis Code with Detailed Binary-Evolution Simulations
Authors:
Tassos Fragos,
Jeff J. Andrews,
Simone S. Bavera,
Christopher P. L. Berry,
Scott Coughlin,
Aaron Dotter,
Prabin Giri,
Vicky Kalogera,
Aggelos Katsaggelos,
Konstantinos Kovlakas,
Shamal Lalvani,
Devina Misra,
Philipp M. Srivastava,
Ying Qin,
Kyle A. Rocha,
Jaime Roman-Garza,
Juan Gabriel Serra,
Petter Stahle,
Meng Sun,
Xu Teng,
Goce Trajcevski,
Nam Hai Tran,
Zepei Xing,
Emmanouil Zapartas,
Michael Zevin
Abstract:
Most massive stars are members of a binary or a higher-order stellar systems, where the presence of a binary companion can decisively alter their evolution via binary interactions. Interacting binaries are also important astrophysical laboratories for the study of compact objects. Binary population synthesis studies have been used extensively over the last two decades to interpret observations of…
▽ More
Most massive stars are members of a binary or a higher-order stellar systems, where the presence of a binary companion can decisively alter their evolution via binary interactions. Interacting binaries are also important astrophysical laboratories for the study of compact objects. Binary population synthesis studies have been used extensively over the last two decades to interpret observations of compact-object binaries and to decipher the physical processes that lead to their formation. Here, we present POSYDON, a novel, binary population synthesis code that incorporates full stellar-structure and binary-evolution modeling, using the MESA code, throughout the whole evolution of the binaries. The use of POSYDON enables the self-consistent treatment of physical processes in stellar and binary evolution, including: realistic mass-transfer calculations and assessment of stability, internal angular-momentum transport and tides, stellar core sizes, mass-transfer rates and orbital periods. This paper describes the detailed methodology and implementation of POSYDON, including the assumed physics of stellar- and binary-evolution, the extensive grids of detailed single- and binary-star models, the post-processing, classification and interpolation methods we developed for use with the grids, and the treatment of evolutionary phases that are not based on pre-calculated grids. The first version of POSYDON targets binaries with massive primary stars (potential progenitors of neutron stars or black holes) at solar metallicity.
△ Less
Submitted 7 August, 2022; v1 submitted 11 February, 2022;
originally announced February 2022.
-
Stochastic gravitational-wave background as a tool to investigate multi-channel astrophysical and primordial black-hole mergers
Authors:
Simone S. Bavera,
Gabriele Franciolini,
Giulia Cusin,
Antonio Riotto,
Michael Zevin,
Tassos Fragos
Abstract:
The formation of merging binary black holes can occur through multiple astrophysical channels such as, e.g., isolated binary evolution and dynamical formation or, alternatively, have a primordial origin. Increasingly large gravitational-wave catalogs of binary black-hole mergers have allowed for the first model selection studies between different theoretical predictions to constrain some of their…
▽ More
The formation of merging binary black holes can occur through multiple astrophysical channels such as, e.g., isolated binary evolution and dynamical formation or, alternatively, have a primordial origin. Increasingly large gravitational-wave catalogs of binary black-hole mergers have allowed for the first model selection studies between different theoretical predictions to constrain some of their model uncertainties and branching ratios. In this work, we show how one could add an additional and independent constraint to model selection by using the stochastic gravitational-wave background. In contrast to model selection analyses that have discriminating power only up to the gravitational-wave detector horizons (currently at redshifts $z\lesssim 1$ for LIGO-Virgo), the stochastic gravitational-wave background accounts for the redshift integration of all gravitational-wave signals in the Universe. As a working example, we consider the branching ratio results from a model selection study that includes potential contribution from astrophysical and primordial channels. We renormalize the relative contribution of each channel to the detected event rate to compute the total stochastic gravitational-wave background energy density. The predicted amplitude lies below the current observational upper limits of GWTC-2 by LIGO-Virgo, indicating that the results of the model selection analysis are not ruled out by current background limits. Furthermore, given the set of population models and inferred branching ratios, we find that, even though the predicted background will not be detectable by current generation gravitational-wave detectors, it will be accessible by third-generation detectors such as the Einstein Telescope and space-based detectors such as LISA.
△ Less
Submitted 21 December, 2021; v1 submitted 13 September, 2021;
originally announced September 2021.
-
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.
-
Revisiting the explodability of single massive star progenitors of stripped-envelope supernovae
Authors:
E. Zapartas,
M. Renzo,
T. Fragos,
A. Dotter,
J. J. Andrews,
S. S. Bavera,
S. Coughlin,
D. Misra,
K. Kovlakas,
J. Román-Garza,
J. G. Serra,
Y. Qin,
K. A. Rocha,
N. H. Tran,
Z. P. Xing
Abstract:
Stripped-envelope supernovae (Types IIb, Ib, and Ic) that show little or no hydrogen comprise roughly one-third of the observed explosions of massive stars. Their origin and the evolution of their progenitors are not yet fully understood. Very massive single stars stripped by their own winds ($\gtrsim 25-30 M_{\odot}$ at solar metallicity) are considered viable progenitors of these events. However…
▽ More
Stripped-envelope supernovae (Types IIb, Ib, and Ic) that show little or no hydrogen comprise roughly one-third of the observed explosions of massive stars. Their origin and the evolution of their progenitors are not yet fully understood. Very massive single stars stripped by their own winds ($\gtrsim 25-30 M_{\odot}$ at solar metallicity) are considered viable progenitors of these events. However, recent 1D core-collapse simulations show that some massive stars may collapse directly into black holes after a failed explosion, with a weak or no visible transient. In this letter, we estimate the effect of direct collapse into a black hole on the rates of stripped-envelope supernovae that arise from single stars. For this, we compute single-star MESA models at solar metallicity and map their final state to their core-collapse outcome following prescriptions commonly used in population synthesis. According to our models, no single stars that have lost their entire hydrogen-rich envelope are able to explode, and only a fraction of progenitors left with a thin hydrogen envelope do (IIb progenitor candidates), unless we use a prescription that takes the effect of turbulence into account or invoke increased wind mass-loss rates. This result increases the existing tension between the single-star paradigm to explain most stripped-envelope supernovae and their observed rates and properties. At face value, our results point toward an even higher contribution of binary progenitors to stripped-envelope supernovae. Alternatively, they may suggest inconsistencies in the common practice of mapping different stellar models to core-collapse outcomes and/or higher overall mass loss in massive stars.
△ Less
Submitted 17 December, 2021; v1 submitted 9 June, 2021;
originally announced June 2021.
-
Approximations to the spin of close Black-hole-Wolf-Rayet binaries
Authors:
Simone S. Bavera,
Michael Zevin,
Tassos Fragos
Abstract:
Population synthesis studies of binary black-hole mergers often lack robust black-hole spin estimates as they cannot accurately follow tidal spin-up during the late black-hole-Wolf-Rayet evolutionary phase. We provide an analytical approximation of the dimensionless second-born black-hole spin given the binary orbital period and Wolf-Rayet stellar mass at helium depletion or carbon depletion. Thes…
▽ More
Population synthesis studies of binary black-hole mergers often lack robust black-hole spin estimates as they cannot accurately follow tidal spin-up during the late black-hole-Wolf-Rayet evolutionary phase. We provide an analytical approximation of the dimensionless second-born black-hole spin given the binary orbital period and Wolf-Rayet stellar mass at helium depletion or carbon depletion. These approximations are obtained from fitting a sample of around $10^5$ detailed MESA simulations that follow the evolution and spin up of close black-hole--Wolf-Rayet systems with metallicities in the range $[10^{-4},1.5Z_\odot]$. Following the potential spin up of the Wolf-Rayet progenitor, the second-born black-hole spin is calculated using up-to-date core collapse prescriptions that account for any potential disk formation in the collapsing Wolf-Rayet star. The fits for second-born black hole spin provided in this work can be readily applied to any astrophysical modeling that relies on rapid population synthesis, and will be useful for the interpretation of gravitational-wave sources using such models.
△ Less
Submitted 19 May, 2021;
originally announced May 2021.
-
The role of core-collapse physics in the observability of black-hole neutron-star mergers as multi-messenger sources
Authors:
Jaime Román-Garza,
Simone S. Bavera,
Tassos Fragos,
Emmanouil Zapartas,
Devina Misra,
Jeff Andrews,
Scotty Coughlin,
Aaron Dotter,
Konstantinos Kovlakas,
Juan Gabriel Serra,
Ying Qin,
Kyle A. Rocha,
Nam Hai Tran
Abstract:
Recent detailed 1D core-collapse simulations have brought new insights on the final fate of massive stars, which are in contrast to commonly used parametric prescriptions. In this work, we explore the implications of these results to the formation of coalescing black-hole (BH) - neutron-star (NS) binaries, such as the candidate event GW190426_152155 reported in GWTC-2. Furthermore, we investigate…
▽ More
Recent detailed 1D core-collapse simulations have brought new insights on the final fate of massive stars, which are in contrast to commonly used parametric prescriptions. In this work, we explore the implications of these results to the formation of coalescing black-hole (BH) - neutron-star (NS) binaries, such as the candidate event GW190426_152155 reported in GWTC-2. Furthermore, we investigate the effects of natal kicks and the NS's radius on the synthesis of such systems and potential electromagnetic counterparts linked to them. Synthetic models based on detailed core-collapse simulations result in an increased merger detection rate of BH-NS systems ($\sim 2.3$ yr$^{-1}$), 5 to 10 times larger than the predictions of "standard" parametric prescriptions. This is primarily due to the formation of low-mass BH via direct collapse, and hence no natal kicks, favored by the detailed simulations. The fraction of observed systems that will produce an electromagnetic counterpart, with the detailed supernova engine, ranges from $2$-$25$%, depending on uncertainties in the NS equation of state. Notably, in most merging systems with electromagnetic counterparts, the NS is the first-born compact object, as long as the NS's radius is $\lesssim 12\,\mathrm{km}$. Furthermore, core-collapse models that predict the formation of low-mass BHs with negligible natal kicks increase the detection rate of GW190426_152155-like events to $\sim 0.6 \, $yr$^{-1}$; with an associated probability of electromagnetic counterpart $\leq 10$% for all supernova engines. However, increasing the production of direct-collapse low-mass BHs also increases the synthesis of binary BHs, over-predicting their measured local merger density rate. In all cases, models based on detailed core-collapse simulation predict a ratio of BH-NSs to binary BHs merger rate density that is at least twice as high as other prescriptions.
△ Less
Submitted 3 December, 2020;
originally announced December 2020.
-
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.
-
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.
-
The origin of spin in binary black holes: Predicting the distributions of the main observables of Advanced LIGO
Authors:
Simone S. Bavera,
Tassos Fragos,
Ying Qin,
Emmanouil Zapartas,
Coenraad J. Neijssel,
Ilya Mandel,
Aldo Batta,
Sebastian M. Gaebel,
Chase Kimball,
Simon Stevenson
Abstract:
We study the formation of coalescing binary black holes via the evolution of isolated field binaries that go through the common envelope phase in order to obtain the combined distributions of the main observables of Advanced LIGO. We used a hybrid technique that combines the parametric binary population synthesis code COMPAS with detailed binary evolution simulations performed with the MESA code.…
▽ More
We study the formation of coalescing binary black holes via the evolution of isolated field binaries that go through the common envelope phase in order to obtain the combined distributions of the main observables of Advanced LIGO. We used a hybrid technique that combines the parametric binary population synthesis code COMPAS with detailed binary evolution simulations performed with the MESA code. We then convolved our binary evolution calculations with the redshift- and metallicity-dependent star-formation rate and the selection effects of gravitational-wave detectors to obtain predictions of observable properties. By assuming efficient angular momentum transport, we are able to present a model that is capable of simultaneously predicting the three main gravitational-wave observables: the effective inspiral spin parameter $χ_{eff}$, the chirp mass $M_{chirp}$, and the cosmological redshift of merger $z_{merger}$. We find an excellent agreement between our model and the ten events from the first two advanced detector observing runs. We make predictions for the third observing run O3 and for Advanced LIGO design sensitivity. We expect around 80% of events with $χ_{eff} < 0.1$, while the remaining 20% of events with $χ_{eff} \ge 0.1$ are split into ~10% with $M_{chirp} < 15$ M$_\odot$ and ~10% with $M_{chirp} \ge 15$ M$_\odot$. In conclusion, the favorable comparison of the existing LIGO/Virgo observations with our model predictions gives support to the idea that the majority, if not all of the observed mergers, originate from the evolution of isolated binaries. The first-born black hole has negligible spin because it lost its envelope after it expanded to become a giant star, while the spin of the second-born black hole is determined by the tidal spin up of its naked helium star progenitor by the first-born black hole companion after the binary finished the common-envelope phase.
△ Less
Submitted 24 January, 2020; v1 submitted 28 June, 2019;
originally announced June 2019.