计算化学公社

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

[Lammps] conpute_vacf命令计算速度自相关函数求助

[复制链接 Copy URL]

24

帖子

0

威望

117

eV
积分
141

Level 2 能力者

   请问一下各位老师,我用lammps的conpute_vacf命令计算的VACF,然后根据数据作图,发现并不一样,请问一下各位我应该怎样改正,另外计算vacf的命令如下:compute VACF  all vacffix VACF all ave/time 1 12 ${ndump} c_VACF[1] c_VACF[2] c_VACF[3] c_VACF[4] file vacf.lammps



QQ图片20200914205633.png (118.67 KB, 下载次数 Times of downloads: 19)

QQ图片20200914205633.png

QQ图片20200914205633.png (118.67 KB, 下载次数 Times of downloads: 21)

左图是文献,右图是自己画的

左图是文献,右图是自己画的

24

帖子

0

威望

117

eV
积分
141

Level 2 能力者

2#
 楼主 Author| 发表于 Post on 2020-9-14 21:01:22 | 只看该作者 Only view this author
compute VACF  all vacf
fix VACF all ave/time 1 12 ${ndump} c_VACF[1] c_VACF[2] c_VACF[3] c_VACF[4] file vacf.lammps
这是命令行,发现发的串行了

33

帖子

0

威望

1980

eV
积分
2013

Level 5 (御坂)

3#
发表于 Post on 2020-9-15 12:15:55 | 只看该作者 Only view this author
txh 发表于 2020-9-14 21:01
compute VACF  all vacf
fix VACF all ave/time 1 12 ${ndump} c_VACF[1] c_VACF[2] c_VACF[3] c_VACF[4]  ...

你取的时间太长了,文献2ps相关性就没有了。
在就是,你的帧数够多吗?

24

帖子

0

威望

117

eV
积分
141

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2020-9-15 15:28:05 | 只看该作者 Only view this author
dummy1 发表于 2020-9-15 12:15
你取的时间太长了,文献2ps相关性就没有了。
在就是,你的帧数够多吗?

好像确实不够,文献是0.005ps一步,我取得有点大了,另外是不是还要进行归一化处理?

33

帖子

0

威望

1980

eV
积分
2013

Level 5 (御坂)

5#
发表于 Post on 2020-9-16 01:43:20 | 只看该作者 Only view this author
txh 发表于 2020-9-15 15:28
好像确实不够,文献是0.005ps一步,我取得有点大了,另外是不是还要进行归一化处理?

一般需要归一化, 但是应该不影响你的线形

24

帖子

0

威望

117

eV
积分
141

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2020-9-16 08:40:15 | 只看该作者 Only view this author
dummy1 发表于 2020-9-16 01:43
一般需要归一化, 但是应该不影响你的线形

我按照文献,0.005 ps取了一个值,又算了一遍但是还是算不出文献那种图形,虽然体系不一样,但是我看VACF文献中都是那种形式的曲线,我这个为什么还是不对啊,还请老师解答一下。

QQ图片20200916083755.png (10.12 KB, 下载次数 Times of downloads: 24)

QQ图片20200916083755.png

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 18:00 , Processed in 0.168798 second(s), 23 queries , Gzip On.

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