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

It has recently been shown that the SOiCI approach [J. Phys.: Condens. Matter 34 (2022) 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 [J. Chem. Phys. 138 (2013) 104113] but with the iCISCF algorithm [J. Chem. Theory Comput. 17 (2021) 7545] for large active spaces. The resulting SOiCISCF employs both double group and time reversal symmetries for computational efficiency and assignment of target states. The fine structures of $p$-block elements are taken as showcases to reveal the efficacy of SOiCISCF.