iVI‐TD‐DFT: An iterative vector interaction method for exterior/interior roots of TD‐DFT

Unified construction of relativistic Hamiltonians

It is shown that the four-component (4C), quasi-four-component (Q4C), and exact two-component (X2C) relativistic Hartree-Fock equations can be implemented in a unified manner by making use of the atomic nature of the small components of molecular 4-spinors. A model density matrix approximation can first be invoked for the small-component charge/current density functions, which gives rise to a static, pre-molecular mean field to be combined with the one-electron term. As a result, only the nonrelativistic-like two-electron term of the 4C/Q4C/X2C Fock matrix needs to be updated during the iterations. A "one-center small-component" approximation can then be invoked in the evaluation of relativistic integrals, that is, all atom-centered small-component basis functions are regarded as extremely localized near the position of the atom to which they belong such that they have vanishing overlaps with all small- or large-component functions centered at other nuclei. Under these approximations, the 4C, Q4C, and X2C mean-field and many-electron Hamiltonians share precisely the same structure and accuracy. Beyond these is the effective quantum electrodynamics Hamiltonian that can be constructed in the same way. Such approximations lead to errors that are orders of magnitude smaller than other sources of errors (e.g., truncation errors in the one- and many-particle bases as well as uncertainties of experimental measurements) and are, hence, safe to use for whatever purposes. The quaternion forms of the 4C, Q4C, and X2C equations are also presented in the most general way, based on which the corresponding Kramers-restricted open-shell variants are formulated for "high-spin" open-shell systems.

SOiCISCF: Combining SOiCI and iCISCF for Variational Treatment of Spin-Orbit Coupling

It has recently been shown that the SOiCI approach [Zhang, N.; J. Phys.: Condens. Matter 2022, 34, 224007], in conjunction with the spin-separated exact two-component relativistic Hamiltonian, can provide very accurate fine structures of systems containing heavy elements by treating electron correlation and spin-orbit coupling (SOC) on an equal footing. Nonetheless, orbital relaxations/polarizations induced by SOC are not yet fully accounted for due to the use of scalar relativistic orbitals. This issue can be resolved by further optimizing the still real-valued orbitals self-consistently in the presence of SOC, as done in the spin-orbit coupled CASSCF approach [Ganyushin, D.; et al. J. Chem. Phys. 2013, 138, 104113] but with the iCISCF algorithm [Guo, Y.; J. Chem. Theory Comput. 2021, 17, 7545-7561] for large active spaces. The resulting SOiCISCF employs both double group and time reversal symmetries for computational efficiency and the assignment of target states. The fine structures of p-block elements are taken as showcases to reveal the efficacy of SOiCISCF.

Self-Adaptive Real-Time Time-Dependent Density Functional Theory for X-ray Absorptions

Real-time time-dependent density functional theory (RT-TDDFT) can in principle access the whole absorption spectrum of a many-electron system exposed to a narrow pulse. However, this requires an accurate and efficient propagator for the numerical integration of the time-dependent Kohn-Sham equation. While a low-order time propagator is already sufficient for the low-lying valence absorption spectra, it is no longer the case for the X-ray absorption spectra (XAS) of systems composed even only of light elements, for which the use of a high-order propagator is indispensable. It is then crucial to choose a largest possible time step and a shortest possible simulation time, so as to minimize the computational cost. To this end, we propose here a robust AutoPST approach to determine automatically (Auto) the propagator (P), step (S), and time (T) for relativistic RT-TDDFT simulations of XAS.

Computational approaches for XANES, VtC-XES, and RIXS using linear-response time-dependent density functional theory based methods

The emergence of state-of-the-art X-ray light sources has paved the way for novel spectroscopies that take advantage of their atomic specificity to shed light on fundamental physical, chemical, and biological processes both in the static and time domains. The success of these experiments hinges on the ability to interpret and predict core-level spectra, which has opened avenues for theory to play a key role. Over the last two decades, linear-response time-dependent density functional theory (LR-TDDFT), despite various theoretical challenges, has become a computationally attractive and versatile framework to study excited-state spectra including X-ray spectroscopies. In this context, we focus our discussion on LR-TDDFT approaches for the computation of X-ray Near-Edge Structure (XANES), Valence-to-Core X-ray Emission (VtC-XES), and Resonant Inelastic X-ray Scattering (RIXS) spectroscopies in molecular systems with an emphasis on Gaussian basis set implementations. We illustrate these approaches with applications and provide a brief outlook of possible new directions.

SOiCI and iCISO: combining iterative configuration interaction with spin-orbit coupling in two ways

The near-exact iCIPT2 approach for strongly correlated systems of electrons, which stems from the combination of iterative configuration interaction (iCI, an exact solver of full CI) with configuration selection for static correlation and second-order perturbation theory (PT2) for dynamic correlation, is extended to the relativistic domain. In the spirit of spin separation, relativistic effects are treated in two steps: scalar relativity is treated by the infinite-order, spin-free part of the exact two-component (X2C) relativistic Hamiltonian, whereas spin-orbit coupling (SOC) is treated by the first-order, Douglas-Kroll-Hess-like SOC operator derived from the same X2C Hamiltonian. Two possible combinations of iCIPT2 with SOC are considered, i.e., SOiCI and iCISO. The former treats SOC and electron correlation on an equal footing, whereas the latter treats SOC in the spirit of state interaction, by constructing and diagonalizing an effective spin-orbit Hamiltonian matrix in a small number of correlated scalar states. Both double group and time reversal symmetries are incorporated to simplify the computation. Pilot applications reveal that SOiCI is very accurate for the spin-orbit splitting (SOS) of heavy atoms, whereas the computationally very cheap iCISO can safely be applied to the SOS of light atoms and even of systems containing heavy atoms when SOC is largely quenched by ligand fields.

Accurate X-ray Absorption Spectra near L- and M-Edges from Relativistic Four-Component Damped Response Time-Dependent Density Functional Theory

The simulation of X-ray absorption spectra requires both scalar and spin-orbit (SO) relativistic effects to be taken into account, particularly near L- and M-edges where the SO splitting of core p and d orbitals dominates. Four-component Dirac-Coulomb Hamiltonian-based linear damped response time-dependent density functional theory (4c-DR-TDDFT) calculates spectra directly for a selected frequency region while including the relativistic effects variationally, making the method well suited for X-ray applications. In this work, we show that accurate X-ray absorption spectra near L2,3- and M4,5-edges of closed-shell transition metal and actinide compounds with different central atoms, ligands, and oxidation states can be obtained by means of 4c-DR-TDDFT. While the main absorption lines do not change noticeably with the basis set and geometry, the exchange-correlation functional has a strong influence with hybrid functionals performing the best. The energy shift compared to the experiment is shown to depend linearly on the amount of Hartee-Fock exchange with the optimal value being 60% for spectral regions above 1000 eV, providing relative errors below 0.2% and 2% for edge energies and SO splittings, respectively. Finally, the methodology calibrated in this work is used to reproduce the experimental L2,3-edge X-ray absorption spectra of [RuCl2(DMSO)2(Im)2] and [WCl4(PMePh2)2], and resolve the broad bands into separated lines, allowing an interpretation based on ligand field theory and double point groups. These results support 4c-DR-TDDFT as a reliable method for calculating and analyzing X-ray absorption spectra of chemically interesting systems, advance the accuracy of state-of-the art relativistic DFT approaches, and provide a reference for benchmarking more approximate techniques.

iCISCF: An Iterative Configuration Interaction-Based Multiconfigurational Self-Consistent Field Theory for Large Active Spaces

An iterative configuration interaction (iCI)-based multiconfigurational self-consistent field (SCF) theory, iCISCF, is proposed to handle systems that require large active spaces. The success of iCISCF stems from three ingredients: (1) efficient selection of individual configuration state functions spanning the active space while maintaining full spin symmetry; (2) the use of Jacobi rotation for optimization of the active orbitals in conjunction with a quasi-Newton algorithm for the core/active-virtual and core-active orbital rotations; (3) a second-order perturbative treatment of the residual space left over by the selection procedure (i.e., iCISCF(2)). Several examples that go beyond the capability of CASSCF are taken as showcases to reveal the efficacy of iCISCF and iCISCF(2), facilitated by iCAS for imposed automatic selection and localization of active orbitals.

The Static-Dynamic-Static Family of Methods for Strongly Correlated Electrons: Methodology and Benchmarking

A series of methods (SDSCI, SDSPT2, iCI, iCIPT2, iCISCF(2), iVI, and iCAS) is introduced to accurately describe strongly correlated systems of electrons. Born from the (restricted) static-dynamic-static (SDS) framework for designing many-electron wave functions, SDSCI is a minimal multireference (MR) configuration interaction (CI) approach that constructs and diagonalizes a [Formula: see text] matrix for [Formula: see text] states, regardless of the numbers of orbitals and electrons to be correlated. If the full molecular Hamiltonian H in the QHQ block (which describes couplings between functions of the first-order interaction space Q) of the SDSCI CI matrix is replaced with a zeroth-order Hamiltonian [Formula: see text] before the diagonalization is taken, we obtain SDSPT2, a CI-like second-order perturbation theory (PT2). Unlike most variants of MRPT2, SDSPT2 treats single and multiple states in the same way and is particularly advantageous in the presence of near degeneracy. On the other hand, if the SDSCI procedure is repeated until convergence, we will have iterative CI (iCI), which can converge quickly from the above to the exact solutions (full CI) even when starting with a poor guess. When further combined with the selection of important configurations followed by a PT2 treatment of dynamic correlation, iCI becomes iCIPT2, which is a near-exact theory for medium-sized systems. The microiterations of iCI for relaxing the coefficients of contracted many-electron functions can be generalized to an iterative vector interaction (iVI) approach for finding exterior or interior roots of a given matrix, in which the dimension of the search subspace is fixed by either the number of target roots or the user-specified energy window. Naturally, iCIPT2 can be employed as the active space solver of the complete active space (CAS) self-consistent field, leading to iCISCF(2), which can further be combined with iCAS for automated selection of active orbitals and assurance of the same CAS for all states and all geometries. The methods are calibrated by taking the Thiel set of benchmark systems as examples. Results for the corresponding cations, a new set of benchmark systems, are also reported.

Dynamic-then-Static Approach for Core Excitations of Open-Shell Molecules

Delta self-consistent-field methods are widely used in studies of electronically excited states. However, the nonaufbau determinants are generally spin-contaminated. Here, we describe a general approach for spin-coupling interactions of open-shell molecules, making use of multistate density functional theory (MSDFT). In particular, the effective exchange integrals that determine spin coupling are obtained by enforcing the multiplet degeneracy of the S+1 state in the M_S = S manifold. Consequently, they are consistent with the energy of the high-spin state that is adequately treated by Kohn-Sham density functional theory (DFT) and, thereby, free of double counting of correlation. The method was applied to core excitations of open-shell molecules and compared with those by spin-adapted time-dependent DFT. An excellent agreement with experiment was found employing the BLYP functional and aug-cc-pCVQZ basis set. Overall, MSDFT provides an effective combination of the strengths of DFT and wave function theory to achieve efficiency and accuracy.

iCAS: Imposed Automatic Selection and Localization of Complete Active Spaces

It is shown that in the spirit of "from fragments to molecule" for localizing molecular orbitals [J. Chem. Theory Comput. 2011, 7, 3643], a prechosen set of occupied/virtual valence/core atomic/fragmental orbitals can be transformed to an equivalent set of localized occupied/virtual pre-localized molecular orbitals (pre-LMO), which can then be taken as probes to select the same number of maximally matching localized occupied/virtual Hartree-Fock (HF) or restricted open-shell HF (ROHF) molecular orbitals as the initial local orbitals spanning the desired complete active space (CAS). In each cycle of the self-consistent field (SCF) calculation, the CASSCF orbitals can be localized by means of the noniterative "top-down least-change" algorithm for localizing ROHF orbitals [J. Chem. Phys. 2017, 146, 104104] such that the maximum matching between the orbitals of two adjacent iterations can readily be monitored, leading finally to converged localized CASSCF orbitals that overlap most the guess orbitals. Such an approach is to be dubbed as "imposed CASSCF" (iCASSCF or simply iCAS in short) for good reasons: (1) it has been assumed that only those electronic states that have largest projections onto the active space defined by the prechosen atomic/fragmental orbitals are to be targeted. This is certainly an imposed constraint but has wide applications in organic and transition metal chemistry where valence (or core) atomic/fragmental orbitals can readily be identified. (2) The selection of both initial and optimized local active orbitals is imposed from the very beginning by the pre-LMOs (which span the same space as the prechosen atomic/fragmental orbitals). Apart from the (imposed) automation and localization, iCAS has two additional merits: (a) the guess orbitals are guaranteed to be the same for all geometries, for the pre-LMOs do not change in character with geometry and (b) the use of localized orbitals facilitates the SCF convergence, particularly for large active spaces. Both organic molecules and transition-metal complexes are taken as showcases to reveal the efficacy of iCAS.

Relativistic EOM-CCSD for Core-Excited and Core-Ionized State Energies Based on the Four-Component Dirac-Coulomb(-Gaunt) Hamiltonian

We report an implementation of the core-valence separation approach to the four-component relativistic Hamiltonian-based equation-of-motion coupled-cluster with singles and doubles theory (CVS-EOM-CCSD) for the calculation of relativistic core-ionization potentials and core-excitation energies. With this implementation, which is capable of exploiting double group symmetry, we investigate the effects of the different CVS-EOM-CCSD variants and the use of different Hamiltonians based on the exact two-component (X2C) framework on the energies of different core-ionized and -excited states in halogen- (CH₃I, HX, and X^-, X = Cl-At) and xenon-containing (Xe, XeF₂) species. Our results show that the X2C molecular mean-field approach [Sikkema, J.; J. Chem. Phys. 2009, 131, 124116], based on four-component Dirac-Coulomb mean-field calculations (²DC^M), is capable of providing core excitations and ionization energies that are nearly indistinguishable from the reference four-component energies for up to and including fifth-row elements. We observe that two-electron integrals over the small-component basis sets lead to non-negligible contributions to core binding energies for the K and L edges for atoms such as iodine or astatine and that the approach based on Dirac-Coulomb-Gaunt mean-field calculations (²DCG^M) are significantly more accurate than X2C calculations for which screened two-electron spin-orbit interactions are included via atomic mean-field integrals.

Further Development of iCIPT2 for Strongly Correlated Electrons

The efficiency of the recently proposed iCIPT2 [iterative configuration interaction (iCI) with selection and second-order perturbation theory (PT2); J. Chem. Theory Comput. 2020, 16, 2296] for strongly correlated electrons is further enhanced (by up to 20×) by using (1) a new ranking criterion for configuration selection, (2) a new particle-hole algorithm for Hamiltonian construction over randomly selected configuration state functions (CSF), and (3) a new data structure for the quick sorting of the variational and first-order interaction spaces. Meanwhile, the memory requirement is also significantly reduced. As a result, this improved implementation of iCIPT2 can handle 1 order of magnitude more CSFs than the previous version, as revealed by taking the chromium dimer and an iron-sulfur cluster, [Fe₂S₂(SCH₃)]₄^2-, as examples.