The molecular origins of the two-dimensional Raman spectrum of an atomic liquid. II. Instantaneous-normal-mode theory

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.

Multidimensional Measures of Response and Fluctuations in Stochastic Dynamical Systems

A new class of experiments is proposed which involve multiple measurements combined with multiple perturbations of a nonlinear classical complex system. A family of multipoint n + m - 1 dimensional measures R^(n,m) that provide complimentary information on complex systems is obtained by combining m non linear stimuli and n measurements. They represent the combined effect of causal response and noncausal spontaneous fluctuations. The proposed signals can be measured either in the frequency domain or in the time domain. Generalized fluctuation-dissipation relations which hold in the nonlinear regime restrict the number of independent techniques. Two dimensional correlation plots derived from such measurements can reveal various types of couplings among collective modes.

Coherent multidimensional vibrational spectroscopy of biomolecules: concepts, simulations, and challenges

The response of complex molecules to sequences of femtosecond infrared pulses provides a unique window into their structure, dynamics, and fluctuating environments. Herein we survey the basic principles of modern two-dimensional infrared (2DIR) spectroscopy, which analogous to those of multidimensional NMR spectroscopy. The perturbative approach for computing the nonlinear optical response of coupled localized chromophores is introduced and applied to the amide backbone transitions of proteins, liquid water, membrane lipids, and amyloid fibrils. The signals are analyzed using classical molecular dynamics simulations combined with an effective fluctuating Hamiltonian for coupled localized anharmonic vibrations whose dependence on the local electrostatic environment is parameterized by an ab initio map. Several simulation methods, (cumulant expansion of Gaussian fluctuation, quasiparticle scattering, the stochastic Liouville equations, direct numerical propagation) are surveyed. Chirality-induced techniques which dramatically enhance the resolution are demonstrated. Signatures of conformational and hydrogen-bonding fluctuations, protein folding, and chemical-exchange processes are discussed.

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.

Fifth-order two-dimensional Raman spectroscopy of liquid water, crystalline ice Ih and amorphous ices: Sensitivity to anharmonic dynamics and local hydrogen bond network structure

The theoretical study of off-resonant fifth-order two-dimensional (2D)-Raman spectroscopy is made to analyze the intermolecular dynamics of liquid and solid water. The 2D-Raman spectroscopy is susceptible to the nonlinear anharmonic dynamics and local hydrogen bond structure in water. It is found that the distinct 2D-Raman response appears as the negative signal near the t(2) axis. The origin of this negative signal for t(2)<15 fs is from the nonlinear polarizability in the librational motions, whereas that for 30 fs<t(2)<150 fs is attributed to the anharmonic translational motions. It is found that the mechanical anharmonicity and nonlinear polarizability couplings among modes clearly can be observed as the sum- and difference-frequency peaks in the 2D-Raman spectrum (i.e., Fourier transforms of the response). The 2D-Raman spectroscopies of ice Ih and amorphous ices, i.e., low density, high density, and very high density amorphous ices, are also investigated. It is found that the 2D-Raman spectroscopy is very sensitive to the anisotropy of the structure of ice Ih. The strong hydrogen bond stretching band is seen in the 2D-Raman spectroscopy of the polarization directions parallel to the c axis, whereas the contributions of the librational motion can be also seen in the spectrum with the polarization directions parallel to the a axis. The 2D-Raman spectroscopy is also found to be also very sensitive to the differences in local hydrogen bond network structures in various amorphous phases.

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.

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

A practical method to evaluate the contributions of the nonlinear polarizability and anharmonicity of potentials from the experimental and simulation data by using double Fourier transformation is presented. In a Lennard-Jones potential system, an approximated expression of the fifth-order response function using the ratio between nonlinear polarizability and anharmonicity exhibits a good agreement with the results of the molecular dynamics simulation. In a soft-core case, the fifth-order Raman signal indicates that the system consists of the delocalized and localized modes, and only the delocalized mode affects the dramatic change of the fifth-order Raman response functions between solid and liquid phases through nonlinear polarizability.

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.

Why Does the Intermolecular Dynamics of Liquid Biphenyl so Closely Resemble that of Liquid Benzene? Molecular Dynamics Simulation of the Optical-Kerr-Effect Spectra

The combination of optical-Kerr-effect (OKE) spectroscopy and molecular dynamics simulations has provided us with a newfound ability to delve into the librational dynamics of liquids, revealing, in the process, some surprising commonalities among aromatic liquids. Benzene and biphenyl, for example, have remarkably similar OKE spectra despite marked differences in their shapes, sizes, and moments of inertia--and even more chemically distinct aromatics tend to have noticeable similarities in their spectra. We explore this universality by using a molecular dynamics simulation to investigate the librational dynamics of molten biphenyl and to predict its OKE spectrum, comparing the results with our previous calculations for liquid benzene. We suggest that the impressive level of quantitative agreement between these two liquids is largely a reflection of the fact that librations in these and other aromatic liquids act as torsional oscillations with oscillator frequencies selected from the liquid's librational bands. Since these bands are centered about the librational Einstein frequencies, the quantitative similarities between the liquids are essentially reflections of the near identities of their Einstein frequencies. Why then are the Einstein frequencies themselves so insensitive to molecular details? We show that, for nearly planar molecules, mean-square torques and moments of inertia tend to scale with molecular dimensions in much the same way. We demonstrate that this near cancellation provides both a quantitative explanation of the close relationship between benzene and biphenyl and a more general perspective on the similarities seen in the ultrafast dynamics of aromatic liquids.

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.

Qualitative features of the two-dimensional Raman spectrum in liquids

The theory presented earlier [J. Kim and T. Keyes, Phys. Rev. E 66, 051110 (2002)] is analyzed to determine the information available from the two-dimensional Raman spectrum R((5))(t(2),t(1)) in liquids. The known spectra are well represented by the sum of two products of ordinary time correlations predicted by the theory. The shape of R((5)) is related in general to the values of simple same-time averages and concepts amenable to physical intuition. Using standard models for the time correlations entering the theory, specific analytic expressions for the spectrum are obtained depending on two parameters and a time scale, and the behavior of the spectrum is mapped out in the parameter space.

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.

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.