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.
