Quantum Wave Packet Study of Nonadiabatic Effects in O(1D) + H2 → OH + H

Stephen K. Gray, Carlo Petrongolo, Karen Drukker, George C Schatz

Research output: Contribution to journalArticle

68 Citations (Scopus)

Abstract

We develop a wave packet approach to treating the electronically nonadiabatic reaction dynamics of O(1D) + H2 → OH + H, allowing for the 11A′ and 21A′ potential energy surfaces and couplings, as well as the three internal nuclear coordinates. Two different systems of coupled potential energy surfaces are considered, a semiempirical diatomics-in-molecules (DIM) system due to Kuntz, Niefer, and Sloan, and a recently developed ab initio system due to Dobbyn and Knowles (DK). Nonadiabatic quantum results, with total angular momentum J = 0, are obtained and discussed. Several single surface calculations are carried out for comparison with the nonadiabatic results. Comparisons with trajectory surface hopping (TSH) calculations, and with approximate quantum calculations, are also included. The electrostatic coupling produces strong interactions between the 11A′ and 21A′ states at short range (where these states have a conical intersection) and weak but, interestingly, nonnegligible interactions between these states at longer range. Our wave packet results show that if the initial state is chosen to be effectively the 1A′ state (for which insertion to form products occurs on the adiabatic surface), then there is very little difference between the adiabatic and coupled surface results. In either case the reaction probability is a relatively flat function of energy, except for resonant oscillations. However, the 2A′ reaction, dynamics (which involves a collinear transition state) is strongly perturbed by nonadiabatic effects in two distinct ways. At energies above the transition state barrier, the diabatic limit is dominant, and the 2A′ reaction probability is similar to that for 1A″, which has no coupling with the other surfaces. At energies below the barrier, we find a significant component of the reaction probability from long range electronic coupling that effectively allows the wave packet to avoid having to tunnel through the barrier. This effect, which is observed on both the DIM and DK surfaces, is estimated to cause a 10% contribution to the room temperature rate constant from nonadiabatic effects. Similar results are obtained from the TSH and approximate quantum calculations.

Original languageEnglish
Pages (from-to)9448-9459
Number of pages12
JournalJournal of Physical Chemistry A
Volume103
Issue number47
Publication statusPublished - Nov 25 1999

Fingerprint

Wave packets
wave packets
Potential energy surfaces
Trajectories
Molecules
potential energy
Angular momentum
trajectories
Electrostatics
Rate constants
Tunnels
intersections
tunnels
energy
molecules
insertion
angular momentum
electrostatics
oscillations
causes

ASJC Scopus subject areas

  • Physical and Theoretical Chemistry

Cite this

Quantum Wave Packet Study of Nonadiabatic Effects in O(1D) + H2 → OH + H. / Gray, Stephen K.; Petrongolo, Carlo; Drukker, Karen; Schatz, George C.

In: Journal of Physical Chemistry A, Vol. 103, No. 47, 25.11.1999, p. 9448-9459.

Research output: Contribution to journalArticle

Gray, SK, Petrongolo, C, Drukker, K & Schatz, GC 1999, 'Quantum Wave Packet Study of Nonadiabatic Effects in O(1D) + H2 → OH + H', Journal of Physical Chemistry A, vol. 103, no. 47, pp. 9448-9459.
Gray, Stephen K. ; Petrongolo, Carlo ; Drukker, Karen ; Schatz, George C. / Quantum Wave Packet Study of Nonadiabatic Effects in O(1D) + H2 → OH + H. In: Journal of Physical Chemistry A. 1999 ; Vol. 103, No. 47. pp. 9448-9459.
@article{bcb5a1ac40cf41568c21accc215cf88d,
title = "Quantum Wave Packet Study of Nonadiabatic Effects in O(1D) + H2 → OH + H",
abstract = "We develop a wave packet approach to treating the electronically nonadiabatic reaction dynamics of O(1D) + H2 → OH + H, allowing for the 11A′ and 21A′ potential energy surfaces and couplings, as well as the three internal nuclear coordinates. Two different systems of coupled potential energy surfaces are considered, a semiempirical diatomics-in-molecules (DIM) system due to Kuntz, Niefer, and Sloan, and a recently developed ab initio system due to Dobbyn and Knowles (DK). Nonadiabatic quantum results, with total angular momentum J = 0, are obtained and discussed. Several single surface calculations are carried out for comparison with the nonadiabatic results. Comparisons with trajectory surface hopping (TSH) calculations, and with approximate quantum calculations, are also included. The electrostatic coupling produces strong interactions between the 11A′ and 21A′ states at short range (where these states have a conical intersection) and weak but, interestingly, nonnegligible interactions between these states at longer range. Our wave packet results show that if the initial state is chosen to be effectively the 1A′ state (for which insertion to form products occurs on the adiabatic surface), then there is very little difference between the adiabatic and coupled surface results. In either case the reaction probability is a relatively flat function of energy, except for resonant oscillations. However, the 2A′ reaction, dynamics (which involves a collinear transition state) is strongly perturbed by nonadiabatic effects in two distinct ways. At energies above the transition state barrier, the diabatic limit is dominant, and the 2A′ reaction probability is similar to that for 1A″, which has no coupling with the other surfaces. At energies below the barrier, we find a significant component of the reaction probability from long range electronic coupling that effectively allows the wave packet to avoid having to tunnel through the barrier. This effect, which is observed on both the DIM and DK surfaces, is estimated to cause a 10{\%} contribution to the room temperature rate constant from nonadiabatic effects. Similar results are obtained from the TSH and approximate quantum calculations.",
author = "Gray, {Stephen K.} and Carlo Petrongolo and Karen Drukker and Schatz, {George C}",
year = "1999",
month = "11",
day = "25",
language = "English",
volume = "103",
pages = "9448--9459",
journal = "Journal of Physical Chemistry A",
issn = "1089-5639",
publisher = "American Chemical Society",
number = "47",

}

TY - JOUR

T1 - Quantum Wave Packet Study of Nonadiabatic Effects in O(1D) + H2 → OH + H

AU - Gray, Stephen K.

AU - Petrongolo, Carlo

AU - Drukker, Karen

AU - Schatz, George C

PY - 1999/11/25

Y1 - 1999/11/25

N2 - We develop a wave packet approach to treating the electronically nonadiabatic reaction dynamics of O(1D) + H2 → OH + H, allowing for the 11A′ and 21A′ potential energy surfaces and couplings, as well as the three internal nuclear coordinates. Two different systems of coupled potential energy surfaces are considered, a semiempirical diatomics-in-molecules (DIM) system due to Kuntz, Niefer, and Sloan, and a recently developed ab initio system due to Dobbyn and Knowles (DK). Nonadiabatic quantum results, with total angular momentum J = 0, are obtained and discussed. Several single surface calculations are carried out for comparison with the nonadiabatic results. Comparisons with trajectory surface hopping (TSH) calculations, and with approximate quantum calculations, are also included. The electrostatic coupling produces strong interactions between the 11A′ and 21A′ states at short range (where these states have a conical intersection) and weak but, interestingly, nonnegligible interactions between these states at longer range. Our wave packet results show that if the initial state is chosen to be effectively the 1A′ state (for which insertion to form products occurs on the adiabatic surface), then there is very little difference between the adiabatic and coupled surface results. In either case the reaction probability is a relatively flat function of energy, except for resonant oscillations. However, the 2A′ reaction, dynamics (which involves a collinear transition state) is strongly perturbed by nonadiabatic effects in two distinct ways. At energies above the transition state barrier, the diabatic limit is dominant, and the 2A′ reaction probability is similar to that for 1A″, which has no coupling with the other surfaces. At energies below the barrier, we find a significant component of the reaction probability from long range electronic coupling that effectively allows the wave packet to avoid having to tunnel through the barrier. This effect, which is observed on both the DIM and DK surfaces, is estimated to cause a 10% contribution to the room temperature rate constant from nonadiabatic effects. Similar results are obtained from the TSH and approximate quantum calculations.

AB - We develop a wave packet approach to treating the electronically nonadiabatic reaction dynamics of O(1D) + H2 → OH + H, allowing for the 11A′ and 21A′ potential energy surfaces and couplings, as well as the three internal nuclear coordinates. Two different systems of coupled potential energy surfaces are considered, a semiempirical diatomics-in-molecules (DIM) system due to Kuntz, Niefer, and Sloan, and a recently developed ab initio system due to Dobbyn and Knowles (DK). Nonadiabatic quantum results, with total angular momentum J = 0, are obtained and discussed. Several single surface calculations are carried out for comparison with the nonadiabatic results. Comparisons with trajectory surface hopping (TSH) calculations, and with approximate quantum calculations, are also included. The electrostatic coupling produces strong interactions between the 11A′ and 21A′ states at short range (where these states have a conical intersection) and weak but, interestingly, nonnegligible interactions between these states at longer range. Our wave packet results show that if the initial state is chosen to be effectively the 1A′ state (for which insertion to form products occurs on the adiabatic surface), then there is very little difference between the adiabatic and coupled surface results. In either case the reaction probability is a relatively flat function of energy, except for resonant oscillations. However, the 2A′ reaction, dynamics (which involves a collinear transition state) is strongly perturbed by nonadiabatic effects in two distinct ways. At energies above the transition state barrier, the diabatic limit is dominant, and the 2A′ reaction probability is similar to that for 1A″, which has no coupling with the other surfaces. At energies below the barrier, we find a significant component of the reaction probability from long range electronic coupling that effectively allows the wave packet to avoid having to tunnel through the barrier. This effect, which is observed on both the DIM and DK surfaces, is estimated to cause a 10% contribution to the room temperature rate constant from nonadiabatic effects. Similar results are obtained from the TSH and approximate quantum calculations.

UR - http://www.scopus.com/inward/record.url?scp=0001124375&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0001124375&partnerID=8YFLogxK

M3 - Article

AN - SCOPUS:0001124375

VL - 103

SP - 9448

EP - 9459

JO - Journal of Physical Chemistry A

JF - Journal of Physical Chemistry A

SN - 1089-5639

IS - 47

ER -