计算化学公社

 找回密码 Forget password
 注册 Register
Views: 811|回复 Reply: 2
打印 Print 上一主题 Last thread 下一主题 Next thread

[CP2K] 求助:当不存在ADMM时,包含镧系元素体系的杂化泛函计算如何进行

[复制链接 Copy URL]

43

帖子

0

威望

1567

eV
积分
1610

Level 5 (御坂)

各位老师好,我在计算包含镧系元素的体系时,遇到了一个问题,相应的元素(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文件中不存在的元素自行定义?

202308011504037261..png (6.58 KB, 下载次数 Times of downloads: 9)

镧系元素缺少ADMM相关定义

镧系元素缺少ADMM相关定义

202308011536528813..png (75.63 KB, 下载次数 Times of downloads: 8)

ALL_POTENTIALS文件中相应参数的说明

ALL_POTENTIALS文件中相应参数的说明

5万

帖子

99

威望

5万

eV
积分
112390

管理员

公社社长

2#
发表于 Post on 2023-8-1 15:41:18 | 只看该作者 Only view this author
如果Lu只是体系中几个别原子,可以用其主基组当做ADMM辅助基组,此时对于Lu的描述既不会节约时间又不会降低精度

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

评分 Rate

参与人数
Participants 2
eV +6 收起 理由
Reason
njust-lbc + 1 好物!
qinjiu + 5 谢谢

查看全部评分 View all ratings

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

86

帖子

0

威望

521

eV
积分
607

Level 4 (黑子)

3#
发表于 Post on 2023-11-13 10:47:16 | 只看该作者 Only view this author
请问你是如何设置的镧系原子关于ADMM辅助基组?我可以看一下您的语句吗?

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-25 17:52 , Processed in 0.194936 second(s), 24 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list