Relativistic electronic structure theory

Accurate Spin-Orbit Coupling by Relativistic Mixed-Reference Spin-Flip-TDDFT

Relativistic mixed-reference spin-flip (MRSF)-TDDFT is developed considering the spin-orbit coupling (SOC) within the mean-field approximation. The resulting SOC-MRSF faithfully reproduces the experiments with very high accuracy, which is also consistent with the values by four-component (4c) relativistic CASSCF and 4c-CASPT2 in the spin-orbit-energy splitting calculations of the C, Si, and Ge atoms. Even for the fifth-row element Sn, the SOC-MRSF yielded accurate splittings (∼ 3 % error). In the SOC calculations of the molecular 4-thiothymine with a third-row element, SOC-MRSF values are in excellent agreement with those of the SO-GMC-QDPT2 level, regardless of geometries and exchange-correlation functionals. The same SOC-MRSF predicted the anticipated chance of S₁ (nπ*) → T₁ (ππ*) intersystem crossing, even in thymine with only second-row elements. With its accuracy and practicality, thus, SOC-MRSF is a promising electronic structure protocol in challenging situations such as nonadiabatic molecular dynamics (NAMD) incorporating both internal conversions and intersystem crossings in large systems.

Four-Component Relativistic State-Specific Multireference Perturbation Theory with a Simplified Treatment of Static Correlation

The relativistic multireference (MR) perturbative approach is one of the most successful tools for the description of computationally demanding molecular systems of heavy elements. We present here the ground state dissociation energy surfaces, equilibrium bond lengths, harmonic frequencies, and dissociation energies of Ag₂, Cu₂, Au₂, and I₂ computed using the four-component (4c) relativistic spinors based state-specific MR perturbation theory (SSMRPT) with improved virtual orbital complete active space configuration interaction (IVO-CASCI) functions. The IVO-CASCI method is a simple, robust, useful and lower cost alternative to the complete active space self-consistent field approach for treating quasidegenerate situations. The redeeming features of the resulting method, termed as 4c-IVO-SSMRPT, lies in (i) manifestly size-extensivity, (ii) exemption from intruder problems, (iii) the freedom of convenient multipartitionings of the Hamiltonian, (iv) flexibility of the relaxed and unrelaxed descriptions of the reference coefficients, and (v) manageable cost/accuracy ratio. The present method delivers accurate descriptions of dissociation processes of heavy element systems. Close agreement with reference values has been found for the calculated molecular constants indicating that our 4c-IVOSSMRPT provides a robust and economic protocol for determining the structural properties for the ground state of heavy element molecules with eloquent MR character as it treats correlation and relativity on equal footing.

Implementation of Analytical Energy Gradient of Spin-Dependent General Hartree-Fock Method Based on the Infinite-Order Douglas-Kroll-Hess Relativistic Hamiltonian with Local Unitary Transformation

An analytical energy gradient for the spin-dependent general Hartree-Fock method based on the infinite-order Douglas-Kroll-Hess (IODKH) method was developed. To treat realistic systems, the local unitary transformation (LUT) scheme was employed both in energy and energy gradient calculations. The present energy gradient method was numerically assessed to investigate the accuracy in several diatomic molecules containing fifth- and sixth-period elements and to examine the efficiency in one-, two-, and three-dimensional silver clusters. To arrive at a practical calculation, we also determined the geometrical parameters of fac-tris(2-phenylpyridine)iridium and investigated the efficiency. The numerical results confirmed that the present method describes a highly accurate relativistic effect with high efficiency. The present method can be a powerful scheme for determining geometries of large molecules, including heavy-element atoms.

Potential energy curves for the ground and low-lying excited states of CuAg

The ground and low-lying excited states of heteronuclear diatomic CuAg are examined by multi-reference configuration interaction (MRCI) method. Relativistic effects were treated and probed in two steps. Scalar terms were considered using the spin-free DKH Hamiltonian as a priori and spin-orbit coupling was calculated perturbatively via the spin-orbit terms of the Breit-Pauli Hamiltonian based on MRCI wavefunctions. Potential energy curves of the spin-free states and their corresponding Ω components correlating with the separated atom limits (2)S(Cu) + (2)S(Ag) and (2)D(Cu) + (2)S(Ag) are obtained. The results are in fine agreement with the experimental measurements and tentative conclusions for the ion-pair B0(+) state are confirmed by our theoretical calculations. Illustrative results are presented to reveal the relative importance and magnitude of the scalar and spin-orbit effects on the spectroscopic properties of this molecule. Time dependent density functional theory calculations, using the LDA, BLYP, B3LYP, and SAOP functionals have been carried out for CuAg and the accuracy of TD-DFT has been compared with ab initio results.

RECENT ADVANCES IN ELECTRONIC STRUCTURE THEORY

Accurate quantum computational chemistry has evolved dramatically. The size of molecular systems, which can be studied accurately using molecular theory is increasing very rapidly. Theoretical chemistry has opened up a world of new possibilities. It can treat real systems with predictable accuracy. Computational chemistry is becoming an integral part of chemistry research. Theory can now make very significant contribution to chemistry.

This review will focus on our recent developments in the theoretical and computational methodology for the study of molecular structure and molecular interactions. We are aiming at developing accurate molecular theory on systems containing hundreds of atoms. We continue our research in the following three directions: (i) development of new ab initio theory, particularly multireference-based perturbation theory, (ii) development of exchange and correlation functionals in density functional theory, and (iii) development of molecular theory including relativistic effects.

We have enjoyed good progress in each of the above areas. We are very excited about our discoveries of new theory and new algorithms and would like to share this enthusiasm with readers.

Calculation of zero-field splitting parameters: Comparison of a two-component noncolinear spin-density-functional method and a one-component perturbational approach

Two different sets of approaches for the density-functional calculation of the spin-orbit contributions to zero-field splitting (ZFS) parameters of high-spin systems have been implemented within the same quantum chemistry code ReSpect and have been validated and compared for a series of model systems. The first approach includes spin-orbit coupling variationally in a two-component calculation, using either an all-electron Douglas-Kroll-Hess ansatz or two-component relativistic pseudopotentials. The ZFS parameters are computed directly from energy differences between different relativistic states. Additionally, an approximate second-order perturbation theory approach has been implemented, based on nonrelativistic or scalar relativistic wave functions. For a series of group 16 triplet diatomics and for the octet GdH3 molecules, two-component density functional calculations underestimate the zero-field splitting D systematically by a factor of 2. This may be rationalized readily by the incomplete description of states with absolute value MJ < J by a single-determinantal wave function built from two-component spinors. In the case of two 3d transition metal complexes and for GdH3, the results depend furthermore sensitively on exchange-correlation functional. Results of the alternative one-component approach agree strikingly with the two-component data for systems with small spin-orbit effects and start to deviate from them only for heavier systems with large spin-orbit effects. These results have fundamental implications for the achievable accuracy of one-component density-functional approaches used widely to compute ZFS parameters in the field of molecular magnetism. Possible refinements of both one-and two-component approaches are discussed.