-
Real-time gravitational-wave inference for binary neutron stars using machine learning
Authors:
Maximilian Dax,
Stephen R. Green,
Jonathan Gair,
Nihar Gupte,
Michael Pürrer,
Vivien Raymond,
Jonas Wildberger,
Jakob H. Macke,
Alessandra Buonanno,
Bernhard Schölkopf
Abstract:
Mergers of binary neutron stars (BNSs) emit signals in both the gravitational-wave (GW) and electromagnetic (EM) spectra. Famously, the 2017 multi-messenger observation of GW170817 led to scientific discoveries across cosmology, nuclear physics, and gravity. Central to these results were the sky localization and distance obtained from GW data, which, in the case of GW170817, helped to identify the…
▽ More
Mergers of binary neutron stars (BNSs) emit signals in both the gravitational-wave (GW) and electromagnetic (EM) spectra. Famously, the 2017 multi-messenger observation of GW170817 led to scientific discoveries across cosmology, nuclear physics, and gravity. Central to these results were the sky localization and distance obtained from GW data, which, in the case of GW170817, helped to identify the associated EM transient, AT 2017gfo, 11 hours after the GW signal. Fast analysis of GW data is critical for directing time-sensitive EM observations; however, due to challenges arising from the length and complexity of signals, it is often necessary to make approximations that sacrifice accuracy. Here, we develop a machine learning approach that performs complete BNS inference in just one second without making any such approximations. This is enabled by a new method for explicit integration of physical domain knowledge into neural networks. Our approach enhances multi-messenger observations by providing (i) accurate localization even before the merger; (ii) improved localization precision by $\sim30\%$ compared to approximate low-latency methods; and (iii) detailed information on luminosity distance, inclination, and masses, which can be used to prioritize expensive telescope time. Additionally, the flexibility and reduced cost of our method open new opportunities for equation-of-state and waveform systematics studies. Finally, we demonstrate that our method scales to extremely long signals, up to an hour in length, thus serving as a blueprint for data analysis for next-generation ground- and space-based detectors.
△ Less
Submitted 12 July, 2024;
originally announced July 2024.
-
Evidence for eccentricity in the population of binary black holes observed by LIGO-Virgo-KAGRA
Authors:
Nihar Gupte,
Antoni Ramos-Buades,
Alessandra Buonanno,
Jonathan Gair,
M. Coleman Miller,
Maximilian Dax,
Stephen R. Green,
Michael Pürrer,
Jonas Wildberger,
Jakob Macke,
Bernhard Schölkopf
Abstract:
Binary black holes (BBHs) in eccentric orbits produce distinct modulations the emitted gravitational waves (GWs). The measurement of orbital eccentricity can provide robust evidence for dynamical binary formation channels. We analyze 57 GW events from the first, second and third observing runs of the LIGO-Virgo-KAGRA (LVK) Collaboration using a multipolar aligned-spin inspiral-merger-ringdown wave…
▽ More
Binary black holes (BBHs) in eccentric orbits produce distinct modulations the emitted gravitational waves (GWs). The measurement of orbital eccentricity can provide robust evidence for dynamical binary formation channels. We analyze 57 GW events from the first, second and third observing runs of the LIGO-Virgo-KAGRA (LVK) Collaboration using a multipolar aligned-spin inspiral-merger-ringdown waveform model with two eccentric parameters: eccentricity and relativistic anomaly. This is made computationally feasible with the machine-learning code DINGO which accelerates inference by 2-3 orders of magnitude compared to traditional inference. First, we find eccentric aligned-spin versus quasi-circular aligned-spin $\log_{10}$ Bayes factors of 1.84 to 4.75 (depending on the glitch mitigation) for GW200129, 3.0 for GW190701 and 1.77 for GW200208_22. We measure $e_{\text{gw}, 10Hz}$ to be $0.27_{-0.12}^{+0.10}$ to $0.17_{-0.13}^{+0.14}$ for GW200129, $0.35_{-0.11}^{+0.32}$ for GW190701 and $0.35_{-0.21}^{+0.18}$ for GW200208_22. Second, we find $\log_{10}$ Bayes factors between the eccentric aligned-spin versus quasi-circular precessing-spin hypothesis between 1.43 and 4.92 for GW200129, 2.61 for GW190701 and 1.23 for GW200208_22. Third, our analysis does not show evidence for eccentricity in GW190521, which has an eccentric aligned-spin against quasi-circular aligned-spin $\log_{10}$ Bayes factor of 0.04. Fourth, we estimate that if we neglect the spin-precession and use an astrophysical prior, the probability of one out of the 57 events being eccentric is greater than 99.5% or $(100 - 8.4 \times 10^{-4})$% (depending on the glitch mitigation). Fifth, we study the impact on parameter estimation when neglecting either eccentricity or higher modes in eccentric models. These results underscore the importance of including eccentric parameters in the characterization of BBHs for GW detectors.
△ Less
Submitted 22 April, 2024;
originally announced April 2024.
-
Ultralight vector dark matter search using data from the KAGRA O3GK run
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
the KAGRA Collaboration,
A. G. Abac,
R. Abbott,
H. Abe,
I. Abouelfettouh,
F. Acernese,
K. Ackley,
C. Adamcewicz,
S. Adhicary,
N. Adhikari,
R. X. Adhikari,
V. K. Adkins,
V. B. Adya,
C. Affeldt,
D. Agarwal,
M. Agathos,
O. D. Aguiar,
I. Aguilar,
L. Aiello,
A. Ain,
P. Ajith,
T. Akutsu,
S. Albanesi
, et al. (1778 additional authors not shown)
Abstract:
Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we prese…
▽ More
Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we present the result of a search for $U(1)_{B-L}$ gauge boson DM using the KAGRA data from auxiliary length channels during the first joint observation run together with GEO600. By applying our search pipeline, which takes into account the stochastic nature of ultralight DM, upper bounds on the coupling strength between the $U(1)_{B-L}$ gauge boson and ordinary matter are obtained for a range of DM masses. While our constraints are less stringent than those derived from previous experiments, this study demonstrates the applicability of our method to the lower-mass vector DM search, which is made difficult in this measurement by the short observation time compared to the auto-correlation time scale of DM.
△ Less
Submitted 5 March, 2024;
originally announced March 2024.
-
Gravitational wave populations and cosmology with neural posterior estimation
Authors:
Konstantin Leyde,
Stephen R. Green,
Alexandre Toubiana,
Jonathan Gair
Abstract:
We apply neural posterior estimation for fast-and-accurate hierarchical Bayesian inference of gravitational wave populations. We use a normalizing flow to estimate directly the population hyper-parameters from a collection of individual source observations. This approach provides complete freedom in event representation, automatic inclusion of selection effects, and (in contrast to likelihood esti…
▽ More
We apply neural posterior estimation for fast-and-accurate hierarchical Bayesian inference of gravitational wave populations. We use a normalizing flow to estimate directly the population hyper-parameters from a collection of individual source observations. This approach provides complete freedom in event representation, automatic inclusion of selection effects, and (in contrast to likelihood estimation) without the need for stochastic samplers to obtain posterior samples. Since the number of events may be unknown when the network is trained, we split into sub-population analyses that we later recombine; this allows for fast sequential analyses as additional events are observed. We demonstrate our method on a toy problem of dark siren cosmology, and show that inference takes just a few minutes and scales to $\sim 600$ events before performance degrades. We argue that neural posterior estimation therefore represents a promising avenue for population inference with large numbers of events.
△ Less
Submitted 6 April, 2024; v1 submitted 20 November, 2023;
originally announced November 2023.
-
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.
-
Search for Eccentric Black Hole Coalescences during the Third Observing Run of LIGO and Virgo
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
the KAGRA Collaboration,
A. G. Abac,
R. Abbott,
H. Abe,
F. Acernese,
K. Ackley,
C. Adamcewicz,
S. Adhicary,
N. Adhikari,
R. X. Adhikari,
V. K. Adkins,
V. B. Adya,
C. Affeldt,
D. Agarwal,
M. Agathos,
O. D. Aguiar,
I. Aguilar,
L. Aiello,
A. Ain,
P. Ajith,
T. Akutsu,
S. Albanesi,
R. A. Alfaidi
, et al. (1750 additional authors not shown)
Abstract:
Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effect…
▽ More
Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effects of eccentricity. Here, we present observational results for a waveform-independent search sensitive to eccentric black hole coalescences, covering the third observing run (O3) of the LIGO and Virgo detectors. We identified no new high-significance candidates beyond those that were already identified with searches focusing on quasi-circular binaries. We determine the sensitivity of our search to high-mass (total mass $M>70$ $M_\odot$) binaries covering eccentricities up to 0.3 at 15 Hz orbital frequency, and use this to compare model predictions to search results. Assuming all detections are indeed quasi-circular, for our fiducial population model, we place an upper limit for the merger rate density of high-mass binaries with eccentricities $0 < e \leq 0.3$ at $0.33$ Gpc$^{-3}$ yr$^{-1}$ at 90\% confidence level.
△ Less
Submitted 7 August, 2023;
originally announced August 2023.
-
Adapting to noise distribution shifts in flow-based gravitational-wave inference
Authors:
Jonas Wildberger,
Maximilian Dax,
Stephen R. Green,
Jonathan Gair,
Michael Pürrer,
Jakob H. Macke,
Alessandra Buonanno,
Bernhard Schölkopf
Abstract:
Deep learning techniques for gravitational-wave parameter estimation have emerged as a fast alternative to standard samplers $\unicode{x2013}$ producing results of comparable accuracy. These approaches (e.g., DINGO) enable amortized inference by training a normalizing flow to represent the Bayesian posterior conditional on observed data. By conditioning also on the noise power spectral density (PS…
▽ More
Deep learning techniques for gravitational-wave parameter estimation have emerged as a fast alternative to standard samplers $\unicode{x2013}$ producing results of comparable accuracy. These approaches (e.g., DINGO) enable amortized inference by training a normalizing flow to represent the Bayesian posterior conditional on observed data. By conditioning also on the noise power spectral density (PSD) they can even account for changing detector characteristics. However, training such networks requires knowing in advance the distribution of PSDs expected to be observed, and therefore can only take place once all data to be analyzed have been gathered. Here, we develop a probabilistic model to forecast future PSDs, greatly increasing the temporal scope of DINGO networks. Using PSDs from the second LIGO-Virgo observing run (O2) $\unicode{x2013}$ plus just a single PSD from the beginning of the third (O3) $\unicode{x2013}$ we show that we can train a DINGO network to perform accurate inference throughout O3 (on 37 real events). We therefore expect this approach to be a key component to enable the use of deep learning techniques for low-latency analyses of gravitational waves.
△ Less
Submitted 16 November, 2022;
originally announced November 2022.
-
Neural Importance Sampling for Rapid and Reliable Gravitational-Wave Inference
Authors:
Maximilian Dax,
Stephen R. Green,
Jonathan Gair,
Michael Pürrer,
Jonas Wildberger,
Jakob H. Macke,
Alessandra Buonanno,
Bernhard Schölkopf
Abstract:
We combine amortized neural posterior estimation with importance sampling for fast and accurate gravitational-wave inference. We first generate a rapid proposal for the Bayesian posterior using neural networks, and then attach importance weights based on the underlying likelihood and prior. This provides (1) a corrected posterior free from network inaccuracies, (2) a performance diagnostic (the sa…
▽ More
We combine amortized neural posterior estimation with importance sampling for fast and accurate gravitational-wave inference. We first generate a rapid proposal for the Bayesian posterior using neural networks, and then attach importance weights based on the underlying likelihood and prior. This provides (1) a corrected posterior free from network inaccuracies, (2) a performance diagnostic (the sample efficiency) for assessing the proposal and identifying failure cases, and (3) an unbiased estimate of the Bayesian evidence. By establishing this independent verification and correction mechanism we address some of the most frequent criticisms against deep learning for scientific inference. We carry out a large study analyzing 42 binary black hole mergers observed by LIGO and Virgo with the SEOBNRv4PHM and IMRPhenomXPHM waveform models. This shows a median sample efficiency of $\approx 10\%$ (two orders-of-magnitude better than standard samplers) as well as a ten-fold reduction in the statistical uncertainty in the log evidence. Given these advantages, we expect a significant impact on gravitational-wave inference, and for this approach to serve as a paradigm for harnessing deep learning methods in scientific applications.
△ Less
Submitted 30 May, 2023; v1 submitted 11 October, 2022;
originally announced October 2022.
-
Using machine learning to parametrize postmerger signals from binary neutron stars
Authors:
Tim Whittaker,
William E. East,
Stephen R. Green,
Luis Lehner,
Huan Yang
Abstract:
There is growing interest in the detection and characterization of gravitational waves from postmerger oscillations of binary neutron stars. These signals contain information about the nature of the remnant and the high-density and out-of-equilibrium physics of the postmerger processes, which would complement any electromagnetic signal. However, the construction of binary neutron star postmerger w…
▽ More
There is growing interest in the detection and characterization of gravitational waves from postmerger oscillations of binary neutron stars. These signals contain information about the nature of the remnant and the high-density and out-of-equilibrium physics of the postmerger processes, which would complement any electromagnetic signal. However, the construction of binary neutron star postmerger waveforms is much more complicated than for binary black holes: (i) there are theoretical uncertainties in the neutron-star equation of state and other aspects of the high-density physics, (ii) numerical simulations are expensive and available ones only cover a small fraction of the parameter space with limited numerical accuracy, and (iii) it is unclear how to parametrize the theoretical uncertainties and interpolate across parameter space. In this work, we describe the use of a machine-learning method called a conditional variational autoencoder (CVAE) to construct postmerger models for hyper/massive neutron star remnant signals based on numerical-relativity simulations. The CVAE provides a probabilistic model, which encodes uncertainties in the training data within a set of latent parameters. We estimate that training such a model will ultimately require $\sim 10^4$ waveforms. However, using synthetic training waveforms as a proof-of-principle, we show that the CVAE can be used as an accurate generative model and that it encodes the equation of state in a useful latent representation.
△ Less
Submitted 17 January, 2022;
originally announced January 2022.
-
Nonlinear effects in the black hole ringdown: absorption-induced mode excitation
Authors:
Laura Sberna,
Pablo Bosch,
William E. East,
Stephen R. Green,
Luis Lehner
Abstract:
Gravitational-wave observations of black hole ringdowns are commonly used to characterize binary merger remnants and to test general relativity. These analyses assume linear black hole perturbation theory, in particular that the ringdown can be described in terms of quasinormal modes even for times approaching the merger. Here we investigate a nonlinear effect during the ringdown, namely how a mod…
▽ More
Gravitational-wave observations of black hole ringdowns are commonly used to characterize binary merger remnants and to test general relativity. These analyses assume linear black hole perturbation theory, in particular that the ringdown can be described in terms of quasinormal modes even for times approaching the merger. Here we investigate a nonlinear effect during the ringdown, namely how a mode excited at early times can excite additional modes as it is absorbed by the black hole. This is a third-order secular effect: the change in the black-hole mass causes a shift in the mode spectrum, so that the original mode is projected onto the new ones. Using nonlinear simulations, we study the ringdown of a spherically-symmetric scalar field around an asymptotically anti-de Sitter black hole, and we find that this "absorption-induced mode excitation" (AIME) is the dominant nonlinear effect. We show that this effect takes place well within the nonadiabatic regime, so we can analytically estimate it using a sudden mass-change approximation. Adapting our estimation technique to asymptotically-flat Schwarzschild black holes, we expect AIME to play a role in the analysis and interpretation of current and future gravitational wave observations.
△ Less
Submitted 21 December, 2021;
originally announced December 2021.
-
Group equivariant neural posterior estimation
Authors:
Maximilian Dax,
Stephen R. Green,
Jonathan Gair,
Michael Deistler,
Bernhard Schölkopf,
Jakob H. Macke
Abstract:
Simulation-based inference with conditional neural density estimators is a powerful approach to solving inverse problems in science. However, these methods typically treat the underlying forward model as a black box, with no way to exploit geometric properties such as equivariances. Equivariances are common in scientific models, however integrating them directly into expressive inference networks…
▽ More
Simulation-based inference with conditional neural density estimators is a powerful approach to solving inverse problems in science. However, these methods typically treat the underlying forward model as a black box, with no way to exploit geometric properties such as equivariances. Equivariances are common in scientific models, however integrating them directly into expressive inference networks (such as normalizing flows) is not straightforward. We here describe an alternative method to incorporate equivariances under joint transformations of parameters and data. Our method -- called group equivariant neural posterior estimation (GNPE) -- is based on self-consistently standardizing the "pose" of the data while estimating the posterior over parameters. It is architecture-independent, and applies both to exact and approximate equivariances. As a real-world application, we use GNPE for amortized inference of astrophysical binary black hole systems from gravitational-wave observations. We show that GNPE achieves state-of-the-art accuracy while reducing inference times by three orders of magnitude.
△ Less
Submitted 30 May, 2023; v1 submitted 25 November, 2021;
originally announced November 2021.
-
Gamma-radiation sky maps from compact binaries
Authors:
Néstor Ortiz,
Federico Carrasco,
Stephen R. Green,
Luis Lehner,
Steven L. Liebling,
John Ryan Westernacher-Schneider
Abstract:
We study sky maps and light curves of gamma-ray emission from neutron stars in compact binaries, and in isolation. We briefly review some gamma-ray emission models, and reproduce sky maps from a standard isolated pulsar in the Separatrix Layer model. We consider isolated pulsars with several variations of a dipole magnetic field, including superpositions, and predict their gamma-ray emission. Our…
▽ More
We study sky maps and light curves of gamma-ray emission from neutron stars in compact binaries, and in isolation. We briefly review some gamma-ray emission models, and reproduce sky maps from a standard isolated pulsar in the Separatrix Layer model. We consider isolated pulsars with several variations of a dipole magnetic field, including superpositions, and predict their gamma-ray emission. Our results provide new heuristics on what can and cannot be inferred about the magnetic field configuration of pulsars from high-energy observations. We find that typical double-peak light curves can be produced by pulsars with significant multipole structure beyond a single dipole. For binary systems, we also present a simple approximation that is useful for rapid explorations of binary magnetic field structure. Finally, we predict the gamma-ray emission pattern from a compact black hole-neutron star binary moments before merger by applying the Separatrix Layer model to data simulated in full general relativity; we find that face-on observers receive little emission, equatorial observers see one broad peak, and more generic observers typically see two peaks.
△ Less
Submitted 17 February, 2022; v1 submitted 14 July, 2021;
originally announced July 2021.
-
Real-time gravitational-wave science with neural posterior estimation
Authors:
Maximilian Dax,
Stephen R. Green,
Jonathan Gair,
Jakob H. Macke,
Alessandra Buonanno,
Bernhard Schölkopf
Abstract:
We demonstrate unprecedented accuracy for rapid gravitational-wave parameter estimation with deep learning. Using neural networks as surrogates for Bayesian posterior distributions, we analyze eight gravitational-wave events from the first LIGO-Virgo Gravitational-Wave Transient Catalog and find very close quantitative agreement with standard inference codes, but with inference times reduced from…
▽ More
We demonstrate unprecedented accuracy for rapid gravitational-wave parameter estimation with deep learning. Using neural networks as surrogates for Bayesian posterior distributions, we analyze eight gravitational-wave events from the first LIGO-Virgo Gravitational-Wave Transient Catalog and find very close quantitative agreement with standard inference codes, but with inference times reduced from O(day) to a minute per event. Our networks are trained using simulated data, including an estimate of the detector-noise characteristics near the event. This encodes the signal and noise models within millions of neural-network parameters, and enables inference for any observed data consistent with the training distribution, accounting for noise nonstationarity from event to event. Our algorithm -- called "DINGO" -- sets a new standard in fast-and-accurate inference of physical parameters of detected gravitational-wave events, which should enable real-time data analysis without sacrificing accuracy.
△ Less
Submitted 30 May, 2023; v1 submitted 23 June, 2021;
originally announced June 2021.
-
Complete parameter inference for GW150914 using deep learning
Authors:
Stephen R. Green,
Jonathan Gair
Abstract:
The LIGO and Virgo gravitational-wave observatories have detected many exciting events over the past five years. As the rate of detections grows with detector sensitivity, this poses a growing computational challenge for data analysis. With this in mind, in this work we apply deep learning techniques to perform fast likelihood-free Bayesian inference for gravitational waves. We train a neural-netw…
▽ More
The LIGO and Virgo gravitational-wave observatories have detected many exciting events over the past five years. As the rate of detections grows with detector sensitivity, this poses a growing computational challenge for data analysis. With this in mind, in this work we apply deep learning techniques to perform fast likelihood-free Bayesian inference for gravitational waves. We train a neural-network conditional density estimator to model posterior probability distributions over the full 15-dimensional space of binary black hole system parameters, given detector strain data from multiple detectors. We use the method of normalizing flows---specifically, a neural spline normalizing flow---which allows for rapid sampling and density estimation. Training the network is likelihood-free, requiring samples from the data generative process, but no likelihood evaluations. Through training, the network learns a global set of posteriors: it can generate thousands of independent posterior samples per second for any strain data consistent with the prior and detector noise characteristics used for training. By training with the detector noise power spectral density estimated at the time of GW150914, and conditioning on the event strain data, we use the neural network to generate accurate posterior samples consistent with analyses using conventional sampling techniques.
△ Less
Submitted 7 August, 2020;
originally announced August 2020.
-
Gravitational-wave parameter estimation with autoregressive neural network flows
Authors:
Stephen R. Green,
Christine Simpson,
Jonathan Gair
Abstract:
We introduce the use of autoregressive normalizing flows for rapid likelihood-free inference of binary black hole system parameters from gravitational-wave data with deep neural networks. A normalizing flow is an invertible mapping on a sample space that can be used to induce a transformation from a simple probability distribution to a more complex one: if the simple distribution can be rapidly sa…
▽ More
We introduce the use of autoregressive normalizing flows for rapid likelihood-free inference of binary black hole system parameters from gravitational-wave data with deep neural networks. A normalizing flow is an invertible mapping on a sample space that can be used to induce a transformation from a simple probability distribution to a more complex one: if the simple distribution can be rapidly sampled and its density evaluated, then so can the complex distribution. Our first application to gravitational waves uses an autoregressive flow, conditioned on detector strain data, to map a multivariate standard normal distribution into the posterior distribution over system parameters. We train the model on artificial strain data consisting of IMRPhenomPv2 waveforms drawn from a five-parameter $(m_1, m_2, φ_0, t_c, d_L)$ prior and stationary Gaussian noise realizations with a fixed power spectral density. This gives performance comparable to current best deep-learning approaches to gravitational-wave parameter estimation. We then build a more powerful latent variable model by incorporating autoregressive flows within the variational autoencoder framework. This model has performance comparable to Markov chain Monte Carlo and, in particular, successfully models the multimodal $φ_0$ posterior. Finally, we train the autoregressive latent variable model on an expanded parameter space, including also aligned spins $(χ_{1z}, χ_{2z})$ and binary inclination $θ_{JN}$, and show that all parameters and degeneracies are well-recovered. In all cases, sampling is extremely fast, requiring less than two seconds to draw $10^4$ posterior samples.
△ Less
Submitted 18 February, 2020;
originally announced February 2020.
-
Extended I-Love relations for slowly rotating neutron stars
Authors:
Jérémie Gagnon-Bischoff,
Stephen R. Green,
Philippe Landry,
Néstor Ortiz
Abstract:
Observations of gravitational waves from inspiralling neutron star binaries---such as GW170817---can be used to constrain the nuclear equation of state by placing bounds on stellar tidal deformability. For slowly rotating neutron stars, the response to a weak quadrupolar tidal field is characterized by four internal-structure-dependent constants called "Love numbers." The tidal Love numbers…
▽ More
Observations of gravitational waves from inspiralling neutron star binaries---such as GW170817---can be used to constrain the nuclear equation of state by placing bounds on stellar tidal deformability. For slowly rotating neutron stars, the response to a weak quadrupolar tidal field is characterized by four internal-structure-dependent constants called "Love numbers." The tidal Love numbers $k_2^\text{el}$ and $k_2^\text{mag}$ measure the tides raised by the gravitoelectric and gravitomagnetic components of the applied field, and the rotational-tidal Love numbers $\mathfrak{f}^\text{o}$ and $\mathfrak{k}^\text{o}$ measure those raised by couplings between the applied field and the neutron star spin. In this work we compute these four Love numbers for perfect fluid neutron stars with realistic equations of state. We discover (nearly) equation-of-state independent relations between the rotational-tidal Love numbers and the moment of inertia, thereby extending the scope of I-Love-Q universality. We find that similar relations hold among the tidal and rotational-tidal Love numbers. These relations extend the applications of I-Love universality in gravitational-wave astronomy. As our findings differ from those reported in the literature, we derive general formulas for the rotational-tidal Love numbers in post-Newtonian theory and confirm numerically that they agree with our general-relativistic computations in the weak-field limit.
△ Less
Submitted 2 April, 2018; v1 submitted 15 November, 2017;
originally announced November 2017.
-
A Simple, Heuristic Derivation of our "No Backreaction" Results
Authors:
Stephen R. Green,
Robert M. Wald
Abstract:
We provide a simple discussion of our results on the backreaction effects of density inhomogeneities in cosmology, without mentioning one-parameter families or weak limits. Emphasis is placed on the manner in which "averaging" is done and the fact that one is solving Einstein's equation. The key assumptions and results that we rigorously derived within our original mathematical framework are there…
▽ More
We provide a simple discussion of our results on the backreaction effects of density inhomogeneities in cosmology, without mentioning one-parameter families or weak limits. Emphasis is placed on the manner in which "averaging" is done and the fact that one is solving Einstein's equation. The key assumptions and results that we rigorously derived within our original mathematical framework are thereby explained in a heuristic way.
△ Less
Submitted 25 January, 2016;
originally announced January 2016.
-
Nonlinear Evolution and Final Fate of Charged Anti-de Sitter Black Hole Superradiant Instability
Authors:
Pablo Bosch,
Stephen R. Green,
Luis Lehner
Abstract:
We describe the full nonlinear development of the superradiant instability for a charged massless scalar field, coupled to general relativity and electromagnetism, in the vicinity of a Reissner-Nordstrom-AdS black hole. The presence of the negative cosmological constant provides a natural context for considering perfectly reflecting boundary conditions and studying the dynamics as the scalar field…
▽ More
We describe the full nonlinear development of the superradiant instability for a charged massless scalar field, coupled to general relativity and electromagnetism, in the vicinity of a Reissner-Nordstrom-AdS black hole. The presence of the negative cosmological constant provides a natural context for considering perfectly reflecting boundary conditions and studying the dynamics as the scalar field interacts repeatedly with the black hole. At early times, small superradiant perturbations grow as expected from linearized studies. Backreaction then causes the black hole to lose charge and mass until the perturbation becomes nonsuperradiant, with the final state described by a stable hairy black hole. For large gauge coupling, the instability extracts a large amount of charge per unit mass, resulting in greater entropy increase. We discuss the implications of the observed behavior for the general problem of superradiance in black hole spacetimes.
△ Less
Submitted 30 March, 2016; v1 submitted 6 January, 2016;
originally announced January 2016.
-
Comments on Backreaction
Authors:
Stephen R. Green,
Robert M. Wald
Abstract:
We respond to the criticisms of a recent paper of Buchert et al. [arXiv:1505.07800]
We respond to the criticisms of a recent paper of Buchert et al. [arXiv:1505.07800]
△ Less
Submitted 29 October, 2015; v1 submitted 21 June, 2015;
originally announced June 2015.
-
How well is our universe described by an FLRW model?
Authors:
Stephen R. Green,
Robert M. Wald
Abstract:
Extremely well! In the $Λ$CDM model, the spacetime metric, $g_{ab}$, of our universe is approximated by an FLRW metric, $g_{ab}^{(0)}$, to about 1 part in $10^4$ or better on both large and small scales, except in the immediate vicinity of very strong field objects, such as black holes. However, derivatives of $g_{ab}$ are not close to derivatives of $g_{ab}^{(0)}$, so there can be significant dif…
▽ More
Extremely well! In the $Λ$CDM model, the spacetime metric, $g_{ab}$, of our universe is approximated by an FLRW metric, $g_{ab}^{(0)}$, to about 1 part in $10^4$ or better on both large and small scales, except in the immediate vicinity of very strong field objects, such as black holes. However, derivatives of $g_{ab}$ are not close to derivatives of $g_{ab}^{(0)}$, so there can be significant differences in the behavior of geodesics and huge differences in curvature. Consequently, observable quantities in the actual universe may differ significantly from the corresponding observables in the FLRW model. Nevertheless, as we shall review here, we have proven general results showing that---within the framework of our approach to treating backreaction---the large matter inhomogeneities that occur on small scales cannot produce significant effects on large scales, so $g_{ab}^{(0)}$ satisfies Einstein's equation with the averaged stress-energy tensor of matter as its source. We discuss the flaws in some other approaches that have suggested that large backreaction effects may occur. As we also will review here, with a suitable "dictionary," Newtonian cosmologies provide excellent approximations to cosmological solutions to Einstein's equation (with dust and a cosmological constant) on all scales. Our results thereby provide strong justification for the mathematical consistency and validity of the $Λ$CDM model within the context of general relativistic cosmology.
△ Less
Submitted 16 October, 2014; v1 submitted 30 July, 2014;
originally announced July 2014.
-
A Holographic Path to the Turbulent Side of Gravity
Authors:
Stephen R. Green,
Federico Carrasco,
Luis Lehner
Abstract:
We study the dynamics of a 2+1 dimensional relativistic viscous conformal fluid in Minkowski spacetime. Such fluid solutions arise as duals, under the "gravity/fluid correspondence", to 3+1 dimensional asymptotically anti-de Sitter (AAdS) black brane solutions to the Einstein equation. We examine stability properties of shear flows, which correspond to hydrodynamic quasinormal modes of the black b…
▽ More
We study the dynamics of a 2+1 dimensional relativistic viscous conformal fluid in Minkowski spacetime. Such fluid solutions arise as duals, under the "gravity/fluid correspondence", to 3+1 dimensional asymptotically anti-de Sitter (AAdS) black brane solutions to the Einstein equation. We examine stability properties of shear flows, which correspond to hydrodynamic quasinormal modes of the black brane. We find that, for sufficiently high Reynolds number, the solution undergoes an inverse turbulent cascade to long wavelength modes. We then map this fluid solution, via the gravity/fluid duality, into a bulk metric. This suggests a new and interesting feature of the behavior of perturbed AAdS black holes and black branes, which is not readily captured by a standard quasinormal mode analysis. Namely, for sufficiently large perturbed black objects (with long-lived quasinormal modes), nonlinear effects transfer energy from short to long wavelength modes via a turbulent cascade within the metric perturbation. As long wavelength modes have slower decay, this lengthens the overall lifetime of the perturbation. We also discuss various implications of this behavior, including expectations for higher dimensions, and the possibility of predicting turbulence in more general gravitational scenarios.
△ Less
Submitted 7 October, 2013; v1 submitted 30 September, 2013;
originally announced September 2013.
-
Dynamic and Thermodynamic Stability of Relativistic, Perfect Fluid Stars
Authors:
Stephen R. Green,
Joshua S. Schiffrin,
Robert M. Wald
Abstract:
We consider perfect fluid bodies (stars) in general relativity, characterized by particle number density and entropy per particle. A star is said to be in dynamic equilibrium if it is a stationary, axisymmetric solution to the Einstein-fluid equations with circular flow. We prove that for a star in dynamic equilibrium, the necessary and sufficient condition for thermodynamic equilibrium (extremum…
▽ More
We consider perfect fluid bodies (stars) in general relativity, characterized by particle number density and entropy per particle. A star is said to be in dynamic equilibrium if it is a stationary, axisymmetric solution to the Einstein-fluid equations with circular flow. We prove that for a star in dynamic equilibrium, the necessary and sufficient condition for thermodynamic equilibrium (extremum of total entropy S) is constancy of angular velocity (Ω), redshifted temperature, and redshifted chemical potential. Friedman previously identified positivity of canonical energy, E, as a criterion for dynamic stability and argued that all rotating stars are dynamically unstable to sufficiently nonaxisymmetric perturbations (the CFS instability), so our main focus is on axisymmetric stability. We show that for a star in dynamic equilibrium, mode stability holds with respect to all axisymmetric perturbations if E is positive on a certain subspace, V, of axisymmetric Lagrangian perturbations that, in particular, have no Lagrangian change in angular momentum density. Conversely, if E fails to be positive on V, then there exist perturbations that can't become asymptotically stationary at late times. We further show that for a star in thermodynamic equilibrium, the canonical energy in the rotating frame, E_r, is related to second order changes in ADM mass, M, and angular momentum, J, by E_r = δ^2 M - Ωδ^2 J. Thus, positivity of E_r for perturbations with δJ = 0 is a necessary condition for thermodynamic stability (local maximum of S). For axisymmetric perturbations we have E = E_r, so a necessary condition for thermodynamic stability with respect to axisymmetric perturbations is positivity of E on all perturbations with δJ = 0, not merely on the perturbations in V. Many of our results are in close parallel to results of Hollands and Wald for the theory of black holes.
△ Less
Submitted 22 December, 2013; v1 submitted 31 August, 2013;
originally announced September 2013.
-
Examples of backreaction of small scale inhomogeneities in cosmology
Authors:
Stephen R. Green,
Robert M. Wald
Abstract:
In previous work, we introduced a new framework to treat large scale backreaction effects due to small scale inhomogeneities in general relativity. We considered one-parameter families of spacetimes for which such backreaction effects can occur, and we proved that, provided the weak energy condition on matter is satisfied, the leading effect of small scale inhomogeneities on large scale dynamics i…
▽ More
In previous work, we introduced a new framework to treat large scale backreaction effects due to small scale inhomogeneities in general relativity. We considered one-parameter families of spacetimes for which such backreaction effects can occur, and we proved that, provided the weak energy condition on matter is satisfied, the leading effect of small scale inhomogeneities on large scale dynamics is to produce a traceless effective stress-energy tensor that itself satisfies the weak energy condition. In this work, we illustrate the nature of our framework by providing two explicit examples of one-parameter families with backreaction. The first, based on previous work of Berger, is a family of polarized vacuum Gowdy spacetimes on a torus, which satisfies all of the assumptions of our framework. As the parameter approaches its limiting value, the metric uniformly approaches a smooth background metric, but spacetime derivatives of the deviation of the metric from the background metric do not converge uniformly to zero. The limiting metric has nontrivial backreaction from the small scale inhomogeneities, with an effective stress-energy that is traceless and satisfies the weak energy condition, in accord with our theorems. Our second one-parameter family consists of metrics which have a uniform Friedmann-Lemaitre-Robertson-Walker limit. This family satisfies all of our assumptions with the exception of the weak energy condition for matter. In this case, the limiting metric has an effective stress-energy tensor which is not traceless. We emphasize the importance of imposing energy conditions on matter in studies of backreaction.
△ Less
Submitted 19 June, 2013; v1 submitted 8 April, 2013;
originally announced April 2013.
-
Newtonian and Relativistic Cosmologies
Authors:
Stephen R. Green,
Robert M. Wald
Abstract:
Cosmological N-body simulations are now being performed using Newtonian gravity on scales larger than the Hubble radius. It is well known that a uniformly expanding, homogeneous ball of dust in Newtonian gravity satisfies the same equations as arise in relativistic FLRW cosmology, and it also is known that a correspondence between Newtonian and relativistic dust cosmologies continues to hold in li…
▽ More
Cosmological N-body simulations are now being performed using Newtonian gravity on scales larger than the Hubble radius. It is well known that a uniformly expanding, homogeneous ball of dust in Newtonian gravity satisfies the same equations as arise in relativistic FLRW cosmology, and it also is known that a correspondence between Newtonian and relativistic dust cosmologies continues to hold in linearized perturbation theory in the marginally bound/spatially flat case. Nevertheless, it is far from obvious that Newtonian gravity can provide a good global description of an inhomogeneous cosmology when there is significant nonlinear dynamical behavior at small scales. We investigate this issue in the light of a perturbative framework that we have recently developed, which allows for such nonlinearity at small scales. We propose a relatively straightforward "dictionary"---which is exact at the linearized level---that maps Newtonian dust cosmologies into general relativistic dust cosmologies, and we use our "ordering scheme" to determine the degree to which the resulting metric and matter distribution solve Einstein's equation. We find that Einstein's equation fails to hold at "order 1" at small scales and at "order $ε$" at large scales. We then find the additional corrections to the metric and matter distribution needed to satisfy Einstein's equation to these orders. While these corrections are of some interest in their own right, our main purpose in calculating them is that their smallness should provide a criterion for the validity of the original dictionary (as well as simplified versions of this dictionary). We expect that, in realistic Newtonian cosmologies, these additional corrections will be very small; if so, this should provide strong justification for the use of Newtonian simulations to describe relativistic cosmologies, even on scales larger than the Hubble radius.
△ Less
Submitted 7 February, 2012; v1 submitted 13 November, 2011;
originally announced November 2011.
-
Constraints on String Cosmology
Authors:
Stephen R. Green,
Emil J. Martinec,
Callum Quigley,
Savdeep Sethi
Abstract:
String theory contains sources like orientifold planes that support higher derivative interactions. These interactions make possible static flux compactifications which are forbidden in supergravity. They can also lead to violations of the strong energy condition (SEC) which is needed for an accelerating universe. We examine how large a violation is possible in the context of the heterotic string…
▽ More
String theory contains sources like orientifold planes that support higher derivative interactions. These interactions make possible static flux compactifications which are forbidden in supergravity. They can also lead to violations of the strong energy condition (SEC) which is needed for an accelerating universe. We examine how large a violation is possible in the context of the heterotic string compactified to four dimensions. We find that de Sitter solutions are still not possible but that classically forbidden anti-de Sitter solutions are possible.
△ Less
Submitted 29 February, 2012; v1 submitted 3 October, 2011;
originally announced October 2011.
-
Numerical Parameter Survey of Nonradiative Black Hole Accretion -- Flow Structure and Variability of the Rotation Measure
Authors:
Bijia Pang,
Ue-Li Pen,
Christopher D. Matzner,
Stephen R. Green,
Matthias Liebendörfer
Abstract:
We conduct a survey of numerical simulations to probe the structure and appearance of non-radiative black hole accretion flows like the supermassive black hole at the Galactic centre. We find a generic set of solutions, and make specific predictions for currently feasible rotation measure (RM) observations, which are accessible to current instruments including the EVLA, GMRT and ALMA. The slow tim…
▽ More
We conduct a survey of numerical simulations to probe the structure and appearance of non-radiative black hole accretion flows like the supermassive black hole at the Galactic centre. We find a generic set of solutions, and make specific predictions for currently feasible rotation measure (RM) observations, which are accessible to current instruments including the EVLA, GMRT and ALMA. The slow time variability of the RM is a key quantitative signature of this accretion flow. The time variability of RM can be used to quantitatively measure the nature of the accretion flow, and to differentiate models. Sensitive measurements of RM can be achieved using RM synthesis or using pulsars.
Our energy conserving ideal magneto-hydrodynamical simulations, which achieve high dynamical range by means of a deformed-mesh algorithm, stretch from several Bondi radii to about one thousandth of that radius, and continue for tens of Bondi times. Magnetized flows which lack outward convection possess density slopes around -1, almost independent of physical parameters, and are more consistent with observational constraints than are strongly convective flows We observe no tendency for the flows to become rotationally supported in their centres, or to develop steady outflow.
We support these conclusions with formulae which encapsulate our findings in terms of physical and numerical parameters. We discuss the relation of these solutions to other approaches. The main potential uncertainties are the validity of ideal MHD and the absence of a fully relativistic inner boundary condition. The RM variability predictions are testable with current and future telescopes.
△ Less
Submitted 24 November, 2010; v1 submitted 24 November, 2010;
originally announced November 2010.
-
A new framework for analyzing the effects of small scale inhomogeneities in cosmology
Authors:
Stephen R. Green,
Robert M. Wald
Abstract:
We develop a new, mathematically precise framework for treating the effects of nonlinear phenomena occurring on small scales in general relativity. Our approach is an adaptation of Burnett's formulation of the "shortwave approximation", which we generalize to analyze the effects of matter inhomogeneities as well as gravitational radiation. Our framework requires the metric to be close to a "backgr…
▽ More
We develop a new, mathematically precise framework for treating the effects of nonlinear phenomena occurring on small scales in general relativity. Our approach is an adaptation of Burnett's formulation of the "shortwave approximation", which we generalize to analyze the effects of matter inhomogeneities as well as gravitational radiation. Our framework requires the metric to be close to a "background metric", but allows arbitrarily large stress-energy fluctuations on small scales. We prove that, within our framework, if the matter stress-energy tensor satisfies the weak energy condition (i.e., positivity of energy density in all frames), then the only effect that small scale inhomogeneities can have on the dynamics of the background metric is to provide an "effective stress-energy tensor" that is traceless and has positive energy density---corresponding to the presence of gravitational radiation. In particular, nonlinear effects produced by small scale inhomogeneities cannot mimic the effects of dark energy. We also develop "perturbation theory" off of the background metric. We derive an equation for the "long-wavelength part" of the leading order deviation of the metric from the background metric, which contains the usual terms occurring in linearized perturbation theory plus additional contributions from the small-scale inhomogeneities. Under various assumptions concerning the absence of gravitational radiation and the non-relativistic behavior of the matter, we argue that the "short wavelength" deviations of the metric from the background metric near a point $x$ should be accurately described by Newtonian gravity, taking into account only the matter lying within a "homogeneity lengthscale" of $x$. Finally, we argue that our framework should provide an accurate description of the actual universe.
△ Less
Submitted 7 April, 2011; v1 submitted 22 November, 2010;
originally announced November 2010.