Methods of calculation of potential energy curves or surfaces, including dissociation energies, bond distances, and vibration frequencies, are discussed as well as recently obtained results for several molecules. The ab initio relativistic methods involve the derivation of “shape‐consistent” effective potentials from Dirac–Fock atomic calculations. These effective potentials are averaged and differenced with respect to spin with the differences, p3/2 – p1/2, etc., yielding spin‐orbit operators. The molecular calculations are then set up in a familiar manner through the SCF stage using spin‐averaged effective potentials. The final stage is a configuration‐interaction calculation including the spin‐orbit terms as well as the electron repulsion terms. Calculations that have been made for several low‐lying excited states as well as the ground state for Au2, TlH, Tl2, Sn2, and Pb2 are reviewed. Good agreement is obtained with spectroscopic data and a number of interesting predictions are made.