这里会显示出您选择的修订版和当前版本之间的差别。
后一修订版 | 前一修订版 | ||
adf:bsse2020 [2020/11/21 12:03] – 创建 liu.jun | adf:bsse2020 [2022/05/12 14:41] (当前版本) – [如何消除BSSE] liu.jun | ||
---|---|---|---|
行 1: | 行 1: | ||
- | ====== 如何计算BSSE ====== | + | ====== |
=====BSSE的含义===== | =====BSSE的含义===== | ||
简单的说,考虑BSSE的原因,就是计算结合能的时候,需要统一基函数,这样结合能定义上才更严格。 | 简单的说,考虑BSSE的原因,就是计算结合能的时候,需要统一基函数,这样结合能定义上才更严格。 | ||
行 15: | 行 15: | ||
我们用H2和苯分子之间的BSSE来举例: | 我们用H2和苯分子之间的BSSE来举例: | ||
- | 1,首先需要对苯+$H_2$的复合体系进行结构优化,几何结构优化的使用可以参见[[adf: | + | 1,首先需要对苯+H< |
2,基于优化好的复合体系的分子结构,先计算该符合体系的单点能(*logfile尾部的bond energy或*.out中的total bonding energy),参数设置参见[[adf: | 2,基于优化好的复合体系的分子结构,先计算该符合体系的单点能(*logfile尾部的bond energy或*.out中的total bonding energy),参数设置参见[[adf: | ||
- | 3,如下图所示,计算$H_2$单体的能量(需要将苯设置为Ghost): | + | 3,如下图所示,计算H< |
{{ adf: | {{ adf: | ||
- | 然后参数设置(如基组、泛函等等)与2一致,计算单点能。这样得到的单点能,就是$H_2$的单点能。因为基函数包含了$H_2$的基函数,也包含苯的基函数,但没有苯的电子、原子核。 | + | 然后参数设置(如基组、泛函等等)与2一致,计算单点能。这样得到的单点能,就是H< |
4,类似地,计算苯单体的单点能能量(需要将H2设置为Ghost) | 4,类似地,计算苯单体的单点能能量(需要将H2设置为Ghost) | ||
行 29: | 行 29: | ||
{{ adf: | {{ adf: | ||
- | 然后参数设置(如基组、泛函等等)与2一致,计算单点能。这样得到的单点能,就是苯的单点能。因为基函数包含了苯的基函数,也包含$H_2$的基函数,但没有$H_2$的电子、原子核。 | + | 然后参数设置(如基组、泛函等等)与2一致,计算单点能。这样得到的单点能,就是苯的单点能。因为基函数包含了苯的基函数,也包含H< |
- | 5,如此得到的复合体系的总能量(第2步的单点能)减去单体的能量(即3, | + | 5,如此得到的复合体系的总能量(第2步的单点能)减去单体的能量(即3, |
<color lightgrey> | <color lightgrey> | ||
- | AMS软件提供**免费试用**(一般为一个月),试用申请方式参考:[[adf: | + | AMS软件提供**免费试用**,试用申请方式参考:[[adf: |