Ab Initio Real-Time Gravitational-Wave Parameter Estimation

We present a specialised GPU-native nested sampling kernel targeting rapid parameter estimation for gravitational wave inference problems. Building upon a Slice-within-Gibbs (SwiG) structure for rapid mixing, we investigate how far we can push baseline stochastic sampling techniques on modern GPU hardware. We demonstrate that for typical long-duration binary neutron star signals observed by the LIGO and Virgo detectors, we can achieve well calibrated posterior inference on the full uncompressed data of a three detector network in a median of twelve minutes on a single GPU. This falls to five minutes when sharded across four devices. Utilising heterodyning to compress the data reduces the median wall time across an injection campaign to 89 seconds -- less than the length of the segment itself -- and enables inference with precessing spin, tidal waveforms on GW170817 in around two minutes. This pushes stochastic sampling techniques using full physical waveform calculations, launched from an uninformed prior state, towards real-time gravitational wave parameter estimation.