Source author record

Michael Coughlin

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

29works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

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

29 published item(s)

preprint2022arXiv

Physics Community Needs, Tools, and Resources for Machine Learning

Machine learning (ML) is becoming an increasingly important component of cutting-edge physics research, but its computational requirements present significant challenges. In this white paper, we discuss the needs of the physics community regarding ML across latency and throughput regimes, the tools and resources that offer the possibility of addressing these needs, and how these can be best utilized and accessed in the coming years.

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.

preprint2021arXiv

Initial Characterization of Active Transitioning Centaur, P/2019 LD2 (ATLAS), using Hubble, Spitzer, ZTF, Keck, APO and GROWTH Visible & Infrared Imaging and Spectroscopy

We present visible and mid-infrared imagery and photometry of temporary Jovian co-orbital comet P/2019 LD$_2$ taken with HST/WFC3, Spitzer/IRAC, the GROWTH telescope network, visible spectroscopy from Keck/LRIS and archival ZTF observations taken between 2019 April and 2020 August. Our observations indicate that the nucleus of LD$_2$ has a radius between 0.2-1.8 km assuming a 0.08 albedo and a coma dominated by $\sim$100$μ$ m-scale dust ejected at $\sim$1 m/s speeds with a $\sim$1'' jet pointing in the SW direction. LD$_2$ experienced a total dust mass loss of $\sim$10$^8$ kg at a loss rate of $\sim$6 kg/s with Af$ρ$/cross-section varying between $\sim$85 cm/125 km$^2$ and $\sim$200 cm/310 km$^2$ from 2019 April 9 to 2019 Nov 8. If the increase in Af$ρ$/cross-section remained constant, it implies LD$_2$'s activity began $\sim$2018 November when within 4.8 au of the Sun, implying the onset of H$_2$O sublimation. We measure CO/CO$_2$ gas production of $\lesssim$10$^{27}$ mol/s /$\lesssim$10$^{26}$ mol/s from our 4.5 $μ$m Spitzer observations, $g$-$r$ = 0.59$\pm$0.03, $r$-$i$ = 0.18$\pm$0.05, $i$-$z$ = 0.01$\pm$0.07 from GROWTH observations, H$_2$O gas production of $\lesssim$80 kg/s scaling from our estimated $C_2$ production of $Q_{C_2}\lesssim$7.5$\times10^{24}$ mol/s from Keck/LRIS spectroscopy. We determine that the long-term orbit of LD$_2$ is similar to Jupiter family comets having close encounters with Jupiter within $\sim$0.5 Hill radius in the last $\sim$3 y, within 0.8 Hill radius in $\sim$9 y. Additionally, 78.8$\%$ of our orbital clones are ejected from the Solar System within $1 \times 10^{6}$ years having a dynamical half-life of 3.4 $\times 10^5$ years.

preprint2020arXiv

Cataclysmic Variables in the First Year of the Zwicky Transient Facility

Using selection criteria based on amplitude, time and color, we have identified 329 objects as known or candidate cataclysmic variable (CVs) during the first year of testing and operation of the Zwicky Transient Facility (ZTF). Of these, 90 are previously confirmed CVs, 218 are strong candidates based on the shape and color of their light curves obtained during 3-562 days of observations, and the remaining 21 are possible CVs but with too few data points to be listed as good candidates. Almost half the strong candidates are within 10 deg of the galactic plane, in contrast to most other large surveys which have avoided crowded fields. The available Gaia parallaxes are consistent with sampling the low mass transfer CVs, as predicted by population models. Our followup spectra have confirmed Balmer/helium emission lines in 27 objects, with four showing high excitation HeII emission, including candidates for an AM CVn, a polar and an intermediate polar. Our results demonstrate that a complete survey of the galactic plane is needed to accomplish an accurate determination of the number of CVs existing in the Milky Way.

preprint2020arXiv

Characterization of Temporarily-Captured Minimoon 2020 CD$_3$ by Keck Time-resolved Spectrophotometry

We present time-resolved visible spectrophotometry of minimoon 2020 CD$_3$, the second asteroid known to become temporarily captured by the Earth-Moon system's gravitational field. The spectrophotometry was taken with Keck I/LRIS between wavelengths 434 nm and 912 nm in $B$, $g$, $V$, $R$, $I$ and RG850 filters as it was leaving the Earth-Moon system on 2020 March 23 UTC. The spectrophotometry of 2020 CD$_3$ most closely resembles the spectra of V-type asteroids and some Lunar rock samples with a reddish slope of ~18$\%$/100 nm between 434 nm and 761 nm corresponding to colors of $g$-$r$ = 0.62$\pm$0.08, $r$-$i$ = 0.21 $\pm$ 0.06 and an absorption band at ~900 nm corresponding to $i$-$z$ = -0.54$\pm$0.10. Combining our measured 31.9$\pm$0.1 absolute magnitude with an albedo of 0.35 typical for V-type asteroids, we determine 2020 CD$_3$'s diameter to be ~0.9$\pm$0.1 m making it the first minimoon and one of the smallest asteroids to be spectrally studied. We use our time-series photometry to detect periodic lightcurve variations with a $<$10$^{-4}$ false alarm probability corresponding to a lightcurve period of ~573 s and a lightcurve amplitude of ~1 mag implying 2020 CD$_3$ possesses a $b/a$ axial ratio of ~2.5. In addition, we extend the observational arc of 2020 CD$_3$ to 37 days between 2020 February 15 UTC and 2020 March 23 UTC. From the improved orbital solution for 2020 CD$_3$, we estimate its likely duration of its capture to be ~2 y, and we measure the non-gravitation perturbation on its orbit due to radiation pressure with an area-to-mass ratio of 6.9$\pm$2.4$\times$10$^{-4}$ m$^2$/kg implying a density of 2.3$\pm$0.8 g/cm$^3$, broadly compatible with the densities of other meter-scale asteroids and Lunar rock. We searched for pre-discovery detections of 2020 CD$_3$ in the ZTF archive as far back as 2018 October, but were unable to locate any positive detections.

preprint2020arXiv

Characterization of the Nucleus, Morphology and Activity of Interstellar Comet 2I/Borisov by Optical and Near-Infrared GROWTH, Apache Point, IRTF, ZTF and Keck Observations

We present visible and near-infrared photometric and spectroscopic observations of interstellar object 2I/Borisov taken from 2019 September 10 to 2019 November 29 using the GROWTH, the APO ARC 3.5 m and the NASA/IRTF 3.0 m combined with post and pre-discovery observations of 2I obtained by ZTF from 2019 March 17 to 2019 May 5. Comparison with imaging of distant Solar System comets shows an object very similar to mildly active Solar System comets with an out-gassing rate of $\sim$10$^{27}$ mol/sec. The photometry, taken in filters spanning the visible and NIR range shows a gradual brightening trend of $\sim0.03$ mags/day since 2019 September 10 UTC for a reddish object becoming neutral in the NIR. The lightcurve from recent and pre-discovery data reveals a brightness trend suggesting the recent onset of significant H$_2$O sublimation with the comet being active with super volatiles such as CO at heliocentric distances $>$6 au consistent with its extended morphology. Using the advanced capability to significantly reduce the scattered light from the coma enabled by high-resolution NIR images from Keck adaptive optics taken on 2019 October 04, we estimate a diameter of 2I's nucleus of $\lesssim$1.4 km. We use the size estimates of 1I/'Oumuamua and 2I/Borisov to roughly estimate the slope of the ISO size-distribution resulting in a slope of $\sim$3.4$\pm$1.2, similar to Solar System comets and bodies produced from collisional equilibrium.

preprint2020arXiv

Enhancing Gravitational-Wave Science with Machine Learning

Machine learning has emerged as a popular and powerful approach for solving problems in astrophysics. We review applications of machine learning techniques for the analysis of ground-based gravitational-wave detector data. Examples include techniques for improving the sensitivity of Advanced LIGO and Advanced Virgo gravitational-wave searches, methods for fast measurements of the astrophysical parameters of gravitational-wave sources, and algorithms for reduction and characterization of non-astrophysical detector noise. These applications demonstrate how machine learning techniques may be harnessed to enhance the science that is possible with current and future gravitational-wave detectors.

preprint2020arXiv

Noise Reduction in Gravitational-wave Data via Deep Learning

With the advent of gravitational wave astronomy, techniques to extend the reach of gravitational wave detectors are desired. In addition to the stellar-mass black hole and neutron star mergers already detected, many more are below the surface of the noise, available for detection if the noise is reduced enough. Our method (DeepClean) applies machine learning algorithms to gravitational wave detector data and data from on-site sensors monitoring the instrument to reduce the noise in the time-series due to instrumental artifacts and environmental contamination. This framework is generic enough to subtract linear, non-linear, and non-stationary coupling mechanisms. It may also provide handles in learning about the mechanisms which are not currently understood to be limiting detector sensitivities. The robustness of the noise reduction technique in its ability to efficiently remove noise with no unintended effects on gravitational-wave signals is also addressed through software signal injection and parameter estimation of the recovered signal. It is shown that the optimal SNR ratio of the injected signal is enhanced by $\sim 21.6\%$ and the recovered parameters are consistent with the injected set. We present the performance of this algorithm on linear and non-linear noise sources and discuss its impact on astrophysical searches by gravitational wave detectors.

preprint2020arXiv

Observation of a potential future sensitivity limitation from ground motion at LIGO Hanford

A first detection of terrestrial gravity noise in gravitational-wave detectors is a formidable challenge. With the help of environmental sensors, it can in principle be achieved before the noise becomes dominant by estimating correlations between environmental sensors and the detector. The main complication is to disentangle different coupling mechanisms between the environment and the detector. In this paper, we analyze the relations between physical couplings and correlations that involve ground motion and LIGO strain data h(t) recorded during its second science run in 2016 and 2017. We find that all noise correlated with ground motion was more than an order of magnitude lower than dominant low-frequency instrument noise, and the dominant coupling over part of the spectrum between ground and h(t) was residual coupling through the seismic-isolation system. We also present the most accurate gravitational coupling model so far based on a detailed analysis of data from a seismic array. Despite our best efforts, we were not able to unambiguously identify gravitational coupling in the data, but our improved models confirm previous predictions that gravitational coupling might already dominate linear ground-to-h(t) coupling over parts of the low-frequency, gravitational-wave observation band.

preprint2016arXiv

Towards a first design of a Newtonian-noise cancellation system for Advanced LIGO

Newtonian gravitational noise from seismic fields is predicted to be a limiting noise source at low frequency for second generation gravitational-wave detectors. Mitigation of this noise will be achieved by Wiener filtering using arrays of seismometers deployed in the vicinity of all test masses. In this work, we present optimized configurations of seismometer arrays using a variety of simplified models of the seismic field based on seismic observations at LIGO Hanford. The model that best fits the seismic measurements leads to noise reduction limited predominantly by seismometer self-noise. A first simplified design of seismic arrays for Newtonian-noise cancellation at the LIGO sites is presented, which suggests that it will be sufficient to monitor surface displacement inside the buildings.

preprint2015arXiv

A Daytime Measurement of the Lunar Contribution to the Night Sky Brightness in LSST's ugrizy Bands-- Initial Results

We report measurements from which we determine the spatial structure of the lunar contribution to night sky brightness, taken at the LSST site on Cerro Pachon in Chile. We use an array of six photodiodes with filters that approximate the Large Synoptic Survey Telescope's {\it u, g, r, i, z,} and {\it y} bands. We use the sun as a proxy for the moon, and measure sky brightness as a function of zenith angle of the point on sky, zenith angle of the sun, and angular distance between the sun and the point on sky. We make a correction for the difference between the illumination spectrum of the sun and the moon. Since scattered sunlight totally dominates the daytime sky brightness, this technique allows us to cleanly determine the contribution to the (cloudless) night sky from backscattered moonlight, without contamination from other sources of night sky brightness. We estimate our uncertainty in the relative lunar night sky brightness vs. zenith and lunar angle to be 10\,\%. This information is useful in planning the optimal execution of the LSST survey, and perhaps for other astronomical observations as well. Although our primary objective is to map out the angular structure and spectrum of the scattered light from the atmosphere and particulates, we also make an estimate of the expected number of scattered lunar photons per pixel per second in LSST, and find values that are in overall agreement with previous estimates.

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.

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.

preprint2015arXiv

Constraining the r-mode saturation amplitude from a hypothetical detection of r-mode gravitational waves from a newborn neutron star - sensitivity study

This paper consists of two related parts: In the first part we derive an expression of the moment of inertia (MOI) of a neutron star as a function of observables from a hypothetical r-mode gravitational wave detection. For a given r-mode detection we show how the value of the MOI of a neutron star constrains the equation of state (EOS) of the matter in the core of the neutron star. Subsequently, for each candidate EOS, we derive a possible value of the saturation amplitude, α, of the r-mode oscillations on the neutron star. Additionally, we argue that a r-mode detection will provide clues about the cooling rate mechanism of the neutron star. The above physics that can be derived from a hypothetical r-mode detection constitute our motivation for the second part of the paper. In that part we present a detection strategy to efficiently search for r-modes in gravitational-wave data. R-mode signals were injected into simulated noise colored with the advanced LIGO (aLIGO) and Einstein Telescope (ET) sensitivity curves. The r-mode waveforms used are those predicted by early theories based on a polytropic equation of state (EOS) neutron star matter. In our best case scenario αof order 10^{-1}, the maximum detection distance when using the aLIGO sensitivity curve is 1 Mpc (supernova event rate of 3-4 per century) while the maximum detection distance when using the ET sensitivity curve is 10 Mpc (supernova event rate of 1-2 per year).

preprint2015arXiv

Detecting gravitational-wave transients at five sigma: a hierarchical approach

As second-generation gravitational-wave detectors prepare to analyze data at unprecedented sensitivity, there is great interest in searches for unmodeled transients, commonly called bursts. Significant effort has yielded a variety of techniques to identify and characterize such transient signals, and many of these methods have been applied to produce astrophysical results using data from first-generation detectors. However, the computational cost of background estimation remains a challenging problem; it is difficult to claim a 5σ detection with reasonable computational resources without paying for efficiency with reduced sensitivity. We demonstrate a hierarchical approach to gravitational-wave transient detection, focusing on long-lived signals, which can be used to detect transients with significance in excess of 5σ using modest computational resources. In particular, we show how previously developed seedless clustering techniques can be applied to large datasets to identify high-significance candidates without having to trade sensitivity for speed.

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.

preprint2014arXiv

Constraining the gravitational wave energy density of the Universe in the Range 0.1 Hz to 1 Hz using the Apollo Seismic Array

In this paper, we describe an analysis of Apollo era lunar seismic data that places an upper limit on an isotropic stochastic gravitational-wave background integrated over a year in the frequency range 0.1\,Hz -- 1\,Hz. We find that because the Moon's ambient noise background is much quieter than that of the Earth, significant improvements over an Earth based analysis were made. We find an upper limit of $Ω_{\rm GW}<1.2\times 10^{5}$, which is three orders of magnitude smaller than a similar analysis of a global network of broadband seismometers on Earth and the best limits in this band to date. We also discuss the benefits of a potential Earth-Moon correlation search and compute the time-dependent overlap reduction function required for such an analysis. For this search, we find an upper limit an order of magnitude larger than the Moon-Moon search.

preprint2014arXiv

Constraining the gravitational wave energy density of the Universe using Earth's ring

The search for gravitational waves is one of today's major scientific endeavors. A gravitational wave can interact with matter by exciting vibrations of elastic bodies. Earth itself is a large elastic body whose so-called normal-mode oscillations ring up when a gravitational wave passes. Therefore, precise measurement of vibration amplitudes can be used to search for the elusive gravitational-wave signals. Earth's free oscillations that can be observed after high-magnitude earthquakes have been studied extensively with gravimeters and low-frequency seismometers over many decades leading to invaluable insight into Earth's structure. Making use of our detailed understanding of Earth's normal modes, numerical models are employed for the first time to accurately calculate Earth's gravitational-wave response, and thereby turn a network of sensors that so far has served to improve our understanding of Earth, into an astrophysical observatory exploring our Universe. In this article, we constrain the energy density of gravitational waves to values in the range 0.035 - 0.15 normalized by the critical energy density of the Universe at frequencies between 0.3mHz and 5mHz, using 10 years of data from the gravimeter network of the Global Geodynamics Project that continuously monitors Earth's oscillations. This work is the first step towards a systematic investigation of the sensitivity of gravimeter networks to gravitational waves. Further advance in gravimeter technology could improve sensitivity of these networks and possibly lead to gravitational-wave detection.

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.

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.

preprint2014arXiv

Real-time earthquake warning for astronomical observatories

Early earthquake warning is a rapidly developing capability that has significant ramifications for many fields, including astronomical observatories. In this work, we describe the susceptibility of astronomical facilities to seismic events, including large telescopes as well as second-generation ground-based gravitational-wave interferometers. We describe the potential warning times for observatories from current seismic networks and propose locations for future seismometers to maximize warning times.

preprint2014arXiv

Seedless clustering in all-sky searches for gravitational-wave transients

The problem of searching for unmodeled gravitational-wave bursts can be thought of as a pattern recognition problem: how to find statistically significant clusters in spectrograms of strain power when the precise signal morphology is unknown. In a previous publication, we showed how "seedless clustering" can be used to dramatically improve the sensitivity of searches for long-lived gravitational-wave transients. In order to manage the computational costs, this initial analysis focused on externally triggered searches where the source location and emission time are both known to some degree of precision. In this paper, we show how the principle of seedless clustering can be extended to facilitate computationally-feasible, all-sky searches where the direction and emission time of the source are entirely unknown. We further demonstrate that it is possible to achieve a considerable reduction in computation time by using graphical processor units (GPUs), thereby facilitating more sensitive searches.

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.

preprint2014arXiv

Upper Limit on a Stochastic Background of Gravitational Waves from Seismic Measurements in the Range 0.05 Hz to 1 Hz

In this paper, we present an upper limit of $Ω_{\rm GW}<1.2\times 10^{8}$ on an isotropic stochastic gravitational-wave (GW) background integrated over a year in the frequency range 0.05 Hz - 1 Hz, which improves current upper limits from high-precision laboratory experiments by about 9 orders of magnitude. The limit is obtained using the response of Earth itself to GWs via a free-surface effect described more than 40 years ago by F.J. Dyson. The response was measured by a global network of broadband seismometers selected to maximize the sensitivity.

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.

preprint2013arXiv

Fortifying the characterization of binary mergers in LIGO data

The study of compact binary in-spirals and mergers with gravitational wave observatories amounts to optimizing a theoretical description of the data to best reproduce the true detector output. While most of the research effort in gravitational wave data modeling focuses on the gravitational wave- forms themselves, here we will begin to improve our model of the instrument noise by introducing parameters which allow us to determine the background instrumental power spectrum while simul- taneously characterizing the astrophysical signal. We use data from the fifth LIGO science run and simulated gravitational wave signals to demonstrate how the introduction of noise parameters results in resilience of the signal characterization to variations in an initial estimation of the noise power spectral density. We find substantial improvement in the consistency of Bayes factor calculations when we are able to marginalize over uncertainty in the instrument noise level.

preprint2013arXiv

Searching for gravitational-wave transients with a qualitative signal model: seedless clustering strategies

Gravitational-wave bursts are observable as bright clusters of pixels in spectrograms of strain power. Clustering algorithms can be used to identify candidate gravitational-wave events. Clusters are often identified by grouping together seed pixels in which the power exceeds some threshold. If the gravitational-wave signal is long-lived, however, the excess power may be spread out over many pixels, none of which are bright enough to become seeds. Without seeds, the problem of detection through clustering becomes more complicated. In this paper we investigate seedless clustering algorithms in searches for long-lived narrowband gravitational-wave bursts. Using four astrophysically motivated test waveforms, we compare a seedless clustering algorithm to two algorithms using seeds. We find that the seedless algorithm can detect gravitational-wave signals (at fixed false-alarm and false-dismissal rate) at distances between 150-200% greater than those achieved with the seed-based clustering algorithms, corresponding to significantly increased detection volumes: 420-740%. This improvement in sensitivity may extend the reach of second-generation detectors such as Advanced LIGO and Advanced Virgo deeper into astrophysically interesting distances.

preprint2012arXiv

Seismic topographic scattering in the context of GW detector site selection

In this paper, we present a calculation of seismic scattering from irregular surface topography in the Born approximation. Based on US-wide topographic data, we investigate topographic scattering at specific sites to demonstrate its impact on Newtonian-noise estimation and subtraction for future gravitational-wave detectors. We find that topographic scattering at a comparatively flat site in Oregon would not pose any problems, whereas scattering at a second site in Montana leads to significant broadening of wave amplitudes in wavenumber space that would make Newtonian-noise subtraction very challenging. Therefore, it is shown that topographic scattering should be included as criterion in the site-selection process of future low-frequency gravitational-wave detectors.

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.