Ab Initio Calculation of Electric-Field-Gradient Tensors Of
Total Page:16
File Type:pdf, Size:1020Kb
American Mineralogist, Volume 81, pages 545-549, 1996 Ab initio calculation of electric-field-gradient tensorsof forsterite B.rdnN WrNxlrnrr Pnrnn BllHn,2.q,No KARLHTTNzScuwl.nz2 'Mineralogisch-PetrographischesInstitut der Christian-Albrechts Universitiit, Olshausenstrasse40, D 24098 Kiel, Germany 'zInstitutfiir TechnischeElektrochemie. Technische Universitiit Wien. A-1060 Vienna, Austna Ansrn-lcr Ab initio band-structure calculations basedon the density functional theory have been performed for forsterite to obtain, with a parameter-freemodel, electric-field gradientsfor all nuclei. Calculations basedon the generalizedgradient approximation yield the ratio of the largest components of the calculated electric-field gradients within l0loof the experi- mental value for 25Mgand within 2oloof the experimental value for '7O. The absolute values differ by about 5ol0,depending on which nuclear quadrupole moment is used in the conversion. The asymmetry parameters are also in good agreement with experimental data. Values obtained with a gradient-correctedexchange-correlation potential are better than those based on the standard local-density approximation. The calculated anglesbe- tween the principal axes of the quadrupole coupling tensor and the crystallographicaxes agreevery well with the experimental data. For the Ml site the maximum deviation is 1.8",and for the M2 site the maximum deviation is 0.6". The current calculationsallow an evaluation of three sets of experimental data for therTO electric-field gradient. They also confirm a proposed assignmentof the measuredelectric-field gradient tensors to spe- cific O atoms. INrnooucrroN V,", which varies between0 (axial symmetry) and 1. Here we follow the sameconventions as in the work of Schwarz In computational mineralogy, predictive microscopic et al. (1990). models are increasingly used to interpret experimental However, QCC and 4 cannot be derived from NMR data for which conflicting phenomenological explana- spectradirectly but must be determined by spectral sim- tions exist. Although reliable quantum mechanicalmodel ulation, which can be problematic for broad, overlapping, calculations require large computational resources,they or weak features. Conventionally, EFGs have been ex- offer the possibility of calculating rather subtle physical plained by point-chargemodels, which contain adjustable properties, including those due to electronic effects. parameters, such as formal chargesor the Sternheimer Therefore, they can be applied to a large variety ofprob- antishielding factors. The predictive power of such par- lems and yield better results than conventional static lat- ameterizedcalculations is limited, and the results are of- tice-energy-minimization calculations. Here we use pa- ten ambiguous. Hence, attempts have been made to de- rameter-free energy-band-structurecalculations to velop parameter-freecalculations. For example, Tossell compute electric-field gradient (EFG) tensors,which are and coworkers [e.g., Tossell (1992) and referencescited experimentally determined by nuclear magnetic reso- therein] usedHartree-Fock calculationson clustersto ob- nance (NMR), nuclear quadrupole resonance,or, M<iss- tain NMR parameters. These calculations required the bauer spectroscopy. use of several approximations owing to limitations in NMR has proven to be an invaluable tool for the in- available computing power. With a different method, vestigation of local structural details, and, because of which could copewith periodic lattices,Blaha et al. (1985) progressin the experimental technique, many nuclei rel- showed that full-potential, Iinearized, augmented plane- evant for the earth sciencescan now be used.An instruc- wave (FP-LAPW) calculations can be successfullyused tive review of applications of NMR spectroscopyin the to obtain EFGs. The crucial quantity that determines an earth scienceswas given by Kirkpatrick (1988). Structural EFG is the nonspherical chargedistribution close to the information can be obtained from NMR measurements nucleus,which itself is directly related to chemical bond- in which the interaction betweenthe electric quadrupole ing and different occupations of the valence electrons. moment, eQ, and the local electric-field gradient, eq : During the last decade,this method was used to inves- V,", is determined in terms of the quadrupole-coupling tigate severalhcp metals (Blaha et al. 1988) and many constant, QCC, givenby QCC : eV""Q/h. The EFG de- other substances,and its feasibility for complex struc- pendson the nonsphericalcharge distribution around the tures, for example, was demonstrated by successfulcal- nuclei, and its deviation from axial symmetry is defined culations for a high-2" cuprate superconductor(Schwarz by the asymmetryparameter, 4, given by ,t: (4" - V-)/ et al. 1990). From these calculations it can be deduced, 0003-004x/96l0506-0545$05.00 545 546 WINKLER ET AL.: ELECTRIC.FIELD GRADIENTS OF FORSTERITE TABLE1. Atomic coordinates used in the calculations and the proximation (LDA) or the generalizedgradient approxi- corresponding remaining forces mation (GGA). A recent in-depth summary can be found Fujino in the book by Singh (1994).For such calculations,one Hazen Forces Brown Forces et al. Forces Forces of the most accurate schemesis the LAPW method, in (1e76) (GGA) (1e70) (GGA) (1981) (LDA) (GGA) which the unit cell is divided into spherescentered at the Mg2x 0.9915 0.6 0.9896 1.0 0.99169 0.0 0.4 atomic positions and an interstitial region. In the latter -0.1 -0.4 -0.3 -0.4 Mg2y 0.2774 0.2776 0.27739 the basis set consistsof plane waves that are augmented Si x 0.4262 6.1 0 4226 48 6 0.42645 -1.3 2.2 Si y 0.0940 -0.9 0.0945 -4.5 0.09403 -1.9 1.4 by atomic-like solutions (numerical radial functions mul- 01 x 0.7657 18.6 0.7667 -3.4 0.76594 17.0 17.4 tiplied with sphericalharmonics) inside the spheres.The -0.4 01 y 0.0913 2.2 0.0918 1.7 0.09156 1.0 potential within the sphereis not restrictedto have spher- 02 x 0.2215 6.6 0.2202 11 3 0.22164 41 4.7 02 y 0.4474 -21.6 0 4477 -25.9 0.44705 -5.4 -14.8 ical symmetry, as in the older "muffin-tin" calculations 03 x 0.2777 -7.4 0.2781 -13 5 0.27751 -4.5 -5.9 but is allowedto be general,i.e., a full potential without y 03 0.1628 11 .9 0.1633 14.4 0.16310 2.5 8.5 It is well known that 03 z 0.0331 -15.6 0.0337 -23.9 0.03304 -6.0 -12.9 any shape approximation is used. DFT calculations in general give structural parameters given : Nofer Forces in mRyd/au;1 Ryd 13.6 eV, 1 au: 0.529A. within 1-20loof the experimental data. In the presentstudy, the full-potential, linearized, aug- mentedplane-wave package, WIEN95 (Blahaet al. 1995) that in generalFP-LAPW calculations yield EFGs within was used.[WIEN95 is an improved and updatedversion approximately 100/oand asymmetry parameters within ofthe original copyrightedcode by Blaha et al. (1990).1 -0. I of the experimentalvalues. WIEN95 is maintained at the Technical University of Consideringthe widespreaduse of NMR spectroscopy Vienna and includes the capability to employ local or- and the recent progressin quantum mechanical calcula- bitals (Singh 1991) and to calculateresidual forces on tions it is worthwhile to compute electric-field-gradient atoms (Yu et al. l99l: Kohler et al. 1996).The calcula- tensorsfor silicates,which are characterizedby their broad tion of the EFGs is based on the work of Blaha et al. range of bond types and bond strengths. We note that (1985); a more detailed description can be found in NMR spectroscopyis a local technique, and therefore it Schwarzet al. (1990).The precisionand accuracyofthe is often used to investigateshort-range order or local dis- calculationsare controlled by only a few parameters.For tortions. Ab initio calculations, however, exploit the pe- the expansion ofthe chargedensity (potential) in the in- riodicity of a perfectly ordered structure, so that solid terstitialregion a Fourier serieswith lG | = t Z wasused. solutions and the effectsofshort-range order are currently For the wave functions inside the atomic spheresangular approximated rather crudely. Also, the present calcula- momentum componentsup to /: 12 wereincluded, and tions are restricted to ground-state(0 K) properties. up to 18 k points in the irreducible wedgeof the Brillouin For fully ordered structures,theoretical studiesare pos- zonewere used. This correspondsto nearly 150 sampling sible and can be used for the interpretation of measured points in the whole Brillouin zone. The cut-off for the spectra,which is sometimes difficult using experimental plane-wavebasis set was chosenas k-u* : 4.5 ait and results only. Reliable calculations allow evaluations of was checkedby additional calculations,which proved that proposedexperiments, which canbe costlywhen isotopic the results presentedhere are well converged. enrichment is required, becausethey can predict NMR parametersand specify the necessaryresolution. For the Rnsur,rs interpretation of spectraof silicates,correlations between structural and NMR parameters have been established Structure (Sherritretal. l99l; Smith et al. 1983).These studies rely The all-electron ab initio techniques currently in use on available data, which may be scarcefor unusual iso- do not yet allow a constant-pressurerelaxation of low- topes. In theoretical studies, however, structures can be symmetry structures.In fact, even at constant volume a easily deformed, and thus NMR parametersas a function relaxation of an orthorhombic structure containing 28 of structural parameters can be calculated. Also, theory atoms using an all-electron approach would at present can help in evaluating sets of experimental data if the consumea prohibitive amount of computer time. Hence, latter are not in good agreement. in the present study all structural parameterswere fixed Forsterite,MgrSiOo, is a model systemin which the at their experimental values. To investigatethe influence same isotopes occupy different structural environments of small structural distortions on the calculated physical 2sMg and in which the EFGs for both and "O have been propertieswe usedthree setsof internal structural param- studied in detail. The small differencesin the Mg and O eters in our calculations.