计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 量子化学 (Quantum Chemistry) › 求助势能面扫描质子转移曲线中势垒计算

求助势能面扫描质子转移曲线中势垒计算

查看数: 5095 | 评论数: 4 | 收藏 Add to favorites 3

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

zhb

发布时间: 2021-9-2 20:58

正文摘要:

老师您好，第一次接触势能面扫描，想问一下，像这种按键长扫描后的势能曲线，质子转移势垒和逆转移势垒是怎么算出来呢？可能问题有些低级，希望老师见谅，谢谢老师

回复 Reply

sobereva 发表于 Post on 2021-9-3 01:31:48

如果没有特殊情况，不建议做扫描，而应当找质子转移的过渡态，然后跑完整的IRC。

当年TDDFT在Gaussian里不支持解析Hessian，做ESIPT的人不得不靠柔性扫描来说事，而如今G16 TDDFT都有解析Hessian了，对激发态找TS、走IRC和基态一个流程，就再也没理由勉强通过柔性扫描来说事了

wzkchem5 发表于 Post on 2021-9-2 23:27:20

zhb 发表于 2021-9-2 15:28
谢谢老师回复，老师，那这个某点的能量是选例如这个Scan图下方的这个Total Energy (Hartree) 能量值吗？

对

wzkchem5 发表于 Post on 2021-9-2 22:08:28

作者专门把那3个点画大，看图猜也能猜出来了。就是找到最高点，代表过渡态，再在键长更短和更长的地方分别找到一个局部极小值点，代表反应物和产物（哪个是反应物哪个是产物要从分子结构和键长判断），然后用过渡态能量分别减去两个极小值点的能量就行了。
不过需要注意的是，文献这个计算方法是近似的，只有在扫描步长无穷小时才是严格的。当然步长无穷小是做不到的，所以严格的做法需要优化过渡态和极小值点的结构，再减。再者还有一个小问题，除非用的是精度极高的特殊计算方法，否则kcal/mol为单位的能量一般应该保留到小数点后1位，而不是2位

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register