A shell-core model is applied to investigate the effects of Zr doping on the surface energies and surface structures of the three low Miller index surfaces in UO2 using the molecular dynamics (MD) technique. The surface energies and atomic structures of the Zr-doped and undoped UO2 (100), (110) and (111) surfaces are compared. Simulation results indicate that (i) the surface energy of (U1-y,Zry)O-2 depend on the crystallographic orientation, as well as of undoped UO2. The (100) surface exhibits the highest surface energy, followed by the (110) surface, and the (111) surface; (ii) Zr doping will significantly increase the surface energy of UO2 by approximately 20% on (100) surface, 10% on (110) surface and 15% on (111) surface with the ZrO2 contents ranging from 0 to 12.5 mol%, respectively; (iii) the surface energies of the three low Miller index surfaces decrease with increasing temperature both in undoped UO2 and in (U1-y,Zry)O-2; (iv) the addition of Zr induces a severe distortion of the (U1-y,Zry)O-2 surface structure, and the outermost top layer exhibits the strongest rumpling; (v) the considerable reconstructions can be observed in the two top layers of Zr-doped and undoped UO2 surfaces when the temperature is elevated to 900-1500K. (C) 2015 Elsevier B.V. All rights reserved.