计算化学公社

标题: 为什么Gaussian中IRC产物与结构优化最低点构型没有完全一致？ [打印本页]

作者
Author: yaoyuan0711 时间: 2021-11-10 14:53
标题: 为什么Gaussian中IRC产物与结构优化最低点构型没有完全一致？
本帖最后由 yaoyuan0711 于 2021-11-10 14:56 编辑

各位老师好！我用Gaussian(M062X/def2TZVP int=ultrafine)算出两个极小点结构： (, 下载次数 Times of downloads: 23) 构型1（两个甲基在异侧）； (, 下载次数 Times of downloads: 45) 构型2（两个甲基在同侧）。
然后计算构型转变的能垒，即甲基和环丙基旋转的能垒。用M062X/def2SVP opt(TS,calcfc,noeigen) int=ultrafine找到了一个“过渡态”。再用M062X/def2SVP IRC=(calcfc,maxpoint=50) int=ultrafine得到IRC曲线 (, 下载次数 Times of downloads: 24) 。
但是IRC的正向和反向产物与两个极小点只是相似，并没有完全一致（甲基在同侧和异侧，但其他原子结构相差较大,如C=N-N=C未像极小点结构一样在一个平面内）。这种情况是因为计算基组不同导致的？能不能说明该“过渡态”基本正确，或者该如何处理呢？ (, 下载次数 Times of downloads: 24) (, 下载次数 Times of downloads: 11)

作者
Author: hebrewsnabla 时间: 2021-11-10 15:38
本来就只应该相似，不是完全一致。IRC收敛的判据和结构优化不一样。

作者
Author: yaoyuan0711 时间: 2021-11-10 17:56
本帖最后由 yaoyuan0711 于 2021-11-10 17:58 编辑

hebrewsnabla 发表于 2021-11-10 15:38
本来就只应该相似，不是完全一致。IRC收敛的判据和结构优化不一样。

那这个IRC结果可以认为是过渡态？我把maxpoint调到200，结构更加相似，但C=N-N=C仍然不在一个平面

作者
Author: hebrewsnabla 时间: 2021-11-10 18:10

yaoyuan0711 发表于 2021-11-10 17:56
那这个IRC结果可以认为是过渡态？我把maxpoint调到200，结构更加相似，但C=N-N=C仍然不在一个平面

可以吧。
或者只要你拿这个irc末端结构去优化，能很快收敛就行

作者
Author: yaoyuan0711 时间: 2021-11-11 10:14

hebrewsnabla 发表于 2021-11-10 18:10
可以吧。
或者只要你拿这个irc末端结构去优化，能很快收敛就行

好的，谢谢您！我再去试试

欢迎光临计算化学公社 (http://bbs.keinsci.com/)