计算化学公社

标题: 求助：li_conductivity.py脚本处理XDATCAR得到的write_msd中的时间和总时间对不上 [打印本页]

作者
Author: 我爱学习d 时间: 2024-4-15 15:44
标题: 求助：li_conductivity.py脚本处理XDATCAR得到的write_msd中的时间和总时间对不上
# 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设置错误。
有没有大佬解答一下，谢谢！
(, 下载次数 Times of downloads: 16)
(, 下载次数 Times of downloads: 14)
(, 下载次数 Times of downloads: 18)

作者
Author: dreamflying 时间: 2024-5-9 09:05
应该timesteps=NSW，还有计算的时候在diff = DiffusionAnalyzer.from_structures(traj,'Li',300,1,1,smoothed=False)这样就可以保留所有的数据，还有就是step_skip我现在也不太清楚是等于NBLOCK还是KBLOCK还是这两个的乘积，这里也没弄明白

作者
Author: 我爱学习d 时间: 2024-6-3 15:34

dreamflying 发表于 2024-5-9 09:05
应该timesteps=NSW，还有计算的时候在diff = DiffusionAnalyzer.from_structures(traj,'Li',300,1,1,smooth ...

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

作者
Author: 我爱学习d 时间: 2024-9-9 10:50

我爱学习d 发表于 2024-6-3 15:34
谢谢，加上smoothed=false之后能够保留原有的数据了

已经解决了谢谢

欢迎光临计算化学公社 (http://bbs.keinsci.com/)