 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 9 | 1/1 | 返回列表 |
| 查看: 4277 | 回復(fù): 8 | ||||
[求助]
從LAMMPS得到的原子速度，如何計算速度自相關(guān)函數(shù)（VAF）？求助 已有3人參與
|
| 比如說我的系統(tǒng)里有4000個原子，在某個狀態(tài)下跑了100步，每一步的時間間隔是確定的（2fs），每一步時每個原子的三維速度都會輸出。我想計算這個系統(tǒng)的速度自相關(guān)函數(shù)VAF，計算公式是<v(0)·v(t)> 。我不太明白這個公式的含義是計算什么，用FORTRAN或MATLAB應(yīng)該如何編程實現(xiàn)呢？緊急求助，不勝感激。 |
回復(fù)此樓» 猜你喜歡
312求調(diào)劑
已經(jīng)有7人回復(fù)
-
化學(xué)工程321分求調(diào)劑
已經(jīng)有16人回復(fù)
-
280求調(diào)劑
已經(jīng)有3人回復(fù)
-
266求調(diào)劑
已經(jīng)有7人回復(fù)
-
工科材料085601 279求調(diào)劑
已經(jīng)有5人回復(fù)
-
環(huán)境工程調(diào)劑
已經(jīng)有8人回復(fù)
-
302求調(diào)劑
已經(jīng)有9人回復(fù)
-
297求調(diào)劑
已經(jīng)有4人回復(fù)
-
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑
已經(jīng)有4人回復(fù)
-
焦慮
已經(jīng)有10人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- lammps中計算彈性常數(shù)的例子是基于什么原理？
已經(jīng)有1人回復(fù)
- 從lammps模擬使用的勢函數(shù)畫出勢能曲線
已經(jīng)有0人回復(fù)
- 急求解答：lammps安裝問題
已經(jīng)有10人回復(fù)
- 對原子數(shù)過多的晶格進行結(jié)構(gòu)優(yōu)化
已經(jīng)有3人回復(fù)
- 【轉(zhuǎn)】（已搜索無重復(fù)）LAMMPS中fix ave spatial空間平均命令中文版
已經(jīng)有26人回復(fù)
- 求助？怎么從劍橋晶體庫得到物質(zhì)的結(jié)構(gòu)數(shù)據(jù)
已經(jīng)有5人回復(fù)
- lammps 時間關(guān)聯(lián)函數(shù)求助
已經(jīng)有1人回復(fù)
- LAMMPS溫度－步長曲線震蕩
已經(jīng)有11人回復(fù)
- 有用lammps模擬過聚合物的相關(guān)性能嗎？
已經(jīng)有12人回復(fù)
- python之蟒開啟理論計算化學(xué)的新時代
已經(jīng)有22人回復(fù)
- 關(guān)于towhee 回答zhang_jaj的問題
已經(jīng)有44人回復(fù)
- 怎樣用lammps得到無定形Si的結(jié)構(gòu)
已經(jīng)有3人回復(fù)
- 站里關(guān)于MS Analysis計算材料的應(yīng)力應(yīng)變與Cubic以及Trigonal cell彈性常數(shù)的討論
已經(jīng)有34人回復(fù)
chuntao118
新蟲 (著名寫手)
- 應(yīng)助: 93 (初中生)
- 金幣: 2590.9
- 散金: 1402
- 紅花: 25
- 帖子: 1611
- 在線: 241.2小時
- 蟲號: 1791567
- 注冊: 2012-05-02
- 專業(yè): 計算數(shù)學(xué)與科學(xué)工程計算
2樓2016-03-22 14:32:44
challenger9
銀蟲 (初入文壇)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 278.1
- 紅花: 1
- 帖子: 50
- 在線: 55.1小時
- 蟲號: 2382504
- 注冊: 2013-03-27
- 專業(yè): 電磁場與波
3樓2016-06-24 10:06:08
brucefan
專家顧問 (著名寫手)
-
專家經(jīng)驗: +218 - 模擬EPI: 4
- 應(yīng)助: 422 (碩士)
- 金幣: 2294.8
- 散金: 1720
- 紅花: 88
- 帖子: 1423
- 在線: 670.2小時
- 蟲號: 3048965
- 注冊: 2014-03-14
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
- 管轄: 分子模擬
【答案】應(yīng)助回帖
|
% a matlab script for calculating VACF from velocity data clear; load v.txt; % assume your velocity data are in the above file and assume the format is (N is the number of atoms, M is the number of time points): % vx_1 vy_1 vz_1 % time point 1 % vx_2 vy_2 vz_2 % time point 1 % ... % time point 1 % vx_N, vy_N, vz_N % time point 1 % vx_1 vy_1 vz_1 % time point 2 % vx_2 vy_2 vz_2 % time point 2 % ... % time point 2 % vx_N, vy_N, vz_N % time point 2 % ... % vx_1 vy_1 vz_1 % time point M % vx_2 vy_2 vz_2 % time point M % ... % time point M % vx_N, vy_N, vz_N % time point M N = xxx; % number of atoms in your system M = length(v)/N; % number of time points for your velocity data dt = xxx; % the time interval between two set of velocities (in some unit) Nt = xxx; % maximum length of the correlation your want (usually Nt = M/10 is a good choice) time = dt*(0:Nt-1); M = M-Nt; % you have to waste a small portion of data vacf=zeros(Nt,1); for nt=0:Nt-1 for m=1:M vacf(nt+1,  =sum(sum(v_all((m-1)*N+1:m*N,.*v_all((m+nt-1)*N+1 m+nt)*N,));end end %vacf=vacf/M; % you can also normalize it by using [vacf=vacf(1);] if you want % now you can plot the result: close all; figure; plot(time, vacf,'o-'); xlabel('time (some unit)'); ylabel('VACF (some unit)'); |
4樓2016-06-27 18:45:07
brucefan
專家顧問 (著名寫手)
-
專家經(jīng)驗: +218 - 模擬EPI: 4
- 應(yīng)助: 422 (碩士)
- 金幣: 2294.8
- 散金: 1720
- 紅花: 88
- 帖子: 1423
- 在線: 670.2小時
- 蟲號: 3048965
- 注冊: 2014-03-14
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
- 管轄: 分子模擬
5樓2016-06-27 18:52:22
cuso4852cuoh
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 223.4
- 散金: 20
- 帖子: 170
- 在線: 49.6小時
- 蟲號: 1481571
- 注冊: 2011-11-08
- 專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域
6樓2018-11-26 19:15:36
looyady
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1238.3
- 帖子: 120
- 在線: 15.8小時
- 蟲號: 2707548
- 注冊: 2013-10-08
- 專業(yè): 工程熱物理與能源利用
7樓2018-12-14 15:21:29
8樓2018-12-19 16:27:11
9樓2019-04-10 15:39:51
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 安安安子 的主題更新
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 304求調(diào)劑 +9 | 小熊joy 2026-03-14 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +3 | 脫穎而出 2026-03-16 | 6/300 |
|
|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-16 | 8/400 |
|
|
|
[考研] 085600材料與化工求調(diào)劑 +5 | 緒幸與子 2026-03-17 | 5/250 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研]
 11408 一志愿西電，277分求調(diào)劑
+3
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
|
[考研] 286求調(diào)劑 +3 | lemonzzn 2026-03-16 | 5/250 |
|
|
|
[考研] 0703一志愿211 285分求調(diào)劑 +5 | ly3471z 2026-03-13 | 5/250 |
|
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對?/a> 2026-03-13 | 5/250 |
|
|
|
[考研] 0703，333分求調(diào)劑 一志愿鄭州大學(xué)-物理化學(xué) +3 | 李魔女斗篷 2026-03-11 | 3/150 |
|
|
|
[考研] 329求調(diào)劑 +3 | miaodesi 2026-03-12 | 4/200 |
|
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
|
[考研] 314求調(diào)劑 +7 | 無懈可擊的巨人 2026-03-12 | 7/350 |
|
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級已過 本科山東雙非 求調(diào)劑！ +3 | 不想理你 2026-03-12 | 3/150 |
|
|
|
[論文投稿]
投稿問題
 5+4
|
星光燦爛xt 2026-03-12 | 6/300 |
|
|
|
[考研] 277求調(diào)劑 +4 | anchor17 2026-03-12 | 4/200 |
|
|
|
[考研] 化工學(xué)碩306求調(diào)劑 +9 | 42838695 2026-03-12 | 9/450 |
|
|
|
[考博] 2026年博士申請 +3 | QwQwQW10 2026-03-11 | 3/150 |
|