Li Qian Shu, Zhang Yue, Zhang Shaowen
School of Science, Beijing Institute of Technology, 100081 Beijing, People's Republic of China.
J Chem Phys. 2004 Nov 15;121(19):9474-80. doi: 10.1063/1.1807391.
We present a direct ab initio dynamics study of thermal rate constants of the hydrogen abstraction reaction of CH(3)O+H-->CH(2)O+H(2). The unrestricted Becke's half-and-half hybrid functional using the Lee-Yang-Parr correlation functional with Dunning's correlation consistent polarized valence double-zeta basis set, the unrestricted quadratic configuration interaction calculation including single and double substitutions with Dunning's correlation consistent polarized valence double-zeta basis set, and the unrestricted quadratic configuration interaction calculation including single and double substitutions with a triples contribution with Dunning's correlation consistent polarized valence triple-zeta basis set methods were employed to optimize the structures and to calculate frequencies for all stationary points. Minimum energy paths were obtained by the unrestricted Becke's half-and-half hybrid functional using the Lee-Yang-Parr correlation functional and the unrestricted quadratic configuration interaction calculation including single and double substitutions with the same Dunning's correlation consistent polarized valence double-zeta basis set levels of theory. No barrier is found at the unrestricted Becke's half-and-half hybrid functional using the Lee-Yang-Parr correlation functional with Dunning's correlation consistent polarized valence double-zeta basis set level of theory in contrast to a small barrier of 1.43 kcal mol(-1) at the unrestricted quadratic configuration interaction calculation including single and double substitutions with Dunning's correlation consistent polarized valence double-zeta basis set level of theory. In particular, the barrier vanishes as the energies along the minimum energy path MEP are refined at the unrestricted quadratic configuration interaction calculation including single and double substitutions with a triples contribution with Dunning's correlation consistent polarized valence triple-zeta basis set level of theory. Smaller barriers of 0.47 and 0.17 kcal mol(-1) were obtained at the unrestricted quadratic configuration interaction calculation including single and double substitutions with a triples contribution with Dunning's correlation consistent polarized valence triple-zeta basis set and the unrestricted quadratic configuration interaction calculation including single and double substitutions with a triples contribution with Dunning's correlation consistent polarized valence triple-zeta basis set based on the geometries at the unrestricted quadratic configuration interaction calculation including single and double substitutions with Dunning's correlation consistent polarized valence triple-zeta basis set levels of theory, respectively. The forward rate constants are evaluated with the canonical variational transition state theory in the temperature range of 300-2500 K. The calculated forward rate constants at the unrestricted quadratic configuration interaction calculation including single and double substitutions with a triples contribution with Dunning's correlation consistent polarized valence triple-zeta basis set based on the geometries at the unrestricted quadratic configuration interaction calculation including single and double substitutions with Dunning's correlation consistent polarized valence double-zeta basis set level of theory are in good agreement with the available experimental data. The kinetic isotope effects are estimated.
我们对CH(3)O + H→CH(2)O + H(2)氢提取反应的热速率常数进行了直接的从头算动力学研究。使用Lee-Yang-Parr相关泛函和Dunning相关一致极化价双ζ基组的无限制Becke半混合泛函、包括单双取代的无限制二次组态相互作用计算(使用Dunning相关一致极化价双ζ基组)以及包括单双取代且带有三重贡献的无限制二次组态相互作用计算(使用Dunning相关一致极化价三重ζ基组)方法,用于优化所有驻点的结构并计算频率。通过使用Lee-Yang-Parr相关泛函的无限制Becke半混合泛函以及包括单双取代的无限制二次组态相互作用计算(使用相同的Dunning相关一致极化价双ζ基组理论水平)获得最小能量路径。在使用Lee-Yang-Parr相关泛函和Dunning相关一致极化价双ζ基组理论水平的无限制Becke半混合泛函中未发现势垒,而在包括单双取代的无限制二次组态相互作用计算(使用Dunning相关一致极化价双ζ基组理论水平)中存在1.43 kcal mol(-1)的小势垒。特别是,在包括单双取代且带有三重贡献的无限制二次组态相互作用计算(使用Dunning相关一致极化价三重ζ基组理论水平)中,随着沿最小能量路径MEP的能量被细化,势垒消失。在包括单双取代且带有三重贡献的无限制二次组态相互作用计算(使用Dunning相关一致极化价三重ζ基组)以及基于包括单双取代的无限制二次组态相互作用计算(使用Dunning相关一致极化价三重ζ基组理论水平)的几何结构的包括单双取代且带有三重贡献的无限制二次组态相互作用计算(使用Dunning相关一致极化价三重ζ基组)中,分别获得了0.47和0.17 kcal mol(-1)的较小势垒。在300 - 2500 K的温度范围内,使用正则变分过渡态理论评估正向速率常数。在基于包括单双取代的无限制二次组态相互作用计算(使用Dunning相关一致极化价双ζ基组理论水平)的几何结构的包括单双取代且带有三重贡献的无限制二次组态相互作用计算(使用Dunning相关一致极化价三重ζ基组)中计算得到的正向速率常数与现有实验数据吻合良好。还估计了动力学同位素效应。