Relativistic two-component calculations of electronic g-tensors that include spin polarization

The first two-component relativistic density-functional approach for the calculation of electronic g-tensors is reported that includes spin polarization using noncollinear spin-density functionals. The method is based on the relativistic Douglas-Kroll-Hess Hamiltonian and has been implemented into the ReSpect program package. Using three self-consistent-field calculations with orthogonal orientations of total magnetization J, the full g-matrix may be obtained. In contrast to previous spin-restricted two-component treatments, results with the new approach agree excellently with spin-polarized one-component calculations for light-atom radicals. Additionally, unlike one-component approaches, the method also reproduces successfully the negative deltag(parallel)-values of heavy-atom 2sigma radicals and the negative deltag(perpendicular) components in cysteinyl. The new method removes effectively the dilemma existing up to now regarding the simultaneous inclusion of spin polarization and higher-order spin-orbit effects in g-tensor calculations. It is straightforwardly applicable to higher than doublet spin multiplicities and has been implemented with hybrid functionals.     

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 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.     

Electronic structure of molecules using one-component wave functions and relativistic effective core potentials

A one-component approach to molecular electronic structure is discussed that includes the dominant relativistic effects on valence electrons and yet allows the use of the traditional quantum-chemistry techniques. The approach starts with one-component Cowan–Griffin relativistic orbitals that successfully incorporate the effects of the mass-velocity and Darwin terms present in more complicated wave functions such as the Dirac–Hartree–Fock. The approach then constructs “relativistic” effective core potentials (RECPS ) from these orbitals, and uses these to bring the relativistic effects into the molecular electronic calculations. The use of effective one-electron spin-orbit operators in conjunction with these one-component wave functions to include the effects of spin-orbit coupling is discussed. Applications to molecular systems involving heavy atoms and comparisons with available spectroscopic data on molecular geometries and excitation energies are presented. Finally, a new approach to the construction of RECPS encompassing the Hamiltonian and shapeconsistent approach is presented together with a novel analysis of the long-range behavior of the RECPS .     

A fully relativistic method for calculation of nuclear magnetic shielding tensors with a restricted magnetically balanced basis in the framework of the matrix Dirac-Kohn-Sham equation

A new relativistic four-component density functional approach for calculations of NMR shielding tensors has been developed and implemented. It is founded on the matrix formulation of the Dirac-Kohn-Sham (DKS) method. Initially, unperturbed equations are solved with the use of a restricted kinetically balanced basis set for the small component. The second-order coupled perturbed DKS method is then based on the use of restricted magnetically balanced basis sets for the small component. Benchmark relativistic calculations have been carried out for the (1)H and heavy-atom nuclear shielding tensors of the HX series (X=F,Cl,Br,I), where spin-orbit effects are known to be very pronounced. The restricted magnetically balanced basis set allows us to avoid additional approximations and/or strong basis set dependence which arises in some related approaches. The method provides an attractive alternative to existing approximate two-component methods with transformed Hamiltonians for relativistic calculations of chemical shifts and spin-spin coupling constants of heavy-atom systems. In particular, no picture-change effects arise in property calculations.     

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.     

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.     

Solvent effects on heavy atom nuclear spin-spin coupling constants: a theoretical study of Hg-C and Pt-P couplings.

The computation of indirect nuclear spin-spin coupling constants, based on the relativistic two-component zeroth order regular approximate Hamiltonian, has been recently implemented by us into the Amsterdam Density Functional program. Applications of the code for the calculation of one-bond metal-ligand couplings of coordinatively unsaturated compounds containing (195)Pt and (199)Hg, including spin-orbit coupling or coordination effects by solvent molecules, show that relativistic density functional calculations are able to reproduce the experimental findings with good accuracy for the systems under investigation. Spin-orbit effects are rather small for these cases, while coordination of the heavy atoms by solvent molecules has a great impact on the calculated couplings. Experimental trends for different solvents are reproduced. An orbital-based analysis of the solvent effect is presented. The scalar relativistic increase of the coupling constants is of the same order of magnitude as the nonrelativistically obtained values, making a relativistic treatment essential for obtaining quantitatively correct results. Solvent effects can be of similar importance.     

A simplified relativistic time-dependent density-functional theory formalism for the calculations of excitation energies including spin-orbit coupling effect

In the present work we have proposed an approximate time-dependent density-functional theory (TDDFT) formalism to deal with the influence of spin-orbit coupling effect on the excitation energies for closed-shell systems. In this formalism scalar relativistic TDDFT calculations are first performed to determine the lowest single-group excited states and the spin-orbit coupling operator is applied to these single-group excited states to obtain the excitation energies with spin-orbit coupling effects included. The computational effort of the present method is much smaller than that of the two-component TDDFT formalism and this method can be applied to medium-size systems containing heavy elements. The compositions of the double-group excited states in terms of single-group singlet and triplet excited states are obtained automatically from the calculations. The calculated excitation energies based on the present formalism show that this formalism affords reasonable excitation energies for transitions not involving 5p and 6p orbitals. For transitions involving 5p orbitals, one can still obtain acceptable results for excitations with a small truncation error, while the formalism will fail for transitions involving 6p orbitals, especially 6p1/2 spinors.     

The calculation of excitation energies based on the relativistic two-component zeroth-order regular approximation and time-dependent density-functional with full use of symmetry

In the present work, we propose a relativistic time-dependent density-functional theory (TDDFT) based on the two-component zeroth-order regular approximation and a noncollinear exchange-correlation (XC) functional. This two-component TDDFT formalism has the correct nonrelativistic limit and affords the correct threefold degeneracy of triplet excitations. The relativistic TDDFT formalism is implemented into the AMSTERDAM DENSITY FUNCTIONAL program package for closed-shell systems with full use of double-group symmetry to reduce the computational effort and facilitate the assignments. The performance of the formalism is tested on some closed-shell atoms, ions, and a few diatomic molecules containing heavy elements. The results show that the fine structure of the excited states for most atoms and ions studied here can be accurately accounted for with a proper XC potential. For the excitation energies of the molecules studied here, the present formalism shows promise and the error encountered is comparable to that of nonrelativistic TDDFT calculations on light elements.     

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.     

First principle simulation of the temperature dependent magnetic circular dichroism of a trinuclear copper complex in the presence of zero field splitting

We present a test of a recently developed density functional theory (DFT) based methodology for the calculation of magnetic circular dichroism (MCD) spectra in the presence of zero-field splitting (ZFS). The absorption and MCD spectra of the trinuclear copper complex μ(3)O ([Cu(3)(L)(μ(3)-O)](4+)), which models the native intermediate produced in the catalytic cycle of the multicopper oxidases, have been simulated from first principle within the framework of adiabatic time dependent density functional theory. The effects of the ZFS of the quartet (4)A(2) ground state on the theoretical MCD spectrum of μ(3)O have been analyzed. The simulated spectra are consistent with the experimental ones. The theoretical assignments of the MCD spectra are based on direct simulation as well as a detailed analysis of the molecular orbitals in μ(3)O. Some of the assignments differ from those given in previous studies. The ZFS effects in the presence of a strong external magnetic field (7 T) prove negligible. The change of the sign of the ZFS changes systematically the intensity of the MCD bands of the z-polarized excitations. The effect of the ZFS on the x,y-polarized excitations is not uniform.     

Density functional calculations of NMR shielding tensors for paramagnetic systems with arbitrary spin multiplicity: validation on 3d metallocenes

The calculation of nuclear shieldings for paramagnetic molecules has been implemented in the ReSpect program, which allows the use of modern density functional methods with accurate treatments of spin-orbit effects for all relevant terms up to order Omicron(alpha4) in the fine structure constant. Compared to previous implementations, the methodology has been extended to compounds of arbitrary spin multiplicity. Effects of zero-field splittings in high-spin systems are approximately accounted for. Validation of the new implementation is carried out for the 13C and 1H NMR signal shifts of the 3d metallocenes 4VCp2, 3CrCp2, 2MnCp2, 6MnCp2, 2CoCp2, and 3NiCp2. Zero-field splitting effects on isotropic shifts tend to be small or negligible. Agreement with experimental isotropic shifts is already good with the BP86 gradient-corrected functional and is further improved by admixture of Hartree-Fock exchange in hybrid functionals. Decomposition of the shieldings confirms the dominant importance of the Fermi-contact shifts, but contributions from spin-orbit dependent terms are frequently also non-negligible. Agreement with 13C NMR shift tensors from solid-state experiments is of similar quality as for isotropic shifts.     

Ab initio and DFT studies of the spin-orbit and spin-spin contributions to the zero-field splitting tensors of triplet nitrenes with aryl scaffolds

Spin-orbit and spin-spin contributions to the zero-field splitting (ZFS) tensors (D tensors) of spin-triplet phenyl-, naphthyl-, and anthryl-nitrenes in their ground state are investigated by quantum chemical calculations, focusing on the effects of the ring size and substituted position of nitrene on the D tensor. A hybrid CASSCF/MRMP2 approach to the spin-orbit term of the D tensor (D(SO) tensor), which was recently proposed by us, has shown that the spin-orbit contribution to the entire D value, termed the ZFS parameter or fine-structure constant, is about 10% in all the arylnitrenes under study and less depends on the size and connectivity of the aryl groups. Order of the absolute values for D(SO) can be explained by the perturbation on the energy level and spatial distributions of π-SOMO through the orbital interaction between SOMO of the nitrene moiety and frontier orbitals of the aryl scaffolds. Spin-spin contribution to the D tensor (D(SS) tensor) has been calculated in terms of the McWeeny-Mizuno equation with the DFT/EPR-II spin densities. The D(SS) value calculated with the RO-B3LYP spin density agrees well with the D(Exptl) -D(SO) reference value in phenylnitrene, but agreement with the reference value gradually becomes worse as the D value decreases. Exchange-correlation functional dependence on the D(SS) tensor has been explored with standard 23 exchange-correlation functionals in both RO- and U-DFT methodologies, and the RO-HCTH/407 method gives the best agreement with the D(Exptl) -D(SO) reference value. Significant exchange-correlation functional dependence is observed in spin-delocalized systems such as 9-anthrylnitrene (6). By employing the hybrid CASSCF/MRMP2 approach and the McWeeny-Mizuno equation combined with the RO-HCTH/407/EPR-II//U-HCTH/407/6-31G* spin densities for D(SO) and D(SS), respectively, a quantitative agreement with the experiment is achieved with errors less than 10% in all the arylnitrenes under study. Guidelines to the putative approaches to D(SS) tensor calculations are given.     

Theoretical study of the electronic spectra of square-planar platinum (II) complexes based on the two-component relativistic time-dependent density-functional theory

In the present work the electronic spectra of [PtCl(4)](2-), [PtBr(4)](2-), and [Pt(CN)(4)](2-) are studied with a recently proposed relativistic time-dependent density-functional theory (TDDFT) based on the two-component zeroth-order regular approximation and a noncollinear exchange-correlation (XC) functional. The contribution to the double group excited states in terms of singlet and triplet single group excited states is estimated through the inner product of the transition density matrix obtained from two-component and scalar relativistic TDDFT calculations to better understand the double group excited states. Spin-orbital coupling effects are found to be very important in order to simulate the electronic spectra of these complexes. The results show that the two-component TDDFT formalism can afford excitation energies with high accuracy for the transition-metal systems studied here when use is made of a proper XC potential.     

Localization scheme for relativistic spinors

A new method to determine localized complex-valued one-electron functions in the occupied space is presented. The approach allows the calculation of localized orbitals regardless of their structure and of the entries in the spinor coefficient matrix, i.e., one-, two-, and four-component Kramers-restricted or unrestricted one-electron functions with real or complex expansion coefficients. The method is applicable to localization schemes that maximize (or minimize) a functional of the occupied spinors and that use a localization operator for which a matrix representation is available. The approach relies on the approximate joint diagonalization (AJD) of several Hermitian (symmetric) matrices which is utilized in electronic signal processing. The use of AJD in this approach has the advantage that it allows a reformulation of the localization criterion on an iterative 2 × 2 pair rotating basis in an analytical closed form which has not yet been described in the literature for multi-component (complex-valued) spinors. For the one-component case, the approach delivers the same Foster-Boys or Pipek-Mezey localized orbitals that one obtains from standard quantum chemical software, whereas in the multi-component case complex-valued spinors satisfying the selected localization criterion are obtained. These localized spinors allow the formulation of local correlation methods in a multi-component relativistic framework, which was not yet available. As an example, several heavy and super-heavy element systems are calculated using a Kramers-restricted self-consistent field and relativistic two-component pseudopotentials in order to investigate the effect of spin-orbit coupling on localization.     

Optical and magnetic properties of disubstituted benzophenones with lowest 3nπ* states

The phosphorescence and optically detected magnetic resonance (ODMR) spectra of 4,4′-dimethoxy-, 4,4′-dimethyl-,3,3′-dibromo- and 2,2′-diazaben have been measured. For 4,4′-dimethoxybenzophenone two X-traps have been found. The zero-field splitting (ZFS) parameters and the phosphorescence dec indicate that the lowest triplet states have mainly nπ* character, The effects of the substituents on the ZFS parameters cannot be explained via spin coupling with only nearby excited states.     

Approximate Relativistic calculations within the one-center approximation for the series CH4 TO PbH4

Using one-component relativistic theory combined with an Xαβ local exchange, we have carried out calculations for the series CH₄ to PbH₄ within the spherically symmetric one-center expansion approximation. Bond lenghts, breathing force constant, ionization potentials and spin-orbit splitting of the filled valence p shell thus obtained, are found to be in excellent agreement with Dirac-Hartre-Fock results.