We theoretically analyze two-dimensional coherent spectroscopy (2DCS) signals for disordered superconductors in two limits: One is the narrow-band limit with sinusoidal pulse waves, and the other is the broad-band limit with delta-function pulses. While the 2DCS signal in the narrow-band limit is related to the third-order nonlinear susceptibilities
χ(3)(3Ω;Ω,Ω,Ω) (third harmonic generation) and
χ(3)(Ω;Ω,Ω,−Ω) (ac Kerr effect), we find that in the broad-band limit the signal along the diagonal and horizontal lines in the two-dimensional frequency space is related to another nonlinear susceptibility
χ(3)(Ω;Ω,0,0) (dc Kerr effect). We numerically evaluate those susceptibilities for a lattice model of superconductors based on the BCS mean-field theory and self-consistent Born approximation for impurities. The 2DCS signals in the narrow-band and broad-band limits show threshold and resonance behaviors at the superconducting-gap frequency, respectively, whose physical origin is discussed in light of quasiparticle and Higgs-mode excitations.
Multidimensional coherent spectroscopy (MDCS) has been established in quantum chemistry as a powerful tool for studying the nonlinear response and nonequilibrium dynamics of molecular systems. More recently, the technique has also been applied to correlated electron materials, where the interplay of localized and itinerant states makes the interpretation of the spectra more challenging. Here we use the Keldysh contour representation of effective models and nonequilibrium dynamical mean field theory to systematically study the MDCS signals of prototypical correlated lattice systems. By analyzing the current induced by sequences of ultrashort laser pulses we demonstrate the usefulness of MDCS as a diagnostic tool for excitation pathways and coherent processes in correlated solids. We also show that this technique allows to extract detailed information on the nature and evolution of photo-excited nonequilibrium states.
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.
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.
Multidimensional spectroscopy has become an established tool for disentangling microscopic processes that remain unresolved in conventional linear spectroscopy. Its extension to the terahertz regime now enables analogous questions to be addressed for low-energy collective excitations in quantum materials and synthetic many-body systems. Unlike isolated few-level systems, however, these collective modes generally cannot be described within the standard framework of nonlinear spectroscopy based on a few-level description. Their response is instead governed by their extended nature, interplay of interactions, fluctuations, disorder, dissipation, and coupling to additional degrees of freedom. This thesis develops a theoretical framework for computing and interpreting multidimensional nonlinear spectroscopy in quantum many-body systems. Using real-time nonequilibrium field theory, we formulate nonlinear susceptibilities diagrammatically in terms of propagators, interaction vertices, and self-energies, with particular emphasis on the third-order susceptibility relevant to infrared-active collective excitations. This approach places the established concepts and capabilities of multidimensional spectroscopy within a systematic many-body field-theoretical framework. Furthermore, this formulation makes it possible to move beyond the established probing capabilities, extending multidimensional spectroscopy beyond few-level descriptions and revealing new signatures and functionalities specific to interacting many-body systems. At the same time, it provides a controlled framework for constructing approximations that preserve the underlying symmetries and conservation laws of the system while remaining amenable to both analytical and numerical calculations. By incorporating quantum and thermal fluctuations within the random-phase approximation, we demonstrate that two-dimensional spectra can reveal the growth of fluctuations, identify field-theoretical bound states, and characterize their coupling to many-body continua. We further develop a diagrammatic description of rephasing in extended bosonic systems in the presence of static disorder. We identify the disorder vertex corrections responsible for echo formation and show how echo spectroscopy separates disorder-induced broadening from intrinsic relaxation. This framework is extended to coupled modes and strongly interacting regimes, where interactions can either enhance or suppress rephasing. A further result is the identification of a nonrephasing four-wave-mixing signal that is critically enhanced near a phase transition, in contrast to simpler nonlinear protocols like third-harmonic generation. We apply this framework to electromagnons in NiI
2, Josephson plasmons in La
2−xSr
xCuO
4, and the low-energy collective excitations in coupled one-dimensional Bose-Einstein condensates. Across these systems, we show how multidimensional spectroscopy can disentangle disorder from intrinsic relaxation, reveal the effects of interactions and fluctuations, identify coupling between collective modes and continua, and extract microscopic information that remains obscured in conventional spectroscopy. Taken together, the results of this thesis provide a diagrammatic framework for calculating and interpreting multidimensional nonlinear spectroscopic signals in many-body systems, while establishing this technique as a selective probe of collective dynamics across a broad range of quantum materials and synthetic platforms.
Detecting Kitaev interactions in real materials remains challenge, as conventional experimental techniques often have difficulty distinguishing fractionalized excitations from other normal contributions. Terahertz two-dimensional coherent spectroscopy (2DCS) offers a novel approach for probing many-body phenomena, such as exotic excitations in quantum magnets. Motivated by recent experiments on CoNb
2O
6 and the development of the terahertz spectroscopy in Kitaev quantum spin liquid, we proposed a twisted Kitaev model for CoNb
2O
6 and determined the precise twist angle according to experimental specific-heat phase diagram. With this calibrated model, we found that non-rephasing diagonal and rephasing anti-diagonal signals appear in the 2DCS nonlinear response. The
x and
y components of the spin superexchange interactions split the rephasing signals into a grid of discrete peaks. We further demonstrate that the diagonal and the discrete rephasing signals primarily originate from two-spinon and four-spinon excitation processes based on numerical projection method. These findings indicate that even weak Kitaev interactions in quantum materials can be effectively detected via two-dimensional coherent spectroscopy .