The bare and hydrogen-covered diamond (100) surfaces were investigated through pseudopotential density-functional calculations within the local-density approximation. Di erent coverages, ranging from one to two, were considered. These corresponded to di erent structures including 1 1, 2 1, and 3 1, and di erent hydrogen-carbon arrangements including monohydride, dihydride, and con gurations in between. Assuming the system was in equilibrium with a hydrogen reservoir, the formation energy of each phase was expressed as a function of hydrogen chemical potential. As the chemical potential increased, the stable phase successively changed from bare 2 1 to (2 1):H, to (3 1):1.33H, and nally to the canted (1 1):2H. Setting the chemical potential at the energy per hydrogen in H2 and in a free atom gave the (3 1):1.33H and the canted (1 1):2H phase as the most stable one, respectively. However, after comparing with the formation energy of CH4, only the (2 1):H and (3 1):1.33H phases were stable against spontaneous formation of CH4. The former existed over a chemical potential range ten times larger than the latter, which may explain why the latter, despite of having a low energy, has not been observed so far. Finally, the vibrational energies of the C H stretch mode were calculated for the (2 1):H phase.