| 5 | 1/1 | 返回列表 |
| 查看: 6959 | 回復(fù): 30 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||||
敏sunflower金蟲 (小有名氣)
|
[求助]
計(jì)算擴(kuò)散系數(shù) 已有2人參與
|
|||
|
初學(xué)模擬,急需各位大俠的幫助! 我最近在模擬A、B兩分子在聚合物中的擴(kuò)散系數(shù),但是遇到了幾個(gè)問題,想向各位請(qǐng)教一下 模擬過程如下:1.建Cell——2.能量最小化——3.NPT——4.NVT平衡結(jié)構(gòu)——5.NVT用于數(shù)據(jù)分析 (1)做完第5步NVT動(dòng)力學(xué),分析得到MSD-t曲線,從而得到擴(kuò)散系數(shù)。但是發(fā)現(xiàn)對(duì)第4步結(jié)束時(shí)的最后一幀結(jié)構(gòu)平行做兩次NVT動(dòng)力學(xué),分析得到的擴(kuò)散系數(shù)有時(shí)相差很大,這是什么原因那? (2)有人說是擴(kuò)散系數(shù)對(duì)初始結(jié)構(gòu)很敏感,但重新建模之后進(jìn)行上述1-5步模擬,發(fā)現(xiàn)得到的擴(kuò)散系數(shù)也相差很大 面對(duì)上述情況,怎樣取擴(kuò)散系數(shù)的值那? ![]() |
擴(kuò)散系數(shù) |
|
本帖內(nèi)容被屏蔽 |
|
本帖內(nèi)容被屏蔽 |
金蟲 (小有名氣)
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗(yàn): +49 |
|
聚合物中如果有孔洞的話,嚴(yán)格意義上來講粒子在孔洞中的擴(kuò)散是受限的。所以會(huì)出現(xiàn)粒子長(zhǎng)時(shí)間只在某個(gè)區(qū)域的孔洞中運(yùn)動(dòng)的情況,這樣的話如果你從不同的初始條件開始而采樣時(shí)間不夠長(zhǎng)的話,采集到的MSD只是粒子在某個(gè)區(qū)域里的MSD,擬合得到的擴(kuò)散系數(shù)自然會(huì)相差很大。即便你從相同的初始條件出發(fā),由于MD算法自身的誤差以及NVT動(dòng)力學(xué)的一些特點(diǎn)(這個(gè)有些復(fù)雜,這里不講),也不可能得到兩條相同的軌跡或是MSD。 一般對(duì)于粒子在這種受限條件下的擴(kuò)散,采樣時(shí)間都在納秒級(jí)別以保證充分采樣。雖然你用300ps就可以達(dá)到平衡,但采樣同樣是要花很長(zhǎng)時(shí)間的。如果你不想跑長(zhǎng)的simulation,可以考慮以不同的初始條件跑幾個(gè)600ps的simulation,然后把求到的幾條MSD取平均。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +9 | 辰沐5211314 2026-03-26 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 各位老師好,我的一志愿為北京科技大學(xué)085601材料專碩 +7 | Koxui 2026-03-28 | 7/350 |
|
|
[考研] 305求調(diào)劑 +7 | RuiFairyrui 2026-03-28 | 7/350 |
|
|
[考研] 332求調(diào)劑 +6 | 蕉蕉123 2026-03-28 | 6/300 |
|
|
[考研] 085600,專業(yè)課化工原理,321分求調(diào)劑 +4 | 大饞小子 2026-03-28 | 4/200 |
|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑 +7 | 芋泥小鈴鐺 2026-03-28 | 7/350 |
|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
[考研] 本科新能源科學(xué)與工程,一志愿華理能動(dòng)285求調(diào)劑 +7 | AZMK 2026-03-28 | 11/550 |
|
|
[考研] 321求調(diào)劑 +6 | 璞玉~~ 2026-03-25 | 7/350 |
|
|
[考研] 275求調(diào)劑 +10 | Micky11223 2026-03-25 | 14/700 |
|
|
[考研] 081200-314 +3 | LILIQQ 2026-03-27 | 4/200 |
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|
|
[考研] 324求調(diào)劑 +5 | hanamiko 2026-03-26 | 5/250 |
|
|
[考研] 一志愿北化求調(diào)劑 +3 | Jsman 2026-03-22 | 3/150 |
|
|
[考研] 327求調(diào)劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
[考研]
材料調(diào)劑
5+4
|
想要一壺桃花水 2026-03-25 | 10/500 |
|
|
[考研]
|
平樂樂樂 2026-03-26 | 4/200 |
|
|
[考研]
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
[考研] 333求調(diào)劑 +3 | ALULU4408 2026-03-23 | 3/150 |
|