计算化学公社

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

[CP2K] 单斜盒子NPT模拟出现温度异常突变

[复制链接 Copy URL]

37

帖子

0

威望

203

eV
积分
240

Level 3 能力者

大佬们,我最近最NPT模拟单斜晶胞的NPT分子动力学模拟,第二步温度直接飙升到了4000K,是因为必须使用正交盒子吗,还是我的inp文件有问题啊,感谢大佬解答,之前用相同的设置跑正交盒子就没有遇到问题。

cp2k.inp

6.37 KB, 下载次数 Times of downloads: 9

cp2k-1.ener

58.34 KB, 下载次数 Times of downloads: 2

5万

帖子

99

威望

5万

eV
积分
112499

管理员

公社社长

2#
发表于 Post on 2024-5-3 20:13:36 | 只看该作者 Only view this author
跟单斜没必然关系
这么小晶胞做MD根本没意义,会等效施加严重的虚假的周期性约束
北京科音自然科学研究中心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!

37

帖子

0

威望

203

eV
积分
240

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2024-5-4 16:57:37 | 只看该作者 Only view this author
sobereva 发表于 2024-5-3 20:13
跟单斜没必然关系
这么小晶胞做MD根本没意义,会等效施加严重的虚假的周期性约束

sob老师,我这个已经扩胞3x3x3了啊,500多个原子了,还是小吗

5万

帖子

99

威望

5万

eV
积分
112499

管理员

公社社长

4#
发表于 Post on 2024-5-4 17:32:10 | 只看该作者 Only view this author
wjybit 发表于 2024-5-4 16:57
sob老师,我这个已经扩胞3x3x3了啊,500多个原子了,还是小吗

我没注意到MULTIPLE_UNIT_CELL  3  3  3,当前可以

尝试先用NVT跑,或者改用DFT跑。
北京科音自然科学研究中心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!

37

帖子

0

威望

203

eV
积分
240

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2024-5-4 18:01:40 | 只看该作者 Only view this author
sobereva 发表于 2024-5-4 17:32
我没注意到MULTIPLE_UNIT_CELL  3  3  3,当前可以

尝试先用NVT跑,或者改用DFT跑。

好的sob老师

37

帖子

0

威望

203

eV
积分
240

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2024-5-4 20:05:26 | 只看该作者 Only view this author
本帖最后由 wjybit 于 2024-5-4 20:25 编辑
sobereva 发表于 2024-5-4 17:32
我没注意到MULTIPLE_UNIT_CELL  3  3  3,当前可以

尝试先用NVT跑,或者改用DFT跑。

sob老师,我如果把使用multiwfn对这个单斜晶胞进行晶胞变换,就是a和b之间的变换,就会导致NPT模拟的温度突变,是multiwfn无法对单斜晶胞坐标轴交换吗

5万

帖子

99

威望

5万

eV
积分
112499

管理员

公社社长

7#
发表于 Post on 2024-5-5 23:43:15 | 只看该作者 Only view this author
wjybit 发表于 2024-5-4 20:05
sob老师,我如果把使用multiwfn对这个单斜晶胞进行晶胞变换,就是a和b之间的变换,就会导致NPT模拟的温度 ...

那个功能目前仅限正交晶胞。
变换完了是否合理,直接在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!

11

帖子

0

威望

89

eV
积分
100

Level 2 能力者

8#
发表于 Post on 2024-9-19 14:55:47 | 只看该作者 Only view this author

同学,你这个最后搞定了吗?我也是NPT出现了和你类似的情况,第二步开始就飙升到6000-7000K。换了poisson solver也是一样的问题,没法搞定。还是说这种情况只能跑NVT了?

本版积分规则 Credits rule

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

GMT+8, 2024-11-27 12:37 , Processed in 0.225413 second(s), 24 queries , Gzip On.

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