计算化学公社

标题: 求助：当不存在ADMM时，包含镧系元素体系的杂化泛函计算如何进行 [打印本页]

作者
Author: qinjiu 时间: 2023-8-1 15:37
标题: 求助：当不存在ADMM时，包含镧系元素体系的杂化泛函计算如何进行
各位老师好，我在计算包含镧系元素的体系时，遇到了一个问题，相应的元素（Lu）没有基组和赝势的文件，无法进行杂化泛函级别的计算。

我做了以下的尝试：
1. 在cp2k 2023.1 /data/目录下，grep Lu *，找到Lu 相关的基组有：BASIS_MOLOPT_LnPP2:Lu DZVP-MOLOPT-GTH-q11 系列； BASIS_MOLOPT_UCL:Lu SZV-MOLOPT-SR-GTH SZV-MOLOPT-SR-GTH-q25 系列；赝势有： GTH-PBE、GTH-GGA、GTH-SCAN、GTH-MGGA、GTH-PBE0、GTH-HYB；不存在辅助基组。
2. 我研究的体系不包含镧系元素时，121核并行杂化泛函+DZVP-MOLOPT-GTH-SR基组+ADMM+GPW计算分子动力学，平均20s一步；包含镧系元素之后，使用杂化泛函+GTH赝势基组，由于libint积分库没有对这种情况优化，缺少ADMM时计算效率极低（几万秒才能一步），所以几乎不能计算。
3. 于是我考虑使用GAPW+全电子基组进行计算，我在BSE上找到了想要元素的基组def2svp (cp2k格式)，在 &DFT 模块中定义了BASIS_SET_FILE_NAME，但是ALL_POTENTIALS文件中没有Lu的定义，计算报错The requested atomic potential <ALL> for element <LU> was not found in the potential file <POTENTIAL> 。

请教各位老师，包含镧系元素体系的杂化泛函计算，是否可以在CP2K中完成？如果可行的话，我应该修改哪些模块？是否可以对ALL_POTENTIALS文件中不存在的元素自行定义？

作者
Author: sobereva 时间: 2023-8-1 15:41
如果Lu只是体系中几个别原子，可以用其主基组当做ADMM辅助基组，此时对于Lu的描述既不会节约时间又不会降低精度

def2-svp在CP2K里没法用于Lu，对Lu是赝势基组，搭配的是SDD赝势，然而CP2K不支持SDD赝势

作者
Author: halouhapily 时间: 2023-11-13 10:47
请问你是如何设置的镧系原子关于ADMM辅助基组？我可以看一下您的语句吗？

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