计算化学公社

标题: 求助势能面扫描质子转移曲线中势垒计算 [打印本页]

作者
Author: zhb 时间: 2021-9-2 20:58
标题: 求助势能面扫描质子转移曲线中势垒计算
老师您好，第一次接触势能面扫描，想问一下，像这种按键长扫描后的势能曲线，质子转移势垒和逆转移势垒是怎么算出来呢？可能问题有些低级，希望老师见谅，谢谢老师

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

作者
Author: zhb 时间: 2021-9-2 22:28

wzkchem5 发表于 2021-9-2 22:08
作者专门把那3个点画大，看图猜也能猜出来了。就是找到最高点，代表过渡态，再在键长更短和更长的地方分别 ...

谢谢老师回复，老师，那这个某点的能量是选例如这个Scan图下方的这个Total Energy (Hartree) 能量值吗？

作者
Author: wzkchem5 时间: 2021-9-2 23:27

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

对

作者
Author: sobereva 时间: 2021-9-3 01:31
如果没有特殊情况，不建议做扫描，而应当找质子转移的过渡态，然后跑完整的IRC。

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

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