A time correlation function theory for the fifth order Raman response function with applications to liquid CS2

2D spectroscopies from condensed phase dynamics: Accessing third-order response properties from equilibrium multi-time correlation functions

The third-order response lies at the heart of simulating and interpreting nonlinear spectroscopies ranging from two dimensional infrared (2D-IR) to 2D electronic (2D-ES), and 2D sum frequency generation (2D-SFG). The extra time and frequency dimensions in these spectroscopies provides access to rich information on the electronic and vibrational states present, the coupling between them, and the resulting rates at which they exchange energy that are obscured in linear spectroscopy, particularly for condensed phase systems that usually contain many overlapping features. While the exact quantum expression for the third-order response is well established it is incompatible with the methods that are practical for calculating the atomistic dynamics of large condensed phase systems. These methods, which include both classical mechanics and quantum dynamics methods that retain quantum statistical properties while obeying the symmetries of classical dynamics, such as LSC-IVR, Centroid Molecular Dynamics (CMD) and Ring Polymer Molecular Dynamics (RPMD) naturally provide short-time approximations to the multi-time symmetrized Kubo transformed correlation function. Here, we show how the third-order response can be formulated in terms of equilibrium symmetrized Kubo transformed correlation functions. We demonstrate the utility and accuracy of our approach by showing how it can be used to obtain the third-order response of a series of model systems using both classical dynamics and RPMD. In particular, we show that this approach captures features such as anharmonically induced vertical splittings and peak shifts while providing a physically transparent framework for understanding multidimensional spectroscopies.

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.

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.

Inclusion of nuclear quantum effects for simulations of nonlinear spectroscopy

The computation and interpretation of nonlinear vibrational spectroscopy is of vital importance for understanding a wide range of dynamical processes in molecular systems. Here, we introduce an approach to evaluate multi-time response functions in terms of multi-time double symmetrized Kubo transformed thermal correlation functions. Furthermore, we introduce a multi-time extension of ring polymer molecular dynamics to evaluate these Kubo transforms. Benchmark calculations show that the approximations are useful for short times even for nonlinear operators, providing a consistent improvement over classical simulations of multi-time correlation functions. The introduced methodology thus provides a practical way of including nuclear quantum effects in multi-time response functions of non-linear optical spectroscopy.

A theoretical study of the sum frequency vibrational spectroscopy of the carbon tetrachloride/water interface

Theoretical approximations to the sum frequency vibrational spectroscopy (SFVS) of the carbon tetrachloride/water interface are constructed using the quantum-corrected time correlation functions (TCF) to aid in interpretation of experimental data and to predict novel vibrational modes. Instantaneous normal mode (INM) methods are used to characterize the observed modes leading to the TCF signal, thus providing molecular resolution of the vibrational lineshapes. Detailed comparisons of the theoretical signals are made with those obtained experimentally and show excellent agreement for the spectral peaks in the O-H stretching region of water. An intermolecular mode, unique to the interface, at 848 cm(-1) is also identifiable, similar to the one seen for the water/vapor interface. INM analysis reveals the resonance is due to a wagging mode (hindered rotation) that was previously identified (Perry et al 2005 J. Chem. Phys. 123 144705) as localized on a single water molecule with both hydrogens displaced normal to the interface-generally it is found that the symmetry breaking at the interface leads to hindered translations and rotations at hydrophilic/hydrophobic interfaces that assume finite vibrational frequencies due to anchoring at the aqueous interface. Additionally, examination of the real and imaginary parts of the theoretical SFVS spectra reveal the spectroscopic species attributed the resonances and possible subspecies in the O-H region; these results are consistent with extant experimental data and associated analysis.

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.

A time correlation function theory describing static field enhanced third order optical effects at interfaces

Sum vibrational frequency spectroscopy, a second order optical process, is interface specific in the dipole approximation. At charged interfaces, there exists a static field, and as a direct consequence, the experimentally detected signal is a combination of enhanced second and static field induced third order contributions. There is significant evidence in the literature of the importance/relative magnitude of this third order contribution, but no previous molecularly detailed approach existed to separately calculate the second and third order contributions. Thus, for the first time, a molecularly detailed time correlation function theory is derived here that allows for the second and third order contributions to sum frequency vibrational spectra to be individually determined. Further, a practical, molecular dynamics based, implementation procedure for the derived correlation functions that describe the third order phenomenon is also presented. This approach includes a novel generalization of point atomic polarizability models to calculate the hyperpolarizability of a molecular system. The full system hyperpolarizability appears in the time correlation functions responsible for third order contributions in the presence of a static field.

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.

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.

Theoretical Investigation of the Temperature Dependence of the Fifth-Order Raman Response Function of Fluid and Liquid Xenon

The temperature dependence of the fifth-order Raman response function, R(5)(t1,t2), is calculated for fluid xenon by employing a recently developed time-correlation function (TCF) theory. The TCF theory expresses the two-dimensional (2D) Raman quantum response function in terms of a two-time, computationally tractable, classical TCF. The theory was shown to be in excellent agreement with existing exact classical MD calculations for liquid xenon as well as reproducing line shape characteristics predicted by earlier theoretical work. It is applied here to investigate the temperature dependence of the fifth-order Raman response function in fluid xenon. In general, the characteristic line shapes are preserved over the temperature range investigated (for the reduced temperature points T* = 0.5, 1.0, and 2.0); differences in the signal decay times and a large decline in intensity with decreasing temperature (and associated anharmonicity) are observed. In addition, there are some signature features that were not observed in earlier results for T* = 1. The most dramatic difference in line shape is observed for the polarization condition, xxzzxx, that shows a vibrational echo peak. In contrast, the fully polarized signal changes mainly in magnitude.

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 theoretical description of the polarization dependence of the sum frequency generation spectroscopy of the water/vapor interface

An improved time correlation function (TCF) description of sum frequency generation (SFG) spectroscopy was developed and applied to theoretically describing the spectroscopy of the ambient water/vapor interface. A more general TCF expression than was published previously is presented-it is valid over the entire vibrational spectrum for both the real and imaginary parts of the signal. Computationally, earlier time correlation function approaches were limited to short correlation times that made signal processing challenging. Here, this limitation is overcome, and well-averaged spectra are presented for the three independent polarization conditions that are possible for electronically nonresonant SFG. The theoretical spectra compare quite favorably in shape and relative magnitude to extant experimental results in the O-H stretching region of water for all polarization geometries. The methodological improvements also allow the calculation of intermolecular SFG spectra. While the intermolecular spectrum of bulk water shows relatively little structure, the interfacial spectra (for polarizations that are sensitive to dipole derivatives normal to the interface--SSP and PPP) show a well-defined intermolecular mode at 875 cm(-1) that is comparable in intensity to the rest of the intermolecular structure, and has an intensity that is approximately one-sixth of the magnitude of the intense free O-H stretching peak. Using instantaneous normal mode methods, the resonance is shown to be due to a wagging mode localized on a single water molecule, almost parallel to the interface, with two hydrogens displaced normal to the interface, and the oxygen anchored in the interface. We have also uncovered the origin of another intermolecular mode at 95 cm(-1) for the SSP and PPP spectra, and at 220 cm(-1) for the SPS spectra. These resonances are due to hindered translations perpendicular to the interface for the SSP and PPP spectra, and translations parallel to the interface for the SPS spectra. Further, by examining the real and imaginary parts of the SFG signal, several resonances are shown to be due to a single spectroscopic species while the "donor" O-H region is shown to consist of three distinct species-consistent with an earlier experimental analysis.

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.

Tractable theory of nonlinear response and multidimensional nonlinear spectroscopy

Nonlinear spectroscopy provides insights into dynamics, but the response functions required for its interpretation pose a challenge to theorists. We proposed an approach in which the fifth-order response function [R5( t1, t2)] was expressed as a two-time classical time correlation function (TCF). Here, we present TCF theory results for R5( t1, t2) in liquid xenon. Using a first-order dipole-induced dipole polarizability model, the result is compared to an exact numerical calculation showing remarkable agreement. In addition, R5( t1, t2) is calculated using the exactly solved polarizability model, yielding different results and predicting an echo signal.

A time correlation function theory of two-dimensional infrared spectroscopy with applications to liquid water

A theory describing the third-order response function R((3))(t(1),t(2),t(3)), which is associated with two-dimensional infrared (2DIR) spectroscopy, has been developed. R((3)) can be written as sums and differences of four distinct quantum mechanical dipole (multi)time correlation functions (TCF's), each with the same classical limit; the combination of TCF's has a leading contribution of order variant Planck's over 2pi (3) and thus there is no obvious classical limit that can be written in terms of a TCF. In order to calculate the response function in a form amenable to classical mechanical simulation techniques, it is rewritten approximately in terms of a single classical TCF, B(R)(t(1),t(2),t(3))=micro(j)(t(2)+t(1))micro(i)(t(3)+t(2)+t(1))micro(k)(t(1))micro(l)(0), where the subscripts denote the Cartesian dipole directions. The response function is then given, in the frequency domain, as the Fourier transform of a classical TCF multiplied by frequency factors. This classical expression can then further be quantum corrected to approximate the true response function, although for low frequency spectroscopy no correction is needed. In the classical limit, R((3)) becomes the sum of multidimensional time derivatives of B(R)(t(1),t(2),t(3)). To construct the theory, the response function's four TCF's are rewritten in terms of a single TCF: first, two TCF's are eliminated from R((3)) using frequency domain detailed balance relationships, and next, two more are removed by relating the remaining TCF's to each other within a harmonic oscillator approximation; the theory invokes a harmonic approximation only in relating the TCF's and applications of theory involve fully anharmonic, atomistically detailed molecular dynamics (MD). Writing the response function as a single TCF thus yields a form amenable to calculation using classical MD methods along with a suitable spectroscopic model. To demonstrate the theory, the response function is obtained for liquid water with emphasis on the OH stretching portion of the spectrum. This approach to evaluating R((3)) can easily be applied to chemically interesting systems currently being explored experimentally by 2DIR and to help understand the information content of the emerging multidimensional spectroscopy.