Simone S. Bavera

With over two dozen detections in the Milky Way, double neutron stars (DNSs) provide a unique window into massive binary evolution. We use the POSYDON binary population synthesis code to model DNS populations and compare them to the observed Galactic sample. By tracing their origins to underlying single and binary star physics, we place constraints on the detailed evolutionary stages leading to DNS formation. Our study reveals a bifurcation within the well-known common envelope (CE) formation channel for DNSs, which naturally explains an observed split in the orbital periods of the Galactic systems. The two subchannels are defined by whether the donor star has a helium core (Case B mass transfer) or a carbon-oxygen core (Case C) at the onset of the CE, with only the helium core systems eventually merging due to gravitational-wave-modulated orbital decay. We find that across different treatments of the CE phase, the formation of DNSs through both subchannels requires either a generous core definition of ≃30% H-fraction or a high CE ejection efficiency of α _CE ≳ 1.2. By testing different supernova kick velocity models, we find that galactic DNSs are best reproduced using a prescription that favors low velocity kicks (≲50 km s ⁻¹ ), in agreement with previous studies. Furthermore, our models indicate that merging DNSs are born from a stripped progenitor with a median pre-supernova envelope mass ∼0.2 M _⊙ . Our results highlight the value of detailed evolutionary models for improving our understanding of exotic binary star formation.

Whether considering rare astrophysical events on cosmological scales or unresolved stellar populations, accurate models must account for the integrated contribution from the entire history of star formation upon which that population is built. Here, we describe the second version of POSYDON, an open-source binary population synthesis code based on extensive grids of detailed binary evolution models computed using the MESA code, which follows both stars’ structures as a binary system evolves through its complete evolution from the zero-age main sequence, through multiple phases of mass transfer and supernovae, to their death as compact objects. To generate synthetic binary populations, POSYDON uses advanced methods to interpolate between our large, densely spaced grids of simulated binaries. In our updated version of POSYDON, we account for the evolution of stellar binaries across a cosmological range of metallicities, extending from 10⁻⁴ to 2 Z _⊙, including grids specifically focused on the Small and Large Magellanic Clouds (0.2 and 0.45 Z _⊙). In addition to describing our model grids and detailing our methodology, we outline several improvements to POSYDON. These include the incorporation of single stars in stellar populations, a treatment for stellar mergers, and a careful modeling of “reverse-mass transferring” binaries in which a once-accreting star later becomes a donor star. Our simulations are focused on binaries with at least one high-mass component, such as those that host neutron stars and black holes, and we provide postprocessing methods to account for the cosmological evolution of metallicity and star formation as well as rate calculations for transient events.

Context. Long gamma-ray bursts (LGRBs) are generally observed in low-metallicity environments. However, 10% to 20% of LGRBs at redshift z < 2 are associated with near-solar to super-solar metallicity environments, remaining unexplained by traditional LGRB formation pathways that favor low metallicity progenitors. Aims. In this work, we propose a novel formation channel for LGRBs that is dominant at high metallicities. We explore how a stripped primary star in a binary can be spun up by a second stable reverse-mass-transfer phase, initiated by the companion star. Methods. We used POSYDON, a state-of-the-art population synthesis code that incorporates detailed single- and binary-star mode grids, to investigate the metallicity dependence of the stable reverse-mass-transfer LGRB formation channel. We determine the available energy to power an LGRB from the rotational profile and internal structure of a collapsing star and investigated how the predicted rate density of the proposed channel changes with different star formation histories and criteria for defining a successful LGRB. Results. Stable reverse mass transfer can produce rapidly rotating, stripped stars at collapse. These stars retain enough angular momentum to account for approximately 10%–20% of the observed local LGRB rate density, under a reasonable assumption for the definition of a successful LGRB. However, the local rate density of LGRBs from stable reverse mass transfer can vary significantly, between 1 and 100 Gpc⁻³ yr⁻¹, due to strong dependencies on cosmic star formation rate and metallicity evolution, as well as the assumed criteria for successful LGRBs.

Modeling of large populations of binary stellar systems is an integral part of many areas of astrophysics, from radio pulsars and supernovae to X-ray binaries, gamma-ray bursts, and gravitational-wave mergers. Binary population synthesis codes that employ self-consistently the most advanced physics treatment available for stellar interiors and their evolution and are at the same time computationally tractable have started to emerge only recently. One element that is still missing from these codes is the ability to generate the complete time evolution of binaries with arbitrary initial conditions using precomputed three-dimensional grids of binary sequences. Here, we present a highly interpretable method, from binary evolution track interpolation. Our method implements simulation generation from irregularly sampled time series. Our results indicate that this method is appropriate for applications within binary population synthesis and computational astrophysics with time-dependent simulations in general. Furthermore, we point out and offer solutions to the difficulty surrounding evaluating the performance of signals exhibiting extreme morphologies akin to discontinuities.

We investigate Roche lobe overflow mass transfer (MT) in eccentric binary systems between stars and compact objects (COs), modeling the coupled evolution of both the star and the orbit due to eccentric MT (eMT) in a self-consistent framework. We implement the analytic expressions for secular rates of change of the orbital semimajor axis and eccentricity, assuming a delta function MT at periapse, in the binary stellar evolution code MESA. Two scenarios are examined: (1) a simplified model isolating the effects of eMT on stellar and orbital evolution, and (2) realistic binary configurations that include angular momentum exchange (e.g., tides, mass loss, spin–orbit coupling, and gravitational-wave (GW) radiation). Unlike the ad hoc approach of instant circularization that is often employed, explicit modeling of eMT reveals that a large fraction of binaries can remain eccentric post-MT. Even binaries that naturally circularize during eMT have different properties (donor mass and orbital size) compared to predictions from instant circularization, with some showing fundamentally different evolutionary outcomes (e.g., stable vs. unstable MT). We demonstrate that a binary’s initial mass ratio and eccentricity are predictive of whether it will remain eccentric or circularize after eMT. These findings underscore the importance of eMT in understanding CO-hosting binary populations, including X-ray binaries, GW sources, and other high-energy transients.

Binaries harboring a millisecond pulsar (MSP) and a black hole (BH) are a key observing target for current and upcoming pulsar surveys. We model the formation and evolution of such binaries in isolation at solar metallicity using the next-generation binary population synthesis code POSYDON. We examine neutron star (NS)–BH binaries where the NS forms first (labeled NSBH), as the NS must be able to spin up to MSP rotation periods before the BH forms in these systems. We find that NSBHs are very rare and have a birth rate −1 for a Milky Way–like galaxy in our typical models. The NSBH formation rate is 2–3 orders of magnitude smaller than that for NS–BH binaries where the BH forms first (labeled BHNS). These rates are also sensitive to model assumptions about the supernova (SN) remnant masses, natal kicks, metallicity, and common-envelope (CE) evolution parameters. We find that 100% of NSBHs undergo a mass ratio reversal before the first SN and up to 52% of NSBHs undergo a double CE phase after the mass ratio reversal occurs. Most importantly, no NSBH binaries in our populations undergo a mass transfer phase, either stable or unstable, after the first SN. This implies that there is no possibility of pulsar spin-up via accretion, and thus MSP–BH binaries cannot form. Thus, dynamical environments and processes may provide the only formation channels for such MSP–BH binaries.

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_⊙ 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.

Binaries harboring a millisecond pulsar (MSP) and a black hole (BH) are a key observing target for current and upcoming pulsar surveys. We model the formation and evolution of such binaries in isolation at solar metallicity using the next-generation binary population synthesis code POSYDON. We examine neutron star (NS)-BH binaries where the NS forms first (labeled NSBH), as the NS must be able to spin-up to MSP rotation periods before the BH forms in these systems. We find that NSBHs are very rare and have a birth rate < 1 Myr$^{-1}$ for a Milky Way-like galaxy in our typical models. The NSBH birth rate is 2-3 orders of magnitude smaller than that for NS-BHs where the BH forms first (labeled BHNS). These rates are also sensitive to model assumptions about the supernova (SN) remnant masses, natal kicks, and common-envelope efficiency. We find that 100% of NSBHs undergo a mass ratio reversal before the first SN and up to 64% of NSBHs undergo a double common envelope phase after the mass ratio reversal occurs. Most importantly, no NSBH binaries in our populations undergo a mass transfer phase, either stable or unstable, after the first SN. This implies that there is no possibility of pulsar spin-up via accretion, and thus MSP-BH binaries cannot form. Thus, dynamical environments and processes may provide the only formation channels for such MSP-BH binaries.