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 DiracHartree-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 2 Cl 2 ), 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 massindependent fractionations for several tin-bearing crystals are calculated from 119 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.Mössbauer spectroscopy | nuclear field shift | mass independent fractionation A mong the various causes of mass-independent isotope fractionation discovered over the past 40 years, the nuclear field shift effect (1) is unique in that it is an equilibrium thermodynamic phenomenon (and may also impart an identifiable and distinct signature in disequilibrium conditions) (2, 3). It has also proven amenable to prediction by electronic structure methods, at least in simple molecular systems (e.g., refs. 4-7). The goal of the present study is to extend the scope of systems that can be modeled to include condensed phases and crystals, and to speed up calculations for complex materials while retaining enough accuracy to be useful.Nuclear field shift effects result from the finite volume and sometimes nonspherical shapes of atomic nuclei, which are slightly different from one isotope to another. These differences in size and shape affect the coulomb potential felt by bound electrons, and thus the electronic structures of atoms and molecules. Nuclear field shifts have lo...