Numerical Relativity Simulations Research Articles

Mapping eccentricity evolutions between numerical relativity and effective-one-body gravitational waveforms

Orbital eccentricity in compact binaries is considered to be a key tracer of their astrophysical origin, and can be inferred from gravitational-wave observations due to its imprint on the emitted signal. For a robust measurement, accurate waveform models are needed. However, ambiguities in the definition of eccentricity can obfuscate the physical meaning and result in seemingly discrepant measurements. In this work we present a suite of 28 new numerical relativity simulations of eccentric, aligned-spin binary black holes with mass ratios between 1 and 6 and initial post-Newtonian eccentricities between 0.05 and 0.3. We then develop a robust pipeline for measuring the eccentricity evolution as a function of frequency from gravitational-wave observables that is applicable even to signals that span at least ≳7 orbits. We assess the reliability of our procedure and quantify its robustness under different assumptions on the data. Using the eccentricity measured at the first apastron, we initialize effective-one-body waveforms and quantify how the precision in the eccentricity measurement, and therefore the choice of the initial conditions, impacts the agreement with the numerical data. We find that even small deviations in the initial eccentricity can lead to non-negligible differences in the phase and amplitude of the waveforms. However, we demonstrate that we can reliably map the eccentricities between the simulation data and analytic models, which is crucial for robustly building eccentric hybrid waveforms, and to improve the accuracy of eccentric waveform models in the strong-field regime. Published by the American Physical Society 2024

A Review of Gravitational Memory and BMS Frame Fixing in Numerical Relativity

Abstract Gravitational memory effects and the BMS freedoms exhibited at future null infinity have recently been resolved and utilized in numerical relativity simulations. With this, gravitational wave models and our understanding of the fundamental nature of general relativity have been vastly improved. In this paper, we review the history and intuition behind memory effects and BMS symmetries, how they manifest in gravitational waves, and how controlling the infinite number of BMS freedoms of numerical relativity simulations can crucially improve the waveform models that are used by gravitational wave detectors. We reiterate the fact that, with memory effects and BMS symmetries, not only can these next-generation numerical waveforms be used to observe never-before-seen physics, but they can also be used to test GR and learn new astrophysical information about our universe.

Self-interacting scalar dark matter around binary black holes

Gravitational waves can provide crucial insights about the environments in which black holes live. In this work, we use numerical relativity simulations to study the behavior of self-interacting scalar (wavelike) dark matter clouds accreting onto isolated and binary black holes. We find that repulsive self-interactions smoothen the “spike” of an isolated black hole and saturate the density. Attractive self-interactions enhance the growth and result in more cuspy profiles, but can become unstable and undergo explosions akin to the superradiant bosenova that reduce the local cloud density. We quantify the impact of self-interactions on an equal-mass black hole merger by computing the dephasing of the gravitational-wave signal for a range of couplings. We find that repulsive self-interactions saturate the density of the cloud, thereby reducing the dephasing. For attractive self-interactions, the dephasing may be larger, but if these interactions dominate prior to the merger, the dark matter can undergo bosenova during the inspiral phase, disrupting the cloud and subsequently reducing the dephasing. Published by the American Physical Society 2024

Gravitational-Wave Data Analysis with High-Precision Numerical Relativity Simulations of Boson Star Mergers.

Gravitational-wave signals detected to date are commonly interpreted under the paradigm that they originate from pairs of black holes or neutron stars. Here, we explore the alternative scenario of boson-star signals being present in the data stream. We perform accurate and long (∼20 orbits) numerical simulations of boson-star binaries and inject the resulting strain into LIGO noise. Our Bayesian inference reveals that some boson-star signals are degenerate with current approximants, albeit with biased parameters, while others exhibit smoking-gun signatures leaving behind conspicuous residuals.

Radioactive Gamma-Ray Lines from Long-lived Neutron Star Merger Remnants

The observation of a kilonova AT2017gfo associated with the gravitational wave event GW170817 provides the first strong evidence that neutron star mergers are dominant contributors to the production of heavy r-process elements. Radioactive gamma-ray lines emitted from neutron star merger remnants provide a unique probe for investigating the nuclide composition and tracking its evolution. In this work, we studied the gamma-ray line features arising from the radioactive decay of heavy nuclei in the merger remnants based on the r-process nuclear reaction network and the astrophysical inputs derived from numerical relativity simulations. The decay chain of 50126 Sn (T 1/2 = 230 kyr) → 51126 Sb (T 1/2 = 12.35 days) → 52126 Te (stable) produces several bright gamma-ray lines with energies of 415, 667, and 695 keV, making it the most promising decay chain during the remnant phase. The photon fluxes of these bright gamma-ray lines reach ∼10−5 γ cm−2 s−1 for Galactic merger remnants with ages less than 100 kyr, which can be detected by the high energy resolution MeV gamma-ray detectors like the MASS mission.

The Impact of Anisotropic Sky Sampling on the Hubble Constant in Numerical Relativity

We study the impact of nearby inhomogeneities on an observer’s inference of the Hubble constant. Large-scale structures induce a dependence of cosmological parameters on observer position as well as an anisotropic variance of those parameters across an observer’s sky. While the former has been explored quite thoroughly, the latter has not. Incomplete sampling of an anisotropic sky could introduce a bias in our cosmological inference if we assume an isotropic expansion law. In this work, we use numerical relativity simulations of large-scale structure combined with ray tracing to produce synthetic catalogs mimicking the low-redshift Pantheon supernova (SN) data set. Our data contain all general-relativistic contributions to fluctuations in the distances and redshifts along geodesics in the simulation. We use these synthetic observations to constrain H 0 for a set of randomly positioned observers. We study both the dependence on observer position as well as the impact of rotating the sample of SNe on the observer’s sky. We find a 1%–2% variance in H 0 between observers when they use an isotropic sample of objects. However, we find the inferred value of H 0 can change by up to 3%–6% when observers simply rotate their Pantheon data set on the sky. While the variances we find are below the level of the “Hubble tension,” our results may suggest a reduction in the significance of the tension if anisotropy of expansion can be correctly accounted for.

RisingTides: An Analytic Modeling Code of Tidal Effects in Binary Neutron Star Mergers

Gravitational waves produced by binary neutron star mergers offer a unique window into matter behavior under extreme conditions. In this context, we analytically model the effect of matter on gravitational waves from binary neutron star mergers. We start with a binary black hole system, leveraging the post-Newtonian formalism for the inspiral and the Backwards-one-Body model for the merger. We combine the two methods to generate a baseline waveform and we validate our results against numerical relativity simulations. Next, we integrate tidal effects in phase and amplitude to account for matter and spacetime interaction using the NRTidal model and test its accuracy against numerical relativity predictions for two equations of state, finding a mismatch around the merger. Subsequently, we lift the restriction on the coefficients to be independent of the tidal deformability and recalibrate them using the numerical relativity predictions. We obtain better fits for phase and amplitude around the merger and are able to extend the phase modeling beyond the merger. We implement our method in new open-source, user-friendly Python code, steered by a Jupyter Notebook, named RisingTides. Our research offers new perspectives on analytically modeling the effect of tides on the gravitational waves from binary neutron star mergers.

Numerical relativity simulations of black hole and relativistic jet formation

ABSTRACT We investigate impacts of stellar rotation and magnetic fields on black hole (BH) formation and its subsequent explosive activities, by conducting axisymmetric radiation-magnetohydrodynamics simulations of gravitational collapse of a 70 $\mathrm{M}_\odot$ star with two-moment multi energy neutrino transport in full general relativity for the first time. Due to its dense stellar structure, all models cannot avoid the eventual BH formation even though a strongly magnetized model experiences the so-called magnetorotational explosion prior to the BH formation. One intriguing phenomenon observed in the strongly magnetized model is the formation of a relativistic jet in the post-BH formation. The relativistic jet is the outcome of a combination of strong magnetic fields and low-density materials above the BH. The jet further enhances the explosion energy beyond $\sim 10^{52}$ erg, which is well exceeding the gravitational overburden ahead of the shock. Our self-consistent supernova models demonstrate that rotating magnetized massive stars at the high-mass end of supernova progenitors could be a potential candidate of hypernova and long gamma-ray burst progenitors.

On primordial universe in anti-de Sitter landscape

How the spacetime evolved non-perturbatively in a landscape with multiple anti-de Sitter (AdS) vacua, which is theoretically well-motivated, has always been a matter of concern. As a step towards this issue, we perform (3+1)D numerical relativity simulations for the inhomogeneous universe in an AdS landscape, and find that large inhomogeneity of scalar field can develop into not only sphere-like bubbles, but also novel tube-like structures. It is observed that the bubble or tube wall (across the potential barrier) likely inflates as a quasi-dS space, while the different regions separated by the walls are in different AdS vacua and will collapse towards singularity.

Probabilistic model for the gravitational wave signal from merging black holes

Parametrized models that predict the gravitational-wave (GW) signal from merging black holes are used to extract source properties from GW observations. The majority of research in this area has focused on developing methods capable of producing highly accurate, point estimate, predictions for the GW signal. A key element missing from every model used in the analysis of GW data is an estimate for how confident the model is in its prediction. This omission increases the risk of biased parameter estimation of source properties. Current strategies include running analyses with multiple models to measure systematic bias however, this fails to accurately reflect the true uncertainty in the models. In this work we develop a probabilistic extension to the phenomenological modeling workflow for nonspinning black holes and demonstrate that the model not only produces accurate point estimates for the GW signal but can be used to provide well-calibrated local estimates for its uncertainty. Our analysis highlights that there is a lack of numerical relativity (NR) simulations available at multiple resolutions which can be used to estimate their numerical error and implore the NR community to continue to improve their estimates for the error in NR solutions published. Waveform models that are not only accurate in their point-estimate predictions but also in their error estimates are a potential way to mitigate bias in GW parameter estimation of compact binaries due to unconfident waveform model extrapolations. Published by the American Physical Society 2024

MRI turbulence in vertically stratified accretion discs at large magnetic Prandtl numbers

ABSTRACT The discovery of the first binary neutron star merger, GW170817, has spawned a plethora of global numerical relativity simulations. These simulations are often ideal (with dissipation determined by the grid) and/or axisymmetric (invoking ad hoc mean-field dynamos). However, binary neutron star mergers (similar to X-ray binaries and active galactic nuclei inner discs) are characterized by large magnetic Prandtl numbers, $\rm Pm$, (the ratio of viscosity to resistivity). $\rm Pm$ is a key parameter determining dynamo action and dissipation but it is ill-defined (and likely of order unity) in ideal simulations. To bridge this gap, we investigate the magnetorotational instability (MRI) and associated dynamo at large magnetic Prandtl numbers using fully compressible, three-dimensional, vertically stratified, isothermal simulations of a local patch of a disc. We find that, within the bulk of the disc (z ≲ 2H, where H is the scale-height), the turbulent intensity (parametrized by the stress-to-thermal-pressure ratio α), and the saturated magnetic field energy density, Emag, produced by the MRI dynamo, both scale as a power with Pm at moderate Pm (4 ≲ Pm ≲ 32): Emag ∼ Pm0.74 and α ∼ Pm0.71, respectively. At larger Pm (≳ 32), we find deviations from power-law scaling and the onset of a plateau. Compared to our recent unstratified study, this scaling with Pm becomes weaker further away from the disc mid-plane, where the Parker instability dominates. We perform a thorough spectral analysis to understand the underlying dynamics of small-scale MRI-driven turbulence in the mid-plane and of large-scale Parker-unstable structures in the atmosphere.

Comparison of neural network architectures for feature extraction from binary black hole merger waveforms

We evaluate several neural-network architectures, both convolutional and recurrent, for gravitational-wave time-series feature extraction by performing point parameter estimation on noisy waveforms from binary-black-hole mergers. We build datasets of 100 000 elements for each of four different waveform models (or approximants) in order to test how approximant choice affects feature extraction. Our choices include SEOBNRv4P and IMRPhenomPv3, which contain only the dominant quadrupole emission mode, alongside IMRPhenomPv3HM and NRHybSur3dq8, which also account for high-order modes. Each dataset element is injected into detector noise corresponding to the third observing run of the LIGO-Virgo-KAGRA (LVK) collaboration. We identify the temporal convolutional network architecture as the overall best performer in terms of training and validation losses and absence of overfitting to data. Comparison of results between datasets shows that the choice of waveform approximant for the creation of a dataset conditions the feature extraction ability of a trained network. Hence, care should be taken when building a dataset for the training of neural networks, as certain approximants may result in better network convergence of evaluation metrics. However, this performance does not necessarily translate to data which is more faithful to numerical relativity simulations. We also apply this network on actual signals from LVK runs, finding that its feature-extracting performance can be effective on real data.

Catalog of precessing black-hole-binary numerical-relativity simulations

Catalog of precessing black-hole-binary numerical-relativity simulations

What if GW190425 did not produce a black hole promptly?

ABSTRACT It is widely believed that the binary neutron star merger GW190425 produced a black hole promptly upon merger. Motivated by the potential association with the fast radio burst FRB 20190425A, which took place 2.5 h after the merger, we revisit the question of the outcome of GW190425 by means of numerical relativity simulations. We show that current laboratory and astrophysical constraints on the equation of state of dense matter do not rule out the formation of a long-lived remnant. However, the formation of a stable remnant would have produced a bright kilonova, in tension with upper limits by ZTF at the location and time of FRB 20190425A. Moreover, the ejecta would have been optically thick to radio emission for days to months, preventing a putative FRB from propagating out. The predicted dispersion measure is also several orders of magnitude larger than that observed for FRB 20190425A. Our results indicate that FRB 20190425A and GW190425 are not associated. However, we cannot completely rule out the formation of a long-lived remnant, due to the incomplete coverage of the relevant sky regions. More observations of GW190425-like events, including potential upper limit, have the potential to constrain nuclear physics. To this aim, it is important that follow-up observational campaigns of gravitational wave events are informed by the properties of the source, such as their chirp mass, and we urge the LIGO-Virgo-KAGRA collaboration to promptly release them publicly.

The skewness of the distance-redshift relation in ΛCDM

Starting from a recently proposed framework for the evaluation of the cosmological averages, we evaluate the higher-order moments for the distribution of a given observable. Then, we explicitly discuss the case of the Hubble-Lemaître diagram and evaluate its skewness at the leading order in the cosmological perturbative expansion of the gravitational potential. In particular, we focus on perturbations of the luminosity distance due to gravitational lensing. Finally, we discuss our findings in view of recent numerical relativistic simulations, confirming that the skewness in the Hubble-Lemaître diagram primarily originates from the late-time matter bispectrum, with other line-of-sight projection effects being sub-dominant.

What Can We Learn about the Unstable Equation-of-state Branch from Neutron Star Mergers?

The equation of state (EOS) of dense strongly interacting matter can be probed by astrophysical observations of neutron stars (NS), such as X-ray detections of pulsars or the measurement of the tidal deformability of NSs during the inspiral stage of NS mergers. These observations constrain the EOS at most up to the density of the maximum-mass configuration, n TOV, which is the highest density that can be explored by stable NSs for a given EOS. However, under the right circumstances, binary neutron star (BNS) mergers can create a postmerger remnant that explores densities above n TOV. In this work, we explore whether the EOS above n TOV can be measured from gravitational-wave or electromagnetic observations of the postmerger remnant. We perform a total of 25 numerical-relativity simulations of BNS mergers for a range of EOSs and find no case in which different descriptions of the matter above n TOV have a detectable impact on postmerger observables. Hence, we conclude that the EOS above n TOV can likely not be probed through BNS merger observations for the current and next generation of detectors.

Effects of onset of phase transition on binary neutron star mergers

ABSTRACT Quantum Chromodynamics predicts phase transition from hadronic matter to quark matter at high density, which is highly probable in astrophysical systems like binary neutron star mergers. To explore the critical density where such phase transition can occur, we performed numerical relativity simulations of binary neutron star mergers with various masses (equal and unequal binaries). We aim to understand the effect of the onset of phase transition on the merger dynamics and gravitational wave spectra. We generated a set of equations of states by agnostically changing the onset of phase transition, having the hadronic matter part and quark matter part fixed. This particular arrangement of the equation of states explores the scenario of mergers where mixed phases of matter are achieved before or during the merger. Under these circumstances, if the matter properties with hadronic and quark degrees differ significantly, it is reflected in the stability of the final merger product for the intermediate mass binary. We performed a case study on mixed species merger, where one of the binary companions is hybrid star. If quark matter appears at low densities, we observe significant change in post-merger gravitational wave analysis in terms of higher peak frequencies and post-merger frequencies in power spectral density. We report indications expressed as spikes in phase difference plots at merger time for mixed mergers. We found that the expression of phase transition in post-merger gravitational wave signals is more significant for unequal mass binary than for equal mass binary having the same total baryonic mass.

Collapse of Rotating Massive Stars Leading to Black Hole Formation and Energetic Supernovae

We explore a possible explosion scenario resulting from core collapses of rotating massive stars that leave a black hole by performing radiation-viscous-hydrodynamics simulations in numerical relativity. We take moderately and rapidly rotating compact pre-collapse stellar models with zero-age main-sequence masses of 9M ⊙ and 20M ⊙ based on stellar evolution calculations as the initial conditions. We find that viscous heating in the disk formed around the central black hole is the power source for an outflow. The moderately rotating models predict a small ejecta mass of the order of 0.1M ⊙ and an explosion energy of ≲1051 erg. Due to the small ejecta mass, these models may predict a short-timescale transient with a rise time of 3–5 days. This can lead to a bright (∼1044 erg s−1) transient, like superluminous supernovae in the presence of a dense massive circumstellar medium. For hypothetically rapidly rotating models that have a high mass-infall rate onto the disk, the explosion energy is ≳3 × 1051 erg, which is comparable to or larger than that of typical stripped-envelope supernovae, indicating that a fraction of such supernovae may be explosions powered by black hole accretion disks. The explosion energy is still increasing at the end of the simulations with a rate of >1050 erg s−1, and thus, it may reach ∼1052 erg. A nucleosynthesis calculation shows that the mass of 56Ni amounts to ≳0.1M ⊙, which, together with the high explosion energy, may satisfy the required amount for broad-lined type Ic supernovae. Irrespective of the models, the lowest value of the electron fraction of the ejecta is ≳0.4; thus, synthesis of heavy r-process elements is not found in our models.

Starting inflation from inhomogeneous initial conditions with momentum

We investigate the circumstances under which cosmic inflation can arise from very inhomogeneous initial conditions using numerical relativity simulations. Previous studies have not considered cases with non-zero momentum density due to technical challenges with solving the coupled Einstein constraint equations. Here we address these, introducing and comparing several different ways of constructing cosmological initial conditions with inhomogeneous scalar field and time derivative profiles. We evolve such initial conditions with large inhomogeneities in both single- and two-field inflationary models. We study cases where the initial gradient and kinetic energy are much larger than the inflationary energy scale, and black holes can form, as well as cases where the initial scalar potential energy is comparable, as in scenarios where inflation occurs at nearly Planckian densities, finding large-field inflation to be generally robust. We consider examples of initial conditions where a large scalar field velocity towards non-inflationary values would prevent inflation from occurring in the homogeneous case, finding that the addition of large gradients in the scalar field can actually dilute this effect, with the increased expansion and non-vanishing restoring force leading to inflation.

Hours-long Near-UV/Optical Emission from Mildly Relativistic Outflows in Black Hole–Neutron Star Mergers

The ongoing LIGO–Virgo–KAGRA observing run O4 provides an opportunity to discover new multimessenger events, including binary neutron star (BNS) mergers such as GW170817 and the highly anticipated first detection of a multimessenger black hole–neutron star (BH–NS) merger. While BNS mergers were predicted to exhibit early optical emission from mildly relativistic outflows, it has remained uncertain whether the BH–NS merger ejecta provides the conditions for similar signals to emerge. We present the first modeling of early near-ultraviolet/optical emission from mildly relativistic outflows in BH–NS mergers. Adopting optimal binary properties, a mass ratio of q = 2, and a rapidly rotating BH, we utilize numerical relativity and general relativistic magnetohydrodynamic (GRMHD) simulations to follow the binary’s evolution from premerger to homologous expansion. We use an M1 neutrino transport GRMHD simulation to self-consistently estimate the opacity distribution in the outflows and find a bright near-ultraviolet/optical signal that emerges due to jet-powered cocoon cooling emission, outshining the kilonova emission at early time. The signal peaks at an absolute magnitude of ∼−15 a few hours after the merger, longer than previous estimates, which did not consider the first principles–based jet launching. By late 2024, the Rubin Observatory will have the capability to track the entire signal evolution or detect its peak up to distances of ≳1 Gpc. In 2026, ULTRASAT will conduct all-sky surveys within minutes, detecting some of these events within ∼200 Mpc. The BH–NS mergers with higher mass ratios or lower BH spins would produce shorter and fainter signals.