Nonlinear Spectroscopy as a Magnon Breakdown Diagnosis and its Efficient Simulation

TL;DR

This paper introduces an efficient computational method for simulating nonlinear spectroscopic responses in quantum magnets, enabling the identification of magnon breakdown and fractionalized excitations through 2DCS measurements.

Contribution

It develops a Lanczos-based approach for calculating second-order susceptibilities directly in frequency domain and applies it to extended Kitaev models, revealing distinct nonlinear responses.

Findings

Different nonlinear responses at various magnetic field strengths in Kitaev models.

Derivation of universal 2DCS response for conventional magnons within linear spin-wave theory.

Deviations from the universal form indicate non-magnon excitations, aiding experimental identification.

Abstract

Identifying quantum spin liquids, magnon breakdown, or fractionalized excitations in quantum magnets is an ongoing challenge due to the ambiguity of possible origins of excitation continua occurring in linear response probes. Recently, it was proposed that techniques measuring higher-order response, such as two-dimensional coherent spectroscopy (2DCS), could resolve such ambiguities. Numerically simulating nonlinear response functions can, however, be computationally very demanding. We present an efficient Lanczos-based method to compute second-order susceptibilities $\chi^{2}\omega_t,\omega_\tau)$ directly in the frequency domain. Applying this to extended Kitaev models describing $\alpha$-RuCl$_3$, we find qualitatively different nonlinear responses between intermediate magnetic field strengths and the high-field regime. To put these results into context, we derive the general 2DCS response of partially-polarized magnets within the linear spin-wave approximation, establishing that $\chi^2(\omega_t,\omega_\tau)$ is restricted to a distinct universal form if the excitations are conventional magnons. Deviations from this form, as predicted in our (Lanczos-based) simulations for $\alpha$-RuCl$_3$, can hence serve in 2DCS experiments as direct criteria to determine whether an observed excitation continuum is of conventional two-magnon type or of different nature.