Recent Development of Relativistic Molecular Theory

Today it is common knowledge that relativistic effects are important in the heavy-element chemistry. The continuing development of the relativistic molecular theory is opening up rows of the periodic table that are impossible to treat with the non-relativistic approach. The most straightforward way to treat relativistic effects on heavy-element systems is to use the four-component Dirac-Hartree-Fock approach and its electron-correlation methods based on the Dirac-Coulomb(-Breit) Hamiltonian. The Dirac-Hartree-Fock (DHF) or Dirac-Kohn-Sham (DKS) equation with the four-component spinors composed of the large- and small-components demands severe computational efforts to solve, and its applications to molecules including heavy elements have been limited to small- to medium-size systems. Recently, we have developed a very efficient algorithm for the four-component DHF and DKS approaches. As an alternative approach, several quasi-relativistic approximations have also been proposed instead of explicitly solving the four-component relativistic equation. We have developed the relativistic elimination of small components (RESC) and higher-order Douglas-Kroll (DK) Hamiltonians within the framework of the two-component quasi-relativistic approach. The developing four-component relativistic and approximate quasi-relativistic methods have been implemented into a program suite named REL4D.In this article, we will introduce the efficient relativistic molecular theories to treat heavy-atomic molecular systems accurately via the four-component relativistic and the two-component quasi-relativistic approaches. We will also show several chemical applications including heavy-element systems with our relativistic molecular approaches.     

Relativistic electronic structure theory

The theoretical and technical foundations are presented for the efficient relativistic electronic structure theories to treat heavy-atomic molecular systems. This review contains two surveys of four-component and two-component quasi-relativistic approaches. First, we review our highly efficient computational scheme for four-component relativistic ab initio molecular orbital (MO) methods over generally contracted spherical harmonic Gaussian-type spinors (GTSs). Illustrative calculations, which are performed with a new four-component relativistic ab initio molecular orbital program package REL4D, clearly show the efficiency of our computational scheme by the Dirac-Hartree-Fock (DHF) and Dirac-Hartree-Fock (DKS) methods. Next, in the two-component quasi-relativistic framework, two relativistic Hamiltonians, RESC and higher order Douglas-Kroll (DK) Hamiltonians, are introduced, and several illustrative calculations are shown. Numerical results for several systems show that good accuracy can be obtained with our third-order DK (DK3) Hamiltonian.     

Quasi-relativistic approximation in the SCF-MO-LCAO method

A quasi-relativistic approach to the MO-LCAO method is formulated taking into account the relativistic effects with an accuracy up to (v/c)² terms, the relativistic part of the electronic interaction in the Hamiltonian being neglected. In the framework of this approximation a set of SCF equations of the Roothaan form is derived; here only the relativistic analogue to the closed shell systems with one-determinant wave functions is considered. In so doing three types of relativistic corrections arise which are quite similar to those of the Pauli equation for one-electron atoms. The new matrix elements appearing due to these corrections can be reduced to some common integrals, which have to be calculated with relativistic radial atomic functions. The method allows a semi-empirical approach to the problem and does not require the Dirac four-component atomic functions (unknown in the most cases), thus making possible approximate quasi-relativistic electronic structure calculations of heavy-atom compounds.     

A direct relativistic four-component multiconfiguration self-consistent-field method for molecules

A new direct relativistic four-component Kramers-restricted multiconfiguration self-consistent-field (KR-MCSCF) code for molecules has been implemented. The program is based upon Kramers-paired spinors and a full implementation of the binary double groups (D(2h)(*) and subgroups). The underlying quaternion algebra for one-electron operators was extended to treat two-electron integrals and density matrices in an efficient and nonredundant way. The iterative procedure is direct with respect to both configurational and spinor variational parameters; this permits the use of large configuration expansions and many basis functions. The relativistic minimum-maximum principle is implemented in a second-order restricted-step optimization algorithm, which provides sharp and well-controlled convergence. This paper focuses on the necessary modifications of nonrelativistic MCSCF methodology to obtain a fully variational KR-MCSCF implementation. The general implementation also allows for the use of molecular integrals from a two-component relativistic Hamiltonian as, for example, the Douglas-Kroll-Hess variants. Several sample applications concern the determination of spectroscopic properties of heavy-element atoms and molecules, demonstrating the influence of spin-orbit coupling in MCSCF approaches to such systems and showing the potential of the new method.     

Accurate and efficient treatment of two-electron contributions in quasirelativistic high-order Douglas-Kroll density-functional calculations

Two-component quasirelativistic approaches are in principle capable of reproducing results from fully relativistic calculations based on the four-component Dirac equation (with fixed particle number). For one-electron systems, this also holds in practice, but in many-electron systems one has to transform the two-electron interaction, which is necessary because a picture change occurs when going from the Dirac equation to a two-component method. For one-electron properties, one can take full account of picture change in a manageable way, but for the electron interaction, this would spoil the computational advantages which are the main reason to perform quasirelativistic calculations. Exploiting those picture change effects are largest in the atomic cores, which in molecular applications do not differ too much from the cores of isolated neutral atoms, we propose an elegant, efficient, and accurate approximation to the two-electron picture change problem. The new approach, called the "model potential" approach because it makes use of atomic (four- and two-component) data to estimate picture change effects in molecules, shares with the nuclear-only approach that the Douglas-Kroll operator needs to be constructed only once (not in each self-consistent-field iteration) and that no time-consuming multicenter relativistic two-electron integrals need to be calculated. The new approach correctly describes the screening of both the nearest nucleus and distant nuclei, for the scalar-relativistic as well as the spin-orbit parts of the Hamiltonian. The approach is tested on atomic and molecular-orbital energies as well as spectroscopic constants of the lead dimer.     

Relativistic perturbation theory of chemical properties

Summary Double perturbation theory is developed for the case where relativity is one perturbation and the other perturbation describes a chemically interesting observable such as molecular structure, force constant or polarizability. Relativity is treated according to Rutkowski's nonsingular perturbation approach. Expressions for four-component and two-component wave-functions and for the Hartree-Fock approximation are given. The method is applied analytically to the relativistic corrections of the electric polarizability of the H atom, and algebraically to the potential curve of the H ₂ ⁺ molecule. Second and third order double perturbation interchange relations are numerically verified. In the present formalism, terms up to third order are needed to qualitatively understand the relativistic corrections of chemical observables.     

Relativistic effects on the Fukui function

The extent of relativistic effects on the Fukui function, which describes local reactivity trends within conceptual density functional theory (DFT), and frontier orbital densities has been analysed on the basis of three benchmark molecules containing the heavy elements: Au, Pb, and Bi. Various approximate relativistic approaches have been tested and compared with the four-component fully relativistic reference. Scalar relativistic effects, as described by the scalar zeroth-order regular approximation methodology and effective core potential calculations, already provide a large part of the relativistic corrections. Inclusion of spin–orbit coupling effects improves the results, especially for the heavy p-block compounds. We thus expect that future conceptual DFT-based reactivity studies on heavy-element molecules can rely on one of the approximate relativistic methodologies.     

Resolution of identity Dirac-Kohn-Sham method using the large component only: Calculations of g-tensor and hyperfine tensor

A new relativistic two-component density functional approach, based on the Dirac-Kohn-Sham method and an extensive use of the technique of resolution of identity (RI), has been developed and is termed the DKS2-RI method. It has been applied to relativistic calculations of g and hyperfine tensors of coinage-metal atoms and some mercury complexes. The DKS2-RI method solves the Dirac-Kohn-Sham equations in a two-component framework using explicitly a basis for the large component only, but it retains all contributions coming from the small component. The DKS2-RI results converge to those of the four-component Dirac-Kohn-Sham with an increasing basis set since the error associated with the use of RI will approach zero. The RI approximation provides a basis for a very efficient implementation by avoiding problems associated with complicated integrals otherwise arising from the elimination of the small component. The approach has been implemented in an unrestricted noncollinear two-component density functional framework. DKS2-RI is related to Dyall's [J. Chem. Phys. 106, 9618 (1997)] unnormalized elimination of the small component method (which was formulated at the Hartree-Fock level and applied to one-electron systems only), but it takes advantage of the local Kohn-Sham exchange-correlation operators (as, e.g., arising from local or gradient-corrected functionals). The DKS2-RI method provides an attractive alternative to existing approximate two-component methods with transformed Hamiltonians (such as Douglas-Kroll-Hess [Ann. Phys. 82, 89 (1974); Phys. Rev. A 33, 3742 (1986)] method, zero-order regular approximation, or related approaches) for relativistic calculations of the structure and properties of heavy-atom systems. In particular, no picture-change effects arise in the property calculations.     

Economical treatments of relativistic effects and electron correlation in WH6

The equilibrium geometries and relative stabilities of several structural isomers of tungsten hexahydride, WH₆, have been obtained at different levels of quantum chemical calculations. The performance of various strategies to (i) include electron correlation, viz. density functional theory based approaches, Møller/Plesset perturbation and coupled cluster theory, and to (ii) account for scalar relativistic effects, viz. various relativistic effective core potentials, first order perturbation theory, a quasi-relativistic treatment employing a Pauli Hamiltonian, and use of the Douglas/Kroll operator, are compared to the best theoretical data available. It is shown that relativistic and electron correlation effects are most important for the high-symmetry species, that these effects give rise to opposite trends in relative energies, and that overall the relativistic effects dominate. The most efficient way to incorporate relativistic effects appears to be via the use of relativistic effective core potentials, while the correlation energies are best taken account of using a conventional method such as CCSD(T). © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1604–1611, 1998     

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.     

Perturbation energy expansions based on two-component relativistic Hamiltonians

The approximate elimination of the small-component approach provides ansätze for the relativistic wave function. The assumed form of the small component of the wave function in combination with the Dirac equation define transformed but exact Dirac equations. The present derivation yields a family of two-component relativistic Hamiltonians which can be used as zeroth-order approximation to the Dirac equation. The operator difference between the Dirac and the two-component relativistic Hamiltonians can be used as a perturbation operator. The first-order perturbation energy corrections have been obtained from a direct perturbation theory scheme based on these two-component relativistic Hamiltonians. At the two-component relativistic level, the errors of the relativistic correction to the energies are proportional to ⁴ Z ⁴, whereas for the relativistic energy corrections including the first-order perturbation theory contributions, the errors are of the order of ⁶ Z ⁶–⁸ Z ⁸ depending on the zeroth-order Hamiltonian.Contribution to the Björn Roos Honorary Issue     

Electronic structure,spectra, and thermodynamical functions of molybdenum pentachloride

The SCF-X -SW method in non-relativistic and quasi-relativistic versions has been used to calculate the electronic structure, ionization potentials, energies and oscillator strengths of the optical transitions in MoCl₅. The electronic absorption spectrum of the gaseous MoCl₅ has been measured. The interpretation of the photoelectron and optical spectra of MoCl₅ is given. Spinpolarization effects and relativistic corrections are discussed. The thermodynamical functions of MoCl₅ (gas) are calculated.     

A closed-shell coupled-cluster treatment of the Breit-Pauli first-order relativistic energy correction

First-order relativistic corrections to the energy of closed-shell molecular systems are calculated, using all terms in the two-component Breit-Pauli Hamiltonian. In particular, we present the first implementation of the two-electron Breit orbit-orbit integrals, thus completing the first-order relativistic corrections within the two-component Pauli approximation. Calculations of these corrections are presented for a series of small and light molecules, at the Hartree-Fock and coupled-cluster levels of theory. Comparisons with four-component Dirac-Coulomb-Breit calculations demonstrate that the full Breit-Pauli energy corrections represent an accurate approximation to a fully relativistic treatment of such systems. The Breit interaction is dominated by the spin-spin interaction, the orbit-orbit interaction contributing only about 10% to the total two-electron relativistic correction in molecules consisting of light atoms. However, the relative importance of the orbit-orbit interaction increases with increasing nuclear charge, contributing more than 20% in H(2)S.     

Biomarkers of Aspergillus spores: Strain typing and protein identification

We applied both matrix-assisted laser desorption/ionization time of flight (MALDI-TOF) mass spectrometric and 1D sodium dodecylsulfate polyacrylamide gel electrophoretic (1D-PAGE) approaches for direct analysis of intact fungal spores of twenty four Aspergillus species. In parallel, we optimized various protocols for protein extraction from Aspergillus spores using acidic conditions, step organic gradient and variable sonication treatment. The MALDI-TOF mass spectra obtained from optimally prepared samples provided a reproducible fingerprint demonstrating the capability of the MALDI-TOF approach to type and characterize different fungal strains within the Aspergillus genus. Mass spectra of intact fungal spores provided signals mostly below 20 kDa. The minimum material amount represented 0.3 μg (10,000 spores). Proteins with higher molecular weight were detected by 1D-PAGE. Eleven proteins were identified from three selected strains in the range 5–25 kDa by the proteomic approach. Hemolysin and hydrophobin have the highest relevance in host–pathogen interactions.     

Mechanism of formation and monomolecular decomposition of aci-nitromethanes: a quantum-chemical study

A quantum-chemical study of the reactions of formation of aci-nitromethane (aci-NM) and aci-dinitromethane (aci-DNM) and their decomposition with elimination of water was carried out. The methods employed were the ab initio RHF method with inclusion of electron correlation at the MP2 level of theory and the Dunning—Hay double zeta basis set augmented with polarization d-functions on heavy-element atoms, the DFT approach at the B3LYP level, and the semiempirical PM3 method. The formation of aci-NM and aci-DNM was found to be the limiting stage of the mechanism under study. For DNM, the barrier to reaction is substantially lower than for NM. The estimates of the heights of the barriers to formation found from density functional calculations at the B3LYP/6-311++G(df,p) level (258 kJ mol^–1 for aci-NM and 218.5 kJ mol^–1 for aci-DNM) are thought to be the most reliable.     

Systematic Sequences of Geometric Relativistic Basis Sets. I:s- and p-Block Elements up to Xe

In this work a scheme for constructing systematic sequences of relativistic SCF basis sets at a reasonable computational cost is presented and applied to atoms of the s- and p-block up to Xe. This scheme, which couples simplex optimization and the use of geometric series given by four-term polynomial expressions for the logarithm of the exponents, allows for the construction of basis sets that exhibit very regular patterns of convergence to the numerical reference values of atomic total energies, spinor energies and radial expectation values. This regularity, together with the broad range of basis set sizes presented, enables these sets to be used as building blocks for basis sets applicable in both routine and benchmark relativistic calculations on atomic and molecular systems.Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Density functional calculations of molecular parity-violating effects within the zeroth-order regular approximation

A (quasirelativistic) two-component density functional theory (DFT) approach to the computation of parity-violating energy differences between enantiomers is presented which is based on the zeroth-order regular approximation (ZORA). This approach is employed herein to compute parity-violating energy differences between several P and M conformations of dihydrogen dichalcogenides (H2X2 with X=O, S, Se, Te, Po), of which some compounds have recently been suggested as potential molecular candidates for the first experimental measurement of parity-violating effects in chiral molecules. The DFT ZORA results obtained in this work with "pure" density functionals are anticipated to deviate by well less than 1% from data that would be computed within related (relativistic) four-component Dirac-Kohn-Sham-Coulomb schemes. In our implementation of the ZORA slightly larger relative deviations are expected for hybrid functionals, depending on the amount of "exact" exchange. For B3LYP (20% exact exchange) differences are estimated to amount to at most 3% in hydrogen peroxide, 2% in disulfane, and 1% or less for the heavier homologs. Thus, the present two-component approach is expected to perform excellently when compared to four-component density functional schemes while being at the same time computationally more efficient. The ZORA approach will therefore be of particular interest for the prediction of parity-violating vibrational frequency shifts, for instance, in isotopomers of H(2)Se(2) and H(2)Te(2).