| 8 | 1/1 | 返回列表 |
| 查看: 6331 | 回復(fù): 7 | ||
fengjingyi鐵蟲 (初入文壇)
|
[求助]
lammps新手用tip4p模型模擬水
|
|
現(xiàn)在想用tip4p模型模擬液態(tài)水,計算時總是出現(xiàn)錯誤ERROR on proc:TIP4P hydrogen has incorrect atom type 一開始用pair_style lj/cut/tip4p/long 1 2 1 1 0.15 8.5 8.5就會出錯,試著換了個勢函數(shù) pair_style lj/cut 8.5 就沒出錯了,所以我覺得可能是勢函數(shù)的問題。不過我比較想用tip4p。有沒有什么辦法呢?本人是非常新的新手,希望能講的詳細(xì)一點點。 下面是我的in文件: units metal atom_style full dimension 3 boundary p p p newton on read_data ye.data pair_style lj/cut/tip4p/long 1 2 1 1 0.15 8.5 8.5 pair_coeff 1 1 0.1521 3.154000 8.5 pair_coeff * 2 0.000000 0.000000 8.5 bond_style harmonic bond_coeff 1 450 0.9572 angle_style harmonic angle_coeff 1 55 104.52 kspace_style pppm/tip4p 0.0001 neigh_modify delay 0 every 1 neighbor 2.0 bin timestep 0.0002 velocity all create 400 72489 mom yes rot yes fix 5 all box/relax iso 10 vmax 0.001 min_style cg minimize 1.0e-12 1.0e-12 10000 10000 unfix 5 reset_timestep 0 fix 2 all npt temp 300 300 0.02 iso 10 10 0.2 pchain 8 drag 0.5 thermo 500 #(每200步輸出熱力學(xué)結(jié)果) thermo_style custom step temp pe etotal press vol density enthalpy xlo xhi ylo yhi zlo zhi dump 1 all cfg 1000 dump2.temp_*.cfg mass type xs ys zs x y z run 10000 |
鐵蟲 (初入文壇)
銀蟲 (小有名氣)
新蟲 (小有名氣)
新蟲 (小有名氣)
|
帖子是很久以前的了,不知道樓主后來怎么辦。分享一點我的經(jīng)驗。 根據(jù)mannual里面的語句介紹 pair_style lj/cut/tip4p/long 1 2 1 1 0.15 8.5 8.5 中第一個“1”表示O原子atom type =1, 后面的“2”表示H原子atom type=2,這個要與輸入的data文件對應(yīng),可能樓主的data文件里是反的,H的atom type=1了。 在 pair_style lj/cut 8.5語句中,沒有關(guān)于這類參數(shù)的定義,只設(shè)定截斷距離,因此不存在與data文件對應(yīng)問題,所以不報錯吧。 希望能給以后遇到同樣問題的木蟲一點小小的幫助~ |
銀蟲 (小有名氣)
銀蟲 (小有名氣)
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 考研調(diào)劑 +4 | 來好運來來來 2026-03-21 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 286求調(diào)劑 +10 | Faune 2026-03-21 | 10/500 |
|
|
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|
[考研] 化學(xué)調(diào)劑 +5 | yzysaa 2026-03-21 | 5/250 |
|
|
[基金申請]
學(xué)校已經(jīng)提交到NSFC,還能修改嗎?
40+4
|
babangida 2026-03-19 | 9/450 |
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級已過 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-16 | 3/150 |
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩,總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂29 2026-03-18 | 5/250 |
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
[考研] 一志愿華中農(nóng)業(yè)071010,總分320求調(diào)劑 +3 | 困困困困坤坤 2026-03-20 | 3/150 |
|
|
[考研] 環(huán)境工程調(diào)劑 +9 | 大可digkids 2026-03-16 | 9/450 |
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|
[考研] 085601專碩,總分342求調(diào)劑,地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|