Jolien D. E. Creighton

FINDCHIRP: an algorithm for detection of gravitational waves from inspiraling compact binaries

Matched-filter searches for gravitational waves from coalescing compact binaries by the LIGO nScientific Collaboration use the FINDCHIRP algorithm: an implementation of the optimal filter with ninnovations to account for unknown signal parameters and to improve performance on detector data nthat has nonstationary and non-Gaussian artifacts. We provide details on the FINDCHIRP algorithm as used nin the search for subsolar mass binaries, binary neutron stars, neutron starblack hole binaries, and binary nblack holes.

Gravitational-Wave Stochastic Background from Cosmic Strings

We consider the stochastic background of gravitational waves produced by a network of cosmic strings and assess their accessibility to current and planned gravitational wave detectors, as well as to big bang nucleosynthesis (BBN), cosmic microwave background (CMB), and pulsar timing constraints. We find that current data from interferometric gravitational wave detectors, such as Laser Interferometer Gravitational Wave Observatory (LIGO), are sensitive to areas of parameter space of cosmic string models complementary to those accessible to pulsar, BBN, and CMB bounds. Future more sensitive LIGO runs and interferometers such as Advanced LIGO and Laser Interferometer Space Antenna (LISA) will be able to explore substantial parts of the parameter space.

Measuring the neutron star equation of state with gravitational wave observations

We report the results of a first study that uses numerical simulations to estimate the accuracy with which one can use gravitational wave observations of double neutron-star inspiral to measure parameters of the neutron-star equation of state. The simulations use the evolution and initial-data codes of Shibata and Uryu to compute the last several orbits and the merger of neutron stars, with matter described by a parametrized equation of state. Previous work suggested the use of an effective cutoff frequency to place constraints on the equation of state. We find, however, that greater accuracy is obtained by measuring departures from the point-particle limit of the gravitational waveform produced during the late inspiral. As the stars approach their final plunge and merger, the gravitational wave phase accumulates more rapidly for smaller values of the neutron-star compactness (the ratio of the mass of the neutron-star to its radius). We estimate that realistic equations of state will lead to gravitational waveforms that are distinguishable from point-particle inspirals at an effective distance (the distance to an optimally oriented and located system that would produce an equivalent waveform amplitude) of 100 Mpc or less. As Lattimer and Prakash observed, neutron-star radius is closely tied tomorexa0» the pressure at density not far above nuclear. Our results suggest that broadband gravitational wave observations at frequencies between 500 and 1000 Hz will constrain this pressure, and we estimate the accuracy with which it can be measured. Related first estimates of radius measurability show that the radius can be determined to an accuracy of {delta}R{approx}1 km at 100 Mpc.«xa0less

Gravitational wave bursts from cosmic (super)strings: Quantitative analysis and constraints

We discuss data analysis techniques that can be used in the search for gravitational wave bursts from cosmic strings. When data from multiple interferometers are available, we describe consistency checks that can be used to greatly reduce the false alarm rates. We construct an expression for the rate of bursts for arbitrary cosmic string loop distributions and apply it to simple known solutions. The cosmology is solved exactly and includes the effects of a late-time acceleration. We find substantially lower burst rates than previous estimates suggest and explain the disagreement. Initial LIGO is unlikely to detect field-theoretic cosmic strings with the usual loop sizes, though it may detect cosmic superstrings as well as cosmic strings and superstrings with nonstandard loop sizes (which may be more realistic). In the absence of a detection, we show how to set upper limits based on the loudest event. Using Initial LIGO sensitivity curves, we show that these upper limits may result in interesting constraints on the parameter space of theories that lead to the production of cosmic strings.

Matter effects on binary neutron star waveforms

Using an extended set of equations of state and a multiple-group multiple-code collaborative effort to generate waveforms, we improve numerical-relativity-based data-analysis estimates of the measurability of matter effects in neutron-star binaries. We vary two parameters of a parameterized piecewise-polytropic equation of state (EOS) to analyze the measurability of EOS properties, via a parameter {Lambda} that characterizes the quadrupole deformability of an isolated neutron star. We find that, to within the accuracy of the simulations, the departure of the waveform from point-particle (or spinless double black-hole binary) inspiral increases monotonically with {Lambda}, and changes in the EOS that did not change {Lambda} are not measurable. We estimate with two methods the minimal and expected measurability of {Lambda} in second- and third- generation gravitational-wave detectors. The first estimate, using numerical waveforms alone, shows two EOS which vary in radius by 1.3km are distinguishable in mergers at 100Mpc. The second estimate relies on the construction of hybrid waveforms by matching to post-Newtonian inspiral, and estimates that the same EOS are distinguishable in mergers at 300Mpc. We calculate systematic errors arising from numerical uncertainties and hybrid construction, and we estimate the frequency at which such effects would interfere with template-based searches.

Optimal strategies for gravitational wave stochastic background searches in pulsar timing data

A low frequency stochastic background of gravitational waves may be detected by pulsar timing experiments in the next 5 to 10 yr. Using methods developed to analyze interferometric gravitational wave data, in this paper we lay out the optimal techniques to detect a background of gravitational waves using a pulsar timing array. We show that for pulsar distances and gravitational wave frequencies typical of pulsar timing experiments, neglecting the effect of the metric perturbation at the pulsar does not result in a significant deviation from optimality. We discuss methods for setting upper limits using the optimal statistic, show how to construct skymaps using the pulsar timing array, and consider several issues associated with realistic analysis of pulsar timing data.

OPTIMAL STRATEGIES FOR CONTINUOUS GRAVITATIONAL WAVE DETECTION IN PULSAR TIMING ARRAYS

Supermassive black hole binaries (SMBHBs) are expected to emit a continuous gravitational wave signal in the pulsar timing array (PTA) frequency band (10–9 to 10–7 Hz). The development of data analysis techniques aimed at efficient detection and characterization of these signals is critical to the gravitational wave detection effort. In this paper, we leverage methods developed for LIGO continuous wave gravitational searches and explore the use of the -statistic for such searches in pulsar timing data. Babak & Sesana have used this approach in the context of PTAs to show that one can resolve multiple SMBHB sources in the sky. Our work improves on several aspects of prior continuous wave search methods developed for PTA data analysis. The algorithm is implemented fully in the time domain, which naturally deals with the irregular sampling typical of PTA data and avoids spectral leakage problems associated with frequency domain methods. We take into account the fitting of the timing model and have generalized our approach to deal with both correlated and uncorrelated colored noise sources. We also develop an incoherent detection statistic that maximizes over all pulsar-dependent contributions to the likelihood. To test the effectiveness and sensitivity of our detection statistics, we perform a number of Monte Carlo simulations. We produce sensitivity curves for PTAs of various configurations and outline an implementation of a fully functional data analysis pipeline. Finally, we present a derivation of the likelihood maximized over the gravitational wave phases at the pulsar locations, which results in a vast reduction of the search parameter space.

Analysis Framework for the Prompt Discovery of Compact Binary Mergers in Gravitational-wave Data

We describe a stream-based analysis pipeline to detect gravitational waves from the merger of binary neutron stars, binary black holes, and neutron-star–black-hole binaries within ∼1 min of the arrival of the merger signal at Earth. Such low-latency detection is crucial for the prompt response by electromagnetic facilities in order to observe any fading electromagnetic counterparts that might be produced by mergers involving at least one neutron star. Even for systems expected not to produce counterparts, low-latency analysis of the data is useful for deciding when not to point telescopes, and as feedback to observatory operations. Analysts using this pipeline were the first to identify GW151226, the second gravitational-wave event ever detected. The pipeline also operates in an offline mode, in which it incorporates more refined information about data quality and employs acausal methods that are inapplicable to the online mode. The pipeline’s offline mode was used in the detection of the first two gravitational-wave events, GW150914 and GW151226, as well as the identification of a third candidate, LVT151012.

Systematic and statistical errors in a Bayesian approach to the estimation of the neutron-star equation of state using advanced gravitational wave detectors

Advanced ground-based gravitational-wave detectors are capable of measuring tidal influences in binary neutron-star systems. In this work, we report on the statistical uncertainties in measuring tidal deformability with a full Bayesian parameter estimation implementation. We show how simultaneous measurements of chirp mass and tidal deformability can be used to constrain the neutron-star equation of state. We also study the effects of waveform modeling bias and individual instances of detector noise on these measurements. We notably find that systematic error between post-Newtonian waveform families can significantly bias the estimation of tidal parameters, thus motivating the continued development of waveform models that are more reliable at high frequencies.

Learning about compact binary merger: The interplay between numerical relativity and gravitational-wave astronomy

Activities in data analysis and numerical simulation of gravitational waves have to date largely proceeded independently. In this work we study how waveforms obtained from numerical simulations could be effectively used within the data analysis effort to search for gravitational waves from black hole binaries. To this end we analyze the cross-correlation between different numerical waveforms weighted by the detectors noise. This allow us to propose measures to quantify the accuracy of numerical waveforms for the purpose of data analysis, study how sensitive the analysis is to errors in the waveforms, and propose a way to efficiently encode the waveforms information for its use as a member of the template bank. We estimate that {approx}100 templates (and {approx}10 simulations with different mass ratios) are needed to detect waves from nonspinning binary black holes with total masses in the range 100M{sub {center_dot}}{<=}M{<=}400M{sub {center_dot}} using initial LIGO. Of course, many more simulation runs will be needed to confirm that the correct physics is captured in the numerical evolutions. From this perspective, we also discuss sources of systematic errors in numerical waveform extraction and provide order of magnitude estimates for the computational cost of simulations that could be used to estimate the costmorexa0» of parameter space surveys. Finally, we discuss what information from near-future numerical simulations of compact binary systems would be most useful for enhancing the detectability of such events with contemporary gravitational-wave detectors and emphasize the role of numerical simulations for the interpretation of eventual gravitational-wave observations.«xa0less