Simultaneous Treatment of Electron Correlation and Spin-Orbit Coupling with State-Averaged Driven Similarity Renormalization Group Second-Order Perturbation Theory.

We report the state-averaged spin-orbit driven similarity renormalization group second-order perturbation theory (SA-SO-DSRG-PT2) for efficient computations of spin-orbit coupling (SOC) effects in multireference systems. This approach adopts a set of multiconfigurational reference states obtained with a spin-free relativistic Hamiltonian and introduces SOC during the perturbative treatment of dynamical correlation. Specifically, the zeroth-order Hamiltonian is defined from the spin-free Hamiltonian, while the remaining spin-free contributions and the spin-dependent Hamiltonian enter as first-order perturbations in the DSRG transformation. This formulation generates an SOC-containing SA-DSRG-PT2 effective Hamiltonian while largely preserving the computational cost of the conventional spin-free theory. Two spin-dependent Hamiltonians are examined: the first-order Douglas-Kroll-Hess spin-orbit Hamiltonian with mean-field approximated two-electron contributions and the spin-dependent part of the one-electron X2C Hamiltonian with the screened-nuclear spin-orbit approximation. Benchmarks on main-group atoms and diatomic molecules, transition-metal elements, trivalent lanthanide cations, and actinide dioxide cations show consistent accuracy across the periodic table for both spin-dependent Hamiltonians. Applications to Co(II) single-ion magnets, including systems described with up to 1790 basis functions, further highlight the promise of SA-SO-DSRG-PT2 for efficient SOC calculations in large open-shell molecules.