Selecting Quantum-Chemical Methods for Lanthanide-Containing Molecules: A Balance between Accuracy and Efficiency

An analysis of how different density functionals, basis sets, and relativistic approximations affect the computed properties of lanthanide-containing molecules allows one to determine which method provides the highest accuracy. Historically, many different density functional methods have been employed to perform calculations on lanthanide complexes and so herein is a detailed analysis of how different methodological combinations change the computed properties of three different families of lanthanide-bearing species: lanthanide diatomic molecules (fluorides and oxides) and their dissociation energies; larger, molecular complexes and their geometries; and lanthanide bis(2-ethylhexyl)phosphate structures and their separation free energies among the lanthanide series. The B3LYP/Sapporo/Douglas-Kroll-Hess (DKH) method was shown to most accurately reproduce dissociation energies calculated at the CCSDT(Q) level of theory with a mean absolute deviation of 1.3 kcal/mol. For the calculations of larger, molecular complexes, the TPSSh/Sapporo/DKH method led to the smallest deviation from experimentally refined crystal structures. Finally, this same method led to calculated separation factors for lanthanide bis(2-ethylhexyl)phosphate structures that matched very closely with experimental values.