BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Lawrence E. Kidder

Lawrence E. Kidder contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
gr-qcastro-ph.HEastro-ph.COhep-thphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
21works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

21 published item(s)

preprint2024arXiv

Nonlinear Effects In Black Hole Ringdown From Scattering Experiments I: spin and initial data dependence of quadratic mode coupling

We investigate quadratic quasinormal mode coupling in black hole spacetime through numerical simulations of single perturbed black holes using both numerical relativity and second-order black hole perturbation theory. Focusing on the dominant $\ell=|m|=2$ quadrupolar modes, we find good agreement (within $\sim10\%$) between these approaches, with discrepancies attributed to truncation error and uncertainties from mode fitting. Our results align with earlier studies extracting the coupling coefficients from select binary black hole merger simulations, showing consistency for the same remnant spins. Notably, the coupling coefficient is insensitive to a diverse range of initial data, including configurations that led to a significant (up to $5\%$) increase in the remnant black hole mass. These findings present opportunities for testing the nonlinear dynamics of general relativity with ground-based gravitational wave observatories. Lastly, we provide evidence of a bifurcation in coupling coefficients between counter-rotating and co-rotating quasinormal modes as black hole spin increases.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Numerical relativity surrogate model with memory effects and post-Newtonian hybridization

Numerical relativity simulations provide the most precise templates for the gravitational waves produced by binary black hole mergers. However, many of these simulations use an incomplete waveform extraction technique -- extrapolation -- that fails to capture important physics, such as gravitational memory effects. Cauchy-characteristic evolution (CCE), by contrast, is a much more physically accurate extraction procedure that fully evolves Einstein's equations to future null infinity and accurately captures the expected physics. In this work, we present a new surrogate model, NRHybSur3dq8$\_$CCE, built from CCE waveforms that have been mapped to the post-Newtonian (PN) BMS frame and then hybridized with PN and effective one-body (EOB) waveforms. This model is trained on 102 waveforms with mass ratios $q\leq8$ and aligned spins $χ_{1z}, \, χ_{2z} \in \left[-0.8, 0.8\right]$. The model spans the entire LIGO-Virgo-KAGRA (LVK) frequency band (with $f_{\text{low}}=20\text{Hz}$) for total masses $M\gtrsim2.25M_{\odot}$ and includes the $\ell\leq4$ and $(\ell,m)=(5,5)$ spin-weight $-2$ spherical harmonic modes, but not the $(3,1)$, $(4,2)$ or $(4,1)$ modes. We find that NRHybSur3dq8$\_$CCE can accurately reproduce the training waveforms with mismatches $\lesssim2\times10^{-4}$ for total masses $2.25M_{\odot}\leq M\leq300M_{\odot}$ and can, for a modest degree of extrapolation, capably model outside of its training region. Most importantly, unlike previous waveform models, the new surrogate model successfully captures memory effects.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A scalable elliptic solver with task-based parallelism for the SpECTRE numerical relativity code

Elliptic partial differential equations must be solved numerically for many problems in numerical relativity, such as initial data for every simulation of merging black holes and neutron stars. Existing elliptic solvers can take multiple days to solve these problems at high resolution and when matter is involved, because they are either hard to parallelize or require a large amount of computational resources. Here we present a new solver for linear and nonlinear elliptic problems that is designed to scale with resolution and to parallelize on computing clusters. To achieve this we employ a discontinuous Galerkin discretization, an iterative multigrid-Schwarz preconditioned Newton-Krylov algorithm, and a task-based parallelism paradigm. To accelerate convergence of the elliptic solver we have developed novel subdomain-preconditioning techniques. We find that our multigrid-Schwarz preconditioned elliptic solves achieve iteration counts that are independent of resolution, and our task-based parallel programs scale over 200 million degrees of freedom to at least a few thousand cores. Our new code solves a classic initial data problem for binary black holes faster than the spectral code SpEC when distributed to only eight cores, and in a fraction of the time on more cores. It is publicly accessible in the next-generation SpECTRE numerical relativity code. Our results pave the way for highly parallel elliptic solves in numerical relativity and beyond.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Eccentric binary black holes: Comparing numerical relativity and small mass-ratio perturbation theory

The modelling of unequal mass binary black hole systems is of high importance to detect and estimate parameters from these systems. Numerical relativity (NR) is well suited to study systems with comparable component masses, $m_1\sim m_2$, whereas small mass ratio (SMR) perturbation theory applies to binaries where $q=m_2/m_1<< 1$. This work investigates the applicability for NR and SMR as a function of mass ratio for eccentric non-spinning binary black holes. We produce $52$ NR simulations with mass ratios between $1:10$ and $1:1$ and initial eccentricities up to $0.7$. From these we extract quantities like gravitational wave energy and angular momentum fluxes and periastron advance, and assess their accuracy. To facilitate comparison, we develop tools to map between NR and SMR inspiral evolutions of eccentric binary black holes. We derive post-Newtonian accurate relations between different definitions of eccentricity. Based on these analyses, we introduce a new definition of eccentricity based on the (2,2)-mode of the gravitational radiation, which reduces to the Newtonian definition of eccentricity in the Newtonian limit. From the comparison between NR simulations and SMR results, we quantify the unknown next-to-leading order SMR contributions to the gravitational energy and angular momentum fluxes, and periastron advance. We show that in the comparable mass regime these contributions are subdominant and higher order SMR contributions are negligible.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Fixing the BMS Frame of Numerical Relativity Waveforms

Understanding the Bondi-Metzner-Sachs (BMS) frame of the gravitational waves produced by numerical relativity is crucial for ensuring that analyses on such waveforms are performed properly. It is also important that models are built from waveforms in the same BMS frame. Up until now, however, the BMS frame of numerical waveforms has not been thoroughly examined, largely because the necessary tools have not existed. In this paper, we show how to analyze and map to a suitable BMS frame for numerical waveforms calculated with the Spectral Einstein Code (SpEC). However, the methods and tools that we present are general and can be applied to any numerical waveforms. We present an extensive study of 13 binary black hole systems that broadly span parameter space. From these simulations, we extract the strain and also the Weyl scalars using both SpECTRE&#39;s Cauchy-characteristic extraction module and also the standard extrapolation procedure with a displacement memory correction applied during postprocessing. First, we show that the current center-of-mass correction used to map these waveforms to the center-of-mass frame is not as effective as previously thought. Consequently, we also develop an improved correction that utilizes asymptotic Poincaré charges instead of a Newtonian center-of-mass trajectory. Next, we map our waveforms to the post-Newtonian (PN) BMS frame using a PN strain waveform. This helps us find the unique BMS transformation that minimizes the $L^{2}$ norm of the difference between the numerical and PN strain waveforms during the early inspiral phase. We find that once the waveforms are mapped to the PN BMS frame, they can be hybridized with a PN strain waveform much more effectively than if one used any of the previous alignment schemes, which only utilize the Poincaré transformations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Gravitational-wave echoes from numerical-relativity waveforms via space-time construction near merging compact objects

We propose a new approach toward reconstructing the late-time near-horizon geometry of merging binary black holes, and toward computing gravitational-wave echoes from exotic compact objects. A binary black-hole merger spacetime can be divided by a time-like hypersurface into a Black-Hole Perturbation (BHP) region, in which the space-time geometry can be approximated by homogeneous linear perturbations of the final Kerr black hole, and a nonlinear region. At late times, the boundary between the two regions is an infalling shell. The BHP region contains late-time gravitational-waves emitted toward the future horizon, as well as those emitted toward future null infinity. In this region, by imposing no-ingoing wave conditions at past null infinity, and matching out-going waves at future null infinity with waveforms computed from numerical relativity, we can obtain waves that travel toward the future horizon. In particular, the Newman-Penrose $ψ_0$ associated with the in-going wave on the horizon is related to tidal deformations measured by fiducial observers floating above the horizon. We further determine the boundary of the BHP region on the future horizon by imposing that $ψ_0$ inside the BHP region can be faithfully represented by quasi-normal modes. Using a physically-motivated way to impose boundary conditions near the horizon, and applying the so-called Boltzmann reflectivity, we compute the quasi-normal modes of non-rotating ECOs, as well as gravitational-wave echoes. We also investigate the detectability of these echoes in current and future detectors, and prospects for parameter estimation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

High Precision Ringdown Modeling: Multimode Fits and BMS Frames

Quasi-normal mode (QNM) modeling is an invaluable tool for characterizing remnant black holes, studying strong gravity, and testing GR. Only recently have QNM studies begun to focus on multimode fitting to numerical relativity (NR) strain waveforms. As GW observatories become even more sensitive they will be able to resolve higher-order modes. Consequently, multimode QNM fits will be critically important, and in turn require a more thorough treatment of the asymptotic frame at $\mathscr{I}^+$. The first main result of this work is a method for systematically fitting a QNM model containing many modes to a numerical waveform produced using Cauchy-characteristic extraction (CCE), an extraction technique which is known to resolve memory effects. We choose the modes to model based on their power contribution to the residual between numerical and model waveforms. We show that the all-mode strain mismatch improves by a factor of $\sim10^5$ when using multimode fitting as opposed to only fitting the $(2,\pm2,n)$ modes. Our most significant result addresses a critical point that has been overlooked in the QNM literature: the importance of matching the Bondi-van der Burg-Metzner-Sachs (BMS) frame of the numerical waveform to that of the QNM model. We show that by mapping the numerical waveforms$-$which exhibit the memory effect$-$to a BMS frame known as the super rest frame, there is an improvement of $\sim10^5$ in the all-mode strain mismatch compared to using a strain waveform whose BMS frame is not fixed. Furthermore, we find that by mapping CCE waveforms to the super rest frame, we can obtain all-mode mismatches that are, on average, a factor of $\sim4$ better than using the publicly-available extrapolated waveforms. We illustrate the effectiveness of these modeling enhancements by applying them to families of waveforms produced by NR and comparing our results to previous QNM studies.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Multipole moments on the common horizon in a binary-black-hole simulation

We construct the covariantly defined multipole moments on the common horizon of an equal-mass, non-spinning, quasicircular binary-black-hole system. We see a strong correlation between these multipole moments and the gravitational waveform. We find that the multipole moments are well described by the fundamental quasinormal modes at sufficiently late times. For each multipole moment, at least two fundamental modes of different $\ell$ are detectable in the best model. These models provide faithful estimates of the true mass and spin of the remnant black hole. We also show that by including overtones, the $\ell=m=2$ mass multipole moment admits an excellent quasinormal-mode description at all times after the merger. This demonstrates the perhaps surprising power of perturbation theory near the merger.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Simulating magnetized neutron stars with discontinuous Galerkin methods

Discontinuous Galerkin methods are popular because they can achieve high order where the solution is smooth, because they can capture shocks while needing only nearest-neighbor communication, and because they are relatively easy to formulate on complex meshes. We perform a detailed comparison of various limiting strategies presented in the literature applied to the equations of general relativistic magnetohydrodynamics. We compare the standard minmod/$ΛΠ^N$ limiter, the hierarchical limiter of Krivodonova, the simple WENO limiter, the HWENO limiter, and a discontinuous Galerkin-finite-difference hybrid method. The ultimate goal is to understand what limiting strategies are able to robustly simulate magnetized TOV stars without any fine-tuning of parameters. Among the limiters explored here, the only limiting strategy we can endorse is a discontinuous Galerkin-finite-difference hybrid method.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Targeted large mass ratio numerical relativity surrogate waveform model for GW190814

Gravitational wave observations of large mass ratio compact binary mergers like GW190814 highlight the need for reliable, high-accuracy waveform templates for such systems. We present NRHybSur2dq15, a new surrogate model trained on hybridized numerical relativity (NR) waveforms with mass ratios $q\leq15$, and aligned spins $|χ_{1z}|\leq0.5$ and $χ_{2z}=0$. We target the parameter space of GW190814-like events as large mass ratio NR simulations are very expensive. The model includes the (2,2), (2,1), (3,3), (4,4), and (5,5) spin-weighted spherical harmonic modes, and spans the entire LIGO bandwidth (with $f_{\mathrm{low}}=20$ Hz) for total masses $M \gtrsim 9.5 \, M_{\odot}$. NRHybSur2dq15 accurately reproduces the hybrid waveforms, with mismatches below $\sim 2 \times 10^{-3}$ for total masses $10 \, M_{\odot} \leq M \leq 300 \, M_{\odot}$. This is at least an order of magnitude improvement over existing semi-analytical models for GW190814-like systems. Finally, we reanalyze GW190814 with the new model and obtain source parameter constraints consistent with previous work.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Adding Gravitational Memory to Waveform Catalogs using BMS Balance Laws

Accurate models of gravitational waves from merging binary black holes are crucial for detectors to measure events and extract new science. One important feature that is currently missing from the Simulating eXtreme Spacetimes (SXS) Collaboration&#39;s catalog of waveforms for merging black holes, and other waveform catalogs, is the gravitational memory effect: a persistent, physical change to spacetime that is induced by the passage of transient radiation. We find, however, that by exploiting the Bondi-Metzner-Sachs (BMS) balance laws, which come from the extended BMS transformations, we can correct the strain waveforms in the SXS catalog to include the missing displacement memory. Our results show that these corrected waveforms satisfy the BMS balance laws to a much higher degree of accuracy. Furthermore, we find that these corrected strain waveforms coincide especially well with the waveforms obtained from Cauchy-characteristic extraction (CCE) that already exhibit memory effects. These corrected strain waveforms also evade the transient junk effects that are currently present in CCE waveforms. Lastly, we make our code for computing these contributions to the BMS balance laws and memory publicly available as a part of the python package $\texttt{sxs}$, thus enabling anyone to evaluate the expected memory effects and violation of the BMS balance laws.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

An aligned-spin neutron-star--black-hole waveform model based on the effective-one-body approach and numerical-relativity simulations

After the discovery of gravitational waves from binary black holes (BBHs) and binary neutron stars (BNSs) with the LIGO and Virgo detectors, neutron-star--black-holes (NSBHs) are the natural next class of binary systems to be observed. In this work, we develop a waveform model for aligned-spin neutron-star--black-holes (NSBHs) combining a BBH baseline waveform (available in the effective-one-body approach) with a phenomenological description of tidal effects (extracted from numerical-relativity simulations), and correcting the amplitude during the late inspiral, merger and ringdown to account for the NS tidal disruption. We calibrate the amplitude corrections using NSBH waveforms obtained with the SpEC and the SACRA codes. The model was calibrated using simulations with NS masses in the range $1.2-1.4 M_\odot$, tidal deformabilities up to $4200$ (for a 1.2 $M_\odot$ NS), and dimensionless BH spin magnitude up to 0.9. Based on the simulations used, and on checking that sensible waveforms are produced, we recommend our model to be employed with NS mass in the range $1\mbox{--}3 M_\odot$, tidal deformability $0\mbox{--}5000$, and BH spin magnitude up to $0.9$. We also validate our model against two new, highly accurate NSBH waveforms with BH spin 0.9 and mass ratios 3 and 4, characterized by tidal disruption, produced with SpEC, and find very good agreement. We find that it will be challenging for the advanced LIGO-Virgo--detector network at design sensitivity to distinguish different source classes. We perform parameter-estimation on a synthetic numerical-relativity signal in zero noise to study parameter biases. Finally, we reanalyze GW170817, with the hypothesis that it is a NSBH. We do not find evidence to distinguish the BNS and NSBH hypotheses, however the posterior for the mass ratio is shifted to less equal masses under the NSBH hypothesis. [Abstract abridged for arxiv].

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Computation of Displacement and Spin Gravitational Memory in Numerical Relativity

We present the first numerical relativity waveforms for binary black hole mergers produced using spectral methods that show both the displacement and the spin memory effects. Explicitly, we use the SXS Collaboration&#39;s $\texttt{SpEC}$ code to run a Cauchy evolution of a binary black hole merger and then extract the gravitational wave strain using $\texttt{SpECTRE}$&#39;s version of a Cauchy-characteristic extraction. We find that we can accurately resolve the strain&#39;s traditional $m=0$ memory modes and some of the $m\not=0$ oscillatory memory modes that have previously only been theorized. We also perform a separate calculation of the memory using equations for the Bondi-Metzner-Sachs charges as well as the energy and angular momentum fluxes at asymptotic infinity. Our new calculation uses only the gravitational wave strain and two of the Weyl scalars at infinity. Also, this computation shows that the memory modes can be understood as a combination of a memory signal throughout the binary&#39;s inspiral and merger phases, and a quasinormal mode signal near the ringdown phase. Additionally, we find that the magnetic memory, up to numerical error, is indeed zero as previously conjectured. Lastly, we find that signal-to-noise ratios of memory for LIGO, the Einstein Telescope (ET), and the Laser Interferometer Space Antenna (LISA) with these new waveforms and new memory calculation are larger than previous expectations based on post-Newtonian or Minimal Waveform models.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Extending Gravitational Wave Extraction Using Weyl Characteristic Fields

We present a detailed methodology for extracting the full set of Newman-Penrose Weyl scalars from numerically generated spacetimes without requiring a tetrad that is completely orthonormal or perfectly aligned to the principal null directions. We also describe how to implement an extrapolation technique for computing the Weyl scalars&#39; contribution at asymptotic null infinity in postprocessing. These methods have been used to produce $Ψ_4$ and $h$ waveforms for the Simulating eXtreme Spacetimes (SXS) waveform catalog and now have been expanded to produce the entire set of Weyl scalars. These new waveform quantities are critical for the future of gravitational wave astronomy in order to understand the finite-amplitude gauge differences that can occur in numerical waveforms. We also present a new analysis of the accuracy of waveforms produced by the Spectral Einstein Code. While ultimately we expect Cauchy characteristic extraction to yield more accurate waveforms, the extraction techniques described here are far easier to implement and have already proven to be a viable way to produce production-level waveforms that can meet the demands of current gravitational-wave detectors.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Up-down instability of binary black holes in numerical relativity

Binary black holes with spins that are aligned with the orbital angular momentum do not precess. However, post-Newtonian calculations predict that &#34;up-down&#34; binaries, in which the spin of the heavier (lighter) black hole is aligned (antialigned) with the orbital angular momentum, are unstable when the spins are slightly perturbed from perfect alignment. This instability provides a possible mechanism for the formation of precessing binaries in environments where sources are preferentially formed with (anti) aligned spins. In this paper, we present the first full numerical relativity simulations capturing this instability. These simulations span $\sim 100$ orbits and $\sim 3$-$5$ precession cycles before merger, making them some of the longest numerical relativity simulations to date. Initialized with a small perturbation of $1^{\circ}$-$10^{\circ}$, the instability causes a dramatic growth of the spin misalignments, which can reach $\sim 90^{\circ}$ near merger. We show that this leaves a strong imprint on the subdominant modes of the gravitational wave signal, which can potentially be used to distinguish up-down binaries from other sources. Finally, we show that post-Newtonian and effective-one-body approximants are able to reproduce the unstable dynamics of up-down binaries extracted from numerical relativity.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Constraining the parameters of GW150914 & GW170104 with numerical relativity surrogates

Gravitational-wave detectors have begun to observe coalescences of heavy black holes at a consistent pace for the past few years. Accurate models of gravitational waveforms are essential for unbiased and precise estimation of source parameters. Recently developed surrogate models based on high-accuracy numerical relativity (NR) simulations are ideal for constraining physical parameters of heavy black hole merger events. In this paper, we first demonstrate the viability of these multi-modal surrogates as reliable parameter estimation tools. We show that NR surrogates can extract additional information from GW data that is inaccessible to traditional models, by analyzing a set of synthetic signals with the NR surrogate and other approximants. We also consider the case of two of the earliest binary black holes detected by the LIGO observatories: GW150914 and GW170104. We reanalyze their data with fully-precessing NR-surrogate templates and freely provide the resulting posterior samples as supplemental material. We find that our refined analysis is able to extract information from sub-dominant GW harmonics in data, and therefore better resolve the degeneracy in measuring source luminosity distance and orbital inclination for both events. We estimate the sources of both events to be 20-25% further away than was previously estimated. Our analyses also constrain their orbital orientation more tightly around face-on or face-off configurations than before. Additionally, for GW150914 we constrain the effective inspiral spin more tightly around zero. This work is one of the first to unambiguously extract sub-dominant GW mode information from real events. It is also a first step toward eliminating the approximations used in semi-analytic waveform models from GW parameter estimation. It also motivates that NR surrogates be extended to cover more of the binary black hole parameter space.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Distinguishing the nature of comparable-mass neutron star binary systems with multimessenger observations: GW170817 case study

The discovery of GW170817 with gravitational waves (GWs) and electromagnetic (EM) radiation is prompting new questions in strong-gravity astrophysics. Importantly, it remains unknown whether the progenitor of the merger comprised two neutron stars (NSs), or a NS and a black hole (BH). Using new numerical-relativity simulations and incorporating modeling uncertainties we produce novel GW and EM observables for NS-BH mergers with similar masses. A joint analysis of GW and EM measurements reveals that if GW170817 is a NS-BH merger, <40% of the binary parameters consistent with the GW data are compatible with EM observations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

High-accuracy waveforms for black hole-neutron star systems with spinning black holes

The availability of accurate numerical waveforms is an important requirement for the creation and calibration of reliable waveform models for gravitational wave astrophysics. For black hole-neutron star binaries, very few accurate waveforms are however publicly available. Most recent models are calibrated to a large number of older simulations with good parameter space coverage for low-spin non-precessing binaries but limited accuracy, and a much smaller number of longer, more recent simulations limited to non-spinning black holes. In this paper, we present long, accurate numerical waveforms for three new systems that include rapidly spinning black holes, and one precessing configuration. We study in detail the accuracy of the simulations, and in particular perform for the first time in the context of BHNS binaries a detailed comparison of waveform extrapolation methods to the results of Cauchy Characteristic Extraction. The new waveforms have $<0.1\,{\rm rad}$ phase errors during inspiral, rising to $\sim (0.2-0.4)\,{\rm rad}$ errors at merger, and $\lesssim 1\%$ error in their amplitude. We compute the faithfulness of recent analytical models to these numerical results, and find that models specifically designed for BHNS binaries perform well ($F>0.99$) for binaries seen face-on. For edge-on observations, particularly for precessing systems, disagreements between models and simulations increase, and models that include precession and/or higher-order modes start to perform better than BHNS models that currently lack these features.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Impact of subdominant modes on the interpretation of gravitational-wave signals from heavy binary black hole systems

Over the past year, a handful of new gravitational wave models have been developed to include multiple harmonic modes thereby enabling for the first time fully Bayesian inference studies including higher modes to be performed. Using one recently-developed numerical relativity surrogate model, NRHybSur3dq8, we investigate the importance of higher modes on parameter inference of coalescing massive binary black holes. We focus on examples relevant to the current three-detector network of observatories, with a detector-frame mass set to $120 M_\odot$ and with signal amplitude values that are consistent with plausible candidates for the next few observing runs. We show that for such systems the higher mode content will be important for interpreting coalescing binary black holes, reducing systematic bias, and computing properties of the remnant object. Even for comparable-mass binaries and at low signal amplitude, the omission of higher modes can influence posterior probability distributions. We discuss the impact of our results on source population inference and self-consistency tests of general relativity. Our work can be used to better understand asymmetric binary black hole merger events, such as GW190412. Higher modes are critical for such systems, and their omission usually produces substantial parameter biases.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Multipolar Effective-One-Body Waveforms for Precessing Binary Black Holes: Construction and Validation

As gravitational-wave detectors become more sensitive, we will access a greater variety of signals emitted by compact binary systems, shedding light on their astrophysical origin and environment. A key physical effect that can distinguish among formation scenarios is the misalignment of the spins with the orbital angular momentum, causing the spins and the binary&#39;s orbital plane to precess. To accurately model such systems, it is crucial to include multipoles beyond the dominant quadrupole. Here, we develop the first multipolar precessing waveform model in the effective-one-body (EOB) formalism for the inspiral, merger and ringdown (IMR) of binary black holes: SEOBNRv4PHM. In the nonprecessing limit, the model reduces to SEOBNRv4HM, which was calibrated to numerical-relativity (NR) simulations, and waveforms from perturbation theory. We validate SEOBNRv4PHM by comparing it to the public catalog of 1405 precessing NR waveforms of the Simulating eXtreme Spacetimes (SXS) collaboration, and also to new 118 precessing NR waveforms, which span mass ratios 1-4 and spins up to 0.9. We stress that SEOBNRv4PHM is not calibrated to NR simulations in the precessing sector. We compute the unfaithfulness against the 1523 SXS precessing NR waveforms, and find that, for $94\%$ ($57\%$) of the cases, the maximum value, in the total mass range $20-200 M_\odot$, is below $3\%$ ($1\%$). Those numbers become $83\%$ ($20\%$) when using the IMR, multipolar, precessing phenomenological model IMRPhenomPv3HM. We investigate the impact of such unfaithfulness values with two parameter-estimation studies on synthetic signals. We also compute the unfaithfulness between those waveform models and identify in which part of the parameter space they differ the most. We validate them also against the multipolar, precessing NR surrogate model NRSur7dq4, and find that the SEOBNRv4PHM model outperforms IMRPhenomPv3HM.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Unequal Mass Binary Neutron Star Simulations with Neutrino Transport: Ejecta and Neutrino Emission

We present twelve new simulations of unequal mass neutron star mergers. The simulations were preformed with the SpEC code, and utilize nuclear-theory based equations of state and a two-moment gray neutrino transport scheme with an improved energy estimate based on evolving the number density. We model the neutron stars with the SFHo, LS220 and DD2 equations of state (EOS) and we study the neutrino and matter emission of all twelve models to search for robust trends between binary parameters and emission characteristics. We find that the total mass of the dynamical ejecta exceeds $0.01M_\odot$ only for SFHo with weak dependence on the mass-ratio across all models. We find that the ejecta have a broad electron fraction ($Y_e$) distribution ($\approx 0.06-0.48$), with mean $0.2$. $Y_e$ increases with neutrino irradiation over time, but decreases with increasing binary asymmetry. We also find that the models have ejecta with a broad asymptotic velocity distribution ($\approx 0.05-0.7c$). The average velocity lies in the range $0.2c - 0.3c$ and decreases with binary asymmetry. Furthermore, we find that disk mass increases with binary asymmetry and stiffness of the EOS. The $Y_e$ of the disk increases with softness of the EOS. The strongest neutrino emission occurs for the models with soft EOS. For (anti) electron neutrinos we find no significant dependence of the magnitude or angular distribution or neutrino luminosity with mass-ratio. The heavier neutrino species have a luminosity dependence on mass-ratio but an angular distribution which does not change with mass-ratio.

0 citations0 reviewsNo code linkedOpen access