Source author record

Renate Meyer

Renate Meyer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

gr-qc physics.data-an astro-ph.IM Methodology Applications astro-ph Computation

Catalog footprint

What is connected

11works

7topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Ability of LISA to detect a gravitational-wave background of cosmological origin: the cosmic string case

We investigate the ability of the Laser Interferometer Space Antenna (LISA) to detect a stochastic gravitational-wave background (GWB) produced by cosmic strings, and to subsequently estimate the string tension $Gμ$ in the presence of instrument noise, an astrophysical background from compact binaries, and the galactic foreground from white dwarf binaries. Fisher Information and Markov Chain Monte Carlo methods provide estimates of the LISA noise and the parameters for the different signal sources. We demonstrate the ability of LISA to simultaneously estimate the galactic foreground, as well as the astrophysical and cosmic string produced backgrounds. Considering the expected astrophysical background and a galactic foreground, a cosmic string tension in the $G μ\approx 10^{-16}$ to $Gμ\approx 10^{-15}$ range or bigger could be measured by LISA, with the galactic foreground affecting this limit more than the astrophysical background. The parameter estimation methods presented here can be applied to other cosmological backgrounds in the LISA observation band.

preprint2022arXiv

Parameter estimation with gravitational waves

The new era of gravitational wave astronomy truly began on September 14, 2015 with the detection of GW150914, the sensational first direct observation of gravitational waves from the inspiral and merger of two black holes by the two Advanced LIGO detectors. In the subsequent first three observing runs of the LIGO/Virgo network, gravitational waves from $\sim 50$ compact binary mergers have been announced, with more results to come. The events have mostly been produced by binary black holes, but two binary neutron star mergers have so far been observed, as well as the mergers of two neutron star - black hole systems. Furthermore, gravitational waves emitted by core-collapse supernovae, pulsars and the stochastic gravitational wave background are within the LIGO/Virgo/KAGRA sensitivity band and are likely to be observed in future observation runs. Beyond signal detection, a major challenge has been the development of statistical and computational methodology for estimating the physical waveform parameters and quantifying their uncertainties in order to accurately characterise the emitting system. These methods depend on the sources of the gravitational waves and the gravitational waveform model that is used. This article reviews the main waveform models and parameter estimation methods used to extract physical parameters from gravitational wave signals detected to date by LIGO and Virgo and from those expected to be observed in the future, which will include KAGRA, and how these methods interface with various aspects of LIGO/Virgo/KAGRA science. Also presented are the statistical methods used by LIGO and Virgo to estimate detector noise, test general relativity, and draw conclusions about the rates of compact binary mergers in the universe. Furthermore, a summary of major publicly available gravitational wave parameter estimation software packages is given.

preprint2021arXiv

Bayesian spectral density estimation using P-splines with quantile-based knot placement

This article proposes a Bayesian approach to estimating the spectral density of a stationary time series using a prior based on a mixture of P-spline distributions. Our proposal is motivated by the B-spline Dirichlet process prior of Edwards et al. (2019) in combination with Whittle's likelihood and aims at reducing the high computational complexity of its posterior computations. The strength of the B-spline Dirichlet process prior over the Bernstein-Dirichlet process prior of Choudhuri et al. (2004) lies in its ability to estimate spectral densities with sharp peaks and abrupt changes due to the flexibility of B-splines with variable number and location of knots. Here, we suggest to use P-splines of Eilers and Marx (1996) that combine a B-spline basis with a discrete penalty on the basis coefficients. In addition to equidistant knots, a novel strategy for a more expedient placement of knots is proposed that makes use of the information provided by the periodogram about the steepness of the spectral power distribution. We demonstrate in a simulation study and two real case studies that this approach retains the flexibility of the B-splines, achieves similar ability to accurately estimate peaks due to the new data-driven knot allocation scheme but significantly reduces the computational costs.

preprint2020arXiv

Computational Techniques for Parameter Estimation of Gravitational Wave Signals

Since the very first detection of gravitational waves from the coalescence of two black holes in 2015, Bayesian statistical methods have been routinely applied by LIGO and Virgo to extract the signal out of noisy interferometric measurements, obtain point estimates of the physical parameters responsible for producing the signal, and rigorously quantify their uncertainties. Different computational techniques have been devised depending on the source of the gravitational radiation and the gravitational waveform model used. Prominent sources of gravitational waves are binary black hole or neutron star mergers, the only objects that have been observed by detectors to date. But also gravitational waves from core collapse supernovae, rapidly rotating neutron stars, and the stochastic gravitational wave background are in the sensitivity band of the ground-based interferometers and expected to be observable in future observation runs. As nonlinearities of the complex waveforms and the high-dimensional parameter spaces preclude analytic evaluation of the posterior distribution, posterior inference for all these sources relies on computer-intensive simulation techniques such as Markov chain Monte Carlo methods. A review of state-of-the-art Bayesian statistical parameter estimation methods will be given for researchers in this cross-disciplinary area of gravitational wave data analysis.

preprint2015arXiv

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

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.

preprint2014arXiv

Bayesian parameter estimation of core collapse supernovae using gravitational wave simulations

Using the latest numerical simulations of rotating stellar core collapse, we present a Bayesian framework to extract the physical information encoded in noisy gravitational wave signals. We fit Bayesian principal component regression models with known and unknown signal arrival times to reconstruct gravitational wave signals, and subsequently fit known astrophysical parameters on the posterior means of the principal component coefficients using a linear model. We predict the ratio of rotational kinetic energy to gravitational energy of the inner core at bounce by sampling from the posterior predictive distribution, and find that these predictions are generally very close to the true parameter values, with $90\%$ credible intervals $\sim 0.04$ and $\sim 0.06$ wide for the known and unknown arrival time models respectively. Two supervised machine learning methods are implemented to classify precollapse differential rotation, and we find that these methods discriminate rapidly rotating progenitors particularly well. We also introduce a constrained optimization approach to model selection to find an optimal number of principal components in the signal reconstruction step. Using this approach, we select 14 principal components as the most parsimonious model.

preprint2013arXiv

Bayesian inference on EMRI signals using low frequency approximations

Extreme mass ratio inspirals (EMRIs) are thought to be one of the most exciting gravitational wave sources to be detected with LISA. Due to their complicated nature and weak amplitudes the detection and parameter estimation of such sources is a challenging task. In this paper we present a statistical methodology based on Bayesian inference in which the estimation of parameters is carried out by advanced Markov chain Monte Carlo (MCMC) algorithms such as parallel tempering MCMC. We analysed high and medium mass EMRI systems that fall well inside the low frequency range of LISA. In the context of the Mock LISA Data Challenges, our investigation and results are also the first instance in which a fully Markovian algorithm is applied for EMRI searches. Results show that our algorithm worked well in recovering EMRI signals from different (simulated) LISA data sets having single and multiple EMRI sources and holds great promise for posterior computation under more realistic conditions. The search and estimation methods presented in this paper are general in their nature, and can be applied in any other scenario such as AdLIGO, AdVIRGO and Einstein Telescope with their respective response functions.

preprint2010arXiv

Modelling coloured residual noise in gravitational-wave signal processing

We introduce a signal processing model for signals in non-white noise, where the exact noise spectrum is a priori unknown. The model is based on a Student's t distribution and constitutes a natural generalization of the widely used normal (Gaussian) model. This way, it allows for uncertainty in the noise spectrum, or more generally is also able to accommodate outliers (heavy-tailed noise) in the data. Examples are given pertaining to data from gravitational wave detectors.

preprint2006arXiv

Bayesian inference on compact binary inspiral gravitational radiation signals in interferometric data

Presented is a description of a Markov chain Monte Carlo (MCMC) parameter estimation routine for use with interferometric gravitational radiational data in searches for binary neutron star inspiral signals. Five parameters associated with the inspiral can be estimated, and summary statistics are produced. Advanced MCMC methods were implemented, including importance resampling and prior distributions based on detection probability, in order to increase the efficiency of the code. An example is presented from an application using realistic, albeit fictitious, data.

preprint2004arXiv

Detecting gravitational radiation from neutron stars using a six-parameter adaptive MCMC method

We present a Markov chain Monte Carlo technique for detecting gravitational radiation from a neutron star in laser interferometer data. The algorithm can estimate up to six unknown parameters of the target, including the rotation frequency and frequency derivative, using reparametrization, delayed rejection and simulated annealing. We highlight how a simple extension of the method, distributed over multiple computer processors, will allow for a search over a narrow frequency band. The ultimate goal of this research is to search for sources at a known locations, but uncertain spin parameters, such as may be found in SN1987A.

preprint2001arXiv

II: Bayesian Methods for Cosmological Parameter Estimation from Cosmic Microwave Background Measurements

We present a strategy for a statistically rigorous Bayesian approach to the problem of determining cosmological parameters from the results of observations of anisotropies in the cosmic microwave background. Our strategy relies on Markov chain Monte Carlo methods, specifically the Metropolis-Hastings algorithm, to perform the necessary high-dimensional integrals. We describe the Metropolis-Hastings algorithm in detail and discuss the results of our test on simulated data.