 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4284 | 回復: 8 | ||||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||||
[求助]
從LAMMPS得到的原子速度，如何計算速度自相關函數(shù)（VAF）？求助 已有3人參與
|
||||
| 比如說我的系統(tǒng)里有4000個原子，在某個狀態(tài)下跑了100步，每一步的時間間隔是確定的（2fs），每一步時每個原子的三維速度都會輸出。我想計算這個系統(tǒng)的速度自相關函數(shù)VAF，計算公式是<v(0)·v(t)> 。我不太明白這個公式的含義是計算什么，用FORTRAN或MATLAB應該如何編程實現(xiàn)呢？緊急求助，不勝感激。 |
回復此樓» 猜你喜歡
271材料工程求調劑
已經(jīng)有5人回復
-
281求調劑（0805）
已經(jīng)有16人回復
-
304求調劑
已經(jīng)有6人回復
-
材料工程專碩調劑
已經(jīng)有6人回復
-
一志愿天大材料與化工（085600）總分338
已經(jīng)有4人回復
-
085700資源與環(huán)境308求調劑
已經(jīng)有3人回復
-
求材料調劑
已經(jīng)有8人回復
-
294求調劑材料與化工專碩
已經(jīng)有5人回復
-
一志愿華中科技大學，080502，354分求調劑
已經(jīng)有4人回復
-
一志愿吉林大學材料學碩321求調劑
已經(jīng)有6人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- lammps中計算彈性常數(shù)的例子是基于什么原理？
已經(jīng)有1人回復
- 從lammps模擬使用的勢函數(shù)畫出勢能曲線
已經(jīng)有0人回復
- 急求解答：lammps安裝問題
已經(jīng)有10人回復
- 對原子數(shù)過多的晶格進行結構優(yōu)化
已經(jīng)有3人回復
- 【轉】（已搜索無重復）LAMMPS中fix ave spatial空間平均命令中文版
已經(jīng)有26人回復
- 求助？怎么從劍橋晶體庫得到物質的結構數(shù)據(jù)
已經(jīng)有5人回復
- lammps 時間關聯(lián)函數(shù)求助
已經(jīng)有1人回復
- LAMMPS溫度－步長曲線震蕩
已經(jīng)有11人回復
- 有用lammps模擬過聚合物的相關性能嗎？
已經(jīng)有12人回復
- python之蟒開啟理論計算化學的新時代
已經(jīng)有22人回復
- 關于towhee 回答zhang_jaj的問題
已經(jīng)有44人回復
- 怎樣用lammps得到無定形Si的結構
已經(jīng)有3人回復
- 站里關于MS Analysis計算材料的應力應變與Cubic以及Trigonal cell彈性常數(shù)的討論
已經(jīng)有34人回復
challenger9
銀蟲 (初入文壇)
- 應助: 1 (幼兒園)
- 金幣: 278.1
- 紅花: 1
- 帖子: 50
- 在線: 55.1小時
- 蟲號: 2382504
- 注冊: 2013-03-27
- 專業(yè): 電磁場與波
3樓2016-06-24 10:06:08
chuntao118
新蟲 (著名寫手)
- 應助: 93 (初中生)
- 金幣: 2590.9
- 散金: 1402
- 紅花: 25
- 帖子: 1611
- 在線: 241.2小時
- 蟲號: 1791567
- 注冊: 2012-05-02
- 專業(yè): 計算數(shù)學與科學工程計算
2樓2016-03-22 14:32:44
brucefan
專家顧問 (著名寫手)
-
專家經(jīng)驗: +218 - 模擬EPI: 4
- 應助: 422 (碩士)
- 金幣: 2294.8
- 散金: 1720
- 紅花: 88
- 帖子: 1423
- 在線: 670.2小時
- 蟲號: 3048965
- 注冊: 2014-03-14
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
- 管轄: 分子模擬
【答案】應助回帖
|
% 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
- 應助: 422 (碩士)
- 金幣: 2294.8
- 散金: 1720
- 紅花: 88
- 帖子: 1423
- 在線: 670.2小時
- 蟲號: 3048965
- 注冊: 2014-03-14
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
- 管轄: 分子模擬
5樓2016-06-27 18:52:22
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 330求調劑 +3 | 小材化本科 2026-03-18 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿中海洋材料工程專碩330分求調劑 +5 | 小材化本科 2026-03-18 | 5/250 |
|
|
|
[考研] 354求調劑 +4 | Tyoumou 2026-03-18 | 7/350 |
|
|
|
[考研] 274求調劑 +5 | S.H1 2026-03-18 | 5/250 |
|
|
|
[考研] 267一志愿南京工業(yè)大學0817化工求調劑 +8 | SUICHILD 2026-03-12 | 8/400 |
|
|
|
[考研] 085600材料與化工 +5 | 安全上岸！ 2026-03-16 | 5/250 |
|
|
|
[考研] 314求調劑 +8 | 無懈可擊的巨人 2026-03-12 | 8/400 |
|
|
|
[考研] 331求調劑（0703有機化學 +7 | ZY-05 2026-03-13 | 8/400 |
|
|
|
[考研] 288求調劑，一志愿華南理工大學071005 +4 | ioodiiij 2026-03-17 | 4/200 |
|
|
|
[考研] 278求調劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
|
[考研] 277調劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +8 | 一鴨鴨喲 2026-03-14 | 10/500 |
|
|
|
[考研] 302求調劑 +4 | 小賈同學123 2026-03-15 | 8/400 |
|
|
|
[考研] 一志愿211 0703方向310分求調劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
|
[考研] 326求調劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
|
[考研] 294求調劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
|
[考研] 材料080500調劑求收留 +3 | 一顆meteor 2026-03-13 | 3/150 |
|
|
|
[考研] ［0860］321分求調劑，ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
|
[考研] 329求調劑 +3 | miaodesi 2026-03-12 | 4/200 |
|
|
|
[考研] 308求調劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|