Multireference configuration interaction calculations of electronicg-tensors for NO2, H2O+, and CO+

State-Interaction Approach for Evaluating g-Tensors within EOM-CC and RAS-CI Frameworks: Theory and Benchmarks

Among various techniques designed for studying open-shell species, electron paramagnetic resonance (EPR) spectroscopy plays an important role. The key quantity measured by EPR is the g-tensor, describing the coupling between an external magnetic field and molecular electronic spin. One theoretical framework for quantum chemistry calculations of g-tensors is based on response theory, which involves substantial developments that are specific to the underlying electronic structure models. A simplified and easier-to-implement approach is based on the state-interaction scheme, in which perturbation is included by considering a small number of states. We describe and benchmark the state-interaction approach using equation-of-motion coupled-cluster and restricted-active-space configuration interaction wave functions. The analysis confirms that this approach can deliver accurate results and highlights caveats of applying it, such as a choice of the reference state, convergence with respect to the number of states used in calculations, etc. The analysis also contributes toward a better understanding of challenges in calculations of higher-order properties using approximate wave functions.

A relativistic configuration interaction method with general expansions and initial applications to electronic g-factors

A new implementation of the orbital-based two-component relativistic configuration interaction approach is reported and applied to calculations of the electronic g-shifts of three diatomic radicals: AlO, HgF, and PdH. The new implementation augments efficient routines for the calculation of nonrelativistic Hamiltonians with new vectorized routines for the calculation of the action of the one-electron spin-orbit operator and allows efficient calculations for the expansion of generalized active space type. The program makes full use of double group as well as time-reversal symmetry. Particle-hole reorganization of the operators is used to improve the efficiency for expansions with nearly fully occupied orbital spaces. The flexibility of the algorithm and program is used to investigate the convergence of electronic g-shifts for the three diatomic radicals as functions of the active space, states included in the orbital optimization, and excitation levels. It was possible to converge to the valence limits within a few percent using expansions containing up to quadruple excitations. However, when excitations from the core orbitals were added, it was not possible to demonstrate convergence to within a few percent with expansions containing at most 10 × 109 determinants.

Application of multireference linearized coupled-cluster theory to atomic and molecular systems

Applications of the multireference linearized coupled-cluster single-doubles (MRLCCSD) to atomic and molecular systems have been carried out. MRLCCSD is exploited to calculate the ground-state energies of HF, H2O, NH3, CH4, N2, BF, and C2with basis sets, cc-pVDZ, cc-pVTZ and cc-pVQZ. The equilibrium bond lengths and vibration frequencies of HF, HCl, Li2, LiH, LiF, LiBr, BH, and AlF are computed with MRLCCSD and compared with the experimental data. The electron affinities of F and CH as well as the proton affinities of H2O and NH3are also calculated with MRLCCSD. These results are compared with the results produced with second-order perturbation theory, linearized coupled-cluster doubles (LCCD), coupled-cluster doubles (CCD), coupled-cluster singles and doubles (CCSD), CCSD with perturbative triples correction (CCSD(T)). It is shown that all results obtained with MRLCCSD are reliable and accurate.

Benchmark coupled-cluster g-tensor calculations with full inclusion of the two-particle spin-orbit contributions

We present a parallel implementation to compute electron spin resonance g-tensors at the coupled-cluster singles and doubles (CCSD) level which employs the ACES III domain-specific software tools for scalable parallel programming, i.e., the super instruction architecture language and processor (SIAL and SIP), respectively. A unique feature of the present implementation is the exact (not approximated) inclusion of the five one- and two-particle contributions to the g-tensor [i.e., the mass correction, one- and two-particle paramagnetic spin-orbit, and one- and two-particle diamagnetic spin-orbit terms]. Like in a previous implementation with effective one-electron operators [J. Gauss et al., J. Phys. Chem. A 113, 11541-11549 (2009)], our implementation utilizes analytic CC second derivatives and, therefore, classifies as a true CC linear-response treatment. Therefore, our implementation can unambiguously appraise the accuracy of less costly effective one-particle schemes and provide a rationale for their widespread use. We have considered a large selection of radicals used previously for benchmarking purposes including those studied in earlier work and conclude that at the CCSD level, the effective one-particle scheme satisfactorily captures the two-particle effects less costly than the rigorous two-particle scheme. With respect to the performance of density functional theory (DFT), we note that results obtained with the B3LYP functional exhibit the best agreement with our CCSD results. However, in general, the CCSD results agree better with the experimental data than the best DFT/B3LYP results, although in most cases within the rather large experimental error bars.

QM/MM investigations of organic chemistry oriented questions

About 35 years after its first suggestion, QM/MM became the standard theoretical approach to investigate enzymatic structures and processes. The success is due to the ability of QM/MM to provide an accurate atomistic picture of enzymes and related processes. This picture can even be turned into a movie if nuclei-dynamics is taken into account to describe enzymatic processes. In the field of organic chemistry, QM/MM methods are used to a much lesser extent although almost all relevant processes happen in condensed matter or are influenced by complicated interactions between substrate and catalyst. There is less importance for theoretical organic chemistry since the influence of nonpolar solvents is rather weak and the effect of polar solvents can often be accurately described by continuum approaches. Catalytic processes (homogeneous and heterogeneous) can often be reduced to truncated model systems, which are so small that pure quantum-mechanical approaches can be employed. However, since QM/MM becomes more and more efficient due to the success in software and hardware developments, it is more and more used in theoretical organic chemistry to study effects which result from the molecular nature of the environment. It is shown by many examples discussed in this review that the influence can be tremendous, even for nonpolar reactions. The importance of environmental effects in theoretical spectroscopy was already known. Due to its benefits, QM/MM can be expected to experience ongoing growth for the next decade.In the present chapter we give an overview of QM/MM developments and their importance in theoretical organic chemistry, and review applications which give impressions of the possibilities and the importance of the relevant effects. Since there is already a bunch of excellent reviews dealing with QM/MM, we will discuss fundamental ingredients and developments of QM/MM very briefly with a focus on very recent progress. For the applications we follow a similar strategy.

First-principles calculation of parameters of electron paramagnetic resonance spectroscopy in solids

The hyperfine A-tensor and Zeeman g-tensor parameterize the interaction of an 'effective' electron spin with the magnetic field due to the nuclear spin and the homogeneous external magnetic field, respectively. The A- and g-tensors are the quantities of primary interest in electron paramagnetic resonance (EPR) spectroscopy. In this paper, we review our work [E.S. Kadantsev, T. Ziegler, J. Phys. Chem. A 2008, 112, 4521; E. S. Kadantsev, T. Ziegler, J. Phys. Chem. A 2009, 113, 1327] on the calculation of these EPR parameters under periodic boundary conditions (PBC) from first-principles. Our methodology is based on the Kohn-Sham DFT (KS DFT), explicit usage of Bloch basis set made up of numerical and Slater-type atomic orbitals (NAOs/STOs), and is implemented in the 'full potential' program BAND. Our implementation does not rely on the frozen core approximation. The NAOs/STOs basis is well suited for the accurate representation of the electron density near the nuclei, a prerequisite for the calculation of highly accurate hyperfine parameters. In the case of g-tensor, our implementation is based on the method of Van Lenthe et al. [E. van Lenthe, P. E. S. Wormer, A. van der Avoird, J. Chem. Phys. 1997, 107, 2488] in which the spin-orbital coupling is taken into account variationally. We demonstrate the viability of our scheme by calculating EPR parameters of paramagnetic defects in solids. We consider the A-tensor of 'normal' and 'anomalous' muonium defect in IIIA-VA semiconductors as well as the S2 anion radical in KCl host crystal lattice.

Multireference ab initio calculations of g tensors for trinuclear copper clusters in multicopper oxidases

EPR spectroscopy has proven to be an indispensable tool in elucidating the structure of metal sites in proteins. In recent years, experimental EPR data have been complemented by theoretical calculations, which have become a standard tool of many quantum chemical packages. However, there have only been a few attempts to calculate EPR g tensors for exchange-coupled systems with more than two spins. In this work, we present a quantum chemical study of structural, electronic, and magnetic properties of intermediates in the reaction cycle of multicopper oxidases and of their inorganic models. All these systems contain three copper(II) ions bridged by hydroxide or O(2-) anions and their ground states are antiferromagnetically coupled doublets. We demonstrate that only multireference methods, such as CASSCF/CASPT2 or MRCI can yield qualitatively correct results (compared to the experimental values) and consider the accuracy of the calculated EPR g tensors as the current benchmark of quantum chemical methods. By decomposing the calculated g tensors into terms arising from interactions of the ground state with the various excited states, the origin of the zero-field splitting is explained. The results of the study demonstrate that a truly quantitative prediction of the g tensors of exchange-coupled systems is a great challenge to contemporary theory. The predictions strongly depend on small energy differences that are difficult to predict with sufficient accuracy by any quantum chemical method that is applicable to systems of the size of our target systems.

An Alternative Approach to the g-Matrix: Theory and Applications

Starting from the formula proposed by Gerloch and McMeeking in 1975, the electronic g-matrix is expressed as a sum of two matrices called Lambda and Sigma describing the orbital and spin contributions respectively. This approach is applied on benchmark diatomic and triatomic molecules, and on TiF3 and Cu(NH3)4(2+) using either CASPT2 or CCSD(T) methods to calculate the spin-free states and SO-RASSI to calculate spin-orbit coupling. Results compare very well to experimental data and to previous theoretical work; and, for each molecule, the anisotropy of the g-matrix is modeled by the mean of a few parameters.

Theoretical Studies on Dications and Trications of FH, ClH, and BrH. Properties of the Bound 15Σ- States. Electron-Spin g-Factors and Fine/Hyperfine Constants of the Metastable X3Σ- States of ClH2+ and BrH2+

This theoretical study reports calculations on the fine and hyperfine structure parameters of the metastable X(3)Sigma(-)(sigma(2)pi(2)) state of ClH(2+) and BrH(2+). Data on the repulsive FH(2+) system are also included for comparison purposes. The hyperfine structure (hfs) coupling constants for magnetic (A(iso), A(dip)) and quadrupole (eQq) interactions are evaluated using B3LYP, MP4SDQ, CCSD, and QCISD methods and several basis sets. The fine structure (fs) constants (zero-field splitting lambda and spin-rotation coupling gamma) and electron-spin magnetic moments (g-factor) are evaluated in 2nd-order perturbation theory using multireference CI (MRCI) wave functions. Our calculations find for (35)Cl of ClH(2+) A(iso)/A(dip) = 110/-86 MHz; eQq(0) = -59 MHz; 2lambda = 20.4 cm(-1); g( perpendicular)(v = 0) = 2.02217; and gamma = -0.31 cm(-1) (to be compared with the available experimental A(iso)/A(dip)= 162/-30 MHz). For (79)BrH(2+), the corresponding values are 300/-400 MHz; 368 MHz; 362.6 cm(-1); 2.07302; and -0.98 cm(-1) (experimental 2lambda = 445(+/-80) cm(-1)). We find g( perpendicular)(ClH(2+)) to increase by about 0.0054 between v = 0 and 2, whereas the experimental effective g( perpendicular) changes drastically with vibrational excitation. Nuclear quadrupole coupling constants for halogen atoms X are found to be as large as corresponding A(dip)(X)'s, indicating that both terms may have to be included in the Hamiltonian used to interpret XH(2+) hyperfine spectra. A novel finding relates to the bound character of the 1(5)Sigma(-)(sigmapi(2)sigma) state in FH(2+), as already known for ClH(2+) and BrH(2+), but having a deeper potential well D(e) approximately 4,000 cm(-1) (versus 1,000 cm(-1) in the heavier radicals). Vertical ionization potentials for formation of XH(3+) trications are also discussed.

Rare Gas Effects on Hyperfine Coupling Constants of BO, AlO, and GaO

Using density functional theory methods and large basis sets, we calculated hyperfine coupling constants (HFCCs) for the (11)B, (17)O, (27)Al, and (69)Ga nuclei of the radicals BO, AlO, and GaO (XO), embedded in 2-14 rare gas (Rg) Ne and Ar atoms. Kr atoms were included for AlO. The distance of the Rg atoms from XO was varied from 4 to 12 bohr. Matrix effects cause A(iso)(X) to increase, accompanied by decreases in A(dip)(X) and A(dip)(O), while A(iso)(O) remains close to zero. Changes are largest for AlO, slightly smaller for GaO, and very small for BO, in line with the molecular polarizabilities. Observed changes of A(iso)(X) and A(dip)(X) for BO in Ne matrixes and for AlO in Ne, Ar, and Kr matrixes are reproduced in complexes with 12 Rg atoms at distances of 5-6 bohr or 14 Rg atoms at distances of 6-7 bohr. For GaO, experimental data are available only in Ne matrixes. Theoretical results obtained for HFCCs of (17)O could not be verified due to insufficient experimental information. Estimates of HFCCs in matrixes not yet experimentally studied and for GaO in the gas phase have been made. Due to the interaction with rare gas atoms, p-spin density on the X and O atoms of XO is converted into s-spin density on X, thereby causing an increase (in magnitude) of A(iso)(X), accompanied by decreases in A(dip) of X and O. The higher polarizability of XO along the bond axis is reflected in complexes that have axial Rg atoms showing larger changes in HFCCs than comparable complexes without axial Rg atoms.

Electron spin resonance g tensors for complexes of Ne and Ar with AlO: theoretical studies related to the large matrix effect observed for AlO

For Ne(n)-AlO (n=2, 4, 6, 8, 10) and Ar(n)-AlO clusters (n=2, 4, 6, 8), the perpendicular (relative to AlO) component of the g tensor was calculated by second-order perturbation theory, using multireference configuration-interaction wave functions. The rare-gas (Rg) atoms were placed axially and/or off axially (one or two rings of four Rg atoms each), and the distance of the Rg atoms from the Al and O atoms, or from the AlO axis, was varied from 4 to 12 bohrs. Rg atoms placed axially mostly increase g(perpendicular), whereas off-axially placed ones lower it below the gas-phase value of AlO. The largest deviations from g(perpendicular) of isolated AlO occur at Ne-Al,O distances of 5-6 bohrs, and Ar-Al,O distances of 6-9 bohrs, with maximal lowerings of about 1600 ppm for Ne and about 2200 ppm (estimated) for Ar in the case of two axial and eight off-axial Rg atoms. Electron spin resonance studies by Knight and Weltner found large matrix effects for AlO, with downshifts of g(perpendicular) observed to be about 450 and 1150 ppm in Ne and Ar matrices, respectively.