| 5 | 1/1 | 返回列表 |
| 查看: 2730 | 回復(fù): 10 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||||
astringent銅蟲 (著名寫手)
|
[交流]
【求助】均方位移 已有4人參與
|
|||
|
最近準(zhǔn)備寫均方位移(MSD)的程序,但是有一點(diǎn)不太明白。均方位移指的是分子質(zhì)心的均方位移,還是采用動(dòng)力學(xué)模擬后得到的軌跡直接算就可以?如果是指質(zhì)心的均方位移,,那么在計(jì)算時(shí)是不是首先要考慮質(zhì)心的移動(dòng)?另外我看文獻(xiàn)上說 unfold each molecule,這個(gè)unfold又是什么意思呢? 請(qǐng)大家不吝賜教,先謝謝了。 |
xuexijisuan |
銅蟲 (著名寫手)
|
不好意思,我是新手,dcd文件是二進(jìn)制,無法讀。你能把如何看軌跡文件里任何時(shí)刻的每一個(gè)粒子坐標(biāo)是不是都沒有超過盒子邊界具體如何操作說的再具體些嗎?另外,如果方便的話,我想看看你的msd程序,我的郵箱是yanbochen123@gmail.com 多謝了。 |

|
首先,看你用的程序輸出的是真實(shí)坐標(biāo)還是應(yīng)用周期化邊界以后的坐標(biāo),如果是前者那很好辦,直接求分子質(zhì)心的均方位移就行。如果是后者,就沒轍,你只能取分析的這段時(shí)間區(qū)間內(nèi)體系中運(yùn)動(dòng)沒有超過盒子邊界的分子來做msd。我自己寫的程序就是為了算msd就把輸出坐標(biāo)改成真實(shí)坐標(biāo)了。 至于unfold each molecule,沒有上下文,我猜是將那些由于應(yīng)用了周期性邊界的而發(fā)生了坐標(biāo)分裂分子(比如說一個(gè)分子的一部分坐標(biāo)在盒子左邊界,一部分在盒子右邊界)去除周期性邊界條件將他們還原到盒子的一側(cè)來。 |
銅蟲 (著名寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +14 | 李李不服輸 2026-03-26 | 14/700 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué)321分求調(diào)劑 +6 | 三dd. 2026-03-30 | 7/350 |
|
|
[考研] 322求調(diào)劑 +9 | 宋明欣 2026-03-27 | 9/450 |
|
|
[考研] 070300化學(xué)354求調(diào)劑 +7 | 101次希望 2026-03-28 | 7/350 |
|
|
[考研] 085602 化學(xué)工程專碩 340分求調(diào)劑 +4 | qianbai11 2026-03-29 | 4/200 |
|
|
[考研] 296求調(diào)劑 +10 | 彼岸t 2026-03-29 | 10/500 |
|
|
[考研] 318求調(diào)劑 +7 | 陳晨79 2026-03-30 | 7/350 |
|
|
[考研] 285求調(diào)劑 +5 | AZMK 2026-03-29 | 8/400 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-26 | 7/350 |
|
|
[考研] 一志愿北京工業(yè)大學(xué),324分求調(diào)劑 +6 | 零八# 2026-03-28 | 6/300 |
|
|
[考博] 26申博自薦 +6 | whh869393 2026-03-24 | 6/300 |
|
|
[考研] 298求調(diào)劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
[考研] 330分求調(diào)劑 +5 | qzenlc 2026-03-29 | 5/250 |
|
|
[考研] 295材料工程專碩求調(diào)劑 +7 | 1428151015 2026-03-27 | 7/350 |
|
|
[考研] 283求調(diào)劑 +3 | A child 2026-03-28 | 3/150 |
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分,英語已過四六級(jí),六級(jí)521,求調(diào)劑 +4 | 美味可樂雞翅 2026-03-26 | 4/200 |
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
[考研] 化工專碩求調(diào)劑 +3 | question挽風(fēng) 2026-03-24 | 3/150 |
|