1
|
Sayer T, Montoya-Castillo A. Generalized quantum master equations can improve the accuracy of semiclassical predictions of multitime correlation functions. J Chem Phys 2024; 161:011101. [PMID: 38949578 DOI: 10.1063/5.0219205] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2024] [Accepted: 06/12/2024] [Indexed: 07/02/2024] Open
Abstract
Multitime quantum correlation functions are central objects in physical science, offering a direct link between the experimental observables and the dynamics of an underlying model. While experiments such as 2D spectroscopy and quantum control can now measure such quantities, the accurate simulation of such responses remains computationally expensive and sometimes impossible, depending on the system's complexity. A natural tool to employ is the generalized quantum master equation (GQME), which can offer computational savings by extending reference dynamics at a comparatively trivial cost. However, dynamical methods that can tackle chemical systems with atomistic resolution, such as those in the semiclassical hierarchy, often suffer from poor accuracy, limiting the credence one might lend to their results. By combining work on the accuracy-boosting formulation of semiclassical memory kernels with recent work on the multitime GQME, here we show for the first time that one can exploit a multitime semiclassical GQME to dramatically improve both the accuracy of coarse mean-field Ehrenfest dynamics and obtain orders of magnitude efficiency gains.
Collapse
Affiliation(s)
- Thomas Sayer
- Department of Chemistry, University of Colorado Boulder, Boulder, Colorado 80309, USA
| | | |
Collapse
|
2
|
Wiethorn ZR, Hunter KE, Zuehlsdorff TJ, Montoya-Castillo A. Beyond the Condon limit: Condensed phase optical spectra from atomistic simulations. J Chem Phys 2023; 159:244114. [PMID: 38153146 DOI: 10.1063/5.0180405] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2023] [Accepted: 12/06/2023] [Indexed: 12/29/2023] Open
Abstract
While dark transitions made bright by molecular motions determine the optoelectronic properties of many materials, simulating such non-Condon effects in condensed phase spectroscopy remains a fundamental challenge. We derive a Gaussian theory to predict and analyze condensed phase optical spectra beyond the Condon limit. Our theory introduces novel quantities that encode how nuclear motions modulate the energy gap and transition dipole of electronic transitions in the form of spectral densities. By formulating the theory through a statistical framework of thermal averages and fluctuations, we circumvent the limitations of widely used microscopically harmonic theories, allowing us to tackle systems with generally anharmonic atomistic interactions and non-Condon fluctuations of arbitrary strength. We show how to calculate these spectral densities using first-principles simulations, capturing realistic molecular interactions and incorporating finite-temperature, disorder, and dynamical effects. Our theory accurately predicts the spectra of systems known to exhibit strong non-Condon effects (phenolate in various solvents) and reveals distinct mechanisms for electronic peak splitting: timescale separation of modes that tune non-Condon effects and spectral interference from correlated energy gap and transition dipole fluctuations. We further introduce analysis tools to identify how intramolecular vibrations, solute-solvent interactions, and environmental polarization effects impact dark transitions. Moreover, we prove an upper bound on the strength of cross correlated energy gap and transition dipole fluctuations, thereby elucidating a simple condition that a system must follow for our theory to accurately predict its spectrum.
Collapse
Affiliation(s)
- Zachary R Wiethorn
- Department of Chemistry, University of Colorado Boulder, Boulder, Colorado 80309, USA
| | - Kye E Hunter
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, USA
| | - Tim J Zuehlsdorff
- Department of Chemistry, Oregon State University, Corvallis, Oregon 97331, USA
| | | |
Collapse
|
3
|
Seliya P, Bonn M, Grechko M. Experimental Access to Mode-Specific Coupling between Quantum Molecular Vibrations and Classical Bath Modes. J Phys Chem Lett 2023; 14:8630-8637. [PMID: 37728562 PMCID: PMC10544034 DOI: 10.1021/acs.jpclett.3c01974] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2023] [Accepted: 09/06/2023] [Indexed: 09/21/2023]
Abstract
The interaction of quantum-mechanical systems with a fluctuating thermal environment (bath) is fundamental to molecular mechanics and energy transport/dissipation. Its complete picture requires mode-specific measurements of this interaction and an understanding of its nature. Here, we present a combined experimental and theoretical study providing detailed insights into the coupling between a high-frequency vibrational two-level system and thermally excited terahertz modes. Experimentally, two-dimensional terahertz-infrared-visible spectroscopy reports directly on the coupling between quantum oscillators represented by CH3 streching vibrations in liquid dimethyl sulfoxide and distinct low-frequency modes. Theoretically, we present a mixed quantum-classical formalism of the sample response to enable the simultaneous quantum description of high-frequency oscillators and a classical description of the bath. We derive the strength and nature of interaction and find different coupling between CH3 stretch and low-frequency modes. This general approach enables quantitative and mode-specific analysis of coupled quantum and classical dynamics in complex chemical systems.
Collapse
Affiliation(s)
- Pankaj Seliya
- Department of Molecular Spectroscopy, Max Planck Institute for Polymer Research, Ackermannweg 10, D-55128 Mainz, Germany
| | - Mischa Bonn
- Department of Molecular Spectroscopy, Max Planck Institute for Polymer Research, Ackermannweg 10, D-55128 Mainz, Germany
| | - Maksim Grechko
- Department of Molecular Spectroscopy, Max Planck Institute for Polymer Research, Ackermannweg 10, D-55128 Mainz, Germany
| |
Collapse
|
4
|
Lai Y, Geva E. Electronic Absorption Spectra from Off-Diagonal Quantum Master Equations. J Chem Phys 2022; 157:104115. [DOI: 10.1063/5.0106888] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Quantum master equations (QMEs) provide a general framework for describing electronic dynamics within a complex molecular system. Off-diagonal QMEs (OD-QMEs) correspond to a family of QMEs that describe the electronic dynamics in the interaction picture based on treating the off-diagonal coupling terms between electronic states as a small perturbation within the framework of second-order perturbation theory. The fact that OD-QMEs are given in terms of the interaction picture makes it non-trivial to obtain Schrodinger picture electronic coherences from them. A key experimental quantity that relies on the ability to obtain accurate Schrodinger picture electronic coherences is the absorption spectrum. In this paper, we propose using a recently introduced procedure for extracting Schrodinger picture electronic coherences from interaction picture inputs to calculate electronic absorption spectra from electronic dynamics generated by OD-QMEs. The accuracy of the absorption spectra obtained in this way is studied in the context of a biexciton benchmark model, by comparing spectra calculated based on time-local and time-nonlocal OD-QMEs to spectra calculated based on a Redfield-type QME and the non-perturbative and quantum-mechanically exact hierarchical equations of motion (HEOM) method.
Collapse
Affiliation(s)
- Yifan Lai
- Chemistry, University of Michigan, United States of America
| | - Eitan Geva
- Department of Chemistry, University of Michigan Department of Chemistry, United States of America
| |
Collapse
|
5
|
Gao X, Lai Y, Geva E. Simulating Absorption Spectra of Multiexcitonic Systems via Quasiclassical Mapping Hamiltonian Methods. J Chem Theory Comput 2020; 16:6465-6480. [DOI: 10.1021/acs.jctc.0c00709] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Xing Gao
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Yifan Lai
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
6
|
Gao X, Geva E. A Nonperturbative Methodology for Simulating Multidimensional Spectra of Multiexcitonic Molecular Systems via Quasiclassical Mapping Hamiltonian Methods. J Chem Theory Comput 2020; 16:6491-6502. [DOI: 10.1021/acs.jctc.0c00843] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Xing Gao
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
7
|
Kananenka AA, Hsieh CY, Cao J, Geva E. Nonadiabatic Dynamics via the Symmetrical Quasi-Classical Method in the Presence of Anharmonicity. J Phys Chem Lett 2018; 9:319-326. [PMID: 29239614 DOI: 10.1021/acs.jpclett.7b03002] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
The symmetrical quasi-classical (SQC) method recently proposed by Miller and Cotton allows one to simulate nonadiabatic dynamics based on an algorithm with classical-like scaling with respect to system size. This is made possible by casting the electronic degrees of freedom in terms of mapping variables that can be propagated in a classical-like manner. While SQC was shown to be rather accurate when applied to benchmark models with harmonic electronic potential energy surfaces, it was also found to become inaccurate and to suffer numerical instabilities when applied to anharmonic systems. In this paper, we propose an extended SQC (E-SQC) methodology for overcoming those discrepancies by describing the anharmonic nuclear modes, which are coupled to the electronic degrees of freedom, in terms of classical-like mapping variables. The accuracy of E-SQC relative to standard SQC is demonstrated on benchmark models with quartic and Morse potential energy surfaces.
Collapse
Affiliation(s)
- Alexei A Kananenka
- Department of Chemistry, University of Michigan , Ann Arbor, Michigan 48109, United States
| | - Chang-Yu Hsieh
- Department of Chemistry, Massachusetts Institute of Technology , Cambridge, Massachusetts 02319, United States
- Singapore-MIT Alliance for Research and Technology (SMART) Center , Singapore 487373
| | - Jianshu Cao
- Department of Chemistry, Massachusetts Institute of Technology , Cambridge, Massachusetts 02319, United States
- Singapore-MIT Alliance for Research and Technology (SMART) Center , Singapore 487373
| | | |
Collapse
|
8
|
Straight SC, Paesani F. Exploring Electrostatic Effects on the Hydrogen Bond Network of Liquid Water through Many-Body Molecular Dynamics. J Phys Chem B 2016; 120:8539-46. [DOI: 10.1021/acs.jpcb.6b02366] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Shelby C. Straight
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, United States
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California, San Diego, La Jolla, California 92093, United States
| |
Collapse
|
9
|
Terenziani F, Painelli A. Two-dimensional electronic-vibrational spectra: modeling correlated electronic and nuclear motion. Phys Chem Chem Phys 2015; 17:13074-81. [DOI: 10.1039/c5cp01485a] [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/28/2022]
Abstract
2D-EV spectra of solvated organic dyes modeled in terms of a reduced set of electronic diabatic states.
Collapse
Affiliation(s)
- F. Terenziani
- Dipartimento di Chimica
- Università di Parma & INSTM UdR Parma
- 43124 Parma
- Italy
| | - A. Painelli
- Dipartimento di Chimica
- Università di Parma & INSTM UdR Parma
- 43124 Parma
- Italy
| |
Collapse
|
10
|
Li J, Li X, Iyengar SS. Vibrational Properties of Hydrogen-Bonded Systems Using the Multireference Generalization to the “On-the-Fly” Electronic Structure within Quantum Wavepacket ab Initio Molecular Dynamics (QWAIMD). J Chem Theory Comput 2014; 10:2265-80. [DOI: 10.1021/ct5002347] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Junjie Li
- Department of Chemistry and
Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Xiaohu Li
- Department of Chemistry and
Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Srinivasan S. Iyengar
- Department of Chemistry and
Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|
11
|
Tomkins J, Hanna G. Signatures of nanoconfinement on the linear and nonlinear vibrational spectroscopy of a model hydrogen-bonded complex dissolved in a polar solvent. J Phys Chem B 2013; 117:13619-30. [PMID: 24079369 DOI: 10.1021/jp407469f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The one-dimensional IR (1D-IR) absorption and IR pump-probe spectra of a hydrogen stretch in a model hydrogen-bonded complex dissolved in a polar solvent confined in spherical hydrophobic cavities of different sizes were simulated using ground-state mixed quantum-classical dynamics. Due to a thorough analysis of key properties of the complex and solvent from equilibrium trajectory data, we were able to gain insight into the microscopic details underlying the spectra. Both the 1D-IR and IR pump-probe spectra manifested the effects of confinement on the relative stabilities of the covalent and ionic forms of the complex through pronounced changes in their peak intensities and numbers. However, in contrast to the 1D-IR spectra, the time-resolved pump-probe spectra were found to be uniquely sensitive to the changes in the molecular dynamics as the cavity size is varied. In particular, it was found that the variations in the time evolutions of the peak intensities in the pump-probe spectra reflect the differences in the solvation dynamics associated with the various forms of the complex in different locations within the cavities. The ability to detect these differences underscores the advantage of using pump-probe spectroscopy for studying nanoconfined systems.
Collapse
Affiliation(s)
- Joseph Tomkins
- Department of Chemistry, University of Alberta , Edmonton, Alberta, Canada T6G 2G2
| | | |
Collapse
|
12
|
Sow CS, Tomkins J, Hanna G. Computational study of the one- and two-dimensional infrared spectra of a proton-transfer mode in a hydrogen-bonded complex dissolved in a polar nanocluster. Chemphyschem 2013; 14:3309-18. [PMID: 23946271 DOI: 10.1002/cphc.201300610] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2013] [Indexed: 11/07/2022]
Abstract
The signatures of nanosolvation on the one- and two-dimensional (1D and 2D) IR spectra of a proton-transfer mode in a hydrogen-bonded complex dissolved in polar solvent molecule nanoclusters of varying size are elucidated by using mixed quantum-classical molecular dynamics simulations. For this particular system, increasing the number of solvent molecules successively from N=7 to N=9 initiates the transition of the system from a cluster state to a bulk-like state. Both the 1D and 2D IR spectra reflect this transition through pronounced changes in their peak intensities and numbers, but the time-resolved 2D IR spectra also manifest spectral features that uniquely identify the onset of the cluster-to-bulk transition. In particular, it is observed that in the 1D IR spectra, the relative intensities of the peaks change such that the number of peaks decreases from three to two as the size of the cluster increases from N=7 to N=9. In the 2D IR spectra, off-diagonal peaks are observed in the N=7 and N=8 cases at zero waiting time, but not in the N=9 case. It is known that there are no off-diagonal peaks in the 2D IR spectrum of the bulk version of this system at zero waiting time, so the disappearance of these peaks is a unique signature of the onset of bulk-like behavior. Through an examination of the trajectories of various properties of the complex and solvent, it is possible to relate the emergence of these off-diagonal peaks to an interplay between the vibrations of the complex and the solvent polarization dynamics.
Collapse
Affiliation(s)
- Chia Shen Sow
- Department of Chemistry, University of Alberta, Edmonton, AB T6G 2G2 (Canada)
| | | | | |
Collapse
|
13
|
Martinez F, Rekik N, Hanna G. Simulation of nonlinear optical signals via approximate solutions of the quantum–classical Liouville equation: Application to the pump–probe spectroscopy of a condensed phase electron transfer reaction. Chem Phys Lett 2013. [DOI: 10.1016/j.cplett.2013.04.018] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
14
|
Kwac K, Geva E. A mixed quantum-classical molecular dynamics study of the hydroxyl stretch in methanol/carbon tetrachloride mixtures: equilibrium hydrogen-bond structure and dynamics at the ground state and the infrared absorption spectrum. J Phys Chem B 2011; 115:9184-94. [PMID: 21675789 DOI: 10.1021/jp204245z] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a mixed quantum-classical molecular dynamics study of the structure and dynamics of the hydroxyl stretch in methanol/carbon tetrachloride mixtures. One of the methanol molecules is tagged, and its hydroxyl stretch is treated quantum-mechanically, while the remaining degrees of freedom are treated classically. The adiabatic Hamiltonian of the quantum-mechanical hydroxyl is diagonalized on-the-fly to obtain the corresponding adiabatic energy levels and wave functions which depend parametrically on the instantaneous configuration of the classical degrees of freedom. The dynamics of the classical degrees of freedom are in turn affected by the quantum-mechanical state of the tagged hydroxyl stretch via the corresponding Hellmann-Feynman forces. The ability of five different force-field combinations to reproduce the experimental absorption infrared spectrum of the hydroxyl stretch is examined for different isotopomers and on a wide range of compositions. It is found that, in addition to accounting for the anharmonic nature of the hydroxyl stretch, one also has to employ polarizable force fields and account for the damping of the polarizability at short distances. The equilibrium ground-state hydrogen-bonding structure and dynamics is analyzed, and its signature on the absorption infrared spectrum of the hydroxyl stretch is investigated in detail. Five different hydroxyl stretch subpopulations are identified and spectrally assigned: monomers (α), hydrogen-bond acceptors (β), hydrogen-bond donors (γ), simultaneous hydrogen-bond donors and acceptors (δ), and simultaneous hydrogen-bond donors and double-acceptors (ε). The fundamental transition frequencies of the α and β subpopulations are found to be narrowly distributed and to overlap, thereby giving rise to a single narrow band whose intensity is significantly diminished by rotational relaxation. The fundamental transition frequency distributions of the γ, δ, and ε subpopulations are found to be broader and to partially overlap, thereby giving rise to a single broad band which is red-shifted relative to the αβ band. The γδε band is also found to be inhomogeneously broadened and unaffected by rotational relaxation. The exchange rates between the different subpopulations and corresponding branching ratios are reported and explained. Finally, nonlinear mapping relations between the hydroxyl transition frequency and bond length and the electric field along the hydroxyl bond axis are established, which can be used to reduce the computational cost of the mixed quantum-classical treatment to that of a purely classical molecular dynamics simulation.
Collapse
Affiliation(s)
- Kijeong Kwac
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109-1055, USA
| | | |
Collapse
|
15
|
Sakurai A, Tanimura Y. Does ℏ Play a Role in Multidimensional Spectroscopy? Reduced Hierarchy Equations of Motion Approach to Molecular Vibrations. J Phys Chem A 2011; 115:4009-22. [DOI: 10.1021/jp1095618] [Citation(s) in RCA: 56] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Atsunori Sakurai
- Department of Chemistry, Graduate School of Science, Kyoto Universiy, Kyoto 606-8502, Japan
| | - Yoshitaka Tanimura
- Department of Chemistry, Graduate School of Science, Kyoto Universiy, Kyoto 606-8502, Japan
| |
Collapse
|
16
|
Hanna G, Geva E. Signature of Nonadiabatic Transitions on the Pump−Probe Infrared Spectra of a Hydrogen-Bonded Complex Dissolved in a Polar Solvent: A Computational Study. J Phys Chem B 2010; 115:5191-200. [DOI: 10.1021/jp1061495] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Gabriel Hanna
- Department of Chemistry, University of Alberta, Edmonton, Alberta, Canada
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109
| |
Collapse
|
17
|
Hanna G, Geva E. Computational study of the signature of hydrogen-bond strength on the infrared spectra of a hydrogen-bonded complex dissolved in a polar liquid. Chem Phys 2010. [DOI: 10.1016/j.chemphys.2010.01.013] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
18
|
McRobbie PL, Hanna G, Shi Q, Geva E. Signatures of nonequilibrium solvation dynamics on multidimensional spectra. Acc Chem Res 2009; 42:1299-309. [PMID: 19552404 DOI: 10.1021/ar800280s] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Abstract
Multidimensional electronic and vibrational spectroscopies have established themselves over the last decade as uniquely detailed probes of intramolecular structure and dynamics. However, these spectroscopies can also provide powerful tools for probing solute-solvent interactions and the solvation dynamics that they give rise to. To this end, it should be noted that multidimensional spectra can be expressed in terms of optical response functions that differ with respect to the chromophore's quantum state during the various time intervals separating light-matter interactions. The dynamics of the photoinactive degrees of freedom during those time intervals (that is, between pulses) is dictated by potential energy surfaces that depend on the corresponding state of the chromophore. One therefore expects the system to hop between potential surfaces in a manner dictated by the optical response functions. Thus, the corresponding spectra should reflect the system's dynamics during the resulting sequence of nonequilibrium solvation processes. However, the interpretation of multidimensional spectra is often based on the assumption that they reflect the equilibrium dynamics of the photoinactive degrees of freedom on the potential surface that corresponds to the chromophore's ground state. In this Account, we present a systematic analysis of the signature of nonequilibrium solvation dynamics on multidimensional spectra and the ability of various computational methods to capture it. The analysis is performed in the context of the following three model systems: (A) a two-state chromophore with shifted harmonic potential surfaces that differ in frequency, (B) a two-state atomic chromophore in an atomic liquid, and (C) the hydrogen stretch of a moderately strong hydrogen-bonded complex in a dipolar liquid. The following computational methods are employed and compared: (1) exact quantum dynamics (model A only), (2) the semiclassical forward-backward initial value representation (FB-IVR) method (models A and B only), (3) the linearized semiclassical (LSC) method (all three models), and (4) the standard ground-state equilibrium dynamics approach (all three models). The results demonstrate how multidimensional spectra can be used to probe nonequilibrium solvation dynamics in real time and with an unprecedented level of detail. We also show that, unlike the standard method, the LSC and FB-IVR methods can accurately capture the signature of solvation dynamics on the spectra. Our results also suggest that LSC and FB-IVR yield similar results in the presence of rapid dephasing, which is typical in complex condensed-phase systems. This observation gives credence to the use of the LSC method for modeling spectra in complex systems for which an exact or even FB-IVR-based calculation is prohibitively expensive.
Collapse
Affiliation(s)
- Porscha L. McRobbie
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109-1055
| | - Gabriel Hanna
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109-1055
| | - Qiang Shi
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, Institute of Chemistry, Chinese Academy of Sciences (CAS), Zhongguancun, Beijing 100190, P. R. China
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109-1055
| |
Collapse
|
19
|
Tanimura Y, Ishizaki A. Modeling, calculating, and analyzing multidimensional vibrational spectroscopies. Acc Chem Res 2009; 42:1270-9. [PMID: 19441802 DOI: 10.1021/ar9000444] [Citation(s) in RCA: 77] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
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.
Collapse
Affiliation(s)
- Yoshitaka Tanimura
- Department of Chemistry, Graduate School of Science, Kyoto University Kitashirakawa, Sakyoku, Kyoto 606-8502, Japan
| | - Akihito Ishizaki
- Department of Chemistry, Graduate School of Science, Kyoto University Kitashirakawa, Sakyoku, Kyoto 606-8502, Japan
| |
Collapse
|
20
|
Roberts ST, Petersen PB, Ramasesha K, Tokmakoff A, Ufimtsev IS, Martinez TJ. Observation of a Zundel-like transition state during proton transfer in aqueous hydroxide solutions. Proc Natl Acad Sci U S A 2009; 106:15154-9. [PMID: 19666493 PMCID: PMC2741221 DOI: 10.1073/pnas.0901571106] [Citation(s) in RCA: 103] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2009] [Indexed: 11/18/2022] Open
Abstract
It is generally accepted that the anomalous diffusion of the aqueous hydroxide ion results from its ability to accept a proton from a neighboring water molecule; yet, many questions exist concerning the mechanism for this process. What is the solvation structure of the hydroxide ion? In what way do water hydrogen bond dynamics influence the transfer of a proton to the ion? We present the results of femtosecond pump-probe and 2D infrared experiments that probe the O-H stretching vibration of a solution of dilute HOD dissolved in NaOD/D(2)O. Upon the addition of NaOD, measured pump-probe transients and 2D IR spectra show a new feature that decays with a 110-fs time scale. The calculation of 2D IR spectra from an empirical valence bond molecular dynamics simulation of a single NaOH molecule in a bath of H(2)O indicates that this fast feature is due to an overtone transition of Zundel-like H(3)O(2)(-) states, wherein a proton is significantly shared between a water molecule and the hydroxide ion. Given the frequency of vibration of shared protons, the observations indicate the shared proton state persists for 2-3 vibrational periods before the proton localizes on a hydroxide. Calculations based on the EVB-MD model argue that the collective electric field in the proton transfer direction is the appropriate coordinate to describe the creation and relaxation of these Zundel-like transition states.
Collapse
Affiliation(s)
- Sean T. Roberts
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139; and
| | - Poul B. Petersen
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139; and
| | - Krupa Ramasesha
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139; and
| | - Andrei Tokmakoff
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA 02139; and
| | - Ivan S. Ufimtsev
- Department of Chemistry, University of Illinois at Urbana–Champaign, Urbana, IL 61801
| | - Todd J. Martinez
- Department of Chemistry, University of Illinois at Urbana–Champaign, Urbana, IL 61801
| |
Collapse
|
21
|
Hanna G, Geva E. Multidimensional Spectra via the Mixed Quantum-Classical Liouville Method: Signatures of Nonequilibrium Dynamics. J Phys Chem B 2009; 113:9278-88. [DOI: 10.1021/jp902797z] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Gabriel Hanna
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109-1055
| | - Eitan Geva
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109-1055
| |
Collapse
|