Rotational glitches in some rotation-powered pulsars display power-law size and exponential waiting time distributions. These statistics are consistent with a state-dependent Poisson process, where the glitch rate is an increasing function of a global stress variable (e.g. crust-superfluid angular velocity lag), diverges at a threshold stress, increases smoothly while the star spins down, and decreases step-wise at each glitch. A minimal, seven-parameter, maximum likelihood model is calculated for PSR J1740-3015, PSR J0534+2200, and PSR J0631+1036, the three objects with the largest samples whose glitch activity is Poisson-like. The estimated parameters have theoretically reasonable values and contain useful information about the glitch microphysics. It is shown that the maximum likelihood, state-dependent Poisson model is a marginally (23-27 per cent) better post factum "predictor" of historical glitch epochs than a homogeneous Poisson process for PSR J1740-3015 and PSR J0631+1036 and a comparable predictor for PSR J0534+2200. Monte Carlo simulations imply that > 50 glitches are needed to test reliably whether one model outperforms the other. It is predicted that the next glitch will occur at Modified Julian Date (MJD) 57784 +/- 256.8, 60713 +/- 1935, and 57406 +/- 1444 for the above three objects respectively. The analysis does not apply to quasiperiodic glitchers like PSR J0537-6910 and PSR J0835-4510, which are not described accurately by the state-dependent Poisson model in its original form.

The equilibrium configurations of neutron superfluid vortices interacting with proton superconductor flux tubes in a rotating, harmonic trap are non-trivial in general, when the magnetorotational symmetry is broken. A non-zero angle $\theta$ between the magnetic and rotation axes leads to tangled vorticity due to competition between vortex-vortex repulsion and vortex-flux-tube pinning. Here we investigate the far-from-equilibrium behaviour of the vortices, as the trap decelerates, by solving the time-dependent, stochastic, Gross-Pitaevskii equation numerically in three dimensions. The numerical simulations reveal new vortex behaviours. Key geometrical attributes of the evolving vortex tangle are characterised, as is the degree to which pinning impedes the deceleration of the neutron condensate as a function of $\eta$, the pinning strength, and $\theta$. The simulated system is a partial analogue of the outer core of a decelerating neutron star, albeit in a very different parameter regime.

Three-dimensional, Gross-Pitaevskii equation (GPE) simulations are presented of the interaction between neutron superfluid vortices and proton superconductor flux tubes in a rotating, harmonic trap, representing an idealised model of the outer core of a neutron star. Low-energy states of the neutron condensate are calculated by evolving the GPE in imaginary time in the presence of a prescribed, static, rectilinear flux tube array. The calculations are carried out as a function of the angle between the global magnetic and rotation axes, and the amplitude and sign of the current-current and density couplings between the neutron and proton condensates. It is found that the system is frustrated by the competition between vortex-vortex repulsion and vortex-flux-tube attraction (pinning), leading to the formation of vortex tangles and "glassy" behaviour characterized by multiple metastable states spaced closely in energy. The dimensionless parameters in the simulations are ordered as one expects in a neutron star, but the dynamic range is many orders of magnitude smaller than in reality, so caution must be exercised when assessing the astrophysical implications. Nevertheless the results suggest that tangled vorticity may be endemic in neutron star outer cores.