The electrical resistivity of the accreted mountain in a millisecond pulsar is limited by the observed spin-down rate of binary radio millisecond pulsars (BRMSPs) and the spins and X-ray fluxes of accreting millisecond pulsars (AMSPs). We find $\eta \ge 10^{-28}\,\mathrm{s}\, (\tau_\mathrm{SD}/1\,\mathrm{Gyr})^{-0.8}$ (where $\tau_\mathrm{SD}$ is the spin-down age) for BRMSPs and $\eta \ge 10^{-25}\,\mathrm{s}\,(\dot{M}_\mathrm{a}/\dot{M}_\mathrm{E})^{0.6}$ (where $\dot{M}_\mathrm{a}$ and $\dot{M}_\mathrm{E}$ are the actual and Eddington accretion rates) for AMSPs. These limits are inferred assuming that the mountain attains a steady state, where matter diffuses resistively across magnetic flux surfaces but is replenished at an equal rate by infalling material. The mountain then relaxes further resistively after accretion ceases. The BRMSP spin-down limit approaches the theoretical electron-impurity resistivity at temperatures $\ga 10^5$ K for an impurity concentration of $\sim 0.1$, while the AMSP stalling limit falls two orders of magnitude below the theoretical electron-phonon resistivity for temperatures above $10^8$ K. Hence BRMSP observations are already challenging theoretical resistivity calculations in a useful way. Next-generation gravitational-wave interferometers will constrain $\eta$ at a level that will be competitive with electromagnetic observations.

We perform ideal-magnetohydrodynamic axisymmetric simulations of magnetically confined mountains on an accreting neutron star, with masses less than ~0.12 solar masses. We consider two scenarios, in which the mountain sits atop a hard surface or sinks into a soft, fluid base. We find that the ellipticity of the star, due to a mountain grown on a hard surface, approaches ~2e-4 for accreted masses greater than ~1.2e-3 solar masses, and that sinking reduces the ellipticity by between 25% and 60%. The consequences for gravitational radiation from low-mass x-ray binaries are discussed.

Three-dimensional numerical magnetohydrodynamic (MHD) simulations are performed to investigate how a magnetically confined mountain on an accreting neutron star relaxes resistively. No evidence is found for non-ideal MHD instabilities on a short time-scale, such as the resistive ballooning mode or the tearing mode. Instead, the mountain relaxes gradually as matter is transported across magnetic surfaces on the diffusion time-scale, which evaluates to $\tau_\mathrm{I} \sim 10^5 - 10^8$ yr (depending on the conductivity of the neutron star crust) for an accreted mass of $M_a = 1.2 \times 10^{-4} M_\odot$. The magnetic dipole moment simultaneously reemerges as the screening currents dissipate over $\tau_\mathrm{I}$. For nonaxisymmetric mountains, ohmic dissipation tends to restore axisymmetry by magnetic reconnection at a filamentary neutral sheet in the equatorial plane. Ideal-MHD oscillations on the Alfv\'{e}n time-scale, which can be excited by external influences, such as variations in the accretion torque, compress the magnetic field and hence decrease $\tau_\mathrm{I}$ by \change{one order of magnitude} relative to its standard value (as computed for the static configuration). The implications of long-lived mountains for gravitational wave emission from low-mass X-ray binaries are briefly explored.

We give an improved estimate of the detectability of gravitational waves from magnetically confined mountains on accreting neutron stars. The improved estimate includes the following effects for the first time: three-dimensional hydromagnetic ("fast") relaxation, three-dimensional resistive ("slow") relaxation, realistic accreted masses $M_a \la 2 \times 10^{-3} M_\odot$, (where the mountain is grown ab initio by injection), and verification of the curvature rescaling transformation employed in previous work. Typically, a mountain does not relax appreciably over the lifetime of a low-mass X-ray binary. The ellipticity reaches $\epsilon \approx 2 \times 10^{-5}$ for $M_a=2\times 10^{-3} M_\odot$. The gravitational wave spectrum for triaxial equilibria contains an additional line, which, although weak, provides valuable information about the mountain shape. We evaluate the detectability of magnetic mountains with Initial and Advanced LIGO. For a standard, coherent matched filter search, we find a signal-to-noise ratio of $d = 28 (M_a/10^{-4} M_\odot) (1+5.5 M_a/10^{-4} M_\odot)^{-1} (D/10 \mathrm{kpc})^{-1} (T_0/14 \mathrm{d})^{1/2}$ for Initial LIGO, where $D$ is the distance and $T_0$ is the observation time. From the nondetection of gravitational waves from low-mass X-ray binaries to date, and the wave strain limits implied by the spin frequency distribution of these objects (due to gravitational wave braking), we conclude that there are other, as yet unmodelled, physical effects that further reduce he quadrupole moment of a magnetic mountain, most notably sinking into the crust.

We compute the continuous part of the ideal-magnetohydrodynamic (ideal-MHD) frequency spectrum of a polar mountain produced by magnetic burial on an accreting neutron star. Applying the formalism developed by Hellsten & Spies (1979), extended to include gravity, we solve the singular eigenvalue problem subject to line-tying boundary conditions. This spectrum divides into an Alfv\'{e}n part and a cusp part. The eigenfunctions are chirped and anharmonic with an exponential envelope, and the eigenfrequencies cover the whole spectrum above a minimum $\omega_\mathrm{low}$. For equilibria with accreted mass $1.2 \times 10^{-6} \la M_a/M_\odot \la 1.7 \times 10^{-4}$ and surface magnetic fields $10^{11} \la B_\ast/\mathrm{G} \la 10^{13}$, $\omega_\mathrm{low}$ is approximately independent of $B_\ast$, and increases with $M_a$. The results are consistent with the Alfv\'{e}n spectrum excited in numerical simulations with the \textsc{zeus-mp} solver. The spectrum is modified substantially by the Coriolis force in neutron stars spinning faster than $\sim 100$ Hz. The implications for gravitational wave searches for low-mass X-ray binaries are considered briefly.

It is widely assumed that the observed reduction of the magnetic field of millisecond pulsars can be connected to the accretion phase during which the pulsar is spun up by mass accretion from a companion. A wide variety of reduction mechanisms have been proposed, including the burial of the field by a magnetic mountain, formed when the accreted matter is confined to the poles by the tension of the stellar magnetic field. A magnetic mountain effectively screens the magnetic dipole moment. On the other hand, observational data suggests that accreting neutron stars are sources of gravitational waves, and magnetic mountains are a natural source of a time-dependent quadrupole moment. We show that the emission is sufficiently strong to be detectable by current and next generation long-baseline interferometers. Preliminary results from fully three-dimensional magnetohydrodynamic (MHD) simulations are presented. We find that the initial axisymmetric state relaxes into a nearly axisymmetric configuration via toroidal magnetic modes. A substantial quadrupole moment is still present in the final state, which is stable (in ideal MHD) yet highly distorted.

We examine the hydromagnetic stability of magnetically confined mountains, which arise when material accumulates at the magnetic poles of an accreting neutron star. We extend a previous axisymmetric stability analysis by performing three-dimensional simulations using the ideal-magnetohydrodynamic (ideal-MHD) code \textsc{zeus-mp}, investigating the role played by boundary conditions, accreted mass, stellar curvature, and (briefly) toroidal magnetic field strength. We find that axisymmetric equilibria are susceptible to the undular sub-mode of the Parker instability but are not disrupted. The line-tying boundary condition at the stellar surface is crucial in stabilizing the mountain. The nonlinear three-dimensional saturation state of the instability is characterized by a small degree of nonaxisymmetry ($\la 0.1$ per cent) and a mass ellipticity of $\epsilon \sim 10^{-5}$ for an accreted mass of $M_a = 10^{-5} M_\odot$. Hence there is a good prospect of detecting gravitational waves from accreting millisecond pulsars with long-baseline interferometers such as Advanced LIGO. We also investigate the ideal-MHD spectrum of the system, finding that long-wavelength poloidal modes are suppressed in favour of toroidal modes in the nonaxisymmetric saturation state.

We present three-dimensional, nonrelativistic, hydrodynamic simulations of bow shocks in pulsar wind nebulae. The simulations are performed for a range of initial and boundary conditions to quantify the degree of asymmetry produced by latitudinal variations in the momentum flux of the pulsar wind, radiative cooling in the postshock flow, and density gradients in the interstellar medium (ISM). We find that the bow shock is stable even when travelling through a strong ISM gradient. We demonstrate how the shape of the bow shock changes when the pulsar encounters density variations in the ISM. We show that a density wall can account for the peculiar bow shock shapes of the nebulae around PSR J2124-3358 and PSR B0740-28. A wall produces kinks in the shock, whereas a smooth ISM density gradient tilts the shock. We conclude that the anisotropy of the wind momentum flux alone cannot explain the observed bow shock morphologies but it is instead necessary to take into account external effects. We show that the analytic (single layer, thin shell) solution is a good approximation when the momentum flux is anisotropic, fails for a steep ISM density gradient, and ap- proaches the numerical solution for efficient cooling. We provide analytic expressions for the latitudinal dependence of a vacuum-dipole wind and the associated shock shape, and compare the results to a split-monopole wind. We find that we are unable to distinguish between these two wind models purely from the bow shock morphology.