Probing Majorana Wave Functions in Kitaev Honeycomb Spin Liquids with Second-Order Two-Dimensional Spectroscopy

Two-dimensional coherent terahertz spectroscopy (2DCS) emerges as a valuable tool to probe the nature, couplings, and lifetimes of excitations in quantum materials. It thus promises to identify unique signatures of spin liquid states in quantum magnets by directly probing properties of their exotic fractionalized excitations. Here, we calculate the second-order 2DCS of the Kitaev honeycomb model and demonstrate that distinct spin liquid fingerprints appear already in this lowest-order nonlinear response χ_{yzx}^{(2)}(ω_{1},ω_{2}) when using crossed light polarizations. We further relate the off-diagonal 2DCS peaks to the localized nature of the matter Majorana excitations trapped by Z_{2} flux excitations and show that 2DCS thus directly probes the inverse participation ratio of Majorana wave functions. By providing experimentally observable features of spin liquid states in the 2D spectrum, our Letter can guide future 2DCS experiments on Kitaev magnets.