Unambiguous identification of fractionalized excitations in quantum spin liquids has been a long-standing issue in correlated topological phases. Conventional spectroscopic probes, such as the dynamical spin structure factor, can only detect composites of fractionalized excitations, leading to a broad continuum in energy. Lacking a clear signature in conventional probes has been the biggest obstacle in the field. In this work, we theoretically investigate what kinds of distinctive signatures of fractionalized excitations can be probed in two-dimensional nonlinear spectroscopy by considering the exactly solvable Kitaev spin liquids. We demonstrate the existence of a number of salient features of the Majorana fermions and fluxes in two-dimensional nonlinear spectroscopy, which provide crucial information about such excitations.
Long-range entanglement in quantum spin liquids (QSLs) leads to novel low-energy excitations with fractionalized quantum numbers and (in two dimensions) statistics. Experimental detection and manipulation of these excitations present a challenge particularly in view of diverse candidate magnets. A promising probe of fractionalization is their coupling to phonons. Here, we present Raman scattering results for the
S=1/2 honeycomb iridate
Cu2IrO3, a candidate Kitaev QSL with fractionalized Majorana fermions and Ising flux excitations. We observe anomalous low-temperature frequency shift and linewidth broadening of the Raman intensities in addition to a broad magnetic continuum, both of which, as we derive, are naturally attributed to the phonon decaying into itinerant Majoranas. The dynamic Raman susceptibility marks a crossover from the QSL to a thermal paramagnet at $\ensuremath{\sim}120$ K. The phonon anomalies below this temperature demonstrate a strong phonon-Majorana coupling. These results provide evidence of spin fractionalization in
Cu2IrO3.
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
Z2 flux excitations and show that 2DCS thus directly probes the inverse participation ratio of Majorana wavefunctions. By providing experimentally observable features of spin liquid states in the 2D spectrum, our work can guide future 2DCS experiments on Kitaev magnets.