Source author record

E. A. Huerta

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

ResearcherUnclaimed source record

gr-qc astro-ph.HE astro-ph.IM astro-ph.CO Artificial Intelligence Machine Learning hep-ex physics.comp-ph astro-ph.GA Databases Distributed, Parallel, and Cluster Computing eess.SP hep-ph Neural and Evolutionary Computing physics.data-an physics.flu-dyn physics.ins-det physics.med-ph

Catalog footprint

What is connected

69works

18topics

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.

preprint2023arXiv

FAIR AI Models in High Energy Physics

The findable, accessible, interoperable, and reusable (FAIR) data principles provide a framework for examining, evaluating, and improving how data is shared to facilitate scientific discovery. Generalizing these principles to research software and other digital products is an active area of research. Machine learning (ML) models -- algorithms that have been trained on data without being explicitly programmed -- and more generally, artificial intelligence (AI) models, are an important target for this because of the ever-increasing pace with which AI is transforming scientific domains, such as experimental high energy physics (HEP). In this paper, we propose a practical definition of FAIR principles for AI models in HEP and describe a template for the application of these principles. We demonstrate the template's use with an example AI model applied to HEP, in which a graph neural network is used to identify Higgs bosons decaying to two bottom quarks. We report on the robustness of this FAIR AI model, its portability across hardware architectures and software frameworks, and its interpretability.

preprint2022arXiv

A FAIR and AI-ready Higgs boson decay dataset

To enable the reusability of massive scientific datasets by humans and machines, researchers aim to adhere to the principles of findability, accessibility, interoperability, and reusability (FAIR) for data and artificial intelligence (AI) models. This article provides a domain-agnostic, step-by-step assessment guide to evaluate whether or not a given dataset meets these principles. We demonstrate how to use this guide to evaluate the FAIRness of an open simulated dataset produced by the CMS Collaboration at the CERN Large Hadron Collider. This dataset consists of Higgs boson decays and quark and gluon background, and is available through the CERN Open Data Portal. We use additional available tools to assess the FAIRness of this dataset, and incorporate feedback from members of the FAIR community to validate our results. This article is accompanied by a Jupyter notebook to visualize and explore this dataset. This study marks the first in a planned series of articles that will guide scientists in the creation of FAIR AI models and datasets in high energy particle physics.

preprint2022arXiv

Constraints from LIGO O3 data on gravitational-wave emission due to r-modes in the glitching pulsar PSR J0537-6910

We present a search for continuous gravitational-wave emission due to r-modes in the pulsar PSR J0537-6910 using data from the LIGO-Virgo Collaboration observing run O3. PSR J0537-6910 is a young energetic X-ray pulsar and is the most frequent glitcher known. The inter-glitch braking index of the pulsar suggests that gravitational-wave emission due to r-mode oscillations may play an important role in the spin evolution of this pulsar. Theoretical models confirm this possibility and predict emission at a level that can be probed by ground-based detectors. In order to explore this scenario, we search for r-mode emission in the epochs between glitches by using a contemporaneous timing ephemeris obtained from NICER data. We do not detect any signals in the theoretically expected band of 86-97 Hz, and report upper limits on the amplitude of the gravitational waves. Our results improve on previous amplitude upper limits from r-modes in J0537-6910 by a factor of up to 3 and place stringent constraints on theoretical models for r-mode driven spin-down in PSR J0537-6910, especially for higher frequencies at which our results reach below the spin-down limit defined by energy conservation.

preprint2022arXiv

Inference-optimized AI and high performance computing for gravitational wave detection at scale

We introduce an ensemble of artificial intelligence models for gravitational wave detection that we trained in the Summit supercomputer using 32 nodes, equivalent to 192 NVIDIA V100 GPUs, within 2 hours. Once fully trained, we optimized these models for accelerated inference using NVIDIA TensorRT. We deployed our inference-optimized AI ensemble in the ThetaGPU supercomputer at Argonne Leadership Computer Facility to conduct distributed inference. Using the entire ThetaGPU supercomputer, consisting of 20 nodes each of which has 8 NVIDIA A100 Tensor Core GPUs and 2 AMD Rome CPUs, our NVIDIA TensorRT-optimized AI ensemble processed an entire month of advanced LIGO data (including Hanford and Livingston data streams) within 50 seconds. Our inference-optimized AI ensemble retains the same sensitivity of traditional AI models, namely, it identifies all known binary black hole mergers previously identified in this advanced LIGO dataset and reports no misclassifications, while also providing a 3X inference speedup compared to traditional artificial intelligence models. We used time slides to quantify the performance of our AI ensemble to process up to 5 years worth of advanced LIGO data. In this synthetically enhanced dataset, our AI ensemble reports an average of one misclassification for every month of searched advanced LIGO data. We also present the receiver operating characteristic curve of our AI ensemble using this 5 year long advanced LIGO dataset. This approach provides the required tools to conduct accelerated, AI-driven gravitational wave detection at scale.

preprint2022arXiv

Interpretable AI forecasting for numerical relativity waveforms of quasi-circular, spinning, non-precessing binary black hole mergers

We present a deep-learning artificial intelligence model that is capable of learning and forecasting the late-inspiral, merger and ringdown of numerical relativity waveforms that describe quasi-circular, spinning, non-precessing binary black hole mergers. We used the NRHybSur3dq8 surrogate model to produce train, validation and test sets of $\ell=|m|=2$ waveforms that cover the parameter space of binary black hole mergers with mass-ratios $q\leq8$ and individual spins $|s^z_{\{1,2\}}| \leq 0.8$. These waveforms cover the time range $t\in[-5000\textrm{M}, 130\textrm{M}]$, where $t=0M$ marks the merger event, defined as the maximum value of the waveform amplitude. We harnessed the ThetaGPU supercomputer at the Argonne Leadership Computing Facility to train our AI model using a training set of 1.5 million waveforms. We used 16 NVIDIA DGX A100 nodes, each consisting of 8 NVIDIA A100 Tensor Core GPUs and 2 AMD Rome CPUs, to fully train our model within 3.5 hours. Our findings show that artificial intelligence can accurately forecast the dynamical evolution of numerical relativity waveforms in the time range $t\in[-100\textrm{M}, 130\textrm{M}]$. Sampling a test set of 190,000 waveforms, we find that the average overlap between target and predicted waveforms is $\gtrsim99\%$ over the entire parameter space under consideration. We also combined scientific visualization and accelerated computing to identify what components of our model take in knowledge from the early and late-time waveform evolution to accurately forecast the latter part of numerical relativity waveforms. This work aims to accelerate the creation of scalable, computationally efficient and interpretable artificial intelligence models for gravitational wave astrophysics.

preprint2022arXiv

Interpreting a Machine Learning Model for Detecting Gravitational Waves

We describe a case study of translational research, applying interpretability techniques developed for computer vision to machine learning models used to search for and find gravitational waves. The models we study are trained to detect black hole merger events in non-Gaussian and non-stationary advanced Laser Interferometer Gravitational-wave Observatory (LIGO) data. We produced visualizations of the response of machine learning models when they process advanced LIGO data that contains real gravitational wave signals, noise anomalies, and pure advanced LIGO noise. Our findings shed light on the responses of individual neurons in these machine learning models. Further analysis suggests that different parts of the network appear to specialize in local versus global features, and that this difference appears to be rooted in the branched architecture of the network as well as noise characteristics of the LIGO detectors. We believe efforts to whiten these "black box" models can suggest future avenues for research and help inform the design of interpretable machine learning models for gravitational wave astrophysics.

preprint2022arXiv

Search for anisotropic gravitational-wave backgrounds using data from Advanced LIGO and Advanced Virgo's first three observing runs

We report results from searches for anisotropic stochastic gravitational-wave backgrounds using data from the first three observing runs of the Advanced LIGO and Advanced Virgo detectors. For the first time, we include Virgo data in our analysis and run our search with a new efficient pipeline called {\tt PyStoch} on data folded over one sidereal day. We use gravitational-wave radiometry (broadband and narrow band) to produce sky maps of stochastic gravitational-wave backgrounds and to search for gravitational waves from point sources. A spherical harmonic decomposition method is employed to look for gravitational-wave emission from spatially-extended sources. Neither technique found evidence of gravitational-wave signals. Hence we derive 95\% confidence-level upper limit sky maps on the gravitational-wave energy flux from broadband point sources, ranging from $F_{α, Θ} < {\rm (0.013 - 7.6)} \times 10^{-8} {\rm erg \, cm^{-2} \, s^{-1} \, Hz^{-1}},$ and on the (normalized) gravitational-wave energy density spectrum from extended sources, ranging from $Ω_{α, Θ} < {\rm (0.57 - 9.3)} \times 10^{-9} \, {\rm sr^{-1}}$, depending on direction ($Θ$) and spectral index ($α$). These limits improve upon previous limits by factors of $2.9 - 3.5$. We also set 95\% confidence level upper limits on the frequency-dependent strain amplitudes of quasimonochromatic gravitational waves coming from three interesting targets, Scorpius X-1, SN 1987A and the Galactic Center, with best upper limits range from $h_0 < {\rm (1.7-2.1)} \times 10^{-25},$ a factor of $\geq 2.0$ improvement compared to previous stochastic radiometer searches.

preprint2021arXiv

Advances in Machine and Deep Learning for Modeling and Real-time Detection of Multi-Messenger Sources

We live in momentous times. The science community is empowered with an arsenal of cosmic messengers to study the Universe in unprecedented detail. Gravitational waves, electromagnetic waves, neutrinos and cosmic rays cover a wide range of wavelengths and time scales. Combining and processing these datasets that vary in volume, speed and dimensionality requires new modes of instrument coordination, funding and international collaboration with a specialized human and technological infrastructure. In tandem with the advent of large-scale scientific facilities, the last decade has experienced an unprecedented transformation in computing and signal processing algorithms. The combination of graphics processing units, deep learning, and the availability of open source, high-quality datasets, have powered the rise of artificial intelligence. This digital revolution now powers a multi-billion dollar industry, with far-reaching implications in technology and society. In this chapter we describe pioneering efforts to adapt artificial intelligence algorithms to address computational grand challenges in Multi-Messenger Astrophysics. We review the rapid evolution of these disruptive algorithms, from the first class of algorithms introduced in early 2017, to the sophisticated algorithms that now incorporate domain expertise in their architectural design and optimization schemes. We discuss the importance of scientific visualization and extreme-scale computing in reducing time-to-insight and obtaining new knowledge from the interplay between models and data.

preprint2021arXiv

Open data from the first and second observing runs of Advanced LIGO and Advanced Virgo

Advanced LIGO and Advanced Virgo are actively monitoring the sky and collecting gravitational-wave strain data with sufficient sensitivity to detect signals routinely. In this paper we describe the data recorded by these instruments during their first and second observing runs. The main data products are the gravitational-wave strain arrays, released as time series sampled at 16384 Hz. The datasets that include this strain measurement can be freely accessed through the Gravitational Wave Open Science Center at http://gw-openscience.org, together with data-quality information essential for the analysis of LIGO and Virgo data, documentation, tutorials, and supporting software.

preprint2021arXiv

Search for lensing signatures in the gravitational-wave observations from the first half of LIGO-Virgo's third observing run

We search for signatures of gravitational lensing in the gravitational-wave signals from compact binary coalescences detected by Advanced LIGO and Advanced Virgo during O3a, the first half of their third observing run. We study: 1) the expected rate of lensing at current detector sensitivity and the implications of a non-observation of strong lensing or a stochastic gravitational-wave background on the merger-rate density at high redshift; 2) how the interpretation of individual high-mass events would change if they were found to be lensed; 3) the possibility of multiple images due to strong lensing by galaxies or galaxy clusters; and 4) possible wave-optics effects due to point-mass microlenses. Several pairs of signals in the multiple-image analysis show similar parameters and, in this sense, are nominally consistent with the strong lensing hypothesis. However, taking into account population priors, selection effects, and the prior odds against lensing, these events do not provide sufficient evidence for lensing. Overall, we find no compelling evidence for lensing in the observed gravitational-wave signals from any of these analyses.

preprint2020arXiv

A guide to LIGO-Virgo detector noise and extraction of transient gravitational-wave signals

The LIGO Scientific Collaboration and the Virgo Collaboration have cataloged eleven confidently detected gravitational-wave events during the first two observing runs of the advanced detector era. All eleven events were consistent with being from well-modeled mergers between compact stellar-mass objects: black holes or neutron stars. The data around the time of each of these events have been made publicly available through the gravitational-wave open science center. The entirety of the gravitational-wave strain data from the first and second observing runs have also now been made publicly available. There is considerable interest among the broad scientific community in understanding the data and methods used in the analyses. In this paper, we provide an overview of the detector noise properties and the data analysis techniques used to detect gravitational-wave signals and infer the source properties. We describe some of the checks that are performed to validate the analyses and results from the observations of gravitational-wave events. We also address concerns that have been raised about various properties of LIGO-Virgo detector noise and the correctness of our analyses as applied to the resulting data.

preprint2020arXiv

A Joint Fermi-GBM and LIGO/Virgo Analysis of Compact Binary Mergers From the First and Second Gravitational-wave Observing Runs

We present results from offline searches of Fermi Gamma-ray Burst Monitor (GBM) data for gamma-ray transients coincident with the compact binary coalescences observed by the gravitational-wave (GW) detectors Advanced LIGO and Advanced Virgo during their first and second observing runs. In particular, we perform follow-up for both confirmed events and low significance candidates reported in the LIGO/Virgo catalog GWTC-1. We search for temporal coincidences between these GW signals and GBM triggered gamma-ray bursts (GRBs). We also use the GBM Untargeted and Targeted subthreshold searches to find coincident gamma-rays below the on-board triggering threshold. This work implements a refined statistical approach by incorporating GW astrophysical source probabilities and GBM visibilities of LIGO/Virgo sky localizations to search for cumulative signatures of coincident subthreshold gamma-rays. All search methods recover the short gamma-ray burst GRB 170817A occurring ~1.7 s after the binary neutron star merger GW170817. We also present results from a new search seeking GBM counterparts to LIGO single-interferometer triggers. This search finds a candidate joint event, but given the nature of the GBM signal and localization, as well as the high joint false alarm rate of $1.1 \times 10^{-6}$ Hz, we do not consider it an astrophysical association. We find no additional coincidences.

preprint2020arXiv

Artificial neural network subgrid models of 2-D compressible magnetohydrodynamic turbulence

We explore the suitability of deep learning to capture the physics of subgrid-scale ideal magnetohydrodynamics turbulence of 2-D simulations of the magnetized Kelvin-Helmholtz instability. We produce simulations at different resolutions to systematically quantify the performance of neural network models to reproduce the physics of these complex simulations. We compare the performance of our neural networks with gradient models, which are extensively used in the extensively in the magnetohydrodynamic literature. Our findings indicate that neural networks significantly outperform gradient models at reproducing the effects of magnetohydrodynamics turbulence. To the best of our knowledge, this is the first exploratory study on the use of deep learning to learn and reproduce the physics of magnetohydrodynamics turbulence.

preprint2020arXiv

Binary Black Hole Population Properties Inferred from the First and Second Observing Runs of Advanced LIGO and Advanced Virgo

We present results on the mass, spin, and redshift distributions with phenomenological population models using the ten binary black hole mergers detected in the first and second observing runs completed by Advanced LIGO and Advanced Virgo. We constrain properties of the binary black hole (BBH) mass spectrum using models with a range of parameterizations of the BBH mass and spin distributions. We find that the mass distribution of the more massive black hole in such binaries is well approximated by models with no more than 1% of black holes more massive than $45\,M_\odot$, and a power law index of $α= {1.3}^{+1.4}_{-1.7}$ (90% credibility). We also show that BBHs are unlikely to be composed of black holes with large spins aligned to the orbital angular momentum. Modelling the evolution of the BBH merger rate with redshift, we show that it is at or increasing with redshift with 93% probability. Marginalizing over uncertainties in the BBH population, we find robust estimates of the BBH merger rate density of $R = {53.2}^{+55.8}_{-28.2}$ Gpc$^{-3}$ yr$^{-1}$ (90% credibility). As the BBH catalog grows in future observing runs, we expect that uncertainties in the population model parameters will shrink, potentially providing insights into the formation of black holes via supernovae, binary interactions of massive stars, stellar cluster dynamics, and the formation history of black holes across cosmic time.

preprint2020arXiv

Deep Learning for Cardiologist-level Myocardial Infarction Detection in Electrocardiograms

Myocardial infarction is the leading cause of death worldwide. In this paper, we design domain-inspired neural network models to detect myocardial infarction. First, we study the contribution of various leads. This systematic analysis, first of its kind in the literature, indicates that out of 15 ECG leads, data from the v6, vz, and ii leads are critical to correctly identify myocardial infarction. Second, we use this finding and adapt the ConvNetQuake neural network model--originally designed to identify earthquakes--to attain state-of-the-art classification results for myocardial infarction, achieving $99.43\%$ classification accuracy on a record-wise split, and $97.83\%$ classification accuracy on a patient-wise split. These two results represent cardiologist-level performance level for myocardial infarction detection after feeding only 10 seconds of raw ECG data into our model. Third, we show that our multi-ECG-channel neural network achieves cardiologist-level performance without the need of any kind of manual feature extraction or data pre-processing.

preprint2020arXiv

Deep learning for gravitational wave forecasting of neutron star mergers

We introduce deep learning time-series forecasting for gravitational wave detection of binary neutron star mergers. This method enables the identification of these signals in real advanced LIGO data up to 30 seconds before merger. When applied to GW170817, our deep learning forecasting method identifies the presence of this gravitational wave signal 10 seconds before merger. This novel approach requires a single GPU for inference, and may be used as part of an early warning system for time-sensitive multi-messenger searches.

preprint2020arXiv

Deep Transfer Learning for Star Cluster Classification: I. Application to the PHANGS-HST Survey

We present the results of a proof-of-concept experiment which demonstrates that deep learning can successfully be used for production-scale classification of compact star clusters detected in HST UV-optical imaging of nearby spiral galaxies (D < 20 Mpc) in the PHANGS-HST survey. Given the relatively small nature of existing, human-labelled star cluster samples, we transfer the knowledge of state-of-the-art neural network models for real-object recognition to classify star clusters candidates into four morphological classes. We perform a series of experiments to determine the dependence of classification performance on: neural network architecture (ResNet18 and VGG19-BN); training data sets curated by either a single expert or three astronomers; and the size of the images used for training. We find that the overall classification accuracies are not significantly affected by these choices. The networks are used to classify star cluster candidates in the PHANGS-HST galaxy NGC 1559, which was not included in the training samples. The resulting prediction accuracies are 70%, 40%, 40-50%, 50-70% for class 1, 2, 3 star clusters, and class 4 non-clusters respectively. This performance is competitive with consistency achieved in previously published human and automated quantitative classification of star cluster candidate samples (70-80%, 40-50%, 40-50%, and 60-70%). The methods introduced herein lay the foundations to automate classification for star clusters at scale, and exhibit the need to prepare a standardized dataset of human-labelled star cluster classifications, agreed upon by a full range of experts in the field, to further improve the performance of the networks introduced in this study.

preprint2020arXiv

GW190412: Observation of a Binary-Black-Hole Coalescence with Asymmetric Masses

We report the observation of gravitational waves from a binary-black-hole coalescence during the first two weeks of LIGO's and Virgo's third observing run. The signal was recorded on April 12, 2019 at 05:30:44 UTC with a network signal-to-noise ratio of 19. The binary is different from observations during the first two observing runs most notably due to its asymmetric masses: a ~30 solar mass black hole merged with a ~8 solar mass black hole companion. The more massive black hole rotated with a dimensionless spin magnitude between 0.22 and 0.60 (90% probability). Asymmetric systems are predicted to emit gravitational waves with stronger contributions from higher multipoles, and indeed we find strong evidence for gravitational radiation beyond the leading quadrupolar order in the observed signal. A suite of tests performed on GW190412 indicates consistency with Einstein's general theory of relativity. While the mass ratio of this system differs from all previous detections, we show that it is consistent with the population model of stellar binary black holes inferred from the first two observing runs.

preprint2020arXiv

GW190425: Observation of a Compact Binary Coalescence with Total Mass $\sim 3.4 M_{\odot}$

On 2019 April 25, the LIGO Livingston detector observed a compact binary coalescence with signal-to-noise ratio 12.9. The Virgo detector was also taking data that did not contribute to detection due to a low signal-to-noise ratio, but were used for subsequent parameter estimation. The 90% credible intervals for the component masses range from 1.12 to 2.52 $M_{\odot}$ (1.45 to 1.88 $M_{\odot}$ if we restrict the dimensionless component spin magnitudes to be smaller than 0.05). These mass parameters are consistent with the individual binary components being neutron stars. However, both the source-frame chirp mass $1.44^{+0.02}_{-0.02} M_{\odot}$ and the total mass $3.4^{+0.3}_{-0.1}\,M_{\odot}$ of this system are significantly larger than those of any other known binary neutron star system. The possibility that one or both binary components of the system are black holes cannot be ruled out from gravitational-wave data. We discuss possible origins of the system based on its inconsistency with the known Galactic binary neutron star population. Under the assumption that the signal was produced by a binary neutron star coalescence, the local rate of neutron star mergers is updated to $250-2810 \text{Gpc}^{-3}\text{yr}^{-1}$.

preprint2020arXiv

GW190521: A Binary Black Hole Merger with a Total Mass of $150 ~ M_{\odot}$

On May 21, 2019 at 03:02:29 UTC Advanced LIGO and Advanced Virgo observed a short duration gravitational-wave signal, GW190521, with a three-detector network signal-to-noise ratio of 14.7, and an estimated false-alarm rate of 1 in 4900 yr using a search sensitive to generic transients. If GW190521 is from a quasicircular binary inspiral, then the detected signal is consistent with the merger of two black holes with masses of $85^{+21}_{-14} M_{\odot}$ and $66^{+17}_{-18} M_{\odot}$ (90 % credible intervals). We infer that the primary black hole mass lies within the gap produced by (pulsational) pair-instability supernova processes, and has only a 0.32 % probability of being below $65 M_{\odot}$. We calculate the mass of the remnant to be $142^{+28}_{-16} M_{\odot}$, which can be considered an intermediate mass black hole (IMBH). The luminosity distance of the source is $5.3^{+2.4}_{-2.6}$ Gpc, corresponding to a redshift of $0.82^{+0.28}_{-0.34}$. The inferred rate of mergers similar to GW190521 is $0.13^{+0.30}_{-0.11}\,\mathrm{Gpc}^{-3}\,\mathrm{yr}^{-1}$.

preprint2020arXiv

GW190814: Gravitational Waves from the Coalescence of a 23 M$_\odot$ Black Hole with a 2.6 M$_\odot$ Compact Object

We report the observation of a compact binary coalescence involving a 22.2 - 24.3 $M_{\odot}$ black hole and a compact object with a mass of 2.50 - 2.67 $M_{\odot}$ (all measurements quoted at the 90$\%$ credible level). The gravitational-wave signal, GW190814, was observed during LIGO's and Virgo's third observing run on August 14, 2019 at 21:10:39 UTC and has a signal-to-noise ratio of 25 in the three-detector network. The source was localized to 18.5 deg$^2$ at a distance of $241^{+41}_{-45}$ Mpc; no electromagnetic counterpart has been confirmed to date. The source has the most unequal mass ratio yet measured with gravitational waves, $0.112^{+0.008}_{-0.009}$, and its secondary component is either the lightest black hole or the heaviest neutron star ever discovered in a double compact-object system. The dimensionless spin of the primary black hole is tightly constrained to $\leq 0.07$. Tests of general relativity reveal no measurable deviations from the theory, and its prediction of higher-multipole emission is confirmed at high confidence. We estimate a merger rate density of 1-23 Gpc$^{-3}$ yr$^{-1}$ for the new class of binary coalescence sources that GW190814 represents. Astrophysical models predict that binaries with mass ratios similar to this event can form through several channels, but are unlikely to have formed in globular clusters. However, the combination of mass ratio, component masses, and the inferred merger rate for this event challenges all current models for the formation and mass distribution of compact-object binaries.

preprint2020arXiv

Model comparison from LIGO-Virgo data on GW170817's binary components and consequences for the merger remnant

GW170817 is the very first observation of gravitational waves originating from the coalescence of two compact objects in the mass range of neutron stars, accompanied by electromagnetic counterparts, and offers an opportunity to directly probe the internal structure of neutron stars. We perform Bayesian model selection on a wide range of theoretical predictions for the neutron star equation of state. For the binary neutron star hypothesis, we find that we cannot rule out the majority of theoretical models considered. In addition, the gravitational-wave data alone does not rule out the possibility that one or both objects were low-mass black holes. We discuss the possible outcomes in the case of a binary neutron star merger, finding that all scenarios from prompt collapse to long-lived or even stable remnants are possible. For long-lived remnants, we place an upper limit of 1.9 kHz on the rotation rate. If a black hole was formed any time after merger and the coalescing stars were slowly rotating, then the maximum baryonic mass of non-rotating neutron stars is at most 3.05 $M_\odot$, and three equations of state considered here can be ruled out. We obtain a tighter limit of 2.67 $M_\odot$ for the case that the merger results in a hypermassive neutron star.

preprint2020arXiv

Modeling the uncertainties of solar-system ephemerides for robust gravitational-wave searches with pulsar timing arrays

The regularity of pulsar emissions becomes apparent once we reference the pulses' times of arrivals to the inertial rest frame of the solar system. It follows that errors in the determination of Earth's position with respect to the solar-system barycenter can appear as a time-correlated bias in pulsar-timing residual time series, affecting the searches for low-frequency gravitational waves performed with pulsar timing arrays. Indeed, recent array datasets yield different gravitational-wave background upper limits and detection statistics when analyzed with different solar-system ephemerides. Crucially, the ephemerides do not generally provide usable error representations. In this article we describe the motivation, construction, and application of a physical model of solar-system ephemeris uncertainties, which focuses on the degrees of freedom (Jupiter's orbital elements) most relevant to gravitational-wave searches with pulsar timing arrays. This model, BayesEphem, was used to derive ephemeris-robust results in NANOGrav's 11-yr stochastic-background search, and it provides a foundation for future searches by NANOGrav and other consortia. The analysis and simulations reported here suggest that ephemeris modeling reduces the gravitational-wave sensitivity of the 11-yr dataset; and that this degeneracy will vanish with improved ephemerides and with the longer pulsar timing datasets that will become available in the near future.

preprint2020arXiv

Physics-inspired deep learning to characterize the signal manifold of quasi-circular, spinning, non-precessing binary black hole mergers

The spin distribution of binary black hole mergers contains key information concerning the formation channels of these objects, and the astrophysical environments where they form, evolve and coalesce. To quantify the suitability of deep learning to characterize the signal manifold of quasi-circular, spinning, non-precessing binary black hole mergers, we introduce a modified version of WaveNet trained with a novel optimization scheme that incorporates general relativistic constraints of the spin properties of astrophysical black holes. The neural network model is trained, validated and tested with 1.5 million $\ell=|m|=2$ waveforms generated within the regime of validity of NRHybSur3dq8, i.e., mass-ratios $q\leq8$ and individual black hole spins $ | s^z_{\{1,\,2\}} | \leq 0.8$. Using this neural network model, we quantify how accurately we can infer the astrophysical parameters of black hole mergers in the absence of noise. We do this by computing the overlap between waveforms in the testing data set and the corresponding signals whose mass-ratio and individual spins are predicted by our neural network. We find that the convergence of high performance computing and physics-inspired optimization algorithms enable an accurate reconstruction of the mass-ratio and individual spins of binary black hole mergers across the parameter space under consideration. This is a significant step towards an informed utilization of physics-inspired deep learning models to reconstruct the spin distribution of binary black hole mergers in realistic detection scenarios.

preprint2020arXiv

Properties and astrophysical implications of the 150 Msun binary black hole merger GW190521

The gravitational-wave signal GW190521 is consistent with a binary black hole merger source at redshift 0.8 with unusually high component masses, $85^{+21}_{-14}\,M_{\odot}$ and $66^{+17}_{-18}\,M_{\odot}$, compared to previously reported events, and shows mild evidence for spin-induced orbital precession. The primary falls in the mass gap predicted by (pulsational) pair-instability supernova theory, in the approximate range $65 - 120\,M_{\odot}$. The probability that at least one of the black holes in GW190521 is in that range is 99.0%. The final mass of the merger $(142^{+28}_{-16}\,M_{\odot})$ classifies it as an intermediate-mass black hole. Under the assumption of a quasi-circular binary black hole coalescence, we detail the physical properties of GW190521's source binary and its post-merger remnant, including component masses and spin vectors. Three different waveform models, as well as direct comparison to numerical solutions of general relativity, yield consistent estimates of these properties. Tests of strong-field general relativity targeting the merger-ringdown stages of coalescence indicate consistency of the observed signal with theoretical predictions. We estimate the merger rate of similar systems to be $0.13^{+0.30}_{-0.11}\,{\rm Gpc}^{-3}\,\rm{yr}^{-1}$. We discuss the astrophysical implications of GW190521 for stellar collapse, and for the possible formation of black holes in the pair-instability mass gap through various channels: via (multiple) stellar coalescence, or via hierarchical merger of lower-mass black holes in star clusters or in active galactic nuclei. We find it to be unlikely that GW190521 is a strongly lensed signal of a lower-mass black hole binary merger. We also discuss more exotic possible sources for GW190521, including a highly eccentric black hole binary, or a primordial black hole binary.

preprint2020arXiv

Search for Eccentric Binary Black Hole Mergers with Advanced LIGO and Advanced Virgo during their First and Second Observing Runs

When formed through dynamical interactions, stellar-mass binary black holes may retain eccentric orbits ($e>0.1$ at 10 Hz) detectable by ground-based gravitational-wave detectors. Eccentricity can therefore be used to differentiate dynamically-formed binaries from isolated binary black hole mergers. Current template-based gravitational-wave searches do not use waveform models associated to eccentric orbits, rendering the search less efficient to eccentric binary systems. Here we present results of a search for binary black hole mergers that inspiral in eccentric orbits using data from the first and second observing runs (O1 and O2) of Advanced LIGO and Advanced Virgo. The search uses minimal assumptions on the morphology of the transient gravitational waveform. We show that it is sensitive to binary mergers with a detection range that is weakly dependent on eccentricity for all bound systems. Our search did not identify any new binary merger candidates. We interpret these results in light of eccentric binary formation models.

preprint2020arXiv

Search for intermediate mass black hole binaries in the first and second observing runs of the Advanced LIGO and Virgo network

Gravitational wave astronomy has been firmly established with the detection of gravitational waves from the merger of ten stellar mass binary black holes and a neutron star binary. This paper reports on the all-sky search for gravitational waves from intermediate mass black hole binaries in the first and second observing runs of the Advanced LIGO and Virgo network. The search uses three independent algorithms: two based on matched filtering of the data with waveform templates of gravitational wave signals from compact binaries, and a third, model-independent algorithm that employs no signal model for the incoming signal. No intermediate mass black hole binary event was detected in this search. Consequently, we place upper limits on the merger rate density for a family of intermediate mass black hole binaries. In particular, we choose sources with total masses $M=m_1+m_2\in[120,800]$M$_\odot$ and mass ratios $q = m_2/m_1 \in[0.1,1.0]$. For the first time, this calculation is done using numerical relativity waveforms (which include higher modes) as models of the real emitted signal. We place a most stringent upper limit of $0.20$~Gpc$^{-3}$yr$^{-1}$ (in co-moving units at the 90% confidence level) for equal-mass binaries with individual masses $m_{1,2}=100$M$_\odot$ and dimensionless spins $χ_{1,2}= 0.8$ aligned with the orbital angular momentum of the binary. This improves by a factor of $\sim 5$ that reported after Advanced LIGO's first observing run.

preprint2020arXiv

Searches for Gravitational Waves from Known Pulsars at Two Harmonics in 2015-2017 LIGO Data

We present a search for gravitational waves from 221 pulsars with rotation frequencies $\gtrsim 10$ Hz. We use advanced LIGO data from its first and second observing runs spanning 2015-2017, which provides the highest-sensitivity gravitational-wave data so far obtained. In this search we target emission from both the $l = m = 2$ mass quadrupole mode, with a frequency at twice that of the pulsar's rotation, and from the $l = 2$, $m = 1$ mode, with a frequency at the pulsar rotation frequency. The search finds no evidence for gravitational-wave emission from any pulsar at either frequency. For the $l = m = 2$ mode search, we provide updated upper limits on the gravitational-wave amplitude, mass quadrupole moment, and fiducial ellipticity for 167 pulsars, and the first such limits for a further 55. For 20 young pulsars these results give limits that are below those inferred from the pulsars' spin-down. For the Crab and Vela pulsars our results constrain gravitational-wave emission to account for less than 0.017% and 0.18% of the spin-down luminosity, respectively. For the recycled millisecond pulsar J0711-6830 our limits are only a factor of 1.3 above the spin-down limit, assuming the canonical value of $10^{38}$ kg m$^2$ for the star's moment of inertia, and imply a gravitational-wave-derived upper limit on the star's ellipticity of $1.2\!\times\!10^{-8}$. We also place new limits on the emission amplitude at the rotation frequency of the pulsars.

preprint2019arXiv

An Optically Targeted Search for Gravitational Waves emitted by Core-Collapse Supernovae during the First and Second Observing Runs of Advanced LIGO and Advanced Virgo

We present the results from a search for gravitational-wave transients associated with core-collapse supernovae observed within a source distance of approximately 20 Mpc during the first and second observing runs of Advanced LIGO and Advanced Virgo. No significant gravitational-wave candidate was detected. We report the detection efficiencies as a function of the distance for waveforms derived from multidimensional numerical simulations and phenomenological extreme emission models. For neutrino-driven explosions the distance at which we reach 50% detection efficiency is approaching 5 kpc, and for magnetorotationally-driven explosions is up to 54 kpc. However, waveforms for extreme emission models are detectable up to 28 Mpc. For the first time, the gravitational-wave data enabled us to exclude part of the parameter spaces of two extreme emission models with confidence up to 83%, limited by coincident data coverage. Besides, using ad hoc harmonic signals windowed with Gaussian envelopes we constrained the gravitational-wave energy emitted during core-collapse at the levels of $4.27\times 10^{-4}\,M_\odot c^2$ and $1.28\times 10^{-1}\,M_\odot c^2$ for emissions at 235 Hz and 1304 Hz respectively. These constraints are two orders of magnitude more stringent than previously derived in the corresponding analysis using initial LIGO, initial Virgo and GEO 600 data.

preprint2019arXiv

Search for sub-solar mass ultracompact binaries in Advanced LIGO's second observing run

We present an Advanced LIGO and Advanced Virgo search for sub-solar mass ultracompact objects in data obtained during Advanced LIGO's second observing run. In contrast to a previous search of Advanced LIGO data from the first observing run, this search includes the effects of component spin on the gravitational waveform. We identify no viable gravitational wave candidates consistent with sub-solar mass ultracompact binaries with at least one component between 0.2 - 1.0 solar masses. We use the null result to constrain the binary merger rate of (0.2 solar mass, 0.2 solar mass) binaries to be less than 3.7 x 10^5 Gpc^-3 yr^-1 and the binary merger rate of (1.0 solar mass, 1.0 solar mass) binaries to be less than 5.2 x 10^3 Gpc^-3 yr^-1. Sub-solar mass ultracompact objects are not expected to form via known stellar evolution channels, though it has been suggested that primordial density fluctuations or particle dark matter with cooling mechanisms and/or nuclear interactions could form black holes with sub-solar masses. Assuming a particular primordial black hole formation model, we constrain a population of merging 0.2 solar mass black holes to account for less than 16% of the dark matter density and a population of merging 1.0 solar mass black holes to account for less than 2% of the dark matter density. We discuss how constraints on the merger rate and dark matter fraction may be extended to arbitrary black hole population models that predict sub-solar mass binaries.

preprint2019arXiv

The NANOGrav 11-Year Data Set: Evolution of Gravitational Wave Background Statistics

An ensemble of inspiraling supermassive black hole binaries should produce a stochastic background of very low frequency gravitational waves. This stochastic background is predicted to be a power law, with a spectral index of -2/3, and it should be detectable by a network of precisely timed millisecond pulsars, widely distributed on the sky. This paper reports a new "time slicing" analysis of the 11-year data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) using 34 millisecond pulsars. Methods to flag potential "false positive" signatures are developed, including techniques to identify responsible pulsars. Mitigation strategies are then presented. We demonstrate how an incorrect noise model can lead to spurious signals, and show how independently modeling noise across 30 Fourier components, spanning NANOGrav's frequency range, effectively diagnoses and absorbs the excess power in gravitational-wave searches. This results in a nominal, and expected, progression of our gravitational-wave statistics. Additionally we show that the first interstellar medium event in PSR J1713+0747 pollutes the common red noise process with low-spectral index noise, and use a tailored noise model to remove these effects.

preprint2019arXiv

The NANOGrav 11-Year Data Set: Limits on Gravitational Wave Memory

The mergers of supermassive black hole binaries (SMBHBs) promise to be incredible sources of gravitational waves (GWs). While the oscillatory part of the merger gravitational waveform will be outside the frequency sensitivity range of pulsar timing arrays (PTAs), the non-oscillatory GW memory effect is detectable. Further, any burst of gravitational waves will produce GW memory, making memory a useful probe of unmodeled exotic sources and new physics. We searched the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GW memory. This dataset is sensitive to very low frequency GWs of $\sim3$ to $400$ nHz (periods of $\sim11$ yr $-$ $1$ mon). Finding no evidence for GWs, we placed limits on the strain amplitude of GW memory events during the observation period. We then used the strain upper limits to place limits on the rate of GW memory causing events. At a strain of $2.5\times10^{-14}$, corresponding to the median upper limit as a function of source sky position, we set a limit on the rate of GW memory events at $<0.4$ yr$^{-1}$. That strain corresponds to a SMBHB merger with reduced mass of $ηM \sim 2\times10^{10}M_\odot$ and inclination of $ι=π/3$ at a distance of 1 Gpc. As a test of our analysis, we analyzed the NANOGrav 9-year data set as well. This analysis found an anomolous signal, which does not appear in the 11-year data set. This signal is not a GW, and its origin remains unknown.

preprint2018arXiv

Search for Multi-messenger Sources of Gravitational Waves and High-energy Neutrinos with Advanced LIGO during its first Observing Run, ANTARES and IceCube

Astrophysical sources of gravitational waves, such as binary neutron star and black hole mergers or core-collapse supernovae, can drive relativistic outflows, giving rise to non-thermal high-energy emission. High-energy neutrinos are signatures of such outflows. The detection of gravitational waves and high-energy neutrinos from common sources could help establish the connection between the dynamics of the progenitor and the properties of the outflow. We searched for associated emission of gravitational waves and high-energy neutrinos from astrophysical transients with minimal assumptions using data from Advanced LIGO from its first observing run O1, and data from the ANTARES and IceCube neutrino observatories from the same time period. We focused on candidate events whose astrophysical origin could not be determined from a single messenger. We found no significant coincident candidate, which we used to constrain the rate density of astrophysical sources dependent on their gravitational wave and neutrino emission processes.

preprint2017arXiv

Denoising Gravitational Waves using Deep Learning with Recurrent Denoising Autoencoders

Gravitational wave astronomy is a rapidly growing field of modern astrophysics, with observations being made frequently by the LIGO detectors. Gravitational wave signals are often extremely weak and the data from the detectors, such as LIGO, is contaminated with non-Gaussian and non-stationary noise, often containing transient disturbances which can obscure real signals. Traditional denoising methods, such as principal component analysis and dictionary learning, are not optimal for dealing with this non-Gaussian noise, especially for low signal-to-noise ratio gravitational wave signals. Furthermore, these methods are computationally expensive on large datasets. To overcome these issues, we apply state-of-the-art signal processing techniques, based on recent groundbreaking advancements in deep learning, to denoise gravitational wave signals embedded either in Gaussian noise or in real LIGO noise. We introduce SMTDAE, a Staired Multi-Timestep Denoising Autoencoder, based on sequence-to-sequence bi-directional Long-Short-Term-Memory recurrent neural networks. We demonstrate the advantages of using our unsupervised deep learning approach and show that, after training only using simulated Gaussian noise, SMTDAE achieves superior recovery performance for gravitational wave signals embedded in real non-Gaussian LIGO noise.

preprint2017arXiv

Search for High-energy Neutrinos from Binary Neutron Star Merger GW170817 with ANTARES, IceCube, and the Pierre Auger Observatory

The Advanced LIGO and Advanced Virgo observatories recently discovered gravitational waves from a binary neutron star inspiral. A short gamma-ray burst (GRB) that followed the merger of this binary was also recorded by the Fermi Gamma-ray Burst Monitor (Fermi-GBM), and the Anticoincidence Shield for the Spectrometer for the International Gamma-Ray Astrophysics Laboratory (INTEGRAL), indicating particle acceleration by the source. The precise location of the event was determined by optical detections of emission following the merger. We searched for high-energy neutrinos from the merger in the GeV--EeV energy range using the ANTARES, IceCube, and Pierre Auger Observatories. No neutrinos directionally coincident with the source were detected within $\pm500$ s around the merger time. Additionally, no MeV neutrino burst signal was detected coincident with the merger. We further carried out an extended search in the direction of the source for high-energy neutrinos within the 14-day period following the merger, but found no evidence of emission. We used these results to probe dissipation mechanisms in relativistic outflows driven by the binary neutron star merger. The non-detection is consistent with model predictions of short GRBs observed at a large off-axis angle.

preprint2016arXiv

A First Targeted Search for Gravitational-Wave Bursts from Core-Collapse Supernovae in Data of First-Generation Laser Interferometer Detectors

We present results from a search for gravitational-wave bursts coincident with a set of two core-collapse supernovae observed between 2007 and 2011. We employ data from the Laser Interferometer Gravitational-wave Observatory (LIGO), the Virgo gravitational-wave observatory, and the GEO 600 gravitational-wave observatory. The targeted core-collapse supernovae were selected on the basis of (1) proximity (within approximately 15 Mpc), (2) tightness of observational constraints on the time of core collapse that defines the gravitational-wave search window, and (3) coincident operation of at least two interferometers at the time of core collapse. We find no plausible gravitational-wave candidates. We present the probability of detecting signals from both astrophysically well-motivated and more speculative gravitational-wave emission mechanisms as a function of distance from Earth, and discuss the implications for the detection of gravitational waves from core-collapse supernovae by the upgraded Advanced LIGO and Virgo detectors.

preprint2016arXiv

All-sky search for long-duration gravitational wave transients with LIGO

We present the results of a search for long-duration gravitational wave transients in two sets of data collected by the LIGO Hanford and LIGO Livingston detectors between November 5, 2005 and September 30, 2007, and July 7, 2009 and October 20, 2010, with a total observational time of 283.0 days and 132.9 days, respectively. The search targets gravitational wave transients of duration 10 - 500 s in a frequency band of 40 - 1000 Hz, with minimal assumptions about the signal waveform, polarization, source direction, or time of occurrence. All candidate triggers were consistent with the expected background; as a result we set 90% confidence upper limits on the rate of long-duration gravitational wave transients for different types of gravitational wave signals. For signals from black hole accretion disk instabilities, we set upper limits on the source rate density between $3.4 \times 10^{-5}$ - $9.4 \times 10^{-4}$ Mpc$^{-3}$ yr$^{-1}$ at 90% confidence. These are the first results from an all-sky search for unmodeled long-duration transient gravitational waves.

preprint2016arXiv

An improved analysis of GW150914 using a fully spin-precessing waveform model

This paper presents updated estimates of source parameters for GW150914, a binary black-hole coalescence event detected by the Laser Interferometer Gravitational-wave Observatory (LIGO) on September 14, 2015 [1]. Reference presented parameter estimation [2] of the source using a 13-dimensional, phenomenological precessing-spin model (precessing IMRPhenom) and a 11-dimensional nonprecessing effective-one-body (EOB) model calibrated to numerical-relativity simulations, which forces spin alignment (nonprecessing EOBNR). Here we present new results that include a 15-dimensional precessing-spin waveform model (precessing EOBNR) developed within the EOB formalism. We find good agreement with the parameters estimated previously [2], and we quote updated component masses of $35^{+5}_{-3}\mathrm{M}_\odot$ and $30^{+3}_{-4}\mathrm{M}_\odot$ (where errors correspond to 90% symmetric credible intervals). We also present slightly tighter constraints on the dimensionless spin magnitudes of the two black holes, with a primary spin estimate $0.65$ and a secondary spin estimate $0.75$ at 90% probability. Reference [2] estimated the systematic parameter-extraction errors due to waveform-model uncertainty by combining the posterior probability densities of precessing IMRPhenom and nonprecessing EOBNR. Here we find that the two precessing-spin models are in closer agreement, suggesting that these systematic errors are smaller than previously quoted.

preprint2016arXiv

Comprehensive All-sky Search for Periodic Gravitational Waves in the Sixth Science Run LIGO Data

We report on a comprehensive all-sky search for periodic gravitational waves in the frequency band 100-1500 Hz and with a frequency time derivative in the range of $[-1.18, +1.00]\times 10^{-8}$ Hz/s. Such a signal could be produced by a nearby spinning and slightly non-axisymmetric isolated neutron star in our galaxy. This search uses the data from the Initial LIGO sixth science run and covers a larger parameter space with respect to any past search. A Loosely Coherent detection pipeline was applied to follow up weak outliers in both Gaussian (95% recovery rate) and non-Gaussian (75% recovery rate) bands. No gravitational wave signals were observed, and upper limits were placed on their strength. Our smallest upper limit on worst-case (linearly polarized) strain amplitude $h_0$ is ${9.7}\times 10^{-25}$ near 169 Hz, while at the high end of our frequency range we achieve a worst-case upper limit of ${5.5}\times 10^{-24}$. Both cases refer to all sky locations and entire range of frequency derivative values.

preprint2016arXiv

Directly comparing GW150914 with numerical solutions of Einstein's equations for binary black hole coalescence

We compare GW150914 directly to simulations of coalescing binary black holes in full general relativity, accounting for all the spin-weighted quadrupolar modes, and separately accounting for all the quadrupolar and octopolar modes. Consistent with the posterior distributions reported in LVC_PE[1] (at 90% confidence), we find the data are compatible with a wide range of nonprecessing and precessing simulations. Followup simulations performed using previously-estimated binary parameters most resemble the data. Comparisons including only the quadrupolar modes constrain the total redshifted mass Mz \in [64 - 82M_\odot], mass ratio q = m2/m1 \in [0.6,1], and effective aligned spin χ_eff \in [-0.3, 0.2], where χ_{eff} = (S1/m1 + S2/m2) \cdot\hat{L} /M. Including both quadrupolar and octopolar modes, we find the mass ratio is even more tightly constrained. Simulations with extreme mass ratios and effective spins are highly inconsistent with the data, at any mass. Several nonprecessing and precessing simulations with similar mass ratio and χ_{eff} are consistent with the data. Though correlated, the components' spins (both in magnitude and directions) are not significantly constrained by the data. For nonprecessing binaries, interpolating between simulations, we reconstruct a posterior distribution consistent with previous results. The final black hole's redshifted mass is consistent with Mf,z between 64.0 - 73.5M_\odot and the final black hole's dimensionless spin parameter is consistent with af = 0.62 - 0.73. As our approach invokes no intermediate approximations to general relativity and can strongly reject binaries whose radiation is inconsistent with the data, our analysis provides a valuable complement to LVC_PE[1].

preprint2016arXiv

GW150914: First results from the search for binary black hole coalescence with Advanced LIGO

On September 14, 2015 at 09:50:45 UTC the two detectors of the Laser Interferometer Gravitational-wave Observatory (LIGO) simultaneously observed the binary black hole merger GW150914. We report the results of a matched-filter search using relativistic models of compact-object binaries that recovered GW150914 as the most significant event during the coincident observations between the two LIGO detectors from September 12 to October 20, 2015. GW150914 was observed with a matched filter signal-to-noise ratio of 24 and a false alarm rate estimated to be less than 1 event per 203 000 years, equivalent to a significance greater than 5.1 σ.

preprint2016arXiv

Localization and broadband follow-up of the gravitational-wave transient GW150914

A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared with 63 teams of observers covering radio, optical, near-infrared, X-ray, and gamma-ray wavelengths with ground- and space-based facilities. In this Letter we describe the low-latency analysis of the GW data and present the sky localization of the first observed compact binary merger. We summarize the follow-up observations reported by 25 teams via private Gamma-ray Coordinates Network circulars, giving an overview of the participating facilities, the GW sky localization coverage, the timeline and depth of the observations. As this event turned out to be a binary black hole merger, there is little expectation of a detectable electromagnetic (EM) signature. Nevertheless, this first broadband campaign to search for a counterpart of an Advanced LIGO source represents a milestone and highlights the broad capabilities of the transient astronomy community and the observing strategies that have been developed to pursue neutron star binary merger events. Detailed investigations of the EM data and results of the EM follow-up campaign are being disseminated in papers by the individual teams.

preprint2016arXiv

Properties of the Binary Black Hole Merger GW150914

On September 14, 2015, the Laser Interferometer Gravitational-wave Observatory (LIGO) detected a gravitational-wave transient (GW150914); we characterize the properties of the source and its parameters. The data around the time of the event were analyzed coherently across the LIGO network using a suite of accurate waveform models that describe gravitational waves from a compact binary system in general relativity. GW150914 was produced by a nearly equal mass binary black hole of $36^{+5}_{-4} M_\odot$ and $29^{+4}_{-4} M_\odot$; for each parameter we report the median value and the range of the 90% credible interval. The dimensionless spin magnitude of the more massive black hole is bound to be $<0.7$ (at 90% probability). The luminosity distance to the source is $410^{+160}_{-180}$ Mpc, corresponding to a redshift $0.09^{+0.03}_{-0.04}$ assuming standard cosmology. The source location is constrained to an annulus section of $610$ deg$^2$, primarily in the southern hemisphere. The binary merges into a black hole of $62^{+4}_{-4} M_\odot$ and spin $0.67^{+0.05}_{-0.07}$. This black hole is significantly more massive than any other inferred from electromagnetic observations in the stellar-mass regime.

preprint2016arXiv

Search for transient gravitational waves in coincidence with short duration radio transients during 2007-2013

We present an archival search for transient gravitational-wave bursts in coincidence with 27 single pulse triggers from Green Bank Telescope pulsar surveys, using the LIGO, Virgo and GEO interferometer network. We also discuss a check for gravitational-wave signals in coincidence with Parkes Fast Radio Bursts using similar methods. Data analyzed in these searches were collected between 2007 and 2013. Possible sources of emission of both short-duration radio signals and transient gravitational-wave emission include starquakes on neutron stars, binary coalescence of neutron stars, and cosmic string cusps. While no evidence for gravitational-wave emission in coincidence with these radio transients was found, the current analysis serves as a prototype for similar future searches using more sensitive second-generation interferometers.

preprint2016arXiv

Supplement: Localization and broadband follow-up of the gravitational-wave transient GW150914

This Supplement provides supporting material for arXiv:1602.08492 . We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the different bands.

preprint2016arXiv

Supplement: The Rate of Binary Black Hole Mergers Inferred from Advanced LIGO Observations Surrounding GW150914

Supplemental information for a Letter reporting the rate of binary black hole (BBH) coalescences inferred from 16 days of coincident Advanced LIGO observations surrounding the transient gravitational wave signal GW150914. In that work we reported various rate estimates whose 90\% credible intervals fell in the range $2$--$600 \, \mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$. Here we give details of our method and computations, including information about our search pipelines, a derivation of our likelihood function for the analysis, a description of the astrophysical search trigger distribution expected from merging BBHs, details on our computational methods, a description of the effects and our model for calibration uncertainty, and an analytic method of estimating our detector sensitivity that is calibrated to our measurements.

preprint2016arXiv

The Rate of Binary Black Hole Mergers Inferred from Advanced LIGO Observations Surrounding GW150914

A transient gravitational-wave signal, GW150914, was identified in the twin Advanced LIGO detectors on September 14, 2015 at 09:50:45 UTC. To assess the implications of this discovery, the detectors remained in operation with unchanged configurations over a period of 39 d around the time of the signal. At the detection statistic threshold corresponding to that observed for GW150914, our search of the 16 days of simultaneous two-detector observational data is estimated to have a false alarm rate (FAR) of $< 4.9 \times 10^{-6} \, \mathrm{yr}^{-1}$, yielding a $p$-value for GW150914 of $< 2 \times 10^{-7}$. Parameter estimation followup on this trigger identifies its source as a binary black hole (BBH) merger with component masses $(m_1, m_2) = \left(36^{+5}_{-4},29^{+4}_{-4}\right) \, M_\odot$ at redshift $z = 0.09^{+0.03}_{-0.04}$ (median and 90\% credible range). Here we report on the constraints these observations place on the rate of BBH coalescences. Considering only GW150914, assuming that all BBHs in the Universe have the same masses and spins as this event, imposing a search FAR threshold of 1 per 100 years, and assuming that the BBH merger rate is constant in the comoving frame, we infer a 90% credible range of merger rates between $2$--$53 \, \mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$ (comoving frame). Incorporating all search triggers that pass a much lower threshold while accounting for the uncertainty in the astrophysical origin of each trigger, we estimate a higher rate, ranging from $13$--$600 \, \mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$ depending on assumptions about the BBH mass distribution. All together, our various rate estimates fall in the conservative range $2$--$600 \, \mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$.

preprint2016arXiv

Upper limits on the rates of binary neutron star and neutron-star--black-hole mergers from Advanced LIGO's first observing run

We report here the non-detection of gravitational waves from the merger of binary neutron star systems and neutron-star--black-hole systems during the first observing run of Advanced LIGO. In particular we searched for gravitational wave signals from binary neutron star systems with component masses $\in [1,3] M_{\odot}$ and component dimensionless spins $< 0.05$. We also searched for neutron-star--black-hole systems with the same neutron star parameters, black hole mass $\in [2,99] M_{\odot}$ and no restriction on the black hole spin magnitude. We assess the sensitivity of the two LIGO detectors to these systems, and find that they could have detected the merger of binary neutron star systems with component mass distributions of $1.35\pm0.13 M_{\odot}$ at a volume-weighted average distance of $\sim$ 70Mpc, and for neutron-star--black-hole systems with neutron star masses of $1.4M_\odot$ and black hole masses of at least $5M_\odot$, a volume-weighted average distance of at least $\sim$ 110Mpc. From this we constrain with 90% confidence the merger rate to be less than 12,600 Gpc$^{-3}$yr$^{-1}$ for binary-neutron star systems and less than 3,600 Gpc$^{-3}$yr$^{-1}$ for neutron-star--black-hole systems. We find that if no detection of neutron-star binary mergers is made in the next two Advanced LIGO and Advanced Virgo observing runs we would place significant constraints on the merger rates. Finally, assuming a rate of $10^{+20}_{-7}$Gpc$^{-3}$yr$^{-1}$ short gamma ray bursts beamed towards the Earth and assuming that all short gamma-ray bursts have binary-neutron-star (neutron-star--black-hole) progenitors we can use our 90% confidence rate upper limits to constrain the beaming angle of the gamma-ray burst to be greater than ${2.3^{+1.7}_{-1.1}}^{\circ}$ (${4.3^{+3.1}_{-1.9}}^{\circ}$).

preprint2015arXiv

A search of the Orion spur for continuous gravitational waves using a "loosely coherent" algorithm on data from LIGO interferometers

We report results of a wideband search for periodic gravitational waves from isolated neutron stars within the Orion spur towards both the inner and outer regions of our Galaxy. As gravitational waves interact very weakly with matter, the search is unimpeded by dust and concentrations of stars. One search disk (A) is $6.87^\circ$ in diameter and centered on $20^\textrm{h}10^\textrm{m}54.71^\textrm{s}+33^\circ33'25.29"$, and the other (B) is $7.45^\circ$ in diameter and centered on $8^\textrm{h}35^\textrm{m}20.61^\textrm{s}-46^\circ49'25.151"$. We explored the frequency range of 50-1500 Hz and frequency derivative from $0$ to $-5\times 10^{-9}$ Hz/s. A multi-stage, loosely coherent search program allowed probing more deeply than before in these two regions, while increasing coherence length with every stage. Rigorous followup parameters have winnowed initial coincidence set to only 70 candidates, to be examined manually. None of those 70 candidates proved to be consistent with an isolated gravitational wave emitter, and 95% confidence level upper limits were placed on continuous-wave strain amplitudes. Near $169$ Hz we achieve our lowest 95% CL upper limit on worst-case linearly polarized strain amplitude $h_0$ of $6.3\times 10^{-25}$, while at the high end of our frequency range we achieve a worst-case upper limit of $3.4\times 10^{-24}$ for all polarizations and sky locations.

preprint2015arXiv

Detecting eccentric supermassive black hole binaries with pulsar timing arrays: Resolvable source strategies

The couplings between supermassive black-hole binaries and their environments within galactic nuclei have been well studied as part of the search for solutions to the final parsec problem. The scattering of stars by the binary or the interaction with a circumbinary disk may efficiently drive the system to sub-parsec separations, allowing the binary to enter a regime where the emission of gravitational waves can drive it to merger within a Hubble time. However, these interactions can also affect the orbital parameters of the binary. In particular, they may drive an increase in binary eccentricity which survives until the system's gravitational-wave signal enters the pulsar-timing array band. Therefore, if we can measure the eccentricity from observed signals, we can potentially deduce some of the properties of the binary environment. To this end, we build on previous techniques to present a general Bayesian pipeline with which we can detect and estimate the parameters of an eccentric supermassive black-hole binary system with pulsar-timing arrays. Additionally, we generalize the pulsar-timing array $\mathcal{F}_e$-statistic to eccentric systems, and show that both this statistic and the Bayesian pipeline are robust when studying circular or arbitrarily eccentric systems. We explore how eccentricity influences the detection prospects of single gravitational-wave sources, as well as the detection penalty incurred by employing a circular waveform template to search for eccentric signals, and conclude by identifying important avenues for future study.

preprint2015arXiv

Detection of eccentric supermassive black hole binaries with pulsar timing arrays: Signal-to-noise ratio calculations

We present a detailed analysis of the expected signal-to-noise ratios of supermassive black hole binaries on eccentric orbits observed by pulsar timing arrays. We derive several analytical relations that extend the results of Peters and Mathews [Phys. Rev. D 131, 435 (1963)] to quantify the impact of eccentricity in the detection of single resolvable binaries in the pulsar timing array band. We present ready-to-use expressions to compute the increase/loss in signal-to-noise ratio of eccentric single resolvable sources whose dominant harmonic is located in the low/high frequency sensitivity regime of pulsar timing arrays. Building upon the work of Phinney (arXiv:astro-ph/0108028) and Enoki and Nagashima [Prog. Theor. Phys. 117, 241 (2007)], we present an analytical framework that enables the construction of rapid spectra for a stochastic gravitational-wave background generated by a cosmological population of eccentric sources. We confirm previous findings which indicate that, relative to a population of quasicircular binaries, the strain of a stochastic, isotropic gravitational-wave background generated by a cosmological population of eccentric binaries will be suppressed in the frequency band of pulsar timing arrays. We quantify this effect in terms of signal-to-noise ratios in a pulsar timing array.

preprint2014arXiv

Accurate and efficient waveforms for compact binaries on eccentric orbits

Compact binaries that emit gravitational waves in the sensitivity band of ground-based detectors can have non-negligible eccentricities just prior to merger, depending on the formation scenario. We develop a purely analytic, frequency-domain model for gravitational waves emitted by compact binaries on orbits with small eccentricity, which reduces to the quasi-circular post-Newtonian approximant TaylorF2 at zero eccentricity and to the post-circular approximation of Yunes et al. (2009) at small eccentricity. Our model uses a spectral approximation to the (post-Newtonian) Kepler problem to model the orbital phase as a function of frequency, accounting for eccentricity effects up to ${\cal{O}}(e^8)$ at each post-Newtonian order. Our approach accurately reproduces an alternative time-domain eccentric waveform model for eccentricities $e\in [0, 0.4]$ and binaries with total mass less than 12 solar masses. As an application, we evaluate the signal amplitude that eccentric binaries produce in different networks of existing and forthcoming gravitational waves detectors. Assuming a population of eccentric systems containing black holes and neutron stars that is uniformly distributed in co-moving volume, we estimate that second generation detectors like Advanced LIGO could detect approximately 0.1-10 events per year out to redshift $z\sim 0.2$, while an array of Einstein Telescope detectors could detect hundreds of events per year to redshift $z \sim 2.3$.

preprint2014arXiv

Application of a Hough search for continuous gravitational waves on data from the 5th LIGO science run

We report on an all-sky search for periodic gravitational waves in the frequency range $\mathrm{50-1000 Hz}$ with the first derivative of frequency in the range $-8.9 \times 10^{-10}$ Hz/s to zero in two years of data collected during LIGO's fifth science run. Our results employ a Hough transform technique, introducing a $χ^2$ test and analysis of coincidences between the signal levels in years 1 and 2 of observations that offers a significant improvement in the product of strain sensitivity with compute cycles per data sample compared to previously published searches. Since our search yields no surviving candidates, we present results taking the form of frequency dependent, 95$%$ confidence upper limits on the strain amplitude $h_0$. The most stringent upper limit from year 1 is $1.0\times 10^{-24}$ in the $\mathrm{158.00-158.25 Hz}$ band. In year 2, the most stringent upper limit is $\mathrm{8.9\times10^{-25}}$ in the $\mathrm{146.50-146.75 Hz}$ band. This improved detection pipeline, which is computationally efficient by at least two orders of magnitude better than our flagship Einstein$@$Home search, will be important for "quick-look" searches in the Advanced LIGO and Virgo detector era.

preprint2014arXiv

Characterization of the LIGO detectors during their sixth science run

In 2009-2010, the Laser Interferometer Gravitational-wave Observa- tory (LIGO) operated together with international partners Virgo and GEO600 as a network to search for gravitational waves of astrophysical origin. The sensitiv- ity of these detectors was limited by a combination of noise sources inherent to the instrumental design and its environment, often localized in time or frequency, that couple into the gravitational-wave readout. Here we review the performance of the LIGO instruments during this epoch, the work done to characterize the de- tectors and their data, and the effect that transient and continuous noise artefacts have on the sensitivity of LIGO to a variety of astrophysical sources.

preprint2014arXiv

Constraints on cosmic strings from the LIGO-Virgo gravitational-wave detectors

Cosmic strings can give rise to a large variety of interesting astrophysical phenomena. Among them, powerful bursts of gravitational waves (GWs) produced by cusps are a promising observational signature. In this Letter we present a search for GWs from cosmic string cusps in data collected by the LIGO and Virgo gravitational wave detectors between 2005 and 2010, with over 625 days of live time. We find no evidence of GW signals from cosmic strings. From this result, we derive new constraints on cosmic string parameters, which complement and improve existing limits from previous searches for a stochastic background of GWs from cosmic microwave background measurements and pulsar timing data. In particular, if the size of loops is given by the gravitational backreaction scale, we place upper limits on the string tension $Gμ$ below $10^{-8}$ in some regions of the cosmic string parameter space.

preprint2014arXiv

Gravitational waves from known pulsars: results from the initial detector era

We present the results of searches for gravitational waves from a large selection of pulsars using data from the most recent science runs (S6, VSR2 and VSR4) of the initial generation of interferometric gravitational wave detectors LIGO (Laser Interferometric Gravitational-wave Observatory) and Virgo. We do not see evidence for gravitational wave emission from any of the targeted sources but produce upper limits on the emission amplitude. We highlight the results from seven young pulsars with large spin-down luminosities. We reach within a factor of five of the canonical spin-down limit for all seven of these, whilst for the Crab and Vela pulsars we further surpass their spin-down limits. We present new or updated limits for 172 other pulsars (including both young and millisecond pulsars). Now that the detectors are undergoing major upgrades, and, for completeness, we bring together all of the most up-to-date results from all pulsars searched for during the operations of the first-generation LIGO, Virgo and GEO600 detectors. This gives a total of 195 pulsars including the most recent results described in this paper.

preprint2014arXiv

Searching for stochastic gravitational waves using data from the two co-located LIGO Hanford detectors

Searches for a stochastic gravitational-wave background (SGWB) using terrestrial detectors typically involve cross-correlating data from pairs of detectors. The sensitivity of such cross-correlation analyses depends, among other things, on the separation between the two detectors: the smaller the separation, the better the sensitivity. Hence, a co-located detector pair is more sensitive to a gravitational-wave background than a non-co-located detector pair. However, co-located detectors are also expected to suffer from correlated noise from instrumental and environmental effects that could contaminate the measurement of the background. Hence, methods to identify and mitigate the effects of correlated noise are necessary to achieve the potential increase in sensitivity of co-located detectors. Here we report on the first SGWB analysis using the two LIGO Hanford detectors and address the complications arising from correlated environmental noise. We apply correlated noise identification and mitigation techniques to data taken by the two LIGO Hanford detectors, H1 and H2, during LIGO's fifth science run. At low frequencies, 40 - 460 Hz, we are unable to sufficiently mitigate the correlated noise to a level where we may confidently measure or bound the stochastic gravitational-wave signal. However, at high frequencies, 460-1000 Hz, these techniques are sufficient to set a $95%$ confidence level (C.L.) upper limit on the gravitational-wave energy density of Ω(f)<7.7 x 10^{-4} (f/ 900 Hz)^3, which improves on the previous upper limit by a factor of $\sim 180$. In doing so, we demonstrate techniques that will be useful for future searches using advanced detectors, where correlated noise (e.g., from global magnetic fields) may affect even widely separated detectors.

preprint2014arXiv

Self-forced evolutions of an implicit rotating source: A natural framework to model comparable and intermediate mass-ratio systems from inspiral through ringdown

We develop a waveform model to describe the inspiral, merger and ringdown of binary systems with comparable and intermediate mass-ratios. This model incorporates first-order conservative self-force corrections to the energy and angular momentum, which are valid in the strong-field regime [1]. We model the radiative part of the self-force by deriving second-order radiative corrections to the energy flux. These corrections are obtained by minimizing the phase discrepancy between our self-force model and the effective one body model [2, 3] for a variety of mass-ratios. We show that our model performs substantially better than post-Newtonian approximants currently used to model neutron star-black hole mergers from early inspiral to the innermost stable circular orbit. In order to match the late inspiral evolution onto the plunge regime, we extend the 'transition phase' developed by Ori and Thorne [4] by including finite mass-ratio corrections and modelling the orbital phase evolution using an implicit rotating source [5]. We explicitly show that the implicit rotating source approach provides a natural transition from late-time radiation to ringdown that is equivalent to ringdown waveform modelling based on a sum of quasinormal modes.

preprint2013arXiv

A directed search for continuous Gravitational Waves from the Galactic Center

We present the results of a directed search for continuous gravitational waves from unknown, isolated neutron stars in the Galactic Center region, performed on two years of data from LIGO's fifth science run from two LIGO detectors. The search uses a semi-coherent approach, analyzing coherently 630 segments, each spanning 11.5 hours, and then incoherently combining the results of the single segments. It covers gravitational wave frequencies in a range from 78 to 496 Hz and a frequency-dependent range of first order spindown values down to -7.86 x 10^-8 Hz/s at the highest frequency. No gravitational waves were detected. We place 90% confidence upper limits on the gravitational wave amplitude of sources at the Galactic Center. Placing 90% confidence upper limits on the gravitational wave amplitude of sources at the Galactic Center, we reach ~3.35x10^-25 for frequencies near 150 Hz. These upper limits are the most constraining to date for a large-parameter-space search for continuous gravitational wave signals.

preprint2013arXiv

Effect of eccentricity on binary neutron star searches in Advanced LIGO

Binary neutron stars (BNSs) are the primary source of gravitational waves for the Laser Interferometer Gravitational-Wave Observatory (LIGO) and its international partners Virgo and KAGRA. Current BNS searches target field binaries whose orbits will have circularized by radiation reaction before their gravitational waves enter the Advanced LIGO sensitive band at 15 Hz. It has been suggested that a population of BNSs may form by n-body interactions near supermassive black holes or in globular clusters and that these systems may have non-negligible eccentricity in the Advanced LIGO band. We show that for BNS systems with total mass of 2.4 (6.0) solar masses, the effect of eccentricity e < 0.02 (0.05) is negligible and a circular search is effectual for these binaries. For eccentricities up to e = 0.4, we investigate the selection bias caused by neglecting eccentricity in BNS searches. If such high eccentricity systems exist, searches that specifically target eccentric binaries will be needed in Advanced LIGO and Virgo.

preprint2013arXiv

First Searches for Optical Counterparts to Gravitational-wave Candidate Events

During the LIGO and Virgo joint science runs in 2009-2010, gravitational wave (GW) data from three interferometer detectors were analyzed within minutes to select GW candidate events and infer their apparent sky positions. Target coordinates were transmitted to several telescopes for follow-up observations aimed at the detection of an associated optical transient. Images were obtained for eight such GW candidates. We present the methods used to analyze the image data as well as the transient search results. No optical transient was identified with a convincing association with any of these candidates, and none of the GW triggers showed strong evidence for being astrophysical in nature. We compare the sensitivities of these observations to several model light curves from possible sources of interest, and discuss prospects for future joint GW-optical observations of this type.

preprint2013arXiv

Search for long-lived gravitational-wave transients coincident with long gamma-ray bursts

Long gamma-ray bursts (GRBs) have been linked to extreme core-collapse supernovae from massive stars. Gravitational waves (GW) offer a probe of the physics behind long GRBs. We investigate models of long-lived (~10-1000s) GW emission associated with the accretion disk of a collapsed star or with its protoneutron star remnant. Using data from LIGO's fifth science run, and GRB triggers from the swift experiment, we perform a search for unmodeled long-lived GW transients. Finding no evidence of GW emission, we place 90% confidence level upper limits on the GW fluence at Earth from long GRBs for three waveforms inspired by a model of GWs from accretion disk instabilities. These limits range from F<3.5 ergs cm^-2 to $F<1200 ergs cm^-2, depending on the GRB and on the model, allowing us to probe optimistic scenarios of GW production out to distances as far as ~33 Mpc. Advanced detectors are expected to achieve strain sensitivities 10x better than initial LIGO, potentially allowing us to probe the engines of the nearest long GRBs.

preprint2012arXiv

Accurate modeling of intermediate-mass-ratio inspirals: Exploring the form of the self-force in the intermediate-mass-ratio regime

In this paper we develop a waveform model that accurately reproduces the dynamical evolution of intermediate-mass-ratio inspirals, as predicted by the effective-one-body (EOB) model introduced in [1], and which enables us to shed some light on the form of the self-force for events with mass-ratio 1:6, 1:10 and 1:100. To complement this study, we make use of self-force results in the extreme-mass-ratio regime, and of predictions of the EOB model introduced in [1], to derive a prescription for the shift of the orbital frequency at the innermost stable circular orbit which consistently captures predictions from the extreme, intermediate and comparable mass-ratio regimes.

preprint2012arXiv

Importance of including small body spin effects in the modelling of intermediate mass-ratio inspirals. II Accurate parameter extraction of strong sources using higher-order spin effects

We improve the numerical kludge waveform model introduced in [1] in two ways. We extend the equations of motion for spinning black hole binaries derived by Saijo et al. [2] using spin-orbit and spin-spin couplings taken from perturbative and post-Newtonian (PN) calculations at the highest order available. We also include first-order conservative self-force corrections for spin-orbit and spin-spin couplings, which are derived by comparison to PN results. We generate the inspiral evolution using fluxes that include the most recent calculations of small body spin corrections, spin-spin and spin-orbit couplings and higher-order fits to solutions of the Teukolsky equation. Using a simplified version of this model in [1], we found that small body spin effects could be measured through gravitational wave observations from intermediate-mass ratio inspirals (IMRIs) with mass ratio eta ~ 0.001, when both binary components are rapidly rotating. In this paper we study in detail how the spin of the small/big body affects parameter measurement using a variety of mass and spin combinations for typical IMRIs sources. We find that for IMRI events of a moderately rotating intermediate mass black hole (IMBH) of ten thousand solar masses, and a rapidly rotating central supermassive black hole (SMBH) of one million solar masses, gravitational wave observations made with LISA at a fixed signal-to-noise ratio (SNR) of 1000 will be able to determine the inspiralling IMBH mass, the central SMBH mass, the SMBH spin magnitude, and the IMBH spin magnitude to within fractional errors of ~0.001, 0.001, 0.0001, and 9%, respectively. LISA can also determine the location of the source in the sky and the SMBH spin orientation to within ~0.0001 steradians. We show that by including conservative corrections up to 2.5PN order, systematic errors no longer dominate over statistical errors for IMRIs with typical SNR ~1000.

preprint2011arXiv

Importance of including small body spin effects in the modelling of extreme and intermediate mass-ratio inspirals

We explore the ability of future low-frequency gravitational wave detectors to measure the spin of stellar mass and intermediate mass black holes that inspiral onto super-massive Kerr black holes (SMBHs). We develop a kludge waveform model based on the equations of motion derived by Saijo et al. [Phys Rev D 58, 064005, 1998] for spinning BH binaries, augmented with spin-orbit and spin-spin couplings taken from perturbative and post-Newtonian (PN) calculations, and the associated conservative self-force corrections, derived by comparison to PN results. We model the inspiral phase using accurate fluxes which include perturbative corrections for the spin of the inspiralling body, spin-spin couplings and higher-order fits to solutions of the Teukolsky equation. We present results of Monte Carlo simulations of parameter estimation errors and of the model errors that arise when we omit conservative corrections from the waveform template. For a source 5000+10^6 solar mass observed with an SNR of 1000, LISA will be able to determine the two masses to within a fractional error of ~0.001, measure the SMBH spin magnitude, q, and the spin magnitude of the inspiralling BH to 0.0001, 10%, respectively, and determine the location of the source in the sky and the SMBH spin orientation to within 0.0001 steradians. For a 10+10^6 solar mass system observed with SNR of 30, LISA will not be able to determine the spin magnitude of the inspiralling BH, although the measurement of the other waveform parameters is not significantly degraded by the presence of spin. The model errors which arise from ignoring conservative corrections become significant for mass-ratios above 0.0001, but including these corrections up to 2PN order may be sufficient to reduce these systematic errors to an acceptable level.

preprint2011arXiv

Influence of conservative corrections on parameter estimation for extreme-mass-ratio inspirals

We present an improved numerical kludge waveform model for circular, equatorial EMRIs. The model is based on true Kerr geodesics, augmented by radiative self-force corrections derived from perturbative calculations, and in this paper for the first time we include conservative self-force corrections that we derive by comparison to post-Newtonian results. We present results of a Monte Carlo simulation of parameter estimation errors computed using the Fisher Matrix and also assess the theoretical errors that would arise form omitting the conservative correction terms we include here. We present results for three different types of system, namely the inspirals of black holes, neutron stars or white dwarfs into a supermassive black hole (SMBH). The analysis shows that for a typical source (a 10 solar mass compact object captured by a one million solar mass SMBH at signal to noise ratio of 30) we expect to determine the two masses to within a fractional error of ~0.0001, measure the spin parameter q to ~0.00003 and determine the location of the source on the sky and the spin orientation to within 0.001 steradians. We show that, for this kludge model, omitting the conservative corrections leads to a small error over much of the parameter space, i.e., the ratio R of the theoretical model error to the Fisher Matrix error is R<1 for all ten parameters in the model. For the few systems with larger errors typically R<3 and hence the conservative corrections can be marginally ignored. In addition, we use our model and first order self-force results for Schwarzschild black holes to estimate the error that arises from omitting the second-order radiative piece of the self-force. This indicates that it may not be necessary to go beyond first order to recover accurate parameter estimates.

preprint2011arXiv

Intermediate-mass-ratio-inspirals in the Einstein Telescope: I. Signal-to-noise ratio calculations

The Einstein Telescope (ET) is a proposed third generation ground-based interferometer, for which the target is a sensitivity that is a factor of ten better than Advanced LIGO and a frequency range that extends down to about 1Hz. ET will provide opportunities to test Einstein's theory of relativity in the strong field and will realize precision gravitational wave astronomy with a thousandfold increase in the expected number of events over the advanced ground-based detectors. A design study for ET is currently underway, so it is timely to assess the science that could be done with such an instrument. This paper is the first in a series that will carry out a detailed study of intermediate-mass-ratio inspirals (IMRIs) for ET. In the context of ET, an IMRI is the inspiral of a neutron star or stellar-mass black hole into an intermediate mass black hole (IMBH). In this paper we focus on the development of IMRI waveform models for circular and equatorial inspirals. We consider two approximations for the waveforms, which both incorporate the inspiral, merger and ringdown phases in a consistent way. One approximation, valid for IMBHs of arbitrary spin, uses the transition model of Ori and Thorne [1] to describe the merger, and this is then matched smoothly onto a ringdown waveform. The second approximation uses the Effective One Body (EOB) approach to model the merger phase of the waveform and is valid for non-spinning IMBHs. In this paper, we use both waveform models to compute signal-to-noise ratios (SNRs) for IMRI sources detectable by ET. At a redshift of z=1, we find typical SNRs for IMRI systems with masses 1.4+100 solar masses, 10+100 solar masses, 1.4+500 solar masses and 10+500 solar masses of about 10-25, 40-80, 3-15 and 10-60, respectively. We also find that the two models make predictions for non-spinning inspirals that are consistent to about ten percent.

preprint2011arXiv

Intermediate-mass-ratio-inspirals in the Einstein Telescope. II. Parameter estimation errors

We explore the precision with which the Einstein Telescope (ET) will be able to measure the parameters of intermediate-mass-ratio inspirals (IMRIs). We calculate the parameter estimation errors using the Fisher Matrix formalism and present results of a Monte Carlo simulation of these errors over choices for the extrinsic parameters of the source. These results are obtained using two different models for the gravitational waveform which were introduced in paper I of this series. These two waveform models include the inspiral, merger and ringdown phases in a consistent way. One of the models, based on the transition scheme of Ori & Thorne [1], is valid for IMBHs of arbitrary spin, whereas the second model, based on the Effective One Body (EOB) approach, has been developed to cross-check our results in the non-spinning limit. In paper I of this series, we demonstrated the excellent agreement in both phase and amplitude between these two models for non-spinning black holes, and that their predictions for signal-to-noise ratios (SNRs) are consistent to within ten percent. We now use these models to estimate parameter estimation errors for binary systems with masses 1.4+100, 10+100, 1.4+500 and 10+500 solar masses (SMs), and various choices for the spin of the central intermediate-mass black hole (IMBH). Assuming a detector network of three ETs, the analysis shows that for a 10 SM compact object (CO) inspiralling into a 100 SM IMBH with spin q=0.3, detected with an SNR of 30, we should be able to determine the CO and IMBH masses, and the IMBH spin magnitude to fractional accuracies of 0.001, 0.0003, and 0.001, respectively. We also expect to determine the location of the source in the sky and the luminosity distance to within 0.003 steradians, and 10%, respectively. We also assess how the precision of parameter determination depends on the network configuration.

preprint2010arXiv

On the relevance of gravitational self-force corrections on parameter estimation errors for extreme-mass-ratio inspirals

It is not currently clear how important it will be to include conservative self-force (SF) corrections in the models for extreme-mass-ratio inspiral (EMRI) waveforms that will be used to detect such signals in LISA (Laser Interferometer Space Antenna) data. These proceedings will address this issue for circular-equatorial inspirals using an approximate EMRI model that includes conservative corrections at leading post-Newtonian order. We will present estimates of the magnitude of the parameter estimation errors that would result from omitting conservative corrections, and compare these to the errors that will arise from noise fluctuations in the detector. We will also use this model to explore the relative importance of the second-order radiative piece of the SF, which is not presently known.

E. A. Huerta

What is connected

Connect this record

See the researcher in context

Building this map preview

69 published item(s)

FAIR AI Models in High Energy Physics

A FAIR and AI-ready Higgs boson decay dataset

Constraints from LIGO O3 data on gravitational-wave emission due to r-modes in the glitching pulsar PSR J0537-6910

Inference-optimized AI and high performance computing for gravitational wave detection at scale

Interpretable AI forecasting for numerical relativity waveforms of quasi-circular, spinning, non-precessing binary black hole mergers

Interpreting a Machine Learning Model for Detecting Gravitational Waves

Search for anisotropic gravitational-wave backgrounds using data from Advanced LIGO and Advanced Virgo's first three observing runs

Advances in Machine and Deep Learning for Modeling and Real-time Detection of Multi-Messenger Sources

Open data from the first and second observing runs of Advanced LIGO and Advanced Virgo

Search for lensing signatures in the gravitational-wave observations from the first half of LIGO-Virgo's third observing run

A guide to LIGO-Virgo detector noise and extraction of transient gravitational-wave signals

A Joint Fermi-GBM and LIGO/Virgo Analysis of Compact Binary Mergers From the First and Second Gravitational-wave Observing Runs

Artificial neural network subgrid models of 2-D compressible magnetohydrodynamic turbulence

Binary Black Hole Population Properties Inferred from the First and Second Observing Runs of Advanced LIGO and Advanced Virgo

Deep Learning for Cardiologist-level Myocardial Infarction Detection in Electrocardiograms

Deep learning for gravitational wave forecasting of neutron star mergers

Deep Transfer Learning for Star Cluster Classification: I. Application to the PHANGS-HST Survey

GW190412: Observation of a Binary-Black-Hole Coalescence with Asymmetric Masses

GW190425: Observation of a Compact Binary Coalescence with Total Mass $\sim 3.4 M_{\odot}$

GW190521: A Binary Black Hole Merger with a Total Mass of $150 ~ M_{\odot}$

GW190814: Gravitational Waves from the Coalescence of a 23 M$_\odot$ Black Hole with a 2.6 M$_\odot$ Compact Object

Model comparison from LIGO-Virgo data on GW170817's binary components and consequences for the merger remnant

Modeling the uncertainties of solar-system ephemerides for robust gravitational-wave searches with pulsar timing arrays

Physics-inspired deep learning to characterize the signal manifold of quasi-circular, spinning, non-precessing binary black hole mergers

Properties and astrophysical implications of the 150 Msun binary black hole merger GW190521

Search for Eccentric Binary Black Hole Mergers with Advanced LIGO and Advanced Virgo during their First and Second Observing Runs

Search for intermediate mass black hole binaries in the first and second observing runs of the Advanced LIGO and Virgo network

Searches for Gravitational Waves from Known Pulsars at Two Harmonics in 2015-2017 LIGO Data

An Optically Targeted Search for Gravitational Waves emitted by Core-Collapse Supernovae during the First and Second Observing Runs of Advanced LIGO and Advanced Virgo

Search for sub-solar mass ultracompact binaries in Advanced LIGO's second observing run

The NANOGrav 11-Year Data Set: Evolution of Gravitational Wave Background Statistics

The NANOGrav 11-Year Data Set: Limits on Gravitational Wave Memory

Search for Multi-messenger Sources of Gravitational Waves and High-energy Neutrinos with Advanced LIGO during its first Observing Run, ANTARES and IceCube

Denoising Gravitational Waves using Deep Learning with Recurrent Denoising Autoencoders

Search for High-energy Neutrinos from Binary Neutron Star Merger GW170817 with ANTARES, IceCube, and the Pierre Auger Observatory

A First Targeted Search for Gravitational-Wave Bursts from Core-Collapse Supernovae in Data of First-Generation Laser Interferometer Detectors

All-sky search for long-duration gravitational wave transients with LIGO

An improved analysis of GW150914 using a fully spin-precessing waveform model

Comprehensive All-sky Search for Periodic Gravitational Waves in the Sixth Science Run LIGO Data

Directly comparing GW150914 with numerical solutions of Einstein's equations for binary black hole coalescence

GW150914: First results from the search for binary black hole coalescence with Advanced LIGO

Localization and broadband follow-up of the gravitational-wave transient GW150914

Properties of the Binary Black Hole Merger GW150914

Search for transient gravitational waves in coincidence with short duration radio transients during 2007-2013

Supplement: Localization and broadband follow-up of the gravitational-wave transient GW150914

Supplement: The Rate of Binary Black Hole Mergers Inferred from Advanced LIGO Observations Surrounding GW150914

The Rate of Binary Black Hole Mergers Inferred from Advanced LIGO Observations Surrounding GW150914

Upper limits on the rates of binary neutron star and neutron-star--black-hole mergers from Advanced LIGO's first observing run

A search of the Orion spur for continuous gravitational waves using a "loosely coherent" algorithm on data from LIGO interferometers

Detecting eccentric supermassive black hole binaries with pulsar timing arrays: Resolvable source strategies

Detection of eccentric supermassive black hole binaries with pulsar timing arrays: Signal-to-noise ratio calculations

Accurate and efficient waveforms for compact binaries on eccentric orbits

Application of a Hough search for continuous gravitational waves on data from the 5th LIGO science run

Characterization of the LIGO detectors during their sixth science run

Constraints on cosmic strings from the LIGO-Virgo gravitational-wave detectors

Gravitational waves from known pulsars: results from the initial detector era

Searching for stochastic gravitational waves using data from the two co-located LIGO Hanford detectors

Self-forced evolutions of an implicit rotating source: A natural framework to model comparable and intermediate mass-ratio systems from inspiral through ringdown

A directed search for continuous Gravitational Waves from the Galactic Center

Effect of eccentricity on binary neutron star searches in Advanced LIGO

First Searches for Optical Counterparts to Gravitational-wave Candidate Events

Search for long-lived gravitational-wave transients coincident with long gamma-ray bursts

Accurate modeling of intermediate-mass-ratio inspirals: Exploring the form of the self-force in the intermediate-mass-ratio regime

Importance of including small body spin effects in the modelling of intermediate mass-ratio inspirals. II Accurate parameter extraction of strong sources using higher-order spin effects

Importance of including small body spin effects in the modelling of extreme and intermediate mass-ratio inspirals

Influence of conservative corrections on parameter estimation for extreme-mass-ratio inspirals

Intermediate-mass-ratio-inspirals in the Einstein Telescope: I. Signal-to-noise ratio calculations

Intermediate-mass-ratio-inspirals in the Einstein Telescope. II. Parameter estimation errors

On the relevance of gravitational self-force corrections on parameter estimation errors for extreme-mass-ratio inspirals