1
|
Amati G, Runeson JE, Richardson JO. On detailed balance in nonadiabatic dynamics: From spin spheres to equilibrium ellipsoids. J Chem Phys 2023; 158:064113. [PMID: 36792511 DOI: 10.1063/5.0137828] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023] Open
Abstract
Trajectory-based methods that propagate classical nuclei on multiple quantum electronic states are often used to simulate nonadiabatic processes in the condensed phase. A long-standing problem of these methods is their lack of detailed balance, meaning that they do not conserve the equilibrium distribution. In this article, we investigate ideas for restoring detailed balance in mixed quantum-classical systems by tailoring the previously proposed spin-mapping approach to thermal equilibrium. We find that adapting the spin magnitude can recover the correct long-time populations but is insufficient to conserve the full equilibrium distribution. The latter can however be achieved by a more flexible mapping of the spin onto an ellipsoid, which is constructed to fulfill detailed balance for arbitrary potentials. This ellipsoid approach solves the problem of negative populations that has plagued previous mapping approaches and can therefore be applied also to strongly asymmetric and anharmonic systems. Because it conserves the thermal distribution, the method can also exploit efficient sampling schemes used in standard molecular dynamics, which drastically reduces the number of trajectories needed for convergence. The dynamics does however still have mean-field character, as is observed most clearly by evaluating reaction rates in the golden-rule limit. This implies that although the ellipsoid mapping provides a rigorous framework, further work is required to find an accurate classical-trajectory approximation that captures more properties of the true quantum dynamics.
Collapse
Affiliation(s)
- Graziano Amati
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | - Johan E Runeson
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | | |
Collapse
|
2
|
Roosta S, Galami F, Elstner M, Xie W. Efficient Surface Hopping Approach for Modeling Charge Transport in Organic Semiconductors. J Chem Theory Comput 2022; 18:1264-1274. [PMID: 35179894 DOI: 10.1021/acs.jctc.1c00944] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The trajectory surface hopping (TSH) method is nowadays widely applied to study the charge/exciton transport process in organic semiconductors (OSCs). In the present study, we systematically examine the performance of two approximations in the fewest switched surface hopping (FSSH) simulations for charge transport (CT) in several representative OSCs. These approximations include (i) the substitution of the nuclear velocity scaling along the nonadiabatic coupling vector (NCV) by rescaling the hopping probability with the Boltzmann factor (Boltzmann correction (BC)) and (ii) a phenomenological approach to treat the quantum feedback from the electronic system to the nuclear system (implicit charge relaxation (IR)) in the OSCs. We find that charge mobilities computed by FSSH-BC-IR are in very good agreement with the mobilities obtained by standard FSSH simulations with explicit charge relaxation (FSSH-ER), however, at reduced computational cost. A key parameter determining the charge carrier mobility is the reorganization energy, which is sensitively dependent on DFT functionals applied. By employing the IR approximation, the FSSH method allows systematic investigation of the effect of the reorganization energies obtained by different DFT functionals like B3LYP or ωB97XD on CT in OSCs. In comparison to the experiments, FSSH-BC-IR using ωB97XD reorganization energy underestimates mobilities in the low-coupling regime, which may indicate the lack of nuclear quantum effects (e.g., zero point energy (ZPE)) in the simulations. The mobilities obtained by FSSH-BC-IR using the B3LYP reorganization energy agree well with experimental values in 3 orders of magnitude. The accidental agreement may be the consequence of the underestimation of the reorganization energy by the B3LYP functional, which compensates for the neglect of nuclear ZPE in the simulations.
Collapse
Affiliation(s)
- Sara Roosta
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Farhad Galami
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Marcus Elstner
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany.,Institute of Biological Interfaces (IBG-2), Karlsruhe Institute of Technology, 76131 Karlsruhe, Germany
| | - Weiwei Xie
- Key Laboratory of Advanced Energy Materials Chemistry (Ministry of Education), Renewable Energy Conversion and Storage Center, College of Chemistry, Nankai University, Tianjin 300071, China
| |
Collapse
|
3
|
Amadei A, Aschi M. Theoretical-Computational Modelling of the Vibrational Relaxation of Small Inorganic Species in Condensed Phase. RESULTS IN CHEMISTRY 2022. [DOI: 10.1016/j.rechem.2022.100305] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022] Open
|
4
|
Nematiaram T, Troisi A. Modeling charge transport in high-mobility molecular semiconductors: Balancing electronic structure and quantum dynamics methods with the help of experiments. J Chem Phys 2020; 152:190902. [DOI: 10.1063/5.0008357] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Affiliation(s)
- Tahereh Nematiaram
- Department of Chemistry and Materials Innovation Factory, University of Liverpool, Liverpool L69 7ZD, United Kingdom
| | - Alessandro Troisi
- Department of Chemistry and Materials Innovation Factory, University of Liverpool, Liverpool L69 7ZD, United Kingdom
| |
Collapse
|
5
|
Xie W, Holub D, Kubař T, Elstner M. Performance of Mixed Quantum-Classical Approaches on Modeling the Crossover from Hopping to Bandlike Charge Transport in Organic Semiconductors. J Chem Theory Comput 2020; 16:2071-2084. [PMID: 32176844 DOI: 10.1021/acs.jctc.9b01271] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
In the present study, several mixed quantum-classical (MQC) methods are applied to on-the-fly nonadiabatic molecular dynamics simulations of hole transport in molecular organic semiconductors (OSCs). The tested MQC methods contain the mean-field Ehrenfest (MFE), trajectory surface hopping (TSH) approaches based on Tully's fewest switches surface hopping (FSSH) and the global flux surface hopping (GFSH), the latter in the diabatic/adiabatic representation, and a Landau-Zener type trajectory surface hopping (LZSH). We also tested several correction schemes which were proposed to identify trivial crossings and to remove unphysical long-range charge transfers due to decoherence corrections. In addition, several cost-effective approaches for the nuclear velocity adjustment after an energy-allowed/energy-forbidden hop are investigated with respect to detailed balance and internal consistency conditions. To model a broad spectrum of OSCs with different charge transport characteristics, we derived from the anthracene structural model the construction of two additional models by uniformly scaling down the electronic couplings by the factors of 0.1 and 0.5. Anthracene shows a bandlike charge transport mechanism, characterized by slightly delocalized charge carriers 'diffusing' through the crystal. For smaller couplings, the mechanism changes to a hopping type, characteristically differing in the charge delocalization and temperature dependence. The MFE and corrected adiabatic TSH approaches are able to quantitatively reproduce the expected behavior, while the diabatic LZSH method fails for large couplings, as do approaches which are based on the hopping of localized charge between neighboring sites. Moreover, we find that while the hole mobility of the anthracene crystal simulated using the celebrated Marcus theory is in good agreement with the experimental value, its agreement has to be regarded as an accident due to the overestimation of the prefactor in the Marcus rate equation.
Collapse
Affiliation(s)
- Weiwei Xie
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Daniel Holub
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Tomáš Kubař
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Marcus Elstner
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany.,Institute of Biological Interfaces (IBG-2), Karlsruhe Institute of Technology, 76131 Karlsruhe, Germany
| |
Collapse
|
6
|
Smith B, Akimov AV. Modeling nonadiabatic dynamics in condensed matter materials: some recent advances and applications. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2020; 32:073001. [PMID: 31661681 DOI: 10.1088/1361-648x/ab5246] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
This review focuses on recent developments in the field of nonadiabatic molecular dynamics (NA-MD), with particular attention given to condensed-matter systems. NA-MD simulations for small molecular systems can be performed using high-level electronic structure (ES) calculations, methods accounting for the quantization of nuclear motion, and using fewer approximations in the dynamical methodology itself. Modeling condensed-matter systems imposes many limitations on various aspects of NA-MD computations, requiring approximations at various levels of theory-from the ES, to the ways in which the coupling of electrons and nuclei are accounted for. Nonetheless, the approximate treatment of NA-MD in condensed-phase materials has gained a spin lately in many applied studies. A number of advancements of the methodology and computational tools have been undertaken, including general-purpose methods, as well as those tailored to nanoscale and condensed matter systems. This review summarizes such methodological and software developments, puts them into the broader context of existing approaches, and highlights some of the challenges that remain to be solved.
Collapse
Affiliation(s)
- Brendan Smith
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260-3000, United States of America
| | | |
Collapse
|
7
|
Smith B, Akimov AV. A comparative analysis of surface hopping acceptance and decoherence algorithms within the neglect of back-reaction approximation. J Chem Phys 2019; 151:124107. [DOI: 10.1063/1.5122770] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Affiliation(s)
- Brendan Smith
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260-3000, USA
| | - Alexey V. Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260-3000, USA
| |
Collapse
|
8
|
Gu B, Franco I. When can quantum decoherence be mimicked by classical noise? J Chem Phys 2019; 151:014109. [DOI: 10.1063/1.5099499] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Bing Gu
- Department of Chemistry, University of Rochester, Rochester, New York 14627, USA
| | - Ignacio Franco
- Department of Chemistry, University of Rochester, Rochester, New York 14627, USA
- Department of Physics, University of Rochester, Rochester, New York 14627, USA
| |
Collapse
|
9
|
Chitsazi R, Wagner AF. Pressure effects on the vibrational and rotational relaxation of vibrationally excited OH (ν, J) in an argon bath. J Chem Phys 2019; 150:114303. [PMID: 30902000 DOI: 10.1063/1.5063923] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Quasi-classical molecular dynamics simulations were used to study the energy relaxation of an initially non-rotating, vibrationally excited (ν = 4) hydroxyl radical (OH) in an Ar bath at 300 K and at high pressures from 50 atm to 400 atm. A Morse oscillator potential represented the OH, and two sets of interaction potentials were used based on whether the Ar-H potential was a Buckingham (Exp6) or a Lennard-Jones (LJ) potential. The vibrational and rotational energies were monitored for 25 000-90 000 ps for Exp6 trajectories and 5000 ps for LJ trajectories. Comparisons to measured vibrational relaxation rates show that Exp6 rates are superior. Simulated initial vibrational relaxation rates are linearly proportional to pressure, implying no effect of high-pressure breakdown in the isolated binary collision approximation. The vibrational decay curves upward from single-exponential decay. A model based on transition rates that exponentially depend on the anharmonic energy gap between vibrational levels fits the vibrational decay well at all pressures, suggesting that anharmonicity is a major cause of the curvature. Due to the competition of vibration-to-rotation energy transfer and bath gas relaxation, the rotational energy overshoots and then relaxes to its thermal value. Approximate models with adjustable rates for this competition successfully reproduced the rotational results. These models show that a large fraction of the vibrational energy loss is initially converted to rotational energy but that fraction decreases rapidly as the vibrational energy content of OH decreases. While simulated rates change dramatically between Exp6 and LJ potentials, the mechanisms remain the same.
Collapse
Affiliation(s)
- Rezvan Chitsazi
- Department of Chemistry, University of Missouri-Columbia, Columbia, Missouri 65211-7600, USA
| | - Albert F Wagner
- Argonne National Laboratory, Chemical Sciences and Engineering Division, Argonne, Illinois 60439, USA
| |
Collapse
|
10
|
Zauleck JPP, Peschel MT, Rott F, Thallmair S, de Vivie-Riedle R. Ultrafast Reactive Quantum Dynamics Coupled to Classical Solvent Dynamics Using an Ehrenfest Approach. J Phys Chem A 2018; 122:2849-2857. [PMID: 29498853 DOI: 10.1021/acs.jpca.7b10372] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
The inclusion of solvent effects in the theoretical analysis of molecular processes becomes increasingly important. Currently, it is not feasible to directly include the solvent on the quantum level. We use an Ehrenfest approach to study the coupled time evolution of quantum dynamically treated solutes and classical solvents system. The classical dynamics of the solvent is coupled to the wavepacket dynamics of the solute and rotational and translational degrees of freedom of the solute are included classically. This allows quantum dynamics simulations for ultrafast processes that are decided by environment interactions without explicit separation of time scales. We show the application to the dissociation of ICN in liquid Ar as a proof of principal system and to the more applied example of uracil in water.
Collapse
Affiliation(s)
- Julius P P Zauleck
- Department of Chemistry , Ludwig-Maximilians-Universität München , D-81377 München , Germany
| | - Martin T Peschel
- Department of Chemistry , Ludwig-Maximilians-Universität München , D-81377 München , Germany
| | - Florian Rott
- Department of Chemistry , Ludwig-Maximilians-Universität München , D-81377 München , Germany
| | - Sebastian Thallmair
- Department of Chemistry , Ludwig-Maximilians-Universität München , D-81377 München , Germany
| | - Regina de Vivie-Riedle
- Department of Chemistry , Ludwig-Maximilians-Universität München , D-81377 München , Germany
| |
Collapse
|
11
|
Amin M, Askerka M, Batista VS, Brudvig GW, Gunner MR. X-ray Free Electron Laser Radiation Damage through the S-State Cycle of the Oxygen-Evolving Complex of Photosystem II. J Phys Chem B 2017; 121:9382-9388. [DOI: 10.1021/acs.jpcb.7b08371] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Affiliation(s)
- Muhamed Amin
- Center
for Photonics and Smart Materials, Zewail City of Science and Technology, Sheikh
Zayed District, 6th of October City, 12588 Giza, Egypt
- Center for Free-Electron Laser Science, Deutsches Elektronen-Synchrotron DESY, Notkestrasse 85, 22607 Hamburg, Germany
| | - Mikhail Askerka
- Department
of Chemistry, Yale University, New Haven, Connecticut 06520-8107, United States
| | - Victor S. Batista
- Department
of Chemistry, Yale University, New Haven, Connecticut 06520-8107, United States
| | - Gary W. Brudvig
- Department
of Chemistry, Yale University, New Haven, Connecticut 06520-8107, United States
| | - M. R. Gunner
- Department
of Physics, City College of New York, New York, New York 10031, United States
| |
Collapse
|
12
|
Kranz JJ, Elstner M. Simulation of Singlet Exciton Diffusion in Bulk Organic Materials. J Chem Theory Comput 2016; 12:4209-21. [DOI: 10.1021/acs.jctc.6b00235] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Julian J. Kranz
- Institute
of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany
| | - Marcus Elstner
- Institute
of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany
| |
Collapse
|
13
|
Nguyen TS, Nanguneri R, Parkhill J. How electronic dynamics with Pauli exclusion produces Fermi-Dirac statistics. J Chem Phys 2015; 142:134113. [DOI: 10.1063/1.4916822] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Triet S. Nguyen
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - Ravindra Nanguneri
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - John Parkhill
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| |
Collapse
|
14
|
Sherman MC, Corcelli SA. Thermal equilibrium properties of surface hopping with an implicit Langevin bath. J Chem Phys 2015; 142:024110. [DOI: 10.1063/1.4905253] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Affiliation(s)
- M. C. Sherman
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - S. A. Corcelli
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| |
Collapse
|
15
|
Bai S, Xie W, Shi Q. A new trajectory branching approximation to propagate the mixed quantum-classical Liouville equation. J Phys Chem A 2014; 118:9262-71. [PMID: 24964189 DOI: 10.1021/jp503522g] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Starting from the mixed quantum-classical Liouville (MQCL) equation, we derive a new trajectory branching method as a modification to the conventional mean field approximation. In the new method, the mean field approximation is used to propagate the mixed quantum-classical dynamics for short times. When the mean field description becomes invalid, new trajectories are added in the simulation by branching the single trajectory into multiple ones. To achieve this, a new set of variables are defined to monitor the deviations of the dynamics on different potential energy surfaces from the reference mean field trajectory, and their equations of motion are derived from the MQCL equation based on the method of first moment expansion. The new method is tested on several one-dimensional two surface problems and is shown to correctly solve the problem of the mean field approximation in several cases.
Collapse
Affiliation(s)
- Shuming Bai
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, Institute of Chemistry, Chinese Academy of Sciences , Zhongguancun, Beijing 100190, China
| | | | | |
Collapse
|
16
|
Talapatra S, Geva E. Vibrational Lifetimes of Cyanide Ion in Aqueous Solution from Molecular Dynamics Simulations: Intermolecular vs Intramolecular Accepting Modes. J Phys Chem B 2014; 118:7395-404. [DOI: 10.1021/jp504460q] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Surma Talapatra
- Department
of Chemistry, University of Michigan, 930 North University Avenue, Ann Arbor, Michigan 48109, United States
| | - Eitan Geva
- Department
of Chemistry, University of Michigan, 930 North University Avenue, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
17
|
Yao Y, Zhao Y. A variational surface hopping algorithm for the sub-Ohmic spin-boson model. J Chem Phys 2013; 139:014102. [DOI: 10.1063/1.4811777] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
18
|
|
19
|
Aghtar M, Liebers J, Strümpfer J, Schulten K, Kleinekathöfer U. Juxtaposing density matrix and classical path-based wave packet dynamics. J Chem Phys 2012; 136:214101. [PMID: 22697524 DOI: 10.1063/1.4723669] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Abstract
In many physical, chemical, and biological systems energy and charge transfer processes are of utmost importance. To determine the influence of the environment on these transport processes, equilibrium molecular dynamics simulations become more and more popular. From these simulations, one usually determines the thermal fluctuations of certain energy gaps, which are then either used to perform ensemble-averaged wave packet simulations, also called Ehrenfest dynamics, or to employ a density matrix approach via spectral densities. These two approaches are analyzed through energy gap fluctuations that are generated to correspond to a predetermined spectral density. Subsequently, density matrix and wave packet simulations are compared through population dynamics and absorption spectra for different parameter regimes. Furthermore, a previously proposed approach to enforce the correct long-time behavior in the wave packet simulations is probed and an improvement is proposed.
Collapse
Affiliation(s)
- Mortaza Aghtar
- School of Engineering and Science, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | | | | | | | | |
Collapse
|
20
|
Alonso JL, Clemente-Gallardo J, Cuchí JC, Echenique P, Falceto F. Ehrenfest dynamics is purity non-preserving: A necessary ingredient for decoherence. J Chem Phys 2012; 137:054106. [DOI: 10.1063/1.4737861] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
21
|
Bastida A, Zúñiga J, Requena A, Miguel B. Molecular dynamics with quantum transitions study of the vibrational relaxation of the HOD bend fundamental in liquid D2O. J Chem Phys 2012; 136:234507. [DOI: 10.1063/1.4729251] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
22
|
Jiang R, Sibert EL. Surface hopping simulation of vibrational predissociation of methanol dimer. J Chem Phys 2012; 136:224104. [DOI: 10.1063/1.4724219] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
23
|
Bastida A, Soler MA, Zúñiga J, Requena A, Kalstein A, Fernández-Alberti S. Hybrid quantum/classical simulations of the vibrational relaxation of the amide I mode of N-methylacetamide in D2O solution. J Phys Chem B 2012; 116:2969-80. [PMID: 22304000 DOI: 10.1021/jp210727u] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Hybrid quantum/classical molecular dynamics (MD) is applied to simulate the vibrational relaxation (VR) of the amide I mode of deuterated N-methylacetamide (NMAD) in aqueous (D(2)O) solution. A novel version of the vibrational molecular dynamics with quantum transitions (MDQT) treatment is developed in which the amide I mode is treated quantum mechanically while the remaining degrees of freedom are treated classically. The instantaneous normal modes of the initially excited NMAD molecule (INM(0)) are used as internal coordinates since they provide a proper initial partition of the system in quantum and classical subsystems. The evolution in time of the energy stored in each individual normal mode is subsequently quantified using the hybrid quantum-classical instantaneous normal modes (INM(t)). The identities of both the INM(0)s and the INM(t)s are tracked using the equilibrium normal modes (ENMs) as templates. The results extracted from the hybrid MDQT simulations show that the quantum treatment of the amide I mode accelerates the whole VR process versus pure classical simulations and gives better agreement with experiments. The relaxation of the amide I mode is found to be essentially an intramolecular vibrational redistribution (IVR) process with little contribution from the solvent, in agreement with previous theoretical and experimental studies. Two well-defined relaxation mechanisms are identified. The faster one accounts for ≈40% of the total vibrational energy that flows through the NMAD molecule and involves the participation of the lowest frequency vibrations as short-life intermediate modes. The second and slower mechanism accounts for the remaining ≈60% of the energy released and is associated to the energy flow through specific mid-range and high-frequency modes.
Collapse
Affiliation(s)
- Adolfo Bastida
- Departamento de Química Física, Universidad de Murcia, 30100 Murcia, Spain.
| | | | | | | | | | | |
Collapse
|
24
|
Bastida A, Soler MA, Zúñiga J, Requena A, Kalstein A, Fernández-Alberti S. Instantaneous normal modes, resonances, and decay channels in the vibrational relaxation of the amide I mode of N-methylacetamide-D in liquid deuterated water. J Chem Phys 2010; 132:224501. [DOI: 10.1063/1.3435212] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
25
|
Applications of Mixed-Quantum/Classical Trajectories to the Study of Nuclear Quantum Effects in Chemical Reactions and Vibrational Relaxation Processes. ADVANCES IN QUANTUM CHEMISTRY 2010. [DOI: 10.1016/s0065-3276(10)59008-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/08/2023]
|
26
|
Bastida A, Zúñiga J, Requena A, Miguel B. Hybrid quantum/classical simulation of the vibrational relaxation of the bend fundamental in liquid water. J Chem Phys 2009; 131:204505. [DOI: 10.1063/1.3266834] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
27
|
Andrade X, Castro A, Zueco D, Alonso JL, Echenique P, Falceto F, Rubio Á. Modified Ehrenfest Formalism for Efficient Large-Scale ab initio Molecular Dynamics. J Chem Theory Comput 2009; 5:728-42. [DOI: 10.1021/ct800518j] [Citation(s) in RCA: 106] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Xavier Andrade
- Nano-bio Spectroscopy Group and European Theoretical Spectroscopy Facility (ETSF), Departamento de Física de Materiales, Universidad del País Vasco UPV/EHU, Centro Mixto CSIC-UPV, and DIPC, Edificio Korta, Av. Tolosa 72, E-20018 San Sebastián, Spain, Institut für Theoretisch Physik, Freie Universität Berlin, Arnimallee, 14, Berlin 14195, Deutschland, Institut für Physik, Universität Augsburg, Universitätsstraβe 1, D-86135 Augsburg, Germany, Instituto de Biocomputación y Física de Sistemas Complejos (BIFI
| | - Alberto Castro
- Nano-bio Spectroscopy Group and European Theoretical Spectroscopy Facility (ETSF), Departamento de Física de Materiales, Universidad del País Vasco UPV/EHU, Centro Mixto CSIC-UPV, and DIPC, Edificio Korta, Av. Tolosa 72, E-20018 San Sebastián, Spain, Institut für Theoretisch Physik, Freie Universität Berlin, Arnimallee, 14, Berlin 14195, Deutschland, Institut für Physik, Universität Augsburg, Universitätsstraβe 1, D-86135 Augsburg, Germany, Instituto de Biocomputación y Física de Sistemas Complejos (BIFI
| | - David Zueco
- Nano-bio Spectroscopy Group and European Theoretical Spectroscopy Facility (ETSF), Departamento de Física de Materiales, Universidad del País Vasco UPV/EHU, Centro Mixto CSIC-UPV, and DIPC, Edificio Korta, Av. Tolosa 72, E-20018 San Sebastián, Spain, Institut für Theoretisch Physik, Freie Universität Berlin, Arnimallee, 14, Berlin 14195, Deutschland, Institut für Physik, Universität Augsburg, Universitätsstraβe 1, D-86135 Augsburg, Germany, Instituto de Biocomputación y Física de Sistemas Complejos (BIFI
| | - J. L. Alonso
- Nano-bio Spectroscopy Group and European Theoretical Spectroscopy Facility (ETSF), Departamento de Física de Materiales, Universidad del País Vasco UPV/EHU, Centro Mixto CSIC-UPV, and DIPC, Edificio Korta, Av. Tolosa 72, E-20018 San Sebastián, Spain, Institut für Theoretisch Physik, Freie Universität Berlin, Arnimallee, 14, Berlin 14195, Deutschland, Institut für Physik, Universität Augsburg, Universitätsstraβe 1, D-86135 Augsburg, Germany, Instituto de Biocomputación y Física de Sistemas Complejos (BIFI
| | - Pablo Echenique
- Nano-bio Spectroscopy Group and European Theoretical Spectroscopy Facility (ETSF), Departamento de Física de Materiales, Universidad del País Vasco UPV/EHU, Centro Mixto CSIC-UPV, and DIPC, Edificio Korta, Av. Tolosa 72, E-20018 San Sebastián, Spain, Institut für Theoretisch Physik, Freie Universität Berlin, Arnimallee, 14, Berlin 14195, Deutschland, Institut für Physik, Universität Augsburg, Universitätsstraβe 1, D-86135 Augsburg, Germany, Instituto de Biocomputación y Física de Sistemas Complejos (BIFI
| | - Fernando Falceto
- Nano-bio Spectroscopy Group and European Theoretical Spectroscopy Facility (ETSF), Departamento de Física de Materiales, Universidad del País Vasco UPV/EHU, Centro Mixto CSIC-UPV, and DIPC, Edificio Korta, Av. Tolosa 72, E-20018 San Sebastián, Spain, Institut für Theoretisch Physik, Freie Universität Berlin, Arnimallee, 14, Berlin 14195, Deutschland, Institut für Physik, Universität Augsburg, Universitätsstraβe 1, D-86135 Augsburg, Germany, Instituto de Biocomputación y Física de Sistemas Complejos (BIFI
| | - Ángel Rubio
- Nano-bio Spectroscopy Group and European Theoretical Spectroscopy Facility (ETSF), Departamento de Física de Materiales, Universidad del País Vasco UPV/EHU, Centro Mixto CSIC-UPV, and DIPC, Edificio Korta, Av. Tolosa 72, E-20018 San Sebastián, Spain, Institut für Theoretisch Physik, Freie Universität Berlin, Arnimallee, 14, Berlin 14195, Deutschland, Institut für Physik, Universität Augsburg, Universitätsstraβe 1, D-86135 Augsburg, Germany, Instituto de Biocomputación y Física de Sistemas Complejos (BIFI
| |
Collapse
|
28
|
|
29
|
Bastida A, Zúñiga J, Requena A, Miguel B. Full quantum vibrational simulation of the relaxation of the cyanide ion in water using the Ehrenfest method with quantum corrections. J Chem Phys 2008; 129:154501. [DOI: 10.1063/1.2992617] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
30
|
Sibert EL, Ramesh SG, Gulmen TS. Vibrational Relaxation of OH and CH Fundamentals of Polar and Nonpolar Molecules in the Condensed Phase. J Phys Chem A 2008; 112:11291-305. [DOI: 10.1021/jp8068442] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Edwin L. Sibert
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin—Madison, Madison, Wisconsin 53706
| | - Sai G. Ramesh
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin—Madison, Madison, Wisconsin 53706
| | - Tolga S. Gulmen
- Department of Chemistry and Theoretical Chemistry Institute, University of Wisconsin—Madison, Madison, Wisconsin 53706
| |
Collapse
|
31
|
Schmidt JR, Parandekar PV, Tully JC. Mixed quantum-classical equilibrium: Surface hopping. J Chem Phys 2008; 129:044104. [DOI: 10.1063/1.2955564] [Citation(s) in RCA: 132] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|