This paper starts from a toy model for inflation in a class of modified theories of gravity in the metric formalism. Instead of the standard procedure -- assuming a non-linear Lagrangian $f(R)$ in the Jordan frame -- we start from a simple $\phi^2$ potential in the Einstein frame and investigate the corresponding $f(R)$ in the former picture. The addition of an ad-hoc Cosmological Constant in the Einstein frame leads to a Thermodynamical interpretation of this physical system, which allows further insight on its (meta)stability and evolution.

We will make a comparison between the dynamics of spherical gravitational collapse for a perturbed FLRW universe to first order in the context of general relativity, with the corresponding results obtained for the gravitational collapse under theories of modified gravity f(R). This work is aimed at obtaining an analytical model of explanation a source of large scale structures in the universe presenting an approximation of model spherical gravitational collapse under theories of modified gravity f(R).

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.

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.

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.

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 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 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.