计算化学公社

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

作者
Author: 战泽天 时间: 2021-5-26 15:35
标题: 求助cp2k的赝势和基组选择
各位老师好，最近刚学cp2k，有以下几点求助各位老师：1、该如何基于自己的体系选择基组和赝势呢？对于不同基组的精度和介绍可以在哪看呢

2、我用DZVP-MOLOPT-SR-GTH基组和GTH-PBE赝势计算了以下Fe2SiO4的结构优化（包含晶格参数优化），发现计算得到的晶格参数比实验值小很多，随着优化步数的增加，晶格参数越来越小，是我的设置有什么问题吗？麻烦大家帮忙看一下

作者
Author: sobereva 时间: 2021-5-27 05:09
以后发帖时标题里不要带叹号，给你改了

1 一般情况优化用DZVP-MOLOPT-SR-GTH就可以（赝势自然是结合标配的GTH赝势。如果有专门对你用的泛函优化的GTH赝势就用专用的，其它情况用PBE的凑合）。能量可以用更好的TZV2P-MOLOPT-GTH（但不是所有元素都有）。不同的MOLOPT系列基组的精度对比看MOLOPT基组原文J. Chem. Phys., 127, 114105 (2007)。对于结合multiwfn做波函数分析，pob-TZVP比较适合，见比如《使用Multiwfn结合CP2K通过NCI和IGM方法图形化考察固体和表面的弱相互作用》（http://sobereva.com/588）里的说明。

2 要贴就贴带完整结构信息的输入文件。应当考虑k点，单胞跑的话至少2 4 5。当前计算也不算是夸张的越来越小，也不要指望PBE优化得到的晶胞参数能和实验特别相符。还要注意对磁性的考虑。

作者
Author: 战泽天 时间: 2021-5-27 10:33

sobereva 发表于 2021-5-27 05:09
以后发帖时标题里不要带叹号，给你改了

1 一般情况优化用DZVP-MOLOPT-SR-GTH就可以（赝势自然是结合标配 ...

谢谢sob老师的解答，我用的基组和赝势跟您推荐的一样，我试试设置K点计算

作者
Author: qq870624836 时间: 2022-3-8 16:24
本帖最后由 qq870624836 于 2022-3-8 19:07 编辑

sobereva 发表于 2021-5-27 05:09
以后发帖时标题里不要带叹号，给你改了

1 一般情况优化用DZVP-MOLOPT-SR-GTH就可以（赝势自然是结合标配 ...

C:\Users\Administrator\Desktop\QQ图片20220308162033,老师好，def2tzvp这个基组应该是可以计算Pb元素的把为什么在cp2k这里会报错呀。我该怎么解决呢

作者
Author: sobereva 时间: 2022-3-9 04:54

qq870624836 发表于 2022-3-8 16:24
,老师好，def2tzvp这个基组应该是可以计算Pb元素的把为什么在cp2k这里会报错呀。我该怎么解决呢

CP2K里没法用def2-TZVP赝势基组（def2系列对第五周期开始是赝势基组）

作者
Author: hydrogen2333 时间: 6 day ago

sobereva 发表于 2021-5-27 05:09
以后发帖时标题里不要带叹号，给你改了

1 一般情况优化用DZVP-MOLOPT-SR-GTH就可以（赝势自然是结合标配 ...

社长，请问周期性晶体进行变胞优化的时候也可以使用DZVP-MOLOPT-SR-GTH基组和GTH赝势吗，就是Multiwfn默认设置里的组合

作者
Author: sobereva 时间: 5 day ago

hydrogen2333 发表于 2025-8-16 10:53
社长，请问周期性晶体进行变胞优化的时候也可以使用DZVP-MOLOPT-SR-GTH基组和GTH赝势吗，就是Multiwfn默 ...

可以

作者
Author: sobereva 时间: 5 day ago

sobereva 发表于 2022-3-9 04:54
CP2K里没法用def2-TZVP赝势基组（def2系列对第五周期开始是赝势基组）

偶然看到了我自己之前的回答，这里更新一下，如今CP2K已经支持用Gaussian等量子化学程序用的SDD赝势和赝势基组了。较新版本的北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）课程里讲 “能量的计算及相关问题” 部分给了具体使用例子

作者
Author: hydrogen2333 时间: 5 day ago

sobereva 发表于 2025-8-17 00:21
可以

谢谢社长。我的遇到cp2k的SCF不收敛的问题了，我看过http://sobereva.com/665这篇博文。我的体系是Li6PS5I周期性体系，在PBE泛函、DZVP-MOLOPT-SR-GTH基组下进行过cell_opt并收敛，接着进行单点计算，由于TZV2P-MOLOPT-GTH缺少Li元素所以我换成了pob-TZVP-rev2，所有的输入文件都是用Multiwfn生成的，在128个电子步内依然没有收敛且呈现振荡趋势，请问这种情况应该怎么办呢？我想尝试用更小的的基组，收敛后用其波函数做初猜，但是博文里提到

注意两次的基组特征必须定性一致，比如之前计算用的大核赝势，而现在用小核赝势或全电子基组，这时候读之前的波函数就没意义了

请问一下我还可以使用哪些更小的基组呢？或者说还有其他的方法吗

作者
Author: sobereva 时间: 4 day ago

hydrogen2333 发表于 2025-8-17 11:19
谢谢社长。我的遇到cp2k的SCF不收敛的问题了，我看过http://sobereva.com/665这篇博文。我的体系是Li6PS5 ...

TZV2P-MOLOPT-GTH缺少Li就换其它的呗，又不是所有元素用的基组必须相同，诸如若TZVP-MOLOPT-GTH有就用TZVP-MOLOPT-GTH的，没有就再找更低的，DZ或带SR版的。并且不要只盯着一个基组文件，MOLOPT系列基组的文件有好几个，北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）里我讲到了：

(, 下载次数 Times of downloads: 0)

作者
Author: hydrogen2333 时间: 4 day ago

sobereva 发表于 2025-8-18 03:00
TZV2P-MOLOPT-GTH缺少Li就换其它的呗，又不是所有元素用的基组必须相同，诸如若TZVP-MOLOPT-GTH有就用TZV ...

好的谢谢社长，我错误地以为每个元素都要用相同的赝势。现在换成了TZV2P-MOLOPT-PBE-GTH，这个基组元素覆盖比较全面，泛函用的GTH-PBE，算能量时收敛趋势非常明显了

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