Relativistic two-component projection-based quantum embedding for open-shell systems

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.

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.

On the accuracy of orbital based multi-level approaches for closed-shell transition metal chemistry

In this work, we investigate the accuracy of the local molecular orbital molecular orbital (LMOMO) scheme and projection-based wave function-in-density functional theory (WF-in-DFT) embedding for the prediction of reaction energies and barriers of typical reactions involving transition metals. To analyze the dependence of the accuracy on the system partitioning, we apply a manual orbital selection for LMOMO as well as the so-called direct orbital selection (DOS) for both approaches. We benchmark these methods on 30 closed shell reactions involving 16 different transition metals. This allows us to devise guidelines for the manual selection as well as settings for the DOS that provide accurate results within an error of 2 kcal mol^-1 compared to local coupled cluster. To reach this accuracy, on average 55% of the occupied orbitals have to be correlated with coupled cluster for the current test set. Furthermore, we find that LMOMO gives more reliable relative energies for small embedded regions than WF-in-DFT embedding.

Relativistic Kramers-Unrestricted Exact-Two-Component Density Matrix Renormalization Group

In this work we develop a variational relativistic density matrix renormalization group (DMRG) approach within the exact-two-component (X2C) framework (X2C-DMRG), using spinor orbitals optimized with the two-component relativistic complete active space self-consistent field. We investigate fine-structure splittings of p- (Ga, In, Tl) and d-block (Sc, Y, La) atoms and excitation energies of monohydride molecules (GeH, SnH, and TlH) with X2C-DMRG calculations using an all-electron relativistic Hamiltonian in a Kramers-unrestricted basis. We find that X2C-DMRG yields accurate ²P and ²D splittings compared to multireference configuration interaction with singles and doubles (MRCISD). We also investigated the degree of symmetry breaking in the atomic multiplets and convergence of electron correlation in the total energies. Symmetry breaking can be large in some cases (∼30 meV); however, increasing the number of renormalized block states m for the DMRG optimization recovers the symmetry breaking by several orders of magnitude. Encouragingly, we find the convergence of electron correlation to be close to MRCISDTQ5 quality. Relativistic X2C-DMRG approaches are important for cases where spin-orbit coupling is significant and the underlying reference wave function requires a large determinantal space. We are able to obtain quantitatively correct fine-structure splittings for systems up to 10¹⁹ number of determinants with traditional CI approaches, which are currently unfeasible to converge for the field.

Relativistic nonorthogonal configuration interaction: application to L2,3-edge X-ray spectroscopy

In this article, we develop a relativistic exact-two-component nonorthogonal configuration interaction (X2C-NOCI) for computing L-edge X-ray spectra. This article to our knowledge is the first time NOCI has been used for relativistic wave functions. A set of molecular complexes, including SF₆, SiCl₄ and [FeCl₆]^3-, are used to demonstrate the accuracy and computational scaling of the X2C-NOCI method. Our results suggest that X2C-NOCI is able to satisfactorily capture the main features of the L_2,3-edge X-ray absorption spectra. Excitations from the core require a large amount of orbital relaxation to yield reasonable energies and X2C-NOCI allows us to treat orbital optimization explicitly. However, the cost of computing the nonorthogonal coupling is higher than in conventional CI. Here, we propose an improved integral screening using overlap-scaled density combined with a continuous measure of the generalized Slater-Condon rules that allows us to estimate if an element is zero before attempting a two-electron integral contraction.

Exact-Two-Component Multiconfiguration Pair-Density Functional Theory

Molecules containing late-row elements exhibit large relativistic effects. To account for both relativistic effects and electron correlation in a computationally inexpensive way, we derived a formulation of multiconfiguration pair-density functional theory with the relativistic exact-two-component Hamiltonian (X2C-MC-PDFT). In this new method, relativistic effects are included during variational optimization of a reference wave function by exact-two-component complete active-space self-consistent-field (X2C-CASSCF) theory, followed by an energy evaluation using pair-density functional theory. Benchmark studies of excited-state and ground-state fine-structure splitting of atomic species show that X2C-MC-PDFT can significantly improve the X2C-CASSCF results by introducing additional state-specific electron correlation.

Relativistic Effects in Modeling the Ligand K-Edge X-ray Absorption Near-Edge Structure of Uranium Complexes

Accurate modeling of the complex electronic structure of actinide complexes requires full inclusion of relativistic effects. In this study, we examine the effect of explicit inclusion of spin-orbit coupling (SOC) versus scalar relativistic effects on the predicted spectra for heavy-element complexes. In this study, we employ a relativistic two-component Hamiltonian in the X2C form with all of the electrons in the system being considered explicitly to compare and contrast with previous studies that included the relativistic effects by means of relativistic effective core potentials (RECPs). A few uranium complexes are chosen as model systems. Comparison of the computed Cl K-edge X-ray absorption spectra with experimental data shows significantly improved agreement when a variational relativistic treatment of SOC is performed. In particular, we note the importance of SOC terms to obtain not only correct transition energies but also correct intensities for these heavy-element complexes because of the redistribution of ligand bonding character among the valence MOs. While RECPs generally agree well with all-electron scalar relativistic calculations, there are some differences in the predicted spectra of open-shell systems. These methods are still suitable for broad application to analyze the qualitative nature of transitions in X-ray absorption spectra, but caution is recommended for quantitative analysis, as SOC can be non-negligible for both open- and closed-shell heavy-element systems.

Direct orbital selection within the domain-based local pair natural orbital coupled-cluster method

Domain-based local pair natural orbital coupled cluster (DLPNO-CC) has become increasingly popular to calculate relative energies (e.g., reaction energies and reaction barriers). It can be applied within a multi-level DLPNO-CC-in-DLPNO-CC ansatz to reduce the computational cost and focus the available computational resources on a specific subset of the occupied orbitals. We demonstrate how this multi-level DLPNO-CC ansatz can be combined with our direct orbital selection (DOS) approach [M. Bensberg and J. Neugebauer, J. Chem. Phys. 150, 214106 (2019)] to automatically select orbital sets for any multi-level calculation. We find that the parameters for the DOS procedure can be chosen conservatively such that they are transferable between reactions. The resulting automatic multi-level DLPNO-CC method requires no user input and is extremely robust and accurate. The computational cost is easily reduced by a factor of 3 without sacrificing accuracy. We demonstrate the accuracy of the method for a total of 61 reactions containing up to 174 atoms and use it to predict the relative stability of conformers of a Ru-based catalyst.

Exact-two-component block-localized wave function: A simple scheme for the automatic computation of relativistic ΔSCF

Block-localized wave function is a useful method for optimizing constrained determinants. In this article, we extend the generalized block-localized wave function technique to a relativistic two-component framework. Optimization of excited state determinants for two-component wave functions presents a unique challenge because the excited state manifold is often quite dense with degenerate states. Furthermore, we test the degree to which certain symmetries result naturally from the ΔSCF optimization such as time-reversal symmetry and symmetry with respect to the total angular momentum operator on a series of atomic systems. Variational optimizations may often break the symmetry in order to lower the overall energy, just as unrestricted Hartree-Fock breaks spin symmetry. Overall, we demonstrate that time-reversal symmetry is roughly maintained when using Hartree-Fock, but less so when using Kohn-Sham density functional theory. Additionally, maintaining total angular momentum symmetry appears to be system dependent and not guaranteed. Finally, we were able to trace the breaking of total angular momentum symmetry to the relaxation of core electrons.

Density functional theory based embedding approaches for transition-metal complexes

Transition metal species are commonly discussed by considering the metal atom embedded in a ligand environment. This apparently makes them interesting targets for modern embedding strategies based on Kohn-Sham density functional theory (DFT), which aim at modelling accurate predictions for large systems by combining different quantum chemical methods. In this perspective, we will focus on subsystem density functional theory and projection-based embedding. We review the developments in the field for transition metal species, demonstrate benefits, drawbacks and analyse error sources of the different strategies using the example of chromium hexacarbonyle, before giving a perspective where the field is currently heading.