A quantum mechanical calculation was performed to study the hydrogen atom transfer of collinear OH…O/OD…O system, for which Delves' coordinates and R-matrix propagation method were applied in a Melius-Blint potential energy surface. The calculation result showed that the state-state H atom transfer probability comported strong oscillation phenomena and collision delay time of the title system was in the fs-ps time scale. The kinetic isotope effect was calculated in this work too.
According to the calculation results of the intrapair and interpair correlation energy for the title systems, it has been found that the intrapair correlation energy of K shell of Cl is almost a constant and both the intrashell and intershell correlation energy of K and L shell changes little. It has also been found that in MCI series compounds the value of Cl correlation energy contribution depends on the ionicity of MCI compounds, i.e., the Cl correlation energy contribution increases with the increase of the ionic bond strength of the compound and this value is always less than the correlation energy of Cl" anion but always larger than that of Cl atom. These rules are helpful for the estimation of the correlation energy of ionic compounds and the energy changes of chemical reactions.
The pair correlation energy of bonding electrons is used and analyzed in the cal- culation of CH and CY (Y = F, O, N) bonding electron pairs in CH3X (X = F, OH, NH2) isoelec- tronic systems based on intra- and interpair correlation energy results at both MP2-OPT2/6- 311++G(d) and MP2-OPT2/cc-pVtz levels with MELD program. Comparison of two set results shows that cc-pVtz and 6-311++G(d) give more correlation energy of valence electrons and innermost core electron pairs, respectively in these systems, resulting that the total correlation energy with cc-pVtz basis of each system is larger than that with 6-311++G(d). Investigations of pair correlation energy show that with the decrease of electronegativity of X atom and the increase of H atoms in these CH3X (X = F, OH, NH2) systems, the pair correlation energy of 1sC2 of the C atoms is transferable, and the correlation energy of CH bonding electron pair with little changes is of approximate transferability, while those of CY (CF, CO, CN) bonding electron pair decrease in a large extent from CH3F through CH3OH to CH3NH2 molecules. It is suggested that the study of pair correlation energy of bonding electrons will further deepen the understanding of electron corre- lation effect from traditional chemical bonding concept.
The basic theory of symplectic algorithm was introduced. A comparison between Runge-Kutta method and symplectic integration method was preformed in the simulation of the long time behavior of H + H2 system on BKMP potential energy surface. Our results reveal a dis-sipative behavior in the integral of ordinary differential equation by the fourth order Runge-Kutta method, which causes incorrect simulation results in QCT calculations. However, when the symplectic integration method is applied, the dissipative behavior is not found in the same system. When the initial state is the same, the energy deviation of fourth order symplectic integral method is almost one percent of that of fourth order Runge-Kutta method in a 60000-step simulation, and that of sixth order symplectic integral method is much less. These results show that the symplectic integral methods are always the better choice in the integral calculation of the long time behavior in maintaining energy conservation.
