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 used broad-band imaging data for 10 cool-core brightest cluster galaxies (BCGs) and conducted a Bayesian analysis using stellar population synthesis to determine the likely properties of the constituent stellar populations. Determination of ongoing star formation rates (SFRs), in particular, has a direct impact on our understanding of the cooling of the intracluster medium (ICM), star formation and AGN-regulated feedback. Our model consists of an old stellar population and a series of young stellar components. We calculated marginalized posterior probability distributions for various model parameters and obtained 68% plausible intervals from them. The 68% plausible interval on the SFRs is broad, owing to a wide range of models that are capable of fitting the data, which also explains the wide dispersion in the star formation rates available in the literature. The ranges of possible SFRs are robust and highlight the strength in such a Bayesian analysis. The SFRs are correlated with the X-ray mass deposition rates (the former are factors of 4 to 50 lower than the latter), implying a picture where the cooling of the ICM is a contributing factor to star formation in cool-core BCGs. We find that 9 out of 10 BCGs have been experiencing starbursts since 6 Gyr ago. While four out of 9 BCGs seem to require continuous SFRs, 5 out of 9 seem to require periodic star formation on intervals ranging from 20 Myr to 200 Myr. This time scale is similar to the cooling-time of the ICM in the central (< 5 kpc) regions.
The parameter space for continuous gravitational waves (GWs) can be divided into amplitude parameters (signal amplitude, inclination and polarization angles describing the orientation of the source, and an initial phase) and phase-evolution parameters. The division is useful in part because the Jaranowski-Krolak-Schutz (JKS) coordinates on the four-dimensional amplitude parameter space allow the GW signal to be written as a linear combination of four template waveforms with the JKS coordinates as coefficients. We define a new set of coordinates on the amplitude parameter space, with the same properties, which is more closely connected to the physical amplitude parameters. These naturally divide into two pairs of Cartesian-like coordinates on two-dimensional subspaces, one corresponding to left- and the other to right-circular polarization. We thus refer to these as CPF (circular polarization factored) coordinates. The corresponding two sets of polar coordinates (known as CPF-polar) can be related in a simple way to the physical parameters. We illustrate some simplifying applications for these various coordinate systems, such as: a calculation of Jacobians between various coordinate systems; an illustration of the signal coordinate singularities associated with left- and right-circular polarization, which correspond to the origins of the two two-dimensional subspaces; and an elucidation of the form of the log-likelihood ratio between hypotheses of Gaussian noise with and without a continuous GW signal. These are used to illustrate some of the prospects for approximate evaluation of a Bayesian detection statistic defined by marginalization over the physical parameter space. Additionally, in the presence of simplifying assumptions about the observing geometry, we are able to explicitly evaluate the integral for the Bayesian detection statistic, and compare it to the approximate results.
Uncertainty in the calibration of gravitational-wave (GW) detector data leads to systematic errors which must be accounted for in setting limits on the strength of GW signals. When cross-correlation measurements are made using data from a pair of instruments, as in searches for a stochastic GW background, the calibration uncertainties of the individual instruments can be combined into an uncertainty associated with the pair. With the advent of multi-baseline GW observation (e.g., networks consisting of multiple detectors such as the LIGO observatories and Virgo), a more sophisticated treatment is called for. We describe how the correlations between calibration factors associated with different pairs can be taken into account by marginalizing over the uncertainty associated with each instrument.
We describe an F-statistic search for continuous gravitational waves from galactic white-dwarf binaries in simulated LISA Data. Our search method employs a hierarchical template-grid based exploration of the parameter space. In the first stage, candidate sources are identified in searches using different simulated laser signal combinations (known as TDI variables). Since each source generates a primary maximum near its true Doppler parameters (intrinsic frequency and sky position) as well as numerous secondary maxima of the F-statistic in Doppler parameter space, a search for multiple sources needs to distinguish between true signals and secondary maxima associated with other, louder signals. Our method does this by applying a coincidence test to reject candidates which are not found at nearby parameter space positions in searches using each of the three TDI variables. For signals surviving the coincidence test, we perform a fully coherent search over a refined parameter grid to provide an accurate parameter estimation for the final candidates. Suitably tuned, the pipeline is able to extract 1989 true signals with only 5 false alarms. The use of the rigid adiabatic approximation allows recovery of signal parameters with errors comparable to statistical expectations, although there is still some systematic excess with respect to statistical errors expected from Gaussian noise. An experimental iterative pipeline with seven rounds of signal subtraction and re-analysis of the residuals allows us to increase the number of signals recovered to a total of 3419 with 29 false alarms.
We report on our F-statistic search for white-dwarf binary signals in the Mock LISA Data Challenge 1B (MLDC1B). We focus in particular on the improvements in our search pipeline since MLDC1, namely refinements in the search pipeline and the use of a more accurate detector response (rigid adiabatic approximation). The search method employs a hierarchical template-grid based exploration of the parameter space, using a coincidence step to distinguish between primary (``true) and secondary maxima, followed by a final (multi-TDI) ``zoom stage to provide an accurate parameter estimation of the final candidates.
The F-statistic is an optimal detection statistic for continuous gravitational waves, i.e., long-duration (quasi-)monochromatic signals with slowly-varying intrinsic frequency. This method was originally developed in the context of ground-based detectors, but it is equally applicable to LISA where many signals fall into this class of signals. We report on the application of a LIGO/GEO F-statistic code to LISA data-analysis using the long-wavelength limit (LWL), and we present results of our search for white-dwarf binary signals in the first Mock LISA Data Challenge. Somewhat surprisingly, the LWL is found to be sufficient -- even at high frequencies -- for detection of signals and their accurate localization on the sky and in frequency, while a more accurate modelling of the TDI response only seems necessary to correctly estimate the four amplitude parameters.
