Merger Time Research Articles

Merger-induced galaxy transformations in the artemis simulations

ABSTRACT Using the artemis set of 45 high-resolution cosmological simulations, we investigate a range of merger-induced dynamical transformations of Milky Way-like galaxies. We first identify populations of accreted stars on highly radial orbits, similar to the ‘Gaia Sausage’ in the Milky Way. We show that ≈1/3 of the artemis galaxies contain a similar feature, and confirm that they usually comprise stellar debris from the most massive accreted satellite. Selecting 15 galaxies with discs at the present-day, we study their changes around the times of the GS-like mergers. Dark matter haloes of many of these exhibit global changes in shape and orientation, with almost half becoming significantly more spherical when the mergers occur. Focusing on the galaxies themselves, we find that 4/15 have stellar discs which experience large changes in the orientation of their angular momentum (AM) axes, at rates of up to ∼60 degrees Gyr−1. By calculating the orbital angular momentum axes of the satellites before they are accreted, we show that there is a tendency for the disc’s AM to become more aligned with this axis after the merger. We also investigate the origin of in situ retrograde stars, analogous to the ‘Splash’ in the Milky Way. Tracing them back to earlier snapshots, we demonstrate that they were often disrupted on to their extreme orbits by multiple early mergers. We also find that the total mass of these stars outside the central regions positively correlates with the total accreted stellar mass.

Emergence of galactic morphologies at cosmic dawn: input from numerical modelling

ABSTRACTWe employ a series of high-resolution zoom-in cosmological simulations to analyse the emerging morphology of main galaxies in dark matter haloes at z ≳ 2. We choose haloes of similar masses, ${\rm log}\, M_{\rm vir}/{\rm M_\odot }\sim 11.65\pm 0.05$, at the target zf = 6, 4, and 2. The rationale for this choice allows us to analyse how the different growth rate in these haloes propagates down to galaxy scales, affecting their basic parameters. Halos were embedded in high/low overdensity regions, and two versions of a galactic wind feedback were employed. Our main results are: (1) Although our galaxies evolve in different epochs, their global parameters remain within narrow range. Their morphology, kinematics, and stellar populations differ substantially, yet all host sub-kpc stellar bars; (2) The star formation rates appear higher for larger zf; (3) Bulges and stellar spheroids were separated by stellar kinematics, discy bulges were revealed using the Sersic method and photometry.The bulge-to-total mass ratios appear independent of the last merger time for all zf. The spheroid-to-total mas ratios lie within ∼0.5–0.8; (4) The synthetic redshifted, pixelized, and PSF-degraded JWST images allow detection of stellar discs at all zf. (5) Based on the kinematic decomposition, rotational support in discs depends on the feedback type, but increases with decreasing zf; (6) Finally, the ALMA images detect discs at all zf, but spiral structure is detectable in zf = 2 galaxies. Moreover, galaxies follow the Tully–Fisher relation, being separated only by the galactic wind feedback.

Stellar mergers as the origin of the blue main-sequence band in young star clusters

Recent high-quality Hubble Space Telescope (HST) photometry shows that the main sequences (MS) stars of young star clusters form two discrete components in the color-magnitude diagram (CMD). Based on their distribution in the CMD, we show that stars of the blue MS component can be understood as slow rotators originating from stellar mergers. We derive the masses of the blue MS stars, and find that they follow a nearly flat mass function, which supports their unusual formation path. Our results imply that the cluster stars gain their mass in two different ways, by disk accretion leading to rapid rotation, contributing to the red MS, or by binary merger leading to slow rotation and populating the blue MS. We also derive the approximate merger time of the individual stars of the blue MS component, and find a strong early peak in the merger rate, with a lower level merger activity prevailing for tens of Myr. This supports recent binary formation models, and explains new velocity dispersion measurements for members of young star clusters. Our findings shed new light on the origin of the bi-modal mass, spin, and magnetic field distributions of main-sequence stars.

An accurate comprehensive approach to substructure: III. Masses and formation times of the host haloes

ABSTRACT With this paper, we complete a comprehensive study of substructure in dark matter haloes. In Paper I, we derived the radial distribution and mass function (MF) of accreted subhaloes (scaled to the radius and mass of the host halo) and showed that they are essentially universal. This is not the case, however, for those of stripped subhaloes, which depend on halo mass and assembly history. In Paper II, we derived these latter properties in the simplest case of purely accreting haloes. Here, we extend the study to ordinary haloes having suffered major mergers. After showing that all the properties of substructure are encoded in the mean truncated-to-original subhalo mass ratio profile, we demonstrate that the dependence of the subhalo MF on halo mass arises from their mass-dependent concentration, while the shape of the subhalo radial distribution depends on the time of the last major merger of the host halo. In this sense, the latter property is a better probe of halo formation time than the former. Unfortunately, this is not the case for the radial distribution of satellites as this profile is essentially disconnected from subhalo stripping and the properties of accreted subhaloes are independent of the halo formation time.

The Dependence of the Fraction of Radio Luminous Quasars on Redshift and its Theoretical Implications

Abstract While radio emission in quasars can be contributed to by a variety of processes (involving star-forming regions, accretion disk coronas and winds, and jets), the powering of the radio loudest quasars must involve very strong jets, presumably launched by the Blandford–Znajek mechanism incorporating the magnetically arrested disk (MAD) scenario. We focus on the latter and investigate the dependence of their fraction on redshift. We also examine the dependence of the radio-loud fraction (RLF) on BH mass (M BH) and Eddington ratio (λ Edd), while excluding the redshift bias by narrowing its range. In both of these investigations, we remove the bias associated with: (1) the diversity of source selection by constructing two well-defined, homogeneous samples of quasars (first within 0.7 ≤ z ≤ 1.9, second within 0.5 ≤ z ≤ 0.7); and (2) a strong drop in the RLF of quasars at smaller BH masses by choosing those with BH masses larger than 108.5 M ⊙. We confirm some of the previous results showing the increase in the fraction of radio-loud quasars with cosmic time and that this trend can be even steeper if we account for the bias introduced by the dependence of the RLF on BH mass, whereas the bias introduced by the dependence of the RLF on Eddington ratio is shown to be negligible. Assuming that quasar activities are triggered by galaxy mergers, we argue that such an increase can result from the slower drop with cosmic time of mixed mergers than of wet mergers.

Binary Black Hole Formation with Detailed Modeling: Stable Mass Transfer Leads to Lower Merger Rates

Abstract Rapid binary population synthesis codes are often used to investigate the evolution of compact-object binaries. They typically rely on analytical fits of single-star evolutionary tracks and parameterized models for interactive phases of evolution (e.g., mass transfer on a thermal timescale, determination of dynamical instability, and common envelope) that are crucial to predict the fate of binaries. These processes can be more carefully implemented in stellar structure and evolution codes such as MESA. To assess the impact of such improvements, we compare binary black hole mergers as predicted in models with the rapid binary population synthesis code COSMIC to models ran with MESA simulations through mass transfer and common-envelope treatment. We find that results significantly differ in terms of formation paths, the orbital periods and mass ratios of merging binary black holes, and consequently merger rates. While common-envelope evolution is the dominant formation channel in COSMIC, stable mass transfer dominates in our MESA models. Depending upon the black hole donor mass, and mass-transfer and common-envelope physics, at subsolar metallicity, COSMIC overproduces the number of binary black hole mergers by factors of 2–35 with a significant fraction of them having merger times orders of magnitude shorter than the binary black holes formed when using detailed MESA models. Therefore we find that some binary black hole merger rate predictions from rapid population syntheses of isolated binaries may be overestimated by factors of ∼ 5–500. We conclude that the interpretation of gravitational-wave observations requires the use of detailed treatment of these interactive binary phases.

Semianalytical approach for sky localization of gravitational waves

Rapid sky localization of gravitational wave sources is crucial to enable prompt electromagnetic follow-ups. In this article, we present a novel semianalytical approach for sky localization of gravitational waves from compact binary coalescences. We use the Bayesian framework with an analytical approximation to the prior distributions for a given astrophysical model. We derive a semianalytical solution to the posterior distribution of source directions. This method only requires one-fold numerical integral that marginalizes over the merger time, compared to the five-fold numerical integration otherwise needed in the Bayesian localization method. The performance of the method is demonstrated using a set of binary neutron stars (BNS) injections on Gaussian noise using LIGO-Virgo's design and O2 sensitivity. We find the median of 90% confidence area in O2 sensitivity to be $\mathcal{O}(10^2) ~\mathrm{deg}^2$, comparable to that of the existing LIGO-Virgo online localization method Bayestar and parameter estimation toolkit LALInference. In the end, we apply this method to localize the BNS event GW170817 and find the 50% (90%) confidence region of 11 $\mathrm{deg}^2$ (50 $\mathrm{deg}^2$). The detected optical counterpart of GW170817 resides within our 50% confidence area.

Parameter estimation bias from overlapping binary black hole events in second generation interferometers

Since the initial detection of Gravitational Waves in 2015, 50 candidate events have been reported by the LIGO-Virgo-KAGRA collaboration. As the current generation of detectors move towards their design sensitivity the rate of these detections will increase. The next generation of detectors are likely to have high enough sensitivities that multiple merging binaries will be visible at the same time. In this paper we show that this is likely to happen before the end of the decade, with the move to the LIGO-Voyager detector. We investigate the situation of overlapping Binary-Black-Hole mergers in these detectors. We find that current parameter estimation techniques are capable of distinguishing the louder of two merging BBH events, without significant bias, when their merger times are not less than $\sim0.1$ seconds apart and when the ratio of the signal-to-noise ratios of the systems is uneven. This region of overlapping parameter space is dependent upon the sky locations of the signals and the relation of those locations to the light travel time between detectors. We also find that, if two signals are highly overlapping, then the recovered set of parameters often show strong evidence of precession. Finally we show that bias can occur even when the signal causing the bias is below the detection threshold.

Extremely efficient mergers of intermediate-mass black hole binaries in nucleated dwarf galaxies

ABSTRACT Gravitational waves emitted by merging black holes between ∼104 and 107 M⊙ will be detectable by the Laser Interferometer Space Antenna (LISA) with signal-to-noise ratios of several hundred out to redshift 20. Supermassive black hole (107–1010 M⊙) binary formation, coalescence, and merger within massive galaxies is well-studied. However, intermediate-mass black holes (IMBHs) between ∼104 and 106 M⊙ are hosted by low mass and dwarf galaxies; it is not trivial to extrapolate black hole merger time-scales to this IMBH binary regime, due to the starkly different host galaxy structure, kinematics, and morphology compared to massive galaxy hosts. We perform ultrahigh resolution N-body simulations to study IMBH dynamics in nucleated dwarf galaxies whose structural parameters are obtained from observations of nearby dwarf galaxies. Starting from 50 pc, an IMBH quickly forms a binary. Thereafter, the binary orbit shrinks rapidly due to the high central stellar densities furnished by nuclear star clusters (NSCs). We find high eccentricities (e ∼ 0.4–0.99) in our suite of IMBH binaries, and residual eccentricity may persist to the LISA regime. IMBH merger times are typically a few hundred million years, with a few exceptionally short merger times for high eccentricities. We find that IMBH-stellar encounters originate pre-dominantly from NSCs, if the NSC-to-IMBH binary mass ratio is greater than 10; otherwise, bulge stars contribute significantly. As the IMBH binary ejects stars, however, the NSCs are disrupted. We conclude that comparable-mass IMBHs merge very efficiently in nucleated dwarf galaxies, making them promising LISA sources, as well as a channel for IMBH growth.

Biases in parameter estimation from overlapping gravitational-wave signals in the third-generation detector era

In the past few years, the detection of gravitational waves from compact binary coalescences with the Advanced LIGO and Advanced Virgo detectors has become routine. Future observatories will detect even larger numbers of gravitational-wave signals, which will also spend a longer time in the detectors' sensitive band. This will eventually lead to overlapping signals, especially in the case of Einstein Telescope (ET) and Cosmic Explorer (CE). Using realistic distributions for the merger rate as a function of redshift as well as for component masses in binary neutron star and binary black hole coalescences, we map out how often signal overlaps of various types will occur in an ET-CE network over the course of a year. We find that a binary neutron star signal will typically have tens of overlapping binary black hole and binary neutron star signals. Moreover, it will happen up to tens of thousands of times per year that two signals will have their end times within seconds of each other. In order to understand to what extent this would lead to measurement biases with current parameter estimation methodology, we perform injection studies with overlapping signals from binary black hole and/or binary neutron star coalescences. Varying the signal-to-noise ratios, the durations of overlap, and the kinds of overlapping signals, we find that in most scenarios the intrinsic parameters can be recovered with negligible bias. However, we find large offsets for a short binary black hole or a quieter binary neutron star signal overlapping with a long and louder binary neutron star event when the merger times are sufficiently close. Although based on a limited number of simulations, our studies may be an indicator of where improvements are required to ensure reliable estimation of source parameters for all detected compact binary signals as we go from second-generation to third-generation detectors.

The impact of turbulent mixing on the galactic r-process enrichment by binary neutron star mergers

ABSTRACT We study the enrichment of the interstellar medium (ISM) with rapid neutron capture (r-process) elements produced in binary neutron star (BNS) mergers. We use a semi-analytic model to describe galactic evolution, with merger rates and time delay distributions of BNS mergers consistent with the latest population synthesis models. In order to study the dispersion of the relative abundances of r-process elements and iron, we applied a turbulent mixing scheme, where the freshly synthesized elements are gradually dispersed in the ISM. We show that within our model the abundances observed in Milky Way stars, in particular the scatter at low metallicities, can be entirely explained by BNS mergers. Our results suggest that BNS mergers could be the dominant source of r-process elements in the Galaxy.

Progenitors of low-mass binary black-hole mergers in the isolated binary evolution scenario

Context. The formation history, progenitor properties, and expected rates of the binary black holes discovered by the LIGO-Virgo collaboration via the gravitational-wave emission during their coalescence are a topic of active research. Aims. We aim to study the progenitor properties and expected rates of the two lowest-mass binary black hole mergers, GW151226 and GW170608, detected within the first two Advanced LIGO-Virgo observing runs, in the context of the classical isolated binary-evolution scenario. Methods. We used the publicly available 1D-hydrodynamic stellar-evolution code MESA, which we adapted to include the black-hole formation and the unstable mass transfer developed during the so-called common-envelope phase. Using more than 60 000 binary simulations, we explored a wide parameter space for initial stellar masses, separations, metallicities, and mass-transfer efficiencies. We obtained the expected distributions for the chirp mass, mass ratio, and merger time delay by accounting for the initial stellar binary distributions. We predicted the expected merger rates and compared them with those of the detected gravitational-wave events. We studied the dependence of our predictions with respect to the (as yet) unconstrained parameters inherent to binary stellar evolution. Results. Our simulations for both events show that while the progenitors we obtained are compatible over the entire range of explored metallicities, they show a strong dependence on the initial masses of the stars, according to stellar winds. All the progenitors we found follow a similar evolutionary path, starting from binaries with initial separations in the 30−200 R⊙ range experiencing a stable mass transfer interaction before the formation of the first black hole, followed by a second unstable mass-transfer episode leading to a common-envelope ejection that occurs either when the secondary star crosses the Hertzsprung gap or when it is burning He in its core. The common-envelope phase plays a fundamental role in the considered low-mass range: only progenitors experiencing such an unstable mass-transfer phase are able to merge in less than a Hubble time. Conclusions. We find integrated merger-rate densities in the range 0.2–5.0 yr−1 Gpc−3 in the Local Universe for the highest mass-transfer efficiencies explored here. The highest rate densities lead to detection rates of 1.2–3.3 yr−1, which are compatible with the observed rates. The common-envelope efficiency αCE has a strong impact on the progenitor populations. A high-efficiency scenario with αCE = 2.0 is favoured when comparing the expected rates with the observations.

Are we observing an NSC in course of formation in the NGC 4654 galaxy?

ABSTRACT We use direct N-body simulations to explore some possible scenarios for the future evolution of two massive clusters observed towards the centre of NGC 4654, a spiral galaxy with mass similar to that of the Milky Way. Using archival HST data, we obtain the photometric masses of the two clusters, M = 3 × 105 M⊙ and M = 1.7 × 106 M⊙, their half-light radii, Reff ∼ 4 pc and Reff ∼ 6 pc, and their projected distances from the photometric centre of the galaxy (both <22 pc). The knowledge of the structure and separation of these two clusters (∼24 pc) provides a unique view for studying the dynamics of a galactic central zone hosting massive clusters. Varying some of the unknown cluster orbital parameters, we carry out several N-body simulations showing that the future evolution of these clusters will inevitably result in their merger. We find that, mainly depending on the shape of their relative orbit, they will merge into the galactic centre in less than 82 Myr. In addition to the tidal interaction, a proper consideration of the dynamical friction braking would shorten the merging times up to few Myr. We also investigate the possibility to form a massive nuclear star cluster (NSC) in the centre of the galaxy by this process. Our analysis suggests that for low-eccentricity orbits, and relatively long merger times, the final merged cluster is spherical in shape, with an effective radius of few parsecs and a mass within the effective radius of the order of $10^5\, \mathrm{M_{\odot }}$. Because the central density of such a cluster is higher than that of the host galaxy, it is likely that this merger remnant could be the likely embryo of a future NSC.

Massive black hole binary inspiral and spin evolution in a cosmological framework

ABSTRACT Massive black hole (MBH) binary inspiral time-scales are uncertain, and their spins are even more poorly constrained. Spin misalignment introduces asymmetry in the gravitational radiation, which imparts a recoil kick to the merged MBH. Understanding how MBH binary spins evolve is crucial for determining their recoil velocities, their gravitational wave (GW) waveforms detectable with Laser Interferometer Space Antenna, and their retention rate in galaxies. Here, we introduce a sub-resolution model for gas- and gravitational wave (GW)-driven MBH binary spin evolution using accreting MBHs from the Illustris cosmological hydrodynamic simulations. We also model binary inspiral via dynamical friction, stellar scattering, viscous gas drag, and GW emission. Our model assumes that the circumbinary disc always removes angular momentum from the binary. It also assumes differential accretion, which causes greater alignment of the secondary MBH spin in unequal-mass mergers. We find that 47 per cent of the MBHs in our population merge by z = 0. Of these, 19 per cent have misaligned primaries and 10 per cent have misaligned secondaries at the time of merger in our fiducial model with initial eccentricity of 0.6 and accretion rates from Illustris. The MBH misalignment fraction depends strongly on the accretion disc parameters, however. Reducing accretion rates by a factor of 100, in a thicker disc, yields 79 and 42 per cent misalignment for primaries and secondaries, respectively. Even in the more conservative fiducial model, more than 12 per cent of binaries experience recoils of >500 km s−1, which could displace them at least temporarily from galactic nuclei. We additionally find that a significant number of systems experience strong precession.

The Stellar Merger Scenario for Black Holes in the Pair-instability Gap

Abstract The recent detection of GW190521 stimulated ideas on how to populate the predicted black hole (BH) pair-instability (PI) mass gap. One proposal is the dynamical merger of two stars below the PI regime forming a star with a small core and an oversized envelope. We outline the main challenges this scenario faces to form one BH in the gap. In particular, the core needs to avoid growing during the merger, and the merger product needs to retain enough mass, including in the subsequent evolution, and at core collapse (CC). We explore this scenario with detailed stellar evolution calculations, starting with ad hoc initial conditions enforcing no core growth during the merger. We find that these massive merger products are likely to be helium-rich and spend most of their remaining lifetime within regions of instabilities in the Herzsprung–Russell diagram, such as luminous blue variable eruptions. An energetic estimate of the amount of mass loss neglecting the back reaction of the star suggests that the total amount of mass that can be removed at low metallicity is ≲1 M ⊙. This is small enough that at CC our models are retaining sufficient mass to form BHs in the PI gap similar to the recent ones detected by LIGO/Virgo. However, mass loss at the time of merger, the resulting core structure, and the mass loss at CC still need to be quantified for these models to confirm the viability of this scenario.

Gravitational and electromagnetic radiation from binary black holes with electric and magnetic charges: Circular orbits on a cone

We derive the equations of motion of dyonic black hole binaries with electric and magnetic charges and explore features of static orbits. By using a Newtonian method with the inclusion of radiation reaction, we calculate the total emission rate of energy and angular momentum due to gravitational radiation and electromagnetic radiation for circular orbits. Moreover, we obtain the evolutions of orbits and calculate merger times of dyonic binaries. We show that electric and magnetic charges significantly suppress the merger times of dyonic binaries. Our results provide rich information about dyonic binaries and can be used to test black holes with magnetic charges.

The Milky Way’s Shell Structure Reveals the Time of a Radial Collision

Abstract We identify shell structures in the Milky Way for the first time. We find two shells in the Virgo Overdensity region and two shells in the Hercules Aquila Cloud region using Sloan Digital Sky Survey, Gaia, and LAMOST data. These shell stars are a subset of the substructure previously identified as the Virgo Radial Merger (VRM). Timing arguments for these shells indicate that their progenitor dwarf galaxy passed through the Galactic center 2.7 ± 0.2 Gyr ago. Based on the time of collision, it is also possible that the VRM is related to the phenomenon that created phase-space spirals in the vertical motion of the disk and/or the Splash and could have caused a burst of star formation in the inner disk. We analyze phase mixing in a collection of radial merger N-body simulations and find that shell structure similar to that observed in Milky Way data disappears by 5 Gyr after collision with the Galactic center. The method used to calculate the merger time of the VRM was able to reliably recover the correct merger times for these simulations. Previous work supports the idea that the VRM and the Gaia Sausage/Gaia–Enceladus Merger are the same. However, the Gaia Sausage is widely believed to be 8–11 Gyr old. The disparate ages could be reconciled if the larger age is associated with an infall time when the progenitor crossed the virial radius; we do not constrain the time at which the progenitor became bound to the Milky Way. Alternatively, the Gaia Sausage could be younger than previously thought.

Chirp mass and spin of binary black holes from first star remnants

ABSTRACT We performed Population III (Pop III) binary evolution using population synthesis simulations for seven different models. We found that Pop III binaries tend to be binary black holes (BBHs) with chirp mass Mchirp ∼ 30 M⊙ and they can merge in the present day, due to a long merger time. The merger rate densities of Pop III BBHs at z = 0 are in the range 3.34–21.2 $\rm yr^{-1}\,Gpc^{-3}$ which is consistent with the Advanced Laser Interferometer Gravitational Wave Observatory (aLIGO)/Advanced Virgo (aVIRGO) result of 9.7–101 $\rm yr^{-1}\,Gpc^{-3}$. These Pop III binaries might contribute some portion of the massive BBH gravitational wave (GW) sources detected by aLIGO/aVIRGO. We also calculated the redshift dependence of Pop III BBH mergers. We found that Pop III low-spin BBHs tend to merge at low redshift, while Pop III high-spin BBHs merge at high redshift, which can be confirmed by future GW detectors such as Einstein Telescope (ET), Cosmic Explorer (CE), and DECi-hertz Interferometer Gravitational wave Observatory (DECIGO). These detectors can also check the redshift dependence of the BBH merger rate and spin distribution. Our results show that, except for one model, the mean effective spin 〈χeff〉 at z = 0 lies in the range 0.02–0.3, while at z = 10 it is 0.16–0.64. Therefore, massive stellar-mass BBH detection by GWs will be key for stellar evolution study in the early Universe.

Joint galaxy–galaxy lensing and clustering constraints on galaxy formation

ABSTRACT We compare predictions for galaxy–galaxy lensing profiles and clustering from the Henriques et al. public version of the Munich semi-analytical model (SAM) of galaxy formation and the IllustrisTNG suite, primarily TNG300, with observations from KiDS + GAMA and SDSS-DR7 using four different selection functions for the lenses (stellar mass, stellar mass and group membership, stellar mass and isolation criteria, and stellar mass and colour). We find that this version of the SAM does not agree well with the current data for stellar mass-only lenses with $M_\ast \gt 10^{11}\, \mathrm{ M}_\odot$. By decreasing the merger time for satellite galaxies as well as reducing the radio-mode active galactic nucleus accretion efficiency in the SAM, we obtain better agreement, both for the lensing and the clustering, at the high-mass end. We show that the new model is consistent with the signals for central galaxies presented in Velliscig et al. Turning to the hydrodynamical simulation, TNG300 produces good lensing predictions, both for stellar mass-only (χ2 = 1.81 compared to χ2 = 7.79 for the SAM) and locally brightest galaxy samples (χ2 = 3.80 compared to χ2 = 5.01). With added dust corrections to the colours it matches the SDSS clustering signal well for red low-mass galaxies. We find that both the SAMs and TNG300 predict $\sim 50\, {{\ \rm per\ cent}}$ excessive lensing signals for intermediate-mass red galaxies with 10.2 < log10M*[M⊙] < 11.2 at $r \approx 0.6\, h^{-1}\, \text{Mpc}$, which require further theoretical development.

Mass Ratios of Merging Double Neutron Stars as Implied by the Milky Way Population

Abstract Of the seven known double neutron stars (DNSs) with precisely measure masses in the Milky Way that will merge within a Hubble time, all but one has a mass ratio, q, close to unity. Recently, precise measurements of three post-Keplerian parameters in the DNS J1913+1102 constrain this system to have a significantly non-unity mass ratio of 0.78 ± 0.03. One may be tempted to conclude that approximately one out of seven (14%) DNS mergers detected by gravitational-wave observatories will have mass ratios significantly different from unity. However J1913+1102 has a relatively long lifetime (pulsar characteristic age plus the system’s merger time due to gravitational-wave radiation) of ≈3 Gyr. We show that when system lifetimes and observational biases are taken into account, the population of Galactic DNSs implies that ≃98% of all merging DNSs will have q > 0.9. We then apply two separate fitting formulas informed by 3D hydrodynamic simulations of DNS mergers to our results on Galactic DNS masses, finding that either ≃0.004 or ≃0.009 of material will be ejected at merger, depending on which formula is used. These ejecta masses have implications for both the peak bolometric luminosities of electromagnetic counterparts (which we find to be ∼1041 erg s−1) as well as the r-process enrichment of the Milky Way.