Analysis of 2D THz-Raman spectroscopy using a non-Markovian Brownian oscillator model with nonlinear system-bath interactions

Classical and quantum thermodynamics described as a system-bath model: The dimensionless minimum work principle

We formulate a thermodynamic theory applicable to both classical and quantum systems. These systems are depicted as thermodynamic system-bath models capable of handling isothermal, isentropic, thermostatic, and entropic processes. Our approach is based on the use of a dimensionless thermodynamic potential expressed as a function of the intensive and extensive thermodynamic variables. Using the principles of dimensionless minimum work and dimensionless maximum entropy derived from quasi-static changes of external perturbations and temperature, we obtain the Massieu-Planck potentials as entropic potentials and the Helmholtz-Gibbs potentials as free energy. These potentials can be interconverted through time-dependent Legendre transformations. Our results are verified numerically for an anharmonic Brownian system described in phase space using the low-temperature quantum Fokker-Planck equations in the quantum case and the Kramers equation in the classical case, both developed for the thermodynamic system-bath model. Thus, we clarify the conditions for thermodynamics to be valid even for small systems described by Hamiltonians and establish a basis for extending thermodynamics to non-equilibrium conditions.

On selection rules in two-dimensional terahertz-infrared-visible spectroscopy

Two-dimensional terahertz-infrared-visible (2D TIRV) spectroscopy directly measures the coupling between quantum high-frequency vibrations and classical low-frequency modes of molecular motion. In addition to coupling strength, the signal intensity in 2D TIRV spectroscopy can also depend on the selection rules of the excited transitions. Here, we explore the selection rules in 2D TIRV spectroscopy by studying the coupling between the high-frequency CH3 stretching and low-frequency vibrations of liquid dimethyl sulfoxide (DMSO). Different excitation pathways are addressed using variations in laser pulse timing and different polarizations of exciting pulses and detected signals. The DMSO signals generated via different excitation pathways can be readily distinguished in the spectrum. The intensities of different excitation pathways vary unequally with changes in polarization. We explain how this difference stems from the intensities of polarized and depolarized Raman and hyper-Raman spectra of high-frequency modes. These results apply to various systems and will help design and interpret new 2D TIRV spectroscopy experiments.

Experimental Access to Mode-Specific Coupling between Quantum Molecular Vibrations and Classical Bath Modes

The interaction of quantum-mechanical systems with a fluctuating thermal environment (bath) is fundamental to molecular mechanics and energy transport/dissipation. Its complete picture requires mode-specific measurements of this interaction and an understanding of its nature. Here, we present a combined experimental and theoretical study providing detailed insights into the coupling between a high-frequency vibrational two-level system and thermally excited terahertz modes. Experimentally, two-dimensional terahertz-infrared-visible spectroscopy reports directly on the coupling between quantum oscillators represented by CH3 streching vibrations in liquid dimethyl sulfoxide and distinct low-frequency modes. Theoretically, we present a mixed quantum-classical formalism of the sample response to enable the simultaneous quantum description of high-frequency oscillators and a classical description of the bath. We derive the strength and nature of interaction and find different coupling between CH3 stretch and low-frequency modes. This general approach enables quantitative and mode-specific analysis of coupled quantum and classical dynamics in complex chemical systems.

Discretized hierarchical equations of motion in mixed Liouville-Wigner space for two-dimensional vibrational spectroscopies of liquid water

A model of a bulk water system describing the vibrational motion of intramolecular and intermolecular modes is constructed, enabling analysis of its linear and nonlinear vibrational spectra as well as the energy transfer processes between the vibrational modes. The model is described as a system of four interacting anharmonic oscillators nonlinearly coupled to their respective heat baths. To perform a rigorous numerical investigation of the non-Markovian and nonperturbative quantum dissipative dynamics of the model, we derive discretized hierarchical equations of motion in mixed Liouville-Wigner space, with Lagrange-Hermite mesh discretization being employed in the Liouville space of the intramolecular modes and Lagrange-Hermite mesh discretization and Hermite discretization in the Wigner space of the intermolecular modes. One-dimensional infrared and Raman spectra and two-dimensional terahertz-infrared-visible and infrared-infrared-Raman spectra are computed as demonstrations of the quantum dissipative description provided by our model.

Low-frequency anharmonic couplings in crystalline bromoform: Theory

Theoretical calculations of the low-frequency anharmonic couplings of the β-phase of crystalline bromoform are presented based on density functional theory quantum chemistry calculations. The electrical and mechanical anharmonicities between intra- and intermolecular modes are calculated, revealing that the electrical anharmonicity dominates the cross-peak intensities in the 2D Raman-THz response and crystalline, as well as liquid, bromoform. Furthermore, the experimentally observed difference in relative cross-peak intensities between the two intramolecular modes of bromoform and the intermolecular modes can be explained by the C_3v-symmetry of bromoform in combination with orientational averaging. The good agreement with the experimental results provides further evidence for our interpretation that the 2D Raman-THz response of bromoform is, indeed, related to the anharmonic coupling between the intra- and intermolecular modes.

Equilibrium-nonequilibrium ring-polymer molecular dynamics for nonlinear spectroscopy

Two-dimensional Raman and hybrid terahertz-Raman spectroscopic techniques provide invaluable insight into molecular structures and dynamics of condensed-phase systems. However, corroborating experimental results with theory is difficult due to the high computational cost of incorporating quantum-mechanical effects in the simulations. Here, we present the equilibrium-nonequilibrium ring-polymer molecular dynamics (RPMD), a practical computational method that can account for nuclear quantum effects on the two-time response function of nonlinear optical spectroscopy. Unlike a recently developed approach based on the double Kubo transformed (DKT) correlation function, our method is exact in the classical limit, where it reduces to the established equilibrium-nonequilibrium classical molecular dynamics method. Using benchmark model calculations, we demonstrate the advantages of the equilibrium-nonequilibrium RPMD over classical and DKT-based approaches. Importantly, its derivation, which is based on the nonequilibrium RPMD, obviates the need for identifying an appropriate Kubo transformed correlation function and paves the way for applying real-time path-integral techniques to multidimensional spectroscopy.

Rate of Insulin Dimer Dissociation: Interplay between Memory Effects and Higher Dimensionality

We calculate the rate of dissociation of an insulin dimer into two monomers in water. The rate of this complex reaction is determined by multiple factors that are elucidated. By employing advanced sampling techniques, we first obtain the reaction free energy surface for the dimer dissociation as a function of two order parameters, namely, the distance between the center-of-mass of two monomers (R) and the number of cross-contacts (Q) among the backbone C_α atoms of two monomers. We then construct an orthogonal 2D reaction energy surface by introducing the reaction coordinate X to denote the minimum energy pathway and a conjugate coordinate Y that spans the orthogonal direction. The free energy landscape is rugged with multiple maxima and minima. We calculate the rate by employing not only the non-Markovian multidimensional rate theory but also several other theoretical approaches. The necessary reaction frequencies and the frictions are calculated from the time correlation function formalism. Our best estimate of the rate is 0.4 μs^-1. Our study reveals interesting opposite influences of dimensionality and memory in determining the rate constant of the reaction. We gain interesting insights into the dimer dissociation process by looking directly at the trajectories obtained from molecular dynamics simulation.

Open Quantum System Response from the Hierarchy of Pure States

The spectral properties of a quantum system are essential when probing theoretical predictions against experimental data. For an open quantum system strongly interacting with its environment, spectral features are challenging to calculate. Here we demonstrate that the stochastic Hierarchy of Pure States (HOPS) approach is well suited to calculate the response of an open quantum system to a, possibly strong, coherent probe driving. For weak driving, where Kubo's linear response theory is applicable, it turns out that the HOPS method is highly efficient since fluctuations inherent to the stochastic dynamics cancel for the response function and, thus, allow us to obtain the susceptibility easily. Our results are in agreement with experimental data for a strongly damped spin system showing that the transition from oscillatory to overdamped motion is also reflected by the transmission spectrum. As a further application we demonstrate that the susceptibility, quantifying the amplitude of the response, as a function of temperature exhibits a maximum which is the hallmark of stochastic resonance. Beyond the linear regime, the exact open system dynamics shows the asymptotic Floquet state. We use the topic of probe driving and response to present the HOPS approach in a novel and self-contained way. This includes the importance sampling scheme which yields the nonlinear HOPS as well as the stochastic treatment of a thermal initial environmental state within the zero temperature formalism. Special attention is given to the exponential representation of the algebraic Ohmic bath correlation function and the truncation condition for the hierarchy.

Full molecular dynamics simulations of molecular liquids for single-beam spectrally controlled two-dimensional Raman spectroscopy

Single-beam spectrally controlled (SBSC) two-dimensional (2D) Raman spectroscopy is a unique 2D vibrational measurement technique utilizing trains of short pulses that are generated from a single broadband pulse by pulse shaping. This approach overcomes the difficulty of 2D Raman spectroscopy in dealing with small-signal extraction and avoids complicated low-order cascading effects, thus providing a new possibility for measuring the intramolecular and intermolecular modes of molecular liquids using fifth-order 2D Raman spectroscopy. Recently, for quantitatively investigating the mode-mode coupling mechanism, Hurwitz et al. [Opt. Express 28, 3803 (2020)] have developed a new pulse design for this measurement to separate the contributions of the fifth- and third-order polarizations, which are often overlapped in the original single-beam measurements. Here, we describe a method for simulating these original measurements and the new 2D Raman measurements on the basis of a second-order response function approach. We carry out full molecular dynamics simulations for carbon tetrachloride and liquid water using an equilibrium-nonequilibrium hybrid algorithm, with the aim of explaining the key features of the SBSC 2D Raman spectroscopic method from a theoretical point of view. The predicted signal profiles and intensities provide valuable information that can be applied to 2D spectroscopy experiments, allowing them to be carried out more efficiently.

Two-dimensional terahertz spectroscopy of condensed-phase molecular systems

Nonlinear terahertz (THz) spectroscopy relies on the interaction of matter with few-cycle THz pulses of electric field amplitudes up to megavolts/centimeter (MV/cm). In condensed-phase molecular systems, both resonant interactions with elementary excitations at low frequencies such as intra- and intermolecular vibrations and nonresonant field-driven processes are relevant. Two-dimensional THz (2D-THz) spectroscopy is a key method for following nonequilibrium processes and dynamics of excitations to decipher the underlying interactions and molecular couplings. This article addresses the state of the art in 2D-THz spectroscopy by discussing the main concepts and illustrating them with recent results. The latter include the response of vibrational excitations in molecular crystals up to the nonperturbative regime of light-matter interaction and field-driven ionization processes and electron transport in liquid water.

Numerically "exact" approach to open quantum dynamics: The hierarchical equations of motion (HEOM)

An open quantum system refers to a system that is further coupled to a bath system consisting of surrounding radiation fields, atoms, molecules, or proteins. The bath system is typically modeled by an infinite number of harmonic oscillators. This system-bath model can describe the time-irreversible dynamics through which the system evolves toward a thermal equilibrium state at finite temperature. In nuclear magnetic resonance and atomic spectroscopy, dynamics can be studied easily by using simple quantum master equations under the assumption that the system-bath interaction is weak (perturbative approximation) and the bath fluctuations are very fast (Markovian approximation). However, such approximations cannot be applied in chemical physics and biochemical physics problems, where environmental materials are complex and strongly coupled with environments. The hierarchical equations of motion (HEOM) can describe the numerically "exact" dynamics of a reduced system under nonperturbative and non-Markovian system-bath interactions, which has been verified on the basis of exact analytical solutions (non-Markovian tests) with any desired numerical accuracy. The HEOM theory has been used to treat systems of practical interest, in particular, to account for various linear and nonlinear spectra in molecular and solid state materials, to evaluate charge and exciton transfer rates in biological systems, to simulate resonant tunneling and quantum ratchet processes in nanodevices, and to explore quantum entanglement states in quantum information theories. This article presents an overview of the HEOM theory, focusing on its theoretical background and applications, to help further the development of the study of open quantum dynamics.

A Feynman diagram description of the 2D-Raman-THz response of amorphous ice

The 2D-Raman-THz response in all possible time-orderings (Raman-THz-THz, THz-Raman-THz, and THz-THz-Raman) of amorphous water ice is calculated in two ways: from atomistic molecular dynamics simulations and with the help of a Feynman diagram model, the latter of which power-expands the potential energy surface and the dipole and polarizability surfaces up to leading order. Comparing both results allows one to dissect the 2D-Raman-THz response into contributions from mechanical anharmonicity, as well as electrical dipole and polarizability anharmonicities. Mechanical anharmonicity dominates the 2D-Raman-THz response of the hydrogen-bond stretching and hydrogen-bond bending bands of water, and dipole anharmonicity dominates that of the librational band, while the contribution of polarizability anharmonicity is comparably weak. A distinct echo of the hydrogen-bond stretching band is observed for the THz-Raman-THz pulse sequence, again dominated by mechanical anharmonicity. A peculiar mechanism is discussed, which is based on the coupling between the many normal modes within the hydrogen-bond stretching band and which will inevitably generate such an echo for an amorphous structure.

Distinguishing Nonlinear Terahertz Excitation Pathways with Two-Dimensional Spectroscopy

High-field terahertz (THz) spectroscopy is enabling the ultrafast study and control of matter in new and exciting ways. However, when intense electromagnetic pulses are used in any kind of pump-probe spectroscopy, several nonlinear excitation pathways can result, leading to scenarios that required the development of multidimensional spectroscopies to illuminate the observed dynamics. Here we demonstrate a clear example where two-dimensional (2D) THz vibrational spectroscopy is needed to distinguish between nonlinear-excitation pathways in CdWO_{4}. We nonlinearly excite a set of Raman-active vibrational modes in CdWO_{4} with broadband THz pulses, and 2D spectroscopy allows us to determine the dominant excitation pathway. We provide a general framework for 2D THz and multi-THz nonlinear phonon spectroscopy in solid systems, which has important implications in contributing needed clarity to the nascent field of nonlinear phononics.

Reconsideration of the relaxational and vibrational line shapes of liquid water based on ultrabroadband dielectric spectroscopy

The Debye relaxation function is widely used to describe the large dielectric dispersion of ambient water around 20 GHz. However, from a theoretical point of view, this function is supposed to give incorrect predictions at high frequencies owing to the inappropriate assumption that inertial effects and intermolecular interactions do not affect the relaxation dynamics. Our ultrabroadband spectroscopy investigation of liquid water ranging from 500 MHz to 400 THz did demonstrate that the Debye function is inaccurate far above the microwave region. As an alternative, we tried a stochastic frequency modulation (SFM) model assuming instantaneous modification of the line shapes by the correlation with the surrounding system. The SFM relaxation model reproduced the experimental dielectric spectra up to 400 THz, showing that the hydrogen-bond dynamics are associated with the inertial effect that causes the non-exponential relaxation behaviour in a very short time (typically 25 fs). Within the framework of this relaxation model, the hindered translation modes are able to be approximated as fast modulation (homogeneous) line shapes because the interaction time with frequency modulation is too short. Compared with them, the libation mode is found to have a relatively slow modulation (inhomogeneous) origin, where disturbance of water hydrogen bonds induced by the hindered translations leads to fluctuations in the libration frequency.

Operational Markov Condition for Quantum Processes

We derive a necessary and sufficient condition for a quantum process to be Markovian which coincides with the classical one in the relevant limit. Our condition unifies all previously known definitions for quantum Markov processes by accounting for all potentially detectable memory effects. We then derive a family of measures of non-Markovianity with clear operational interpretations, such as the size of the memory required to simulate a process or the experimental falsifiability of a Markovian hypothesis.

Intermolecular charge fluxes and far-infrared spectral intensities of liquid formamide

Intermolecular charge fluxes induced by hydrogen-bond length modulations occurring upon molecular librations lead to intensity enhancement of the far-infrared spectrum.

Terahertz echoes reveal the inhomogeneity of aqueous salt solutions

The structural and dynamical properties of water are known to be affected by ion solvation. However, a consistent molecular picture that describes how and to what extent ions perturb the water structure is still missing. Here we apply 2D Raman-terahertz spectroscopy to investigate the impact of monatomic cations on the relaxation dynamics of the hydrogen-bond network in aqueous salt solutions. The inherent ability of multidimensional spectroscopy to deconvolute heterogeneous relaxation dynamics is used to reveal the correlation between the inhomogeneity of the collective intermolecular hydrogen-bond modes and the viscosity of a salt solution. Specifically, we demonstrate that the relaxation time along the echo direction t₁= t₂ correlates with the capability of a given cation to 'structure' water. Moreover, we provide evidence that the echo originates from the water-water modes, and not the water-cation modes, which implies that cations can structure the hydrogen-bond network to a certain extent.

Effects of Intermolecular Charge Transfer in Liquid Water on Raman Spectra

The low-frequency vibrational spectrum of liquid water is composed of contributions from the intermolecular librational and translation modes. The existence of these two modes introduces difficulty into the simulation of experimentally obtained Raman spectra. We constructed a polarizability function for a water model that includes intramolecular charge flow (CF) effects, intermolecular charge transfer (CT) effects, and intermolecular dipole-induced-dipole (DID) effects. We computed the one-dimensional (1D) Raman and terahertz (THz) spectra with all of these effects included (CFCT-DID) and compared with experimental spectra. We find that the CFCT-DID function provides a better description of the experimental results, because the CT effects reduce the polarizability only for translational motion to which parallelly polarized (VV) and perpendicularly polarized (VH) Raman spectra are sensitive. In our calculations of two-dimensional (2D) Raman and THz-Raman spectra, we observe the enhancement of echo signals in both cases. The details of the CFCT-DID function, along with its source code, are provided in the Supporting Information.

Coherent two-dimensional terahertz-terahertz-Raman spectroscopy

We present 2D terahertz-terahertz-Raman (2D TTR) spectroscopy, the first technique, to our knowledge, to interrogate a liquid with multiple pulses of terahertz (THz) light. This hybrid approach isolates nonlinear signatures in isotropic media, and is sensitive to the coupling and anharmonicity of thermally activated THz modes that play a central role in liquid-phase chemistry. Specifically, by varying the timing between two intense THz pulses, we control the orientational alignment of molecules in a liquid, and nonlinearly excite vibrational coherences. A comparison of experimental and simulated 2D TTR spectra of bromoform (CHBr3), carbon tetrachloride (CCl4), and dibromodichloromethane (CBr2Cl2) shows previously unobserved off-diagonal anharmonic coupling between thermally populated vibrational modes.

Nonlinear terahertz coherent excitation of vibrational modes of liquids

We report the first coherent excitation of intramolecular vibrational modes via the nonlinear interaction of a TeraHertz (THz) light field with molecular liquids. A terahertz-terahertz-Raman pulse sequence prepares the coherences with a broadband, high-energy, (sub)picosecond terahertz pulse, that are then measured in a terahertz Kerr effect spectrometer via phase-sensitive, heterodyne detection with an optical pulse. The spectrometer reported here has broader terahertz frequency coverage, and an increased sensitivity relative to previously reported terahertz Kerr effect experiments. Vibrational coherences are observed in liquid diiodomethane at 3.66 THz (122 cm(-1)), and in carbon tetrachloride at 6.50 THz (217 cm(-1)), in exact agreement with literature values of those intramolecular modes. This work opens the door to 2D spectroscopies, nonlinear in terahertz field, that can study the dynamics of condensed-phase molecular systems, as well as coherent control at terahertz frequencies.

Onsets of hierarchy truncation and self-consistent Born approximation with quantum mechanics prescriptions invariance

The issue of efficient hierarchy truncation is related to many approximate theories. In this paper, we revisit this issue from both the numerical efficiency and quantum mechanics prescription invariance aspects. The latter requires that the truncation approximation made in Schrödinger picture, such as the quantum master equations and their self-consistent-Born-approximation improvements, should be transferable to their Heisenberg-picture correspondences, without further approximations. We address this issue with the dissipaton equation of motion (DEOM), which is a unique theory for the dynamics of not only reduced systems but also hybrid bath environments. We also highlight the DEOM theory is not only about how its dynamical variables evolve in time, but also the underlying dissipaton algebra. We demonstrate this unique feature of DEOM with model systems and report some intriguing nonlinear Fano interferences characteristics that are experimentally measurable.

Notes on simulating two-dimensional Raman and terahertz-Raman signals with a full molecular dynamics simulation approach

Recent developments in two-dimensional (2D) THz-Raman and 2D Raman spectroscopies have created the possibility for quantitatively investigating the role of many dynamic and structural aspects of the molecular system. We explain the significant points for properly simulating 2D vibrational spectroscopic studies of intermolecular modes using the full molecular dynamics approach, in particular, regarding the system size, the treatment of the thermostat, and inclusion of an Ewald summation for the induced polarizability. Moreover, using the simulation results for water employing various polarization functions, we elucidate the roles of permanent and induced optical properties in determining the 2D profiles of the signal.