 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4315 | 回復(fù): 8 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||
[求助]
從LAMMPS得到的原子速度，如何計(jì)算速度自相關(guān)函數(shù)（VAF）？求助 已有3人參與
|
||||
| 比如說我的系統(tǒng)里有4000個(gè)原子，在某個(gè)狀態(tài)下跑了100步，每一步的時(shí)間間隔是確定的（2fs），每一步時(shí)每個(gè)原子的三維速度都會(huì)輸出。我想計(jì)算這個(gè)系統(tǒng)的速度自相關(guān)函數(shù)VAF，計(jì)算公式是<v(0)·v(t)> 。我不太明白這個(gè)公式的含義是計(jì)算什么，用FORTRAN或MATLAB應(yīng)該如何編程實(shí)現(xiàn)呢？緊急求助，不勝感激。 |
回復(fù)此樓» 猜你喜歡
323求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑
已經(jīng)有4人回復(fù)
-
352求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿東華大學(xué)化學(xué)070300，求調(diào)劑
已經(jīng)有8人回復(fù)
-
277材料科學(xué)與工程080500求調(diào)劑
已經(jīng)有7人回復(fù)
-
317求調(diào)劑
已經(jīng)有18人回復(fù)
-
293求調(diào)劑
已經(jīng)有5人回復(fù)
-
280分求調(diào)劑 一志愿085802
已經(jīng)有7人回復(fù)
-
0854電子信息求調(diào)劑
已經(jīng)有3人回復(fù)
-
263求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- lammps中計(jì)算彈性常數(shù)的例子是基于什么原理？
已經(jīng)有1人回復(fù)
- 從lammps模擬使用的勢(shì)函數(shù)畫出勢(shì)能曲線
已經(jīng)有0人回復(fù)
- 急求解答：lammps安裝問題
已經(jīng)有10人回復(fù)
- 對(duì)原子數(shù)過多的晶格進(jìn)行結(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 時(shí)間關(guān)聯(lián)函數(shù)求助
已經(jīng)有1人回復(fù)
- LAMMPS溫度－步長(zhǎng)曲線震蕩
已經(jīng)有11人回復(fù)
- 有用lammps模擬過聚合物的相關(guān)性能嗎？
已經(jīng)有12人回復(fù)
- python之蟒開啟理論計(jì)算化學(xué)的新時(shí)代
已經(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計(jì)算材料的應(yīng)力應(yīng)變與Cubic以及Trigonal cell彈性常數(shù)的討論
已經(jīng)有34人回復(fù)
brucefan
專家顧問 (著名寫手)
-
專家經(jīng)驗(yàn): +218 - 模擬EPI: 4
- 應(yīng)助: 422 (碩士)
- 金幣: 2294.8
- 散金: 1720
- 紅花: 88
- 帖子: 1423
- 在線: 670.2小時(shí)
- 蟲號(hào): 3048965
- 注冊(cè): 2014-03-14
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
- 管轄: 分子模擬
5樓2016-06-27 18:52:22
chuntao118
新蟲 (著名寫手)
- 應(yīng)助: 93 (初中生)
- 金幣: 2590.9
- 散金: 1402
- 紅花: 25
- 帖子: 1611
- 在線: 241.2小時(shí)
- 蟲號(hào): 1791567
- 注冊(cè): 2012-05-02
- 專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算
2樓2016-03-22 14:32:44
challenger9
銀蟲 (初入文壇)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 278.1
- 紅花: 1
- 帖子: 50
- 在線: 55.1小時(shí)
- 蟲號(hào): 2382504
- 注冊(cè): 2013-03-27
- 專業(yè): 電磁場(chǎng)與波
3樓2016-06-24 10:06:08
brucefan
專家顧問 (著名寫手)
-
專家經(jīng)驗(yàn): +218 - 模擬EPI: 4
- 應(yīng)助: 422 (碩士)
- 金幣: 2294.8
- 散金: 1720
- 紅花: 88
- 帖子: 1423
- 在線: 670.2小時(shí)
- 蟲號(hào): 3048965
- 注冊(cè): 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
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿武理材料工程348求調(diào)劑 +5 | ￣^￣゜汗 2026-03-19 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑 +3 | jiajunX 2026-03-22 | 3/150 |
|
|
|
[考研] 324求調(diào)劑 +6 | lucky呀呀呀鴨 2026-03-20 | 6/300 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3 | 晨昏線與星海 2026-03-20 | 3/150 |
|
|
|
[考研] 305分求調(diào)劑（食品工程） +4 | Sxy112 2026-03-21 | 6/300 |
|
|
|
[考研] 生物學(xué)一志愿985，分?jǐn)?shù)349求調(diào)劑 +4 | zxts12 2026-03-21 | 7/350 |
|
|
|
[考研] 化學(xué)調(diào)劑 +5 | yzysaa 2026-03-21 | 5/250 |
|
|
|
[考研] 材料求調(diào)劑 +5 | @taotao 2026-03-21 | 5/250 |
|
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級(jí)已過 本科山東雙非 求調(diào)劑！ +3 | 不想理你 2026-03-16 | 3/150 |
|
|
|
[考研] 一志愿中國石油大學(xué)（華東） 本科齊魯工業(yè)大學(xué) +3 | 石能偉 2026-03-17 | 3/150 |
|
|
|
[考研]
 271材料工程求調(diào)劑
+8
|
.6lL 2026-03-18 | 8/400 |
|
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
|
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
|
[考研]
 295復(fù)試調(diào)劑
+8
|
簡(jiǎn)木ChuFront 2026-03-19 | 8/400 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +7 | heming3743 2026-03-16 | 7/350 |
|
|
|
[論文投稿]
申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒變，給編輯又寫郵件了，沒回復(fù)
 10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|