Analyzing atomic liquids and solids by means of two-dimensional Raman spectra in frequency domain

Low-Frequency Spectra of Hydrated Ionic Liquids with Kosmotropic and Chaotropic Anions

In this study, we investigated the water concentration dependence of the intermolecular vibrations of two hydrated ionic liquids (ILs), cholinium dihydrogen phosphate ([ch][dhp]) and cholinium bromide ([ch]Br), using femtosecond Raman-induced Kerr effect spectroscopy (fs-RIKES). The anions of the former and latter hydrated ILs are kosmotropic and chaotropic, respectively. We found that the spectral peak of ∼50 cm-1 shifted to the low-frequency side in hydrated [ch][dhp], indicating the weakening of its intermolecular interactions. In contrast, no change in the peak frequency of the low-frequency band at ∼50 cm-1 was observed with increasing water concentration in hydrated [ch]Br. The vibrational density of states (VDOS) spectra generated from molecular dynamics (MD) simulations were in qualitative agreement with the experimental results. Decomposition analysis of the VDOS spectra for each component revealed that the red shift of the low-frequency band in the hydrated [ch][dhp] upon water addition was essentially due to the contributions of anions and water rather than that of the cholinium cation. We also found from the low-frequency spectra of the two hydrated ILs that they differed in the concentration dependence of the 180 cm-1 band, which is assigned as a hindered translational motion of water molecules combined to form O···O stretching motions. From the relationship between the peak frequency of the low-frequency band and the bulk parameter, which is the square root of the surface tension divided by the density, we found that the peak frequency in the hydrated IL with kosmotropic [dhp]- depends on the bulk parameter, similar to the case for an aqueous solution of the typical deep eutectic solvent reline. However, the peak frequency of the hydrated IL with chaotropic Br- is constant with the bulk parameter.

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.

Phase-space algorithm for simulating quantum nonlinear response functions of bosons using stochastic classical trajectories

Using the positive P-representation of the density matrix, we develop an algorithm for calculating the quantum many-body nonlinear response functions of a system of bosons driven impulsively by external fields. The formalism maps the quantum time evolution of N boson degrees of freedom into a stochastic dynamics of 4N classical degrees of freedom. The first- and the third-order response functions are calculated by propagating the parameters of the P-representation using a set of coupled Langevin equations with multiplicative noise. These parameters serve as classical variables. Two classical ways for computing the response functions are presented. In the nonequilibrium method, an observable is calculated for weak impulsive pulses, and the response functions are obtained by taking its derivatives with respect to the pulse amplitudes. In the alternative, equilibrium simulation, the response functions are expressed in terms of time-correlation functions involving the P-representation parameters and stability matrices representing the perturbation of the trajectories. The stability matrices can be propagated simultaneously with the Langevin equations for the parameters. The formalism is generalized for a many-body boson system coupled to a harmonic bath.

Polarization Selectivity of Third-Order and Fifth-Order Raman Spectroscopies in Liquids and Solids

Polarization selectivity of third-order and fifth-order Raman spectroscopies is examined for both isotropic liquids and periodic lattices. Our approach directly applies the symmetry property of the probed system to decompose the polarization tensor elements into independent components. The polarization selectivity predicted by symmetry analysis is rigorous and applicable to higher-order Raman spectroscopy. The different polarization selectivities of isotropic systems and periodic lattices can be used as a signature of the liquid-solid phase transition.

Two-dimensional infrared surface spectroscopy for CO on Cu(100): Detection of intermolecular coupling of adsorbates

Surface-specific infrared signals obtained by subjecting the system to two infrared laser pulses are calculated for an admixture of CO and isotopic CO on Cu(100) by using molecular dynamics simulation based on a stability matrix formalism. The two-dimensional profiles of the signals in the frequency domain show both diagonal and cross peaks. The former peaks mainly arise from the overtones of the CO and isotopic CO, while the latter represent the couplings between those. As temperature is increased, the phases of cross peaks in a second-order infrared response function change significantly, while those of diagonal peaks are unchanged. The authors show that the phase shifts are originated from the potential anharmonicities due to the electronic interaction between adsorbates. Using a model with two dipole moments, they find that the frustrated rotational mode activated with temperature has effects on the anharmonicity. These results indicate that two-dimensional infrared surface spectroscopy reveals the anharmonic couplings between adsorbates and surface atoms or between adsorbates which cannot be observed in first-order spectroscopy.

Time correlation function and finite field approaches to the calculation of the fifth order Raman response in liquid xenon

The fifth order, two-dimensional Raman response in liquid xenon is calculated via a time correlation function (TCF) theory and the numerically exact finite field method. Both employ classical molecular dynamics simulations. The results are shown to be in excellent agreement, suggesting the efficacy of the TCF approach, in which the response function is written approximately in terms of a single classical multitime TCF.

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.