Source author record

Matthew Pitkin

Matthew Pitkin appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2026arXiv

Monitoring of Continuous-Wave Hardware Injections in LIGO Interferometers during the O4 Observing Run

Although there have now been hundreds of transient gravitational-wave detections of merging compact stars by the LIGO-Virgo-KAGRA (LVK) detector network, no continuous-wave (CW) signals have yet been discovered. To ensure that such signals, expected to be exceedingly weak, can be detected in the ongoing O4 observing run by coherent integration over years, simulated waveforms ('hardware injections') are injected directly into the LIGO data by continuously modulating the positions of the interferometer mirrors so as to mimic nearly sinusoidal signals from fast-spinning galactic neutron stars. A set of 18 such simulated CW sources are injected with signal frequencies spanning much of the LIGO detection band and with varying sky locations. By verifying the successful recovery of the simulated signals, including preservation of absolute phase over as many as 10^{11} signal cycles, we validate our understanding of detector response and end-to-end search pipelines, including data cleaning. Daily and weekly monitoring of the signal reconstruction is meant to catch any unexpected sudden changes in interferometer response, to verify that signal-to-noise ratio increases as expected and to verify that simulated source parameters are recovered correctly. We describe three methods of monitoring: 1) a highly templated matched filter to extract signal amplitude and phase precisely; 2) a frequentist Fstatistic evaluation that marginalizes over amplitude, phase and orientation of the star; and 3) a Bayesian reconstruction of the source parameters together with noise characterization. Results from each method are shown, with emphasis on the new templated method, which yields precise measurement of the critical phase offset parameter and therefore validates understanding of absolute timing delays in the detector response and data stream.

preprint2022arXiv

Nested sampling for physical scientists

We review Skilling's nested sampling (NS) algorithm for Bayesian inference and more broadly multi-dimensional integration. After recapitulating the principles of NS, we survey developments in implementing efficient NS algorithms in practice in high-dimensions, including methods for sampling from the so-called constrained prior. We outline the ways in which NS may be applied and describe the application of NS in three scientific fields in which the algorithm has proved to be useful: cosmology, gravitational-wave astronomy, and materials science. We close by making recommendations for best practice when using NS and by summarizing potential limitations and optimizations of NS.

preprint2020arXiv

Establishing the significance of continuous gravitational-wave detections from known pulsars

We present a method for assigning a statistical significance to detection candidates in targeted searches for continuous gravitational waves from known pulsars, without assuming the detector noise is Gaussian and stationary. We take advantage of the expected Doppler phase modulation of the signal induced by Earth's orbital motion, as well as the amplitude modulation induced by Earth's spin, to effectively blind the search to real astrophysical signals from a given location in the sky. We use this "sky-shifting" to produce a large number of noise-only data realizations to empirically estimate the background of a search and assign detection significances, in a similar fashion to the use of timeslides in searches for compact binaries. We demonstrate the potential of this approach by means of simulated signals, as well as hardware injections into real detector data. In a study of simulated signals in non-Gaussian noise, we find that our method outperforms another common strategy for evaluating detection significance. We thus demonstrate that this and similar techniques have the potential to enable a first confident detection of continuous gravitational waves.

preprint2016arXiv

A comparison of methods for the detection of gravitational waves from unknown neutron stars

Rapidly rotating neutron stars are promising sources of continuous gravitational wave radiation for the LIGO and Virgo interferometers. The majority of neutron stars in our galaxy have not been identified with electromagnetic observations. All-sky searches for isolated neutron stars offer the potential to detect gravitational waves from these unidentified sources. The parameter space of these blind all-sky searches, which also cover a large range of frequencies and frequency derivatives, presents a significant computational challenge. Different methods have been designed to perform these searches within acceptable computational limits. Here we describe the first benchmark in a project to compare the search methods currently available for the detection of unknown isolated neutron stars. We employ a mock data challenge to compare the ability of each search method to recover signals simulated assuming a standard signal model. We find similar performance among the short duration search methods, while the long duration search method achieves up to a factor of two higher sensitivity. We find the absence of second derivative frequency in the search parameter space does not degrade search sensivity for signals with physically plausible second derivative frequencies. We also report on the parameter estimation accuracy of each search method, and the stability of the sensitivity in frequency, frequency derivative and in the presence of detector noise.

preprint2016arXiv

A targeted spectral interpolation algorithm for the detection of continuous gravitational waves

We present an improved method of targeting continuous gravitational-wave signals in data from the LIGO and Virgo detectors with a higher efficiency than the time-domain Bayesian pipeline used in many previous searches. Our spectral interpolation algorithm, SplInter, removes the intrinsic phase evolution of the signal from source rotation and relative detector motion. We do this in the frequency domain and generate a time series containing only variations in the signal due to the antenna pattern. Although less flexible than the classic heterodyne approach, SplInter allows for rapid analysis of putative signals from isolated (and some binary) pulsars, and efficient follow-up searches for candidate signals generated by other search methods. The computational saving over the heterodyne approach can be many orders of magnitude, up to a factor of around fifty thousand in some cases, with a minimal impact on overall sensitivity for most targets.

preprint2015arXiv

Detecting Beyond-Einstein Polarizations of Continuous Gravitational Waves

The direct detection of gravitational waves with the next generation detectors, like Advanced LIGO, provides the opportunity to measure deviations from the predictions of General Relativity. One such departure would be the existence of alternative polarizations. To measure these, we study a single detector measurement of a continuous gravitational wave from a triaxial pulsar source. We develop methods to detect signals of any polarization content and distinguish between them in a model independent way. We present LIGO S5 sensitivity estimates for 115 pulsars.

preprint2015arXiv

Robust parameter estimation for compact binaries with ground-based gravitational-wave observations using the LALInference software library

The Advanced LIGO and Advanced Virgo gravitational wave (GW) detectors will begin operation in the coming years, with compact binary coalescence events a likely source for the first detections. The gravitational waveforms emitted directly encode information about the sources, including the masses and spins of the compact objects. Recovering the physical parameters of the sources from the GW observations is a key analysis task. This work describes the LALInference software library for Bayesian parameter estimation of compact binary signals, which builds on several previous methods to provide a well-tested toolkit which has already been used for several studies. We show that our implementation is able to correctly recover the parameters of compact binary signals from simulated data from the advanced GW detectors. We demonstrate this with a detailed comparison on three compact binary systems: a binary neutron star, a neutron star black hole binary and a binary black hole, where we show a cross-comparison of results obtained using three independent sampling algorithms. These systems were analysed with non-spinning, aligned spin and generic spin configurations respectively, showing that consistent results can be obtained even with the full 15-dimensional parameter space of the generic spin configurations. We also demonstrate statistically that the Bayesian credible intervals we recover correspond to frequentist confidence intervals under correct prior assumptions by analysing a set of 100 signals drawn from the prior. We discuss the computational cost of these algorithms, and describe the general and problem-specific sampling techniques we have used to improve the efficiency of sampling the compact binary coalescence parameter space.

preprint2014arXiv

Glitch or Anti-Glitch: a Bayesian View

The sudden spin-down in the rotation of magnetar 1E 2259+586 observed by Archibald et al. (2013) was a rare event. However this particular event, referred to as an anti-glitch, was followed by another event which Archibald et al. (2013) suggested could either be a conventional glitch or another anti-glitch. Although there is no accompanied radiation activity or pulse profile change, there is decisive evidence for the existence of the second timing event, judging from the timing data. We apply Bayesian Model Selection to quantitatively determine which of these possibilities better explains the observed data. We show that the observed data strongly supports the presence of two successive anti-glitches with a Bayes Factor, often called the odds ratio, greater than 40. Furthermore, we show that the second anti-gtlich has an associated frequency change $Δν$ of $-8.2 \times 10^{-8}$ Hz. We discuss the implications of these results for possible physical mechanisms behind this anti-glitch.

preprint2013arXiv

Reconstructing the sky location of gravitational-wave detected compact binary systems: methodology for testing and comparison

The problem of reconstructing the sky position of compact binary coalescences detected via gravitational waves is a central one for future observations with the ground-based network of gravitational-wave laser interferometers, such as Advanced LIGO and Advanced Virgo. Different techniques for sky localisation have been independently developed. They can be divided in two broad categories: fully coherent Bayesian techniques, which are high-latency and aimed at in-depth studies of all the parameters of a source, including sky position, and "triangulation-based" techniques, which exploit the data products from the search stage of the analysis to provide an almost real-time approximation of the posterior probability density function of the sky location of a detection candidate. These techniques have previously been applied to data collected during the last science runs of gravitational-wave detectors operating in the so-called initial configuration. Here, we develop and analyse methods for assessing the self-consistency of parameter estimation methods and carrying out fair comparisons between different algorithms, addressing issues of efficiency and optimality. These methods are general, and can be applied to parameter estimation problems other than sky localisation. We apply these methods to two existing sky localisation techniques representing the two above-mentioned categories, using a set of simulated inspiral-only signals from compact binary systems with total mass $\le 20\,M_\odot$ and non-spinning components. We compare the relative advantages and costs of the two techniques and show that sky location uncertainties are on average a factor $\approx 20$ smaller for fully coherent techniques than for the specific variant of the "triangulation-based" technique used during the last science runs, at the expense of a factor $\approx 1000$ longer processing time.

preprint2012arXiv

Extending gravitational wave burst searches with pulsar timing arrays

Pulsar timing arrays (PTAs) are being used to search for very low frequency gravitational waves. A gravitational wave signal appears in pulsar timing residuals through two components: one independent of and one dependent on the pulsar's distance, called the 'Earth term' (ET) and 'pulsar term' (PT), respectively. The signal of a burst (or transient) gravitational wave source in pulsars' residuals will in general have the Earth and pulsar terms separated by times of the order of the time of flight from the pulsar to the Earth. Therefore, both terms are not observable over a realistic observation span, but the ETs observed in many pulsars should be correlated. We show that pairs (or more) of pulsars can be aligned in such a way that the PTs caused by a source at certain sky locations can arrive at Earth within a time window short enough to be captured during a realistic observation span. We find that for the pulsars within the International Pulsar Timing Array (IPTA) ~67 per cent of the sky produces such alignments for pulsars terms separated by less than 10 years. We compare estimates of the source event rate that would be required to observe one signal in the IPTA if searching for the correlated ETs, or in searching via the PTs, and find that event rates would need to be about two orders of magnitude higher to observe an event with the PTs than the ETs. We also find that an array of hundreds of thousands of pulsars would be required to achieve similar numbers of observable events in PT or ET searches. This disfavours PTs being used for all-sky searches, but they could potentially be used target specific sources and be complementary to ET only searches.

preprint2011arXiv

Gravitational Wave Detection by Interferometry (Ground and Space)

Significant progress has been made in recent years on the development of gravitational wave detectors. Sources such as coalescing compact binary systems, neutron stars in low-mass X-ray binaries, stellar collapses and pulsars are all possible candidates for detection. The most promising design of gravitational wave detector uses test masses a long distance apart and freely suspended as pendulums on Earth or in drag-free craft in space. The main theme of this review is a discussion of the mechanical and optical principles used in the various long baseline systems in operation around the world - LIGO (USA), Virgo (Italy/France), TAMA300 and LCGT (Japan), and GEO600 (Germany/U.K.) - and in LISA, a proposed space-borne interferometer. A review of recent science runs from the current generation of ground-based detectors will be discussed, in addition to highlighting the astrophysical results gained thus far. Looking to the future, the major upgrades to LIGO (Advanced LIGO), Virgo (Advanced Virgo), LCGT and GEO600 (GEO-HF) will be completed over the coming years, which will create a network of detectors with significantly improved sensitivity required to detect gravitational waves. Beyond this, the concept and design of possible future "third generation" gravitational wave detectors, such as the Einstein Telescope (ET), will be discussed.

preprint2011arXiv

Prospects of observing continuous gravitational waves from known pulsars

Several past searches for gravitational waves from a selection of known pulsars have been performed with data from the science runs of the Laser Inferometer Gravitational-wave Observatory (LIGO) gravitational wave detectors. So far these have lead to no detection, but upper limits on the gravitational wave amplitudes have been set. Here we study our intrinsic ability to detect, and estimate the gravitational wave amplitude for non-accreting pulsars. Using spin-down limits on emission as a guide we examine amplitudes that would be required to observe known pulsars with future detectors (Advanced LIGO, Advanced Virgo and the Einstein Telescope), assuming that they are triaxial stars emitting at precisely twice the known rotation frequency. Maximum allowed amplitudes depend on the stars' equation of state (e.g. a normal neutron star, a quark star, a hybrid star) and the theoretical mass quadrupoles that they can sustain. We study what range of quadrupoles, and therefore equations of state, would be consistent with being able to detect these sources. For globular cluster pulsars, with spin-downs masked by accelerations within the cluster, we examine what spin-down values gravitational wave observations would be able to set. For all pulsars we also alternatively examine what internal magnetic fields they would need to sustain observable ellipticities.