BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Nelson Christensen

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

ResearcherUnclaimed source record
gr-qcastro-ph.IMastro-ph.HEastro-ph.COphysics.data-anastro-ph.SRphysics.geo-phApplicationsastro-phastro-ph.GAComputationMethodologyphysics.optics

Catalog footprint

What is connected

40works
13topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

40 published item(s)

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Impact of correlated seismic and correlated Newtonian noise on the Einstein Telescope

Correlated noise could impact the search for the gravitational wave background at future Earth-based gravitational-wave detectors. Due to the small distance ($\sim$ 400 m) between the different interferometers of the Einstein Telescope, correlated seismic noise could have a significant effect. To this extent, we study the seismic correlations at the Earth's surface, as well as underground, between seismometers and geophones separated by several hundreds of meters, in the frequency range 0.05 Hz - 50 Hz. Based on these correlated seismic fields we predict the levels of correlated Newtonian noise (NN). We construct upper limits on the allowed seismic coupling function such that correlated seismic noise does not affect the search for an isotropic gravitational wave background. Assuming a facility located 300 m below the surface, the impact on the search for a gravitational wave background of correlated NN from Rayleigh waves are found to be problematic up to $\sim$ 5 Hz. The NN from body waves, however, constitutes a serious threat to the search of a gravitational wave background. Correlated NN from body waves could be up to five to seven orders of magnitude above the planned sensitivity at $\sim$ 3 Hz and it could impede any search for a gravitational wave background below 40 Hz. With a factor 10 of NN reduction via NN cancellation in each interferometer, the effects of the NN on the stochastic search could be eliminated above 30 Hz.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Impact of Schumann resonances on the Einstein Telescope and projections for the magnetic coupling function

Correlated magnetic noise in the form of Schumann resonances could introduce limitations to the gravitational-wave background searches of future Earth-based gravitational-wave detectors. We consider recorded magnetic activity at a candidate site for the Einstein Telescope, and forecast the necessary measures to ensure that magnetic contamination will not pose a threat to the science goals of this third-generation detector. In addition to global magnetic effects, we study local magnetic noise and the impact it might have on co-located interferometers. We express our results as upper limits on the coupling function of magnetic fields to the interferometer arms, implying that any larger values of magnetic coupling into the strain channel would lead to a reduction in the detectors' sensitivity. For gravitational-wave background searches below $\sim 30$ Hz it will be necessary for the Einstein Telescope magnetic isolation coupling to be two to four orders of magnitude better than that measured in the current Advanced LIGO and Virgo detectors.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Searches for Modulated γ-Ray Precursors to Compact Binary Mergers in Fermi-GBM Data

GW170817 is the only gravitational-wave (GW) event, for which a confirmed γ-ray counterpart, GRB 170817A, has been detected. Here we present a method to search for another type of γ-ray signal, a γ-ray burst precursor, associated with a compact binary merger. If emitted shortly before the coalescence, a high-energy electromagnetic (EM) flash travels through a highly dynamical and relativistic environment, created by the two compact objects orbiting each other. Thus, the EM signal arriving at an Earth observer could present a somewhat predictable time-dependent modulation. We describe a targeted search method for lightcurves exhibiting such a modulation, parameterized by the observer-frame component masses and binary merger time, using Fermi-GBM data. The sensitivity of the method is assessed based on simulated signals added to GBM data. The method is then applied to a selection of potentially interesting compact binary mergers detected during the second (O2) and third (O3) observing runs of Advanced LIGO and Advanced Virgo. We find no significant modulated γ-ray precursor signal associated with any of the considered events.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Gravitational-Wave Geodesy: Defining False Alarm Probabilities with Respect to Correlated Noise

Future searches for a gravitational-wave background using Earth-based gravitational-wave detectors might be impacted by correlated noise sources. A well known example are the Schumann resonances, which are extensively studied in the context of searches for a gravitational-wave background. Earlier work has shown that a technique termed "gravitational-wave geodesy" can be used to generically differentiate observations of a gravitational-wave background from signals due to correlated terrestrial effects, requiring true observations to be consistent with the known geometry of our detector network. The key result of this test is a Bayes factor between the hypotheses that a candidate signal is astrophysical or terrestrial in origin. Here, we further formalize the geodesy test, mapping distributions of false-alarm and false-acceptance probabilities to quantify the degree with which a given Bayes factor will boost or diminish our confidence in an apparent detection of the gravitational-wave background. To define the false alarm probability of a given Bayes factor, we must have knowledge of our null hypothesis: the space of all possible correlated terrestrial signals. Since we do not have this knowledge we instead construct a generic space of smooth functions in the frequency domain using Gaussian processes, which we tailor to be conservative. This enables us to use draws from our Gaussian processes as a proxy for all possible non-astrophysical signals. As a demonstration, we apply the tool to the SNR = 1.25 excess observed for a 2/3-power law by the LIGO and Virgo collaborations during their second observing run.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Measuring the Hubble Constant with a sample of kilonovae

Kilonovae produced by the coalescence of compact binaries with at least one neutron star are promising standard sirens for an independent measurement of the Hubble constant ($H_0$). Through their detection via follow-up of gravitational-wave (GW), short gamma-ray bursts (sGRBs) or optical surveys, a large sample of kilonovae (even without GW data) can be used for $H_0$ contraints. Here, we show measurement of $H_0$ using light curves associated with four sGRBs, assuming these are attributable to kilonovae, combined with GW170817. Including a systematic uncertainty on the models that is as large as the statistical ones, we find $H_0 = 73.8^{+6.3}_{-5.8}$\,$\mathrm{km}$ $\mathrm{s}^{-1}$ $\mathrm{Mpc}^{-1}$ and $H_0 = 71.2^{+3.2}_{-3.1}$\,$\mathrm{km}$ $\mathrm{s}^{-1}$ $\mathrm{Mpc}^{-1}$ for two different kilonova models that are consistent with the local and inverse-distance ladder measurements. For a given model, this measurement is about a factor of 2-3 more precise than the standard-siren measurement for GW170817 using only GWs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Numerical solutions for phase noise due to pointing jitter with the LISA telescope

The aim of the Laser Interferometer Space Antenna (LISA) is to detect gravitational waves through a phase modulation in long (2.5 Mkm) laser light links between spacecraft. Among other noise sources to be addressed are the phase fluctuations caused by a possible angular jitter of the emitted beam. The present paper follows our preceding one (Vinet et al 2019 Class. Quant. Grav. 36, 205 003) based on an analytical study of the far field phase. We address here a numerical treatment of the phase, to first order in the emitted wavefront aberrations, but without any assumptions on the static bias term. We verify that, in the phase change, the higher order terms in the static mispointing are consistent with the results found in our preceding paper.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On standardizing kilonovae and their use as standard candles to measure the Hubble constant

The detection of GW170817 is revolutionizing many areas of astrophysics with the joint observation of gravitational waves and electromagnetic emissions. These multi-messenger events provide a new approach to determine the Hubble constant, thus, they are a promising candidate for mitigating the tension between measurements of Type Ia supernovae via the local distance ladder and the Cosmic Microwave Background. In addition to the "standard siren" provided by the gravitational-wave measurement, the kilonova itself has characteristics that allow to improve existing measurements or to perform yet another, independent measurement of the Hubble constant without gravitational-wave information. Here, we employ standardization techniques borrowed from the type-Ia community and apply them to kilonovae, not using any information from the gravitational-wave signal. We use two versions of this technique, one derived from direct observables measured from the lightcurve, and the other based on inferred ejecta parameters, e.g., mass, velocity, and composition, for two different models. These lead to Hubble Constant measurements of $H_0 = 109^{+49}_{-35}$\,km $\mathrm{s}^{-1}$ $\mathrm{Mpc}^{-1}$ for the measured analysis, and $H_0 = 85^{+22}_{-17}$\,km $\mathrm{s}^{-1}$ $\mathrm{Mpc}^{-1}$ and $H_0 = 79^{+23}_{-15}$\,km $\mathrm{s}^{-1}$ $\mathrm{Mpc}^{-1}$ for the inferred analyses. This measurement has error bars within ~$\sim$\,2 to the gravitational-wave measurements ($H_0=74^{+16}_{-8}$\,km $\mathrm{s}^{-1}$ $\mathrm{Mpc}^{-1}$), showing its promise as an independent constraint on $H_0$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Simultaneous estimation of astrophysical and cosmological stochastic gravitational-wave backgrounds with terrestrial detectors

The recent Advanced LIGO and Advanced Virgo joint observing runs have not claimed a stochastic gravitational-wave background detection, but one expects this to change as the sensitivity of the detectors improves. The challenge of claiming a true detection will be immediately succeeded by the difficulty of relating the signal to the sources that contribute to it. In this paper, we consider backgrounds that comprise compact binary coalescences and additional cosmological sources, and we set simultaneous upper limits on these backgrounds. We find that the Advanced LIGO, Advanced Virgo network, operating at design sensitivity, will not allow for separation of the sources we consider. Third generation detectors, sensitive to most individual compact binary mergers, can reduce the astrophysical signal via subtraction of individual sources, and potentially reveal a cosmological background. Our Bayesian analysis shows that, assuming a detector network containing Cosmic Explorer and Einstein Telescope and reasonable levels of individual source subtraction, we can detect cosmological signals $Ω_{\rm{CS}} (25\,\rm{Hz})=4.5 \times 10^{-13}$ for cosmic strings, and $Ω_{\rm BPL}(25\,\rm{Hz})= 2.2 \times 10^{-13}$ for a broken power law model of an early universe phase transition.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Using machine learning for transient classification in searches for gravitational-wave counterparts

The large sky localization regions offered by the gravitational-wave interferometers require efficient follow-up of the many counterpart candidates identified by the wide field-of-view telescopes. Given the restricted telescope time, the creation of prioritized lists of the many identified candidates becomes mandatory. Towards this end, we use \text{\astrorapid}, a multi-band photometric lightcurve classifier, to differentiate between kilonovae, supernovae, and other possible transients. We demonstrate our method on the photometric observations of real events. In addtion, the classification performance is tested on simulated lightcurves, both ideally and realistically sampled. We show that after only a few days of observations of an astronomical object, it is possible to rule out candidates as supernovae and other known transients.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Mock Data and Science Challenge for Detecting an Astrophysical Stochastic Gravitational-Wave Background with Advanced LIGO and Advanced Virgo

The purpose of this mock data and science challenge is to prepare the data analysis and science interpretation for the second generation of gravitational-wave experiments Advanced LIGO-Virgo in the search for a stochastic gravitational-wave background signal of astrophysical origin. Here we present a series of signal and data challenges, with increasing complexity, whose aim is to test the ability of current data analysis pipelines at detecting an astrophysically produced gravitational-wave background, test parameter estimation methods and interpret the results. We introduce the production of these mock data sets that includes a realistic observing scenario data set where we account for different sensitivities of the advanced detectors as they are continuously upgraded toward their design sensitivity. After analysing these with the standard isotropic cross-correlation pipeline we find that we are able to recover the injected gravitational-wave background energy density to within $2σ$ for all of the data sets and present the results from the parameter estimation. The results from this mock data and science challenge show that advanced LIGO and Virgo will be ready and able to make a detection of an astrophysical gravitational-wave background within a few years of operations of the advanced detectors, given a high enough rate of compact binary coalescing events.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Proposed Search for the Detection of Gravitational Waves from Eccentric Binary Black Holes

Most of compact binary systems are expected to circularize before the frequency of emitted gravitational waves (GWs) enters the sensitivity band of the ground based interferometric detectors. However, several mechanisms have been proposed for the formation of binary systems, which retain eccentricity throughout their lifetimes. Since no matched-filtering algorithm has been developed to extract continuous GW signals from compact binaries on orbits with low to moderate values of eccentricity, and available algorithms to detect binaries on quasi-circular orbits are sub-optimal to recover these events, in this paper we propose a search method for detection of gravitational waves produced from the coalescences of eccentric binary black holes (eBBH). We study the search sensitivity and the false alarm rates on a segment of data from the second joint science run of LIGO and Virgo detectors, and discuss the implications of the eccentric binary search for the advanced GW detectors.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

All-sky, narrowband, gravitational-wave radiometry with folded data

Gravitational-wave radiometry is a powerful tool by which weak signals with unknown signal morphologies are recovered through a process of cross correlation. Radiometry has been used, e.g., to search for persistent signals from known neutron stars such as Scorpius X-1. In this paper, we demonstrate how a more ambitious search--for persistent signals from unknown neutron stars--can be efficiently carried out using folded data, in which an entire ~year-long observing run is represented as a single sidereal day. The all-sky, narrowband radiometer search described here will provide a computationally tractable means to uncover gravitational-wave signals from unknown, nearby neutron stars in binary systems, which can have modulation depths of ~0.1-2 Hz. It will simultaneously provide a sensitive search algorithm for other persistent, narrowband signals from unexpected sources.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Chasing 5-sigma: Prospects for searches for long-duration gravitational-waves without time slides

The detection of unmodeled gravitational-wave bursts by ground-based interferometric gravitational-wave detectors is a major goal for the advanced detector era. These searches are commonly cast as pattern recognition problems, where the goal is to identify statistically significant clusters in spectrograms of strain power when the precise signal morphology is unknown. In previous work, we have introduced a clustering algorithm referred to as "seedless clustering," and shown that it is a powerful tool for detecting weak long-lived (10-1000s) signals in background. However, as the algorithm is currently conceived, in order to carry out an all-sky search on a $\approx$ year of data, significant computational resources may be required in order to carry out background estimation. Alternatively, some of the sensitivity of the search must be sacrificed to control computational costs. The sensitivity of the algorithm is limited by the amount of computing resources due to the requirement of performing background studies to assign significance in gravitational-wave searches. In this paper, we present an analytic method for estimating the background generated by the seedless clustering algorithm and compare the performance to both Monte Carlo Gaussian noise and time-shifted gravitational-wave data from a week of LIGO's 5th Science Run. We demonstrate qualitative agreement between the model and measured distributions and argue that the approximation will be useful to supplement conventional background estimation techniques for advanced detector searches for long-duration gravitational-wave transients.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Detecting very long-lived gravitational-wave transients lasting hours to weeks

We explore the possibility of very long-lived gravitational-wave transients (and detector artifacts) lasting hours to weeks. Such very long signals are both interesting in their own right and as a potential source of systematic error in searches for persistent signals, e.g., from a stochastic gravitational-wave background. We review possible mechanisms for emission on these time scales and discuss computational challenges associated with their detection: namely, the substantial volume of data involved in a search for very long transients can require vast computer memory and processing time. These computational difficulties can be addressed through a form of data compression known as coarse-graining, in which information about short time spans is discarded in order to reduce the computational requirements of a search. Using data compression, we demonstrate an efficient radiometer (cross-correlation) algorithm for the detection of very long transients. In the process, we identify features of a very long transient search (related to the rotation of the Earth) that make it more complicated than a search for shorter transient signals. We implement suitable solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

High-energy electromagnetic offline follow-up of LIGO-Virgo gravitational-wave binary coalescence candidate events

We present two different search methods for electromagnetic counterparts to gravitational-wave (GW) events from ground-based detectors using archival NASA high-energy data from the Fermi-GBM and RXTE-ASM instruments. To demonstrate the methods, we use a limited number of representative GW background noise events produced by a search for binary neutron star coalescence over the last two months of the LIGO-Virgo S6/VSR3 joint science run. Time and sky location provided by the GW data trigger a targeted search in the high-energy photon data. We use two custom pipelines: one to search for prompt gamma-ray counterparts in GBM, and the other to search for a variety of X-ray afterglow model signals in ASM. We measure the efficiency of the joint pipelines to weak gamma-ray burst counterparts, and a family of model X-ray afterglows. By requiring a detectable signal in either electromagnetic instrument coincident with a GW event, we are able to reject a large majority of GW candidates. This reduces the signal-to-noise of the loudest surviving GW background event by around 15-20%.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Improved constraint on the primordial gravitational-wave density using recent cosmological data and its impact on cosmic string models

The production of a primordial stochastic gravitational-wave background by processes occuring in the early Universe is expected in a broad range of models. Observing this background would open a unique window onto the Universe's evolutionary history. Probes like the Cosmic Microwave Background (CMB) or the Baryon Acoustic Oscillations (BAO) can be used to set upper limits on the stochastic gravitational-wave background energy density $Ω_{GW}$ for frequencies above $10^{-15}$ Hz. We perform a profile likelihood analysis of the Planck CMB temperature anisotropies and gravitational lensing data combined with WMAP low-$\ell$ polarization, BAO, South Pole Telescope and Atacama Cosmology Telescope data. We find that $Ω_{GW}h_{0}^{2} < 3.8 \times 10^{-6}$ at 95\% confidence level for adiabatic initial conditions which improves over the previous limit by a factor 2.3. Assuming that the primordial gravitational waves have been produced by a network of cosmic strings, we have derived exclusion limits in the cosmic string parameter space. If the size of the loops is determined by gravitational back-reaction, string tension values greater than $\sim 4 \times 10^{-9}$ are excluded for a reconnection probability of $10^{-3}$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Planned search for LIGO-GBM coincidence in the first advanced LIGO data run

In the fall of 2015 the first scientific observing run (O1) of the advanced LIGO detectors will be conducted. Based on the recent commissioning progress at the LIGO Hanford and Livingston sites, the gravitational wave detector range for a neutron star binary inspiral is expected to be of order 60 Mpc. We describe here our planning for an O1 search for coincidence between a LIGO gravitational wave detection and a gamma-ray signal from the Fermi Gamma-ray Burst Monitor. Such a coincidence would constitute measurement of an electromagnetic counterpart to a gravitational wave signal, with significant corresponding scientific benefits, including revealing the central engine powering the gamma-ray burst, enhanced confidence in the event as a genuine astrophysical detection, and a determination of the relative speed of the photon and graviton.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Robust parameter estimation for compact binaries with ground-based gravitational-wave observations using the LALInference software library

The Advanced LIGO and Advanced Virgo gravitational wave (GW) detectors will begin operation in the coming years, with compact binary coalescence events a likely source for the first detections. The gravitational waveforms emitted directly encode information about the sources, including the masses and spins of the compact objects. Recovering the physical parameters of the sources from the GW observations is a key analysis task. This work describes the LALInference software library for Bayesian parameter estimation of compact binary signals, which builds on several previous methods to provide a well-tested toolkit which has already been used for several studies. We show that our implementation is able to correctly recover the parameters of compact binary signals from simulated data from the advanced GW detectors. We demonstrate this with a detailed comparison on three compact binary systems: a binary neutron star, a neutron star black hole binary and a binary black hole, where we show a cross-comparison of results obtained using three independent sampling algorithms. These systems were analysed with non-spinning, aligned spin and generic spin configurations respectively, showing that consistent results can be obtained even with the full 15-dimensional parameter space of the generic spin configurations. We also demonstrate statistically that the Bayesian credible intervals we recover correspond to frequentist confidence intervals under correct prior assumptions by analysing a set of 100 signals drawn from the prior. We discuss the computational cost of these algorithms, and describe the general and problem-specific sampling techniques we have used to improve the efficiency of sampling the compact binary coalescence parameter space.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Correlated noise in networks of gravitational-wave detectors: subtraction and mitigation

One of the key science goals of advanced gravitational-wave detectors is to observe a stochastic gravitational-wave background. However, recent work demonstrates that correlated magnetic fields from Schumann resonances can produce correlated strain noise over global distances, potentially limiting the sensitivity of stochastic background searches with advanced detectors. In this paper, we estimate the correlated noise budget for the worldwide Advanced LIGO network and conclude that correlated noise may affect upcoming measurements. We investigate the possibility of a Wiener filtering scheme to subtract correlated noise from Advanced LIGO searches, and estimate the required specifications. We also consider the possibility that residual correlated noise remains following subtraction, and we devise an optimal strategy for measuring astronomical parameters in the presence of correlated noise. Using this new formalism, we estimate the loss of sensitivity for a broadband, isotropic stochastic background search using 1 yr of LIGO data at design sensitivity. Given our current noise budget, the uncertainty with which LIGO can estimate energy density will likely increase by a factor of ~4--if it is impossible to achieve significant subtraction. Additionally, narrowband cross-correlation searches may be severely affected at low frequencies f < 45 Hz without effective subtraction.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Detecting compact binary coalescences with seedless clustering

Compact binary coalescences are a promising source of gravitational waves for second-generation interferometric gravitational-wave detectors. Although matched filtering is the optimal search method for well-modeled systems, alternative detection strategies can be used to guard against theoretical errors (e.g., involving new physics and/or assumptions about spin/eccentricity) while providing a measure of redundancy. In previous work, we showed how "seedless clustering" can be used to detect long-lived gravitational-wave transients in both targeted and all-sky searches. In this paper, we apply seedless clustering to the problem of low-mass ($M_\text{total}\leq10M_\odot$) compact binary coalescences for both spinning and eccentric systems. We show that seedless clustering provides a robust and computationally efficient method for detecting low-mass compact binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Instrumental vetoes for transient gravitational-wave triggers using noise-coupling models: The bilinear-coupling veto

LIGO and Virgo recently completed searches for gravitational waves at their initial target sensitivities, and soon Advanced LIGO and Advanced Virgo will commence observations with even better capabilities. In the search for short duration signals, such as coalescing compact binary inspirals or "burst" events, noise transients can be problematic. Interferometric gravitational-wave detectors are highly complex instruments, and, based on the experience from the past, the data often contain a large number of noise transients that are not easily distinguishable from possible gravitational-wave signals. In order to perform a sensitive search for short-duration gravitational-wave signals it is important to identify these noise artifacts, and to "veto" them. Here we describe such a veto, the bilinear-coupling veto, that makes use of an empirical model of the coupling of instrumental noise to the output strain channel of the interferometric gravitational-wave detector. In this method, we check whether the data from the output strain channel at the time of an apparent signal is consistent with the data from a bilinear combination of auxiliary channels. We discuss the results of the application of this veto on recent LIGO data, and its possible utility when used with data from Advanced LIGO and Advanced Virgo.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Method for estimation of gravitational-wave transient model parameters in frequency-time maps

A common technique for detection of gravitational-wave signals is searching for excess power in frequency-time maps of gravitational-wave detector data. In the event of a detection, model selection and parameter estimation will be performed in order to explore the properties of the source. In this paper, we develop a Bayesian statistical method for extracting model-dependent parameters from observed gravitational-wave signals in frequency-time maps. We demonstrate the method by recovering the parameters of model gravitational-wave signals added to simulated advanced LIGO noise. We also characterize the performance of the method and discuss prospects for future work.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The detectability of eccentric compact binary coalescences with advanced gravitational-wave detectors

Compact binary coalescences are a promising source of gravitational waves for second-generation interferometric gravitational-wave detectors such as advanced LIGO and advanced Virgo. While most binaries are expected to possess circular orbits, some may be eccentric, for example, if they are formed through dynamical capture. Eccentric orbits can create difficulty for matched filtering searches due to the challenges of creating effective template banks to detect these signals. In previous work, we showed how seedless clustering can be used to detect low-mass ($M_\text{total}\leq10M_\odot$) compact binary coalescences for both spinning and eccentric systems, assuming a circular post-Newtonian expansion. Here, we describe a parameterization that is designed to maximize sensitivity to low-eccentricity ($0\leqε\leq0.6$) systems, derived from the analytic equations. We show that this parameterization provides a robust and computationally efficient method for detecting eccentric low-mass compact binaries. Based on these results, we conclude that advanced detectors will have a chance of detecting eccentric binaries if optimistic models prove true. However, a null observation is unlikely to firmly rule out models of eccentric binary populations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Wiener filtering with a seismic underground array at the Sanford Underground Research Facility

A seismic array has been deployed at the Sanford Underground Research Facility in the former Homestake mine, South Dakota, to study the underground seismic environment. This includes exploring the advantages of constructing a third-generation gravitational-wave detector underground. A major noise source for these detectors would be Newtonian noise, which is induced by fluctuations in the local gravitational field. The hope is that a combination of a low-noise seismic environment and coherent noise subtraction using seismometers in the vicinity of the detector could suppress the Newtonian noise to below the projected noise floor for future gravitational-wave detectors. In this paper, we use Wiener filtering techniques to subtract coherent noise in a seismic array in the frequency band 0.05 -- 1\,Hz. This achieves more than an order of magnitude noise cancellation over a majority of this band. We show how this subtraction would benefit proposed future low-frequency gravitational wave detectors. The variation in the Wiener filter coefficients over the course of the day, including how local activities impact the filter, is analyzed. We also study the variation in coefficients over the course of a month, showing the stability of the filter with time. How varying the filter order affects the subtraction performance is also explored. It is shown that optimizing filter order can significantly improve subtraction of seismic noise, which gives hope for future gravitational-wave detectors to address Newtonian noise.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Correlated magnetic noise in global networks of gravitational-wave interferometers: observations and implications

One of the most ambitious goals of gravitational-wave astronomy is to observe the stochastic gravitational-wave background. Correlated noise in two or more detectors can introduce a systematic error, which limits the sensitivity of stochastic searches. We report on measurements of correlated magnetic noise from Schumann resonances at the widely separated LIGO and Virgo detectors. We investigate the effect of this noise on a global network of interferometers and derive a constraint on the allowable coupling of environmental magnetic fields to test mass motion in gravitational-wave detectors. We find that while correlated noise from global electromagnetic fields could be safely ignored for initial LIGO stochastic searches, it could severely impact Advanced LIGO and third-generation detectors.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Reconstructing the sky location of gravitational-wave detected compact binary systems: methodology for testing and comparison

The problem of reconstructing the sky position of compact binary coalescences detected via gravitational waves is a central one for future observations with the ground-based network of gravitational-wave laser interferometers, such as Advanced LIGO and Advanced Virgo. Different techniques for sky localisation have been independently developed. They can be divided in two broad categories: fully coherent Bayesian techniques, which are high-latency and aimed at in-depth studies of all the parameters of a source, including sky position, and "triangulation-based" techniques, which exploit the data products from the search stage of the analysis to provide an almost real-time approximation of the posterior probability density function of the sky location of a detection candidate. These techniques have previously been applied to data collected during the last science runs of gravitational-wave detectors operating in the so-called initial configuration. Here, we develop and analyse methods for assessing the self-consistency of parameter estimation methods and carrying out fair comparisons between different algorithms, addressing issues of efficiency and optimality. These methods are general, and can be applied to parameter estimation problems other than sky localisation. We apply these methods to two existing sky localisation techniques representing the two above-mentioned categories, using a set of simulated inspiral-only signals from compact binary systems with total mass $\le 20\,M_\odot$ and non-spinning components. We compare the relative advantages and costs of the two techniques and show that sky location uncertainties are on average a factor $\approx 20$ smaller for fully coherent techniques than for the specific variant of the "triangulation-based" technique used during the last science runs, at the expense of a factor $\approx 1000$ longer processing time.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Techniques for targeted Fermi-GBM follow-up of gravitational-wave events

The Advanced LIGO and Advanced Virgo ground-based gravitational-wave detectors are projected to come online 2015-2016, reaching a final sensitivity sufficient to observe dozens of binary neutron star mergers per year by 2018. We present a fully-automated, targeted search strategy for prompt gamma-ray counterparts in offline Fermi-GBM data. The multi-detector method makes use of a detailed model response of the instrument, and benefits from time and sky location information derived from the gravitational-wave signal.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Characterization of the seismic environment at the Sanford Underground Laboratory, South Dakota

An array of seismometers is being developed at the Sanford Underground Laboratory, the former Homestake mine, in South Dakota to study the properties of underground seismic fields and Newtonian noise, and to investigate the possible advantages of constructing a third-generation gravitational-wave detector underground. Seismic data were analyzed to characterize seismic noise and disturbances. External databases were used to identify sources of seismic waves: ocean-wave data to identify sources of oceanic microseisms, and surface wind-speed data to investigate correlations with seismic motion as a function of depth. In addition, sources of events contributing to the spectrum at higher frequencies are characterized by studying the variation of event rates over the course of a day. Long-term observations of spectral variations provide further insight into the nature of seismic sources. Seismic spectra at three different depths are compared, establishing the 4100-ft level as a world-class low seismic-noise environment.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Status of NINJA: the Numerical INJection Analysis project

The 2008 NRDA conference introduced the Numerical INJection Analysis project (NINJA), a new collaborative effort between the numerical relativity community and the data analysis community. NINJA focuses on modeling and searching for gravitational wave signatures from the coalescence of binary system of compact objects. We review the scope of this collaboration and the components of the first NINJA project, where numerical relativity groups shared waveforms and data analysis teams applied various techniques to detect them when embedded in colored Gaussian noise.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Testing gravitational-wave searches with numerical relativity waveforms: Results from the first Numerical INJection Analysis (NINJA) project

The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave data analysis communities. The purpose of NINJA is to study the sensitivity of existing gravitational-wave search algorithms using numerically generated waveforms and to foster closer collaboration between the numerical relativity and data analysis communities. We describe the results of the first NINJA analysis which focused on gravitational waveforms from binary black hole coalescence. Ten numerical relativity groups contributed numerical data which were used to generate a set of gravitational-wave signals. These signals were injected into a simulated data set, designed to mimic the response of the Initial LIGO and Virgo gravitational-wave detectors. Nine groups analysed this data using search and parameter-estimation pipelines. Matched filter algorithms, un-modelled-burst searches and Bayesian parameter-estimation and model-selection algorithms were applied to the data. We report the efficiency of these search methods in detecting the numerical waveforms and measuring their parameters. We describe preliminary comparisons between the different search methods and suggest improvements for future NINJA analyses.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access