-
We present the results of a LIGO search for gravitational waves (GWs)
associated with GRB 051103, a short-duration hard-spectrum gamma-ray burst
(GRB) whose electromagnetically determined sky position is coincident with the
spiral galaxy M81, which is 3.6 Mpc from Earth. Possible progenitors for
short-hard GRBs include compact object mergers and soft gamma repeater (SGR)
giant flares. A merger progenitor would produce a characteristic GW signal that
should be detectable at the distance of M81, while GW emission from an SGR is
not expected to be detectable at that distance. We found no evidence of a GW
signal associated with GRB 051103. Assuming weakly beamed gamma-ray emission
with a jet semi-angle of 30 deg we exclude a binary neutron star merger in M81
as the progenitor with a confidence of 98%. Neutron star-black hole mergers are
excluded with > 99% confidence. If the event occurred in M81 our findings
support the the hypothesis that GRB 051103 was due to an SGR giant flare,
making it the most distant extragalactic magnetar observed to date.
-
We report on a search for gravitational waves from coalescing compact
binaries using LIGO and Virgo observations between July 7, 2009 and October 20,
2010. We searched for signals from binaries with total mass between 2 and 25
solar masses; this includes binary neutron stars, binary black holes, and
binaries consisting of a black hole and neutron star. The detectors were
sensitive to systems up to 40 Mpc distant for binary neutron stars, and further
for higher mass systems. No gravitational-wave signals were detected. We report
upper limits on the rate of compact binary coalescence as a function of total
mass, including the results from previous LIGO and Virgo observations. The
cumulative 90%-confidence rate upper limits of the binary coalescence of binary
neutron star, neutron star- black hole and binary black hole systems are 1.3 x
10^{-4}, 3.1 x 10^{-5} and 6.4 x 10^{-6} Mpc^{-3}yr^{-1}, respectively. These
upper limits are up to a factor 1.4 lower than previously derived limits. We
also report on results from a blind injection challenge.
-
Aims. A transient astrophysical event observed in both gravitational wave
(GW) and electromagnetic (EM) channels would yield rich scientific rewards. A
first program initiating EM follow-ups to possible transient GW events has been
developed and exercised by the LIGO and Virgo community in association with
several partners. In this paper, we describe and evaluate the methods used to
promptly identify and localize GW event candidates and to request images of
targeted sky locations.
Methods. During two observing periods (Dec 17 2009 to Jan 8 2010 and Sep 2 to
Oct 20 2010), a low-latency analysis pipeline was used to identify GW event
candidates and to reconstruct maps of possible sky locations. A catalog of
nearby galaxies and Milky Way globular clusters was used to select the most
promising sky positions to be imaged, and this directional information was
delivered to EM observatories with time lags of about thirty minutes. A Monte
Carlo simulation has been used to evaluate the low-latency GW pipeline's
ability to reconstruct source positions correctly.
Results. For signals near the detection threshold, our low-latency algorithms
often localized simulated GW burst signals to tens of square degrees, while
neutron star/neutron star inspirals and neutron star/black hole inspirals were
localized to a few hundred square degrees. Localization precision improves for
moderately stronger signals. The correct sky location of signals well above
threshold and originating from nearby galaxies may be observed with ~50% or
better probability with a few pointings of wide-field telescopes.
-
We report on an all-sky search for periodic gravitational waves in the
frequency band 50-800 Hz and with the frequency time derivative in the range of
0 through -6e-9 Hz/s. Such a signal could be produced by a nearby spinning and
slightly non-axisymmetric isolated neutron star in our galaxy. After recent
improvements in the search program that yielded a 10x increase in computational
efficiency, we have searched in two years of data collected during LIGO's fifth
science run and have obtained the most sensitive all-sky upper limits on
gravitational wave strain to date. Near 150 Hz our upper limit on worst-case
linearly polarized strain amplitude $h_0$ is 1e-24, while at the high end of
our frequency range we achieve a worst-case upper limit of 3.8e-24 for all
polarizations and sky locations. These results constitute a factor of two
improvement upon previously published data. A new detection pipeline utilizing
a Loosely Coherent algorithm was able to follow up weaker outliers, increasing
the volume of space where signals can be detected by a factor of 10, but has
not revealed any gravitational wave signals. The pipeline has been tested for
robustness with respect to deviations from the model of an isolated neutron
star, such as caused by a low-mass or long-period binary companion.
-
The gravitational-wave (GW) sky may include nearby pointlike sources as well
as astrophysical and cosmological stochastic backgrounds. Since the relative
strength and angular distribution of the many possible sources of GWs are not
well constrained, searches for GW signals must be performed in a
model-independent way. To that end we perform two directional searches for
persistent GWs using data from the LIGO S5 science run: one optimized for
pointlike sources and one for arbitrary extended sources. The latter result is
the first of its kind. Finding no evidence to support the detection of GWs, we
present 90% confidence level (CL) upper-limit maps of GW strain power with
typical values between 2-20x10^-50 strain^2 Hz^-1 and 5-35x10^-49 strain^2
Hz^-1 sr^-1 for pointlike and extended sources respectively. The limits on
pointlike sources constitute a factor of 30 improvement over the previous best
limits. We also set 90% CL limits on the narrow-band root-mean-square GW strain
from interesting targets including Sco X-1, SN1987A and the Galactic Center as
low as ~7x10^-25 in the most sensitive frequency range near 160 Hz. These
limits are the most constraining to date and constitute a factor of 5
improvement over the previous best limits.
-
We present direct upper limits on continuous gravitational wave emission from
the Vela pulsar using data from the Virgo detector's second science run. These
upper limits have been obtained using three independent methods that assume the
gravitational wave emission follows the radio timing. Two of the methods
produce frequentist upper limits for an assumed known orientation of the star's
spin axis and value of the wave polarization angle of, respectively,
$1.9\ee{-24}$ and $2.2\ee{-24}$, with 95% confidence. The third method, under
the same hypothesis, produces a Bayesian upper limit of $2.1\ee{-24}$, with 95%
degree of belief. These limits are below the indirect {\it spin-down limit} of
$3.3\ee{-24}$ for the Vela pulsar, defined by the energy loss rate inferred
from observed decrease in Vela's spin frequency, and correspond to a limit on
the star ellipticity of $\sim 10^{-3}$. Slightly less stringent results, but
still well below the spin-down limit, are obtained assuming the star's spin
axis inclination and the wave polarization angles are unknown.
-
Soft gamma repeaters (SGRs) and anomalous X-ray pulsars (AXPs) are thought to
be magnetars: neutron stars powered by extreme magnetic fields. These rare
objects are characterized by repeated and sometimes spectacular gamma-ray
bursts. The burst mechanism might involve crustal fractures and excitation of
non-radial modes which would emit gravitational waves (GWs). We present the
results of a search for GW bursts from six galactic magnetars that is sensitive
to neutron star f-modes, thought to be the most efficient GW emitting
oscillatory modes in compact stars. One of them, SGR 0501+4516, is likely ~1
kpc from Earth, an order of magnitude closer than magnetars targeted in
previous GW searches. A second, AXP 1E 1547.0-5408, gave a burst with an
estimated isotropic energy >10^{44} erg which is comparable to the giant
flares. We find no evidence of GWs associated with a sample of 1279
electromagnetic triggers from six magnetars occurring between November 2006 and
June 2009, in GW data from the LIGO, Virgo, and GEO600 detectors. Our lowest
model-dependent GW emission energy upper limits for band- and time-limited
white noise bursts in the detector sensitive band, and for f-mode ringdowns (at
1090 Hz), are 3.0x10^{44} d_1^2 erg and 1.4x10^{47} d_1^2 erg respectively,
where d_1 = d_{0501} / 1 kpc and d_{0501} is the distance to SGR 0501+4516.
These limits on GW emission from f-modes are an order of magnitude lower than
any previous, and approach the range of electromagnetic energies seen in SGR
giant flares for the first time.
-
We present the first modeled search for gravitational waves using the
complete binary black hole gravitational waveform from inspiral through the
merger and ringdown for binaries with negligible component spin. We searched
approximately 2 years of LIGO data taken between November 2005 and September
2007 for systems with component masses of 1-99 solar masses and total masses of
25-100 solar masses. We did not detect any plausible gravitational-wave signals
but we do place upper limits on the merger rate of binary black holes as a
function of the component masses in this range. We constrain the rate of
mergers for binary black hole systems with component masses between 19 and 28
solar masses and negligible spin to be no more than 2.0 per Mpc^3 per Myr at
90% confidence.
-
The physical mechanisms responsible for pulsar timing glitches are thought to
excite quasi-normal mode oscillations in their parent neutron star that couple
to gravitational wave emission. In August 2006, a timing glitch was observed in
the radio emission of PSR B0833-45, the Vela pulsar. At the time of the glitch,
the two co-located Hanford gravitational wave detectors of the Laser
Interferometer Gravitational-wave observatory (LIGO) were operational and
taking data as part of the fifth LIGO science run (S5). We present the first
direct search for the gravitational wave emission associated with oscillations
of the fundamental quadrupole mode excited by a pulsar timing glitch. No
gravitational wave detection candidate was found. We place Bayesian 90%
confidence upper limits of 6.3e-21 to 1.4e-20 on the peak intrinsic strain
amplitude of gravitational wave ring-down signals, depending on which spherical
harmonic mode is excited. The corresponding range of energy upper limits is
5.0e44 to 1.3e45 erg.
-
We apply the postquasistatic approximation, an iterative method for the
evolution of self-gravitating spheres of matter, to study the evolution of
dissipative and electrically charged distributions in General Relativity. We
evolve nonadiabatic distributions assuming an equation of state that accounts
for the anisotropy induced by the electric charge. Dissipation is described by
streaming out or diffusion approximations. We match the interior solution, in
noncomoving coordinates, with the Vaidya-Reissner-Nordstr\"om exterior
solution. Two models are considered: i) a Schwarzschild-like shell in the
diffusion limit; ii) a Schwarzschild-like interior in the free streaming limit.
These toy models tell us something about the nature of the dissipative and
electrically charged collapse. Diffusion stabilizes the gravitational collapse
producing a spherical shell whose contraction is halted in a short
characteristic hydrodynamic time. The streaming out radiation provides a more
efficient mechanism for emission of energy, redistributing the electric charge
on the whole sphere, while the distribution collapses indefinitely with a
longer hydrodynamic time scale.
-
We present a search for periodic gravitational waves from the neutron star in
the supernova remnant Cassiopeia A. The search coherently analyzes data in a
12-day interval taken from the fifth science run of the Laser Interferometer
Gravitational-Wave Observatory. It searches gravitational wave frequencies from
100 to 300 Hz, and covers a wide range of first and second frequency
derivatives appropriate for the age of the remnant and for different spin-down
mechanisms. No gravitational wave signal was detected. Within the range of
search frequencies, we set 95% confidence upper limits of 0.7--1.2e-24 on the
intrinsic gravitational wave strain, 0.4--4e-4 on the equatorial ellipticity of
the neutron star, and 0.005--0.14 on the amplitude of r-mode oscillations of
the neutron star. These direct upper limits beat indirect limits derived from
energy conservation and enter the range of theoretical predictions involving
crystalline exotic matter or runaway r-modes. This is the first gravitational
wave search to present upper limits on r-modes.
-
We apply the postquasistatic approximation to study the evolution of
spherically symmetric fluid distributions undergoing dissipation in the form of
radial heat flow. For a model which corresponds to an incompressible fluid
departing from the static equilibrium, it is not possible to go far from the
initial state after the emission of a small amount of energy. Initially
collapsing distributions of matter are not permitted. Emission of energy can be
considered as a mechanism to avoid the collapse. If the distribution collapses
initially and emits one hundredth of the initial mass only the outermost layers
evolve. For a model which corresponds to a highly compressed Fermi gas, only
the outermost shell can evolve with a shorter hydrodynamic time scale.
-
The Laser Interferometer Gravitational Wave Observatory (LIGO) is a network
of three detectors built to detect local perturbations in the space-time metric
from astrophysical sources. These detectors, two in Hanford, WA and one in
Livingston, LA, are power-recycled Fabry-Perot Michelson interferometers. In
their fifth science run (S5), between November 2005 and October 2007, these
detectors accumulated one year of triple coincident data while operating at
their designed sensitivity. In this paper, we describe the calibration of the
instruments in the S5 data set, including measurement techniques and
uncertainty estimation.
-
We summarize the sensitivity achieved by the LIGO and Virgo gravitational
wave detectors for compact binary coalescence (CBC) searches during LIGO's
fifth science run and Virgo's first science run. We present noise spectral
density curves for each of the four detectors that operated during these
science runs which are representative of the typical performance achieved by
the detectors for CBC searches. These spectra are intended for release to the
public as a summary of detector performance for CBC searches during these
science runs.
-
We apply the post-quasi--static approximation, an iterative method for the
evolution of self-gravitating spheres of matter, to study the evolution of
anisotropic non-adiabatic radiating and dissipative distributions in General
Relativity. Dissipation is described by viscosity and free-streaming radiation,
assuming an equation of state to model anisotropy induced by the shear
viscosity. We match the interior solution, in non-comoving coordinates, with
the Vaidya exterior solution. Two simple models are presented, based on the
Schwarzschild and Tolman VI solutions, in the non--adiabatic and adiabatic
limit. In both cases the eventual collapse or expansion of the distribution is
mainly controlled by the anisotropy induced by the viscosity.
-
We present an up-to-date, comprehensive summary of the rates for all types of
compact binary coalescence sources detectable by the Initial and Advanced
versions of the ground-based gravitational-wave detectors LIGO and Virgo.
Astrophysical estimates for compact-binary coalescence rates depend on a number
of assumptions and unknown model parameters, and are still uncertain. The most
confident among these estimates are the rate predictions for coalescing binary
neutron stars which are based on extrapolations from observed binary pulsars in
our Galaxy. These yield a likely coalescence rate of 100 per Myr per Milky Way
Equivalent Galaxy (MWEG), although the rate could plausibly range from 1 per
Myr per MWEG to 1000 per Myr per MWEG. We convert coalescence rates into
detection rates based on data from the LIGO S5 and Virgo VSR2 science runs and
projected sensitivities for our Advanced detectors. Using the detector
sensitivities derived from these data, we find a likely detection rate of 0.02
per year for Initial LIGO-Virgo interferometers, with a plausible range between
0.0002 and 0.2 per year. The likely binary neutron-star detection rate for the
Advanced LIGO-Virgo network increases to 40 events per year, with a range
between 0.4 and 400 per year.
-
Progenitor scenarios for short gamma-ray bursts (short GRBs) include
coalescenses of two neutron stars or a neutron star and black hole, which would
necessarily be accompanied by the emission of strong gravitational waves. We
present a search for these known gravitational-wave signatures in temporal and
directional coincidence with 22 GRBs that had sufficient gravitational-wave
data available in multiple instruments during LIGO's fifth science run, S5, and
Virgo's first science run, VSR1. We find no statistically significant
gravitational-wave candidates within a [-5, +1) s window around the trigger
time of any GRB. Using the Wilcoxon-Mann-Whitney U test, we find no evidence
for an excess of weak gravitational-wave signals in our sample of GRBs. We
exclude neutron star-black hole progenitors to a median 90% CL exclusion
distance of 6.7 Mpc.
-
(Abridged.) The mean-square current quadrupole moment associated with
vorticity fluctuations in high-Reynolds-number turbulence in a differentially
rotating neutron star is calculated analytically, as are the amplitude and
decoherence time of the resulting, stochastic gravitational wave signal. The
calculation resolves the subtle question of whether the signal is dominated by
the smallest or largest turbulent eddies: for the Kolmogorov-like power
spectrum observed in superfluid spherical Couette simulations, the wave strain
is controlled by the largest eddies, and the decoherence time approximately
equals the maximum eddy turnover time. For a neutron star with spin frequency
$\nu_s$ and Rossby number $Ro$, at a distance $d$ from Earth, the
root-mean-square wave strain reaches $h_{RMS} \approx 3\times 10^{-24} Ro^3
(\nu_s / 30 Hz)^3 (d/1 kpc)^{-1}$. A cross-correlation search can detect such a
source in principle, because the signal decoheres over the time-scale $\tau_c
\approx 10^{-3} Ro^{-1} (\nu_s / 30 Hz)^{-1} s$, which is adequately sampled by
existing long-baseline interferometers. Hence hydrodynamic turbulence imposes a
fundamental noise floor on gravitational wave observations of neutron stars,
although its polluting effect may be muted by partial decoherence in the
hectohertz band, where current continuous-wave searches are concentrated, for
the highest frequency (and hence most powerful) sources.
-
Experimental and numerical evidence is reviewed for the existence of a
Stewartson layer in spherical Couette flow at small Ekman and Rossby numbers
($\Ek \lsim 10^{-3}$, $\Ro \lsim 10^{-2}$), the relevant hydrodynamic regime in
the superfluid outer core of a neutron star. Numerical simulations of a
superfluid Stewartson layer are presented for the first time, showing how the
layer is disrupted by nonaxisymmetric instabilities. The unstable ranges of
$\Ek$ and $\Ro$ are compared with estimates of these quantities in radio
pulsars that exhibit glitches. It is found that glitching pulsars lie on the
stable side of the instability boundary, allowing differential rotation to
build up before a glitch.
-
We solve numerically for the first time the two-fluid,
Hall--Vinen--Bekarevich--Khalatnikov (HVBK) equations for a He-II-like
superfluid contained in a differentially rotating, spherical shell,
generalizing previous simulations of viscous spherical Couette flow (SCF) and
superfluid Taylor--Couette flow. In axisymmetric superfluid SCF, the number of
meridional circulation cells multiplies as $\Rey$ increases, and their shapes
become more complex, especially in the superfluid component, with multiple
secondary cells arising for $\Rey > 10^3$. The torque exerted by the normal
component is approximately three times greater in a superfluid with anisotropic
Hall--Vinen (HV) mutual friction than in a classical viscous fluid or a
superfluid with isotropic Gorter-Mellink (GM) mutual friction. HV mutual
friction also tends to "pinch" meridional circulation cells more than GM mutual
friction. The boundary condition on the superfluid component, whether no slip
or perfect slip, does not affect the large-scale structure of the flow
appreciably, but it does alter the cores of the circulation cells, especially
at lower $\Rey$. As $\Rey$ increases, and after initial transients die away,
the mutual friction force dominates the vortex tension, and the streamlines of
the superfluid and normal fluid components increasingly resemble each other. In
nonaxisymmetric superfluid SCF, three-dimensional vortex structures are
classified according to topological invariants.
-
Experimental evidence is reviewed for the existence of superfluid turbulence
in a differentially rotating, spherical shell at high Reynolds numbers
($\Rey\gsim 10^3$), such as the outer core of a neutron star. It is shown that
torque variability increases with $\Rey$, suggesting that glitch activity in
radio pulsars may be a function of $\Rey$ as well. The $\Rey$ distribution of
the 67 glitching radio pulsars with characteristic ages $\tau_c \leq 10^6$ {\rm
yr} is constructed from radio timing data and cooling curves and compared with
the $\Rey$ distribution of all 348 known pulsars with $\tau_c \leq 10^6$ {\rm
yr}. The two distributions are different, with a Kolmogorov-Smirnov probability
$\geq 1 - 3.9 \times 10^{-3}$. The conclusion holds for (modified) Urca and
nonstandard cooling, and for Newtonian and superfluid viscosities.
-
We test statistically the hypothesis that radio pulsar glitches result from
an avalanche process, in which angular momentum is transferred erratically from
the flywheel-like superfluid in the star to the slowly decelerating, solid
crust via spatially connected chains of local, impulsive, threshold-activated
events, so that the system fluctuates around a self-organised critical state.
Analysis of the glitch population (currently 285 events from 101 pulsars)
demonstrates that the size distribution in individual pulsars is consistent
with being scale invariant, as expected for an avalanche process. The
waiting-time distribution is consistent with being exponential in seven out of
nine pulsars where it can be measured reliably, after adjusting for
observational limits on the minimum waiting time, as for a constant-rate
Poisson process. PSR J0537$-$6910 and PSR J0835$-$4510 are the exceptions;
their waiting-time distributions show evidence of quasiperiodicity. In each
object, stationarity requires that the rate $\lambda$ equals $- \epsilon
\dot{\nu} / <\Delta\nu>$, where $\dot{\nu}$ is the angular acceleration of the
crust, $<\Delta\nu>$ is the mean glitch size, and $\epsilon\dot{\nu}$ is the
relative angular acceleration of the crust and superfluid. There is no evidence
that $\lambda$ changes monotonically with spin-down age. The rate distribution
itself is fitted reasonably well by an exponential for $\lambda \geq 0.25 {\rm
yr^{-1}}$. For $\lambda < 0.25 {\rm yr^{-1}}$, its exact form is unknown; the
exponential overestimates the number of glitching pulsars observed at low
$\lambda$, where the limited total observation time exercises a selection bias.
-
We investigate the global transition from a turbulent state of superfluid
vorticity to a laminar state, and vice versa, in the outer core of a neutron
star. By solving numerically the hydrodynamic Hall-Vinen-Bekarevich-Khalatnikov
equations for a rotating superfluid in a differentially rotating spherical
shell, we find that the meridional counterflow driven by Ekman pumping exceeds
the Donnelly-Glaberson threshold throughout most of the outer core, exciting
unstable Kelvin waves which disrupt the rectilinear vortex array, creating a
vortex tangle. In the turbulent state, the torque exerted on the crust
oscillates, and the crust-core coupling is weaker than in the laminar state.
This leads to a new scenario for the rotational glitches observed in radio
pulsars: a vortex tangle is sustained in the differentially rotating outer core
by the meridional counterflow, a sudden spin-up event brings the crust and core
into corotation, the vortex tangle relaxes back to a rectilinear vortex array,
then the crust spins down electromagnetically until enough meridional
counterflow builds up to reform a vortex tangle. The turbulent-laminar
transition can occur uniformly or in patches; the associated time-scales are
estimated from vortex filament theory. We calculate numerically the global
structure of the flow with and without an inviscid superfluid component, for
Hall-Vinen and Gorter-Mellink forms of the mutual friction. We also calculate
the post-glitch evolution of the angular velocity of the crust and its time
derivative, and compare the results with radio pulse timing data, predicting a
correlation between glitch activity and Reynolds number.
-
The gravitational wave signal generated by global, nonaxisymmetric shear
flows in a neutron star is calculated numerically by integrating the
incompressible Navier--Stokes equation in a spherical, differentially rotating
shell. At Reynolds numbers $\Rey \gsim 3 \times 10^{3}$, the laminar Stokes
flow is unstable and helical, oscillating Taylor--G\"ortler vortices develop.
The gravitational wave strain generated by the resulting kinetic-energy
fluctuations is computed in both $+$ and $\times$ polarizations as a function
of time. It is found that the signal-to-noise ratio for a coherent,
$10^{8}$-{\rm s} integration with LIGO II scales as $ 6.5 (\Omega_*/10^{4} {\rm
rad} {\rm s}^{-1})^{7/2}$ for a star at 1 {\rm kpc} with angular velocity
$\Omega_*$. This should be regarded as a lower limit: it excludes pressure
fluctuations, herringbone flows, Stuart vortices, and fully developed
turbulence (for $\Rey \gsim 10^{6}$).
-
We integrate for the first time the hydrodynamic
Hall-Vinen-Bekarevich-Khalatnikov equations of motion of a $^{1}S_{0}$-paired
neutron superfluid in a rotating spherical shell, using a pseudospectral
collocation algorithm coupled with a time-split fractional scheme. Numerical
instabilities are smoothed by spectral filtering. Three numerical experiments
are conducted, with the following results. (i) When the inner and outer spheres
are put into steady differential rotation, the viscous torque exerted on the
spheres oscillates quasiperiodically and persistently (after an initial
transient). The fractional oscillation amplitude ($\sim 10^{-2}$) increases
with the angular shear and decreases with the gap width. (ii) When the outer
sphere is accelerated impulsively after an interval of steady differential
rotation, the torque increases suddenly, relaxes exponentially, then oscillates
persistently as in (i). The relaxation time-scale is determined principally by
the angular velocity jump, whereas the oscillation amplitude is determined
principally by the gap width. (iii) When the mutual friction force changes
suddenly from Hall-Vinen to Gorter-Mellink form, as happens when a rectilinear
array of quantized Feynman-Onsager vortices is destabilized by a counterflow to
form a reconnecting vortex tangle, the relaxation time-scale is reduced by a
factor of $\sim 3$ compared to (ii), and the system reaches a stationary state
where the torque oscillates with fractional amplitude $\sim 10^{-3}$ about a
constant mean value. Preliminary scalings are computed for observable
quantities like angular velocity and acceleration as functions of Reynolds
number, angular shear, and gap width. The results are applied to the timing
irregularities (e.g., glitches and timing noise) observed in radio pulsars.
We present the results of a LIGO search for gravitational waves (GWs) associated with GRB 051103, a short-duration hard-spectrum gamma-ray burst (GRB) whose electromagnetically determined sky position is coincident with the spiral galaxy M81, which is 3.6 Mpc from Earth. Possible progenitors for short-hard GRBs include compact object mergers and soft gamma repeater (SGR) giant flares. A merger progenitor would produce a characteristic GW signal that should be detectable at the distance of M81, while GW emission from an SGR is not expected to be detectable at that distance. We found no evidence of a GW signal associated with GRB 051103. Assuming weakly beamed gamma-ray emission with a jet semi-angle of 30 deg we exclude a binary neutron star merger in M81 as the progenitor with a confidence of 98%. Neutron star-black hole mergers are excluded with > 99% confidence. If the event occurred in M81 our findings support the the hypothesis that GRB 051103 was due to an SGR giant flare, making it the most distant extragalactic magnetar observed to date.
We report on a search for gravitational waves from coalescing compact binaries using LIGO and Virgo observations between July 7, 2009 and October 20, 2010. We searched for signals from binaries with total mass between 2 and 25 solar masses; this includes binary neutron stars, binary black holes, and binaries consisting of a black hole and neutron star. The detectors were sensitive to systems up to 40 Mpc distant for binary neutron stars, and further for higher mass systems. No gravitational-wave signals were detected. We report upper limits on the rate of compact binary coalescence as a function of total mass, including the results from previous LIGO and Virgo observations. The cumulative 90%-confidence rate upper limits of the binary coalescence of binary neutron star, neutron star- black hole and binary black hole systems are 1.3 x 10^{-4}, 3.1 x 10^{-5} and 6.4 x 10^{-6} Mpc^{-3}yr^{-1}, respectively. These upper limits are up to a factor 1.4 lower than previously derived limits. We also report on results from a blind injection challenge.
Aims. A transient astrophysical event observed in both gravitational wave (GW) and electromagnetic (EM) channels would yield rich scientific rewards. A first program initiating EM follow-ups to possible transient GW events has been developed and exercised by the LIGO and Virgo community in association with several partners. In this paper, we describe and evaluate the methods used to promptly identify and localize GW event candidates and to request images of targeted sky locations. Methods. During two observing periods (Dec 17 2009 to Jan 8 2010 and Sep 2 to Oct 20 2010), a low-latency analysis pipeline was used to identify GW event candidates and to reconstruct maps of possible sky locations. A catalog of nearby galaxies and Milky Way globular clusters was used to select the most promising sky positions to be imaged, and this directional information was delivered to EM observatories with time lags of about thirty minutes. A Monte Carlo simulation has been used to evaluate the low-latency GW pipeline's ability to reconstruct source positions correctly. Results. For signals near the detection threshold, our low-latency algorithms often localized simulated GW burst signals to tens of square degrees, while neutron star/neutron star inspirals and neutron star/black hole inspirals were localized to a few hundred square degrees. Localization precision improves for moderately stronger signals. The correct sky location of signals well above threshold and originating from nearby galaxies may be observed with ~50% or better probability with a few pointings of wide-field telescopes.
We report on an all-sky search for periodic gravitational waves in the frequency band 50-800 Hz and with the frequency time derivative in the range of 0 through -6e-9 Hz/s. Such a signal could be produced by a nearby spinning and slightly non-axisymmetric isolated neutron star in our galaxy. After recent improvements in the search program that yielded a 10x increase in computational efficiency, we have searched in two years of data collected during LIGO's fifth science run and have obtained the most sensitive all-sky upper limits on gravitational wave strain to date. Near 150 Hz our upper limit on worst-case linearly polarized strain amplitude $h_0$ is 1e-24, while at the high end of our frequency range we achieve a worst-case upper limit of 3.8e-24 for all polarizations and sky locations. These results constitute a factor of two improvement upon previously published data. A new detection pipeline utilizing a Loosely Coherent algorithm was able to follow up weaker outliers, increasing the volume of space where signals can be detected by a factor of 10, but has not revealed any gravitational wave signals. The pipeline has been tested for robustness with respect to deviations from the model of an isolated neutron star, such as caused by a low-mass or long-period binary companion.
The gravitational-wave (GW) sky may include nearby pointlike sources as well as astrophysical and cosmological stochastic backgrounds. Since the relative strength and angular distribution of the many possible sources of GWs are not well constrained, searches for GW signals must be performed in a model-independent way. To that end we perform two directional searches for persistent GWs using data from the LIGO S5 science run: one optimized for pointlike sources and one for arbitrary extended sources. The latter result is the first of its kind. Finding no evidence to support the detection of GWs, we present 90% confidence level (CL) upper-limit maps of GW strain power with typical values between 2-20x10^-50 strain^2 Hz^-1 and 5-35x10^-49 strain^2 Hz^-1 sr^-1 for pointlike and extended sources respectively. The limits on pointlike sources constitute a factor of 30 improvement over the previous best limits. We also set 90% CL limits on the narrow-band root-mean-square GW strain from interesting targets including Sco X-1, SN1987A and the Galactic Center as low as ~7x10^-25 in the most sensitive frequency range near 160 Hz. These limits are the most constraining to date and constitute a factor of 5 improvement over the previous best limits.
We present direct upper limits on continuous gravitational wave emission from the Vela pulsar using data from the Virgo detector's second science run. These upper limits have been obtained using three independent methods that assume the gravitational wave emission follows the radio timing. Two of the methods produce frequentist upper limits for an assumed known orientation of the star's spin axis and value of the wave polarization angle of, respectively, $1.9\ee{-24}$ and $2.2\ee{-24}$, with 95% confidence. The third method, under the same hypothesis, produces a Bayesian upper limit of $2.1\ee{-24}$, with 95% degree of belief. These limits are below the indirect {\it spin-down limit} of $3.3\ee{-24}$ for the Vela pulsar, defined by the energy loss rate inferred from observed decrease in Vela's spin frequency, and correspond to a limit on the star ellipticity of $\sim 10^{-3}$. Slightly less stringent results, but still well below the spin-down limit, are obtained assuming the star's spin axis inclination and the wave polarization angles are unknown.
Soft gamma repeaters (SGRs) and anomalous X-ray pulsars (AXPs) are thought to be magnetars: neutron stars powered by extreme magnetic fields. These rare objects are characterized by repeated and sometimes spectacular gamma-ray bursts. The burst mechanism might involve crustal fractures and excitation of non-radial modes which would emit gravitational waves (GWs). We present the results of a search for GW bursts from six galactic magnetars that is sensitive to neutron star f-modes, thought to be the most efficient GW emitting oscillatory modes in compact stars. One of them, SGR 0501+4516, is likely ~1 kpc from Earth, an order of magnitude closer than magnetars targeted in previous GW searches. A second, AXP 1E 1547.0-5408, gave a burst with an estimated isotropic energy >10^{44} erg which is comparable to the giant flares. We find no evidence of GWs associated with a sample of 1279 electromagnetic triggers from six magnetars occurring between November 2006 and June 2009, in GW data from the LIGO, Virgo, and GEO600 detectors. Our lowest model-dependent GW emission energy upper limits for band- and time-limited white noise bursts in the detector sensitive band, and for f-mode ringdowns (at 1090 Hz), are 3.0x10^{44} d_1^2 erg and 1.4x10^{47} d_1^2 erg respectively, where d_1 = d_{0501} / 1 kpc and d_{0501} is the distance to SGR 0501+4516. These limits on GW emission from f-modes are an order of magnitude lower than any previous, and approach the range of electromagnetic energies seen in SGR giant flares for the first time.
We present the first modeled search for gravitational waves using the complete binary black hole gravitational waveform from inspiral through the merger and ringdown for binaries with negligible component spin. We searched approximately 2 years of LIGO data taken between November 2005 and September 2007 for systems with component masses of 1-99 solar masses and total masses of 25-100 solar masses. We did not detect any plausible gravitational-wave signals but we do place upper limits on the merger rate of binary black holes as a function of the component masses in this range. We constrain the rate of mergers for binary black hole systems with component masses between 19 and 28 solar masses and negligible spin to be no more than 2.0 per Mpc^3 per Myr at 90% confidence.
The physical mechanisms responsible for pulsar timing glitches are thought to excite quasi-normal mode oscillations in their parent neutron star that couple to gravitational wave emission. In August 2006, a timing glitch was observed in the radio emission of PSR B0833-45, the Vela pulsar. At the time of the glitch, the two co-located Hanford gravitational wave detectors of the Laser Interferometer Gravitational-wave observatory (LIGO) were operational and taking data as part of the fifth LIGO science run (S5). We present the first direct search for the gravitational wave emission associated with oscillations of the fundamental quadrupole mode excited by a pulsar timing glitch. No gravitational wave detection candidate was found. We place Bayesian 90% confidence upper limits of 6.3e-21 to 1.4e-20 on the peak intrinsic strain amplitude of gravitational wave ring-down signals, depending on which spherical harmonic mode is excited. The corresponding range of energy upper limits is 5.0e44 to 1.3e45 erg.
We apply the postquasistatic approximation, an iterative method for the evolution of self-gravitating spheres of matter, to study the evolution of dissipative and electrically charged distributions in General Relativity. We evolve nonadiabatic distributions assuming an equation of state that accounts for the anisotropy induced by the electric charge. Dissipation is described by streaming out or diffusion approximations. We match the interior solution, in noncomoving coordinates, with the Vaidya-Reissner-Nordstr\"om exterior solution. Two models are considered: i) a Schwarzschild-like shell in the diffusion limit; ii) a Schwarzschild-like interior in the free streaming limit. These toy models tell us something about the nature of the dissipative and electrically charged collapse. Diffusion stabilizes the gravitational collapse producing a spherical shell whose contraction is halted in a short characteristic hydrodynamic time. The streaming out radiation provides a more efficient mechanism for emission of energy, redistributing the electric charge on the whole sphere, while the distribution collapses indefinitely with a longer hydrodynamic time scale.
We present a search for periodic gravitational waves from the neutron star in the supernova remnant Cassiopeia A. The search coherently analyzes data in a 12-day interval taken from the fifth science run of the Laser Interferometer Gravitational-Wave Observatory. It searches gravitational wave frequencies from 100 to 300 Hz, and covers a wide range of first and second frequency derivatives appropriate for the age of the remnant and for different spin-down mechanisms. No gravitational wave signal was detected. Within the range of search frequencies, we set 95% confidence upper limits of 0.7--1.2e-24 on the intrinsic gravitational wave strain, 0.4--4e-4 on the equatorial ellipticity of the neutron star, and 0.005--0.14 on the amplitude of r-mode oscillations of the neutron star. These direct upper limits beat indirect limits derived from energy conservation and enter the range of theoretical predictions involving crystalline exotic matter or runaway r-modes. This is the first gravitational wave search to present upper limits on r-modes.
We apply the postquasistatic approximation to study the evolution of spherically symmetric fluid distributions undergoing dissipation in the form of radial heat flow. For a model which corresponds to an incompressible fluid departing from the static equilibrium, it is not possible to go far from the initial state after the emission of a small amount of energy. Initially collapsing distributions of matter are not permitted. Emission of energy can be considered as a mechanism to avoid the collapse. If the distribution collapses initially and emits one hundredth of the initial mass only the outermost layers evolve. For a model which corresponds to a highly compressed Fermi gas, only the outermost shell can evolve with a shorter hydrodynamic time scale.
The Laser Interferometer Gravitational Wave Observatory (LIGO) is a network of three detectors built to detect local perturbations in the space-time metric from astrophysical sources. These detectors, two in Hanford, WA and one in Livingston, LA, are power-recycled Fabry-Perot Michelson interferometers. In their fifth science run (S5), between November 2005 and October 2007, these detectors accumulated one year of triple coincident data while operating at their designed sensitivity. In this paper, we describe the calibration of the instruments in the S5 data set, including measurement techniques and uncertainty estimation.
We summarize the sensitivity achieved by the LIGO and Virgo gravitational wave detectors for compact binary coalescence (CBC) searches during LIGO's fifth science run and Virgo's first science run. We present noise spectral density curves for each of the four detectors that operated during these science runs which are representative of the typical performance achieved by the detectors for CBC searches. These spectra are intended for release to the public as a summary of detector performance for CBC searches during these science runs.
We apply the post-quasi--static approximation, an iterative method for the evolution of self-gravitating spheres of matter, to study the evolution of anisotropic non-adiabatic radiating and dissipative distributions in General Relativity. Dissipation is described by viscosity and free-streaming radiation, assuming an equation of state to model anisotropy induced by the shear viscosity. We match the interior solution, in non-comoving coordinates, with the Vaidya exterior solution. Two simple models are presented, based on the Schwarzschild and Tolman VI solutions, in the non--adiabatic and adiabatic limit. In both cases the eventual collapse or expansion of the distribution is mainly controlled by the anisotropy induced by the viscosity.
We present an up-to-date, comprehensive summary of the rates for all types of compact binary coalescence sources detectable by the Initial and Advanced versions of the ground-based gravitational-wave detectors LIGO and Virgo. Astrophysical estimates for compact-binary coalescence rates depend on a number of assumptions and unknown model parameters, and are still uncertain. The most confident among these estimates are the rate predictions for coalescing binary neutron stars which are based on extrapolations from observed binary pulsars in our Galaxy. These yield a likely coalescence rate of 100 per Myr per Milky Way Equivalent Galaxy (MWEG), although the rate could plausibly range from 1 per Myr per MWEG to 1000 per Myr per MWEG. We convert coalescence rates into detection rates based on data from the LIGO S5 and Virgo VSR2 science runs and projected sensitivities for our Advanced detectors. Using the detector sensitivities derived from these data, we find a likely detection rate of 0.02 per year for Initial LIGO-Virgo interferometers, with a plausible range between 0.0002 and 0.2 per year. The likely binary neutron-star detection rate for the Advanced LIGO-Virgo network increases to 40 events per year, with a range between 0.4 and 400 per year.
Progenitor scenarios for short gamma-ray bursts (short GRBs) include coalescenses of two neutron stars or a neutron star and black hole, which would necessarily be accompanied by the emission of strong gravitational waves. We present a search for these known gravitational-wave signatures in temporal and directional coincidence with 22 GRBs that had sufficient gravitational-wave data available in multiple instruments during LIGO's fifth science run, S5, and Virgo's first science run, VSR1. We find no statistically significant gravitational-wave candidates within a [-5, +1) s window around the trigger time of any GRB. Using the Wilcoxon-Mann-Whitney U test, we find no evidence for an excess of weak gravitational-wave signals in our sample of GRBs. We exclude neutron star-black hole progenitors to a median 90% CL exclusion distance of 6.7 Mpc.
(Abridged.) The mean-square current quadrupole moment associated with vorticity fluctuations in high-Reynolds-number turbulence in a differentially rotating neutron star is calculated analytically, as are the amplitude and decoherence time of the resulting, stochastic gravitational wave signal. The calculation resolves the subtle question of whether the signal is dominated by the smallest or largest turbulent eddies: for the Kolmogorov-like power spectrum observed in superfluid spherical Couette simulations, the wave strain is controlled by the largest eddies, and the decoherence time approximately equals the maximum eddy turnover time. For a neutron star with spin frequency $\nu_s$ and Rossby number $Ro$, at a distance $d$ from Earth, the root-mean-square wave strain reaches $h_{RMS} \approx 3\times 10^{-24} Ro^3 (\nu_s / 30 Hz)^3 (d/1 kpc)^{-1}$. A cross-correlation search can detect such a source in principle, because the signal decoheres over the time-scale $\tau_c \approx 10^{-3} Ro^{-1} (\nu_s / 30 Hz)^{-1} s$, which is adequately sampled by existing long-baseline interferometers. Hence hydrodynamic turbulence imposes a fundamental noise floor on gravitational wave observations of neutron stars, although its polluting effect may be muted by partial decoherence in the hectohertz band, where current continuous-wave searches are concentrated, for the highest frequency (and hence most powerful) sources.
Experimental and numerical evidence is reviewed for the existence of a Stewartson layer in spherical Couette flow at small Ekman and Rossby numbers ($\Ek \lsim 10^{-3}$, $\Ro \lsim 10^{-2}$), the relevant hydrodynamic regime in the superfluid outer core of a neutron star. Numerical simulations of a superfluid Stewartson layer are presented for the first time, showing how the layer is disrupted by nonaxisymmetric instabilities. The unstable ranges of $\Ek$ and $\Ro$ are compared with estimates of these quantities in radio pulsars that exhibit glitches. It is found that glitching pulsars lie on the stable side of the instability boundary, allowing differential rotation to build up before a glitch.
We solve numerically for the first time the two-fluid, Hall--Vinen--Bekarevich--Khalatnikov (HVBK) equations for a He-II-like superfluid contained in a differentially rotating, spherical shell, generalizing previous simulations of viscous spherical Couette flow (SCF) and superfluid Taylor--Couette flow. In axisymmetric superfluid SCF, the number of meridional circulation cells multiplies as $\Rey$ increases, and their shapes become more complex, especially in the superfluid component, with multiple secondary cells arising for $\Rey > 10^3$. The torque exerted by the normal component is approximately three times greater in a superfluid with anisotropic Hall--Vinen (HV) mutual friction than in a classical viscous fluid or a superfluid with isotropic Gorter-Mellink (GM) mutual friction. HV mutual friction also tends to "pinch" meridional circulation cells more than GM mutual friction. The boundary condition on the superfluid component, whether no slip or perfect slip, does not affect the large-scale structure of the flow appreciably, but it does alter the cores of the circulation cells, especially at lower $\Rey$. As $\Rey$ increases, and after initial transients die away, the mutual friction force dominates the vortex tension, and the streamlines of the superfluid and normal fluid components increasingly resemble each other. In nonaxisymmetric superfluid SCF, three-dimensional vortex structures are classified according to topological invariants.
Experimental evidence is reviewed for the existence of superfluid turbulence in a differentially rotating, spherical shell at high Reynolds numbers ($\Rey\gsim 10^3$), such as the outer core of a neutron star. It is shown that torque variability increases with $\Rey$, suggesting that glitch activity in radio pulsars may be a function of $\Rey$ as well. The $\Rey$ distribution of the 67 glitching radio pulsars with characteristic ages $\tau_c \leq 10^6$ {\rm yr} is constructed from radio timing data and cooling curves and compared with the $\Rey$ distribution of all 348 known pulsars with $\tau_c \leq 10^6$ {\rm yr}. The two distributions are different, with a Kolmogorov-Smirnov probability $\geq 1 - 3.9 \times 10^{-3}$. The conclusion holds for (modified) Urca and nonstandard cooling, and for Newtonian and superfluid viscosities.
We test statistically the hypothesis that radio pulsar glitches result from an avalanche process, in which angular momentum is transferred erratically from the flywheel-like superfluid in the star to the slowly decelerating, solid crust via spatially connected chains of local, impulsive, threshold-activated events, so that the system fluctuates around a self-organised critical state. Analysis of the glitch population (currently 285 events from 101 pulsars) demonstrates that the size distribution in individual pulsars is consistent with being scale invariant, as expected for an avalanche process. The waiting-time distribution is consistent with being exponential in seven out of nine pulsars where it can be measured reliably, after adjusting for observational limits on the minimum waiting time, as for a constant-rate Poisson process. PSR J0537$-$6910 and PSR J0835$-$4510 are the exceptions; their waiting-time distributions show evidence of quasiperiodicity. In each object, stationarity requires that the rate $\lambda$ equals $- \epsilon \dot{\nu} / <\Delta\nu>$, where $\dot{\nu}$ is the angular acceleration of the crust, $<\Delta\nu>$ is the mean glitch size, and $\epsilon\dot{\nu}$ is the relative angular acceleration of the crust and superfluid. There is no evidence that $\lambda$ changes monotonically with spin-down age. The rate distribution itself is fitted reasonably well by an exponential for $\lambda \geq 0.25 {\rm yr^{-1}}$. For $\lambda < 0.25 {\rm yr^{-1}}$, its exact form is unknown; the exponential overestimates the number of glitching pulsars observed at low $\lambda$, where the limited total observation time exercises a selection bias.