| 5 | 1/1 | 返回列表 |
| 查看: 6958 | 回復(fù): 30 | ||||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||||
敏sunflower金蟲 (小有名氣)
|
[求助]
計算擴散系數(shù) 已有2人參與
|
|||
|
初學模擬,急需各位大俠的幫助! 我最近在模擬A、B兩分子在聚合物中的擴散系數(shù),但是遇到了幾個問題,想向各位請教一下 模擬過程如下:1.建Cell——2.能量最小化——3.NPT——4.NVT平衡結(jié)構(gòu)——5.NVT用于數(shù)據(jù)分析 (1)做完第5步NVT動力學,分析得到MSD-t曲線,從而得到擴散系數(shù)。但是發(fā)現(xiàn)對第4步結(jié)束時的最后一幀結(jié)構(gòu)平行做兩次NVT動力學,分析得到的擴散系數(shù)有時相差很大,這是什么原因那? (2)有人說是擴散系數(shù)對初始結(jié)構(gòu)很敏感,但重新建模之后進行上述1-5步模擬,發(fā)現(xiàn)得到的擴散系數(shù)也相差很大 面對上述情況,怎樣取擴散系數(shù)的值那? ![]() |
擴散系數(shù) |
金蟲 (小有名氣)
|
本帖內(nèi)容被屏蔽 |
|
本帖內(nèi)容被屏蔽 |
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗: +49 |
|
聚合物中如果有孔洞的話,嚴格意義上來講粒子在孔洞中的擴散是受限的。所以會出現(xiàn)粒子長時間只在某個區(qū)域的孔洞中運動的情況,這樣的話如果你從不同的初始條件開始而采樣時間不夠長的話,采集到的MSD只是粒子在某個區(qū)域里的MSD,擬合得到的擴散系數(shù)自然會相差很大。即便你從相同的初始條件出發(fā),由于MD算法自身的誤差以及NVT動力學的一些特點(這個有些復(fù)雜,這里不講),也不可能得到兩條相同的軌跡或是MSD。 一般對于粒子在這種受限條件下的擴散,采樣時間都在納秒級別以保證充分采樣。雖然你用300ps就可以達到平衡,但采樣同樣是要花很長時間的。如果你不想跑長的simulation,可以考慮以不同的初始條件跑幾個600ps的simulation,然后把求到的幾條MSD取平均。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 348求調(diào)劑 +3 | 小懶蟲不懶了 2026-03-28 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 生物學學碩,一志愿湖南大學,初試成績338 +6 | YYYYYNNNNN 2026-03-26 | 7/350 |
|
|
[考研] 求化學調(diào)劑 +4 | wulanna 2026-03-28 | 4/200 |
|
|
[考研] 322求調(diào)劑 +5 | 舊吢 2026-03-24 | 5/250 |
|
|
[考研] 0703一志愿9,初試成績:338,四六級已過,有科研經(jīng)歷,求調(diào)劑! +4 | Zuhui0306 2026-03-25 | 4/200 |
|
|
[考研] 339求調(diào)劑,想調(diào)回江蘇 +6 | 烤麥芽 2026-03-27 | 8/400 |
|
|
[考研] 286求調(diào)劑 +4 | 丟掉懶惰 2026-03-27 | 7/350 |
|
|
[考研] 331環(huán)境科學與工程求調(diào)劑 +3 | 熠然好運氣 2026-03-27 | 3/150 |
|
|
[考研] 322求調(diào)劑 +4 | 我真的很想學習 2026-03-23 | 4/200 |
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學大學 ,080500材料科學與工程學碩 +4 | @taotao 2026-03-26 | 5/250 |
|
|
[考研] 333求調(diào)劑 +6 | wfh030413@ 2026-03-23 | 6/300 |
|
|
[考研]
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
[考研] 材料科學與工程 317求調(diào)劑 +4 | JKSOIID 2026-03-26 | 4/200 |
|
|
[考研] 081700 調(diào)劑 267分 +11 | 迷人的哈哈 2026-03-23 | 11/550 |
|
|
[考研] 一志愿中南大學化學學碩0703總分337求調(diào)劑 +7 | niko- 2026-03-22 | 7/350 |
|
|
[考研] 290分調(diào)劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
[考研] 344求調(diào)劑 +3 | desto 2026-03-24 | 3/150 |
|
|
[考研] 361求調(diào)劑 +3 | Glack 2026-03-22 | 3/150 |
|