Ultrafast 2D spectroscopy uses correlated multiple light−matter interactions for retrieving dynamic features that may otherwise be hidden under the linear spectrum; its extension to the terahertz regime of the electromagnetic spectrum, where a rich variety of material degrees of freedom reside, remains an experimental challenge. We report a demonstration of ultrafast 2D terahertz spectroscopy of gas-phase molecular rotors at room temperature. Using time-delayed terahertz pulse pairs, we observe photon echoes and other nonlinear signals resulting from molecular dipole orientation induced by multiple terahertz field−dipole interactions. The nonlinear time domain orientation signals are mapped into the frequency domain in 2D rotational spectra that reveal J-state-resolved nonlinear rotational dynamics. The approach enables direct observation of correlated rotational transitions and may reveal rotational coupling and relaxation pathways in the ground electronic and vibrational state.ultrafast spectroscopy | multidimensional coherent spectroscopy | terahertz | rotational dynamics R ecent years have witnessed increasing interest in 2D infrared (IR) vibrational spectroscopy techniques for studying structural dynamics and correlations between coupled molecular motions in biological systems such as water, proteins, and DNA (1-3). Multidimensional optical spectroscopies were applied to probe the high-order correlations of excitons in quantum wells and organic complexes (4, 5). However, it remains an experimental challenge to extend multidimensional spectroscopies into the terahertz (THz) frequency range where a rich variety of material degrees of freedom, including gas molecular rotations, lattice vibrations in solids, magnetization dynamics in magnetically ordered materials, and many others (6), find fundamental and technological importance. Despite challenges involved, there have been examples of 2D THz and 2D THz−Raman spectroscopies in the studies of electronic nonlinearities in solids (7,8), coupled vibrational modes in liquid molecules (9), and hydrogen bond dynamics in water (10).Rotations of gas molecules have been the subjects of intensive recent efforts in coherent spectroscopy and coherent control, motivated by interests in rotational angular momentum and energy relaxation processes, high-order optical interactions with multilevel quantum systems, and applications of molecular alignment and orientation for various forms of orbital tomography, for example using high-harmonic generation to probe the molecular orbital structures of rotating molecules (11,12). Strong optical fields have been used to drive high-order rotational coherences and to populate high-lying rotational levels (13), to control molecular alignment dynamics (14), and to induce rotational alignment echoes (15, 16). Strong microwave fields can drive high-order rotational responses and have been used for double-resonance and 2D rotational spectroscopies (17-22), mostly at low rotational temperatures. In recent work, we demonstrated second-order int...