老八蜜汁小憨包 发表于 2020-12-17 11:03 vaspkit做的没问题啊 |
|
您好，请问如何得到均方位移啊，用什么软件呢？我用的vaspkit报错了 |
k64_cc 发表于 2020-11-10 16:02 对，体系一大就很费劲了，唉，我先把20ps的跑出来我们瞅瞅，要是还不行就得想办法扩胞试试，就是太费时间了，有点耗不起，感谢大佬回答  |
maokun_wu 发表于 2020-11-10 15:48 cell太小，周期性影响很大，就很危险。你的体系我也不太了解，不过加尺寸会不会负担不起？ |
k64_cc 发表于 2020-11-10 15:27 我也觉得很危险，跑出来都不知道对不对，我模拟的是一个原子在4*4的二维材料表面上的迁移，体系一共40多个原子，我看网上跑MD要100多个原子，但是也有文献拿五六十个原子来跑的，我的体系是不是太小了？刚才您说的这个体系尺寸设计什么意思啊？方便解释一下么？感谢感谢 |
maokun_wu 发表于 2020-11-10 14:37 听起来有点危险……你确定你体系尺寸设计的没问题吗？ 等会，我也不知道你跑的是啥体系啊。 |
k64_cc 发表于 2020-11-10 13:27 我刚看了下10ps下跑的结构，基本过了5ps骨架就发生了较大变化，但是本征骨架下也就是没有扩散离子的时候结构是没问题的，想请问一下，这种情况是MD跑的问题还是说可能有了扩散粒子骨架不稳定？ |
k64_cc 发表于 2020-11-10 13:27 感谢您的回答，看到过莫一非老师关于AIMD的文献，基本时间超过70%由于统计误差就不准了，我先算一下20ps的，然后把结果贴出来大家看看，感谢感谢 |
本帖最后由 k64_cc 于 2020-11-10 13:28 编辑 wzkchem5 发表于 2020-11-10 13:18 dt>8 ps的部分其实只是开头2 ps和最后2 ps的位移差，data也过分少了，万一他体系里有个什么行为需要2 ps完成，那MSD长这样就完全不奇怪；或者他体系是个cage，粒子在8ps的时候撞墙弹回来，也可能搞成这样。说实话拿10 ps算MSD本身就不太稳当，小分子弛豫时间搞不好也就是这个量级，说不定10 ps全跑完都没能independent，多跑点吧还是…… |
k64_cc 发表于 2020-11-10 13:15 那请问一下，MSD过了1/2 t以后只是均值漂移，噪音并没有增大，也是正常的吗？ 如果只是因为可用的数据点变少，噪音应该也变大才对吧 |
| 10 ps的data，5 ps之后每个delta t能用到的data是递减的，越往后越不准。过了1/2 t之后的data就不用看了，你这个msd没毛病。 |
wzkchem5 发表于 2020-11-10 12:14 多谢多谢，我试试，然后和您说一下结果 |
本帖最后由 wzkchem5 于 2020-11-10 13:50 编辑 maokun_wu 发表于 2020-11-10 12:07 这样，你试一下，如果跑20ps，MSD仍然是8ps的时候最大吗？ |
wzkchem5 发表于 2020-11-10 12:00 应该不会，我用的vaspkit，重复过LiGPS的文献，是没有问题的，但是换到自己的体系MSD就这样了，唉，太惆怅了 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-18 09:09 , Processed in 0.200832 second(s), 31 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites