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.     

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.     

Spin ground state and magnetic properties of cobalt(II): relativistic DFT calculations guided by EPR measurements of bis(2,4-acetylacetonate)cobalt(II)-based complexes

The spin ground state of the core ion and structure of the bis(2,4-acetylacetonate)cobalt(II) model complex and its synthetic aqua and ethanol derivatives, Co(acac)(2)L(n), (L = EtOH, H(2)O), were examined by means of density functional theory (DFT) calculations supported by electron paramagnetic resonance (EPR) measurements. Geometry optimizations were carried out for low-spin (doublet) and high-spin (quartet) states. For the Co(acac)(2) complex two possible conformations, a square-planar and a tetrahedral one, were taken into account. For all structures relative energies were calculated with both "pure" and hybrid functionals. The calculated data were complemented with the results of the EPR investigations carried out at liquid helium temperature, allowing for definite assignment of the high-spin state for the Co(acac)(2)(EtOH)(2) complex. However, because of the unresolved spectral features, only effective g-values could be assessed, whereas the zero-field splitting parameters (ZFS) were calculated by means of the spin-orbit mean field (SOMF) relativistic DFT method for which direct spin-spin (SS) and spin-orbit coupling (SOC) contributions were quantified.     

First-principles calculations of zero-field splitting parameters

In this work, an implementation of an approach to calculate the zero-field splitting (ZFS) constants in the framework of ab initio methods such as complete active space self-consistent field, multireference configuration interaction, or spectroscopy oriented configuration interaction is reported. The spin-orbit coupling (SOC) contribution to ZFSs is computed using an accurate multicenter mean-field approximation for the Breit-Pauli Hamiltonian. The SOC parts of ZFS constants are obtained directly after diagonalization of the SOC operator in the basis of a preselected number of roots of the spin-free Hamiltonian. This corresponds to an infinite order treatment of the SOC in terms of perturbation theory. The spin-spin (SS) part is presently estimated in a mean-field fashion and appears to yield results close to the more complete treatments available in the literature. Test calculations for the first- and second-row atoms as well as first-row transition metal atoms and a set of diatomic molecules show accurate results for the SOC part of ZFSs. SS contributions have been found to be relatively small but not negligible (exceeding 1 cm(-1) for oxygen molecule). At least for the systems studied in this work, it is demonstrated that the presented method provides much more accurate estimations for the SOC part of ZFS constants than the emerging density functional theory approaches.     

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.     

Conformational analysis of singlet-triplet state mixing in Paternò-Büchi diradicals

Conformational dependence of spin-orbit coupling (SOC) in flexible Paternò-Büchi (PB) diradicals has been studied with high-level ab initio methods using both (i) one-electron spin-orbit Hamiltonian with parametrized (effective) nuclear charges in conjunction with a state-averaged MCSCF wave function as implemented by Robb in Gaussian 98 and (ii) complete one- and two-electron SOC with a fully optimized MCSCF triplet wave function and frozen core singlet as implemented by Furlani in the GAMESS computational package. The ab initio results revealed two distinct areas of elevated SOC values, one corresponding to the region whereby a cisoid conformation in the C-C-O-C fragment brings the two odd-electron orbitals closer to each other, and the other area corresponding to the partially eclipsed conformation lacking direct overlap between the spin centers. In this second region the 1,4-electronic communication is mediated by the oxygen's 2p-lone pair, which is suitably oriented to play the role of a "relay-antenna". The other critical factor affecting the rate of intersystem crossing (ISC)--singlet-triplet energy separation--was computed utilizing a multireference CASSCF-MP2 method to include dynamic correlation effects. The largest singlet-triplet energy gap, approximately 2 kcal/mol, was found for a gauche conformer (also a minimum SOC conformation). Rotation about the central C-O bond either toward the fully eclipsed (0 degrees ) or the partially eclipsed (120 degrees ) conformations decreases the singlet-triplet gap while increasing the value of the SOC matrix element. These computational findings support the Griesbeck model for stereochemistry of triplet PB reactions and provide a rigorous basis for predicting the probability of ISC in diradicals separated by a partially conjugated spacer.     

TD-DFT Accuracy in Determining Excited-state Structures and Fluorescence Spectra of Firefly Emitter

We analyzed the excited-state structures and emission spectra of firefly emitter, the anionic keto form of firefly oxyluciferin(keto-1), determined by the time dependent-density functional theory(TD-DFT) approach. The analysis is based on a direct comparison with the highly correlated CASSCF(MS-CASPT2) ab initio approach. 49 DFT functionals were considered and applied to the study. Among the tested functionals, mPW3PBE, B3PW91 and B3P86 give the best performance for ground-state geometry, absorption spectrum, excited-state geometry and emission spectrum.     

Detailed ab initio first-principles study of the magnetic anisotropy in a family of trigonal pyramidal iron(II) pyrrolide complexes

A theoretical, computational, and conceptual framework for the interpretation and prediction of the magnetic anisotropy of transition metal complexes with orbitally degenerate or orbitally nearly degenerate ground states is explored. The treatment is based on complete active space self-consistent field (CASSCF) wave functions in conjunction with N-electron valence perturbation theory (NEVPT2) and quasidegenerate perturbation theory (QDPT) for treatment of magnetic field- and spin-dependent relativistic effects. The methodology is applied to a series of Fe(II) complexes in ligand fields of almost trigonal pyramidal symmetry as provided by several variants of the tris-pyrrolylmethyl amine ligand (tpa). These systems have recently attracted much attention as mononuclear single-molecule magnet (SMM) complexes. This study aims to establish how the ligand field can be fine tuned in order to maximize the magnetic anisotropy barrier. In trigonal ligand fields high-spin Fe(II) complexes adopt an orbitally degenerate (5)E ground state with strong in-state spin-orbit coupling (SOC). We study the competing effects of SOC and the (5)E?ε multimode Jahn-Teller effect as a function of the peripheral substituents on the tpa ligand. These subtle distortions were found to have a significant effect on the magnetic anisotropy. Using a rigorous treatment of all spin multiplets arising from the triplet and quintet states in the d(6) configuration the parameters of the effective spin-Hamiltonian (SH) approach were predicted from first principles. Being based on a nonperturbative approach we investigate under which conditions the SH approach is valid and what terms need to be retained. It is demonstrated that already tiny geometric distortions observed in the crystal structures of four structurally and magnetically well-documented systems, reported recently, i.e., [Fe(tpa(R))](-) (R = tert-butyl, Tbu (1), mesityl, Mes (2), phenyl, Ph (3), and 2,6-difluorophenyl, Dfp (4), are enough to lead to five lowest and thermally accessible spin sublevels described sufficiently well by S = 2 SH provided that it is extended with one fourth order anisotropy term. Using this most elementary parametrization that is consistent with the actual physics, the reported magnetization data for the target systems were reinterpreted and found to be in good agreement with the ab initio results. The multiplet energies from the ab initio calculations have been fitted with remarkable consistency using a ligand field (angular overlap) model (ab initio ligand field, AILFT). This allows for determination of bonding parameters and quantitatively demonstrates the correlation between increasingly negative D values and changes in the σ-bond strength induced by the peripheral ligands. In fact, the sigma-bonding capacity (and hence the Lewis basicity) of the ligand decreases along the series 1 > 2 > 3 > 4.     

Comparing ab initio density-functional and wave function theories: the impact of correlation on the electronic density and the role of the correlation potential

The framework of ab initio density-functional theory (DFT) has been introduced as a way to provide a seamless connection between the Kohn-Sham (KS) formulation of DFT and wave-function based ab initio approaches [R. J. Bartlett, I. Grabowski, S. Hirata, and S. Ivanov, J. Chem. Phys. 122, 034104 (2005)]. Recently, an analysis of the impact of dynamical correlation effects on the density of the neon atom was presented [K. Jankowski, K. Nowakowski, I. Grabowski, and J. Wasilewski, J. Chem. Phys. 130, 164102 (2009)], contrasting the behaviour for a variety of standard density functionals with that of ab initio approaches based on second-order M?ller-Plesset (MP2) and coupled cluster theories at the singles-doubles (CCSD) and singles-doubles perturbative triples [CCSD(T)] levels. In the present work, we consider ab initio density functionals based on second-order many-body perturbation theory and coupled cluster perturbation theory in a similar manner, for a range of small atomic and molecular systems. For comparison, we also consider results obtained from MP2, CCSD, and CCSD(T) calculations. In addition to this density based analysis, we determine the KS correlation potentials corresponding to these densities and compare them with those obtained for a range of ab initio density functionals via the optimized effective potential method. The correlation energies, densities, and potentials calculated using ab initio DFT display a similar systematic behaviour to those derived from electronic densities calculated using ab initio wave function theories. In contrast, typical explicit density functionals for the correlation energy, such as VWN5 and LYP, do not show behaviour consistent with this picture of dynamical correlation, although they may provide some degree of correction for already erroneous explicitly density-dependent exchange-only functionals. The results presented here using orbital dependent ab initio density functionals show that they provide a treatment of exchange and correlation contributions within the KS framework that is more consistent with traditional ab initio wave function based methods.     

Ab initio study of the ground and lower-lying excited electronic states of NiX2 and FeX2 (X=F,Cl, Br,I) molecules

The ground and lower-lying excited electronic states of FeX₂ and NiX₂ (X=F, Cl, Br, I) molecules are systematically investigated by ab initio method at the complete active space self-consistent field (CASSCF) and multiconfigurational quasi-degenerate second-order perturbation (MCQDPT2) levels of theory. It is concluded that the dynamic electron correlation has to be taken into account in the prediction of the properties for such kind of molecules. The equilibrium bond lengths r_e(M–X), force constants and harmonic vibrational frequencies are calculated for the ground and lower-lying excited electronic states. The spin-orbit coupling (SOC) effects are analysed.     

Ab initio calculation of anisotropic magnetic properties of complexes. I. Unique definition of pseudospin Hamiltonians and their derivation

A methodology for the rigorous nonperturbative derivation of magnetic pseudospin Hamiltonians of mononuclear complexes and fragments based on ab initio calculations of their electronic structure is described. It is supposed that the spin-orbit coupling and other relativistic effects are already taken fully into account at the stage of quantum chemistry calculations of complexes. The methodology is based on the establishment of the correspondence between the ab initio wave functions of the chosen manifold of multielectronic states and the pseudospin eigenfunctions, which allows to define the pseudospin Hamiltonians in the unique way. Working expressions are derived for the pseudospin Zeeman and zero-field splitting Hamiltonian corresponding to arbitrary pseudospins. The proposed calculation methodology, already implemented in the SINGLE_ANISO module of the MOLCAS-7.6 quantum chemistry package, is applied for a first-principles evaluation of pseudospin Hamiltonians of several complexes exhibiting weak, moderate, and very strong spin-orbit coupling effects.     

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.     

Matrix isolation ESR spectroscopy and quantum chemical calculations on 5-methylhexa-1,2,4-triene-1,3-diyl, a highly delocalized triplet "hybrid" carbene

The ESR spectrum of 5-methylhexa-1,2,4-triene-1,3-diyl (1) was recorded in an argon matrix at 15 K. The derived zero-field splitting (ZFS) parameters (D = 0.5054 ± 0.0006 cm(-1) and E = 0.0045 ± 0.0002 cm(-1)) fall between those determined previously for propargylene (2) and vinylcarbene (3). DFT and ab initio (CAS and MRCI) quantum-chemical calculations of the ZFS parameters of 1, 2, and 3 were performed. These calculations indicate that multireference methods are needed to successfully predict ZFS parameters of delocalized carbenes/biradicals such as 1-3. The calculated singly occupied MOs and spin density distributions show that the spin is more delocalized in 1 than in 2 and 3, indicating that 1 is a "hybrid" of the constituent ethynyl- and vinylcarbenes, 2 and 3, respectively. The dominant contribution to the D-value in 1 and 2 is found to result from spin-spin interactions on the C atoms of the propylidene moiety, which is strongly affected by spin polarization. Accurate values for the D-parameter are also predicted for other types of delocalized triplet carbenes such as HC(5)H and HCCN.     

Kinetic study of the reaction H2O2+H→H2O+OH by ab initio and density functional theory calculations

Theoretical investigations on the kinetics of the elementary reaction H₂O₂+H→H₂O+OH were performed using the transition state theory (TST). Ab initio (MP2//CASSCF) and density functional theory (B3LYP) methods were used with large basis set to predict the kinetic parameters; the classical barrier height and the pre-exponential factor. The ZPE and BSSE corrected value of the classical barrier height was predicted to be 4.1 kcal mol⁻¹ for MP2//CASSCF and 4.3 kcal mol⁻¹ for B3LYP calculations. The experimental value fitted from Arrhenius expressions ranges from 3.6 to 3.9 kcal mol⁻¹. Thermal rate constants of the title reaction, based on the ab initio and DFT calculations, was evaluated for temperature ranging from 200 to 2500 K assuming a direct reaction mechanism. The modeled ab initio-TST and DFT–TST rate constants calculated without tunneling were found to be in reasonable agreement with the observed ones indicating that the contribution of the tunneling effect to the reaction was predicted to be unimportant at ambient temperature.     

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.     

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.     

Computational methods for analysis of an unsaturated carbocycle: heptafulvene

With respect to geometric optimizations, harmonic vibrational frequencies and single point conformational energies, various computational methods [HF, MP2, CCSD(T), BD(T), CASSCF, CASPT2, and DFT] were evaluated for their suitability to describe the heptafulvene system. We found that a significant number of basis sets lead to wrong predictions of folded minima, when ab initio methods including dynamic electron correlation are used. Possible explanations for these inconsistencies, such as wave function instabilities, near linear dependences of the basis sets and inadequate inclusion of polarization functions in the basis set, are discussed. Such concerns are likewise important for other classes of π-conjugated compounds, such that the results are expected to be of interest not only for heptafulvenes.     

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.