Understanding the evolution of the angle $\chi$ between a magnetar's rotation and magnetic axes sheds light on the star's birth properties. This evolution is coupled with that of the stellar rotation $\Omega$, and depends on the competing effects of internal viscous dissipation and external torques. We study this coupled evolution for a model magnetar with a strong internal toroidal field, extending previous work by modelling -- for the first time in this context -- the strong proto-magnetar wind acting shortly after birth. We also account for the effect of buoyancy forces on viscous dissipation at late times. Typically we find that $\chi\to 90^\circ$ shortly after birth, then decreases towards $0^\circ$ over hundreds of years. From observational indications that magnetars typically have small $\chi$, we infer that these stars are subject to a stronger average exterior torque than radio pulsars, and that they were born spinning faster than $\sim 100-300$ Hz. Our results allow us to make quantitative predictions for the gravitational and electromagnetic signals from a newborn rotating magnetar. We also comment briefly on the possible connection with periodic Fast Radio Burst sources.

In the current work we explore the applicability of standard theoretical models of accretion to the observed properties of M51 ULX-7. The spin-up rate and observed X-ray luminosity are evidence of a neutron star with a surface magnetic field of $2-7\times10^{13}$ G, rotating near equilibrium. Analysis of the X-ray light-curve of the system (Swift/XRT data) reveals the presence of a $\sim$39 d super-orbital period. We argue that the super-orbital periodicity is due to disc precession, and that material is accreted onto the neutron star at a constant rate throughout it. Moreover, by attributing this modulation to the free precession of the neutron star we estimate a surface magnetic field strength of $3-4\times10^{13}$ G. The agreement of these two independent estimates provide strong constraints on the surface polar magnetic field strength of the neutron star.

Under normal conditions in a neutron-star crust, ions are locked in place in the crustal lattice and only electrons are mobile, and magnetic-field evolution is thus directly related to the electron velocity. The evolution, however, builds magnetic stresses that can become sufficiently large for the crust to exceed its elastic limit, and to flow plastically. We consider the nature of this plastic flow and the back-reaction on the crustal magnetic field evolution. We formulate a plane-parallel model for the local failure, showing that surface motions are inevitable once the crust yields, in the absence of extra dissipative mechanisms. We perform numerical evolutions of the crustal magnetic field under the joint effect of Hall drift and Ohmic decay, tracking the build-up of magnetic stresses, and diagnosing crustal failure with the von Mises criterion. Beyond this point we solve for the coupled evolution of the plastic velocity and magnetic field. Our results suggest that to have a coexistence of a magnetar corona with small-scale magnetic features, the viscosity of the plastic flow must be roughly $10^{36}-10^{37}\ \textrm{g cm}^{-1}\textrm{s}^{-1}$. We find significant motion at the surface at a rate of $10-100$ centimetres per year, and that the localised magnetic field is weaker than in evolutions without plastic flow. We discuss astrophysical implications, and how our local simulations could be used to build a global model of field evolution in the neutron-star crust.

A rotating fluid star, endowed with a magnetic field, can undergo a form of precessional motion: a sum of rigid-body free precession and a non-rigid response. On secular timescales this motion is dissipated by bulk and shear viscous processes in the stellar interior and magnetospheric braking in the exterior, changing the inclination angle between the rotation and magnetic axes. Using our recent solutions for the non-rigid precessional dynamics, and viscous dissipation integrals derived in this paper, we make the only self-consistent calculation to date of these dissipation rates. We present the first results for the full coupled evolution of spindown and inclination angle for a model of a late-stage proto-neutron star with a strong toroidal magnetic field, allowing for both electromagnetic torques and internal dissipation when evolving the inclination angle. We explore this coupled evolution for a range of initial inclination angles, rotation rates and magnetic field strengths. For fixed initial inclination angle, our results indicate that the neutron-star population naturally evolves into two classes: near-aligned and near-orthogonal rotators -- with typical pulsars falling into the latter category. Millisecond magnetars can evolve into the near-aligned rotators which mature magnetars appear to be, but only for small initial inclination angle and internal toroidal fields stronger than roughly $5\times 10^{15}$ G. Once any model has evolved to either an aligned or orthogonal state, there appears to be no further evolution away from this state at later times.

In the standard form of the relativistic heat equation used in astrophysics, information propagates instantaneously, rather than being limited by the speed of light as demanded by relativity. We show how this equation nonetheless follows from a more general, causal theory of heat propagation in which the entropy plays the role of a fluid. In deriving this result, however, we see that it is necessary to make some assumptions which are not universally valid: the dynamical timescales of the process must be long compared with the explicitly causal physics of the theory, the heat flow must be sufficiently steady, and the spacetime static. Generalising the heat equation (e.g. restoring causality) would thus entail retaining some of the terms we neglected. As a first extension, we derive the heat equation for the spacetime associated with a slowly-rotating star or black hole, showing that it only differs from the static result by an additional advection term due to the rotation, and as a consequence demonstrate that a hotspot on a neutron star will be seen to be modulated at the rotation frequency by a distant observer.

Stars are, generically, rotating and magnetised objects with a misalignment between their magnetic and rotation axes. Since a magnetic field induces a permanent distortion to its host, it provides effective rigidity even to a fluid star, leading to bulk stellar motion which resembles free precession. This bulk motion is however accompanied by induced interior velocity and magnetic field perturbations, which are oscillatory on the precession timescale. Extending previous work, we show that these quantities are described by a set of second-order perturbation equations featuring cross-terms scaling with the product of the magnetic and centrifugal distortions to the star. For the case of a background toroidal field, we reduce these to a set of differential equations in radial functions, and find a method for their solution. The resulting magnetic-field and velocity perturbations show complex multipolar structure and are strongest towards the centre of the star.

The activity of magnetars is believed to be powered by colossal magnetic energy reservoirs. We sketch an evolutionary picture in which internal field evolution in magnetars generates a twisted corona, form which energy may be released suddenly in a single giant flare, or more gradually through smaller outbursts and persistent emission. Given the ages of magnetars and the energy of their giant flares, we suggest that their evolution is driven by a novel mechanism: magnetic flux transport/decay due to persistent plastic flow in the crust, which would invalidate the common assumption that the crustal lattice is static and evolves only under Hall drift and Ohmic decay. We estimate the field strength required to induce plastic flow as a function of crustal depth, and the viscosity of the plastic phase. The star's superconducting core may also play a role in magnetar field evolution, depending on the star's spindown history and how rotational vortices and magnetic fluxtubes interact.

We construct equilibrium configurations of magnetized, two-fluid neutron stars using an iterative numerical method. Working in Newtonian framework we assume that the neutron star has two regions: the core, which is modelled as a two-component fluid consisting of type-II superconducting protons and superfluid neutrons, and the crust, a region composed of normal matter. Taking a new step towards more complete equilibrium models, we include the effect of entrainment, which implies that a magnetic force acts on neutrons, too. We consider purely poloidal field cases and present improvements to an earlier numerical scheme for solving equilibrium equations, by introducing new convergence criteria. We find that entrainment results in qualitative differences in the structure of field lines along the magnetic axis.

The presence of superconducting and superfluid components in the core of mature neutron stars calls for the rethinking of a number of key magnetohydrodynamical notions like resistivity, the induction equation, magnetic energy and flux-freezing. Using a multi-fluid magnetohydrodynamics formalism, we investigate how the magnetic field evolution is modified when neutron star matter is composed of superfluid neutrons, type-II superconducting protons and relativistic electrons. As an application of this framework, we derive an induction equation where the resistive coupling originates from the mutual friction between the electrons and the vortex/fluxtube arrays of the neutron and proton condensates. The resulting induction equation allows the identification of two timescales that are significantly different from those of standard magnetohydrodynamics. The astrophysical implications of these results are briefly discussed.

Crustquake events may be connected with both rapid spin-up `glitches' within the regular slowdown of neutron stars, and high-energy magnetar flares. We argue that magnetic field decay builds up stresses in a neutron star's crust, as the elastic shear force resists the Lorentz force's desire to rearrange the global magnetic-field equilibrium. We derive a criterion for crust-breaking induced by a changing magnetic-field configuration, and use this to investigate strain patterns in a neutron star's crust for a variety of different magnetic-field models. Universally, we find that the crust is most liable to break if the magnetic field has a strong toroidal component, in which case the epicentre of the crustquake is around the equator. We calculate the energy released in a crustquake as a function of the fracture depth, finding that it is independent of field strength. Crust-breaking is, however, associated with a characteristic local field strength of $2.4\times 10^{14}$ G for a breaking strain of $0.001$, or $2.4\times 10^{15}$ G at a breaking strain of $0.1$. We find that even the most luminous magnetar giant flare could have been powered by crustal energy release alone.

We study equilibrium magnetic field configurations in a neutron star (NS) whose core has type-II superconducting protons. Unlike the equations for normal matter, which feature no special field strength, those for superconductors contain the lower critical field, of order $10^{15}$ G. We find that the ratio of this critical field to the smooth-averaged stellar field at the crust-core boundary is the key feature dictating the field geometry. Our results suggest that pulsar and magnetar-strength fields have notably different configurations. Field decay for NSs with $B_\textrm{pole}\sim 10^{14}$ G could thus result in substantial internal rearrangements, pushing the toroidal field component out of the core; this may be related to observed magnetar activity. In addition, we calculate the magnetically-induced ellipticities of our models.

We study the time-evolution of axisymmetric oscillations of superfluid magnetars with a poloidal magnetic field and an elastic crust, working in Newtonian gravity. Extending earlier models, we study the effects of composition gradients and entrainment on the magneto-elastic wave spectrum and on the potential identification of the observed quasi periodic oscillations (QPOs). We use two-fluid polytropic equations of state to construct our stellar models, which mimic realistic composition gradient configurations. The basic features of the axial axisymmetric spectrum of normal fluid stars are reproduced by our results and in addition we find several magneto-elastic waves with a mixed character. In the core, these oscillations mimic the shear mode pattern of the crust as a result of the strong dynamical coupling between these two regions. Incorporating the most recent entrainment configurations in our models, we find that they have a double effect on the spectrum: the magnetic oscillations of the core have a frequency enhancement, while the mixed magneto-elastic waves originating in the crust are moved towards the frequencies of the single-fluid case. The distribution of lower-frequency magneto-elastic oscillations for our models is qualitatively similar to the observed magnetar QPOs with $\nu < 155 $Hz. \emph{In particular, some of these QPOs could represent mixed magneto-elastic oscillations with frequencies not greatly different from the crustal modes of an unmagnetised star.} We find that many QPOs could even be accounted for using a model with a relatively weak polar field of $B_{p} \simeq 3\cdot 10^{14}$G, because of the superfluid enhancement of magnetic oscillations. Finally, we discuss the possible identification of 625 and 1837~Hz QPOs either with non-axisymmetric modes or with high-frequency axisymmetric QPOs excited by crustal mode overtones.

We construct hydromagnetic neutron star equilibria which allow for a non-zero electric current distribution in the exterior. The novelty of our models is that the neutron star's interior field is in equilibrium with its magnetosphere, thus bridging the gap between previous work in this area which either solves for the interior assuming a vacuum exterior or solves for the magnetosphere without modelling the star itself. We consider only non-rotating stars in this work, so our solutions are most immediately applicable to slowly-rotating systems such as magnetars. Nonetheless, we demonstrate that magnetospheres qualitatively resembling those expected for both magnetars and pulsars are possible within our framework. The "inside-out" approach taken in this paper should be more generally applicable to rotating neutron stars, where the interior and exterior regions are again not independent but evolve together.

The interior of a neutron star is likely to be predominantly a mixture of superfluid neutrons and superconducting protons. This results in the quantisation of the star's magnetic field into an array of thin fluxtubes, producing a macroscopic force very different from the Lorentz force of normal matter. We show that in an axisymmetric superconducting equilibrium the behaviour of a magnetic field is governed by a single differential equation. Solving this, we present the first self-consistent superconducting neutron star equilibria with poloidal and mixed poloidal-toroidal fields, also giving the first quantitative results for the corresponding magnetically-induced distortions to the star. The poloidal component is dominant in all our configurations. We suggest that the transition from normal to superconducting matter in a young neutron star may cause a large-scale field rearrangement.

The violent giant flares of magnetars excite QPOs which persist for hundreds of seconds, as seen in the X-ray tail following the initial burst. Recent studies, based on single-fluid barotropic magnetar models, have suggested that the lower-frequency QPOs correspond to magneto-elastic oscillations of the star. The higher frequencies, however, in particular the strong 625 Hz peak, have proved harder to explain, except as high mode multipoles. In this work we study the time evolutions of non-axisymmetric oscillations of two-fluid Newtonian magnetars with no crust. We consider models with superfluid neutrons and normal protons, and poloidal and toroidal background field configurations. We show that multi-fluid physics (composition-gradient stratification, entrainment) tends to increase Alfven mode frequencies significantly from their values in a single-fluid barotropic model. The higher-frequency magnetar QPOs may then be naturally interpreted as Alfven oscillations of the multi-fluid stellar core. The lower-frequency QPOs are less easily explained within our purely fluid core model, but we discuss the possibility that these are crustal modes.

We construct barotropic stellar equilibria, containing magnetic fields with both poloidal and toroidal field components. We extend earlier results by exploring the effect of different magnetic field and current distributions. Our results suggest that the boundary treatment plays a major role in whether the poloidal or toroidal field component is globally dominant. Using time evolutions we provide the first stability test for mixed poloidal-toroidal fields in barotropic stars, finding that all these fields suffer instabilities due to one of the field components: these are localised around the pole for toroidal-dominated equilibria and in the closed-field line region for poloidal-dominated equilibria. Rotation provides only partial stabilisation. There appears to be very limited scope for the existence of stable magnetic fields in barotropic stars. We discuss what additional physics from real stars may allow for stable fields.

We construct two-fluid equilibrium configurations for neutron stars with magnetic fields, using a self-consistent and nonlinear numerical approach. The two-fluid approach - likely to be valid for large regions of all but the youngest NSs - provides us with a straightforward way to introduce stratification and allows for more realistic models than the ubiquitous barotropic assumption. In all our models the neutrons are modelled as a superfluid, whilst for the protons we consider two cases: one where they are a normal fluid and another where they form a type-II superconductor. We consider a variety of field configurations in the normal-proton case and purely toroidal fields in the superconducting case. We find that stratification allows for a stronger toroidal component in mixed-field configurations, though the poloidal component remains the largest in all our models. We provide quantitative results for magnetic ellipticities of NSs, both in the normal- and superconducting-proton cases.

Traditionally, the subject of hydromagnetic equilibrium in neutron stars has been addressed in the context of standard magnetohydrodynamics, with matter obeying a barotropic equation of state. In this paper we take a step towards a more realistic treatment of the problem by considering neutron stars with interior superfluid components. In this multifluid model stratification associated with a varying matter composition (the relative proton to neutron density fraction) enters as a natural ingredient, leading to a non-barotropic system. After formulating the hydromagnetic equilibrium of superfluid/superconducting neutron stars as a perturbation problem, we focus on the particular case of a three-fluid system consisting of superfluid neutrons and normal protons and electrons. We determine the equilibrium structure of dipolar magnetic fields with a mixed poloidal-toroidal composition. We find that, with respect to barotropic models, stratification has the generic effect of leading to equilibria with a higher fraction of magnetic energy stored in the toroidal component. However, even in models with strong stratification the poloidal and toroidal components are comparable, with the former contributing the bulk of the magnetic energy.

The problem of the stability of magnetic fields in stars has a long history and has been investigated in detail in perturbation theory. Here we consider the nonlinear evolution of a non-rotating neutron star with a purely poloidal magnetic field, in general relativity. We find that an instability develops in the region of the closed magnetic field lines and over an Alfven timescale, as predicted by perturbation theory. After the initial unstable growth, our evolutions show that a toroidal magnetic field component is generated, which increases until it is locally comparable in strength with the poloidal one. On longer timescales the system relaxes to a new non-axisymmetric configuration with a reorganization of the stellar structure and large-amplitude oscillations, mostly in the fundamental mode. We discuss the energies involved in the instability and the impact they may have on the phenomenology of magnetar flares and on their detectability through gravitational-wave emission.

We study the time evolution of non-axisymmetric linear perturbations of a rotating magnetised neutron star, whose magnetic field is multipolar and purely poloidal. The background stellar configurations are generated self-consistently, allowing for distortions to the density distribution from rotational and magnetic forces. We find that the behaviour of axial-led perturbations is dominated by an instability generic to poloidal fields, which is localised around the `neutral line' where the background field vanishes. Rotation acts to reduce the effect of this instability. Polar-led perturbations do not appear to be unstable and in this case we find global Alfv\'en modes, whose restoring force is the magnetic field. In a rotating magnetised star there are no pure Alfv\'en modes or pure inertial modes, but hybrids of these. We discuss the nature of magnetic instabilities and oscillations in magnetars and pulsars, finding the dominant Alfv\'en mode has a frequency comparable with observed magnetar QPOs.

We study $m=1$ oscillations and instabilities of magnetised neutron stars, by numerical time-evolution of linear perturbations of the system. The background stars are stationary equilibrium configurations with purely toroidal magnetic fields. We find that an $m=1$ instability of toroidal magnetic fields, already known from local analyses, may also be found in our relatively low-resolution global study. We present quantitative results for the instability growth rate and its suppression by rotation. The instability is discussed as a possible trigger mechanism for Soft Gamma Repeater (SGR) flares. Although our primary focus is evolutions of magnetised stars, we also consider perturbations about unmagnetised background stars in order to study $m=1$ inertial modes. We track these modes up to break-up frequency, extending known slow-rotation results.

We investigate the oscillation spectrum of rotating Newtonian neutron stars endowed with purely toroidal magnetic fields, using a time evolution code to evolve linear perturbations in the Cowling approximation. The background star is generated by numerically solving the MHD equilibrium equations and may be nonspherical by virtue of both rotation and magnetic effects; hence our perturbations and background are fully consistent. Whilst the background field is purely toroidal, the perturbed field is mixed poloidal-toroidal. From Fourier analysis of the perturbations we are able to identify a number of magnetically-restored Alfv\'en (or $a$-) modes. We show that in a rotating star pure inertial and $a$-modes are replaced by hybrid magneto-inertial modes, which reduce to $a$-modes in the nonrotating limit and inertial modes in the nonmagnetic limit. We show that the $r$-mode instability is suppressed by magnetic fields in sufficiently slowly rotating stars. In addition, we determine magnetic frequency shifts in the $f$-mode. We discuss the astrophysical relevance of our results, in particular for magnetar oscillations.

We derive general equations for axisymmetric Newtonian MHD and use these as the basis of a code for calculating equilibrium configurations of rotating magnetised neutron stars in a stationary state. We investigate the field configurations that result from our formalism, which include purely poloidal, purely toroidal and mixed fields. For the mixed-field formalism the toroidal component appears to be bounded at less than 7%. We calculate distortions induced both by magnetic fields and by rotation. From our non-linear work we are able to look at the realm of validity of perturbative work: we find for our results that perturbative-regime formulae for magnetic distortions agree to within 10% of the nonlinear results if the ellipticity is less than 0.15 or the average field strength is less than $10^{17}$ G. We also consider how magnetised equilibrium structures vary for different polytropic indices.