Cross-correlation method for intermediate-duration gravitational wave searches associated with gamma-ray bursts
CCross-correlation method for intermediate-duration gravitational wave searchesassociated with gamma-ray bursts
Robert Coyne,
1, 2, ∗ Alessandra Corsi, and Benjamin J. Owen
1, 3 Department of Physics, Texas Tech University, Lubbock, TX 79409-1051 (USA) Department of Physics, George Washington University, Washington, DC 20052 (USA) Department of Physics, The Pennsylvania State University, University Park, Pennsylvania 16802-6300 (USA) (compiled 2016 August 24)Several models of gamma-ray burst progenitors suggest that the gamma-ray event may be followedby gravitational wave signals of 10 –10 seconds duration (possibly accompanying the so-called X-rayafterglow “plateaus”). We term these signals “intermediate-duration” because they are shorter thancontinuous wave signals but longer than signals traditionally considered as gravitational wave bursts,and are difficult to detect with most burst and continuous wave methods. The cross-correlationtechnique proposed by [S. Dhurandhar et al., Phys. Rev. D , 082001 (2008)], which so far has beenused only on continuous wave signals, in principle unifies both burst and continuous wave (as well asmatched filtering and stochastic background) methods, reducing them to different choices of whichdata to correlate on which time scales. Here we perform the first tuning of this cross-correlationtechnique to intermediate-duration signals. We derive theoretical estimates of sensitivity in Gaussiannoise in different limits of the cross-correlation formalism, and compare them to the performanceof a prototype search code on simulated Gaussian-noise data. We estimate that the code is likelyable to detect some classes of intermediate-duration signals (such as the ones described in [A. Corsi& P. M´esz´aros, Astrophys. J., , 1171 (2009)]) from sources located at astrophysically-relevantdistances of several tens of Mpc. I. INTRODUCTION
Over the last decade, the LIGO and Virgo gravita-tional wave (GW) detectors have carried out triggered(or targeted) GW searches in coincidence with Gamma-Ray Bursts (GRBs) and other electromagnetic transients[1–15] as well as persistent electromagnetic sources [16–28]. These searches have traditionally been optimizedto detect well-modeled “chirp” signals from neutron star(NS)-NS and/or black-hole (BH)-NS binary inspirals, un-modeled short ( (cid:46) −
10 s) duration bursts of GWs in as-sociation with electromagnetic transients, and persistent(continuous) GWs from nearby rotating NSs. Searchesbased on methods for a stochastic background have alsobeen adapted to continuous wave targets [23, 29].Methods targeting an as of yet largely unexplored classof “intermediate duration” GW signals have also been de-veloped [30–32] and two so far have led to a search on realdata [13, 33]. Intermediate duration GWs are of specialinterest in several astrophysical scenarios (e.g., [13, 35–43]), and their detectability over a large parameter spaceremains mostly unexplored compared to the more tradi-tional inspiral, burst, or continuous wave signals.In this work, we focus on the possibility of detecting10 − s duration GWs in coincidence with GRBs. Our ∗ [email protected] Those works use “long” to refer to signals of O (10 ) s duration,because these durations are long compared to the O ( (cid:46)
1) s du-ration signals traditionally targeted in burst data analyses. Theterm “very long duration” signals has also been adopted to referto GWs lasting from hours to weeks, e.g. [34]. Here, we use “in-termediate” to put the discussion in the broader context, whichincludes the substantially longer continuous wave signals. study is motivated by the need for a data analysis tech-nique that is optimized to probe some of the long-livedprogenitor scenarios for (long and short) GRBs, such asthe so-called “magnetar model”. The magnetized NS(magnetar) scenario has been invoked to explain X-ray“plateaus” (10 − s-long periods of relatively constantemission) observed in (cid:38)
50% of long, and in several short,GRB afterglows [44–51]. Gravitational collapse leadingto the formation of a NS, in turn, has long been consid-ered an observable source of GWs. In a rotating, newlyborn NS, non-axisymmetric instabilities such as the secu-lar Chandrasekhar-Friedman-Schutz [CFS, 52, 53] insta-bilities, can yield GW emission with high efficiency [54].If the newly born GRB-magnetar emits GWs over theplateau timescale ( ∼ s), GW detectors such as the ad-vanced LIGO (aLIGO) and Virgo detectors may be ableto directly probe the source of the observed prolonged en-ergy injection, and clarify one of the key open questionson the nature of GRB central engines [38, 55].Detecting intermediate-duration GW signals, such asthe ones possibly associated with GRB plateaus, requiressearch techniques that can bridge the gap (both in termsof science reach and signal detection strategies) betweentraditional inspiral/burst searches, and continuous waveor stochastic ones. Traditional short duration inspiraland long duration continuous wave searches make use ofhighly sensitive coherent (and computationally limitedsemi-coherent) techniques that leverage accurate knowl-edge of the expected GW waveform (as a function of a setof physical parameters). Traditional burst and stochas-tic searches, on the other hand, assume little a prioriknowledge of the signal and depend respectively on ex-cess signal power (above the background noise) and cross-correlation of power between interferometers for detec- a r X i v : . [ g r- q c ] A ug tion.Here we address the problem of searching forintermediate-duration, large frequency bandwidth sig-nals by adapting the cross-correlation method of [56].While originally developed in the context of continu-ous waves, the method by [56] encompasses all of theaforementioned traditional search techniques when vari-ous parameters are taken to the appropriate limits, andit shows how to make best use of the information avail-able about each type of signal. (A Bayesian framework ofsimilarly broad relevance was developed later in Cornishand Romano [57], but here like Dhurandhar et al. [56] wepresent an essentially frequentist analysis.) We correctsome small errors in the original formalism of [56], andapply it for the first time to intermediate-duration signalsby developing a code, the performance of which we testedon simulated data. We restrict ourselves to intermediate-duration signals with large frequency bandwidth (such asthe ones described in [38]), since intermediate-durationnarrow band signals have different astrophysical originsand are treated with adaptations of continuous wavesearches (see e.g. [58]).Our paper is organized as follows. In Sec. II we mo-tivate the application of Dhurandhar et al.’s [56] cross-correlation technique to intermediate-duration GWs. InSec. III we describe our notation and assumptions. InSec. IV we briefly re-derive the general statistical be-havior of the cross-correlation method, discuss explic-itly its limits and intermediate regimes, and show howseveral assumptions made in [56] need to be modifiedfor the search of non well-modeled GW transients evolv-ing on 10 − s timescales. In Sec. V we apply thecross-correlation technique to the model of secularly un-stable GRB-magnetars described in [38], thus provid-ing an example of applicability to astrophysically mo-tivated waveforms of intermediate-duration. Finally, inSection VI, we compare our results with other data anal-ysis techniques that have been proposed to search forintermediate-duration GW signals, and give our conclu-sions. II. MOTIVATION FOR ACROSS-CORRELATION SEARCH
GWs signals are typically predicted to have strengthsso close to the level of noise in the detectors that it is nec-essary to filter the interferometer data streams to detectthe real GW events amongst spurious noise events. Whenthe functional form of the predicted GW signal is verywell known (as a function of a set of physical parameters),matched filtering with template waveforms is the optimalstrategy (e.g., [59, 60]). Matched filtering involves com-puting the cross-correlation between the interferometeroutput and a template waveform, weighted inversely bythe noise spectrum of the detector. The signal-to-noiseratio (SNR) is defined as the cross-correlation of the tem-plate with a particular stretch of data divided by the root-mean-squared (rms) value of the cross-correlationof the template with pure detector noise.Usually, a family of templates spanning the possiblerange of parameter values (a so-called template bank)is used in real data analyses. A template bank addsto the search statistics a trial factor, which has to betaken into account when estimating the detection sensi-tivity. A template bank also involves more computationalcost since each template must be cross-correlated withthe data. While the parameters describing the searchtemplates typically vary continuously throughout a fi-nite range of values, a realistic template bank is com-posed of templates, the parameter values of which varyin discrete steps within the allowed range. The “mis-match” between the signal and nearest of the discretetemplates causes some reduction in the expected matchedfilter SNR. Thus, the number of templates to be used in asearch is a compromise between the maximum computa-tional cost one can sustain, and the maximum mismatchthat one is willing to tolerate (e.g., [61–64]).When the maximum sustainable computational costimplies a mismatch such that the loss in SNR reduces thesensitivity of the search to a very limited portion of theparameter space, modifications to the matched filteringstrategy toward sub-optimal techniques are mandatory.In addition, in many cases, the GW signal waveform isnot known well enough for matched filtering. Indeed,even if a very finely spaced discrete template bank isused, a search may fail to detected a signal if the tem-plates do not represent with sufficient accuracy the rele-vant physics. In other words, a realistic search is affectednot only by the mismatch but also by the so-called “fit-ting factor” [65–68], the fractional loss in SNR caused bythe fact that even the best template in a family is only a“fit” to a hypothetical exact gravitational waveform. Inthe context of GWs from compact binaries, where numer-ical relativity can be used to quantify the fitting factorof phenomenological waveforms used to construct tem-plate banks for matched filter searches (e.g., [69]), it hasbeen estimated that fitting factors <
3% are needed toachieve detection efficiencies >
90% (see e.g. [65, 70]).Indeed, matched filtering is by construction highly likelyto miss a signal even for moderately bad fitting factors.On the other hand, sub-optimal (less sensitive) detectiontechniques are more robust against the intrinsic uncer-tainties in the underlying physics [71–73].In the case of secular bar-mode GW signals from GRBafterglow plateaus, given the uncertainties related to thephysics of GRB central engines, the derived gravitationalwaveforms are to be considered as simplified phenomeno-logical models. Thus, a more robust (when comparedto matched filtering) search is necessary. A very ro-bust approach against signal uncertainties consists of us-ing the cross-correlation between the output of different,non-colocated detectors. This approach (which, differ-ently from matched filtering, requires no a-priori knowl-edge of the signal waveform and its properties) is typi-cally used for stochastic GW background searches (e.g.,[29, 74, 75]). The cross-correlation between different,non-colocated detectors, only relies on the fact that, inthe presence of a GW signal, the output from distinctdetectors (at the same times, after correcting for thelight-travel time between detectors) should be correlated,while pure noise would remain uncorrelated. Of course,this technique also implies a poor resolution in the pa-rameter space, and more expensive follow-ups to verifypossible detections [56].It is important to note that the cross-correlation is atthe basis of two opposite search strategies: the (highlysensitive) matched filtering (cross-correlation of the datawith a template), and the (very robust) “stochasticsearch” (cross-correlation of different detectors’ output).Indeed, by noticing this fundamental fact, Dhurandhar etal. 2008 [56] have provided an elegant formulation of thecross-correlation statistic for periodic GW searches suchthat, depending on the maximum duration over whichone believes phase coherence is preserved by the signal,the statistic can be tuned to go from a “stochastic-type”search using data from distinct detectors, to the semi-coherent time-frequency methods with increasing coher-ent time baselines (e.g., [62]), and all the way to a fullycoherent search (nearly recovering the matched filteringstatistic).Dhurandar et al.’s formulation of the cross-correlationstatistic [56] leads to a unified framework that can beused to make informed trade-offs between computationalcost, sensitivity, and robustness against signal uncertain-ties. Studies based on the cross-correlation statistic asformulated by [56] have focused on continuous GW emis-sion from Supernova 1987a and Scorpius X-1 [76, 77], anda number of refinements to the cross-correlation methodhave also been published in recent years, particularly forthe treatment of spectral leakage [77, 78]. In what fol-lows, we present a strategy tuned for the detection of intermediate-duration ( (cid:46) s) quasi-periodic GW sig-nals, and discuss its application to the case of secularlyunstable GRB magnetars (Section V). III. NOTATION AND ASSUMPTIONSA. The Short-time Fourier Transform
The Short-time Fourier Transform (SFT) is a usefultool when examining a signal in which frequency con-tent is evolving with time. The time-domain output ofLIGO/Virgo detectors, x ( t ), can be represented as thelinear combination of a GW signal h ( t ), and backgroundnoise n ( t ): x ( t ) = h ( t ) + n ( t ) . (3.1)The SFT of the detector output is constructed by divid-ing the time-series x ( t ) into N SFT segments of duration∆ T SFT (generally speaking, these segments may or maynot overlap), and by taking the Discrete Fourier Trans- form (DFT) of each of these segments:˜ x I [ f k ] = 1 f s N bin − (cid:88) l =0 x [ t l ]e − πif k ( t l − T I +∆ T SFT / , (3.2)where f s is the sampling frequency (typically f s =16 ,
384 Hz for the LIGO detectors); N bin = ∆ T SFT × f s is the number of frequency bins of each SFT; and f k isthe frequency corresponding to the k -th frequency bin: f k = k ∆ T SFT for k = 0 , ..., N bin / − , (3.3) f k = ( k − N bin )∆ T SFT for k = N bin / , ..., N bin − . (3.4)Note that t l in Eq. (3.2) corresponds to the l -th timesample i.e., t l = T I − ∆ T SFT / l/f s . For each I =0 , , ...T obs / ∆ T SFT (where T obs is the total duration ofthe signal) and l = 0 , , ..., N bin , t l spans the time inter-val T I − ∆ T SFT / ≤ t l ≤ T I + ∆ T SFT . Note also thatwe distinguish between continuous time series x ( ... ) andtheir associated discretely-sampled time series x [ ... ] byusing square brackets.To reduce spectral leakage, a windowing function w [ t l ]is often applied to the DFT [79]:˜ x I [ f k ] = N bin − (cid:88) l =0 w [ t l ] x [ t l ]e − πif k ( t l − T I +∆ T SFT / . (3.5)For simplicity, and following [56], hereafter we neglect thewindow function (but discuss some of the related issuesin Section IV D). B. Detector noise and its PSD
In this Section we consider the detector output in theabsence of a signal. In the continuum limit of Eq. (3.1),the frequency ( f ) content of the detector noise can bedescribed by its Fourier transform:˜ n ( f ) = ∞ (cid:90) −∞ dt n ( t )e − πift . (3.6)The single-sided ( f (cid:38)
0) Power Spectral Density (PSD)of the noise, S n ( f ), is defined as: S n ( f ) := 2 ∞ (cid:90) −∞ dτ (cid:104) n ( t ) n ( t + τ ) (cid:105) e − πifτ , (3.7)where (cid:104) n ( t ) n ( t + τ ) (cid:105) is the autocorrelation function ofthe noise, and the expectation value (cid:104)·(cid:105) represents anaverage over an ensemble of noise realizations. The noiseautocorrelation function thus forms a Fourier transformpair with its PSD. Note that hereafter we assume thenoise is stationary and Gaussian (with zero mean), thusits autocorrelation function is independent of t .From Eq. (3.6), it follows that (see also [80]): (cid:104) ˜ n ∗ ( f (cid:48) )˜ n ( f ) (cid:105) = (cid:42) ∞ (cid:90) −∞ dt (cid:48) n ∗ ( t (cid:48) )e πif (cid:48) t (cid:48) ∞ (cid:90) −∞ dt n ( t )e − πift (cid:43) . (3.8)This product of independent integrals can be recast as: (cid:104) ˜ n ∗ ( f (cid:48) )˜ n ( f ) (cid:105) = (cid:42) ∞ (cid:90) −∞ dt (cid:48) ∞ (cid:90) −∞ dt n ∗ ( t (cid:48) ) n ( t )e πif (cid:48) t (cid:48) e − πift (cid:43) . (3.9)Noting that real detector output implies n ∗ ( t ) = n ( t ),and given the linearity and limited multiplicativity of theexpectation value, we have: (cid:104) ˜ n ∗ ( f (cid:48) )˜ n ( f ) (cid:105) = ∞ (cid:90) −∞ dt (cid:48) ∞ (cid:90) −∞ dt (cid:104) n ( t (cid:48) ) n ( t ) (cid:105) e πif (cid:48) t (cid:48) e − πift . (3.10)Setting t = t (cid:48) + τ , yields: (cid:104) ˜ n ∗ ( f (cid:48) )˜ n ( f ) (cid:105) = ∞ (cid:90) −∞ dt (cid:48) e − πi ( f − f (cid:48) ) t (cid:48) ∞ (cid:90) −∞ dτ (cid:104) n ( t (cid:48) ) n ( t (cid:48) + τ ) (cid:105) e − πifτ . (3.11)Then, using Eq. (3.7), we replace the integral over dτ with the PSD, (cid:104) ˜ n ∗ ( f (cid:48) )˜ n ( f ) (cid:105) = S n ( f )2 ∞ (cid:90) −∞ dt (cid:48) e − πi ( f − f (cid:48) ) t (cid:48) . (3.12)The remaining integral over dt (cid:48) is simply a delta function, (cid:104) ˜ n ∗ ( f (cid:48) )˜ n ( f ) (cid:105) = 12 δ ( f − f (cid:48) ) S n ( f ) , (3.13)and using the finite time approximation of the delta func-tion: δ ∆ T SFT ( f ) = sin( πf ∆ T SFT ) πf , (3.14)which reduces to ∆ T SFT in the limit of f →
0, we canrelate the variance of the Fourier transformed detectoroutput to the PSD: (cid:104)| ˜ n I [ f k ] | (cid:105) ≈ ∆ T SFT S n [ f k ] . (3.15) The expectation value (cid:104) XY (cid:105) of random variables X, Y is multi-plicative if Cov(
X, Y ) = 0. That is, only if X and Y are statis-tically independent. C. Short-duration Fourier Transform of the signal
We make the hypothesis that the GW signal h ( t ) isquasi-periodic (by taking a sufficiently small time in-terval the signal in such an interval can be consideredmonochromatic), and assume that its time-frequencyevolution is described with sufficient physical accuracy,for a time interval of T coh , via some known function ofa given set of parameters (although this function maynot have a closed form expression). By definition, this“coherence timescale” is less than or equal to the totalobservation time T obs over which the signal is expectedto last (e.g. T coh (cid:46) T obs (cid:46) s for the type of signals ofinterest in the context of GRB afterglow plateaus).Since the signal is quasi-periodic, we can define an SFTbaseline ∆ T SFT ≤ T coh such that, within the baseline,all of the signal power is concentrated in a single SFTbin. More specifically, around each time T I we can ap-proximate the signal received by the detector in the timeinterval T I − ∆ T SFT (cid:46) t (cid:46) T I + ∆ T SFT , as: h ( t ) ≈ h ( T I ) A + F + cos(Φ( T I ) + 2 πf ( T I )( t − T I ))+ h ( T I ) A × F × sin(Φ( T I ) + 2 πf ( T I )( t − T I )) , (3.16)where A + , A × are amplitude factors dependent on thephysical system’s inclination angle ι (for on-axis GRBs, ι is the angle between the jet axis and the line of sight): A + = 1 + cos ι , (3.17) A × = cos ι, (3.18)and F + , F × are the antenna factors that quantify the de-tector’s sensitivity to each polarization state. Note thatfor triggered searches targeting GRBs (as is the case inSec. V), the line of sight is expected to be nearly alignedwith the jet axis, thus ι ≈ A + ≈ A × ≈ T obs (cid:46) s so that, for a given GW detector, F + and F × can be treated as constants as a functionof time (see e.g. [56]).2. If ˙ f ( t ) is the time derivative of the signal fre-quency at a given time t , then the effects of ˙ f ( t )on the signal phase should be negligible during thetime interval ∆ T SFT . Using the quarter-cycle crite-rion, this leads to 2 π | ˙ f ( T I ) | (cid:0) ∆ T SFT (cid:1) < π . Thus,∆ T SFT < / (cid:113) | ˙ f ( T I ) | .3. ∆ T SFT is small enough that h ( t ) ≈ h ( T I ) (con-stant amplitude approximation) in the interval That is, the line of sight is within the jet-opening angle, whichis expected to be of the order 5 −
20 deg for long GRBs [81, 82]. T I − ∆ T SFT / (cid:46) t (cid:46) T I + ∆ T SFT /
2. We considerthis condition satisfied if (cid:12)(cid:12)(cid:12) ˙ h ( T I ) (cid:12)(cid:12)(cid:12) ∆ T SFT /h ( T I ) (cid:46) ∼
10% [83]; thus, any change of signal am-plitude below 10% is not expected to significantlyaffect the goodness of this approximation).In addition, hereafter we assume that ∆ T SFT is largeenough that the corresponding frequency resolution,(∆ T SFT ) − , still enables one to track the time-frequencyevolution of the signal.Using Eq. (3.2), we can calculate the DFT of the signalin Eq. (3.16) (see also Eq. (2.25) in [56]):˜ h I [ f k ] = h ( T I ) e iπf k,I ∆ T SFT × [ e i Φ( T I ) A + F + ,I − i A × F × ,I δ ∆ T SFT ( f k − f k,I )+ e − i Φ( T I ) A + F + ,I + i A × F × ,I δ ∆ T SFT ( f k + f k,I )] , (3.19)or, equivalently,˜ h I [ f k ] = (cid:113) A F ,I + A × F × ,I h ( T I ) e iπf k,I ∆ T SFT × (cid:104) e i Φ( T I ) e iϕ I δ ∆ T SFT ( f k − f k,I )+ e − i Φ( T I ) e − iϕ I δ ∆ T SFT ( f k + f k,I ) (cid:105) , (3.20)where we have set: A + F + ,I ± i A × F × ,I = (cid:113) A F ,I + A × F × ,I e ∓ iϕ I , (3.21)and ϕ I = arctan( −A × F × ,I / A + F + ,I ) . (3.22)Note that, while in our limit of intermediate-durationGW signals the antenna response from one detector canbe considered constant over the observed duration of thesignal, for the multiple detector case the antenna re-sponses refer to the specific GW detector from whoseoutput the I -th SFT is taken. IV. THE CROSS-CORRELATION STATISTIC
Following [56], we define the raw cross-correlationstatistic as: Y IJ = ˜ x ∗ I [ f k,I ]˜ x J [ f k (cid:48) ,J ]∆ T SF T , (4.1)where the frequency f k,I is the frequency at which allof the signal power is concentrated during the I th timeinterval (see Eq. (3.19)), and is related to the frequency f k (cid:48) ,J at which all of the signal power is concentrated dur-ing the J th time interval via the relation: f k (cid:48) ,J = f k,I − ∆ f IJ . (4.2)In the above relation, ∆ f IJ is the frequency differencepredicted by the model’s time-frequency evolution (inthis analysis the signal time-frequency evolution is as-sumed to be known to some level of accuracy; see SectionIII C). Note that, because for any I -th SFT the associ-ated frequency bin k is fixed by the model’s predictions,we omit the indexes k, k (cid:48) from Y IJ for simplicity.For a signal embedded in stationary Gaussian noisewith zero mean, the {Y IJ } are themselves random vari-ables with mean and variance given by µ IJ = h ( T I ) h ( T J ) ˜ G IJ , (4.3) σ IJ = 14∆ T S n [ f k,I ] S n [ f k (cid:48) ,J ] , (4.4)where we have used Eqs. (3.15) and (3.19), and the factthat:˜ h ∗ I [ f k ]˜ h J [ f k +∆ f IJ ] = h ( T I ) h ( T J ) ˜ G IJ δ T SFT ( f k − f k,I ) . (4.5)In the above equations, ˜ G IJ is the signal cross-correlationfunction, defined here as˜ G IJ = (cid:113) A F ,I + A × F × ,I (cid:113) A F ,J + A × F × ,J e − i ∆ θ IJ , (4.6)with ∆ θ IJ = θ I − θ J = π ∆ f IJ ∆ T SFT + ∆Φ IJ + ∆ ϕ IJ . Ingeneral, the subscripts ( I ) , ( J ) in the antenna responsesrefer to the specific GW detector from whose output the I -th (or J -th) SFT is taken. Indeed, in the definitionof the {Y IJ } , there is total freedom to correlate pairsfrom one single detector or from an arbitrary number ofdetectors.Note that the e − iπ ∆ f IJ ∆ T SFT term that arises from∆ θ IJ in Eq. (4.6) is absent from the definition of thesignal-cross-correlation function given in [56]. This dis-crepancy was first noted in [76], and is discussed therein detail. This term proves essential to properly trackingthe frequency evolution of a given signal across SFTs, sowe call attention to it here.When cross-correlation pairs are only taken from theoutput of a single detector over timescales of T obs (cid:46) s,then F + , × ,I = F + , × ,J = F + , × . This simplifies Eq. (4.6)considerably:˜ G IJ = A F + A × F × e − iπ ∆ f IJ ∆ T SFT e − i ∆Φ IJ . (4.7)For two or more detectors, such as LIGO Hanford (H)and LIGO Livingston (L), the indexes I and J are free torange over SFTs from either detector, and so the abovesimplification does not generally apply (even if the an-tenna factors for each detector are approximately con-stant within the considered time interval).Following [56], our detection statistic is then con-structed as a weighted sum of the Y IJ ρ = (cid:88) IJ ( u IJ Y IJ + u ∗ IJ Y ∗ IJ ) , (4.8)with nearly optimal weights u IJ = ˜ G ∗ IJ σ IJ . (4.9)For stationary Gaussian distributed white noise (seeEq. 4.4), σ IJ does not depend on frequency nor on time,but it might still depend on the detector. Thus: σ IJ = 14∆ T S n , (4.10)for IJ pairs from a single detector (or identical detec-tors), or: σ IJ = 14∆ T S Hn S Ln , (4.11)for e.g. a LIGO Hanford-Livingston IJ pair. Thus, usingthe above equations and Eq. (4.6), we have in general: u IJ = (cid:113) ( A F ,I + A × F × ,I )( A F ,J + A × F × ,J )∆ T − e − i ∆ θ IJ S n [ f k,I ] S n [ f k,J ] , (4.12)where, again, the antenna responses and detector’s noiserefer to the specific GW detector from whose output the I -th (or J -th) SFT is taken.As we describe in more detail in what follows, themean and variance of ρ , as well its statistical distribu-tion, depend on the choice of which SFT pairs are cross-correlated. Because of the freedom in choosing whichdata-segment pairs to correlate, we can naturally con-sider one single detector or an arbitrary number of de-tectors (with no need to modify our statistic), and wecan work in one of the following limits [56]:1. We can choose to correlate only data segmentstaken from distinct detectors at the same times (af-ter correcting for the light travel time between dif-ferent detectors; Section IV A). This limit is anal-ogous in spirit to the methods of stochastic GWsearches, such as [84–87], and we hence refer to itas the “stochastic limit”. In this case, the compu-tational cost of the search is small and the searchis very robust against signal uncertainties. But thesensitivity is the poorest, as is the resolution in pa-rameter space. Strictly speaking, these weights are only optimal when self-pairsare excluded, as in [56]. For sufficiently small amplitude signals,these weights remain optimal, to first order, even when self-pairsare considered. For situations where this may not be the case,we refer the reader to the discussion in the Appendix of [56].
2. At the other extreme, we can correlate all possibleSFT segments (Section IV B): This (nearly) corre-sponds to a full matched filter statistic described forcoalescing compact binaries and continuous wavesin e.g. [60, 62, 88, 89]. The parameter space resolu-tion becomes very fine and while this is ideally themost sensitive method, is it also the most computa-tionally expensive (prohibitive for wide parameterspace searches) and the least robust against signaluncertainties.3. In intermediate regimes, we can correlate data seg-ments separated by a maximum coherence time T coh (cid:46) T obs (Section IV C). This “semi-coherent”approach is similar to several methods used forcontinuous waves [71, 90–93] (though on signaltimescales much longer than what considered inthis work). Because in this limit the sensitivityand robustness of the search can be tuned to the ex-pected accuracy of a given model, this is the regimeof greatest interest to us.4. Finally, one can consider all pairs except self-correlations. This was the main focus of the anal-ysis presented in [56] (see their Section IV). Here,we do not focus on this limit because we considerit a special case of the ones above (with no par-ticular advantages for the detection of the type ofsignals considered in our study and with some com-plications added to the statistical properties of ρ ).However, in what follows, we do discuss the maindifferences of (1)-(3) above with respect to this case(see also Section IV of [56]).In discussing the above limits, it is useful to note thatwe can re-write Eq. (4.8) in terms of Eq. (4.1) as: ρ = 1∆ T (cid:88) IJ u IJ ˜ x ∗ I [ f k,I ]˜ x J [ f k (cid:48) ,J ] + u ∗ IJ ˜ x I [ f k,I ]˜ x ∗ J [ f k (cid:48) ,J ] , (4.13)which is equivalent to: ρ = 2∆ T (cid:88) IJ (cid:60) { u IJ ˜ x ∗ I [ f k,I ]˜ x J [ f k (cid:48) ,J ] } . (4.14) A. Stochastic limit (independent pairs only)
Consider the output of two different detectors, ˜ x H and ˜ x L . Each detector’s output can be divided into T obs / ∆ T SFT = N SFT segments. Of the (2 N SFT ) possi-ble SFT pairs that can contribute to ρ we correlate onlypairs of SFTs from different detectors at the same time(after correcting for the light travel time between detec-tors), so that N pairs = N SFT . In this limit, Eq. (4.13)becomes: ρ = 2 (cid:88) I (cid:60) { u II Y II } , (4.15)where the weights are described by e.g. Eq. (4.11). Writ-ten explicitly, this becomes ρ = 2∆ T (cid:88) I (cid:60) (cid:8) u II ˜ x ∗ H I [ f k,I ]˜ x L I [ f k (cid:48) ,I ] (cid:9) , (4.16)i.e., a weighted sum of completely independent randomvariables that are each the product of two Gaussian vari- ables with mean and variance given by Eqs. (4.3) and(4.4). Thus, ρ converges to a Gaussian distribution (bythe Central Limit Theorem) with mean (see Eqs. (4.3),(4.6), (4.12), and [56]) and variance (see also Eq. (4.4)and Dhurandhar et al. [56]): µ ρ = ( A F ,H + A × F × ,H )( A F ,L + A × F × ,L ) ∆ T (cid:88) I h ( T I ) S Hn [ f k,I ] S Ln [ f k,I ] , (4.17) σ ρ = ( A F ,H + A × F × ,H )( A F ,L + A × F × ,L ) ∆ T (cid:88) I S Hn [ f k,I ] S Ln [ f k,I ] . (4.18)The detection threshold is easily derived in terms ofthe Cumulative Distribution Function (CDF) of a normaldistribution, F N ( ρ ) = 12 (cid:34) − erfc (cid:32) ρ − µ ρ σ ρ √ (cid:33)(cid:35) , (4.19)and its inverse (see also [56]), where erfc is the comple-mentary error function. For a False Alarm Probability(FAP) α , the associated threshold is simply 1 − α = F N ( ρ th ), thus: ρ th = √ σ ρ erfc − (2 α ) , (4.20)where we have used the fact that the background distri-bution is considered in the absence of a signal ( µ ρ = 0).When a signal is present, the detection probability γ , or,equivalently, the False Dismissal Probability (FDP) 1 − γ ,is given by γ = F N ( ρ th ), i.e.: γ = 12 erfc (cid:32) ρ th − µ ρ σ ρ √ (cid:33) . (4.21)Thus, the detectability condition reads: µ ρ σ ρ (cid:38) √ S , (4.22)where S = erfc − (2 α ) − erfc − (2 γ ). In the case of whiteGaussian noise, using Eqs. (4.17)–(4.18), the detectabil-ity condition implies: h rms (cid:38) √ S / ∆ T − / N − / ( S Hn S Ln ) / [( A F ,H + A × F × ,H )( A F ,L + A × F × ,L )] / , (4.23)which generalizes Eq. (4.15) in Dhurandhar et al. [56]to the case of a non-constant signal amplitude for which(see also Eq. (3.16)): h rms = (cid:113) (cid:104) h ( T I ) (cid:105) I = (cid:115) (cid:80) I h ( T I ) N SFT . (4.24) › ρ fi N o r m a li z e d p r o b a b ili t y
1e 48
Signal (Stochastic Limit)
Predicted NoisePredicted SignalSimulated NoiseSimulated Signal
FIG. 1. Comparison between the simulated and predicteddistribution of ρ in the stochastic limit, for 2048 s of simulatedwhite Gaussian noise sampled at a rate of f s = 2048 Hz, fromtwo detector’s outputs x H [ t ] , x L [ t ]. We have used an SFTbaseline of ∆ T SFT = 2 s and, for simplicity, we assumed twoidealized, colocated, and optimally oriented detectors withaLIGO-equivalent PSDs S n ≈ . × − Hz − , see Fig.2. The simulated signal is a line of constant frequency f =128 Hz and constant amplitude h ≈ − . In Fig. 1 we show the distribution of ρ in the absenceof a signal for simulated Gaussian white noise, and inthe presence of a GW signal of constant amplitude h and constant frequency f . (A signal with constant fre-quency represents the simplest time-frequency evolutionto which the technique here presented can be applied andis particularly useful for illustrative purposes.)We stress that the independence of the pairs that areadded in ρ is essential for the validity of the conclusionregarding the Gaussianity of ρ , and for the validity ofEqs. (4.23)-(4.24). While pairs are truly independent inthe stochastic limit analyzed in this Section, this is notstrictly true for the combination of pairs considered in Frequency (Hz) -48 -47 -46 -45 -44 -43 P o w e r Sp e c t r a l D e n s i t y ( H z − ) Gaussian Noise Approximation aLIGO NoiseGaussian ApproximationSignal Frequency Range
FIG. 2. Method for generating Gaussian-distributed noisewith aLIGO PSD. A frequency range is defined with respectto the maximum and minimum value of an injected signal’sfrequency (blue dotted lines). The constant Gaussian PSD S n (red dashed line) is calculated so that the area underneath it(shaded red) is equal to the area underneath the aLIGO PSD(shaded blue). For signals of constant frequency f (e.g. Figs.1, 3, and 4), the red area is taken between 0 . f and 1 . f . Section IV of [56] ( ρ includes all possible pairs but selfones - see also case 4 in Section IV) and in the Appendixof [56] ( ρ includes all possible SFT pairs - see also case2 in Section IV). In these cases, ρ is a sum of productsthat are not all independent. Thus, while the expres-sions for the mean and variance of ρ presented in SectionIV of [56] (or, equivalently, Eqs. (4.17) and (4.18) here)remain valid, we caution the reader that the lack of inde-pendence affects the shape of the background distribution ,and in some limits, results in a distribution that cannot be reduced to a Gaussian. Thus, the detection thresholdneeds to be modified accordingly. Some brief discussionof these corrections to [56] is also presented in AppendixB of [77]. In what follows, our in depth discussion ofcases 2–3 (Section IV) shows explicitly that the correc-tions to [56] are crucial for the detection of the family ofintermediate-duration GW signals that we target in thisanalysis. B. Matched filter limit (all pairs)
In this limit, we choose to correlated all possible SFTsegments (from one or multiple detectors). Starting fromEq. (4.14), we replace the weights with their explicit formgiven by Eq. (4.12), ρ = 2 (cid:60) N pairs (cid:88) I,J (cid:113) ( A F ,I + A × F × ,I ) S n [ f k,I ] (cid:113) ( A F ,J + A × F × ,J ) S n [ f k,J ] ˜ x ∗ I [ f k,I ]˜ x J [ f k,J ] e i ∆ θ IJ , (4.25)where N pairs = N and N SFT = N det T obs / ∆ T SFT , with N det being the number of detectors from which data aretaken. Under the change of variable˜ x (cid:48) I [ f k,I ] = (cid:113) ( A F ,I + A × F × ,I ) S n [ f k,I ] ˜ x I [ f k,I ] e − iθ I , (4.26)Eq. (4.25) simplifies to: ρ = 2 (cid:40) N SFT (cid:88) I | ˜ x (cid:48) I [ f k,I ] | + 2 N SFT (cid:88)
I>J (cid:60) [˜ x (cid:48)∗ I [ f k,I ]˜ x (cid:48) J [ f k (cid:48) ,J ]] (cid:41) . (4.27)It then follows that, ρ = 2 (cid:12)(cid:12)(cid:12)(cid:12)(cid:12) N SFT (cid:88) I ˜ x (cid:48) I [ f k,I ] (cid:12)(cid:12)(cid:12)(cid:12)(cid:12) . (4.28)Or alternatively, ρ = 2 (cid:32) N SFT (cid:88) I (cid:60) (˜ x (cid:48) I [ f k,I ]) (cid:33) + (cid:32) N SFT (cid:88) I (cid:61) (˜ x (cid:48) I [ f k,I ]) (cid:33) . (4.29) For stationary Gaussian noise with zero mean, ˜ x I [ f k ]follows a complex normal distribution. We note that thescaling and complex rotation applied in Eq. (4.26) haveno effect on the shape of the distribution of the ˜ x (cid:48) I whencompared to the ˜ x I (but they do change the mean andvariance of the distribution). Thus, the real and imagi-nary parts of ˜ x (cid:48) I are still Gaussian distributed as the ˜ x I ,and so are their sums. Indeed, in the absence of a signal,the sums of the real and imaginary parts of the ˜ x (cid:48) I areGaussian variables with zero mean and variance (see Eq.(3.15)): σ = N SFT (cid:88) I (cid:34) ∆ T SFT ( A F ,I + A × F × ,I )4 S n [ f k,I ] (cid:35) . (4.30)So we can re-write the expression for ρ as: ρ = C χ × (cid:32) (cid:80) N SFT I (cid:60) (˜ x (cid:48) I [ f k,I ]) σ Σ (cid:33) + (cid:32) (cid:80) N SFT I (cid:61) (˜ x (cid:48) I [ f k,I ]) σ Σ (cid:33) , (4.31) › ρ fi N o r m a li z e d p r o b a b ili t y
1e 50
Matched Filter Limit Distributions
Predicted SignalPredicted NoiseSimulated SignalSimulated Noise
FIG. 3. We simulate 2048 s of white Gaussian noise for a sin-gle optimally-oriented detector with aLIGO-equivalent noisePSD given by S n ≈ . × − Hz − . We used a sam-pling frequency of f s = 2048 Hz and an SFT baseline of∆ T SFT = 2 s. All possible pairs are included in the cross-correlation statistic ρ which is thus distributed as a scaled χ distribution with 2 degrees of freedom. The simulated signalwas a line of constant frequency f = 128 Hz and constantamplitude h ≈ . × − . which is the sum of the squares of two normally dis-tributed variables, scaled by a factor: C χ = 2 σ = N SFT (cid:88) I (cid:34) ∆ T SFT ( A F ,I + A × F × ,I )2 S n [ f k,I ] (cid:35) . (4.32)Thus, the resulting ρ statistic is distributed as a χ with2 degrees of freedom (Fig. 3; see also [77]).Continuing from Eq. (4.31), in the absence of a signal,the variance of ρ is simply, σ ρ = 4 C χ = 2 N SFT (cid:88) I (cid:34) ∆ T SFT ( A F ,I + A × F × ,I ) S n [ f k,I ] (cid:35) . (4.33)In the presence of a signal, the distribution of ρ inEq. (4.31) becomes a non-central χ with two degrees offreedom, χ nc (2; λ ), of mean: µ ρ = C χ (2 + λ ) . (4.34)The non-centrality parameter can be derived using theabove relation, and noting that µ ρ can be easily calcu-lated using Eqs. (4.5), (4.6), and (4.25). This yields (seealso Eq. (4.24) and Fig. 3): λ = N SFT (cid:88) I h ( T I ) (cid:34) ∆ T SFT ( A F ,I + A × F × ,I ) S n [ f k,I ] (cid:35) . (4.35)Note that in this limit the number of SFT pairs onlyaffects the variance (and mean) of the two Gaussian variables (cid:80) N SFT I (cid:60) (˜ x (cid:48) I [ f k,I ]) and (cid:80) N SFT I (cid:61) (˜ x (cid:48) I [ f k,I ]). Itdoes not affect the number of degrees of freedom in ρ ,which remains two independently of the number of SFTs.Thus, as N SFT increases, the distribution of ρ does not approach a Gaussian. This is a critical distinction tomake, since it changes the (false alarm and false dis-missal) thresholds of ρ significantly from the ones thatwere adopted in the appendix of [56], where a Gaussiandistribution was incorrectly assumed for ρ .In the case in which all pairs come from a single detec-tor (or from colocated, equally oriented detectors, withidentical S n ), the variance of ρ simplifies substantially to: σ ρ = 4 C χ = 2 T obs (cid:20) ( A F + A × F × ) S n (cid:21) , (4.36)where we have used T obs = N SFT ∆ T SFT . The non-centrality parameter likewise simplifies, yielding, λ = h T obs ( A F + A × F × ) S n , (4.37)where we have used Eq. (4.24).In either case, the corresponding detection thresholdfor a given false alarm and detection rate is now sub-stantially different than in the stochastic limit: ρ th = C χ F − χ (1 − α ; 2) , (4.38) γ = F nc χ ( ρ th /C χ ; 2 , λ ) . (4.39)The CDF for the χ (2) is known in closed form (and iseven invertible), while the CDF for the non-central casecan be calculated numerically, with results as shown inFig. 5.In this limit, the sensitivity approaches that ofmatched filtering. However there is one significant er-ror in the description in [56]: the limit approached isthat of filtering with an unknown overall phase constant,which is commonly handled by summing the squares oftwo matched filters a quarter cycle out of phase witheach other—e.g., [61]. Hence the resulting statistic is dis-tributed as a χ with 2 degrees of freedom rather thana Gaussian. Under idealized circumstances, this reducesthe sensitivity by approximately 13% with respect to aGaussian distribution (with FAP=0.1% and FDP=50%). C. Semi-coherent regime
As discussed in Section II, the semi-coherent regime isthe most relevant for an astrophysically motivated searchwhere the expected GW signal is known to limited accu-racy. In this regime, the total observation time T obs isbroken up into N coh coherent segments, each of dura-tion T coh . The coherence time ( T coh ) is once again de-fined as the length of time wherein the signal is expectedto maintain phase coherence (and therefore good agree-ment) with the model predictions. All possible SFT-pairswithin each coherent time segment are cross-correlated,0 › ρ fi N o r m a li z e d p r o b a b ili t y
1e 50
Matched Filter Limit Distributions
Predicted SignalPredicted NoiseSimulated SignalSimulated Noise
FIG. 4. Comparison between the simulated and predicted dis-tribution of ρ in the semi-coherent limit, for 1024 s of simu-lated white Gaussian noise sampled at a rate of f s = 2048 Hz,from one detector’s output x ( t ). We have used an SFT base-line of ∆ T SFT = 2 s and we assumed an optimally orienteddetector with PSD S n ≈ . × − Hz − . The coherencetime is T coh = 256 s for a total of N coh = 4 coherent seg-ments. The simulated signal was a line of constant frequency f = 128 Hz and constant amplitude h ≈ . × − . and the results for each segment are then combined in-coherently. In order for the resulting sum of χ (2) distributed vari-ables to add to a χ (2 N coh ) distributed detection statis-tic, it is essential that all coherent segments have identi-cal scale factors. This condition is satisfied for a detec-tor network of arbitrary size only if the detectors havesimilar antenna factors for the given sky location of theevent, and each detector has (stationary) white Gaussiannoise (although the frequency independent S n of the de-tectors need not be identical). In the case of colorednoise, the scale factors will vary between coherence seg-ments (since the frequency of the signal is evolving withtime, which causes S n [ f k,I ] to change from segment tosegment). Thus, in the presence of colored noise, whiten-ing the data over the signal bandwidth prior to analysisis desirable.Changes in the antenna factors F + , F × over the du-ration of a signal in a non-idealized search i.e., devia-tions from assumption 1 in Section III C, can also af-fect the statistic. For the GRB X-ray plateaus of inter- An alternative, but equivalent, description is to define a “coher-ence window” of duration T coh which is then stepped across theSFT according to a given spacing criterion. All segments in eachstep are cross-correlated then combined incoherently. In general, for random χ variables X i , their linear combina-tion Y = (cid:80) i C i X i is itself a χ variable if and only if the scalecoefficients C i are identical (or 0). However, if the normalized co-efficients C i / (cid:104)C i (cid:105) are close to unity, Y is reasonably approximatedby a χ distribution. est to Section V, >
50% of events with sufficiently shal-low plateau decays have plateau durations (cid:46) s [94].For circularly-polarized signals of this duration, we ten-tatively estimate that time-varying antenna factors willcause fluctuations of ≈
15% in amplitude sensitivity, com-parable to LIGO amplitude calibration uncertainties [83].We leave to future work a more in depth examination ofdeviations from this assumption.When all coherent segments have identical scale fac-tors, ρ is an incoherent sum of N coh independent vari-ables, each distributed as a scaled χ (2) distribution.The scale parameter for each coherent segment is givenby Eq. (4.32) but now with N SFT = T coh / ∆ T SFT , sothat: C SC χ = C χ N coh . (4.40)The variance of the semi-coherent ρ then reads: σ ρ, SC = 2 C SC χ (2 N coh ) = 4 C χ , (4.41)which is identical to the variance in matched-filter limit,see Eq. (4.33).When a signal is present, the non-centrality parameterfor each coherent segment will, in general, vary from onesemi-coherent chunk to the other due to the time-varyingamplitude of the signal in Eq. (4.35). But, since thetotal λ for the semi-coherent regime is additive across all coherent segments, the total non-centrality parameter forthe semi-coherent ρ likewise remains unchanged from thematched filter limit. The resulting distribution thus hasmean: µ ρ, SC = C SC χ (2 N coh + λ ) . (4.42)The above Equation reduces to (4.34) in the limit of N coh = 1 (matched-filter limit). The detection thresh-old for a given signal will differ from the matched filterlimit due to the higher number of degrees of freedom ofthe χ distribution of the semi-coherent ρ , and can becalculated numerically as shown in Fig. 5 (along withother limits).In the limit of large N coh , the χ (2 N coh ) distributiontends toward a Gaussian. Continuous wave searches us-ing this cross-correlation technique, e.g. [77] can have N coh of order 10 , and hence can set their thresh-olds based on Gaussian statistics as assumed by [56].However, searches for intermediate-duration GW signals(such as those of interest to this paper) can have N coh smaller by 1–2 orders of magnitude, so it is essential tocorrect the corresponding detection thresholds to accountfor non-Gaussianity. In particular, the Central LimitTheorem reduces the skew of the χ (2 N coh ) distributionrelatively slowly as N coh grows. We consider specifically Type IIa GRBs as described in [94] withplateau power law decay indexes of magnitudes (cid:46) -5 -4 -3 -2 -1 False Alarm Probability h m i n
1e 24
Min Amplitude vs FAP (FDR=50%)
X-Cor ( N coh = 1 )X-Cor ( N coh = 64 )X-Cor ( N coh = N SFT = 512 ) X-Cor (Stoch)DhurandharMatched Filter (shaded)
FIG. 5. The smallest detectable GW amplitude h min is plot-ted versus FAP with a set FDP of 1 − γ = 50%. A matchedfilter with known initial phase (black dotted line, with grayshading) is the idealized optimal search, and it provides anabsolute limit on the sensitivity of any real search. Hence,the shaded gray area is forbidden. The matched-filter limitof the cross-correlation method (dashed blue line) is expectedto approach (but not converge with) the black-dotted line.The semi-coherent limit (dash/dotted purple lines) becomesless sensitive for increasing N coh , eventually approaching thestochastic limit (dotted red). The N coh = N SFT = 512 limitof the cross-correlation method (red dash/double-dot red line)differs from the stochastic limit in that it includes self-pairs(autocorrelations). As discussed in the text, the assumptionsof Gaussian statistics and known phase constant in [56] yieldincorrect results as the resulting sensitivity (green solid line)does better than the optimal matched filter for sufficientlysmall FAP.
D. Spectral Leakage Effects
Several of the assumptions made in the previous Sec-tions are expected to lead to some amount of spectralleakage. These include the finite-time approximation ofthe delta function in Eq. (3.14), the quarter-cycle cri-terion, and SFT windowing effects (that is, the simpli-fication of using a simple rectangular window). A fulltreatment of the effects of spectral leakage is outside thescope of this paper but we mention some of its effectshere.As shown in Fig. 6, spectral leakage is an issue anytime the signal frequency does not precisely correspondto the center of one of the SFT frequency bins. In the simplest case of a constant frequency periodic sig-nal, spectral leakage can cause a reduction of up to 50%in the SNR ( µ ρ, signal /σ ρ, noise ) for the ρ statistic in eachof the fully coherent segments. This effect is worsenedwhen one considers time-varying frequencies: while thequarter-cycle criterion restricts the leakage from first or-der terms ( ˙ f ), higher order components of the frequencyevolution ( ¨ f , ... f , etc) can lead to additional leakage. Thenet result is that, on average, neglecting spectral leakagewill result in reduced SNR that is roughly 75% of theidealized case, see Figures 6a and 6b, and also [78].The typical solution for this problem is to introducea windowing function for the SFT, but this is not with-out tradeoffs. Each windowing function (of which thereare many) has different strengths and weaknesses. Thecommonly used Hann window is well equipped to handlespectral leakage and maintains good frequency resolu-tion, but suffers in amplitude accuracy [78]. SFT win-dows must then be overlapped in an attempt to regainsome of the lost amplitude information, increasing com-putational cost. The Tukey window, commonly used incontinuous wave searches, is – by contrast – not as goodat diminishing the effects of spectral leakage but retainsmore of the original power. Recent work within the cross-correlation framework has examined the effects of differ-ent windowing functions in detail [77, 78].Other methods can also be used to reduce spectralleakage. These include over-resolving each SFT by zero-padding (although this can still lead to some spectralleakage for signals in which frequency varies continuouslywith time), sinc-interpolating between SFT bins (thusleveraging the sampling theorem [95]), or simply addingcontributions from neighboring SFT bins. Including justthe two adjacent SFT bins when cross-correlating can im-prove recovery of the expected SNR from ≈
77% to ≈ ≈
75% of theidealized value for the ρ statistic (i.e. up to a factor of √ . ≈
87% in signal amplitude and/or distance reachfor cases in which ¨ f and higher terms may not be neg-ligible). This is consistent with the estimate of 77 . f via the quarter-cycle criterion (assumption √ . ≈
71% in amplitude sensitivity, see Fig. 6a).
V. GRB PLATEAU SEARCH SENSITIVITY
In this Section we apply the cross-correlation statisticto the specific model of intermediate-duration GW sig-nals described in [38]. This model describes the scenarioof a secularly unstable GRB-magnetar possibly associ-ated with a GRB afterglow plateau (see also Section I).2 f (Hz) S N R SNR variation w/line injection
Expected SNRRecovered SNR (a) SNR variation with respect to the location of f . The SNR for aline ( f ( t ) = f ) is at a maximum when f lies at the center of one ofthe bins, and at a minimum when it lies on an edge. In the case ofthe latter, the SNR is reduced by a factor of 2, as the leakage leadsto approximately half of the signal power leaking into each adjacentbin. Thus, the maximum loss in amplitude sensitivity from spectralleakage should be no more than √ ≈
71% for a signal satisfyingthe quarter-cycle criterion (see Sec. III C). The grey shaded regioncorresponds to the 3 σ error region resulting from 512 independentsimulations. f (Hz) S N R SNR variation w/pulsar-like injection
Expected SNRRecovered SNR (b) SNR variation with respect to the location of f for apulsar-like evolution, f ( t ) = f + f t , where f = − / / ∆ T SFT
Hz) over the entire duration of the signal ( T obs ), wellwithin the quarter-cycle criterion. Here, the factor of f guaranteesthat the signal never precisely corresponds to the center of the bin,which averages out to roughly 75% of the SNR (or √ ≈
87% inamplitude) for any given f The grey shaded region corresponds tothe 3 σ error region resulting from 512 independent simulations. FIG. 6. The effects of spectral leakage on signals of the form h ( t ) = h sin Φ( t ) with Φ( t ) = 2 π (cid:82) f ( t ) dt and h = 10 − injectedinto Gaussian noise with zero mean and S n ≈ . × − Hz − . Two frequency evolutions are considered: a line feature ofconstant frequency f ( t ) = f (left); a pulsar-like evolution of the form f ( t ) = f − f t (here f = 1 / / s, right). The SFTbaseline used to calculate the SNR is ∆ T SFT = 2 s, resulting in SFT bin widths of 1 / f = 128 . .
25 and 127 .
75, red vertical dashed lines). In all cases the total duration of the signal is T obs = 512 s. In the Newtonian limit, the l = m = 2 f -mode be-comes secularly unstable when the ratio β = T / | W | ofthe rotational kinetic energy T to the gravitational bind-ing energy | W | is between 0 .
14 and 0 .
27. This modehas the shortest growth time of all polar fluid modes,1 s (cid:46) τ GW (cid:46) × s for 0 . (cid:38) β (cid:38) .
15 [54] and maybe an important source of GWs. Under the hypothesisthat a secular bar-mode instability does indeed set in fora magnetar left over after a GRB explosion, Corsi andM´esz´aros [38] have followed the NS quasi-static evolu-tion under the effect of gravitational radiation accordingto the analytical formulation given by [54]. Since τ GW is generally much longer than the dynamical time of thestar, the evolution is quasi-static, i.e., the star evolvesalong an equilibrium sequence of Riemann-S ellipsoids.Differently from what was done by Lai and Shapiro [54],Corsi and M´esz´aros [38] added into the evolution energylosses due to magnetic dipole radiation, assuming thatthose will not substantially modify the dynamics, butwill act to speed up the overall evolution along the samesequence of Riemann-S ellipsoids that the NS would havefollowed in the absence of radiative losses.In the model proposed by Corsi & M´esz´aros 2009 [38],the resulting quasi-periodic GW signal depends on fiveparameters: β , the initial kinetic-to-gravitational po-tential energy ratio of the magnetized NS [54]; n , the NS polytropic index; M , the NS mass; R , the unper-turbed NS radius; and B , the initial dipolar magneticfield strength at poles. For a typical parameter choiceof M = 1 . (cid:12) , R = 20 km, n = 1, B = 10 G, and β = 0 .
20 (Fig. 7, red), [38] have estimated a distancereach (assuming a matched filter search) of ≈
100 Mpcfor the aLIGO-Virgo detectors (for FAP ≈ × − andFDP=50%).We have tested the detectability of this class of sig-nals using the adaptation of the cross-correlation statisticdescribed in the previous Sections, and assuming Gaus-sian noise with S n ( f ) ≈ . × − Hz − approximatelyequal to that of whitened aLIGO-Virgo noise in the sig-nal’s frequency band (Fig. 2). The results are reportedin Table I for an optimally oriented GRB .A matched filter analysis yields the highest sensitivity,and thus the largest horizon distance limits. For a typicalchoice of model parameters (e.g. β = 0 . B = 10 G, M = 1 . M (cid:12) , R = 20 km), if we assume that the initial Here, “optimally oriented” is taken to mean that the GRB jetis aligned with the line of sight (so that ι = 0 and the GW iscircularly polarized, i.e. A + = A × = 1, see Sec. III C) and theGRB sky location is such that the line of sight is orthogonal tothe plane containing the detector (so that F + F × = 1). -2 -1 Time (s) F r e q u e n c y ( H z ) Magnetar Frequency Evolutions β =0 . β =0 . FIG. 7. Frequency evolution for two representative signalsgenerated via the Corsi and M´esz´aros model. The signals arefor a typical choice of parameters M = 1 . (cid:12) , R = 20 km, n = 1, B = 10 G, and different β ( β = 0 .
20 is in the centerof the allowed range of 0 . < β < . phase is known as in [38], we obtain a distance limit of ≈
139 Mpc for a FAP of 0.1% and a FDP of 50% usingdata from a single detector. This is consistent with theestimate of ≈
100 Mpc reported in [38] (which assumed asmaller FAP). A real matched filter search will have anoverall unknown phase constant (see the last paragraphin Sec. IV B), which reduces the horizon distance to ≈
118 Mpc (for the same model parameters). Our cross-correlation matched filter limit yields an horizon distanceof ≈
103 Mpc, or 103 Mpc / √ ≈
119 Mpc when correct-ing for spectral leakage (see Fig. 6), in agreement withthe (real) matched filter.We finally note that signals with faster frequency evo-lutions are affected more by spectral leakage, for a fixedchoice of ∆ T SFT (satisfying the quarter cycle criterion).For example, the GW signal from a magnetar with β =0 .
26 would have a faster frequency evolution than thatfrom a source with β = 0 .
20 (with other source parame-ters unchanged; Fig. 7). The distance horizon we achievein the cross-correlation matched filter limit for β = 0 . T SFT equal to the one used for the β = 0 . ≈
216 Mpc. This is ≈ √
57% of the expectedmatched filter horizon of 287 Mpc (see Table I), worsethan what we would have expected for average spectralleakage losses of √ √
50% (see Fig. 6a). Such extremallosses are consistent with baselines very near to (but notexceeding) the maximum value set by the quarter-cyclecriterion. These losses can be improved by optimizing thesize of the baseline, given each frequency evolution, andis planned for future work. These results are summarizedin Table I.While a detailed study of the parameter space of themodel by [38] is beyond the scope of this paper, we also
TABLE I. Single-detector distance horizons for simulations inwhich the search is performed on the “correct” frequency-timetrack for with B = 10 G, M = 1 . M (cid:12) , R = 20 km andvarying values of β using the model proposed by [38]. Thesearch techniques used are matched filtering with unknownphase (MF), the cross-correlation matched filter limit ( χ MF,see Sec. IV B), and the cross-correlation stochastic limit (seeSec. IV A). β Value Distance Horizon (Mpc)MF χ MF Stochastic0.20 118 103 200.26 287 216 40 carried out several simulations to demonstrate the effec-tiveness of a semi-coherent approach in: (i) enhancing therobustness of the search against signal uncertainties whencompared to a matched-filter limit; and (ii) enhancing thesensitivity of the search when compared to a “stochasticapproach”. We do so by calculating the distance horizonsfor situations in which the assumed time-frequency trackdiffers from the actual signal by some amount. This dif-ference is quantified by an error ( δM, δR, δB ) on thevalues of the true signal parameters (
M, R, B ). Thesizes of these errors help determine the parameter spaceresolution for an effective search. The results of thesetests are summarized in Tables II and III.Because an error in signal parameters implies a mis-match between the true signal time-frequency evolutionand the time-frequency track adopted for the calculationof the ρ statistic, we expect the cross-correlation search tocompletely miss the signal in the limit of large coherencetimescales, T coh → T obs (approaching the matched-filterlimit, which is not robust against such deviations). Onthe other hand, in the limit of small coherence timescales, T coh → ∆ T SFT , while the search is expected to be robustagainst signal uncertainties, the sensitivity is significantlylower than the matched-filter case. Thus, for a given pa-rameter space resolution, one can define an optimal co-herence timescale, which can then be used to quantifythe distance reach of the semi-coherent regime (for givenFAP and FDP).We obtain the optimal coherence time ( T opt ) by calcu-lating the detection efficiency for given FAP (here, 0.1%)as a function of T coh , for a signal at a fixed distance. The T coh that is associated with the maximal detection effi-ciency is then used for a series of injections of varyingdistance, but fixed T coh . The distance that is associatedwith an efficiency of 50% (which is equivalent to a FDPof 50%) is then taken to be the distance horizon for thatstep size. The step sizes taken for each model parame-ter informs the size of the parameter space that a semi-coherent cross-correlation search should cover. We ransimulations with two classes of step size: “large” stepsthat correspond to a coarse grid in the parameter space,4and “small” steps that correspond to finer (and subse-quently, more computationally intensive) grid in the pa-rameter space.The results for the large steps are shown in Fig. 8,and summarized in Table II. Optimal coherence times,see Fig. 8 (left), are of O (1) s, which lead to maximaldetection distances around 20 −
30 Mpc (recovering only ≈
25% of the matched filter limit), see Fig. 8 (right).In the case where all three parameters are stepped si-multaneously ( δ All), the optimal coherence time is onlytwice the SFT baseline of ∆ T SFT = 0 .
25 s and providesno significant gains over the stochastic limit, see TableII.
TABLE II. Single-detector distance horizons for large steps ineach of the model parameters with, β = 0 . B = 10 G + δ B , M = 1 . M (cid:12) + δM , R = 20 km + δ R . The resultingdistance horizons are approximately 20-30 Mpc, which is up toa 50% improvement over the stochastic limit, a but only ≈ O (1) Mpc.Parameter Step size T opt (sec) Distance Horizon (Mpc)Semicoh Stochastic δB G 1 22 20 δM × − M (cid:12) δR
20 m 2 29 20 δ All As above 0.5 20 20 a A factor of 1.5 in distance horizon increases the expecteddetection rate by a factor of 1 . ≈ The small step sizes produce optimal coherence timesof as high as 256 s, see Fig. 9 (left), which lead to max-imal detection distances of ≈ ≈ VI. DISCUSSION AND CONCLUSION
We have explored the application of the cross-correlation technique described in [56] to a new classof intermediate duration GW signals of duration T obs (cid:46) TABLE III. Single-detector distance horizons for small stepsin each of the model parameters with, β = 0 . B =10 G + δ B , M = 1 . M (cid:12) + δM , R = 20 km + δ R . Thesimulation used T obs = 1024 s and ∆ T SFT = 0 .
25 s. The re-sulting distance horizons are approximately 60-80 Mpc, up tofour times as large as the stochastic limit, a and (cid:38)
75% of thematched filter limit. All errors of order O (1) Mpc.Parameter Step size T opt (sec) Distance Horizon (Mpc)Semicoh Stochastic δB G 64 61 20 δM × − M (cid:12)
256 73 20 δR . δ All As above 64 58 20 a A factor of 4 in distance horizon increases the expecteddetection rate by a factor of 4 = 64. s, specifically the bar mode instability model for mil-lisecond magnetars developed in [38]. In doing so, wehave corrected the statistical properties of the cross-correlation statistic reported in [56] for both the semi-coherent, and fully-coherent matched-filter limits. In ad-dition, we have done a cursory exploration of the param-eter space for this model.There are several parallels between limits of the cross-correlation method and other search techniques used forLIGO data analysis. Natural examples are the tech-niques derived from efforts to quantify the stochastic GWbackground. Two such methods are the Stochastic Tran-sient Analysis Multi-detector Pipeline (STAMP), a cross-power statistic widely used for LIGO all-sky searches[30, 98], and stochtrack, a seedless clustering algorithmthat has been tested on signal models comparable in du-ration to those considered here [31]. Both these methodsare similar (in spirit, if not necessarily implementation)to the stochastic limit of the cross-correlation approach.Because of their significant robustness against signaluncertainties (and relatively low computational costs)stochastic-inspired methods (as the two described above)are attractive for many search regimes, and especially asa first pass when searching for viable GW candidates withwide parameter spaces. On the other hand, the improve-ment in sensitivity (and therefore distance reach) enabledby the semi-coherent limit of the cross-correlation ap-proach lends itself to deeper searches. A potential wayto leverage the strengths of both regimes is to develop aframework in which a stochastic-inspired search is usedfor discovery, with semi-coherent cross-correlation fol-lowup for parameter estimation and refinement. Thiscould be done entirely within cross-correlation methoddescribed in this work, or by using an established stochas-tic technique (e.g. STAMP) for discovery and cross-correlation for follow-up.Overall, the results of our study are encouraging:5The tunable robustness versus sensitivity of the cross-correlation technique is well suited for intermediate-duration GW signals that evolve on timescales of 10 -10 s, and can reach astrophysically relevant distancehorizons with the expected noise characteristics of GWdetectors such as aLIGO and Virgo. However, a full pa-rameter space exploration is yet to be completed, as istesting on real instrument noise. Additionally, the trialsfactor for a full parameter space search will reduce, tosome extent, the idealized horizon distances calculatedhere. We intend to explore these aspects of the analyses in future work. ACKNOWLEDGMENTS
This work is supported by NSF grant PHY-1456447(PI: Corsi). B.O. acknowledges support from NSF grantsPHY-1206027, PHY-1544295, and PHY-1506311. A.C.and B.O. thank P. Meszaros for useful discussions in theearly stages of this work. A.C. also thanks C. Palombafor early discussions regarding continuous wave searches.This paper has been assigned LIGO Document No, ligo-p1500226. [1] B. Abbott, R. Abbott, R. Adhikari, A. Ageev, B. Allen,R. Amin, S. B. Anderson, W. G. Anderson, M. Araya,H. Armandula, et al.,
Search for gravitational waves as-sociated with the gamma ray burst GRB030329 usingthe LIGO detectors , Phys. Rev. D , 042002 (2005),arXiv:gr-qc/0501068.[2] B. Abbott, R. Abbott, R. Adhikari, J. Agresti, P. Ajith,B. Allen, R. Amin, S. B. Anderson, W. G. Anderson,M. Arain, et al., Search for gravitational wave radiationassociated with the pulsating tail of the SGR 1806-20 hy-perflare of 27 December 2004 using LIGO , Phys. Rev. D , 062003 (2007), arXiv:astro-ph/0703419.[3] F. Acernese, P. Amico, M. Alshourbagy, F. Antonucci,S. Aoudia, P. Astone, S. Avino, D. Babusci, G. Ballardin,F. Barone, et al., Gravitational waves by gamma-raybursts and the Virgo detector: the case of GRB 050915a ,Class. Quant. Grav. , S671 (2007).[4] B. Abbott, R. Abbott, R. Adhikari, J. Agresti, P. Ajith,B. Allen, R. Amin, S. B. Anderson, W. G. Anderson,M. Arain, et al., Search for gravitational waves associatedwith 39 gamma-ray bursts using data from the second,third, and fourth LIGO runs , Phys. Rev. D , 062004(2008), 0709.0766.[5] B. Abbott, R. Abbott, R. Adhikari, J. Agresti, P. Ajith,B. Allen, R. Amin, S. B. Anderson, W. G. Anderson,M. Arain, et al., Implications for the Origin of GRB070201 from LIGO Observations , Astrophys. J. ,1419 (2008), 0711.1163.[6] F. Acernese, M. Alshourbagy, P. Amico, F. Antonucci,S. Aoudia, K. G. Arun, P. Astone, S. Avino, L. Baggio,G. Ballardin, et al.,
Search for gravitational waves asso-ciated with GRB 050915a using the Virgo detector , Class.Quant. Grav. , 225001 (2008), 0803.0376.[7] B. P. Abbott, R. Abbott, R. Adhikari, P. Ajith, B. Allen,G. Allen, R. S. Amin, S. B. Anderson, W. G. Ander-son, M. A. Arain, et al., Stacked Search for GravitationalWaves from the 2006 SGR 1900+14 Storm , Astrophys.J. Lett. , L68 (2009), 0905.0005.[8] B. P. Abbott, R. Abbott, F. Acernese, R. Adhikari,P. Ajith, B. Allen, G. Allen, M. Alshourbagy, R. S. Amin,S. B. Anderson, et al.,
Search For Gravitational-waveBursts Associated with Gamma-ray Bursts using Datafrom LIGO Science Run 5 and Virgo Science Run 1 , As-trophys. J. , 1438 (2010), 0908.3824. [9] J. Abadie, B. P. Abbott, R. Abbott, T. Accadia, F. Ac-ernese, R. Adhikari, P. Ajith, B. Allen, G. Allen,E. Amador Ceron, et al.,
Search for Gravitational-waveInspiral Signals Associated with Short Gamma-ray BurstsDuring LIGO’s Fifth and Virgo’s First Science Run , As-trophys. J. , 1453 (2010), 1001.0165.[10] J. Abadie, B. P. Abbott, R. Abbott, R. Adhikari,P. Ajith, B. Allen, G. Allen, E. Amador Ceron, R. S.Amin, S. B. Anderson, et al.,
Search for gravitationalwaves associated with the August 2006 timing glitch of theVela pulsar , Phys. Rev. D , 042001 (2011), 1011.1357.[11] J. Abadie, B. P. Abbott, R. Abbott, M. Abernathy,T. Accadia, F. Acernese, C. Adams, R. Adhikari, C. Af-feldt, B. Allen, et al., Search for Gravitational WaveBursts from Six Magnetars , Astrophys. J. Lett. , L35(2011), 1011.4079.[12] J. Abadie, B. P. Abbott, T. D. Abbott, R. Abbott,M. Abernathy, C. Adams, R. Adhikari, C. Affeldt,P. Ajith, et al.,
Implications For The Origin Of GRB051103 From LIGO Observations , ArXiv e-prints (2012).[13] J. Aasi, J. Abadie, B. P. Abbott, R. Abbott, T. Abbott,M. R. Abernathy, T. Accadia, F. Acernese, C. Adams,T. Adams, et al.,
Search for long-lived gravitational-wavetransients coincident with long gamma-ray bursts , Phys.Rev. D , 122004 (2013), 1309.6160.[14] J. Aasi et al. (VIRGO, LIGO Scientific), Methods andresults of a search for gravitational waves associated withgamma-ray bursts using the GEO600, LIGO, and Virgodetectors , Phys. Rev. D
D89 , 122004 (2014), 1405.1053.[15] J. Aasi et al. (VIRGO, IPN, LIGO Scientific),
Searchfor gravitational waves associated with γ -ray bursts de-tected by the Interplanetary Network , Phys. Rev. Lett. , 011102 (2014), 1403.6639.[16] B. Abbott et al. (LIGO Scientific), Setting upper limitson the strength of periodic gravitational waves using thefirst science data from the GEO 600 and LIGO detectors ,Phys. Rev. D
D69 , 082004 (2004), gr-qc/0308050.[17] B. Abbott et al. (LIGO Scientific),
Limits on gravita-tional wave emission from selected pulsars using LIGOdata , Phys. Rev. Lett. , 181103 (2005), gr-qc/0410007.[18] B. Abbott, R. Abbott, R. Adhikari, J. Agresti, P. Ajith,B. Allen, R. Amin, S. B. Anderson, W. G. Anderson,M. Arain, et al., Upper limits on gravitational wave emis-sion from 78 radio pulsars , Phys. Rev. D , 042001(2007), arXiv:gr-qc/0702039. E ff i c i e n c y d inj =22 T opt =1 .
10 20 30 40 500.00.20.40.60.81.01.2 s t e p B d max =22 ± E ff i c i e n c y d inj =28 T opt =2 .
10 20 30 40 50 600.00.20.40.60.81.01.2 s t e p M d max =28 ± E ff i c i e n c y d inj =29 T opt =2 .
10 20 30 40 50 600.00.20.40.60.81.01.2 s t e p R d max =29 ± Coherence Time (s) E ff i c i e n c y d inj =20 T opt =0 .
10 20 30 40 50
Distance (Mpc) s t e p A ll d max =20 ± Coherence Time / Distance Optimizations (Large Steps)
FIG. 8. Efficiency (1-FDP) plots for large steps δB = 10 G (blue), δM = 5 × − M (cid:12) (green), δR = 20 m (red) and allthree combined (purple). All plots assume FAP=0.1% and distances are extracted using FDP=50% (black dotted line and grayshaded area). On the left, optimal coherence time plots. The signal is injected at a constant distance, T coh is then varied tofind the value that maximizes detection efficiency ( T opt ). On the right, T coh is fixed at the optimum value for each step, andthen distance is varied. The result is fit by an asymmetric sigmoid of the form sig(x) = [1 + exp(p { x − p } )] − / p (where p , p , p are constants to be fit), which is then used to interpolate and determine the max distance ( d max ).[19] B. Abbott, R. Abbott, R. Adhikari, J. Agresti, P. Ajith,B. Allen, R. Amin, S. B. Anderson, W. G. Anderson,M. Arain, et al., Searches for periodic gravitational wavesfrom unknown isolated sources and Scorpius X-1: Resultsfrom the second LIGO science run , Phys. Rev. D , 082001 (2007), arXiv:gr-qc/0605028.[20] B. Abbott, R. Abbott, R. Adhikari, P. Ajith, B. Allen,G. Allen, R. Amin, S. B. Anderson, W. G. Anderson,M. A. Arain, et al., Beating the Spin-Down Limit onGravitational Wave Emission from the Crab Pulsar , As- E ff i c i e n c y d inj =61 T opt =64 .
30 40 50 60 70 80 900.00.20.40.60.81.01.2 s t e p B _ s m a ll d max =61 ± E ff i c i e n c y d inj =73 T opt =256 .
40 50 60 70 80 90 1000.00.20.40.60.81.01.2 s t e p M _ s m a ll d max =73 ± E ff i c i e n c y d inj =76 T opt =256 .
40 50 60 70 80 90 1000.00.20.40.60.81.01.2 s t e p R _ s m a ll d max =76 ± Coherence Time (s) E ff i c i e n c y d inj =58 T opt =64 .
30 40 50 60 70 80 90
Distance (Mpc) s t e p A ll _ s m a ll d max =58 ± Coherence Time / Distance Optimizations (Small Steps)
FIG. 9. Efficiency (1-FDP) plots for small steps in δB = 10 G (blue), δM = 5 × − M (cid:12) (green), δR = 0 . T coh is then variedto find the value that maximizes detection efficiency ( T opt ). On the right, T coh is fixed at the optimum value for each step,and then distance is varied. The result is fit by an asymmetric sigmoid of the form sig(x) = [1 + exp(p { x − p } )] − / p (where p , p , p are constants to be fit), which is then used to interpolate and determine the max distance ( d max ). trophys. J. Lett. , L45 (2008), 0805.4758.[21] J. Abadie, B. P. Abbott, R. Abbott, M. Abernathy,C. Adams, R. Adhikari, P. Ajith, B. Allen, G. Allen,E. Amador Ceron, et al., First Search for GravitationalWaves from the Youngest Known Neutron Star , Astro- phys. J. , 1504 (2010), 1006.2535.[22] J. Abadie, B. P. Abbott, R. Abbott, M. Abernathy,T. Accadia, F. Acernese, C. Adams, R. Adhikari, C. Af-feldt, B. Allen, et al.,
Beating the Spin-down Limit onGravitational Wave Emission from the Vela Pulsar , As- trophys. J. , 93 (2011), 1104.2712.[23] J. Abadie, B. P. Abbott, R. Abbott, M. Abernathy,T. Accadia, F. Acernese, C. Adams, R. Adhikari,P. Ajith, B. Allen, et al., Directional Limits on Persis-tent Gravitational Waves Using LIGO S5 Science Data ,Phys. Rev. Lett. , 271102 (2011), 1109.1809.[24] J. Aasi, J. Abadie, B. P. Abbott, R. Abbott, T. Abbott,M. R. Abernathy, T. Accadia, F. Acernese, C. Adams,T. Adams, et al.,
Directed search for continuous gravita-tional waves from the Galactic center , Phys. Rev. D ,102002 (2013), 1309.6221.[25] J. Aasi et al. (LIGO Scientific), Gravitational waves fromknown pulsars: results from the initial detector era , As-trophys. J. , 119 (2014), 1309.4027.[26] J. Aasi, B. P. Abbott, R. Abbott, T. Abbott, M. R. Aber-nathy, F. Acernese, K. Ackley, C. Adams, T. Adams,T. Adams, et al.,
Narrow-band search of continuousgravitational-wave signals from Crab and Vela pulsarsin Virgo VSR4 data , Phys. Rev. D , 022004 (2015),1410.8310.[27] J. Aasi et al. (VIRGO, LIGO Scientific), Directed searchfor gravitational waves from Scorpius X-1 with ini-tial LIGO data , Phys. Rev. D
D91 , 062008 (2015),1412.0605.[28] J. Aasi et al.,
Searches for continuous gravitational wavesfrom nine young supernova remnants (2014), 1412.5942.[29] B. Abbott, R. Abbott, R. Adhikari, J. Agresti, P. Ajith,B. Allen, R. Amin, S. B. Anderson, W. G. Anderson,M. Arain, et al.,
Upper limit map of a background ofgravitational waves , Phys. Rev. D , 082003 (2007),arXiv:astro-ph/0703234.[30] E. Thrane, S. Kandhasamy, C. D. Ott, W. G. An-derson, N. L. Christensen, M. W. Coughlin, S. Dor-sher, S. Giampanis, V. Mandic, A. Mytidis, et al., Long gravitational-wave transients and associated detec-tion strategies for a network of terrestrial interferome-ters , Phys. Rev. D , 083004 (2011), 1012.2150.[31] E. Thrane and M. Coughlin, Seedless clustering in all-skysearches for gravitational-wave transients , Phys. Rev. D (2014).[32] E. Thrane, V. Mandic, and N. Christensen, Detectingvery long-lived gravitational-wave transients lasting hoursto weeks , Phys. Rev. D
D91 , 104021 (2015), 1501.06648.[33] The LIGO Scientific Collaboration, the Virgo Collab-oration, B. P. Abbott, R. Abbott, T. D. Abbott,M. R. Abernathy, F. Acernese, K. Ackley, C. Adams,T. Adams, et al.,
An all-sky search for long-durationgravitational wave transients with LIGO , ArXiv e-prints(2015), 1511.04398.[34] E. Thrane, V. Mandic, and N. Christensen,
Detectingvery long-lived gravitational-wave transients lasting hoursto weeks , Phys. Rev. D , 104021 (2015), 1501.06648.[35] M. H. van Putten, Proposed Source of Gravitational Ra-diation from a Torus around a Black Hole , Phys. Rev.Lett. , 091101 (2001), astro-ph/0107007.[36] S. Dall’Osso and L. Stella, Newborn magnetars as sourcesof gravitational radiation: constraints from high energyobservations of magnetar candidates , Astrophys. SpaceSci. , 119 (2007), astro-ph/0702075.[37] A. L. Piro and E. Pfahl,
Fragmentation of CollapsarDisks and the Production of Gravitational Waves , As-trophys. J. , 1173 (2007), astro-ph/0610696.[38] A. Corsi and P. M´esz´aros,
Gamma-ray Burst AfterglowPlateaus and Gravitational Waves: Multi-messenger Sig- nature of a Millisecond Magnetar? , Astrophys. J. ,1171 (2009), 0907.2290.[39] M. H. P. M. van Putten,
Gravitational Waveforms ofKerr Black Holes Interacting with High-Density Matter , Astrophys. J. Lett. , L91 (2008).[40] C. D. Ott,
TOPICAL REVIEW: The gravitational-wave signature of core-collapse supernovae , Class. Quant.Grav. , 063001 (2009), 0809.0695.[41] A. L. Piro and C. D. Ott, Supernova Fallback onto Mag-netars and Propeller-powered Supernovae , Astrophys. J. , 108 (2011), 1104.0252.[42] A. L. Piro and E. Thrane,
Gravitational Waves from Fall-back Accretion onto Neutron Stars , Astrophys. J. , 63(2012), 1207.3805.[43] D. D. Doneva, K. D. Kokkotas, and P. Pnigouras,
A grav-itational wave afterglow in binary neutron star mergers ,ArXiv e-prints (2015).[44] J. A. Nousek, C. Kouveliotou, D. Grupe, K. L. Page,J. Granot, E. Ramirez-Ruiz, S. K. Patel, D. N. Burrows,V. Mangano, S. Barthelmy, et al.,
Evidence for a Canoni-cal Gamma-Ray Burst Afterglow Light Curve in the SwiftXRT Data , Astrophys. J. , 389 (2006), arXiv:astro-ph/0508332.[45] B. Zhang, Y. Z. Fan, J. Dyks, S. Kobayashi, P. M´esz´aros,D. N. Burrows, J. A. Nousek, and N. Gehrels,
PhysicalProcesses Shaping Gamma-Ray Burst X-Ray AfterglowLight Curves: Theoretical Implications from the SwiftX-Ray Telescope Observations , Astrophys. J. , 354(2006), arXiv:astro-ph/0508321.[46] E. Liang, B. Zhang, and B. Zhang, Astrophys. J. ,565 (2007).[47] R. Starling et al.,
Swift captures the spectrally evolvingprompt emission of grb070616 , Mon. Not. R. Astron. Soc. , 504S (2008).[48] M. G. Bernardini, R. Margutti, J. Mao, E. Zaninoni,and G. Chincarini,
The X-ray light curve of gamma-raybursts: clues to the central engine , Astron. Astrophys. , A3 (2012), 1112.1058.[49] B. P. Gompertz, P. T. O’Brien, G. A. Wynn, andA. Rowlinson,
Can magnetar spin-down power extendedemission in some short GRBs? , Mon. Not. R. Astron.Soc. , 1745 (2013), 1302.3643.[50] A. Rowlinson, P. T. O’Brien, B. D. Metzger, N. R. Tan-vir, and A. J. Levan,
Signatures of magnetar central en-gines in short GRB light curves , Mon. Not. R. Astron.Soc. , 1061 (2013), 1301.0629.[51] S. X. Yi, Z. G. Dai, X. F. Wu, and F. Y. Wang,
X-RayAfterglow Plateaus of Long Gamma-Ray Bursts: Fur-ther Evidence for Millisecond Magnetars , ArXiv e-prints(2014), 1401.1601.[52] S. Chandrasekhar,
Solutions of Two Problems in the The-ory of Gravitational Radiation , Phys. Rev. Lett. , 611(1970).[53] J. L. Friedman and B. F. Schutz, Secular instability ofrotating Newtonian stars , Astrophys. J. , 281 (1978).[54] D. Lai and S. L. Shapiro,
Gravitational radiation fromrapidly rotating nascent neutron stars , Astrophys. J. ,259 (1995), arXiv:astro-ph/9408053.[55] B. Zhang and P. M´esz´aros,
Gamma-Ray Burst Afterglowwith Continuous Energy Injection: Signature of a HighlyMagnetized Millisecond Pulsar , Astrophys. J. Lett. ,L35 (2001), arXiv:astro-ph/0011133.[56] S. Dhurandhar, B. Krishnan, H. Mukhopadhyay, andJ. T. Whelan,
Cross-correlation search for periodic gravitational waves , Phys. Rev. D , 082001 (2008),0712.1578.[57] N. J. Cornish and J. D. Romano, Towards a unified treat-ment of gravitational-wave data analysis , Phys. Rev. D
D87 , 122003 (2013), 1305.2934.[58] R. Prix, S. Giampanis, and C. Messenger,
Search methodfor long-duration gravitational-wave transients from neu-tron stars , Phys. Rev. D
D84 , 023007 (2011), 1104.1704.[59] B. F. Schutz,
The Detection of Gravitational Waves (Cambridge University press, 1991).[60] P. Jaranowski, A. Kr´olak, and B. F. Schutz,
Data anal-ysis of gravitational-wave signals from spinning neutronstars: The signal and its detection , Phys. Rev. D ,063001 (1998), arXiv:gr-qc/9804014.[61] B. J. Owen, Search templates for gravitational waves frominspiraling binaries: Choice of template spacing , Phys.Rev. D , 6749 (1996), arXiv:gr-qc/9511032.[62] P. R. Brady, T. Creighton, C. Cutler, and B. F. Schutz, Searching for periodic sources with LIGO , Phys. Rev. D , 2101 (1998), arXiv:gr-qc/9702050.[63] B. J. Owen and B. S. Sathyaprakash, Matched filteringof gravitational waves from inspiraling compact binaries:Computational cost and template placement , Phys. Rev.D , 022002 (1999), arXiv:gr-qc/9808076.[64] K. Wette, B. J. Owen, B. Allen, M. Ashley, J. Betzwieser,N. Christensen, T. D. Creighton, V. Dergachev, I. Gho-lami, E. Goetz, et al., Searching for gravitational wavesfrom Cassiopeia A with LIGO , Class. Quant. Grav. ,235011 (2008), 0802.3332.[65] T. A. Apostolatos, Search templates for gravitationalwaves from precessing, inspiraling binaries , Phys. Rev.D , 605 (1995), URL http://link.aps.org/doi/10.1103/PhysRevD.52.605 .[66] L. Lindblom, B. J. Owen, and D. A. Brown, Model wave-form accuracy standards for gravitational wave data anal-ysis , Phys. Rev. D , 124020 (2008), 0809.3844.[67] L. Lindblom, J. G. Baker, and B. J. Owen, Improvedtime-domain accuracy standards for model gravitationalwaveforms , Phys. Rev. D , 084020 (2010), 1008.1803.[68] T. Damour, A. Nagar, and M. Trias, Accuracy and effec-tualness of closed-form, frequency-domain waveforms fornonspinning black hole binaries , Phys. Rev. D , 024006(2011), 1009.5998.[69] I. Hinder, A. Buonanno, M. Boyle, Z. B. Etienne,J. Healy, N. K. Johnson-McDaniel, A. Nagar, H. Nakano,Y. Pan, H. P. Pfeiffer, et al., Error-analysis and compar-ison to analytical models of numerical waveforms pro-duced by the NRAR Collaboration , Classical and Quan-tum Gravity , 025012 (2013), 1307.5307.[70] L. Sampson, N. Cornish, and N. Yunes, Mismodelingin gravitational-wave astronomy: The trouble with tem-plates , Phys. Rev. D , 064037 (2014), 1311.4898.[71] P. R. Brady and T. Creighton, Searching for periodicsources with LIGO. II. Hierarchical searches , Phys. Rev.D , 082001 (2000), arXiv:gr-qc/9812014.[72] B. Allen, χ time-frequency discriminator for gravita-tional wave detection , Phys. Rev. D , 062001 (2005),arXiv:gr-qc/0405045.[73] C. Cutler, I. Gholami, and B. Krishnan, Improved stack-slide searches for gravitational-wave pulsars , Phys. Rev.D , 042004 (2005), arXiv:gr-qc/0505082.[74] B. P. Abbott, R. Abbott, F. Acernese, R. Adhikari,P. Ajith, B. Allen, G. Allen, M. Alshourbagy, R. S. Amin,S. B. Anderson, et al., An upper limit on the stochastic gravitational-wave background of cosmological origin , Na-ture (London) , 990 (2009), 0910.5772.[75] The LIGO Scientific Collaboration, the Virgo Collabora-tion, J. Aasi, J. Abadie, B. P. Abbott, R. Abbott, T. Ab-bott, M. R. Abernathy, T. Accadia, F. Acernese, et al.,
Searching for stochastic gravitational waves using datafrom the two co-located LIGO Hanford detectors , ArXive-prints (2014), 1410.6211.[76] C. T. Y. Chung, A. Melatos, B. Krishnan, and J. T. Whe-lan,
Designing a cross-correlation search for continuous-wave gravitational radiation from a neutron star in thesupernova remnant snr 1987a , Mon. Not. R. Astron. Soc. , 2650 (2011).[77] J. T. Whelan, S. Sundaresan, Y. Zhang, and P. Peiris,
Model-based cross-correlation search for gravitationalwaves from scorpius x-1 , Phys. Rev. D (2015).[78] S. Sundaresan and J. T. Whelan, Window-ing and leakage in the cross-correlation searchfor periodic gravitational waves , Tech. Rep.,LIGO Technical Document (2012), URL http://dcc.ligo.org/LIGO-T1200431-v1/public .[79] D. B. Percival and A. T. Walden,
Spectral Analysisfor Physical Applications (Cambridge University Press,1993).[80] J. D. Creighton and W. Anderson,
Gravitational-WavePhysics and Astronomy (Wiley-VCH, 2011).[81] D. A. Frail, S. R. Kulkarni, R. Sari, S. G. Djorgovski,J. S. Bloom, T. J. Galama, D. E. Reichart, E. Berger,F. A. Harrison, P. A. Price, et al.,
Beaming in Gamma-Ray Bursts: Evidence for a Standard Energy Reservoir , Astrophys. J. Lett. , L55 (2001), astro-ph/0102282.[82] A. Panaitescu and P. Kumar,
Fundamental Physical Pa-rameters of Collimated Gamma-Ray Burst Afterglows , Astrophys. J. Lett. , L49 (2001), astro-ph/0108045.[83] J. Abadie et al. (LIGO Scientific),
Calibration of theLIGO Gravitational Wave Detectors in the Fifth Sci-ence Run , Nucl. Instrum. Methods
A624 , 223 (2010),1007.3973.[84] E. E. Flanagan,
The Sensitivity of the laser interferome-ter gravitational wave observatory (LIGO) to a stochasticbackground, and its dependence on the detector orienta-tions , Phys. Rev. D
D48 , 2389 (1993), astro-ph/9305029.[85] B. Allen and J. D. Romano,
Detecting a stochasticbackground of gravitational radiation: Signal processingstrategies and sensitivities , Phys. Rev. D
D59 , 102001(1999), gr-qc/9710117.[86] N. J. Cornish,
Mapping the gravitational wave back-ground , Class. Quant. Grav. , 4277 (2001), astro-ph/0105374.[87] S. W. Ballmer, A Radiometer for stochastic gravita-tional waves , Class. Quant. Grav. , S179 (2006), gr-qc/0510096.[88] L. S. Finn, Detection, measurement and gravitational ra-diation , Phys. Rev. D
D46 , 5236 (1992), gr-qc/9209010.[89] P. Astone, S. D’Antonio, S. Frasca, and C. Palomba,
A method for detection of known sources of continuousgravitational wave signals in non-stationary data , Class.Quant. Grav. , 194016 (2010).[90] B. Krishnan, A. M. Sintes, M. A. Papa, B. F. Schutz,S. Frasca, and C. Palomba, The Hough transform searchfor continuous gravitational waves , Phys. Rev. D
D70 ,082001 (2004), gr-qc/0407001.[91] V. Dergachev,
Description of powerflux algorithms andimplementation (2005), lIGO-T050186, see dcc.ligo.org. [92] V. Dergachev, Loosely coherent searches for sets of well-modeled signals , Phys. Rev. D
D85 , 062003 (2012),1110.3297.[93] E. Goetz and K. Riles,
Coherently combining short datasegments for all-sky semi-coherent continuous gravita-tional wave searches (2015), 1510.06820.[94] R. Margutti, E. Zaninoni, M. G. Bernardini, G. Chincar-ini, F. Pasotti, C. Guidorzi, L. Angelini, D. N. Burrows,M. Capalbi, P. A. Evans, et al.,
The prompt-afterglowconnection in gamma-ray bursts: a comprehensive statis-tical analysis of Swift X-ray light curves , Mon. Not. R.Astron. Soc. , 729 (2013), 1203.1059.[95] C. E. Shannon,
Communication in the presence of noise ,Proceedings of the IEEE , 447 (1998). [96] C. Tinney, R. Stathakis, R. Cannon, T. Galama,M. Wieringa, D. A. Frail, S. R. Kulkarni, J. L. Higdon,R. Wark, J. S. Bloom, et al., GRB 980425 , Int. Astron.Union , 3 (1998).[97] T. J. Galama, P. M. Vreeswijk, J. van Paradijs, C. Kou-veliotou, T. Augusteijn, H. B¨ohnhardt, J. P. Brewer,V. Doublier, J.-F. Gonzalez, B. Leibundgut, et al.,
Anunusual supernova in the error box of the γ -ray burst of25 April 1998 , Nature (London) , 670 (1998), astro-ph/9806175.[98] M. W. Coughlin, LIGO Scientific Collaboration, andVirgo Collaboration, Identification of long-duration noisetransients in LIGO and Virgo , Class. Quant. Grav.28