24小时热门版块排行榜    

查看: 2249  |  回复: 3
【悬赏金币】回答本帖问题,作者志远T将赠送您 25 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

志远T

新虫 (初入文坛)

[求助] lammps计算速度自相关函数 已有1人参与

我用lammps的compute vacf 计算vacf然后输出得到的图如下,但是跟文献并不一样,是我没有进行归一化处理还是需要在加别的命令进行计算啊,请求各位老师指点一下,下图左边是文献,右边是自己画的图

lammps计算速度自相关函数


@我想更懂你29 @甄士隐207 @jiagle @jiagle 发自小木虫IOS客户端
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

志远T

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by mbchen at 2020-09-22 17:50:00
补充:我上面的答复都是讨论在已经做了归一化处理后的行为,不然,就一定不会vacf(t=0)=1。

你好,这是已经做了归一化的图形,左图和右图的体系是不一样的,我只是想做出她那种形式的图形,但是我感觉我这个图的表现形式跟作图不太一样
lammps计算速度自相关函数-1
QQͼƬ20200916083755.png

4楼2020-10-13 09:01:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

mbchen

专家顾问 (著名写手)

【答案】应助回帖

我不知道你处理的数据来源。但是,无论如何速度自相关函数vacf(t)随时间t的变化需要满足如下条件:

(1)先假设你这里的速度自相关函数vacf(t)是指两个不同时间的速度向量标积(不是其他乘积)的系综平均值。

(2)因为处理的都是随机行为(无论理论模型还是实验测量),涨落是无法避免的。越是粒子质量重的体系,运动时间尺度越长。
平动时间尺度约在10^(-15)即fs量级,那么右图涨落很明显,更像vacf(t)图。而左图的涨落很不明显,看起来时间尺度很长,近于或大于ps量级,那是什么运动呢?不像转动,更像频率很低的振动模。或是像很重的粒子团之间的相对运动?

(3)无论是复杂体系还是简单体系,总有
t=0时,vacf(t=0)=1;而且1是自相关的极大值,t不等于0时,不可能大于1。
那么,
左图符合此条件。
右图那个t=0的点应该远高于纵坐标上端,不在视线内。

先作此讨论,供参考。
2楼2020-09-22 17:45:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mbchen

专家顾问 (著名写手)

【答案】应助回帖

补充:我上面的答复都是讨论在已经做了归一化处理后的行为,不然,就一定不会vacf(t=0)=1。
3楼2020-09-22 17:50:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿天津大学材料与化工275求调剂 10+4 穿只靴子 2026-03-07 18/900 2026-03-10 18:44 by chixmc
[考研] 求调剂材料专硕293 +6 段_(:з」∠)_ 2026-03-10 6/300 2026-03-10 18:22 by ms629
[考研] 0703求调剂 +3 jtyq001 2026-03-10 3/150 2026-03-10 16:45 by ztnimte
[考研] 材料工程专硕调剂自荐信,初试323 +6 李白26 2026-03-07 6/300 2026-03-10 16:16 by 无际的草原
[考研] 083000环境科学与工程调剂 +6 mingmingry 2026-03-09 7/350 2026-03-10 14:37 by kkkdtdde
[考研] 材料专硕调剂 +7 慕辰123 2026-03-05 10/500 2026-03-10 14:22 by peike
[考研] 307求调剂 +5 超级伊昂大王 2026-03-10 5/250 2026-03-10 14:01 by xxrvvv
[考研] 085600材料与化工 326 求调剂 +4 热爱生活ing 2026-03-09 4/200 2026-03-10 08:25 by Linda Hu
[考研] 083000环境科学与工程调剂 +5 加油呀fxy 2026-03-07 6/300 2026-03-08 21:07 by sdgdsdgd
[考研] 0701-322 求调剂 +3 jiliuxian 2026-03-06 8/400 2026-03-08 19:31 by jiliuxian
[教师之家] 交大前校长王树国:现在最先进的科技并不在大学实验室,而是在企业研究院 +4 zju2000 2026-03-08 6/300 2026-03-08 19:15 by zju2000
[考研] 085701环境工程317分求调剂 +9 6汆尼9 2026-03-07 9/450 2026-03-08 06:41 by 刘兵
[考研] 材料科学(0805)338 求调剂 +7 xiaokang3286 2026-03-07 7/350 2026-03-08 03:42 by lfhuang
[考研] 276求调剂 本科天津工业大学 一志愿北京化工大学0855 +4 路lyh123 2026-03-06 6/300 2026-03-07 11:37 by shdliugang
[考研] 308求调剂 +7 倘若起风了呢 2026-03-05 9/450 2026-03-06 22:59 by 倘若起风了呢
[考研] 材料与化工304求调剂 +7 邱gl 2026-03-05 10/500 2026-03-06 01:11 by 卷柏卷柏
[考研] 材料调剂 +5 ws 上岸鸭 2026-03-05 5/250 2026-03-05 20:06 by houyaoxu
[考研] 纺织、生物、化学、材料等专业 +3 Eember. 2026-03-05 7/350 2026-03-05 16:35 by Eember.
[考研] 322分 085600求调剂,有互联网+国金及主持省级大创经历 +6 熊境喆 2026-03-04 6/300 2026-03-04 20:32 by kakakapanpan
[考研] 材料专硕346求调剂 +3 旺一下 2026-03-04 3/150 2026-03-04 16:26 by sslc1985
信息提示
请填处理意见