| 8 | 1/1 | 返回列表 |
| 查看: 6330 | 回復(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。有沒有什么辦法呢?本人是非常新的新手,希望能講的詳細一點點。 下面是我的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ā)表 | |
|---|---|---|---|---|
|
[考研] 305分求調(diào)劑(食品工程) +3 | Sxy112 2026-03-21 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 286分人工智能專業(yè)請求調(diào)劑愿意跨考! +4 | lemonzzn 2026-03-17 | 8/400 |
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
[考研] 一志愿深大,0703化學(xué),總分302,求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
[考研] 278求調(diào)劑 +9 | 煙火先于春 2026-03-17 | 9/450 |
|
|
[考研] 302求調(diào)劑 +12 | 呼呼呼。。。。 2026-03-17 | 12/600 |
|
|
[考研] 336求調(diào)劑 +5 | rmc8866 2026-03-21 | 5/250 |
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +5 | 紙魚ly 2026-03-21 | 5/250 |
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
[考研] 材料與化工(0856)304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
[考研] 316求調(diào)劑 +6 | 梁茜雯 2026-03-19 | 6/300 |
|
|
[考研] 化學(xué)求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
[考研] 265求調(diào)劑 +9 | 梁梁校校 2026-03-17 | 9/450 |
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +6 | 想上岸的鯉魚 2026-03-18 | 7/350 |
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生(接收計算機專業(yè)調(diào)劑) +3 | 123阿標 2026-03-17 | 3/150 |
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|