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 e
xplored 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 c
ondensates. 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} [(Deltaomega/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, $Deltaomega/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 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 domai
ns 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 $-2leq 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 powe
r 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.
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.3leq a leq -2.0$ and $-5.5leq bleq -2.2$, and mean glitching rates in the range $0.0023leqlambdaleq0.13$ in units of inverse time.
