Time resolved four‐ and six‐wave mixing in liquids. I. Theory

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 Raman spectroscopy of Lennard-Jones liquids via ring-polymer molecular dynamics

The simulation of multidimensional vibrational spectroscopy of condensed-phase systems including nuclear quantum effects is challenging since full quantum-mechanical calculations are still intractable for large systems comprising many degrees of freedom. Here, we apply the recently developed double Kubo transform (DKT) methodology in combination with ring-polymer molecular dynamics (RPMD) for evaluating multi-time correlation functions [K. A. Jung et al., J. Chem. Phys. 148, 244105 (2018)], providing a practical method for incorporating nuclear quantum effects in nonlinear spectroscopy of condensed-phase systems. We showcase the DKT approach in the simulation of the fifth-order two-dimensional (2D) Raman spectroscopy of Lennard-Jones liquids as a prototypical example, which involves nontrivial nonlinear spectroscopic observables of systems described by anharmonic potentials. Our results show that the DKT can faithfully reproduce the 2D Raman response of liquid xenon at high temperatures, where the system behaves classically. In contrast, liquid neon at low temperatures exhibits moderate but discernible nuclear quantum effects in the 2D Raman response compared to the responses obtained with classical molecular dynamics approaches. Thus, the DKT formalism in combination with RPMD simulations enables simulations of multidimensional optical spectroscopy of condensed-phase systems that partially account for nuclear quantum effects.

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.

Inter- and intramolecular hydrogen bonding in phenol derivatives: a model system for poly-L-tyrosine

The ultrafast dynamics of solutions of phenol and two phenol derivatives--hydroquinone (1,4-benzenediol) and pyrocatechol (1,2-benzenediol)--have been studied with Optically Heterodyne-Detected Optical Kerr-Effect (OHD-OKE) spectroscopy. The solvents, methanol and acetonitrile, were selected to provide strong and weak solvent-solute hydrogen-bonding interactions, respectively, while pyrocatechol features an intramolecular hydrogen bond. Together these provide a series of model systems for polypeptides such as polytyrosine, which facilitate the direct study of inter- and intramolecular hydrogen bonding. A broad contribution to the Raman spectral density of the methanol solutions at frequencies between 150 and 300 cm(-1) has been observed that is absent in acetonitrile. This contribution has been assigned to solvent-solute hydrogen-bond stretching vibrations. The OHD-OKE response of poly-L-tyrosine has been measured and was found to contain a similar contribution. Density functional theory geometry optimizations and normal mode calculations have been performed using the B3LYP hybrid functional and 6-311++G** basis set. These have yielded a complete assignment of the low-frequency Raman and far-infrared spectra of pyrocatechol for the first time, which has provided information on the nature of the intramolecular hydrogen bond of pyrocatechol.

Ultrafast dynamics in complex fluids observed through the ultrafast optically-heterodyne-detected optical-Kerr-effect (OHD-OKE)

The ultrafast molecular dynamics of complex fluids have been recorded using the optically-heterodyne-detected optical-Kerr-effect (OHD-OKE). The OHD-OKE method is reviewed and some recent refinements to the method are described. Applications to a range of complex fluids, including microemulsions, polymer melts and solutions, liquid crystal and ionic liquids are surveyed. The level of detail attainable with the OHD-OKE method in these complex fluids is discussed. The prospects for future experiments are discussed.

Molecular dynamics investigation with the time resolved optical Kerr effect on the CS2–C6H6 mixtures

An investigation of the molecular dynamics in pure liquids and in mixtures through the technique of time resolved optical Kerr effect is performed. The samples studied were the mixtures of carbon disulfide (CS(2)) with benzene (C(6)H(6)). The molecular dynamics of the pure liquids is briefly discussed while the main results are obtained for the mixtures. A slow dynamics component is observed for the optical heterodyne detected optical Kerr effect transient decaying exponentially with time constants on picosecond time scale. The fast subpicosend time relaxations are analyzed in terms of the nondiffusive component of the spectral response that is associated with the molecular dynamics. The modifications of the spectrum are quantified, and the explanation of the observed changes is given in terms of the structural interaction configurations that produced changes in the intermolecular potential within which the molecules execute librational motions.

Fifth-Order Raman Spectroscopy of Liquid Benzene: Experiment and Theory

The heterodyned fifth-order Raman response of liquid benzene has been measured and characterized by exploiting the passive-phase stabilization of diffractive optics. This result builds on our previous work with liquid carbon disulfide and extends the spectroscopy to a new liquid for the first time. The all-parallel and Dutch Cross polarization tensor elements are presented for both the experimental results and a finite-field molecular dynamics simulation. The overall response characteristics are similar to those of liquid carbon disulfide: a complete lack of signal along the pump delay, an elongated signal along the probe delay, and a short-lived signal along the time diagonal. Of particular interest is the change in phase between the nuclear and electronic response along the probe delay and diagonal which is not seen in CS2. Good agreement is achieved between the experiment and the finite-field molecular dynamics simulation. The measurement of the low-frequency Raman two-time delay correlation function indicates the intermolecular modes of liquid benzene to be primarily homogeneously broadened and that the liquid loses its nuclear rephasing ability within 300 fs. This rapid loss of nuclear correlations indicates a lack of modal character in the low-frequency motions of liquid benzene. This result is a validation of the general nature of the technique and represents an important step forward with respect to the use of nonlinear spectroscopy to directly access information on the anharmonic motions of liquids.

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.

Temperature- and solvation-dependent dynamics of liquid sulfur dioxide studied through the ultrafast optical Kerr effect

The ultrafast dynamics of liquid sulphur dioxide have been studied over a wide temperature range and in solution. The optically heterodyne-detected and spatially masked optical Kerr effect (OKE) has been used to record the anisotropic and isotropic third-order responses, respectively. Analysis of the anisotropic response reveals two components, an ultrafast nonexponential relaxation and a slower exponential relaxation. The slower component is well described by the Stokes-Einstein-Debye equation for diffusive orientational relaxation. The simple form of the temperature dependence and the agreement between collective (OKE) and single molecule (e.g., NMR) measurements of the orientational relaxation time suggests that orientational pair correlation is not significant in this liquid. The relative contributions of intermolecular interaction-induced and single-molecule orientational dynamics to the ultrafast part of the spectral density are discussed. Single-molecule librational-orientational dynamics appear to dominate the ultrafast OKE response of liquid SO2. The temperature-dependent OKE data are transformed to the frequency domain to yield the Raman spectral density for the low-frequency intermolecular modes. These are bimodal with the lowest-frequency component arising from diffusive orientational relaxation and a higher-frequency component connected with the ultrafast time-domain response. This component is characterized by a shift to higher frequency at lower temperature. This result is analyzed in terms of a harmonic librational oscillator model, which describes the data accurately. The observed spectral shifts with temperature are ascribed to increasing intermolecular interactions with increasing liquid density. Overall, the dynamics of liquid SO2 are found to be well described in terms of molecular orientational relaxation which is controlled over every relevant time range by intermolecular interactions.

Applications of a time correlation function theory for the fifth-order Raman response function I: Atomic liquids

Multidimensional spectroscopy has the ability to provide great insight into the complex dynamics and time-resolved structure of liquids. Theoretically describing these experiments requires calculating the nonlinear-response function, which is a combination of quantum-mechanical time correlation functions R5(t1,t2) was expressed with a two-time, computationally tractable, classical TCF. Writing the response function in terms of classical TCFs brings the full power of atomistically detailed molecular dynamics to the problem. In this paper, the new TCF theory is employed to calculate the fifth-order Raman response function for liquid xenon and investigate several of the polarization conditions for which experiments can be performed on an isotropic system. The theory is shown to reproduce line-shape characteristics predicted by earlier theoretical work.

A new ultrafast technique for measuring the terahertz dynamics of chiral molecules: the theory of optical heterodyne-detected Raman-induced Kerr optical activity

Optical heterodyne-detected Raman-induced Kerr optical activity (OHD-RIKOA) is a nonresonant ultrafast chiroptical technique for measuring the terahertz-frequency Raman spectrum of chirally active modes in liquids, solutions, and glasses of chiral molecules. OHD-RIKOA has the potential to provide much more information on the structure of molecules and the symmetries of librational and vibrational modes than the well-known nonchirally sensitive technique optical heterodyne-detected Raman-induced Kerr-effect spectroscopy (OHD-RIKES). The theory of OHD-RIKOA is presented and possible practical ways of performing the experiments are analyzed.

Orientational and interaction induced dynamics in the isotropic phase of a liquid crystal: Polarization resolved ultrafast optical Kerr effect spectroscopy

The ultrafast dynamics of the isotropic phase of a liquid crystal 4'-pentyl-4-p-biphenylcarbonitrile (5CB) have been investigated using polarization resolved optical Kerr effect spectroscopy. Measurements were made as a function of both temperature and dilution in nonpolar solvents. To separate single molecule and interaction induced components to the relaxation of the induced birefringence, measurements of both the anisotropic and isotropic response were made. The isotropic response was found to be dominated by a damped low-frequency mode of intramolecular origin. There is a minor additional component assigned to an interaction induced contribution. There is at most an extremely weak isotropic signal beyond 1 ps, showing that the picosecond time scale dynamics of 5CB are dominated by orientational relaxation. The isotropic response is independent of temperature in the range studied (0.2-50 K above the nematic to isotropic phase-transition temperature). The anisotropic response exhibits relaxation dynamics on time scales spanning subpicosecond to several hundred picoseconds and beyond. The fastest components are dominated by a librational response, but there are smaller contributions from three low-frequency intramolecular modes, and a contribution from interaction induced effects. The low-frequency spectral density extracted from these data are independent of temperature in the range studied, 0.2-30 K above the phase-transition temperature, but shift to lower frequency on dilution in alkane solvents. In neat 5CB the picosecond time scale orientational dynamics are dominated by temperature-independent reorientation within the pseudonematic domains, while in solution these are disrupted, and the orientational response becomes faster and temperature dependent.

Two-dimensional Raman and infrared vibrational spectroscopy for a harmonic oscillator system nonlinearly coupled with a colored noise bath

Multidimensional vibrational response functions of a harmonic oscillator are reconsidered by assuming nonlinear system-bath couplings. In addition to a standard linear-linear (LL) system-bath interaction, we consider a square-linear (SL) interaction. The LL interaction causes the vibrational energy relaxation, while the SL interaction is mainly responsible for the vibrational phase relaxation. The dynamics of the relevant system are investigated by the numerical integration of the Gaussian-Markovian Fokker-Planck equation under the condition of strong couplings with a colored noise bath, where the conventional perturbative approach cannot be applied. The response functions for the fifth-order nonresonant Raman and the third-order infrared (or equivalently the second-order infrared and the seventh-order nonresonant Raman) spectra are calculated under the various combinations of the LL and the SL coupling strengths. Calculated two-dimensional response functions demonstrate that those spectroscopic techniques are very sensitive to the mechanism of the system-bath couplings and the correlation time of the bath fluctuation. We discuss the primary optical transition pathways involved to elucidate the corresponding spectroscopic features and to relate them to the microscopic sources of the vibrational nonlinearity induced by the system-bath interactions. Optical pathways for the fifth-order Raman spectroscopies from an "anisotropic" medium were newly found in this study, which were not predicted by the weak system-bath coupling theory or the standard Brownian harmonic oscillator model.

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.

Heterodyne-detected fifth-order nonresonant Raman scattering from room temperature CS2

Actively phase-locked heterodyne-detected fifth-order nonresonant Raman scattering from room temperature CS2 has been measured. The experimental signals have similar magnitudes, shapes, and sign changes as calculated responses obtained via molecular dynamics simulations [S. Saito and I. Ohmine, Phys. Rev. Lett. 88, 207401 (2002)]. The measured signals contain sign changes that appear to be associated with the coupling of rotational motions both to each other and to translational motions.