计算化学公社

标题: 关于BSSE+ZPE的双重矫正问题 [打印本页]

作者
Author: wangxuran 时间: 2024-5-6 17:32
标题: 关于BSSE+ZPE的双重矫正问题
本帖最后由 wangxuran 于 2024-5-7 21:42 编辑

各位老师好，学生最近读了Ng@B40的文献，在文中作者计算了BSSE和ZPE校正后的解离能。
一、想问老师是否可以以下几步操作是否正确？
①先优化（opt+freq）整体Ng@B40结构，后利用counterpoise=2计算的解离能Ediss1=E_AB-E_A-E_B，读取了complexation energy后的数值作为BSSE校正后的解离能Ediss1。
②ZPE校正：利用在整体结构中的A,B分结构，只加freq得到A,B的零点矫正能，不进行单独优化
③Ediss总=Ediss1+整体结构的ZPE-A结构的ZPE-B结构的ZPE,得到双重矫正下的解离能

二、学生计算后还发现一个问题，在ZPE校正中，从整体优化后的结构抽出的未opt的A,B结构存在虚频，与单独经过opt的自由A,B单体来比较存在一定的变形能，请问这个矫正是否可以把相应的变形能加进最终的计算中？

附件中为相应的文献。恳请各位老师指教

作者
Author: wzkchem5 时间: 2024-5-7 23:22
这个问题不属于波函数分析，应该发在量子化学版面
ZPE校正只有在平衡结构才有意义，所以只有优化了A和B，才能算ZPE。BSSE校正里的A和B不优化，得到BSSE校正（注意我指的是校正值本身，不是校正后的complexation energy），然后优化A和B，算ZPE校正的解离能，再用你已经得到的BSSE校正，去校正解离能

作者
Author: wangxuran 时间: 2024-5-8 15:04

wzkchem5 发表于 2024-5-7 23:22
这个问题不属于波函数分析，应该发在量子化学版面
ZPE校正只有在平衡结构才有意义，所以只有优化了A和B， ...

明白啦，非常感谢老师！

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