Small-core Effective Core Potentials Research Articles

Improving the theoretical description of Ln(III)/An(III) separation with phosphinic acid ligands: a benchmarking study of structure and selectivity.

The efficient separation of trivalent lanthanides from minor actinides with soft-donor ligands, while showing experimental promise, has theorists continuing to search for suitable approaches for describing and interpreting selectivity. To remedy this, we employ several computational methods in describing the structure of model M(H2PX2)3 complexes, with M = Eu and Am, and X = O, S, Se, and Te, and predicting the selectivity of model phosphinic acid ligands in Eu(III)/Am(III) separation. After first establishing a set of MP2 and CCSD(T)-DKH3 results as benchmarks, we evaluate several density functionals and descriptions of valence shells for their accuracy with respect to metal-ligand bonding and selectivity. We find that commonly employed functionals with a 0-27% range of exact exchange used with small-core effective core potentials or with an explicit treatment of the relativistic effects (DKH2) incorrectly predict a decrease in the metal-ligand bond distance in going from Eu(III) to Am(III) and completely fail to track a selectivity trend, even giving a wrong sign for some or all ligands. Surprisingly, when these functionals are used in conjunction with an f-in-core description of metal ions, the correct trend in selectivity is recovered, though Am-X distances are overestimated in relation to Eu-X. Functionals with high components of exact exchange (50%) and double-hybrid functionals are reasonably aligned with benchmark results, pointing to the problems of DFT with small exact exchange fractions to handle f-electrons. Natural bond orbital analyses reveal that these poorly performing functionals disproportionately overpopulate outer f orbitals in the model complexes. We anticipate that recommendations resulting from this work will lead to more accurate theoretical descriptions of lanthanide/actinide selectivity with soft-donor chalcogen-based ligands in the future.

Polarizabilities, dispersion coefficients, and retardation functions at the complete basis set CCSD limit: From Be to Ba plus Yb.

Static and dynamic polarizabilities of alkaline earth metal atoms Be-Ba and of the Yb atom, as well as dispersion coefficients and retardation functions for their long-range interactions, are used as a benchmark for the restricted coupled cluster method with singles and doubles (CCSD) and noniterative triples added [CCSD(T)] and related polarization propagator CCSD(3) methods at the complete basis set limit. The latter is attained through the sequence of the augmented correlation-consistent polarized weighted core valence n-zeta basis sets with the exact 2-component approximation for the scalar relativistic effects and with the small-core effective core potentials (for Ca, Sr, and Ba). At the converged level of core correlation treatment, the finite-field CCSD(T) method reproduces the best available data for the static dipole and quadrupole polarizabilities better than 1% and 4%, respectively. Systematic cancelation of the triple contribution in the CCSD(3) calculations of the dynamic polarizabilities of alkaline earth metal atoms makes their dispersion coefficients accurate within 3%. The retardation functions are computed and used for the analysis of the long-range interactions in the homonuclear dimers. Implications to accurate ab initio calculations of the global interaction potentials are discussed.

The Relativistic Effects on the Carbon-Carbon Coupling Constants Mediated by a Heavy Atom.

The (2)JCC, (3)JCC, and (4)JCC spin-spin coupling constants in the systems with a heavy atom (Cd, In, Sn, Sb, Te, Hg, Tl, Pb, Bi, and Po) in the coupling path have been calculated by means of density functional theory. The main goal was to estimate the relativistic effects on spin-spin coupling constants and to explore the factors which may influence them, including the nature of the heavy atom and carbon hybridization. The methods applied range, in order of reduced complexity, from the Dirac-Kohn-Sham (DKS) method (density functional theory with four-component Dirac-Coulomb Hamiltonian), through DFT with two- and one-component zeroth-order regular approximation (ZORA) Hamiltonians, to scalar effective core potentials (ECPs) with the nonrelativistic Hamiltonian. The use of DKS and ZORA methods leads to very similar results, and small-core ECPs of the MDF and MWB variety reproduce correctly the scalar relativistic effects. Scalar relativistic effects usually are larger than the spin-orbit coupling effects. The latter tend to influence the most the coupling constants of the sp(3)-hybridized carbon atoms and in compounds of the p-block heavy atoms. Large spin-orbit coupling contributions for the Po compounds are probably connected with the inverse of the lowest triplet excitation energy.

Europium-doped silicon clusters EuSi n (n = 3–11) and their anions: structures, thermochemistry, electron affinities, and magnetic moments

The structures, electron affinities, and dissociation energies of EuSi n (n = 3–11) and their anions have been examined by means of four hybrid and pure density functional theory (DFT) methods. Basis sets used in this work are of segmented (SEG) Gaussian valence basis sets and relativistic small-core effective core potentials (ECP) with additional diffuse 2pdfg functions, denoted aug-SEG/ECP for Eu atoms and aug-cc-pVTZ for Si atoms. The geometries are fully optimized with each DFT method independently. The ground-state structures for all of these species are found to be substitutional type, which can be regarded as being derived from the ground-state structure of Si n+1 (and/or Si +1 − ) by replacing a Si atom with a Eu atom. The theoretical adiabatic electron affinities (AEAs) of EuSi n predicted by the four DFT schemes are in excellent agreement with the experimental data, especially the AEAs of TPSSh and B2PLYP. The average absolute deviations from experiment are by 0.10, 0.06, 0.07, and 0.05 eV, and the largest deviations are 0.16, 0.12, 0.18, and 0.10 eV at the B3LYP, TPSSh, PBE, and B2PLYP levels, respectively. The AEA of EuSi n (n = 3–11) is less than that of Si n . With the increase in silicon cluster size, the AEA of EuSi n may be close to that of Si n , but cannot be larger than that of Si n . The Eu atom acts as an electron donor, and the bonding between Eu and silicon clusters is ionic in nature. The bond between Eu and silicon clusters of neutral EuSi n (n = 3–11) is stronger than that of the anions. The total magnetic moments of EuSi n /EuSi (n = 3–11) and the magnetic moments on the Eu atom do not quench, and the total magnetic moments are contributed by Eu atom. The dissociation energies of Eu atom from EuSi n and their anions have also been calculated to examine relative stabilities.

The Predicted Spectrum and Singlet–Triplet Interaction of the Hypermetallic Molecule SrOSr

In accordance with previous studies in our group on Be, Mg, and Ca hypermetallic oxides, we find that SrOSr has a linear X̃1Σ(g)+ ground electronic state and a very low lying first excited ã3Σ(u)+ triplet electronic state. No gas-phase spectrum of this molecule has been assigned yet, and to encourage and assist in its discovery we present a complete ab initio simulation, with absolute intensities, of the infrared absorption spectrum for both electronic states. The three-dimensional potential energy surfaces and the electric dipole moment surfaces of the X̃1Σ(g)+ and ã3Σ(u)+ electronic states are calculated using a multireference configuration interaction (MRCISD) approach in combination with internally contracted multireference perturbation theory (RS2C) based on complete active space self-consistent field (CASSCF) wave functions applying a Sadlej pVTZ basis set for both O and Sr and the Stuttgart relativistic small-core effective core potential for Sr. The infrared spectra are simulated using the MORBID program system. We also calculate vertical excitation energies and transition moments for several excited singlet and triplet electronic states in order to predict the positions and intensities of the most prominent singlet and triplet electronic absorption bands. Finally, for this heavy molecule, we calculate the singlet–triplet interaction matrix elements between close-lying vibronic levels of the X̃ and ã electronic states and find them to be very small.

Theoretical investigation of the He–I2(E3Πg) ion-pair state: Ab initio intermolecular potential and vibrational levels

We present a theoretical study on the potential energy surface and vibrational bound states of the E electronic excited state of the HeI(2) van der Waals system. The interaction energies are computed using accurate ab initio methods and large basis sets. Relativistic small-core effective core potentials in conjunction with a quintuple-zeta quality basis set are employed for the heavy iodine atoms in multireference configuration interaction calculations for the (3)A' and (3)A" states. For the representation of the potential energy surface we used a general interpolation technique for constructing potential surfaces from ab initio data based on the reproducing kernel Hilbert space method. The surface presents global and local minima for T-shaped configurations with well-depths of 33.2 and 4.6 cm(-1), respectively. Vibrational energies and states are computed through variational quantum mechanical calculations. We found that the binding energy of the HeI(2)(E) T-shaped isomer is 16.85 cm(-1), in excellent agreement with recent experimental measurements. In lieu of more experimental data we also report our predictions on higher vibrational levels and we analyze the influence of the underlying surface on them. This is the first attempt to represent the potential surface of such a highly excited electronic state of a van der Waals complex, and it demonstrates the capability of the ab initio technology to provide accurate results for carrying out reliable studies to model experimental data.

Stannylenes: Structures, Electron Affinities, Ionization Energies, and Singlet–Triplet Gaps of SnX2/SnXY and XSnR/SnR2/RSnR′ Species (X; Y = H, F, Cl, Br, I, and R; R′ = CH3, SiH3, GeH3, SnH3)

Systematic computational studies of stannylene derivatives SnX(2)/SnXY and XSnR/SnR(2)/RSnR' were carried out using density functional theory. The basis sets used for H, F, Cl, Br, C, Si, and Ge atoms are of double-ζ plus polarization quality with additional s- and p-type diffuse functions, denoted DZP++. For the iodine and tin atoms, the Stuttgart-Dresden basis sets, with relativistic small-core effective core potentials (ECP), are used. All geometries are fully optimized with three functionals (BHLYP, BLYP, and B3LYP). Harmonic vibrational wavenumber analyses are performed to evaluate zero-point energy corrections and to determine the nature of the stationary points located. Predicted are four types of neutral-anion separations, plus adiabatic ionization energies (E(IE)) and singlet-triplet energy gaps (ΔE(S-T)). The dependence of all three energetic properties upon choice of substituent is remarkably strong. The EA(ad(ZPVE)) values (eV) obtained with the B3LYP functional range from 0.70 eV [Sn(CH(3))(2)] to 2.36 eV [SnI(2)]. The computed E(IE) values lie between 7.33 eV [Sn(SnH(3))(2)] and 11.15 eV [SnF(2)], while the singlet-triplet splittings range from 0.60 eV [Sn(SnH(3))(2)] to 3.40 eV [SnF(2)]. The geometries and energetics compare satisfactorily with the few available experiments, while most of these species are investigated for the first time. Some unusual structures are encountered for the SnXI(+) (X = F, Cl, and Br) cations. The structural parameters and energetics are discussed and compared with the carbene, silylene, and germylene analogues.

The enthalpies of formation of AsXn molecules, where X = H, F or Cl, and n = 1, 2 or 3, by RCCSD(T) and UCCSD(T)-F12x calculations

RCCSD(T) and UCCSD(T)-F12x calculations were performed on AsX(n) molecules, where X = H, F or Cl, and n = 1, 2 or 3, and related species, in order to evaluate their enthalpies of formation (ΔH(f)(Ø)). The recommended ΔH(f)(Ø) values obtained from the present investigation are AsH, 57.7(2); AsF, -7.9(3); AsCl, 27.2(4); AsH(2), 39.8(4); AsF(2), -96.6(9); AsCl(2), -17.8(10); AsH(3), 17.1(4); AsF(3)-196.0(5) and AsCl(3), -59.1(27) kcal mole(-1). These values are anchored only on one thermodynamic quantity, namely, ΔH(f)(Ø)(As) (= 70.3 kcal mole(-1)). In the calculations, the fully-relativistic small-core effective core potential (ECP10MDF) was used for As. Contributions from outer core correlation of As 3d(10) electrons were computed explicitly in both RCCSD(T) and UCCSD(T)-F12 calculations with additional tight basis functions designed for As 3d(10) electrons. Basis sets of up to augmented correlation-consistent polarized valence quintuple-zeta (aug-cc-pV5Z) quality were used in RCCSD(T) calculations and computed relative electronic energies were extrapolated to the complete basis set (CBS) limit. For the simplified, explicitly correlated UCCSD(T)-F12x calculations, basis sets of up to quadruple-zeta (QZ) quality were employed. Based on the RCCSD(T)/CBS benchmark values, the reliability of available theoretical and experimental values have been assessed.

Ab initio calculations on the X̃ B21 and à A21 states of AsH2, and Franck–Condon simulation, including anharmonicity, of the Ã(,,)-X̃ single vibronic level emission spectrum of AsH2

Restricted-spin coupled-cluster single-double plus perturbative triple excitation {RCCSD(T)} calculations were carried out on the X (2)B(1) and A (2)A(1) states of AsH(2) employing the fully relativistic small-core effective core potential (ECP10MDF) for As and basis sets of up to the augmented correlation-consistent polarized valence quintuple-zeta (aug-cc-pV5Z) quality. Minimum-energy geometrical parameters and relative electronic energies were evaluated, including contributions from extrapolation to the complete basis set limit and from outer core correlation of the As 3d(10) electrons employing additional tight 4d3f2g2h functions designed for As. In addition, simplified, explicitly correlated CCSD(T)-F12 calculations were also performed employing different atomic orbital basis sets of up to aug-cc-pVQZ quality, and associated complementary auxiliary and density-fitting basis sets. The best theoretical estimate of the relative electronic energy of the A (2)A(1) state of AsH(2) relative to the X (2)B(1) state including zero-point energy correction (T(0)) is 19,954(32) cm(-1), which agrees very well with available experimental T(0) values of 19,909.4531(18) and 19,909.4910(17) cm(-1) obtained from recent laser induced fluorescence and cavity ringdown absorption spectroscopic studies. In addition, potential energy functions (PEFs) of the X (2)B(1) and A (2)A(1) states of AsH(2) were computed at different RCCSD(T) and CCSD(T)-F12 levels. These PEFs were used in variational calculations of anharmonic vibrational wave functions, which were then utilized to calculate Franck-Condon factors (FCFs) between these two states, using a method which includes allowance for anharmonicity and Duschinsky rotation. The A(0,0,0)-X single vibronic level (SVL) emission spectrum of AsH(2) was simulated using these computed FCFs. Comparison between simulated and available experimental vibrationally resolved spectra of the A(0,0,0)-X SVL emission of AsH(2), which consist essentially of the bending (2(n)) series, suggests that there is a significant loss in intensity in the low emission energy region of the experimental spectrum.

Franck–Condon simulation of the photoelectron spectrum of AsF2 and the photodetachment spectrum of AsF2− using ab initio calculations: Ionization energy and electron affinity of AsF2

RCCSD(T) and/or CASSCF/MRCI calculations were carried out on the X(2)B(1) state of AsF(2), the X(1)A(1), ã(3)B(1) and A(1)B(1) states of AsF(2)(+), and the X(1)A(1) state of AsF(2)(-) employing the fully-relativistic small-core effective core potential (ECP10MDF) for As and basis sets of up to augmented correlation-consistent polarized valence quintuple-zeta (aug-cc-pV5Z) quality. Minimum-energy geometrical parameters and relative electronic energies were evaluated, including contributions from extrapolation to the complete basis set limit and from outer core correlation of the As 3d(10) electrons. In addition, simplified, explicitly correlated RHF/UCCSD(T)-F12x calculations were also performed employing different atomic orbital basis sets, and associated complementary auxiliary and density-fitting basis sets. The best theoretical estimates of the adiabatic ionization energies (AIE(0)) of AsF(2)(X(2)B(1)) to the X(1)A(1) and ã(3)B(1) states of AsF(2)(+), including corrections for zero-point vibrational energy (DeltaZPE), are 9.099(8) and 13.290(22) eV respectively. The best estimated electron affinity (EA(0)) of AsF(2) is 1.182(16) eV, also including Delta(ZPE). These are currently the most reliable AIE(0) and EA(0) values for AsF(2). Potential energy functions (PEFs) of the X(2)B(1) state of AsF(2), the X(1)A(1) and ã(3)B(1) states of AsF(2)(+) and X(1)A(1) state of AsF(2)(-) were computed at RCCSD(T)/aug-cc-pV5Z, RCCSD(T)/aug-cc-pCV5Z and RHF/UCCSD(T)-F12a/aug-cc-pCVTZ levels. These PEFs were employed in variational calculations of anharmonic vibrational wavefunctions, which were then utilised to calculate Franck-Condon factors (FCFs), using a method which includes allowance for anharmonicity and Duschinsky rotation. The computed FCFs were used to simulate the first two bands in the photoelectron spectrum of AsF(2) and the first band in the photodetachment spectrum of AsF(2)(-), both yet to be recorded. The simulated spectra obtained using different sets of PEFs were found to be almost identical, suggesting that the simplified explicitly correlated UCCSD(T)-F12x method with a relatively small basis set can be a reliable alternative to the conventional RCCSD(T) correlation methods with a relatively large basis set, but at a significantly lower cost.

Theoretical Actinide Molecular Science

Interest in the chemistry of the early actinide elements (notably uranium through americium) usually results either from the nuclear waste problem or the unique chemistry of these elements that result from 5f contributions to bonding. Computational actinide chemistry provides one useful tool for studying these processes. Theoretical actinide chemistry is challenging because three principal axes of approximation have to be optimized. These are the model chemistry (the choice of approximate electron-electron correlation method and basis sets), the approximate relativistic method, and a method for modeling solvent (condensed phase) effects. In this Account, we arrange these approximations in a three-dimensional diagram, implying that they are relatively independent of each other. A fourth level of approximation concerns the choice of suitable models for situations too complex to treat in their entirety. We discuss test cases for each of these approximations. Gas-phase data for uranium fluorides and oxofluorides such as UF(6) and UO(2)F(2) show that GGA functionals provide accurate geometries and frequencies while hybrid density functional theory (DFT) functionals are superior for energetics. MP2 is seen to be somewhat erratic for this set of compounds, and CCSD(T) gives the most accurate results. Three different relativistic methods, small-core effective core potentials (SC-ECP), ZORA, and all-electron scalar, provide comparable results. The older large-core ECP (LC-ECP) approach is consistently worse and should not be used. We confirmed these conclusions through studies of the actinyl aquo complexes [AnO(2)(OH(2))(5)](n+), (An = U, Np, or Pu and n = 1 or 2) that are also used to test solvation models. As long as the first coordination sphere of the metal is included explicitly, continuum solvation models are reliable, and we found no clear advantage for the (costly) explicit treatment of the second coordination sphere. Spin-orbit effects must be included to reproduce the correct trend in An(VI)/An(V) reduction potentials. We propose a multistep mechanism for the experimentally observed oxygen exchange of UO(2)(2+) cations in highly alkaline solutions present in tank wastes. This process involves an equilibrium between [UO(2)(OH)(4)](2-) and [UO(2)(OH)(5)](3-), followed by formation of the stable [UO(3)(OH)(3)](3-) intermediate that forms from [UO(2)(OH)(5)](3-) through intramolecular water elimination. The [UO(3)(OH)(3)](3-) intermediate facilitates oxygen exchange through proton shuttling. We explain the experimentally observed stabilization of the pentavalent oxidation state of actinyl ions by macrocyclic ligands (such as 18-crown-6) as an effect of solvation: the large macrocycle screens the positive charge of the ion from the polarizable solvent. Alkyl-substituted isoamethyrin complexes are bent despite being aromatic because of steric factors, rather than fit/misfit criteria regarding the actinyl ion. By application of an efficient DFT code, actinide molecules with more than 100 atoms can now be studied routinely. "Real" chemical questions can be answered as long as we take great care to apply methods that are accurate with respect to the three axes of approximation described above. While the exclusive focus of this Account has been on the early actinide elements, these conclusions also apply elsewhere in the periodic table.

Crown Ether Inclusion Complexes of the Early Actinide Elements, [AnO2(18-crown-6)]n+, An = U, Np, Pu andn= 1, 2: A Relativistic Density Functional Study

The title compounds, [AnO2(18-crown-6)]n+, An = U, Np, and Pu and n = 1 and 2, as well as the related (experimentally observed) complex [UO2(dicyclohexyl-18-crown-6)]2+ are studied using relativistic density functional theory (DFT). Different relativistic methods (large-core and small-core effective core potentials, all-electron scalar four-component) and two flavors of approximate DFT (B3LYP and PBE) are used. Calculated bond lengths agree well with the available experimental data for the NpV complex, while larger differences for the UVI complexes appear to be related to the large uncertainties in the experimental data. The axial AnO bonds are found to be weaker and longer than in the corresponding penta-aquo complexes, though still of partial triple-bond character. The AnO bond lengths and strengths decrease along the actinide series, consistent with the actinide contraction. Gas-phase binding energies calculated for the penta-aquo complexes and crown-ether complexes of the actinides studied, as well as ligand-exchange energies, show that there is no intrinsic preference, or "better fit", for actinyl(V) cations as compared to actinyl(VI) ones. Rather, the ability of NpO2+ (NpV) to form in-cavity 18-crown-6 complexes in water, which is impossible for UO22+, is traced to solvation effects in polar solvents. Thus, the experimentally observed stabilization of the pentavalent oxidation state as compared to the hexavalent one is due to the effective screening of the charge provided by the macrocycle, and this leads to destabilization of the AnVI crown complexes relative to their AnV counterparts.

Diffuse basis functions for small-core relativistic pseudopotential basis sets and static dipole polarizabilities of selected lanthanides La, Sm, Eu, Tm and Yb

The compact sets of the diffuse functions are suggested for better description of the polarization properties of the lanthanide atoms using the small-core effective core potentials and the corresponding basis sets. The pd and 2pdfg augmentations with exponents 0.028 and 0.015 (p), 0.032 (d) and 0.05 (f,g) are recommended for all lanthanides. Scalar and tensor dipole polarizabilities calculated using the multireference averaged quadratic coupled cluster method agree well with the literature data.

Optimized accurate auxiliary basis sets for RI-MP2 and RI-CC2 calculations for the atoms Rb to Rn

The introduction of the resolution-of-the-identity (RI) approximation for electron repulsion integrals in quantum chemical calculations requires in addition to the orbital basis so-called auxiliary or fitting basis sets. We report here such auxiliary basis sets optimized for second-order Moller–Plesset perturbation theory for the recently published (Weigend and Ahlrichs Phys Chem Chem Phys, 2005, 7, 3297–3305) segmented contracted Gaussian basis sets of split, triple-ζ and quadruple-ζ valence quality for the atoms Rb–Rn (except lanthanides). These basis sets are designed for use in connection with small-core effective core potentials including scalar relativistic corrections. Hereby accurate resolution-of-the-identity calculations with second-order Moller–Plesset perturbation theory (MP2) and related methods can now be performed for molecules containing elements from H to Rn. The error of the RI approximation has been evaluated for a test set of 385 small and medium sized molecules, which represent the common oxidation states of each element, and is compared with the one-electron basis set error, estimated based on highly accurate explicitly correlated MP2–R12 calculations. With the reported auxiliary basis sets the RI error for MP2 correlation energies is typically two orders of magnitude smaller than the one-electron basis set error, independent on the position of the atoms in the periodic table.

Ab initio study of a Bi3+ impurity in Cs2NaYCl6 and Y2O3: Comparison of perturbative and variational electron correlation methods

Ab initio study of excitation energies and oscillator strengths for absorption towards the (3)P(1) and (1)P(1) states of the Bi(3+) ion has been performed for the Bi(3+) ion in gas phase and as a dopant of the cubic elpasolite Cs(2)NaYCl(6) and the yttria Y(2)O(3) crystal using the ab initio embedded-cluster method. The ground and excited states were computed with a relativistic spin-orbit configuration interaction approach suited for heavy elements. Electron correlation was treated in the scalar relativistic scheme with perturbative, variational, and coupled-cluster methods. Intermediate coupling is included via an effective-Hamiltonian based spin-orbit configuration interaction approach. Small-core (60 electrons) and large-core (78 electrons) relativistic effective core potentials (ECPs) have been used to describe the bismuth ion. The best match with experiment was obtained with the small-core ECP. The accuracy of excitation energies strongly depends on the electron correlation method used. The agreement between experimental data and the results obtained using second-order multiconfigurational perturbation theory is greatly improved with the shifted zeroth-order Hamiltonian proposed by Ghido et al. [Chem. Phys. Lett. 396, 142 (2004)]. Although quite time consuming, coupled-cluster and variational methods yield good agreement with experimental data. The first absorption band recorded for the doped elpasolite crystal is positioned with an excellent accuracy while the computed energy of the second absorbing manifold is in poorer agreement with experimental data. This suggests that interactions with neglected close-lying excited states with a ligand-to-metal charge transfer character may be significant. Calculations of the spectrum of Bi(3+) doping yttria in both the S(6) and C(2) site symmetries indicate that the absorbing manifold arises from electronic excitations localized on the Bi(3+) doping ion with main triplet 6s6p character. Our results predict the first absorbing peak to lie about 0.5 eV lower for the S(6) sites than for the C(2) site, thus attributing the violet and the green emission wavelengths to the S(6) and C(2) sites, respectively. A subsequent study of Stokes shift and emission wavelength should hopefully lead to a final assignment of the measured excitation spectra.

Density Functional Studies of Actinyl Aquo Complexes Studied Using Small-Core Effective Core Potentials and a Scalar Four-Component Relativistic Method

ADVERTISEMENT RETURN TO ISSUEPREVAddition/CorrectionORIGINAL ARTICLEThis notice is a correctionDensity Functional Studies of Actinyl Aquo Complexes Studied Using Small-Core Effective Core Potentials and a Scalar Four-Component Relativistic MethodGrigory A. Shamov and Georg SchreckenbachCite this: J. Phys. Chem. A 2006, 110, 43, 12072Publication Date (Web):October 11, 2006Publication History Published online11 October 2006Published inissue 1 November 2006https://doi.org/10.1021/jp0662855Copyright © 2006 American Chemical SocietyRIGHTS & PERMISSIONSArticle Views523Altmetric-Citations43LEARN ABOUT THESE METRICSArticle Views are the COUNTER-compliant sum of full text article downloads since November 2008 (both PDF and HTML) across all institutions and individuals. These metrics are regularly updated to reflect usage leading up to the last few days.Citations are the number of other articles citing this article, calculated by Crossref and updated daily. Find more information about Crossref citation counts.The Altmetric Attention Score is a quantitative measure of the attention that a research article has received online. Clicking on the donut icon will load a page at altmetric.com with additional details about the score and the social media presence for the given article. Find more information on the Altmetric Attention Score and how the score is calculated. Share Add toView InAdd Full Text with ReferenceAdd Description ExportRISCitationCitation and abstractCitation and referencesMore Options Share onFacebookTwitterWechatLinked InReddit PDF (13 KB) Get e-Alerts Get e-Alerts

On the Structure and Spin States of Fe(III)-EDDHA Complexes

DFT methods are suitable for predicting both the geometries and spin states of EDDHA-Fe(III) complexes. Thus, extensive DFT computational studies have shown that the racemic-Fe(III) EDDHA complex is more stable than the meso isomer, regardless of the spin state of the central iron atom. A comparison of the energy values obtained for the complexes under study has also shown that high-spin (S = 5/2) complexes are more stable than low-spin (S = 1/2) ones. These computational results matched the experimental results of the magnetic susceptibility values of both isomers. In both cases, their behavior has been fitted as being due to isolated high-spin Fe(III) in a distorted octahedral environment. The study of the correlation diagram also confirms the high-spin iron in complex 2b. The geometry optimization of these complexes performed with the standard 3-21G* basis set for hydrogen, carbon, oxygen, and nitrogen and the Hay-Wadt small-core effective core potential (ECP) including a double-xi valence basis set for iron, followed by single-point energy refinement with the 6-31G* basis set, is suitable for predicting both the geometries and the spin-states of EDDHA-Fe(III) complexes. The presence of a high-spin iron in Fe(III)-EDDHA complexes could be the key to understanding their lack of reactivity in electron-transfer processes, either chemically or electrochemically induced, and their resistance to photodegradation.

Density Functional Studies of Actinyl Aquo Complexes Studied Using Small-Core Effective Core Potentials and a Scalar Four-Component Relativistic Method

The title compounds, [AnO(2)(H(2)O)(5)](n)(+), n = 1 or 2 and An = U, Np, and Pu, are studied using relativistic density functional theory (DFT). Three rather different relativistic methods are used, small-core effective core potentials (SC-ECP), a scalar four-component all-electron relativistic method, and the zeroeth-order regular approximation. The methods provide similar results for a variety of properties, giving confidence in their accuracy. Spin-orbit and multiplet corrections to the An(VI)/An(V) reduction potential are added in an approximate fashion but are found to be essential. Bulk solvation effects are modeled with continuum solvation models (CPCM, COSMO). These models are tested by comparing explicit (cluster), continuum, and mixed cluster/continuum solvation models as applied to various properties. The continuum solvation models are shown to accurately account for the effects of the solvent, provided that at least the first coordination sphere is included. Reoptimizing the structures in the presence of the bulk solvent is seen to be important for the equatorial bond lengths but less relevant for energetics. Explicit inclusion of waters in the second coordination sphere has a modest influence on the energetics. For the first time, free energies of solvation are calculated for all six [AnO(2)(H(2)O)(5)](n)(+) species. The calculated numbers are within the experimental error margins, and the experimental trend is reproduced correctly. By comparison of different relativistic methods, it is shown that an accurate relativistic description leads to marked improvements over the older large-core ECP (LC-ECP) method for bond lengths, vibrational frequencies, and, in particular, the An(VI)/An(V) reduction potential. Two approximate DFT methods are compared, B3LYP, a hybrid DFT method, and PBE, a generalized gradient approximation. Either method yields An(VI)/An(V) reduction potentials of comparable quality. Overall, the experimental reduction potentials are accurately reproduced by the calculations.