1
|
Reppert M, Reppert D. Equivalence of quantum and classical third order response for weakly anharmonic coupled oscillators. J Chem Phys 2023; 158:114114. [PMID: 36948800 DOI: 10.1063/5.0135260] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/25/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Mike Reppert
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| | - Deborah Reppert
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| |
Collapse
|
2
|
Abstract
Multidimensional optical spectra are measured from the response of a material system to a sequence of laser pulses and have the capacity to elucidate specific molecular interactions and dynamics whose influences are absent or obscured in a conventional linear absorption spectrum. Interpretation of complex spectra is supported by theoretical modeling of the spectroscopic observable, requiring implementation of quantum dynamics for coupled electrons and nuclei. Performing numerically correct quantum dynamics in this context may pose computational challenges, particularly in the condensed phase. Semiclassical methods based on calculating classical trajectories offer a practical alternative. Here I review the recent application of some semiclassical, trajectory-based methods to nonlinear molecular vibrational and electronic spectra. Expected final online publication date for the Annual Review of Physical Chemistry, Volume 73 is April 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.
Collapse
Affiliation(s)
- Roger F. Loring
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York, USA
| |
Collapse
|
3
|
Polley K, Loring RF. Spectroscopic response theory with classical mapping Hamiltonians. J Chem Phys 2020; 153:204103. [PMID: 33261495 DOI: 10.1063/5.0029231] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Exact quantum dynamics with a time-independent Hamiltonian in a discrete state space can be computed using classical mechanics through the classical Meyer-Miller-Stock-Thoss mapping Hamiltonian. In order to compute quantum response functions from classical dynamics, we extend this mapping to a quantum Hamiltonian with time-dependence arising from a classical field. This generalization requires attention to time-ordering in quantum and classical propagators. Quantum response theory with the original quantum Hamiltonian is equivalent to classical response theory with the classical mapping Hamiltonian. We elucidate the structure of classical response theory with the mapping Hamiltonian, thereby generating classical versions of the two-sided quantum density operator diagrams conventionally used to describe spectroscopic processes. This formal development can provide a foundation for new semiclassical approximations to spectroscopic observables for models in which classical nuclear degrees of freedom are introduced into a mapping Hamiltonian describing electronic states. Calculations of the temperature-dependence of two-dimensional electronic spectra for an exciton dimer using two semiclassical approaches are compared with benchmark calculations using the hierarchical equations of motion method.
Collapse
Affiliation(s)
- Kritanjan Polley
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, New York 14853, USA
| | - Roger F Loring
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
4
|
Polley K, Loring RF. Two-dimensional vibronic spectra from classical trajectories. J Chem Phys 2019; 150:164114. [DOI: 10.1063/1.5093911] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Affiliation(s)
- Kritanjan Polley
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, New York 14853, USA
| | - Roger F. Loring
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
5
|
Reppert M, Brumer P. Classical coherent two-dimensional vibrational spectroscopy. J Chem Phys 2018; 148:064101. [DOI: 10.1063/1.5017985] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Mike Reppert
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
| | - Paul Brumer
- Chemical Physics Theory Group, Department of Chemistry, and Center for Quantum Information and Quantum Control, University of Toronto, Toronto, Ontario M5S 3H6, Canada
| |
Collapse
|
6
|
Loring RF. Mean-trajectory approximation for electronic and vibrational-electronic nonlinear spectroscopy. J Chem Phys 2017; 146:144106. [DOI: 10.1063/1.4979621] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Roger F. Loring
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
7
|
Ito H, Jo JY, Tanimura Y. Notes on simulating two-dimensional Raman and terahertz-Raman signals with a full molecular dynamics simulation approach. STRUCTURAL DYNAMICS (MELVILLE, N.Y.) 2015; 2:054102. [PMID: 26798823 PMCID: PMC4711663 DOI: 10.1063/1.4932597] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/29/2015] [Accepted: 09/25/2015] [Indexed: 06/05/2023]
Abstract
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.
Collapse
Affiliation(s)
- Hironobu Ito
- Department of Chemistry, Graduate School of Science, Kyoto University , Sakyoku, Kyoto 606-8502, Japan
| | - Ju-Yeon Jo
- Department of Chemistry, Graduate School of Science, Kyoto University , Sakyoku, Kyoto 606-8502, Japan
| | - Yoshitaka Tanimura
- Department of Chemistry, Graduate School of Science, Kyoto University , Sakyoku, Kyoto 606-8502, Japan
| |
Collapse
|
8
|
Alemi M, Loring RF. Two-Dimensional Vibrational Spectroscopy of a Dissipative System with the Optimized Mean-Trajectory Approximation. J Phys Chem B 2015; 119:8950-9. [PMID: 25275943 PMCID: PMC4383732 DOI: 10.1021/jp5076884] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2014] [Revised: 09/10/2014] [Indexed: 11/30/2022]
Abstract
The optimized mean-trajectory (OMT) approximation is a semiclassical method for computing vibrational response functions from action-quantized classical trajectories connected by discrete transitions representing radiation-matter interactions. Here we apply this method to an anharmonic chromophore coupled to a harmonic bath. A forward-backward trajectory implementation of the OMT method is described that addresses the numerical challenges of applying the OMT to large systems with disparate frequency scales. The OMT is shown to well reproduce line shapes and waiting time dynamics in the pure dephasing limit of weak coupling to an off-resonant bath. The OMT is also shown to describe a case where energy transfer is the predominant source of line broadening.
Collapse
Affiliation(s)
- Mallory Alemi
- Department
of Chemistry and
Chemical Biology, Baker Laboratory, Cornell
University, Ithaca, New York 14853, United
States
| | - Roger F. Loring
- Department
of Chemistry and
Chemical Biology, Baker Laboratory, Cornell
University, Ithaca, New York 14853, United
States
| |
Collapse
|
9
|
Alemi M, Loring RF. Vibrational coherence and energy transfer in two-dimensional spectra with the optimized mean-trajectory approximation. J Chem Phys 2015; 142:212417. [PMID: 26049437 DOI: 10.1063/1.4916644] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
The optimized mean-trajectory (OMT) approximation is a semiclassical method for computing vibrational response functions from action-quantized classical trajectories connected by discrete transitions that represent radiation-matter interactions. Here, we extend the OMT to include additional vibrational coherence and energy transfer processes. This generalized approximation is applied to a pair of anharmonic chromophores coupled to a bath. The resulting 2D spectra are shown to reflect coherence transfer between normal modes.
Collapse
Affiliation(s)
- Mallory Alemi
- Baker Laboratory, Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Roger F Loring
- Baker Laboratory, Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
10
|
Ito H, Hasegawa T, Tanimura Y. Calculating two-dimensional THz-Raman-THz and Raman-THz-THz signals for various molecular liquids: The samplers. J Chem Phys 2014; 141:124503. [PMID: 25273447 DOI: 10.1063/1.4895908] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Hironobu Ito
- Department of Chemistry, Graduate School of Science, Kyoto University, Sakyoku, Kyoto 606-8502, Japan
| | - Taisuke Hasegawa
- Department of Chemistry, Graduate School of Science, Kyoto University, Sakyoku, Kyoto 606-8502, Japan
- Department of Physics, University of Hamburg, Centre for Free Electron Laser Science, DESY, Notkestrasse 85, D-22607 Hamburg, Germany
| | - Yoshitaka Tanimura
- Department of Chemistry, Graduate School of Science, Kyoto University, Sakyoku, Kyoto 606-8502, Japan
| |
Collapse
|
11
|
Gerace M, Loring RF. Two-dimensional spectroscopy of coupled vibrations with the optimized mean-trajectory approximation. J Phys Chem B 2013; 117:15452-61. [PMID: 23924378 PMCID: PMC3865215 DOI: 10.1021/jp405225g] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The optimized mean-trajectory (OMT) approximation is a semiclassical representation of the nonlinear vibrational response function used to compute multidimensional infrared spectra. In this method, response functions are calculated from a sequence of classical trajectories linked by discontinuities representing the effects of radiation-matter interactions, thus providing an approximation to quantum dynamics using classical inputs. This approach was previously formulated and assessed numerically for a single anharmonic degree of freedom. Our previous work is generalized here in two respects. First, the derivation of the OMT is extended to any number of coupled anharmonic vibrations by determining semiclassical approximations for pairs of double-sided Feynman diagrams. Second, an efficient numerical procedure is developed for calculating two-dimensional infrared spectra of coupled anharmonic vibrations in the OMT approximation. The OMT approximation is shown to reproduce the fundamental features of the quantum response function including both coherence and population dynamics.
Collapse
Affiliation(s)
- Mallory Gerace
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, New York, 14853, USA
| | - Roger F. Loring
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, New York, 14853, USA
| |
Collapse
|
12
|
Gerace M, Loring RF. An optimized semiclassical approximation for vibrational response functions. J Chem Phys 2013; 138:124104. [PMID: 23556706 DOI: 10.1063/1.4795941] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Mallory Gerace
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, New York 14853, USA
| | | |
Collapse
|
13
|
Gruenbaum SM, Loring RF. Semiclassical quantization in Liouville space for vibrational dynamics. J Phys Chem B 2011; 115:5148-56. [PMID: 21375227 DOI: 10.1021/jp104872r] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Semiclassical approximations to quantum mechanics can include quantum coherence effects in dynamical calculations based on classical mechanics. The Herman-Kluk (HK) semiclassical propagator has been demonstrated to reproduce quantum effects in nonlinear vibrational response functions of anharmonic oscillators but does not provide a practical numerical route to calculations for multiple degrees of freedom. In an HK calculation of a response function, quantum coherence effects enter through interference between pairs of classical trajectories. We have previously elucidated the mechanism by which the HK approximation reproduces quantum effects in response functions in the regime of quasiperiodic dynamics. We have applied this understanding to significantly simplify the semiclassical calculation of response functions in this dynamical regime. The phase space difference between trajectories is treated perturbatively in anharmonicity, allowing integration over these differences to be performed analytically and leaving integration over mean trajectories to be performed numerically. This mean-trajectory (MT) approximation has been applied to linear and nonlinear vibrational response functions for isolated and coupled anharmonic motions. Here, we derive an MT approximation for the Liouville space time evolution operator or superoperator that propagates the density operator. This analysis provides a form of the MT approximation that is readily applicable to other dynamical quantities besides response functions and clarifies the connection between semiclassical quantization of propagators for the wave function and for the density operator.
Collapse
Affiliation(s)
- Scott M Gruenbaum
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, New York 14853, USA
| | | |
Collapse
|
14
|
Palmieri B, Nagata Y, Mukamel S. Phase-space algorithm for simulating quantum nonlinear response functions of bosons using stochastic classical trajectories. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2010; 82:046706. [PMID: 21230411 PMCID: PMC3709579 DOI: 10.1103/physreve.82.046706] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/24/2010] [Indexed: 11/07/2022]
Abstract
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.
Collapse
Affiliation(s)
- Benoit Palmieri
- Department of Chemistry, University of California, Irvine, California 92697-2025, USA
| | - Yuki Nagata
- Department of Chemistry, University of California, Irvine, California 92697-2025, USA
| | - Shaul Mukamel
- Department of Chemistry, University of California, Irvine, California 92697-2025, USA
| |
Collapse
|
15
|
Sahrapour MM, Makri N. Multitime response functions and nonlinear spectra for model quantum dissipative systems. J Chem Phys 2010; 132:134506. [DOI: 10.1063/1.3336463] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
|
16
|
Gruenbaum SM, Loring RF. Semiclassical nonlinear response functions for coupled anharmonic vibrations. J Chem Phys 2009; 131:204504. [PMID: 19947691 DOI: 10.1063/1.3266566] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Scott M Gruenbaum
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, New York 14853, USA
| | | |
Collapse
|
17
|
Gruenbaum SM, Loring RF. Semiclassical mean-trajectory approximation for nonlinear spectroscopic response functions. J Chem Phys 2008; 129:124510. [DOI: 10.1063/1.2978167] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
18
|
Malinin SV, Chernyak VY. Classical nonlinear response of a chaotic system. I. Collective resonances. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2008; 77:056201. [PMID: 18643136 DOI: 10.1103/physreve.77.056201] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/30/2008] [Indexed: 05/26/2023]
Abstract
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.
Collapse
Affiliation(s)
- Sergey V Malinin
- Department of Chemistry, Wayne State University, 5101 Cass Avenue, Detroit, Michigan 48202, USA
| | | |
Collapse
|
19
|
Malinin SV, Chernyak VY. Classical nonlinear response of a chaotic system. II. Langevin dynamics and spectral decomposition. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2008; 77:056202. [PMID: 18643137 DOI: 10.1103/physreve.77.056202] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/30/2008] [Indexed: 05/26/2023]
Abstract
The spectrum of a strongly chaotic system consists of discrete complex Ruelle-Pollicott (RP) resonances. We interpret the RP resonances as eigenstates and eigenvalues of the Fokker-Planck operator obtained by adding an infinitesimal diffusion term to the first-order Liouville operator. We demonstrate how the deterministic expression for the linear response is reproduced in the limit of vanishing noise. For the second-order response function we establish an equivalence of the spectral decomposition in the limit of vanishing noise and the long-time asymptotic expansion in the deterministic case.
Collapse
Affiliation(s)
- Sergey V Malinin
- Department of Chemistry, Wayne State University, 5101 Cass Avenue, Detroit, Michigan 48202, USA
| | | |
Collapse
|
20
|
Gruenbaum SM, Loring RF. Interference and quantization in semiclassical response functions. J Chem Phys 2008; 128:124106. [DOI: 10.1063/1.2841943] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
21
|
Malinin SV, Chernyak VY. Collective oscillations in the classical nonlinear response of a chaotic system. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2008; 77:025201. [PMID: 18352079 DOI: 10.1103/physreve.77.025201] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/16/2007] [Revised: 09/26/2007] [Indexed: 05/26/2023]
Abstract
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.
Collapse
Affiliation(s)
- Sergey V Malinin
- Department of Chemistry, Wayne State University, 5101 Cass Avenue, Detroit, Michigan 48202, USA
| | | |
Collapse
|
22
|
Noid WG, Loring RF. Classical and quantum mechanical infrared echoes from resonantly coupled molecular vibrations. J Chem Phys 2007; 122:174507. [PMID: 15910045 DOI: 10.1063/1.1888485] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- W G Noid
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, New York 14853, USA
| | | |
Collapse
|
23
|
Saito S, Ohmine I. 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. J Chem Phys 2006; 125:084506. [PMID: 16965028 DOI: 10.1063/1.2232254] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Shinji Saito
- Department of Computational Molecular Science, Institute for Molecular Science, Myodaiji, Okazaki, Aichi 444-8585, Japan.
| | | |
Collapse
|
24
|
Hasegawa T, Tanimura Y. Calculating fifth-order Raman signals for various molecular liquids by equilibrium and nonequilibrium hybrid molecular dynamics simulation algorithms. J Chem Phys 2006; 125:074512. [PMID: 16942356 DOI: 10.1063/1.2217947] [Citation(s) in RCA: 53] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Taisuke Hasegawa
- Department of Chemistry, Graduate School of Science, Kyoto University, Sakyo, Kyoto 606-8502, Japan.
| | | |
Collapse
|
25
|
Goj A, Loring RF. Effect of noise on the classical and quantum mechanical nonlinear response of resonantly coupled anharmonic oscillators. J Chem Phys 2006; 124:194101. [PMID: 16729797 DOI: 10.1063/1.2198203] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Anne Goj
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, NY 14853, USA
| | | |
Collapse
|
26
|
|
27
|
Noid WG, Loring * RF. Nonlinear spectroscopy of resonantly coupled classical mechanical molecular vibrations. Mol Phys 2005. [DOI: 10.1080/00268970500245999] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
28
|
Kryvohuz M, Cao J. Quantum-classical correspondence in response theory. PHYSICAL REVIEW LETTERS 2005; 95:180405. [PMID: 16383881 DOI: 10.1103/physrevlett.95.180405] [Citation(s) in RCA: 43] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/20/2005] [Indexed: 05/05/2023]
Abstract
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.
Collapse
Affiliation(s)
- Maksym Kryvohuz
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | | |
Collapse
|