Vibrational Echoes: Dephasing, Rephasing, and the Stability of Classical Trajectories

Equivalence of quantum and classical third order response for weakly anharmonic coupled oscillators

Two-dimensional (2D) infrared (IR) spectra are commonly interpreted using a quantum diagrammatic expansion that describes the changes to the density matrix of quantum systems in response to light-matter interactions. Although classical response functions (based on Newtonian dynamics) have shown promise in computational 2D IR modeling studies, a simple diagrammatic description has so far been lacking. Recently, we introduced a diagrammatic representation for the 2D IR response functions of a single, weakly anharmonic oscillator and showed that the classical and quantum 2D IR response functions for this system are identical. Here, we extend this result to systems with an arbitrary number of bilinearly coupled, weakly anharmonic oscillators. As in the single-oscillator case, quantum and classical response functions are found to be identical in the weakly anharmonic limit or, in experimental terms, when the anharmonicity is small relative to the optical linewidth. The final form of the weakly anharmonic response function is surprisingly simple and offers potential computational advantages for application to large, multi-oscillator systems.

A theory of average response to large jump perturbations

A key feature of the classical Fluctuation Dissipation theorem is its ability to approximate the average response of a dynamical system to a sufficiently small external perturbation from an appropriate time correlation function of the unperturbed dynamics of this system. In the present work, we examine the situation where the state of a nonlinear dynamical system is perturbed by a finitely large, instantaneous external perturbation (jump), for example, the Earth climate perturbed by an extinction level event. Such jump can be either deterministic or stochastic, and in the case of a stochastic jump its randomness can be spatial, or temporal, or both. We show that, even for large instantaneous jumps, the average response of the system can be expressed in the form of a suitable time correlation function of the corresponding unperturbed dynamics. For stochastic jumps, we consider two situations: one where a single spatially random jump of a system state occurs at a predetermined time, and the other where jumps occur randomly in time with small space-time dependent statistical intensity. For all studied configurations, we compute the corresponding average response formulas in the form of suitable time correlation functions of the unperturbed dynamics. Some efficiently computable approximations are derived for practical modeling scenarios.

Echo dephasing and heat capacity from constrained and unconstrained dynamics of triiodothyronine nuclear receptor protein

The objective of this study is to observe the echo feature curves, vibrational dephasing, and heat capacity of a protein-hormone system taking thyroid hormone receptor-beta (THR-β) as an example. Constrained and unconstrained molecular dynamics simulations are performed by implementing the theory of velocity reassignments to probe the phase coherent state in terms of echo pulses. The constrained vibrations are incorporated by adjusting rigid bonds to all hydrogen atoms with an integrator parameter of 2 fs/step in order to reduce the degrees of freedom whereas 1 fs/step is used in the free vibrations of the atomic cluster. The nature of temperature auto-correlation functions changes so that echo feature curves also show a distinct nature in the cases of constrained and unconstrained vibrations. There is a large variation in kinetic temperature and internal potential energy in the echo time zone. The temperature rate of change of internal potential energy is the main contributor to the heat capacity of the native state protein-hormone system. The heat capacity of proteins estimated from this technique is in good agreement with the values from experiments. This study shows that triiodothyronine (T3) hormone makes some differences in heat capacity upon binding to the THR-β ligand binding domain (LBD). The physical properties of unliganded THR-β and T3-bound THR-β LBD in the cases of constrained and unconstrained dynamics are observed distinctly under the effect of anharmonicity on the phase coherent state of normal modes and the dephasing time lies in a range of 0.6-0.8 ps when the systems are perturbed suddenly.

Nonlinear response theory in chemical kinetics

A theory of nonlinear response of chemical kinetics, in which multiple perturbations are used to probe the time evolution of nonlinear chemical systems, is developed. Expressions for nonlinear chemical response functions and susceptibilities, which can serve as multidimensional measures of the kinetic pathways and rates, are derived. A new class of multidimensional measures that combine multiple perturbations and measurements is also introduced. Nonlinear fluctuation-dissipation relations for steady-state chemical systems, which replace operations of concentration measurement and perturbations, are proposed. Several applications to the analysis of complex reaction mechanisms are provided.

An optimized semiclassical approximation for vibrational response functions

The observables of multidimensional infrared spectroscopy may be calculated from nonlinear vibrational response functions. Fully quantum dynamical calculations of vibrational response functions are generally impractical, while completely classical calculations are qualitatively incorrect at long times. These challenges motivate the development of semiclassical approximations to quantum mechanics, which use classical mechanical information to reconstruct quantum effects. The mean-trajectory (MT) approximation is a semiclassical approach to quantum vibrational response functions employing classical trajectories linked by deterministic transitions representing the effects of the radiation-matter interaction. Previous application of the MT approximation to the third-order response function R(3)(t3, t2, t1) demonstrated that the method quantitatively describes the coherence dynamics of the t3 and t1 evolution times, but is qualitatively incorrect for the waiting-time t2 period. Here we develop an optimized version of the MT approximation by elucidating the connection between this semiclassical approach and the double-sided Feynman diagrams (2FD) that represent the quantum response. Establishing the direct connection between 2FD and semiclassical paths motivates a systematic derivation of an optimized MT approximation (OMT). The OMT uses classical mechanical inputs to accurately reproduce quantum dynamics associated with all three propagation times of the third-order vibrational response function.

Phase-space algorithm for simulating quantum nonlinear response functions of bosons using stochastic classical trajectories

Using the positive P-representation of the density matrix, we develop an algorithm for calculating the quantum many-body nonlinear response functions of a system of bosons driven impulsively by external fields. The formalism maps the quantum time evolution of N boson degrees of freedom into a stochastic dynamics of 4N classical degrees of freedom. The first- and the third-order response functions are calculated by propagating the parameters of the P-representation using a set of coupled Langevin equations with multiplicative noise. These parameters serve as classical variables. Two classical ways for computing the response functions are presented. In the nonequilibrium method, an observable is calculated for weak impulsive pulses, and the response functions are obtained by taking its derivatives with respect to the pulse amplitudes. In the alternative, equilibrium simulation, the response functions are expressed in terms of time-correlation functions involving the P-representation parameters and stability matrices representing the perturbation of the trajectories. The stability matrices can be propagated simultaneously with the Langevin equations for the parameters. The formalism is generalized for a many-body boson system coupled to a harmonic bath.

Semiclassical nonlinear response functions for coupled anharmonic vibrations

Observables in linear and nonlinear infrared spectroscopy may be computed from vibrational response functions describing nuclear dynamics on a single electronic surface. We demonstrate that the Herman-Kluk (HK) semiclassical approximation to the quantum propagator yields an accurate representation of quantum coherence effects in linear and nonlinear response functions for coupled anharmonic oscillators. A considerable numerical price is paid for this accuracy; the calculation requires a multidimensional integral over a highly oscillatory integrand that also grows without bound as a function of evolution times. The interference among classical trajectories in the HK approximation produces quantization of good action variables. By treating this interference analytically, we develop a mean-trajectory (MT) approximation that requires only the propagation of classical trajectories linked by transitions in action. The MT approximation accurately reproduces coherence effects in response functions of coupled anharmonic oscillators in a regime in which the observables are strongly influenced by these interactions among vibrations.

Computational vibrational spectroscopy of peptides and proteins in one and two dimensions

Vibrational spectroscopy provides direct information on molecular environment and motions but, its interpretation is often hampered by band broadening. Over the past decade, two-dimensional (2D) vibrational spectroscopy has emerged as a promising technique to overcome a number of difficulties associated with linear spectroscopy and provided significantly detailed information on the structure and dynamics of complex molecules in condensed phases. This Account reviews recently developed computational methods used to simulate 1D and 2D vibrational spectra. The central quantity to calculate in computational spectroscopy is the spectroscopic response function, which is the product of many contributing factors such as vibrational transition energies, transition moments, and their modulations by fluctuating local environment around a solute. Accurate calculations of such linear and nonlinear responses thus require a concerted effort employing a wide range of methods including electronic structure calculation (ESC) and molecular dynamics (MD) simulation. The electronic structure calculation can provide fundamental quantities such as normal-mode frequencies and transition multipole strengths. However, since the treatable system size is limited with this method, classical MD simulation has also been used to account for the dynamics of the solvent environment. To achieve chemical accuracy, these two results are combined to generate time series of fluctuating transition frequencies and transition moments with the distributed multipole analysis, and this particular approach has been known as the hybrid ESC/MD method. For coupled multichromophore systems, vibrational properties of each chromophore such as a peptide are individually calculated by electronic structure methods and the Hessian matrix reconstruction scheme was used to obtain local mode frequencies and couplings of constituting anharmonic oscillators. The spectra thus obtained, especially for biomolecules including polypeptides and proteins, have proven to be reliable and in good agreement with experimental spectra. An alternative to the hybrid method has also been developed, where the classical limit of the vibrational response function was considered. Its main attraction is the capability to obtain the spectra directly from a set of MD trajectories. A novel development along this direction has been achieved by using quantum mechanical/molecular mechanical (QM/MM) force fields for the accurate description of vibrational anharmonicity and chromophore polarization effects. The latter aspects are critical in the 2D case because classical force fields employing harmonic intramolecular potential cannot produce reliable 2D signal. We anticipate that the computational methods presented here will continue to evolve along with experimental advancements and will be of use to further elucidate ultrafast dynamics of chemical and biological systems.

Classical nonlinear response of a chaotic system. I. Collective resonances

We develop a general semiquantitative picture of nonlinear classical response in strongly chaotic systems. In contrast to behavior in integrable or almost integrable systems, the nonlinear classical response in chaotic systems vanishes at long times. The exponential decay of the response functions in the case of strong chaos is attributed to both exponentially decaying and growing elements in the stability matrices. We calculate the linear and second-order response in one of the simplest chaotic systems: free classical motion on a compact surface of constant negative curvature. The response reveals certain features of collective resonances which do not correspond to any periodic classical trajectories. We demonstrate the relevance of the model for the interpretation of spectroscopic experiments.

Collective oscillations in the classical nonlinear response of a chaotic system

We establish a general semiquantitative phase-space picture of the classical nonlinear response in a strongly chaotic system. As opposed to the case of stable dynamics, the response functions decay exponentially at long times. Damped oscillations in response functions are attributed to collective resonances which do not correspond to any periodic classical motions. We calculate analytically the second-order response in a simple chaotic system and demonstrate the relevance of the concept for the interpretation of spectroscopic data.

Classical and quantum mechanical infrared echoes from resonantly coupled molecular vibrations

The nonlinear response function associated with the infrared vibrational echo is calculated for a quantum mechanical model of resonantly coupled, anharmonic oscillators at zero temperature. The classical mechanical response function is determined from the quantum response function by setting variant Planck's over 2pi-->0, permitting the comparison of the effects of resonant vibrational coupling among an arbitrary number of anharmonic oscillators on quantum and classical vibrational echoes. The quantum response function displays a time dependence that reflects both anharmonicity and resonant coupling, while the classical response function depends on anharmonicity only through a time-independent amplitude, and shows a time dependence controlled only by the resonant coupling. In addition, the classical response function grows without bound in time, a phenomenon associated with the nonlinearity of classical mechanics, and absent in quantum mechanics. This unbounded growth was previously identified in the response function for a system without resonant vibrational energy transfer, and is observed to persist in the presence of resonant coupling among vibrations. Quantitative agreement between classical and quantum response functions is limited to a time scale of duration inversely proportional to the anharmonicity.

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.

Effect of noise on the classical and quantum mechanical nonlinear response of resonantly coupled anharmonic oscillators

Multidimensional infrared spectroscopy probes coupled molecular vibrations in complex, condensed phase systems. Recent theoretical studies have focused on the analytic structure of the nonlinear response functions required to calculate experimental observables in a perturbative treatment of the radiation-matter interaction. Classical mechanical nonlinear response functions have been shown to exhibit unbounded growth for anharmonic, integrable systems, as a consequence of the nonlinearity of classical mechanics, a feature that is absent in a quantum mechanical treatment. We explore the analytic structure of the third-order vibrational response function for an exactly solvable quantum mechanical model that includes some of the important and theoretically challenging aspects of realistic models of condensed phase systems: anharmonicity, resonant coupling, fluctuations, and a well-defined classical mechanical limit.

Classical divergence of nonlinear response functions

The time divergence of classical nonlinear response functions reveals the fundamental difficulty of dynamic perturbation based on classical mechanics. The nature of the divergence is established for systems in regular motions using asymptotic decomposition of Fourier integrals. The asymptotic analysis shows that the divergence cannot be removed by phase-space averaging such as the Boltzmann distribution function. The implications of this study are discussed in the context of the conceptual development of quantum-classical correspondence in dynamic response.

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.

Multidimensional infrared spectroscopy for molecular vibrational modes with dipolar interactions, anharmonicity, and nonlinearity of dipole moments and polarizability

We present an analytical expression for the linear and nonlinear infrared spectra of interacting molecular vibrational motions. Each of the molecular modes is explicitly represented by a classical damped oscillator on an anharmonic multidimensional potential-energy surface. The two essential interactions, the dipole-dipole (DD) and the dipole-induced-dipole (DID) interactions, are taken into account, and each dipole moment and polarizability are expanded to nonlinear order with respect to the nuclear vibrational coordinate. Our analytical treatment leads to expressions for the contributions of anharmonicity, DD and DID interactions, and the nonlinearity of dipole moments and polarizability elements to the one-, two-, and three-dimensional spectra as separated terms, which allows us to discuss the relative importance of these respective contributions. We can calculate multidimensional signals for various configurations of molecules interacting through DD and DID interactions for different material parameters over the whole range of frequencies. We demonstrate that contributions from the DD and DID interactions and anharmonicity are separately detectable through the third-order three-dimensional IR spectroscopy, whereas they cannot be distinguished from each other in either the linear or the second-order IR spectroscopies. The possibility of obtaining the intra- or intermolecular structural information from multidimensional spectra is also discussed.

Quantum-classical correspondence in response theory

The correspondence principle between the quantum commutator [A, B] and the classical Poisson brackets iota h{A, B} is examined in the context of response theory. The classical response function is obtained as the leading term of the expansion of the phase space representation of the response function in terms of Weyl-Wigner transformations and is shown to increase without bound at long times as a result of ignoring divergent higher-order contributions. Systematical inclusion of higher-order contributions improves the accuracy of the h expansion at finite times. Resummation of all the higher-order terms establishes the classical-quantum correspondence <v + n|alpha(t)v> <--> alpha n e iota n omega t|Jv + nh/2. The time interval of the validity of the simple classical limit [A(t), B(0)] --> iota h{A(t), B(0)} is estimated for quasiperiodic dynamics and is shown to be inversely proportional to anharmonicity.

Nondivergent classical response functions from uncertainty principle: Quasiperiodic systems

Time-divergence in linear and nonlinear classical response functions can be removed by taking a phase-space average within the quantized uncertainty volume O(hn) around the microcanonical energy surface. For a quasiperiodic system, the replacement of the microcanonical distribution density in the classical response function with the quantized uniform distribution density results in agreement of quantum and classical expressions through Heisenberg's correspondence principle: each matrix element (u/alpha(t)/v) corresponds to the (u-v)th Fourier component of alpha(t) evaluated along the classical trajectory with mean action (Ju+Jv)/2. Numerical calculations for one- and two-dimensional systems show good agreement between quantum and classical results. The generalization to the case of N degrees of freedom is made. Thus, phase-space averaging within the quantized uncertainty volume provides a useful way to establish the classical-quantum correspondence for the linear and nonlinear response functions of a quasiperiodic system.

Interpreting nonlinear vibrational spectroscopy with the classical mechanical analogs of double-sided Feynman diagrams

Observables in coherent, multiple-pulse infrared spectroscopy may be computed from a vibrational nonlinear response function. This response function is conventionally calculated quantum-mechanically, but the challenges in applying quantum mechanics to large, anharmonic systems motivate the examination of classical mechanical vibrational nonlinear response functions. We present an approximate formulation of the classical mechanical third-order vibrational response function for an anharmonic solute oscillator interacting with a harmonic solvent, which establishes a clear connection between classical and quantum mechanical treatments. This formalism permits the identification of the classical mechanical analog of the pure dephasing of a quantum mechanical degree of freedom, and suggests the construction of classical mechanical analogs of the double-sided Feynman diagrams of quantum mechanics, which are widely applied to nonlinear spectroscopy. Application of a rotating wave approximation permits the analytic extraction of signals obeying particular spatial phase matching conditions from a classical-mechanical response function. Calculations of the third-order response function for an anharmonic oscillator coupled to a harmonic solvent are compared to numerically correct classical mechanical results.