计算化学公社

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

[Multiwfn使用咨询] 请问在做弱相互作用NCI分析时,cp2k中后周期过渡金属元素的基组该选哪个呢?

[复制链接 Copy URL]

29

帖子

0

威望

183

eV
积分
212

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
sob老师教程里提到在做NCI分析周期性结构时,cp2k计算要选用6-311G**之类的基组,可是后周期元素例如Mo并没有这些基组,请问这些元素该用哪个基组计算呢?
我的体系包含270个原子(MoS2),用MOLOPT基组和GTH-PBE赝势得到的molden文件在用Multiwfn分析时非常慢,几个小时都分析不完。

5万

帖子

99

威望

5万

eV
积分
112554

管理员

公社社长

2#
发表于 Post on 2022-11-6 09:43:50 | 只看该作者 Only view this author
说具体当前用的什么MOLOPT基组,以及计算条件
另外,恰当设置计算的盒子范围,只包含感兴趣的区域,能大幅节约时间
北京科音自然科学研究中心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!

29

帖子

0

威望

183

eV
积分
212

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2022-11-6 10:26:16 | 只看该作者 Only view this author
sobereva 发表于 2022-11-6 09:43
说具体当前用的什么MOLOPT基组,以及计算条件
另外,恰当设置计算的盒子范围,只包含感兴趣的区域,能大幅 ...

sob老师您好,目前S用的是DZVP-MOLOPT-SR-GTH-q6和GTH-PBE-q6,Mo用的是DZVP-MOLOPT-SR-GTH-q14和GTH-PBE-q14。我研究的是块体,整个盒子都要研究的。用的是48核的服务器,cp2k任务几分钟就可以完成,但是Multiwfn的NCI分析要5个小时以上。

5万

帖子

99

威望

5万

eV
积分
112554

管理员

公社社长

4#
发表于 Post on 2022-11-6 10:57:54 | 只看该作者 Only view this author
yufeng 发表于 2022-11-6 10:26
sob老师您好,目前S用的是DZVP-MOLOPT-SR-GTH-q6和GTH-PBE-q6,Mo用的是DZVP-MOLOPT-SR-GTH-q14和GTH-PBE ...

改用TZVP-GTH

恰当设置盒子范围
北京科音自然科学研究中心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!

29

帖子

0

威望

183

eV
积分
212

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2022-11-6 11:09:11 | 只看该作者 Only view this author
本帖最后由 yufeng 于 2022-11-6 12:03 编辑
sobereva 发表于 2022-11-6 10:57
改用TZVP-GTH

恰当设置盒子范围

谢谢您的回复!TZVP精度更高,不是会导致计算更慢吗?是Multiwfn对于更好的波函数分析时间可以减少吗?

我现在S原子用TZVP-MOLOPT-PBE-GTH-q6,Mo原子用TZVP-MOLOPT-PBE-GTH-q14基组,然后格子大小是16*16*4 埃,格点间距用0.2,Multiwfn分析速度还是非常慢,半个小时只能分析不到8%。

GTH_BASIS_SETS中没有Mo的TZVP-GTH基组,我只好用BASIS_MOLOPT_UZH中的TZVP-MOLOPT-PBE-GTH基组。



5万

帖子

99

威望

5万

eV
积分
112554

管理员

公社社长

6#
发表于 Post on 2022-11-6 13:47:37 | 只看该作者 Only view this author
yufeng 发表于 2022-11-6 11:09
谢谢您的回复!TZVP精度更高,不是会导致计算更慢吗?是Multiwfn对于更好的波函数分析时间可以减少吗?
...

TZVP-GTH和TZVP-MOLOPT-GTH完全是两码事,S完全没必要用昂贵的TZVP-MOLOPT-GTH
也可以对S用def2-SVP
并且确保Multiwfn利用所有物理核心数并行了

北京科音自然科学研究中心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!

5万

帖子

99

威望

5万

eV
积分
112554

管理员

公社社长

7#
发表于 Post on 2022-11-6 14:09:16 | 只看该作者 Only view this author
还有个节约耗时的办法,把settings.ini里的expcutoff_PBC改为-10,会牺牲一些精度
北京科音自然科学研究中心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!

29

帖子

0

威望

183

eV
积分
212

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2022-11-6 20:36:31 | 只看该作者 Only view this author
本帖最后由 yufeng 于 2022-11-6 20:37 编辑
sobereva 发表于 2022-11-6 13:47
TZVP-GTH和TZVP-MOLOPT-GTH完全是两码事,S完全没必要用昂贵的TZVP-MOLOPT-GTH
也可以对S用def2-SVP
并 ...

好的,谢谢sob老师。我现在用的S原子是TZVP-GTH,Mo原子是DZVP-MOLOPT-SR-GTH,应该不需要改了吧?
然后Multiwfn并行用的-nt 48。
16*16*4 埃的格子,不损失精度的情况下,Multiwfn大概40分钟能分析完,这样正常了吗?

5万

帖子

99

威望

5万

eV
积分
112554

管理员

公社社长

9#
发表于 Post on 2022-11-7 10:12:28 | 只看该作者 Only view this author
yufeng 发表于 2022-11-6 20:36
好的,谢谢sob老师。我现在用的S原子是TZVP-GTH,Mo原子是DZVP-MOLOPT-SR-GTH,应该不需要改了吧?
然后 ...

正常
北京科音自然科学研究中心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!

本版积分规则 Credits rule

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

GMT+8, 2024-11-27 22:25 , Processed in 4.906841 second(s), 21 queries , Gzip On.

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