 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4279 | 回復(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ù)此樓» 猜你喜歡
0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有4人回復(fù)
-
0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國(guó)級(jí)證書40+數(shù)一英一287
已經(jīng)有5人回復(fù)
-
070300化學(xué)319求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿西南交大，求調(diào)劑
已經(jīng)有3人回復(fù)
-
286求調(diào)劑
已經(jīng)有6人回復(fù)
-
焦慮
已經(jīng)有11人回復(fù)
-
331求調(diào)劑（0703有機(jī)化學(xué)
已經(jīng)有7人回復(fù)
-
298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑
已經(jīng)有5人回復(fù)
-
297求調(diào)劑
已經(jīng)有5人回復(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ù)
- 求助？怎么從劍橋晶體庫(kù)得到物質(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ù)
9樓2019-04-10 15:39:51
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ā)表 | |
|---|---|---|---|---|
|
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +6 | ZY-05 2026-03-13 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +5 | 手機(jī)用戶 2026-03-17 | 5/250 |
|
|
[考博] 環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
|
[考研] 307求調(diào)劑 +3 | 冷笙123 2026-03-17 | 3/150 |
|
|
|
[考研] 265求調(diào)劑 +3 | 梁梁校校 2026-03-17 | 3/150 |
|
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
|
[考研] 328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-16 | 8/400 |
|
|
|
[考研] 求調(diào)劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學(xué)。調(diào)劑到任何專業(yè)都可以 +4 | 小丁想進(jìn)步 2026-03-11 | 5/250 |
|
|
|
[考研] 311求調(diào)劑 +5 | 26研0 2026-03-15 | 5/250 |
|
|
|
[考研] 085600材料與化工 求調(diào)劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
|
[考研] 一志愿華中師范071000，325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
|
[考研] 學(xué)碩285求調(diào)劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
|
[考研] 307求調(diào)劑 +5 | 超級(jí)伊昂大王 2026-03-12 | 5/250 |
|
|
|
[考研] 工科278分求調(diào)劑 +5 | 周慢熱啊 2026-03-12 | 7/350 |
|
|
|
[考研] 314求調(diào)劑 +7 | 無懈可擊的巨人 2026-03-12 | 7/350 |
|
|
|
[考研] 290求調(diào)劑 +3 | ADT 2026-03-13 | 3/150 |
|
|
|
[考研] 283求調(diào)劑，材料、化工皆可 +8 | 蘇打水7777 2026-03-11 | 10/500 |
|