A framework is presented for a statistical theory of neutron star glitches, motivated by the results emerging from recent Gross-Pitaevskii simulations of pinned, decelerating quantum condensates. It is shown that the observed glitch size distributions cannot be reproduced if superfluid vortices unpin independently via a Poisson process; the central limit theorem yields a narrow Gaussian for the size distribution, instead of the broad, power-law tail observed. This conclusion is not altered fundamentally when a range of pinning potentials is included, which leads to excavation of the potential distribution of occupied sites, vortex accumulation at strong pinning sites, and hence the occasional, abnormally large glitch. Knock-on processes are therefore needed to make the unpinning rate of a vortex conditional on the pinning state of its near and/or remote neighbours, so that the Gaussian size distributions resulting generically from the central limit theorem are avoided. At least two knock-on processes, nearest- neighbour proximity knock-on and remote acoustic knock-on, are clearly evident in the Gross-Pitaevskii simulation output. It is shown that scale-invariant (i.e. power-law) vortex avalanches occur when knock-on is included, provided that two specific relations hold between the temperature and spin-down torque. This fine tuning is unlikely in an astronomical setting, leaving the overall problem partly unsolved. A state-dependent Poisson formalism is presented which will form the basis of future studies in this area.

The pinning and collective unpinning of superfluid vortices in a decelerating container is a key element of the canonical model of neutron star glitches and laboratory spin-down experiments with helium II. Here the dynamics of vortex (un)pinning is explored using numerical Gross-Pitaevskii calculations, with a view to understanding the triggers for catastrophic unpinning events (vortex avalanches) that lead to rotational glitches. We explicitly identify three triggers: rotational shear between the bulk condensate and the pinned vortices, a vortex proximity effect driven by the repulsive vortex-vortex interaction, and sound waves emitted by moving and repinning vortices. So long as dissipation is low, sound waves emitted by a repinning vortex are found to be sufficiently strong to unpin a nearby vortex. For both ballistic and forced vortex motion, the maximum inter-vortex separation required to unpin scales inversely with pinning strength.

The current-quadrupole gravitational-wave signal emitted during the spin-up phase of a pulsar glitch is calculated from first principles by modeling the vortex dynamics observed in recent Gross-Pitaevskii simulations of pinned, decelerating quantum condensates. Homogeneous and inhomogeneous unpinning geometries, representing creep- and avalanche-like glitches, provide lower and upper bounds on the gravitational wave signal strength respectively. The signal arising from homogeneous glitches is found to scale with the square root of glitch size, whereas the signal from inhomogeneous glitches scales proportional to glitch size. The signal is also computed as a function of vortex travel distance and stellar angular velocity. Convenient amplitude scalings are derived as functions of these parameters. For the typical astrophysical situation, where the glitch duration (in units of the spin period) is large compared to the vortex travel distance (in units of the stellar radius), an individual glitch from an object $1\,\rm{kpc}$ from Earth generates a wave strain of $10^{-24} [(\Delta\omega/\omega) / 10^{-7}] (\omega/10^2 \rm{rad s}^{-1})^3 (\Delta r / 10^{-2} \rm{m})^{-1}$, where $\Delta r$ is the average distance travelled by a vortex during a glitch, $\Delta\omega/\omega$ is the fractional glitch size, and $\omega$ is the pulsar angular velocity. The non-detection of a signal from the 2006 Vela glitch in data from the fifth science run conducted by the Laser Interferometer Gravitational-Wave Observatory implies that the glitch duration exceeds $\sim 10^{-4}\,\rm{ms}$. This represents the first observational lower bound on glitch duration to be obtained.

The first large-scale quantum mechanical simulations of pulsar glitches are presented, using a Gross-Pitaevskii model of the crust-superfluid system in the presence of pinning. Power-law distributions of simulated glitch sizes are obtained, in accord with astronomical observations, with exponents ranging from -0.55 to -1.26. Examples of large-scale simulations, containing $\sim 200$ vortices, reveal that these statistics persist in the many-vortex limit. Waiting-time distributions are also constructed. These and other statistics support the hypothesis that catastrophic unpinning mediated by collective vortex motion produces glitches; indeed, such collective events are seen in time-lapse movies of superfluid density. Three principal trends are observed. (1) The glitch rate scales proportional to the electromagnetic spin-down torque. (2) A strong positive correlation is found between the strength of vortex pinning and mean glitch size. (3) The spin-down dynamics depend less on the pinning site abundance once the latter exceeds one site per vortex, suggesting that unpinned vortices travel a distance comparable to the inter-vortex spacing before repinning.

The scale-invariant glitch statistics observed in individual pulsars (exponential waiting-time and power-law size distributions) are consistent with a critical self-organization process, wherein superfluid vortices pin metastably in macroscopic domains and unpin collectively via nearest-neighbor avalanches. Macroscopic inhomogeneity emerges naturally if pinning occurs at crustal faults. If, instead, pinning occurs at lattice sites and defects, which are macroscopically homogeneous, we show that an alternative, noncritical self-organization process operates, termed coherent noise, wherein the global Magnus force acts uniformly on vortices trapped in a range of pinning potentials and undergoing thermal creep. It is found that vortices again unpin collectively, but not via nearest-neighbor avalanches, and that, counterintuitively, the resulting glitch sizes are scale invariant, in accord with observational data. A mean-field analytic theory of the coherent noise process, supported by Monte-Carlo simulations, yields a power-law size distribution, between the smallest and largest glitch, with exponent $a$ in the range $-2\leq a \leq 0$. When the theory is fitted to data from the nine most active pulsars, including the two quasiperiodic glitchers PSR J0537$-$6910 and PSR J0835$-$4510, it directly constrains the distribution of pinning potentials in the star, leading to two conclusions: (i) the potentials are broadly distributed, with the mean comparable to the standard deviation; and (ii) the mean potential decreases with characteristic age. An observational test is proposed to discriminate between nearest-neighbor avalanches and coherent noise.

We incorporate a contribution to reionization from X-rays within analytic and semi-numerical simulations of the 21-cm signal arising from neutral hydrogen during the epoch of reionization. We explore the impact that X-ray ionizations have on the power spectrum (PS) of 21-cm fluctuations by varying both the average X-ray MFP and the fractional contribution of X-rays to reionization. In general, prior to the epoch when the intergalactic medium is dominated by ionized regions (H {\sevensize II} regions), X-ray-induced ionization enhances fluctuations on spatial scales smaller than the X-ray MFP, provided that X-ray heating does not strongly supress galaxy formation. Conversely, at later times when \H2 regions dominate, small-scale fluctuations in the 21-cm signal are suppressed by X-ray ionization. Our modelling also shows that the modification of the 21-cm signal due to the presence of X-rays is sensitive to the relative scales of the X-ray MFP, and the characteristic size of \H2 regions. We therefore find that X-rays imprint an epoch and scale-dependent signature on the 21-cm PS, whose prominence depends on fractional X-ray contribution. The degree of X-ray heating of the IGM also determines the extent to which these features can be discerned. We show that the MWA will have sufficient sensitivity to detect this modification of the PS, so long as the X-ray photon MFP falls within the range of scales over which the array is most sensitive ($\sim0.1$ Mpc$^{-1}$). In cases in which this MFP takes a much smaller value, an array with larger collecting area would be required.

We investigate the impact of neutral hydrogen (HI) in galaxies on the statistics of 21-cm fluctuations using analytic and semi-numerical modelling. Following the reionisation of hydrogen the HI content of the Universe is dominated by damped absorption systems (DLAs), with a cosmic density in HI that is observed to be constant at a level equal to ~2% of the cosmic baryon density from z~1 to z~5. We show that extrapolation of this constant fraction into the reionisation epoch results in a reduction of 10-20% in the amplitude of 21-cm fluctuations over a range of spatial scales. The assumption of a different percentage during the reionisation era results in a proportional change in the 21-cm fluctuation amplitude. We find that consideration of HI in galaxies/DLAs reduces the prominence of the HII region induced shoulder in the 21-cm power spectrum (PS), and hence modifies the scale dependence of 21-cm fluctuations. We also estimate the 21cm-galaxy cross PS, and show that the cross PS changes sign on scales corresponding to the HII regions. From consideration of the sensitivity for forthcoming low-frequency arrays we find that the effects of HI in galaxies/DLAs on the statistics of 21-cm fluctuations will be significant with respect to the precision of a PS or cross PS measurement. In addition, since overdense regions are reionised first we demonstrate that the cross-correlation between galaxies and 21-cm emission changes sign at the end of the reionisation era, providing an alternative avenue to pinpoint the end of reionisation. The sum of our analysis indicates that the HI content of the galaxies that reionise the universe will need to be considered in detailed modelling of the 21-cm intensity PS in order to correctly interpret measurements from forthcoming low-frequency arrays.

A cellular automaton model of pulsar glitches is described, based on the superfluid vortex unpinning paradigm. Recent analyses of pulsar glitch data suggest that glitches result from scale-invariant avalanches \citep{Melatos07a}, which are consistent with a self-organized critical system (SOCS). A cellular automaton provides a computationally efficient means of modelling the collective behaviour of up to $10^{16}$ vortices in the pulsar interior, whilst ensuring that the dominant aspects of the microphysics are not lost. The automaton generates avalanche distributions that are qualitatively consistent with a SOCS and with glitch data. The probability density functions of glitch sizes and durations are power laws, and the probability density function of waiting times between successive glitches is Poissonian, consistent with statistically independent events. The output of the model depends on the physical and computational paramters used. The fitted power law exponents for the glitch sizes ($a$) and durations ($b$) decreases as the strength of the vortex pinning increases. Similarly the exponents increase as the fraction of vortices that are pinned decreases. For the physical and computational parameters considered, one finds $-4.3\leq a \leq -2.0$ and $-5.5\leq b\leq -2.2$, and mean glitching rates in the range $0.0023\leq\lambda\leq0.13\$ in units of inverse time.