We present and benchmark a quantum computing approach to calculate the two-dimensional coherent spectrum (2DCS) of high-spin models. Our approach is based on simulating their real-time dynamics in the presence of several magnetic field pulses, which are spaced in time. We utilize the adaptive variational quantum dynamics simulation (AVQDS) algorithm for the study due to its compact circuits, which enables simulations over sufficiently long times to achieve the required resolution in frequency space. Specifically, we consider an antiferromagnetic quantum spin model that incorporates Dzyaloshinskii-Moriya interactions and single-ion anisotropy. The obtained 2DCS spectra exhibit distinct peaks at multiples of the magnon frequency, arising from transitions between different eigenstates of the unperturbed Hamiltonian. By comparing the one-dimensional coherent spectrum with 2DCS, we demonstrate that 2DCS provides a higher resolution of the energy spectrum. We further investigate how the quantum resources scale with the magnitude of the spin using two different binary encodings of the high-spin operators: the standard binary encoding and the Gray code. At low magnetic fields both encodings require comparable quantum resources, but at larger field strengths the Gray code is advantageous. Numerical simulations for spin models with increasing number of sites indicate a polynomial system-size scaling for quantum resources. Lastly, we compare the numerical 2DCS with experimental results on a rare-earth orthoferrite system. The observed strength of the magnonic high-harmonic generation signals in the 2DCS of the quantum high-spin model aligns well with the experimental data, showing significant improvement over the corresponding mean-field results.
We present and benchmark a quantum computing approach to calculate the two-dimensional coherent spectrum (2DCS) of high-spin models. Our approach is based on simulating their real-time dynamics in the presence of several magnetic field pulses, which are spaced in time. We utilize the adaptive variational quantum dynamics simulation (AVQDS) algorithm for the study due to its compact circuits, which enables simulations over sufficiently long times to achieve the required resolution in frequency space. Specifically, we consider an antiferromagnetic quantum spin model that incorporates Dzyaloshinskii-Moriya interactions and single-ion anisotropy. The obtained 2DCS spectra exhibit distinct peaks at multiples of the magnon frequency, arising from transitions between different eigenstates of the unperturbed Hamiltonian. By comparing the one-dimensional coherent spectrum with 2DCS, we demonstrate that 2DCS provides a higher resolution of the energy spectrum. We further investigate how the quantum resources scale with the magnitude of the spin using two different binary encodings of the high-spin operators: the standard binary encoding and the Gray code. At low magnetic fields both encodings require comparable quantum resources, but at larger field strengths the Gray code is advantageous. Numerical simulations for spin models with increasing number of sites indicate a polynomial system-size scaling for quantum resources. Lastly, we compare the numerical 2DCS with experimental results on a rare-earth orthoferrite system. The observed strength of the magnonic high-harmonic generation signals in the 2DCS of the quantum high-spin model aligns well with the experimental data, showing significant improvement over the corresponding mean-field results.
Abstract We present and benchmark a quantum computing approach to calculate the two-dimensional coherent spectrum (2DCS) of high-spin models. Our approach is based on simulating their real-time dynamics in the presence of several magnetic field pulses, which are spaced in time. We utilize the adaptive variational quantum dynamics simulation algorithm for the study due to its compact circuits, which enables simulations over sufficiently long times to achieve the required resolution in frequency space. Specifically, we consider an antiferromagnetic quantum spin model that incorporates Dzyaloshinskii-Moriya interactions and single-ion anisotropy. The obtained 2DCS spectra exhibit distinct peaks at multiples of the magnon frequency, arising from transitions between different eigenstates of the unperturbed Hamiltonian. By comparing the one-dimensional coherent spectrum with 2DCS, we demonstrate that 2DCS provides a higher resolution of the energy spectrum. We further investigate how the quantum resources scale with the magnitude of the spin using two different binary encodings of the high-spin operators: the standard binary encoding and the Gray code. At low magnetic fields both encodings require comparable quantum resources, but at larger field strengths the Gray code is advantageous. Numerical simulations for spin models with increasing number of sites indicate a polynomial system-size scaling for quantum resources. Lastly, we compare the numerical 2DCS with experimental results on a rare-earth orthoferrite system. The observed strength of the magnonic high-harmonic generation signals in the 2DCS of the quantum high-spin model aligns well with the experimental data, showing significant improvement over the corresponding mean-field results.
Two-dimensional coherent spectroscopy (2DCS) with terahertz radiation offers a promising new avenue for the exploration of many-body phenomena in quantum magnets. This includes the potential diagnosis of fractionalized excitations, for which linear response often struggles due to the indistinguishability of a continuum of fractional excitations from that caused by disorders or impurities. However, the interpretation of the complex results produced by 2DCS remains a challenge. Here, we explore a numerical approach based on exact diagonalization (ED) to help push forward our understanding of 2DCS and, as an example, use it to study the impact of confinement in the one-dimensional transverse field Ising (TFIM) model. We first validate our approach by comparing numerical ED and exact analytical results for the spectroscopic signatures of fractionalized spinons in the model, which exhibits a sharp spinon-echo signal in 2DCS, and show how to understand the finite-size effects inevitable with the inherently small system sizes in ED calculations. By augmenting the model with a longitudinal field, which breaks integrability and introduces confinement of spinon pairs, we observe significant changes to the 2DCS spectrum, such as nonsymmetric broadening of the spinon-echo signal. To further elucidate these changes, we introduce a ``four-kink'' approximation, an effective model that captures the interactions between multiple (confined) spinon pairs. Comparing ED data to these four-kink results provides understanding of the multiple features of the 2DCS spectra, allowing us to interpolate between finite-system size data and the thermodynamic limit. One advantage of our ED approach is its possible extension to finite temperatures, which we explore using thermally pure quantum states and demonstrate how the intensity and spectroscopic patterns of 2DCS change when going beyond the integrable model. Our numerically exact results provide a benchmark for future experiments and theoretical studies relying on approximation methods and pave the way for the exploration of fractionalized excitations in quantum magnets.
Two-dimensional coherent spectroscopy (2DCS) provides insights into the nonlinear response of correlated lattice systems. We simulate multipulse excitations in the Hubbard model using nonequilibrium dynamical mean-field theory to extract the 2DCS signal of weak-coupling antiferromagnets with and without local potential disorder. By comparing calculations with static and dynamic Hartree terms, and analyzing the waiting-time dependence of the signal, we identify the contribution of the collective amplitude mode to the spectroscopic features and the relevant underlying processes. With broadband pulses, the rephasing and nonrephasing peaks at the gap energy are found to be of predominant amplitude mode character. Using narrow-band pulses, we also demonstrate a strong enhancement of these amplitude mode-related signals at a pulse frequency of half the gap size.
Identifying quantum spin liquids, magnon breakdown, or fractionalized excitations in quantum magnets is an ongoing challenge due to the ambiguity in interpreting 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 χ 2 ( ω t , ω τ ) directly in the frequency domain. Applying this to extended Kitaev models describing α − 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 χ 2 ( ω t , ω τ ) 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 α − 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.
We investigate the quantum dynamics of the spins on different Heisenberg antiferromagnetic spin lattice systems. Firstly, we applied the coupled-cluster method to the spin-1/2 antiferromagnetic \emph{XXZ} model on a square lattice by employing an approximation which contains two-body long-range correlations and high-order four-body local correlations. Improvement is found for the ground-state energy, sublattice magnetization, and the critical anisotropy when comparing with the approximation including the two-body correlations alone. We also obtain the full excitation spectrum which is in good agreement with the quantum Monte Carlo results and the high-order spin-wave theory.Secondly, we study the longitudinal excitations of quantum antiferromagnets on a triangular lattice by a recently proposed microscopic many-body approach based on magnon-density waves. We calculate the full longitudinal excitation spectra of the antiferromagnetic Heisenberg model for a general spin quantum number in the isotropic limit. Similar to the square lattice model, we find that, at the center of the first hexagonal Brillouin zone Γ(q=0) and at the magnetic ordering wavevectors ±[Q=(4π/3,0)], the excitation spectra become gapless in the thermodynamic limit, due to the slow, logarithmic divergence of the structure factor. However, these longitudinal modes on two-dimensional models may be considered as quasi-gapped, as any finite-size effect or small anisotropy will induce a large energy gap, when compared with the counterpart of the transverse spin-wave excitations. We have also investigated the excited states of the quasi-one-dimensional quantum antiferromagnets on hexagonal lattices, including the longitudinal modes based on the magnon-density waves. A model Hamiltonian with a uniaxial single-ion anisotropy is first studied by a spin-wave theory based on the one-boson method; the ground state thus obtained is employed for the study of the longitudinal modes. The full energy spectra of both the transverse modes (i.e., magnons) and the longitudinal modes are obtained as functions of the nearest-neighbor coupling and the anisotropy constants. We have found two longitudinal modes due to the non-collinear nature of the triangular antiferromagnetic order, similar to that of the phenomenological field theory approach by Affleck. The excitation energy gaps due to the anisotropy and the energy gaps of the longitudinal modes without anisotropy are then investigated. We then compares our results for the longitudinal energy gaps at the magnetic wavevectors with the experimental results for several antiferromagnetic compounds with both integer and non-integer spin quantum numbers, and we find good agreements after the higher-order contributions are included in our calculations.
Inelastic neutron scattering (INS) has traditionally been one of the primary methods for investigating quantum magnets, particularly in identifying a continuum of excitations as a hallmark of spin fractionalization in quantum spin liquids (QSLs). However, INS faces severe limitations due to its inability to distinguish between such QSL signatures and similar excitation continua arising from highly frustrated magnetic orders with large unit cells or classical spin liquids. In contrast, two-dimensional coherent spectroscopy (2DCS) has emerged as a powerful tool to probe nonlinear excitation dynamics, offering insights into the underlying mechanisms behind these broad spectral features. In this paper, we utilize classical molecular dynamics (MD) techniques to explore the 2DCS responses of frustrated magnets with dominant Kitaev interactions. Comparing the classical and quantum versions of the pure Kitaev model our results indicate both clear similarities, in the form of sharp line features, and clear distinctions, in the locations of these features and in selection rules. Moreover, in the extended KΓΓ′ model, we show that the 2DCS response of the Kitaev spin liquid is completely distinct from that of large unit cell magnetic orders, despite both generating a broad continuum in INS. Additionally, we demonstrate the extreme sensitivity of classical 2DCS to thermal fluctuations and discuss the potential significance of quantum coherence in experimental settings. Overall, our work illustrates the potential of 2DCS in resolving the complex physics underlying ambiguous spin excitation continua, thereby enhancing our understanding of the dynamics in these frustrated systems.