Comparison of Gravitational Wave Detector Network Sky Localization Approximations

Gravitational waves emitted during compact binary coalescences are a promising source for gravitational-wave detector networks. The accuracy with which the location of the source on the sky can be inferred from gravitational wave data is a limiting factor for several potential scientific goals of gravitational-wave astronomy, including multi-messenger observations. Various methods have been used to estimate the ability of a proposed network to localize sources. Here we compare two techniques for predicting the uncertainty of sky localization -- timing triangulation and the Fisher information matrix approximations -- with Bayesian inference on the full, coherent data set. We find that timing triangulation alone tends to over-estimate the uncertainty in sky localization by a median factor of $4$ for a set of signals from non-spinning compact object binaries ranging up to a total mass of $20 M_\odot$, and the over-estimation increases with the mass of the system. We find that average predictions can be brought to better agreement by the inclusion of phase consistency information in timing-triangulation techniques. However, even after corrections, these techniques can yield significantly different results to the full analysis on specific mock signals. Thus, while the approximate techniques may be useful in providing rapid, large scale estimates of network localization capability, the fully coherent Bayesian analysis gives more robust results for individual signals, particularly in the presence of detector noise.Comment: 11 pages, 7 Figure

Observing the dynamics of super-massive black hole binaries with Pulsar Timing Arrays

Pulsar Timing Arrays are a prime tool to study unexplored astrophysical regimes with gravitational waves. Here we show that the detection of gravitational radiation from individually resolvable super-massive black hole binary systems can yield direct information about the masses and spins of the black holes, provided that the gravitational-wave induced timing fluctuations both at the pulsar and at the Earth are detected. This in turn provides a map of the non-linear dynamics of the gravitational field and a new avenue to tackle open problems in astrophysics connected to the formation and evolution of super-massive black holes. We discuss the potential, the challenges and the limitations of these observations.Comment: 5 pages, 1 figur

Early Advanced LIGO binary neutron-star sky localization and parameter estimation

2015 will see the first observations of Advanced LIGO and the start of the gravitational-wave (GW) advanced-detector era. One of the most promising sources for ground-based GW detectors are binary neutron-star (BNS) coalescences. In order to use any detections for astrophysics, we must understand the capabilities of our parameter-estimation analysis. By simulating the GWs from an astrophysically motivated population of BNSs, we examine the accuracy of parameter inferences in the early advanced-detector era. We find that sky location, which is important for electromagnetic follow-up, can be determined rapidly (~5 s), but that sky areas may be hundreds of square degrees. The degeneracy between component mass and spin means there is significant uncertainty for measurements of the individual masses and spins; however, the chirp mass is well measured (typically better than 0.1%).Comment: 4 pages, 2 figures. Published in the proceedings of Amaldi 1

Towards a generic test of the strong field dynamics of general relativity using compact binary coalescence: Further investigations

In this paper we elaborate on earlier work by the same authors in which a novel Bayesian inference framework for testing the strong-field dynamics of General Relativity using coalescing compact binaries was proposed. Unlike methods that were used previously, our technique addresses the question whether one or more 'testing coefficients' (e.g. in the phase) parameterizing deviations from GR are non-zero, rather than all of them differing from zero at the same time. The framework is well-adapted to a scenario where most sources have low signal-to-noise ratio, and information from multiple sources as seen in multiple detectors can readily be combined. In our previous work, we conjectured that this framework can detect generic deviations from GR that can in principle not be accomodated by our model waveforms, on condition that the change in phase near frequencies where the detectors are the most sensitive is comparable to that induced by simple shifts in the lower-order phase coefficients of more than a few percent ($\sim 5$ radians at 150 Hz). To further support this claim, we perform additional numerical experiments in Gaussian and stationary noise according to the expected Advanced LIGO/Virgo noise curves, and coherently injecting signals into the network whose phasing differs structurally from the predictions of GR, but with the magnitude of the deviation still being small. We find that even then, a violation of GR can be established with good confidence.Comment: 15 pages, 7 figures, Amaldi 9 proceeding

The Effect of Quantized Magnetic Flux Lines on the Dynamics of Superfluid Neutron Star Cores

We investigate dynamical coupling timescales of a neutron star's superfluid core, taking into account the interactions of quantized neutron vortices with quantized flux lines of the proton superconductor in addition to the previously considered scattering of the charged components against the spontaneous magnetization of the neutron vortex line. We compare the cases where vortex motion is constrained in different ways by the array of magnetic flux tubes associated with superconducting protons. This includes absolute pinning to and creep across a uniform array of flux lines. The effect of a toroidal arrangement of flux lines is also considered. The inclusion of a uniform array of flux tubes in the neutron star core significantly decreases the timescale of coupling between the neutron and proton fluid constituents in all cases. For the toroidal component, creep response similar to that of the inner crust superfluid is possible.Comment: Submitted to MNRA

Enhancing gravitational wave astronomy with galaxy catalogues

Joint gravitational wave (GW) and electromagnetic (EM) observations, as a key research direction in multi-messenger astronomy, will provide deep insight into the astrophysics of a vast range of astronomical phenomena. Uncertainties in the source sky location estimate from gravitational wave observations mean follow-up observatories must scan large portions of the sky for a potential companion signal. A general frame of joint GW-EM observations is presented by a multi-messenger observational triangle. Using a Bayesian approach to multi-messenger astronomy, we investigate the use of galaxy catalogue and host galaxy information to reduce the sky region over which follow-up observatories must scan, as well as study its use for improving the inclination angle estimates for coalescing binary compact objects. We demonstrate our method using a simulated neutron stars inspiral signal injected into simulated Advanced detectors noise and estimate the injected signal sky location and inclination angle using the Gravitational Wave Galaxy Catalogue. In this case study, the top three candidates in rank have $72\%$, $15\%$ and $8\%$ posterior probability of being the host galaxy, receptively. The standard deviation of cosine inclination angle (0.001) of the neutron stars binary using gravitational wave-galaxy information is much smaller than that (0.02) using only gravitational wave posterior samples.Comment: Proceedings of the Sant Cugat Forum on Astrophysics. 2014 Session on 'Gravitational Wave Astrophysics

Modelling pulsar glitches with realistic pinning forces: a hydrodynamical approach

Although pulsars are one of the most stable clocks in the universe, many of them are observed to 'glitch', i.e. to suddenly increase their spin frequency (\nu) with fractional increases that range from \Delta\nu/\nu \approx 10^{-11} to 10^{-5}. In this paper we focus on the 'giant' glitches, i.e. glitches with fractional increases in the spin rate of the order of \Delta\nu/{\nu} \approx 10^{-6}, that are observed in a sub class of pulsars including the Vela. We show that giant glitches can be modelled with a two-fluid hydrodynamical approach. The model is based on the formalism for superfluid neutron stars of Andersson and Comer (2006) and on the realistic pinning forces of Grill and Pizzochero (2011). We show that all stages of Vela glitches, from the rise to the post-glitch relaxation, can be reproduced with a set of physically reasonable parameters and that the sizes and waiting times between giant glitches in other pulsars are also consistent with our model.Comment: submitted to MNRA

Towards a generic test of the strong field dynamics of general relativity using compact binary coalescence

Coalescences of binary neutron stars and/or black holes are amongst the most likely gravitational-wave signals to be observed in ground based interferometric detectors. Apart from the astrophysical importance of their detection, they will also provide us with our very first empirical access to the genuinely strong-field dynamics of General Relativity (GR). We present a new framework based on Bayesian model selection aimed at detecting deviations from GR, subject to the constraints of the Advanced Virgo and LIGO detectors. The method tests the consistency of coefficients appearing in the waveform with the predictions made by GR, without relying on any specific alternative theory of gravity. The framework is suitable for low signal-to-noise ratio events through the construction of multiple subtests, most of which involve only a limited number of coefficients. It also naturally allows for the combination of information from multiple sources to increase one's confidence in GR or a violation thereof. We expect it to be capable of finding a wide range of possible deviations from GR, including ones which in principle cannot be accommodated by the model waveforms, on condition that the induced change in phase at frequencies where the detectors are the most sensitive is comparable to the effect of a few percent change in one or more of the low-order post-Newtonian phase coefficients. In principle the framework can be used with any GR waveform approximant, with arbitrary parameterized deformations, to serve as model waveforms. In order to illustrate the workings of the method, we perform a range of numerical experiments in which simulated gravitational waves modeled in the restricted post-Newtonian, stationary phase approximation are added to Gaussian and stationary noise that follows the expected Advanced LIGO/Virgo noise curves.Comment: 26 pages, 23 figures, Accepted by PR

Entropy entrainment and dissipation in superfluid Helium

Building on a general variational framework for multi-fluid dynamics, we discuss finite temperature effects in superfluids. The main aim is to provide insight into the modelling of more complex finite temperature superfluid systems, like the mixed neutron superfluid/proton superconductor that is expected in the outer core of a neutron star. Our final results can also (to a certain extent) be used to describe colour-flavour locked quark superconductors that may be present at the extreme densities in the deep neutron star core. As a demonstration of the validity of the model, which is based on treating the excitations in the system as a massless ``entropy'' fluid, we show that it is formally equivalent to the traditional two-fluid approach for superfluid Helium. In particular, we highlight the fact that the entropy entrainment encodes the ``normal fluid density'' of the traditional approach. We also show how the superfluid constraint of irrotationality reduces the number of dissipation coefficients in the system. This analysis provides insight into the more general problem when vortices are present in the superfluid, and we discuss how the so-called mutual friction force can be accounted for in our framework. The end product is a hydrodynamic formalism for finite temperature effects in a single superfluid condensate. This framework can readily be extended to more complex situations.Comment: revised version, clarifies points regarding entrainment in different context

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