Fast LISA likelihood approximations by downsampling

Jethro Linley
2024-02-03
Abstract:The LISA gravitational wave observatory is due to launch in 2035. Researchers are currently working out the data analysis picture, including determining requirements of, and prospects for, parameterised gravitational wave models used in Bayesian inference, of which early stage inspirals of low and intermediate mass black hole binaries constitute a significant and important part. With datasets consisting of up to hundreds of millions of datapoints, their likelihood functions can be extremely expensive to compute. We present an approximation procedure for accurately reproducing the likelihood function (from time-domain data) for such signals in simulated environments, significantly reducing this cost. The method is simply to discard most of the datapoints and define a new inner product operator in terms of the remaining data, which closely resembles the original inner product on the model parameter space. Highly accurate reproductions of the likelihood can be achieved using just a few hundred to a few thousand datapoints, which can reduce parameter estimation (PE) time by factors of up to tens of thousands (compared to estimated frequency domain PE time), however, the end user must be responsible for ensuring convergence. The time-domain implementation is particularly useful for modelling waveform modifications arising from non-trivial astrophysical environments. The contents of this article provide the theoretical basis of the software package dolfen.
General Relativity and Quantum Cosmology,High Energy Astrophysical Phenomena,Instrumentation and Methods for Astrophysics
What problem does this paper attempt to address?