Ab initio local-density- functional- theory calculations of formation energies, surface stress, and multilayer relaxations are reported for the (111), (100), and (110) surfaces of Rh. The study is performed using ultrasoft pseudopotentials and plane waves in a parallel implementation.