In a previous paper, we reported simulations of the evolution of the magnetic field in neutron star cores through ambipolar diffusion, taking the neutrons as a motionless uniform background. However, in real neutron stars, neutrons are free to move, and a strong composition gradient leads to stable stratification (stability against convective motions) both of which might impact on the time-scales of evolution. Here we address these issues by providing the first long-term two-fluid simulations of the evolution of an axially symmetric magnetic field in a neutron star core composed of neutrons, protons, and electrons with density and composition gradients. Again, we find that the magnetic field evolves towards barotropic ``Grad-Shafranov equilibria'', in which the magnetic force is balanced by the degeneracy pressure gradient and gravitational force of the charged particles. However, the evolution is found to be faster than in the case of motionless neutrons, as the movement of charged particles (which are coupled to the magnetic field, but are also limited by the collisional drag forces exerted by neutrons) is less constrained, since neutrons are now allowed to move. The possible impact of non-axisymmetric instabilities on these equilibria, as well as beta decays, proton superconductivity, and neutron superfluidity, are left for future work.

As another step towards understanding the long-term evolution of the magnetic field in neutron stars, we provide the first simulations of ambipolar diffusion in a spherical star. Restricting ourselves to axial symmetry, we consider a charged-particle fluid of protons and electrons carrying the magnetic flux through a motionless, uniform background of neutrons that exerts a collisional drag force on the former. We also ignore the possible impact of beta decays, proton superconductivity, and neutron superfluidity. All initial magnetic field configurations considered are found to evolve on the analytically expected time-scales towards "barotropic equilibria" satisfying the "Grad-Shafranov equation", in which the magnetic force is balanced by the degeneracy pressure gradient, so ambipolar diffusion is choked. These equilibria are so-called "twisted torus" configurations, which include poloidal and toroidal components, the latter restricted to the toroidal volumes in which the poloidal field lines close inside the star. In axial symmetry, they appear to be stable, although they are likely to undergo non-axially symmetric instabilities.

The Adaptive Optics Lucky Imager (AOLI) is a new instrument designed to combine adaptive optics (AO) and lucky imaging (LI) techniques to deliver high spatial resolution in the visible, 20 mas, from ground-based telescopes. Here we present details of the integration and verification phases explaining the defiance that we have faced and the innovative and versatile solution of modular integration for each of its subsystems that we have developed. Modularity seems a clue key for opto-mechanical integration success in the extremely-big telescopes era. We present here the very fresh preliminary results after its first fully-working observing run on the WHT.

AOLI (Adaptive Optics Lucky Imager) is a state-of-art instrument that combines adaptive optics (AO) and lucky imaging (LI) with the objective of obtaining diffraction limited images in visible wavelength at mid- and big-size ground-based telescopes. The key innovation of AOLI is the development and use of the new TP3-WFS (Two Pupil Plane Positions Wavefront Sensor). The TP3-WFS, working in visible band, represents an advance over classical wavefront sensors such as the Shack-Hartmann WFS (SH-WFS) because it can theoretically use fainter natural reference stars, which would ultimately provide better sky coverages to AO instruments using this newer sensor. This paper describes the software, algorithms and procedures that enabled AOLI to become the first astronomical instrument performing real-time adaptive optics corrections in a telescope with this new type of WFS, including the first control-related results at the William Herschel Telescope (WHT).

The Adaptive Optics Lucky Imager, AOLI, is an instrument developed to deliver the highest spatial resolution ever obtained in the visible, 20 mas, from ground-based telescopes. In AOLI a new philosophy of instrumental prototyping has been applied, based on the modularization of the subsystems. This modular concept offers maximum flexibility regarding the instrument, telescope or the addition of future developments.

Upper main sequence stars, white dwarfs and neutron stars are known to possess stable, large-scale magnetic fields. Numerical works have confirmed that stable MHD equilibria can exist in non-barotropic, stably stratified stars. On the other hand, it is unclear whether stable equilibria are possible in barotropic stars, although the existing evidence suggests that they are all unstable. This work aims to construct barotropic equilibria in order to study their properties, as a first step to test their stability. We have assumed that the star is a perfectly conducting, axially symmetric fluid, allowing for both poloidal and toroidal components of the magnetic field. In addition, we made the astrophysically justified assumption that the magnetic force has a negligible influence on the fluid structure, in which case the equilibrium is governed by the Grad-Shafranov equation, involving two arbitrary functions of the poloidal flux. We built a numerical code to solve this equation, allowing for an arbitrary prescription for these functions. Taking particularly simple, but physically reasonable choices for these functions with a couple of adjustable parameters, all of the equilibria found present only a small ($\lesssim 10\%{}$) fraction of the magnetic energy stored in the toroidal component, confirming previous results. We developed an analytical model in order to study in more detail the behavior of the magnetic energy over the full range of parameters. The model confirms that the toroidal fraction of the energy and the ratio of toroidal to poloidal flux are bounded from above for the whole range of parameters.

In stably stratified stars, numerical magneto-hydrodynamics simulations have shown that arbitrary initial magnetic fields evolve into stable equilibrium configurations, usually containing nearly axisymmetric, linked poloidal and toroidal fields that stabilize each other. In this work, we test the hypothesis that stable stratification is a requirement for the existence of such stable equilibria. For this purpose, we follow numerically the evolution of magnetic fields in barotropic (and thus neutrally stable) stars, starting from two different types of initial conditions, namely random disordered magnetic fields, as well as linked poloidal-toroidal configurations resembling the previously found equilibria. With many trials, we always find a decay of the magnetic field over a few Alfv\'en times, never a stable equilibrium. This strongly suggests that there are no stable equilibria in barotropic stars, thus clearly invalidating the assumption of barotropic equations of state often imposed on the search of magnetic equilibria. It also supports the hypothesis that, as dissipative processes erode the stable stratification, they might destabilize previously stable magnetic field configurations, leading to their decay.

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.

We present solutions for Hall equilibria applicable to neutron star crusts. Such magnetic configurations satisfy a Grad-Shafranov-type equation, which is solved analytically and numerically. The solutions presented cover a variety of configurations, from purely poloidal fields connected to an external dipole to poloidal-toroidal fields connected to an external vacuum field, or fully confined within the star. We find that a dipole external field should be supported by a uniformly rotating electron fluid. The energy of the toroidal magnetic field is generally found to be a few percent of the total magnetic field energy for the fields with an external component. We discuss the evolution due to Ohmic dissipation which leads to slowing down of the electron fluid. We also find that the transition from an MHD equilibrium to a state governed by Hall effect, generates spontaneously an additional toroidal field in regions where the electron fraction changes.

Although barotropic matter does not constitute a realistic model for magnetic stars, it would be interesting to confirm a recent conjecture that states that magnetized stars with a barotropic equation of state would be dynamically unstable (Reisenegger 2009). In this work we construct a set of barotropic equilibria, which can eventually be tested using a stability criterion. A general description of the ideal MHD equations governing these equilibria is summarized, allowing for both poloidal and toroidal magnetic field components. A new finite-difference numerical code is developed in order to solve the so-called Grad-Shafranov equation describing the equilibrium of these configurations, and some properties of the equilibria obtained are briefly discussed.

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.