Comparing models of the periodic variations in spin-down and beam-width for PSR B1828-11

We build a framework using tools from Bayesian data analysis to evaluate models explaining the periodic variations in spin-down and beam-width of PSR B1828-11. The available data consists of the time averaged spin-down rate, which displays a distinctive double-peaked modulation, and measurements of the beam-width. Two concepts exist in the literature that are capable of explaining these variations; we formulate predictive models from these and quantitatively compare them. The first concept is phenomenological and stipulates that the magnetosphere undergoes periodic switching between two meta-stable states as first suggested by Lyne et al. The second concept, precession, was first considered as a candidate for the modulation of B1828-11 by Stairs et al.. We quantitatively compare models built from these concepts using a Bayesian odds-ratio. Because the phenomenological switching model itself was informed by this data in the first place, it is difficult to specify appropriate parameter-space priors that can be trusted for an unbiased model comparison. Therefore we first perform a parameter estimation using the spin-down data, and then use the resulting posterior distributions as priors for model comparison on the beam-width data. We find that a precession model with a simple circular Gaussian beam geometry fails to appropriately describe the data, while allowing for a more general beam geometry provides a good fit to the data. The resulting odds between the precession model (with a general beam geometry) and the switching model are estimated as $10^{2.7 \pm 0.5}$ in favour of the precession model.Comment: 20 pages, 15 figures; removed incorrect factor of (2\pi) from equation (15), allowed for arbitrary braking index, and revised prior ranges; overall conclusions unchange

Stationary structure of relativistic superfluid neutron stars

We describe recent progress in the numerical study of the structure of rapidly rotating superfluid neutron star models in full general relativity. The superfluid neutron star is described by a model of two interpenetrating and interacting fluids, one representing the superfluid neutrons and the second consisting of the remaining charged particles (protons, electrons, muons). We consider general stationary configurations where the two fluids can have different rotation rates around a common rotation axis. The previously discovered existence of configurations with one fluid in a prolate shape is confirmed.Comment: 5 pages, 2 figures. Conference proceedings for the 26th Spanish Relativity Meeting (ERE 2002), Menorca, Spain, 22-24 Sept. 200

The effect of timing noise on targeted and narrow-band coherent searches for continuous gravitational waves from pulsars

Most searches for continuous gravitational-waves from pulsars use Taylor expansions in the phase to model the spin-down of neutron stars. Studies of pulsars demonstrate that their electromagnetic (EM) emissions suffer from \emph{timing noise}, small deviations in the phase from Taylor expansion models. How the mechanism producing EM emission is related to any continuous gravitational-wave (CW) emission is unknown; if they either interact or are locked in phase then the CW will also experience timing noise. Any disparity between the signal and the search template used in matched filtering methods will result in a loss of signal-to-noise ratio (SNR), referred to as `mismatch'. In this work we assume the CW suffers a similar level of timing noise to its EM counterpart. We inject and recover fake CW signals, which include timing noise generated from observational data on the Crab pulsar. Measuring the mismatch over durations of order $\sim 10$ months, the effect is for the most part found to be small. This suggests recent so-called `narrow-band' searches which placed upper limits on the signals from the Crab and Vela pulsars will not be significantly affected. At a fixed observation time, we find the mismatch depends upon the observation epoch. Considering the averaged mismatch as a function of observation time, we find that it increases as a power law with time, and so may become relevant in long baseline searches.Comment: 9 pages, 5 figure

Global parameter-space correlations of coherent searches for continuous gravitational waves

The space of phase-parameters (sky-position, frequency, spindowns) of a coherent matched-filtering search for continuous gravitational waves from isolated neutron stars shows strong global correlations (``circles in the sky''). In the local limit this can be analysed in terms of a parameter-space metric, but the global properties are less well studied. In this work we report on our recent progress in understanding these global correlations analytically for short to intermediate (less than a month, say) observation times and neglecting spindowns. The location of these correlation-circles in parameter-space is found to be determined mostly by the orbital velocity of the earth, while the spin-motion of the detector and the antenna-patterns only contribute significantly to the amplitude of the detection statistic along these circles.Comment: 10 pages, 6 figures; contribution to GWDAW9, submitted to CQ

Improved all-sky search method for continuous gravitational waves from unknown neutron stars in binary systems

Continuous gravitational waves from spinning deformed neutron stars have not been detected yet, and are one of the most promising signals for future detection. All-sky searches for continuous gravitational waves from unknown neutron stars in binary systems are the most computationally challenging search type. Consequently, very few search algorithms and implementations exist for these sources, and only a handful of such searches have been performed so far. In this paper, we present a new all-sky binary search method, BinarySkyHou$\mathcal{F}$, which extends and improves upon the earlier BinarySkyHough method, and which was the basis for a recent search (Covas et al. [1]). We compare the sensitivity and computational cost to previous methods, showing that it is both more sensitive and computationally efficient, which allows for broader and more sensitive searches. <br

Improved short-segment detection statistic for continuous gravitational waves

Continuous gravitational waves represent one of the long-sought types of signals that have yet to be detected. Due to their small amplitude, long observational datasets (months-years) have to be analyzed together, thereby vastly increasing the computational cost of these searches. All-sky searches face the most severe computational obstacles, especially searches for sources in unknown binary systems, which need to break the data into very short segments in order to be computationally feasible. In this paper, we present a new detection statistic that improves sensitivity by up to 19% compared to the standard $\mathcal{F}$-statistic for segments shorter than a few hours

F-statistic search for white-dwarf binaries in the first Mock LISA Data Challenge

The F-statistic is an optimal detection statistic for continuous gravitational waves, i.e. long-duration (quasi-)monochromatic signals with slowlyvarying 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 Fstatistic 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 ¯rst Mock LISA Data Challenge. Somewhat surprisingly, the LWL is found to be su±cient { 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

Searching for Galactic White Dwarf Binaries in Mock LISA Data using an F-Statistic Template Bank

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.Comment: 29 pages, 11 figures; submitted to Classical and Quantum Gravit