The effect of small inter-pulsar distance variations in stochastic gravitational wave background searches with Pulsar Timing Arrays

One of the primary objectives for Pulsar Timing Arrays (PTAs) is to detect a stochastic background generated by the incoherent superposition of gravitational waves (GWs), in particular from the cosmic population of supermassive black hole binaries. Current stochastic background searches assume that pulsars in a PTA are separated from each other and the Earth by many GW wavelengths. As more millisecond pulsars are discovered and added to PTAs, some may be separated by only a few radiation wavelengths or less, resulting in correlated GW phase changes between close pulsars in the array. Here we investigate how PTA overlap reduction functions (ORFs), up to quadrupole order, are affected by these additional correlated phase changes, and how they are in turn affected by relaxing the assumption that all pulsars are equidistant from the solar system barycenter. We find that in the low frequency GW background limit of $f\sim10^{-9}$~Hz, and for pulsars at varying distances from the Earth, that these additional correlations only affect the ORFs by a few percent for pulsar pairs at large angular separations, as expected. However when nearby (order 100 pc) pulsars are separated by less than a few degrees, the correlated phase changes can introduce variations of a few tens of percent in the magnitude of the isotropic ORF, and much larger fractional differences in the anisotropic ORFs-- up to 188 in the $m=0$, $l=2$ ORF for equidistant pulsars separated by 3 degrees. In fact, the magnitude of most of the anisotropic ORFs is largest at small, but non-zero, pulsar separations. Finally, we write down a small angle approximation for the correlated phase changes which can easily be implemented in search pipelines, and for completeness, examine the behavior of the ORFs for pulsars which lie at a radiation wavelength from the Earth.Comment: 16 pages, 8 figures, submitted to PR

Characterising gravitational wave stochastic background anisotropy with Pulsar Timing Arrays

Detecting a stochastic gravitational wave background, particularly radiation from individually unresolvable super-massive black hole binary systems, is one of the primary targets for Pulsar Timing Arrays. Increasingly more stringent upper limits are being set on these signals under the assumption that the background radiation is isotropic. However, some level of anisotropy may be present and the characterisation of the power at different angular scales carries important information. We show that the standard analysis for isotropic backgrounds can be generalised in a conceptually straightforward way to the case of generic anisotropic background radiation by decomposing the angular distribution of the gravitational wave power on the sky into multipole moments. We introduce the concept of generalised overlap reduction functions which characterise the effect of the anisotropy multipoles on the correlation of the timing residuals from the pulsars timed by a Pulsar Timing Array. In a search for a signal characterised by a generic anisotropy, the generalised overlap reduction functions play the role of the so-called Hellings and Downs curve used for isotropic radiation. We compute the generalised overlap reduction functions for a generic level of anisotropy and Pulsar Timing Array configuration. We also provide an order of magnitude estimate of the level of anisotropy that can be expected in the background generated by super-massive black hole binary systems.Comment: 12 pages plus 5 page Appendix. Accepted to PR

Effect of small interpulsar distances in stochastic gravitational wave background searches with pulsar timing arrays

One of the primary objectives for pulsar timing arrays (PTAs) is to detect a stochastic background generated by the incoherent superposition of gravitational waves (GWs), in particular from the cosmic population of supermassive black hole binaries. Current stochastic background searches assume that pulsars in a PTA are separated from each other and the Earth by many GW wavelengths. As more millisecond pulsars are discovered and added to PTAs, some may be separated by only a few radiation wavelengths or less, resulting in correlated GW phase changes between close pulsars in the array. Here we investigate how PTA overlap reduction functions (ORFs), up to quadrupole order, are affected by these additional correlated phase changes, and how they are in turn affected by relaxing the assumption that all pulsars are equidistant from the solar system barycenter. We find that in the low-frequency GW background limit of f ∼ 10^(−9) Hz, and for pulsars at varying distances from the Earth, these additional correlations only affect the ORFs by a few percent for pulsar pairs at large angular separations, as expected. However, when nearby (order 100 pc) pulsars are separated by less than a few degrees, the correlated phase changes can introduce variations of a few tens of percent in the magnitude of the isotropic ORF, and much larger fractional differences in the anisotropic ORFs—up to 188 in the m = 0, l = 2 ORF for equidistant pulsars separated by 3°. In fact, the magnitude of most of the anisotropic ORFs is largest at small, but nonzero, pulsar separations. Finally, we write down a small angle approximation for the correlated phase changes which can easily be implemented in search pipelines, and for completeness, examine the behavior of the ORFs for pulsars which lie at a radiation wavelength from the Earth

The First Two Years of Electromagnetic Follow-Up with Advanced LIGO and Virgo

We anticipate the first direct detections of gravitational waves (GWs) with Advanced LIGO and Virgo later this decade. Though this groundbreaking technical achievement will be its own reward, a still greater prize could be observations of compact binary mergers in both gravitational and electromagnetic channels simultaneously. During Advanced LIGO and Virgo's first two years of operation, 2015 through 2016, we expect the global GW detector array to improve in sensitivity and livetime and expand from two to three detectors. We model the detection rate and the sky localization accuracy for binary neutron star (BNS) mergers across this transition. We have analyzed a large, astrophysically motivated source population using real-time detection and sky localization codes and higher-latency parameter estimation codes that have been expressly built for operation in the Advanced LIGO/Virgo era. We show that for most BNS events the rapid sky localization, available about a minute after a detection, is as accurate as the full parameter estimation. We demonstrate that Advanced Virgo will play an important role in sky localization, even though it is anticipated to come online with only one-third as much sensitivity as the Advanced LIGO detectors. We find that the median 90% confidence region shrinks from ~500 square degrees in 2015 to ~200 square degrees in 2016. A few distinct scenarios for the first LIGO/Virgo detections emerge from our simulations.Comment: 17 pages, 11 figures, 5 tables. For accompanying data, see http://www.ligo.org/scientists/first2year

Parameter estimation on gravitational waves from neutron-star binaries with spinning components

Inspiraling binary neutron stars are expected to be one of the most significant sources of gravitational-wave signals for the new generation of advanced ground-based detectors. We investigate how well we could hope to measure properties of these binaries using the Advanced LIGO detectors, which began operation in September 2015. We study an astrophysically motivated population of sources (binary components with masses $1.2~\mathrm{M}_\odot$--$1.6~\mathrm{M}_\odot$ and spins of less than $0.05$) using the full LIGO analysis pipeline. While this simulated population covers the observed range of potential binary neutron-star sources, we do not exclude the possibility of sources with parameters outside these ranges; given the existing uncertainty in distributions of mass and spin, it is critical that analyses account for the full range of possible mass and spin configurations. We find that conservative prior assumptions on neutron-star mass and spin lead to average fractional uncertainties in component masses of $\sim 16\%$, with little constraint on spins (the median $90\%$ upper limit on the spin of the more massive component is $\sim 0.7$). Stronger prior constraints on neutron-star spins can further constrain mass estimates, but only marginally. However, we find that the sky position and luminosity distance for these sources are not influenced by the inclusion of spin; therefore, if LIGO detects a low-spin population of BNS sources, less computationally expensive results calculated neglecting spin will be sufficient for guiding electromagnetic follow-up.Comment: 10 pages, 9 figure

Parameter Estimation for Binary Neutron-star Coalescences with Realistic Noise during the Advanced LIGO Era

Advanced ground-based gravitational-wave (GW) detectors begin operation imminently. Their intended goal is not only to make the first direct detection of GWs, but also to make inferences about the source systems. Binary neutron-star mergers are among the most promising sources. We investigate the performance of the parameter-estimation (PE) pipeline that will be used during the first observing run of the Advanced Laser Interferometer Gravitational-wave Observatory (aLIGO) in 2015: we concentrate on the ability to reconstruct the source location on the sky, but also consider the ability to measure masses and the distance. Accurate, rapid sky localization is necessary to alert electromagnetic (EM) observatories so that they can perform follow-up searches for counterpart transient events. We consider PE accuracy in the presence of non-stationary, non-Gaussian noise. We find that the character of the noise makes negligible difference to the PE performance at a given signal-to-noise ratio. The source luminosity distance can only be poorly constrained, since the median 90% (50%) credible interval scaled with respect to the true distance is 0.85 (0.38). However, the chirp mass is well measured. Our chirp-mass estimates are subject to systematic error because we used gravitational-waveform templates without component spin to carry out inference on signals with moderate spins, but the total error is typically less than 10^(-3) M_☉. The median 90% (50%) credible region for sky localization is ~ 600 deg^2 (~150 deg^2), with 3% (30%) of detected events localized within 100 deg^2. Early aLIGO, with only two detectors, will have a sky-localization accuracy for binary neutron stars of hundreds of square degrees; this makes EM follow-up challenging, but not impossible

Superfluid neutron star dynamics, mutual friction and turbulence

This thesis investigates the role of superfluidity in neutron stars and associated phenomena. We model the internal fluid of a neutron star as a two-component system: one of charged particles and one of superfluid neutrons. We derive a set of multi-constituent hydrodynamic equations that allows for a mutual friction between the constituents. We show that when a velocity difference exists between the two constituents the momentum of each constituent is modified by an entrainment parameter. Throughout all of this work we take direction from both theoretical and experimental work on superfluid Helium. This suggests that a force due to vortex lines in the superfluid acts between the two constituents. The hydrodynamic equations are on a scale at which the effect of vortices can be averaged over. The form of the mutual friction between the two constituents depends on the configuration of the vortices. Firstly, we concentrate on an array of vortices. The mutual friction is calculated both for a straight array, and then extended to a ‘moderately’ curved array. We also investigate a turbulent model for the superfluid neutrons in which the vortices are in a tangle. To include rotation in our model we use a phenomenological approach to construct the mutual friction for a polarised tangle. The hydrodynamic equations are used to investigate how entrainment and mutual friction affect plane waves. We show that there are conditions in which the waves are unstable and discuss how this may lead to turbulence. As a first step in considering the neutron star crust we consider how oscillations in the fluid are dissipated on a boundary. As before, we concentrate on the effects of entrainment and mutual friction. Finally, we consider a simple global model of the glitch phenomenon seen in neutron stars in which the important process is a reconfiguration of the vortex array. We use this model to consider how the observational data may constrain parameters