我用lammps的compute vacf 计算vacf然后输出得到的图如下,但是跟文献并不一样,是我没有进行归一化处理还是需要在加别的命令进行计算啊,请求各位老师指点一下,下图左边是文献,右边是自己画的图 返回小木虫查看更多
我不知道你处理的数据来源。但是,无论如何速度自相关函数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的点应该远高于纵坐标上端,不在视线内。 先作此讨论,供参考。
补充:我上面的答复都是讨论在已经做了归一化处理后的行为,不然,就一定不会vacf(t=0)=1。
我不知道你处理的数据来源。但是,无论如何速度自相关函数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的点应该远高于纵坐标上端,不在视线内。
先作此讨论,供参考。
补充:我上面的答复都是讨论在已经做了归一化处理后的行为,不然,就一定不会vacf(t=0)=1。
你好,这是已经做了归一化的图形,左图和右图的体系是不一样的,我只是想做出她那种形式的图形,但是我感觉我这个图的表现形式跟作图不太一样
QQͼƬ20200916083755.png
,