Calculation of Chemical-Shift Tensors of Heavy Nuclei: A DFT/ZORA Investigation of 199 Hg Chemical-Shift Tensors in Solids, and the Effects of Cluster Size and Electronic-State Approximations

Calculations of the nuclear magnetic resonance chemical-shielding tensors of a suite of mercury-containing materials using various cluster models for the structures provide a stringent test of the procedures for forming models and for calculation with various methods. The inclusion of higher co-ordination shells in the molecular clusters permits quantum chemical calculations of 199 Hg chemical-shielding tensor elements within 3% of the experimental values. We show that it is possible to reduce the size of computationally expensive molecular-cluster calculations with limited effect on calculated NMR parameters by carefully introducing the frozen core approximation. The importance of the relativistic Hamiltonian for accurate predictions of chemical-shielding values is demonstrated within the molecular cluster approach. The results demonstrate that careful design of a cluster to represent the solid-state structure, inclusion of relativistic components in the Hamiltonian at least at the spin-orbit level, and judicious use of approximations are essential to obtain good agreement with experimental results.


Introduction
Nuclear magnetic resonance (NMR) chemical shieldings and quadrupolar couplings provide a means of examining the electronic state of a material. The full tensor interaction, available in the solid, gives a description of the three-dimensional electronic structure. 1,2 The connection between the electronic state and these parameters is embodied in complex equations that rely on accurate models of the electronic wave function. In general, very good approximate computational predictions of isotropic NMR parameters for light nuclei such as 13 C and 1 H in molecules have become routine nowadays, and such calculations of the parameters are often a part of the assignment of spectra. 3,4 Prediction of NMR parameters of heavier nuclei such as 199 Hg, 205 Tl, and 207 Pb is usually limited to prediction of properties of isolated structures in the gas phase [5][6][7][8] and in the solution phase. [8][9][10][11] The experimentally determined ranges of chemical shieldings (chemical shifts) for heavy nuclei like 199 Hg and 207 Pb are generally much larger than the ranges for lighter nuclei such as 13 C and 1 H. For example, the range of experimental principal components for 199 Hg in the various materials studied is almost 4500 ppm. The range of experimental principal components for 207 Pb is almost 11 700 ppm. 12 Calculations of NMR parameters of heavy nuclei in a solid have two obstacles not generally encountered in ''standard'' quantum chemical calculations on molecules in solution or the gas phase. First, some, if not all, of the electrons must be treated as relativistic particles, in part because of the high nuclear charge. 13,14 Second, the local intermolecular environment with its periodic structure has geometric effects that are not averaged by rapid motion. 15 To address the latter, two methods have been proposed to calculate NMR parameters of nuclei in solids: (1) use of periodic boundary conditions based on the full crystal symmetry, and (2) treatment of model molecular clusters that define the local environment.
In the last ten years, periodic-boundary calculations have become popular for estimating NMR parameters in solids of known structure. An important method which employs periodic boundary conditions was developed by Sebastiani et al. 16,17 with the use of localized Wannier orbitals as implemented in the program CPMD. 18 Another popular formalism is the Gauge Included Projected Augmented Wave (GIPAW) method of Pickard and Mauri. 19 The general theory and applications to solid-state NMR of the GIPAW method have been reviewed recently. 20 The GIPAW method is currently not very suitable for calculation of the NMR parameters of heavy nuclei for the following two reasons. First, such calculations are computationally very expensive for unit cells with large numbers of electrons, which is usually the case for systems containing heavy nuclei like 207 Pb and 199 Hg. Second, relativistic effects, which can be very important for heavy nuclei, are treated in a limited manner (or not at all) in codes that employ periodic-boundary conditions.
Calculations using model clusters that simulate the local solid-state environment are an alternative to the periodicboundary approach. In early studies, Tossell et al. investigated the effects of next-nearest neighbors on NMR chemical shielding for 15 N and 23 Na by employing molecular clusters to mimic the solid-state structure. [21][22][23] Similarly, Valerio et al. employed molecular-cluster-based calculations for NMR chemical shielding of 29 Si and 27 Al in zeolite structures. 24,25 In these cases, it was shown that more accurate chemical-shift and quadrupolar parameters result when the cluster size is sufficient to reflect the nature of local structure on the electronic state. More recently, 19 F chemical shifts have been computed for several different fluorine sites with the cluster approach. 26 Another attempt at using cluster models to mimic the solid-state structure involves embedding a cluster in a point-charge field. [27][28][29] The cluster-based approach has been employed in some cases to obtain solid-state NMR parameters of heavy nuclei such as 207 Pb and 199 Hg. [30][31][32][33][34][35][36][37][38][39] In general, NMR calculations on clusters have centered on the effects of nearest neighbors, principally due to the large number of electrons in these systems. The agreement with experimental data is usually qualitative.
There are important issues in employing molecular clusters to simulate the solid-state environment. For example, there is frequently a large number of parameters (charge of the cluster, dangling bonds, etc.) which must be optimized to give reasonable results. One important challenge in a cluster-based approach is to define the co-ordination sphere of the NMR nucleus and the appropriate size of the cluster. Some differences between cluster and periodic-boundary calculations for determining NMR parameters have been summarized in a recent article. 29 In this contribution, we report calculations of NMR chemicalshielding tensors for a series of solid mercury compounds, in which we use the cluster approach. We systematically evaluate the effect of molecular-cluster size on NMR parameters of a suite of 199 Hg-containing materials. We examine the effect of application of the frozen core approximation (FCA) in calculating NMR chemical shielding using molecular clusters, and we examine the effect of neglect of various parts of the relativistic Hamiltonian in NMR calculations by treating the molecular clusters at different levels of theory. Our aim is to provide examples of how various approximations affect the calculation of NMR parameters in a variety of solids.

Calculational protocols
The formal definition of the chemical-shielding interaction is given by the components of a rank-2 three-dimensional Cartesian tensor, r: In eqn (1), E represents the total energy of the system, l is the magnetic moment of the nucleus of interest, and B ext is the external magnetic field. The chemical-shielding tensor can be represented as a 3 Â 3 matrix in a Cartesian co-ordinate system. (2) Although the chemical-shielding tensor may have nine independent elements in a given co-ordinate system (as shown in eqn (2)), 40 it is often considered that the antisymmetric components are small or zero, which reduces the number of independent elements to six. In one specific co-ordinate system, the principal-axis system (PAS), the chemical-shielding tensor is diagonal. These principal components are conventionally assigned as: s 11 r s 22 r s 33 , which is known as the frequency-ordered convention. The chemical-shielding tensor in other co-ordinate systems may be derived from this tensor in the PAS by rotational transformation between the two co-ordinate systems.
The absolute chemical shielding is not measured in a NMR experiment. Instead, the chemical shift, the shielding relative to some reference material, is generally reported. The principal components of the chemical shielding and the observable principal chemical-shift components are related by eqn (3).
An approximation of eqn (3)  There are several other descriptions in use. In one convention, the tensor is expressed in terms of its irreducible spherical-tensor components. 40 Another convenient description for spectra of powders is given by the Maryland convention, a set of three NMR parameters describing the shape of the spectrum of a random powder, the isotropic chemical shift (d iso ), the span (O), and the skew (k). 43 These three parameters are related to the principal components of the chemical-shift tensor by the following relations: All computations were performed using the ADF (ADF v2010 and ADF v2013) suite of programs. [44][45][46] Calculations were carried out at the DFT/BP86 level. 47,48 For the results in Sections 3.1 and 3.2, the TZ2P/all-electron (AE) basis set is used for all atoms in the cluster. Relativistic effects are treated with the ZORA Hamiltonian at the spin-orbit level. [49][50][51][52] For the calculations in Section 3.3, the TZ2P/AE basis set is employed for atoms in the near vicinity of the NMR nucleus of interest, including the nearest neighbors, and the frozen core approximation 53,54 (FCA) is used with TZPand DZ-type basis sets for the remainder of the cluster. For the calculations in Section 3.4, we examine how inclusion of various relativistic terms affects the NMR parameters by comparison to results that neglect these effects. In all cases, NMR parameters are computed within the GIAO formalism. 11,55,56 An important issue in calculations of NMR parameters of molecular clusters is the size of the cluster and how it affects predicted parameters. All clusters examined are formed from experimental X-ray geometries found in the literature ( Table 1). The set includes materials with relatively small first-coordination shells, and does not contain examples of large mercury-containing complexes. For the solid systems with hydrogen atoms [Hg(acetate) 2 and Hg 2 (NO 3 ) 2 Á2H 2 O], an optimization was carried out on hydrogen positions of small clusters at the ZORA/scalar level with the BP86 density functional and TZ2P/AE basis set. For HgF 2 and HgI 2 , the clusters are terminated with protons to reduce the charge on the cluster and obtain SCF convergence.
To specify the agreement or lack of agreement between calculated shifts and experimental observations, we report the residual, the square root of the sum of the squares of the deviations of the calculated principal components from the experimental values. The smaller the residual, the better a model agrees with experiment.

The effect of size and charge
Although NMR chemical shielding is known to be primarily a local property, the electronic state from which one derives chemical shielding is influenced strongly by long-range interactions, especially in solids. To check for the importance of this long-range effect, we evaluate chemical shielding for two types of molecular clusters (called, respectively, small and large clusters). Small clusters consist of only the atoms with NMRactive nuclei ( 199 Hg) and the nearest neighbors (ligands). Large clusters include the extended co-ordination shell, next-nearest neighbors, and the rest of the molecule of which the next-nearest neighbors and extended co-ordination shell atoms are a part. Pictorial representations of the small and large clusters of HgCl 2 are given in Fig. 1.
We investigated the nature of cluster size on the NMR chemical shielding of Hg 2 Cl 2 described by various clusters formed from the X-ray crystal structure of the compound (Fig. 2). All calculations were performed within the ZORA/SO framework while employing a TZ2P/AE basis set. For convenience, the predicted NMR chemical-shift principal components were referenced to a single value, the calculated chemical shielding of the isolated dimethylmercury (DMM) molecule at the same level of theory, 7965 ppm. The calculated NMR parameters are tabulated in Table 2, along with the experimental results. For cluster I, which consists of the isolated Hg 2 Cl 2 unit, the calculated isotropic shift is about 580 ppm more negative than the experimental value, whereas the span is in error by almost 2000 ppm. This error in the span mostly arises from an error in the d 33 component of the chemical-shielding tensor, whereas the other two components are within 100 ppm of the respective experimental values.
Cluster II includes the extended co-ordination shell of the NMR nucleus of interest. The added atoms (compared to cluster I)  The small cluster consists of NMR-active nuclei (pink cross-hatched) and the first co-ordination shell which are nearest neighbors to NMR-active nuclei (green and gray beach ball). A large cluster consists of the central molecule with its NMR-active nuclei and the first co-ordination shell, as well as the extended co-ordination shell (orange and gray), nearest neighbors of the first co-ordination shell (blue striped), and atoms to complete the molecules (red). are the four Cl À ions. The distance from the mercury nucleus to these additional chlorine centers is 0.321 nm, compared to the Hg-Cl bond distance of 0.243 nm. The addition of these four chlorine centers significantly improves the calculated value of d 33 relative to the neglect of these centers in cluster I. The introduction of these centers, however, causes significant errors in the values of d 11 and d 22 that were not seen for cluster I. The calculated span is much smaller than the experimental span. Cluster II does not have certain symmetry elements that cluster I does have, in particular the mirror plane between the two Hg centers. It appears that this lack of symmetry is one possible cause of the difference between the values of these two components for cluster I and cluster II. Cluster III adds elements to include the mirror symmetry that was lost in the definition of cluster II. Compared to cluster II, cluster III has four additional chlorine centers present. Calculations on cluster III give an isotropic shift that is 800 ppm more negative than the experimental value, whereas the span is underestimated by about 700 ppm. In general, all the calculated principal components of cluster III are more shielded than the experimental values. This observation is often the fingerprint of excess charge on the molecular cluster, which is the case for both cluster II and cluster III, having charges of À4 and À8, respectively.
Cluster IV is designed to compensate charge by termination with hydrogen atoms. These hydrogen atoms were inserted along the mercury-chlorine bond axis at 0.127 nm from the chlorine. This position maintains the symmetry at the mercury site whose NMR parameters we calculate. The compensation of charge, while maintaining the symmetry, decreases the isotropic shielding of the Hg center compared to that of cluster III. Additionally, residuals of the principal components are much  smaller than those of cluster III. Saturation of the dangling bonds to decrease charge on the cluster is important in using clusters to model the local structure in solids for calculating NMR parameters, as these calculations show. However, in such treatments, the calculated chemical-shielding tensor is usually dependent on the positions of the hydrogen atoms.
In cluster V, we include only the Hg 2 Cl 2 units which contain the four next-nearest Hg atoms to the NMR nucleus of interest. The Hg-Hg distances between the Hg 2 Cl 2 units are 0.448 nm. As expected, inclusion of next-nearest Hg centers has a smaller effect on NMR parameters than inclusion of the next-nearest Cl centers, since the Hg-Hg distance is larger than the Hg-Cl distance. Moreover, since the next-nearest chlorines are missing in this cluster, the agreement between experiment and calculation is not good. Compared to cluster I, there is a considerable difference between the principal components, which indicates that next-nearest Hg atoms should be included in the extended co-ordination shell in designing the cluster.
In cluster VI, not only is the first co-ordination shell present, but the extended co-ordination shell and next-nearest neighbors are also included. Additionally, the structure is completed by adding atoms to produce uncharged Hg 2 Cl 2 units, as shown in Fig. 1. Inclusion of these additional centers improves the calculated values of the chemical-shift tensor components, so that the residual is the smallest of any of the clusters examined. We have found that, with available computational resources to create clusters, a residual of 200 ppm or less for these kinds of heavy atoms indicates a reasonable model cluster. With this accuracy and the known wide ranges of shifts for 199 Hg, calculation with these cluster models allows assignment of shifts to nuclei in unique structural motifs.

Calculations on clusters in the all-electron basis
We report NMR parameters for a series of mercury compounds, using large-cluster models and small-cluster models (Fig. 3). The calculated NMR chemical-shift parameters are tabulated in Table 3, along with reported experimental values. (The calculated chemical shieldings from which these shifts are derived are in Table S1 in the ESI. †) s ref for determining chemical shifts is taken as the intersection of the best-fit linear correlation for each type of cluster ( Fig. 4a and b). These values are 7853 ppm and 8043 ppm for small and large clusters, respectively.
In general, calculations on small clusters predict isotropic chemical shifts (d iso ) to within 300-1000 ppm of the experimental values. As seen in Table 3, the deviations of the principal components of the chemical-shift tensor calculated for small clusters are often more severely in error (relative to experimental data) than the isotropic shifts. The differences between computed spans (which do not depend on errors in referencing) and experimental spans often exceed 1000 ppm. For this reason, conclusions about the quality of models based on the near agreement of a calculated isotropic chemical shift with an experimental chemical shift can be misleading. More importantly, the principal components of a chemical-shift tensor reflect important properties of the electronic structure that cannot be discerned by concentrating on the isotropic chemical shift alone. For example, the isotropic chemical shift of HgCl 2 calculated with a small cluster is only 27 ppm from the experimental value, but the calculated span deviates by 1457 ppm from the experimental value. The spans calculated for large clusters are closer to the experimental values for all cases we have examined. Additionally, there is better agreement between experimental and predicted isotropic chemical shifts, as well as better agreement between calculated and experimental principal components of the shift tensor, as can be discerned from the smaller residuals for the large clusters. With the exception of the hydrated mercurous nitrate, the residuals of all the large-cluster-model calculations are near or below 200 ppm, which seems to be some limit for clusters of this size for these systems. But, even for the hydrated mercurous nitrate, the large-cluster model gives better agreement with experiment than the small-cluster model.
The performance of small-and large-cluster models is shown graphically in Fig. 4a and b. In Fig. 4a, the correlation between calculated chemical-shielding tensor components and experimental chemical-shift tensor components is shown for small clusters. In Fig. 4b, the same correlation is shown for large clusters. The correlation for the large clusters has a slope that deviates by only 2.3% from the ideal slope of À1, whereas the linear correlation for the small clusters deviates by 36% from this ideal value.
From the intersection of the best-fit linear correlation lines in Fig. 4a and b, the absolute shielding of DMM can be estimated within the two-component ZORA/SO framework for large and small clusters. Within this framework, the absolute shielding of DMM is predicted to be 7853 ppm by the small clusters and 8043 ppm by the large clusters. In separate calculations on an isolated DMM molecule at the same level of theory and with either (a) a fully optimized geometry or (b) the experimental geometry 63 with optimized hydrogen positions, the predicted isotropic shielding in case (a) was 8120 ppm, and in case (b) it was 7965 ppm. In a previous study, Taylor et al. 33 showed (with the same level of theory) that the absolute shielding of DMM was between 7929 and 8095 ppm.
The shielding constant of DMM has also been calculated by means of four-component relativistic methods. 64 64 showed that 199 Hg shielding constants calculated with ZORA/SO and with four-component DFT follow a similar trend, and the chemical shifts calculated with these two methods are in agreement within 60 ppm. Autschbach shows that the valence-shell properties such as chemical shift and J coupling are well described in the ZORA formalism. 67,68 This observation is supported by our findings that calculations using ZORA produce results in agreement with experiment.
In Ramsey's formulation, 69 the shielding is evaluated as an integral of operators which vary as 1/r 3 , where r is the distance from the electron to the nucleus of interest. One expects that major contributions to NMR shielding are predominantly from orbitals that place the electron near the nucleus. For solids, the positions of nearby nuclei are exceedingly important in determining the general structure of orbitals. Comparison of the results for small and large clusters shows that medium-to-long-range effects must be taken into account to determine meaningful NMR tensor parameters.
There are two factors that contribute to the difference between the results for small clusters and for large clusters. First, missing atoms in the near region in the small cluster for a 199 Hg-containing system may affect the chemical shielding.
A secondary effect results from the fact that the nearest neighbors of the ligands are not included in the small clusters. As a result, the molecular orbitals (MOs) forming Hg-L bonds are largely localized, which would not adequately represent the structural effects on NMR chemical shielding. Both effects depend on the crystal structure, and they contribute differently to the 199 Hg shielding in the two models.

The use of the frozen core approximation
The frozen core approximation (FCA) is commonly used in many applications to trim the computational time. It is generally thought that the deep core electrons are not strongly influenced by changes that may influence the valence electrons. We have investigated whether use of the FCA for these clusters yields reliable NMR parameters by the following procedure. For the large clusters of Fig. 3, the electrons of the first co-ordination shell near the nucleus of interest were treated with an all-electron (AE) basis set, whereas the rest of the cluster was treated with the FCA. TZP and DZ basis sets with a frozen large core are used in the part of the cluster treated with the FCA. (The details of the frozen core basis sets are given in Table S2 of the ESI. †) The results of these calculations are given in Table 4 for various mercury-containing solids. (The corresponding chemical-shielding  values are given in Table S3, ESI. †) Fig. 5a and b show the correlation of the calculated FCA/AE principal components with the experimental chemical-shift principal components. From the correlation, s ref is found to be 8056 ppm or 8087 ppm, respectively, for FCA(TZP)/AE and FCA(DZ)/AE basis sets. The results determined with FCA(TZP)/AE are closer to the all-electron results, as expected. The difference between principal components calculated with FCA(TZP)/AE and the AE basis sets does not exceed 25 ppm for most cases. The largest differences are 113 ppm for d 33 of HgI 2 and 69 ppm for d iso (= d 11 = d 22 = d 33 for this case) of HgF 2 . The differences between principal components determined with FCA(DZ)/AE compared to the all-electron basis set are slightly larger. Nevertheless, given the accuracy of the NMR calculations on heavy nuclei, calculations using the FCA for the remote atoms and calculations with the full all-electron basis set produce almost the same NMR parameters. In fact, Fig. 5a and b show that the correlation for results using the FCA is slightly closer to the ideal case than for results using the full all-electron basis set for the large clusters. This difference is within the uncertainty of the experimental measurements. The two calculations essentially give the same results within experimental error.
Although there are very small differences in the computed NMR parameters, there is a substantial difference in computational time that favors using the FCA. In Table 5, we compare the computational time and number of Cartesian functions employed for selected calculations, when using the all-electron basis set versus the FCA(TZP)/AE and FCA(DZ)/AE basis sets. The calculations with the triple-zeta basis set generally take longer than those with the double-zeta basis set, but they are both substantially shorter than a calculation on the same cluster using the all-electron wave function. For example, the all-electron calculation for a large cluster of Hg 2 Cl 2 took nearly 2.5 days, whereas the FCA(TZP)/AE treatment of the same cluster required 3.6 hours and the FCA(DZ)/AE treatment took only 35 minutes. These results indicate that the FCA can be applied to the peripheral portions of large clusters in solid-state NMR calculations with minimal introduction of significant errors.
An important criterion in the cluster-based approaches is the convergence of NMR parameters with increasing cluster size. Since the computational requirements are relatively low at the FCA level of theory, it is possible to investigate whether calculated parameters are converged with these larger clusters. For this reason, we calculated the shielding values for Hg(SCN) 2 and Hg 2 Cl 2 at the FCA(TZP)/AE level of theory, and with molecular clusters which are extended another co-ordination shell from the 'large clusters' (which we call extended clusters). (Fig. 6) Table 6 shows the comparison of chemical shielding components for large clusters and extended clusters of Hg(SCN) 2 À990 À1014 À4047 À2017 3032 489 FCA(DZ)/AE Hg(2) À820 À1084 À4048 À1984 3201 461 FCA(TZP)/AE Hg(1) À993 À1069 À4070 À2044 3053 516 FCA(TZP)/AE Hg(2) À820 À1090 À4066 À1992 3220 468 All-electron Hg(1) À988 À1078 À4058 À2042 3045 515 All-electron Hg (2) À813 À1092 À4056 À1987 3217 464 and Hg 2 Cl 2 . For Hg(SCN) 2 , the large cluster and the extended cluster produce similar results. The differences between s 33 and s 11 are 22 ppm and 9 ppm, whereas s 22 is the same for the large cluster and the extended cluster. The difference in s iso is only 4 ppm. For Hg 2 Cl 2 , the differences are slightly larger than those of Hg(SCN) 2 . The largest difference is for the s 33 . (79 ppm) On the other hand, the differences for s iso are 5 ppm and 13 ppm for Hg (1) and Hg (2) respectively. For all calculated principal components, the variations between the two clusters are less than 1%, which would be a valid convergence criterion for heavy-nucleus shieldings.

Effect of level of relativistic theory on NMR parameters
An important consideration in the calculation of chemical shielding of heavy nuclei such as 199 Hg, 207 Pb, and 119 Sn is the inclusion of relativistic effects. 7 Including relativistic terms in the Hamiltonian increases the computational time, which makes it of some practical importance to determine how significantly these affect the calculated chemical shielding. We investigated these effects by comparing results of calculations that (a) neglect all relativistic terms (nonrelativistic results) with (b) results of calculations that include only scalar relativistic corrections, and with (c) results of calculations that include the full spin-orbit interaction for large clusters. Fig. 7a shows the correlation between calculated chemical shieldings and experimental chemical shifts for all tensor components of the large clusters, when the system is treated with only nonrelativistic terms. The linear correlation has a slope of À0.6274 and an intercept of 6232 ppm. Fig. 7b shows a similar correlation for the same systems when only scalar relativistic corrections are considered. The slope of the linear correlation is À0.6697 for this system and the intercept is 5976 ppm. The inclusion of the scalar relativistic terms improves agreement with the experimental data marginally over the nonrelativistic results.
The data of Fig. 5b show the result of inclusion of relativistic effects at the spin-orbit level. The slope of the linear correlation   line is À1.0234 and the intercept is 8043 ppm. Thus, within less than 3%, the slope of this linear trend is correct. In addition, comparison of Fig. 7 and 4b shows that scatter about the ideal line is much smaller when the ZORA spin-orbit Hamiltonian is used. Significantly, there is a systematic error of approximately 2000 ppm in the predicted value of the chemical shielding of DMM when the spin-orbit terms are neglected. From these calculations, it is clear that one must include spin-orbit relativistic terms in the Hamiltonian when calculating chemical shieldings of mercury-containing solids.

Conclusions
We have used the cluster-model approach to predict 199 Hg NMR parameters of eleven Hg-containing solids. The comparison of results on small clusters (often a single isolated unit) with large clusters demonstrates that to predict the NMR chemical shielding accurately, one must take into account the effects of the surrounding environment, and specifically one must treat the problem with a large cluster that maintains symmetry at the site of the nucleus of interest. We have found that, in some cases, the isotropic shift may be fortuitously predicted for an isolated unit, but that the tensor components may be in substantial error. The cluster model must be chosen so that the local site symmetry at the nucleus of interest is maintained, which determines the minimal size. Care must be taken to ensure that the charge on a cluster (which is sometimes necessary) is not too large. The situation for large organomercury complexes may involve the higher coordination shells to a lesser degree than for the small molecules investigated here, in which case a single molecule model may be adequate. Of concern in calculation of chemical shielding in these heavynucleus-containing materials is the computation time that a large-cluster model may require. We examined the use of the frozen core approximation (FCA) for atomic centers removed from the locality of interest as a means to shorten this calculation time. The calculations show that using the FCA on more remote centers makes a minimal difference in computed NMR parameters, compared to calculations with models using the full all-electron basis set. The use of the FCA for these more remote centers presents a substantial savings in computational time when compared with the all-electron approach, and we suggest that using the FCA in this manner for computation of NMR parameters may allow one to specify the environment with larger clusters that better define the effects of structure, without sacrificing much accuracy. The impact of relativistic terms in the Hamiltonian is significant for these heavy nuclei. In particular, inclusion of spin-orbit terms is essential to get chemical shieldings that reflect the experimental data. Total neglect of relativistic terms is the poorest approximation, inclusion of scalar relativistic terms improves the calculation slightly, but for the highest accuracy, one must include spin-orbit components in the Hamiltonian.