Hybrid equilibrium-nonequilibrium molecular dynamics approach for two-dimensional solute-pump/solvent-probe spectroscopy

Two-dimensional infrared-Raman spectroscopy as a probe of water's tetrahedrality

Two-dimensional spectroscopic techniques combining terahertz (THz), infrared (IR), and visible pulses offer a wealth of information about coupling among vibrational modes in molecular liquids, thus providing a promising probe of their local structure. However, the capabilities of these spectroscopies are still largely unexplored due to experimental limitations and inherently weak nonlinear signals. Here, through a combination of equilibrium-nonequilibrium molecular dynamics (MD) and a tailored spectrum decomposition scheme, we identify a relationship between the tetrahedral order of liquid water and its two-dimensional IR-IR-Raman (IIR) spectrum. The structure-spectrum relationship can explain the temperature dependence of the spectral features corresponding to the anharmonic coupling between low-frequency intermolecular and high-frequency intramolecular vibrational modes of water. In light of these results, we propose new experiments and discuss the implications for the study of tetrahedrality of liquid water.

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.

Linear-Response and Nonlinear-Response Formulations of the Instantaneous Marcus Theory for Nonequilibrium Photoinduced Charge Transfer

Instantaneous Marcus theory (IMT) offers a way for capturing the time-dependent charge transfer (CT) rate coefficient in nonequilibrium photoinduced CT processes, where the system was photoexcited from its equilibrated ground state vertically to the excitonic state, followed by an electronic transition to the CT state. As derived from the linearized semiclassical nonequilibrium Fermi's golden rule (LSC NE-FGR), the original IMT requires expensive all-atom nonequilibrium molecular dynamics (NEMD) simulations. In this work, we propose computationally efficient linear-response and nonlinear-response formulations for IMT rate calculations, which only require equilibrium molecular dynamics simulations. The linear- and nonlinear-response IMT methods were tested to predict the transient behavior in the photoinduced CT dynamics of the carotenoid-porphyrin-C₆₀ molecular triad solvated in explicit tetrahydrofuran. Our result demonstrated that the nonlinear-response IMT is in excellent agreement with the benchmark NEMD for all cases investigated here, whereas the linear-response IMT predicts the correct trend for all cases but overestimates the transient CT rate in one case involving a significant nonequilibrium relaxation. This mild breakdown of linear-response IMT is due to neglecting the higher-order terms in the exact nonlinear-response IMT. Taking advantage of time translational symmetry, the linear- and nonlinear-response approaches were demonstrated to be able to reduce the computational cost by 80% and 60% compared with NEMD simulations, respectively. Thus, we highly recommend the readily applicable and accurate nonlinear-response IMT approach for simulating nonequilibrium CT processes in complex molecular systems in the condensed phase.

Solvation Dynamics in Water. 4. On the Initial Regime of Solvation Relaxation

It is shown, by means of numerical and analytic work, that initial molecular momenta play little significant role in the initial fast solvation relaxation that follows electronic excitation of, and charge creation for, a standard model system of a solute in water. Instead, the nonequilibrium dynamics are predominantly described by noninertial "steering" by the torques directly generated by the newly created charge distribution. It is this process that largely overcomes inertia and drives the relaxation dynamics on a time scale of a few tens of femtoseconds in the key initial regime of the dynamics. These results are discussed in the context of commonly employed descriptions such as inertial, Gaussian, and underdamped dynamical behavior.

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.