请选择 进入手机版 | 继续访问电脑版
第9届北京科音分子动力学与GROMACS培训班将于4月17~20日于北京举办,请点击此链接查看培训详情,欢迎参加和相互转告!

计算化学公社

 找回密码
 现在注册!
查看: 1255|回复: 19

[CP2K] 请教:CP2K做过渡态、电子结构等静态计算时的可靠性

[复制链接]

125

帖子

0

威望

735

eV
积分
860

Level 4 (黑子)

发表于 2021-2-18 17:59:31 | 显示全部楼层 |阅读模式
本帖最后由 djjj148 于 2021-2-18 17:59 编辑

听说CP2K运行速度快,但是擅长做AIMD的它做静态计算的精度没有VASP和QE准。于是测试了自己的研究体系:同一个初猜结构下,CP2K在做结构优化速度十分的快,而且收敛性非常好(监控每一个离子步的能量变化发现几乎都是逐渐降低),反观VASP在结构优化时,最大受力总是会不断振荡才缓慢收敛。70多原子体系(晶格参数:16.7308 12.7134 25)的结构优化任务,CP2K(OT算法111K点)可以在1-3小时做完,而VASP(K点:1 2 1)需要1-2天,就算改为111 K点搭配vasp_gam版本也大概需要一天多。结构优化尚且如此,如果做过渡态计算或更大体系的话用时应该会更加悬殊。
然而,对于优化收敛的结构,分别做bader电荷分析发现,后处理VASP的输出文件生成的ACF.dat文件中比较正常,而CP2K的ACF.dat竟然出现了这样的情况:大部分C原子的CHARGE都是接近赝势中的价电子,没有差很多,然而小部分的C原子的CHARGE却是0,这似乎已经到了离谱的状态。相应的CP2K的inp文件已上传,请教一下大家:
1. 是否是因为inp中的设置的精度不够才导致如此“离谱”的结果?
2. 我在ACS.cata这样的期刊上看到的纯模拟的文章都是用VASP做的静态计算(如过渡态),CP2K做AIMD或者全部用VASP。然而,由于自己能用的计算硬件不够多和硬,还是眼馋CP2K的速度,不知道在恰当设置inp后,质量较好的期刊能否接受"CP2K做结构优化和过渡态计算得到的能量"?
3. 如果用CP2K优化极小点或过渡态结构后,用VASP算单点能量和进行下一步电子结构分析,应该可以提高结果的可靠性吧?

cp2k.inp

3.17 KB, 下载次数: 15

CP2K输入文件

3万

帖子

25

威望

3万

eV
积分
65563

管理员

公社社长+计算化学玩家

发表于 2021-2-18 19:47:42 | 显示全部楼层
1 不是精度的事。没有具体结构,也不好去判断。用CP2K没必要用bader电荷(AIM电荷)说事,本身CP2K自己就能算诸如Hirshfled-I、REPEAT这样通常比很垃圾的AIM电荷有意义得多的原子电荷,还可以结合Multiwfn算CM5电荷。而且CP2K用medium的print级别直接就输出Mulliken电荷,也能定性说明问题。

2 完全可以。TZV2PX-MOLOPT-GTH结合足够高的截断能已经能达到很高的精度。

3 没这个必要。CP2K结合恰当的设置产生的波函数做电子结构分析已经很好了。

顺带推荐个我之前审的一篇算沸石的文章,Materials Today Communications 26 (2021) 102028(https://doi.org/10.1016/j.mtcomm.2021.102028),全都是CP2K算的,还用Multiwfn做了IGM分析直观考察了沸石对芳香分子的吸附。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班。这些培训是计算化学快速入门以及全面系统性提升研究水平的最佳途径,培训各种相关信息见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395(已满),2号:466017436(已满)。3号:764390338(可加),合计8000人,讨论范畴相同
思想家公社的门口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!

125

帖子

0

威望

735

eV
积分
860

Level 4 (黑子)

 楼主| 发表于 2021-2-18 21:45:38 | 显示全部楼层
sobereva 发表于 2021-2-18 19:47
1 不是精度的事。没有具体结构,也不好去判断。用CP2K没必要用bader电荷(AIM电荷)说事,本身CP2K自己就能 ...

谢谢sob老师指点。还想请教一下,因为TZV2PX-MOLOPT-GTH是常用的BASIS_MOLOPT中精度最高的基组,我刚才测试了一下,其他参数不变,用TZV2PX-MOLOPT-GTH的耗时:
  Step     Update method      Time    Convergence         Total energy    Change
  ------------------------------------------------------------------------------
     1 OT CG       0.12E+00    9.1     0.01652627      -457.3860222786 -4.57E+02
     2 OT LS       0.17E+00    8.5                     -470.3961998629
     3 OT LS       0.22E+00    8.5                     -472.8613114815
     4 OT CG       0.18E+00   15.7     0.01742506      -473.7860931582 -1.64E+01

用DZVP-MOLOPT-SR-GTH:
     1 OT CG       0.12E+00    2.5     0.02632822      -459.3078393090 -4.59E+02
     2 OT LS       0.17E+00    1.7                     -471.9480597834
     3 OT LS       0.16E+00    1.6                     -472.8442457990
     4 OT CG       0.13E+00    3.1     0.03062221      -472.8461111027 -1.35E+01
耗时大概变成了原来3.6-5倍,幅度还是不小的。我想请教一下,像过渡态这样的优化任务用DZVP或TZVP-MOLOPT-GTH得到结构,然后用TZV2PX-MOLOPT-GTH和更高的截断能算单点和分析波函数可以吗?

2265

帖子

3

威望

8626

eV
积分
10951

Level 6 (一方通行)

Ab Initio Amateur

发表于 2021-2-18 21:50:02 | 显示全部楼层
本帖最后由 卡开发发 于 2021-2-18 22:00 编辑

1、CP2K的GPW是结合赝势的,所以赝波函数只有价层的且是smooth的,有可能不能正确识别NCP,也许得考虑缀加芯密度甚至是额外价层的赝波函数变换到全电子价层密度才行。实在不行也可以试试GAPW的方式。
3、并没有办法保证CP2K算得的结构能量差(反应能和能垒)能够和VASP一致,如果出现这种不一致可能会比较难解释(比如DOI: 10.1021/acs.jpcc.9b03554)。

如果只讨论精度,CP2K其实做了deltatest(https://www.cp2k.org/_media/events:2019_ghent:automation.pdf;其他程序的deltatest测试可以看https://molmod.ugent.be/deltacodesdft),可靠性方面我不做评价,可自行比对。
满招损,谦受益。热衷于理论和方法研究水平不高但欢迎讨论。

3万

帖子

25

威望

3万

eV
积分
65563

管理员

公社社长+计算化学玩家

发表于 2021-2-19 07:22:19 | 显示全部楼层
djjj148 发表于 2021-2-18 21:45
谢谢sob老师指点。还想请教一下,因为TZV2PX-MOLOPT-GTH是常用的BASIS_MOLOPT中精度最高的基组,我刚才测 ...

可以

另外,对于大多数波函数分析的目的,比如RDG、IRI、ELF、AIM分析等,其实TZVP-MOLOPT-GTH就够了,明显没有能量对基组敏感性那么高。

PS:对于在Multiwfn中做波函数分析,我建议额外用pob-TZVP基组做个GAPW单点计算得到molden文件。一方面pob-TZVP是基于片段收缩的,在Multiwfn里分析耗时低得多得多(比DZVP-MOLOPT-SR-GTH等价的GTF数能少好几倍),另一方面这样产生的波函数直接描述了内核电子,便于涉及到内核电子结构的分析。

评分

参与人数 1eV +5 收起 理由
djjj148 + 5 谢谢

查看全部评分

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班。这些培训是计算化学快速入门以及全面系统性提升研究水平的最佳途径,培训各种相关信息见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395(已满),2号:466017436(已满)。3号:764390338(可加),合计8000人,讨论范畴相同
思想家公社的门口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!

125

帖子

0

威望

735

eV
积分
860

Level 4 (黑子)

 楼主| 发表于 2021-2-19 09:47:18 | 显示全部楼层
sobereva 发表于 2021-2-19 07:22
可以

另外,对于大多数波函数分析的目的,比如RDG、IRI、ELF、AIM分析等,其实TZVP-MOLOPT-GTH就够了 ...

谢谢sob老师!请问现在Multiwfn是可以分析CP2K产生的基于周期性结构的波函数了吗?我看了一下好像还没有相应的博文。还是目前按照分析Gaussian和ORCA的fchk、molden文件一样直接载入Multiwfn中分析呢?如果是这样,晶格边界的电子性质是不是不可用呢?

125

帖子

0

威望

735

eV
积分
860

Level 4 (黑子)

 楼主| 发表于 2021-2-19 09:48:11 | 显示全部楼层
卡开发发 发表于 2021-2-18 21:50
1、CP2K的GPW是结合赝势的,所以赝波函数只有价层的且是smooth的,有可能不能正确识别NCP,也许得考虑缀加 ...

谢谢卡卡老师~!

2265

帖子

3

威望

8626

eV
积分
10951

Level 6 (一方通行)

Ab Initio Amateur

发表于 2021-2-19 10:09:37 | 显示全部楼层
djjj148 发表于 2021-2-19 09:47
谢谢sob老师!请问现在Multiwfn是可以分析CP2K产生的基于周期性结构的波函数了吗?我看了一下好像还没有 ...

最新的Multiwfn手册2.9.2节有提到。

评分

参与人数 1eV +3 收起 理由
djjj148 + 3 谢谢

查看全部评分

满招损,谦受益。热衷于理论和方法研究水平不高但欢迎讨论。

104

帖子

0

威望

346

eV
积分
450

Level 3 能力者

发表于 2021-2-19 12:20:31 | 显示全部楼层
可以用 QZVP开头的基组和适合的截断能获得cp2k能达到的最高精度

评分

参与人数 2eV +2 收起 理由
djjj148 + 1 谢谢
卡开发发 + 1 GAPW的情况下理论上能用到更大的。

查看全部评分

330

帖子

0

威望

1413

eV
积分
1743

Level 5 (御坂)

发表于 2021-2-19 12:53:59 | 显示全部楼层
CP2K的能量精度很可以的。我以前测过几个有机反应体系,在调好截断能的情况下计算两个构象间能量差。同个泛函、相同水平基组下,CP2K和Gaussian的dE差距甚小。

MOLOPT基组的使用细节我记得前几天有帖子。

2265

帖子

3

威望

8626

eV
积分
10951

Level 6 (一方通行)

Ab Initio Amateur

发表于 2021-2-19 13:35:11 | 显示全部楼层
k64_cc 发表于 2021-2-19 12:53
CP2K的能量精度很可以的。我以前测过几个有机反应体系,在调好截断能的情况下计算两个构象间能量差。同个泛 ...

如果只是CHO体系,从delta test的结果反映来看误差也确实不大,但到了含过渡金属或者其他元素的情况就不是很有保障了,比如前面提到的CO@ZnO的情况。

评分

参与人数 1eV +1 收起 理由
k64_cc + 1 那就很艰难了……

查看全部评分

满招损,谦受益。热衷于理论和方法研究水平不高但欢迎讨论。

330

帖子

0

威望

1413

eV
积分
1743

Level 5 (御坂)

发表于 2021-2-20 00:51:56 | 显示全部楼层
本帖最后由 k64_cc 于 2021-2-20 00:53 编辑
卡开发发 发表于 2021-2-19 13:35
如果只是CHO体系,从delta test的结果反映来看误差也确实不大,但到了含过渡金属或者其他元素的情况就不 ...

啊这个,我看他输入文件里也没有过渡金属,最多就一个B,就觉得应该还行吧……

评分

参与人数 1eV +1 收起 理由
卡开发发 + 1 被坑过23333

查看全部评分

3万

帖子

25

威望

3万

eV
积分
65563

管理员

公社社长+计算化学玩家

发表于 2021-2-20 15:23:14 | 显示全部楼层
djjj148 发表于 2021-2-19 09:47
谢谢sob老师!请问现在Multiwfn是可以分析CP2K产生的基于周期性结构的波函数了吗?我看了一下好像还没有 ...

最近会写一系列帖子。前阵子太忙没空写
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班。这些培训是计算化学快速入门以及全面系统性提升研究水平的最佳途径,培训各种相关信息见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395(已满),2号:466017436(已满)。3号:764390338(可加),合计8000人,讨论范畴相同
思想家公社的门口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!

3万

帖子

25

威望

3万

eV
积分
65563

管理员

公社社长+计算化学玩家

发表于 2021-2-20 15:34:43 | 显示全部楼层
顺带发一个最近看到的Chem. Sci., 2017, 8, 3554–3565的补充材料的对比,可见CP2K可以和量子化学程序吻合得相当好,即便是用GTH赝势做GPW也是如此

Clipboard01.png

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班。这些培训是计算化学快速入门以及全面系统性提升研究水平的最佳途径,培训各种相关信息见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395(已满),2号:466017436(已满)。3号:764390338(可加),合计8000人,讨论范畴相同
思想家公社的门口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!

125

帖子

0

威望

735

eV
积分
860

Level 4 (黑子)

 楼主| 发表于 2021-2-20 23:46:09 | 显示全部楼层
sobereva 发表于 2021-2-20 15:23
最近会写一系列帖子。前阵子太忙没空写

给卢老师打call,Multiwfn支持对CP2K的分析让这把瑞士军刀更为锋利了!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )

GMT+8, 2021-4-15 07:03 , Processed in 0.230855 second(s), 30 queries .

快速回复 返回顶部 返回列表