 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4318 | 回復(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ù)此樓» 猜你喜歡
【雙一流院校新能源、環(huán)境材料，材料加工與模擬招收大量調(diào)劑】
已經(jīng)有7人回復(fù)
-
08工學(xué)調(diào)劑
已經(jīng)有12人回復(fù)
-
材料與化工328分調(diào)劑
已經(jīng)有4人回復(fù)
-
279分求調(diào)劑 一志愿211
已經(jīng)有20人回復(fù)
-
293求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿華東理工大學(xué)081700，初試分?jǐn)?shù)271
已經(jīng)有6人回復(fù)
-
085404電子信息284分求調(diào)劑
已經(jīng)有3人回復(fù)
-
306求0703調(diào)劑一志愿華中師范
已經(jīng)有9人回復(fù)
-
一志愿吉大化學(xué)322求調(diào)劑
已經(jīng)有6人回復(fù)
-
求調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- lammps中計(jì)算彈性常數(shù)的例子是基于什么原理？
已經(jīng)有1人回復(fù)
- 從lammps模擬使用的勢函數(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溫度－步長曲線震蕩
已經(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ù)
8樓2018-12-19 16:27:11
chuntao118
新蟲 (著名寫手)
- 應(yīng)助: 93 (初中生)
- 金幣: 2590.9
- 散金: 1402
- 紅花: 25
- 帖子: 1611
- 在線: 241.2小時(shí)
- 蟲號(hào): 1791567
- 注冊: 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
- 注冊: 2013-03-27
- 專業(yè): 電磁場與波
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
- 注冊: 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ā)表 | |
|---|---|---|---|---|
|
[考研] 335求調(diào)劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿國科過程所081700，274求調(diào)劑 +3 | 三水研0水立方 2026-03-23 | 3/150 |
|
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|
|
|
[考研]
 材料與化工考研調(diào)劑
+4
|
孅華 2026-03-22 | 4/200 |
|
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +8 | tcxiaoxx 2026-03-20 | 9/450 |
|
|
|
[考研] 298求調(diào)劑 +8 | 上岸6666@ 2026-03-20 | 8/400 |
|
|
|
[考研] 352求調(diào)劑 +3 | 大米飯！ 2026-03-22 | 3/150 |
|
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-18 | 6/300 |
|
|
|
[考研] 307求調(diào)劑 +11 | 冷笙123 2026-03-17 | 11/550 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3 | 晨昏線與星海 2026-03-20 | 3/150 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-20 | 7/350 |
|
|
|
[考研] 311求調(diào)劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
|
[考研] 332求調(diào)劑 +4 | ydfyh 2026-03-17 | 4/200 |
|
|
|
[考研] 324分 085600材料化工求調(diào)劑 +4 | llllkkkhh 2026-03-18 | 4/200 |
|
|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
|
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|