Extension and acceleration of relativistic density functional theory based on transformed density operator

We report an extension of relativistic density functional theory (RDFT) within one-component or two-component expressions that relies on a unitary-transformed density operator as well as a unitary-transformed Hamiltonian [Oyama et al., Chem. Phys. Lett. 680, 37 (2017)]. The transformed density operator is introduced to avoid the picture-change effect in the electron density, density gradient, kinetic energy density, and exchange-correlation potential. We confirmed that the implementation based on the spin-free infinite-order Douglas-Kroll-Hess method gives total, orbital, and excitation energies close to the reference values given by four-component RDFT calculations. To reduce the computational cost due to the transformed density operator, the local unitary transformation was also implemented. Numerical assessments revealed that the present scheme enabled the RDFT calculation of polyatomic systems with negligibly small picture-change effect.