Quasiclassical trajectory study of the reaction NH(X~3∑~-)+H→N(~4S)+H_2

The dynamics of the NH + H→N+H2 reaction has been investigated by means of the 3D quasiclassical trajectory approach by using the LEPS potential energy surface.The calculated rate coefficient is in good agreement with the experimental value.The reaction was found to occur via a direct channel.The product H2 has a cold excitation of rotational state,but has a reverse distribution of the vibrational state with a peak at v=1.Based on the potential energy surface and the trajectory analysis,the reaction mechanism has been explained successfully.     

Stereo‐dynamics of the F + HCl → HF + Cl reaction

We present a quasi‐classical trajectory (QCT) study on product polarization for the reaction F(²P) + HCl(v = 0, j = 0) → HF + Cl(²P) on a recently computed 1² A′ ground‐state surface reported by Deskevich et al. J Chem Phys, 2006, 124, 224303. Four polarization dependent generalized differential cross‐sections (2π/σ)(dσ₀₀/dω_t), (2π/σ)(dσ₂₀/dω_t), (2π/σ)(dσ₂₂₊/dω_t), and (2π/σ)(dσ_21?/dω_t) were calculated in the center‐of‐mass frame at four different collision energies. The obtained P(θ_r), P(?_r), and P(θ_r, ?_r), which denote respectively the distribution of angles between k and j′, the distribution of dihedral angle denoting k‐k′‐j′ correlation and the angular distribution of product rotational vectors in the form of polar plots, indicate that the degree of rotational alignment of the product HF molecule is strong and the degree of the rotational alignment decreases as collision energy increases. The product rotational angular momentum vector j′ is not only aligned, but also oriented along the y‐axis, and the molecular rotation of the product prefers an in‐plane reaction mechanism rather than the out‐of‐plane mechanism. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Dynamical studies of the reaction H~++He+He→HHe~++He

The classical trajectory method is employed to calculate the rate coefficient kr for the reaction H++He+He-HHe++He at the temperatures ranging 200-350 K, based on an ab initio potential energy surface. The results show that kr is strongly dependent on the temperature, which can be well fitted by the function kr=ATDr-3 with A=4.192x10-31 cm6/s and the reaction dimension Dr=2.706. The product molecules HHe+ are found in high vibrational states.     

Energy resonance and inverse population of the product vibrational states for the reaction F + H2(v = 0) → HF(v′) + H,LCAC‐SW theoretical quantum scattering study

LCAC‐SW (linear combination of arrangement channel‐scattering wavefunction) method was used to calculate collinear state‐to‐state reaction probabilities for the reaction F + H₂(v = 0) → HF(v′) + H on the 6SEC potential energy surface. The results show that reaction probabilities P₀₂ and P₀₃ [i. e., v′ = 2,3 for reaction F + H₂ (v = 0) + HF(v′) + H] are primary, the population of product vibrational states is inverse and the reaction probabilities are oscillatory with collision energies, i.e., there is energy resonance in this reaction, which agrees with a new experiment.     

A quantum scattering study of collinear state-testate reaction probabilities for the F + H2(v) → HF(v') + H system

A new quantum scattering approach (linear combination of arrangement channels-scattering wavefunction, LCAC-SW) proposed by Deng and his co-workers is used to calculate collinear state-to-state reaction probabilities for the F + H₂(v) → HF(v') + H system. Several interesting problems such M threshold energy, compound states and enhance by translational energy of the reactants and the vibration excitation of products are discussed and they are compared with other theoretical investigations reported in the literature. It is shown that the LCAC-SW approach is the successful one of quantum scattering methods.     

Topological studies on IRC paths of X+H2→XH+H reactions

Topological properties of potential energy and electronic density distribution on five reaction paths X+H₂→XH+H (X=H, N, HN, H₂C, NC) are investigated at the level of UMP2/6–311G(d,p). It has been found that in the region of the reaction paths studied, where B(r_c)|_s>0 [B(r_c)|_s is the product of ρ(r_c) and ∇²ρ(r_c) at the point of reaction process, i.e., B(r_c)|_s=ρ(r_c)∇² ρ(r_c)] is basically the same as the region of V′(s)<0[V′(s) is the second derivative of potential energy with respect to the reaction coordinate, i.e., V′(s)=d²V/ds²], and the point with maximum B(r_c)|_s is almost coincident with the point of minimum V′(s). It can be concluded from the calculated results that there is a good correlation between the topological properties of potential energy and electronic density distribution along the reaction path. The structure transition state of such collinear reactions may be determined by topological analysis of electronic density. © 1997 John Wiley & Sons, Inc. J Comput Chem 18: 1167–1174     

The photoassociation and photodissociation in laser‐assisted collision reaction H + Cl+

The collision reaction H + Cl⁺ assisted by the ultra‐short laser pulse is investigated using the time‐dependent quantum wave packet method. The probability of dissociation depends on the yield ratio of association product HCl⁺. The greater the laser frequency is, the lower the vibrational level of HCl⁺ is. With lowering laser frequency, the probabilities of photoassociation and photodissociation increase, and the ratio of products H⁺ + Cl(²P⁰) to H(²S) + Cl⁺(¹D) also increases. The kinetic energy spectra of the dissociated fragments at low frequency are wider than those at high frequency. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Quantum wave packet calculation of reaction probabilities,cross sections,and rate constants for H+ + LiH → Li + H reaction

The H⁺ + LiH → Li + H reactive scattering has been studied using a quantum real wave packet method. The state‐to‐state and state‐to‐all reaction probabilities for the entitled collision have been calculated at zero total angular momentum. The probabilities for J > 0 are estimated from the J = 0 results by using J‐shifting approximation based on the Capture model. The integral cross sections and thermal rate constants are then calculated. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Temperature dependence of the rate constant of ionic analogs of H + H2 → H2 + H reaction by quantum instanton method

Rate constant ratios k(T)/k(1,500K) for two symmetrical reactions H^? + H₂ → H₂ + H^? and H⁺ + H₂ → H₂ + H⁺ are reported. Direct method based on quantum instanton approximation for evaluation of the temperature dependence of the quantum‐mechanical reaction rate constant is used. Implementation of the theory involves thermodynamic integration and path integral Monte Carlo method. Results of anionic case shows resemblance to neutral case, whereas cationic case is significantly different and below 1,000K rate constant shows strong deviation form linearity of Arrhenius plot due to high activation barrier. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Real wave packet and flux analysis studies of the H + F2 → HF + F reaction

The H + F₂ → HF + F reaction on ground state potential energy surface is investigated using the quantum mechanical real wave packet and Flux analysis method based on centrifugal sudden approximation. The initial state selected reaction probabilities for total angular momentum J = 0 have been calculated by both methods while the probabilities for J > 0 have been calculated by Flux analysis method. The initial state selected reaction probabilities, integral cross sections and rate coefficients have been calculated for a broad range of collision energy. The results show a large rotational enhancement of the reaction probability. Some resonances were seen in the state‐to‐state reaction probabilities while state‐to‐all reaction probabilities and the reaction cross section do not manifest any oscillations and the initial state selected reaction rate constants are sensitive to the temperature. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012     

Nonadiabatic collision model calculations of ion-pair formation cross sections of Cs+O2→Cs++O

This paper has improved Hickman's nonadiabatic collision model by substituting Hickman's constant velocity classical straight line trajectory approximation with the solution of motion equation mR=?dV(R)/dR, and has calculated the cross sections of ion-pair formation Cs+O₂→Cs⁺+O^?₂ with the improved nonadiabatic collision model (INCM). A comparison of our results with other theoretical and experimental results has been made.     

Quantum trajectory method calculation of collinear collisions in the hydrogen exchange reaction H + ClH′ → HCl + H′

The quantum trajectory method was used to study the collinear reaction H + ClH′ → HCl + H′. The potential energy surface was calculated on the QCISD(T)/6-311++G(3df,3pd) level. The reaction probabilities are in good accord with the results obtained by solving the Schroedinger equation using the finite difference method. Translated from Teoreticheskaya i éksperimental’naya Khimiya, Vol. 45, No. 3, pp. 156-159, May-June, 2009.     

An investigation on the reaction mechanism of the F2+Cl2→2ClF using the B3LYP method

The reaction mechanism of F₂+Cl₂→2ClF has been investigated with the density functional theory at the B3LYP/6‐311G* level. Six transition states have been found for the three possible reaction paths and verified by the normal mode vibrational and IRC analyses. Ab initio MP2/6‐311G* geometry optimizations and CCSD(T)/6‐311G(2df)//MP2/6‐311G* single‐point energy calculations have been performed for comparison. It is found that when the F₂ (or Cl₂) molecule decomposes into atoms first and then the F (or Cl) atom reacts with the molecule Cl₂ (or F₂) nearly along the molecular axis, the energy barrier is very low. The calculated energy barrier of F attacking Cl₂ is zero and that of Cl attacking F₂ is only 15.57 kJ?mol^?1 at the B3LYP level. However, the calculated dissociation energies of F₂ and Cl₂ are as high as 145.40 and 192.48 kJ?mol^?1, respectively. When the reaction proceeds through a bimolecular reaction mechanism, two four‐center transition states are obtained and the lower energy barrier is 218.69 kJ?mol^?1. Therefore, the title reaction F₂+Cl₂→2ClF is most probably initiated from the atomization of the F₂ molecule and terminated by the reaction of F attacking Cl₂ nearly along the Cl? Cl bond. MP2 calculations lead to the same conclusion, but the geometry of TS and the energy barrier are somewhat different. © 2002 John Wiley & Sons, Inc. Int J Quantum Chem, 2002     

Rate coefficients of the CF3CHFCF3 + H → CF3CFCF3 + H2 reaction at different temperatures calculated by transition state theory with ab initio and DFT reaction paths

The minimum energy path (MEP) of the reaction, CF₃CHFCF₃ + H → transition state (TS) → CF₃CFCF₃ + H₂, has been computed at different ab initio levels and with density functional theory (DFT) using different functionals. The computed B3LYP/6‐31++G**, BH&HLYP/cc‐pVDZ, BMK/6‐31++G**, M05/6‐31+G**, M05‐2X/6‐31+G**, UMP2/6‐31++G**, PUMP2/6‐31++G**//UMP2/6‐31++G**, RCCSD(T)/aug‐cc‐pVDZ//UMP2/6‐31++G**, RCCSD(T)/aug‐cc‐pVTZ(spd,sp)//UMP2//6‐31++G**, RCCSD(T)/CBS//M05/6‐31+G**, and RCCSD(T)/CBS//UMP2/6‐31++G** MEPs, and associated gradients and Hessians, were used in reaction rate coefficient calculations based on the transition state theory (TST). Reaction rate coefficients were computed between 300 and 1500 K at various levels of TST, which include conventional TST, canonical variational TST (CVT) and improved CVT (ICVT), and with different tunneling corrections, namely, Wigner, zero‐curvature, and small‐curvature (SCT). The computed rate coefficients obtained at different ab initio, DFT and TST levels are compared with experimental values available in the 1000–1200 K temperature range. Based on the rate coefficients computed at the ICVT/SCT level, the highest TST level used in this study, the BH&HLYP functional performs best among all the functionals used, while the RCCSD(T)/CBS//MP2/6‐31++G** level is the best among all the ab initio levels used. Comparing computed reaction rate coefficients obtained at different levels of theory shows that, the computed barrier height has the strongest effect on the computed reaction rate coefficients as expected. Variational effects on the computed rate coefficients are found to be negligibly small. Although tunneling effects are relatively small at high temperatures (～1500 K), SCT corrections are significant at low temperatures (～300 K), and both barrier heights and the magnitudes of the imaginary frequencies affect SCT corrections. © 2012 Wiley Periodicals, Inc.     

Exploring the reaction dynamics of O(3P)+ (X2 )→OH+(X3Σ−)+ H(2S) reaction with time‐dependent wave packet method

The O(³P)+ reaction has been investigated by employing time‐dependent quantum wave packet with split operator method on potential energy surface of the doublet ground‐state H₂O⁺(1²A″). The reaction probabilities and integral cross sections are calculated using centrifugal sudden approximation, which basically agree with the quasi‐classical results of Paniagua et al. [Phys. Chem. Chem. Phys. 2014, 16, 23594]. Moreover, the effect of vibrational and rotational excitation of reactant is investigated. The results show that the vibrational and rotational excitation effects on the integral cross section are not obvious. The little differences between Coriolis coupling results and centrifugal sudden approximation ones show that the cheaper centrifugal sudden calculations here reported are effective for this reaction.     

Z and E rotamers of N‐formyl‐1‐bromo‐4‐hydroxy‐3‐methoxymorphinan‐6‐one and their interconversion as studied by 1H/13C NMR spectroscopy and quantum chemical calculations

N‐Formyl‐1‐bromo‐4‐hydroxy‐3‐methoxymorphinan‐6‐one (compound 2 ), an important intermediate in the NIH Opiate Total Synthesis, presumably exists as a mixture of two rotamers (Z and E) in both CHCl₃ and DMSO at room temperature due to the hindered rotation of its N‐C18 bond in the amide moiety. By comparing the experimental ¹H and ¹³C chemical shifts of a single rotamer and the mixture of compound 2 in CDCl₃ with the calculated chemical shifts of the geometry optimized Z and E rotamers utilizing density functional theory, the crystalline rotamer of compound 2 was characterized as having the E configuration. The energy barrier between the two rotamers was also determined with the temperature dependence of ¹H and ¹³C NMR coalescence experiments, and then compared with that from the reaction path for the interconversion of the two rotamers calculated at the level of B3LYP/6‐31G*. Detailed geometry of the ground state and the transition states of both rotamers are given and discussed. Copyright © 2012 This article is a US Government work and is in the public domain in the USA.     

Solid‐state reaction as a mechanism of 1T ↔ 2H transformation in MoS2 monolayers

Monolayers of molybdenum disulfide MoS₂ are considered to be prospective materials for nanoelectronics and various catalytic processes. Since in certain conditions they undergo 1T ? 2H phase transitions, studying these phase changes is an urgent task. We present a DFT research of these transitions to show that they can proceed as a solid‐state reaction. Two transition states were discovered with energy barriers 1.03 and 1.40 eV. Sulfur atoms in the transition states are shown to be displaced relative to molybdenum atoms so that a tendency of one structural modification to transform into the other modification is seen. This kind of displacements agrees with electron microscopy data reported earlier. The energy parameters indicate that 1T → 2H reactions are exothermic for both transition states and can possibly proceed in a self‐sustained manner when initially activated by some external energy impact. © 2015 Wiley Periodicals, Inc.     

Global X2A′ potential energy surface of Li2H and quantum dynamics of H + Li2 (X1Σg+) → Li + LiH (X1Σ+) reaction

A global potential energy surface (PES) for the electronic ground state of Li₂H system is constructed over a large configuration space. About 30 000 ab initio energy points have been calculated by MRCI‐F12 method with aug‐cc‐pVTZ basis set. The neural network method is applied to fit the PES and the root mean square error of the current PES is only 1.296 meV. The reaction dynamics of the title reaction has been carried out by employing time‐dependent wave packet approach with second order split operator on the new PES. The reaction probability, integral cross section and thermal rate constant are obtained from the dynamics calculation. In most of the collision energy regions, the integral cross sections are in well agreement with the results reported by Gao et al. The rate constant calculated from the new PES increases in the temperature range of present investigation.