计算化学公社

标题: 石墨烯负载金属原子体系使用TZVP级别的基组算结构优化后的能量，可以忽略BSSE设置吗 [打印本页]

作者
Author: 1764360701 时间: 2024-3-23 10:25
标题: 石墨烯负载金属原子体系使用TZVP级别的基组算结构优化后的能量，可以忽略BSSE设置吗
本帖最后由 1764360701 于 2024-3-27 09:17 编辑

最近在做一个负载两个金属原子的slab，用来吸附一个N2分子。CP2K软件包算的，打算用几何结构优化后的能量来计算这个slab对N2分子的吸附能，体系共有五十个原子。由于计算经验不足，导致我目前计算出现了一些疑问，希望老师可以给一些建议。
输入文件都是用的Multiwfn产生的，然后自己稍加修改，基组是用的TZVP-MOLOPT-GTH的，这种体系周期性我设置的是XY，真空层设置的是15 ai。最近在考虑结构优化出来的能量的准确性问题，结构优化能跑，用的对角化算法配合SMEAR k点取的 2，2，1。并且收敛效果理想。
1.请问用这样的基组算几何结构优化后的能量还要考虑BSSE的误差吗？
2.只设置XY方向的周期性是否合理？
3.关于K点的选取，X和Y方向的晶胞边长分别是10 ai和9 ai，K点选取2，2，1是否合理？Z方向应该不需要考虑设置k点吧。
真的很感谢能帮助到我的每一个老师，谢谢

作者
Author: sobereva 时间: 2024-3-23 10:39
基组名必须完整写全，否则TZVP谁知道你是指pob-TZVP还是TZVP-GTH还是TZVP-MOLOPT-GTH还是Ahlrichs-def2-TZVP

作者
Author: 1764360701 时间: 2024-3-23 10:55

sobereva 发表于 2024-3-23 10:39
基组名必须完整写全，否则TZVP谁知道你是指pob-TZVP还是TZVP-GTH还是TZVP-MOLOPT-GTH还是Ahlrichs-def2-TZV ...

老师，已修改

作者
Author: sobereva 时间: 2024-3-23 12:19
1 不用。但当前体系不大，可以用更好的TZV2P-MOLOPT-GTH
2 合理
3 埃不要写成ai，看着难受死了
3,3,1更好

作者
Author: 1764360701 时间: 2024-3-23 17:31

sobereva 发表于 2024-3-23 12:19
1 不用。但当前体系不大，可以用更好的TZV2P-MOLOPT-GTH
2 合理
3 埃不要写成ai，看着难受死了

谢谢老师的指导

作者
Author: 1764360701 时间: 2024-3-23 18:00

sobereva 发表于 2024-3-23 12:19
1 不用。但当前体系不大，可以用更好的TZV2P-MOLOPT-GTH
2 合理
3 埃不要写成ai，看着难受死了

老师您好，在.inp文件里Ti这个元素报错没有找到TZV2P-MOLOPT-GTH这个基组，然后我发现Ti元素有TZV2P-MOLOPT-PBE-GTH这个基组，然后我用TZV2P-MOLOPT-PBE-GTH-q12，也就是混合基组，造成的结果误差影响大吗？谢谢老师。

作者
Author: sobereva 时间: 2024-3-23 21:31

1764360701 发表于 2024-3-23 18:00
老师您好，在.inp文件里Ti这个元素报错没有找到TZV2P-MOLOPT-GTH这个基组，然后我发现Ti元素有TZV2P-MOLO ...

MOLOPT基组尽量都用同一个基组文件里的，不要诸如原始的MOLOPT、UZH的MOLOPT、UCL的MOLOPT基组文件里定义的混用

作者
Author: 1764360701 时间: 2024-3-23 22:59

sobereva 发表于 2024-3-23 21:31
MOLOPT基组尽量都用同一个基组文件里的，不要诸如原始的MOLOPT、UZH的MOLOPT、UCL的MOLOPT基组文件里定义 ...

就是比如说我某一个元素用的BASIS_MOLOPT_UZH文件夹里的基组，那么其他元素也尽量要用这个文件夹里的基组，是这个意思吗老师

作者
Author: sobereva 时间: 2024-3-24 11:05

1764360701 发表于 2024-3-23 22:59
就是比如说我某一个元素用的BASIS_MOLOPT_UZH文件夹里的基组，那么其他元素也尽量要用这个文件夹里的基组 ...

虽然没有绝对要求，但最好如此，否则显得莫名其妙

作者
Author: 1764360701 时间: 2024-3-25 09:32

sobereva 发表于 2024-3-24 11:05
虽然没有绝对要求，但最好如此，否则显得莫名其妙

好的，感谢老师指导

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