The (111), (110), and (001) surfaces properties of PuO 2 are studied by using density-functional theory+U method. The total-energy static calculations determine the relative order of stability for low-index PuO 2 surfaces, namely, O-terminated (111) > (110) > defective (001) > polar (001).The effect of thickness is shown to modestly modulate the surface stability and chemical activity of the (110) (111) is found to be the most stable surface. Whereas under O-reducing conditions, the on-surface O-vacancy of C V = 1/9 is stable, and for high reducing conditions, the (111) surface with nearly one monolayer subsurface oxygen removed (C V = 8/9) becomes most stable.