A combined DFT and restricted open-shell configuration interaction method including spin-orbit coupling: application to transition metal L-edge X-ray absorption spectroscopy.

A novel restricted-open-shell configuration interaction with singles (ROCIS) approach for the calculation of transition metal L-edge X-ray absorption spectra is introduced. In this method, one first calculates the ground state and a number of excited states of the non-relativistic Hamiltonian. By construction, the total spin is a good quantum number in each of these states. For a ground state with total spin S excited states with spin S' = S, S - 1, and S + 1 are constructed. Using Wigner-Eckart algebra, all magnetic sublevels with MS = S,..., -S for each multiplet of spin S are obtained. The spin-orbit operator is represented by a mean-field approximation to the full Breit-Pauli spin-orbit operator and is diagonalized over this N-particle basis. This is equivalent to a quasi-degenerate treatment of the spin-orbit interaction to all orders. Importantly, the excitation space spans all of the molecular multiplets that arise from the atomic Russell-Saunders terms. Hence, the method represents a rigorous first-principles approach to the complicated low-symmetry molecular multiplet problem met in L-edge X-ray absorption spectroscopy. In order to gain computational efficiency, as well as additional accuracy, the excitation space is restricted to single excitations and the configuration interaction matrix is slightly parameterized in order to account for dynamic correlation effects in an average way. To this end, it is advantageous to employ Kohn-Sham rather than Hartree-Fock orbitals thus defining the density functional theory∕ROCIS method. However, the method can also be used in an entirely non-empirical fashion. Only three global empirical parameters are introduced and have been determined here for future application of the method to any system containing any transition metal. The three parameters were carefully calibrated using the L-edge X-ray absorption spectroscopy spectra of a test set of coordination complexes containing first row transition metals. These parameters are universal and transferable. Hence, there are no adjustable parameters that are used to fit experimental X-ray absorption spectra. Thus, the new approach classifies as a predictive first-principles method rather than an analysis tool. A series of calculations on transition metal compounds containing Cu, Ti, Fe, and Ni in various oxidation and spin states is investigated and a detailed comparison to experimental data is reported. In most cases, the approach yields good to excellent agreement with experiment. In addition, the origin of the observed spectral features is discussed in terms of the electronic structure of the investigated compounds.

[1]  F. Neese,et al.  First principles calculations of the structure and V L-edge X-ray absorption spectra of V2O5 using local pair natural orbital coupled cluster theory and spin-orbit coupled configuration interaction approaches. , 2013, Physical chemistry chemical physics : PCCP.

[2]  Frank Neese,et al.  Excited states of large open-shell molecules: an efficient, general, and spin-adapted approach based on a restricted open-shell ground state wave function. , 2013, The journal of physical chemistry. A.

[3]  Michael Odelius,et al.  Ab Initio Calculations of X-ray Spectra: Atomic Multiplet and Molecular Orbital Effects in a Multiconfigurational SCF Approach to the L-Edge Spectra of Transition Metal Complexes. , 2012, The journal of physical chemistry letters.

[4]  F. Neese,et al.  One-electron contributions to the g-tensor for second-order Douglas-Kroll-Hess theory. , 2012, The Journal of chemical physics.

[5]  F. Neese,et al.  Manganese K-edge X-ray absorption spectroscopy as a probe of the metal-ligand interactions in coordination compounds. , 2012, Inorganic chemistry.

[6]  M. Zwierzycki,et al.  Multiplet ligand-field theory using Wannier orbitals , 2011, 1111.4940.

[7]  Frank Neese,et al.  Prediction of high-valent iron K-edge absorption spectra by time-dependent density functional theory. , 2011, Dalton transactions.

[8]  Isao Tanaka,et al.  Ab initio charge transfer multiplet calculations on the L_{2,3} XANES and ELNES of 3d transition metal oxides , 2011 .

[9]  F. D. de Groot,et al.  The CTM4XAS program for EELS and XAS spectral shape analysis of transition metal L edges. , 2010, Micron.

[10]  F. Neese,et al.  Accurate modeling of spin-state energetics in spin-crossover systems with modern density functional theory. , 2010, Inorganic chemistry.

[11]  F. Neese Prediction of molecular properties and molecular spectroscopy with density functional theory: From fundamental theory to exchange-coupling , 2009 .

[12]  Eli Stavitski,et al.  Multiplet calculations of L2,3 x-ray absorption near-edge structures for 3d transition-metal compounds , 2009, Journal of physics. Condensed matter : an Institute of Physics journal.

[13]  T. Möller,et al.  Size dependence of L2,3 branching ratio and 2p core-hole screening in x-ray absorption of metal clusters. , 2008, Physical review letters.

[14]  Frank Neese,et al.  Prediction of iron K-edge absorption spectra using time-dependent density functional theory. , 2008, The journal of physical chemistry. A.

[15]  Frank Neese,et al.  All-Electron Scalar Relativistic Basis Sets for Third-Row Transition Metal Atoms. , 2008, Journal of chemical theory and computation.

[16]  H. Freund,et al.  A new analysis of X-ray adsorption branching ratios: Use of Russell-Saunders coupling , 2008 .

[17]  F. D. Groot,et al.  Ligand and metal X-ray absorption in transition metal complexes , 2008 .

[18]  Frank Neese,et al.  Description of the ground-state covalencies of the bis(dithiolato) transition-metal complexes from X-ray absorption spectroscopy and time-dependent density-functional calculations. , 2007, Chemistry.

[19]  Frank Neese,et al.  Advanced aspects of ab initio theoretical optical spectroscopy of transition metal complexes: Multiplets, spin-orbit coupling and resonance Raman intensities , 2007 .

[20]  Ruiqin Q. Zhang,et al.  Effective simulation of biological systems: Choice of density functional and basis set for heme-containing complexes , 2007 .

[21]  Shubin Liu,et al.  Towards understanding performance differences between approximate density functionals for spin states of iron complexes. , 2006, The Journal of chemical physics.

[22]  K. Pierloot,et al.  Relative energy of the high-(5T2g) and low-(1A1g) spin states of [Fe(H2O)6]2+, [Fe(NH3)6]2+, and [Fe(bpy)3]2+: CASPT2 versus density functional theory. , 2006, The Journal of chemical physics.

[23]  K. Hodgson,et al.  Fe L-edge XAS studies of K4[Fe(CN)6] and K3[Fe(CN)6]: a direct probe of back-bonding. , 2006, Journal of the American Chemical Society.

[24]  Frank Neese,et al.  First-principles calculations of zero-field splitting parameters. , 2006, The Journal of chemical physics.

[25]  F. Weigend Accurate Coulomb-fitting basis sets for H to Rn. , 2006, Physical chemistry chemical physics : PCCP.

[26]  F. Weigend,et al.  Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. , 2005, Physical chemistry chemical physics : PCCP.

[27]  Yukinori Koyama,et al.  First-principles multielectron calculations of Ni L 2 , 3 NEXAFS and ELNES for Li Ni O 2 and related compounds , 2005 .

[28]  E. Solomon,et al.  Metal and ligand K-edge XAS of organotitanium complexes: metal 4p and 3d contributions to pre-edge intensity and their contributions to bonding. , 2005, Journal of the American Chemical Society.

[29]  Frank Neese,et al.  Efficient and accurate approximations to the molecular spin-orbit coupling operator and their use in molecular g-tensor calculations. , 2005, The Journal of chemical physics.

[30]  Frank Neese,et al.  A spectroscopy oriented configuration interaction procedure , 2003 .

[31]  K. Hodgson,et al.  L-edge X-ray absorption spectroscopy of non-heme iron sites: experimental determination of differential orbital covalency. , 2003, Journal of the American Chemical Society.

[32]  F. Neese Metal and ligand hyperfine couplings in transition metal complexes: The effect of spin-orbit coupling as studied by coupled perturbed Kohn-Sham theory , 2003 .

[33]  Abhik Ghosh,et al.  High-level ab initio calculations on the energetics of low-lying spin states of biologically relevant transition metal complexes: a first progress report. , 2003, Current opinion in chemical biology.

[34]  Markus Reiher,et al.  The generalized Douglas–Kroll transformation , 2002 .

[35]  Daniel Maynau,et al.  Metal-ligand delocalization in magnetic orbitals of binuclear complexes , 2002 .

[36]  J. Conesa,et al.  Density functional calculations for modeling the oxidized states of the active site of nickel-iron hydrogenases. 1. Verification of the method with paramagnetic Ni and Co complexes. , 2002, Inorganic chemistry.

[37]  Y. Koyama,et al.  Relativistic cluster calculation of ligand-field multiplet effects on cation L-2,L-3 x-ray-absorption edges of SrTiO3, NiO, and CaF2 , 2001 .

[38]  K. Hodgson,et al.  Relationship between the Dipole Strength of Ligand Pre-Edge Transitions and Metal-Ligand Covalency. , 1999, Inorganic chemistry.

[39]  S. Grimme,et al.  A COMBINATION OF KOHN-SHAM DENSITY FUNCTIONAL THEORY AND MULTI-REFERENCE CONFIGURATION INTERACTION METHODS , 1999 .

[40]  T. Hatsui,et al.  Polarized Ni K- and L-edge and S K-edge XANES study of [Ni(III)(mnt)2]1-. , 1999, Journal of synchrotron radiation.

[41]  T. Hatsui,et al.  Strong metal-to-ligand charge transfer bands observed in Ni K- and L-edge XANES of planar Ni complexes. , 1999, Journal of synchrotron radiation.

[42]  F. Neese,et al.  Calculation of Zero-Field Splittings, g-Values, and the Relativistic Nephelauxetic Effect in Transition Metal Complexes. Application to High-Spin Ferric Complexes. , 1998, Inorganic chemistry.

[43]  Christoph van Wüllen,et al.  Molecular density functional calculations in the regular relativistic approximation: Method, application to coinage metal diatomics, hydrides, fluorides and chlorides, and comparison with first-order relativistic calculations , 1998 .

[44]  P. Wormer,et al.  Density functional calculations of molecular hyperfine interactions in the zero order regular approximation for relativistic effects , 1998 .

[45]  B. Hoffman,et al.  Pulsed ENDOR and ESEEM Study of [Bis(maleonitriledithiolato)nickel](-): An Investigation into the Ligand Electronic Structure. , 1998, Inorganic chemistry.

[46]  Florian Weigend,et al.  Auxiliary basis sets for main row atoms and transition metals and their use to approximate Coulomb potentials , 1997 .

[47]  K. Hodgson,et al.  Cu K-Edge XAS Study of the [Cu2(μ-O)2] Core: Direct Experimental Evidence for the Presence of Cu(III) , 1997 .

[48]  J. Rehr,et al.  TI K-EDGE XANES STUDIES OF TI COORDINATION AND DISORDER IN OXIDE COMPOUNDS: COMPARISON BETWEEN THEORY AND EXPERIMENT , 1997 .

[49]  K. Hodgson,et al.  A Multiplet Analysis of Fe K-Edge 1s → 3d Pre-Edge Features of Iron Complexes , 1997 .

[50]  Stefan Grimme,et al.  Density functional calculations with configuration interaction for the excited states of molecules , 1996 .

[51]  Christel M. Marian,et al.  A mean-field spin-orbit method applicable to correlated wavefunctions , 1996 .

[52]  U. Wahlgren,et al.  A new mean-field and ECP-based spin-orbit method. Applications to Pt and PtH , 1996 .

[53]  Marco Häser,et al.  Auxiliary basis sets to approximate Coulomb potentials , 1995 .

[54]  K. Hodgson,et al.  Ligand K-edge X-ray absorption spectroscopic studies. Metal-ligand covalency in a series of transition metal tetrachlorides , 1995 .

[55]  K. Hodgson,et al.  Ligand K-edge x-ray absorption spectroscopy as a probe of ligand-metal bonding: Charge donation and covalency in copper-chloride systems , 1994 .

[56]  J. Almlöf,et al.  Integral approximations for LCAO-SCF calculations , 1993 .

[57]  Evert Jan Baerends,et al.  Relativistic regular two‐component Hamiltonians , 1993 .

[58]  John F. Stanton,et al.  The equation of motion coupled‐cluster method. A systematic biorthogonal approach to molecular excitation energies, transition probabilities, and excited state properties , 1993 .

[59]  Jean-Paul Malrieu,et al.  Specific CI calculation of energy differences: Transition energies and bond energies , 1993 .

[60]  E. Solomon,et al.  Copper L-edge spectral studies. A direct experimental probe of the ground-state covalency in the blue copper site in plastocyanin , 1993 .

[61]  A. Becke Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[62]  A. Becke A New Mixing of Hartree-Fock and Local Density-Functional Theories , 1993 .

[63]  R. Caballol,et al.  Variational calculation of small energy differences. The singlet-triplet gap in [Cu2Cl6]2− , 1992 .

[64]  A. Kotani,et al.  CHARGE-TRANSFER SATELLITES AND MULTIPLET SPLITTING IN X-RAY PHOTOEMISSION SPECTRA OF LATE TRANSITION-METAL HALIDES , 1992 .

[65]  Thole,et al.  2p x-ray absorption of 3d transition-metal compounds: An atomic multiplet description including the crystal field. , 1990, Physical review. B, Condensed matter.

[66]  C. Breneman,et al.  Determining atom‐centered monopoles from molecular electrostatic potentials. The need for high sampling density in formamide conformational analysis , 1990 .

[67]  Rodney J. Bartlett,et al.  The equation-of-motion coupled-cluster method: Excitation energies of Be and CO , 1989 .

[68]  A. Becke,et al.  Density-functional exchange-energy approximation with correct asymptotic behavior. , 1988, Physical review. A, General physics.

[69]  Thole,et al.  Branching ratio in x-ray absorption spectroscopy. , 1988, Physical review. B, Condensed matter.

[70]  Thole,et al.  Linear relation between x-ray absorption branching ratio and valence-band spin-orbit expectation value. , 1988, Physical review. A, General physics.

[71]  Thole,et al.  Multiplet structure in the L2,3 x-ray-absorption spectra: A fingerprint for high- and low-spin Ni2+ compounds. , 1988, Physical review. B, Condensed matter.

[72]  Ingvar Lindgren,et al.  Diagonalisation of the Dirac Hamiltonian as a basis for a relativistic many-body procedure , 1986 .

[73]  Hess,et al.  Relativistic electronic-structure calculations employing a two-component no-pair formalism with external-field projection operators. , 1986, Physical review. A, General physics.

[74]  Fink,et al.  2p absorption spectra of the 3d elements. , 1985, Physical review. B, Condensed matter.

[75]  K. Wieghardt,et al.  Crystallographic study of the low-spin iron(II) and iron(III) bis complexes of 1,4,7-triazacyclononane , 1985 .

[76]  P. Siegbahn,et al.  The effect of electron correlation on the metal-ligand interaction in iron pentacarbonyl (Fe(CO)5) , 1985 .

[77]  P. Pianetta,et al.  CopperL2,3near-edge structure inCu2O , 1984 .

[78]  P. Fuentealba,et al.  On the reliability of semi-empirical pseudopotentials: simulation of Hartree-Fock and Dirac-Fock results , 1983 .

[79]  S. Suga,et al.  Inner-core excitation spectra of transition-metal compounds: II. p-d absorption spectra , 1982 .

[80]  C. Moore,et al.  Crystal and molecular structures of rinneite, sodium tripotassium hexachloroferrate(II), and hexaamminecobalt(III) hexachloroferrate(III). Comparison of iron-chloride distances in hexachloroferrate(II) and hexachloroferrate(III) , 1982 .

[81]  Shik Shin,et al.  Final State Interactions of Inner Core Absorption Excitation in Transition Metal Halides , 1982 .

[82]  B. Teo Chemical applications of extended x-ray absorption fine structure (EXAFS) spectroscopy , 1980 .

[83]  John R. Sabin,et al.  On some approximations in applications of Xα theory , 1979 .

[84]  M. Krause,et al.  Natural widths of atomic K and L levels, Kα X‐ray lines and several KLL Auger lines , 1979 .

[85]  A. Kobayashi,et al.  One-dimensional System of Square-planar Bis(1,2-dicyanovinylene-1,2-dithiolato)metal Complexes. I. The Crystal Structure of [(C4H9)4N]2[Ni(mnt)2] and [(C2H5)4N][Ni(mnt)2] , 1977 .

[86]  S. Sugano,et al.  Satellites in X-Ray Photoelectron Spectra of Transition-Metal Compounds , 1976 .

[87]  P. Eisenberger,et al.  Observations and interpretation of x-ray absorption edges in iron compounds and proteins. , 1976, Proceedings of the National Academy of Sciences of the United States of America.

[88]  Robert J. Buenker,et al.  Energy extrapolation in CI calculations , 1975 .

[89]  J. Ibers,et al.  Structure of tetramethylammonium tetrachloroferrate(II). Comparison of iron(II) and iron(III) bond lengths in high-spin tetrahedral environments , 1975 .

[90]  Robert J. Buenker,et al.  Individualized configuration selection in CI calculations with subsequent energy extrapolation , 1974 .

[91]  Rajendra P. Gupta,et al.  Calculation of multiplet structure of core p -vacancy levels. II , 1974 .

[92]  Evert Jan Baerends,et al.  Self-consistent molecular Hartree—Fock—Slater calculations I. The computational procedure , 1973 .

[93]  A. Maki,et al.  Electronic ground state of bis (maleonitrile-dithiolene)nickel monoanion. Sulfur-33 hyperfine interaction , 1968 .

[94]  R. Mcweeny On the Origin of Spin‐Hamiltonian Parameters , 1965 .

[95]  N. Edelstein,et al.  Electron Paramagnetic Resonance Studies of the Electronic Structures of Bis(maleonitriledithiolato)copper(II), -nickel(III), -cobalt(II), and -rhodium(II) Complexes , 1964 .

[96]  N. Edelstein,et al.  E.s.r. Studies of Four-Coördinate Complexes of Nickel, Palladium and Platinum Related by Electron Transfer Reactions , 1963 .

[97]  P. Löwdin On the Non‐Orthogonality Problem Connected with the Use of Atomic Wave Functions in the Theory of Molecules and Crystals , 1950 .

[98]  Frank Neese,et al.  The ORCA program system , 2012 .

[99]  F. Neese,et al.  A Modern First-Principles View on Ligand Field Theory Through the Eyes of Correlated Multireference Wavefunctions , 2011 .

[100]  Frank M. F. de Groot,et al.  Multiplet effects in X-ray spectroscopy , 2005 .

[101]  E. Baerends,et al.  g- and a-tensor calculations in the zero-order approximation for relativistic effects of Ni complexes Ni(mnt)(2)(-) and Ni(CO)(3)H as model complexes for the active center of [NiFe]-hydrogenase , 2001 .

[102]  G. Thornton,et al.  X-ray absorption fine structure study of the bound state electronic transitions at the vanadium K and L edges in low symmetry, molecular, vanadium-(IV) and -(V) complexes with oxyoxime and oxyoximate ligands , 1998 .

[103]  Markus P. Fülscher,et al.  Multiconfigurational perturbation theory: Applications in electronic spectroscopy , 1996 .

[104]  J. Sinfelt,et al.  X-ray absorption edge studies of the electronic structure of metal catalysts , 1993 .

[105]  A. Klamt,et al.  COSMO : a new approach to dielectric screening in solvents with explicit expressions for the screening energy and its gradient , 1993 .

[106]  N. Mårtensson,et al.  Core-level binding energies in metals , 1980 .

[107]  E. Davidson The iterative calculation of a few of the lowest eigenvalues and corresponding eigenvectors of large real-symmetric matrices , 1975 .

[108]  Marvin Douglas,et al.  Quantum electrodynamical corrections to the fine structure of helium , 1971 .