Issue 31, 2024

Rate coefficients for the O + H2 and O + D2 reactions: how well ring polymer molecular dynamics accounts for tunelling

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.

Graphical abstract: Rate coefficients for the O + H2 and O + D2 reactions: how well ring polymer molecular dynamics accounts for tunelling

Supplementary files

Article information

Article type
Paper
Submitted
25 Apr 2024
Accepted
04 Jul 2024
First published
05 Jul 2024
This article is Open Access
Creative Commons BY-NC license

Phys. Chem. Chem. Phys., 2024,26, 20947-20961

Rate coefficients for the O + H2 and O + D2 reactions: how well ring polymer molecular dynamics accounts for tunelling

M. Menéndez, A. Veselinova, A. Zanchet, P. G. Jambrina and F. J. Aoiz, Phys. Chem. Chem. Phys., 2024, 26, 20947 DOI: 10.1039/D4CP01711K

This article is licensed under a Creative Commons Attribution-NonCommercial 3.0 Unported Licence. You can use material from this article in other publications, without requesting further permission from the RSC, provided that the correct acknowledgement is given and it is not used for commercial purposes.

To request permission to reproduce material from this article in a commercial publication, please go to the Copyright Clearance Center request page.

If you are an author contributing to an RSC publication, you do not need to request permission provided correct acknowledgement is given.

If you are the author of this article, you do not need to request permission to reproduce figures and diagrams provided correct acknowledgement is given. If you want to reproduce the whole article in a third-party commercial publication (excluding your thesis/dissertation for which permission is not required) please go to the Copyright Clearance Center request page.

Read more about how to correctly acknowledge RSC content.

Social activity

Spotlight

Advertisements