Molecular dynamics simulations of reduced CeO2: bulk and surfaces