The method of extrapolation by intrinsic scaling, recently introduced to obtain correlation energies, is generalized to multiconfigurational reference functions and used to calculate the binding energies of the diatomic molecules C-2, N-2, O-2, and F-2. First, accurate approximations to the full configuration interaction energies of the individual molecules and their constituent atoms are determined, employing Dunning's correlation consistent double-, triple- and quadruple basis sets. Then, these energies are extrapolated to their full basis set limits. Chemical accuracy is attained for the binding energies of all molecules. (c) 2005 American Institute of Physics.