We present three-dimensional force-free electrodynamics simulations of magnetar magnetospheres that demonstrate the instability of certain degenerate, high energy equilibrium solutions of the Grad-Shafranov equation. This result indicates the existence of an unstable branch of twisted magnetospheric solutions and allows to formulate an instability criterion. The rearrangement of magnetic field lines as a consequence of this instability triggers the dissipation of up to 30% of the magnetospheric energy on a thin layer above the magnetar surface. During this process, we predict an increase of the mechanical stresses onto the stellar crust, which can potentially result in a global mechanical failure of a significant fraction of it. We find that the estimated energy release and the emission properties are compatible with the observed giant flare events. The newly identified instability is a candidate for recurrent energy dissipation, which could explain part of the phenomenology observed in magnetars.

We present a detailed analysis of the properties of twisted, force-free magnetospheres of non-rotating neutron stars, which are of interest in the modelling of magnetar properties and evolution. In our models the magnetic field smoothly matches to a current-free (vacuum) solution at some large external radius, and they are specifically built to avoid pathological surface currents at any of the interfaces. By exploring a large range of parameters, we find a few remarkable general trends. We find that the total dipolar moment can be increased by up to $40\%$ with respect to a vacuum model with the same surface magnetic field, due to the contribution of magnetospheric currents to the global magnetic field. Thus, estimates of the surface magnetic field based on the large-scale dipolar braking torque are slightly overestimating the surface value by the same amount. Consistently, there is a moderate increase in the total energy of the model with respect to the vacuum solution of up to $25\%$, which would be the available energy budget in the event of a fast, global magnetospheric reorganization commonly associated with magnetar flares. We have also found the interesting result of the existence of a critical twist ($\varphi_{\rm max} \lesssim 1.5$ rad), beyond which we cannot find any more numerical solutions. Combining the models considered in this paper with the evolution of the interior of neutron stars will allow us to study the influence of the magnetosphere on the long-term magnetic, thermal, and rotational evolution.

We extend our previous study of equilibrium solutions of non-rotating force-free magnetospheres of neutron stars. We show that multiple solutions exist for the same sets of parameters, implying that the solutions are degenerate. We are able to obtain configurations with disconnected field lines, however, in nearly all cases these correspond to degenerate higher energy solutions. We carry out a wide parametric search in order to understand the properties of the solutions. We confirm our previous results that the lower energy solutions have up to $\sim 25\%$ more energy than the vacuum case, helicity of the order of $\sim 5$ (in some defined units), maximum twist of $\sim 1.5$ rad, and a dipole strength that is up to $\sim 40\%$ larger than the vacuum dipole. Including the degenerate higher energy solutions allows for larger theoretical limits of up to $\sim 80\%$ more energy with respect to the vacuum case, helicity of the order of $\sim 8$, and a dipole strength that can now be up to four times that of the vacuum dipole, while the twist can be significantly larger and even diverge for configurations with disconnected domains. The higher energy solutions are probably unstable, therefore, it is unlikely that such magnetospheres exist under normal conditions in magnetars and high magnetic field pulsars.

We study the long-term quasi-steady evolution of the force-free magnetosphere of a magnetar coupled to its internal magnetic field. We find that magnetospheric currents can be maintained on long timescales of the order of thousands of years. Meanwhile, the energy, helicity and twist stored in the magnetosphere all gradually increase over the course of this evolution, until a critical point is reached, beyond which a force-free magnetosphere cannot be constructed. At this point, some large-scale magnetospheric rearrangement, possibly resulting in an outburst or a flare, must occur, releasing a large fraction of the stored energy, helicity and twist. After that, the quasi-steady evolution should continue in a similar manner from the new initial conditions. The timescale for reaching this critical point depends on the overall magnetic field strength and on the relative fraction of the toroidal field. The energy stored in the force-free magnetosphere is found to be up to $\sim 30\%$ larger than the corresponding vacuum energy. This implies that for a $10^{14}$ G field at the pole, the energy budget available for fast magnetospheric events is of the order of a few $10^{44}$ erg. The spindown rate is estimated to increase by up to $\sim 60\%$, since the dipole content in the magnetosphere is enhanced by the currents present there. A rough estimate of the braking index $n$ reveals that it is systematically $n < 3$ for the most part of the evolution, consistent with actual measurements for pulsars and early estimates for several magnetars.

We study the coupling of the force-free magnetosphere to the long-term internal evolution of a magnetar. We allow the relation between the poloidal and toroidal stream functions - that characterizes the magnetosphere - to evolve freely without constraining its particular form. We find that, on time-scales of the order of kyr, the energy stored in the magnetosphere gradually increases, as the toroidal region grows and the field lines expand outwards. This continues until a critical point is reached beyond which force-free solutions for the magnetosphere can no longer be constructed, likely leading to some large-scale magnetospheric reorganization. The energy budget available for such events can be as high as several $10^{45}\,$erg for fields of $10^{14}\,$G. Subsequently, starting from the new initial conditions, the evolution proceeds in a similar manner. The time-scale to reach the critical point scales inversely with the magnetic field amplitude. Allowing currents to pass through the last few meters below the surface, where the magnetic diffusivity is orders of magnitude larger than in the crust, should give rise to a considerable amount of energy deposition through Joule heating. We estimate that the effective surface temperature could increase locally from $\sim 0.1\,$keV to $\sim 0.3 - 0.6\,$keV, in good agreement with observations. Similarly, the power input from the interior into the magnetosphere could be as high as $10^{35} - 10^{36}\,$erg/s, which is consistent with peak luminosities observed during magnetar outbursts. Therefore, a detailed treatment of currents flowing through the envelope may be needed to explain the thermal properties of magnetars.

We revisit the various approximations employed to study the long-term evolution of the magnetic field in neutron star cores and discuss their limitations and possible improvements. A recent controversy on the correct form of the induction equation and the relevant evolution timescale in superconducting neutron star cores is addressed and clarified. We show that this ambiguity in the estimation of timescales arises as a consequence of nominally large terms that appear in the induction equation, but which are, in fact, mostly irrotational. This subtlety leads to a discrepancy by many orders of magnitude when velocity fields are absent or ignored. Even when internal velocity fields are accounted for, only the solenoidal part of the electric field contributes to the induction equation, which can be substantially smaller than the irrotational part. We also argue that stationary velocity fields must be incorporated in the slow evolution of the magnetic field as the next level of approximation.

We study ambipolar diffusion in strongly magnetised neutron stars, with special focus on the effects of neutrino reaction rates and the impact of a superfluid/superconducting transition in the neutron star core. For axisymmetric magnetic field configurations, we determine the deviation from $\beta-$equilibrium induced by the magnetic force and calculate the velocity of the slow, quasi-stationary, ambipolar drift. We study the temperature dependence of the velocity pattern and clearly identify the transition to a predominantly solenoidal flow. For stars without superconducting/superfluid constituents and with a mixed poloidal-toroidal magnetic field of typical magnetar strength, we find that ambipolar diffusion proceeds fast enough to have a significant impact on the magnetic field evolution only at low core temperatures, $T \lesssim 1-2\times10^8$ K. The ambipolar diffusion timescale becomes appreciably shorter when fast neutrino reactions are present, because the possibility to balance part of the magnetic force with pressure gradients is reduced. We also find short ambipolar diffusion timescales in the case of superconducting cores for $T \lesssim 10^9$ K, due to the reduced interaction between protons and neutrons. In the most favourable scenario, with fast neutrino reactions and superconducting cores, ambipolar diffusion results in advection velocities of several km/kyr. This velocity can substantially reorganize magnetic fields in magnetar cores, in a way that can only be confirmed by dynamical simulations.

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.

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.

In 1977, Flowers and Ruderman described a perturbation that destabilises a purely dipolar magnetic field in a fluid star. They considered the effect of cutting the star in half along a plane containing the symmetry axis and rotating each half by $90\degr$ in opposite directions, which would cause the energy of the magnetic field in the exterior of the star to be greatly reduced, just as it happens with a pair of aligned magnets. We formally solve for the energy of the external magnetic field and check that it decreases monotonously along the entire rotation. We also describe the instability using perturbation theory, and see that it happens due to the work done by the interaction of the magnetic field with surface currents. Finally, we consider the stabilising effect of adding a toroidal field by studying the potential energy perturbation when the rotation is not done along a sharp cut, but with a continuous displacement field that switches the direction of rotation across a region of small but finite width. Using these results, we estimate the relative strengths of the toroidal and poloidal field needed to make the star stable to this displacement and see that the energy of the toroidal field required for stabilisation is much smaller than the energy of the poloidal field. We also show that, contrary to a common argument, the Flowers-Ruderman instability cannot be applied many times in a row to reduce the external magnetic energy indefinitely.

We review the interaction in intermediate and high mass stars between their evolution and magnetic and chemical properties. We describe the theory of Ap-star `fossil' fields, before touching on the expected secular diffusive processes which give rise to evolution of the field. We then present recent results from a spectropolarimetric survey of Herbig Ae/Be stars, showing that magnetic fields of the kind seen on the main-sequence already exist during the pre-main sequence phase, in agreement with fossil field theory, and that the origin of the slow rotation of Ap/Bp stars also lies early in the pre-main sequence evolution; we also present results confirming a lack of stars with fields below a few hundred gauss. We then seek which macroscopic motions compete with atomic diffusion in determining the surface abundances of AmFm stars. While turbulent transport and mass loss, in competition with atomic diffusion, are both able to explain observed surface abundances, the interior abundance distribution is different enough to potentially lead to a test using asterosismology. Finally we review progress on the turbulence-driving and mixing processes in stellar radiative zones.