计算化学公社

标题: cp2k的赝势和基组选择问题 [打印本页]

作者
Author: 1739140886 时间: 2022-5-6 17:25
标题: cp2k的赝势和基组选择问题
刚开始接触CP2K，我计算的体系主要是有机分子（前两周期元素，原子数大概500左右），Gaussian有些吃力，也在发帖询问后了解到可以使用ORCA，后续会对ORCA继续了解。

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

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

作者
Author: yjb 时间: 2022-5-6 18:03
&KIND 查看元素所有的可选基组/赝势用： grep Si（元素）BASIS_MOLOPT

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

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

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

作者
Author: sobereva 时间: 2022-5-6 18:35
TPSS结合TZVP-MOLOPT-GTH
对于很大的分子，特别是柔性的，偶极矩会对构象有显著的依赖性，注意想清楚计算偶极矩到底有没有实际意义。

作者
Author: 1739140886 时间: 2022-5-6 20:59

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

感谢sob老师！！

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