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 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, \phi_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 $\phi_0$ posterior. Finally, we train the autoregressive latent variable model on an expanded parameter space, including also aligned spins $(\chi_{1z}, \chi_{2z})$ and binary inclination $\theta_{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.

Highlights

  • The task of gravitational-wave parameter estimation is to determine the system parameters that gave rise to observed detector strain data

  • In this subsection we review the concept of a masked autoregressive flow, a type of normalizing flow that we use in our work to map simple distributions into more complex ones

  • In this work we introduced the use of masked autoregressive flows for improving deep learning of gravitationalwave posterior distributions

Read more

Summary

INTRODUCTION

The task of gravitational-wave parameter estimation is to determine the system parameters that gave rise to observed detector strain data. Since conditional distributions can be parametrized by functions, neural networks can be used to model gravitational-wave posteriors. We use the method of normalizing flows [9] to define very flexible model distributions with improved ability to approximate complex gravitational-wave posteriors. To describe a conditional distribution, such as a gravitational-wave posterior, a normalizing flow can be made to depend on additional variables; for our application, we take f 1⁄4 fyðuÞ to depend on detector strain data y. Since the density can be evaluated directly and differentiated with respect to the neural network parameters defining the flow, the model can be trained using stochastic gradient descent to maximize the likelihood that the training data [ðx; yÞ pairs] came from the model.

PRELIMINARIES
Neural network models of gravitational-wave posteriors
Variational autoencoders
Masked autoregressive flows
Combined models
Encoder with inverse autoregressive flow
Decoder with masked autoregressive flow
Prior with masked autoregressive flow
NONSPINNING BINARIES
Training data
INCLINED BINARIES WITH ALIGNED SPINS
CONCLUSIONS
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call