计算化学公社

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

[其它程序] 求助:li_conductivity.py脚本处理XDATCAR得到的write_msd中的时间和总时间对不上

[复制链接 Copy URL]

13

帖子

0

威望

85

eV
积分
98

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
# dt = timesteps * self.time_step * self.step_skip
设置了NBLOCK = 100 ,NSW = 30000,timesteps应该等于30000/100 = 300, self.time_step = POTIM = 2 ,所以self.step_skip应该=NBLOCK = 100(个人理解)
但是在 write_msd 中的时长是47000多fs,也不是一个整数,按理应该是60000fs,可能是时间间隔即self.step_skip设置错误。
有没有大佬解答一下,谢谢!




1

帖子

0

威望

7

eV
积分
8

Level 1 能力者

2#
发表于 Post on 2024-5-9 09:05:48 | 只看该作者 Only view this author
应该timesteps=NSW,还有计算的时候在diff = DiffusionAnalyzer.from_structures(traj,'Li',300,1,1,smoothed=False)这样就可以保留所有的数据,还有就是step_skip我现在也不太清楚是等于NBLOCK还是KBLOCK还是这两个的乘积,这里也没弄明白

13

帖子

0

威望

85

eV
积分
98

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2024-6-3 15:34:20 | 只看该作者 Only view this author
dreamflying 发表于 2024-5-9 09:05
应该timesteps=NSW,还有计算的时候在diff = DiffusionAnalyzer.from_structures(traj,'Li',300,1,1,smooth ...

谢谢,加上smoothed=false之后能够保留原有的数据了

13

帖子

0

威望

85

eV
积分
98

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2024-9-9 10:50:00 | 只看该作者 Only view this author
我爱学习d 发表于 2024-6-3 15:34
谢谢,加上smoothed=false之后能够保留原有的数据了

已经解决了 谢谢

本版积分规则 Credits rule

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

GMT+8, 2024-11-23 09:44 , Processed in 0.190804 second(s), 29 queries , Gzip On.

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