|
|
本帖最后由 ysc19951114 于 2023-4-25 22:17 编辑
本人目前在学习Gaussian程序的使用,并对实验中苯酚在对位被氯正离子取代氢这个亲电取代反应进行了过渡态计算。根据本人查到的相关资料,该反应分为两步,首先是氯正离子进攻苯环,形成σ络合物,然后σ络合物再脱去一个氢离子生成一个氯苯酚和一个氢离子。因此我将该反应分成了两歩进行计算。其中第一步很快就计算出了看起来很正常的结果(如压缩包step1内所示)。但是第二步脱去氢离子的过程我试了几个星期修改初猜,都只能找到如下两种过渡态构型(a,b)(文件见附件step2)。其中a是氢离子往氯原子方向移动,b是氢离子从对位C往间位C上移动。因此我选择a作为更加可能的反应方向,在a的构型上将H-Cl的距离适当拉大作为新的初猜进行了第三次过渡态搜索,这次的结果为c,氢离子明显远离了苯环,虚频震动也是氢离子的震动,因此我将该过渡态进行了IRC运算(文件见附件step3).IRC端点能优化至我预想的结构,但IRC曲线有点奇怪,显示反应过程没有能垒,且氯苯酚与氢离子的能量高于σ络合物。这是为何?
|
-
a.jpg
(99.74 KB, 下载次数 Times of downloads: 14)
-
b.jpg
(109.41 KB, 下载次数 Times of downloads: 12)
-
c.jpg
(78.26 KB, 下载次数 Times of downloads: 12)
-
-
step1.zip
396 KB, 下载次数 Times of downloads: 12
-
-
step2.zip
353.75 KB, 下载次数 Times of downloads: 5
-
-
step3.zip
211.59 KB, 下载次数 Times of downloads: 5
|