计算化学公社

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

[CP2K] CP2K跑MD盒子大小的问题

[复制链接 Copy URL]

1016

帖子

0

威望

2096

eV
积分
3112

Level 5 (御坂)

天上地下三山亖海五湖六合八荒九州之泛粤大典編輯委員會編外委員

跳转到指定楼层 Go to specific reply
楼主
各位老师,我要用CP2K跑几个MD,要加两个方向或三个方向的周期性边界条件。网上一些资料说CP2K设盒子不能太小,否则基函数会超出盒子边界,导致计算出问题。请问盒子至少要多大才行?如果开OT(orbital transformation),好像也要对盒子大小有要求,翻了OT的文章,没有提及盒子大小的问题。发帖求助。请问是否有相关资料。

我知道cutoff和relative cutoff可以自己测试一下,如果要测试盒子大小的影响,如果操作?



ORCA大法好!CP2K大法好!
嶺南粤音 泛粤典 https://jyutjam.org/

山河自落蒼梧月,風雨猶驅草木兵。
san huɔɞ ki lɔk tʰɔŋ ŋ ᵑgut,fʊŋ ʝi ʝiu kʰui tʰɔu ᵐbʊk pɪŋ

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

2#
发表于 Post on 2022-7-26 03:01:04 | 只看该作者 Only view this author
真空区最小需要多大,和体系特征(尤其是基函数指数、电子的弥散程度)有很大关系
可以做个能量随真空层厚度的收敛性测试,取能令能量收敛的厚度
北京科音自然科学研究中心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

1016

帖子

0

威望

2096

eV
积分
3112

Level 5 (御坂)

天上地下三山亖海五湖六合八荒九州之泛粤大典編輯委員會編外委員

3#
 楼主 Author| 发表于 Post on 2022-7-26 08:41:11 | 只看该作者 Only view this author
sobereva 发表于 2022-7-26 03:01
真空区最小需要多大,和体系特征(尤其是基函数指数、电子的弥散程度)有很大关系
可以做个能量随真空层厚 ...

谢谢!如果只有一个方向有真空区,其他方向是纯液态PBC,或者三个方向都是液态PBC,如何测试?
ORCA大法好!CP2K大法好!
嶺南粤音 泛粤典 https://jyutjam.org/

山河自落蒼梧月,風雨猶驅草木兵。
san huɔɞ ki lɔk tʰɔŋ ŋ ᵑgut,fʊŋ ʝi ʝiu kʰui tʰɔu ᵐbʊk pɪŋ

553

帖子

0

威望

3324

eV
积分
3877

Level 5 (御坂)

4#
发表于 Post on 2022-7-26 09:18:00 | 只看该作者 Only view this author
带周期性的部分没有边界问题,不带的才要考虑,具体自己调吧,有边界问题的时候SCF会给warning

1016

帖子

0

威望

2096

eV
积分
3112

Level 5 (御坂)

天上地下三山亖海五湖六合八荒九州之泛粤大典編輯委員會編外委員

5#
 楼主 Author| 发表于 Post on 2022-7-26 09:20:09 | 只看该作者 Only view this author
k64_cc 发表于 2022-7-26 09:18
带周期性的部分没有边界问题,不带的才要考虑,具体自己调吧,有边界问题的时候SCF会给warning

谢谢,我是否可以理解为,不带PBC的方向,在盒子边缘波函数要接近零,所以才有边界问题?
ORCA大法好!CP2K大法好!
嶺南粤音 泛粤典 https://jyutjam.org/

山河自落蒼梧月,風雨猶驅草木兵。
san huɔɞ ki lɔk tʰɔŋ ŋ ᵑgut,fʊŋ ʝi ʝiu kʰui tʰɔu ᵐbʊk pɪŋ

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

6#
发表于 Post on 2022-7-26 15:46:45 | 只看该作者 Only view this author
chands 发表于 2022-7-26 08:41
谢谢!如果只有一个方向有真空区,其他方向是纯液态PBC,或者三个方向都是液态PBC,如何测试?

仅需要对不考虑对称性但又设了真空区以定义平面波计算所用的盒子的那个方向,才需要做测试
北京科音自然科学研究中心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

1016

帖子

0

威望

2096

eV
积分
3112

Level 5 (御坂)

天上地下三山亖海五湖六合八荒九州之泛粤大典編輯委員會編外委員

7#
 楼主 Author| 发表于 Post on 2022-7-26 15:57:18 | 只看该作者 Only view this author
sobereva 发表于 2022-7-26 15:46
仅需要对不考虑对称性但又设了真空区以定义平面波计算所用的盒子的那个方向,才需要做测试

谢谢。我知道不考虑PBC的方向要加真空区,否则出问题。请问加真空区的原因什么能详细说一下吗?
ORCA大法好!CP2K大法好!
嶺南粤音 泛粤典 https://jyutjam.org/

山河自落蒼梧月,風雨猶驅草木兵。
san huɔɞ ki lɔk tʰɔŋ ŋ ᵑgut,fʊŋ ʝi ʝiu kʰui tʰɔu ᵐbʊk pɪŋ

20

帖子

0

威望

954

eV
积分
974

Level 4 (黑子)

8#
发表于 Post on 2022-7-26 17:35:49 | 只看该作者 Only view this author
老师构建的晶体模型做MD一定得建真空层吗?

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

9#
发表于 Post on 2022-7-26 18:03:24 | 只看该作者 Only view this author
chands 发表于 2022-7-26 15:57
谢谢。我知道不考虑PBC的方向要加真空区,否则出问题。请问加真空区的原因什么能详细说一下吗?

因为CP2K用平面波当辅助基组,注定得设置盒子才能用平面波基函数。在非周期性方向,真空区没有或者太小的话,都等于错误地截断了电子密度
北京科音自然科学研究中心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万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

10#
发表于 Post on 2022-7-26 18:03:55 | 只看该作者 Only view this author
LaplaceL 发表于 2022-7-26 17:35
老师构建的晶体模型做MD一定得建真空层吗?

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

20

帖子

0

威望

954

eV
积分
974

Level 4 (黑子)

11#
发表于 Post on 2022-7-26 20:08:39 | 只看该作者 Only view this author

好的,看了老师回复其他人的我也明白了

236

帖子

0

威望

1227

eV
积分
1463

Level 4 (黑子)

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

12#
发表于 Post on 2025-3-5 20:09:49 | 只看该作者 Only view this author
sobereva 发表于 2022-7-26 15:46
仅需要对不考虑对称性但又设了真空区以定义平面波计算所用的盒子的那个方向,才需要做测试

老师,看您前面的表述:1.如果做三维周期性的FPMD(pbe搭配MOLOPT系列基组),是不是就不需要考虑盒子大小测试问题。只有对非周期性方向才做盒子大小测试,或者说真空区高度测试。
2.此外,FPMD中还需要进行截断能测试吗,一般大体系(550-1000+)的FPMD比较耗时,很少有人做FPMD的CUTOFF测试,我看都是计算能量时才进行CUTOFF测试。
目前专攻:
基于DFT、MD模拟的自由能计算
基于过渡态理论的自由能垒和反应速率常数计算
基于MD模拟的交联聚合物计算

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

13#
发表于 Post on 2025-3-6 02:17:05 | 只看该作者 Only view this author
JCenter 发表于 2025-3-5 20:09
老师,看您前面的表述:1.如果做三维周期性的FPMD(pbe搭配MOLOPT系列基组),是不是就不需要考虑盒子大 ...

1 看具体测试指什么目的。真空区厚度的收敛性测试显然对于本质是周期性方向上是无意义的

2 一般根据基组的最大高斯函数指数判断就够了,一般不需要测试
北京科音自然科学研究中心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模拟爱好者

14#
发表于 Post on 2025-3-6 13:47:27 | 只看该作者 Only view this author
sobereva 发表于 2025-3-6 02:17
1 看具体测试指什么目的。真空区厚度的收敛性测试显然对于本质是周期性方向上是无意义的

2 一般根据基 ...

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

本版积分规则 Credits rule

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

GMT+8, 2025-8-15 19:41 , Processed in 0.478535 second(s), 20 queries , Gzip On.

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