An improved potential energy surface and multi-temperature quasiclassical trajectory calculations of N2 + N2 dissociation reactions

Jason D. Bender, Paolo Valentini, Ioannis Nompelis, Yuliya Paukku, Z. Varga, Donald G. Truhlar, Thomas Schwartzentruber, Graham V. Candler

Research output: Contribution to journalArticle

82 Citations (Scopus)

Abstract

Accurate modeling of high-temperature hypersonic flows in the atmosphere requires consideration of collision-induced dissociation of molecular species and energy transfer between the translational and internal modes of the gas molecules. Here, we describe a study of the N2 + N2→N2 + 2N and N2 + N2→4N nitrogen dissociation reactions using the quasiclassical trajectory (QCT) method. The simulations used a new potential energy surface for the N4 system; the surface is an improved version of one that was presented previously. In the QCT calculations, initial conditions were determined based on a two-temperature model that approximately separates the translational-rotational temperature from the vibrational temperature of the N2 diatoms. Five values from 8000 K to 30 000 K were considered for each of the two temperatures. Over 2.4 × 109 trajectories were calculated. We present results for ensemble-averaged dissociation rate constants as functions of the translational-rotational temperature T and the vibrational temperature Tv. The rate constant depends more strongly on T when Tv is low, and it depends more strongly on Tv when T is low. Quasibound reactant states contribute significantly to the rate constants, as do exchange processes at higher temperatures. We discuss two sets of runs in detail: an equilibrium test set in which T = Tv and a nonequilibrium test set in which Tv <T. In the equilibrium test set, high-v and moderately-low-j molecules contribute most significantly to the overall dissociation rate, and this state specificity becomes stronger as the temperature decreases. Dissociating trajectories tend to result in a major loss of vibrational energy and a minor loss of rotational energy. In the nonequilibrium test set, as Tv decreases while T is fixed, higher-j molecules contribute more significantly to the dissociation rate, dissociating trajectories tend to result in a greater rotational energy loss, and the dissociation probability's dependence on v weakens. In this way, as Tv decreases, rotational energy appears to compensate for the decline in average vibrational energy in promoting dissociation. In both the equilibrium and nonequilibrium test sets, in every case, the average total internal energy loss in the dissociating trajectories is between 10.2 and 11.0 eV, slightly larger than the equilibrium potential energy change of N2 dissociation.

Original languageEnglish
Article number054304
JournalThe Journal of Chemical Physics
Volume143
Issue number5
DOIs
Publication statusPublished - Aug 7 2015

Fingerprint

Potential energy surfaces
potential energy
Trajectories
trajectories
dissociation
Temperature
temperature
Rate constants
Molecules
energy dissipation
Energy dissipation
molecules
hypersonic flow
energy
Hypersonic flow
algae
internal energy
Potential energy
Energy transfer
energy transfer

ASJC Scopus subject areas

  • Physics and Astronomy(all)
  • Physical and Theoretical Chemistry

Cite this

An improved potential energy surface and multi-temperature quasiclassical trajectory calculations of N2 + N2 dissociation reactions. / Bender, Jason D.; Valentini, Paolo; Nompelis, Ioannis; Paukku, Yuliya; Varga, Z.; Truhlar, Donald G.; Schwartzentruber, Thomas; Candler, Graham V.

In: The Journal of Chemical Physics, Vol. 143, No. 5, 054304, 07.08.2015.

Research output: Contribution to journalArticle

Bender, Jason D. ; Valentini, Paolo ; Nompelis, Ioannis ; Paukku, Yuliya ; Varga, Z. ; Truhlar, Donald G. ; Schwartzentruber, Thomas ; Candler, Graham V. / An improved potential energy surface and multi-temperature quasiclassical trajectory calculations of N2 + N2 dissociation reactions. In: The Journal of Chemical Physics. 2015 ; Vol. 143, No. 5.
@article{6f9a625badec49438e7bf2e1eaeb07ae,
title = "An improved potential energy surface and multi-temperature quasiclassical trajectory calculations of N2 + N2 dissociation reactions",
abstract = "Accurate modeling of high-temperature hypersonic flows in the atmosphere requires consideration of collision-induced dissociation of molecular species and energy transfer between the translational and internal modes of the gas molecules. Here, we describe a study of the N2 + N2→N2 + 2N and N2 + N2→4N nitrogen dissociation reactions using the quasiclassical trajectory (QCT) method. The simulations used a new potential energy surface for the N4 system; the surface is an improved version of one that was presented previously. In the QCT calculations, initial conditions were determined based on a two-temperature model that approximately separates the translational-rotational temperature from the vibrational temperature of the N2 diatoms. Five values from 8000 K to 30 000 K were considered for each of the two temperatures. Over 2.4 × 109 trajectories were calculated. We present results for ensemble-averaged dissociation rate constants as functions of the translational-rotational temperature T and the vibrational temperature Tv. The rate constant depends more strongly on T when Tv is low, and it depends more strongly on Tv when T is low. Quasibound reactant states contribute significantly to the rate constants, as do exchange processes at higher temperatures. We discuss two sets of runs in detail: an equilibrium test set in which T = Tv and a nonequilibrium test set in which Tv <T. In the equilibrium test set, high-v and moderately-low-j molecules contribute most significantly to the overall dissociation rate, and this state specificity becomes stronger as the temperature decreases. Dissociating trajectories tend to result in a major loss of vibrational energy and a minor loss of rotational energy. In the nonequilibrium test set, as Tv decreases while T is fixed, higher-j molecules contribute more significantly to the dissociation rate, dissociating trajectories tend to result in a greater rotational energy loss, and the dissociation probability's dependence on v weakens. In this way, as Tv decreases, rotational energy appears to compensate for the decline in average vibrational energy in promoting dissociation. In both the equilibrium and nonequilibrium test sets, in every case, the average total internal energy loss in the dissociating trajectories is between 10.2 and 11.0 eV, slightly larger than the equilibrium potential energy change of N2 dissociation.",
author = "Bender, {Jason D.} and Paolo Valentini and Ioannis Nompelis and Yuliya Paukku and Z. Varga and Truhlar, {Donald G.} and Thomas Schwartzentruber and Candler, {Graham V.}",
year = "2015",
month = "8",
day = "7",
doi = "10.1063/1.4927571",
language = "English",
volume = "143",
journal = "Journal of Chemical Physics",
issn = "0021-9606",
publisher = "American Institute of Physics Publising LLC",
number = "5",

}

TY - JOUR

T1 - An improved potential energy surface and multi-temperature quasiclassical trajectory calculations of N2 + N2 dissociation reactions

AU - Bender, Jason D.

AU - Valentini, Paolo

AU - Nompelis, Ioannis

AU - Paukku, Yuliya

AU - Varga, Z.

AU - Truhlar, Donald G.

AU - Schwartzentruber, Thomas

AU - Candler, Graham V.

PY - 2015/8/7

Y1 - 2015/8/7

N2 - Accurate modeling of high-temperature hypersonic flows in the atmosphere requires consideration of collision-induced dissociation of molecular species and energy transfer between the translational and internal modes of the gas molecules. Here, we describe a study of the N2 + N2→N2 + 2N and N2 + N2→4N nitrogen dissociation reactions using the quasiclassical trajectory (QCT) method. The simulations used a new potential energy surface for the N4 system; the surface is an improved version of one that was presented previously. In the QCT calculations, initial conditions were determined based on a two-temperature model that approximately separates the translational-rotational temperature from the vibrational temperature of the N2 diatoms. Five values from 8000 K to 30 000 K were considered for each of the two temperatures. Over 2.4 × 109 trajectories were calculated. We present results for ensemble-averaged dissociation rate constants as functions of the translational-rotational temperature T and the vibrational temperature Tv. The rate constant depends more strongly on T when Tv is low, and it depends more strongly on Tv when T is low. Quasibound reactant states contribute significantly to the rate constants, as do exchange processes at higher temperatures. We discuss two sets of runs in detail: an equilibrium test set in which T = Tv and a nonequilibrium test set in which Tv <T. In the equilibrium test set, high-v and moderately-low-j molecules contribute most significantly to the overall dissociation rate, and this state specificity becomes stronger as the temperature decreases. Dissociating trajectories tend to result in a major loss of vibrational energy and a minor loss of rotational energy. In the nonequilibrium test set, as Tv decreases while T is fixed, higher-j molecules contribute more significantly to the dissociation rate, dissociating trajectories tend to result in a greater rotational energy loss, and the dissociation probability's dependence on v weakens. In this way, as Tv decreases, rotational energy appears to compensate for the decline in average vibrational energy in promoting dissociation. In both the equilibrium and nonequilibrium test sets, in every case, the average total internal energy loss in the dissociating trajectories is between 10.2 and 11.0 eV, slightly larger than the equilibrium potential energy change of N2 dissociation.

AB - Accurate modeling of high-temperature hypersonic flows in the atmosphere requires consideration of collision-induced dissociation of molecular species and energy transfer between the translational and internal modes of the gas molecules. Here, we describe a study of the N2 + N2→N2 + 2N and N2 + N2→4N nitrogen dissociation reactions using the quasiclassical trajectory (QCT) method. The simulations used a new potential energy surface for the N4 system; the surface is an improved version of one that was presented previously. In the QCT calculations, initial conditions were determined based on a two-temperature model that approximately separates the translational-rotational temperature from the vibrational temperature of the N2 diatoms. Five values from 8000 K to 30 000 K were considered for each of the two temperatures. Over 2.4 × 109 trajectories were calculated. We present results for ensemble-averaged dissociation rate constants as functions of the translational-rotational temperature T and the vibrational temperature Tv. The rate constant depends more strongly on T when Tv is low, and it depends more strongly on Tv when T is low. Quasibound reactant states contribute significantly to the rate constants, as do exchange processes at higher temperatures. We discuss two sets of runs in detail: an equilibrium test set in which T = Tv and a nonequilibrium test set in which Tv <T. In the equilibrium test set, high-v and moderately-low-j molecules contribute most significantly to the overall dissociation rate, and this state specificity becomes stronger as the temperature decreases. Dissociating trajectories tend to result in a major loss of vibrational energy and a minor loss of rotational energy. In the nonequilibrium test set, as Tv decreases while T is fixed, higher-j molecules contribute more significantly to the dissociation rate, dissociating trajectories tend to result in a greater rotational energy loss, and the dissociation probability's dependence on v weakens. In this way, as Tv decreases, rotational energy appears to compensate for the decline in average vibrational energy in promoting dissociation. In both the equilibrium and nonequilibrium test sets, in every case, the average total internal energy loss in the dissociating trajectories is between 10.2 and 11.0 eV, slightly larger than the equilibrium potential energy change of N2 dissociation.

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

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

U2 - 10.1063/1.4927571

DO - 10.1063/1.4927571

M3 - Article

AN - SCOPUS:84938808775

VL - 143

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 5

M1 - 054304

ER -