In the solid crusts of neutron stars, the advection of the magnetic field by the current-carrying electrons, an effect known as Hall drift, should play a very important role as the ions remain essentially fixed (as long as the solid does not break). Although Hall drift preserves the magnetic field energy, it has been argued that it may drive a turbulent cascade to scales at which Ohmic dissipation becomes effective, allowing a much faster decay in objects with very strong fields. On the other hand, it has been found that there are "Hall equilibria", i.e., field configurations that are unaffected by Hall drift. Here, we address the crucial question of the stability of these equilibria through axially symmetric (2D) numerical simulations of Hall drift and Ohmic diffusion, with the simplifying assumption of uniform electron density and conductivity. We demonstrate the 2D-stability of a purely poloidal equilibrium, for which Ohmic dissipation makes the field evolve towards an attractor state through adjacent stable configurations, around which damped oscillations occur. For this field, the decay scales with the Ohmic timescale. We also study the case of an unstable equilibrium consisting of both poloidal and toroidal field components that are confined within the crust. This field evolves into a stable configuration, which undergoes damped oscillations superimposed on a slow evolution towards an attractor, just as the purely poloidal one.

Magnetic protection of potentially habitable planets plays a central role in determining their actual habitability and/or the chances of detecting atmospheric biosignatures. We develop here a thermal evolution model of potentially habitable Earth-like planets and super-Earths. Using up-to-date dynamo scaling laws we predict the properties of core dynamo magnetic fields and study the influence of thermal evolution on their properties. The level of magnetic protection of tidally locked and unlocked planets is estimated by combining simplified models of the planetary magnetosphere and a phenomenological description of the stellar wind. Thermal evolution introduces a strong dependence of magnetic protection on planetary mass and rotation rate. Tidally locked terrestrial planets with an Earth-like composition would have early dayside magnetospause distances between 1.5 and 4.0 Rp, larger than previously estimated. Unlocked planets with periods of rotation ~1 day are protected by magnetospheres extending between 3 and 8 Rp. Our results are robust against variations in planetary bulk composition and uncertainties in other critical model parameters. For illustration purposes the thermal evolution and magnetic protection of the potentially habitable super-Earths GL 581d, GJ 667Cc and HD 40307g were also studied. Assuming an Earth-like composition we found that the dynamos of these planets are already extinct or close to being shut down. While GL 581d is the best protected, the protection of HD 40307g cannot be reliably estimated. GJ 667Cc, even under optimistic conditions, seems to be severely exposed to the stellar wind and, under the conditions of our model, has probably suffered massive atmospheric losses.

We present a model for the evolution of the magnetic properties of habitable terrestrial planets and their effects on the protection of planetary atmosphere against the erosive action of stellar wind. Using up-to-date thermal evolution models and dynamo scaling laws we predict the evolution of the planetary dipole moment as a function of planetary mass and rotation rate. Combining these results with models for the evolution of the stellar wind, stellar XUV fluxes and planetary exosphere characteristics, we determine the properties of the magnetosphere and the exobase radius in order to estimate the level of atmospheric mass losses. We use this model to evaluate the magnetic protection of the potentially habitable super-Earths GJ 667Cc, Gl 581d and HD 85512b. We confirm that Earth-like planets, even under the highest attainable magnetic field strengths, will lose a significant fraction of their atmospheric volatiles if they are tidally locked in the habitable zone of dM stars, or even if having N/O-rich atmospheres they are in habitable zones closer than $\sim$ 0.8 AU. Similar mass-dependent inner limits have been found for super-Earths $M_p\gtrsim 3 M_\oplus$ that in any case seem to have better chances of preserving their atmospheres even if they are tidally locked. We predict that the atmosphere of GJ 667Cc has probably already been obliterated and it is presently uninhabitable. On the other hand, our model predicts that the atmospheres of Gl 581d and HD 85512b would be well protected by intrinsic magnetic fields, even under the worst expected conditions of stellar aggression. (abrigded abstract).

We study the effect of the non-linear process of ambipolar diffusion (joint transport of magnetic flux and charged particles relative to neutral particles) on the long-term behavior of a non-uniform magnetic field in a one-dimensional geometry. Our main focus is the dissipation of magnetic energy inside neutron stars(particularly magnetars), but our results have a wider application, particularly to the interstellar medium and the loss of magnetic flux from collapsing molecular cloud cores. Our system is a weakly ionized plasma in which neutral and charged particles can be converted into each other through nuclear beta decays (or ionization-recombination processes). In the "weak-coupling" limit of infrequent inter-particle interactions, the evolution of the magnetic field is controlled by the beta decay rate and can be described by a non-linear partial integro-differential equation. In the opposite, "strong-coupling" regime, the evolution is controlled by the inter-particle collisions and can be modelled through a non-linear diffusion equation. We show numerically that, in both regimes, ambipolar diffusion tends to spread out the magnetic flux, but, contrary to the normal Ohmic diffusion, it produces sharp magnetic field gradients with associated current sheets around those regions where the magnetic field is weak.

This paper is the first in a series aimed at understanding the long-term evolution of neutron star magnetic fields. We model the stellar matter as an electrically neutral and lightly ionized plasma composed of three moving particle species: neutrons, protons, and electrons, which can be converted into each other by weak interactions (beta decays), suffer binary collisions, and be affected by each other's macroscopic electromagnetic fields. Since the evolution of the magnetic field occurs over thousands of years or more, compared to dynamical time scales (sound and Alfv\'en) of milliseconds to seconds, we use a slow-motion approximation in which we neglect the inertial terms in the equations of motion for the particles. We restrict ourselves to a one-dimensional geometry in which the magnetic field points in one Cartesian direction but varies only along an orthogonal direction. We study the evolution of the system in three different ways: (i) estimating time scales directly from the equations, guided by physical intuition; (ii) a normal-mode analysis in the limit of a nearly uniform system; and (iii) a finite-difference numerical integration of the equations of motion. We find good agreement between our analytical normal-mode solutions and the numerical simulations. We show that the magnetic field and the particles evolve through successive quasi-equilibrium states, on time scales that can be understood by physical arguments. Depending of the parameter values the magnetic field can evolve by ohmic diffusion or by ambipolar diffusion, the latter being limited either by interparticle collisions or by relaxation to chemical equilibrium through beta decays. The numerical simulations are further validated by verifying that they satisfy the known conservation laws also in highly non-linear situations.

Using a numerical simulation, we study the effects of ambipolar diffusion and ohmic diffusion on the magnetic field evolution in the interior of an isolated neutron star. We are interested in the behavior of the magnetic field on a long time scale, over which all Alfven and sound waves have been damped. We model the stellar interior as an electrically neutral plasma composed of neutrons, protons and electrons, which can interact with each other through collisions and electromagnetic forces. Weak interactions convert neutrons and charged particles into each other, erasing chemical imbalances. As a first step, we assume that the magnetic field points in one fixed Cartesian direction but can vary along an orthogonal direction. We start with a uniform-density background threaded by a homogeneous magnetic field and study the evolution of a magnetic perturbation as well as the density fluctuations it induces in the particles. We show that the system evolves through different quasi-equilibrium states and estimate the characteristic time scales on which these quasi-equilibria occur.