BZPEERPapers, people, evidence and research paths
Menu

Discover

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

Research tools

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

Network

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

Opportunities

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

Account

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

Source author record

Eric Thrane

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

ResearcherUnclaimed source record
gr-qcastro-ph.HEastro-ph.IMastro-ph.COastro-ph.GAastro-ph.SRphysics.ins-det

Catalog footprint

What is connected

54works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

54 published item(s)

preprint2022arXiv

A Rosetta Stone for eccentric gravitational waveform models

Orbital eccentricity is a key signature of dynamical binary black hole formation. The gravitational waves from a coalescing binary contain information about its orbital eccentricity, which may be measured if the binary retains sufficient eccentricity near merger. Dedicated waveforms are required to measure eccentricity. Several models have been put forward, and show good agreement with numerical relativity at the level of a few percent or better. However, there are multiple ways to define eccentricity for inspiralling systems, and different models internally use different definitions of eccentricity, making it difficult to directly compare eccentricity measurements. In this work, we systematically compare two eccentric waveform models, $\texttt{SEOBNRE}$ and $\texttt{TEOBResumS}$, by developing a framework to translate between different definitions of eccentricity. This mapping is constructed by minimizing the relative mismatch between the two models over eccentricity and reference frequency, before evolving the eccentricity of one model to the same reference frequency as the other model. We show that for a given value of eccentricity passed to $\texttt{SEOBNRE}$, one must input a $20$-$50\%$ smaller value of eccentricity to $\texttt{TEOBResumS}$ in order to obtain a waveform with the same empirical eccentricity. We verify this mapping by repeating our analysis for eccentric numerical relativity simulations, demonstrating that $\texttt{TEOBResumS}$ reports a correspondingly smaller value of eccentricity than $\texttt{SEOBNRE}$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Building better spin models for merging binary black holes: Evidence for non-spinning and rapidly spinning nearly aligned sub-populations

Recent work paints a conflicting portrait of the distribution of black hole spins in merging binaries measured with gravitational waves. Some analyses find that a significant fraction of merging binaries contain at least one black hole with a spin tilt $>90^\circ$ with respect to the orbital angular momentum vector, which has been interpreted as a signature for dynamical assembly. Other analyses find the data are consistent with a bimodal population in which some binaries contain black holes with negligible spin while the rest contain black holes with spin vectors preferentially aligned with the orbital angular momentum vector. In this work, we scrutinize models for the distribution of black hole spins to pinpoint possible failure modes in which the model yields a faulty conclusion. We reanalyze data from the second LIGO--Virgo gravitational-wave transient catalog (GWTC-2) using a revised spin model, which allows for a sub-population of black holes with negligible spins. In agreement with recent results by Roulet et al., we show that the GWTC-2 detections are consistent with two distinct sub-populations. We estimate that $29-75\%$ (90\% credible interval) of merging binaries contain black holes with negligible spin $χ\approx 0$. The remaining binaries are part of a second sub-population in which the spin vectors are preferentially (but not exactly) aligned to the orbital angular momentum. The black holes in this second sub-population are characterized by spins of $χ\sim0.45$. We suggest that the inferred spin distribution is consistent with the hypothesis that all merging binaries form via the field formation scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Consistency of the Parkes Pulsar Timing Array Signal with a Nanohertz Gravitational-wave Background

Pulsar timing array experiments have recently reported strong evidence for a common-spectrum stochastic process with a strain spectral index consistent with that expected of a nanohertz-frequency gravitational-wave background, but with negligible yet non-zero evidence for spatial correlations required for a definitive detection. However, it was pointed out by the Parkes Pulsar Timing Array (PPTA) collaboration that the same models used in recent analyses resulted in strong evidence for a common-spectrum process in simulations where none is present. In this work, we introduce a methodology to distinguish pulsar power spectra with the same amplitude from noise power spectra of similar but distinct amplitudes. The former is the signature of a spatially uncorrelated pulsar term of a nanohertz gravitational-wave background, whereas the latter could represent ensemble pulsar noise properties. We test the methodology on simulated data sets. We find that the reported common process in PPTA pulsars is indeed consistent with the spectral feature of a pulsar term. We recommend this methodology as one of the validity tests that the real astrophysical and cosmological backgrounds should pass, as well as for inferences about the spatially uncorrelated component of the background.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Fast, flexible, and accurate evaluation of gravitational-wave Malmquist bias with machine learning

Many astronomical surveys are limited by the brightness of the sources, and gravitational-wave searches are no exception. The detectability of gravitational waves from merging binaries is affected by the mass and spin of the constituent compact objects. To perform unbiased inference on the distribution of compact binaries, it is necessary to account for this selection effect, which is known as Malmquist bias. Since systematic error from selection effects grows with the number of events, it will be increasingly important over the coming years to accurately estimate the observational selection function for gravitational-wave astronomy. We employ density estimation methods to accurately and efficiently compute the compact binary coalescence selection function. We introduce a simple pre-processing method, which significantly reduces the complexity of the required machine learning models. We demonstrate that our method has smaller statistical errors at comparable computational cost than the method currently most widely used allowing us to probe narrower distributions of spin magnitudes. The currently used method leaves $10-50\%$ of the interesting black hole spin models inaccessible; our new method can probe $>99\%$ of the models and has a lower uncertainty for $>80\%$ of the models.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Linking the rates of neutron star binaries and short gamma-ray bursts

Short gamma-ray bursts are believed to be produced by both binary neutron star (BNS) and neutron star-black hole (NSBH) mergers. We use current estimates for the BNS and NSBH merger rates to calculate the fraction of observable short gamma-ray bursts produced through each channel. This allows us to constrain merger rates of BNS to $\mathcal{R}_{\rm{BNS}}=384^{+431}_{-213}{\rm{Gpc}^{-3} \rm{yr}^{-1}}$ ($90\%$ credible interval), a $16\%$ decrease in the rate uncertainties from the second LIGO--Virgo Gravitational-Wave Transient Catalog, GWTC-2. Assuming a top-hat emission profile with a large Lorentz factor, we constrain the average opening angle of gamma-ray burst jets produced in BNS mergers to $\approx 15^\circ$. We also measure the fraction of BNS and NSBH mergers that produce an observable short gamma-ray burst to be $0.02^{+0.02}_{-0.01}$ and $0.01 \pm 0.01$, respectively and find that $\gtrsim 40\%$ of BNS mergers launch jets (90\% confidence). We forecast constraints for future gravitational-wave detections given different modelling assumptions, including the possibility that BNS and NSBH jets are different. With $24$ BNS and $55$ NSBH observations, expected within six months of the LIGO-Virgo-KAGRA network operating at design sensitivity, it will be possible to constrain the fraction of BNS and NSBH mergers that launch jets with $10\%$ precision. Within a year of observations, we can determine whether the jets launched in NSBH mergers have a different structure than those launched in BNS mergers and rule out whether $\gtrsim 80\%$ of binary neutron star mergers launch jets. We discuss the implications of future constraints on understanding the physics of short gamma-ray bursts and binary evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Orbital dynamics and extreme scattering event properties from long-term scintillation observations of PSR J1603-7202

We model long-term variations in the scintillation of binary pulsar PSR J1603$-$7202, observed by the 64 m Parkes radio telescope (Murriyang) between 2004 and 2016. We find that the time variation in the scintillation arc curvature is well-modelled by scattering from an anisotropic thin screen of plasma between the Earth and the pulsar. Using our scintillation model, we measure the inclination angle and longitude of ascending node of the orbit, yielding a significant improvement over the constraints from pulsar timing. From our measurement of the inclination angle, we place a lower bound on the mass of J1603$-$7202's companion of $\gtrsim 0.5\,\text{M}_\odot$ assuming a pulsar mass of $\gtrsim1.2\,\text{M}_\odot$. We find that the scintillation arcs are most pronounced when the electron column density along the line of sight is increased, and that arcs are present during a known extreme scattering event. We measure the distance to the interstellar plasma and its velocity, and we discuss some structures seen in individual scintillation arcs within the context of our model.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

When models fail: an introduction to posterior predictive checks and model misspecification in gravitational-wave astronomy

Bayesian inference is a powerful tool in gravitational-wave astronomy. It enables us to deduce the properties of merging compact-object binaries and to determine how these mergers are distributed as a population according to mass, spin, and redshift. As key results are increasingly derived using Bayesian inference, there is increasing scrutiny on Bayesian methods. In this review, we discuss the phenomenon of \textit{model misspecification}, in which results obtained with Bayesian inference are misleading because of deficiencies in the assumed model(s). Such deficiencies can impede our inferences of the true parameters describing physical systems. They can also reduce our ability to distinguish the "best fitting" model: it can be misleading to say that Model~A is preferred over Model~B if both models are manifestly poor descriptions of reality. Broadly speaking, there are two ways in which models fail: models that fail to adequately describe the data (either the signal or the noise) have misspecified likelihoods. Population models -- designed, for example, to describe the distribution of black hole masses -- may fail to adequately describe the true population due to a misspecified prior. We recommend tests and checks that are useful for spotting misspecified models using examples inspired by gravitational-wave astronomy. We include companion python notebooks to illustrate essential concepts.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Black-hole spectroscopy, the no-hair theorem and GW150914: Kerr vs. Occam

The "no-hair" theorem states that astrophysical black holes are fully characterised by just two numbers: their mass and spin. The gravitational-wave emission from a perturbed black-hole consists of a superposition of damped sinusoids, known as \textit{quasi-normal modes}. Quasi-normal modes are specified by three integers $(\ell,m,n)$: the $(\ell, m)$ integers describe the angular properties and $(n)$ specifies the (over)tone. If the no-hair theorem holds, the frequencies and damping times of quasi-normal modes are determined uniquely by the mass and spin of the black hole, while phases and amplitudes depend on the particular perturbation. Current tests of the no-hair theorem, attempt to identify these modes in a semi-agnostic way, without imposing priors on the source of the perturbation. This is usually known as \textit{black-hole spectroscopy}. Applying this framework to GW150914, the measurement of the first overtone led to the confirmation of the theorem to $20\%$ level. We show, however, that such semi-agnostic tests cannot provide strong evidence in favour of the no-hair theorem, even for extremely loud signals, given the increasing number of overtones (and free parameters) needed to fit the data. This can be solved by imposing prior assumptions on the origin of the perturbed black hole that can further constrain the explored parameters: in particular, our knowledge that the ringdown is sourced by a binary black hole merger. Applying this strategy to GW150914 we find a natural log Bayes factor of $\sim 6.5$ in favour of the Kerr nature of its remnant, indicating that the hairy object hypothesis is disfavoured with $<1:600$ with respect to the Kerr black-hole one.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Probing Extremal Gravitational-Wave Events with Coarse-Grained Likelihoods

As catalogs of gravitational-wave transients grow, new records are set for the most extreme systems observed to date. The most massive observed black holes probe the physics of pair instability supernovae while providing clues about the environments in which binary black hole systems are assembled. The least massive black holes, meanwhile, allow us to investigate the purported neutron star-black hole mass gap, and binaries with unusually asymmetric mass ratios or large spins inform our understanding of binary and stellar evolution. Existing outlier tests generally implement leave-one-out analyses, but these do not account for the fact that the event being left out was by definition an extreme member of the population. This results in a bias in the evaluation of outliers. We correct for this bias by introducing a coarse-graining framework to investigate whether these extremal events are true outliers or whether they are consistent with the rest of the observed population. Our method enables us to study extremal events while testing for population model misspecification. We show that this ameliorates biases present in the leave-one-out analyses commonly used within the gravitational-wave community. Applying our method to results from the second LIGO--Virgo transient catalog, we find qualitative agreement with the conclusions of Abbott et al, ApJL 913 L7 (2021). GW190814 is an outlier because of its small secondary mass. We find that neither GW190412 nor GW190521 are outliers.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Temperature dependent appearance of exotic matter makes nascent neutron stars spin faster

Neutron stars offer the opportunity to study the behaviour of matter at densities and temperatures inaccessible to terrestrial experiments. Gravitational-wave observations of binary neutron star coalescences can constrain the neutron-star equation of state before and after merger. After the neutron star binary merges, hyperons can form in the remnant, changing the behaviour of the neutron-star equation of state. In this study, we use finite-entropy equations of state to show that a post-merger remnant can spin up due to cooling. The magnitude of the spin-up depends on the neutron-star equation of state. If hyperons are present, the post-merger spin-up changes the peak gravitational-wave frequency by $\sim 540$ Hz, when the entropy per baryon drops from $s=2$ $k_B$ to $s=0$ $k_B$. If hyperons are not present, the post-merger spin-up changes by $\sim 360$ Hz, providing a gravitational-wave signature for exotic matter. We expect the same qualitative behaviour whenever temperature dependent phase transitions are triggered.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The imprint of superradiance on hierarchical black hole mergers

Ultralight bosons are a proposed solution to outstanding problems in cosmology and particle physics: they provide a dark-matter candidate while potentially explaining the strong charge-parity problem. If they exist, ultralight bosons can interact with black holes through the superradiant instability. In this work we explore the consequences of this instability on the evolution of hierarchical black holes within dense stellar clusters. By reducing the spin of individual black holes, superradiance reduce the recoil velocity of merging binary black holes, which, in turn, increases the retention fraction of hierarchical merger remnants. We show that the existence of ultralight bosons with mass $ 2\times10^{-14}\lesssim μ/\textrm{eV} \lesssim2\times10^{-13}$ would lead to an increased rate of hierarchical black hole mergers in nuclear star clusters. An ultralight boson in this energy range would result in up to $\approx60\%$ more present-day nuclear star clusters supporting hierarchical growth. The presence of an ultralight boson can also double the rate of intermediate mass black hole mergers to $\approx0.08$\,Gpc$^{-3}$\,yr$^{-1}$ in the local Universe. These results imply that a select range of ultralight boson mass can have far-reaching consequences for the population of black holes in dense stellar environments. Future studies into black hole cluster populations and the spin distribution of hierarchically formed black holes will test this scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Constraining short gamma-ray burst jet properties with gravitational waves and gamma rays

Gamma-ray burst (GRB) prompt emission is highly beamed, and understanding the jet geometry and beaming configuration can provide information on the poorly understood central engine and circum-burst environment. Prior to the advent of gravitational-wave astronomy, astronomers relied on observations of jet breaks in the multi-wavelength afterglow to determine the GRB opening angle, since the observer's viewing angle relative to the system cannot be determined from the electromagnetic data alone. Gravitational-wave observations, however, provide an independent measurement of the viewing angle. We describe a Bayesian method for determining the geometry of short GRBs using coincident electromagnetic and gravitational-wave observations. We demonstrate how an ensemble of multi-messenger detections can be used to measure the distributions of the jet energy, opening angle, Lorentz factor, and angular profile of short GRBs; we find that for a population of 100 such observations, we can constrain the mean of the opening angle distribution to within 10 degrees regardless of the angular emission profile. Conversely, the constraint on the energy distribution depends on the shape of the profile, which can be distinguished.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Constraining temperature distribution inside LIGO test masses from frequencies of their vibrational modes

Thermal distortion of test masses, as well as thermal drift of their vibrational mode frequencies, present a major challenge for operation of the Advanced LIGO and Advanced VIRGO interferometers, reducing optical efficiency, which limits sensitivity and potentially causing instabilities which reduce duty-cycle. In this paper, we demonstrate that test-mass vibrational mode frequency data can be used to overcome some of these difficulties. First, we derive a general expression for the change in a mode frequency as a function of temperature distribution inside the test mass. Then we show how the mode frequency dependence on temperature distribution can be used to identify the wavefunction of observed vibrational modes. We then show how monitoring the frequencies of multiple vibrational modes allows the temperature distribution inside the test mass to be strongly constrained. Finally, we demonstrate using simulations, the potential to improve the thermal model of the test mass, providing independent and improved estimates of important parameters such as the coating absorption coefficient and the location of point absorbers.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Constraining the Gravitational-Wave Afterglow From a Binary Neutron Star Coalescence

Binary neutron star mergers are rich laboratories for physics, accessible with ground-based interferometric gravitational-wave detectors such as the Advanced LIGO and Advanced Virgo. If a neutron star remnant survives the merger, it can emit gravitational waves that might be detectable with the current or next generation detectors. The physics of the long-lived post-merger phase is not well understood and makes modelling difficult. In particular the phase of the gravitational-wave signal is not well modelled. In this paper, we explore methods for using long duration post-merger gravitational-wave signals to constrain the parameters and the properties of the remnant. We develop a phase-agnostic likelihood model that uses only the spectral content for parameter estimation and demonstrate the calculation of a Bayesian upper limit in the absence of a signal. With the millisecond magnetar model, we show that for an event like GW170817, the ellipticity of a long-lived remnant can be constrained to less than about 0.5 in the parameter space used.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Gravitational-wave astronomy with a physical calibration model

We carry out astrophysical inference for compact binary merger events in LIGO-Virgo's first gravitational-wave transient catalog (GWTC-1) using a physically motivated calibration model. We demonstrate that importance sampling can be used to reduce the cost of what would otherwise be a computationally challenging analysis. We show that including the physical estimate for the calibration error distribution has negligible impact on the inference of parameters for the events in GWTC-1. Studying a simulated signal with matched filter signal-to-noise ratio $\text{SNR}=200$, we project that a calibration error estimate typical of GWTC-1 is likely to be negligible for the current generation of gravitational-wave detectors. We argue that other sources of systematic error---from waveforms, prior distributions, and noise modelling---are likely to be more important. Finally, using the events in GWTC-1 as standard sirens, we infer an astrophysically-informed improvement on the estimate of the calibration error in the LIGO interferometers.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Gravitational-wave astronomy with an uncertain noise power spectral density

In order to extract information about the properties of compact binaries, we must estimate the noise power spectral density of gravitational-wave data, which depends on the properties of the gravitational-wave detector. In practice, it is not possible to know this perfectly, only to estimate it from the data. Multiple estimation methods are commonly used and each has a corresponding statistical uncertainty. However, this uncertainty is widely ignored when measuring the physical parameters describing compact binary coalescences, and the appropriate likelihoods which account for the uncertainty are not well known. In order to perform increasingly precise astrophysical inference and model selection, it will be essential to account for this uncertainty. In this work, we derive the correct likelihood for one of the most widely used estimation methods in gravitational-wave transient analysis, the median average. We demonstrate that simulated Gaussian noise follows the predicted distributions. We then examine real gravitational-wave data at and around the time of GW151012, a relatively low-significance binary black hole merger event. We show that the data are well described by stationary-Gaussian noise and explore the impact of different noise power spectral density estimation methods on the astrophysical inferences we draw about GW151012.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Inferring the population properties of binary black holes from unresolved gravitational waves

The vast majority of compact binary mergers in the Universe produce gravitational waves that are too weak to yield unambiguous detections; they are unresolved. We present a method to infer the population properties of compact binaries -- such as their merger rates, mass spectrum, and spin distribution -- using both resolved and unresolved gravitational waves. By eliminating entirely the distinction between resolved and unresolved signals, we eliminate bias from selection effects. To demonstrate this method, we carry out a Monte Carlo study using an astrophysically motivated population of binary black holes. We show that some population properties of compact binaries are well constrained by unresolved signals after about one week of observation with Advanced LIGO at design sensitivity.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Is there a spectral turnover in the spin noise of millisecond pulsars?

Pulsar timing arrays provide a unique means to detect nanohertz gravitational waves through long-term measurements of pulse arrival times from an ensemble of millisecond pulsars. After years of observations, some timing array pulsars have been shown to be dominated by low-frequency red noise, including spin noise that might be associated with pulsar rotational irregularities. The power spectral density of pulsar timing red noise is usually modeled with a power law or a power law with a turnover frequency below which the noise power spectrum plateaus. If there is a turnover in the spin noise of millisecond pulsars, residing within the observation band of current and/or future pulsar timing measurements, it may be easier than projected to resolve the gravitational-wave background from supermassive binary black holes. Additionally, the spectral turnover can provide valuable insights on neutron star physics. In the recent study by Melatos and Link, the authors provided a derivation of the model for power spectral density of spin noise from superfluid turbulence in the core of a neutron star, from first principles. The model features a spectral turnover, which depends on the dynamical response time of the superfluid and the steady-state angular velocity lag between the crust and the core of the star. In this work, we search for a spectral turnover in spin noise using the first data release of the International Pulsar Timing Array. Through Bayesian model selection, we find no evidence of a spectral turnover. Our analysis also shows that data from pulsars J1939+2134, J1024$-$0719 and J1713+0747 prefers the power-law model to the superfluid turbulence model.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Measuring the neutron star equation of state with gravitational waves: the first forty binary neutron star mergers

Gravitational waves from binary neutron star coalescences contain rich information about matter at supranuclear densities encoded by the neutron star equation of state. We can measure the equation of state by analyzing the tidal interactions between neutron stars, which is quantified by the tidal deformability. Multiple merger events are required to probe the equation of state over a range of neutron star masses. The more events included in the analysis, the stronger the constraints on the equation of state. In this paper, we build on previous work to explore the constraints that LIGO and Virgo are likely to place on the neutron star equation of state by combining the first forty binary neutron star detections, a milestone we project to be reached during the first year of accumulated design-sensitivity data. We carry out Bayesian inference on a realistic mock dataset of binaries to obtain posterior distributions for neutron star tidal parameters. In order to combine posterior samples from multiple observations, we employ a random forest regressor, which allows us to efficiently interpolate the likelihood distribution. Assuming a merger rate of 1540 Gpc$^{-3}$ yr$^{-1}$ and a LIGO-Virgo detector network operating for one year at the sensitivity of the third-observation run, plus an additional eight months of design sensitivity, we find that the radius of a 1.4 $M_\odot$ neutron star can be constrained to $\sim 10$% at 90% confidence. At the same time, the pressure at twice the nuclear saturation density can be constrained to $\sim 45$ % at 90% confidence. Finally, we add an appendix following publication of the paper in the journal, showing the posterior distribution of the maximum neutron star mass allowed by the equation of state. We find that the maximum mass can be constrained to $\sim 0.3$ $M_\odot$ at 90% confidence.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Memory Effect or Cosmic String? Classifying Gravitational-Wave Bursts with Bayesian Inference

A variety of gravitational-wave transient sources can be modeled in the Fourier domain using a power law. This simple power-law model provides a reasonable approximation for gravitational-wave bursts from cosmic string cusps, cosmic string kinks, and the memory effect. Each of these sources is described using a different spectral index. In this work, we simulate interferometer strain data with injections of power-law and memory bursts to demonstrate parameter estimation, signal detection, and model selection. We show how Bayesian inference can be used to measure the power-law spectral index, thereby distinguishing between different astrophysical scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the origin of GW190425

The LIGO/Virgo collaborations recently announced the detection of a likely binary neutron star merger, GW190425. The total mass of GW190425 is significantly larger than the masses of Galactic double neutron stars known through radio astronomy. This suggests that GW190425 formed differently from Galactic double neutron stars. We hypothesize that GW190425 formed via unstable "case BB" mass transfer. According to this hypothesis, the progenitor of GW190425 was a binary consisting of a neutron star and a ${\sim} 4-5 {M_\odot}$ helium star, which underwent a common-envelope process. Following the supernovae of the helium star core, a tight, eccentric, double neutron star was formed, which merged in ${\lesssim}$ 10 Myr. The helium star progenitor may explain the unusually large mass of GW190425, while the short time to merger may explain why we do not see similar systems in radio. In order to test this hypothesis, we measure the eccentricity of GW190425 using publicly available LIGO/Virgo data. We constrain the eccentricity at 10 Hz to be $e \leq 0.007$ with $90\%$ confidence. This result provides no evidence for or against the unstable mass transfer scenario because the binary is likely to have circularized to $e\lesssim10^{-4}$ by the time it entered the LIGO/Virgo band. Future detectors operating in lower frequency bands will enable us to discern the formation channel of mergers similar to GW190425 using eccentricity measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Standard-siren cosmology using gravitational waves from binary black holes

Gravitational-wave astronomy provides a unique new way to study the expansion history of the Universe. In this work, we investigate the impact future gravitational-wave observatories will have on cosmology. Third-generation observatories like the Einstein Telescope and Cosmic Explorer will be sensitive to essentially all of the binary black hole coalescence events in the Universe. Recent work by \cite{farr2019future} points out that features in the stellar-mass black hole population break the mass-redshift degeneracy, facilitating precise determination of the Hubble parameter without electromagnetic counterparts or host galaxy catalogues. Using a hierarchical Bayesian inference model, we show that with one year of observations by the Einstein Telescope, the Hubble constant will be measured to $\lesssim 1\%$. We also show that this method can be used to perform Bayesian model selection between cosmological models. As an illustrative example, we find that a decisive statement can be made comparing the $Λ$CDM and RHCT cosmological models using two weeks of data from the Einstein Telescope.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The astrophysical odds of GW151216

The gravitational-wave candidate GW151216 is a proposed binary black hole event from the first observing run of the Advanced LIGO detectors. Not identified as a bona fide signal by the LIGO--Virgo collaboration, there is disagreement as to its authenticity, which is quantified by $p_\text{astro}$, the probability that the event is astrophysical in origin. Previous estimates of $p_\text{astro}$ from different groups range from 0.18 to 0.71, making it unclear whether this event should be included in population analyses, which typically require $p_\text{astro}>0.5$. Whether GW151216 is an astrophysical signal or not has implications for the population properties of stellar-mass black holes and hence the evolution of massive stars. Using the astrophysical odds, a Bayesian method which uses the signal coherence between detectors and a parameterised model of non-astrophysical detector noise, we find that $p_\text{astro}=0.03$, suggesting that GW151216 is unlikely to be a genuine signal. We also analyse GW150914 (the first gravitational-wave detection) and GW151012 (initially considered to be an ambiguous detection) and find $p_\text{astro}$ values of 1 and 0.997 respectively. We argue that the astrophysical odds presented here improve upon traditional methods for distinguishing signals from noise.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Toward the unambiguous identification of supermassive binary black holes through Bayesian inference

Supermassive binary black holes at sub-parsec orbital separations have yet to be discovered, with the possible exception of blazar OJ~287. In parallel to the global hunt for nanohertz gravitational waves from supermassive binaries using pulsar timing arrays, there has been a growing sample of candidates reported from electromagnetic surveys, particularly searches for periodic variations in optical light curves of quasars. However, the periodicity search is prone to false positives from quasar red noise and quasi-periodic oscillations from the accretion disc of a single supermassive black hole---especially when the data span fewer than a few signal cycles. We present a Bayesian method for the detection of quasar (quasi-)periodicity in the presence of red noise. We apply this method to the binary candidate PG1302$-$102, and show that a) there is very strong support (Bayes factor $>10^6$) for quasi-periodicity, and b) the data slightly favour a quasi-periodic oscillation over a sinusoidal signal, which we interpret as modest evidence against the binary black hole hypothesis. We also find that the prevalent damped random walk red-noise model is disfavored with more than 99.9\% credibility. Finally, we outline future work that may enable the unambiguous identification of supermassive binary black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Searching for Eccentricity: Signatures of Dynamical Formation in the First Gravitational-Wave Transient Catalogue of LIGO and Virgo

Binary black holes are thought to form primarily via two channels: isolated evolution and dynamical formation. The component masses, spins, and eccentricity of a binary black hole system provide clues to its formation history. We focus on eccentricity, which can be a signature of dynamical formation. Employing the spin-aligned eccentric waveform model SEOBNRE, we perform Bayesian inference to measure the eccentricity of binary black hole merger events in the first Gravitational-Wave Transient Catalogue of LIGO and Virgo. We find that all of these events are consistent with zero eccentricity. We set upper limits on eccentricity ranging from 0.02 to 0.05 with 90% confidence at a reference frequency of 10 Hz. These upper limits do not significantly constrain the fraction of LIGO-Virgo events formed dynamically in globular clusters, because only ~ 5% are expected to merge with measurable eccentricity. However, with the Gravitational-Wave Transient Catalogue set to expand dramatically over the coming months, it may soon be possible to significantly constrain the fraction of mergers taking place in globular clusters using eccentricity measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Thanks for the memory: measuring gravitational-wave memory in the first LIGO/Virgo gravitational-wave transient catalog

Gravitational-wave memory, a strong-field effect of general relativity, manifests itself as a permanent displacement in spacetime. We develop a Bayesian framework to detect gravitational-wave memory with the Advanced LIGO/Virgo detector network. We apply this algorithm on the ten binary black hole mergers in LIGO/Virgo's first transient gravitational-wave catalog. We find no evidence of memory, which is consistent with expectations. In order to estimate when memory will be detected, we use the best current population estimates to construct a realistic sample of binary black hole observations for LIGO/Virgo at design sensitivity. We show that an ensemble of $\mathcal{O}(2000)$ binary black hole observations can be used to find definitive evidence for gravitational-wave memory. We conclude that memory is likely to be detected in the early A+/Virgo+ era.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Ultra-relativistic astrophysics using multi-messenger observations of double neutron stars with LISA and the SKA

Recent work highlights that tens of Galactic double neutron stars are likely to be detectable in the millihertz band of the space-based gravitational-wave observatory, LISA. Kyutoku and Nishino point out that some of these binaries might be detectable as radio pulsars using the Square Kilometer Array (SKA). We point out that the joint LISA+SKA detection of a $f_\text{gw}\gtrsim$1 mHz binary, corresponding to a binary period of $\lesssim$400 s, would enable precision measurements of ultra-relativistic phenomena. We show that, given plausible assumptions, multi-messenger observations of ultra-relativistic binaries can be used to constrain the neutron star equation of state with remarkable fidelity. It may be possible to measure the mass-radius relation with a precision of $\approx$0.2% after 10 yr of observations with the SKA. Such a measurement would be roughly an order of magnitude more precise than possible with other proposed observations. We summarize some of the other remarkable science made possible with multi-messenger observations of millihertz binaries, and discuss the prospects for the detection of such objects.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Detectability of Gravitational Waves from High-Redshift Binaries

Recent non-detection of gravitational-wave backgrounds from pulsar timing arrays casts further uncertainty on the evolution of supermassive black hole binaries. We study the capabilities of current gravitational-wave observatories to detect individual binaries and demonstrate that, contrary to conventional wisdom, some are in principle detectable throughout the Universe. In particular, a binary with rest-frame mass $\gtrsim10^{10}\,M_\odot$ can be detected by current timing arrays at arbitrarily high redshifts. The same claim will apply for less massive binaries with more sensitive future arrays. As a consequence, future searches for nanohertz gravitational waves could be expanded to target evolving high-redshift binaries. We calculate the maximum distance at which binaries can be observed with pulsar timing arrays and other detectors, properly accounting for redshift and using realistic binary waveforms.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Detecting gravitational-wave memory with LIGO: implications of GW150914

It may soon be possible for Advanced LIGO to detect hundreds of binary black hole mergers per year. We show how the accumulation of many such measurements will allow for the detection of gravitational-wave memory: a permanent displacement of spacetime that comes from strong-field, general relativistic effects. We estimate that Advanced LIGO operating at design sensitivity may be able to make a signal-to-noise ratio 3(5) detection of memory with ~35 (90) events with masses and distance similar to GW150914. We highlight the importance of incorporating higher-order gravitational-wave modes for parameter estimation of binary black hole mergers, and describe how our methods can also be used to detect higher-order modes themselves before Advanced LIGO reaches design sensitivity.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Gravitational-wave cosmology across 29 decades in frequency

Quantum fluctuations of the gravitational field in the early Universe, amplified by inflation, produce a primordial gravitational-wave background across a broad frequency band. We derive constraints on the spectrum of this gravitational radiation, and hence on theories of the early Universe, by combining experiments that cover 29 orders of magnitude in frequency. These include Planck observations of cosmic microwave background temperature and polarization power spectra and lensing, together with baryon acoustic oscillations and big bang nucleosynthesis measurements, as well as new pulsar timing array and ground-based interferometer limits. While individual experiments constrain the gravitational-wave energy density in specific frequency bands, the combination of experiments allows us to constrain cosmological parameters, including the inflationary spectral index, $n_t$, and the tensor-to-scalar ratio, $r$. Results from individual experiments include the most stringent nanohertz limit of the primordial background to date from the Parkes Pulsar Timing Array, $Ω_{\rm gw}(f)<2.3\times10^{-10}$. Observations of the cosmic microwave background alone limit the gravitational-wave spectral index at 95\% confidence to $n_t\lesssim5$ for a tensor-to-scalar ratio of $r = 0.11$. However, the combination of all the above experiments limits $n_t<0.36$. Future Advanced LIGO observations are expected to further constrain $n_t<0.34$ by 2020. When cosmic microwave background experiments detect a non-zero $r$, our results will imply even more stringent constraints on $n_t$ and hence theories of the early Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Limits of Astrophysics with Gravitational-Wave Backgrounds

The recent Advanced LIGO detection of gravitational waves from the binary black hole GW150914 suggests there exists a large population of merging binary black holes in the Universe. Although most are too distant to be individually resolved by advanced detectors, the superposition of gravitational waves from many unresolvable binaries is expected to create an astrophysical stochastic background. Recent results from the LIGO and Virgo collaborations show that this astrophysical background is within reach of Advanced LIGO. In principle, the binary black hole background encodes interesting astrophysical properties, such as the mass distribution and redshift distribution of distant binaries. However, we show that this information will be difficult to extract with the current configuration of advanced detectors (and using current data analysis tools). Additionally, the binary black hole background also constitutes a foreground that limits the ability of advanced detectors to observe other interesting stochastic background signals, for example from cosmic strings or phase transitions in the early Universe. We quantify this effect.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Subtraction of correlated noise in global networks of gravitational-wave interferometers

The recent discovery of merging black holes suggests that a stochastic gravitational-wave background is within reach of the advanced detector network operating at design sensitivity. However, correlated magnetic noise from Schumann resonances threatens to contaminate observation of a stochastic background. In this paper, we report on the first effort to eliminate intercontinental correlated noise from Schumann resonances using Wiener filtering. Using magnetometers as proxies for gravitational-wave detectors, we demonstrate as much as a factor of two reduction in the coherence between magnetometers on different continents. While much work remains to be done, our results constitute a proof-of-principle and motivate follow-up studies with a dedicated array of magnetometers.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Detecting compact binary coalescences with seedless clustering

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Measuring neutron-star ellipticity with measurements of the stochastic gravitational-wave background

Galactic neutron stars are a promising source of gravitational waves in the analysis band of detectors such as LIGO and Virgo. Previous searches for gravitational waves from neutron stars have focused on the detection of individual neutron stars, which are either nearby or highly non-spherical. Here we consider the stochastic gravitational-wave signal arising from the ensemble of Galactic neutron stars. Using a population synthesis model, we estimate the single-sigma sensitivity of current and planned gravitational-wave observatories to average neutron star ellipticity $ε$ as a function of the number of in-band Galactic neutron stars $N_\text{tot}$. For the plausible case of $N_\text{tot}\approx 53000$, and assuming one year of observation time with colocated initial LIGO detectors, we find it to be $σ_ε=2.1\times10^{-7}$, which is comparable to current bounds on some nearby neutron stars. (The current best $95\%$ upper limits are $ε\lesssim7\times10^{-8}.$) It is unclear if Advanced LIGO can significantly improve on this sensitivity using spatially separated detectors. For the proposed Einstein Telescope, we estimate that $σε=5.6\times10^{-10}$. Finally, we show that stochastic measurements can be combined with measurements of individual neutron stars in order to estimate the number of in-band Galactic neutron stars. In this way, measurements of stochastic gravitational waves provide a complementary tool for studying Galactic neutron stars.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Statistical properties of astrophysical gravitational-wave backgrounds

We investigate how a stochastic gravitational wave background, produced from a discrete set of astrophysical sources, differs from an idealised model consisting of an isotropic, unpolarised, and Gaussian background. We focus, in particular, on the different signatures produced from these two cases, as observed in a cross-correlation search. We show that averaged over many realisations of an astrophysical background, the cross-correlation measurement of an astrophysical background is identical to that of an idealised background. However, any one realisation of an astrophysical background can produce a different signature. Using a model consisting of an ensemble of binary neutron star coalescences, we quantify the typical difference between the signal from individual realisations of the astrophysical background and the idealised case. For advanced detectors, we find that, using a cross-correlation analysis, astrophysical backgrounds from many discrete sources are probably indistinguishable from an idealised background.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The detectability of cosmological gravitational-wave backgrounds: a rule of thumb

The recent claim by BICEP2 of evidence for primordial gravitational waves from inflation has focused interest on the potential for early-Universe cosmology using observations of gravitational waves. In addition to cosmic microwave background detectors, efforts are underway to carry out gravitational-wave astronomy over a wide range of frequencies including pulsar timing arrays (nHz), space-based detectors (mHz), and terrestrial detectors ($\sim$10-2000 Hz). This multiband effort will probe a wide range of times in the early Universe (each corresponding to a different energy scale), during which gravitational-wave backgrounds may have been produced through processes such as phase transitions or preheating. In this letter, we derive a rule of thumb (not quite so strong as an upper limit) governing the maximum energy density of cosmological backgrounds. For most cosmological scenarios, we expect the energy density spectrum to peak at values of $Ω_\text{gw}(f)\lesssim10^{-12\pm2}$. We discuss the applicability of this rule of thumb and the implications for gravitational-wave astronomy.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Measuring the non-Gaussian stochastic gravitational-wave background: a method for realistic interferometer data

A stochastic gravitational-wave background (SGWB) can arise from the superposition of many independent events. If the rate of events per unit time is sufficiently high, the resulting background is Gaussian, which is to say that it is characterized only by a gravitational-wave strain power spectrum. Alternatively, if the event rate is low, we expect a non-Gaussian background, characterized by intermittent sub-threshold bursts. Many experimentally accessible models of the SGWB, such as the SGWB arising from compact binary coalescences, are expected to be of this non-Gaussian variety. Primordial backgrounds from the early universe, on the other hand, are more likely to be Gaussian. Measuring the Gaussianity of the SGWB can therefore provide additional information about its origin. In this paper we introduce a novel maximum likelihood estimator that can be used to estimate the non-Gaussian component of an SGWB signature measured in a network of interferometers. This method can be robustly applied to spatially separated interferometers with colored, non-Gaussian noise. Furthermore, it can be cast as a generalization of the widely used stochastic radiometer algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Sensitivity curves for searches for gravitational-wave backgrounds

We propose a graphical representation of detector sensitivity curves for stochastic gravitational-wave backgrounds that takes into account the increase in sensitivity that comes from integrating over frequency in addition to integrating over time. This method is valid for backgrounds that have a power-law spectrum in the analysis band. We call these graphs "power-law integrated curves." For simplicity, we consider cross-correlation searches for unpolarized and isotropic stochastic backgrounds using two or more detectors. We apply our method to construct power-law integrated sensitivity curves for second-generation ground-based detectors such as Advanced LIGO, space-based detectors such as LISA and the Big Bang Observer, and timing residuals from a pulsar timing array. The code used to produce these plots is available at https://dcc.ligo.org/LIGO-P1300115/public for researchers interested in constructing similar sensitivity curves.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gravitational Waves from Fallback Accretion onto Neutron Stars

Massive stars generally end their lives as neutron stars (NSs) or black holes (BHs), with NS formation typically occurring at the low mass end and collapse to a BH more likely at the high mass end. In an intermediate regime, with a mass range that depends on the uncertain details of rotation and mass loss during the star's life, a NS is initially formed which then experiences fallback accretion and collapse to a BH. The electromagnetic consequence of such an event is not clear. Depending on the progenitor's structure, possibilities range from a long gamma-ray burst to a Type II supernova (that may or may not be jet-powered) to a collapse with a weak electromagnetic signature. Gravitational waves (GWs) provide the exciting opportunity to peer through the envelope of a dying massive star and directly probe what is occurring inside. We explore whether fallback onto young NSs can be detected by ground-based interferometers. When the incoming material has sufficient angular momentum to form a disk, the accretion spins up the NS sufficiently to produce non-axisymmetric instabilities and gravitational radiation at frequencies of ~700-2400 Hz for ~30-3000 s until collapse to a BH occurs. Using a realistic excess cross-power search algorithm, we show that such events are detectable by Advanced LIGO out to ~17 Mpc. From the rate of nearby core-collapse supernovae, we estimate that there will be ~1-2 events each year that are worth checking for fallback GWs. The observation of these unique GW signatures coincident with electromagnetic detections would identify the transient events that are associated with this channel of BH formation, while providing information about the protoneutron star progenitor.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Identification of noise artifacts in searches for long-duration gravitational-wave transients

We present an algorithm for the identification of transient noise artifacts (glitches) in cross-correlation searches for long O(10s) gravitational-wave transients. The algorithm utilizes the auto-power in each detector as a discriminator between well-behaved Gaussian noise (possibly including a gravitational-wave signal) and glitches. We test the algorithm with both Monte Carlo noise and time-shifted data from the LIGO S5 science run and find that it is effective at removing a significant fraction of glitches while keeping the vast majority (99.6%) of the data. Using an accretion disk instability signal model, we estimate that the algorithm is accidentally triggered at a rate of less than 10^-5% by realistic signals, and less than 3% even for exceptionally loud signals. We conclude that the algorithm is a safe and effective method for cleaning the cross-correlation data used in searches for long gravitational-wave transients.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Parameter Estimation in Searches for the Stochastic Gravitational-Wave Background

The stochastic gravitational-wave background (SGWB) is expected to arise from the superposition of many independent and unresolved gravitational-wave signals of either cosmological or astrophysical origin. The spectral content of the SGWB carries signatures of the physics that generated it. We present a Bayesian framework for estimating the parameters associated with different SGWB models using data from gravitational-wave detectors. We apply this technique to recent results from LIGO to produce the first simultaneous 95% confidence level limits on multiple parameters in generic power-law SGWB models and in SGWB models of compact binary coalescences. We also estimate the sensitivity of the upcoming second-generation detectors such as Advanced LIGO/Virgo to these models and demonstrate how SGWB measurements can be combined and compared with observations of individual compact binary coalescences in order to build confidence in the origin of an observed SGWB signal. In doing so, we demonstrate a novel means of differentiating between different sources of the SGWB.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Long gravitational-wave transients and associated detection strategies for a network of terrestrial interferometers

Searches for gravitational waves (GWs) traditionally focus on persistent sources (e.g., pulsars or the stochastic background) or on transients sources (e.g., compact binary inspirals or core-collapse supernovae), which last for timescales of milliseconds to seconds. We explore the possibility of long GW transients with unknown waveforms lasting from many seconds to weeks. We propose a novel analysis technique to bridge the gap between short O(s) burst analyses and persistent stochastic analyses. Our technique utilizes frequency-time maps of GW strain cross-power between two spatially separated terrestrial GW detectors. The application of our cross-power statistic to searches for GW transients is framed as a pattern recognition problem, and we discuss several pattern-recognition techniques. We demonstrate these techniques by recovering simulated GW signals in simulated detector noise. We also recover environmental noise artifacts, thereby demonstrating a novel technique for the identification of such artifacts in GW interferometers. We compare the efficiency of this framework to other techniques such as matched filtering.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Probing the anisotropies of a stochastic gravitational-wave background using a network of ground-based laser interferometers

We present a maximum-likelihood analysis for estimating the angular distribution of power in an anisotropic stochastic gravitational-wave background using ground-based laser interferometers. The standard isotropic and gravitational-wave radiometer searches (optimal for point sources) are recovered as special limiting cases. The angular distribution can be decomposed with respect to any set of basis functions on the sky, and the single-baseline, cross-correlation analysis is easily extended to a network of three or more detectors-that is, to multiple baselines. A spherical harmonic decomposition, which provides maximum-likelihood estimates of the multipole moments of the gravitational-wave sky, is described in detail. We also discuss: (i) the covariance matrix of the estimators and its relationship to the detector response of a network of interferometers, (ii) a singular-value decomposition method for regularizing the deconvolution of the detector response from the measured sky map, (iii) the expected increase in sensitivity obtained by including multiple baselines, and (iv) the numerical results of this method when applied to simulated data consisting of both point-like and di#use sources. Comparisions between this general method and the standard isotropic and radiometer searches are given throughout, to make contact with the existing literature on stochastic background searches.

0 citations0 reviewsNo code linkedOpen access