Data analysis strategies for the detection of gravitational waves in non-Gaussian noise

In order to analyze data produced by the kilometer-scale gravitational wave detectors that will begin operation early next century, one needs to develop robust statistical tools capable of extracting weak signals from the detector noise. This noise will likely have non-stationary and non-Gaussian components. To facilitate the construction of robust detection techniques, I present a simple two-component noise model that consists of a background of Gaussian noise as well as stochastic noise bursts. The optimal detection statistic obtained for such a noise model incorporates a natural veto which suppresses spurious events that would be caused by the noise bursts. When two detectors are present, I show that the optimal statistic for the non-Gaussian noise model can be approximated by a simple coincidence detection strategy. For simulated detector noise containing noise bursts, I compare the operating characteristics of (i) a locally optimal detection statistic (which has nearly-optimal behavior for small signal amplitudes) for the non-Gaussian noise model, (ii) a standard coincidence-style detection strategy, and (iii) the optimal statistic for Gaussian noise.Comment: 5 pages RevTeX, 4 figure

Is the squeezing of relic gravitational waves produced by inflation detectable?

Grishchuk has shown that the stochastic background of gravitational waves produced by an inflationary phase in the early Universe has an unusual property: it is not a stationary Gaussian random process. Due to squeezing, the phases of the different waves are correlated in a deterministic way, arising from the process of parametric amplification that created them. The resulting random process is Gaussian but non-stationary. This provides a unique signature that could in principle distinguish a background created by inflation from stationary stochastic backgrounds created by other types of processes. We address the question: could this signature be observed with a gravitational wave detector? Sadly, the answer appears to be "no": an experiment which could distinguish the non-stationary behavior would have to last approximately the age of the Universe at the time of measurement. This rules out direct detection by ground and space based gravitational wave detectors, but not indirect detections via the electromagnetic Cosmic Microwave Background Radiation (CMBR).Comment: 17 pages, 4 Postscript figures, uses revtex, psfig, to be submitted to PRD, minor revisions - appendix B clarified, corrected typos, added reference

Measuring gravitational waves from binary black hole coalescences: II. the waves' information and its extraction, with and without templates

We discuss the extraction of information from detected binary black hole (BBH) coalescence gravitational waves, focusing on the merger phase that occurs after the gradual inspiral and before the ringdown. Our results are: (1) If numerical relativity simulations have not produced template merger waveforms before BBH detections by LIGO/VIRGO, one can band-pass filter the merger waves. For BBHs smaller than about 40 solar masses detected via their inspiral waves, the band pass filtering signal to noise ratio indicates that the merger waves should typically be just barely visible in the noise for initial and advanced LIGO interferometers. (2) We derive an optimized (maximum likelihood) method for extracting a best-fit merger waveform from the noisy detector output; one "perpendicularly projects" this output onto a function space (specified using wavelets) that incorporates our prior knowledge of the waveforms. An extension of the method allows one to extract the BBH's two independent waveforms from outputs of several interferometers. (3) If numerical relativists produce codes for generating merger templates but running the codes is too expensive to allow an extensive survey of the merger parameter space, then a coarse survey of this parameter space, to determine the ranges of the several key parameters and to explore several qualitative issues which we describe, would be useful for data analysis purposes. (4) A complete set of templates could be used to test the nonlinear dynamics of general relativity and to measure some of the binary parameters. We estimate the number of bits of information obtainable from the merger waves (about 10 to 60 for LIGO/VIRGO, up to 200 for LISA), estimate the information loss due to template numerical errors or sparseness in the template grid, and infer approximate requirements on template accuracy and spacing.Comment: 33 pages, Rextex 3.1 macros, no figures, submitted to Phys Rev

Fluctuations of the vacuum energy density of quantum fields in curved spacetime via generalized zeta functions

For quantum fields on a curved spacetime with an Euclidean section, we derive a general expression for the stress energy tensor two-point function in terms of the effective action. The renormalized two-point function is given in terms of the second variation of the Mellin transform of the trace of the heat kernel for the quantum fields. For systems for which a spectral decomposition of the wave opearator is possible, we give an exact expression for this two-point function. Explicit examples of the variance to the mean ratio $\Delta' = (-^2)/(^2)$ of the vacuum energy density $\rho$ of a massless scalar field are computed for the spatial topologies of $R^d\times S^1$ and $S^3$, with results of $\Delta'(R^d\times S^1) =(d+1)(d+2)/2$, and $\Delta'(S^3) = 111$ respectively. The large variance signifies the importance of quantum fluctuations and has important implications for the validity of semiclassical gravity theories at sub-Planckian scales. The method presented here can facilitate the calculation of stress-energy fluctuations for quantum fields useful for the analysis of fluctuation effects and critical phenomena in problems ranging from atom optics and mesoscopic physics to early universe and black hole physics.Comment: Uses revte

Facilitators and Barriers to Prescribing PreExposure Prophylaxis (PrEP) for the Prevention of HIV

Background: What is PrEP and who gets it? PrEP is the use of medication by individuals to prevent HIV contraction, approved in 2012 after demonstrating safety and efficacy in the iPrEx study and Partners PrEP2 trials. HIV infection risk is 92% lower in patients using PrEP. Truvada®, a combination of tenofovir and emtricitabine taken orally daily, is the only approved PrEP regimen and is intended to compliment other prevention strategies such as condoms. HIV negative-individuals at risk for exposure to HIV have been identified as men who have sex with men (MSM), IV drug users, heterosexuals who have unprotected sex with partners of unknown HIV status, and those in serodiscordant relationships. Barriers to PrEP Implementation PrEP is effective when patients adhere; however, both the medical community and some high-risk populations have been slow to adopt it as an HIV prevention strategy. Surveys have shown clinicians perceived barriers to PrEP such as adverse side effects, viral drug resistance, increased high-risk behavior, cost, and training. HIV in Vermont New diagnoses of HIV among Vermont residents has remained relatively stable over the last twenty years. Vermont CARES, a non-profit, offers free and anonymous HIV tests and in-person risk-reduction counseling. Clients are increasingly asking about PrEP as a prevention strategy, but the response from the medical community is difficult to ascertain.https://scholarworks.uvm.edu/comphp_gallery/1235/thumbnail.jp

Detecting gravitational waves from precessing binaries of spinning compact objects: Adiabatic limit

Black-hole (BH) binaries with single-BH masses m=5--20 Msun, moving on quasicircular orbits, are among the most promising sources for first-generation ground-based gravitational-wave (GW) detectors. Until now, the development of data-analysis techniques to detect GWs from these sources has been focused mostly on nonspinning BHs. The data-analysis problem for the spinning case is complicated by the necessity to model the precession-induced modulations of the GW signal, and by the large number of parameters needed to characterize the system, including the initial directions of the spins, and the position and orientation of the binary with respect to the GW detector. In this paper we consider binaries of maximally spinning BHs, and we work in the adiabatic-inspiral regime to build families of modulated detection templates that (i) are functions of very few physical and phenomenological parameters, (ii) model remarkably well the dynamical and precessional effects on the GW signal, with fitting factors on average >~ 0.97, but (iii) might require increasing the detection thresholds, offsetting at least partially the gains in the fitting factors. Our detection-template families are quite promising also for the case of neutron-star--black-hole binaries, with fitting factors on average ~ 0.93. For these binaries we also suggest (but do not test) a further template family, which would produce essentially exact waveforms written directly in terms of the physical spin parameters.Comment: 38 pages, 16 figures, RevTeX4. Final PRD version. Lingering typos corrected. Small corrections to GW flux terms as per Blanchet et al., PRD 71, 129902(E)-129904(E) (2005

Avalanche dynamics of radio pulsar glitches

We test statistically the hypothesis that radio pulsar glitches result from an avalanche process, in which angular momentum is transferred erratically from the flywheel-like superfluid in the star to the slowly decelerating, solid crust via spatially connected chains of local, impulsive, threshold-activated events, so that the system fluctuates around a self-organised critical state. Analysis of the glitch population (currently 285 events from 101 pulsars) demonstrates that the size distribution in individual pulsars is consistent with being scale invariant, as expected for an avalanche process. The waiting-time distribution is consistent with being exponential in seven out of nine pulsars where it can be measured reliably, after adjusting for observational limits on the minimum waiting time, as for a constant-rate Poisson process. PSR J0537$-$6910 and PSR J0835$-$4510 are the exceptions; their waiting-time distributions show evidence of quasiperiodicity. In each object, stationarity requires that the rate $\lambda$ equals $- \epsilon \dot{\nu} /$, where $\dot{\nu}$ is the angular acceleration of the crust, $$ is the mean glitch size, and $\epsilon\dot{\nu}$ is the relative angular acceleration of the crust and superfluid. There is no evidence that $\lambda$ changes monotonically with spin-down age. The rate distribution itself is fitted reasonably well by an exponential for $\lambda \geq 0.25 {\rm yr^{-1}}$. For $\lambda < 0.25 {\rm yr^{-1}}$, its exact form is unknown; the exponential overestimates the number of glitching pulsars observed at low $\lambda$, where the limited total observation time exercises a selection bias.Comment: Accepted for publication in the Astrophysical Journa

Stochastic Background Search Correlating ALLEGRO with LIGO Engineering Data

We describe the role of correlation measurements between the LIGO interferometer in Livingston, LA, and the ALLEGRO resonant bar detector in Baton Rouge, LA, in searches for a stochastic background of gravitational waves. Such measurements provide a valuable complement to correlations between interferometers at the two LIGO sites, since they are sensitive in a different, higher, frequency band. Additionally, the variable orientation of the ALLEGRO detector provides a means to distinguish gravitational wave correlations from correlated environmental noise. We describe the analysis underway to set a limit on the strength of a stochastic background at frequencies near 900 Hz using ALLEGRO data and data from LIGO's E7 Engineering Run.Comment: 8 pages, 2 encapsulated PostScript figures, uses IOP class files, submitted to the proceedings of the 7th Gravitational Wave Data Analysis Workshop (which will be published in Classical and Quantum Gravity

An excess power statistic for detection of burst sources of gravitational radiation

We examine the properties of an excess power method to detect gravitational waves in interferometric detector data. This method is designed to detect short-duration (< 0.5 s) burst signals of unknown waveform, such as those from supernovae or black hole mergers. If only the bursts' duration and frequency band are known, the method is an optimal detection strategy in both Bayesian and frequentist senses. It consists of summing the data power over the known time interval and frequency band of the burst. If the detector noise is stationary and Gaussian, this sum is distributed as a chi-squared (non-central chi-squared) deviate in the absence (presence) of a signal. One can use these distributions to compute frequentist detection thresholds for the measured power. We derive the method from Bayesian analyses and show how to compute Bayesian thresholds. More generically, when only upper and/or lower bounds on the bursts duration and frequency band are known, one must search for excess power in all concordant durations and bands. Two search schemes are presented and their computational efficiencies are compared. We find that given reasonable constraints on the effective duration and bandwidth of signals, the excess power search can be performed on a single workstation. Furthermore, the method can be almost as efficient as matched filtering when a large template bank is required. Finally, we derive generalizations of the method to a network of several interferometers under the assumption of Gaussian noise.Comment: 22 pages, 6 figure