We propose and discuss an alternative pulsar radio emission mechanism that relies on rotation-driven plasma oscillations, rather than on a beam-driven instability, and suggest that it may be the generic radio emission mechanism for pulsars. We identify these oscillations as superluminal longitudinal waves in the pulsar plasma, and point out that these waves can escape directly in the O~mode. We argue that the frequency of the oscillations is $\omega_0\approx\omega_{\rm p}(2\langle\gamma\rangle)^{1/2}\gamma_{\rm s}$, where $\gamma_{\rm s}$ is the Lorentz factor of bulk streaming motion and $\langle\gamma\rangle$ is the mean Lorentz factor in the rest frame of the plasma. The dependence of the plasma frequency $\omega_{\rm p}$ on radial distance implies a specific frequency-to-radius mapping, $\omega_0\propto r^{-3/2}$. Escape of the energy in these oscillations is possible only if they are generated in overdense, field-aligned regions that we call fibers; the wave energy is initially refracted into underdense regions between the fibers, which act as ducts. Some implications of the model for the interpretation of pulsar radio emission are discussed.

We consider critically the three most widely favored pulsar radio emission mechanisms: coherent curvature emission (CCE), beam-driven relativistic plasma emission (RPE) and anomalous Doppler emission (ADE). We assume that the pulsar plasma is one dimensional (1D), streaming outward with a bulk Lorentz factor $\gamma_{\rm s} \gg \langle\gamma\rangle -1 \gtrsim 1$, where $\langle\gamma\rangle$ is the intrinsic spread in the rest frame of the plasma. We argue that the formation of beams in a multi-cloud model is ineffective in the intrinsically relativistic case for plausible parameters, because the overtaking takes too long. We argue that the default choice for the particle distribution in the rest frame is a J{\"u}ttner distribution and that relativistic streaming should be included by applying a Lorentz transformation to the rest-frame distribution, rather than the widely assumed relativistically streaming Gaussian distribution. We find that beam-driven wave growth is severely restricted by (a) the wave properties in pulsar plasma, (b) a separation condition between beam and background, and (c) the inhomogeneity of the plasma in the pulsar frame. The growth rate for the kinetic instability is much smaller and the bandwidth of the growing waves is much larger for a J{\"u}ttner distribution than for a relativistically streaming Gaussian distribution. No reactive instability occurs at all for a J{\"u}ttner distribution. We conclude that none of CCE, RPE and ADE in tenable as the generic pulsar radio emission mechanism for ``plausible'' assumptions about the pulsar plasma.

Vector magnetogram data are often used as photospheric boundary conditions for force-free coronal magnetic field extrapolations. In general, however, vector magnetogram data are not consistent with the force-free assumption. In this article, we demonstrate a way to deal with inconsistent boundary data, by generalizing the "self-consistency procedure" of Wheatland & Regnier (2009). In that procedure, the inconsistency is resolved by an iterative process of constructing two solutions based on the values of the force-free parameter alpha on the two polarities of the field in the boundary (the P and N polarities), and taking uncertainty-weighted averages of the boundary alpha values in the P and N solutions. When the alpha values in the P and N regions are very different, the self-consistent solution may lose high alpha values from the boundary conditions. We show how, by altering the weighting of the uncertainties in the P or N boundary conditions, we can preserve high alpha values in the self-consistent solution. The weighted self-consistent extrapolation method is demonstrated on an analytic bipole field and applied to vector magnetogram data taken by the Helioseismic and Magnetic Imager (HMI) instrument for NOAA active region AR 12017 on 2014 March 29.

Weak heating events are frequent and ubiquitous in solar corona. They derive their energy from the local magnetic field and form a major source of local heating, signatures of which are seen in EUV and X-ray bands. Their radio emission arise from various plasma instabilities that lead to coherent radiation, making even a weak flare appear very bright. The radio observations hence probe non-equilibrium dynamics providing complementary information on plasma evolution. However, a robust study of radio emission from a weak event among many simultaneous events, requires high dynamic range imaging at sub-second andsub-MHz resolutions owing to their high spectro-temporal variability. Such observations were not possible until recently.This is among the first spatially resolved studies of an active region loop hosting a transient brightening (ARTB) and dynamically linked to a metrewave type-I noise storm. It uses imaging observations at metrewave, EUV and X-ray bands, along with magnetogram data. We believe this is the first spectroscopic radio imaging study of a type-I source, the data for which was obtained using the Murchison Widefield Array. We report the discovery of 30 s quasi-periodic oscillations (QPOs) in the radio light curve, riding on a coherent baseline flux. The strength of the QPOs and the baseline flux enhanced during a mircoflare associated with the ARTB. Our observations suggest a scenario where magnetic stress builds up over an Alfv\'{e}n timescale (30s) across the typical magnetic field braiding scale and then dissipates via a cascade of weak reconnection events.

Beam-driven instabilities are considered in a pulsar plasma assuming that both the background plasma and the beam are relativistic J\"uttner distributions. In the rest frame of the background, the only waves that can satisfy the resonance condition are in a tiny range of slightly subluminal phase speeds. The growth rate for the kinetic (or maser) version of the weak-beam instability is much smaller than has been estimated for a relativistically streaming Gaussian distribution, and the reasons for this are discussed. The growth rate for the reactive version of the weak-beam instability is treated in a conventional way. We compare the results with exact calculations, and find that the approximate solutions are not consistent with the exact results. We conclude that, for plausible parameters, there is no reactive version of the instability. The growth rate in the pulsar frame is smaller than that in the rest frame of the background plasma by a factor $2\gamma_{\rm s}$, where $\gamma_{\rm s} = 10^2 {\rm -} 10^3$ is the Lorentz factor of the bulk motion of the background plasma, placing a further constraint on effective wave growth. Based on these results, we argue that beam-driven wave growth probably plays no role in pulsar radio emission.

Wave dispersion in a pulsar plasma is discussed emphasizing the relevance of different inertial frames, notably the plasma rest frame ${\cal K}$ and the pulsar frame ${\cal K}'$ in which the plasma is streaming with speed $\beta_{\rm s}$. The effect of a Lorentz transformation on both subluminal, $|z|<1$, and superluminal, $|z|>1$, waves is discussed. It is argued that the preferred choice for a relativistically streaming distribution should be a Lorentz-transformed J\"uttner distribution; such a distribution is compared with other choices including a relativistically streaming Gaussian distribution. A Lorentz transformation of the dielectric tensor is written down, and used to derive an explicit relation between the relativistic plasma dispersion functions in ${\cal K}$ and ${\cal K}'$. It is shown that the dispersion equation can be written in an invariant form, implying a one-to-one correspondence between wave modes in any two inertial frames. Although there are only three modes in the plasma rest frame, it is possible for backward-propagating or negative-frequency solutions in ${\cal K}$ to transform into additional forward-propagating, positive-frequency solutions in ${\cal K}'$ that may be regarded as additional modes.

Wave dispersion in a pulsar plasma (a 1D, strongly magnetized, pair plasma streaming highly relativistically with a large spread in Lorentz factors in its rest frame) is discussed, motivated by interest in beam-driven wave turbulence and the pulsar radio emission mechanism. In the rest frame of the pulsar plasma there are three wave modes in the low-frequency, non-gyrotropic approximation. For parallel propagation these are referred to as the X, A and L modes, with the X and A modes having dispersion relation $z=z_A\approx1-1/2\beta_A^2$, where $z=\omega/k_\parallel c$ is the phase speed and $\beta_Ac$ is the Alfven speed. The L mode dispersion relation is determined by a relativistic plasma dispersion function, $z^2W(z)$, which is negative for $ z < z_0 $ and has a sharp maximum at $z=z_m$, with $1-z_m<1-z_0\ll1$. We give numerical estimates for the maximum of $z^2W(z)$ and for $z_m$ and $z_0$ for a 1D Juttner distribution. The L and A modes reconnect, for $z_A>z_0$, to form the O and Alfven modes for oblique propagation ($\theta\neq0$). For $z_A<z_0$ the Alfven and O~mode curves reconnect forming a new mode that exists only for $\tan^2\theta>z_0^2-z_A^2$. The L mode is the nearest counterpart to Langmuir waves in a nonrelativistic plasma, but we argue that there are no `Langmuir-like' waves in pulsar plasma, identifying three features of the L~mode (dispersion relation, ratio of electric to total energy and group speed) that are not Langmuir-like. A beam-driven instability requires a beam speed equal to the phase speed of the wave. This resonance condition can be satisfied for the O mode, but only for an implausibly energetic beam and only for a tiny range of angles for the O~mode around $\theta\approx0$. The resonance is also possible for the Alfven mode but only near a turnover frequency that has no counterpart for Alfven waves in a nonrelativistic plasma.

Rapid and stepwise changes of the magnetic field are often observed during flares but cannot be explained by models yet. Using a 45 min sequence of SDO/HMI 135 s fast-cadence vector magnetograms of the X1 flare on 2014-03-29 we construct, at each timestep, nonlinear force-free models for the coronal magnetic field. Observed flare-related changes in the line-of-sight magnetic field $B_{\rm LOS}$ at the photosphere and chromosphere are compared with changes in the magnetic fields in the models. We find a moderate agreement at the photospheric layer (the basis for the models), but no agreement at chromospheric layers. The observed changes at the photosphere and chromosphere are surprisingly different, and are unlikely to be reproduced by a force-free model. The observed changes are likely to require a change in the magnitude of the field, not just in its direction.

Flares produce sudden and permanent changes in the horizontal photospheric magnetic field. In particular flares generally produce increased magnetic shear in the photospheric field along the neutral line. Recent observations show also that flares can produce sudden photospheric motion. We present a model for the observed changes as the response of the photosphere to a large-amplitude shear Alfv\'{e}n wave propagating down from the corona on either side of the neutral line. The Alfv\'{e}nic front is assumed to impact the photosphere close to the neutral line first, and then successively further away with time, such that the line of impact coincides with the flare ribbon. The wave introduces magnetic shear and velocity shear. The magnetic shear introduced at the photosphere has the same sign on either side of the neutral line, while the velocity shear has the opposite sign. We discuss the possibility that this process is responsible for particle acceleration in flares.

A precondition for the radio emission of pulsars is the existence of strong, small-scale magnetic field structures (`magnetic spots') in the polar cap region. Their creation can proceed via crustal Hall drift out of two qualitatively and quantitatively different initial magnetic field configurations: a field confined completely to the crust and another which penetrates the whole star. The aim of this study is to explore whether these magnetic structures in the crust can deform the star sufficiently to make it an observable source of gravitational waves. We model the evolution of these field configurations, which can develop, within $\sim 10^4$ -- $10^5$ yr, magnetic spots with local surface field strengths $\sim 10^{14}$ G maintained over $\gtrsim 10^6$ yr. Deformations caused by the magnetic forces are calculated. We show that, under favourable initial conditions, a star undergoing crustal Hall drift can have ellipticity $\epsilon\sim 10^{-6}$, even with sub-magnetar polar field strengths, after $\sim 10^5$ yr. A pulsar rotating at $\sim 10^2$ Hz with such $\epsilon$ is a promising gravitational-wave source candidate. Since such large deformations can be caused only by a particular magnetic field configuration that penetrates the whole star and whose maximum magnetic energy is concentrated in the outer core region, gravitational wave emission observed from radio pulsars can thus inform us about the internal field structures of young neutron stars.

The measured spin-down rates in quiescence of the transient accreting millisecond pulsars IGR J00291+5934, XTE J1751-305, SAX J1808.4-3658, and Swift J1756.9-2508 have been used to estimate the magnetic moments of these objects assuming standard magnetic dipole braking. It is shown that this approach leads to an overestimate, if the amount of residual accretion is enough to distort the magnetosphere away from a force-free configuration, through magnetospheric mass loading or crushing, so that the lever arm of the braking torque migrates inside the light cylinder. We derive an alternative spin-down formula and calculate the residual accretion rates where the formula is applicable. As a demonstration, we apply the alternative spin-down formula to produce updated magnetic moment estimates for the four objects above. We note that, based on current uncertain observations of quiescent accretion rates, magnetospheric mass loading and crushing are neither firmly indicated nor ruled out in these four objects. Because quiescent accretion rates are not measured directly (only upper limits are placed), it is impossible to be confident without more data whether the thresholds for magnetospheric mass loading or crushing are reached or not.

Quantum mechanical corrections to the hydromagnetic force balance equation, derived from the microscopic Schr\"{o}dinger-Pauli theory of quantum plasmas, modify the equilibrium structure and hence the mass quadrupole moment of a neutron star. It is shown here that the dominant effect --- spin paramagnetism --- is most significant in a magnetar, where one typically has $\mu_{B}|\boldsymbol{B}|\gtrsim k_B T_e$, where $\mu_{B}$ is the Bohr magneton, $\boldsymbol{B}$ is the magnetic field, and $T_e$ is the electron temperature. The spin paramagnetic deformation of a nonbarotropic magnetar with a linked poloidal-toroidal magnetic field is calculated to be up to ${{\sim 10}}$ times greater than the deformation caused solely by the Lorentz force. It depends on the degree of Pauli blocking by conduction electrons and the propensity to form magnetic domains, processes which are incompletely modelled at magnetar field strengths. The star becomes more oblate, as the toroidal field component strengthens. The result implies that existing classical predictions underestimate the maximum strength of the gravitational wave signal from rapidly spinning magnetars at birth. Turning the argument around, future gravitational-wave upper limits of increasing sensitivity will place ever-stricter constraints on the physics of Pauli blocking and magnetic domain formation under magnetar conditions.

The sudden spin-down event ('anti-glitch') observed in AXP 1E 2259+586 on 2012 April 21 was arguably caused by a decay of its internal toroidal magnetic field component, which turns a stable prolate configuration into an unstable one. We refine previous models of this process by modelling the star's magnetic field self-consistently as a 'twisted torus' configuration in non-barotropic equilibrium (which allows us to explore a greater range of equilibrium configurations). It is shown that, if the star's magnetic field is purely dipolar, the change in the toroidal field strength required to produce an anti-glitch of the observed size can be ~ 10 times larger than previously calculated. If the star has a quadrupolar magnetic field component, then an anti-glitch of similar magnitude can be produced via a decay of the quadrupole component, in addition to a decay of the toroidal component. We show that, if the quadrupole component decays, the minimum initial toroidal field strength and the change in toroidal field strength needed to produce the observed anti-glitch are lower than in the pure dipole twisted torus. In addition, we predict the maximum anti-glitch sizes, assuming that they are caused by a change in ellipticity, in four glitching magnetars and discuss the implications for energetics of accompanying X-ray bursts.

A recipe is presented to construct an analytic, self-consistent model of a non-barotropic neutron star with a poloidal-toroidal field of arbitrary multipole order, whose toroidal component is confined in a torus around the neutral curve inside the star, as in numerical simulations of twisted tori. The recipe takes advantage of magnetic-field-aligned coordinates to ensure continuity of the mass density at the surface of the torus. The density perturbation and ellipticity of such a star are calculated in general and for the special case of a mixed dipole-quadrupole field as a worked example. The calculation generalises previous work restricted to dipolar, poloidal-toroidal and multipolar, poloidal-only configurations. The results are applied, as an example, to magnetars whose observations (e.g., spectral features and pulse modulation) indicate that the internal magnetic fields may be at least one order of magnitude stronger than the external fields, as inferred from their spin downs, and are not purely dipolar.

Certain multi-wavelength observations of neutron stars, such as intermittent radio emissions from rotation-powered pulsars beyond the pair-cascade death line, the pulse profile of the magnetar SGR 1900+14 after its 1998 August 27 giant flare, and X-ray spectral features of PSR J0821-4300 and SGR 0418+5729, suggest that the magnetic fields of non-accreting neutron stars are not purely dipolar and may contain higher-order multipoles. Here, we calculate the ellipticity of a non-barotropic neutron star with (i) a quadrupole poloidal-toroidal field, and (ii) a purely poloidal field containing arbitrary multipoles, deriving the relation between the ellipticity and the multipole amplitudes. We present, as a worked example, a purely poloidal field comprising dipole, quadrupole, and octupole components. We show the correlation between field energy and ellipticity for each multipole, that the l=4 multipole has the lowest energy, and that l=5 has the lowest ellipticity. We show how a mixed multipolar field creates an observationally testable mismatch between the principal axes of inertia (to be inferred from gravitational wave data) and the magnetic inclination angle. Strong quadrupole and octupole components (with amplitudes ~ 10^2 times higher than the dipole) in SGR 0418+5729 still yield ellipticity ~ 10^{-8}, consistent with current gravitational wave upper limits. The existence of higher multipoles in fast-rotating objects (e.g., newborn magnetars) has interesting implications for the braking law and hence phase tracking during coherent gravitational wave searches.

Magnetic fields in upper main-sequence stars, white dwarfs, and neutron stars are known to persist for timescales comparable to their lifetimes. From a theoretical perspective this is problematic, as it can be shown that simple magnetic field configurations are always unstable. In non-barotropic stars, stable stratification allows for a much wider range of magnetic field structures than in barotropic stars, and helps stabilize them by making it harder to induce radial displacements. Recent simulations by Braithwaite and collaborators have shown that, in stably stratified stars, random initial magnetic fields evolve into nearly axisymmetric configurations with both poloidal and toroidal components, which then remain stable for some time. It is desirable to provide an analytic study of the stability of such fields. We write an explicit expression for a plausible equilibrium structure of an axially symmetric magnetic field with both poloidal and toroidal components of adjustable strengths, in a non-barotropic static fluid star, and study its stability using the energy principle. We construct a displacement field that should be a reasonable approximation to the most unstable mode of a toroidal field, and confirm Braithwaite's result that a given toroidal field can be stabilized by a poloidal field containing much less energy than the former. This is consistent with the speculation that the toroidal field is the main reservoir powering magnetar activity. The deformation of a neutron star caused by the hidden toroidal field can also cause emission of gravitational waves.

Recent calculations of the hydromagnetic deformation of a stratified, non-barotropic neutron star are generalized to describe objects with superconducting interiors, whose magnetic permeability \mu is much smaller than the vacuum value \mu_0. It is found that the star remains oblate if the poloidal magnetic field energy is \gtrsim 40% of total magnetic field energy, that the toroidal field is confined to a torus which shrinks as \mu decreases, and that the deformation is much larger (by a factor \sim \mu_0/\mu) than in a non-superconducting object. The results are applied to the latest direct and indirect upper limits on gravitational-wave emission from Laser Interferometer Gravitational Wave Observatory (LIGO) and radio pulse timing (spin-down) observations of 81 millisecond pulsars, to show how one can use these observations to infer the internal field strength. It is found that the indirect spin-down limits already imply astrophysically interesting constraints on the poloidal-toroidal field ratio and diamagnetic shielding factor (by which accretion reduces the observable external magnetic field, e.g. by burial). These constraints will improve following gravitational-wave detections, with implications for accretion-driven magnetic field evolution in recycled pulsars and the hydromagnetic stability of these objects' interiors.

A strong candidate for a source of gravitational waves is a highly magnetised, rapidly rotating neutron star (magnetar) deformed by internal magnetic stresses. We calculate the mass quadrupole moment by perturbing a zeroth-order hydrostatic equilibrium by an axisymmetric magnetic field with a \emph{linked poloidal-toroidal structure}. In this work, we do \emph{not} require the model star to obey a barotropic equation of state (as a realistic neutron star is not barotropic), allowing us to explore the hydromagnetic equilibria with fewer constraints. We derive the relation between the ratio of poloidal-to-total field energy $\Lambda$ and ellipticity $\epsilon$ and briefly compare our results to those obtained using the barotropic assumption. Then, we present some examples of how our results can be applied to astrophysical contexts. First, we show how our formulae, in conjunction with current gravitational wave (non-)detections of the Crab pulsar and the Cassiopeia A central compact object (Cas A CCO), can be used to constrain the strength of the internal toroidal fields of those objects. We find that, for the Crab pulsar (whose canonical equatorial dipole field strength, inferred from spin down, is $4\times 10^8$ T) to emit detectable gravitational radiation, the neutron star must have a strong toroidal field component, with maximum internal toroidal field strength $B_{\mathrm{tm}}=7\times 10^{12}$ T; for gravitational waves to be detected from the Cas A CCO at 300 Hz, $B_{\mathrm{tm}}\sim 10^{13}$ T, whereas detection at 100 Hz would require $B_{\mathrm{tm}}\sim 10^{14}$ T. Using our results, we also show how the gravitational wave signal emitted by a magnetar immediately after its birth (assuming it is born rapidly rotating, with $\Lambda\lesssim 0.2$) makes such a newborn magnetar a stronger candidate for gravitational wave detection than, for example, an SGR giant flare.

The existence of a tachocline in the Sun has been proven by helioseismology. It is unknown whether a similar shear layer, widely regarded as the seat of magnetic dynamo action, also exists in a protoneutron star. Sudden jumps in magnetic diffusivity $\eta$ and turbulent vorticity $\alpha$, for example at the interface between the neutron-finger and convective zones, are known to be capable of enhancing mean-field dynamo effects in a protoneutron star. Here we apply the well-known, plane-parallel, MacGregor-Charbonneau analysis of the Solar interfacial dynamo to the protoneutron star problem and calculate the growth rate analytically under a range of conditions. It is shown that, like the Solar dynamo, it is impossible to achieve self-sustained growth if the discontinuities in $\alpha$, $\eta$, and shear are coincident and the magnetic diffusivity is isotropic. In contrast, when the jumps in $\eta$ and $\alpha$ are situated away from the shear layer, self-sustained growth is possible for $P\lesssim 49.8$ ms (if the velocity shear is located at $0.3R$) or $P\lesssim 83.6$ ms (if the velocity shear is located at $0.6R$). This translates into stronger shear and/or $\alpha$-effect than in the Sun. Self-sustained growth is also possible if the magnetic diffusivity if anisotropic, through the ${\bf{\Omega}}\times{\bf{J}}$ effect, even when the $\alpha$, $\eta$, and shear discontinuities are coincident.

Recent numerical magnetohydrodynamic calculations by Braithwaite and collaborators support the `fossil field' hypothesis regarding the origin of magnetic fields in compact stars and suggest that the resistive evolution of the fossil field can explain the reorganisation and decay of magnetar magnetic fields. Here, these findings are modelled analytically by allowing the stellar magnetic field to relax through a quasistatic sequence of nonaxisymmetric, force-free states, by analogy with spheromak relaxation experiments, starting from a random field. Under the hypothesis that the force-free modes approach energy equipartition in the absence of resistivity, the output of the numerical calculations is semiquantitatively recovered: the field settles down to a linked poloidal-toroidal configuration, which inflates and becomes more toroidal as time passes. A qualitatively similar (but not identical) end state is reached if the magnetic field evolves by exchanging helicity between small and large scales according to an $\alpha$-dynamo-like, mean-field mechanism, arising from the fluctuating electromotive force produced by the initial random field. The impossibility of matching a force-free internal field to a potential exterior field is discussed in the magnetar context.

A recent laboratory experiment (Blaauwgeers et al. 2003) suggests that a Kelvin-Helmholtz (KH) instability at the interface between two superfluids, one rotating and anisotropic, the other stationary and isotropic, may trigger sudden spin-up of the stationary superfluid. This result suggests that a KH instability at the crust-core ($^1S_0$-$^3P_2$-superfluid) boundary of a neutron star may provide a trigger mechanism for pulsar glitches. We calculate the dispersion relation of the KH instability involving two different superfluids including the normal fluid components and their effects on stability, particularly entropy transport. We show that an entropy difference between the core and crust superfluids reduces the threshold differential shear velocity and threshold crust-core density ratio. We evaluate the wavelength of maximum growth of the instability for neutron star parameters and find the resultant circulation transfer to be within the range observed in pulsar glitches.