Semiempirical local spin: Theory and implementation of the ZILSH method for predicting Heisenberg exchange constants of polynuclear transition metal complexes

The local spin formalism ( 3 ) for computing expectation values 〈S_A · S_B〉 that appear in the Heisenberg spin model has been extended to semiempirical single determinant wave functions. An alternative derivation of expectation values in restricted and unrestricted cases is given that takes advantage of the zero differential overlap (ZDO) approximation. A formal connection between single determinant wave functions (which are not in general spin eigenfunctions) and the Heisenberg spin model was established by demonstrating that energies of single determinants that are eigenfunctions of the local spin operators with eigenvalues corresponding to high‐spin radical centers are given by the same Heisenberg coupling constants {J_AB} that describe the true spin states of the system. Unrestricted single determinant wave functions for transition metal complexes are good approximations of local spin eigenfunctions when the metal d orbitals are local in character and all unpaired electrons on each metal have the same spin (although spins on different metals might be reversed). Good approximations of the coupling constants can then be extracted from local spin expectation values 〈S_A · S_B〉 energies of the single determinant wave functions. Once the coupling constants are obtained, diagonalization of the Heisenberg spin Hamiltonian provides predictions of the energies and compositions of the spin states. A computational method is presented for obtaining coupling constants and spin‐state energies in this way for polynuclear transition metal complexes using the intermediate neglect of differential overlap Hamiltonian parameterized for optical spectroscopy (INDO/S) in the ZINDO program. This method is referred to as ZILSH, derived from ZINDO, Davidson's local spin formalism, and the Heisenberg spin model. Coupling constants and spin ground states obtained for 10 iron complexes containing from 2 to 6 metals are found to agree well with experimental results in most cases. In the case of the complex [Fe₆O₃(OAc)₉(OEt)₂(bpy)₂]⁺, a priori predictions of the coupling constants yield a ground‐state spin of zero, in agreement with variable‐temperature magnetization data, and corroborate spin alignments proposed earlier on the basis of structural considerations. This demonstrates the potential of the ZILSH method to aid in understanding magnetic interactions in polynuclear transition metal complexes. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

A thermochemically competitive local hybrid functional without gradient corrections

Following the suggestion of local hybrid functionals with position-dependent exact-exchange admixture [J. Jaramillo, G. E. Scuseria, and M. Ernzerhof, J. Chem. Phys. 118, 1068 (2003)], a functional that mixes only local and exact exchange plus local correlation has been constructed. With a simple local mixing function for the position dependence, this Lh-SVWN functional provides atomization energies for the G2-1 set that are competitive with currently available state-of-the-art functionals like, e.g., B3LYP. This is achieved without generalized gradient approximations for exchange or correlation.     

Broken symmetry density functional study of a mixed‐valence unsymmetrical dinuclear iron complex

Density functional theory (DFT) calculations with different exchange‐correlation functionals were performed for a mixed valence Fe(II)/Fe(III) binuclear complex with μ‐methoxo and two μ‐carboxylate bridging ligands, (1) with geometry optimizations being performed for all possible spin multiplicities (M_S = 2, 4, 6, 8, and 10). Within the exchange‐correlation functionals studied, only the hybrid GGA functionals B3P and B3LYP and also the pure GGA functional RPBE, predicts the geometry with high spin (S = 9/2) to be more stable than the geometry with low spin state (S = 1/2) by 20 kcal/mol, in agreement with the experimental findings. These functionals also predict the same stability order for the different spin states, being M_S = 10>8>6>2>4. The meta‐GGA functionals TPSS and TPSSh and also the pure GGA functionals BLYP and BP86 predict different stability orders. The computed average EPR g‐tensor, g_av, of 2.03, at the B3LYP level, is in good agreement with the experimental findings. Heisenberg exchange coupling constants, J, were calculated within the broken‐symmetry formalism, at the B3LYP level, showing that the two iron centers are antiferromagnetic coupling, with a very weak coupling constant of about ?7 cm^?1, in good agreement with the experimental value. Additionally, the effect of using different multiplicities of the reference geometries on the computed J value is discussed. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Exchange and Delocalization Effects in [Fe6S6]2+ Superclusters

A theoretical study of Heisenberg exchange and double exchange (delocalization) effects in the iron-sulphur supercluster is presented. Such clusters can play important role in biological systems (proteins and enzymes) acting as so-called active centres. The cluster with valence 2+ can be modelled by two Fe(III) and four Fe(II) ions. An idealized structure of double cubane has been considered instead of a more realistic defected double cubane structure of lower symmetry. Energies of the lowest spin states have been calculated numerically depending on the Heisenberg exchange J _i and double exchange b parameters. Possible spin ground states (S=0, 1, 2, 3, 4, 5) have been predicted. The ground state of a given total spin Sis usually achieved for the intermediate spin value of S ₅₆=4 in the case of fully antiferromagnetic as well as partially ferromagnetic spin interactions. In the case of no double exchange, the ground state with the total spin S=3 should always be observed, while a nonzero hopping effect results in narrowing a parameter region of the ground state. If the double exchange is taken into account, then the spin values depend on the Heisenberg integrals. The model results can be applied in order to interpret many structural and magnetic properties of proteins and enzymes possessing the Fe-S active centres.     

Theoretical study of the relationship between the nearest-neighbor exchange coupling interactions and the number of peripheral complexes in the cyano-bridged CrMn6(CN)6 and CrMn2(CN)2 clusters

A theoretical density functional study of the relationship between the nearest-neighbor constants and the number of peripheral complexes in the cyano-bridged [Cr[CNMn(salen)(H2O)]6]3+ and [(5-Brsalen)2(H2O)2Mn2Cr(CN)6] clusters is presented. Two approaches show that the antiferromagnetic coupling interactions between nearest neighbors decrease with the increase of the number of peripheral complexes, although the second approach provides better results using several exchange-correlation functionals. The first approach consisted of evaluating the exchange coupling constant J(ij) between two paramagnetic metal centers i and j in the hexanuclear molecule by calculating the energy differences between the highest and broken-symmetry spin states of a model molecule in which metal atoms except for i and j are substituted by diamagnetic Zn(II) cations, while the second consisted of calculating the different spin-state energies of hexanuclear complexes and using the Heisenberg Hamiltonian to obtain the exchange coupling constants between different metal centers. Moreover, Kahn's qualitative theory succeeded in being applied to interpret the trend.     

Benchmarking approximate density functional theory. I. s/d excitation energies in 3d transition metal cations

The performance of a number of different implementations of density functional theory (DFT) for predicting the s/d interconfigurational energies of the 3d transition metal cations is investigated. Systematic comparisons of computed results with experimental data indicate that gradient corrected correlation functionals, like the LYP GGA, efficiently correct the flaws of the LDA, but reveal shortcomings in the treatment of exchange by currently available GGAs. The admixture of exact exchange in hybrid functionals eventually leads to largely reduced errors. Several basis sets available for the 3d elements are tested in combination with the B3LYP functional. Finally, the influence of variations of the admixture of exact exchange is systematically tested. The results reveal that computed s/d excitation energies obtained for the individual ions depend in markedly different ways on the amount of exact exchange admixture and that there is no single optimal and transferable exchange parameter to create a hybrid functional that yields improved results for all ions alike. Several of the recently devised functionals perform as good as or slightly better than the B3LYP functional in the present study. But given the fact that the B3LYP functional has been identified as the most successful DFT method in an overwhelming number of systematic investigations in very many areas of chemical research, there is no persuasive motivation to recommend its replacement by one of the other functionals, as much less is known about their robustness.     

两个典型氧桥三核铁(Ⅲ)配合物[Fe3O(O2CCH2OC6H5)6(3H2O)]和[Fe3O(TIEO)2(O2CPh)2Cl3]的局域自旋和磁性质的理论探讨

用从头算波函数(UHF或UDFT波函数)代替ZILSH方法中的半经验波函数得到了ABLSH方法，接着用该方法研究了两个典型氧桥三核铁(Ⅲ)配合物[Fe₃O(O₂CCH₂OC₆H₅)₆(3H₂O)]和[Fe₃O(TIEO)₂(O₂CPh)₂Cl₃]的局域自旋和磁性质。通过计算得到的局域自旋结果和前人的具有可比性，同时所得的磁交换耦合常数和实验值很吻合。该方法可作为研究海森堡型磁性系统(HM)的新工具。     

Hybrid exchange correlation functionals and potentials: Concept elaboration

This paper deals with hybrid functionals that contain exact exchange energy and are the most popular and effective functionals in modern density functional theory. Emphasis is laid on generalization of the notion of a “hybrid functional,” which arises from the introduction of the spatial dependence of the exact exchange admixture (local hybrid functionals). Problems inherent in hybrid functionals are considered along with problems inherent in a wider class of so-called orbital-dependent functionals. In particular, the technique for constructing the local and multiplicative potentials, including the optimized effective potential method, is considered in detail. The theoretical approaches under study are illustrated by calculations of atomization molecular energies and magnetic resonance parameters.     

Redox-active ligands as a challenge for electronic structure methods

To improve the catalytic activity of 3d transition metal catalysts, redox-active ligands are a promising tool. These ligands influence the oxidation state of the metal center as well as the ground spin-state and make the experimental determination of both properties challenging. Therefore, first-principles calculations, in particular employing density functional theory with a proper choice of exchange-correlation (xc) functional, are crucial. Common xc functionals were tested on a simple class of metal complexes: homoleptic, octahedral tris(diimine) iron(II) complexes. The spin-state energy splittings for most of these complexes showed the expected linear dependence on the amount of exact exchange included in the xc functionals. Even though varying redox-activity affects the electronic structure of the complexes considerably, the sensitivity of the spin-state energetics to the exact exchange admixture is surprisingly small. For iron(II) complexes with highly redox-active ligands and for a broad range of ligands in the reduced tris(diimine) iron(I) complexes, self-consistent field convergence to local minima was observed, which differ from the global minimum in the redox state of the ligand. This may also result in convergence to a molecular structure that corresponds to an energetically higher-lying local minimum. One criterion to detect such behavior is a change in the sign of the slope for the dependence of the spin-state energy splittings on the amount of exact exchange. We discuss possible protocols for dealing with such artifacts in cases in which a large number of calculations makes checking by hand unfeasible.     

The performance of nonhybrid density functionals for calculating the structures and spin states of Fe(II) and Fe(III) complexes

The local density approximation and a range of nonhybrid gradient corrected density functionals (PW91, BLYP, PBE, revPBE, RPBE) have been assessed with respect to the prediction of geometries and spin-state energy preferences for a range of homoleptic Fe(II)L6 and Fe(III)L6 complexes, where L = Cl-, CN-, NH3, pyridine, imidazole, H2O, O=CH2 and tetrahydrofuran. While the qualitative spin-state energies from in vacuo structure optimizations are reasonable the geometries are relatively poorly treated, especially for [FeCl6]3-/4-. Structural results for all the complexes are significantly improved by including environmental effects. The best compromise between structural and spin-state predictive accuracy was obtained for the RPBE functional in combination with the COSMO solvation approach. This approach systematically overestimates the energetic preference for a low spin state, which is partly due to the well-known effect of the lack of exact exchange in nonhybrid functionals and partly due to the larger solvation stabilization of low-spin complexes that have shorter bond lengths and thus smaller molecular volumes than their high-spin partners. Calculations on low spin [Fe(bipy)3]2+ and [Fe(phen)3]2+ and their ortho methyl substituted analogs, which are high spin at room temperature but cross over to low spin at low temperature, suggest the RPBE/COSMO combination generates low spin states which are too stable by approximately 13 kcal mol(-1).     

Importance of the correlation contribution for local hybrid functionals: range separation and self-interaction corrections

Local hybrid functionals with their position-dependent exact-exchange admixture are a conceptually simple and promising extension of the concept of a hybrid functional. Local hybrids based on a simple mixing of the local spin density approximation (LSDA) with exact exchange have been shown to be successful for thermochemistry, reaction barriers, and a range of other properties. So far, the combination of this generation of local hybrids with an LSDA correlation functional has been found to give the most favorable results for atomization energies, for a range of local mixing functions (LMFs) governing the exact-exchange admixture. Here, we show that the choice of correlation functional to be used with local hybrid exchange crucially influences the parameterization also of the exchange part as well as the overall performance. A novel ansatz for the correlation part of local hybrids is suggested based on (i) range-separation of LSDA correlation into short-range (SR) and long-range (LR) parts, and (ii) partial or full elimination of the one-electron self-correlation from the SR part. It is shown that such modified correlation functionals allow overall larger exact exchange admixture in thermochemically competitive local hybrids than before. This results in improvements for reaction barriers and for other properties crucially influenced by self-interaction errors, as demonstrated by a number of examples. Based on the range-separation approach, a fresh view on the breakdown of the correlation energy into dynamical and non-dynamical parts is suggested.     

DFT calculations on the spin-crossover complex Fe(salen)(NO): a quest for the best functional

DFT calculations on the spin-crossover complex Fe(salen)(NO) provide a striking illustration of the comparative performance of different exchange-correlation functionals vis-à-vis the issue of transition metal spin state energetics. Thus, although the "classic" pure functionals PW91 and BLYP favor the S = 1/2 state by about 10 kcal/mol, relative to the S = 3/2 state, the hybrid functional B3LYP favors the latter state by nearly the same margin. In contrast, the newer pure functionals OLYP and OPBE, based on the OPTX exchange functional, as well as the B3LYP* hybrid functional (which has 15% Hartree-Fock exchange, compared with 20% for B3LYP) predict nearly isoenergetic S = 1/2 and 3/2 states, as required for a spin-crossover complex. Intriguingly, the OLYP and B3LYP* spin density profiles for the S = 1/2 state of Fe(salen)(NO) are substantially dissimilar.     

Bioinorganic chemistry modeled with the TPSSh density functional

In this work, the TPSSh density functional has been benchmarked against a test set of experimental structures and bond energies for 80 transition-metal-containing diatomics. It is found that the TPSSh functional gives structures of the same quality as other commonly used hybrid and nonhybrid functionals such as B3LYP and BP86. TPSSh gives a slope of 0.99 upon linear fitting to experimental bond energies, whereas B3LYP and BP86, representing 20% and 0% exact exchange, respectively, give linear fits with slopes of 0.91 and 1.07. Thus, TPSSh eliminates the large systematic component of the error in other functionals, reducing rms errors from 46-57 to 34 kJ/mol. The nonhybrid version of the functional, TPSS, gives a slope of 1.08, similar to BP86, implying that using 10% exact exchange is the main reason for the success of TPSSh. Typical bioinorganic reactions were then investigated, including spin inversion and electron affinity in iron-sulfur clusters, and breaking or formation of bonds in iron proteins and cobalamins. The results show that differences in reaction energies due to exact exchange can be much larger than the usually cited approximately 20 kJ/mol, sometimes exceeding 100 kJ/mol. The TPSSh functional provides energies approximately halfway between nonhybrids BP86 and TPSS, and 20% exact exchange hybrid B3LYP: Thus, a linear correlation between the amount of exact exchange and the numeric value of the reaction energy is observed in all these cases. For these reasons, TPSSh stands out as a most promising density functional for use and further development within the field of bioinorganic chemistry.     

About the calculation of exchange coupling constants using density-functional theory: the role of the self-interaction error

The effect of the correction of the self-interaction error on the calculation of exchange coupling constants with methods based on density-functional theory has been tested in simple model systems. The inclusion of the self-interaction correction cancels the nondynamical correlation energy contributions simulated by the commonly used functionals. Hence, such correction should be important in the accurate determination of exchange coupling constants. We have also tested several recent functionals to calculate exchange coupling constants in transition-metal complexes, such as meta-GGA functionals or new formulations of hybrid functionals. The influence of the basis set and of the use of pseudopotentials on the calculated J values has also been evaluated for a Fe(III) dinuclear complex in which the paramagnetic centers bear several unpaired electrons.     

Assessment of a Coulomb-attenuated exchange-correlation energy functional

The recently proposed CAM-B3LYP exchange-correlation energy functional, based on a partitioning of the r operator in the exchange interaction into long- and short-range components, is assessed for the determination of molecular thermochemistry, structures, and second order response properties. Rydberg and charge transfer excitation energies and static electronic polarisabilities are notably improved over the standard B3LYP functional; classical reaction barriers also improve. Ionisation potentials, bond lengths, NMR shielding constants and indirect spin-spin coupling constants are comparable with the two functionals. CAM-B3LYP atomisation energies and diatomic harmonic vibrational wavenumbers are less accurate than those of B3LYP. Future research directions are outlined.     

Trinuclear copper(II) complexes derived from Schiff-base ligands based on a 6-amino-6-deoxyglucopyranoside: structural and magnetic characterization

The trinuclear copper(II) complexes ([CuL1)(mu-ac)Cu(mu-ac)CuL1) (1) and ([CuL2)(mu-ac)Cu(mu-ac)CuL2) (2) of the tridentate aminosaccharide-derived Schiff-base ligands H2L1 [6-N-(salicylidene)amino-6-deoxy-1,2,3-tri-O-methyl-alpha-D-glucopyranoside] and H2L2 [6-N-(3,5-di-tert-butylsalicylidene)amino-6-deoxy-1,2,3-tri-O-methyl-alpha-D-glucopyranoside] were synthesized and structurally characterized. The trinuclear complex units can be described as two terminal copper-ligand moieties bridged by a central copper acetate moiety, with the Cu centers arranged in a triangular fashion. IR and UV/vis spectroscopic studies strongly indicate that the trinuclear structure is maintained in a methanolic solution. The temperature dependence of the magnetic susceptibility of both complexes shows a moderate antiferromagnetic coupling and can be well interpreted by applying a symmetric Cua-Cub-Cua' model with linear spin topology. The fit of the magnetic data affords coupling constants J of -34 and -24 cm(-1) for 1 and 2, respectively [H = -J(SaSb + SbSa')]. For mu-alkoxo-mu-acetato-bridged copper(II) complexes with a large dihedral angle between the adjacent coordination planes, as found in 1 and 2, such an antiferromagnetic coupling is unusual. However, density functional theory calculations of 2 using BP86, B3LYP*, and B3LYP density functionals confirmed a symmetric doublet ground state.