We report the results of a series of density functio
nal theory (DFT) calculations aimed at predictingthe
57Fe Mössbauer electric field gradient (EFG) tensors (quadrupole splittings and asymmetry parameters)and their orientations in
S = 0,
1/
2, 1,
3/
2, 2, and
5/
2 metalloproteins and/or model systems. Excellent resultswere found by using a Wachter's all electron basis set for iron, 6-311G* for other heavy atoms, and 6-31G*for hydrogen atoms, BPW91 and B3LYP exchange-correlation functio
nals, and spin-unrestricted methodsfor the paramagnetic systems. For the theory versus experiment correlation, we found
R2 = 0.975, slope= 0.99, intercept = -0.08 mm sec
-1, rmsd = 0.30 mm sec
-1 (
N = 23 points) covering a
EQ range of 5.63mm s
-1 when using the BPW91 functio
nal and
R2 = 0.978, slope = 1.12, intercept = -0.26 mm sec
-1,rmsd = 0.31 mm sec
-1 when using the B3LYP functio
nal.
EQ values in the following systems weresuccessfully predicted: (1) ferric low-spin (
S =
1/
2) systems, including one iron porphyrin with the usual(d
xy)
2(d
xzd
yz)
3 electronic configuration and two iron porphyrins with the more unusual (d
xzd
yz)
4(d
xy)
1 electronicconfiguration; (2) ferrous NO-heme model compounds (
S =
1/
2); (3) ferrous intermediate spin (
S = 1)tetraphenylporphi
nato iron(II); (4) a ferric intermediate spin (
S =
3/
2) iron porphyrin; (5) ferrous high-spin(
S = 2) deoxymyoglobin and deoxyhemoglobin; and (6) ferric high spin (
S =
5/
2) metmyoglobin plus twofive-coordi
nate and one six-coordi
nate iron porphyrins. In addition, seven diamagnetic (
S = 0, d
6 and d
8)systems studied previously were reinvestigated using the same functio
nals and basis set scheme as usedfor the paramagnetic systems. All computed asymmetry parameters were found to be in good agreementwith the available experimental data as were the electric field gradient tensor orientations. In addition, weinvestigated the electronic structures of several systems, including the (d
xy)
2(d
xz,d
yz)
3 and (d
xz,d
yz)
4(d
xy)
1[Fe(III)/porphyri
nate]
+ cations as well as the NO adduct of Fe(II)(octaethylporphi
nate), where interestinginformation on the spin density distributions can be readily obtained from the computed wave functions.