Response functions for electric field induced two-dimensional nonlinear spectroscopy in a Kitaev magnet

We study the dynamical response functions relevant for electric field induced two-dimensional (2D) coherent nonlinear optical spectroscopy in a Kitaev magnet at finite temperature. We show that these response functions are susceptible to both types of fractional quasiparticles of this quantum spin-liquid, i.e. fermions and flux visons. Focusing on the second order response, we find a strong antidiagonal feature in the 2D frequency plane, related to the galvanoelectric effect of the fractional fermions. Perpendicular to the antidiagonal, the width of this feature is set by quasiparticle relaxation rates beyond the bare Kitaev magnet, thereby providing access to single-particle characteristics within a multi-particle continuum. While the 2D spectrum of the response is set by the fermionic quasiparticles and displays Fermi blocking versus temperature, the emergent bond randomness which arises due to thermally populated visons strongly modifies the fermionic spectrum. Therefore also the presence of gauge excitations is manifest in the 2D nonlinear response as the temperature is increased beyond the flux proliferation crossover.