Two-dimensional Raman and infrared vibrational spectroscopy for a harmonic oscillator system nonlinearly coupled with a colored noise bath

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.

Efficient propagation of the hierarchical equations of motion using the Tucker and hierarchical Tucker tensors

We develop new methods to efficiently propagate the hierarchical equations of motion (HEOM) by using the Tucker and hierarchical Tucker (HT) tensors to represent the reduced density operator and auxiliary density operators. We first show that by employing the split operator method, the specific structure of the HEOM allows a simple propagation scheme using the Tucker tensor. When the number of effective modes in the HEOM increases and the Tucker representation becomes intractable, the split operator method is extended to the binary tree structure of the HT representation. It is found that to update the binary tree nodes related to a specific effective mode, we only need to propagate a short matrix product state constructed from these nodes. Numerical results show that by further employing the mode combination technique commonly used in the multi-configuration time-dependent Hartree approaches, the binary tree representation can be applied to study excitation energy transfer dynamics in a fairly large system including over 10⁴ effective modes. The new methods may thus provide a promising tool in simulating quantum dynamics in condensed phases.

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.

Numerically "exact" approach to open quantum dynamics: The hierarchical equations of motion (HEOM)

An open quantum system refers to a system that is further coupled to a bath system consisting of surrounding radiation fields, atoms, molecules, or proteins. The bath system is typically modeled by an infinite number of harmonic oscillators. This system-bath model can describe the time-irreversible dynamics through which the system evolves toward a thermal equilibrium state at finite temperature. In nuclear magnetic resonance and atomic spectroscopy, dynamics can be studied easily by using simple quantum master equations under the assumption that the system-bath interaction is weak (perturbative approximation) and the bath fluctuations are very fast (Markovian approximation). However, such approximations cannot be applied in chemical physics and biochemical physics problems, where environmental materials are complex and strongly coupled with environments. The hierarchical equations of motion (HEOM) can describe the numerically "exact" dynamics of a reduced system under nonperturbative and non-Markovian system-bath interactions, which has been verified on the basis of exact analytical solutions (non-Markovian tests) with any desired numerical accuracy. The HEOM theory has been used to treat systems of practical interest, in particular, to account for various linear and nonlinear spectra in molecular and solid state materials, to evaluate charge and exciton transfer rates in biological systems, to simulate resonant tunneling and quantum ratchet processes in nanodevices, and to explore quantum entanglement states in quantum information theories. This article presents an overview of the HEOM theory, focusing on its theoretical background and applications, to help further the development of the study of open quantum dynamics.

Hierarchical equations of motion method based on Fano spectrum decomposition for low temperature environments

The hierarchical equations of motion (HEOM) method has become one of the most popular methods for the studies of the open quantum system. However, its applicability to systems at ultra-low temperatures is largely restrained by the enormous computational cost, which is caused by the numerous exponential functions required to accurately characterize the non-Markovian memory of the reservoir environment. To overcome this problem, a Fano spectrum decomposition (FSD) scheme has been proposed recently [Cui et al., J. Chem. Phys. 151, 024110 (2019)], which expands the reservoir correlation functions using polynomial-exponential functions and hence greatly reduces the size of the memory basis set. In this work, we explicitly establish the FSD-based HEOM formalisms for both bosonic and fermionic environments. The accuracy and efficiency of the FSD-based HEOM are exemplified by the calculated low-temperature dissipative dynamics of a spin-boson model and the dynamic and static properties of a single-orbital Anderson impurity model in the Kondo regime. The encouraging numerical results highlight the practicality and usefulness of the FSD-based HEOM method for general open systems at ultra-low temperatures.

Notes on simulating two-dimensional Raman and terahertz-Raman signals with a full molecular dynamics simulation approach

Recent developments in two-dimensional (2D) THz-Raman and 2D Raman spectroscopies have created the possibility for quantitatively investigating the role of many dynamic and structural aspects of the molecular system. We explain the significant points for properly simulating 2D vibrational spectroscopic studies of intermolecular modes using the full molecular dynamics approach, in particular, regarding the system size, the treatment of the thermostat, and inclusion of an Ewald summation for the induced polarizability. Moreover, using the simulation results for water employing various polarization functions, we elucidate the roles of permanent and induced optical properties in determining the 2D profiles of the signal.

Investigating the structural origin of trpzip2 temperature dependent unfolding fluorescence line shape based on a Markov state model simulation

Vibrationally resolved fluorescence spectra of the β-hairpin trpzip2 peptide at two temperatures as well as during a T-jump unfolding process are simulated on the basis of a combination of Markov state models and quantum chemistry schemes. The broad asymmetric spectral line shape feature is reproduced by considering the exciton-phonon couplings. The temperature dependent red shift observed in the experiment has been attributed to the state population changes of specific chromophores. Through further theoretical study, it is found that both the environment's electric field and the chromophores' geometry distortions are responsible for tryptophan fluorescence shift.

Realization of a Brownian engine to study transport phenomena: a semiclassical approach

Brownian particles moving in a periodic potential with or without external load are often used as good theoretical models for the phenomenological studies of microscopic heat engines. The model that we propose here, assumes the particle to be moving in a nonequilibrium medium and we have obtained the exact expression for the stationary current density. We have restricted our consideration to the overdamped motion of the Brownian particle. We present here a self-consistent theory based on the system-reservoir coupling model, within a microscopic approach, of fluctuation induced transport in the semiclassical limit for a general system coupled with two heat baths kept at different temperatures. This essentially puts forth an approach to semiclassical state-dependent diffusion. We also explore the possibility of observing a current when the temperature of the two baths are different, and also envisage that our system may act as a Carnot engine even when the bath temperatures are the same. The condition for such a construction has been elucidated.

Stochastic resonance in a generalized quantum Kubo oscillator

We discuss stochastic resonance in a biased linear quantum system that is subject to multiplicative and additive noises. Starting from a microscopic system-reservoir Hamiltonian, we derive a c-number analogue of the generalized Langevin equation. The developed approach puts forth a quantum mechanical generalization of the "Kubo type" oscillator which is a linear system. Such a system is often used in the literature to study various phenomena in nonequilibrium systems via a particular interaction between system and the external noise. Our analytical results proposed here have the ability to reveal the role of external noise and vis-a-vis the mechanisms and detection of subtle underlying signatures of the stochastic resonance behavior in a linear system. In our development, we show that only when the external noise possesses a "finite correlation time" the quantum effect begins to appear. We observe that the quantum effect enhances the resonance in comparison to the classical one.

Quantum transport in a periodic symmetric potential of a driven quantum system

A system-reservoir nonlinear coupling model has been proposed to the situation when the system is driven externally by a random force and the associated bath is kept in thermal equilibrium, in an attempt to put forth a microscopic approach to quantum state-dependent diffusion and multiplicative noises in terms of a quantum Langevin equation in the overdamped limit (quantum Smoluchowski equation). We then obtain the analytical expression for phase induced quantum current in a periodic potential when the external noise has finite correlation time and explore the dependence of the current on various parameters related to the external noise, for example, the noise strength.

Modeling, calculating, and analyzing multidimensional vibrational spectroscopies

Spectral line shapes in a condensed phase contain information from various dynamic processes that modulate the transition energy, such as microscopic dynamics, inter- and intramolecular couplings, and solvent dynamics. Because nonlinear response functions are sensitive to the complex dynamics of chemical processes, multidimensional vibrational spectroscopies can separate these processes. In multidimensional vibrational spectroscopy, the nonlinear response functions of a molecular dipole or polarizability are measured using ultrashort pulses to monitor inter- and intramolecular vibrational motions. Because a complex profile of such signals depends on the many dynamic and structural aspects of a molecular system, researchers would like to have a theoretical understanding of these phenomena. In this Account, we explore and describe the roles of different physical phenomena that arise from the peculiarities of the system-bath coupling in multidimensional spectra. We also present simple analytical expressions for a weakly coupled multimode Brownian system, which we use to analyze the results obtained by the experiments and simulations. To calculate the nonlinear optical response, researchers commonly use a particular form of a system Hamiltonian fit to the experimental results. The optical responses of molecular vibrational motions have been studied in either an oscillator model or a vibration energy state model. In principle, both models should give the same results as long as the energy states are chosen to be the eigenstates of the oscillator model. The energy state model can provide a simple description of nonlinear optical processes because the diagrammatic Liouville space theory that developed in the electronically resonant spectroscopies can easily handle three or four energy states involved in high-frequency vibrations. However, the energy state model breaks down if we include the thermal excitation and relaxation processes in the dynamics to put the system in a thermal equilibrium state. The roles of these excitation and relaxation processes are different and complicated compared with those in the resonant spectroscopy. Observing the effects of such thermal processes is more intuitive with the oscillator model because the bath modes, which cause the fluctuation and dissipation processes, are also described in the coordinate space. This coordinate space system-bath approach complements a realistic full molecular dynamics simulation approach. By comparing the calculated 2D spectra from the coordinate space model and the energy state model, we can examine the role of thermal processes and anharmonic mode-mode couplings in the energy state model. For this purpose, we employed the Brownian oscillator model with the nonlinear system-bath interaction. Using the hierarchy formalism, we could precisely calculate multidimensional spectra for a single and multimode anharmonic system for inter- and intramolecular vibrational modes.

Phase induced current in presence of nonequilibrium bath: A quantum approach

Based on a system-reservoir nonlinear coupling model, where the associated bath is externally driven by a fluctuating force, we present a microscopic approach to quantum state-dependent diffusion and multiplicative noises in terms of a quantum (Markovian) Langevin equation in overdamped limit when the associated bath is in nonequilibrium state. We then explore the possibility of observing a quantum current when the bath is modulated by white noise, the phenomena which is absent in the classical regime.

Directed motion in a periodic potential of a quantum system coupled to a heat bath driven by a colored noise

A system-reservoir nonlinear coupling model is proposed for a quantum system when the associated bath is not in thermal equilibrium but is modulated by an external colored noise, to present a microscopic approach to quantum state-dependent diffusion and multiplicative noise in terms of a quantum Langevin description. Consequently, the Fokker-Planck equation in position space, valid in the overdamped limit, for multiplicative colored noise is constructed to explore the possibility of observing a quantum current and dependence of the current on various parameters of external noise is examined.

Exploring a free energy landscape by means of multidimensional infrared and terahertz spectroscopies

A model for the dipolar crystal system is employed to explore a role of free energy landscape (FEL), in which dipolar molecules are posted on two-dimensional lattice sites with two-state libratinal dynamics. All dipole-dipole interactions are included to have frustrated interactions among the dipoles. For the regular and distorted lattice cases, the FEL is calculated from the interaction energies and the total polarizations for all possible dipolar states at various temperatures. At high temperatures, the shape of the calculated FEL is smooth and parabolic, while it becomes bumpy at low temperatures exhibiting multiple local minima. To study dynamical aspects of the system, the single flip dynamics and the single-double mixed flips dynamics of dipoles are examined from a master equation approach. As the observables of linear absorption and two-dimensional (2D) infrared, far infrared, and terahertz spectroscopies, the first- and third-order response functions of polarization are calculated for different physical conditions characterized by the FEL. While the linear absorption signals decay in time in a similar manner regardless of the FEL profiles, the 2D signals exhibit prominent differences for those profiles. This indicates that we may differentiate the FEL profiles by changing two-time valuables in 2D spectroscopy. As illustrated in the single-double flips case, the FEL study by means of 2D spectroscopy, however, relies on the dynamics which is set independently from the FEL. The Smoluchowski equation is applied to examine the description of the collective dynamics on the microscopically calculated FEL. We found that the one-dimensional and 2D signals calculated from the Smoluchowski equation agree with those from master equation only at temperatures where the FEL becomes parabolic shape.

Multiplicative cross-correlated noise induced escape rate from a metastable state

We present an analytical framework to study the escape rate from a metastable state under the influence of two external multiplicative cross-correlated noise processes. By starting from a phenomenological stationary Langevin description with multiplicative noise processes, we have investigated the Kramers theory for activated rate processes in a nonequilibrium open system (one dimensional in nature) driven by two external cross-correlated noise processes which are Gaussian, stationary, and delta correlated. Based on the Fokker-Planck description in phase space, we then derive the escape rate from a metastable state in the moderate to large friction limit to study the effect of degree of correlation on the same. By employing numerical simulation in the presence of external cross-correlated additive and multiplicative noises, we check the validity of our analytical formalism for constant dissipation, which shows a satisfactory agreement between both the approaches for the specific choice of noise processes. It is evident both from analytical development and the corresponding numerical simulation that the enhancement of rate is possible by increasing the degree of correlation of the external fluctuations.

Dynamics of a Multimode System Coupled to Multiple Heat Baths Probed by Two-Dimensional Infrared Spectroscopy

Reduced equation of motion for a multimode system coupled to multiple heat baths is constructed by extending the quantum Fokker-Planck equation with low-temperature correction terms (J. Phys. Soc. Jpn. 2005, 74, 3131). Unlike such common approaches used to describe intramolecular multimode vibration as a Bloch-Redfield theory and a stochastic theory, the present formalism is defined by the molecular coordinates. To explore the correlation among different modes through baths, we consider two cases of system-bath couplings. One is a correlated case in which two modes are coupled to a single bath, and the other is an uncorrelated case in which each mode is coupled to a different bath. We further classify the correlated case into two cases, the plus- and minus-correlated cases, according to distinct correlation manners. For these, one-dimensional and two-dimensional infrared (2D-IR) spectra are calculated numerically by solving the equation of motion. It is demonstrated that 2D-IR spectroscopy has the ability to analyze the correlation of fluctuation-dissipation processes among different modes.

Vibrational-exciton couplings for the amide I, II, III, and A modes of peptides

The couplings between all amide fundamentals and their overtones and combination vibrational states are calculated. Combined with the level energies reported previously (Hayashi, T.; Zhuang, W.; Mukamel, S. J. Phys. Chem. A 2005, 109, 9747), we obtain a complete effective vibrational Hamiltonian for the entire amide system. Couplings between neighboring peptide units are obtained using the anharmonic vibrational Hamiltonian of glycine dipeptide (GLDP) at the BPW91/6-31G(d,p) level. Electrostatic couplings between non-neighboring units are calculated by the fourth rank transition multipole coupling (TMC) expansion, including 1/R3 (dipole-dipole), 1/R4 (quadrupole-dipole), and 1/R5 (quadrupole-quadrupole and octapole-dipole) interactions. Exciton delocalization length and its variation with frequency in the various amide bands are calculated. The simulated infrared amide I and II absorptions and CD spectra of 24 residue alpha-helical motifs (SPE3) are in good agreement with experiment.

Theoretical calculations of infrared absorption, vibrational circular dichroism, and two-dimensional vibrational spectra of acetylproline in liquids water and chloroform

Infrared absorption, vibrational circular dichroism, and two-dimensional infrared pump-probe and photon echo spectra of acetylproline solutions are theoretically calculated and directly compared with experiments. In order to quantitatively determine interpeptide interaction-induced amide I mode frequency shifts, high-level quantum chemistry calculations were performed. The solvatochromic amide I mode frequency shift and fluctuation were taken into account by carrying out molecular dynamics simulations of acetylproline dissolved in liquids water and chloroform and by using the extrapolation method developed recently. We first studied correlation time scales of the two amide I vibrational frequency fluctuations, cross correlation between the two fluctuating local mode frequencies, ensemble averaged conformations of the acetylproline molecule in liquids water and chloroform. The corresponding conformations of the acetylproline in liquids water and chloroform are close to the ideal 3(10) helix and the C(7) structure, respectively. A few methods proposed to determine the angle between the two transition dipoles associated with the amide I vibrations were tested and their limitations are discussed.

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.

Modeling vibrational dephasing and energy relaxation of intramolecular anharmonic modes for multidimensional infrared spectroscopies

Starting from a system-bath Hamiltonian in a molecular coordinate representation, we examine an applicability of a stochastic multilevel model for vibrational dephasing and energy relaxation in multidimensional infrared spectroscopy. We consider an intramolecular anharmonic mode nonlinearly coupled to a colored noise bath at finite temperature. The system-bath interaction is assumed linear plus square in the system coordinate, but linear in the bath coordinates. The square-linear system-bath interaction leads to dephasing due to the frequency fluctuation of system vibration, while the linear-linear interaction contributes to energy relaxation and a part of dephasing arises from anharmonicity. To clarify the role and origin of vibrational dephasing and energy relaxation in the stochastic model, the system part is then transformed into an energy eigenstate representation without using the rotating wave approximation. Two-dimensional (2D) infrared spectra are then calculated by solving a low-temperature corrected quantum Fokker-Planck (LTC-QFP) equation for a colored noise bath and by the stochastic theory. In motional narrowing regime, the spectra from the stochastic model are quite different from those from the LTC-QFP. In spectral diffusion regime, however, the 2D line shapes from the stochastic model resemble those from the LTC-QFP besides the blueshifts caused by the dissipation from the colored noise bath. The preconditions for validity of the stochastic theory for molecular vibrational motion are also discussed.

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.