-
Waveform Modelling for the Laser Interferometer Space Antenna
Authors:
LISA Consortium Waveform Working Group,
Niayesh Afshordi,
Sarp Akçay,
Pau Amaro Seoane,
Andrea Antonelli,
Josu C. Aurrekoetxea,
Leor Barack,
Enrico Barausse,
Robert Benkel,
Laura Bernard,
Sebastiano Bernuzzi,
Emanuele Berti,
Matteo Bonetti,
Béatrice Bonga,
Gabriele Bozzola,
Richard Brito,
Alessandra Buonanno,
Alejandro Cárdenas-Avendaño,
Marc Casals,
David F. Chernoff,
Alvin J. K. Chua,
Katy Clough,
Marta Colleoni,
Mekhi Dhesi,
Adrien Druart
, et al. (121 additional authors not shown)
Abstract:
LISA, the Laser Interferometer Space Antenna, will usher in a new era in gravitational-wave astronomy. As the first anticipated space-based gravitational-wave detector, it will expand our view to the millihertz gravitational-wave sky, where a spectacular variety of interesting new sources abound: from millions of ultra-compact binaries in our Galaxy, to mergers of massive black holes at cosmologic…
▽ More
LISA, the Laser Interferometer Space Antenna, will usher in a new era in gravitational-wave astronomy. As the first anticipated space-based gravitational-wave detector, it will expand our view to the millihertz gravitational-wave sky, where a spectacular variety of interesting new sources abound: from millions of ultra-compact binaries in our Galaxy, to mergers of massive black holes at cosmological distances; from the beginnings of inspirals that will venture into the ground-based detectors' view to the death spiral of compact objects into massive black holes, and many sources in between. Central to realising LISA's discovery potential are waveform models, the theoretical and phenomenological predictions of the pattern of gravitational waves that these sources emit. This white paper is presented on behalf of the Waveform Working Group for the LISA Consortium. It provides a review of the current state of waveform models for LISA sources, and describes the significant challenges that must yet be overcome.
△ Less
Submitted 20 December, 2023; v1 submitted 2 November, 2023;
originally announced November 2023.
-
Towards exponentially-convergent simulations of extreme-mass-ratio inspirals: A time-domain solver for the scalar Teukolsky equation with singular source terms
Authors:
Manas Vishal,
Scott E. Field,
Katie Rink,
Sigal Gottlieb,
Gaurav Khanna
Abstract:
Gravitational wave signals from extreme mass ratio inspirals are a key target for space-based gravitational wave detectors. These systems are typically modeled as a distributionally-forced Teukolsky equation, where the smaller black hole is treated as a Dirac delta distribution. Time-domain solvers often use regularization approaches that approximate the Dirac distribution that often introduce sma…
▽ More
Gravitational wave signals from extreme mass ratio inspirals are a key target for space-based gravitational wave detectors. These systems are typically modeled as a distributionally-forced Teukolsky equation, where the smaller black hole is treated as a Dirac delta distribution. Time-domain solvers often use regularization approaches that approximate the Dirac distribution that often introduce small length scales and are a source of systematic error, especially near the smaller black hole. We describe a multi-domain discontinuous Galerkin method for solving the distributionally-forced Teukolsky equation that describes scalar fields evolving on a Kerr spacetime. To handle the Dirac delta, we expand the solution in spherical harmonics and recast the sourced Teukolsky equation as a first-order, one-dimensional symmetric hyperbolic system. This allows us to derive the method's numerical flux to correctly account for the Dirac delta. As a result, our method achieves global spectral accuracy even at the source's location. To connect the near field to future null infinity, we use the hyperboloidal layer method, allowing us to supply outer boundary conditions and providing direct access to the far-field waveform. We document several numerical experiments where we test our method, including convergence tests against exact solutions, energy luminosities for circular orbits, the scheme's superconvergence properties at future null infinity, and the late-time tail behavior of the scalar field. We also compare two systems that arise from different choices of the first-order reduction variables, finding that certain choices are numerically problematic in practice. The methods developed here may be beneficial when computing gravitational self-force effects, where the regularization procedure has been developed for the spherical harmonic modes and high accuracy is needed at the Dirac delta's location.
△ Less
Submitted 6 December, 2023; v1 submitted 3 July, 2023;
originally announced July 2023.
-
Remnant black hole properties from numerical-relativity-informed perturbation theory and implications for waveform modelling
Authors:
Tousif Islam,
Scott E. Field,
Gaurav Khanna
Abstract:
During binary black hole (BBH) mergers, energy and momenta are carried away from the binary system as gravitational radiation. Access to the radiated energy and momenta allows us to predict the properties of the remnant black hole. We develop a Python package gw_remnant to compute the remnant mass, remnant spin, peak luminosity, and the final kick imparted on the remnant black hole from the gravit…
▽ More
During binary black hole (BBH) mergers, energy and momenta are carried away from the binary system as gravitational radiation. Access to the radiated energy and momenta allows us to predict the properties of the remnant black hole. We develop a Python package gw_remnant to compute the remnant mass, remnant spin, peak luminosity, and the final kick imparted on the remnant black hole from the gravitational radiation. Using this package, we compute the remnant properties of the final black hole in case of non-spinning BBH mergers with mass ratios ranging from q=2.5 to q=1000 using waveform modes generated from BHPTNRSur1dq1e4, a recently developed numerical-relativity-informed surrogate model based on black-hole perturbation theory framework. We validate our results against the remnant properties estimated from numerical relativity (NR) surrogate models in the comparable mass ratio regime and against recently available high-mass-ratio NR simulations at q=[15,32,64]. We find that our remnant property estimates computed from fluxes at future null infinity closely match the estimates obtained from the NR surrogate model of apparent horizon data. Using Gaussian process regression fitting methods, we train a surrogate model, {BHPTNR_Remnant, for the properties of the remnant black hole arising from BBH mergers with mass ratios from q=2.5 to q=1000. Finally, we discuss potential improvements in the BHPTNRSur1dq1e4 waveform model when including remnant information. We make both the gw_remnant and BHPTNR_Remnant packages publicly available.
△ Less
Submitted 12 February, 2023; v1 submitted 17 January, 2023;
originally announced January 2023.
-
Evidence of large recoil velocity from a black hole merger signal
Authors:
Vijay Varma,
Sylvia Biscoveanu,
Tousif Islam,
Feroz H. Shaik,
Carl-Johan Haster,
Maximiliano Isi,
Will M. Farr,
Scott E. Field,
Salvatore Vitale
Abstract:
The final black hole left behind after a binary black hole merger can attain a recoil velocity, or a "kick", reaching values up to 5000 km/s. This phenomenon has important implications for gravitational wave astronomy, black hole formation scenarios, testing general relativity, and galaxy evolution. We consider the gravitational wave signal from the binary black hole merger GW200129_065458 (hencef…
▽ More
The final black hole left behind after a binary black hole merger can attain a recoil velocity, or a "kick", reaching values up to 5000 km/s. This phenomenon has important implications for gravitational wave astronomy, black hole formation scenarios, testing general relativity, and galaxy evolution. We consider the gravitational wave signal from the binary black hole merger GW200129_065458 (henceforth referred to as GW200129), which has been shown to exhibit strong evidence of orbital precession. Using numerical relativity surrogate models, we constrain the kick velocity of GW200129 to $v_f \sim 1542^{+747}_{-1098}$ km/s or $v_f \gtrsim 698$ km/s (one-sided limit), at 90\% credibility. This marks the first identification of a large kick velocity for an individual gravitational wave event. Given the kick velocity of GW200129, we estimate that there is a less than $0.48\%$ ($7.7\%$) probability that the remnant black hole after the merger would be retained by globular (nuclear star) clusters. Finally, we show that kick effects are not expected to cause biases in ringdown tests of general relativity for this event, although this may change in the future with improved detectors.
△ Less
Submitted 29 April, 2022; v1 submitted 4 January, 2022;
originally announced January 2022.
-
High-precision source characterization of intermediate mass-ratio black hole coalescences with gravitational waves: The importance of higher-order multipoles
Authors:
Tousif Islam,
Scott E. Field,
Carl-Johan Haster,
Rory Smith
Abstract:
Intermediate mass ratio inspiral (IMRI) binaries -- containing stellar-mass black holes coalescing into intermediate-mass black holes ($M>100M_{\odot}$) -- are a highly anticipated source of gravitational waves (GWs) for Advanced LIGO/Virgo. Their detection and source characterization would provide a unique probe of strong-field gravity and stellar evolution. Due to the asymmetric component masses…
▽ More
Intermediate mass ratio inspiral (IMRI) binaries -- containing stellar-mass black holes coalescing into intermediate-mass black holes ($M>100M_{\odot}$) -- are a highly anticipated source of gravitational waves (GWs) for Advanced LIGO/Virgo. Their detection and source characterization would provide a unique probe of strong-field gravity and stellar evolution. Due to the asymmetric component masses and the large primary, these systems generically excite subdominant modes while reducing the importance of the dominant quadrupole mode. Including higher order harmonics can also result in a $10\%-25\%$ increase in signal-to-noise ratio for IMRIs, which may help to detect these systems. We show that by including subdominant GW modes into the analysis we can achieve a precise characterization of IMRI source properties. For example, we find that the source properties for IMRIs can be measured to within $2\%-15\%$ accuracy at a fiducial signal-to-noise ratio of 25 if subdominant modes are included. When subdominant modes are neglected, the accuracy degrades to $9\%-44\%$ and significant biases are seen in chirp mass, mass ratio, primary spin and luminosity distances. We further demonstrate that including subdominant modes in the waveform model can enable an informative measurement of both individual spin components and improve the source localization by a factor of $\sim$10. We discuss some important astrophysical implications of high-precision source characterization enabled by subdominant modes such as constraining the mass gap and probing formation channels.
△ Less
Submitted 19 October, 2021; v1 submitted 10 May, 2021;
originally announced May 2021.
-
Eccentric binary black hole surrogate models for the gravitational waveform and remnant properties: comparable mass, nonspinning case
Authors:
Tousif Islam,
Vijay Varma,
Jackie Lodman,
Scott E. Field,
Gaurav Khanna,
Mark A. Scheel,
Harald P. Pfeiffer,
Davide Gerosa,
Lawrence E. Kidder
Abstract:
We develop new strategies to build numerical relativity surrogate models for eccentric binary black hole systems, which are expected to play an increasingly important role in current and future gravitational-wave detectors. We introduce a new surrogate waveform model, \texttt{NRSur2dq1Ecc}, using 47 nonspinning, equal-mass waveforms with eccentricities up to $0.2$ when measured at a reference time…
▽ More
We develop new strategies to build numerical relativity surrogate models for eccentric binary black hole systems, which are expected to play an increasingly important role in current and future gravitational-wave detectors. We introduce a new surrogate waveform model, \texttt{NRSur2dq1Ecc}, using 47 nonspinning, equal-mass waveforms with eccentricities up to $0.2$ when measured at a reference time of $5500M$ before merger. This is the first waveform model that is directly trained on eccentric numerical relativity simulations and does not require that the binary circularizes before merger. The model includes the $(2,2)$, $(3,2)$, and $(4,4)$ spin-weighted spherical harmonic modes. We also build a final black hole model, \texttt{NRSur2dq1EccRemnant}, which models the mass, and spin of the remnant black hole. We show that our waveform model can accurately predict numerical relativity waveforms with mismatches $\approx 10^{-3}$, while the remnant model can recover the final mass and dimensionless spin with absolute errors smaller than $\approx 5 \times 10^{-4}M$ and $\approx 2 \times10^{-3}$ respectively. We demonstrate that the waveform model can also recover subtle effects like mode-mixing in the ringdown signal without any special ad-hoc modeling steps. Finally, we show that despite being trained only on equal-mass binaries, \texttt{NRSur2dq1Ecc} can be reasonably extended up to mass ratio $q\approx3$ with mismatches $\simeq 10^{-2}$ for eccentricities smaller than $\sim 0.05$ as measured at a reference time of $2000M$ before merger. The methods developed here should prove useful in the building of future eccentric surrogate models over larger regions of the parameter space.
△ Less
Submitted 24 March, 2021; v1 submitted 27 January, 2021;
originally announced January 2021.
-
Improved analysis of GW190412 with a precessing numerical relativity surrogate waveform model
Authors:
Tousif Islam,
Scott E. Field,
Carl-Johan Haster,
Rory Smith
Abstract:
The recent observation of GW190412, the first high-mass ratio binary black-hole (BBH) merger, by the LIGO-Virgo Collaboration (LVC) provides a unique opportunity to probe the impact of subdominant harmonics and precession effects encoded in a gravitational wave signal. We present refined estimates of source parameters for GW190412 using \texttt{NRSur7dq4}, a recently developed numerical relativity…
▽ More
The recent observation of GW190412, the first high-mass ratio binary black-hole (BBH) merger, by the LIGO-Virgo Collaboration (LVC) provides a unique opportunity to probe the impact of subdominant harmonics and precession effects encoded in a gravitational wave signal. We present refined estimates of source parameters for GW190412 using \texttt{NRSur7dq4}, a recently developed numerical relativity waveform surrogate model that includes all $\ell \leq 4$ spin-weighted spherical harmonic modes as well as the full physical effects of precession. We compare our results with two different variants of phenomenological precessing BBH waveform models, \texttt{IMRPhenomPv3HM} and \texttt{IMRPhenomXPHM}, as well as to the LVC results. Our results are broadly in agreement with \texttt{IMRPhenomXPHM} results and the reported LVC analysis compiled with the \texttt{SEOBNRv4PHM} waveform model, but in tension with \texttt{IMRPhenomPv3HM}. Using the \texttt{NRSur7dq4} model, we provide a tighter constraint on the mass-ratio ($0.26^{+0.08}_{-0.06}$) as compared to the LVC estimate of $0.28^{+0.13}_{-0.07}$ (both reported as median values withs 90\% credible intervals). We also constrain the binary to be more face-on, and find a broader posterior for the spin precession parameter. We further find that even though $\ell=4$ harmonic modes have negligible signal-to-noise ratio, omission of these modes will influence the estimated posterior distribution of several source parameters including chirp mass, effective inspiral spin, luminosity distance, and inclination. We also find that commonly used model approximations, such as neglecting the asymmetric modes (which are generically excited during precession), have negligible impact on parameter recovery for moderate SNR-events similar to GW190412.
△ Less
Submitted 20 April, 2021; v1 submitted 9 October, 2020;
originally announced October 2020.
-
A Surrogate Model for Gravitational Wave Signals from Comparable- to Large- Mass-Ratio Black Hole Binaries
Authors:
Nur E. M. Rifat,
Scott E. Field,
Gaurav Khanna,
Vijay Varma
Abstract:
Gravitational wave signals from compact astrophysical sources such as those observed by LIGO and Virgo require a high-accuracy, theory-based waveform model for the analysis of the recorded signal. Current inspiral-merger-ringdown models are calibrated only up to moderate mass ratios, thereby limiting their applicability to signals from high-mass ratio binary systems. We present EMRISur1dq1e4, a re…
▽ More
Gravitational wave signals from compact astrophysical sources such as those observed by LIGO and Virgo require a high-accuracy, theory-based waveform model for the analysis of the recorded signal. Current inspiral-merger-ringdown models are calibrated only up to moderate mass ratios, thereby limiting their applicability to signals from high-mass ratio binary systems. We present EMRISur1dq1e4, a reduced-order surrogate model for gravitational waveforms of 13,500M in duration and including several harmonic modes for non-spinning black hole binary systems with mass-ratios varying from 3 to 10,000 thus vastly expanding the parameter range beyond the current models. This surrogate model is trained on waveform data generated by point-particle black hole perturbation theory (ppBHPT) both for large mass-ratio and comparable mass-ratio binaries. We observe that the gravitational waveforms generated through a simple application of ppBHPT to the comparable mass-ratio cases agree remarkably (and surprisingly) well with those from full numerical relativity after a rescaling of the ppBHPT's total mass parameter. This observation and the EMRISur1dq1e4 surrogate model will enable data analysis studies in the high-mass ratio regime, including potential intermediate mass-ratio signals from LIGO/Virgo and extreme-mass ratio events of interest to the future space-based observatory LISA.
△ Less
Submitted 1 April, 2020; v1 submitted 23 October, 2019;
originally announced October 2019.
-
Surrogate models for precessing binary black hole simulations with unequal masses
Authors:
Vijay Varma,
Scott E. Field,
Mark A. Scheel,
Jonathan Blackman,
Davide Gerosa,
Leo C. Stein,
Lawrence E. Kidder,
Harald P. Pfeiffer
Abstract:
Only numerical relativity simulations can capture the full complexities of binary black hole mergers. These simulations, however, are prohibitively expensive for direct data analysis applications such as parameter estimation. We present two new fast and accurate surrogate models for the outputs of these simulations: the first model, NRSur7dq4, predicts the gravitational waveform and the second mod…
▽ More
Only numerical relativity simulations can capture the full complexities of binary black hole mergers. These simulations, however, are prohibitively expensive for direct data analysis applications such as parameter estimation. We present two new fast and accurate surrogate models for the outputs of these simulations: the first model, NRSur7dq4, predicts the gravitational waveform and the second model, \RemnantModel, predicts the properties of the remnant black hole. These models extend previous 7-dimensional, non-eccentric precessing models to higher mass ratios, and have been trained against 1528 simulations with mass ratios $q\leq4$ and spin magnitudes $χ_1,χ_2 \leq 0.8$, with generic spin directions. The waveform model, NRSur7dq4, which begins about 20 orbits before merger, includes all $\ell \leq 4$ spin-weighted spherical harmonic modes, as well as the precession frame dynamics and spin evolution of the black holes. The final black hole model, \RemnantModel, models the mass, spin, and recoil kick velocity of the remnant black hole. In their training parameter range, both models are shown to be more accurate than existing models by at least an order of magnitude, with errors comparable to the estimated errors in the numerical relativity simulations. We also show that the surrogate models work well even when extrapolated outside their training parameter space range, up to mass ratios $q=6$.
△ Less
Submitted 10 October, 2019; v1 submitted 22 May, 2019;
originally announced May 2019.
-
The SXS Collaboration catalog of binary black hole simulations
Authors:
Michael Boyle,
Daniel Hemberger,
Dante A. B. Iozzo,
Geoffrey Lovelace,
Serguei Ossokine,
Harald P. Pfeiffer,
Mark A. Scheel,
Leo C. Stein,
Charles J. Woodford,
Aaron B. Zimmerman,
Nousha Afshari,
Kevin Barkett,
Jonathan Blackman,
Katerina Chatziioannou,
Tony Chu,
Nicholas Demos,
Nils Deppe,
Scott E. Field,
Nils L. Fischer,
Evan Foley,
Heather Fong,
Alyssa Garcia,
Matthew Giesler,
Francois Hebert,
Ian Hinder
, et al. (15 additional authors not shown)
Abstract:
Accurate models of gravitational waves from merging black holes are necessary for detectors to observe as many events as possible while extracting the maximum science. Near the time of merger, the gravitational waves from merging black holes can be computed only using numerical relativity. In this paper, we present a major update of the Simulating eXtreme Spacetimes (SXS) Collaboration catalog of…
▽ More
Accurate models of gravitational waves from merging black holes are necessary for detectors to observe as many events as possible while extracting the maximum science. Near the time of merger, the gravitational waves from merging black holes can be computed only using numerical relativity. In this paper, we present a major update of the Simulating eXtreme Spacetimes (SXS) Collaboration catalog of numerical simulations for merging black holes. The catalog contains 2018 distinct configurations (a factor of 11 increase compared to the 2013 SXS catalog), including 1426 spin-precessing configurations, with mass ratios between 1 and 10, and spin magnitudes up to 0.998. The median length of a waveform in the catalog is 39 cycles of the dominant $\ell=m=2$ gravitational-wave mode, with the shortest waveform containing 7.0 cycles and the longest 351.3 cycles. We discuss improvements such as correcting for moving centers of mass and extended coverage of the parameter space. We also present a thorough analysis of numerical errors, finding typical truncation errors corresponding to a waveform mismatch of $\sim 10^{-4}$. The simulations provide remnant masses and spins with uncertainties of 0.03% and 0.1% ($90^{\text{th}}$ percentile), about an order of magnitude better than analytical models for remnant properties. The full catalog is publicly available at https://www.black-holes.org/waveforms .
△ Less
Submitted 11 September, 2019; v1 submitted 9 April, 2019;
originally announced April 2019.
-
SpECTRE: A Task-based Discontinuous Galerkin Code for Relativistic Astrophysics
Authors:
Lawrence E. Kidder,
Scott E. Field,
Francois Foucart,
Erik Schnetter,
Saul A. Teukolsky,
Andy Bohn,
Nils Deppe,
Peter Diener,
François Hébert,
Jonas Lippuner,
Jonah Miller,
Christian D. Ott,
Mark A. Scheel,
Trevor Vincent
Abstract:
We introduce a new relativistic astrophysics code, SpECTRE, that combines a discontinuous Galerkin method with a task-based parallelism model. SpECTRE's goal is to achieve more accurate solutions for challenging relativistic astrophysics problems such as core-collapse supernovae and binary neutron star mergers. The robustness of the discontinuous Galerkin method allows for the use of high-resoluti…
▽ More
We introduce a new relativistic astrophysics code, SpECTRE, that combines a discontinuous Galerkin method with a task-based parallelism model. SpECTRE's goal is to achieve more accurate solutions for challenging relativistic astrophysics problems such as core-collapse supernovae and binary neutron star mergers. The robustness of the discontinuous Galerkin method allows for the use of high-resolution shock capturing methods in regions where (relativistic) shocks are found, while exploiting high-order accuracy in smooth regions. A task-based parallelism model allows efficient use of the largest supercomputers for problems with a heterogeneous workload over disparate spatial and temporal scales. We argue that the locality and algorithmic structure of discontinuous Galerkin methods will exhibit good scalability within a task-based parallelism framework. We demonstrate the code on a wide variety of challenging benchmark problems in (non)-relativistic (magneto)-hydrodynamics. We demonstrate the code's scalability including its strong scaling on the NCSA Blue Waters supercomputer up to the machine's full capacity of 22,380 nodes using 671,400 threads.
△ Less
Submitted 21 July, 2017; v1 submitted 31 August, 2016;
originally announced September 2016.
-
Fast and Accurate Inference on Gravitational Waves from Precessing Compact Binaries
Authors:
Rory Smith,
Scott E. Field,
Kent Blackburn,
Carl-Johan Haster,
Michael Pürrer,
Vivien Raymond,
Patricia Schmidt
Abstract:
Inferring astrophysical information from gravitational waves emitted by compact binaries is one of the key science goals of gravitational-wave astronomy. In order to reach the full scientific potential of gravitational-wave experiments we require techniques to mitigate the cost of Bayesian inference, especially as gravitational-wave signal models and analyses become increasingly sophisticated and…
▽ More
Inferring astrophysical information from gravitational waves emitted by compact binaries is one of the key science goals of gravitational-wave astronomy. In order to reach the full scientific potential of gravitational-wave experiments we require techniques to mitigate the cost of Bayesian inference, especially as gravitational-wave signal models and analyses become increasingly sophisticated and detailed. Reduced order models (ROMs) of gravitational waveforms can significantly reduce the computational cost of inference by removing redundant computations. In this paper we construct the first reduced order models of gravitational-wave signals that include the effects of spin-precession, inspiral, merger, and ringdown in compact object binaries, and which are valid for component masses describing binary neutron star, binary black hole and mixed binary systems. This work utilizes the waveform model known as "IMRPhenomPv2". Our ROM enables the use of a fast reduced order quadrature (ROQ) integration rule which allows us to approximate Bayesian probability density functions at a greatly reduced computational cost. We find that the ROQ rule can be used to speed up inference by factors as high as 300 without introducing systematic bias. This corresponds to a reduction in computational time from around half a year to a half a day, for the longest duration/lowest mass signals. The ROM and ROQ rule are available with the main inference library of the LIGO Scientific Collaboration, LALInference.
△ Less
Submitted 29 August, 2016; v1 submitted 27 April, 2016;
originally announced April 2016.
-
Fast and accurate prediction of numerical relativity waveforms from binary black hole coalescences using surrogate models
Authors:
Jonathan Blackman,
Scott E. Field,
Chad R. Galley,
Bela Szilagyi,
Mark A. Scheel,
Manuel Tiglio,
Daniel A. Hemberger
Abstract:
Simulating a binary black hole (BBH) coalescence by solving Einstein's equations is computationally expensive, requiring days to months of supercomputing time. Using reduced order modeling techniques, we construct an accurate surrogate model, which is evaluated in a millisecond to a second, for numerical relativity (NR) waveforms from non-spinning BBH coalescences with mass ratios in $[1, 10]$ and…
▽ More
Simulating a binary black hole (BBH) coalescence by solving Einstein's equations is computationally expensive, requiring days to months of supercomputing time. Using reduced order modeling techniques, we construct an accurate surrogate model, which is evaluated in a millisecond to a second, for numerical relativity (NR) waveforms from non-spinning BBH coalescences with mass ratios in $[1, 10]$ and durations corresponding to about $15$ orbits before merger. We assess the model's uncertainty and show that our modeling strategy predicts NR waveforms {\em not} used for the surrogate's training with errors nearly as small as the numerical error of the NR code. Our model includes all spherical-harmonic ${}_{-2}Y_{\ell m}$ waveform modes resolved by the NR code up to $\ell=8.$ We compare our surrogate model to Effective One Body waveforms from $50$-$300 M_\odot$ for advanced LIGO detectors and find that the surrogate is always more faithful (by at least an order of magnitude in most cases).
△ Less
Submitted 7 October, 2015; v1 submitted 26 February, 2015;
originally announced February 2015.
-
Gravitational wave parameter estimation with compressed likelihood evaluations
Authors:
Priscilla Canizares,
Scott E. Field,
Jonathan R. Gair,
Manuel Tiglio
Abstract:
One of the main bottlenecks in gravitational wave (GW) astronomy is the high cost of performing parameter estimation and GW searches on the fly. We propose a novel technique based on Reduced Order Quadratures (ROQs), an application and data-specific quadrature rule, to perform fast and accurate likelihood evaluations. These are the dominant cost in Markov chain Monte Carlo (MCMC) algorithms, which…
▽ More
One of the main bottlenecks in gravitational wave (GW) astronomy is the high cost of performing parameter estimation and GW searches on the fly. We propose a novel technique based on Reduced Order Quadratures (ROQs), an application and data-specific quadrature rule, to perform fast and accurate likelihood evaluations. These are the dominant cost in Markov chain Monte Carlo (MCMC) algorithms, which are widely employed in parameter estimation studies, and so ROQs offer a new way to accelerate GW parameter estimation. We illustrate our approach using a four dimensional GW burst model embedded in noise. We build an ROQ for this model, and perform four dimensional MCMC searches with both the standard and ROQs quadrature rules, showing that, for this model, the ROQ approach is around 25 times faster than the standard approach with essentially no loss of accuracy. The speed-up from using ROQs is expected to increase for more complex GW signal models and therefore has significant potential to accelerate parameter estimation of GW sources such as compact binary coalescences.
△ Less
Submitted 25 September, 2013; v1 submitted 1 April, 2013;
originally announced April 2013.
-
Reduced Basis representations of multi-mode black hole ringdown gravitational waves
Authors:
Sarah Caudill,
Scott E. Field,
Chad R. Galley,
Frank Herrmann,
Manuel Tiglio
Abstract:
We construct compact and high accuracy Reduced Basis (RB) representations of single and multiple quasinormal modes (QNMs). The RB method determines a hierarchical and relatively small set of the most relevant waveforms. We find that the exponential convergence of the method allows for a dramatic compression of template banks used for ringdown searches. Compressing a catalog with a minimal match…
▽ More
We construct compact and high accuracy Reduced Basis (RB) representations of single and multiple quasinormal modes (QNMs). The RB method determines a hierarchical and relatively small set of the most relevant waveforms. We find that the exponential convergence of the method allows for a dramatic compression of template banks used for ringdown searches. Compressing a catalog with a minimal match $\MMm=0.99$, we find that the selected RB waveforms are able to represent {\em any} QNM, including those not in the original bank, with extremely high accuracy, typically less than $10^{-13}$. We then extend our studies to two-mode QNMs. Inclusion of a second mode is expected to help with detection, and might make it possible to infer details of the progenitor of the final black hole. We find that the number of RB waveforms needed to represent any two-mode ringdown waveform with the above high accuracy is {\em smaller} than the number of metric-based, one-mode templates with $\MMm=0.99$. For unconstrained two-modes, which would allow for consistency tests of General Relativity, our high accuracy RB has around $10^4$ {\em fewer} waveforms than the number of metric-based templates for $\MMm=0.99$. The number of RB elements grows only linearly with the number of multipole modes versus exponentially with the standard approach, resulting in very compact representations even for many multiple modes. The results of this paper open the possibility of searches of multi-mode ringdown gravitational waves.
△ Less
Submitted 18 March, 2012; v1 submitted 26 September, 2011;
originally announced September 2011.
-
Reduced basis catalogs for gravitational wave templates
Authors:
Scott E. Field,
Chad R. Galley,
Frank Herrmann,
Jan S. Hesthaven,
Evan Ochsner,
Manuel Tiglio
Abstract:
We introduce a reduced basis approach as a new paradigm for modeling, representing and searching for gravitational waves. We construct waveform catalogs for non-spinning compact binary coalescences, and we find that for accuracies of 99% and 99.999% the method generates a factor of about $10-10^5$ fewer templates than standard placement methods. The continuum of gravitational waves can be represen…
▽ More
We introduce a reduced basis approach as a new paradigm for modeling, representing and searching for gravitational waves. We construct waveform catalogs for non-spinning compact binary coalescences, and we find that for accuracies of 99% and 99.999% the method generates a factor of about $10-10^5$ fewer templates than standard placement methods. The continuum of gravitational waves can be represented by a finite and comparatively compact basis. The method is robust under variations in the noise of detectors, implying that only a single catalog needs to be generated.
△ Less
Submitted 13 June, 2011; v1 submitted 19 January, 2011;
originally announced January 2011.