  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.
