结构中有La,Cu,O,我在使用cp2k的时候，系统显示BASIS_MOLOPT没有La

liuxp · 发表于 Post on 2024-9-4 10:23:55

The requested basis set <DZVP-MOLOPT-SR-GTH-Q11> for element <LA> *
* | was not found in the basis set files <BASIS_MOLOPT>
我这个结构中有La,Cu,O,我在使用cp2k的时候，系统显示BASIS_MOLOPT，没有La，请问一下对于镧系的元素应该使用哪个基组以及使用的精度设置多大，我这个体系大概在250原子左右，后续我还会加上小分子进行振动分析计算反应的自由能。

Dempey · 发表于 Post on 2024-9-4 10:45:01

使用的基组可以在CP2K的data文件夹下

$ grep La *

复制代码

就能看到每一个基组或赝势文件中含有La元素的基组或赝势了，从这里选就可以。

CUTOFF、REL_CUTOFF、k点的设置严格来说要做收敛性测试

sobereva · 发表于 Post on 2024-9-4 11:57:33

参考北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）这页ppt了解镧系都能用什么

liuxp · 发表于 Post on 2024-9-4 15:45:51

Dempey 发表于 2024-9-4 10:45
使用的基组可以在CP2K的data文件夹下

就能看到每一个基组或赝势文件中含有La元素的基组或赝势了，从这里 ...

好的谢谢大佬我搜到了

liuxp · 发表于 Post on 2024-9-4 15:52:58

sobereva 发表于 2024-9-4 11:57
参考北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）这页ppt了解 ...

sob老师我看BASIS_MOLOPT_UCL没有O元素，那我可以写两行吗像下边那种
&DFT
BASIS_SET_FILE_NAME  BASIS_MOLOPT
BASIS_SET_FILE_NAME  BASIS_MOLOPT_UCL
POTENTIAL_FILE_NAME  POTENTIAL

sobereva · 发表于 Post on 2024-9-5 11:03:40

liuxp 发表于 2024-9-4 15:52
sob老师我看BASIS_MOLOPT_UCL没有O元素，那我可以写两行吗像下边那种
&DFT
BASIS_SET_FILE_NAME ...

可以。就应该这样
如果你用的不同元素的赝势不在同一个文件里，需要把赝势合并到同一个文件中

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

[CP2K] 结构中有La,Cu,O,我在使用cp2k的时候，系统显示BASIS_MOLOPT没有La