1
|
Malpathak S, Ananth N. A Linearized Semiclassical Dynamics Study of the Multiquantum Vibrational Relaxation of NO Scattering from a Au(111) Surface. J Phys Chem Lett 2024; 15:794-801. [PMID: 38232133 DOI: 10.1021/acs.jpclett.3c03041] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2024]
Abstract
The vibrational relaxation of NO molecules scattering from a Au(111) surface has served as the focus of efforts to understand nonadiabatic energy transfer at metal-molecule interfaces. Experimental measurements and previous theoretical efforts suggest that multiquantal NO vibrational energy relaxation occurs via electron-hole pair excitations in the metal. Here, using a linearized semiclassical approach, we accurately predict the vibrational relaxation of NO from the νi = 3 state for different incident translational energies. We also accurately capture the central role of transient electron transfer from the metal to the molecule in mediating the vibrational relaxation process but fall short of quantitatively predicting the full extent of multiquantum relaxation for high incident vibrational excitations (νi = 16).
Collapse
Affiliation(s)
- Shreyas Malpathak
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, New York 14853, United States
| | - Nandini Ananth
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, New York 14853, United States
| |
Collapse
|
2
|
Malpathak S, Ananth N. Non-linear correlation functions and zero-point energy flow in mixed quantum-classical semiclassical dynamics. J Chem Phys 2023; 158:104106. [PMID: 36922136 DOI: 10.1063/5.0133222] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2023] Open
Abstract
Mixed quantum classical (MQC)-initial value representation (IVR) is a recently introduced semiclassical framework that allows for selective quantization of the modes of a complex system. In the quantum limit, MQC reproduces the semiclassical Double Herman-Kluk IVR results, accurately capturing nuclear quantum coherences and conserving zero-point energy. However, in the classical limit, although MQC mimics the Husimi-IVR for real-time correlation functions with linear operators, it is significantly less accurate for non-linear correlation functions with errors even at time zero. Here, we identify the origin of this discrepancy in the MQC formulation and propose a modification. We analytically show that the modified MQC approach is exact for all correlation functions at time zero, and in a study of zero-point energy (ZPE) flow, we numerically demonstrate that it correctly obtains the quantum and classical limits as a function of time. Interestingly, although classical-limit MQC simulations show the expected, unphysical ZPE leakage, we find that it is possible to predict and even modify the direction of ZPE flow through selective quantization of the system, with the quantum-limit modes accepting energy but preserving the minimum quantum mechanically required energy.
Collapse
Affiliation(s)
- Shreyas Malpathak
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, New York 14853, USA
| | - Nandini Ananth
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
3
|
Shiga M. Path integral Brownian chain molecular dynamics: A simple approximation of quantum vibrational dynamics. J Comput Chem 2022; 43:1864-1879. [PMID: 36094104 DOI: 10.1002/jcc.26989] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2022] [Revised: 08/06/2022] [Accepted: 08/09/2022] [Indexed: 11/05/2022]
Abstract
An approximate approach to quantum vibrational dynamics, "Brownian chain molecular dynamics (BCMD)," is proposed to alleviate the chain resonance and curvature problems in the imaginary time-based path integral (PI) simulation. Here the non-centroid velocity is randomized at each step when solving the equation of motion of path integral molecular dynamics. This leads to a combination of the Newton equation and the overdamped Langevin equation for the centroid and non-centroid variables, respectively. BCMD shares the basic properties of other PI approaches such as centroid and ring polymer molecular dynamics: It gives the correct Kubo-transformed correlation function at short times, conserves the time symmetry, has the correct high-temperature/classical limits, gives exactly the position and velocity autocorrelations of harmonic oscillator systems, and does not have the zero-point leakage problem. Numerical tests were done on simple molecular models and liquid water. On-the-fly ab initio BCMD simulations were performed for the protonated water cluster, H 5 O 2 + $$ {\mathrm{H}}_5{\mathrm{O}}_2^{+} $$ , and its isotopologue, D 5 O 2 + $$ {\mathrm{D}}_5{\mathrm{O}}_2^{+} $$ .
Collapse
Affiliation(s)
- Motoyuki Shiga
- Center for Computational Science and e-Systems, Japan Atomic Energy Agency, Chiba, Japan
| |
Collapse
|
4
|
Malpathak S, Church MS, Ananth N. A Semiclassical Framework for Mixed Quantum Classical Dynamics. J Phys Chem A 2022; 126:6359-6375. [PMID: 36070472 DOI: 10.1021/acs.jpca.2c03467] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Semiclassical (SC) approximations for quantum dynamic simulations in complex chemical systems range from rigorously accurate methods that are computationally expensive to methods that exhibit near-classical scaling with system size but are limited in their ability to describe quantum effects. In practical studies of high-dimensional reactions, neither extreme is the best choice: frequently a high-level quantum mechanical description is only required for a handful of modes, while the majority of environment modes that do not play a key role in the reactive event of interest are well served with a lower level of theory. In this feature, we introduce modified Filinov filtration as a powerful tool to construct mixed quantum-classical SC theories where different subsystems can be quantized to different extents without introducing ad hoc intersubsystem interaction terms. We demonstrate that these Filinov-based SC methods can systematically tune between quantum and classical limit SC behavior, offering a practical way forward to accurate and computationally efficient simulations of high-dimensional quantum processes.
Collapse
Affiliation(s)
- Shreyas Malpathak
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University Ithaca, New York 14853, United States
| | - Matthew S Church
- Department of Chemistry, Brown University, Providence, Rhode Island 02906, United States
| | - Nandini Ananth
- Department of Chemistry and Chemical Biology, Baker Laboratory, Cornell University Ithaca, New York 14853, United States
| |
Collapse
|
5
|
Poulsen JA, Nyman G. A Divergence-Free Wigner Transform of the Boltzmann Operator Based on an Effective Frequency Theory. J Phys Chem A 2021; 125:9209-9225. [PMID: 34636560 PMCID: PMC8543443 DOI: 10.1021/acs.jpca.1c05860] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
Abstract
![]()
The centroid effective
frequency representation of path integrals
as developed by Feynman and Kleinert was originally aimed at calculating
partition functions and related quantities in the canonical ensemble.
In its path integral formulation, only closed paths
were relevant. This formulation has been used by the present authors
in order to calculate the many-body Wigner function of the Boltzmann
operator, which includes also open paths. This usage of the theory
outside of the original intention can lead to mathematical divergence
issues for potentials with barriers, particularly at low temperature.
In the present paper, we modify the effective frequency theory of
Feynman and Kleinert by also including open paths in its variational
equations. In this way, a divergence-free approximation to the Boltzmann
operator matrix elements is derived. This generalized version of Feynman
and Kleinert’s formulation is thus more robust and can be applied
to all types of barriers at all temperatures. This new version is
used to calculate the Wigner functions of the Boltzmann operator for
a quartic oscillator and for a double well potential and both static
and dynamic properties are studied at several temperatures. The new
theory is found to be essentially as precise as the original one.
Its advantage is that it will always deliver a well-defined, even
if approximate, Wigner function, which can, for instance, be used
for sampling initial conditions for molecular dynamics simulations.
As will be discussed, the theory can be systematically improved by
including higher-order Fourier modes into the nonquadratic part of
the trial action.
Collapse
Affiliation(s)
- Jens Aage Poulsen
- Department of Chemistry and Molecular Biology, University of Gothenburg, SE 405 30 Gothenburg, Sweden
| | - Gunnar Nyman
- Department of Chemistry and Molecular Biology, University of Gothenburg, SE 405 30 Gothenburg, Sweden
| |
Collapse
|
6
|
Dubey V, Dueby S, Daschakraborty S. Breakdown of the Stokes-Einstein relation in supercooled water: the jump-diffusion perspective. Phys Chem Chem Phys 2021; 23:19964-19986. [PMID: 34515269 DOI: 10.1039/d1cp02202d] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Although water is the most ubiquitous liquid it shows many thermodynamic and dynamic anomalies. Some of the anomalies further intensify in the supercooled regime. While many experimental and theoretical studies have focused on the thermodynamic anomalies of supercooled water, fewer studies explored the dynamical anomalies very extensively. This is due to the intricacy of the experimental measurement of the dynamical properties of supercooled water. Violation of the Stokes-Einstein relation (SER), an important relation connecting the diffusion of particles with the viscosity of the medium, is one of the major dynamical anomalies. In absence of experimentally measured viscosity, researchers used to check the validity of SER indirectly using average translational relaxation time or α-relaxation time. Very recently, the viscosity of supercooled water was accurately measured at a wide range of temperatures and pressures. This allowed direct verification of the SER at different temperature-pressure thermodynamic state points. An increasing breakdown of the SER was observed with decreasing temperature. Increasing pressure reduces the extent of breakdown. Although some well-known theories explained the above breakdown, a detailed molecular mechanism was still elusive. Recently, a translational jump-diffusion (TJD) approach has been able to quantitatively explain the breakdown of the SER in pure supercooled water and an aqueous solution of methanol. The objective of this article is to present a detailed and state-of-the-art analysis of the past and present works on the breakdown of SER in supercooled water with a specific focus on the new TJD approach for explaining the breakdown of the SER.
Collapse
Affiliation(s)
- Vikas Dubey
- Department of Chemistry, Indian Institute of Technology Patna, Bihar 801106, India.
| | - Shivam Dueby
- Department of Chemistry, Indian Institute of Technology Patna, Bihar 801106, India.
| | | |
Collapse
|
7
|
Zheng D, Yuan Y, Wang F. Determining the hydration free energies of selected small molecules with MP2 and local MP2 through adaptive force matching. J Chem Phys 2021; 154:104113. [PMID: 33722038 DOI: 10.1063/5.0044712] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Force fields for seven small solute molecules, ethanol, 2-methyl-1-propanol, 2-butanol, cyclohexene, tetrahydropyran, 1,4-dioxane, and 1,4-butanediol, in dilute aqueous solutions were created with the adaptive force matching (AFM) method using MP2 or local MP2 as reference. The force fields provide a way to predict the hydration free energies (HFEs) of these molecules with only electronic structure calculations as reference. For six of the seven molecules, the predicted HFEs are in very good agreement with experiments. For 1,4-butanediol, the model created by force matching LMP2 provides a HFE that is too positive. Further investigation suggests that LMP2 may not be sufficiently accurate for computing HFEs for alcohols with AFM. Other properties, such as enthalpy of hydration, diffusion constants, and vibrational spectra, are also computed with the force field developed. The force fields developed by AFM provide a bridge for computing ensemble properties of the reference electronic structure method. With MP2 and LMP2 as reference methods, the computed properties of the small molecular solutes are found to be in good agreement with experiments.
Collapse
Affiliation(s)
- Dong Zheng
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, USA
| | - Ying Yuan
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, USA
| | - Feng Wang
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, USA
| |
Collapse
|
8
|
Pieroni C, Marsili E, Lauvergnat D, Agostini F. Relaxation dynamics through a conical intersection: Quantum and quantum-classical studies. J Chem Phys 2021; 154:034104. [PMID: 33499611 DOI: 10.1063/5.0036726] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023] Open
Abstract
We study the relaxation process through a conical intersection of a photo-excited retinal chromophore model. The analysis is based on a two-electronic-state two-dimensional Hamiltonian developed by Hahn and Stock [J. Phys. Chem. B 104 1146 (2000)] to reproduce, with a minimal model, the main features of the 11-cis to all-trans isomerization of the retinal of rhodopsin. In particular, we focus on the performance of various trajectory-based schemes to nonadiabatic dynamics, and we compare quantum-classical results to the numerically exact quantum vibronic wavepacket dynamics. The purpose of this work is to investigate, by analyzing electronic and nuclear observables, how the sampling of initial conditions for the trajectories affects the subsequent dynamics.
Collapse
Affiliation(s)
- Carlotta Pieroni
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| | - Emanuele Marsili
- Department of Chemistry, Durham University, South Road, Durham DH1 3LE, United Kingdom
| | - David Lauvergnat
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| | - Federica Agostini
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405 Orsay, France
| |
Collapse
|
9
|
Frank I. Classical Motion of the Nuclei in a Molecule: A Concept Without Alternatives. ChemistrySelect 2020. [DOI: 10.1002/slct.202000066] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Irmgard Frank
- Universität Hannover, Theoretische Chemie Callinstr. 3 A 30167 Hannover
| |
Collapse
|
10
|
Plé T, Huppert S, Finocchi F, Depondt P, Bonella S. Sampling the thermal Wigner density via a generalized Langevin dynamics. J Chem Phys 2019; 151:114114. [PMID: 31542021 DOI: 10.1063/1.5099246] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The Wigner thermal density is a function of considerable interest in the area of approximate (linearized or semiclassical) quantum dynamics where it is employed to generate initial conditions for the propagation of appropriate sets of classical trajectories. In this paper, we propose an original approach to compute the Wigner density based on a generalized Langevin equation. The stochastic dynamics is nontrivial in that it contains a coordinate-dependent friction coefficient and a generalized force that couples momenta and coordinates. These quantities are, in general, not known analytically and have to be estimated via auxiliary calculations. The performance of the new sampling scheme is tested on standard model systems with highly nonclassical features such as relevant zero point energy effects, correlation between momenta and coordinates, and negative parts of the Wigner density. In its current brute force implementation, the algorithm, whose convergence can be systematically checked, is accurate and has only limited overhead compared to schemes with similar characteristics. We briefly discuss potential ways to further improve its numerical efficiency.
Collapse
Affiliation(s)
- Thomas Plé
- Sorbonne Université, CNRS, Institut des NanoSciences de Paris, INSP, 4 Place Jussieu, F-75005 Paris, France
| | - Simon Huppert
- Sorbonne Université, CNRS, Institut des NanoSciences de Paris, INSP, 4 Place Jussieu, F-75005 Paris, France
| | - Fabio Finocchi
- Sorbonne Université, CNRS, Institut des NanoSciences de Paris, INSP, 4 Place Jussieu, F-75005 Paris, France
| | - Philippe Depondt
- Sorbonne Université, CNRS, Institut des NanoSciences de Paris, INSP, 4 Place Jussieu, F-75005 Paris, France
| | - Sara Bonella
- CECAM Centre Européen de Calcul Atomique et Moléculaire, École Polytechnique Fédérale de Lausanne, Batochimie, Avenue Forel 2, 1015 Lausanne, Switzerland
| |
Collapse
|
11
|
Mangaud E, Huppert S, Plé T, Depondt P, Bonella S, Finocchi F. The Fluctuation–Dissipation Theorem as a Diagnosis and Cure for Zero-Point Energy Leakage in Quantum Thermal Bath Simulations. J Chem Theory Comput 2019; 15:2863-2880. [DOI: 10.1021/acs.jctc.8b01164] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Etienne Mangaud
- Sorbonne Université, CNRS - UMR 7588, Institut des NanoSciences de Paris, INSP, 4 Place Jussieu, F-75005 Paris, France
| | - Simon Huppert
- Sorbonne Université, CNRS - UMR 7588, Institut des NanoSciences de Paris, INSP, 4 Place Jussieu, F-75005 Paris, France
| | - Thomas Plé
- Sorbonne Université, CNRS - UMR 7588, Institut des NanoSciences de Paris, INSP, 4 Place Jussieu, F-75005 Paris, France
| | - Philippe Depondt
- Sorbonne Université, CNRS - UMR 7588, Institut des NanoSciences de Paris, INSP, 4 Place Jussieu, F-75005 Paris, France
| | - Sara Bonella
- CECAM Centre Européen de Calcul Atomique et Moléculaire, École Polytechnique Fédérale de Lausanne, Batochimie, Avenue Forel 2, 1015 Lausanne, Switzerland
| | - Fabio Finocchi
- Sorbonne Université, CNRS - UMR 7588, Institut des NanoSciences de Paris, INSP, 4 Place Jussieu, F-75005 Paris, France
| |
Collapse
|
12
|
Affiliation(s)
- Xinzijian Liu
- Beijing National Laboratory For Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing, China
| | - Jian Liu
- Beijing National Laboratory For Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing, China
| |
Collapse
|
13
|
Marsalek O, Markland TE. Quantum Dynamics and Spectroscopy of Ab Initio Liquid Water: The Interplay of Nuclear and Electronic Quantum Effects. J Phys Chem Lett 2017; 8:1545-1551. [PMID: 28296422 DOI: 10.1021/acs.jpclett.7b00391] [Citation(s) in RCA: 137] [Impact Index Per Article: 19.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Understanding the reactivity and spectroscopy of aqueous solutions at the atomistic level is crucial for the elucidation and design of chemical processes. However, the simulation of these systems requires addressing the formidable challenges of treating the quantum nature of both the electrons and nuclei. Exploiting our recently developed methods that provide acceleration by up to 2 orders of magnitude, we combine path integral simulations with on-the-fly evaluation of the electronic structure at the hybrid density functional theory level to capture the interplay between nuclear quantum effects and the electronic surface. Here we show that this combination provides accurate structure and dynamics, including the full infrared and Raman spectra of liquid water. This allows us to demonstrate and explain the failings of lower-level density functionals for dynamics and vibrational spectroscopy when the nuclei are treated quantum mechanically. These insights thus provide a foundation for the reliable investigation of spectroscopy and reactivity in aqueous environments.
Collapse
Affiliation(s)
- Ondrej Marsalek
- Department of Chemistry, Stanford University , Stanford, California 94305, United States
| | - Thomas E Markland
- Department of Chemistry, Stanford University , Stanford, California 94305, United States
| |
Collapse
|
14
|
Montoya-Castillo A, Reichman DR. Approximate but accurate quantum dynamics from the Mori formalism. II. Equilibrium time correlation functions. J Chem Phys 2017; 146:084110. [DOI: 10.1063/1.4975388] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | - David R. Reichman
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| |
Collapse
|
15
|
Montoya-Castillo A, Reichman DR. Path integral approach to the Wigner representation of canonical density operators for discrete systems coupled to harmonic baths. J Chem Phys 2017; 146:024107. [DOI: 10.1063/1.4973646] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | - David R. Reichman
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| |
Collapse
|
16
|
Lee MK, Huo P, Coker DF. Semiclassical Path Integral Dynamics: Photosynthetic Energy Transfer with Realistic Environment Interactions. Annu Rev Phys Chem 2016; 67:639-68. [DOI: 10.1146/annurev-physchem-040215-112252] [Citation(s) in RCA: 96] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Mi Kyung Lee
- Department of Chemistry, Boston University, Boston, Massachusetts 02215;
| | - Pengfei Huo
- Department of Chemistry, University of Rochester, Rochester, New York 14627;
| | - David F. Coker
- Department of Chemistry, Boston University, Boston, Massachusetts 02215;
| |
Collapse
|
17
|
Montoya-Castillo A, Reichman DR. Approximate but accurate quantum dynamics from the Mori formalism: I. Nonequilibrium dynamics. J Chem Phys 2016; 144:184104. [DOI: 10.1063/1.4948408] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Affiliation(s)
| | - David R. Reichman
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| |
Collapse
|
18
|
Liu J, Zhang Z. Path integral Liouville dynamics: Applications to infrared spectra of OH, water, ammonia, and methane. J Chem Phys 2016; 144:034307. [DOI: 10.1063/1.4939953] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jian Liu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
- State Key Joint Laboratory of Environmental Simulation and Pollution Control, College of Environmental Sciences and Engineering, Peking University, Beijing 100871, China
| | - Zhijun Zhang
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| |
Collapse
|
19
|
Beutier J, Vuilleumier R, Bonella S, Ciccotti G. Gas phase infrared spectra from quasi-classical Kubo time correlation functions. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1064550] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- J. Beutier
- Département de Chimie, École Normale Supérieure-PSL Research University , Paris, France
- Département de Chimie, Sorbonne Universités , UPMC Univ Paris 06, Paris, France
- Département de Chimie, CNRS , UMR 8640 PASTEUR, Paris, France
- Dipartimento di Fisica, Università di Roma ‘La Sapienza’ , Roma, Italy
| | - R. Vuilleumier
- Département de Chimie, École Normale Supérieure-PSL Research University , Paris, France
- Département de Chimie, Sorbonne Universités , UPMC Univ Paris 06, Paris, France
- Département de Chimie, CNRS , UMR 8640 PASTEUR, Paris, France
| | - S. Bonella
- Department of Chemistry, CECAM Centre Européen de Calcul Atomique et Moléculaire, École Polytechnique Fédérale de Lausanne , Batochime, Lausanne, Switzerland
| | - G. Ciccotti
- Dipartimento di Fisica, Università di Roma ‘La Sapienza’ , Roma, Italy
- School of Physics, University College of Dublin UCD , Dublin 4, Ireland
| |
Collapse
|
20
|
Abstract
We show a new imaginary time path integral based method--path integral Liouville dynamics (PILD), which can be derived from the equilibrium Liouville dynamics [J. Liu and W. H. Miller, J. Chem. Phys. 134, 104101 (2011)] in the Wigner phase space. Numerical tests of PILD with the simple (white noise) Langevin thermostat have been made for two strongly anharmonic model problems. Since implementation of PILD does not request any specific form of the potential energy surface, the results suggest that PILD offers a potentially useful approach for general condensed phase molecular systems to have the two important properties: conserves the quantum canonical distribution and recovers exact thermal correlation functions (of even nonlinear operators, i.e., nonlinear functions of position or momentum operators) in the classical, high temperature, and harmonic limits.
Collapse
Affiliation(s)
- Jian Liu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| |
Collapse
|
21
|
Spura T, John C, Habershon S, Kühne TD. Nuclear quantum effects in liquid water from path-integral simulations using anab initioforce-matching approach. Mol Phys 2014. [DOI: 10.1080/00268976.2014.981231] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
22
|
Beutier J, Borgis D, Vuilleumier R, Bonella S. Computing thermal Wigner densities with the phase integration method. J Chem Phys 2014; 141:084102. [DOI: 10.1063/1.4892597] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
|
23
|
|
24
|
|
25
|
Habershon S, Manolopoulos DE, Markland TE, Miller TF. Ring-Polymer Molecular Dynamics: Quantum Effects in Chemical Dynamics from Classical Trajectories in an Extended Phase Space. Annu Rev Phys Chem 2013; 64:387-413. [DOI: 10.1146/annurev-physchem-040412-110122] [Citation(s) in RCA: 448] [Impact Index Per Article: 40.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Scott Habershon
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom;
| | - David E. Manolopoulos
- Physical and Theoretical Chemistry Laboratory, University of Oxford, Oxford OX1 3QZ, United Kingdom;
| | - Thomas E. Markland
- Department of Chemistry, Stanford University, Stanford, California 94305;
| | - Thomas F. Miller
- Department of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125;
| |
Collapse
|
26
|
|
27
|
Monteferrante M, Bonella S, Ciccotti G. Quantum dynamical structure factor of liquid neon via a quasiclassical symmetrized method. J Chem Phys 2013; 138:054118. [DOI: 10.1063/1.4789760] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
28
|
Basire M, Borgis D, Vuilleumier R. Computing Wigner distributions and time correlation functions using the quantum thermal bath method: application to proton transfer spectroscopy. Phys Chem Chem Phys 2013; 15:12591-601. [DOI: 10.1039/c3cp50493j] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
29
|
Khaliullin RZ, Kühne TD. Microscopic properties of liquid water from combined ab initio molecular dynamics and energy decomposition studies. Phys Chem Chem Phys 2013; 15:15746-66. [DOI: 10.1039/c3cp51039e] [Citation(s) in RCA: 51] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
|
30
|
Pinnick ER, Erramilli S, Wang F. Predicting the melting temperature of ice-Ih with only electronic structure information as input. J Chem Phys 2012; 137:014510. [PMID: 22779668 DOI: 10.1063/1.4731693] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023] Open
Abstract
The melting temperature of ice-Ih was calculated with only electronic structure information as input by creating a problem-specific force field. The force field, Water model by AFM for Ice and Liquid (WAIL), was developed with the adaptive force matching (AFM) method by fitting to post-Hartree-Fock quality forces obtained in quantum mechanics∕molecular mechanics calculations. WAIL predicts the ice-Ih melting temperature to be 270 K. The model also predicts the densities of ice and water, the temperature of maximum density of water, the heat of vaporizations, and the radial distribution functions for both ice and water in good agreement with experimental measurements. The non-dissociative WAIL model is very similar to a flexible version of the popular TIP4P potential and has comparable computational cost. By customizing to problem-specific configurations with the AFM approach, the resulting model is remarkably more accurate than any variants of TIP4P for simulating ice-Ih and water in the temperature range from 253 K and 293 K under ambient pressure.
Collapse
Affiliation(s)
- Eric R Pinnick
- Department of Physics, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, USA
| | | | | |
Collapse
|
31
|
Liu J, Miller WH, Fanourgakis GS, Xantheas SS, Imoto S, Saito S. Insights in quantum dynamical effects in the infrared spectroscopy of liquid water from a semiclassical study with an ab initio-based flexible and polarizable force field. J Chem Phys 2012; 135:244503. [PMID: 22225165 DOI: 10.1063/1.3670960] [Citation(s) in RCA: 53] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The dynamical properties of liquid water play an important role in many processes in nature. In this paper, we focus on the infrared (IR) absorption spectrum of liquid water based on the linearized semiclassical initial value representation (LSC-IVR) with the local Gaussian approximation (LGA) [J. Liu and W. H. Miller, J. Chem. Phys. 131, 074113 (2009)] and an ab initio based, flexible, polarizable Thole-type model (TTM3-F) [G. S. Fanourgakis and S. S. Xantheas, J. Chem. Phys. 128, 074506 (2008)]. Although the LSC-IVR (LGA) gives the exact result for the isolated three-dimensional shifted harmonic stretching model, it yields a blueshifted peak position for the more realistic anharmonic stretching potential. By using the short-time information of the LSC-IVR correlation function; however, it is shown how one can obtain more accurate results for the position of the stretching peak. Due to the physical decay in the condensed phase system, the LSC-IVR (LGA) is a good and practical approximate quantum approach for the IR spectrum of liquid water. The present results offer valuable insight into future attempts to improve the accuracy of the TTM3-F potential or other ab initio-based models in reproducing the IR spectrum of liquid water.
Collapse
Affiliation(s)
- Jian Liu
- Department of Chemistry, University of California, Berkeley, California 94720-1460, USA
| | | | | | | | | | | |
Collapse
|
32
|
Huo P, Coker DF. Influence of environment induced correlated fluctuations in electronic coupling on coherent excitation energy transfer dynamics in model photosynthetic systems. J Chem Phys 2012; 136:115102. [DOI: 10.1063/1.3693019] [Citation(s) in RCA: 42] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
33
|
Monteferrante M, Bonella S, Ciccotti G. Linearized symmetrized quantum time correlation functions calculation via phase pre-averaging. Mol Phys 2011. [DOI: 10.1080/00268976.2011.619506] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
34
|
Liu J, Alder BJ, Miller WH. A semiclassical study of the thermal conductivity of low temperature liquids. J Chem Phys 2011; 135:114105. [DOI: 10.1063/1.3639107] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
35
|
Paesani F. Hydrogen bond dynamics in heavy water studied with quantum dynamical simulations. Phys Chem Chem Phys 2011; 13:19865-75. [PMID: 21892511 DOI: 10.1039/c1cp21863h] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
The structure and dynamics of the hydrogen-bond network in heavy water (D(2)O) is studied as a function of the temperature using quantum dynamical simulations. Our approach combines an ab initio-based representation of the water interactions with an explicit quantum treatment of the molecular motion. A direct connection between the calculated linear and nonlinear vibrational spectra and the underlying molecular dynamics is made, which provides new insights into the rearrangement of the hydrogen-bond network in heavy water. A comparison with previous calculations on liquid H(2)O suggests that tunneling does not effectively contribute to the dynamics of the water hydrogen-bond network above the melting point. However, the effects of nuclear quantization are not negligible at all temperatures and become increasingly important near the melting point, which is in agreement with recent experimental analysis of the structural properties of liquid water as well as of the proton momentum distribution in supercooled water.
Collapse
Affiliation(s)
- Francesco Paesani
- Department of Chemistry and Biochemistry, University of California, La Jolla, San Diego, CA 92093, USA.
| |
Collapse
|
36
|
Liu J. Two more approaches for generating trajectory-based dynamics which conserves the canonical distribution in the phase space formulation of quantum mechanics. J Chem Phys 2011; 134:194110. [DOI: 10.1063/1.3589406] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
37
|
Liu J, Miller WH. An approach for generating trajectory-based dynamics which conserves the canonical distribution in the phase space formulation of quantum mechanics. II. Thermal correlation functions. J Chem Phys 2011; 134:104102. [DOI: 10.1063/1.3555274] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
38
|
Andersson S, Arasa C, Yabushita A, Yokoyama M, Hama T, Kawasaki M, Western CM, Ashfold MNR. A theoretical and experimental study on translational and internal energies of H2O and OH from the 157 nm irradiation of amorphous solid water at 90 K. Phys Chem Chem Phys 2011; 13:15810-20. [DOI: 10.1039/c1cp21138b] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
39
|
Czakó G, Kaledin AL, Bowman JM. Zero-point energy constrained quasiclassical, classical, and exact quantum simulations of isomerizations and radial distribution functions of the water trimer using an ab initio potential energy surface. Chem Phys Lett 2010. [DOI: 10.1016/j.cplett.2010.10.015] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
40
|
Bonella S, Monteferrante M, Pierleoni C, Ciccotti G. Path integral based calculations of symmetrized time correlation functions. II. J Chem Phys 2010; 133:164105. [DOI: 10.1063/1.3493449] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
41
|
Shaik MS, Liem SY, Popelier PLA. Properties of liquid water from a systematic refinement of a high-rank multipolar electrostatic potential. J Chem Phys 2010; 132:174504. [PMID: 20459171 DOI: 10.1063/1.3409563] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Abstract
We build on previous work [S. Y. Liem and P. L. A. Popelier, J. Chem. Theory Comput. 4, 353 (2008)], where for the first time, a high-rank multipolar electrostatic potential was used in molecular dynamics simulations of liquid water at a wide range of pressures and temperatures, and using a multipolar Ewald summation. Water is represented as a rigid body, with atomic multipole moments defined by quantum chemical topology partitioning its gas phase electron density. The effect of the level of theory on the local structure of liquid water is systematically addressed. Values for Lennard-Jones (LJ) parameters are optimized, for both oxygen and hydrogen atoms, against bulk properties. The best LJ parameters were then used in a set of simulations at 30 different temperatures (1 atm) and another set at 11 different pressures (at 298 K). Inclusion of the hydrogen LJ parameters significantly increases the self-diffusion coefficient. The behavior of bulk properties was studied and the local water structure analyzed by both radial and spatial distribution functions. Comparisons with familiar point-charge potentials, such as TIP3P, TIP4P, TIP5P, and simple point charge, show the benefits of multipole moments.
Collapse
Affiliation(s)
- Majeed S Shaik
- Manchester Interdisciplinary Biocentre (MIB), University of Manchester, 131 Princess Street, Manchester M1 7DN, United Kingdom
| | | | | |
Collapse
|
42
|
Arasa C, Andersson S, Cuppen HM, van Dishoeck EF, Kroes GJ. Molecular dynamics simulations of the ice temperature dependence of water ice photodesorption. J Chem Phys 2010. [DOI: 10.1063/1.3422213] [Citation(s) in RCA: 58] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
43
|
Linearized approximations for condensed phase non-adiabatic dynamics: Multi-layered baths and Brownian dynamics implementation. Chem Phys 2010. [DOI: 10.1016/j.chemphys.2010.02.021] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
44
|
Czakó G, Kaledin AL, Bowman JM. A practical method to avoid zero-point leak in molecular dynamics calculations: Application to the water dimer. J Chem Phys 2010; 132:164103. [DOI: 10.1063/1.3417999] [Citation(s) in RCA: 65] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
|
45
|
Hyeon-Deuk K, Ando K. Quantum effects of hydrogen atoms on the dynamical rearrangement of hydrogen-bond networks in liquid water. J Chem Phys 2010; 132:164507. [DOI: 10.1063/1.3397809] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
46
|
Ivanov SD, Witt A, Shiga M, Marx D. Communications: On artificial frequency shifts in infrared spectra obtained from centroid molecular dynamics: Quantum liquid water. J Chem Phys 2010; 132:031101. [DOI: 10.1063/1.3290958] [Citation(s) in RCA: 60] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
|
47
|
Bonella S, Ciccotti G, Kapral R. Linearization approximations and Liouville quantum–classical dynamics. Chem Phys Lett 2010. [DOI: 10.1016/j.cplett.2009.11.056] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
48
|
Habershon S, Manolopoulos DE. Zero point energy leakage in condensed phase dynamics: An assessment of quantum simulation methods for liquid water. J Chem Phys 2009; 131:244518. [DOI: 10.1063/1.3276109] [Citation(s) in RCA: 108] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
49
|
Liu J, Miller WH, Paesani F, Zhang W, Case DA. Quantum dynamical effects in liquid water: A semiclassical study on the diffusion and the infrared absorption spectrum. J Chem Phys 2009; 131:164509. [PMID: 19894958 DOI: 10.1063/1.3254372] [Citation(s) in RCA: 71] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jian Liu
- Department of Chemistry and K. S. Pitzer Center for Theoretical Chemistry, University of California, Berkeley, California 94720-1460, USA
| | | | | | | | | |
Collapse
|
50
|
Sun X, Yoo S, Xantheas SS, Dang LX. The reorientation mechanism of hydroxide ions in water: A molecular dynamics study. Chem Phys Lett 2009. [DOI: 10.1016/j.cplett.2009.09.004] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|