The relativistic four-component coupled cluster method for molecules: spectroscopic constants of SnH4

State Interaction for Relativistic Four-Component Methods: Choose the Right Zeroth-Order Hamiltonian for Late-Row Elements

We present several schemes based on the spin-separation of the Dirac-Coulomb-Breit Hamiltonian for the perturbative treatment of relativistic four-component Hamiltonians within the state interaction framework. While state interaction approaches traditionally use zeroth-order scalar-relativistic states, we develop augmented zeroth-order Hamiltonians with increasing accuracy and investigate convergence to the variational limit as a function of the choice of zeroth-order Hamiltonian. The state interaction schemes developed in this work are benchmarked using ground-state fine-structure splitting of late-row atoms and diatomic hyrides. Although the scalar-relativistic zeroth-order Hamiltonian exhibits significant errors in ground-state fine-structure splitting, the predictive accuracy can be improved by augmenting the zeroth-order Hamiltonian with one- and two-electron vector-relativistic operators (e.g., spin-orbit, spin-spin, orbit-orbit). This work lays the theoretical foundation for the development of low-scaling, high-accuracy perturbative relativistic methods suitable for late-row elements.

Relativistic Coupled Cluster with Completely Renormalized and Perturbative Triples Corrections

We have implemented noniterative triples corrections to the energy from coupled-cluster with single and double excitations (CCSD) within the 1-electron exact two-component (1eX2C) relativistic framework. The effectiveness of both the CCSD(T) and the completely renormalized (CR) CC(2,3) approaches are demonstrated by performing all-electron computations of the potential energy curves and spectroscopic constants of copper, silver, and gold dimers in their ground electronic states. Spin-orbit coupling effects captured via the 1eX2C framework are shown to be crucial for recovering the correct shape of the potential energy curves, and the correlation effects due to triples in these systems change the dissociation energies by about 0.1-0.2 eV or about 4-7%. We also demonstrate that relativistic effects and basis set size and contraction scheme are significantly more important in Au2 than in Ag2 or Cu2.

Dirac-Coulomb-Breit Molecular Mean-Field Exact-Two-Component Relativistic Equation-of-Motion Coupled-Cluster Theory

We present a relativistic equation-of-motion coupled-cluster with single and double excitation formalism within the exact two-component framework (X2C-EOM-CCSD), where both scalar relativistic effects and spin-orbit coupling are variationally included at the reference level. Three different molecular mean-field treatments of relativistic corrections, including the one-electron, Dirac-Coulomb, and Dirac-Coulomb-Breit Hamiltonian, are considered in this work. Benchmark calculations include atomic excitations and fine-structure splittings arising from spin-orbit coupling. Comparison with experimental values and relativistic time-dependent density functional theory is also carried out. The computation of the oscillator strength using the relativistic X2C-EOM-CCSD approach allows for studies of spin-orbit-driven processes, such as the spontaneous phosphorescence lifetime.

Comparison of Variational and Perturbative Spin-Orbit Coupling within Two-Component CASSCF

The modeling of spin-orbit coupling (SOC) remains a challenge in computational chemistry due to the high computational cost. With the rising popularity of spin-driven processes and f-block metals in chemistry and materials science, it is incumbent on the community to develop accurate multiconfigurational SOC methods that scale to large systems and understand the limits of different treatments of SOC. Herein, we introduce an implementation of perturbative SOC in scalar-relativistic two-component CASSCF (srX2C-CASSCF-SO). Perspectives on the limitations and accuracy of srX2C-CASSCF-SO are presented via benchmark calculations.

Negative Ion Resonance States: The Fock-Space Coupled-Cluster Way

The negative ion resonance states, which are electron-molecule metastable compound states, play the most important role in free-electron controlled molecular reactions and low-energy free-electron-induced DNA damage. Their electronic structure is often only poorly described but crucial to an understanding of their reaction dynamics. One of the most important challenges to current electronic structure theory is the computation of negative ion resonance states. As a major step forward, coupled-cluster theories, which are well-known for their ability to produce the best approximate bound state electronic eigen solutions, are upgraded to offer the most accurate and effective approximations for negative ion resonance states. The existing Fock-space coupled-cluster (FSCC) and the equation-of-motion coupled-cluster (EOM-CC) approaches that compute bound states are redesigned for the direct and simultaneous determination of both the kinetic energy of the free electron at which the electron-molecule compound states are resonantly formed and the corresponding autodetachment decay rate of the electron from the metastable compound state. This Feature Article reviews the computation of negative ion resonances using the FSCC approach and, in passing, provides the highlights of the equivalent EOM-CC approach.

RELATIVISTIC AND CORRELATION EFFECTS IN ATOMS

This review article deals with some case studies of relativistic and correlation effects in atomic systems. After a brief introduction to relativistic many-electron theory, a number of applications ranging from correlation energy to parity non-conservation in atoms are considered. There is a special emphasis on relativistic coupled-cluster theory as most of the results presented here are based on it.

Is E112 a relatively inert element? Benchmark relativistic correlation study of spectroscopic constants in E112H and its cation

We report the first results of relativistic correlation calculation of the spectroscopic properties for the ground state of E112H and its cation in which spin-orbit interaction is taken into account non-perturbatively. Studying the properties of E112 (eka-Hg) is required for chemical identification of its long-lived isotope, (283)112. It is shown that appropriate accounting for spin-orbit effects leads to dramatic impact on the properties of E112H whereas they are not so important for E112H(+). The calculated equilibrium distance, R(e) (calc)=1.662 Angstrom, in E112H is notably smaller than R(e) (expt)=(1.738+/-0.003) Angstrom and R(e) (calc)=1.738 Angstrom in HgH, whereas the dissociation energy, D(e) (calc)=0.42 eV, in E112H is close to D(e) (expt)=0.46 eV and D(e) (calc)=0.41 eV in HgH. These data are quite different from R(e) (NH)=1.829 Angstrom and D(e) (NH)=0.06 eV obtained for E112H within the scalar-relativistic Douglas-Kroll approximation. Our results indicate that E112 should not be expected to behave like a noble gas in contrast to the results by other authors.