论文部分内容阅读
在液相下的化学反应,实验上早有研究。但是在理论上反应机制的探索却是从上世纪初开始的。相比气相,溶剂的存在会带来更多的反应自由度,更多的统计抽样,直接对反应体系进行单纯的量子力学计算,费时费力。另一方面,已有的分子力学方法在计算上的效率是其最大优势,但不能正确描述反应过程中旧键的断裂,以及新键的生成过程,是其最大的缺陷。70年代,Warshel提出,可以尝试把两种方法混合在一起,建立起了量子力学/分子力学结合的方法,使得液相下和大分子体系的研究发展迅速。本文的工作正是采用量子力学与分子力学方法(QM/MM方法),第一次在理论上研究了水溶液中卤代烃CH3X+Yˉ→CH3Y+Xˉ(X=Br,Cl;Y=OH,CN)在水溶液中的亲核取代反应。计算中引入了三种表象,ESP/MM、DFT/MM、和CCSD(T)/MM,从精度较低但效率高的ESP/MM表象开始,利用表象间能量差值的变换,实现高精度的CCSD(T)层面上的理论成果。 首先将气相下的反应物结构放入水中,水分子考虑的是 SPC(E)模型,选用标准的amber力场。在DFT/MM表象下对体系优化,得到液相下的反应物。模拟键的断裂和生成,找到产物并进行优化。利用 NEB(弹性波带法)建立起反应的初始路径,最高点被怀疑为过渡态,通过数值频率计算,鞍点搜索,找到真实的过渡态结构。在最终确定的反应物和产物基础上,运用过渡态搜索理论,分子动力学模拟及NEB方法等得到了DFT/MM和CCSD(T)/MM两种表象下的最小反应路径和沿反应路径的自由能曲线。其中,在DFT/MM和CCSD(T)/MM两种表象下,CH3Br+OHˉ在水溶液中反应的势垒分别是19.2kcal/mol和22.8kcal/mol,CH3Br+CNˉ在水溶液中的势垒分别是16.5kcal/mol和19.1kcal/mol,CH3Cl+CNˉ在水溶液中的势垒分别是18.3kcal/mol和20.3kcal/mol,与已有的实验结果对比,都非常接近。也可以看出DFT理论对活化能垒的计算偏低,本文中得到的CCSD(T)结果更为精确。液相环境对反应势垒的贡献来源于溶剂化效应和极化效应,其中溶剂化效应是主要的。 本文共分为六章,第一章简要介绍了卤代烃反应的研究进展和基本概念,;第二章介绍研究中用到的理论方法和相关理论知识;第三、第四章和第五章分别阐述液相下的亲核取代反应CH3Br+CNˉ、CH3Br+OHˉ和CH3Cl+CNˉ的反应机制;第六章为总结与展望。