We study the impact of mass-transfer physics on the observable properties of binary black hole populations that formed through isolated binary evolution. We used the
POSYDON
framework to combine ...detailed
MESA
binary simulations with the
COSMIC
population synthesis tool to obtain an accurate estimate of merging binary black hole observables with a specific focus on the spins of the black holes. We investigate the impact of mass-accretion efficiency onto compact objects and common-envelope efficiency on the observed distributions of the effective inspiral spin parameter
χ
eff
, chirp mass
M
chirp
, and binary mass ratio
q
. We find that low common envelope efficiency translates to tighter orbits following the 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-wave detectors as they form and merge at high redshifts (
z
∼ 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 that the local rate density (
z
≃ 0.01) for the common envelope channel is in the range of ∼17–113 Gpc
−3
yr
−1
, considering a range of
α
CE
∈ 0.2, 5.0, while for the stable mass transfer channel the rate density is ∼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 constrained the lower bound of branching fraction from other formation channels in the detected population to be ∼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.
Abstract
Most massive stars are members of a binary or a higher-order stellar system, 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, publicly available, 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 postprocessing, classification, and interpolation methods we developed for use with the grids, and the treatment of evolutionary phases that are not based on precalculated grids. The first version of
POSYDON
targets binaries with massive primary stars (potential progenitors of neutron stars or black holes) at solar metallicity.
A Black Hole Kicked at Birth: MAXI J1305-704 Kimball, Chase; Imperato, Sam; Kalogera, Vicky ...
Astrophysical journal. Letters,
08/2023, Letnik:
952, Številka:
2
Journal Article
Recenzirano
Odprti dostop
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 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 ≳200 km s
−1
. 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
−1
while ejecting less than ≃1
M
⊙
with 95% confidence assuming uninformative priors on mass loss and natal kick velocity.
Abstract
The discovery of periodicity in the arrival times of the fast radio bursts (FRBs) poses a challenge to the oft-studied magnetar scenarios. However, models that postulate that FRBs result ...from magnetized shocks or magnetic reconnection in a relativistic outflow are not specific to magnetar engines; instead, they require only the impulsive injection of relativistic energy into a dense magnetized medium. Motivated thus, we outline a new scenario in which FRBs are powered by short-lived relativistic outflows (“flares”) from accreting black holes or neutron stars, which propagate into the cavity of the pre-existing (“quiescent”) jet. In order to reproduce FRB luminosities and rates, we are driven to consider binaries of stellar-mass compact objects undergoing super-Eddington mass transfer, similar to ultraluminous X-ray (ULX) sources. Indeed, the host galaxies of FRBs, and their spatial offsets within their hosts, show broad similarities with ULXs. Periodicity on timescales of days to years could be attributed to precession (e.g., Lens-Thirring) of the polar accretion funnel, along which the FRB emission is geometrically and relativistically beamed, which sweeps across the observer line of sight. Accounting for the most luminous FRBs via accretion power may require a population of binaries undergoing brief-lived phases of unstable (dynamical-timescale) mass transfer. This will lead to secular evolution in the properties of some repeating FRBs on timescales of months to years, followed by a transient optical/IR counterpart akin to a luminous red nova, or a more luminous accretion-powered optical/X-ray transient. We encourage targeted FRB searches of known ULX sources.
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 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 an 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
NS
,
birth
−
max
. To reflect the observed dearth of low-mass BHs,
COSMIC
and
MESA
models favor
M
NS
,
birth
−
max
≲
2
M
⊙
. 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.
Aims.
We use
N
-body simulations to examine whether a characteristic turnaround radius, as predicted from the spherical collapse model in a ΛCDM Universe, can be meaningfully identified for galaxy ...clusters in the presence of full three-dimensional effects.
Methods.
We use The Dark Sky Simulations and Illustris-TNG dark-matter-only cosmological runs to calculate radial velocity profiles around collapsed structures, extending out to many times the virial radius
R
200
. There, the turnaround radius can be unambiguously identified as the largest nonexpanding scale around a center of gravity.
Results.
We find that: (a) a single turnaround scale can meaningfully describe strongly nonspherical structures. (b) For halos of masses
M
200
> 10
13
M
⊙
, the turnaround radius
R
ta
scales with the enclosed mass
M
ta
as
M
ta
1/3
, as predicted by the spherical collapse model. (c) The deviation of
R
ta
in simulated halos from the spherical collapse model prediction is relatively insensitive to halo asphericity. Rather, it is sensitive to the tidal forces due to massive neighbors when these are present. (d) Halos exhibit a characteristic average density within the turnaround scale. This characteristic density is dependent on cosmology and redshift. For the present cosmic epoch and for concordance cosmological parameters (Ω
m
∼ 0.3; Ω
Λ
∼ 0.7) turnaround structures exhibit a density contrast with the matter density of the background Universe of
δ
∼ 11. Thus,
R
ta
is equivalent to
R
11
– in a way that is analogous to defining the “virial” radius as
R
200
– with the advantage that
R
11
is shown in this work to correspond to a kinematically relevant scale in
N
-body simulations.
In this work, we present detailed constraints on the metallicity dependence of the high-mass X-ray binary (HMXB) X-ray luminosity function (XLF). We analyze 5 Ms of Chandra data for 55 actively ...star-forming galaxies at D 30 Mpc, with gas-phase metallicities spanning 7-9.2. Within the galactic footprints, our sample contains a total of 1311 X-ray point sources, of which 49% are expected to be HMXBs, with the remaining sources likely to be low-mass X-ray binaries (LMXBs; 22%) and unrelated background sources ( 29%). We construct a model that successfully characterizes the average HMXB XLF over the full metallicity range. We demonstrate that the SFR-normalized HMXB XLF shows clear trends with metallicity, showing steadily increasing numbers of luminous and ultraluminous X-ray sources ( (erg s−1) = 38-40.5) with declining metallicity. However, we find that the low-luminosity ( (erg s−1) = 36-38) HMXB XLF appears to show a nearly constant SFR scaling and slope with metallicity. Our model provides a revised scaling relation of integrated LX/SFR versus , and a new characterization of its SFR-dependent stochastic scatter. The general trend of this relation is broadly consistent with past studies based on integrated galaxy emission; however, our model suggests that this relation is driven primarily by the high-luminosity end of the HMXB XLF. Our results have implications for binary population synthesis models, the nature of super-Eddington accreting objects (e.g., ultraluminous X-ray sources), recent efforts to identify active galactic nucleus candidates in dwarf galaxies, and the X-ray radiation fields in the early universe during the epoch of cosmic heating at z 10.
Context.
Many physical processes taking place during the evolution of binary stellar systems remain poorly understood. The ever-expanding observational sample of X-ray binaries (XRBs) makes them ...excellent laboratories for constraining binary evolution theory. Such constraints and useful insights can be obtained by studying the effects of various physical assumptions on synthetic X-ray luminosity functions (XLFs) and comparing them with observed XLFs.
Aims.
In this work we focus on high-mass X-ray binaries (HMXBs) and study the effects on the XLF of various, poorly constrained assumptions regarding physical processes, such as the common-envelope phase, core collapse, and wind-fed accretion.
Methods.
We used the new binary population synthesis code
POSYDON
, which employs extensive precomputed grids of detailed stellar structure and binary evolution models, to simulate the entire evolution of binaries. We generated 96 synthetic XRB populations corresponding to different combinations of model assumptions, including different prescriptions for supernova kicks, supernova remnant masses, common-envelope evolution, circularization at the onset of Roche-lobe overflow, and observable wind-fed accretion.
Results.
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 ∼10
38
erg s
−1
, similar to observed XLFs. However, we also find a general overabundance of XRBs (up to a factor of ∼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 and reduce the discrepancy with observations.
Conclusions.
Our synthetic XLFs do not always agree well with observations, especially at intermediate X-ray luminosities, which is likely due to uncertainties in the adopted physical assumptions. While some model parameters leave distinct imprints on the shape of the synthetic XLFs and can reduce this deviation, others do not have a significant effect overall. Our study reveals the importance of large-scale parameter studies, highlighting the power of XRBs in constraining binary evolution theory.
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 subpopulation 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
B
∼ 0.05 our model predicts a rate density comparable to the observed one, throughout the redshift range, while, at redshift
z
∈ 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.
Abstract
Recent 1D core-collapse simulations indicate a nonmonotonicity of the explodability of massive stars with respect to their precollapse core masses, which is in contrast to commonly used ...prescriptions. In this work, we explore the implications of these results on the formation of coalescing black hole (BH)–neutron star (NS) binaries. Furthermore, we investigate the effects of natal kicks and the NS’s radius on the synthesis of such systems and potential electromagnetic counterparts (EMCs) linked to them. Models based on 1D core-collapse simulations result in a BH–NS merger detection rate ( ∼ 2.3 yr
−1
), 5–10 times larger than the predictions of “standard” prescriptions. This is primarily due to the formation of low-mass BHs via direct collapse, and hence no natal kicks, favored by the 1D simulations. The fraction of observed systems that will produce an EMC, with the supernova engine from 1D simulations, ranges from 2% to 25%, depending on the NS equation of state. Notably, in most merging systems with EMCs, the NS is the first-born compact object, as long as the NS’s radius is ≲ 12 km. Furthermore, models with negligible kicks for low-mass BHs increase the detection rate of GW190426_152155-like events to ∼ 0.6 yr
−1
, with an associated probability of EMC ≤10% for all supernova engines. Finally, models based on 1D core-collapse simulations predict a ratio of BH–NSs to binary BHs’ merger rate density that is at least twice as high as other prescriptions, but at the same time overpredicting the measured local merger density rate of binary black holes.