Interplay of inhomogeneity and anharmonicity in 2D Raman spectroscopy of liquids

Low-frequency anharmonic couplings in bromoform revealed from 2D Raman-THz spectroscopy: from the liquid to the crystalline phase

Two-dimensional (2D) Raman-THz spectroscopy in the frequency up to 7 THz has been applied to study the crystalline beta-phase of bromoform (CHBr3). As for liquid CHBr3, cross peaks are observed, which however sharpen up in the crystalline sample and split into assignable sub-contributions. In the Raman dimension, the frequency positions of these cross peaks coincide with the intramolecular bending modes of the CHBr3 molecules, and in the THz dimension with the IR active lattice modes of the crystal. This work expands the applicability of this new 2D spectroscopic technique to solid samples at cryogenic temperatures. Furthermore, it provides new experimental evidence that the cross peaks indeed originate from the coupling between intra- and intermolecular vibrational modes.

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.

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.

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.

Calculating fifth-order Raman signals for various molecular liquids by equilibrium and nonequilibrium hybrid molecular dynamics simulation algorithms

The fifth-order two-dimensional (2D) Raman signals have been calculated from the equilibrium and nonequilibrium (finite field) molecular dynamics simulations. The equilibrium method evaluates response functions with equilibrium trajectories, while the nonequilibrium method calculates a molecular polarizability from nonequilibrium trajectories for different pulse configurations and sequences. In this paper, we introduce an efficient algorithm which hybridizes the existing two methods to avoid the time-consuming calculations of the stability matrices which are inherent in the equilibrium method. Using nonequilibrium trajectories for a single laser excitation, we are able to dramatically simplify the sampling process. With this approach, the 2D Raman signals for liquid xenon, carbon disulfide, water, acetonitrile, and formamide are calculated and discussed. Intensities of 2D Raman signals are also estimated and the peak strength of formamide is found to be only five times smaller than that of carbon disulfide.

Two-dimensional Raman spectra of atomic solids and liquids

We calculate third- and fifth-order Raman spectra of simple atoms interacting through a soft-core potential by means of molecular-dynamics (MD) simulations. The total polarizability of molecules is treated by the dipole-induced dipole model. Two- and three-body correlation functions of the polarizability at various temperatures are evaluated from equilibrium MD simulations based on a stability matrix formulation. To analyze the processes involved in the spectroscopic measurements, we divide the fifth-order response functions into symmetric and antisymmetric integrated response functions; the symmetric one is written as a simple three-body correlation function, while the antisymmetric one depends on a stability matrix. This analysis leads to a better understanding of the time scales and molecular motions that govern the two-dimensional (2D) signal. The 2D Raman spectra show novel differences between the solid and liquid phases, which are associated with the decay rates of coherent motions. On the other hand, these differences are not observed in the linear Raman spectra.

Higher-order optical correlation spectroscopy in liquids

Linear optical spectroscopies have long been used to study the behavior of liquids. Laser technology has progressed to the point that it has become possible to perform nonlinear optical experiments that probe higher-order correlation functions in liquids, opening a new window into our understanding of the microscopic details of solution-phase processes. Here we review advances that have been made in recent years in employing higher-order electronic and vibrational spectroscopies to study liquid-state dynamics and structure.

Off-resonant fifth-order response function for two-dimensional Raman spectroscopy of liquids CS2 and H2O

The off-resonant fifth-order response functions for two-dimensional (2D) Raman spectroscopy of molecular liquids CS2 and H2O are investigated by using molecular dynamics calculation. This spectroscopy, able to deal with a phase space dynamics, shows the existence of nodal lines in several polarization tensor elements [see L. Kaufman et al., Phys. Rev. Lett. 88, 207402 (2002) for experimental results]. The nodal property is found to arise from the dynamical couplings among rotational modes, not accounted for in a normal mode picture. The effects of anharmonicities and "mode coupling through polarizability" are also investigated by comparing the 2D Raman signal with a constant temperature velocity reassignment echo method.