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

Neil J. Cornish

Neil J. Cornish 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.COastro-ph.IMastro-ph.GAastro-phhep-ph

Catalog footprint

What is connected

50works
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

50 published item(s)

preprint2026arXiv

The NANOGrav 15 yr Data Set: Piecewise Power-Law Reconstruction of the Gravitational-Wave Background

The NANOGrav 15-year (NG15) data set provides evidence for a gravitational-wave background (GWB) signal at nanohertz frequencies, which is expected to originate either from a cosmic population of inspiraling supermassive black-hole binaries or new particle physics in the early Universe. A firm identification of the source of the NG15 signal requires an accurate reconstruction of its frequency spectrum. In this paper, we provide such a spectral characterization of the NG15 signal based on a piecewise power-law (PPL) ansatz that strikes a balance between existing alternatives in the literature. Our PPL reconstruction is more flexible than the standard constant-power-law model, which describes the GWB spectrum in terms of only two parameters: an amplitude A and a spectral index gamma. Concurrently, it better approximates physically realistic GWB spectra -- especially those of cosmological origin -- than the free spectral model, since the latter allows for arbitrary variations in the GWB amplitude from one frequency bin to the next. Our PPL reconstruction of the NG15 signal relies on individual PPL models with a fixed number of internal nodes (i.e., constant power law, broken power law, doubly broken power law, etc.) that are ultimately combined in a Bayesian model average. The data products resulting from our analysis provide the basis for fast refits of spectral GWB models.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Comparison of maximum likelihood mapping methods for gravitational-wave backgrounds

Detection of a stochastic background of gravitational waves is likely to occur in the next few years. Beyond searches for the isotropic component of SGWBs, there have been various mapping methods proposed to target anisotropic backgrounds. Some of these methods have been applied to data taken by the Laser Interferometer Gravitational-wave Observatories (LIGO) and Virgo. Specifically, these directional searches have focused on mapping the intensity of the signal on the sky via maximum likelihood solutions. We compare this intensity mapping approach to a previously proposed, but never employed, amplitude-phase mapping method to understand whether this latter approach may be employed in future searches. We build up our understanding of the differences between these two approaches by analysing simple toy models of time-stream data, and run mock-data mapping tests for the two methods. We find that the amplitude-phase method is only applicable to the case of a background which is phase-coherent on large scales or, at the very least, has an intrinsic coherence scale that is larger than that of the detector. Otherwise, the amplitude-phase mapping method leads to a loss of overall information, with respect to both phase and amplitude. Since we do not expect these phase-coherent properties to hold for any of the gravitational-wave background signals we hope to detect in the near future, we conclude that intensity mapping is the preferred method for such backgrounds.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Fast Bayesian analysis of individual binaries in pulsar timing array data

Searching for gravitational waves in pulsar timing array data is computationally intensive. The data is unevenly sampled, and the noise is heteroscedastic, necessitating the use of a time-domain likelihood function with attendant expensive matrix operations. The computational cost is exacerbated when searching for individual supermassive black hole binaries, which have a large parameter space due to the additional pulsar distance, phase offset and noise model parameters needed for each pulsar. We introduce a new formulation of the likelihood function which can be used to make the Bayesian analysis significantly faster. We divide the parameters into projection and shape parameters. We then accelerate the exploration of the projection parameters by more than four orders of magnitude by precomputing the expensive inner products for each set of shape parameters. The projection parameters include nuisance parameters such as the gravitational wave phase offset at each pulsar. In the new scheme, these troublesome nuisance parameters are efficiently marginalized over using multiple-try Markov chain Monte Carlo sampling as part of a Metropolis-within-Gibbs scheme. The acceleration provided by our method will become increasingly important as pulsar timing datasets rapidly grow. Our method also makes sophisticated analyses more tractable, such as searches for multiple binaries, or binaries with non-negligible eccentricities.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Probing neutron stars with the full pre-merger and post-merger gravitational wave signal from binary coalescences

The gravitational wave signal emitted during the coalescence of two neutron stars carries information about the stars' internal structure. During the long inspiral phase the main matter observable is the tidal interaction between the binary components, an effect that can be parametrically modeled with compact-binary solutions to General Relativity. After the binary merger the main observable is frequency modes of the remnant, most commonly giving rise to a short-duration signal accessible only through numerical simulations. The complicated morphology and the decreasing detector sensitivity in the relevant frequencies currently hinder detection of the post-merger signal and motivate separate analyses for the pre-merger and post-merger data. However, planned and ongoing detector improvements could soon put the post-merger signal within reach. In this study we target the whole pre-merger and post-merger signal without an artificial separation at the binary merger. We construct a hybrid analysis that models the inspiral with templates based on analytical calculations and calibrated to numerical relativity and the post-merger signal with a flexible morphology-independent analysis. Applying this analysis to GW170817 we find, as expected, that the post-merger signal remains undetected. We further study simulated signals and find that we can reconstruct the full signal and simultaneously estimate both the pre-merger tidal deformation and the post-merger signal frequency content. Our analysis allows us to study neutron star physics using all the data available and directly test the pre-merger and post-merger signal for consistency thus probing effects such as the onset of the hadron-quark phase transition.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Searching For Gravitational Waves From Cosmological Phase Transitions With The NANOGrav 12.5-year dataset

We search for a first-order phase transition gravitational wave signal in 45 pulsars from the NANOGrav 12.5 year dataset. We find that the data can be modeled in terms of a strong first order phase transition taking place at temperatures below the electroweak scale. However, we do not observe any strong preference for a phase-transition interpretation of the signal over the standard astrophysical interpretation in terms of supermassive black holes mergers; but we expect to gain additional discriminating power with future datasets, improving the signal to noise ratio and extending the sensitivity window to lower frequencies. An interesting open question is how well gravitational wave observatories could separate such signals.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Low Latency Detection of Massive Black Hole Binaries

The next decade is expected to see the launch of one or more space based gravitational wave detectors: the European lead Laser Interferometer Space Antenna (LISA); and one or more Chinese mission concepts, Taiji and TianQin. One of the primary scientific targets for these missions are the mergers of black holes with masses between $10^3 M_\odot$ and $10^8 M_\odot$. These systems may produce detectable electromagnetic signatures in additional to gravitational waves due to the presence of gas in mini-disks around each black hole, and a circumbinary disk surrounding the system. The electromagnetic emission may occur before, during and after the merger. In order to have the best chance of capturing all phases of the emission it is imperative that the gravitational wave signals can be detected in low latency, and used to produce reliable estimates for the sky location and distance to help guide the search for counterparts. The low latency detection also provides a starting point for the ``global fit'' of the myriad signals that are simultaneously present in the data. Here a low latency analysis pipeline is presented that is capable of analyzing months of data in just a few hours using a laptop from the last decade. The problem of performing a global fit is avoided by whitening out the bright foreground produced by nearby galactic binaries. The performance of the pipeline is illustrated using simulated data from the LISA Data Challenge.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Massive Black Hole Binaries and Where to Find Them with Dual Detector Networks

A single space-based gravitational wave detector will push the boundaries of astronomy and fundamental physics. Having a network of two or more detectors would significantly improve source localization. Here we consider how dual networks of space-based detectors would improve parameter estimation of massive black hole binaries. We consider two scenarios: a network comprised of the Laser Interferometer Space Antenna (LISA) and an additional LISA-like heliocentric detector (e.g. Taiji); and a network comprised of LISA with an an additional geocentric detector (e.g. TianQin). We use Markov chain Monte Carlo techniques and Fisher matrix estimates to explore the impact of a two detector network on sky localization and distance determination. The impact on other source parameters is also studied. With the addition of a Taiji or TianQin, we find orders of magnitude improvements in sky localization for the more massive MBHBs, while also seeing improvements for lower mass systems, and for other source parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The BayesWave analysis pipeline in the era of gravitational wave observations

We describe updates and improvements to the BayesWave gravitational wave transient analysis pipeline, and provide examples of how the algorithm is used to analyze data from ground-based gravitational wave detectors. BayesWave models gravitational wave signals in a morphology-independent manner through a sum of frame functions, such as Morlet-Gabor wavelets or chirplets. BayesWave models the instrument noise using a combination of a parametrized Gaussian noise component and non-stationary and non-Gaussian noise transients. Both the signal model and noise model employ trans-dimensional sampling, with the complexity of the model adapting to the requirements of the data. The flexibility of the algorithm makes it suitable for a variety of analyses, including reconstructing generic unmodeled signals; cross checks against modeled analyses for compact binaries; as well as separating coherent signals from incoherent instrumental noise transients (glitches). The BayesWave model has been extended to account for gravitational wave signals with generic polarization content and the simultaneous presence of signals and glitches in the data. We describe updates in the BayesWave prior distributions, sampling proposals, and burn-in stage that provide significantly improved sampling efficiency. We present standard review checks indicating the robustness and convergence of the BayesWave trans-dimensional sampler.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The NANOGrav 12.5-year Data Set: Search For An Isotropic Stochastic Gravitational-Wave Background

We search for an isotropic stochastic gravitational-wave background (GWB) in the $12.5$-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. Our analysis finds strong evidence of a stochastic process, modeled as a power-law, with common amplitude and spectral slope across pulsars. The Bayesian posterior of the amplitude for an $f^{-2/3}$ power-law spectrum, expressed as the characteristic GW strain, has median $1.92 \times 10^{-15}$ and $5\%$--$95\%$ quantiles of $1.37$--$2.67 \times 10^{-15}$ at a reference frequency of $f_\mathrm{yr} = 1 ~\mathrm{yr}^{-1}$. The Bayes factor in favor of the common-spectrum process versus independent red-noise processes in each pulsar exceeds $10,000$. However, we find no statistically significant evidence that this process has quadrupolar spatial correlations, which we would consider necessary to claim a GWB detection consistent with general relativity. We find that the process has neither monopolar nor dipolar correlations, which may arise from, for example, reference clock or solar system ephemeris systematics, respectively. The amplitude posterior has significant support above previously reported upper limits; we explain this in terms of the Bayesian priors assumed for intrinsic pulsar red noise. We examine potential implications for the supermassive black hole binary population under the hypothesis that the signal is indeed astrophysical in nature.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Black Hole Hunting with LISA

The Laser Interferometer Space Antenna (LISA) will be able to detect massive black hole mergers throughout the visible Universe. These observations will provide unique information about black hole formation and growth, and the role black holes play in galaxy evolution. Here we develop several key building blocks for detecting and characterizing black hole binary mergers with LISA, including fast heterodyned likelihood evaluations, and efficient stochastic search techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Joint search for isolated sources and an unresolved confusion background in pulsar timing array data

Supermassive black hole binaries are the most promising source of gravitational-waves in the frequency band accessible to pulsar timing arrays. Most of these binaries will be too distant to detect individually, but together they will form an approximately stochastic background that can be detected by measuring the correlation pattern induced between pairs of pulsars. A small number of nearby and especially massive systems may stand out from this background and be detected individually. Analyses have previously been developed to search for stochastic signals and isolated signals separately. Here we present BayesHopper, an algorithm capable of jointly searching for both signal components simultaneously using trans-dimensional Bayesian inference. Our implementation uses the Reversible Jump Markov Chain Monte Carlo method for sampling the relevant parameter space with changing dimensionality. We have tested BayesHopper on various simulated datasets. We find that it gives results consistent with fixed-dimensional methods when tested on data with a stochastic background or data with a single binary. For the full problem of analyzing a dataset with both a background and multiple black hole binaries, we find two kinds of interactions between the binary and background components. First, the background effectively increases the noise level, thus making individual binary signals less significant. Second, weak binary signals can be absorbed by the background model due to the natural parsimony of Bayesian inference. Because of its flexible model structure, we anticipate that BayesHopper will outperform existing approaches when applied to realistic data sets produced from population synthesis models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Multi-Messenger Gravitational Wave Searches with Pulsar Timing Arrays: Application to 3C66B Using the NANOGrav 11-year Data Set

When galaxies merge, the supermassive black holes in their centers may form binaries and, during the process of merger, emit low-frequency gravitational radiation in the process. In this paper we consider the galaxy 3C66B, which was used as the target of the first multi-messenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C66B to less than $(1.65\pm0.02) \times 10^9~{M_\odot}$ using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits, and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data to `blind' pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Rapid and reliable sky localization of gravitational wave sources

The first detection of gravitational waves by LIGO from the merger of two compact objects has sparked new interest in detecting electromagnetic counterparts to these violent events. For mergers involving neutron stars, it is thought that prompt high-energy emission in gamma rays and x-rays will be followed days to weeks later by an afterglow in visible light, infrared and radio. Rapid sky localization using the data from a network of gravitational wave detectors is essential to maximize the chances of making a joint detection. Here I describe a new technique that is able to produce accurate, fully Bayesian sky maps in seconds or less. The technique can be applied to spin-precessing compact binaries, and can take into account detector calibration and spectral estimation uncertainties.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

BayesWave: Bayesian Inference for Gravitational Wave Bursts and Instrument Glitches

A central challenge in Gravitational Wave Astronomy is identifying weak signals in the presence of non-stationary and non-Gaussian noise. The separation of gravitational wave signals from noise requires good models for both. When accurate signal models are available, such as for binary Neutron star systems, it is possible to make robust detection statements even when the noise is poorly understood. In contrast, searches for "un-modeled" transient signals are strongly impacted by the methods used to characterize the noise. Here we take a Bayesian approach and introduce a multi-component, variable dimension, parameterized noise model that explicitly accounts for non-stationarity and non-Gaussianity in data from interferometric gravitational wave detectors. Instrumental transients (glitches) and burst sources of gravitational waves are modeled using a Morlet-Gabor continuous wavelet frame. The number and placement of the wavelets is determined by a trans-dimensional Reversible Jump Markov Chain Monte Carlo algorithm. The Gaussian component of the noise and sharp line features in the noise spectrum are modeled using the BayesLine algorithm, which operates in concert with the wavelet model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Constraining the Solution to the Last Parsec Problem with Pulsar Timing

The detection of a stochastic gravitational-wave signal from the superposition of many inspiraling supermassive black holes with pulsar timing arrays (PTAs) is likely to occur within the next decade. With this detection will come the opportunity to learn about the processes that drive black-hole-binary systems toward merger through their effects on the gravitational-wave spectrum. We use Bayesian methods to investigate the extent to which effects other than gravitational-wave emission can be distinguished using PTA observations. We show that, even in the absence of a detection, it is possible to place interesting constraints on these dynamical effects for conservative predictions of the population of tightly bound supermassive black-hole binaries. For instance, if we assume a relatively weak signal consistent with a low number of bound binaries and a low black-hole-mass to galaxy-mass correlation, we still find that a non-detection by a simulated array, with a sensitivity that should be reached in practice within a few years, disfavors gravitational-wave-dominated evolution with an odds ratio of $\sim$30:1. Such a finding would suggest either that all existing astrophysical models for the population of tightly bound binaries are overly optimistic, or else that some dynamical effect other than gravitational-wave emission is actually dominating binary evolution even at the relatively high frequencies/small orbital separations probed by PTAs.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Enabling high confidence detections of gravitational-wave bursts

With the advanced LIGO and Virgo detectors taking observations the detection of gravitational waves is expected within the next few years. Extracting astrophysical information from gravitational wave detections is a well-posed problem and thoroughly studied when detailed models for the waveforms are available. However, one motivation for the field of gravitational wave astronomy is the potential for new discoveries. Recognizing and characterizing unanticipated signals requires data analysis techniques which do not depend on theoretical predictions for the gravitational waveform. Past searches for short-duration un-modeled gravitational wave signals have been hampered by transient noise artifacts, or "glitches," in the detectors. In some cases, even high signal-to-noise simulated astrophysical signals have proven difficult to distinguish from glitches, so that essentially any plausible signal could be detected with at most 2-3 $σ$ level confidence. We have put forth the BayesWave algorithm to differentiate between generic gravitational wave transients and glitches, and to provide robust waveform reconstruction and characterization of the astrophysical signals. Here we study BayesWave's capabilities for rejecting glitches while assigning high confidence to detection candidates through analytic approximations to the Bayesian evidence. Analytic results are tested with numerical experiments by adding simulated gravitational wave transient signals to LIGO data collected between 2009 and 2010 and found to be in good agreement.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fisher vs. Bayes : A comparison of parameter estimation techniques for massive black hole binaries to high redshifts with eLISA

Massive black hole binaries are the primary source of gravitational waves (GW) for the future eLISA observatory. The detection and parameter estimation of these sources to high redshift would provide invaluable information on the formation mechanisms of seed black holes, and on the evolution of massive black holes and their host galaxies through cosmic time. The Fisher information matrix has been the standard tool for GW parameter estimation in the last two decades. However, recent studies have questioned the validity of using the Fisher matrix approach. For example, the Fisher matrix approach sometimes predicts errors of $\geq100\%$ in the estimation of parameters such as the luminosity distance and sky position. With advances in computing power, Bayesian inference is beginning to replace the Fisher matrix approximation in parameter estimation studies. In this work, we conduct a Bayesian inference analysis for 120 sources situated at redshifts of between $0.1\leq z\leq 13.2$, and compare the results with those from a Fisher matrix analysis. The Fisher matrix results suggest that for this particular selection of sources, eLISA would be unable to localize sources at redshifts of $z\lesssim6$. In contrast, Bayesian inference provides finite error estimations for all sources in the study, and shows that we can establish minimum closest distances for all sources. The study further predicts that we should be capable with eLISA, out to a redshift of at least $z\leq13$, of predicting a maximum error in the chirp mass of $\lesssim 1\%$, the reduced mass of $\lesssim20\%$, the time to coalescence of 2 hours, and to a redshift of $z\sim5$, the inclination of the source with a maximum error of $\sim60$ degrees.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Phase-coherent mapping of gravitational-wave backgrounds using ground-based laser interferometers

We extend the formalisms developed in Gair et al. and Cornish and van Haasteren to create maps of gravitational-wave backgrounds using a network of ground-based laser interferometers. We show that in contrast to pulsar timing arrays, which are insensitive to half of the gravitational-wave sky (the curl modes), a network of ground-based interferometers is sensitive to both the gradient and curl components of the background. The spatial separation of a network of interferometers, or of a single interferometer at different times during its rotational and orbital motion around the Sun, allows for recovery of both components. We derive expressions for the response functions of a laser interferometer in the small-antenna limit, and use these expressions to calculate the overlap reduction function for a pair of interferometers. We also construct maximum-likelihood estimates of the + and x-polarization modes of the gravitational-wave sky in terms of the response matrix for a network of ground-based interferometers, evaluated at discrete times during Earth's rotational and orbital motion around the Sun. We demonstrate the feasibility of this approach for some simple simulated backgrounds (a single point source and spatially-extended distributions having only grad or curl components), calculating maximum-likelihood sky maps and uncertainty maps based on the (pseudo)inverse of the response matrix. The distinction between this approach and standard methods for mapping gravitational-wave power is also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Towards Robust Gravitational Wave Detection with Pulsar Timing Arrays

Precision timing of highly stable milli-second pulsars is a promising technique for the detection of very low frequency sources of gravitational waves. In any single pulsar, a stochastic gravitational wave signal appears as an additional source of timing noise that can be absorbed by the noise model, and so it is only by considering the coherent response across a network of pulsars that the signal can be distinguished from other sources of noise. In the limit where there are many gravitational wave sources in the sky, or many pulsars in the array, the signals produce a unique tensor correlation pattern that depends only on the angular separation between each pulsar pair. It is this distinct fingerprint that is used to search for gravitational waves using pulsar timing arrays. Here we consider how the prospects for detection are diminished when the statistical isotropy of the timing array or the gravitational wave signal is broken by having a finite number of pulsars and a finite number of sources. We find the standard tensor-correlation analysis to be remarkably robust, with a mild impact on detectability compared to the isotropic limit. Only when there are very few sources and very few pulsars does the standard analysis begin to fail. Having established that the tensor correlations are a robust signature for detection, we study the use of "sky-scrambles" to break the correlations as a way to increase confidence in a detection. This approach is analogous to the use of "time-slides" in the analysis of data from ground based interferometric detectors.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

When is a gravitational-wave signal stochastic?

We discuss the detection of gravitational-wave backgrounds in the context of Bayesian inference and suggest a practical definition of what it means for a signal to be considered stochastic---namely, that the Bayesian evidence favors a stochastic signal model over a deterministic signal model. A signal can further be classified as Gaussian-stochastic if a Gaussian signal model is favored. In our analysis we use Bayesian model selection to choose between several signal and noise models for simulated data consisting of uncorrelated Gaussian detector noise plus a superposition of sinusoidal signals from an astrophysical population of gravitational-wave sources. For simplicity, we consider co-located and co-aligned detectors with white detector noise, but the method can be extended to more realistic detector configurations and power spectra. The general trend we observe is that a deterministic model is favored for small source numbers, a non-Gaussian stochastic model is preferred for intermediate source numbers, and a Gaussian stochastic model is preferred for large source numbers. However, there is very large variation between individual signal realizations, leading to fuzzy boundaries between the three regimes. We find that a hybrid, trans-dimensional model comprised of a deterministic signal model for individual bright sources and a Gaussian-stochastic signal model for the remaining confusion background outperforms all other models in most instances.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

BayesLine: Bayesian Inference for Spectral Estimation of Gravitational Wave Detector Noise

Gravitational wave data from ground-based detectors is dominated by instrument noise. Signals will be comparatively weak, and our understanding of the noise will influence detection confidence and signal characterization. Mis-modeled noise can produce large systematic biases in both model selection and parameter estimation. Here we introduce a multi-component, variable dimension, parameterized model to describe the Gaussian-noise power spectrum for data from ground-based gravitational wave interferometers. Called BayesLine, the algorithm models the noise power spectral density using cubic splines for smoothly varying broad-band noise and Lorentzians for narrow-band line features in the spectrum. We describe the algorithm and demonstrate its performance on data from the fifth and sixth LIGO science runs. Once fully integrated into LIGO/Virgo data analysis software, BayesLine will produce accurate spectral estimation and provide a means for marginalizing inferences drawn from the data over all plausible noise spectra.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Mapping the nano-Hertz gravitational wave sky

We describe a new method for extracting gravitational wave signals from pulsar timing data. We show that any gravitational wave signal can be decomposed into an orthogonal set of sky maps, with the number of maps equal to the number of pulsars in the timing array. These maps may be used as a basis to construct gravitational wave templates for any type of source, including collections of point sources. A variant of the standard Hellings-Downs correlation analysis is recovered for statistically isotropic signals. The template based approach allows us to probe potential anisotropies in the signal and produce maps of the gravitational wave sky.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Multivariate Classification with Random Forests for Gravitational Wave Searches of Black Hole Binary Coalescence

Searches for gravitational waves produced by coalescing black hole binaries with total masses $\gtrsim25\,$M$_\odot$ use matched filtering with templates of short duration. Non-Gaussian noise bursts in gravitational wave detector data can mimic short signals and limit the sensitivity of these searches. Previous searches have relied on empirically designed statistics incorporating signal-to-noise ratio and signal-based vetoes to separate gravitational wave candidates from noise candidates. We report on sensitivity improvements achieved using a multivariate candidate ranking statistic derived from a supervised machine learning algorithm. We apply the random forest of bagged decision trees technique to two separate searches in the high mass $\left( \gtrsim25\,\mathrm{M}_\odot \right)$ parameter space. For a search which is sensitive to gravitational waves from the inspiral, merger, and ringdown (IMR) of binary black holes with total mass between $25\,$M$_\odot$ and $100\,$M$_\odot$, we find sensitive volume improvements as high as $70_{\pm 13}-109_{\pm 11}$\% when compared to the previously used ranking statistic. For a ringdown-only search which is sensitive to gravitational waves from the resultant perturbed intermediate mass black hole with mass roughly between $10\,$M$_\odot$ and $600\,$M$_\odot$, we find sensitive volume improvements as high as $61_{\pm 4}-241_{\pm 12}$\% when compared to the previously used ranking statistic. We also report how sensitivity improvements can differ depending on mass regime, mass ratio, and available data quality information. Finally, we describe the techniques used to tune and train the random forest classifier that can be generalized to its use in other searches for gravitational waves.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Detecting a Stochastic Gravitational Wave Background in the presence of a Galactic Foreground and Instrument Noise

Detecting a stochastic gravitational wave background requires that we first understand and model any astrophysical foregrounds. In the millihertz frequency band, the predominate foreground signal will be from unresolved white dwarf binaries in the galaxy. We build on our previous work to show that a stochastic gravitational wave background can be detected in the presence of both instrument noise and a galactic confusion foreground. The key to our approach is accurately modeling the spectra for each of the various signal components. We simulate data for a gigameter Laser Interferometer Space Antenna (LISA) operating in the mHz frequency band detector operating with both 6- and 4-links. We obtain posterior distribution functions for the instrument noise parameters, the galaxy level and modulation parameters, and the stochastic background energy density. We find that we are able to detect a scale-invariant stochastic background with energy density as low as Omega_gw = 2e-13 for a 6-link interferometer and Omega_gw = 5e-13 for a 4-link interferometer with one year of data.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fast Fisher Matrices and Lazy Likelihoods

Theoretical studies in gravitational wave astronomy often require the calculation of Fisher Information Matrices and Likelihood functions, which in a direct approach entail the costly step of computing gravitational waveforms. Here I describe an alternative technique that sidesteps the need to compute full waveforms, resulting in significant computational savings. I describe how related techniques can be used to speed up Bayesian inference applied to real gravitational wave data.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Pulsar Timing Array Analysis for Black Hole Backgrounds

An astrophysical population of supermassive black hole binaries is thought to be the strongest source of gravitational waves in the frequency range covered by Pulsar Timing Arrays (PTAs). A potential cause for concern is that the standard cross-correlation method used in PTA data analysis assumes that the signals are isotropically distributed and Gaussian random, while the signals from a black hole population are likely to be anisotropic and deterministic. Here we argue that while the conventional analysis is not optimal, it is not hopeless either, as the standard Hellings-Downs correlation curve turns out to hold for point sources, and the small effective number of signal samples blurs the distinction between Gaussian and deterministic signals. Possible improvements to the standard cross-correlation analysis that account for the anisotropy of the signal are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Towards a unified treatment of gravitational-wave data analysis

We present a unified description of gravitational-wave data analysis that unites the template-based analysis used to detect deterministic signals from well-modeled sources, such as binary-black-hole mergers, with the cross-correlation analysis used to detect stochastic gravitational-wave backgrounds. We also discuss the connection between template-based analyses and those that target poorly-modeled bursts of gravitational waves, and suggest a new approach for detecting burst signals.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Applying Bayesian Inference to the first International Pulsar Timing Array data challenge

This is a very brief summary of the techniques I used to analyze the IPTA challenge 1 data sets. I tried many things, and more failed than succeeded, but in the end I found two approaches that appear to work based on tests done using the open data sets. One approach works directly with the time domain data, and the other works with a specially constructed Fourier transform of the data. The raw data was run through TEMPO2 to produce reduced timing residuals for the analysis. Standard Markov Chain Monte Carlo techniques were used to produce samples from the posterior distribution function for the model parameters. The model parameters include the gravitational wave amplitude and spectral slope, and the white noise amplitude for each pulsar in the array. While red timing noise was only included in Dataset 3, I found that it was necessary to include effective red noise in all the analyses to account for some of the spurious effects introduced by the TEMPO2 timing fit. This added an additional amplitude and slope parameter for each pulsar, so my overall model for the 36 pulsars residuals has 110 parameters. As an alternative to using an effective red noise model, I also tried to simultaneously re-fit the timing model model while looking for the gravitational wave signal, but for reasons that are not yet clear, this approach was not very successful. I comment briefly on ways in which the algorithms could be improved. My best estimates for the gravitational wave amplitudes in the three closed (blind) data sets are: (1) $A=(7.3\pm 1.0)\times 10^{-15}$; (2) $A=(5.7\pm 0.6)\times 10^{-14}$; and (3) $A=(4.6\pm 1.3)\times 10^{-15}$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Constraints on the Topology of the Universe: Extension to General Geometries

We present an update to the search for a non-trivial topology of the universe by searching for matching circle pairs in the cosmic microwave background using the WMAP 7 year data release. We extend the exisiting bounds to encompass a wider range of possible topologies by searching for matching circle pairs with opening angles 10 degree < α< 90 degree and separation angles 11 degree < θ< 180 degree. The extended search reveal two small anomalous regions in the CMB sky. Numerous pairs of well-matched circles are found where both circles pass through one or the other of those regions. As this is not the signature of any known manifold, but is a likely consequence of contamination in those sky regions, we repeat the search excluding circle pairs where both pass through either of the two regions. We then find no statistically significant pairs of matched circles, and so no hints of a non-trivial topology. The absence of matched circles increases the lower limit on the length of the shortest closed null geodesic that self-intersects at our location in the universe (equivalently the injectivity radius at our location) to 98.5% of the diameter of the last scattering surface or approximately 26 Gpc. It extends the limit to any manifolds in which the intersecting arcs of said geodesic form an angle greater than 10^o.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

eLISA: Astrophysics and cosmology in the millihertz regime

This document introduces the exciting and fundamentally new science and astronomy that the European New Gravitational Wave Observatory (NGO) mission (derived from the previous LISA proposal) will deliver. The mission (which we will refer to by its informal name "eLISA") will survey for the first time the low-frequency gravitational wave band (about 0.1 mHz to 1 Hz), with sufficient sensitivity to detect interesting individual astrophysical sources out to z = 15. The eLISA mission will discover and study a variety of cosmic events and systems with high sensitivity: coalescences of massive black holes binaries, brought together by galaxy mergers; mergers of earlier, less-massive black holes during the epoch of hierarchical galaxy and black-hole growth; stellar-mass black holes and compact stars in orbits just skimming the horizons of massive black holes in galactic nuclei of the present era; extremely compact white dwarf binaries in our Galaxy, a rich source of information about binary evolution and about future Type Ia supernovae; and possibly most interesting of all, the uncertain and unpredicted sources, for example relics of inflation and of the symmetry-breaking epoch directly after the Big Bang. eLISA's measurements will allow detailed studies of these signals with high signal-to-noise ratio, addressing most of the key scientific questions raised by ESA's Cosmic Vision programme in the areas of astrophysics and cosmology. They will also provide stringent tests of general relativity in the strong-field dynamical regime, which cannot be probed in any other way. This document not only describes the science but also gives an overview on the mission design and orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gravitational Wave Astronomy: Needle in a Haystack

A world-wide array of highly sensitive interferometers stands poised to usher in a new era in astronomy with the first direct detection of gravitational waves. The data from these instruments will provide a unique perspective on extreme astrophysical phenomena such as neutron stars and black holes, and will allow us to test Einstein's theory of gravity in the strong field, dynamical regime. To fully realize these goals we need to solve some challenging problems in signal processing and inference, such as finding rare and weak signals that are buried in non-stationary and non-Gaussian instrument noise, dealing with high-dimensional model spaces, and locating what are often extremely tight concentrations of posterior mass within the prior volume. Gravitational wave detection using space based detectors and Pulsar Timing Arrays bring with them the additional challenge of having to isolate individual signals that overlap one another in both time and frequency. Promising solutions to these problems will be discussed, along with some of the challenges that remain.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Low-frequency gravitational-wave science with eLISA/NGO

We review the expected science performance of the New Gravitational-Wave Observatory (NGO, a.k.a. eLISA), a mission under study by the European Space Agency for launch in the early 2020s. eLISA will survey the low-frequency gravitational-wave sky (from 0.1 mHz to 1 Hz), detecting and characterizing a broad variety of systems and events throughout the Universe, including the coalescences of massive black holes brought together by galaxy mergers; the inspirals of stellar-mass black holes and compact stars into central galactic black holes; several millions of ultracompact binaries, both detached and mass transferring, in the Galaxy; and possibly unforeseen sources such as the relic gravitational-wave radiation from the early Universe. eLISA's high signal-to-noise measurements will provide new insight into the structure and history of the Universe, and they will test general relativity in its strong-field dynamical regime.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Prospects for observing ultra-compact binaries with space-based gravitational wave interferometers and optical telescopes

Space-based gravitational wave interferometers are sensitive to the galactic population of ultra-compact binaries. An important subset of the ultra-compact binary population are those stars that can be individually resolved by both gravitational wave interferometers and electromagnetic telescopes. The aim of this paper is to quantify the multi-messenger potential of space-based interferometers with arm-lengths between 1 and 5 Gm. The Fisher Information Matrix is used to estimate the number of binaries from a model of the Milky Way which are localized on the sky by the gravitational wave detector to within 1 and 10 square degrees and bright enough to be detected by a magnitude limited survey. We find, depending on the choice of GW detector characteristics, limiting magnitude, and observing strategy, that up to several hundred gravitational wave sources could be detected in electromagnetic follow-up observations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Characterizing Spinning Black Hole Binaries in Eccentric Orbits with LISA

The Laser Interferometer Space Antenna (LISA) is designed to detect gravitational wave signals from astrophysical sources, including those from coalescing binary systems of compact objects such as black holes. Colliding galaxies have central black holes that sink to the center of the merged galaxy and begin to orbit one another and emit gravitational waves. Some galaxy evolution models predict that the binary black hole system will enter the LISA band with significant orbital eccentricity, while other models suggest that the orbits will already have circularized. Using a full seventeen parameter waveform model that includes the effects of orbital eccentricity, spin precession and higher harmonics, we investigate how well the source parameters can be inferred from simulated LISA data. Defining the reference eccentricity as the value one year before merger, we find that for typical LISA sources, it will be possible to measure the eccentricity to an accuracy of parts in a thousand. The accuracy with which the eccentricity can be measured depends only very weakly on the eccentricity, making it possible to distinguish circular orbits from those with very small eccentricities. LISA measurements of the orbital eccentricity can help constraints theories of galaxy mergers in the early universe. Failing to account for the eccentricity in the waveform modeling can lead to a loss of signal power and bias the estimation of parameters such as the black hole masses and spins.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Measuring parameters of massive black hole binaries with partially aligned spins

The future space-based gravitational wave detector LISA will be able to measure parameters of coalescing massive black hole binaries, often to extremely high accuracy. Previous work has demonstrated that the black hole spins can have a strong impact on the accuracy of parameter measurement. Relativistic spin-induced precession modulates the waveform in a manner which can break degeneracies between parameters, in principle significantly improving how well they are measured. Recent studies have indicated, however, that spin precession may be weak for an important subset of astrophysical binary black holes: those in which the spins are aligned due to interactions with gas. In this paper, we examine how well a binary's parameters can be measured when its spins are partially aligned and compare results using waveforms that include higher post-Newtonian harmonics to those that are truncated at leading quadrupole order. We find that the weakened precession can substantially degrade parameter estimation. This degradation is particularly devastating for the extrinsic parameters sky position and distance. Absent higher harmonics, LISA typically localizes the sky position of a nearly aligned binary a factor of $\sim 6$ less accurately than for one in which the spin orientations are random. Our knowledge of a source's sky position will thus be worst for the gas-rich systems which are most likely to produce electromagnetic counterparts. Fortunately, higher harmonics of the waveform can make up for this degradation. By including harmonics beyond the quadrupole in our waveform model, we find that the accuracy with which most of the binary's parameters are measured can be substantially improved. In some cases, parameters can be measured as well in partially aligned binaries as they can be when the binary spins are random.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Computing waveforms for spinning compact binaries in quasi-eccentric orbits

Several scenarios have been proposed in which the orbits of binary black holes enter the band of a gravitational wave detector with significant eccentricity. To avoid missing these signals or biasing parameter estimation it is important that we consider waveform models that account for eccentricity. The ingredients needed to compute post-Newtonian (PN) waveforms produced by spinning black holes inspiralling on quasi-eccentric orbits have been available for almost two decades at 2 PN order, and this work has recently been extended to 2.5 PN order. However, the computational cost of directly implementing these waveforms is high, requiring many steps per orbit to evolve the system of coupled differential equations. Here we employ the standard techniques of a separation of timescales and a generalized Keplerian parameterization of the orbits to produce efficient waveforms describing spinning black hole binaries with arbitrary masses and spins on quasi-eccentric orbits to 1.5 PN order. We separate the fast orbital timescale from the slow spin-orbit precession timescale by solving for the orbital motion in a non-interial frame of reference that follows the orbital precession. We outline a scheme for extending our approach to higher post-Newtonian order.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Discriminating between a Stochastic Gravitational Wave Background and Instrument Noise

The detection of a stochastic background of gravitational waves could significantly impact our understanding of the physical processes that shaped the early Universe. The challenge lies in separating the cosmological signal from other stochastic processes such as instrument noise and astrophysical foregrounds. One approach is to build two or more detectors and cross correlate their output, thereby enhancing the common gravitational wave signal relative to the uncorrelated instrument noise. When only one detector is available, as will likely be the case with the Laser Interferometer Space Antenna (LISA), alternative analysis techniques must be developed. Here we show that models of the noise and signal transfer functions can be used to tease apart the gravitational and instrument noise contributions. We discuss the role of gravitational wave insensitive "null channels" formed from particular combinations of the time delay interferometry, and derive a new combination that maintains this insensitivity for unequal arm length detectors. We show that, in the absence of astrophysical foregrounds, LISA could detect signals with energy densities as low as $Ω_{\rm gw} = 6 \times 10^{-13}$ with just one month of data. We describe an end-to-end Bayesian analysis pipeline that is able to search for, characterize and assign confidence levels for the detection of a stochastic gravitational wave background, and demonstrate the effectiveness of this approach using simulated data from the third round of Mock LISA Data Challenges.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Pulsar Timing Array Observations of Massive Black Hole Binaries

Pulsar timing is a promising technique for detecting low frequency sources of gravitational waves. Historically the focus has been on the detection of diffuse stochastic backgrounds, such as those formed from the superposition of weak signals from a population of binary black holes. More recently, attention has turned to members of the binary population that are nearer and brighter, which stand out from the crowd and can be individually resolved. Here we show that the timing data from an array of pulsars can be used to recover the physical parameters describing an individual black hole binary to good accuracy, even for moderately strong signals. A novel aspect of our analysis is that we include the distance to each pulsar as a search parameter, which allows us to utilize the full gravitational wave signal. This doubles the signal power, improves the sky location determination by an order of magnitude, and allows us to extract the mass and the distance to the black hole binary.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Separating Gravitational Wave Signals from Instrument Artifacts

Central to the gravitational wave detection problem is the challenge of separating features in the data produced by astrophysical sources from features produced by the detector. Matched filtering provides an optimal solution for Gaussian noise, but in practice, transient noise excursions or ``glitches'' complicate the analysis. Detector diagnostics and coincidence tests can be used to veto many glitches which may otherwise be misinterpreted as gravitational wave signals. The glitches that remain can lead to long tails in the matched filter search statistics and drive up the detection threshold. Here we describe a Bayesian approach that incorporates a more realistic model for the instrument noise allowing for fluctuating noise levels that vary independently across frequency bands, and deterministic ``glitch fitting'' using wavelets as ``glitch templates'', the number of which is determined by a trans-dimensional Markov chain Monte Carlo algorithm. We demonstrate the method's effectiveness on simulated data containing low amplitude gravitational wave signals from inspiraling binary black hole systems, and simulated non-stationary and non-Gaussian noise comprised of a Gaussian component with the standard LIGO/Virgo spectrum, and injected glitches of various amplitude, prevalence, and variety. Glitch fitting allows us to detect significantly weaker signals than standard techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Mock LISA Data Challenges: from Challenge 3 to Challenge 4

The Mock LISA Data Challenges are a program to demonstrate LISA data-analysis capabilities and to encourage their development. Each round of challenges consists of one or more datasets containing simulated instrument noise and gravitational waves from sources of undisclosed parameters. Participants analyze the datasets and report best-fit solutions for the source parameters. Here we present the results of the third challenge, issued in Apr 2008, which demonstrated the positive recovery of signals from chirping Galactic binaries, from spinning supermassive--black-hole binaries (with optimal SNRs between ~ 10 and 2000), from simultaneous extreme-mass-ratio inspirals (SNRs of 10-50), from cosmic-string-cusp bursts (SNRs of 10-100), and from a relatively loud isotropic background with Omega_gw(f) ~ 10^-11, slightly below the LISA instrument noise.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Characterizing the Gravitational Wave Signature from Cosmic String Cusps

Cosmic strings are predicted to form kinks and cusps that travel along the string at close to the speed of light. These disturbances are radiated away as highly beamed gravitational waves that produce a burst like pulse as the cone of emission sweeps past an observer. Gravitational wave detectors such as the Laser Interferometer Space Antenna (LISA) and the Laser Interferometer Gravitational wave Observatory (LIGO) will be capable of detecting these bursts for a wide class of string models. Such a detection would illuminate the fields of string theory, cosmology, and relativity. Here we develop template based Markov Chain Monte Carlo (MCMC) techniques that can efficiently detect and characterize the signals from cosmic string cusps. We estimate how well the signal parameters can be recovered by the advanced LIGO-Virgo network and the LISA detector using a combination of MCMC and Fisher matrix techniques. We also consider joint detections by the ground and space based instruments. We show that a parallel tempered MCMC approach can detect and characterize the signals from cosmic string cusps, and we demonstrate the utility of this approach on simulated data from the third round of Mock LISA Data Challenges (MLDCs).

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

The Search for Massive Black Hole Binaries with LISA

In this work we focus on the search and detection of Massive black hole binary (MBHB) systems, including systems at high redshift. As well as expanding on previous works where we used a variant of Markov Chain Monte Carlo (MCMC), called Metropolis-Hastings Monte Carlo, with simulated annealing, we introduce a new search method based on frequency annealing which leads to a more rapid and robust detection. We compare the two search methods on systems where we do and do not see the merger of the black holes. In the non-merger case, we also examine the posterior distribution exploration using a 7-D MCMC algorithm. We demonstrate that this method is effective in dealing with the high correlations between parameters, has a higher acceptance rate than previously proposed methods and produces posterior distribution functions that are close to the prediction from the Fisher Information matrix. Finally, after carrying out searches where there is only one binary in the data stream, we examine the case where two black hole binaries are present in the same data stream. We demonstrate that our search algorithm can accurately recover both binaries, and more importantly showing that we can safely extract the MBHB sources without contaminating the rest of the data stream.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Darwin Meets Einstein: LISA Data Analysis Using Genetic Algorithms

This work presents the first application of the method of Genetic Algorithms (GAs) to data analysis for the Laser Interferometer Space Antenna (LISA). In the low frequency regime of the LISA band there are expected to be tens of thousands galactic binary systems that will be emitting gravitational waves detectable by LISA. The challenge of parameter extraction of such a large number of sources in the LISA data stream requires a search method that can efficiently explore the large parameter spaces involved. As signals of many of these sources will overlap, a global search method is desired. GAs represent such a global search method for parameter extraction of multiple overlapping sources in the LISA data stream. We find that GAs are able to correctly extract source parameters for overlapping sources. Several optimizations of a basic GA are presented with results derived from applications of the GA searches to simulated LISA data.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Detecting the Cosmic Gravitational Wave Background with the Big Bang Observer

The detection of the Cosmic Microwave Background Radiation (CMB) was one of the most important cosmological discoveries of the last century. With the development of interferometric gravitational wave detectors, we may be in a position to detect the gravitational equivalent of the CMB in this century. The Cosmic Gravitational Background (CGB) is likely to be isotropic and stochastic, making it difficult to distinguish from instrument noise. The contribution from the CGB can be isolated by cross-correlating the signals from two or more independent detectors. Here we extend previous studies that considered the cross-correlation of two Michelson channels by calculating the optimal signal to noise ratio that can be achieved by combining the full set of interferometry variables that are available with a six link triangular interferometer. In contrast to the two channel case, we find that the relative orientation of a pair of coplanar detectors does not affect the signal to noise ratio. We apply our results to the detector design described in the Big Bang Observer (BBO) mission concept study and find that BBO could detect a background with $Ω_{gw} > 2.2 \times 10^{-17}$.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

The Effects of Orbital Motion on LISA Time Delay Interferometry

In an effort to eliminate laser phase noise in laser interferometer spaceborne gravitational wave detectors, several combinations of signals have been found that allow the laser noise to be canceled out while gravitational wave signals remain. This process is called time delay interferometry (TDI). In the papers that defined the TDI variables, their performance was evaluated in the limit that the gravitational wave detector is fixed in space. However, the performance depends on certain symmetries in the armlengths that are available if the detector is fixed in space, but that will be broken in the actual rotating and flexing configuration produced by the LISA orbits. In this paper we investigate the performance of these TDI variables for the real LISA orbits. First, addressing the effects of rotation, we verify Daniel Shaddock's result that the Sagnac variables will not cancel out the laser phase noise, and we also find the same result for the symmetric Sagnac variable. The loss of the latter variable would be particularly unfortunate since this variable also cancels out gravitational wave signal, allowing instrument noise in the detector to be isolated and measured. Fortunately, we have found a set of more complicated TDI variables, which we call Delta-Sagnac variables, one of which accomplishes the same goal as the symmetric Sagnac variable to good accuracy. Finally, however, as we investigate the effects of the flexing of the detector arms due to non-circular orbital motion, we show that all variables, including the interferometer variables, which survive the rotation-induced loss of direction symmetry, will not completely cancel laser phase noise when the armlengths are changing with time. This unavoidable problem will place a stringent requirement on laser stability of 5 Hz per root Hz.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

The LISA Response Function

The orbital motion of the Laser Interferometer Space Antenna (LISA) introduces modulations into the observed gravitational wave signal. These modulations can be used to determine the location and orientation of a gravitational wave source. The complete LISA response to an arbitary gravitational wave is derived using a coordinate free approach in the transverse-traceless gauge. The general response function reduces to that found by Cutler (PRD 57, 7089 1998) for low frequency, monochromatic plane waves. Estimates of the noise in the detector are found to be complicated by the time variation of the interferometer arm lengths.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Detecting a stochastic gravitational wave background with the Laser Interferometer Space Antenna

The random superposition of many weak sources will produce a stochastic background of gravitational waves that may dominate the response of the LISA (Laser Interferometer Space Antenna) gravitational wave observatory. Unless something can be done to distinguish between a stochastic background and detector noise, the two will combine to form an effective noise floor for the detector. Two methods have been proposed to solve this problem. The first is to cross-correlate the output of two independent interferometers. The second is an ingenious scheme for monitoring the instrument noise by operating LISA as a Sagnac interferometer. Here we derive the optimal orbital alignment for cross-correlating a pair of LISA detectors, and provide the first analytic derivation of the Sagnac sensitivity curve.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Mapping the gravitational wave background

The gravitational wave sky is expected to have isolated bright sources superimposed on a diffuse gravitational wave background. The background radiation has two components: a confusion limited background from unresolved astrophysical sources; and a cosmological component formed during the birth of the universe. A map of the gravitational wave background can be made by sweeping a gravitational wave detector across the sky. The detector output is a complicated convolution of the sky luminosity distribution, the detector response function and the scan pattern. Here we study the general de-convolution problem, and show how LIGO (Laser Interferometric Gravitational Observatory) and LISA (Laser Interferometer Space Antenna) can be used to detect anisotropies in the gravitational wave background.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Circles in the Sky: Finding Topology with the Microwave Background Radiation

If the universe is finite and smaller than the distance to the surface of last scatter, then the signature of the topology of the universe is writ on the microwave background sky. Previous efforts to search for this topology have focused on one particular model: a toroidal flat universe. We show how both the high degree of spatial symmetry of this topology and the integrability of its geodesics make it unreliable as a paradigmatic example, and discuss why topology on scales significantly smaller than the horizon are not ruled out by previous analyses focussing on this special case. We show that in these small universes the microwave background will be identified at the intersections of the surface of last scattering as seen by different ``copies'' of the observer. Since the surface of last scattering is a sphere, these intersections will be circles, regardless of the background geometry or topology. We therefore propose a statistic that is sensitive to all small finite homogeneous topologies. Here, small means that the distance to the surface of last scatter is smaller than the ``periodicity scale'' of the universe.

0 citations0 reviewsNo code linkedOpen access