1
|
Sah MK, Naskar K, Adhikari S, Smits B, Meyer J, Somers MF. On the quantum dynamical treatment of surface vibrational modes for reactive scattering of H2 from Cu(111) at 925 K. J Chem Phys 2024; 161:014306. [PMID: 38953445 DOI: 10.1063/5.0217639] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2024] [Accepted: 06/13/2024] [Indexed: 07/04/2024] Open
Abstract
We construct the effective Hartree potential for H2 on Cu(111) as introduced in our earlier work [Dutta et al., J. Chem. Phys. 154, 104103 (2021), and Dutta et al., J. Chem. Phys. 157, 194112 (2022)] starting from the same gas-metal interaction potential obtained for 0 K. Unlike in that work, we now explicitly account for surface expansion at 925 K and investigate different models to describe the surface vibrational modes: (i) a cluster model yielding harmonic normal modes at 0 K and (ii) slab models resulting in phonons at 0 and 925 K according to the quasi-harmonic approximation-all consistently calculated at the density functional theory level with the same exchange-correlation potential. While performing dynamical calculations for the H2(v = 0, j = 0)-Cu(111) system employing Hartree potential constructed with 925 K phonons and surface temperature, (i) the calculated chemisorption probabilities are the highest compared to the other approaches over the energy domain and (ii) the threshold for the reaction probability is the lowest, in close agreement with the experiment. Although the survival probabilities (v' = 0) depict the expected trend (lower in magnitude), the excitation probabilities (v' = 1) display a higher magnitude since the 925 K phonons and surface temperature are more effective for the excitation process compared to the phonons/normal modes obtained from the other approaches investigated to describe the surface.
Collapse
Affiliation(s)
- Mantu Kumar Sah
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Koushik Naskar
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Satrajit Adhikari
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Bauke Smits
- Leiden Institute of Chemistry, Gorlaeus Building, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Jörg Meyer
- Leiden Institute of Chemistry, Gorlaeus Building, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Mark F Somers
- Leiden Institute of Chemistry, Gorlaeus Building, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
2
|
Liu XY, Wang SR, Fang WH, Cui G. Nuclear Quantum Effects on Nonadiabatic Dynamics of a Green Fluorescent Protein Chromophore Analogue: Ring-Polymer Surface-Hopping Simulation. J Chem Theory Comput 2024; 20:3426-3439. [PMID: 38656202 DOI: 10.1021/acs.jctc.4c00068] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2024]
Abstract
Herein, we have used the "on-the-fly" ring-polymer surface-hopping simulation method with the centroid approximation (RPSH-CA), in combination with the multireference OM2/MRCI electronic structure calculations to study the photoinduced dynamics of a green fluorescent protein (GFP) chromophore analogue in the gas phase, i.e., o-HBI, at 50, 100, and 300 K with 1, 5, 10, and 15 beads (3600 1 ps trajectories). The electronic structure calculations identified five new minimum-energy conical intersection (MECI) structures, which, together with the previous one, play crucial roles in the excited-state decay dynamics of o-HBI. It is also found that the excited-state intramolecular proton transfer (ESIPT) occurs in an ultrafast manner and is completed within 20 fs in all the simulation conditions because there is no barrier associated with this ESIPT process in the S1 state. However, the other excited-state dynamical results are strongly related to the number of beads. At 50 and 100 K, the nuclear quantum effects (NQEs) are very important; therefore, the excited-state dynamical results change significantly with the bead number. For example, the S1 decay time deduced from time-dependent state populations becomes longer as the bead number increases. Nevertheless, an essentially convergent trend is observed when the bead number is close to 10. In contrast, at 300 K, the NQEs become weaker and the above dynamical results converge very quickly even with 1 bead. Most importantly, the NQEs seriously affect the excited-state decay mechanism of o-HBI. At 50 and 100 K, most trajectories decay to the S0 state via perpendicular keto MECIs, whereas, at 300 K, only twisted keto MECIs are responsible for the excited-state decay. The present work not only comprehensively explores the temperature-dependent photoinduced dynamics of o-HBI, but also demonstrates the importance and necessity of NQEs in nonadiabatic dynamics simulations, especially at relatively low temperatures.
Collapse
Affiliation(s)
- Xiang-Yang Liu
- College of Chemistry and Material Science, Sichuan Normal University, Chengdu 610068, China
| | - Sheng-Rui Wang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
3
|
Gerrits N, Jackson B, Bogaerts A. Accurate Reaction Probabilities for Translational Energies on Both Sides of the Barrier of Dissociative Chemisorption on Metal Surfaces. J Phys Chem Lett 2024; 15:2566-2572. [PMID: 38416779 PMCID: PMC10926167 DOI: 10.1021/acs.jpclett.3c03408] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2023] [Revised: 02/16/2024] [Accepted: 02/26/2024] [Indexed: 03/01/2024]
Abstract
Molecular dynamics simulations are essential for a better understanding of dissociative chemisorption on metal surfaces, which is often the rate-controlling step in heterogeneous and plasma catalysis. The workhorse quasi-classical trajectory approach ubiquitous in molecular dynamics is able to accurately predict reactivity only for high translational and low vibrational energies. In contrast, catalytically relevant conditions generally involve low translational and elevated vibrational energies. Existing quantum dynamics approaches are intractable or approximate as a result of the large number of degrees of freedom present in molecule-metal surface reactions. Here, we extend a ring polymer molecular dynamics approach to fully include, for the first time, the degrees of freedom of a moving metal surface. With this approach, experimental sticking probabilities for the dissociative chemisorption of methane on Pt(111) are reproduced for a large range of translational and vibrational energies by including nuclear quantum effects and employing full-dimensional simulations.
Collapse
Affiliation(s)
- Nick Gerrits
- Leiden
Institute of Chemistry, Gorlaeus Laboratories, Leiden University, Post Office
Box 9502, 2300 RA Leiden, Netherlands
- Research
Group PLASMANT, Department of Chemistry, University of Antwerp, Universiteitsplein 1, BE-2610, Wilrijk, Antwerp, Belgium
| | - Bret Jackson
- Department
of Chemistry, University of Massachusetts
Amherst, Amherst, Massachusetts 01003, United States
| | - Annemie Bogaerts
- Research
Group PLASMANT, Department of Chemistry, University of Antwerp, Universiteitsplein 1, BE-2610, Wilrijk, Antwerp, Belgium
| |
Collapse
|
4
|
Videla PE, Batista VS. Matsubara dynamics approximation for generalized multi-time correlation functions. J Chem Phys 2023; 158:2889027. [PMID: 37154285 DOI: 10.1063/5.0146654] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2023] [Accepted: 04/21/2023] [Indexed: 05/10/2023] Open
Abstract
We introduce a semi-classical approximation for calculating generalized multi-time correlation functions based on Matsubara dynamics, a classical dynamics approach that conserves the quantum Boltzmann distribution. This method is exact for the zero time and harmonic limits and reduces to classical dynamics when only one Matsubara mode is considered (i.e., the centroid). Generalized multi-time correlation functions can be expressed as canonical phase-space integrals, involving classically evolved observables coupled through Poisson brackets in a smooth Matsubara space. Numerical tests on a simple potential show that the Matsubara approximation exhibits better agreement with exact results than classical dynamics, providing a bridge between the purely quantum and classical descriptions of multi-time correlation functions. Despite the phase problem that prevents practical applications of Matsubara dynamics, the reported work provides a benchmark theory for the future development of quantum-Boltzmann-preserving semi-classical approximations for studies of chemical dynamics in condensed phase systems.
Collapse
Affiliation(s)
- Pablo E Videla
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520-8107, USA
| | - Victor S Batista
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520-8107, USA
| |
Collapse
|
5
|
Prada A, Pós ES, Althorpe SC. Comparison of Matsubara dynamics with exact quantum dynamics for an oscillator coupled to a dissipative bath. J Chem Phys 2023; 158:114106. [PMID: 36948794 DOI: 10.1063/5.0138250] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/17/2023] Open
Abstract
We report the first numerical calculations in which converged Matsubara dynamics is compared directly with exact quantum dynamics with no artificial damping of the time-correlation functions (TCFs). The system treated is a Morse oscillator coupled to a harmonic bath. We show that, when the system-bath coupling is sufficiently strong, the Matsubara calculations can be converged by explicitly including up to M = 200 Matsubara modes, with the remaining modes included as a harmonic "tail" correction. The resulting Matsubara TCFs are in near-perfect agreement with the exact quantum TCFs, for non-linear as well as linear operators, at a temperature at which the TCFs are dominated by quantum thermal fluctuations. These results provide compelling evidence that incoherent classical dynamics can arise in the condensed phase at temperatures at which the statistics are dominated by quantum (Boltzmann) effects, as a result of smoothing of imaginary-time Feynman paths. The techniques developed here may also lead to efficient methods for benchmarking system-bath dynamics in the overdamped regime.
Collapse
Affiliation(s)
- Adam Prada
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Eszter S Pós
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Stuart C Althorpe
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
6
|
Murakami T, Iida R, Hashimoto Y, Takahashi Y, Takahashi S, Takayanagi T. Ring-Polymer Molecular Dynamics and Kinetics for the H – + C 2H 2 → H 2 + C 2H – Reaction Using the Full-Dimensional Potential Energy Surface. J Phys Chem A 2022; 126:9244-9258. [DOI: 10.1021/acs.jpca.2c05851] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Tatsuhiro Murakami
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama338-8570, Japan
- Department of Materials & Life Sciences, Faculty of Science & Technology, Sophia University, 7-1 Kioicho, Chiyoda-ku, Tokyo102-8554, Japan
| | - Ryusei Iida
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama338-8570, Japan
| | - Yu Hashimoto
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama338-8570, Japan
| | - Yukinobu Takahashi
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama338-8570, Japan
| | - Soma Takahashi
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama338-8570, Japan
| | - Toshiyuki Takayanagi
- Department of Chemistry, Saitama University, Shimo-Okubo 255, Sakura-ku, Saitama City, Saitama338-8570, Japan
| |
Collapse
|
7
|
García-Alfonso E, Barranco M, Bonhommeau DA, Halberstadt N, Pi M, Calvo F. Clustering, collision, and relaxation dynamics in pure and doped helium nanoclusters: Density- vs particle-based approaches. J Chem Phys 2022; 157:014106. [DOI: 10.1063/5.0091942] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The clustering, collision, and relaxation dynamics of pristine and doped helium nanodroplets is theoretically investigated in cases of pickup and clustering of heliophilic argon, collision of heliophobic cesium atoms, and coalescence of two droplets brought into contact by their mutual long-range van der Waals interaction. Three approaches are used and compared with each other. The He time-dependent density functional theory method considers the droplet as a continuous medium and accounts for its superfluid character. The ring-polymer molecular dynamics method uses a path-integral description of nuclear motion and incorporates zero-point delocalization while bosonic exchange effects are ignored. Finally, the zero-point averaged dynamics approach is a mixed quantum–classical method in which quantum delocalization is described by attaching a frozen wavefunction to each He atom, equivalent to classical dynamics with effective interaction potentials. All three methods predict that the growth of argon clusters is significantly hindered by the helium host droplet due to the impeding shell structure around the dopants and kinematic effects freezing the growing cluster in metastable configurations. The effects of superfluidity are qualitatively manifested by different collision dynamics of the heliophilic atom at high velocities, as well as quadrupole oscillations that are not seen with particle-based methods, for droplets experiencing a collision with cesium atoms or merging with each other.
Collapse
Affiliation(s)
- Ernesto García-Alfonso
- Laboratoire Collisions, Agrégats, Réactivité (LCAR), Université de Toulouse, CNRS, 31062 Toulouse, France
| | - Manuel Barranco
- Laboratoire Collisions, Agrégats, Réactivité (LCAR), Université de Toulouse, CNRS, 31062 Toulouse, France
- Department FQA, Facultat de Física, Universitat de Barcelona, Diagonal 645, 08028 Barcelona, Spain
- Institute of Nanoscience and Nanotechnology (IN2UB), Universitat de Barcelona, Barcelona, Spain
| | - David A. Bonhommeau
- Université de Reims Champagne Ardenne, CNRS, GSMA UMR 7331, 51100 Reims, France
| | - Nadine Halberstadt
- Laboratoire Collisions, Agrégats, Réactivité (LCAR), Université de Toulouse, CNRS, 31062 Toulouse, France
| | - Martí Pi
- Department FQA, Facultat de Física, Universitat de Barcelona, Diagonal 645, 08028 Barcelona, Spain
- Institute of Nanoscience and Nanotechnology (IN2UB), Universitat de Barcelona, Barcelona, Spain
| | - Florent Calvo
- Université Grenoble Alpes, CNRS, LIPHY, F38000 Grenoble, France
| |
Collapse
|
8
|
Töpfer K, Upadhyay M, Meuwly M. Quantitative molecular simulations. Phys Chem Chem Phys 2022; 24:12767-12786. [PMID: 35593769 PMCID: PMC9158373 DOI: 10.1039/d2cp01211a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2022] [Accepted: 04/30/2022] [Indexed: 11/21/2022]
Abstract
All-atom simulations can provide molecular-level insights into the dynamics of gas-phase, condensed-phase and surface processes. One important requirement is a sufficiently realistic and detailed description of the underlying intermolecular interactions. The present perspective provides an overview of the present status of quantitative atomistic simulations from colleagues' and our own efforts for gas- and solution-phase processes and for the dynamics on surfaces. Particular attention is paid to direct comparison with experiment. An outlook discusses present challenges and future extensions to bring such dynamics simulations even closer to reality.
Collapse
Affiliation(s)
- Kai Töpfer
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland.
| | - Meenu Upadhyay
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland.
| | - Markus Meuwly
- Department of Chemistry, University of Basel, Klingelbergstrasse 80, CH-4056 Basel, Switzerland.
| |
Collapse
|
9
|
Begušić T, Tao X, Blake GA, Miller TF. Equilibrium-nonequilibrium ring-polymer molecular dynamics for nonlinear spectroscopy. J Chem Phys 2022; 156:131102. [PMID: 35395895 DOI: 10.1063/5.0087156] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
Two-dimensional Raman and hybrid terahertz-Raman spectroscopic techniques provide invaluable insight into molecular structures and dynamics of condensed-phase systems. However, corroborating experimental results with theory is difficult due to the high computational cost of incorporating quantum-mechanical effects in the simulations. Here, we present the equilibrium-nonequilibrium ring-polymer molecular dynamics (RPMD), a practical computational method that can account for nuclear quantum effects on the two-time response function of nonlinear optical spectroscopy. Unlike a recently developed approach based on the double Kubo transformed (DKT) correlation function, our method is exact in the classical limit, where it reduces to the established equilibrium-nonequilibrium classical molecular dynamics method. Using benchmark model calculations, we demonstrate the advantages of the equilibrium-nonequilibrium RPMD over classical and DKT-based approaches. Importantly, its derivation, which is based on the nonequilibrium RPMD, obviates the need for identifying an appropriate Kubo transformed correlation function and paves the way for applying real-time path-integral techniques to multidimensional spectroscopy.
Collapse
Affiliation(s)
- Tomislav Begušić
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Xuecheng Tao
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Geoffrey A Blake
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Thomas F Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
10
|
Marjollet A, Inhester L, Welsch R. Initial state-selected scattering for the reactions H + CH4/CHD3 and F + CHD3 employing ring polymer molecular dynamics. J Chem Phys 2022; 156:044101. [DOI: 10.1063/5.0076216] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- A. Marjollet
- Center for Free-Electron Laser Science CFEL, Deutsches Elektronen-Synchrotron DESY, Notkestrasse 85, 22607 Hamburg, Germany
- Department of Physics, Universität Hamburg, Notkestr. 9-11, 22607 Hamburg, Germany
| | - L. Inhester
- Center for Free-Electron Laser Science CFEL, Deutsches Elektronen-Synchrotron DESY, Notkestrasse 85, 22607 Hamburg, Germany
| | - R. Welsch
- Center for Free-Electron Laser Science CFEL, Deutsches Elektronen-Synchrotron DESY, Notkestrasse 85, 22607 Hamburg, Germany
| |
Collapse
|
11
|
Ananth N. Path Integrals for Nonadiabatic Dynamics: Multistate Ring Polymer Molecular Dynamics. Annu Rev Phys Chem 2022; 73:299-322. [PMID: 35081325 DOI: 10.1146/annurev-physchem-082620-021809] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
This review focuses on a recent class of path-integral-based methods that simulate nonadiabatic dynamics in the condensed phase using only classical molecular dynamics trajectories in an extended phase space. Specifically, a semiclassical mapping protocol is used to derive an exact, continuous, Cartesian variable path-integral representation for the canonical partition function of a system in which multiple electronic states are coupled to nuclear degrees of freedom. Building on this exact statistical foundation, multistate ring polymer molecular dynamics methods are developed for the approximate calculation of real-time thermal correlation functions. The remarkable promise of these multistate ring polymer methods, their successful applications, and their limitations are discussed in detail.Expected final online publication date for the Annual Review of Physical Chemistry, Volume 73 is April 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.
Collapse
Affiliation(s)
- Nandini Ananth
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York, USA;
| |
Collapse
|
12
|
Haggard C, Sadhasivam VG, Trenins G, Althorpe SC. Testing the quasicentroid molecular dynamics method on gas-phase ammonia. J Chem Phys 2021; 155:174120. [PMID: 34742190 DOI: 10.1063/5.0068250] [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/14/2022] Open
Abstract
Quasicentroid molecular dynamics (QCMD) is a path-integral method for approximating nuclear quantum effects in dynamics simulations, which has given promising results for gas- and condensed-phase water. In this work, by simulating the infrared spectrum of gas-phase ammonia, we test the feasibility of extending QCMD beyond water. Overall, QCMD works as well for ammonia as for water, reducing or eliminating blue shifts from the classical spectrum without introducing the artificial red shifts or broadening associated with other imaginary-time path-integral methods. However, QCMD gives only a modest improvement over the classical spectrum for the position of the symmetric bend mode, which is highly anharmonic (since it correlates with the inversion pathway). We expect QCMD to have similar problems with large-amplitude degrees of freedom in other molecules but otherwise to work as well as for water.
Collapse
Affiliation(s)
- Christopher Haggard
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Vijay Ganesh Sadhasivam
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - George Trenins
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | - Stuart C Althorpe
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
13
|
Li C, Liu Q, Zhang L, Li Y, Jiang B. Ring polymer molecular dynamics in gas-surface reactions: tests on initial sampling and potential energy landscape. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1941367] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Chen Li
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, University of Science and Technology of China, Hefei, People’s Republic of China
| | - Qinghua Liu
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, University of Science and Technology of China, Hefei, People’s Republic of China
| | - Liang Zhang
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, University of Science and Technology of China, Hefei, People’s Republic of China
| | - Yongle Li
- Department of Physics, International Center of Quantum and Molecular Structures and Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University, Shanghai, People’s Republic of China
| | - Bin Jiang
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes, University of Science and Technology of China, Hefei, People’s Republic of China
| |
Collapse
|
14
|
Kroes GJ. Computational approaches to dissociative chemisorption on metals: towards chemical accuracy. Phys Chem Chem Phys 2021; 23:8962-9048. [PMID: 33885053 DOI: 10.1039/d1cp00044f] [Citation(s) in RCA: 33] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We review the state-of-the-art in the theory of dissociative chemisorption (DC) of small gas phase molecules on metal surfaces, which is important to modeling heterogeneous catalysis for practical reasons, and for achieving an understanding of the wealth of experimental information that exists for this topic, for fundamental reasons. We first give a quick overview of the experimental state of the field. Turning to the theory, we address the challenge that barrier heights (Eb, which are not observables) for DC on metals cannot yet be calculated with chemical accuracy, although embedded correlated wave function theory and diffusion Monte-Carlo are moving in this direction. For benchmarking, at present chemically accurate Eb can only be derived from dynamics calculations based on a semi-empirically derived density functional (DF), by computing a sticking curve and demonstrating that it is shifted from the curve measured in a supersonic beam experiment by no more than 1 kcal mol-1. The approach capable of delivering this accuracy is called the specific reaction parameter (SRP) approach to density functional theory (DFT). SRP-DFT relies on DFT and on dynamics calculations, which are most efficiently performed if a potential energy surface (PES) is available. We therefore present a brief review of the DFs that now exist, also considering their performance on databases for Eb for gas phase reactions and DC on metals, and for adsorption to metals. We also consider expressions for SRP-DFs and briefly discuss other electronic structure methods that have addressed the interaction of molecules with metal surfaces. An overview is presented of dynamical models, which make a distinction as to whether or not, and which dissipative channels are modeled, the dissipative channels being surface phonons and electronically non-adiabatic channels such as electron-hole pair excitation. We also discuss the dynamical methods that have been used, such as the quasi-classical trajectory method and quantum dynamical methods like the time-dependent wave packet method and the reaction path Hamiltonian method. Limits on the accuracy of these methods are discussed for DC of diatomic and polyatomic molecules on metal surfaces, paying particular attention to reduced dimensionality approximations that still have to be invoked in wave packet calculations on polyatomic molecules like CH4. We also address the accuracy of fitting methods, such as recent machine learning methods (like neural network methods) and the corrugation reducing procedure. In discussing the calculation of observables we emphasize the importance of modeling the properties of the supersonic beams in simulating the sticking probability curves measured in the associated experiments. We show that chemically accurate barrier heights have now been extracted for DC in 11 molecule-metal surface systems, some of which form the most accurate core of the only existing database of Eb for DC reactions on metal surfaces (SBH10). The SRP-DFs (or candidate SRP-DFs) that have been derived show transferability in many cases, i.e., they have been shown also to yield chemically accurate Eb for chemically related systems. This can in principle be exploited in simulating rates of catalyzed reactions on nano-particles containing facets and edges, as SRP-DFs may be transferable among systems in which a molecule dissociates on low index and stepped surfaces of the same metal. In many instances SRP-DFs have allowed important conclusions regarding the mechanisms underlying observed experimental trends. An important recent observation is that SRP-DFT based on semi-local exchange DFs has so far only been successful for systems for which the difference of the metal work function and the molecule's electron affinity exceeds 7 eV. A main challenge to SRP-DFT is to extend its applicability to the other systems, which involve a range of important DC reactions of e.g. O2, H2O, NH3, CO2, and CH3OH. Recent calculations employing a PES based on a screened hybrid exchange functional suggest that the road to success may be based on using exchange functionals of this category.
Collapse
Affiliation(s)
- Geert-Jan Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, The Netherlands.
| |
Collapse
|
15
|
Chowdhury SN, Huo P. Non-adiabatic Matsubara dynamics and non-adiabatic ring-polymer molecular dynamics. J Chem Phys 2021; 154:124124. [PMID: 33810665 DOI: 10.1063/5.0042136] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Abstract
We present the non-adiabatic Matsubara dynamics, a general framework for computing the time-correlation function (TCF) of electronically non-adiabatic systems. This new formalism is derived based on the generalized Kubo-transformed TCF using the Wigner representation for both the nuclear degrees of freedom and the electronic mapping variables. By dropping the non-Matsubara nuclear normal modes in the quantum Liouvillian and explicitly integrating these modes out from the expression of the TCF, we derived the non-adiabatic Matsubara dynamics approach. Further making the approximation to drop the imaginary part of the Matsubara Liouvillian and enforce the nuclear momentum integral to be real, we arrived at the non-adiabatic ring-polymer molecular dynamics (NRPMD) approach. We have further justified the capability of NRPMD for simulating the non-equilibrium TCF. This work provides the rigorous theoretical foundation for several recently proposed state-dependent RPMD approaches and offers a general framework for developing new non-adiabatic quantum dynamics methods in the future.
Collapse
Affiliation(s)
- Sutirtha N Chowdhury
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| | - Pengfei Huo
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| |
Collapse
|
16
|
Jiang H, Tao X, Kammler M, Ding F, Wodtke AM, Kandratsenka A, Miller TF, Bünermann O. Small Nuclear Quantum Effects in Scattering of H and D from Graphene. J Phys Chem Lett 2021; 12:1991-1996. [PMID: 33596383 DOI: 10.1021/acs.jpclett.0c02933] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
We study nuclear quantum effects in H/D sticking to graphene, comparing scattering experiments at near-zero coverage with classical, quantized, and transition-state calculations. The experiment shows H/D sticking probabilities that are indistinguishable from one another and markedly smaller than those expected from a consideration of zero-point energy shifts of the chemisorption transition state. Inclusion of dynamical effects and vibrational anharmonicity via ring-polymer molecular dynamics (RPMD) yields results that are in good agreement with the experimental results. RPMD also reveals that nuclear quantum effects, while modest, arise primarily from carbon and not from H/D motion, confirming the importance of a C atom rehybridization mechanism associated with H/D sticking on graphene.
Collapse
Affiliation(s)
- Hongyan Jiang
- Department of Dynamics at Surfaces, Max-Planck-Institute for Biophysical Chemistry, Am Faßerg 11, 37077 Göttingen, Germany
| | - Xuecheng Tao
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Marvin Kammler
- Department of Dynamics at Surfaces, Max-Planck-Institute for Biophysical Chemistry, Am Faßerg 11, 37077 Göttingen, Germany
| | - Feizhi Ding
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Alec M Wodtke
- Department of Dynamics at Surfaces, Max-Planck-Institute for Biophysical Chemistry, Am Faßerg 11, 37077 Göttingen, Germany
- Institut für Physikalische Chemie, Georg-Georg-August-Universität Göttingen, Tammanstraße 6, 37077 Göttingen, Germany
- International Center for Advanced Studies of Energy Conversion, Georg-August University of Göttingen, Tammanstraße 6, 37077 Göttingen, Germany
| | - Alexander Kandratsenka
- Department of Dynamics at Surfaces, Max-Planck-Institute for Biophysical Chemistry, Am Faßerg 11, 37077 Göttingen, Germany
| | - Thomas F Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Oliver Bünermann
- Department of Dynamics at Surfaces, Max-Planck-Institute for Biophysical Chemistry, Am Faßerg 11, 37077 Göttingen, Germany
- Institut für Physikalische Chemie, Georg-Georg-August-Universität Göttingen, Tammanstraße 6, 37077 Göttingen, Germany
- International Center for Advanced Studies of Energy Conversion, Georg-August University of Göttingen, Tammanstraße 6, 37077 Göttingen, Germany
| |
Collapse
|
17
|
Brown SE, Shakib FA. Recent progress in approximate quantum dynamics methods for the study of proton-coupled electron transfer reactions. Phys Chem Chem Phys 2021; 23:2535-2556. [PMID: 33367437 DOI: 10.1039/d0cp05166g] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
Proton-coupled electron transfer (PCET) reactions are ubiquitous natural processes at the heart of energy conversion reactions in photosynthesis and respiration, DNA repair, and diverse enzymatic reactions. Theoretical formulation and computational method developments have eyed modeling of thermal and photoinduced PCET for the last three decades. The accumulation of these studies, collected in dozens of reviews, accounts, and perspectives, has firmly established the influence of quantum effects, including non-adiabatic electronic transitions, vibrational relaxation, zero-point energy, and proton tunneling, on the rate and mechanism of PCET reactions. Here, we focus on some recently-developed methods, spanning the last eight years, that can quantitatively capture these effects in the PCET context and provide efficient means for their qualitative description in complex systems. The theoretical background of each method and their accuracy with respect to exact results are discussed and the results of relevant PCET simulations based on each method are presented.
Collapse
Affiliation(s)
- Sandra E Brown
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, CA 90095, USA
| | - Farnaz A Shakib
- Department of Chemistry and Environmental Science, New Jersey Institute of Technology, Newark, NJ 07102, USA.
| |
Collapse
|
18
|
Remsing RC, Bates JE. Effective mass path integral simulations of quasiparticles in condensed phases. J Chem Phys 2020; 153:121104. [PMID: 33003737 DOI: 10.1063/5.0020555] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The quantum many-body problem in condensed phases is often simplified using a quasiparticle description, such as effective mass theory for electron motion in a periodic solid. These approaches are often the basis for understanding many fundamental condensed phase processes, including the molecular mechanisms underlying solar energy harvesting and photocatalysis. Despite the importance of these effective particles, there is still a need for computational methods that can explore their behavior on chemically relevant length and time scales. This is especially true when the interactions between the particles and their environment are important. We introduce an approach for studying quasiparticles in condensed phases by combining effective mass theory with the path integral treatment of quantum particles. This framework incorporates the generally anisotropic electronic band structure of materials into path integral simulation schemes to enable modeling of quasiparticles in quantum confinement, for example. We demonstrate the utility of effective mass path integral simulations by modeling an exciton in solid potassium chloride and electron trapping by a sulfur vacancy in monolayer molybdenum disulfide.
Collapse
Affiliation(s)
- Richard C Remsing
- Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, USA
| | - Jefferson E Bates
- A. R. Smith Department of Chemistry and Fermentation Sciences, Appalachian State University, Boone, North Carolina 28608, USA
| |
Collapse
|
19
|
Jung KA, Videla PE, Batista VS. Ring-polymer, centroid, and mean-field approximations to multi-time Matsubara dynamics. J Chem Phys 2020; 153:124112. [PMID: 33003707 DOI: 10.1063/5.0021843] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Based on a recently developed generalization of Matsubara dynamics to the multi-time realm, we present a formal derivation of multi-time generalizations of ring-polymer molecular dynamics, thermostatted ring-polymer molecular dynamics (TRPMD), centroid molecular dynamics (CMD), and mean-field Matsubara dynamics. Additionally, we analyze the short-time accuracy of each methodology. We find that for multi-time correlation functions of linear operators, (T)RPMD is accurate up to order t3, while CMD is only correct up to t, indicating a degradation in the accuracy of these methodologies with respect to the single-time counterparts. The present work provides a firm justification for the use of path-integral-based approximations for the calculation of multi-time correlation functions.
Collapse
Affiliation(s)
- Kenneth A Jung
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, USA
| | - Pablo E Videla
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, USA
| | - Victor S Batista
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, USA
| |
Collapse
|
20
|
Marjollet A, Welsch R. Nuclear quantum effects in state-selective scattering from ring polymer molecular dynamics. J Chem Phys 2020; 152:194113. [DOI: 10.1063/5.0004179] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Adrien Marjollet
- Center for Free-Electron Laser Science, DESY, Notkestrasse 85, 22607 Hamburg, Germany
- Department of Physics, Universität Hamburg, Jungiusstrasse 9, 20355 Hamburg, Germany
| | - Ralph Welsch
- Center for Free-Electron Laser Science, DESY, Notkestrasse 85, 22607 Hamburg, Germany
- The Hamburg Centre for Ultrafast Imaging, Luruper Chaussee 149, 22761 Hamburg, Germany
| |
Collapse
|
21
|
Tao X, Shushkov P, Miller TF. Microcanonical rates from ring-polymer molecular dynamics: Direct-shooting, stationary-phase, and maximum-entropy approaches. J Chem Phys 2020; 152:124117. [DOI: 10.1063/1.5144307] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Affiliation(s)
- Xuecheng Tao
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Philip Shushkov
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Thomas F. Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
22
|
Amati G, Schilling T. Structural localization in the classical and quantum Fermi-Pasta-Ulam model. CHAOS (WOODBURY, N.Y.) 2020; 30:033116. [PMID: 32237780 DOI: 10.1063/1.5130740] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/09/2019] [Accepted: 02/19/2020] [Indexed: 06/11/2023]
Abstract
We study the statistics and short-time dynamics of the classical and the quantum Fermi-Pasta-Ulam chain in the thermal equilibrium. We analyze the distributions of single-particle configurations by integrating out the rest of the system. At low temperatures, we observe a systematic increase in the mobility of the chain when transitioning from classical to quantum mechanics due to zero-point energy effects. We analyze the consequences of quantum dispersion on the dynamics at short times of configurational correlation functions.
Collapse
Affiliation(s)
- Graziano Amati
- Physikalisches Institut, Albert-Ludwigs-Universität, 79104 Freiburg, Germany
| | - Tanja Schilling
- Physikalisches Institut, Albert-Ludwigs-Universität, 79104 Freiburg, Germany
| |
Collapse
|
23
|
Nelson TR, White AJ, Bjorgaard JA, Sifain AE, Zhang Y, Nebgen B, Fernandez-Alberti S, Mozyrsky D, Roitberg AE, Tretiak S. Non-adiabatic Excited-State Molecular Dynamics: Theory and Applications for Modeling Photophysics in Extended Molecular Materials. Chem Rev 2020; 120:2215-2287. [PMID: 32040312 DOI: 10.1021/acs.chemrev.9b00447] [Citation(s) in RCA: 198] [Impact Index Per Article: 49.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Optically active molecular materials, such as organic conjugated polymers and biological systems, are characterized by strong coupling between electronic and vibrational degrees of freedom. Typically, simulations must go beyond the Born-Oppenheimer approximation to account for non-adiabatic coupling between excited states. Indeed, non-adiabatic dynamics is commonly associated with exciton dynamics and photophysics involving charge and energy transfer, as well as exciton dissociation and charge recombination. Understanding the photoinduced dynamics in such materials is vital to providing an accurate description of exciton formation, evolution, and decay. This interdisciplinary field has matured significantly over the past decades. Formulation of new theoretical frameworks, development of more efficient and accurate computational algorithms, and evolution of high-performance computer hardware has extended these simulations to very large molecular systems with hundreds of atoms, including numerous studies of organic semiconductors and biomolecules. In this Review, we will describe recent theoretical advances including treatment of electronic decoherence in surface-hopping methods, the role of solvent effects, trivial unavoided crossings, analysis of data based on transition densities, and efficient computational implementations of these numerical methods. We also emphasize newly developed semiclassical approaches, based on the Gaussian approximation, which retain phase and width information to account for significant decoherence and interference effects while maintaining the high efficiency of surface-hopping approaches. The above developments have been employed to successfully describe photophysics in a variety of molecular materials.
Collapse
Affiliation(s)
- Tammie R Nelson
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Alexander J White
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Josiah A Bjorgaard
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Andrew E Sifain
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States.,U.S. Army Research Laboratory , Aberdeen Proving Ground , Maryland 21005 , United States
| | - Yu Zhang
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Benjamin Nebgen
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | | | - Dmitry Mozyrsky
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| | - Adrian E Roitberg
- Department of Chemistry , University of Florida , Gainesville , Florida 32611 , United States
| | - Sergei Tretiak
- Theoretical Division , Los Alamos National Laboratory , Los Alamos , New Mexico 87545 , United States
| |
Collapse
|
24
|
Liu Q, Zhang L, Li Y, Jiang B. Ring Polymer Molecular Dynamics in Gas-Surface Reactions: Inclusion of Quantum Effects Made Simple. J Phys Chem Lett 2019; 10:7475-7481. [PMID: 31738557 DOI: 10.1021/acs.jpclett.9b02570] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Accurately modeling gas-surface collision dynamics presents a great challenge for theory, especially in the low-energy (or temperature) regime where quantum effects are important. Here, a path integral-based nonequilibrium ring polymer molecular dynamics (NE-RPMD) approach is adapted to calculate dissociative initial sticking probabilities (S0) of H2 on Cu(111) and D2O on Ni(111), revealing the distinct quantum nature in the two benchmark surface reactions. NE-RPMD successfully captures quantum tunneling in H2 dissociation at very low energies, where the quasi-classical trajectory (QCT) method suddenly fails. Additionally, QCT substantially overestimates S0 of D2O because of severe zero point energy (ZPE) leakage, even at collision energies greater than the ZPE-corrected barrier. Instead, NE-RPMD predicts S0 values of D2O in much improved agreement with reference results obtained by the quantum wavepacket method with reasonable corrections of the thermal contribution. Our results suggest NE-RPMD as a promising approach to model quantum effects in gas-surface reactions.
Collapse
Affiliation(s)
- Qinghua Liu
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes , University of Science and Technology of China , Hefei , Anhui 230026 , China
| | - Liang Zhang
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes , University of Science and Technology of China , Hefei , Anhui 230026 , China
| | - Yongle Li
- Department of Physics, International Center of Quantum and Molecular Structures and Shanghai Key Laboratory of High Temperature Superconductors , Shanghai University , Shanghai 200444 , China
| | - Bin Jiang
- Hefei National Laboratory for Physical Science at the Microscale, Department of Chemical Physics, Key Laboratory of Surface and Interface Chemistry and Energy Catalysis of Anhui Higher Education Institutes , University of Science and Technology of China , Hefei , Anhui 230026 , China
| |
Collapse
|
25
|
Guarini E, Neumann M, Bafile U, Bellissima S, Colognesi D. Dynamical Origin of the Total and Zero-Point Kinetic Energy in a Quantum Fluid. PHYSICAL REVIEW LETTERS 2019; 123:135301. [PMID: 31697543 DOI: 10.1103/physrevlett.123.135301] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/25/2019] [Indexed: 06/10/2023]
Abstract
By applying an exponential mode analysis to ring polymer molecular dynamics simulations of dense fluid parahydrogen, we find that the dynamical processes establishing the time behavior of the Kubo velocity autocorrelation function display the same nature as those already observed in high-density classical fluids. This result permits us to demonstrate that the exponential mode decomposition is a unique tool to identify which dynamical processes lead to one of the most notable properties of quantum fluids: the large value of the mean kinetic energy per particle and the importance of the zero-temperature quantum effects in determining it.
Collapse
Affiliation(s)
- Eleonora Guarini
- Dipartimento di Fisica e Astronomia, Università degli Studi di Firenze, via G. Sansone 1, I-50019 Sesto Fiorentino, Italy
| | - Martin Neumann
- Fakultät für Physik der Universität Wien, Strudlhofgasse 4, A-1090 Wien, Austria
| | - Ubaldo Bafile
- Consiglio Nazionale delle Ricerche, Istituto di Fisica Applicata "Nello Carrara," via Madonna del Piano 10, I-50019 Sesto Fiorentino, Italy
| | - Stefano Bellissima
- Consiglio Nazionale delle Ricerche, Istituto di Fisica Applicata "Nello Carrara," via Madonna del Piano 10, I-50019 Sesto Fiorentino, Italy
| | - Daniele Colognesi
- Consiglio Nazionale delle Ricerche, Istituto di Fisica Applicata "Nello Carrara," via Madonna del Piano 10, I-50019 Sesto Fiorentino, Italy
| |
Collapse
|
26
|
Chowdhury SN, Huo P. State dependent ring polymer molecular dynamics for investigating excited nonadiabatic dynamics. J Chem Phys 2019; 150:244102. [DOI: 10.1063/1.5096276] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Sutirtha N. Chowdhury
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| | - Pengfei Huo
- Department of Chemistry, University of Rochester, 120 Trustee Road, Rochester, New York 14627, USA
| |
Collapse
|
27
|
Jiang H, Kammler M, Ding F, Dorenkamp Y, Manby FR, Wodtke AM, Miller TF, Kandratsenka A, Bünermann O. Imaging covalent bond formation by H atom scattering from graphene. Science 2019; 364:379-382. [DOI: 10.1126/science.aaw6378] [Citation(s) in RCA: 60] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2019] [Accepted: 03/06/2019] [Indexed: 11/02/2022]
Abstract
Viewing the atomic-scale motion and energy dissipation pathways involved in forming a covalent bond is a longstanding challenge for chemistry. We performed scattering experiments of H atoms from graphene and observed a bimodal translational energy loss distribution. Using accurate first-principles dynamics simulations, we show that the quasi-elastic channel involves scattering through the physisorption well where collision sites are near the centers of the six-membered C-rings. The second channel results from transient C–H bond formation, where H atoms lose 1 to 2 electron volts of energy within a 10-femtosecond interaction time. This remarkably rapid form of intramolecular vibrational relaxation results from the C atom’s rehybridization during bond formation and is responsible for an unexpectedly high sticking probability of H on graphene.
Collapse
|
28
|
Kaur R, Welsch R. Probing photodissociation dynamics using ring polymer molecular dynamics. J Chem Phys 2019; 150:114105. [PMID: 30901996 DOI: 10.1063/1.5086218] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The performance of the ring polymer molecular dynamics (RPMD) approach to simulate typical photodissociation processes is assessed. The correct description of photodissociation requires the calculation of correlation functions or expectation values associated with non-equilibrium initial conditions, which was shown to be possible with RPMD very recently [J. Chem. Phys. 145, 204118 (2016)]. This approach is combined with treatment of the nonadiabatic dynamics employing the ring polymer surface hopping approach (RPSH), which is based on Tully's fewest switches surface hopping (FSSH) approach. The performance is tested using one-dimensional photodissociation models. It is found that RPSH with non-equilibrium initial conditions can well reproduce the time-dependent dissociation probability, and adiabatic and diabatic populations for cases where the crossing point is below and above the Franck-Condon point, respectively, while standard FSSH fails to reproduce the exact quantum dynamics in the latter case. Thus, it is shown that RPSH is an efficient and accurate alternative to standard FSSH, which is one of the most widely employed approaches to study photochemistry.
Collapse
Affiliation(s)
- Rajwant Kaur
- Center for Free-Electron Laser Science, DESY, Notkestrasse 85, 22607 Hamburg, Germany
| | - Ralph Welsch
- Center for Free-Electron Laser Science, DESY, Notkestrasse 85, 22607 Hamburg, Germany
| |
Collapse
|
29
|
Bellissima S, Neumann M, Bafile U, Colognesi D, Barocchi F, Guarini E. Density and time scaling effects on the velocity autocorrelation function of quantum and classical dense fluid para-hydrogen. J Chem Phys 2019; 150:074502. [DOI: 10.1063/1.5085202] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Affiliation(s)
- S. Bellissima
- Consiglio Nazionale delle Ricerche, Istituto di Fisica Applicata “Nello Carrara,” via Madonna del Piano 10, I-50019 Sesto Fiorentino, Italy
| | - M. Neumann
- Fakultät für Physik der Universität Wien, Strudlhofgasse 4, A-1090 Wien, Austria
| | - U. Bafile
- Consiglio Nazionale delle Ricerche, Istituto di Fisica Applicata “Nello Carrara,” via Madonna del Piano 10, I-50019 Sesto Fiorentino, Italy
| | - D. Colognesi
- Consiglio Nazionale delle Ricerche, Istituto di Fisica Applicata “Nello Carrara,” via Madonna del Piano 10, I-50019 Sesto Fiorentino, Italy
| | - F. Barocchi
- Dipartimento di Fisica e Astronomia, Università degli Studi di Firenze, via G. Sansone 1, I-50019 Sesto Fiorentino, Italy
| | - E. Guarini
- Dipartimento di Fisica e Astronomia, Università degli Studi di Firenze, via G. Sansone 1, I-50019 Sesto Fiorentino, Italy
| |
Collapse
|
30
|
Welsch R. Kinetic isotope effects in the water forming reaction H2/D2 + OH from rigorous close-coupling quantum dynamics simulations. Phys Chem Chem Phys 2019; 21:17054-17062. [DOI: 10.1039/c9cp02323b] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Rigorous quantum dynamics simulations of thermal rate constants and kinetic isotope effects for the water-forming H2/D2 + OH reaction are presented, which show increased tunneling below 300 K and can serve as benchmarks for approximate methods.
Collapse
Affiliation(s)
- Ralph Welsch
- Center for Free-Electron Laser Science
- DESY
- 22607 Hamburg
- Germany
- The Hamburg Centre for Ultrafast Imaging
| |
Collapse
|
31
|
Ekanayake NT, Garashchuk S. Evaluation of the quantum time-correlation functions employing the Hamilton–Jacobi dynamics framework. Theor Chem Acc 2018. [DOI: 10.1007/s00214-018-2398-6] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
32
|
Welsch R. Low‐Temperature Thermal Rate Constants for the Water Formation Reaction H
2
+OH from Rigorous Quantum Dynamics Calculations. Angew Chem Int Ed Engl 2018; 57:13150-13153. [DOI: 10.1002/anie.201807666] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2018] [Indexed: 11/06/2022]
Affiliation(s)
- Ralph Welsch
- Center for Free-Electron Laser Science, DESY Notkestraße 85 22607 Hamburg Germany
| |
Collapse
|
33
|
Welsch R. Low‐Temperature Thermal Rate Constants for the Water Formation Reaction H
2
+OH from Rigorous Quantum Dynamics Calculations. Angew Chem Int Ed Engl 2018. [DOI: 10.1002/ange.201807666] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Ralph Welsch
- Center for Free-Electron Laser Science, DESY Notkestraße 85 22607 Hamburg Germany
| |
Collapse
|
34
|
Jung KA, Videla PE, Batista VS. Inclusion of nuclear quantum effects for simulations of nonlinear spectroscopy. J Chem Phys 2018; 148:244105. [PMID: 29960352 DOI: 10.1063/1.5036768] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The computation and interpretation of nonlinear vibrational spectroscopy is of vital importance for understanding a wide range of dynamical processes in molecular systems. Here, we introduce an approach to evaluate multi-time response functions in terms of multi-time double symmetrized Kubo transformed thermal correlation functions. Furthermore, we introduce a multi-time extension of ring polymer molecular dynamics to evaluate these Kubo transforms. Benchmark calculations show that the approximations are useful for short times even for nonlinear operators, providing a consistent improvement over classical simulations of multi-time correlation functions. The introduced methodology thus provides a practical way of including nuclear quantum effects in multi-time response functions of non-linear optical spectroscopy.
Collapse
Affiliation(s)
- Kenneth A Jung
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, USA
| | - Pablo E Videla
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, USA
| | - Victor S Batista
- Department of Chemistry, Yale University, P.O. Box 208107, New Haven, Connecticut 06520-8107, USA
| |
Collapse
|
35
|
Mishra PK, Bettaque V, Vendrell O, Santra R, Welsch R. Prospects of Using High-Intensity THz Pulses To Induce Ultrafast Temperature-Jumps in Liquid Water. J Phys Chem A 2018; 122:5211-5222. [DOI: 10.1021/acs.jpca.8b00828] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Pankaj Kr. Mishra
- Center for Free-Electron Laser Science, DESY, Notkestraße 85, D-22607 Hamburg, Germany
- The Hamburg Centre for Ultrafast Imaging, University of Hamburg, Luruper Chaussee 149, D-22761 Hamburg, Germany
| | - Vincent Bettaque
- Department of Physics, University of Hamburg, Jungiusstraße 9, D-20355 Hamburg, Germany
| | - Oriol Vendrell
- Department of Physics and Astronomy, Aarhus University, DK-8000 Aarhus C, Denmark
| | - Robin Santra
- Center for Free-Electron Laser Science, DESY, Notkestraße 85, D-22607 Hamburg, Germany
- The Hamburg Centre for Ultrafast Imaging, University of Hamburg, Luruper Chaussee 149, D-22761 Hamburg, Germany
- Department of Physics, University of Hamburg, Jungiusstraße 9, D-20355 Hamburg, Germany
| | - Ralph Welsch
- Center for Free-Electron Laser Science, DESY, Notkestraße 85, D-22607 Hamburg, Germany
| |
Collapse
|
36
|
Yuan Y, Li J, Li XZ, Wang F. The strengths and limitations of effective centroid force models explored by studying isotopic effects in liquid water. J Chem Phys 2018; 148:184102. [PMID: 29764147 PMCID: PMC5940466 DOI: 10.1063/1.5027433] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2018] [Accepted: 04/20/2018] [Indexed: 12/15/2022] Open
Abstract
The development of effective centroid potentials (ECPs) is explored with both the constrained-centroid and quasi-adiabatic force matching using liquid water as a test system. A trajectory integrated with the ECP is free of statistical noises that would be introduced when the centroid potential is approximated on the fly with a finite number of beads. With the reduced cost of ECP, challenging experimental properties can be studied in the spirit of centroid molecular dynamics. The experimental number density of H2O is 0.38% higher than that of D2O. With the ECP, the H2O number density is predicted to be 0.42% higher, when the dispersion term is not refit. After correction of finite size effects, the diffusion constant of H2O is found to be 21% higher than that of D2O, which is in good agreement with the 29.9% higher diffusivity for H2O observed experimentally. Although the ECP is also able to capture the redshifts of both the OH and OD stretching modes in liquid water, there are a number of properties that a classical simulation with the ECP will not be able to recover. For example, the heat capacities of H2O and D2O are predicted to be almost identical and higher than the experimental values. Such a failure is simply a result of not properly treating quantized vibrational energy levels when the trajectory is propagated with classical mechanics. Several limitations of the ECP based approach without bead population reconstruction are discussed.
Collapse
Affiliation(s)
- Ying Yuan
- State Key Laboratory for Mesoscopic Physics and School of Physics, Peking University, Beijing 100871, People’s Republic of China
| | - Jicun Li
- 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
|
37
|
Affiliation(s)
- Timothy J. H. Hele
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, NY, USA
| |
Collapse
|
38
|
Seki Y, Takayanagi T, Shiga M. Photoexcited Ag ejection from a low-temperature He cluster: a simulation study by nonadiabatic Ehrenfest ring-polymer molecular dynamics. Phys Chem Chem Phys 2017; 19:13798-13806. [DOI: 10.1039/c7cp00888k] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Nonadiabatic ring-polymer molecular dynamics simulations were performed to understand the photoexcitation dynamics of a low-temperature Ag·He500 cluster.
Collapse
Affiliation(s)
- Yusuke Seki
- Department of Chemistry
- Saitama University
- Saitama City
- Japan
| | | | - Motoyuki Shiga
- Center for Computational Science and E-Systems
- Japan Atomic Energy Agency
- Kashiwa
- Japan
| |
Collapse
|
39
|
Fu B, Shan X, Zhang DH, Clary DC. Recent advances in quantum scattering calculations on polyatomic bimolecular reactions. Chem Soc Rev 2017; 46:7625-7649. [DOI: 10.1039/c7cs00526a] [Citation(s) in RCA: 63] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
This review surveys quantum scattering calculations on chemical reactions of polyatomic molecules in the gas phase published in the last ten years.
Collapse
Affiliation(s)
- Bina Fu
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry
- Dalian Institute of Chemical Physics
- Chinese Academy of Sciences
- Dalian 116023
- China
| | - Xiao Shan
- Physical and Theoretical Chemistry Laboratory
- Department of Chemistry
- University of Oxford
- Oxford
- UK
| | - Dong H. Zhang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry
- Dalian Institute of Chemical Physics
- Chinese Academy of Sciences
- Dalian 116023
- China
| | - David C. Clary
- Physical and Theoretical Chemistry Laboratory
- Department of Chemistry
- University of Oxford
- Oxford
- UK
| |
Collapse
|