Implementing Variational Quantum Algorithms to Compute Floquet Quasienergy Bands on Near-Term Hardware
Dane Kalorii
Abstract
Floquet engineering—the periodic driving of quantum systems—offers a powerful pathway to realize exotic nonequilibrium phases, including topological Floquet insulators and discrete time crystals. However, computing the quasienergy band structure of interacting, many-body Floquet systems remains a formidable challenge for classical computers due to the exponential growth of the Hilbert space and the difficulty of diagonalizing the non-Hermitian effective Floquet Hamiltonian. In this work, we introduce and experimentally demonstrate a variational quantum algorithm (VQA) tailored specifically for the near-term computation of Floquet quasienergy bands. Unlike standard variational eigensolvers that minimize static energy expectations, our approach employs a parameterized quantum circuit to prepare variational trial states and directly extracts the quasienergy spectrum via a Hadamard-test-based cost function that measures the complex overlap with the single-period Floquet unitary propagator. By scanning the conserved crystal momentum, we reconstruct the full quasienergy dispersion for both one- and two-dimensional lattice models. We benchmark the algorithm against exact diagonalization for system sizes up to N=8N=8 qubits and subsequently deploy it on superconducting transmon processors, systematically characterizing the impact of decoherence, gate infidelity, and readout errors on the resolved band structure. To mitigate hardware noise, we integrate zero-noise extrapolation and measurement-error mitigation, recovering quasienergy features with fidelities exceeding 90% for shallow circuit depths. Our results establish that variational phase-estimation protocols offer a scalable, hardware-efficient route to exploring driven quantum matter, providing a practical toolkit for simulating Floquet topological phenomena on current-generation quantum devices well before the advent of fault-tolerant quantum computation.