11 research outputs found
Exploring the nature of ultra-luminous X-ray sources across stellar population ages using detailed binary evolution calculations
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.Comment: 22 pages, 12 figures, Accepted by A&A. Fixed typos and updated
references. Referee's comments were addresse
Active Learning for Computationally Efficient Distribution of Binary Evolution Simulations
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
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.Comment: 20 pages (16 main text), 10 figures, submitted to Ap
A Black Hole Kicked At Birth: MAXI J1305-704
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 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. In addition to integrating its trajectory to infer its
velocity upon formation of its black hole, we reconstruct its evolutionary
history, accounting 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 its black hole received a natal kick
of at least 70 km/s with 95\% confidence.Comment: To be submitted; 9 pages, 5 figure
X-ray luminosity function of high-mass X-ray binaries: Studying the signatures of different physical processes using detailed binary evolution calculations
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 erg s, similar to observed XLFs. However, we
find also 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. Our study reveals the importance of large-scale parameter studies,
highlighting the power of XRBs in constraining binary evolution theory.Comment: 31 pages, 32 figures, Accepted by A&A. Fixed typos and updated
references. Referee's comments were addresse
Investigating the Lower Mass Gap with Low Mass X-ray Binary Population Synthesis
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 . To reflect the observed dearth of low-mass BHs, COSMIC and MESA
models favor . 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.Comment: 21 pages, accepted to Ap
POSYDON: A General-Purpose Population Synthesis Code with Detailed Binary-Evolution Simulations
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.Comment: 60 pages, 33 figures, 8 tables, referee's comments addressed. The
code and the accompanying documentations and data products are available at
https:\\posydon.or
The formation of merging black holes with masses beyond 30 M⊙ at solar metallicity
Gravitational-wave astronomy has revealed a population of stellar-mass black holes more massive than observed previously by other means. The maximum mass of black holes formed in isolated binaries is determined by stellar winds, mixing processes and interactions between the binary components. We consider the impact that fully self-consistent, detailed stellar-structure and binary-evolution calculations have on the population synthesis of black-hole binaries at solar metallicity. We find a qualitatively different picture from previous studies employing rapid population-synthesis techniques. Merging binary black holes form with a non-negligible rate (∼4×10−7M−1⊙
) and their progenitor stars with initial masses ≳50 M⊙ do not expand to supergiant radii, thereby largely avoiding substantial dust-driven or luminous blue variable winds. The progenitor stars lose less mass in winds, which results in black holes as massive as ~30 M⊙, and approximately half avoid a mass-transfer episode before forming the first-born black hole. Binaries with initial periods of a few days, some of which may undergo Roche-lobe overflow mass transfer, result in mildly spinning first-born black holes, χBH1 ≲ 0.2, assuming efficient angular-momentum transport
Active learning for computationally efficient distribution of binary evolution simulations
Binary stars undergo a variety of interactions and evolutionary phases, critical for predicting and explaining observations. Binary population synthesis with full simulation of stellar structure and evolution is 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 that are 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 target 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 ∼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 new population synthesis codes such as POSYDON to cover more input parameters while preserving interpolation accuracies
Probing the progenitors of spinning binary black-hole mergers with long gamma-ray bursts
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 fB ∼ 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