 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 1084 | 回復: 2 | ||
[求助]
LAMMPS能量計算值不正確
|
|
大家好，關(guān)于LAMMPS模擬中，能量相關(guān)的部分，有一些問題想要請教各位。 我所采用的in文件及data文件，是moltemplate軟件中的OPLS-AA力場正丁烷范例，條件為260K，1atm。在按照范例運行之后，發(fā)現(xiàn)內(nèi)能/焓的計算值相比于NIST的數(shù)據(jù)小50%左右。我嘗試了其他的溫度和壓力，發(fā)現(xiàn)能量（U及H）的計算值要比NIST數(shù)據(jù)小很多，而且集中在4-5kcal/mol附近，變化很小。 我嘗試修改了special_bonds的參數(shù)，雖然值發(fā)生變化但這應(yīng)該是不正確的。因此想要請教各位是否出現(xiàn)過同樣的問題，以及對應(yīng)的解決方法。 謝謝大家！ |
回復此樓» 猜你喜歡
研究生招生
已經(jīng)有6人回復
-
材料與化工，291，求調(diào)劑
已經(jīng)有12人回復
-
0817一志愿蘇州大學280復試調(diào)劑
已經(jīng)有8人回復
-
求調(diào)劑 一志愿蘇州大學，0856化工323分 | 本科應(yīng)化 | 有專利/競賽/科研助手經(jīng)歷 |
已經(jīng)有9人回復
-
085602化工求調(diào)劑
已經(jīng)有4人回復
-
083000環(huán)境科學與工程調(diào)劑
已經(jīng)有6人回復
-
307求調(diào)劑
已經(jīng)有3人回復
-
求材料讀博院校
已經(jīng)有9人回復
-
哈工大材料324求調(diào)劑
已經(jīng)有4人回復
-
材料調(diào)劑，307分
已經(jīng)有10人回復
|
這里是In文件中的運行部分 # ------------------------------- Run Section ------------------------------- dump dumpeq1 all custom 50 traj_eq1_min.lammpstrj id mol type x y z ix iy iz thermo_style custom step temp press ke etotal enthalpy thermo 50 thermo_modify norm yes velocity all create 900.0 4928459 dist gaussian # -- Equilibration: part 1: initial minimization -- minimize 1.0e-5 1.0e-7 100000 400000 undump dumpeq1 write_data system_after_eq1_min.data # -- Equilibration part 2: reorienting the molecules (NVT) -- timestep 1.0 dump dumpeq2 all custom 200 traj_eq2_reorient.lammpstrj id mol type x y z ix iy iz fix fxlan all langevin 900.0 900.0 120 48279 fix_modify fxlan energy yes fix fxnve all nve run 4000 velocity all scale 900.0 sum yes unfix fxlan unfix fxnve fix fxnvt all nvt temp 900.0 900.0 100.0 fix_modify fxnvt energy yes run 10000 undump dumpeq2 write_data system_after_eq2_reorient.data unfix fxnvt # -- equilibration part 3: Equilibrating the density (NPT) -- dump dumpeq3 all custom 200 traj_eq3_npt.lammpstrj id mol type x y z ix iy iz fix fxnpt all npt temp 900.0 260.0 100.0 iso 500.0 500.0 1000.0 drag 2.0 timestep 1.0 run 20000 unfix fxnpt fix fxnpt all npt temp 260.0 260.0 100.0 iso 500.0 1.0 1000.0 drag 2.0 timestep 1.0 run 100000 write_data system_after_eq3_npt.data |
2樓2019-05-05 10:07:04
3樓2019-05-05 10:08:31
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 Recycled 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿天津大學，英一數(shù)二305分求調(diào)劑，四六級已過 +5 | 小小番的茄 2026-03-09 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 307求調(diào)劑 +3 | 辛仁豆腐 2026-03-08 | 5/250 |
|
|
|
[考研] 求調(diào)劑，數(shù)一英一274分 +4 | 小菲會努力 2026-03-08 | 4/200 |
|
|
|
[考研] 材料工程330分求調(diào)劑，一志愿985 +3 | 小材化本科 2026-03-07 | 3/150 |
|
|
|
[考研] 288求調(diào)劑 +5 | 好好- 2026-03-07 | 5/250 |
|
|
|
[考研] 347求調(diào)劑 +4 | 浮云滿足 2026-03-07 | 4/200 |
|
|
|
[基金申請] 進入個人成果庫好難，一下午都沒進去 +6 | mi_dilee 2026-03-05 | 6/300 |
|
|
|
[考研] 0703化學調(diào)劑 +5 | G212 2026-03-03 | 6/300 |
|
|
|
[考研] 第一志愿上海大學，專業(yè)化學工程與技術(shù)，總分288，求調(diào)劑 +3 | 1829197082 2026-03-07 | 3/150 |
|
|
|
[考研] 一志愿哈爾濱工業(yè)大學0856材料與化工，前三科206，總分283，求調(diào)劑 +7 | 26考研求調(diào)劑 2026-03-06 | 7/350 |
|
|
|
[考研] 一志愿南京大學資源與環(huán)境專業(yè)283分求調(diào)劑 +8 | 追云?? 2026-03-04 | 8/400 |
|
|
|
[考研] 308求調(diào)劑 +7 | 倘若起風了呢 2026-03-05 | 9/450 |
|
|
|
[考研] 282求調(diào)劑 +7 | 夕～日 2026-03-05 | 8/400 |
|
|
|
[考研] 化工282求調(diào)劑一志愿211 +5 | NA0912 2026-03-05 | 6/300 |
|
|
|
[考研] 一志愿武理085601專碩347分求調(diào)劑 +4 | 啊歐歐歐 2026-03-04 | 5/250 |
|
|
|
[考研] 材料085600 303求調(diào)劑 +7 | 1bygone 2026-03-04 | 7/350 |
|
|
|
[考研] 22408-273求調(diào)劑-擔任3個項目的負責人-1篇國際期刊論文（一作）1篇核心期刊論文在投。 +3 | 沒想好取什么名 2026-03-03 | 3/150 |
|
|
|
[考研] 一志愿314求調(diào)劑 +7 | 202111120625 2026-03-03 | 7/350 |
|
|
|
[考研] 312求調(diào)劑 +8 | 醋精華了一下發(fā)?/a> 2026-03-03 | 10/500 |
|
|
|
[考研]
 0703化學學碩 288分
+4
|
旅行中的紫葡萄 2026-03-03 | 4/200 |
|