Dataset Information

On Predicting Mossbauer Parameters of Iron-Containing Molecules with Density-Functional Theory.

ABSTRACT: The performance of six frequently used density functional theory (DFT) methods (RPBE, OLYP, TPSS, B3LYP, B3LYP*, and TPSSh) in the prediction of Mössbauer isomer shifts(?) and quadrupole splittings (?E_Q) is studied for an extended and diverse set of Fe complexes. In addition to the influence of the applied density functional and the type of the basis set, the effect of the environment of the molecule, approximated with the conducting-like screening solvation model (COSMO) on the computed Mössbauer parameters, is also investigated. For the isomer shifts the COSMO-B3LYP method is found to provide accurate ? values for all 66 investigated complexes, with a mean absolute error (MAE) of 0.05 mm s^-1 and a maximum deviation of 0.12 mm s^-1. Obtaining accurate ?E_Q values presents a bigger challenge; however, with the selection of an appropriate DFT method, a reasonable agreement can be achieved between experiment and theory. Identifying the various chemical classes of compounds that need different treatment allowed us to construct a recipe for ?E_Q calculations; the application of this approach yields a MAE of 0.12 mm s^-1 (7% error) and a maximum deviation of 0.55 mm s^-1 (17% error). This accuracy should be sufficient for most chemical problems that concern Fe complexes. Furthermore, the reliability of the DFT approach is verified by extending the investigation to chemically relevant case studies which include geometric isomerism, phase transitions induced by variations of the electronic structure (e.g., spin crossover and inversion of the orbital ground state), and the description of electronically degenerate triplet and quintet states. Finally, the immense and often unexploited potential of utilizing the sign of the ?E_Q in characterizing distortions or in identifying the appropriate electronic state at the assignment of the spectral lines is also shown.

SUBMITTER: Papai M

PROVIDER: S-EPMC4358633 | biostudies-literature | 2013 Nov

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

On Predicting Mössbauer Parameters of Iron-Containing Molecules with Density-Functional Theory.

Pápai Mátyás M Vankó György G

Journal of chemical theory and computation 20131015 11

The performance of six frequently used density functional theory (DFT) methods (RPBE, OLYP, TPSS, B3LYP, B3LYP*, and TPSSh) in the prediction of Mössbauer isomer shifts(δ) and quadrupole splittings (ΔEQ) is studied for an extended and diverse set of Fe complexes. In addition to the influence of the applied density functional and the type of the basis set, the effect of the environment of the molecule, approximated with the conducting-like screening solvation model (COSMO ...[more]

PMID: 25821417

Similar Datasets

Project description:This research investigated the thermodynamic favorability and resulting structures for chemical adsorption of trichloroethylene (TCE) to metallic iron using periodic density functional theory (DFT). Three initial TCE positions having the plane defined by HCC atoms parallel to the iron surface resulted in formation of three different chemisorption complexes between carbon atoms in TCE and the iron surface. The Cl-bridge initial configuration with the HCC plane of TCE perpendicular to the iron surface did not result in C-Fe bond formation. The most energetically favorable complex formed at the C-bridge site where the initial configuration had the C=C bond in TCE at a bridge site between adjacent iron atoms. In the C-bridge complex, one C atom formed two a bonds to different Fe atoms, while the second C atom formed a sigma bond with a second Fe atom. Surface complexation atthe C-bridge site resulted in scission of all three C-Cl bonds and also resulted in a shortening of the C==C bond to a distance intermediate between a double and a triple bond. Initial configurations with the C==C bond adsorbed at top or hollow sites on the iron surface resulted in formation of C-Fe a bonds between a single C and two adjacent Fe atoms, and the scission of only two C==Cl bonds. Bond angles and bond lengths indicated that there were no changes in bond order of the C==C bond for top and hollow adsorption. Chemisorption at the C-bridge site had an activation energy of 49 kJ/mol and an early transition state where all three C-CI bonds were activated. The early transition state and the loss of all three Cl atoms upon chemisorption are consistent with most experimental observations that TCE undergoes complete dechlorination in one interaction with the iron surface. The absence of chemisorption and scission of only two C--Cl bonds at the Cl-bridge site is consistent with experimental observations that trace amounts of chloroacetylene may also be produced from reactions of TCE with iron.

Project description:High-spin Fe(IV)-oxo species are known to be kinetically competent oxidants in non-heme iron enzymes. The properties of these oxidants are not as well understood as the corresponding intermediate-spin oxidants of heme complexes. The present work gives a detailed characterization of the structurally similar complexes [Fe(IV)H(3)buea(O)](-), [Fe(III)H(3)buea(O)](2-), and [Fe(III)H(3)buea(OH)](-) (H(3)buea = tris[(N'-tert-butylureaylato)-N-ethylene]aminato) using Mössbauer and dual-frequency/dual-mode electron paramagnetic resonance (EPR) spectroscopies. The [Fe(IV)H(3)buea(O)](-) complex has a high-spin (S = 2) configuration imposed by the C(3)-symmetric ligand. The EPR spectra of the [Fe(IV)H(3)buea(O)](-) complex presented here represent the first documented examples of an EPR signal from an Fe(IV)-oxo complex, demonstrating the ability to detect and quantify Fe(IV) species with EPR spectroscopy. Quantitative simulations allowed the determination of the zero-field parameter, D = +4.7 cm(-1), and the species concentration. Density functional theory (DFT) calculations of the zero-field parameter were found to be in agreement with the experimental value and indicated that the major contribution to the D value is from spin-orbit coupling of the ground state with an excited S = 1 electronic configuration at 1.2 eV. (17)O isotope enrichment experiments allowed the determination of the hyperfine constants ((17)O)A(z) = 10 MHz for [Fe(IV)H(3)buea(O)](-) and ((17)O)A(y) = 8 MHz, ((17)O)A(z) = 12 MHz for [Fe(III)H(3)buea(OH)](-). The isotropic hyperfine constant (((17)O)A(iso) = -16.8 MHz) was derived from the experimental value to allow a quantitative determination of the spin polarization (ρ(p) = 0.56) of the oxo p orbitals of the Fe-oxo bond in [Fe(IV)H(3)buea(O)](-). This is the first experimental determination for non-heme complexes and indicates significant covalency in the Fe-oxo bond. High-field Mössbauer spectroscopy gave an (57)Fe A(dip) tensor of (+5.6, +5.3, -10.9) MHz and A(iso) = -25.9 MHz for the [Fe(IV)H(3)buea(O)](-) complex, and the results of DFT calculations were in agreement with the nuclear parameters of the complex.

Project description:We tested a battery of density functional theory (DFT) methods ranging from generalized gradient approximation (GGA) via meta-GGA to hybrid meta-GGA schemes as well as Møller-Plesset perturbation theory of the second order and a single and double excitation coupled-cluster (CCSD) theory for their ability to reproduce accurate gas-phase structures of di- and triatomic molecules derived from microwave spectroscopy. We obtained the most accurate molecular structures using the hybrid and hybrid meta-GGA approximations with B3PW91, APF, TPSSh, mPW1PW91, PBE0, mPW1PBE, B972, and B98 functionals, resulting in lowest errors. We recommend using these methods to predict accurate three-dimensional structures of inorganic molecules when intramolecular dispersion interactions play an insignificant role. The structures that the CCSD method predicts are of similar quality although at considerably larger computational cost. The structures that GGA and meta-GGA schemes predict are less accurate with the largest absolute errors detected with BLYP and M11-L, suggesting that these methods should not be used if accurate three-dimensional molecular structures are required. Because of numerical problems related to the integration of the exchange-correlation part of the functional and large scattering of errors, most of the Minnesota models tested, particularly MN12-L, M11, M06-L, SOGGA11, and VSXC, are also not recommended for geometry optimization. When maintaining a low computational budget is essential, the nonseparable gradient functional N12 might work within an acceptable range of error. As expected, the DFT-D3 dispersion correction had a negligible effect on the internuclear distances when combined with the functionals tested on nonweakly bonded di- and triatomic inorganic molecules. By contrast, the dispersion correction for the APF-D functional has been found to shorten the bonds significantly, up to 0.064 Å (AgI), in Ag halides, BaO, BaS, BaF, BaCl, Cu halides, and Li and Na halides and hydrides. These results do not agree well with very accurate structures derived from microwave spectroscopy; we therefore believe that the dispersion correction in the APF-D method should be reconsidered. Finally, we found that inaccurate structures can easily lead to errors of few kcal/mol in single-point energies.

Dataset Information

On Predicting Mossbauer Parameters of Iron-Containing Molecules with Density-Functional Theory.

Publications

On Predicting Mössbauer Parameters of Iron-Containing Molecules with Density-Functional Theory.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets