Modeling nuclear volume isotope effects in crystals

Mass-independent isotope fractionations driven by differences in volumes and shapes of nuclei (the field shift effect) are known in several elements and are likely to be found in more. All-electron relativistic electronic structure calculations can predict this effect but at present are computationally intensive and limited to modeling small gas phase molecules and clusters. Density functional theory, using the projector augmented wave method (DFT-PAW), has advantages in greater speed and compatibility with a three-dimensional periodic boundary condition while preserving information about the effects of chemistry on electron densities within nuclei. These electron density variations determine the volume component of the field shift effect. In this study, DFT-PAW calculations are calibrated against all-electron, relativistic Dirac–Hartree–Fock, and coupled-cluster with single, double (triple) excitation methods for estimating nuclear volume isotope effects. DFT-PAW calculations accurately reproduce changes in electron densities within nuclei in typical molecules, when PAW datasets constructed with finite nuclei are used. Nuclear volume contributions to vapor–crystal isotope fractionation are calculated for elemental cadmium and mercury, showing good agreement with experiments. The nuclear-volume component of mercury and cadmium isotope fractionations between atomic vapor and montroydite (HgO), cinnabar (HgS), calomel (Hg ₂ Cl ₂ ), monteponite (CdO), and the CdS polymorphs hawleyite and greenockite are calculated, indicating preferential incorporation of neutron-rich isotopes in more oxidized, ionically bonded phases. Finally, field shift energies are related to Mössbauer isomer shifts, and equilibrium mass-independent fractionations for several tin-bearing crystals are calculated from ¹¹⁹ Sn spectra. Isomer shift data should simplify calculations of mass-independent isotope fractionations in other elements with Mössbauer isotopes, such as platinum and uranium.