Vibrational spectroscopy via the Caldeira-Leggett model with anharmonic system potentials

General theory of the viscosity of liquids and solids from nonaffine particle motions

A microscopic formula for the viscosity of liquids and solids is derived rigorously from a first-principles (microscopically reversible) Hamiltonian for particle-bath atomistic motion. The derivation is done within the framework of nonaffine linear response theory. This formula may lead to a valid alternative to the Green-Kubo approach to describe the viscosity of condensed matter systems from molecular simulations without having to fit long-time tails. Furthermore, it provides a direct link between the viscosity, the vibrational density of states of the system, and the zero-frequency limit of the memory kernel. Finally, it provides a microscopic solution to Maxwell's interpolation problem of viscoelasticity by naturally recovering Newton's law of viscous flow and Hooke's law of elastic solids in two opposite limits.

Generalized Langevin equation with shear flow and its fluctuation-dissipation theorems derived from a Caldeira-Leggett Hamiltonian

We provide a first-principles derivation of the Langevin equation with shear flow and its corresponding fluctuation-dissipation theorems. Shear flow of simple fluids has been widely investigated by numerical simulations. Most studies postulate a Markovian Langevin equation with a simple shear drag term in the manner of Stokes. However, this choice has never been justified from first principles. We start from a particle-bath system described by a classical Caldeira-Leggett Hamiltonian modified by adding a term proportional to the strain-rate tensor according to Hoover's DOLLS method, and we derive a generalized Langevin equation for the sheared system. We then compute, analytically, the noise time-correlation functions in different regimes. Based on the intensity of the shear rate, we can distinguish between close-to-equilibrium and far-from-equilibrium states. According to the results presented here, the standard, simple, and Markovian form of the Langevin equation with shear flow postulated in the literature is valid only in the limit of extremely weak shear rates compared to the effective vibrational temperature of the bath. For even marginally higher shear rates, the (generalized) Langevin equation is strongly non-Markovian, and nontrivial fluctuation-dissipation theorems are derived.

Discretized hierarchical equations of motion in mixed Liouville-Wigner space for two-dimensional vibrational spectroscopies of liquid water

A model of a bulk water system describing the vibrational motion of intramolecular and intermolecular modes is constructed, enabling analysis of its linear and nonlinear vibrational spectra as well as the energy transfer processes between the vibrational modes. The model is described as a system of four interacting anharmonic oscillators nonlinearly coupled to their respective heat baths. To perform a rigorous numerical investigation of the non-Markovian and nonperturbative quantum dissipative dynamics of the model, we derive discretized hierarchical equations of motion in mixed Liouville-Wigner space, with Lagrange-Hermite mesh discretization being employed in the Liouville space of the intramolecular modes and Lagrange-Hermite mesh discretization and Hermite discretization in the Wigner space of the intermolecular modes. One-dimensional infrared and Raman spectra and two-dimensional terahertz-infrared-visible and infrared-infrared-Raman spectra are computed as demonstrations of the quantum dissipative description provided by our model.

Non-Markovian modeling of protein folding

Protein-folding kinetics is often described as Markovian (i.e., memoryless) diffusion in a one-dimensional free energy landscape, governed by an instantaneous friction coefficient that is fitted to reproduce experimental or simulated folding times. For the α-helix forming polypeptide alanine9 and a specific reaction coordinate that consists of the summed native hydrogen-bond lengths, we demonstrate that the friction extracted from molecular dynamics simulations exhibits significant memory with a decay time that is in the nanosecond range and thus, of the same order as the folding and unfolding times. Our non-Markovian modeling not only reproduces the molecular dynamics simulations accurately but also demonstrates that memory friction effects lead to anomalous and drastically accelerated protein kinetics.

We extract the folding free energy landscape and the time-dependent friction function, the two ingredients of the generalized Langevin equation (GLE), from explicit-water molecular dynamics (MD) simulations of the α-helix forming polypeptide alanine9 for a one-dimensional reaction coordinate based on the sum of the native H-bond distances. Folding and unfolding times from numerical integration of the GLE agree accurately with MD results, which demonstrate the robustness of our GLE-based non-Markovian model. In contrast, Markovian models do not accurately describe the peptide kinetics and in particular, cannot reproduce the folding and unfolding kinetics simultaneously, even if a spatially dependent friction profile is used. Analysis of the GLE demonstrates that memory effects in the friction significantly speed up peptide folding and unfolding kinetics, as predicted by the Grote–Hynes theory, and are the cause of anomalous diffusion in configuration space. Our methods are applicable to any reaction coordinate and in principle, also to experimental trajectories from single-molecule experiments. Our results demonstrate that a consistent description of protein-folding dynamics must account for memory friction effects.

Caldeira-Leggett model vs ab initio potential: A vibrational spectroscopy test of water solvation

We present a semiclassically approximate quantum treatment of solvation with the purpose of investigating the accuracy of the Caldeira-Leggett model. We do that by simulating the vibrational features of water solvation by means of two different approaches. One is entirely based on the adoption of an accurate ab initio potential to describe water clusters of increasing dimensionality. The other one consists of a model made of a central water molecule coupled to a high-dimensional Caldeira-Leggett harmonic bath. We demonstrate the role of quantum effects in the detection of water solvation and show that the computationally cheap approach based on the Caldeira-Leggett bath is only partially effective. The main conclusion of the study is that quantum methods associated with high-level potential energy surfaces are necessary to correctly study solvation features, while simplified models, even if attractive owing to their reduced computational cost, can provide some useful insights but are not able to come up with a comprehensive description of the solvation phenomenon.

Phonon-Induced Optical Dephasing in Single Organic Molecules

We present a joint experiment-theory analysis of the temperature-dependent emission spectra, zero-phonon linewidth, and second-order correlation function of light emitted from a single organic molecule. We observe spectra with a zero-phonon line together with several additional sharp peaks, broad phonon sidebands, and a strongly temperature dependent homogeneous broadening. Our model includes both localized vibrational modes of the molecule and a thermal phonon bath, which we include nonperturbatively, and is able to capture all observed features. For resonant driving we measure Rabi oscillations that become increasingly damped with temperature, which our model naturally reproduces. Our results constitute an essential characterization of the photon coherence of molecules, paving the way to their use in future quantum information applications.

Quantum dynamics and spectroscopy of dihalogens in solid matrices. I. Efficient simulation of the photodynamics of the embedded I2Kr18 cluster using the G-MCTDH method

The molecular dynamics following the electronic BΠu30⁺⟵XΣg+1 photoexcitation of the iodine molecule embedded in solid krypton are studied quantum mechanically using the Gaussian variant of the multiconfigurational time-dependent Hartree method (G-MCTDH). The accuracy of the Gaussian wave packet approximation is validated against numerically exact MCTDH simulations for a fully anharmonic seven-dimensional model of the I₂Kr₁₈ cluster in a crystal Kr cage. The linear absorption spectrum, time-evolving vibrational probability densities, and I₂ energy expectation value are accurately reproduced by the numerically efficient G-MCTDH approach. The reduced density matrix of the chromophore is analyzed in the coordinate, Wigner and energy representations, so as to obtain a multifaceted dynamical view of the guest-host interactions. Vibrational coherences extending over the bond distance range 2.7 Å < R_I-I < 4.0 Å are found to survive for several vibrational periods, despite extensive dissipation. The present results prepare the ground for the simulation of time-resolved coherent Raman spectroscopy of the I₂-krypton system addressed in Paper II.