No Arabic abstract
In this paper we study the use of cross-correlations between multiple gravitational wave (GW) data streams for detecting long-lived periodic signals. Cross-correlation searches between data from multiple detectors have traditionally been used to search for stochastic GW signals, but recently they have also been used in directed searches for periodic GWs. Here we further adapt the cross-correlation statistic for periodic GW searches by taking into account both the non-stationarity and the long term-phase coherence of the signal. We study the statistical properties and sensitivity of this search, its relation to existing periodic wave searches, and describe the precise way in which the cross-correlation statistic interpolates between semi-coherent and fully-coherent methods. Depending on the maximum duration over we wish to preserve phase coherence, the cross-correlation statistic can be tuned to go from a standard cross-correlation statistic using data from distinct detectors, to the semi-coherent time-frequency methods with increasing coherent time baselines, and all the way to a full coherent search. This leads to a unified framework for studying periodic wave searches and can be used to make informed trade-offs between computational cost, sensitivity, and robustness against signal uncertainties.
The cross-correlation search has been previously applied to map the gravitational wave (GW) stochastic background in the sky and also to target GW from rotating neutron stars/pulsars. Here we investigate how the cross-correlation method can be used to target a small region in the sky spanning at most a few pixels, where a pixel in the sky is determined by the diffraction limit which depends on the (i) baseline joining a pair of detectors and (ii) detector bandwidth. Here as one of the promising targets, we consider the Virgo cluster - a hot spot spanning few pixels - which could contain, as estimates suggest $sim 10^{11}$ neutron stars, of which a small fraction would continuously emit GW in the bandwidth of the detectors. For the detector baselines, we consider advanced detector pairs among LCGT, LIGO, Virgo, ET etc. Our results show that sufficient signal to noise can be accumulated with integration times of the order of a year. The results improve for the multibaseline search. This analysis could as well be applied to other likely hot spots in the sky and other possible pairs of detectors.
We consider the cross-correlation search for periodic GWs and its potential application to the LMXB Sco X-1. This method coherently combines data from different detectors at the same time, as well as different times from the same or different detectors. By adjusting the maximum time offset between a pair of data segments to be coherently combined, one can tune the method to trade off sensitivity and computing costs. In particular, the detectable signal amplitude scales as the inverse fourth root of this coherence time. The improvement in amplitude sensitivity for a search with a coherence time of 1hr, compared with a directed stochastic background search with 0.25Hz wide bins is about a factor of 5.4. We show that a search of 1yr of data from Advanced LIGO and Advanced Virgo with a coherence time of 1hr would be able to detect GWs from Sco X-1 at the level predicted by torque balance over a range of signal frequencies from 30-300Hz; if the coherence time could be increased to 10hr, the range would be 20-500Hz. In addition, we consider several technical aspects of the cross-correlation method: We quantify the effects of spectral leakage and show that nearly rectangular windows still lead to the most sensitive search. We produce an explicit parameter-space metric for the cross-correlation search in general and as applied to a neutron star in a circular binary system. We consider the effects of using a signal template averaged over unknown amplitude parameters: the search is sensitive to a combination of the intrinsic signal amplitude and the inclination of the neutron star rotation axis, and the peak of the expected detection statistic is systematically offset from the true signal parameters. Finally, we describe the potential loss of SNR due to unmodelled effects such as signal phase acceleration within the Fourier transform timescale and gradual evolution of the spin frequency.
We describe the application of the lattice covering problem to the placement of templates in a search for continuous gravitational waves from the low-mass X-Ray binary Scorpius X-1. Efficient placement of templates to cover the parameter space at a given maximum mismatch is an application of the sphere covering problem, for which an implementation is available in the LatticeTiling software library. In the case of Sco X-1, potential correlations, in both the prior uncertainty and the mismatch metric, between the orbital period and orbital phase, lead to complications in the efficient construction of the lattice. We define a shearing coordinate transformation which simultaneously minimizes both of these sources of correlation, and allows us to take advantage of the small prior orbital period uncertainty. The resulting lattices have a factor of about 3 fewer templates than the corresponding parameter space grids constructed by the prior straightforward method, allowing a more sensitive search at the same computing cost and maximum mismatch.
We report on an all-sky search for periodic gravitational waves in the frequency band 20-475 Hz and with a frequency time derivative in the range of [-1.0, +0.1]e-8 Hz/s. Such a signal could be produced by a nearby spinning and slightly non-axisymmetric isolated neutron star in our galaxy. This search uses the data from Advanced LIGOs first observational run, O1. No periodic gravitational wave signals were observed, and upper limits were placed on their strengths. The lowest upper limits on worst-case (linearly polarized) strain amplitude h0 are 4e-25 near 170 Hz. For a circularly polarized source (most favorable orientation), the smallest upper limits obtained are 1.5e-25. These upper limits refer to all sky locations and the entire range of frequency derivative values. For a population-averaged ensemble of sky locations and stellar orientations, the lowest upper limits obtained for the strain amplitude are 2.5e-25.
We report on an all-sky search with the LIGO detectors for periodic gravitational waves in the frequency range 50-1000 Hz and with the frequencys time derivative in the range -1.0E-8 Hz/s to zero. Data from the fourth LIGO science run (S4) have been used in this search. Three different semi-coherent methods of transforming and summing strain power from Short Fourier Transforms (SFTs) of the calibrated data have been used. The first, known as StackSlide, averages normalized power from each SFT. A weighted Hough scheme is also developed and used, and which also allows for a multi-interferometer search. The third method, known as PowerFlux, is a variant of the StackSlide method in which the power is weighted before summing. In both the weighted Hough and PowerFlux methods, the weights are chosen according to the noise and detector antenna-pattern to maximize the signal-to-noise ratio. The respective advantages and disadvantages of these methods are discussed. Observing no evidence of periodic gravitational radiation, we report upper limits; we interpret these as limits on this radiation from isolated rotating neutron stars. The best population-based upper limit with 95% confidence on the gravitational-wave strain amplitude, found for simulated sources distributed isotropically across the sky and with isotropically distributed spin-axes, is 4.28E-24 (near 140 Hz). Strict upper limits are also obtained for small patches on the sky for best-case and worst-case inclinations of the spin axes.