cp2k的赝势和基组选择问题

1739140886 · 发表于 Post on 2022-5-6 17:25:31

刚开始接触CP2K，我计算的体系主要是有机分子（前两周期元素，原子数大概500左右），Gaussian有些吃力，也在发帖询问后了解到可以使用ORCA，后续会对ORCA继续了解。

目前CP2K基本的输入可以通过Multiwfn解决（感谢sob老师），但是想了解如何选用合适的基组/赝势/方法（对弥散有要求，需要得到偶极矩信息），谢谢大家的回复！

ps：我也尝试了用B3LYP+BSE下来的def2-SVPD基组在CP2K算正丁烷的偶极矩，发现结果似乎不合理（偶极矩得到16Debye）

yjb · 发表于 Post on 2022-5-6 18:03:38

&KIND 查看元素所有的可选基组/赝势用： grep Si（元素）BASIS_MOLOPT

对一些重元素用：
grep La（重元素） BASIS_MOLOPT_UCL

SZV，DZVP, TZVP...代表基组的劈裂和极化函数情况。对价电子做劈裂，用两个收缩的高斯基函数表示，用的高斯基函数越多，计算越准确，但是计算量越大。对于大体系常用DZVP，如果对精度有高要求用TZVP。

cp2k会有BSSE问题，对计算吸附能、给合能等会有偏负。

sobereva · 发表于 Post on 2022-5-6 18:35:39

TPSS结合TZVP-MOLOPT-GTH
对于很大的分子，特别是柔性的，偶极矩会对构象有显著的依赖性，注意想清楚计算偶极矩到底有没有实际意义。

1739140886 · 发表于 Post on 2022-5-6 20:59:49

sobereva 发表于 2022-5-6 18:35
TPSS结合TZVP-MOLOPT-GTH
对于很大的分子，特别是柔性的，偶极矩会对构象有显著的依赖性，注意想清楚计算 ...

感谢sob老师！！

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] cp2k的赝势和基组选择问题