A systematic density functional study of the zero-field splitting in Mn(II) coordination compounds

This work presents a detailed evaluation of the performance of density functional theory (DFT) for the prediction of zero-field splittings (ZFSs) in Mn(II) coordination complexes. Eighteen experimentally well characterized four-, five-, and six-coordinate complexes of the general formula [Mn(L)nL'2] with L' = Cl, Br, I, NCS, or N3 (L = an oligodentate ligand) are considered. Several DFT-based approaches for the prediction of the ZFSs are compared. For the estimation of the spin-orbit coupling (SOC) part of the ZFS, it was found that the Pederson-Khanna (PK) approach is more successful than the previously proposed quasi-restricted orbitals (QRO)-based method. In either case, accounting for the spin-spin (SS) interaction either with or without the inclusion of the spin-polarization effects improves the results. This argues for the physical necessity of accounting for this important contribution to the ZFS. On average, the SS contribution represents approximately 30% of the axial D parameters. In addition to the SS part, the SOC contributions of d-d spin flip (alphabeta) and ligand-to-metal charge transfer excited states (betabeta) were found to dominate the SOC part of the D parameter; the observed near cancellation between the alphaalpha and betaalpha parts is discussed in the framework of the PK model. The calculations systematically (correlation coefficient approximately 0.99) overestimate the experimental D values by approximately 60%. Comparison of the signs of calculated and measured D values shows that the signs of the calculated axial ZFS parameters are unreliable once E/D > 0.2. Finally, we find that the calculated D and E/D values are highly sensitive to small structural changes. It is observed that the use of theoretically optimized geometries leads to a significant deterioration of the theoretical predictions relative to the experimental geometries derived from X-ray diffraction. The standard deviation of the theoretical predictions for the D values almost doubles from approximately 0.1 to approximately 0.2 cm-1 upon using quantum chemically optimized structures. We do not find any noticeable improvement in considering basis sets larger than standard double- (SVP) or triple-zeta (TZVP) basis sets or using functionals other than the BP functional.     

Spin-spin contributions to the zero-field splitting tensor in organic triplets, carbenes and biradicals-a density functional and ab initio study

An evaluation study for the direct dipolar electron spin-spin (SS) contribution to the zero-field splitting (ZFS) tensor in electron paramagnetic resonance (EPR) spectroscopy is presented. Calculations were performed on a wide variety of organic systems where the SS contribution to the ZFS dominates over the second-order spin-orbit coupling (SOC) contribution. Calculations were performed using (hybrid) density functional theory (DFT), as well as complete active space self-consistent field (CASSCF) wave functions. In the former case, our implementation is an approximation, because we use the two-particle reduced spin-density matrix of the noninteracting reference system. In the latter case, the SS contribution is approximated by a mean-field method which, nevertheless, gives accurate results, compared to the approximation free computation of the SS part in a CASSCF framework. For the case of the triplet dioxygen molecule, it was shown that restricted open-shell density functional theory (RODFT), as well as CASSCF, can provide accurate spin-spin couplings while spin-unrestricted DFT leads to much larger errors. Furthermore, 15 organic radicals, including several 1,3 and 1,5 diradicals, dinitroxide biradicals, and even a chlorophyll a model system, were examined as test cases to demonstrate the accuracy and efficiency of our approach within a DFT framework. Accurate D values with root-mean-square deviations of 0.0035 cm(-1) were obtained. Furthermore, all trends, including those due to substituent effects, were correctly reproduced. In a different set of calculations, the polyacenes benzene, naphthalene, anthracene, and tetracene were studied. Applying DFT, the absolute D values were noticeably underestimated, but it was possible to correctly reproduce the trend to smaller D values with larger size of the systems. Finally, it was demonstrated that our approach is also well-suited for the study of carbenes. The smaller organic radicals of this work were also studied, through the use of CASSCF wave functions. This was a special advantage in the case of the triplet polyacenes, where the CASSCF approach gave better results than the DFT method. In comparing spin-restricted and spin-unrestricted results, it was shown through a natural orbital analysis and comparison to high-level ab initio calculations that even small amounts of spin polarization introduced by the unrestricted calculations lead to large deviations between the unrestricted Kohn-Sham (UKS) and restricted open-shell Kohn-Sham (ROKS) approaches. It is challenging to understand why the ROKS results show much better correlation with the experimental data.     

Calculation of the zero-field splitting tensor on the basis of hybrid density functional and Hartree-Fock theory

The zero-field splitting (ZFS) (expressed in terms of the D tensor) is the leading spin-Hamiltonian parameter for systems with a ground state spin S>12. To first order in perturbation theory, the ZFS arises from the direct spin-spin dipole-dipole interaction. To second order, contributions arise from spin-orbit coupling (SOC). The latter contributions are difficult to treat since the SOC mixes states of different multiplicities. This is an aspect of dominant importance for the correct prediction of the D tensor. In this work, the theory of the D tensor is discussed from the point of view of analytic derivative theory. Starting from a general earlier perturbation treatment [F. Neese and E. I. Soloman, Inorg. Chem. 37, 6568 (1998)], straightforward response equations are derived that are readily transferred to the self-consistent field (SCF) Hartree-Fock (HF) or density functional theory (DFT) framework. The main additional effort in such calculations arises from the solution of nine sets of nonstandard coupled-perturbed SCF equations. These equations have been implemented together with the spin-orbit mean-field representation of the SOC operator and a mean-field treatment of the direct spin-spin interaction into the ORCA electronic structure program. A series of test calculations on diatomic molecules with accurately known zero-field splittings shows that the new approach corrects most of the shortcomings of previous DFT based methods and, on average, leads to predictions within 10% of the experimental values. The slope of the correlation line is essentially unity for the B3LYP and BLYP functionals compared to approximately 0.5 in previous treatments.     

Ab initio and coupled-perturbed density functional theory estimation of zero-field splittings in MnII transition metal complexes

The paper presents a method comparison for the prediction of zero-field splitting (ZFS) parameters in a series of Mn (II) coordination complexes. The test set consists of Mn (II) complexes that are experimentally well-characterized by X-ray diffraction and high-field electron paramagnetic resonance. Their ZFS parameters have been calculated using density functional theory (DFT) as well as complete active space self-consistent field (CASSCF) methods. It is shown that the recently introduced coupled-perturbed spin-orbit coupling (CP-SOC) approach [ Neese, F. J. Chem. Phys. 2007, 127, 164112 ] together with hybrid-DFT functionals leads to a slope of the correlation line (plot of experimental vs calculated D values) that is essentially unity provided that the direct spin-spin interaction is properly included in the treatment. This is different from our previous DFT study on the same series of complexes where a severe overestimation of the D parameter has been found [ Zein, S. ; Duboc, C. ; Lubitz, W. ; Neese, F. Inorg. Chem. 2008, 47, 134 ]. CASSCF methods have been used to evaluate the ZFS in an "ab initio ligand-field" type treatment. The study demonstrates that a substantial part of the relevant physics is lost in such a treatment since only excitations within the manganese d-manifold are accounted for. Thus, a severe underestimation of the D parameter has been found. Because the CASSCF calculations in combination with quasidegenerate perturbation theory treats the SOC to all orders, we have nevertheless verified that second-order perturbation theory is an adequate approximation in the case of the high-spin d (5) configuration.     

Efficient and accurate approximations to the molecular spin-orbit coupling operator and their use in molecular g-tensor calculations

Approximations to the Breit-Pauli form of the spin-orbit coupling (SOC) operator are examined. The focus is on approximations that lead to an effective quasi-one-electron operator which leads to efficient property evaluations. In particular, the accurate spin-orbit mean-field (SOMF) method developed by Hess, Marian, Wahlgren, and Gropen is examined in detail. It is compared in detail with the "effective potential" spin-orbit operator commonly used in density functional theory (DFT) and which has been criticized for not including the spin-other orbit (SOO) contribution. Both operators contain identical one-electron and Coulomb terms since the SOO contribution to the Coulomb term vanishes exactly in the SOMF treatment. Since the DFT correlation functional only contributes negligibly to the SOC the only difference between the two operators is in the exchange part. In the SOMF approximation, the SOO part is equal to two times the spin-same orbit contribution. The DFT exchange contribution is of the wrong sign and numerically shown to be in error by a factor of 2-2.5 in magnitude. The simplest possible improvement in the DFT-SOC treatment [Veff(-2X)-SOC] is to multiply the exchange contribution to the Veff operator by -2. This is verified numerically in calculations of molecular g-tensors and one-electron SOC constants of atoms and ions. Four different ways of handling the computationally critical Coulomb part of the SOMF and Veff operators are discussed and implemented. The resolution of the identity approximation is virtually exact for the SOC with standard auxiliary basis sets which need to be slightly augmented by steep s functions for heavier elements. An almost as efficient seminumerical approximation is equally accurate. The effective nuclear charge model gives results within approximately 10% (on average) of the SOMF treatment. The one-center approximation to the Coulomb and one-electron SOC terms leads to errors on the order of approximately 5%. Small absolute errors are obtained for the one-center approximation to the exchange term which is consequently the method of choice [SOMF(1X)] for large molecules.     

Investigating magnetostructural correlations in the pseudooctahedral trans-[Ni(II){(OPPh2)(EPPh2)N}2(sol)2] complexes (E = S, Se; sol = DMF, THF) by magnetometry, HFEPR, and ab initio quantum chemistry

In this work, magnetometry and high-frequency and -field electron paramagnetic resonance spectroscopy (HFEPR) have been employed in order to determine the spin Hamiltonian (SH) parameters of the non-Kramers, S = 1, pseudooctahedral trans-[Ni(II){(OPPh(2))(EPPh(2))N}(2)(sol)(2)] (E = S, Se; sol = DMF, THF) complexes. X-ray crystallographic studies on these compounds revealed a highly anisotropic NiO(4)E(2) coordination environment, as well as subtle structural differences, owing to the nature of the Ni(II)-coordinated solvent molecule or ligand E atoms. The effects of these structural characteristics on the magnetic properties of the complexes were investigated. The accurately HFEPR-determined SH zero-field-splitting (zfs) D and E parameters, along with the structural data, provided the basis for a systematic density functional theory (DFT) and multiconfigurational ab initio computational analysis, aimed at further elucidating the electronic structure of the complexes. DFT methods yielded only qualitatively useful data. However, already entry level ab initio methods yielded good results for the investigated magnetic properties, provided that the property calculations are taken beyond a second-order treatment of the spin-orbit coupling (SOC) interaction. This was achieved by quasi-degenerate perturbation theory, in conjunction with state-averaged complete active space self-consistent-field calculations. The accuracy in the calculated D parameters improves upon recovering dynamic correlation with multiconfigurational ab initio methods, such as the second-order N-electron valence perturbation theory NEVPT2, the difference dedicated configuration interaction, and the spectroscopy-oriented configuration interaction. The calculations showed that the magnitude of D (～3-7 cm(-1)) in these complexes is mainly dominated by multiple SOC contributions, the origin of which was analyzed in detail. In addition, the observed largely rhombic regime (E/D = 0.16-0.33) is attributed to the highly distorted metal coordination sphere. Of special importance is the insight by this work on the zfs effects of Se coordination to Ni(II). Overall, a combined experimental and theoretical methodology is provided, as a means to probe the electronic structure of octahedral Ni(II) complexes.     

三角对称晶场中 4 A2(3d3 )态离子EPR参量的SS和SOO机制

A CDM/ EPR program has been developed using Visual Basic 6. 0. The spin-spin（SS）and spin-otherorbit（SOO）interactions omitted in published works have also been included in the Hamiltonian. The CDM/EPR program can study not only the EPR parameters but also the CF energy levels and wavefunctions for 4A2 （3d3）states ions in crystals. Utilizing the CDM/ EPR program,the EPR parameters and fine spectra for Ruby and Emerald have been investigated. The theoretical results are in good agreement with the experimental findings. The contributions to the EPR parameters and fine spectra arising from SS and SOO interactions have been studied. The investigation shows:① The EPR parameters are mainly induced by SO coupling interaction;② The contribution to the zerofield-splitting（ZFS）arising from SS interaction is appreciable and cannot be omitted,whereas the contributions to the ZFS parameter D arising from SOO interaction are smaller;③ The contribution to the Zeeman g-factors and spectra arising from SS and SOO interactions is slight.     

Zero-field splittings and local tilting angles tauMn2+ for Mn2+ in ZnGeP2 and CdGeP2 crystals

The electron paramagnetic resonance (EPR) zero-field splittings (ZFSs) D of Mn2+ in ZnGeP2 and CdGeP2 crystals are calculated from both the microscopic spin-orbit coupling mechanism and the empirical superposition model. From the calculations, the ZFS D of ZnGeP2:Mn2+ is reasonably explained by using the local tilting angle tauMn2+ (rather than the corresponding angle tauZn2+ in the host crystal) and the local tilting angle tauMn2+ (which has not been reported) in CdGeP2:Mn2+ is estimated. The intrinsic ZFS parameter b2(R0) approximately -0.052(6)cm(-1) (with R0 approximately 2.43 angstroms) is suggested for Mn2+-P(3-) combination by using the local tilting angles tauMn2+. The value is quite unlike that (approximately 0.4(2)cm(-1)) obtained in the previous paper by combining the crystallographic data of host crystals with the ZFSs for Mn2+ ions in crystals. The reasonableness of these results is discussed.     

Dynamical mean-field theory from a quantum chemical perspective

We investigate the dynamical mean-field theory (DMFT) from a quantum chemical perspective. Dynamical mean-field theory offers a formalism to extend quantum chemical methods for finite systems to infinite periodic problems within a local correlation approximation. In addition, quantum chemical techniques can be used to construct new ab initio Hamiltonians and impurity solvers for DMFT. Here, we explore some ways in which these things may be achieved. First, we present an informal overview of dynamical mean-field theory to connect to quantum chemical language. Next, we describe an implementation of dynamical mean-field theory where we start from an ab initio Hartree-Fock Hamiltonian that avoids double counting issues present in many applications of DMFT. We then explore the use of the configuration interaction hierarchy in DMFT as an approximate solver for the impurity problem. We also investigate some numerical issues of convergence within DMFT. Our studies are carried out in the context of the cubic hydrogen model, a simple but challenging test for correlation methods. Finally, we finish with some conclusions for future directions.     

Theoretical analysis of the spin Hamiltonian parameters in Co(II)S4 complexes, using density functional theory and correlated ab initio methods

A systematic Density Functional Theory (DFT) and multiconfigurational ab initio computational analysis of the Spin Hamiltonian (SH) parameters of tetracoordinate S = 3/2 Co((II))S(4)-containing complexes has been performed. The complexes under study bear either arylthiolato, ArS(-), or dithioimidodiphosphinato, [R(2)P(S)NP(S)R'(2)](-) ligands. These complexes were chosen because accurate structural and spectroscopic data are available, including extensive Electron Paramagnetic Resonance (EPR)/Electron Nuclear Double Resonance (ENDOR) studies. For comparison purposes, the [Co(PPh(3))(2)Cl(2)] complex, which was thoroughly studied in the past by High-Field and Frequency EPR and Variable Temperature, Variable Field Magnetic Circular Dichroism (MCD) spectroscopies, was included in the studied set. The magnitude of the computed axial zero-field splitting parameter D (ZFS), of the Co((II))S(4) systems, was found to be within ~10% of the experimental values, provided that the property calculation is taken beyond the accuracy obtained with a second-order treatment of the spin-orbit coupling interaction. This is achieved by quasi degenerate perturbation theory (QDPT), in conjunction with complete active space configuration interaction (CAS-CI). The accuracy was increased upon recovering dynamic correlation with multiconfigurational ab initio methods. Specifically, spectroscopy oriented configuration interaction (SORCI), and difference dedicated configuration interaction (DDCI) were employed for the calculation of the D-tensor. The sign and magnitude of parameter D was analyzed in the framework of Ligand Field Theory, to reveal the differences in the electronic structures of the investigated Co((II))S(4) systems. For the axial complexes, accurate effective g'-tensors were obtained in the QDPT studies. These provide a diagnostic tool for the adopted ground state configuration (±3/2 or ±1/2) and are hence indicative of the sign of D. On the other hand, for the rhombic complexes, the determination of the sign of D required the SH parameters to be derived along suitably constructed symmetry interconversion pathways. This procedure, which introduces a dynamic perspective into the theoretical investigation, helped to shed some light on unresolved issues of the corresponding experimental studies. The metal hyperfine and ligand super-hyperfine A-tensors of the C(2) [Co{(SPPh(2))(SP(i)Pr(2))N}(2)] complex were estimated by DFT calculations. The theoretical data were shown to be in good agreement with the available experimental data. Decomposition of the metal A-tensor into individual contributions revealed that, despite the large ZFS, the observed significant anisotropy should be largely attributed to spin-dipolar contributions. The analysis of both, metal and ligand A-tensors, is consistent with a highly covalent character of the Co-S bonds.     

Relativistically corrected hyperfine structure constants calculated with the regular approximation applied to correlation corrected ab initio theory

The infinite-order regular approximation (IORA) and IORA with modified metric (IORAmm) is used to develop an algorithm for calculating relativistically corrected isotropic hyperfine structure (HFS) constants. The new method is applied to the calculation of alkali atoms Li-Fr, coinage metal atoms Cu, Ag, and Au, the Hg(+) radical ion, and the mercury containing radicals HgH, HgCH(3), HgCN, and HgF. By stepwise improvement of the level of theory from Hartree-Fock to second-order M?ller-Plesset theory and to quadratic configuration interaction theory with single and double excitations, isotropic HFS constants of high accuracy were obtained for atoms and for molecular radicals. The importance of relativistic corrections is demonstrated.     

A theoretical study of the fine and hyperfine interactions in the NCO and CNO radicals

The geometries, the harmonic vibrational frequencies, and the Renner-Teller parameter have been reported for the NCO(+)(X (3)Sigma(-)), NCO(X (2)Pi,A (2)Sigma(+),B (2)Pi,2 (2)Sigma(+)), NCO(-)(X (1)Sigma(+)), CNO(+)(X), CNO(X (2)Pi,A (2)Sigma(+),B (2)Pi,2 (2)Sigma(+)), and CNO(-)(X (1)Sigma(+)) systems at the full valence-complete active space self-consistent-field (fv-CASSCF) level of theory. The (2)Pi electronic states of the NCO and CNO radicals have two distinct real vibrational frequencies for the bending modes and these states are subject to the type A Renner-Teller effect. The total energy of CNO(+) without zero point energy correction of the linear geometry is approximately 31 cm(-1) higher than the bent geometry at the fv-CASSCF level and the inversion barrier vanishes after the zero point energy correction; therefore, the ground state of the CNO(+) may possess a quasilinear geometry. The spin-orbit coupling constants estimated using atomic mean field Hamiltonian at the fv-CASSCF level of theory are in better agreement with the experimental values. The excitation energies, the electron affinity, and the ionization potential have been computed at the complete active space second order perturbation theory (CASPT2) and the multireference singles and doubles configuration (MRSD-CI) levels of theory. The computed values of the electric hyperfine coupling constants for the (14)N atom in the ground state of the NCO radical agree well with the experimental data. The magnetic hyperfine coupling constants (HFCC's) have been estimated employing the configuration selected MRSD-CI and the multireference singles configuration interaction (MRS-CI) methods using iterative natural orbitals (ino) as one particle basis. Sufficiently accurate value of the isotropic contribution to the HFCC's can be obtained using an MRS-CI-ino procedure.     

Relationships among coarse-grained field theories of fluctuations in polymer liquids

Two closely related field-theoretic approaches have been used in previous work to construct coarse-grained theories of corrections to the random phase approximation for correlations in block copolymer melts and miscible polymer blends. The "auxiliary field" (AF) approach is based on a rigorous expression for the partition function Z of a coarse-grained model as a functional integral of an auxiliary chemical potential field. The "effective Hamiltonian" (EH) approach is instead based on an expression for Z as a functional integral of an observable order parameter field. The exact effective Hamiltonian H(eff) in the EH approach is defined as the free energy of a system with a constrained order parameter field. In practice, however, H(eff) has often been approximated by a mean-field free energy functional, yielding what we call a mean-field effective Hamiltonian (MFEH) approximation. This approximation was the starting point of both the Fredrickson-Helfand analysis of fluctuation effects in diblock copolymers and earlier work on the Ginzburg criterion in polymer blends. A more rigorous EH approach by Holyst and Vilgis used an auxiliary field representation of the exact H(eff) and allowed for Gaussian fluctuations of this field. All applications of both AF and EH approaches have thus far relied upon some form of Gaussian, or "one-loop" approximation for fluctuations of a chemical potential and/or order parameter field about a mean-field saddle-point. The one-loop EH approximation of Holyst and Vilgis and the one-loop AF theory are equivalent to one another, but not to the one-loop MFEH theory. The one-loop AF and MFEH theories are shown to yield predictions for the inverse structure factor S(-1)(q) that (in the absence of further approximations to either theory) differ by a function that is independent of the Flory-Huggins interaction parameter χ. As a result, these theories yield predictions for the peak scattering intensity that exhibit a similar χ-dependence near a spinodal. The Fredrickson-Helfand theory for the structure factor in disordered diblock copolymer melts is an asymptotic approximation to the MFEH one-loop theory that captures the dominant asymptotic behavior of very long, symmetric copolymers very near the order-disorder transition.     

A ULFC method for d(4)(D(2d)) ions and a study of the spin singlets contributions to zero-field splitting of Cr(2+) ions in zinc sulfide crystals

A simple theoretical method is shown to yield a detailed explanation of numerous EPR parameters for a d4 configuration ion in tetragonal ligand field. Using the unified ligand-field-coupling (ULFC) scheme, the formulas relating the microscopic spin Hamiltonian parameters with the crystal structure parameters are derived. On the basis of the theoretical formulas, the 210 x 210 complete energy matrices including all the spin states are constructed within a strong field representation. By diagonalizing the complete energy matrices, the local lattice structure and Jahn-Teller energy of Cr(2+) ions in ZnS:Cr(2+) system have been investigated. It is found that the theoretical results are in good agreement with the experimental values. Moreover, the contributions of the spin singlets to the zero-field splitting (ZFS) parameters of Cr(2+) ions in ZnS crystals are investigated for the first time. The results indicate that the spin singlets contributions to ZFS parameter b(0)(4) is negligible, but the contributions to ZFS parameters b(0)(4) and b(4)(4) cannot be neglected.     

Solvent effects on optically detected magnetic resonance in triplet spin labels

We have calculated solvent effects on the zero-field splitting (ZFS) constants induced by electron spin–spin coupling (SSC) in the low-lying triplet states of azaaromatic molecules in solutions using multiconfiguration self-consistent-field wave functions and the polarizable continuum model. The second-order spin–orbit coupling (SOC) contribution to the splitting of the ^3* states is found to be almost negligible, and the calculations therefore provide a good estimate of the ZFS parameters and their solvent dependence based only on the electron spin–spin coupling expectation values. The correlation between the shift in the ZFS and the phosphorescence frequency that has been observed in optically detected magnetic resonance experiments in low-temperature glasses is supported by our direct SSC calculations without taking SOC into account. This makes it possible to distinguish between the two theories that earlier were proposed to explain the inhomogeneous broadening of triplet state spectra, and discard the one that is exclusively based on the SOC-induced mixing of the singlet and triplet states.Contribution to the Jacopo Tomasi Honorary IssueAcknowledgments. This work was supported (B. M.) by the Swedish Royal Academy of Science (KVA). This work was also supported by the Norwegian Research Council through a grant of computer time from the Program for Supercomputing. We are grateful to B. Schimmelpfenning for his valuable assistance in the computations.     

SPOCK.CI: a multireference spin-orbit configuration interaction method for large molecules

We present SPOCK.CI, a selecting direct multireference spin-orbit configuration interaction (MRSOCI) program based on configuration state functions. It constitutes an extension of the spin-free density functional theory/multireference configuration interaction (DFT/MRCI) code by Grimme and Waletzke [J. Chem. Phys. 111, 5645 (1999)] and includes spin-orbit interaction on the same footing with electron correlation. Key features of SPOCK.CI are a fast determination of coupling coefficients between configuration state functions, the use of a nonempirical effective one-electron spin-orbit atomic mean-field Hamiltonian, the application of a resolution-of-the-identity approximation to computationally expensive spin-free four-index integrals, and the use of an efficient multiroot Davidson diagonalization scheme for the complex Hamiltonian matrix. SPOCK.CI can be run either in ab initio mode or as semiempirical procedure combined with density functional theory (DFT/MRSOCI). The application of these techniques and approximations makes it possible to compute spin-dependent properties of large molecules in ground and electronically excited states efficiently and with high confidence. Second-order properties such as phosphorescence rates are known to converge very slowly when evaluated perturbationally by sum-over-state approaches. We have investigated the performance of SPOCK.CI on these properties in three case studies on 4H-pyran-4-thione, dithiosuccinimide, and free-base porphin. In particular, we have studied the dependence of the computed phosphorescence lifetimes on various technical parameters of the MRSOCI wave function such as the size of the configuration space, selection of single excitations, diagonalization thresholds, etc. The results are compared to the outcome of extensive quasidegenerate perturbation theory (QDPT) calculations as well as experiment. In all three cases, the MRSOCI approach is found to be superior to the QDPT expansion and yields results in very good agreement with experimental findings. For molecules up to the size of free-base porphin, MRSOCI calculations can easily be run on a single-processor personal computer. Total CPU times for the evaluation of the electronic excitation spectrum and the phosphorescence lifetime of this molecule are below 40 h.     

New perspectives in multireference perturbation theory: the n-electron valence state approach

The n-electron valence state perturbation theory (NEVPT) is a form of multireference perturbation theory which is based on a zero order reference wavefunction of CAS-CI type (complete active space configuration interaction) and which is characterized by the utilization of correction functions (zero order wavefunctions external to the CAS) of multireference nature, obtained through the diagonalization of a suitable two-electron model Hamiltonian (Dyall’s Hamiltonian) in some well defined determinant spaces. A review of the NEVPT approach is presented, starting from the original second order state-specific formulation, going through the quasidegenerate multi-state extension and arriving at the recent implementations of the third order in the energy and of the internally contracted configuration interaction. The chief properties of NEVPT—size consistence and absence of intruder states—are analyzed. Finally, an application concerning the calculation of the vertical spectrum of the biologically important free base porphin molecule, is presented.     

The electronic mean-field configuration interaction method. I. Theory and integral formulas

In this article, we introduce a new method for solving the electronic Schrodinger equation. This new method follows the same idea followed by the mean-field configuration interaction method already developed for molecular vibrations; i.e., groups of electronic degrees of freedom are contracted together in the mean field of the other degrees. If the same partition of electronic degrees of freedom is iterated, a self-consistent field method is obtained. Making coarser partitions (i.e., including more degrees in the same groups) and discarding the high energy states, the full configuration interaction limit can be approached. In contrast with the usual group function theory, no strong orthogonality condition is enforced. We have made use of a generalized version of the fundamental formula defining a Hopf algebra structure to derive Hamiltonian and overlap matrix element expressions which respect the group structure of the wave function as well as its fermionic symmetry. These expressions are amenable to a recursive computation.