Potential energy surfaces (PES) for rotation about the N-C(sp(3)) or N-C(aryl) bond and energies of stationary points on PES for rotation about the C(sp(2))-N bond are reported for methylurea, ethylurea, isopropylurea, tert-butylurea, and phenylurea, using the B3LYP/DZVP2 and MP2/aug-cc-pVDZ methods. The analysis of alkylureas reveals cis and (less stable) trans isomers that adopt anti geometries, whereas syn geometries do not correspond to stationary points. In contrast, the analysis of phenylurea reveals that the lowest energy form at the MP2 level is a trans isomer in a syn geometry. The fully optimized geometries are in good agreement with crystal structure data, and PESs are consistent with the experimental dihedral angle distribution. Rotation about the C(Sp(2))-N bond in alkylureas and phenylurea is slightly more hindered (8.6-9.4 kcal/ mol) than the analogous motion in the unsubstituted molecule (8.2 kcal/mol). At the MP2 level of theory, the maximum barriers to rotation for the methyl, ethyl, isopropyl, tert-butyl, and phenyl groups are predicted to be 0.9, 6.2, 6.0, 4.6, and 2.4 kcal/mol, respectively. The results are used to benchmark the performance of the MMFF94 force field. Systematic discrepancies between MMFF94 and MP2 results were improved by modification of several torsional parameters.