| 3 | 1/1 | 返回列表 |
| 查看: 3311 | 回復: 2 | ||
1207150929金蟲 (小有名氣)
|
[求助]
關于用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ā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學666 2026-03-16 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +3 | Mqqqqqq 2026-03-19 | 3/150 |
|
|
[考研] 化學求調(diào)劑 +3 | 臨澤境llllll 2026-03-17 | 4/200 |
|
|
[考研] 一志愿南昌大學,327分,材料與化工085600 +3 | Ncdx123456 2026-03-19 | 3/150 |
|
|
[考研] 0817調(diào)劑 +3 | 沒有答案_ 2026-03-14 | 3/150 |
|
|
[考研] 285化工學碩求調(diào)劑(081700) +11 | 柴郡貓_ 2026-03-12 | 11/550 |
|
|
[考研] 288求調(diào)劑,一志愿華南理工大學071005 +4 | ioodiiij 2026-03-17 | 4/200 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考研] 0703化學求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
[考研] 一志愿南京大學,080500材料科學與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] 機械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
[基金申請] 現(xiàn)在如何回避去年的某一個專家,不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
[考研] 復試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|