P. D. Lasky

Gravitational waves from binary supermassive black holes missing in pulsar observations

Placing bounds on gravitational wave detection Gravitational waves are expected to be generated by the interaction of the massive bodies in black-hole binary systems. As gravitational waves distort spacetime, it should be possible to verify their existence as they interfere with the pulses emitted by millisecond pulsars. However, after monitoring 24 pulsars with the Parkes radio telescope for 12 years, Shannon et al. found no detectable variation in pulsar records. This nondetection result indicates that a new detection strategy for gravitational waves is needed. Science, this issue p. 1522 A lack of observed variations in the timing of pulsars places constraints on the detection of gravitational waves. Gravitational waves are expected to be radiated by supermassive black hole binaries formed during galaxy mergers. A stochastic superposition of gravitational waves from all such binary systems would modulate the arrival times of pulses from radio pulsars. Using observations of millisecond pulsars obtained with the Parkes radio telescope, we constrained the characteristic amplitude of this background, Ac,yr, to be <1.0 × 10−15 with 95% confidence. This limit excludes predicted ranges for Ac,yr from current models with 91 to 99.7% probability. We conclude that binary evolution is either stalled or dramatically accelerated by galactic-center environments and that higher-cadence and shorter-wavelength observations would be more sensitive to gravitational waves.

Uniting old stellar systems: from globular clusters to giant ellipticals

Elliptical galaxies and globular clusters (GCs) have traditionally been regarded as physically distinct entities due to their discontinuous distribution in key scaling diagrams involving size, luminosity and velocity dispersion. Recently this distinctness has been challenged by the discovery of stellar systems with mass intermediate between those of GCs and dwarf ellipticals (such as ultracompact dwarfs and dwarf galaxy transition objects). Here we examine the relationship between the virial and stellar mass for a range of old stellar systems, from GCs to giant ellipticals, and including such intermediate-mass objects (IMOs). Improvements on previous work in this area include the use of (i) near-infrared magnitudes from the Two Micron All Sky Survey (2MASS), (ii) aperture corrections to velocity dispersions, (iii) homogeneous half-light radii and (iv) accounting for the effects of non-homology in galaxies. We find a virial-to-stellar mass relation that ranges from ∼10 4 M ⊙ systems (GCs) to ∼10 12 M⊙ systems (elliptical galaxies). The lack of measured velocity dispersions for dwarf ellipticals with - 16 > M K > - 18 (∼ 10 8 M ⊙ ) currently inhibits our ability to determine how, or indeed if, these galaxies connect continuously with GCs in terms of their virial-to-stellar mass ratios. We find elliptical galaxies to have roughly equal fractions of dark and stellar matter within a virial radius; only in the most massive (greater than 10 12 M ⊙ ) ellipticals does dark matter dominate the virial mass. Although the IMOs reveal slightly higher virial-to-stellar mass ratios than lower mass GCs, this may simply reflect our limited understanding of their initial mass function (and hence their stellar mass-to-light ratios) or structural properties. We argue that most of these IMOs have similar properties to massive GCs, i.e. IMOs are essentially massive star clusters. Only the dwarf spheroidal galaxies exhibit behaviour notably distinct from the other stellar systems examined here, i.e. they display a strongly increasing virial-to-stellar mass ratio (equivalent to higher dark matter fractions) with decreasing stellar mass. The data used in this study are available in electronic format.

Detecting Gravitational-Wave Memory with LIGO: Implications of GW150914

It may soon be possible for Advanced LIGO to detect hundreds of binary black hole mergers per year. We show how the accumulation of many such measurements will allow for the detection of gravitational-wave memory: a permanent displacement of spacetime that comes from strong-field, general relativistic effects. We estimate that Advanced LIGO operating at design sensitivity may be able to make a signal-to-noise ratio 3 (5) detection of memory with ∼35 (90) events with masses and distance similar to GW150914. We highlight the importance of incorporating higher-order gravitational-wave modes for parameter estimation of binary black hole mergers, and describe how our methods can also be used to detect higher-order modes themselves before Advanced LIGO reaches design sensitivity.

Gravitational-wave cosmology across 29 decades in frequency

Quantum fluctuations of the gravitational field in the early Universe, amplified by inflation, produce a primordial gravitational-wave background across a broad frequency band. We derive constraints on the spectrum of this gravitational radiation, and hence on theories of the early Universe, by combining experiments that cover 29 orders of magnitude in frequency. These include Planck observations of cosmic microwave background temperature and polarization power spectra and lensing, together with baryon acoustic oscillations and big bang nucleosynthesis measurements, as well as new pulsar timing array and ground-based interferometer limits. While individual experiments constrain the gravitational-wave energy density in specific frequency bands, the combination of experiments allows us to constrain cosmological parameters, including the inflationary spectral index n_t and the tensor-to-scalar ratio r. Results from individual experiments include the most stringent nanohertz limit of the primordial background to date from the Parkes Pulsar Timing Array, Ω_(GW)(f) < 2.3 × 10^(−10). Observations of the cosmic microwave background alone limit the gravitational-wave spectral index at 95% confidence to n_t ≲ 5 for a tensor-to-scalar ratio of r = 0.11. However, the combination of all the above experiments limits n_t < 0.36. Future Advanced LIGO observations are expected to further constrain n_t < 0.34 by 2020. When cosmic microwave background experiments detect a nonzero r, our results will imply even more stringent constraints on n_t and, hence, theories of the early Universe.

Observationally constraining gravitational wave emission from short gamma-ray burst remnants

Observations of short gamma-ray bursts indicate ongoing energy injection following the prompt emission, with the most likely candidate being the birth of a rapidly rotating, highly magnetised neutron star. We utilise X-ray observations of the burst remnant to constrain properties of the nascent neutron star, including its magnetic field-induced ellipticity and the saturation amplitude of various oscillation modes. Moreover, we derive strict upper limits on the gravitational wave emission from these objects by looking only at the X-ray light curve, showing the burst remnants are unlikely to be detected in the near future using ground-based gravitational wave interferometers such as Advanced LIGO.

From spin noise to systematics: stochastic processes in the first International Pulsar Timing Array data release

We analyse the stochastic properties of the 49 pulsars that comprise the first International Pulsar Timing Array (IPTA) data release. We use Bayesian methodology, performing model selection to determine the optimal description of the stochastic signals present in each pulsar. In addition to spin-noise and dispersion-measure (DM) variations, these models can include timing noise unique to a single observing system, or frequency band. We show the improved radio-frequency coverage and presence of overlapping data from different observing systems in the IPTA data set enables us to separate both system and band-dependent effects with much greater efficacy than in the individual PTA data sets. For example, we show that PSR J1643−1224 has, in addition to DM variations, significant band-dependent noise that is coherent between PTAs which we interpret as coming from time-variable scattering or refraction in the ionised interstellar medium. Failing to model these different contributions appropriately can dramatically alter the astrophysical interpretation of the stochastic signals observed in the residuals. In some cases, the spectral exponent of the spin noise signal can vary from 1.6 to 4 depending upon the model, which has direct implications for the long-term sensitivity of the pulsar to a stochastic gravitational-wave (GW) background. By using a more appropriate model, however, we can greatly improve a pulsars sensitivity to GWs. For example, including system and band-dependent signals in the PSR J0437−4715 data set improves the upper limit on a fiducial GW background by ∼ 60% compared to a model that includes DM variations and spin-noise only.

Are gravitational waves from giant magnetar flares observable

Are giant flares in magnetars viable sources of gravitational radiation? Few theoretical studies have been concerned with this problem, with the small number using either highly idealized models or assuming a magnetic field orders of magnitude beyond what is supported by observations. We perform nonlinear general-relativistic magnetohydrodynamics simulations of large-scale hydromagnetic instabilities in magnetar models. We utilise these models to find gravitational wave emissions over a wide range of energies, from 10^40 to 10^47 erg. This allows us to derive a systematic relationship between the surface field strength and the gravitational wave strain, which we find to be highly nonlinear. In particular, for typical magnetar fields of a few times 10^15 G, we conclude that a direct observation of f-modes excited by global magnetic field reconfigurations is unlikely with present or near-future gravitational wave observatories, though we also discuss the possibility that modes in a low-frequency band up to 100 Hz could be sufficiently excited to be relevant for observation.

Neutron star deformation due to multipolar magnetic fields

Certain multi-wavelength observations of neutron stars, such as intermittent radio emissions from rotation-powered pulsars beyond the pair-cascade death line, the pulse profile of the magnetar SGR 1900+14 after its 1998 August 27 giant flare, and X-ray spectral features of PSR J0821 4300 and SGR 0418+5729, suggest that the magnetic fields of non-accreting neutron stars are not purely dipolar and may contain higherorder multipoles. Here, we calculate the ellipticity of a non-barotropic neutron star with (i) a quadrupole poloidal-toroidal field, and (ii) a purely poloidal field containing arbitrary multipoles, deriving the relation between the ellipticity and the multipole amplitudes. We present, as a worked example, a purely poloidal field comprising dipole, quadrupole, and octupole components. We show the correlation between field energy and ellipticity for each multipole, that the l = 4 multipole has the lowest energy, and that l = 5 has the lowest ellipticity. We show how a mixed multipolar field creates an observationally testable mismatch between the principal axes of inertia (to be inferred from gravitational wave data) and the magnetic inclination angle. Strong quadrupole and octupole components (with amplitudes � 10 2 times higher than the dipole) in SGR 0418+5729 still yield ellipticity � 10 8 , consistent with current gravitational wave upper limits. The existence of higher multipoles in fast-rotating objects (e.g., newborn magnetars) has interesting implications for the braking law and hence phase tracking during coherent gravitational wave searches.

Gravitational waves from Scorpius X-1: A comparison of search methods and prospects for detection with advanced detectors

The low-mass X-ray binary Scorpius X-1 (Sco X-1) is potentially the most luminous source of continuous gravitational-wave radiation for interferometers such as LIGO and Virgo. For low-mass X-ray binaries this radiation would be sustained by active accretion of matter from its binary companion. With the Advanced Detector Era fast approaching, work is underway to develop an array of robust tools for maximizing the science and detection potential of Sco X-1. We describe the plans and progress of a project designed to compare the numerous independent search algorithms currently available. We employ a mock-data challenge in which the search pipelines are tested for their relative proficiencies in parameter estimation, computational efficiency, robustness, and most importantly, search sensitivity. The mock-data challenge data contains an ensemble of 50 Scorpius X-1 (Sco X-1) type signals, simulated within a frequency band of 50–1500 Hz. Simulated detector noise was generated assuming the expected best strain sensitivity of Advanced LIGO [1] and Advanced VIRGO [2] (4×10^(−24)  Hz^(−1/2)). A distribution of signal amplitudes was then chosen so as to allow a useful comparison of search methodologies. A factor of 2 in strain separates the quietest detected signal, at 6.8×10^(−26) strain, from the torque-balance limit at a spin frequency of 300 Hz, although this limit could range from 1.2×10^(−25) (25 Hz) to 2.2×10^(−26) (750 Hz) depending on the unknown frequency of Sco X-1. With future improvements to the search algorithms and using advanced detector data, our expectations for probing below the theoretical torque-balance strain limit are optimistic.

Inhomogeneous cosmology with numerical relativity

We perform three-dimensional numerical relativity simulations of homogeneous and inhomogeneous expanding spacetimes, with a view towards quantifying non-linear effects from cosmological inhomogeneities. We demonstrate fourth-order convergence with errors less than one part in 10^6 in evolving a flat, dust Friedmann-Lemaitre-Roberston-Walker (FLRW) spacetime using the Einstein Toolkit within the Cactus framework. We also demonstrate agreement to within one part in 10^3 between the numerical relativity solution and the linear solution for density, velocity and metric perturbations in the Hubble flow over a factor of ~350 change in scale factor (redshift). We simulate the growth of linear perturbations into the non-linear regime, where effects such as gravitational slip and tensor perturbations appear. We therefore show that numerical relativity is a viable tool for investigating nonlinear effects in cosmology.