Quantum dynamics calculation of reaction probability for H+Cl_2→HC1+Cl

We present in this paper a time-dependent quantum wave packet calculation of the initial state selected reaction probability for H + CI2 based on the GHNS potential energy surface with total angular momentum J= 0. The effects of the translational, vibrational and rotational excitation of CI2 on the reaction probability have been investigated. In a broad region of the translational energy, the rotational excitation enhances the reaction probability while the vibrational excitation depresses the reaction probability. The theoretical results agree well with the fact that it is an early down-hill reaction.     

On the separation of motions in the theory of polyatomic molecular spectra

A method has been suggested for the separation of rotational and vibrational motions of nuclei in polyatomic molecules. The method is based on the introduction of the angular velocity vector associated with the kinetic moment and provides the complete separation of motions in the case of arbitrary molecule deformations. The method requires no a priori assumptions on vibrational amplitudes.     

Quantum dynamics calculation of reaction probability for H + Cl2 → HCl + Cl

We present in this paper a time-dependent quantum wave packet calculation of the initial state selected reaction probability for H + Cl₂ based on the GHNS potential energy surface with total angular momentumJ = 0. The effects of the translational, vibrational and rotational excitation of Cl₂ on the reaction probability have been investigated. In a broad region of the translational energy, the rotational excitation enhances the reaction probability while the vibrational excitation depresses the reaction probability. The theoretical results agree well with the fact that it is an early down-hill reaction.     

Quantum dynamics calculation of reaction probability for H + Cl2 → HCl + Cl

We present in this paper a time-dependent quantum wave packet calculation of the initial state selected reaction probability for H + Cl₂ based on the GHNS potential energy surface with total angular momentumJ = 0. The effects of the translational, vibrational and rotational excitation of Cl₂ on the reaction probability have been investigated. In a broad region of the translational energy, the rotational excitation enhances the reaction probability while the vibrational excitation depresses the reaction probability. The theoretical results agree well with the fact that it is an early down-hill reaction.     

Theoretical study of the thermal unimolecular rearrangement of fluoroethylidenes

The thermal unimolecular isomerization of fluoroethylidenes to the corresponding fluoroethylenes has been studied by the MNDO method. It has been shown that fluorine substitution on the carbene carbon increases the activation energy in comparison with the ethylidene rearrangement. To understand the reason for this increase in the activation energy, the charge-transfer effects have been analyzed. Fluorine substitution at other positions does not significantly affect the activation energies. The thermodynamic parameters for the reaction have been evaluated, using vibrational and rotational spectral data calculated in this work. RRKM calculations have been performed and high-pressure Arrhenius parameters calculated. Hydrogen–deuterium kinetic isotope effects indicate that the reaction rates are altered considerably on isotopic substitution, and the change in reaction rates depends upon the position of deuterium substitution, as well as on the number of hydrogens replaced by deuterium atoms. © 1992 John Wiley & Sons, Inc.     

Phase space geometry of dynamics passing through saddle coupled with spatial rotation

Nonlinear reaction dynamics through a rank-one saddle is investigated for many-particle system with spatial rotation. Based on the recently developed theories of the phase space geometry in the saddle region, we present a theoretical framework to incorporate the spatial rotation which is dynamically coupled with the internal vibrational motions through centrifugal and Coriolis interactions. As an illustrative simple example, we apply it to isomerization reaction of HCN with some nonzero total angular momenta. It is found that no-return transition state (TS) and a set of impenetrable reaction boundaries to separate the "past" and "future" of trajectories can be identified analytically under rovibrational couplings. The three components of the angular momentum are found to have distinct effects on the migration of the "anchor" of the TS and the reaction boundaries through rovibrational couplings and anharmonicities in vibrational degrees of freedom. This method provides new insights in understanding the origin of a wide class of reactions with nonzero angular momentum.     

O(3P)+HBr(DBr)反应的含时量子散射计算

基于LEPS势能面, 用三维含时量子波包法对O(³P)+HBr(DBr)反应进行了准确的动力学计算. 计算的结果表明, 振动激发对这个反应是有效的, 而转动激发在某一能量范围内具有方位效应. 计算得到了该反应的速率常数和反应截面, 速率常数kO+HBr的计算值同实验值符合得很好. 通过对相应结果的对比, 可以发现这个反应具有比较明显的同位素效应.     

Statistical evaporation of rotating clusters. IV. Alignment effects in the dissociation of nonspherical clusters

Unimolecular evaporation in rotating, nonspherical atomic clusters is investigated using phase space theory in its orbiting transition state version. The distributions of the total kinetic energy release epsilon(tr) and the rotational angular momentum J(r) are calculated for oblate top and prolate top main products with an arbitrary degree of deformation. The orientation of the angular momentum of the product cluster with respect to the cluster symmetry axis has also been obtained. This statistical approach is tested in the case of the small eight-atom Lennard-Jones cluster, for which comparison with extensive molecular dynamics simulations is presented. The role of the cluster shape has been systematically studied for larger, model clusters in the harmonic approximation for the vibrational densities of states. We find that the type of deformation (prolate versus oblate) plays little role on the distributions and averages of epsilon(tr) and J(r) except at low initial angular momentum. However, alignment effects between the product angular momentum and the symmetry axis are found to be significant, and maximum at some degree of oblateness. The effects of deformation on the rotational cooling and heating effects are also illustrated.     

O(~3P)+HBr(DBr)反应的含时量子散射计算(英文)

基于LEPS势能面, 用三维含时量子波包法对O(3~P)+HBr(DBr)反应进行了准确的动力学计算. 计算的结果表明, 振动激发对这个反应是有效的, 而转动激发在某一能量范围内具有方位效应. 计算得到了该反应的速率常数和反应截面, 速率常数kO+HBr的计算值同实验值符合得很好. 通过对相应结果的对比, 可以发现这个反应具有比较明显的同位素效应.     

A new method for the derivation of exact vibration-rotational kinetic energy operator in internal coordinates

An efficient angular momentum method is presented and used to derive analytic expressions for the vibration-rotational kinetic energy operator of polyatomic molecules.The vibration-rotational kinetic energy operator is expressed in terms of the total angular momentum operator J,the angular momentum operator J and the momentum operator p conjugate to Z in the molecule-fixed frame Not only the method of derivation is simpler than that in the previous work,but also the expressions ot the kinetic energy operators arc more compact.Particularly,the operator is easily applied to different vibrational or rovibrational problems of the polyatomic molecules by variations of matrix elements Gn of a mass-dependent constant symmetric matrix     

用DVR方法研究Ar-HF和Ar-DF的振转光谱

本文用离散变量表示(DVR)方法研究了Ar-HF和Ar-DF体系的振转光谱。对这两个体系的已观测到的振转能级, 用DVR方法得到的计算值与实验结果十分吻合, 误差一般小于0.1cm^-^1, 最大偏差为0.24cm^-^1对Ar-HF的(1113)态。该研究结果与Hutson等用耦合孔道法得到的计算结果基本一致, 但对(3210)振动能级和对(3002)-(3110)态的能级分裂值, 本文计算结果更接近观测值。     

Bimolecular reactions of vibrationally excited molecules. Roaming atom mechanism at low kinetic energies

Quasiclassical trajectory calculations have been performed for the H + H'X(v) → X + HH' abstraction and H + H'X(v) → XH + H' (X = Cl, F) exchange reactions of the vibrationally excited diatomic reactant at a wide collision energy range extending to ultracold temperatures. Vibrational excitation of the reactant increases the abstraction cross sections significantly. If the vibrational excitation is larger than the height of the potential barrier for reaction, the reactive cross sections diverge at very low collision energies, similarly to capture reactions. The divergence is quenched by rotational excitation but returns if the reactant rotates fast. The thermal rate coefficients for vibrationally excited reactants are very large, approach or exceed the gas kinetic limit because of the capture-type divergence at low collision energies. The Arrhenius activation energies assume small negative values at and below room temperature, if the vibrational quantum number is larger than 1 for HCl and larger than 3 for HF. The exchange reaction also exhibits capture-type divergence, but the rate coefficients are larger. Comparisons are presented between classical and quantum mechanical results at low collision energies. At low collision energies the importance of the exchange reaction is enhanced by a roaming atom mechanism, namely, collisions leading to H atom exchange but bypassing the exchange barrier. Such collisions probably have a large role under ultracold conditions. The calculations indicate that for roaming to occur, long-range attractive interaction and small relative kinetic energy in the chemical reaction at the first encounter are necessary, which ensures that the partners can not leave the attractive well. Large orbital angular momentum of the primary products (equivalent to large rotational excitation in a unimolecular reaction) is favorable for roaming.     

Quasi-classical Trajectory Study of Ba+HI→BaI+H Reaction

The quasi-classical trajectory calculations based on extended London-Eyring-Polanyi-Sato potential energy surface have been used to study the reaction of Ba+HI→BaI+H system. The rotational, vibrational, translational, and angular distributions of the product BaI have been calculated. The calculated results are in good agreement with the experimental ones.     

Role of angular momentum conservation in unimolecular translational energy release: validity of the orbiting transition state theory

The translational kinetic energy release distribution (KERD) for the halogen loss reaction of the bromobenzene and iodobenzene cations has been reinvestigated on the microsecond time scale. Two necessary conditions of validity of the orbiting transition state theory (OTST) for the calculation of kinetic energy release distributions (KERDs) have been formulated. One of them examines the central ion-induced dipole potential approximation. As a second criterion, an adiabatic parameter is derived. The lower the released translational energy and the total angular momentum, the larger the reduced mass, the rotational constant of the molecular fragment, and the polarizability of the released atom, the more valid is the OTST. Only the low-energy dissociation of the iodobenzene ion (E approximately 0.45 eV, where E is the internal energy above the reaction threshold) is found to fulfill the criteria of validity of the OTST. The constraints that act on the dissociation dynamics have been studied by the maximum entropy method. Calculations of entropy deficiencies (which measure the deviation from a microcanonical distribution) show that the pair of fragments does not sample the whole of the phase space that is compatible with the mere specification of the internal energy. The major constraint that results from conservation of angular momentum is related to a reduction of the dimensionality of the dynamics of the translational motion to a two-dimensional space. A second and minor constraint that affects the KERD leads to a suppression of small translational releases, i.e., accounts for threshold behavior. At high internal energies, the effects of curvature of the reaction path and of angular momentum conservation are intricately intermeddled and it is not possible to specify the share of each effect.     

Isotopic replacement in ionic systems: the 4He2+ + 3He-->3He 4He+ + 4He reaction

Full quantum dynamics calculations have been carried out for the ionic reaction (4)He(2) (+)+(3)He and state-to-state reactive probabilities have been obtained using both time-dependent and time-independent approaches. An accurate ab initio potential-energy surface has been employed for the present quantum dynamics and the two sets of results are shown to be in agreement with each other. The results for zero total angular momentum suggest a marked presence of atom exchange (isotopic replacement) reaction with probabilities as high as 60%. The reaction probabilities are only weakly dependent on the initial vibrational state of the reactants, while they are slightly more sensitive to the degree of rotational excitation. A brief discussion of the results for selected higher total angular momentum values is also presented, while the l-shifting approximation [S. K. Gray et al., Phys. Chem. Chem. Phys. 1, 1141 (1999)] has been used to provide estimates of the total reaction rates for the title process. Such rates are found to be large enough to possibly become experimentally accessible.     

Vibrational effects on the reaction of NO(2)(+) with C(2)H(2): effects of bending and bending angular momentum

NO(2)(+) in six different vibrational states was reacted with C(2)H(2) over the center-of-mass energy range from 0.03 to 3.3 eV. The reaction, forming NO(+)+C(2)H(2)O and NO+C(2)H(2)O(+), shows a bimodal dependence on collision energy (E(col)). At low E(col), the reaction is quite inefficient (<2%) despite this being a barrierless, exoergic reaction, and is strongly inhibited by E(col). For E(col)> approximately 0.5 eV, a second mechanism turns on, with an efficiency reaching approximately 27% for E(col)>3 eV. The two reaction channels have nearly identical dependence on E(col) and NO(2)(+) vibrational state, and identical recoil dynamics, leading to the conclusion that they represent a single reaction path throughout most of the collision. All modes of NO(2)(+) vibrational excitation enhance both channels at all E(col), however, the effects of bend (010) and bend overtone (02(0)0) excitation are particularly strong (factor of 4). In contrast, the asymmetric stretch (001), which intuition suggests should be coupled to the reaction coordinate, leads to only a factor of approximately 2 enhancement, as does the symmetric stretch (100). Perhaps the most surprising effect is that of the bending angular momentum, which strongly suppress reaction, even though both the energy and angular momentum involved are tiny compared to the collision energy and angular momentum. The results are interpreted in light of ab initio and Rice-Ramsperger-Kassel-Marcus calculations.     

State-to-state quantum dynamics calculations of the C + OH reaction on the second excited potential energy surface

Accurate three-dimensional quantum-mechanical scattering calculations using a time-indepedent hyperspherical method have been performed for the C((3)P) + OH(X(2)Π) → CO(a(3)Π) + H((2)S) reaction on the second excited potential energy surface of 1(4)A″ symmetry. State-to-state reaction probabilities at a total angular momentum J = 0 have been computed in a wide range of collision energies. Many pronounced resonances have been found, espcially at low energy. The product vibrational distributions are noninverted. The present results therefore suggest that the title reaction proceeds via a long-lived intermediate complex. An approximate quantum-mechanical rate constant has also been calculated, and large differences are observed with the quasi-classical trajectory prediction.     

辛准经典轨迹法研究Cl+H2反应

用辛准经典轨迹法模拟了Cl+H2反应在mBW2势能面上的动力学行为。研究了各种初始条件下的反应碰撞截面,产物的能量分配,角度分布和态分布。另外,我们还比较了反应物的三种能量形式(平动能,转动能和振动能)对反庆的有效性。     

The dynamics of the Cl+n-C4H10-->HCl (v',j') + C4H9 reaction at 0.32 eV

Rotational state resolved center-of-mass angular scattering and kinetic energy release distributions have been determined for the HCl (v' = 0, j' = 0-6) products of the reaction of chlorine with n-butane using the photon-initiated reaction technique, coupled with velocity-map ion imaging. The angular and kinetic energy release distributions derived from the ion images are very similar to those obtained previously for the Cl plus ethane reaction. The angular distributions are found to shift from forward scattering to more isotropic scattering with increasing HCl rotational excitation. The kinetic energy release distributions indicate that around 30% of the available energy is channeled into internal excitation of the butyl radical products. The data analysis also suggests that H-atom abstraction takes place from both primary and secondary carbon atom sites, with the primary site producing rotationally cold, forward scattered HCl (v' = 0) products, and the secondary site yielding more isotropically scattered HCl (v' = 0) possessing higher rotational excitation. The mechanisms leading to these two product channels are discussed in the light of the present findings, and in comparison with studies of other Cl plus alkane reactions.     

Quantum chemistry aspects of the solvent effects on the ene reaction of 1‐Phenyl‐1,3,4‐triazolin‐2,5‐dione and 2‐methyl‐2‐butene

A density functional theory study was used to investigate the quantum aspects of the solvent effects on the kinetic and mechanism of the ene reaction of 1‐phenyl‐1,3,4‐triazolin‐2,5‐dione and 2‐methyl‐2‐butene. Using the B3LYP/6–311++ G(d,p) level of the theory, reaction rates have been calculated in the various solvents and good agreement with the experimental data has been obtained. Natural bond orbital analysis has been applied to calculate the stabilization energy of N18? H19 bond during the reaction. Topological analysis of quantum theory of atom in molecule (QTAIM) studies for the electron charge density in the bond critical point (BCP) of N18? H19 bond of the transition states (TSs) in different solvents shows a linear correlation with the interaction energy. It is also seen form the QTAIM analysis that increase in the electron density in the BCP of N18? H19, raises the corresponding vibrational frequency. Average calculated ratio of 0.37 for kinetic energy density to local potential energy density at the BCPs as functions of N18? H19 bond length in different media confirmed covalent nature of this bond. Using the concepts of the global electrophilicity index, chemical hardness and electronic chemical potentials, some correlations with the rate constants and interaction energy have been established. Mechanism and kinetic studies on 1‐phenyl‐1,3,4‐triazolin‐2,5‐dione and 2‐methyl‐2‐butene ene reaction suggests that the reaction rate will boost with interaction energy enhancement. Interaction energy of the TS depends on the solvent nature and is directly related to electron density of the bonds involved in the reaction proceeding, global electrophilicity index and electronic chemical potential. However, the chemical hardness relationship is reversed. Finally, an interesting and direct correlation between the imaginary vibrational frequency of the N18? H19 critical bond and its electron density at the TS has been obtained. © 2014 Wiley Periodicals, Inc.