| 3 | 1/1 | 返回列表 |
| 查看: 3312 | 回復: 2 | ||
1207150929金蟲 (小有名氣)
|
[求助]
關(guān)于用lammps建立無限延伸聚合物晶體模型的問題
|
|
我現(xiàn)在用lammps(17 Oct 2012 version)用非平衡分子動力學的方法做聚合物晶體聚合物熱導率的計算。模擬建模時,通過拓展polymer的unit cell, 生成supercell,這個supercell即為模擬的box,力場用GAFF,周期性邊界條件加在box的三個方向。 我的具體建模做法如下: 1、用antechamber處理gaussian的輸出文件,生成polymer單胞內(nèi)單體的mol2文件和frcmod文件; 2、用tleap工具將單體的力場信息load到超胞中,生成超胞的prmtop和inpcrd文件; 3、用amber2lammps將超胞的prmtop和inpcrd文件轉(zhuǎn)化成lammps可讀的data文件。 但由于我需要模擬無限延伸的polymer體系,故在box的邊界總會有bond穿越box的邊界。所以我在使用tleap工具將單體的力場信息load到超胞中時,就有意將 polymer鏈最兩端的原子bond到一起。這樣生成的data文件并無錯誤信息。但之后我測試用lammps模擬能量最小化時,在log輸出文件中總會出現(xiàn)錯誤信息: Bond/angle/dihedral extent > half of periodic box length,模擬過程終止。我不知道是不是我建模的過程有問題;我也不確定lammps在模擬的過程中 是否遵循最小鏡像原理,如果遵循的話,lammps應該就不會誤認為將鏈最兩端的原子bond到一起,而是一個box中的最邊緣原子和它近鄰box中的鏡像bond到一起。 |
銅蟲 (小有名氣)

銀蟲 (小有名氣)
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 288求調(diào)劑,一志愿華南理工大學071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 10/500 |
|
|
[考研] 317求調(diào)劑 +3 | 申子申申 2026-03-19 | 6/300 |
|
|
[考研] 化學求調(diào)劑 +3 | 臨澤境llllll 2026-03-17 | 4/200 |
|
|
[考研] 0703化學調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
[考研] 0817調(diào)劑 +3 | 沒有答案_ 2026-03-14 | 3/150 |
|
|
[考研] 304求調(diào)劑 +6 | 司空. 2026-03-18 | 6/300 |
|
|
[考研] 085410人工智能專碩317求調(diào)劑(0854都可以) +3 | xbxudjdn 2026-03-18 | 3/150 |
|
|
[考研] 085600材料與化工 +5 | 安全上岸! 2026-03-16 | 5/250 |
|
|
[考研] 0854可跨調(diào)劑,一作一項核心論文五項專利,省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
[考研] 268求調(diào)劑 +6 | 簡單點0 2026-03-17 | 6/300 |
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學123 2026-03-15 | 8/400 |
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
[考研] 274求調(diào)劑 +5 | 時間點 2026-03-13 | 5/250 |
|
|
[考研] 26考研一志愿中國石油大學(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
[考研] 297求調(diào)劑 +4 | 學海漂泊 2026-03-13 | 4/200 |
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|