  We present predictions of the binding energy per nucleon and the neutron skin
thickness in highly neutron-rich isotopes of Oxygen, Magnesium, and Aluminum.
The calculations are carried out at and below the neutron drip line as
predicted by our model. The nuclear properties are obtained via an energy
functional whose input is the equation of state of isospin-asymmetric infinite
matter. The latter is based on a microscopic derivation of the energy per
particle in neutron matter applying chiral few-nucleon forces together with a
phenomenological model for the equation of state of symmetric nuclear matter.
We highlight the impact of the neutron matter equation of state at different
orders of chiral effective field theory on neutron skins and the binding energy
per particle and quantify the uncertainty carried by our predictions.
