Efficient and accurate approximations to the molecular spin-orbit coupling operator and their use in molecular g-tensor calculations

Approximations to the Breit-Pauli form of the spin-orbit coupling (SOC) operator are examined. The focus is on approximations that lead to an effective quasi-one-electron operator which leads to efficient property evaluations. In particular, the accurate spin-orbit mean-field (SOMF) method developed by Hess, Marian, Wahlgren, and Gropen is examined in detail. It is compared in detail with the "effective potential" spin-orbit operator commonly used in density functional theory (DFT) and which has been criticized for not including the spin-other orbit (SOO) contribution. Both operators contain identical one-electron and Coulomb terms since the SOO contribution to the Coulomb term vanishes exactly in the SOMF treatment. Since the DFT correlation functional only contributes negligibly to the SOC the only difference between the two operators is in the exchange part. In the SOMF approximation, the SOO part is equal to two times the spin-same orbit contribution. The DFT exchange contribution is of the wrong sign and numerically shown to be in error by a factor of 2-2.5 in magnitude. The simplest possible improvement in the DFT-SOC treatment Veff(-2X)-SOC is to multiply the exchange contribution to the Veff operator by -2. This is verified numerically in calculations of molecular g-tensors and one-electron SOC constants of atoms and ions. Four different ways of handling the computationally critical Coulomb part of the SOMF and Veff operators are discussed and implemented. The resolution of the identity approximation is virtually exact for the SOC with standard auxiliary basis sets which need to be slightly augmented by steep s functions for heavier elements. An almost as efficient seminumerical approximation is equally accurate. The effective nuclear charge model gives results within approximately 10% (on average) of the SOMF treatment. The one-center approximation to the Coulomb and one-electron SOC terms leads to errors on the order of approximately 5%. Small absolute errors are obtained for the one-center approximation to the exchange term which is consequently the method of choice SOMF(1X) for large molecules.