ChatGPT scores a bad birdie in counting gravitational-wave chirps

How many gravitational-wave observations from compact object mergers have we seen to date? This seemingly simple question has a surprisingly complex answer that even ChatGPT struggles to answer. To shed light on this, we present a database with the literature's answers to this question. We find values spanning 67-100 for the number of detections from double compact object mergers to date, emphasizing that the exact number of detections is uncertain and depends on the chosen data analysis pipeline and underlying assumptions. We also review the number of gravitational-wave detections expected in the coming decades with future observing runs, finding values up to millions of detections per year in the era of Cosmic Explorer and Einstein Telescope. We present a publicly available code to visualize the detection numbers, highlighting the exponential growth in gravitational-wave observations in the coming decades and the exciting prospects of gravitational-wave astrophysics. See http://www.broekgaarden.nl/floor/wordpress/elementor-967/. We plan to keep this database up-to-date and welcome comments and suggestions for additional references.Comment: 1 April submission, with fun videos for visualizing the landscape of gravitational waves! (they are awesome!) See http://www.broekgaarden.nl/floor/wordpress/elementor-967

Streamlining and standardizing software citations with The Software Citation Station

Software is crucial for the advancement of astronomy especially in the context of rapidly growing datasets that increasingly require algorithm and pipeline development to process the data and produce results. However, software has not always been consistently cited, despite its importance to strengthen support for software development. To encourage, streamline, and standardize the process of citing software in academic work such as publications we introduce 'The Software Citation Station': a publicly available website and tool to quickly find or add software citationsComment: Visit the Software citation station at https://www.tomwagg.com/software-citation-station/ which includes a very easy form to add your own software package! comments welcome

Investigating the Cosmological Rate of Compact Object Mergers from Isolated Massive Binary Stars

Gravitational wave detectors are observing compact object mergers from increasingly far distances, revealing the redshift evolution of the binary black hole (BBH) -- and soon the black hole-neutron star (BHNS) and binary neutron star (BNS) -- merger rate. To help interpret these observations, we investigate the expected redshift evolution of the compact object merger rate from the isolated binary evolution channel. We present a publicly available catalog of compact object mergers and their accompanying cosmological merger rates from population synthesis simulations conducted with the COMPAS software. To explore the impact of uncertainties in stellar and binary evolution, our simulations use two-parameter grids of binary evolution models that vary the common-envelope efficiency with mass transfer accretion efficiency, and supernova remnant mass prescription with supernova natal kick velocity, respectively. We quantify the redshift evolution of our simulated merger rates using the local ($z\sim 0$) rate, the redshift at which the merger rate peaks, and the normalized differential rates (as a proxy for slope). We find that although the local rates span a range of $\sim 10^3$ across our model variations, their redshift-evolutions are remarkably similar for BBHs, BHNSs, and BNSs, with differentials typically within a factor $3$ and peaks of $z\approx 1.2-2.4$ across models. Furthermore, several trends in our simulated rates are correlated with the model parameters we explore. We conclude that future observations of the redshift evolution of the compact object merger rate can help constrain binary models for stellar evolution and gravitational-wave formation channels.Comment: 16 pages, 9 figures, link to code https://github.com/Adam-Boesky/Exploring_Parameter_Space, link to data https://gwlandscape.org.au/compas/publication/Q29tcGFzUHVibGljYXRpb25Ob2RlOjYz

STROOPWAFEL: Simulating rare outcomes from astrophysical populations, with application to gravitational-wave sources

Gravitational-wave observations of double compact object (DCO) mergers are providing new insights into the physics of massive stars and the evolution of binary systems. Making the most of expected near-future observations for understanding stellar physics will rely on comparisons with binary population synthesis models. However, the vast majority of simulated binaries never produce DCOs, which makes calculating such populations computationally inefficient. We present an importance sampling algorithm, STROOPWAFEL, that improves the computational efficiency of population studies of rare events, by focusing the simulation around regions of the initial parameter space found to produce outputs of interest. We implement the algorithm in the binary population synthesis code COMPAS, and compare the efficiency of our implementation to the standard method of Monte Carlo sampling from the birth probability distributions. STROOPWAFEL finds $\sim$25-200 times more DCO mergers than the standard sampling method with the same simulation size, and so speeds up simulations by up to two orders of magnitude. Finding more DCO mergers automatically maps the parameter space with far higher resolution than when using the traditional sampling. This increase in efficiency also leads to a decrease of a factor $\sim$3-10 in statistical sampling uncertainty for the predictions from the simulations. This is particularly notable for the distribution functions of observable quantities such as the black hole and neutron star chirp mass distribution, including in the tails of the distribution functions where predictions using standard sampling can be dominated by sampling noise.Comment: Accepted. Data and scripts to reproduce main results is publicly available. The code for the STROOPWAFEL algorithm will be made publicly available. Early inquiries can be addressed to the lead autho

The Binary Black Hole Merger Rate Deviates From the Cosmic Star Formation Rate: A Tug of War Between Metallicity and Delay Times

Gravitational-wave detectors are now making it possible to investigate how the merger rate of binary black holes (BBHs) evolves with redshift. In this study, we examine whether the BBH merger rate of isolated binaries deviates from a scaled star formation rate density (SFRD) -- a frequently used model in state-of-the-art research. To address this question, we conduct population synthesis simulations using COMPAS with a grid of stellar evolution models, calculate their cosmological merger rates, and compare them to a scaled SFRD. We find that our simulated rates deviate by factors up to $3.5\times$ at $z\sim0$ and $5\times$ at $z\sim 9$ due to two main phenomena: (i) The formation efficiency of BBHs is an order of magnitude higher at low metallicities than at solar metallicity; and (ii) BBHs experience a wide range of delays (from a few Myr to many Gyr) between formation and merger. Deviations are similar when comparing to a $\textit{delayed}$ SFRD, and even larger (up to $\sim 10\times$) when comparing to SFRD-based models scaled to the local merger rate. Interestingly, our simulations find that the BBH delay time distribution is redshift-dependent, increasing the complexity of the redshift distribution of mergers. We find similar results for simulated merger rates of BHNSs and BNSs. We conclude that the rate of BBH, BHNS, and BNS mergers from the isolated channel can significantly deviate from a scaled SFRD, and that future measurements of the merger rate will provide insights into the formation pathways of gravitational-wave sources.Comment: 20 pages, 12 figures, link to code https://github.com/Adam-Boesky/Exploring_Parameter_Space, link to data https://gwlandscape.org.au/compas/publication/Q29tcGFzUHVibGljYXRpb25Ob2RlOjYz

Stellar Black Holes and Compact Stellar Remnants

The recent observations of gravitational waves (GWs) by the LIGO-Virgo-KAGRA collaboration (LVK) have provided a new opportunity for studying our Universe. By detecting several merging events of black holes (BHs), LVK has spurred the astronomical community to improve theoretical models of single, binary, and multiple star evolution in order to better understand the formation of binary black hole (BBH) systems and interpret their observed properties. The final BBH system configuration before the merger depends on several processes, including those related to the evolution of the inner stellar structure and those due to the interaction with the companion and the environment (such as in stellar clusters). This chapter provides a summary of the formation scenarios of stellar BHs in single, binary, and multiple systems. We review all the important physical processes that affect the formation of BHs and discuss the methodologies used to detect these elusive objects and constrain their properties.Comment: To appear in Chapter 1 in the book Black Holes in the Era of Gravitational Wave Astronomy, ed. Arca Sedda, Bortolas, Spera, pub. Elsevier. All authors equally contributed to the chapter. Figures from other publications have been reproduced with permissio

Population synthesis of accreting white dwarfs: Rates and evolutionary pathways of H and He novae

Novae are some of the most commonly detected optical transients and have the potential to provide valuable information about binary evolution. Binary population synthesis codes have emerged as the most effective tool for modelling populations of binary systems, but such codes have traditionally employed greatly simplified nova physics, precluding detailed study. In this work, we implement a model treating H and He novae as individual events into the binary population synthesis code \binaryc. This treatment of novae represents a significant improvement on the `averaging' treatment currently employed in modern population synthesis codes. We discuss the evolutionary pathways leading to these phenomena and present nova event rates and distributions of several important physical parameters. Most novae are produced on massive white dwarfs, with approximately 70 and 55 per cent of nova events occurring on O/Ne white dwarfs for H and He novae respectively. Only 15 per cent of H-nova systems undergo a common-envelope phase, but these systems are responsible for the majority of H nova events. All He-accreting He-nova systems are considered post-common-envelope systems, and almost all will merge with their donor star in a gravitational-wave driven inspiral. We estimate the current annual rate of novae in M31 (Andromeda) to be approximately $41 \pm 4$ for H novae, underpredicting the current observational estimate of $65^{+15}_{-16}$, and $0.14\pm0.015$ for He novae. When varying common-envelope parameters, the H nova rate varies between 20 and 80 events per year.Comment: Accepted, MNRAS. 7 Jun 2020: Minor correction regarding AM CVn masses at period bounce, courtesy of P. Neuteufe

Common-Envelope Episodes that lead to Double Neutron Star formation

Close double neutron stars have been observed as Galactic radio pulsars, while their mergers have been detected as gamma-ray bursts and gravitational-wave sources. They are believed to have experienced at least one common-envelope episode during their evolution prior to double neutron star formation. In the last decades there have been numerous efforts to understand the details of the common-envelope phase, but its computational modelling remains challenging. We present and discuss the properties of the donor and the binary at the onset of the Roche-lobe overflow leading to these common-envelope episodes as predicted by rapid binary population synthesis models. These properties can be used as initial conditions for detailed simulations of the common-envelope phase. There are three distinctive populations, classified by the evolutionary stage of the donor at the moment of the onset of the Roche-lobe overflow: giant donors with fully-convective envelopes, cool donors with partially-convective envelopes, and hot donors with radiative envelopes. We also estimate that, for standard assumptions, tides would not circularise a large fraction of these systems by the onset of Roche-lobe overflow. This makes the study and understanding of eccentric mass-transferring systems relevant for double neutron star populations.Comment: 26 pages, 10 figures. Includes bug fix. Two new figures and an appendix adde

Multi-messenger prospects for black hole - neutron star mergers in the O4 and O5 runs

The existence of merging black hole-neutron star (BHNS) binaries has been ascertained through the observation of their gravitational wave (GW) signals. However, to date, no definitive electromagnetic (EM) emission has been confidently associated with these mergers. Such an association could help unravel crucial information on these systems, for example, their BH spin distribution, the equation of state (EoS) of NS and the rate of heavy element production. We model the multi-messenger (MM) emission from BHNS mergers detectable during the fourth (O4) and fifth (O5) observing runs of the LIGO-Virgo-KAGRA GW detector network, in order to provide detailed predictions that can help enhance the effectiveness of observational efforts and extract the highest possible scientific information from such remarkable events. Our methodology is based on a population synthesis-approach, which includes the modelling of the signal-to-noise ratio of the GW signal in the detectors, the GW-inferred sky localization of the source, the kilonova (KN) optical and near-infrared light curves, the relativistic jet gamma-ray burst (GRB) prompt emission peak photon flux, and the GRB afterglow light curves in the radio, optical and X-ray bands. The resulting prospects for BHNS MM detections during O4 are not promising, with a GW detection rate of $15.0^{+15.4}_{-8.8}$ yr$^{-1}$, but joint MM rates of $\sim 10^{-1}$ yr$^{-1}$ for the KN and $\sim 10^{-2}$ yr$^{-1}$ for the jet-related emission. In O5 we find an overall increase in expected detection rates by around an order of magnitude, owing to both the enhanced sensitivity of the GW detector network, and the coming online of future EM facilities. Finally, we discuss direct searches for the GRB radio afterglow with large-field-of-view instruments as a new possible follow-up strategy in the context of ever-dimming prospects for KN detection.Comment: Submitted to A&A. 17 pages, 11 figures, 2 tables. Comments are welcome

Characterizing Gravitational Wave Detector Networks: From A♯^\sharp to Cosmic Explorer

Gravitational-wave observations by the Laser Interferometer Gravitational-Wave Observatory (LIGO) and Virgo have provided us a new tool to explore the universe on all scales from nuclear physics to the cosmos and have the massive potential to further impact fundamental physics, astrophysics, and cosmology for decades to come. In this paper we have studied the science capabilities of a network of LIGO detectors when they reach their best possible sensitivity, called A#, and a new generation of observatories that are factor of 10 to 100 times more sensitive (depending on the frequency), in particular a pair of L-shaped Cosmic Explorer observatories (one 40 km and one 20 km arm length) in the US and the triangular Einstein Telescope with 10 km arms in Europe. We use a set of science metrics derived from the top priorities of several funding agencies to characterize the science capabilities of different networks. The presence of one or two A# observatories in a network containing two or one next generation observatories, respectively, will provide good localization capabilities for facilitating multimessenger astronomy and precision measurement of the Hubble parameter. A network of two Cosmic Explorer observatories and the Einstein Telescope is critical for accomplishing all the identified science metrics including the nuclear equation of state, cosmological parameters, growth of black holes through cosmic history, and make new discoveries such as the presence of dark matter within or around neutron stars and black holes, continuous gravitational waves from rotating neutron stars, transient signals from supernovae, and the production of stellar-mass black holes in the early universe. For most metrics the triple network of next generation terrestrial observatories are a factor 100 better than what can be accomplished by a network of three A# observatories.Comment: 45 pages, 20 figure