计算化学公社

标题: 求助周期性体系弱相互作用研究基组级别选择 [打印本页]

作者
Author: linlink 时间: 2024-5-19 16:54
标题: 求助周期性体系弱相互作用研究基组级别选择
各位老师好，最近在用Multiwfn生成CP2K输入文件然后对周期性晶胞进行优化，看了sob老师的博文对于基组和理论方法的选择上有两个问题
1.请问老师，研究＞500原子的周期性晶胞的弱相互作用选用怎样的基组和方法比较好，此处优化的时候是否只需对H原子进行优化并且是否只优化结构不优化晶胞（在老师的例子中对H原子优化用的6-31G*，这个是否可用）
2.后续想利用Multiwfn对体系中内部氢键键能，色散，范德华力进行计算，问题1中选用的基组和方法是否还适用
3.从大晶体中挖出簇进行相关分析，是否需要对簇进行优化，优化时选用的基组和方法是否和问题1中一致，是否需要对H原子以外的原子进行优化
谢谢老师

作者
Author: sobereva 时间: 2024-5-20 02:35
1 PBE-D3(BJ)/DZVP-MOLOPT-SR-GTH优化通常够用了。看具体研究哪里的弱相互作用，如果体相内的，不需要优化重原子，除非晶体解析度很差

2 适用

3 同1

作者
Author: linlink 时间: 2024-5-30 08:50

sobereva 发表于 2024-5-20 02:35
1 PBE-D3(BJ)/DZVP-MOLOPT-SR-GTH优化通常够用了。看具体研究哪里的弱相互作用，如果体相内的，不需要优化 ...

老师您好，我后面计算弱相互作用，AIM等挖的团簇并根据您大体系弱相互作用的博文完成了优化和单点计算。但是，我看到您好多博文还有更好的一些方式去计算，虽然对我这个体系有点算不动。不过假如我用算大体系的方法计算并发文章的话，专家会质疑吗？有点纠结，请教一下老师

作者
Author: sobereva 时间: 2024-5-30 08:54

linlink 发表于 2024-5-30 08:50
老师您好，我后面计算弱相互作用，AIM等挖的团簇并根据您大体系弱相互作用的博文完成了优化和单点计算。 ...

看不懂你的意思

作者
Author: linlink 时间: 2024-5-30 09:00

sobereva 发表于 2024-5-30 08:54
看不懂你的意思

不好意思老师，就是我按照大体系弱相互作用的博文计算了弱相互作用和AIM等。但是选择的基组和方法精度不是最高的，只是最适合这个大体系的。想问这样计算发文章，会被审稿人质疑方法基组选取问题吗

作者
Author: reid 时间: 2024-5-30 23:28

linlink 发表于 2024-5-30 09:00
不好意思老师，就是我按照大体系弱相互作用的博文计算了弱相互作用和AIM等。但是选择的基组和方法精度不 ...

结构优化PBE-D3(BJ)/DZVP-MOLOPT-SR-GTH优化通常够用了，如果算能量，再用大一点的基组

作者
Author: sobereva 时间: 2024-5-31 09:34

linlink 发表于 2024-5-30 09:00
不好意思老师，就是我按照大体系弱相互作用的博文计算了弱相互作用和AIM等。但是选择的基组和方法精度不 ...

看具体用什么级别以及体系大小
如果体系不大，用更好的级别也轻易算得动，那就自然可以用更好点的，省得被审稿人嫌弃。比如好几百个原子只能用得动6-31G**，自然没法用更好的。而一百多原子用更好的6-311G**也不费劲。就十几个或者小几十个原子，若追求精度，可以用更好的def2-TZVP，而这时候用6-31G**就显得太low了。

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