### Abstract

We have used the quasiclassical trajectory method and a diatomics-in-molecules potential energy surface to study the state-to-state dynamics of the H_{2}^{+} + H_{2} → H_{3}^{+} + H reaction, with an emphasis on the product H_{3}^{+} vibration/rotation distributions. The product vibration/rotation states were determined by a recently developed Fourier transform method for calculating the H_{3}^{+} good action variables. For H_{2}^{+} (v_{1} = 0, j_{1} = 2) + H_{2} (v_{2} = 0, j_{2} = 1), we find that over 70% of the available energy goes into H_{3}^{+} internal excitation. This result is only weakly dependent on reagent translational energy, but the fraction increases to over 80% when the reagent H_{2}^{+} is excited to v_{1} = 3. Approximately 20% of the available energy goes into product rotation at low translational energies, and larger amounts at higher energies. The distribution of H_{3}^{+} symmetric stretch and degenerate mode (bend plus asymmetric stretch) quantum states is found to be highly nonstatistical, with an average stretch quantum number N_{s} of about 0.6 and an average degenerate quantum number N_{d} of 2.4 for ground-state reagents. The values of N_{s} and N_{d} for the most probable final states are found to be strongly correlated, with N_{d} increasing nonlinearly with increasing N_{s}. The distribution of vibrational angular momentum l is found to show a strong propensity for l = 0, and we correlate this with the most favored orientation of the reagents along the pathway leading to reaction. The product rotational distributions were found to peak near J = 12 at low translational and vibrational energy, with K = 0 being the most likely projection state. This departure from the usual K = J propensity arises because strong vibration/rotation coupling in H_{3}^{+} causes the body-fixed Z projection to be oscillatory in time. Product angular distributions were found to be mainly controlled by a stripping reaction mechanism, but product translational distributions reflect the partitioning of the reaction exoergicity into product degrees of freedom. Where comparison is possible, the product translational, angular, and vibrational distributions agree with experiment. The total reaction cross section and its division into proton-transfer and atom-transfer contributions is also considered and compared with experiment by using a surface hopping model in which the dynamics is assumed to be diabatic up to a certain intermolecular separation and then adiabatic thereafter. We find that our total reactive cross sections are in reasonable agreement with experiment as a function of both reagent translational and vibrational energy. These cross sections are found to arise from roughly equal probabilities of proton and atom transfer.

Original language | English |
---|---|

Pages (from-to) | 2612-2620 |

Number of pages | 9 |

Journal | Journal of Physical Chemistry |

Volume | 89 |

Issue number | 12 |

Publication status | Published - 1985 |

### Fingerprint

### ASJC Scopus subject areas

- Physical and Theoretical Chemistry

### Cite this

_{2}

^{+}+ H

_{2}→ H

_{3}

^{+}+ H reaction.

*Journal of Physical Chemistry*,

*89*(12), 2612-2620.

**A quasiclassical trajectory study of the H _{2}^{+} + H_{2} → H_{3}^{+} + H reaction.** / Eaker, Charles W.; Schatz, George C.

Research output: Contribution to journal › Article

_{2}

^{+}+ H

_{2}→ H

_{3}

^{+}+ H reaction',

*Journal of Physical Chemistry*, vol. 89, no. 12, pp. 2612-2620.

_{2}

^{+}+ H

_{2}→ H

_{3}

^{+}+ H reaction. Journal of Physical Chemistry. 1985;89(12):2612-2620.

}

TY - JOUR

T1 - A quasiclassical trajectory study of the H2+ + H2 → H3+ + H reaction

AU - Eaker, Charles W.

AU - Schatz, George C

PY - 1985

Y1 - 1985

N2 - We have used the quasiclassical trajectory method and a diatomics-in-molecules potential energy surface to study the state-to-state dynamics of the H2+ + H2 → H3+ + H reaction, with an emphasis on the product H3+ vibration/rotation distributions. The product vibration/rotation states were determined by a recently developed Fourier transform method for calculating the H3+ good action variables. For H2+ (v1 = 0, j1 = 2) + H2 (v2 = 0, j2 = 1), we find that over 70% of the available energy goes into H3+ internal excitation. This result is only weakly dependent on reagent translational energy, but the fraction increases to over 80% when the reagent H2+ is excited to v1 = 3. Approximately 20% of the available energy goes into product rotation at low translational energies, and larger amounts at higher energies. The distribution of H3+ symmetric stretch and degenerate mode (bend plus asymmetric stretch) quantum states is found to be highly nonstatistical, with an average stretch quantum number Ns of about 0.6 and an average degenerate quantum number Nd of 2.4 for ground-state reagents. The values of Ns and Nd for the most probable final states are found to be strongly correlated, with Nd increasing nonlinearly with increasing Ns. The distribution of vibrational angular momentum l is found to show a strong propensity for l = 0, and we correlate this with the most favored orientation of the reagents along the pathway leading to reaction. The product rotational distributions were found to peak near J = 12 at low translational and vibrational energy, with K = 0 being the most likely projection state. This departure from the usual K = J propensity arises because strong vibration/rotation coupling in H3+ causes the body-fixed Z projection to be oscillatory in time. Product angular distributions were found to be mainly controlled by a stripping reaction mechanism, but product translational distributions reflect the partitioning of the reaction exoergicity into product degrees of freedom. Where comparison is possible, the product translational, angular, and vibrational distributions agree with experiment. The total reaction cross section and its division into proton-transfer and atom-transfer contributions is also considered and compared with experiment by using a surface hopping model in which the dynamics is assumed to be diabatic up to a certain intermolecular separation and then adiabatic thereafter. We find that our total reactive cross sections are in reasonable agreement with experiment as a function of both reagent translational and vibrational energy. These cross sections are found to arise from roughly equal probabilities of proton and atom transfer.

AB - We have used the quasiclassical trajectory method and a diatomics-in-molecules potential energy surface to study the state-to-state dynamics of the H2+ + H2 → H3+ + H reaction, with an emphasis on the product H3+ vibration/rotation distributions. The product vibration/rotation states were determined by a recently developed Fourier transform method for calculating the H3+ good action variables. For H2+ (v1 = 0, j1 = 2) + H2 (v2 = 0, j2 = 1), we find that over 70% of the available energy goes into H3+ internal excitation. This result is only weakly dependent on reagent translational energy, but the fraction increases to over 80% when the reagent H2+ is excited to v1 = 3. Approximately 20% of the available energy goes into product rotation at low translational energies, and larger amounts at higher energies. The distribution of H3+ symmetric stretch and degenerate mode (bend plus asymmetric stretch) quantum states is found to be highly nonstatistical, with an average stretch quantum number Ns of about 0.6 and an average degenerate quantum number Nd of 2.4 for ground-state reagents. The values of Ns and Nd for the most probable final states are found to be strongly correlated, with Nd increasing nonlinearly with increasing Ns. The distribution of vibrational angular momentum l is found to show a strong propensity for l = 0, and we correlate this with the most favored orientation of the reagents along the pathway leading to reaction. The product rotational distributions were found to peak near J = 12 at low translational and vibrational energy, with K = 0 being the most likely projection state. This departure from the usual K = J propensity arises because strong vibration/rotation coupling in H3+ causes the body-fixed Z projection to be oscillatory in time. Product angular distributions were found to be mainly controlled by a stripping reaction mechanism, but product translational distributions reflect the partitioning of the reaction exoergicity into product degrees of freedom. Where comparison is possible, the product translational, angular, and vibrational distributions agree with experiment. The total reaction cross section and its division into proton-transfer and atom-transfer contributions is also considered and compared with experiment by using a surface hopping model in which the dynamics is assumed to be diabatic up to a certain intermolecular separation and then adiabatic thereafter. We find that our total reactive cross sections are in reasonable agreement with experiment as a function of both reagent translational and vibrational energy. These cross sections are found to arise from roughly equal probabilities of proton and atom transfer.

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

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

M3 - Article

AN - SCOPUS:3342899505

VL - 89

SP - 2612

EP - 2620

JO - Journal of Physical Chemistry

JF - Journal of Physical Chemistry

SN - 0022-3654

IS - 12

ER -