The Deledda project developed, implemented, and validated new machine learning methods for fast and reliable inference in gravitational-wave astronomy. The work addressed three main fronts corresponding to different sources and analysis challenges: compact binary coalescences observed by ground-based detectors, the nanohertz gravitational-wave background probed by pulsar timing arrays, and the general problem of computing Bayesian evidence for model comparison.
The first achievement is the development of Labrador, a simulation-based inference framework that combines neural posterior estimation with domain-specific physical insights. The method compresses detector data through heterodyning against an optimal reference waveform, reparametrizes source parameters to remove degeneracies, and folds the parameter space to eliminate known multimodalities. These design choices make the network approximately equivariant to changes in source parameters, improving both efficiency and interpretability. Labrador achieves state-of-the-art performance with a full end-to-end training time of about one day on a single A100 GPU, representing a major step toward real-time parameter estimation for gravitational-wave events.
The second line of work introduced variational Bayesian inference as a new approach for analyzing pulsar-timing-array datasets. Unlike traditional Markov Chain Monte Carlo techniques, this method optimizes a neural approximation to the posterior distribution using stochastic gradient descent, allowing it to fully exploit the parallelism of modern GPUs. When applied to the NANOGrav 15-year dataset, the approach reduced the analysis time from days to minutes while maintaining statistical accuracy. This breakthrough opens the door to systematic studies of model uncertainties and alternative astrophysical or cosmological scenarios using PTA data.
Finally, the project developed floZ, a general-purpose algorithm to estimate Bayesian evidence directly from posterior samples. Based on normalizing flows, floZ is accurate, robust to sharp posterior features, and scalable to high-dimensional spaces. It provides an efficient alternative to nested sampling and other evidence estimators, and can be integrated with variational or simulation-based inference pipelines.
Together, these results demonstrate the potential of deep learning to transform the analysis of gravitational-wave data, making it faster, more scalable, and physically grounded.