Search

GW230529 is the first compact binary coalescence detected by the LIGO-Virgo-KAGRA collaboration with at least one component mass confidently in the lower mass-gap, corresponding to the range 3-5$M_{\odot}$. If interpreted as a neutron star-black hole merger, this event has the most symmetric mass ratio detected so far and therefore has a relatively high probability of producing electromagnetic (EM) emission. However, no EM counterpart has been reported. At the merger time $t_0$, Swift-BAT and Fermi-GBM together covered 100$\%$ of the sky. Performing a targeted search in a time window $[t_0-20 \text{s},t_0+20 \text{s}]$, we report no detection by the Swift-BAT and the Fermi-GBM instruments. Combining the position-dependent $\gamma-$ray flux upper limits and the gravitational-wave posterior distribution of luminosity distance, sky localization and inclination angle of the binary, we derive constraints on the characteristic luminosity and structure of the jet possibly launched during the merger. Assuming a top-hat jet structure, we exclude at 90$\%$ credibility the presence of a jet which has at the same time an on-axis isotropic luminosity $\gtrsim 10^{48}$ erg s$^{-1}$, in the bolometric band 1 keV-10 MeV, and a jet opening angle $\gtrsim 15$ deg. Similar constraints are derived testing other assumptions about the jet structure profile. Excluding GRB 170817A, the luminosity upper limits derived here are below the luminosity of any GRB observed so far., Comment: 18 pages, 1 table, 11 figures

This study investigates the origins of GW230529, delving into its formation from massive stars within isolated binary systems. Utilizing population synthesis models, we present compelling evidence that the neutron star component forms second. However, the event's low signal-to-noise ratio introduces complexities in identifying the underlying physical mechanisms driving its formation. Augmenting our analysis with insights from numerical relativity, we estimate the final black hole mass and spin to be approximately $5.3 M_\odot$ and $0.53$, respectively. Furthermore, we employ the obtained posterior samples to calculate the ejecta mass and kilonova light curves resulting from r-process nucleosynthesis. We find the ejecta mass to range within $0-0.06 M_{\odot}$, contingent on the neutron star equation of state. The peak brightness of the kilonovae light curves indicates that targeted follow-up observations with a Rubin-like observatory may have detected this emission., Comment: 12 figures, 1 table, 15 pages

General relativity (GR) has proven to be a highly successful theory of gravity since its inception. The theory has thrivingly passed numerous experimental tests, predominantly in weak gravity, low relative speeds, and linear regimes, but also in the strong-field and very low-speed regimes with binary pulsars. Observable gravitational waves (GWs) originate from regions of spacetime where gravity is extremely strong, making them a unique tool for testing GR, in previously inaccessible regions of large curvature, relativistic speeds, and strong gravity. Since their first detection, GWs have been extensively used to test GR, but no deviations have been found so far. Given GR's tremendous success in explaining current astronomical observations and laboratory experiments, accepting any deviation from it requires a very high level of statistical confidence and consistency of the deviation across GW sources. In this paper, we compile a comprehensive list of potential causes that can lead to a false identification of a GR violation in standard tests of GR on data from current and future ground-based GW detectors. These causes include detector noise, signal overlaps, gaps in the data, detector calibration, source model inaccuracy, missing physics in the source and in the underlying environment model, source misidentification, and mismodeling of the astrophysical population. We also provide a rough estimate of when each of these causes will become important for tests of GR for different detector sensitivities. We argue that each of these causes should be thoroughly investigated, quantified, and ruled out before claiming a GR violation in GW observations., Comment: Review article; 1 figure; 1 table; comments welcome

We study the ringdown signal of black holes formed in prompt-collapse binary neutron star mergers. We analyze data from $47$ numerical relativity simulations. We show that the $(\ell=2,m=2)$ and $(\ell=2,m=1)$ multipoles of the gravitational wave signal are well fitted by decaying damped exponentials, as predicted by black-hole perturbation theory. We show that the ratio of the amplitude in the two modes depends on the progenitor binary mass ratio $q$ and reduced tidal parameter $\tilde\Lambda$. Unfortunately, the numerical uncertainty in our data is too large to fully quantify this dependency. If confirmed, these results will enable novel tests of general relativity in the presence of matter with next-generation gravitational-wave observatories., Comment: 15 pages, 7 figures. Accepted version

We study the detectability of postmerger QCD phase transitions in neutron star binaries with next-generation gravitational-wave detectors Cosmic Explorer and Einstein Telescope. We perform numerical relativity simulations of neutron star mergers with equations of state that include a quark deconfinement phase transition through either a Gibbs or Maxwell construction. These are followed by Bayesian parameter estimation of the associated gravitational-wave signals using the $\tt{NRPMw}$ waveform model, with priors inferred from the analysis of the inspiral signal. We assess the ability of the model to measure the postmerger peak frequency $f_2^{\rm peak}$ and identify aspects that should be improved in the model. We show that, even at postmerger signal to noise ratios as low as 10, the model can distinguish (at the 90% level) $f_2^{\rm peak}$ between binaries with and without a phase transition in most cases. Phase-transition induced deviations in the $f_2^{\rm peak}$ from the predictions of equation-of-state insensitive relations can also be detected if they exceed $1.6\,\sigma$. Our results suggest that next-generation gravitational wave detectors can measure phase transition effects in binary neutron star mergers. However, unless the phase transition is ``strong'', disentangling it from other hadronic physics uncertainties will require significant theory improvements.

We present, for the first time, recoil velocity estimates for binary neutron star mergers using data from numerical relativity simulations. We find that binary neutron star merger remnants can have recoil velocity of the order of a few tens of km/s and as high as $150$ km/s in our dataset. These recoils are attained due to equivalent contributions from the anisotropic gravitational wave emission as well as the asymmetric ejection of dynamical matter during the merger. We provide fits for net recoil velocity as well as its ejecta component as a function of the amount of ejected matter, which may be useful when constraints on the ejected matter are obtained through electromagnetic observations. We also estimate the mass and spin of the remnants and find them to be in the range $[2.34, 3.38] M_{\odot}$ and $[0.63, 0.82]$ respectively, for our dataset., Comment: 14 pages, 9 figures

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#, given the infrastructure in which they exist 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. 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. Two Cosmic Explorer observatories are indispensable for achieving precise localization of binary neutron star events, facilitating detection of electromagnetic counterparts and transforming multimessenger astronomy. Their combined operation is even more important in the detection and localization of high-redshift sources, such as binary neutron stars, beyond the star-formation peak, and primordial black hole mergers, which may occur roughly 100 million years after the Big Bang. The addition of the Einstein Telescope to a network of two Cosmic Explorer observatories is critical for accomplishing all the identified science metrics. 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: 48 pages, 20 figures, 14 tables

In this paper, we explore the prospect for improving the measurement accuracy of masses and radii of neutron stars. We consider imminent and long-term upgrades of the Laser Interferometer Gravitational-Wave Observatory (LIGO) and Virgo, as well as next-generation observatories -- the Cosmic Explorer and Einstein Telescope. We find that neutron star radius with single events will be constrained to within roughly 500m with the current generation of detectors and their upgrades. This will improve to 200m, 100m and 50m with a network of observatories that contain one, two or three next-generation observatories, respectively. Combining events in bins of 0.05 solar masses we find that for stiffer (softer) equations-of-state like ALF2 (APR4), a network of three XG observatories will determine the radius to within 30m (100m) over the entire mass range of neutron stars from 1 to 2.0 solar masses (2.2 solar masses), allowed by the respective equations-of-state. Neutron star masses will be measured to within 0.5 percent with three XG observatories irrespective of the actual equation-of-state. Measurement accuracies will be a factor of 4 or 2 worse if the network contains only one or two XG observatories, respectively, and a factor of 10 worse in the case of networks consisting of Advanced LIGO, Virgo KAGRA and their upgrades. Tens to hundreds of high-fidelity events detected by future observatories will allow us to accurately measure the mass-radius curve and hence determine the dense matter equation-of-state to exquisite precision., Comment: Updated version with additional supporting text

Gravitational-wave astronomy has revolutionized humanity's view of the universe, a revolution driven by observations that no other field can make. This white paper describes an observatory that builds on decades of investment by the National Science Foundation and that will drive discovery for decades to come: Cosmic Explorer. Major discoveries in astronomy are driven by three related improvements: better sensitivity, higher precision, and opening new observational windows. Cosmic Explorer promises all three and will deliver an order-of-magnitude greater sensitivity than LIGO. Cosmic Explorer will push the gravitational-wave frontier to almost the edge of the observable universe using technologies that have been proven by LIGO during its development. With the unprecedented sensitivity that only a new facility can deliver, Cosmic Explorer will make discoveries that cannot yet be anticipated, especially since gravitational waves are both synergistic with electromagnetic observations and can reach into regions of the universe that electromagnetic observations cannot explore. With Cosmic Explorer, scientists can use the universe as a laboratory to test the laws of physics and study the nature of matter. Cosmic Explorer allows the United States to continue its leading role in gravitational-wave science and the international network of next-generation observatories. With its extraordinary discovery potential, Cosmic Explorer will deliver revolutionary observations across astronomy, physics, and cosmology including: Black Holes and Neutron Stars Throughout Cosmic Time, Multi-Messenger Astrophysics and Dynamics of Dense Matter, New Probes of Extreme Astrophysics, Fundamental Physics and Precision Cosmology, Dark Matter and the Early Universe.

A direct detection of black hole formation in neutron star mergers would provide invaluable information about matter in neutron star cores and finite temperature effects on the nuclear equation of state. We study black hole formation in neutron star mergers using a set of 190 numerical relativity simulations consisting of long-lived and black-hole-forming remnants. The postmerger gravitational-wave spectrum of a long-lived remnant has greatly reduced power at a frequency $f$ greater than $f_{\rm peak}$, for $f \gtrsim 4\,\rm kHz$, with $f_{\rm peak} \in [2.5, 4]\,\rm kHz$. On the other hand, black-hole-forming remnants exhibit excess power in the same large $f$ region and manifest exponential damping in the time domain characteristic of a quasinormal mode. We demonstrate that the gravitational-wave signal from a collapsed remnant is indeed a quasinormal ringing. We report on the opportunity for direct detections of black hole formation with next-generation gravitational-wave detectors such as Cosmic Explorer and Einstein Telescope and set forth the tantalizing prospect of such observations up to a distance of 100 Mpc for an optimally oriented and located source with an SNR of 4.

Observations by the current generation of gravitational-wave detectors have been pivotal in expanding our understanding of the universe. Although tens of exciting compact binary mergers have been observed, neutron star-black hole (NSBH) mergers remained elusive until they were first confidently detected in 2020. The number of NSBH detections is expected to increase with sensitivity improvements of the current detectors and the proposed construction of new observatories over the next decade. In this work, we explore the NSBH detection and measurement capabilities of these upgraded detectors and new observatories using the following metrics: network detection efficiency and detection rate as a function of redshift, distributions of the signal-to-noise ratios, the measurement accuracy of intrinsic and extrinsic parameters, the accuracy of sky position measurement, and the number of early-warning alerts that can be sent to facilitate the electromagnetic follow-up. Additionally, we evaluate the prospects of performing multi-messenger observations of NSBH systems by reporting the number of expected kilonova detections with the Vera C. Rubin Observatory and the Nancy Grace Roman Space Telescope. We find that as many as $\mathcal{O}(10)$ kilonovae can be detected by these two telescopes every year, depending on the population of the NSBH systems and the equation of state of neutron stars., Comment: 30 pages, 15 figures, 15 tables

Our limited understanding of the physical properties of matter at ultra-high density, high proton/neutron number asymmetry, and low temperature is presently one of the major outstanding problems in physics. As matter in this extreme state is known to only exist stably in the cores of neutron stars (NSs), complementary measurements from electromagnetic and gravitational wave astrophysical observations of NSs, combined with terrestrial laboratory constraints and further theoretical investigations, hold the promise to provide important insight into the properties of matter in a region of the quantum chromodynamics phase space that is otherwise inaccessible. This multidisciplinary endeavor imposes the following requirements for facilities and resources in the upcoming decade and beyond: * A next generation of gravitational wave detectors to uncover more double NS and neutron star-black hole mergers; * Sensitive radio telescopes to find the most massive and fastest spinning NSs; * Large-area, high-time-resolution and/or high angular resolution X-ray telescopes to constrain the NS mass-radius relation; * Suitable laboratory facilities for nuclear physics experiments to constrain the dense matter equation of state; * Funding resources for theoretical studies of matter in this regime; * The availability of modern large-scale high performance computing infrastructure. The same facilities and resources would also enable significant advances in other high-profile fields of inquiry in modern physics such as the nature of dark matter, alternative theories of gravity, nucleon superfluidity and superconductivity, as well as an array of astrophysics, including but not limited to stellar evolution, nucleosynthesis, and primordial black holes., Comment: Submitted to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021) under Cosmic Frontier (CF07: Cosmic probes of fundamental physics); 30 pages, 8 figures

Inference of the equation-of-state (EoS) of dense nuclear matter in neutron-star cores is a principal science goal of X-ray and gravitational-wave observations of neutron stars. In particular, gravitational-wave observations provide an independent probe of the properties of bulk matter in neutron star cores that can then be used to compare with theoretically derived equations of state. In this paper, we quantify the systematic errors arising from the application of EoS-independent \emph{quasi-universal relations} in the estimation of neutron star tidal deformabilities and radii from gravitational-wave measurements and introduce a strategy to correct for the systematic biases in the inferred radii. We apply this method to a simulated population of events expected to be observed by future upgrades of current detectors and the next-generation of ground-based observatories. We show that our approach can accurately correct for the systematic biases arising from approximate universal relations in the mass-radius curves of neutron stars. Using the posterior distributions of the mass and radius for the simulated population we infer the underlying EoS with a good degree of precision. Our method revives the possibility of using the universal relations for rapid Bayesian model selection of dense matter EoS in gravitational-wave observations., Comment: accepted PRD

The origins of type Ia supernovae (SNe Ia) are still debated. Some of the leading scenarios involve a double detonation in double white dwarf (WD) systems. In these scenarios, helium shell detonation occurs on top of a carbon-oxygen (CO) WD, which then drives the detonation of the CO-core, producing a SN Ia. Extensive studies have been done on the possibility of a double helium detonation, following a dynamical helium mass-transfer phase onto a CO-WD. However, 3D self-consistent modeling of the double-WD system, the mass transfer, and the helium shell detonation have been little studied. Here we use 3D hydrodynamical simulations to explore this case in which a helium detonation occurs near the point of Roche lobe overflow of the donor WD and may lead to an SN Ia through the dynamically driven double-degenerate double-detonation (D6) mechanism. We find that the helium layer of the accreting primary WD does undergo a detonation, while the underlying carbon-oxygen core does not, leading to an extremely rapid and faint nova-like transient instead of a luminous SN Ia event. This failed core detonation suggests that D6 SNe Ia may be restricted to the most massive carbon-oxygen primary WDs. We highlight the nucleosynthesis of the long-lived radioisotope $^{44}$Ti during explosive helium burning, which may serve as a hallmark both of successful as well as failed D6 events which subsequently detonate as classical double-degenerate mergers., Comment: 11 pages, 4 figures, 1 table. Published in the Astrophysical Journal Letters

We develop a method to compute synthetic kilonova light curves that combines numerical relativity simulations of neutron star mergers and the $\texttt{SNEC}$ radiation-hydrodynamics code. We describe our implementation of initial and boundary conditions, r-process heating, and opacities for kilonova simulations. We validate our approach by carefully checking that energy conservation is satisfied and by comparing the $\texttt{SNEC}$ results with those of two semi-analytic light curve models. We apply our code to the calculation of color light curves for three binaries having different mass ratios (equal and unequal mass) and different merger outcome (short-lived and long-lived remnants). We study the sensitivity of our results to hydrodynamic effects, nuclear physics uncertainties in the heating rates, and duration of the merger simulations. We find that hydrodynamics effects are typically negligible and that homologous expansion is a good approximation in most cases. However, pressure forces can amplify the impact of uncertainties in the radioactive heating rates. We also study the impact of shocks possibly launched into the outflows by a relativistic jet. None of our models match AT2017gfo, the kilonova in GW170817. This points to possible deficiencies in our merger simulations and kilonova models which neglect non-LTE effects and possible additional energy injection from the merger remnant, and to the need to go beyond the assumption of spherical symmetry adopted in this work.

We determine the threshold mass for prompt (no bounce) black hole formation in equal-mass neutron star (NS) mergers using a new set of 227 numerical relativity simulations. We consider 23 phenomenological and microphysical finite temperature equations of state (EOS), including models with hyperons and first-order phase transitions to deconfined quarks. We confirm the existence of EOS-insensitive relations between the threshold mass, the binary tidal parameter at the threshold ($\Lambda_{th}$), the maximum mass of nonrotating NSs, and the radii of reference mass NSs. We correct the systematic errors in previously reported fitting coefficients that were obtained with approximate general-relativity simulations. We combine the EOS-insensitive relations, phenomenological constraints on NS properties and observational data from GW170817 to derive an improved lower limit on radii of maximum mass and 1.6 M$_\odot$ NS of 9.81 km and 10.90 km, respectively. We also constrain the radius and quadrupolar tidal deformability ($\Lambda$) of a 1.4 $M_\odot$ NS to be larger than 10.74 km and 172, respectively. We consider uncertainties in all independent parameters -- fitting coefficients as well as GW170817 masses while reporting the range of radii constraints. We introduce new methods to constrain the upper as well as lower limit of NS maximum mass using future BNS detections and their identification as prompt or delayed collapse. With future observations it will be possible to derive even tighter constraints on the properties of matter at and above nuclear density using the method proposed in this work., Comment: Accepted for publication in Physical Review D

This Horizon Study describes a next-generation ground-based gravitational-wave observatory: Cosmic Explorer. With ten times the sensitivity of Advanced LIGO, Cosmic Explorer will push gravitational-wave astronomy towards the edge of the observable universe ($z \sim 100$). The goals of this Horizon Study are to describe and evaluate design concepts for Cosmic Explorer; to plan for the United States' leadership in gravitational-wave astronomy; and to envisage the role of Cosmic Explorer in the international effort to build a "Third-Generation" (3G) observatory network that will make discoveries transformative across astronomy, physics, and cosmology., Comment: 173 pages

(abridged) We investigate the quark deconfinement phase transition in the context of binary neutron star (BNS) mergers. We employ a new finite-temperature composition-dependent equation of state (EOS) with a first order phase transition between hadrons and deconfined quarks to perform numerical relativity simulations of BNS mergers. The softening of the EOS due to the phase transition causes the merger remnants to be more compact and to collapse to a black hole (BH) at earlier times. The phase transition is imprinted on the postmerger gravitational wave (GW) signal duration, amplitude, and peak frequency. However, this imprint is only detectable for binaries with sufficiently long-lived remnants. Moreover, the phase transition does not result in significant deviations from quasi-universal relations for the postmerger GW peak frequency. We also study the impact of the phase transition on dynamical ejecta, remnant accretion disk masses, r-process nucleosynthetic yields and associated electromagnetic (EM) counterparts. While there are differences in the EM counterparts and nucleosynthesis yields between the purely hadronic models and the models with phase transitions, these can be primarily ascribed to the difference in remnant collapse time between the two. An exception is the non-thermal afterglow caused by the interaction of the fastest component of the dynamical ejecta and the interstellar medium, which is systematically boosted in the binaries with phase transition as a consequence of the more violent merger they experience., Comment: 28 pages and 16 figures

Extensive dynamical simulations are used to explore the possible existence of sudden sufficiently large energy or rogue fluctuations (RF) at late times and across short time windows in the {\it strongly nonlinear regime} of the $\beta$-Fermi-Pasta-Ulam-Tsingou (FPUT) type systems. Our studies build on a study of RF in the non-dissipative granular chain system and suggest that {\it rare RF may be generic to non-integrable, strongly nonlinear systems at late enough times}. We comment on the role of initial conditions and the intriguing influence of harmonic forces on these strongly nonlinear systems. The RF under focus here are distinct from the well known Peregrine solitons used to describe rogue waves via the weakly nonlinear Schr\"odinger equation., Comment: 16 pages, 5 figures

In this paper, we use large reaction networks to find ignition conditions of single-zone nuclear fuel with compositions typical of white dwarf (WD) matter. The necessary but, possibly not sufficient condition for initiation of detonation is that the nuclear burning proceeds on smaller timescale than the sound-crossing timescale. Under typical white dwarf thermodynamic conditions, the nuclear timescale depends sensitively on the chosen nuclear reaction network, and widely-used analytical formulae are not sufficient to accurately determine unstable ignition conditions the answer to which lies in the fully compressible reactive turbulent simulations. We model the nuclear reactions in the numerical simulation using the astrophysical code TORCH. We report the differences introduced due to the size of the nuclear reaction network. Our findings have implications for numerical results obtained using the popular 19-species network in multidimensional simulations of WD of Chandrasekhar and sub-Chandrasekhar masses. We present a sufficiently accurate general criterion in temperature-density ($T-\rho$) plane for a given initial fuel to ignite in an unstable fashion., Comment: Comments are invited

Binary neutron star mergers have been recently confirmed to be the progenitors of the optical transients kilonovae (KNe). KNe are powered by the radioactive decay of neutron-rich elements (r-process elements) which are believed to be the product of disruption of neutron stars during their merger. KNe exhibit interesting parallels with type Ia supernovae (SNe), whose light curves show specific correlations which allow them to be used as standardizable candles. In this paper, we investigate the possibility of the KN light curves exhibiting similar correlations. While a satisfactory answer to this question can only be provided by future KN observations, employing theoretical models we explore whether there is any ground for harboring such expectations. Using semi-analytic models of KN light curves in conjunction with results from numerical relativity simulations of binary neutron star mergers, we obtain the maximum bolometric luminosity ($L_{\mathrm{Bol}}^{\mathrm{max}}$) and decline in luminosity ($\Delta L_{\mathrm{Bol}}$) for a simulated population of mergers. We find that theoretical light curves of KNe show remarkable correlations despite the complex physics governing their behavior. This presents a possibility of future observations to uncover such correlations in the observed light curves, eventually allowing observers to standardize these light curves and to use them for local distance measurements., Comment: 3 pages; 3 figures; presented in GR22, Valencia