1
|
Menéndez M, Veselinova A, Zanchet A, Jambrina PG, Aoiz FJ. Rate coefficients for the O + H 2 and O + D 2 reactions: how well ring polymer molecular dynamics accounts for tunelling. Phys Chem Chem Phys 2024; 26:20947-20961. [PMID: 39046374 DOI: 10.1039/d4cp01711k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/25/2024]
Abstract
We present here extensive calculations of the O(3P) + H2 and O(3P) + D2 reaction dynamics spanning the temperature range from 200 K to 2500 K. The calculations have been carried out using fully converged time-independent quantum mechanics (TI QM), quasiclassical trajectories (QCT) and ring polymer molecular dynamics (RPMD) on the two lowest lying adiabatic potential energy surfaces (PESs), 13A' and 13A'', calculated by Zanchet et al. [J. Chem. Phys., 2019, 151, 094307]. TI QM rate coefficients were determined using the cumulative reaction probability formalism on each PES including all of the total angular momenta and the Coriolis coupling and can be considered to be essentially exact within the Born-Oppenheimer approximation. The agreement between the rate coefficients calculated by using QM and RPMD is excellent for the reaction with D2 in almost the whole temperature range. For the reaction with H2, although the agreement is very good above 500 K, the deviations are significant at lower temperatures. In contrast, the QCT calculations largely underestimate the rate coefficients for the two isotopic variants due to their inability to account for tunelling. The differences found in the disagreements between RPMD and QM rate coefficients for the reactions for both the isotopologues are indicative of the ability of the RPMD method to accurately describe systems where tunelling plays a relevant role. Considering that both reactions are dominated by tunelling below 500 K, the present results show that RPMD is a very powerful tool for determining rate coefficients. The present QM rate coefficients calculated on adiabatic PESs slightly underestimate the best global fits of the experimental measurements, which we attribute to the intersystem crossing with the singlet 11A' PES.
Collapse
Affiliation(s)
- Marta Menéndez
- Departamento de Química Física, Unidad Asociada CSIC, Universidad Complutense de Madrid, 28040 Madrid, Spain.
| | - Anzhela Veselinova
- Departamento de Química Física, Universidad de Salamanca, 37008 Salamanca, Spain
| | - Alexandre Zanchet
- Instituto de Física Fundamental, CSIC, C/Serrano 121-123, 28006 Madrid, Spain
| | - Pablo G Jambrina
- Departamento de Química Física, Universidad de Salamanca, 37008 Salamanca, Spain
| | - F Javier Aoiz
- Departamento de Química Física, Unidad Asociada CSIC, Universidad Complutense de Madrid, 28040 Madrid, Spain.
| |
Collapse
|
2
|
Espinosa-Garcia J, Bhowmick S. Kinetic study of the CN + C 2H 6 hydrogen abstraction reaction based on an analytical potential energy surface. Phys Chem Chem Phys 2024; 26:8344-8355. [PMID: 38391269 DOI: 10.1039/d3cp05930h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/24/2024]
Abstract
The temperature dependence of the thermal rate constants and kinetic isotope effects (KIE) of the CN + C2H6 gas-phase hydrogen abstraction reaction was theoretically determined within the 25-1000 K temperature range, i.e., from very low- to high-temperature regimes. Based on a recently developed full-dimensional analytical potential energy surface fitted to highly accurate explicitly correlated ab initio calculations, three different kinetic theories were used: canonical variational transition state theory (CVT), quasiclassical trajectory theory (QCT), and ring polymer molecular dynamics (RPMD) method for the computation of rate constants. We found that the thermal rate constants obtained with the three theories show a V-shaped temperature dependence, with a pronounced minimum near 200 K, qualitatively reproducing the experimental measurements. Among the three methods used in this work, the QCT and RPMD methods have the best agreement with the experiment at low and high temperatures, respectively, while the CVT model shows the largest discrepancies. The significant increase in the rate constant at very low temperatures in this very exothermic and practically barrierless reaction could be attributed to the large value of the impact parameter, possibly ruling out the role of the tunneling effect and the intermediate complexes in the entrance channel. The theoretical H/D KIE depicted a "normal" behaviour, i.e., values greater than unity, emulating the experimental measurements and improving previous theoretical results. Finally, the discrepancies between theory and experiments were analysed as a function of several factors, such as limitations of the kinetics theories and the potential energy surface, as well as the uncertainties in the experimental measurements.
Collapse
Affiliation(s)
- Joaquin Espinosa-Garcia
- Departamento de Química Física and Instituto de Computación Científica Avanzada, Universidad de Extremadura, 06071 Badajoz, Spain.
| | - Somnath Bhowmick
- Climate and Atmosphere Research Centre, The Cyprus Institute, Nicosia 2121, Cyprus.
| |
Collapse
|
3
|
Zhang S, Chen Q, Zhang L, Li J, Hu X, Xie D. Dynamics studies for the multi-well and multi-channel reaction of OH with C 2H 2 on a full-dimensional global potential energy surface. Phys Chem Chem Phys 2024; 26:7351-7362. [PMID: 38375620 DOI: 10.1039/d3cp05811e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/21/2024]
Abstract
The C2H2 + OH reaction is an important acetylene oxidation pathway in the combustion process, as well as a typical multi-well and multi-channel reaction. Here, we report an accurate full-dimensional machine learning-based potential energy surface (PES) for the C2H2 + OH reaction at the UCCSD(T)-F12b/cc-pVTZ-F12 level, based on about 475 000 ab initio points. Extensive quasi-classical trajectory (QCT) calculations were performed on the newly developed PES to obtain detailed dynamic data and analyze reaction mechanisms. Below 1000 K, the C2H2 + OH reaction produces H + OCCH2 and CO + CH3. With increasing temperature, the product channels H2O + C2H and H + HCCOH are accessible and the former dominates above 1900 K. It is found that the formation of H2O + C2H is dominated by a direct reaction process, while other channels belong to the indirect mechanism involving long-lived intermediates along the reaction pathways. At low temperatures, the C2H2 + OH reaction behaves like an unimolecular reaction due to the unique PES topographic features, of which the dynamic features are similar to the decomposition of energy-rich complexes formed by C2H2 + OH collision. The classification of trajectories that undergo different reaction pathways to generate each product and their product energy distributions were also reported in this work. This dynamic information may provide a deep understanding of the C2H2 + OH reaction.
Collapse
Affiliation(s)
- Shuwen Zhang
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Qixin Chen
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
| | - Lidong Zhang
- National Synchrotron Radiation Laboratory, University of Science and Technology of China, Hefei 230029, China; State Key Laboratory of Fire Science, University of Science and Technology of China, Hefei 230026, China.
| | - Jun Li
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Chongqing 401331, China.
| | - Xixi Hu
- Kuang Yaming Honors School, Nanjing University, Nanjing 210023, China.
- Hefei National Laboratory, Hefei 230088, China
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210023, China
- Hefei National Laboratory, Hefei 230088, China
| |
Collapse
|
4
|
Steffen J. Caracal: A Versatile Ring Polymer Molecular Dynamics Simulation Package. J Chem Theory Comput 2023; 19:5334-5355. [PMID: 37555628 DOI: 10.1021/acs.jctc.3c00568] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/10/2023]
Abstract
A new open-source program package named Caracal covering simulations of molecular systems with ring polymer molecular dynamics (RPMD) is presented. It combines a powerful RPMD implementation including chemical reaction rate calculations and biased periodic and nonperiodic samplings with a collection of easy to set up potential energy surface (PES) methodologies, thus delivering an all-inclusive approach. Most implemented PESs are based on the QMDFF and EVB-QMDFF methods. Where the quantum mechanically derived force field (QMDFF) can be set up for an arbitrary molecular system in a black-box fashion, the empirical valence bond (EVB)-QMDFF connects two QMDFFs and is able to represent the PES of a chemical reaction. With our previously published flavors of this composite method, PESs for almost arbitrary gas phase thermal ground state reactions can be set up. Given an optimized reaction path, the mechanism of the reaction can be classified and RPMD rate constants can be obtained via umbrella sampling and recrossing calculations on an EVB-QMDFF PES. Further, QMDFFs can be polymerized for the description of liquid systems. In this paper, the internal structure as well as the handling philosophy of Caracal are outlined. Further, examples of the different possible kinds of calculations are given.
Collapse
Affiliation(s)
- Julien Steffen
- Chair of Theoretical Chemistry, Friedrich-Alexander University Erlangen-Nürnberg, Egerlandstraße 3, 91058 Erlangen, Bavaria, Germany
| |
Collapse
|
5
|
Zhang L, Zuo J, Suleimanov YV, Guo H. Ring Polymer Molecular Dynamics Approach to Quantum Dissociative Chemisorption Rates. J Phys Chem Lett 2023; 14:7118-7125. [PMID: 37531595 DOI: 10.1021/acs.jpclett.3c01848] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/04/2023]
Abstract
A ring polymer molecular dynamics (RPMD) method is proposed for the calculation of the dissociative chemisorption rate coefficient on surfaces. The RPMD rate theory is capable of handling quantum effects such as the zero-point energy and tunneling in dissociative chemisorption, while it relies on classical trajectories for the simulation. Applications to H2 dissociative chemisorption are demonstrated. For the highly activated process on Ag(111), strong deviations from Arrhenius behavior are found at low temperatures and attributed to tunneling. On Pt(111), where the dissociation has a barrierless pathway, the RPMD rate coefficient is found to agree with the experimentally derived thermal sticking coefficient within a factor of 2 over a large temperature range. Significant quantum effects are also found.
Collapse
Affiliation(s)
- Liang Zhang
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Junxiang Zuo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Yury V Suleimanov
- American Association for the Advancement of Science, 1200 New York Ave NW, Washington, D.C. 20005, United States
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
6
|
Li C, Li Y, Jiang B. First-principles surface reaction rates by ring polymer molecular dynamics and neural network potential: role of anharmonicity and lattice motion. Chem Sci 2023; 14:5087-5098. [PMID: 37206404 PMCID: PMC10189860 DOI: 10.1039/d2sc06559b] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2022] [Accepted: 04/05/2023] [Indexed: 08/04/2023] Open
Abstract
Elementary gas-surface processes are essential steps in heterogeneous catalysis. A predictive understanding of catalytic mechanisms remains challenging due largely to difficulties in accurately characterizing the kinetics of such steps. Experimentally, thermal rates for elementary surface reactions can now be measured using a novel velocity imaging technique, providing a stringent testing ground for ab initio rate theories. Here, we propose to combine ring polymer molecular dynamics (RPMD) rate theory with state-of-the-art first-principles-determined neural network potential to calculate surface reaction rates. Taking NO desorption from Pd(111) as an example, we show that the harmonic approximation and the neglect of lattice motion in the commonly-used transition state theory overestimates and underestimates the entropy change during the desorption process, respectively, leading to opposite errors in rate coefficient predictions and artificial error cancellations. Including anharmonicity and lattice motion, our results reveal a generally neglected surface entropy change due to significant local structural change during desorption and obtain the right answer for the right reasons. Although quantum effects are found to be less important in this system, the proposed approach establishes a more reliable theoretical benchmark for accurately predicting the kinetics of elementary gas-surface processes.
Collapse
Affiliation(s)
- Chen Li
- Key Laboratory of Precision and Intelligent Chemistry, 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, Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University Shanghai 200444 China
| | - Bin Jiang
- Key Laboratory of Precision and Intelligent Chemistry, 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
|
7
|
Poulsen JA, Nyman G. Computing reaction rates using a Wigner function based on a new effective frequency theory. Mol Phys 2023. [DOI: 10.1080/00268976.2023.2174361] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023]
Affiliation(s)
- Jens Aage Poulsen
- Department of Chemistry and Molecular Biology, University of Gothenburg, Gothenburg, Sweden
| | - Gunnar Nyman
- Department of Chemistry and Molecular Biology, University of Gothenburg, Gothenburg, Sweden
| |
Collapse
|
8
|
Angiolari F, Huppert S, Spezia R. Quantum versus classical unimolecular fragmentation rate constants and activation energies at finite temperature from direct dynamics simulations. Phys Chem Chem Phys 2022; 24:29357-29370. [PMID: 36448557 DOI: 10.1039/d2cp03809a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
In the present work, we investigate how nuclear quantum effects modify the temperature dependent rate constants and, consequently, the activation energies in unimolecular reactions. In the reactions under study, nuclear quantum effects mainly stem from the presence of a large zero point energy. Thus, we investigate the behavior of methods compatible with direct dynamics simulations, the quantum thermal bath (QTB) and ring polymer molecular dynamics (RPMD). To this end, we first compare them with quantum reaction theory for a model Morse potential before extending this comparison to molecular models. Our results show that, in particular in the temperature range comparable with or lower than the zero point energy of the system, the RPMD method is able to correctly capture nuclear quantum effects on rate constants and activation energies. On the other hand, although the QTB provides a good description of equilibrium properties including zero-point energy effects, it largely overestimates the rate constants. The origin of the different behaviours is in the different distance distributions provided by the two methods and in particular how they differently describe the tails of such distributions. The comparison with transition state theory shows that RPMD can be used to study fragmentation of complex systems for which it may be difficult to determine the multiple reaction pathways and associated transition states.
Collapse
Affiliation(s)
- Federica Angiolari
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, 4 Place Jussieu, 75005 Paris, France.
| | - Simon Huppert
- Sorbonne Université, Institut de Nanosciences de Paris, UMR 7588 CNRS, 4 Place Jussieu, 75005 Paris, France
| | - Riccardo Spezia
- Sorbonne Université, Laboratoire de Chimie Théorique, UMR 7616 CNRS, 4 Place Jussieu, 75005 Paris, France.
| |
Collapse
|
9
|
Gui X, Fan W, Sun J, Li Y. New Stable and Fast Ring-Polymer Molecular Dynamics for Calculating Bimolecular Rate Coefficients with an Example of OH + CH 4. J Chem Theory Comput 2022; 18:5203-5212. [PMID: 35983956 DOI: 10.1021/acs.jctc.2c00522] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The accurate and efficient calculation of the rate coefficients of chemical reactions is a key issue in the research of chemical dynamics. In this work, by applying the dimension-free ultrastable Cayley propagator, the thermal rate coefficients of a prototypic high dimensional chemical reaction OH + CH4 → H2O + CH3 in the temperature range of 200 to 1500 K are investigated with ring polymer molecular dynamics (RPMD) on a highly accurate full-dimensional potential energy surface. Kinetic isotope effects (KIEs) for three isotopologues of the title reaction are also studied. The results demonstrate excellent agreement with experimental data, even in the deep tunneling region. Especially, the Cayley propagator shows a high calculation efficiency with little loss of accuracy. The present results confirmed the applicability of the RPMD method, particularly the speed-up using a Cayley propagator, in theoretical calculations of bimolecular reaction rates.
Collapse
Affiliation(s)
- Xiongfei Gui
- Department of Physics, International Center of Quantum and Molecular Structures, and Shanghai Frontiers Science Center of Quantum and Superconducting Matter States, Shanghai 200444, China
| | - Wenbin Fan
- Department of Chemistry, Fudan University, Shanghai 200433, China
| | - Jiace Sun
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Yongle Li
- Department of Physics, International Center of Quantum and Molecular Structures, and Shanghai Frontiers Science Center of Quantum and Superconducting Matter States, Shanghai 200444, China
| |
Collapse
|
10
|
Méndez E, Videla PE, Laria D. Equilibrium and Dynamical Characteristics of Hydrogen Bond Bifurcations in Water-Water and Water-Ammonia Dimers: A Path Integral Molecular Dynamics Study. J Phys Chem A 2022; 126:4721-4733. [PMID: 35834556 DOI: 10.1021/acs.jpca.2c02525] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present path integral molecular dynamics results that describe the effects of nuclear quantum fluctuations on equilibrium and dynamical characteristics pertaining to bifurcation pathways in hydrogen bonded dimers combining water and ammonia, at cryogenic temperatures of the order of 20 K. Along these isomerizations, the hydrogen atoms in the molecules acting as hydrogen-bond donors interchange their original dangling/connective characters. Our results reveal that the resulting quantum transition paths comprise three stages: the initial and final ones involve overall rotations during which the two protons retain their classical-like characteristics. Effects from quantum fluctuation are clearly manifested in the changes operated at the intermediate passages over transition states, as the spatial extents of the protons stretch over typical lengths comparable to the distances between connective and dangling basins of attractions. Consequently, the classical over-the-hill path is replaced by a tunneling controlled mechanism which, within the path integral perspective, can be cast in terms of concerted inter-basin migrations of polymer beads from dangling-to-connective and from connective-to-dangling, at practically no energy costs. We also estimated the characteristic timescales describing such interconversions within the approximate ring polymer rate theory. Effects derived from full and partial deuteration are also discussed.
Collapse
Affiliation(s)
- Emilio Méndez
- Departamento de Química Inorgánica, Analítica y Química-Física and INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
| | - Pablo E Videla
- Department of Chemistry and Energy Sciences Institute, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, United States
| | - Daniel Laria
- Departamento de Física de la Materia Condensada, Comisión Nacional de Energía Atómica, Avenida Libertador 8250, 1429 Buenos Aires, Argentina
| |
Collapse
|
11
|
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
|
12
|
Meng Q, Chen J, Ma J, Zhang X, Chen J. Adiabatic models for the quantum dynamics of surface scattering with lattice effects. Phys Chem Chem Phys 2022; 24:16415-16436. [PMID: 35766107 DOI: 10.1039/d2cp01560a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this contribution, we review models for the lattice effects in quantum dynamics calculations on surface scattering, which is important to modeling heterogeneous catalysis for achieving an interpretation of experimental measurements. Unlike dynamics models for reactions in the gas phase, those for heterogeneous reactions have to include the effects of the surface. For manageable computational costs in calculations, the effects of static surface (SS) are firstly modeled as this is simply and easily implemented. Then, the SS model has to be improved to include the effects of the flexible surface, that is the lattice effects. To do this, various surface models have been designed where the coordinates of the surface atoms are introduced in the Hamiltonian operator, especially those of the top surface atom. Based on this model Hamiltonian operator, extensive multi-dimension quantum dynamics calculations can be performed to recover the lattice effects. Here, we first review an overview of the techniques in constructing the Hamiltonian operator, which is a sum of the kinetic energy operator (KEO) and potential energy surface (PES). Since the PES containing the coordinates of the surface atoms in a cell is still expensive, the SS model is often accepted. We consider a mathematical model, called the coupled harmonic oscillator (CHO) model, to introduce the concepts of adiabatic and diabatic representations for separating the molecule and surface. Under the adiabatic model, we further introduce the expansion model where the potential function is Taylor expanded around the optimized geometry of the surface. By an expansion model truncated at the first and second order, various coupling surface models between the molecule and surface are derived. Moreover, by further and deeply understanding the adiabatic representation, an effective Hamiltonian operator is obtained by optimizing the total wave function in factorized form. By this factorized form of wave function and effective Hamiltonian operator, the geometry phase of the surface wave function is theoretically found. This theoretical prediction may be measured by carefully designing experiments. Finally, discussions on the adiabatic representation, the PES construction, and possibility of the classical-dynamics solutions are given. Based on these discussions, a simple outlook on the dynamics of photocatalytics is finally given.
Collapse
Affiliation(s)
- Qingyong Meng
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China.
| | - Junbo Chen
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China. .,Xi'an Modern Chemistry Research Institute, China North Industries Group Corp., Ltd., East Zhangba Road 168, 710065 Xi'an, China
| | - Jianxing Ma
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China.
| | - Xingyu Zhang
- Department of Chemistry, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi'an, China.
| | - Jun Chen
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Yangqiao Road West 155, 350002 Fuzhou, China.,Fujian Science & Technology Innovation Laboratory for Optoelectronic Information of China, Optoelectronic Industry Base at High-tech Zone, 350108 Fuzhou, China
| |
Collapse
|
13
|
Hoppe H, Manthe U. First-Principles Theory for the Reaction of Chlorine with Methane. J Phys Chem Lett 2022; 13:2563-2566. [PMID: 35285640 DOI: 10.1021/acs.jpclett.2c00407] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
A full-dimensional quantum dynamics simulation of the Cl + CH4 → HCl + CH3 reaction based on first-principles theory is reported. Accurate thermal rate constants are calculated, and perfect agreement with experiment is obtained. Despite the heavy atoms present in both reactants, the passage of the reaction barrier is found to occur within only a few tens of femtoseconds. This surprisingly short time scale results from correlated motion of the transferring hydrogen atom and the hydrogen atoms in the methyl fragment which facilitates irreversible barrier passage without relevant participation of heavy atoms. Resonance effects resulting from the heavy-light-heavy characteristics of the reaction system, which were observed in reactive scattering studies, do not affect the thermal rate constant.
Collapse
Affiliation(s)
- Hannes Hoppe
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstr. 25, D-33615 Bielefeld, Germany
| | - Uwe Manthe
- Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Universitätsstr. 25, D-33615 Bielefeld, Germany
| |
Collapse
|
14
|
Xu X, Chen J, Lu X, Fang W, Liu S, Zhang DH. Strong non-Arrhenius behavior at low temperatures in the OH + HCl → H 2O + Cl reaction due to resonance induced quantum tunneling. Chem Sci 2022; 13:7955-7961. [PMID: 35865883 PMCID: PMC9258319 DOI: 10.1039/d2sc01958b] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2022] [Accepted: 06/11/2022] [Indexed: 11/21/2022] Open
Abstract
The OH + HCl reaction possesses many Feshbach resonances trapped in the hydrogen bond well in the entrance channel, which substantially enhance the reaction rates at low temperatures.
Collapse
Affiliation(s)
- Xin Xu
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, Liaoning, China
| | - Jun Chen
- State Key Laboratory of Structure Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou 350002, Fujian, China
| | - Xiaoxiao Lu
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, Liaoning, China
| | - Wei Fang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, Liaoning, China
| | - Shu Liu
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, Liaoning, China
| | - Dong H. Zhang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, Liaoning, China
| |
Collapse
|
15
|
Cao J, Wu Y, Bian W. Ring polymer molecular dynamics of the C(1D)+H2 reaction on the most recent potential energy surfaces. CHINESE J CHEM PHYS 2021. [DOI: 10.1063/1674-0068/cjcp2110197] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
Affiliation(s)
- Jianwei Cao
- Beijing National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
| | - Yanan Wu
- Beijing National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
- School of Chemical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Wensheng Bian
- Beijing National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
- School of Chemical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| |
Collapse
|
16
|
Bhowmick S, Hernández MI, Campos-Martínez J, Suleimanov YV. Isotopic separation of helium through graphyne membranes: a ring polymer molecular dynamics study. Phys Chem Chem Phys 2021; 23:18547-18557. [PMID: 34612392 DOI: 10.1039/d1cp02121d] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Microscopic-level understanding of the separation mechanism for two-dimensional (2D) membranes is an active area of research due to potential implications of this class of membranes for various technological processes. Helium (He) purification from the natural resources is of particular interest due to the shortfall in its production. In this work, we applied the ring polymer molecular dynamics (RPMD) method to graphdiyne (Gr2) and graphtriyne (Gr3) 2D membranes having variable pore sizes for the separation of He isotopes, and compare for the first time with rigorous quantum calculations. We found that the transmission rate through Gr3 is many orders of magnitude greater than Gr2. The selectivity of either isotope at low temperatures is a consequence of a delicate balance between the zero-point energy effect and tunneling of 4He and 3He. In particular, a remarkable tunneling effect is reported on the Gr2 membrane at 10 K, leading to a much larger permeation of the lighter species as compared to the heavier isotope. RPMD provides an efficient approach for studying the separation of He isotopes, taking into account quantum effects of light nuclei motions at low temperatures, which classical methods fail to capture.
Collapse
Affiliation(s)
- Somnath Bhowmick
- Computation-based Science and Technology Research Center, The Cyprus Institute, 20 Konstantinou Kavafi Street, Nicosia 2121, Cyprus.
| | | | | | | |
Collapse
|
17
|
Fang JH, Fan WB, Yang H, Song JN, Li YL. Rate coefficients and kinetic isotope effects of Cl+XCl→XCl+Cl (X=H, D, Mu) reactions from ring polymer molecular dynamics. CHINESE J CHEM PHYS 2021. [DOI: 10.1063/1674-0068/cjcp2007117] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Affiliation(s)
- Jun-hua Fang
- Department of Chemistry, Shanghai University, Shanghai 200444, China
| | - Wen-bin Fan
- Department of Physics, International Center for Quantum and Molecular Structures, and Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University, Shanghai 200444, China
| | - Hui Yang
- Department of Physics, International Center for Quantum and Molecular Structures, and Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University, Shanghai 200444, China
| | - Jia-ning Song
- Department of Physics, International Center for Quantum and Molecular Structures, and Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University, Shanghai 200444, China
| | - Yong-le Li
- Department of Physics, International Center for Quantum and Molecular Structures, and Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University, Shanghai 200444, China
| |
Collapse
|
18
|
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
|
19
|
Yang H, Fan WB, Fang JH, Song J, Li Y. Rate coefficients of roaming reaction H+MgH using ring polymer molecular dynamics. CHINESE J CHEM PHYS 2021. [DOI: 10.1063/1674-0068/cjcp2007121] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Affiliation(s)
- Hui Yang
- Department of Physics, International Center for Quantum and Molecular Structures, and Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University, Shanghai 200444, China
| | - Wen-bin Fan
- Department of Physics, International Center for Quantum and Molecular Structures, and Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University, Shanghai 200444, China
| | - Jun-hua Fang
- Department of Chemistry, Shanghai University, Shanghai 200444, China
| | - Jianing Song
- Department of Physics, International Center for Quantum and Molecular Structures, and Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University, Shanghai 200444, China
| | - Yongle Li
- Department of Physics, International Center for Quantum and Molecular Structures, and Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University, Shanghai 200444, China
- Devision of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena 91125, USA
| |
Collapse
|
20
|
Cao J, Wu Y, Ma H, Shen Z, Bian W. Dynamics and kinetics of the Si( 1D) + H 2/D 2 reactions on a new global ab initio potential energy surface. Phys Chem Chem Phys 2021; 23:6141-6153. [PMID: 33684184 DOI: 10.1039/d0cp05540a] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Recent studies on the exothermic complex-forming reactions have improved our understanding of complex-forming reactions greatly, however, so far a similar level of study on endothermic ones has been rather limited. In this work, the endothermic complex-forming reaction Si(1D) + H2 → SiH + H and its deuterated isotopic variant are investigated by quantum dynamics (QD) and ring polymer molecular dynamics (RPMD) calculations on a new global ab initio potential energy surface (PES) for the ground electronic state, which is constructed based on 8996 symmetry unique points computed at the icMRCI+Q/aug-cc-pVQZ level. The PES reproduces our ab initio data very well in the dynamically important regions, on which the ro-vibrational energy levels of SiH2 are calculated and general good agreement with experiment is obtained. The integral cross sections and product angular and state distributions are computed in a wide range of collision energies, and interesting dynamics behaviors are revealed. The rate coefficients are also investigated, which display an exponential rise from 2.09 × 10-20 to 6.00 × 10-12 cm3 s-1 for the Si(1D) + H2 reaction as the temperature increases from 300 to 1500 K, in contrast to the nearly temperature-independent behavior of exothermic complex-forming reactions. In addition, the applicability of the RPMD approach is demonstrated, and the kinetic isotope effect is investigated, the ratio of which decreases from 7.89 (300 K) to 1.70 (1500 K). The effects of tunneling and initial rotational excitation are also discussed.
Collapse
Affiliation(s)
- Jianwei Cao
- Beijing National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China.
| | - Yanan Wu
- Beijing National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China. and School of Chemical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Haitao Ma
- Beijing National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China.
| | - Zhitao Shen
- Henan Key Laboratory of Photovoltaic Materials, Henan University, Kaifeng 475004, China
| | - Wensheng Bian
- Beijing National Laboratory for Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China. and School of Chemical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| |
Collapse
|
21
|
Chen Q, Hu X, Guo H, Xie D. Theoretical H + O 3 rate coefficients from ring polymer molecular dynamics on an accurate global potential energy surface: assessing experimental uncertainties. Phys Chem Chem Phys 2021; 23:3300-3310. [PMID: 33506830 DOI: 10.1039/d0cp05771a] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Thermal rate coefficients and kinetic isotope effects have been calculated for an important atmospheric reaction H/D + O3 → OH/OD + O2 based on an accurate permutation invariant polynomial-neural network potential energy surface, using ring polymer molecular dynamics (RPMD), quasi-classical trajectory (QCT) and variational transition-state theory (VTST) with multidimensional tunneling. The RPMD approach yielded results that are generally in better agreement with experimental rate coefficients than the VTST and QCT ones, especially at low temperatures, attributable to its capacity to capture quantum effects such as tunneling and zero-point energy. The theoretical results support one group of existing experiments over the other. In addition, rate coefficients for the D + O3 → OD + O2 reaction are also reported using the same methods, which will allow a stringent assessment of future experimental measurements, thus helping to reduce the uncertainty in the recommended rate coefficients of this reaction.
Collapse
Affiliation(s)
- Qixin Chen
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210093, China
| | - Xixi Hu
- Kuang Yaming Honors School, Nanjing University, Nanjing 210023, China.
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing 210093, China
| |
Collapse
|
22
|
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
|
23
|
Liu Y, Song H, Li J. Kinetic study of the OH + HO 2 → H 2O + O 2 reaction using ring polymer molecular dynamics and quantum dynamics. Phys Chem Chem Phys 2020; 22:23657-23664. [PMID: 33112305 DOI: 10.1039/d0cp04120c] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The reaction OH + HO2 → H2O + O2 is a prototype of radical-radical reactions. It plays an important role in interstellar/atmospheric chemistry and combustion, and considerable attention has thus been dedicated to its kinetics. In our previous work, we reported an accurate full-dimensional potential energy surface for the title reaction on the ground triplet electronic state. The quasi-classical trajectory (QCT) approach was employed to investigate its kinetics. Although the QCT rate coefficients were in good agreement with some experimental and theoretical results, QCT cannot account for the quantum mechanical effects, such as zero-point vibrational energy, recrossing, and tunneling, which may significantly affect the rate coefficients, particularly at low temperatures. In this work, the reduced-dimensional quantum dynamics and ring polymer molecular dynamics calculations were carried out to examine these effects and their impact on rate coefficients over the temperature range of 300-1300 K.
Collapse
Affiliation(s)
- Yang Liu
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Chongqing 401331, China.
| | - Hongwei Song
- State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, Wuhan Institute of Physics and Mathematics, Innovation Academy for Precision Measurement Science and Technology, Chinese Academy of Sciences, Wuhan 430071, China.
| | - Jun Li
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Theoretical and Computational Chemistry, Chongqing University, Chongqing 401331, China.
| |
Collapse
|
24
|
Espinosa-Garcia J, Garcia-Chamorro M, Corchado JC, Bhowmick S, Suleimanov YV. VTST and RPMD kinetics study of the nine-body X + C 2H 6 (X ≡ H, Cl, F) reactions based on analytical potential energy surfaces. Phys Chem Chem Phys 2020; 22:13790-13801. [PMID: 32538410 DOI: 10.1039/d0cp02238a] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Thermal rate constants of nine-atom hydrogen abstraction reactions, X + C2H6 → HX + C2H5 (X ≡ H, Cl, F) with qualitatively different reaction paths, have been investigated using two kinetics approaches - variational transition state theory with multidimensional tunnelling (VTST/MT) and ring polymer molecular dynamics (RPMD) - and full dimensional analytical potential energy surfaces. For the H + C2H6 reaction, which proceeds through a noticeable barrier height of 11.62 kcal mol-1, kinetics approaches showed excellent agreement between them (with differences less than 30%) and with the experiment (with differences less than 60%) in the wide temperature range of 200-2000 K. For X = Cl and F, however, the situation is very different. The barrier height is either low or very low, 2.44 and 0.23 kcal mol-1, respectively, and the presence of van der Waals complexes in the entrance channel leads to a very flat topography and, consequently, imposes theoretical challenges. For the Cl(2P) reaction, VTST/MT underestimates the experimental rate constants (with differences less than 86%), and RPMD demonstrates better agreement (with differences less than 47%), although the temperature dependence is opposite to the experiment at low temperatures. Finally, for the F(2P) reaction, available experimental information shows discrepancies, both in the absolute values of the rate constants and also in the temperature dependence. Unfortunately, kinetics theories did not resolve this discrepancy. Different possible causes of these theory/experiment discrepancies were analyzed.
Collapse
Affiliation(s)
- Joaquin Espinosa-Garcia
- Departamento de Quimica Fisica and Instituto de Computacion Cientifica Avanzada, Universidad de Extremadura, 06071 Badajoz, Spain.
| | - Moises Garcia-Chamorro
- Departamento de Quimica Fisica and Instituto de Computacion Cientifica Avanzada, Universidad de Extremadura, 06071 Badajoz, Spain.
| | - Jose C Corchado
- Departamento de Quimica Fisica and Instituto de Computacion Cientifica Avanzada, Universidad de Extremadura, 06071 Badajoz, Spain.
| | - Somnath Bhowmick
- Computation-based Science and Technology Research Center, The Cyprus Institute, 20 Konstantinou Kavafi Street, Nicosia 2121, Cyprus.
| | - Yury V Suleimanov
- Computation-based Science and Technology Research Center, The Cyprus Institute, 20 Konstantinou Kavafi Street, Nicosia 2121, Cyprus.
| |
Collapse
|
25
|
Song Q, Zhang Q, Meng Q. Revisiting the Gaussian process regression for fitting high-dimensional potential energy surface and its application to the OH + HO2→O2+ H2O reaction. J Chem Phys 2020; 152:134309. [PMID: 32268765 DOI: 10.1063/1.5143544] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Qingfei Song
- School of Chemistry and Chemical Engineering, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi’an, China
- Key Laboratory of Special Functional and Smart Polymer Materials of Ministry of Industry and Information Technology, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi’an, China
| | - Qiuyu Zhang
- School of Chemistry and Chemical Engineering, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi’an, China
- Key Laboratory of Special Functional and Smart Polymer Materials of Ministry of Industry and Information Technology, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi’an, China
| | - Qingyong Meng
- School of Chemistry and Chemical Engineering, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi’an, China
- Key Laboratory of Special Functional and Smart Polymer Materials of Ministry of Industry and Information Technology, Northwestern Polytechnical University, West Youyi Road 127, 710072 Xi’an, China
| |
Collapse
|
26
|
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
|
27
|
Liu Y, Li J. An accurate potential energy surface and ring polymer molecular dynamics study of the Cl + CH4→ HCl + CH3reaction. Phys Chem Chem Phys 2020; 22:344-353. [DOI: 10.1039/c9cp05693a] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Thermal rate coefficients for the Cl + CH4/CD4reactions were studied on a new full-dimensional accurate potential energy surface with the spin–orbit corrections considered in the entrance channel.
Collapse
Affiliation(s)
- Yang Liu
- School of Chemistry and Chemical Engineering
- Chongqing University
- Chongqing 401331
- China
| | - Jun Li
- School of Chemistry and Chemical Engineering
- Chongqing University
- Chongqing 401331
- China
| |
Collapse
|
28
|
Lawrence JE, Manolopoulos DE. Path integral methods for reaction rates in complex systems. Faraday Discuss 2020; 221:9-29. [DOI: 10.1039/c9fd00084d] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We shall use this introduction to the Faraday Discussion on quantum effects in complex systems to review the recent progress that has been made in using imaginary time path integral methods to calculate chemical reaction rates.
Collapse
Affiliation(s)
- Joseph E. Lawrence
- Department of Chemistry
- University of Oxford
- Physical and Theoretical Chemistry Laboratory
- Oxford OX1 3QZ
- UK
| | - David E. Manolopoulos
- Department of Chemistry
- University of Oxford
- Physical and Theoretical Chemistry Laboratory
- Oxford OX1 3QZ
- UK
| |
Collapse
|
29
|
Novikov IS, Shapeev AV, Suleimanov YV. Ring polymer molecular dynamics and active learning of moment tensor potential for gas-phase barrierless reactions: Application to S + H2. J Chem Phys 2019; 151:224105. [DOI: 10.1063/1.5127561] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ivan S. Novikov
- Skolkovo Institute of Science and Technology, Skolkovo Innovation Center, Nobel St. 3, Moscow 143026, Russia
| | - Alexander V. Shapeev
- Skolkovo Institute of Science and Technology, Skolkovo Innovation Center, Nobel St. 3, Moscow 143026, Russia
| | - Yury V. Suleimanov
- Computation-Based Science and Technology Research Center, Cyprus Institute, 20 Kavafi Street, Nicosia 2121, Cyprus
| |
Collapse
|
30
|
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
|
31
|
González-Lezana T, Bossion D, Scribano Y, Bhowmick S, Suleimanov YV. Dynamics of H + HeH +( v = 0, j = 0) → H 2+ + He: Insight on the Possible Complex-Forming Behavior of the Reaction. J Phys Chem A 2019; 123:10480-10489. [PMID: 31725286 DOI: 10.1021/acs.jpca.9b06122] [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/29/2022]
Abstract
The H + HeH+→ He + H2+ reaction has been studied by means of a combination of theoretical approaches: a statistical quantum method (SQM), ring polymer molecular dynamics (RPMD), and the quasiclassical trajectory (QCT) method. Cross sections and rate constants have been calculated in an attempt to investigate the dynamics of the process. The comparison with previous calculations and experimental results reveals that despite the fact that statistical predictions seem to reproduce some of the overall observed features, the analysis at a more detailed state-to-state level shows noticeable deviations from a complex-forming dynamics. We find some differences in cross sections and rate constants obtained in the QCT calculation with a Gaussian binning procedure with respect to previous works in which the standard histogram binning was employed.
Collapse
Affiliation(s)
| | - Duncan Bossion
- Laboratoire Univers et Particules de Montpellier, UMR-CNRS 5299 , Université de Montpellier , 34095 Montpellier Cedex , France
| | - Yohann Scribano
- Laboratoire Univers et Particules de Montpellier, UMR-CNRS 5299 , Université de Montpellier , 34095 Montpellier Cedex , France
| | - Somnath Bhowmick
- Computation-based Science and Technology Research Center , The Cyprus Institute , 20 Konstantinou Kavafi Street , Nicosia 2121 , Cyprus
| | - Yury V Suleimanov
- Computation-based Science and Technology Research Center , The Cyprus Institute , 20 Konstantinou Kavafi Street , Nicosia 2121 , Cyprus
| |
Collapse
|
32
|
Lawrence JE, Fletcher T, Lindoy LP, Manolopoulos DE. On the calculation of quantum mechanical electron transfer rates. J Chem Phys 2019; 151:114119. [DOI: 10.1063/1.5116800] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Joseph E. Lawrence
- Department of Chemistry, University of Oxford, Physical and Theoretical Chemistry Laboratory, South Parks Road, Oxford, OX1 3QZ, United Kingdom
| | - Theo Fletcher
- Department of Chemistry, University of Oxford, Physical and Theoretical Chemistry Laboratory, South Parks Road, Oxford, OX1 3QZ, United Kingdom
| | - Lachlan P. Lindoy
- Department of Chemistry, University of Oxford, Physical and Theoretical Chemistry Laboratory, South Parks Road, Oxford, OX1 3QZ, United Kingdom
| | - David E. Manolopoulos
- Department of Chemistry, University of Oxford, Physical and Theoretical Chemistry Laboratory, South Parks Road, Oxford, OX1 3QZ, United Kingdom
| |
Collapse
|
33
|
Menéndez M, Jambrina PG, Zanchet A, Verdasco E, Suleimanov YV, Aoiz FJ. New Stress Test for Ring Polymer Molecular Dynamics: Rate Coefficients of the O( 3P) + HCl Reaction and Comparison with Quantum Mechanical and Quasiclassical Trajectory Results. J Phys Chem A 2019; 123:7920-7931. [PMID: 31461272 DOI: 10.1021/acs.jpca.9b06695] [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/29/2022]
Abstract
In the past decade, ring polymer molecular dynamics (RPMD) has emerged as a very efficient method to determine thermal rate coefficients for a great variety of chemical reactions. This work presents the application of this methodology to study the O(3P) + HCl reaction, which constitutes a stringent test for any dynamical calculation due to rich resonant structure and other dynamical features. The rate coefficients, calculated on the 3A' and 3A″ potential energy surfaces (PESs) by Ramachandran and Peterson [ J. Chem. Phys. 2003 , 119 , 9590 ], using RPMD and quasiclassical trajectories (QCT) are compared with the existing experimental and the quantum mechanical (QM) results by Xie et al. [ J. Chem. Phys. 2005 122 , 014301 ]. The agreement is very good at T > 600 K, although RPMD underestimates rate coefficients by a factor between 4 and 2 in the 200-500 K interval. The origin of these discrepancies lies in the large contribution from tunneling on the 3A″ PES, which is enhanced by resonances due to quasibound states in the van der Waals wells. Although tunneling is fairly well accounted for by RPMD even below the crossover temperature, the effect of resonances, a long-time effect, is not included in the methodology. At the highest temperatures studied in this work, 2000-3300 K, the RPMD rate coefficients are somewhat larger than the QM ones, but this is shown to be due to limitations in the QM calculations and the RPMD are believed to be more reliable.
Collapse
Affiliation(s)
- M Menéndez
- Departamento de Química Física I, Facultad de Ciencias Químicas , Universidad Complutense de Madrid , 28040 Madrid , Spain
| | - P G Jambrina
- Departamento de Química Física, Facultad de Ciencias Químicas , Universidad de Salamanca , 37008 Salamanca , Spain
| | - A Zanchet
- Departamento de Química Física I, Facultad de Ciencias Químicas , Universidad Complutense de Madrid , 28040 Madrid , Spain
| | - E Verdasco
- Departamento de Química Física I, Facultad de Ciencias Químicas , Universidad Complutense de Madrid , 28040 Madrid , Spain
| | - Y V Suleimanov
- Computation-based Science and Technology Research Center , Cyprus Institute , 20 Kavafi Strasse , Nicosia 2121 , Cyprus.,Department of Chemical Engineering , Massachusetts Institute of Technology , 77 Massachusetts Avenue , Cambridge , Massachusetts 02139 , United States
| | - F J Aoiz
- Departamento de Química Física I, Facultad de Ciencias Químicas , Universidad Complutense de Madrid , 28040 Madrid , Spain
| |
Collapse
|
34
|
Li J. Ring-polymer molecular dynamics studies of thermal rate coefficients for reaction F + H2O → HF + OH. CHINESE J CHEM PHYS 2019. [DOI: 10.1063/1674-0068/cjcp1808186] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
Affiliation(s)
- Jun Li
- School of Chemistry and Chemical Engineering, Chongqing University, Chongqing 400044, China
| |
Collapse
|
35
|
Lu X, Wang X, Fu B, Zhang D. Theoretical Investigations of Rate Coefficients of H + H2O2 → OH + H2O on a Full-Dimensional Potential Energy Surface. J Phys Chem A 2019; 123:3969-3976. [DOI: 10.1021/acs.jpca.9b02526] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Xiaoxiao Lu
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Zhongshan Road 457, Dalian 116023, China
- Department of Chemical Physics, University of Science and Technology of China, Jinzhai Road 96, Hefei 230026, China
| | - Xingan Wang
- Department of Chemical Physics, University of Science and Technology of China, Jinzhai Road 96, Hefei 230026, China
| | - 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, Zhongshan Road 457, Dalian 116023, China
| | - Donghui Zhang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Zhongshan Road 457, Dalian 116023, China
| |
Collapse
|
36
|
Thapa MJ, Fang W, Richardson JO. Nonadiabatic quantum transition-state theory in the golden-rule limit. I. Theory and application to model systems. J Chem Phys 2019; 150:104107. [DOI: 10.1063/1.5081108] [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)
- Manish J. Thapa
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | - Wei Fang
- Laboratory of Physical Chemistry, ETH Zürich, 8093 Zürich, Switzerland
| | | |
Collapse
|
37
|
Tao X, Shushkov P, Miller TF. Simple Flux-Side Formulation of State-Resolved Thermal Reaction Rates for Ring-Polymer Surface Hopping. J Phys Chem A 2019; 123:3013-3020. [DOI: 10.1021/acs.jpca.9b00877] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Xuecheng Tao
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Philip Shushkov
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Thomas F. Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| |
Collapse
|
38
|
Meng Q, Chen J. Ring-polymer molecular dynamics study on rate coefficient of the barrierless OH + CO system at low temperature. J Chem Phys 2019; 150:044307. [PMID: 30709288 DOI: 10.1063/1.5065657] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Qingyong Meng
- Department of Applied Chemistry, Northwestern Polytechnical University, Youyi West Road 127, 710072 Xi’an, China
| | - Jun Chen
- iChEM, College of Chemistry and Chemical Engineering, Xiamen University, Siming South Road 422, 361005 Xiamen, China
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Zhongshan Road 457, 116023 Dalian, China
| |
Collapse
|
39
|
Lu X, Meng Q, Wang X, Fu B, Zhang DH. Rate coefficients of the H + H2O2→ H2+ HO2reaction on an accurate fundamental invariant-neural network potential energy surface. J Chem Phys 2018; 149:174303. [PMID: 30409010 DOI: 10.1063/1.5063613] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Xiaoxiao Lu
- Department of Chemical Physics, University of Science and Technology of China, Jinzhai Road 96, Hefei 230026, China
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical and Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Zhongshan Road 457, Dalian 116023, China
| | - Qingyong Meng
- Department of Applied Chemistry, Northwestern Polytechnical University, Youyi West Road 127, Xi’an 710072, China
| | - Xingan Wang
- Department of Chemical Physics, University of Science and Technology of China, Jinzhai Road 96, Hefei 230026, China
| | - 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, Zhongshan Road 457, Dalian 116023, China
| | - 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, Zhongshan Road 457, Dalian 116023, China
| |
Collapse
|
40
|
Bhowmick S, Bossion D, Scribano Y, Suleimanov YV. The low temperature D + + H 2→ HD + H + reaction rate coefficient: a ring polymer molecular dynamics and quasi-classical trajectory study. Phys Chem Chem Phys 2018; 20:26752-26763. [PMID: 30324962 DOI: 10.1039/c8cp05398g] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The reaction between D+ and H2 plays an important role in astrochemistry at low temperatures and also serves as a prototype for a simple ion-molecule reaction. Its ground X[combining tilde]1A' state has a very small thermodynamic barrier (up to 1.8 × 10-2 eV) and the reaction proceeds through the formation of an intermediate complex lying within the potential well with a depth of at least 0.2 eV, thus representing a challenge for dynamical studies. In the present work, we analyze the title reaction within the temperature range of 20-100 K by means of ring polymer molecular dynamics (RPMD) and quasi-classical trajectory (QCT) methods over the full-dimensional global potential energy surface developed by Aguado et al. [A. Aguado, O. Roncero, C. Tablero, C. Sanz and M. Paniagua, J. Chem. Phys., 2000, 112, 1240]. The computed thermal RPMD and QCT rate coefficients are found to be almost independent of temperature and fall within the range of 1.34-2.01 × 10-9 cm3 s-1. They are also in very good agreement with previous time-independent quantum mechanical and statistical quantum method calculations. Furthermore, we observe that the choice of asymptotic separation distance between the reactants can markedly alter the rate coefficient in the low temperature regime (20-50 K). Therefore it is of utmost importance to correctly assign the value of this parameter for dynamical studies, particularly at very low temperatures of astrochemical importance. We finally conclude that the experimental rate measurements for the title reaction are highly desirable in future.
Collapse
Affiliation(s)
- Somnath Bhowmick
- Computation-based Science and Technology Research Center, The Cyprus Institute, 20 Konstantinou Kavafi Street, Nicosia 2121, Cyprus.
| | - Duncan Bossion
- Laboratoire Univers et Particules de Montpellier, UMR-CNRS 5299, Université de Montpellier, Place Eugène Bataillon, 34095 Montpellier, France
| | - Yohann Scribano
- Laboratoire Univers et Particules de Montpellier, UMR-CNRS 5299, Université de Montpellier, Place Eugène Bataillon, 34095 Montpellier, France
| | - Yury V Suleimanov
- Computation-based Science and Technology Research Center, The Cyprus Institute, 20 Konstantinou Kavafi Street, Nicosia 2121, Cyprus.
| |
Collapse
|
41
|
Meng Q. Ring-Polymer Molecular Dynamics with Coarse-Grained Treatment of the Rate Coefficients of Chlorine Atom Reactions with Methane, Ethane, and Propane. J Phys Chem A 2018; 122:8320-8325. [PMID: 30281311 DOI: 10.1021/acs.jpca.8b08052] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Qingyong Meng
- Department of Applied Chemistry, Northwestern Polytechnical University, Youyi West Road 127, 710072 Xi’an, China
- State Key Laboratory
of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Zhongshan Road 457, 116023 Dalian, China
| |
Collapse
|
42
|
Meng Q. Ring-polymer molecular dynamics study on rate coefficients of hydrogen abstraction of methane: A reduced-dimensional model. Chem Phys Lett 2018. [DOI: 10.1016/j.cplett.2018.06.042] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
43
|
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: 16] [Impact Index Per Article: 2.7] [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
|
44
|
Mazzuca JW, Haut NK. Theoretical description of quantum mechanical permeation of graphene membranes by charged hydrogen isotopes. J Chem Phys 2018; 148:224301. [PMID: 29907032 DOI: 10.1063/1.5027821] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
It has been recently shown that in the presence of an applied voltage, hydrogen and deuterium nuclei can be separated from one another using graphene membranes as a nuclear sieve, resulting in a 10-fold enhancement in the concentration of the lighter isotope. While previous studies, both experimental and theoretical, have attributed this effect mostly to differences in vibrational zero point energy (ZPE) of the various isotopes near the membrane surface, we propose that multi-dimensional quantum mechanical tunneling of nuclei through the graphene membrane influences this proton permeation process in a fundamental way. We perform ring polymer molecular dynamics calculations in which we include both ZPE and tunneling effects of various hydrogen isotopes as they permeate the graphene membrane and compute rate constants across a range of temperatures near 300 K. While capturing the experimentally observed separation factor, our calculations indicate that the transverse motion of the various isotopes across the surface of the graphene membrane is an essential part of this sieving mechanism. An understanding of the multi-dimensional quantum mechanical nature of this process could serve to guide the design of other such isotopic enrichment processes for a variety of atomic and molecular species of interest.
Collapse
Affiliation(s)
- James W Mazzuca
- Chemistry Department, Alma College, Alma, Michigan 48801, USA
| | | |
Collapse
|
45
|
Thompson KM, Gao Y, Marshall P, Wang H, Zhou L, Li Y, Guo H. Experimental and theoretical studies of the reactions of ground-state sulfur atoms with hydrogen and deuterium. J Chem Phys 2018; 147:134302. [PMID: 28987091 DOI: 10.1063/1.4991418] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The gas-phase kinetics of S(3P) atoms with H2 and D2 have been studied via the laser flash photolysis-resonance fluorescence technique. S atoms were generated by pulsed photolysis of CS2 at 193 nm and monitored by time-resolved fluorescence at 181 nm. The rate coefficients for H2 (k1) and D2 (k2), respectively, are summarized as k1(600-1110 K) = 3.0 × 10-9 exp-1.317×105-2.703×107K/T8.314 T/K cm3 molecule-1 s-1 and k2(770-1110 K) = 2.2 × 10-14 (T/298 K)3.55 exp(-5420 K/T) cm3 molecule-1 s-1. Error limits are discussed in the text. The rate coefficients for formation of SH(SD) + H(D) on a newly developed triplet potential energy surface were characterized via ring polymer molecular dynamics and canonical variational transition-state theory. There is excellent agreement above about 1000 K between theory and experiment. At lower temperatures, the experimental rate coefficient is substantially larger than the results computed for the adiabatic reaction, suggesting a significant role for intersystem crossing to the singlet potential energy surface at lower temperatures.
Collapse
Affiliation(s)
| | - Yide Gao
- Department of Chemistry, University of North Texas, Denton, Texas 76203, USA
| | - Paul Marshall
- Department of Chemistry, University of North Texas, Denton, Texas 76203, USA
| | - Han Wang
- Department of Physics, International Center of Quantum and Molecular Structures, and Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University, Shanghai 200444, China
| | - Linsen Zhou
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| | - 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
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, USA
| |
Collapse
|
46
|
Bishop KP, Roy PN. Quantum mechanical free energy profiles with post-quantization restraints: Binding free energy of the water dimer over a broad range of temperatures. J Chem Phys 2018; 148:102303. [DOI: 10.1063/1.4986915] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Kevin P. Bishop
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - Pierre-Nicholas Roy
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| |
Collapse
|
47
|
Nuñez-Reyes D, Hickson KM, Larrégaray P, Bonnet L, González-Lezana T, Suleimanov YV. A combined theoretical and experimental investigation of the kinetics and dynamics of the O( 1D) + D 2 reaction at low temperature. Phys Chem Chem Phys 2018; 20:4404-4414. [PMID: 29372194 DOI: 10.1039/c7cp07843a] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The O(1D) + H2 reaction is a prototype for simple atom-diatom insertion type mechanisms considered to involve deep potential wells. While exact quantum mechanical methods can be applied to describe the dynamics, such calculations are challenging given the numerous bound quantum states involved. Consequently, efforts have been made to develop alternative theoretical strategies to portray accurately the reactive process. Here we report an experimental and theoretical investigation of the O(1D) + D2 reaction over the 50-296 K range. The calculations employ three conceptually different approaches - mean potential phase space theory, the statistical quantum mechanical method and ring polymer molecular dynamics. The calculated rate constants are in excellent agreement over the entire temperature range, exhibiting only weak temperature dependence. The agreement between experiment and theory is also very good, with discrepancies smaller than 26%. Taken together, the present and previous theoretical results validate the hypothesis that long-lived complex formation dominates the reaction dynamics at low temperature.
Collapse
Affiliation(s)
- Dianailys Nuñez-Reyes
- Université de Bordeaux, Institut des Sciences Moléculaires, F-33400 Talence, France.
| | | | | | | | | | | |
Collapse
|
48
|
Laude G, Calderini D, Tew DP, Richardson JO. Ab initio instanton rate theory made efficient using Gaussian process regression. Faraday Discuss 2018; 212:237-258. [PMID: 30230495 DOI: 10.1039/c8fd00085a] [Citation(s) in RCA: 38] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
Ab initio instanton rate theory is a computational method for rigorously including tunnelling effects into the calculations of chemical reaction rates based on a potential-energy surface computed on the fly from electronic-structure theory. This approach is necessary to extend conventional transition-state theory into the deep-tunnelling regime, but it is also more computationally expensive as it requires many more ab initio calculations. We propose an approach which uses Gaussian process regression to fit the potential-energy surface locally around the dominant tunnelling pathway. The method can be converged to give the same result as from an on-the-fly ab initio instanton calculation but it requires far fewer electronic-structure calculations. This makes it a practical approach for obtaining accurate rate constants based on high-level electronic-structure methods. We show fast convergence to reproduce benchmark H + CH4 results and evaluate new low-temperature rates of H + C2H6 in full dimensionality at a UCCSD(T)-F12b/cc-pVTZ-F12 level.
Collapse
Affiliation(s)
- Gabriel Laude
- Laboratory of Physical Chemistry, ETH Zurich, Switzerland. and On exchange from School of Chemistry, University of Edinburgh, UK
| | | | - David P Tew
- Max-Planck-Institut für Festkörperforschung, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | | |
Collapse
|
49
|
Steffen J, Hartke B. Cheap but accurate calculation of chemical reaction rate constants from ab initio data, via system-specific, black-box force fields. J Chem Phys 2017; 147:161701. [DOI: 10.1063/1.4979712] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
50
|
Zuo J, Xie C, Guo H, Xie D. Accurate Determination of Tunneling-Affected Rate Coefficients: Theory Assessing Experiment. J Phys Chem Lett 2017; 8:3392-3397. [PMID: 28685568 DOI: 10.1021/acs.jpclett.7b01296] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
The thermal rate coefficients of a prototypical bimolecular reaction are determined on an accurate ab initio potential energy surface (PES) using ring polymer molecular dynamics (RPMD). It is shown that quantum effects such as tunneling and zero-point energy (ZPE) are of critical importance for the HCl + OH reaction at low temperatures, while the heavier deuterium substitution renders tunneling less facile in the DCl + OH reaction. The calculated RPMD rate coefficients are in excellent agreement with experimental data for the HCl + OH reaction in the entire temperature range of 200-1000 K, confirming the accuracy of the PES. On the other hand, the RPMD rate coefficients for the DCl + OH reaction agree with some, but not all, experimental values. The self-consistency of the theoretical results thus allows a quality assessment of the experimental data.
Collapse
Affiliation(s)
- Junxiang Zuo
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University , Nanjing 210093, China
| | - Changjian Xie
- Department of Chemistry and Chemical Biology, University of New Mexico , Albuquerque, New Mexico 87131, United States
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico , Albuquerque, New Mexico 87131, United States
| | - Daiqian Xie
- Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry, School of Chemistry and Chemical Engineering, Nanjing University , Nanjing 210093, China
- Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China , Hefei, Anhui 230026, China
| |
Collapse
|