24小时热门版块排行榜    

查看: 2461  |  回复: 3
【悬赏金币】回答本帖问题,作者志远T将赠送您 25 个金币

志远T

新虫 (初入文坛)

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

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

lammps计算速度自相关函数


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

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的回帖

志远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的回帖
相关版块跳转 我要订阅楼主 志远T 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 10/500 2026-09-02 11:35 by 大不刘6
[基金申请] 我不理解! +18 Edward_pc 2026-08-26 26/1300 2026-09-02 11:01 by shanzhipu
[基金申请] 科研人应该花精力去思考如何解决问题,而不是去凝练问题 +4 瞬息宇宙 2026-09-01 8/400 2026-09-02 10:35 by 瞬息宇宙
[基金申请] 学科评审组评审是指会评吗? +5 瞬息宇宙 2026-08-31 5/250 2026-09-02 10:13 by 雨冰共舞
[文学芳草园] 梦想 +5 myrtle 2026-08-26 8/400 2026-09-02 02:42 by Leogzhya
[基金申请] 面上合作单位盖章 +6 ssyjh 2026-08-27 9/450 2026-09-01 20:01 by huagongfeihu
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
信息提示
请填处理意见