Adsorption of divalent metal ions (Pb2+,Cd2+, Zn2+,Cu2+) on graphene oxide (GO) was studied using density functional theory (DFT). Adsorption geometries and energies, as well as the nature of the binding energy, were calculated for the interaction of divalent metal ions with oxygen-containing groups on the surface of GO. The configurations of the complexes were modeled by placing the divalent metal ions above the center and perpendicular to the surface. Binding of Cu2+ to the GO sheet was predicted to be much stronger than that for other divalent metal ions. Calculated results show good agreement with experimental observations and provide useful information for environmental pollution cleanup.