计算化学公社

标题: 请教CP2K对La的基组与赝势定义更改与选择 [打印本页]

作者
Author: renzhogn424 时间: 2022-3-31 16:35
标题: 请教CP2K对La的基组与赝势定义更改与选择
本帖最后由 renzhogn424 于 2022-4-1 13:11 编辑

最近学习CP2K，并用其计算一下含镧的体系，试了6-31G和DZVP等，发现都没有定义，然后在SOB老师的指导下，使用“BASIS_MOLOPT_LnPP1”。于是我手动修改 &KIND La部分为：
&KIND La
   ELEMENT La
   BASIS_SET DZV-MOLOPT-GTH
   POTENTIAL GTH-PBE

再往下看到&DFT部分，原本是
&DFT
BASIS_SET_FILE_NAME  BASIS_MOLOPT
POTENTIAL_FILE_NAME  POTENTIAL

我中间插了一行：“ BASIS_SET_FILE_NAME  BASIS_MOLOPT_LnPP1”，可以进行计算了，但是我查了MOLOPT_LnPP1的原文，发现其有自己配套的赝势，它与默认的赝势“POTENTIAL”中La的定义有少许差别。我又发现在cp2k目录里有“LnPP1_POTENTIALS”，并且与文献中一样，于是加在输出文件后面：
&DFT
BASIS_SET_FILE_NAME  BASIS_MOLOPT
BASIS_SET_FILE_NAME  BASIS_MOLOPT_LnPP1
POTENTIAL_FILE_NAME  POTENTIAL
POTENTIAL_FILE_NAME  LnPP1_POTENTIALS
结果报错，说是赝势重复了，这里就想问一下，如果我想用文献里的赝势怎么办？
可否把文件“POTENTIAL”中La的赝势值替换为“LnPP1_POTENTIALS”的？或者其它写法？？？

后来又发现，对La是有基组的，在文件“BASIS_MOLOPT_UCL”里面，还是DZVP-MOLOPT-SR-GTH的，加载这个文件就能直接计算了。
那么问题来了，DZVP-MOLOPT-SR-GTH和之前提到的DZV-MOLOPT-GTH哪个更好？

作者
Author: jiangning198511 时间: 2022-3-31 16:48
估计是因为两个文件里面都有相同名字的势函数，你可以考虑修改成不同名称

作者
Author: renzhogn424 时间: 2022-4-1 13:07

jiangning198511 发表于 2022-3-31 16:48
估计是因为两个文件里面都有相同名字的势函数，你可以考虑修改成不同名称

好的，这应该是个办法。

作者
Author: Phantomzyz 时间: 2024-10-29 16:56
感谢大佬的帖子，间接帮助我解决了问题

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