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.
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 construction of optimal template banks for matched-filtering searches is an example of the sphere covering problem. For parameter spaces with constant-coefficient metrics a (near-) optimal template bank is achieved by the A_n* lattice, which is the best lattice-covering in dimensions n <= 5, and is close to the best covering known for dimensions n <= 16. Generally this provides a substantially more efficient covering than the simpler hyper-cubic lattice. We present an algorithm for generating lattice template banks for constant-coefficient metrics and we illustrate its implementation by generating A_n* template banks in n=2,3,4 dimensions.
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.
