Bayesian semiparametric power spectral density estimation with applications in gravitational wave data analysis

Abstract

The standard noise model in gravitational wave (GW) data analysis assumes detector noise is stationary and Gaussian distributed, with a known power spectral density (PSD) that is usually estimated using clean off-source data. Real GW data often depart from these assumptions, and misspecified parametric models of the PSD could result in misleading inferences. We propose a Bayesian semiparametric approach to improve this. We use a nonparametric Bernstein polynomial prior on the PSD, with weights attained via a Dirichlet process distribution, and update this using the Whittle likelihood. Posterior samples are obtained using a blocked Metropolis-within-Gibbs sampler. We simultaneously estimate the reconstruction parameters of a rotating core collapse supernova GW burst that has been embedded in simulated Advanced LIGO noise. We also discuss an approach to deal with non-stationary data by breaking longer data streams into smaller and locally stationary components.

15 pages, 15 figures

References in corpus (10)

• Advanced Virgo: a 2nd generation interferometric gravitational wave detector
• Advanced LIGO
• BayesWave: Bayesian Inference for Gravitational Wave Bursts and Instrument Glitches
• BayesLine: Bayesian Inference for Spectral Estimation of Gravitational Wave Detector Noise
• Asymptotic spectral theory for nonlinear time series
• Measuring the Angular Momentum Distribution in Core-Collapse Supernova Progenitors with Gravitational Waves
• Rotating stellar core-collapse waveform decompositon: a Principal Component Analysis approach
• Separating Gravitational Wave Signals from Instrument Artifacts
• Data series subtraction with unknown and unmodeled background noise
• Bayesian parameter estimation of core collapse supernovae using gravitational wave simulations

Cited by in corpus (19)

• Parameter estimation with gravitational waves
• Parameter Estimation and Model Selection of Gravitational Wave Signals Contaminated by Transient Detector Noise Glitches
• Gravitational-wave parameter estimation with gaps in LISA: a Bayesian data augmentation method
• GPU-accelerated massive black hole binary parameter estimation with LISA
• Time-Frequency Analysis of Gravitational Wave Data
• Quantifying the Effect of Power Spectral Density Uncertainty on Gravitational-Wave Parameter Estimation for Compact Binary Sources
• Uncovering gravitational-wave backgrounds from noises of unknown shape with LISA
• Accurate modeling and mitigation of overlapping signals and glitches in gravitational-wave data
• Bayesian nonparametric spectral density estimation using B-spline priors
• Stochastic gravitational wave background: methods and Implications
• Identifying and Addressing Nonstationary LISA Noise
• Total-variation methods for gravitational-wave denoising: performance tests on Advanced LIGO data
• Impact of LIGO-Virgo black hole binaries on gravitational wave background searches
• Concurrent estimation of noise and compact-binary signal parameters in gravitational-wave data
• Degeneracy of gravitational waveforms in the context of GW150914
• Heavy-tailed likelihoods for robustness against data outliers: Applications to the analysis of gravitational wave data
• Formalism for power spectral density estimation for non-identical and correlated noise using the null channel in Einstein Telescope
• Computational Techniques for Parameter Estimation of Gravitational Wave Signals
• Beyond Whittle: Nonparametric correction of a parametric likelihood with a focus on Bayesian time series analysis