We present a theoretical approach for computing and analyzing two-dimensional spectroscopy of bosonic collective excitations in disordered many-body systems. Specifically, we employ the Keldysh formalism to derive, within a nonperturbative treatment of disorder effects, the third-order nonlinear response and obtain two-dimensional spectroscopy maps. In the weak nonlinear regime of our formalism, we demonstrate the ability of the echo peak to distinguish between elastic and inelastic scattering processes, in perfect agreement with the intuition developed in isolated two-level systems. Furthermore, we discuss unique many-body effects on the echo peak signature arising from interaction-induced quantum fluctuations. In particular, we show that these quantum fluctuations induce a finite nonrephasable broadening and examine how the echo peak is influenced by the attractive or repulsive nature of the collective excitations.
We present a theoretical approach for computing and analyzing two-dimensional spectroscopy of bosonic collective excitations in disordered many-body systems. Specifically, we employ the Keldysh formalism to derive, within a nonperturbative treatment of disorder effects, the third-order nonlinear response and obtain two-dimensional spectroscopy maps. In the weak nonlinear regime of our formalism, we demonstrate the ability of the echo peak to distinguish between elastic and inelastic scattering processes, in perfect agreement with the intuition developed in isolated two-level systems. Furthermore, we discuss unique many-body effects on the echo peak signature arising from interaction-induced quantum fluctuations. In particular, we show that these quantum fluctuations induce a finite nonrephasable broadening and examine how the echo peak is influenced by the attractive or repulsive nature of the collective excitations.
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.
Emergent collective excitations constitute a hallmark of interacting quantum many-body systems, yet in solid-state platforms their study has been largely limited by the constraints of linear-response probes and by finite momentum resolution. We propose to overcome these limitations by combining the spatial resolution of ultracold atomic systems with the nonlinear probing capabilities of two-dimensional spectroscopy (2DS). As a concrete illustration, we analyze momentum-resolved 2DS of the quantum sine-Gordon model describing the low-energy dynamics of two weakly coupled one-dimensional Bose-Einstein condensates. This approach reveals distinctive many-body signatures, most notably asymmetric cross-peaks reflecting the interplay between isolated (B_{2} breather) and continuum (B_{1} pair) modes. The protocol further enables direct characterization of anharmonicity and disorder, establishing momentum-resolved 2DS as both a powerful diagnostic for quantum simulators and a versatile probe of correlated quantum matter.
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.
The analysis of nonlinear spectroscopy, widely used to study the dynamics and structures of condensed-phase matter, typically employs a perturbative approach noticing the weak interaction between the laser and the matter of interest. However, such a perturbative approach is no longer applicable once the interaction between the laser and the matter is strong. We adapt the method of quantum dynamical evolution into the calculation of signal and present the response formalism of the nonlinear spectroscopy in a nonperturbative approach. In this approach, we demonstrate that, in addition to the third-order term in the perturbative method, the higher-order terms have essential contributions to the nonlinear signal of the two-pulse and three-pulse photon echo (2PPE and 3PPE). The detailed calculations are demonstrated with the example of a three-level V-type system, which is widely used in the studies of quantum optics. We consider the effect of the environment via a pure dephasing mechanism with both the localized modes of each molecule and the shared modes between molecules.
The analysis of nonlinear spectroscopy, widely used to study the dynamics and structures of condensed-phase matter, typically employs a perturbative approach noticing the weak interaction between the laser and the matter of interest. However, such perturbative approach is no longer applicable once the interaction between the laser and the matter is strong. We adapt the method of quantum dynamical evolution into the calculation of signal and present the response formalism of the nonlinear spectroscopy in a nonperturbative approach. In this new approach, we demonstrate that in addition to the third-order term in the perturbative method, the higher-order terms have essential contributions to the nonlinear signal of the two-pulse and three-pulse photon echo (2PPE and 3PPE). The detailed calculations are demonstrated with the example of a three-level V-type system, which is widely used in the studies of quantum optics. We consider the effect of the environment via a pure dephasing mechanism with both the localized modes of each molecule and the shared modes between molecules.