Binary Mergers Research Articles

Compact object populations over cosmic time II. Compact object merger rates and masses over redshift from varying initial conditions

ABSTRACT We perform a first study of the impact of varying two components of the initial conditions in binary population synthesis of compact binary mergers – the initial mass function, which is made metallicity- and star formation rate-dependent, and the orbital parameter (orbital period, mass ratio, and eccentricity) distributions, which are assumed to be correlated – within a larger grid of initial condition models also including alternatives for the primary mass-dependent binary fraction and the metallicity-specific cosmic star formation history. We generate the initial populations with the sampling code bossa and evolve them with the rapid population synthesis code compas. We find strong suggestions that the main role of initial conditions models is to set the relative weights of key features defined by the evolution models. In the two models we compare, black hole–black hole (BHBH) mergers are the most strongly affected, which we connect to a shift from the common envelope to the stable Roche lobe overflow formation channels with decreasing redshift. We also characterize variations in the black hole–neutron star (BHNS) and neutron star–neutron star (NSNS) final parameter distributions. We obtain the merger rate evolution for BHBH, BHNS, and NSNS mergers up to $z=10$, and find a variation by a factor of $\sim 50\textnormal {--}60$ in the local BHBH and BHNS merger rates, suggesting a more important contribution from initial conditions than previously thought, and calling for a complete exploration of the initial conditions model permutations.

Constraining the AGN formation channel for detected black hole binary mergers up to z=1.5 with the Quaia catalogue

Constraining the AGN formation channel for detected black hole binary mergers up to z=1.5 with the Quaia catalogue

Thermal Conductivity and Thermal Hall Effect in Dense Electron-Ion Plasma

In this study, we examine thermal conductivity and the thermal Hall effect in electron-ion plasmas relevant to hot neutron stars, white dwarfs, and binary neutron star mergers, focusing on densities found in the outer crusts of neutron stars and the interiors of white dwarfs. We consider plasma consisting of single species of ions, which could be either iron [56]Fe or carbon [12]C nuclei. The temperature range explored is from the melting temperature of the solid T∼109 K up to 1011 K. This covers both degenerate and non-degenerate electron regimes. We find that thermal conductivity increases with density and temperature for which we provide analytical scaling relations valid in different regimes. The impact of magnetic fields on thermal conductivity is also analyzed, showing anisotropy in low-density regions and the presence of the thermal Hall effect characterized by the Righi–Leduc coefficient. The transition froma degenerate to non-degenerate regime is characterized by a minimum ratio of thermal conductivity to temperature, which is analogous to the minimum observed already in the case of electrical conductivity. We provide also formulas fit to our numerical results, which can be used in dissipative magneto-hydrodynamics simulations of warm compact stars.

Binary neutron star mergers in massive scalar-tensor theory: Properties of postmerger remnants

We investigate the properties of postmerger remnants of binary neutron star mergers in the framework of Damour–Esposito-Farese-type scalar-tensor theory of gravity with a massive scalar field by numerical relativity simulation. It is found that the threshold mass for prompt collapse is raised in the presence of the excited scalar field. Our simulation results also suggest the existence of a long-lived ϕ mode in hypermassive neutron stars due to the presence of the massive scalar field that enhances the quasiradial oscillation in the remnant. We investigate the descalarization condition in hypermassive neutron stars and discover a distinctive signature in postmerger gravitational waves. Published by the American Physical Society 2024

Incidence of afterglow plateaus in gamma-ray bursts associated with binary neutron star mergers

One of the most surprising gamma-ray burst (GRB) features discovered with the Swift X-ray telescope (XRT) is a plateau phase in the early X-ray afterglow light curves. These plateaus are observed in the majority of long GRBs, while their incidence in short GRBs (SGRBs) is still uncertain due to their fainter X-ray afterglow luminosity with respect to long GRBs. An accurate estimate of the fraction of SGRBs with plateaus is of utmost relevance given the implications that the plateau may have for our understanding of the jet structure and possibly of the nature of the binary neutron star (BNS) merger remnant. This work presents the results of an extensive data analysis of the largest and most up-to-date sample of SGRBs observed with the XRT, and for which the redshift has been measured. We find a plateau incidence of 18-37<!PCT!> in SGRBs, which is a significantly lower fraction than that measured in long GRBs (>50<!PCT!>). Although still debated, the plateau phase could be explained as energy injection from the spin-down power of a newly born magnetized neutron star (NS; magnetar). We show that this scenario can nicely reproduce the observed short GRB (SGRBs) plateaus, while at the same time providing a natural explanation for the different plateau fractions between short and long GRBs. In particular, our findings may imply that only a minority of BNS mergers generating SGRBs leave behind a sufficiently stable or long-lived NS to form a plateau. From the probability distribution of the BNS remnant mass, a fraction 18-37<!PCT!> of short GRB plateaus implies a maximum NS mass in the range $ odot

Kinematic Insights into Luminous Blue Variables and B[e] Supergiants

Abstract Recent work suggests that many luminous blue variables (LBVs) and B[e] supergiants (sgB[e]) are isolated, implying that they may be products of massive binaries, kicked by partner supernovae. However, the evidence is somewhat complex and controversial. To test this scenario, we measure the proper-motion velocities for these objects in the LMC and SMC, using Gaia Data Release 3. Our LMC results show that the kinematics, luminosities, and IR properties point to LBVs and sgB[e] stars being distinct classes. We find that Class 1 LBVs, which have dusty nebulae, and sgB[e] stars both show velocity distributions comparable to that of SMC field OBe stars, which are known to have experienced SN kicks. The sgB[e] stars are faster, plausibly due to their lower average masses. However, Class 2 LBVs, which are luminous objects without dusty nebulae, show no signs of acceleration, therefore suggesting that they are single stars, pre-SN binaries, or perhaps binary mergers. The candidate LBV Class 3 stars, which are dominated by hot dust, are all confirmed sgB[e] stars; their luminosities and velocities show that they simply represent the most luminous and massive of the sgB[e] class. There are very few SMC objects, but the sgB[e] stars are faster than their LMC counterparts, which may be consistent with expectations that lower-metallicity binaries are tighter, causing faster ejections. We also examine the distinct class of dust-free, weak-lined sgB[e] stars, finding that the SMC objects have the fastest velocities of the entire sample.

Simulating Short Gamma-Ray Burst Jets in Realistic Late Binary Neutron Star Merger Environments

The electromagnetic emission and the afterglow observations of the binary neutron star merger event GW170817A confirmed the association of the merger with a short gamma-ray burst (GRB) harboring a narrow (5°–10°) and powerful (1049–1050 erg) jet. Using the 1 s long neutrino-radiation general relativistic MHD simulation of coalescing neutron stars of K. Kiuchi et al., and following the semi-analytical estimates of M. Pais et al., we inject a narrow, powerful, unmagnetized jet into the post-merger phase. We explore different opening angles, luminosities, central engine durations, and times after the merger. We explore early (0.1 s following the merger) and late (1 s) jet launches; the latter is consistent with the time delay of ≈1.74 s observed between GW170817 and GRB 170817A. We demonstrate that the semi-analytical estimates correctly predict the jets’ breakout and collimation conditions. When comparing our synthetic afterglow light curves to the observed radio data of GW170807, we find a good agreement for a 3 × 1049 erg jet launched late with an opening angle in the range ≃5°–7°.

Dual Active Galactic Nuclei: Precursors of Binary Supermassive Black Hole Formation and Mergers

The presence of dual active galactic nuclei (AGNs) on scales of a few tens of kiloparsecs can be used to study merger-induced accretion on supermassive black holes (SMBHs) and offer insights about SMBH mergers, using dual AGNs as merger precursors. This study uses the Romulus25 cosmological simulation to investigate the properties and evolution of dual AGNs. We first analyze the properties of AGNs (L bol > 1043 erg s−1) and their neighboring SMBHs (any SMBHs closer than 30 pkpc to an AGN) at z ≤ 2. This is our underlying population. We then applied the luminosity threshold of L bol > 1043 erg s−1 to the neighboring SMBHs thereby identifying dual and multiple AGNs. Our findings indicate an increase in the number of both single and dual AGNs from lower to higher redshifts. We also find that the number of dual AGNs with separations of 0.5–4 kpc is twice the number of duals with separations of 4–30 kpc. All dual AGNs in our sample resulted from major mergers. Compared to single AGNs, duals have a lower black hole-to-halo mass ratio. We found that the properties of dual AGN host halos, including halo mass, stellar mass, star formation rate, and gas mass, are generally consistent with those of single AGN halos, albeit tending toward the higher end of their respective property ranges. Our analysis uncovered a diverse array of evolutionary patterns among dual AGNs, including rapidly evolving systems, slower ones, and instances where SMBH mergers are ineffective.

Unveiling the Progenitors of a Population of Likely Peculiar Gamma-Ray Bursts

Traditionally, gamma-ray bursts (GRBs) are classified as long and short GRBs, with T 90 = 2 s being the threshold duration. Generally, long-duration GRBs (LGRBs; T 90 > 2 s) are associated with the collapse of massive stars, and short-duration (SGRBs; T 90 < 2 s) are associated with compact binary mergers involving at least one neutron star. However, the existence of a population of so-called “peculiar GRBs”—i.e., LGRBs originating from mergers or long Type I GRBs, and SGRBs originating from collapsars or short Type II GRBs—has challenged the traditional paradigm of GRB classification. Finding more peculiar GRBs may help to give more insight into this issue. In this work, we analyze the properties of machine-learning-identified long Type I GRB and short Type II GRB candidates, long GRBs-I and short GRBs-II (the so-called “peculiar GRBs”). We find that long GRBs-I almost always exhibit properties similar to Type I GRBs, which suggests that mergers may indeed produce GRBs with T 90 > 2 s. Furthermore, according to the probability given by the redshift distribution, short GRBs-II almost exhibit properties similar to Type II GRBs. This suggests that the populations of short Type II GRBs are not scarce and that they are hidden in a large number of samples without redshifts, which is unfavorable for the interpretation that the jet progression leads to a missed main emission.

Electromagnetic signatures of black hole clusters in the center of super-Eddington galaxies

Supermassive black holes (SMBHs) at the centers of active galaxies are fed by accretion disks that radiate from the infrared or optical to the X-ray bands. Several types of objects can orbit SMBHs, including massive stars, neutron stars, clouds from the broad- and narrow-line regions, and X-ray binaries. Isolated black holes with a stellar origin (BHs of ∼10 M⊙) should also be present in large numbers within the central parsec of the galaxies. These BHs are expected to form a cluster around the SMBH as a result of the enhanced star formation rate in the inner galactic region and the BH migration caused by gravitational dynamical friction. However, except for occasional microlensing effects on background stars or gravitational waves from binary BH mergers, the presence of a BH population is hard to verify. In this paper, we explore the possibility of detecting electromagnetic signatures of a central cluster of BHs when the accretion rate onto the central SMBH is greater than the Eddington rate. In these supercritical systems, the accretion disk launches powerful winds that interact with the objects orbiting the SMBH. Isolated BHs can capture matter from this dense wind, leading to the formation of small accretion disks around them. If jets are produced in these single microquasars, they could be sites of particle acceleration to relativistic energies. These particles in turn are expected to cool by various radiative processes. Therefore, the wind of the SMBH might illuminate the BHs through the production of both thermal and nonthermal radiation.

Hydrodynamic 3D Simulation of Roche Lobe Overflow in High-mass X-Ray Binaries

Abstract While binary merger events have been an active area of study in both simulations and observational work, the formation channels by which a high-mass star extends from Roche lobe overflow (RLO) in a decaying orbit of a black-hole (BH) companion to a binary black-hole (BBH) system merits further investigation. Variable length-scales must be employed to accurately represent the dynamical fluid transfer and morphological development of the primary star as it conforms to a diminishing Roche lobe under the runaway influence of the proximal BH. We have simulated and evolved binary mass flow under these conditions to better identify the key transitional processes from RLO to BBHs. We demonstrate a new methodology to model RLO systems to unprecedented resolution simultaneously across the envelope, donor wind, tidal stream, and accretion disk regimes without reliance upon previously universal symmetry, mass flux, and angular momentum flux assumptions. We have applied this method to the semidetached high-mass X-ray binary M33 X-7 in order to provide a direct comparison to recent observations of an RLO candidate system at two overflow states of overfilling factors f = 1.01 and f = 1.1. We found extreme overflow (f = 1.1) to be entirely conservative in both mass and angular momentum transport, forming a conical L1 tidal stream of density and deflected angle comparable to existing predictions. This case lies within the unstable mass transfer (MT) regime as recently proposed of M33 X-7. The f = 1.01 case differed in stream geometry, accretion disk size, and efficiency, demonstrating nonconservative stable MT through a ballistic uniform-width stream. The nonconservative and stable nature of the f = 1.01 case MT also suggests that existing assumptions of semidetached binaries undergoing RLO may mischaracterize their role and distribution as progenitors of BBHs and common envelopes.

A dark standard siren measurement of the Hubble constant following LIGO/Virgo/KAGRA O4a and previous runs

Abstract We present a new constraint on the Hubble constant (H0) from the standard dark siren method using a sample of 5 well-covered gravitational waves (GW) alerts reported during the first part of the fourth LIGO/Virgo/KAGRA observing run and with 3 updated standard dark sirens from third observation run in combination the previous constraints from the first three runs. Our methodology relies on the galaxy catalogue method alone. We use a Deep Learning method to derive the full probability density estimation of photometric redshifts using the Legacy Survey catalogues. We add the constraints from well localized Binary Black Hole mergers to the sample of standard dark sirens analysed in our previous work. We combine the H0 posterior for 5 new standard sirens with other 10 previous events (using the most recent available data for the 5 novel events and updated 3 previous posteriors from O3), finding $H_0 = 70.4^{+13.6}_{-11.7}~{\rm km~s^{-1}~Mpc^{-1}}$ (68% Confidence interval) with the catalogue method only. This result represents an improvement of $\sim 23~{{\%}}$ comparing the new 15 dark siren constrain with the previous 10 dark siren constraint and a reduction in uncertainty of $\sim 40~{{\%}}$ from the combination of 15 dark and bright sirens compared with the GW170817 bright siren alone. The combination of dark and bright siren GW170817 with recent jet constraints yields H0 of $68.0^{+4.4}_{-3.8}~{\rm km~s^{-1}~Mpc^{-1}}$, a $\sim 6~{{\%}}$ precision from Standard Sirens, reducing the previous constraint uncertainty by $\sim 10~{{\%}}$ .

A Possible Formation Scenario of the Gaia BH1: Inner Binary Merger in Triple Systems

Based on astrometric measurements and spectral analysis from Gaia DR3, two quiescent black hole (BH) binaries, Gaia BH1 and BH2, have been identified. Their origins remain controversial, particularly for Gaia BH1. By considering a rapidly rotating (ω/ω crit = 0.8) and strongly magnetized (B 0 = 5000 G) merger product, we find that, at typical Galactic metallicity, the merger product can undergo efficient chemically homogeneous evolution. This results in the merger product having a significantly smaller radius during its evolution compared to that of a normally evolving massive star. Under the condition that the initial triple stability is satisfied, we use the Multiple Stellar Evolution code and the MESA code to identify an initial hierarchical triple that can evolve into Gaia BH1. It initially consists of three stars with masses of 9.03 M ⊙, 3.12 M ⊙, and 1 M ⊙, with inner and outer orbital periods of 2.21 days and 121.92 days, and inner and outer eccentricities of 0.41 and 0.45, respectively. This triple initially experiences triple evolution dynamics instability (TEDI) followed by Roche lobe overflow (RLOF). During RLOF, the inner orbit shrinks, and tidal effects gradually suppress the TEDI. Eventually, the inner binary undergoes a merger through contact (or collision). Finally, using models of rapidly rotating and strongly magnetic stars, along with standard core-collapse supernova (SN) or failed supernova (FSN) models, we find that a postmerger binary (PMB) consisting of an 12.11 M ⊙ merger product and a 1 M ⊙ companion star (originally an outer tertiary) can avoid RLOF. After an SN or FSN with a low ejected mass of ∼0.22 M ⊙ and a low kick velocity ( 46−33+25kms−1 or 9−8+16kms−1 ), the PMB can form Gaia BH1 in the Galactic disk.

Binary neutron star merger offsets from their host galaxies. GW 170817 as a case study

The locations of binary neutron star (BNS) mergers within their host galaxies encode the systemic kicks that these systems received in the supernova aftermath. We investigate how the galactic potential and the systemic kicks shape the offset distribution of BNS mergers with a case study of GW 170817 and its host NGC 4993. We derived dynamical constraints on the host potential from integral field spectroscopy with Jeans anisotropic modelling. We evolved the trajectories of synthetic BNSs from the BPASS code in the galactic potential, using two different kick prescriptions to investigate how the observed offsets might differentiate between these two possibilities. The simulation was repeated after swapping the host potential with that of a dwarf galaxy, to test the effect of the potential on the offsets. The location of GW 170817 is entirely consistent with our predictions, regardless of large or small kicks, because the strong potential of NGC 4993 is only diagnostic of very large kicks. In galaxies of similar or greater mass, large offsets can constrain large kicks, while small offsets do not provide much information. In an old dwarf galaxy, on the other hand, small offsets can constrain small kicks, while large offsets would prevent host association.

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

As the sensitivity of ground-based gravitational wave detectors progressively increases, observations of black hole mergers will provide us with the joint distribution of their masses and spins. This will be a critical benchmark to validate different formation scenarios. Merging binary black holes formed through the evolution of isolated binary systems require both components to be stripped of their hydrogen envelopes before core-collapse. The rotation rates of such stripped stars are constrained by the critical rotation limit at their surface, including its deviation from the Keplerian value owing to the outward force provided by radiation. This sets a restriction on their angular momentum content at core-collapse. We aim to determine if this restriction plays a role in the spins of binary black hole mergers. We used detailed calculations of stripped stars with the MESA code at low metallicities ($Z=Z_ $Z_ and $Z_ to determine the dimensionless spins of black holes produced by critically rotating stellar progenitors. To study how such progenitors can arise, we considered their formation through chemically homogeneous evolution (CHE) in binary stars. We used a semi-analytical model to study the physical processes that determine the final angular momentum of CHE binaries, and compared our results against available population synthesis models that rely on detailed binary evolution calculations. We find that above black hole masses of $ 25M_ the dimensionless spin parameter of critically rotating stripped stars ($a=Jc/(GM^2)$) is below unity. This results in an exclusion region at high chirp masses and effective spins that cannot be populated by isolated binary evolution. CHE can produce binaries where both black holes hit this limit, producing a pileup at the boundary of the excluded region. High-spin black holes arise from very low-metallicity CHE systems with short delay times, which merge at higher redshifts. On the other hand, the contribution of CHE to merging binary black holes detected in the third observing run of the LVK collaboration is expected to be dominated by systems with low spins ($ eff &lt;0.5$) that merge near redshift zero. Owing to its higher projected sensitivity and runtime, the fourth observing run of the LVK collaboration can potentially place constraints on the high-spin population and the existence of a limit set by critical rotation.

Off-Axis Color Characteristics of Binary Neutron Star Merger Events: Applications for Space Multi-Band Variable Object Monitor and James Webb Space Telescope

With advancements in gravitational wave detection technology, an increasing number of binary neutron star (BNS) merger events are expected to be detected. Due to the narrow opening angle of jet cores, many BNS merger events occur off-axis, resulting in numerous gamma-ray bursts (GRBs) going undetected. Models suggest that kilonovae, which can be observed off-axis, offer more opportunities to be detected in the optical/near-infrared band as electromagnetic counterparts of BNS merger events. In this study, we calculate kilonova emission using a three-dimensional semi-analytical code and model the GRB afterglow emission with the open-source Python package afterglowpy at various inclination angles. Our results show that it is possible to identify the kilonova signal from the observed color evolution of BNS merger events. We also deduce the optimal observing window for SVOM/VT and JWST/NIRCam, which depends on the viewing angle, jet opening angle, and circumburst density. These parameters can be cross-checked with the multi-band afterglow fitting. We suggest that kilonovae are more likely to be identified at larger inclination angles, which can also help determine whether the observed signals without accompanying GRBs originate from BNS mergers.

Binary neutron star mergers using a discontinuous Galerkin-finite difference hybrid method

Abstract We present a discontinuous Galerkin-finite difference hybrid scheme that allows high-order shock capturing with the discontinuous Galerkin method for general relativistic magnetohydrodynamics in dynamical spacetimes. We present several optimizations and stability improvements to our algorithm that allow the hybrid method to successfully simulate single, rotating, and binary neutron stars. The hybrid method achieves the efficiency of discontinuous Galerkin methods throughout almost the entire spacetime during the inspiral phase, while being able to robustly capture shocks and resolve the stellar surfaces. We also use Cauchy-Characteristic evolution to compute the first gravitational waveforms at future null infinity from binary neutron star mergers. The simulations presented here are the first successful binary neutron star inspiral and merger simulations using discontinuous Galerkin methods.

Nonlinear studies of modifications to general relativity: Comparing different approaches

Studying the dynamical, nonlinear regime of modified theories of gravity remains a theoretical challenge that limits our ability to test general relativity. Here we consider two generally applicable, but approximate, methods for treating modifications to full general relativity that have been used to study binary black hole mergers and other phenomena in this regime, and compare solutions obtained by them to those from solving the full equations of motion. The first method evolves corrections to general relativity order by order in a perturbative expansion, while the second method introduces extra dynamical fields in such a way that strong hyperbolicity is recovered. We use shift-symmetric Einstein-scalar-Gauss-Bonnet gravity as a benchmark theory to illustrate the differences between these methods for several spacetimes of physical interest. We study the formation of scalar hair about initially nonspinning black holes, the collision of black holes with scalar charge, and the inspiral and merger of binary black holes. By directly comparing predictions, we assess the extent to which those from the approximate treatments can be meaningfully confronted with gravitational wave observations. We find that the order-by-order approach cannot faithfully track the solutions when the corrections to general relativity are non-negligible. The second approach, however, can provide consistent solutions, provided the timescale over which the dynamical fields are driven to their target values is made short compared to the physical timescales. Published by the American Physical Society 2024

Prospect of Detecting Magnetic Fields from Strong-magnetized Binary Neutron Stars

Binary neutron star mergers are unique sources of gravitational waves in multi-messenger astronomy. The inspiral phase of binary neutron stars can emit gravitational waves as chirp signals. The present waveform models of gravitational waves only considered the gravitational interaction. In this paper, we derive the waveform of the gravitational wave signal taking into account the presence of magnetic fields. We found that the electromagnetic interaction and radiation can introduce different frequency-dependent power laws for both the amplitude and frequency of the gravitational wave. We show from the results of the Fisher information matrix that the third-generation observation may detect magnetic dipole moments if the magnetic field is ∼1017 G.

Probing kaon meson condensations through gravitational waves during neutron star inspiral phases

Gravitational waves from binary neutron star merger events have shown unparalleled advantages in deciphering the stellar internal structure, yet it remains underexplored to seek clues about kaon meson condensations inside massive neutron stars through these events. To address this gap, this study examines the binary star inspiral processes with kaon meson condensations, employing various kaon meson potential wells to probe unique imprints in the inspiral gravitational waves. Our analysis reveals that the inspiral process of kaon meson condensation binary stars exhibits lower gravitational wave frequencies and shorter retarded times compared to traditional neutron stars. Furthermore, a deeper kaon meson potential well in these systems further shortens the inspiral duration time, resulting in distinctive gravitational waveforms and notable orbital phase deviations. These findings underscore the pivotal role of kaon meson condensations in modulating the inspiral gravitational waves, suggesting the inspiral gravitational waves could serve as a powerful tool for probing potential signals of kaon mesons within neutron stars.