计算化学公社

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

[CP2K] 关于CP2K计算简单离子水溶液体系AIMD的问题

[复制链接 Copy URL]

6

帖子

0

威望

362

eV
积分
368

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
各位老师好!学生最近想使用CP2K计算AIMD来研究简单离子在水溶液中的性质,建模参考sob老师关于使用packmol构建溶剂盒子的博文,AIMD在NVT系综下进行,使用CSVR热浴。有一些问题想要请教各位老师:
1. 我想研究反离子对金属阳离子在水溶液中的性质是否具有一定影响,能否通过对比包含/不包含反离子的溶液体系的AIMD结果实现这一目的?因为对于周期性的溶液盒子,如果不放置反离子来平衡电荷,整个盒子是带电的,我不确定这样做是否正确。如果必须要平衡电荷,请问CP2K中有没有类似加入背景电荷之类的设置来实现这一目的(我没有找到类似的关键词)?
2. 我查阅资料和论坛中的帖子,MD的步长应设置为小于体系振动频率最高的键周期的10%左右,我做了一部分计算(流程为:结构优化-TIMECON=100下的预平衡-TIMECON=200下的较长时间模拟)步长设置为1fs,查阅了一些模拟液体水的文献发现大多数步长设置为0.5fs,1fs的步长对于我所设置的溶液体系来说是否有点太长了?我可以接着使用步长1fs的轨迹进行分析,还是说需要在0.5fs的步长下重新平衡并进行较长时间模拟收集数据?
学生入门时间不长,如果有什么表述的不正确的地方希望各位老师多多指正,非常感谢各位的宝贵意见!

6万

帖子

99

威望

5万

eV
积分
120086

管理员

公社社长

2#
发表于 Post on 2022-11-23 01:49:35 | 只看该作者 Only view this author
1 不适合。抗衡离子总是会出现的

2 常温下,用1fs可以接受。但对于涉及到成键/断键的模拟,1fs就太大了
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6

帖子

0

威望

362

eV
积分
368

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2022-11-23 10:24:16 | 只看该作者 Only view this author
sobereva 发表于 2022-11-23 01:49
1 不适合。抗衡离子总是会出现的

2 常温下,用1fs可以接受。但对于涉及到成键/断键的模拟,1fs就太大了

感谢sob老师回复!学生还有一些问题想请教您和各位老师:
1. 不包含抗衡离子的体系确实实际上不存在,我发现文献JPCA,2018,122,1905中使用CPMD研究Bi(III)的水解时直接研究了Bi(H2O)83+ + 56H2O的体系而没有考虑反离子的存在,想请教一下您直接对一个带电荷的周期性溶液盒子做AIMD模拟理论上是否合理,如果这样做理论上没问题则作为一个虚构体系,也许对反离子存在时的影响能够有一定的参照。
2. 如果我设置更高的温度(比如50/70/90℃),是否需要进一步缩短步长(如0.5fs)?您说的“涉及成键/断键的模拟,1fs太大了”是指1fs的步长太长,观察不到成/断键过程,还是说其他方面的考虑?
3. 关于模拟的步长,我可以通过什么样的标准来判断我选择的步长对一个体系的模拟是合适/不合适的呢?
感谢您和各位老师的指导!

42

帖子

0

威望

926

eV
积分
968

Level 4 (黑子)

4#
发表于 Post on 2022-11-23 22:05:32 | 只看该作者 Only view this author
关于步长的这里,基本上涉及到水溶液如果想要计算精确的化还是需要0.5fs的步长。步长一定程度上与原子震动频率挂钩,如果是其他重一点的原子用1fs都还可以,而氢原子太轻了,震动频率较高,需要用到0.5fs的步长

6万

帖子

99

威望

5万

eV
积分
120086

管理员

公社社长

5#
发表于 Post on 2022-11-24 01:11:30 | 只看该作者 Only view this author
Frobenius 发表于 2022-11-23 10:24
感谢sob老师回复!学生还有一些问题想请教您和各位老师:
1. 不包含抗衡离子的体系确实实际上不存在,我 ...

1 AIMD时候盒子有净电荷不意味着模拟的动力学现象就有什么问题、没实际意义。比如很多文章模拟水合电子,盒子也是带净电荷的。如果你的重点是动力学行为,而不是能量层面的问题,适当的时候还是能对比的。

2 温度越高,原子运动速度越快,步长建议减小以确保精度和动力学稳定性。不过只是相差几十度不至于有多大影响。对于涉及H的成键断键,出现现象时会导致瞬间受力、速度会很大,1fs就太糙了,可能导致明显的错误的动力学行为、产物,甚至模拟崩溃。

3 有很多因素。没经验的话就看近10年的发在像样刊物上的类似研究的文章。如果你怕被审稿人吐槽,计算资源又不是很匮乏的话,宁可浪费一些、用小一点,比如0.5 fs。1 fs属于不涉及成键断键的情况下常温模拟带H的体系的最大可接受步长(而且没用SHAKE之类算法约束住氢相关的键的情况下)。另外,结果随步长做收敛性测试是最能说明问题的,但这种测试本身非常耗时,所以一般不考虑。撑死了用0.025 fs也绝对足够了。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

6

帖子

0

威望

362

eV
积分
368

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2022-11-24 09:56:35 | 只看该作者 Only view this author
asd9998876 发表于 2022-11-23 22:05
关于步长的这里,基本上涉及到水溶液如果想要计算精确的化还是需要0.5fs的步长。步长一定程度上与原子震动 ...

感谢您的回答!

6

帖子

0

威望

362

eV
积分
368

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2022-11-24 09:56:50 | 只看该作者 Only view this author
sobereva 发表于 2022-11-24 01:11
1 AIMD时候盒子有净电荷不意味着模拟的动力学现象就有什么问题、没实际意义。比如很多文章模拟水合电子, ...

感谢sob老师指点!

236

帖子

0

威望

1227

eV
积分
1463

Level 4 (黑子)

实验组内的DFT计算、第一性原理、MD模拟爱好者

8#
发表于 Post on 2025-5-12 12:14:03 | 只看该作者 Only view this author
sobereva 发表于 2022-11-24 01:11
1 AIMD时候盒子有净电荷不意味着模拟的动力学现象就有什么问题、没实际意义。比如很多文章模拟水合电子, ...

老师,请教您一个问题:做fpmd时,体系中如果有钠离子,但其数量也就两三个,使用PBE/DZVP-MOLOPT-SR-GTH,其截断能也需要增加至1000吗,跑MD是不是太奢侈。
目前专攻:
基于DFT、MD模拟的自由能计算
基于过渡态理论的自由能垒和反应速率常数计算
基于MD模拟的交联聚合物计算

882

帖子

3

威望

1668

eV
积分
2610

Level 5 (御坂)

傻傻的木瓜

9#
发表于 Post on 2025-5-12 12:52:07 | 只看该作者 Only view this author
JCenter 发表于 2025-5-12 12:14
老师,请教您一个问题:做fpmd时,体系中如果有钠离子,但其数量也就两三个,使用PBE/DZVP-MOLOPT-SR-GTH ...

看你跑MD的主要目的,如果钠仅用于平衡少量电荷而并非(能量与受力计算,动力学行为等的)重点关注对象,稍微降低点也勉强能接受。但以我这边对NaCl水溶液测试收敛性的经验,其实1000 Ry对钠的受力计算而言仍然有不小误差;要问最低可接受的设置,我觉得800 Ry可能是底限了。

(挖个坑,后面我可能会就CP2K里截断能的问题写写帖子)
√546=23.36664289109

236

帖子

0

威望

1227

eV
积分
1463

Level 4 (黑子)

实验组内的DFT计算、第一性原理、MD模拟爱好者

10#
发表于 Post on 2025-5-12 16:41:49 | 只看该作者 Only view this author
本帖最后由 JCenter 于 2025-5-12 16:45 编辑
Uus/pMeC6H4-/キ 发表于 2025-5-12 12:52
看你跑MD的主要目的,如果钠仅用于平衡少量电荷而并非(能量与受力计算,动力学行为等的)重点关注对象, ...

好嘞,那您之前有尝试过在a.GPW方法下使用赝式基组和高截断能1000 Ry,b.GAPW方法下对钠离子使用全电子基组比如pob-dzvp--rev2,其余部分使用赝式基组,降低截断能至400 Ry,这两种策略下哪种耗时更短、更快得到结果呢,实在不行,就直接换成钾离子,一了百了,不需要那么高的截断能
目前专攻:
基于DFT、MD模拟的自由能计算
基于过渡态理论的自由能垒和反应速率常数计算
基于MD模拟的交联聚合物计算

882

帖子

3

威望

1668

eV
积分
2610

Level 5 (御坂)

傻傻的木瓜

11#
发表于 Post on 2025-5-12 17:10:45 | 只看该作者 Only view this author
JCenter 发表于 2025-5-12 16:41
好嘞,那您之前有尝试过在a.GPW方法下使用赝式基组和高截断能1000 Ry,b.GAPW方法下对钠离子使用全电子基 ...

我觉得普通体系的简单AIMD任务用混合基组还是有点风险,得拿出比“部分元素的赝势和基组比较特别”更充分有力的理由来服众才行。

能否把钠换成钾也得看实际体系和过程特征,如果反应可能发生在碱金属离子第一配位层的配体之间,二者的半径差异是会有影响的。
√546=23.36664289109

6万

帖子

99

威望

5万

eV
积分
120086

管理员

公社社长

12#
发表于 Post on 2025-5-12 19:14:36 | 只看该作者 Only view this author
JCenter 发表于 2025-5-12 12:14
老师,请教您一个问题:做fpmd时,体系中如果有钠离子,但其数量也就两三个,使用PBE/DZVP-MOLOPT-SR-GTH ...

没什么理由不对Na用全电子基组,CUTOFF省多了去了
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

236

帖子

0

威望

1227

eV
积分
1463

Level 4 (黑子)

实验组内的DFT计算、第一性原理、MD模拟爱好者

13#
发表于 Post on 2025-5-12 19:59:46 | 只看该作者 Only view this author
sobereva 发表于 2025-5-12 19:14
没什么理由不对Na用全电子基组,CUTOFF省多了去了

好嘞,感谢卢老师
目前专攻:
基于DFT、MD模拟的自由能计算
基于过渡态理论的自由能垒和反应速率常数计算
基于MD模拟的交联聚合物计算

236

帖子

0

威望

1227

eV
积分
1463

Level 4 (黑子)

实验组内的DFT计算、第一性原理、MD模拟爱好者

14#
发表于 Post on 2025-5-12 20:01:45 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2025-5-12 17:10
我觉得普通体系的简单AIMD任务用混合基组还是有点风险,得拿出比“部分元素的赝势和基组比较特别”更充分 ...

好滴,感谢大佬。我再研究琢磨琢磨
目前专攻:
基于DFT、MD模拟的自由能计算
基于过渡态理论的自由能垒和反应速率常数计算
基于MD模拟的交联聚合物计算

本版积分规则 Credits rule

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

GMT+8, 2025-8-13 15:27 , Processed in 0.161985 second(s), 20 queries , Gzip On.

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