The iterative method recently proposed for determining the internal constitution of the outer crust of a nonaccreted neutron star is extended to magnetars by taking into account the Landau-Rabi quantization of electron motion induced by the presence of a very high magnetic field. It is shown that in the strongly quantizing regime, the method can be efficiently implemented using new analytical solutions for the transitions between adjacent crustal layers. Detailed numerical computations are performed to assess the performance and precision of the method.

Magnetic fields of order $10^{15}$ G have been measured at the surface of some neutron stars, and much stronger magnetic fields are expected to be present in the solid region beneath the surface. The effects of the magnetic field on the equation of state and on the composition of the crust due to Landau-Rabi quantization of electron motion are studied. Both the outer and inner crustal regions are described in a unified and consistent way within the nuclear-energy density functional theory.

Magnetars are neutron stars endowed with surface magnetic fields of the order of $10^{14}-10^{15}$~G, and with presumably much stronger fields in their interior. As a result of Landau quantization of electron motion, the neutron-drip transition in the crust of a magnetar is shifted to either higher or lower densities depending on the magnetic field strength. The impact of nuclear uncertainties is explored considering the recent series of Brussels-Montreal microscopic nuclear mass models. All these models are based on the Hartree-Fock-Bogoliubov method with generalized Skyrme functionals. They differ in their predictions for the symmetry energy coefficient at saturation, and for the stiffness of the neutron-matter equation of state. For comparison, we have also considered the very accurate but more phenomenological model of Duflo and Zuker. Although the equilibrium composition of the crust of a magnetar and the onset of neutron emission are found to be model dependent, the quantum oscillations of the threshold density are essentially universal.

The equilibrium properties of the outer crust of cold nonaccreting magnetars (i.e. neutron stars endowed with very strong magnetic fields) are studied using the latest experimental atomic mass data complemented with a microscopic atomic mass model based on the Hartree-Fock-Bogoliubov method. The Landau quantization of electron motion caused by the strong magnetic field is found to have a significant impact on the composition and the equation of state of crustal matter. It is also shown that the outer crust of magnetars could be much more massive than that of ordinary neutron stars.

We have recently developed a set of equations of state based on the nuclear energy density functional theory providing a unified description of the different regions constituting the interior of neutron stars and magnetars. The nuclear functionals, which were constructed from generalized Skyrme effective nucleon-nucleon interactions, yield not only an excellent fit to essentially all experimental atomic mass data but were also constrained to reproduce the neutron-matter equation of state as obtained from realistic many-body calculations.

The equation of state (EoS) of the outer crust of a cold non-accreting magnetar has been determined using the model of Lai and Shapiro (1991). For this purpose, we have made use of the latest experimental atomic mass data complemented with a Hartree-Fock-Bogoliubov (HFB) mass model. Magnetar crusts are found to be significantly different from the crusts of ordinary neutron stars.