Time-reversal symmetry adaptation in relativistic density matrix renormalization group algorithm

In the nonrelativistic Schrödinger equation, the total spin S and spin projection M are good quantum numbers. In contrast, spin symmetry is lost in the presence of spin-dependent interactions, such as spin-orbit couplings in relativistic Hamiltonians. Therefore, the relativistic density matrix renormalization group algorithm (R-DMRG) only employing particle number symmetry is much more expensive than nonrelativistic DMRG. In addition, artificial breaking of Kramers degeneracy can happen in the treatment of systems with an odd number of electrons. To overcome these issues, we propose time-reversal symmetry adaptation for R-DMRG. Since the time-reversal operator is antiunitary, this cannot be simply achieved in the usual way. We introduce a time-reversal symmetry-adapted renormalized basis and present strategies to maintain the structure of basis functions during the sweep optimization. With time-reversal symmetry adaptation, only half of the renormalized operators are needed, and the computational costs of Hamiltonian-wavefunction multiplication and renormalization are reduced by half. The present construction of the time-reversal symmetry-adapted basis also directly applies to other tensor network states without loops.