arXiv:2007.06112 · Journal of Mathematical Physics
Computing Floquet Hamiltonians with symmetries
Terry A. Loring, Fredy Vides
Abstract
Unitary matrices arise in many ways in physics, in particular as a time evolution operator. For a periodically driven system, one frequently wishes to compute a Floquet Hamiltonian that should be a Hermitian operator H such that e−iTH = U(T), where U(T) is the time evolution operator at time corresponding to the period of the system. That is, we want H to be equal to −i times a matrix logarithm of U(T). If the system has a symmetry, such as time reversal symmetry, one can expect H to have a symmetry beyond being Hermitian. We discuss here practical numerical algorithms on computing matrix logarithms that have certain symmetries, which can be used to compute Floquet Hamiltonians that have appropriate symmetries. Along the way, we prove some results on how a symmetry in the Floquet operator U(T) can lead to a symmetry in a basis of Floquet eigenstates.