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

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.

Exact Two-Component TDDFT with Simple Two-Electron Picture-Change Corrections: X-ray Absorption Spectra Near L- and M-Edges of Four-Component Quality at Two-Component Cost

X-ray absorption spectroscopy (XAS) has gained popularity in recent years as it probes matter with high spatial and elemental sensitivities. However, the theoretical modeling of XAS is a challenging task since XAS spectra feature a fine structure due to scalar (SC) and spin-orbit (SO) relativistic effects, in particular near L and M absorption edges. While full four-component (4c) calculations of XAS are nowadays feasible, there is still interest in developing approximate relativistic methods that enable XAS calculations at the two-component (2c) level while maintaining the accuracy of the parent 4c approach. In this article we present theoretical and numerical insights into two simple yet accurate 2c approaches based on an (extended) atomic mean-field exact two-component Hamiltonian framework, (e)amfX2C, for the calculation of XAS using linear eigenvalue and damped response time-dependent density functional theory (TDDFT). In contrast to the commonly used one-electron X2C (1eX2C) Hamiltonian, both amfX2C and eamfX2C account for the SC and SO two-electron and exchange-correlation picture-change (PC) effects that arise from the X2C transformation. As we demonstrate on L- and M-edge XAS spectra of transition metal and actinide compounds, the absence of PC corrections in the 1eX2C approximation results in a substantial overestimation of SO splittings, whereas (e)amfX2C Hamiltonians reproduce all essential spectral features such as shape, position, and SO splitting of the 4c references in excellent agreement, while offering significant computational savings. Therefore, the (e)amfX2C PC correction models presented here constitute reliable relativistic 2c quantum-chemical approaches for modeling XAS.

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.