Computing Floquet quasienergies in finite and extended systems: Role of electromagnetic and quantum-geometric gauges

We present an approach to compute the Floquet quasienergy spectrum of time-periodic systems. The method allows to characterize the light-matter interaction in finite and extended structures by carefully addressing the resolution of the position operator. In periodic systems we discuss the role of the quantum-geometric gauge freedom of Bloch states and employ a Wannier-based scheme to compute the required matrix elements. As a consequence, the method is accurate and applicable to a broad range of systems, from atoms and molecules to cold atomic gases and materials described by density functional theory, as well as model systems. We demonstrate the applicability of the approach by studying two cases: a particle trapped in a one-dimensional box and the semiconducting material ${\mathrm{BC}}_{2}\mathrm{N}$. We employ the first example to provide a numerical proof of the invariance of the Floquet quasienergy spectrum with respect to the choice of electromagnetic gauge. The analysis of ${\mathrm{BC}}_{2}\mathrm{N}$ then serves to illustrate the physical effects described by the quasienergies, such as multiphoton resonances, and their expected range of occurrence in real materials in terms of external electric field and frequency of the drive pulse.