Three‐dimensional DIM‐3C potential energy surfaces for the reactions H+XY and X+HY (X,Y=F, Cl, Br, I)

New Potential Energy Surface for the H + Cl2 Reaction and Quantum Dynamics Studies

The reaction of H + Cl2 → HCl + Cl plays a crucial role in various fields. However, no previous study has investigated this reaction using accurate quantum mechanical methods. In this paper, we construct a global potential energy surface (PES) using the neural network method with more than 20,000 ab initio energies obtained by the MRCI-F12+Q method with the aug-cc-pV5Z basis and extrapolated to the complete basis set limit. The spin-orbit coupling of the Cl atom has been considered in the PES. With this new PES, product state-resolved quantum dynamics calculations for the H + Cl2 (v0 = 0, j0 = 0-2) → HCl + Cl reaction was carried out. Numerical results show that the initial rotational excitation of the Cl2 has negligible effects on the reactivity. Product state-resolved integral cross sections (ICS) and rate constants reveal that the HCl is most favorably produced in its v' = 2 vibrational state. The calculated product vibrational state-resolved and total reaction rate constants suggest that the new global PES is accurate enough, as compared with the available experimental measurements.

Ring-polymer molecular dynamical calculations for the F + HCl → HF + Cl reaction on the ground 12A' potential energy surface

The reaction kinetics of the heavy-light-heavy abstraction reaction F + HCl → HF + Cl on the ground electronic state potential energy surface (PES) is investigated theoretically by a recently developed ring polymer molecular dynamics (RPMD) approach. First, a new PES is developed by the permutation invariant polynomial neural network (PIP-NN) approach based on 30 620 points sampled over a large configuration space from the latest and most accurate Deskevich-Hayes-Takahashi-Skodje-Nesbitt (DHTSN) PES (J. Chem. Phys., 2006, 124, 224303). Excellent fitting performance was achieved with only 521 parameters. The PIP-NN PES is 11 times faster than the DHTSN PES. Besides, the first analytical derivatives with respect to the coordinates of the atoms have been obtained for the PIP-NN PES. The RPMD rate coefficients on the PIP-NN PES are calculated and compared with available theoretical and experimental results. It is found that the experimental rate coefficients are significantly larger than the theoretical results on the DHTSN PES, due to its overestimated reaction barrier. We conclude that a reliable PES for this important heavy-light-heavy reaction is highly desirable.

Time-dependent wave packet quantum scattering and quasi-classical trajectory calculations of the H + FCl(v=0,j=0) → HF + Cl/HCl + F reaction

The dynamics of the title reaction are investigated using both time-dependent wave packet quantum scattering and quasi-classical trajectory (QCT) methods on adiabatic ground 1(2)A' potential energy surface (PES). Compared with the quantum results of reaction probabilities of H + FCl(J=0) → HF + Cl/HCl + F, the QCT method is proven feasible and further employed to produce integral cross sections and rate constants. Significant resonance structures are observed in the reaction probabilities using the quantum method; however, there are some undulations in the calculated QCT integral cross sections for both product channels. A comparison between the quantum mechanical coupled-channel (CC) calculation and centrifugal sudden approximation calculation reveals the very important role of Coriolis coupling effects in the quantum calculation. Comparisons between the calculated thermal rate constants for both reactions and the previous theoretical and experimental results have been done. HCl product formation is favored over the HF product in the reactive system. Finally, the HF products are found to be mainly forward scattering, and the HCl products are mainly backward scattering.

PRODUCT ROTATIONAL ANGULAR MOMENTUM POLARIZATION IN REACTIONS D + FCl(v = 0,j = 0) → DCl + F, DF + Cl

A comparative quasi-classical trajectory study on the product alignment and orientation of reactions D + FCl (v,j) → DCl + F , DF + Cl are carried out on a recently computed 12A′ ground-state surface reported by Deskevich et al. The reaction probabilities, integral and differential cross-sections as a function of collision energy are presented. The differential cross-sections and 〈P2(j′ ⋅ k)〉 are governed by different reaction mechanisms for these two reactions. P(θr) and P(ϕr) distributions of products DCl and DF at four selected collision energies of 5, 10, 16 and 30 kcal/mol indicate that DCl and DF product molecules are not only aligned, but also oriented along y-axis. By comparing the P(θr) and P(ϕr) distributions of reactions H + FCl (v = 0,j = 0) → HCl + F and H + FCl (v = 0,j = 0) → HF + Cl , mass factor has some influence on stereodynamics of the title reactions, but the features of PES should have major influence on the difference between dynamics of reactions D + FCl (v = 0,j = 0) → DCl + F and D + FCl (v = 0,j = 0) → DF + Cl .

Quantum state resolved scattering dynamics of F+HCl→HF(v,J)+Cl

State-to-state reaction dynamics of the reaction F+HCl-->HF(v,J)+Cl have been studied under single-collision conditions using an intense discharge F atom source in crossed supersonic molecular beams at Ecom=4.3(1.3) kcal/mol. Nascent HF product is monitored by shot-noise limited direct infrared laser absorption, providing quantum state distributions as well as additional information on kinetic energy release from high resolution Dopplerimetry. The vibrational distributions are highly inverted, with 34(4)%, 44(2)%, and 8(1)% of the total population in vHF=1, 2, and 3, respectively, consistent with predominant energy release into the newly formed bond. However, there is a small [14(1)%] but significant formation channel into the vHF=0 ground state, which is directly detectable for the first time via direct absorption methods. Of particular dynamical interest, both the HF(v=2,J) and HF(v=1,J) populations exhibit strongly bimodal J distributions. These results differ significantly from previous flow and arrested-relaxation studies and may signal the presence of microscopic branching in the reaction dynamics.

Barriers of Hydrogen Abstraction vs Halogen Exchange: An Experimental Manifestation of Charge-Shift Bonding

This paper shows that the differences between the barriers of the halogen exchange reactions, in the H + XH systems, and the hydrogen abstraction reactions, in the X + HX systems (X = F, Cl, Br), measure the covalent-ionic resonance energies of the corresponding X-H bonds. These processes are investigated using CCSD(T) calculations as well as the breathing-orbital valence bond (BOVB) method. Thus, the VB analysis shows that (i) at the level of covalent structures the barriers are the same for the two series and (ii) the higher barriers for halogen exchange processes originate solely from the less efficient mixing of the ionic structures into the respective covalent structures. The barrier differences, in the HXH vs XHX series, which decrease as X is varied from F to I, can be estimated as one-quarter of the covalent-ionic resonance energy of the H-X bond. The largest difference (22 kcal/mol) is calculated for X = F in accord with the finding that the H-F bond possesses the largest covalent-ionic resonance energy, 87 kcal/mol, which constitutes the major part of the bonding energy. The H-F bond belongs to the class of "charge-shift" bonds (Shaik, S.; Danovich, D.; Silvi, B.; Lauvergnat, D. L.; Hiberty, P. C. Chem. Eur. J. 2005, 21, 6358), which are all typified by dominant covalent-ionic resonance energies. Since the barrier difference between the two series is an experimental measure of the resonance energy quantity, in the particular case of X = F, the unusually high barrier for the fluorine exchange reaction emerges as an experimental manifestation of charge-shift bonding.

A Simple Picture for the Rotational Enhancement of the Rate for the F + HCl → HF + Cl Reaction: A Dynamical Study Using a New ab initio Potential Energy Surface

Quantum scattering calculations for the reaction F + HCl --> HF + Cl are performed on a new ground-state ab initio potential energy surface. The reagent rotation is found to have a dramatic effect on the reaction probability. Furthermore, the exit channel rotational thresholds leave a strong imprint on the reaction probabilities and even on the cumulative reaction probability. A very simple vibrationally adiabatic model is shown to account for most aspects of the reaction dynamics. In this model, the fast vibrational motion is adiabatically eliminated leaving the key reaction dynamics represented by a reduced atom + rotor collision. The shape of the adiabatic potential surface immediately yields to a simple and intuitive interpretation for the rotational enhancement of the rate. The rotational enhancement is shown to be an effect of the entrance channel dynamics of the atom-rotor problem.