请教CP2K对La的基组与赝势定义更改与选择

renzhogn424 · 发表于 Post on 2022-3-31 16:35:41

本帖最后由 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哪个更好？

jiangning198511 · 发表于 Post on 2022-3-31 16:48:08

估计是因为两个文件里面都有相同名字的势函数，你可以考虑修改成不同名称

renzhogn424 · 发表于 Post on 2022-4-1 13:07:53

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

好的，这应该是个办法。

Phantomzyz · 发表于 Post on 2024-10-29 16:56:34

感谢大佬的帖子，间接帮助我解决了问题

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

[CP2K] 请教CP2K对La的基组与赝势定义更改与选择

评分 Rate