计算化学公社

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

[CP2K] 求助硅酸盐融体的AIMD模拟

[复制链接 Copy URL]

53

帖子

0

威望

578

eV
积分
631

Level 4 (黑子)

各位老师好,我最近利用CP2K跑NaAlSi3O8+H2O熔体的模拟,温度设为1123K,压力为3000bar,为了模拟地壳中硅酸盐的熔体的特征。现在跑了20ps左右,看了一下各元素的点的trajectory, 发现除了氧原子外,其他的原子没有跑满整个盒了,仅在有限的空间内运动。这样可能最终得到的统计结果可能没有啥意义。升高温度能加速粒子的运动,能跑满整个盒子,又与实际温度不符,担心会影响结果。请问各位老师有什么是建议吗?

202503172036403454..png (684.33 KB, 下载次数 Times of downloads: 0)

202503172036403454..png

6万

帖子

99

威望

5万

eV
积分
120117

管理员

公社社长

2#
发表于 Post on 2025-3-18 06:57:35 | 只看该作者 Only view this author
区区20ps只能跑成这样,哪怕你升高温度也没用。要想看起来概率密度分布均匀必须跑长得多的时间,得ns数量级(用合适的力场/势函数)
北京科音自然科学研究中心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

53

帖子

0

威望

578

eV
积分
631

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2025-3-18 08:40:55 | 只看该作者 Only view this author
sobereva 发表于 2025-3-18 06:57
区区20ps只能跑成这样,哪怕你升高温度也没用。要想看起来概率密度分布均匀必须跑长得多的时间,得ns数量级 ...

谢谢sob老师的回答。盒子里有一个过渡金属原子,我用的是DFT方法。我还想问一下,如果要得到金属元素关于O的rdf和cn,必须得要各元素在盒子内密度分布均匀吗?

6万

帖子

99

威望

5万

eV
积分
120117

管理员

公社社长

4#
发表于 Post on 2025-3-19 10:05:44 | 只看该作者 Only view this author
laifeng20151107 发表于 2025-3-18 08:40
谢谢sob老师的回答。盒子里有一个过渡金属原子,我用的是DFT方法。我还想问一下,如果要得到金属元素关于 ...

没有绝对的必要性
北京科音自然科学研究中心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

53

帖子

0

威望

578

eV
积分
631

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2025-3-19 14:21:24 | 只看该作者 Only view this author
sobereva 发表于 2025-3-19 10:05
没有绝对的必要性

明白,谢谢sob老师,那熔体体系只要体系的势能、温度达到平衡,就可以进入production阶段了吧

6万

帖子

99

威望

5万

eV
积分
120117

管理员

公社社长

6#
发表于 Post on 2025-3-19 14:26:56 | 只看该作者 Only view this author
laifeng20151107 发表于 2025-3-19 14:21
明白,谢谢sob老师,那熔体体系只要体系的势能、温度达到平衡,就可以进入production阶段了吧

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

本版积分规则 Credits rule

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

GMT+8, 2025-8-14 21:22 , Processed in 0.168649 second(s), 23 queries , Gzip On.

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