lammps計算速度自相關函數(shù)
我用lammps的compute vacf 計算vacf然后輸出得到的圖如下,但是跟文獻并不一樣,是我沒有進行歸一化處理還是需要在加別的命令進行計算啊,請求各位老師指點一下,下圖左邊是文獻,右邊是自己畫的圖
返回小木蟲查看更多
今日熱帖
我用lammps的compute vacf 計算vacf然后輸出得到的圖如下,但是跟文獻并不一樣,是我沒有進行歸一化處理還是需要在加別的命令進行計算啊,請求各位老師指點一下,下圖左邊是文獻,右邊是自己畫的圖
返回小木蟲查看更多
京公網安備 11010802022153號
我不知道你處理的數(shù)據來源。但是,無論如何速度自相關函數(shù)vacf(t)隨時間t的變化需要滿足如下條件:
(1)先假設你這里的速度自相關函數(shù)vacf(t)是指兩個不同時間的速度向量標積(不是其他乘積)的系綜平均值。
(2)因為處理的都是隨機行為(無論理論模型還是實驗測量),漲落是無法避免的。越是粒子質量重的體系,運動時間尺度越長。
平動時間尺度約在10^(-15)即fs量級,那么右圖漲落很明顯,更像vacf(t)圖。而左圖的漲落很不明顯,看起來時間尺度很長,近于或大于ps量級,那是什么運動呢?不像轉動,更像頻率很低的振動模;蚴窍窈苤氐牧W訄F之間的相對運動?
(3)無論是復雜體系還是簡單體系,總有
t=0時,vacf(t=0)=1;而且1是自相關的極大值,t不等于0時,不可能大于1。
那么,
左圖符合此條件。
右圖那個t=0的點應該遠高于縱坐標上端,不在視線內。
先作此討論,供參考。
補充:我上面的答復都是討論在已經做了歸一化處理后的行為,不然,就一定不會vacf(t=0)=1。
你好,這是已經做了歸一化的圖形,左圖和右圖的體系是不一樣的,我只是想做出她那種形式的圖形,但是我感覺我這個圖的表現(xiàn)形式跟作圖不太一樣
QQͼƬ20200916083755.png
,