计算化学公社

标题: 结构中有La,Cu,O,我在使用cp2k的时候，系统显示BASIS_MOLOPT没有La [打印本页]

作者
Author: liuxp 时间: 2024-9-4 10:23
标题: 结构中有La,Cu,O,我在使用cp2k的时候，系统显示BASIS_MOLOPT没有La
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原子左右，后续我还会加上小分子进行振动分析计算反应的自由能。

作者
Author: Dempey 时间: 2024-9-4 10:45
使用的基组可以在CP2K的data文件夹下

$ grep La *

复制代码

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

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

作者
Author: sobereva 时间: 2024-9-4 11:57
参考北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）这页ppt了解镧系都能用什么

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

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

作者
Author: liuxp 时间: 2024-9-4 15:45

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

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

好的谢谢大佬我搜到了

作者
Author: liuxp 时间: 2024-9-4 15:52

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

作者
Author: sobereva 时间: 2024-9-5 11:03

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

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

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