 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2946 | 回復(fù): 5 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
HitYingph新蟲 (初入文壇)
|
[求助]
Pair style hybrid設(shè)置的問題 已有2人參與
|
|
|
我剛接觸分子動力學(xué)模擬不久，用lammps做碳納米管充氦得仿真，C-C采用airebo力場，C-He 和He-He采用LJ 6 12勢，在進(jìn)行初步的能量最小化弛豫輸出的VDW圖片的時候發(fā)現(xiàn)He原子之間有化學(xué)鍵的存在，我在懷疑是不是我使用的Pair style hybrid設(shè)置的問題，把所有的原子都設(shè)置成了碳原子，導(dǎo)致出現(xiàn)了化學(xué)鍵代碼如下： #atom type 1 C #atom type 2 He pair_style hybrid airebo 3.0 lj/cut 10 pair_coeff * * airebo CH.airebo C C pair_coeff 1 2 lj/cut 0.00165375 3.191 pair_coeff 2 2 lj/cut 0.00093800 2.633 我在這里還很困惑的是，關(guān)于airebo勢里**的定義，按照我這樣的定義，是不是所有的都被設(shè)置成為了airebo勢？我嘗試調(diào)換順序的時候，像這樣： #atom type 1 C #atom type 2 He pair_style hybrid airebo 3.0 lj/cut 10 pair_coeff 1 2 lj/cut 0.00165375 3.191 pair_coeff 2 2 lj/cut 0.00093800 2.633 pair_coeff * * airebo CH.airebo C C 結(jié)果報錯：pair hybrid sub-style is not used 我嘗試將* *替換成1 1的時候還是報錯。那么我這樣的力場寫法有問題么？是不是由于力場設(shè)置問題導(dǎo)致He原子之間有化學(xué)鍵，還是只是VDW顯示的問題？ 附 1.VDW顯示照片已上傳 2.代碼： variable simname index cswnt_he log log.${simname}.txt units metal boundary s s s lattice fcc 1.0 timestep 0.0002 atom_modify map array atom_style full read_data ${simname}.data #atom type 1 C #atom type 2 He pair_style hybrid airebo 3.0 lj/cut 10 pair_coeff * * airebo CH.airebo C C pair_coeff 1 2 lj/cut 0.00165375 3.191 pair_coeff 2 2 lj/cut 0.00093800 2.633 neighbor 2.0 bin neigh_modify delay 10 group cnt type 1 group he type 2 compute stressatom all stress/atom NULL virial compute displaceatom all displace/atom compute prpertyfz all property/atom fz variable fz equal fcm(all,z) variable fzcnt equal fcm(cnt,z) variable fzhe equal fcm(he,z) #compute peatom all pe/atom #compute peatom all reduce sum c_peatom #compute peatomcnt cnt reduce sum c_peatom #compute peatomhe he reduce sum c_peatom thermo 1000 thermo_style custom step temp pe ke etotal evdwl ecoul emol dump 1 all custom 1000 ${simname}.lammpstrj id type x y z min_style cg minimize 1e-25 1e-25 5000 10000 fix fixnvt cnt nvt temp 0.001 0.001 10 run 50000 unfix fixnvt fix fixnve cnt nve run 50000 unfix fixnve write_restart restart.endinitial 1.jpg |
回復(fù)此樓» 猜你喜歡
材料，紡織，生物（0856、0710），化學(xué)招生啦
已經(jīng)有8人回復(fù)
-
0703化學(xué)調(diào)劑 ，六級已過，有科研經(jīng)歷
已經(jīng)有9人回復(fù)
-
工科材料085601 279求調(diào)劑
已經(jīng)有6人回復(fù)
-
070300化學(xué)319求調(diào)劑
已經(jīng)有5人回復(fù)
-
生物學(xué)071000 329分求調(diào)劑
已經(jīng)有3人回復(fù)
-
0703化學(xué)調(diào)劑
已經(jīng)有6人回復(fù)
-
304求調(diào)劑
已經(jīng)有12人回復(fù)
-
288求調(diào)劑，一志愿華南理工大學(xué)071005
已經(jīng)有3人回復(fù)
-
307求調(diào)劑
已經(jīng)有3人回復(fù)
-
312求調(diào)劑
已經(jīng)有8人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 跪求各位大神，幫忙看看小弟的多體勢命令編寫是否合適，求指點(diǎn)��！
已經(jīng)有0人回復(fù)
- Bad Fene bond
已經(jīng)有2人回復(fù)
- 【討論】lammps中利用velocity命令無法控制住速度問題
已經(jīng)有11人回復(fù)
青丘白淺
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1.5
- 帖子: 22
- 在線: 32.2小時
- 蟲號: 3943407
- 注冊: 2015-06-27
- 專業(yè): 高分子材料的加工與成型
5樓2018-05-16 10:43:22
2樓2017-07-18 08:53:10
HitYingph
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 375
- 帖子: 12
- 在線: 13.3小時
- 蟲號: 5894243
- 注冊: 2017-03-09
- 專業(yè): 固體力學(xué)
3樓2017-07-18 10:22:52
匿名 
用戶注銷 (著名寫手)
- 應(yīng)助: 6 (幼兒園)
- 金幣: 2788.2
- 散金: 28
- 紅花: 11
- 帖子: 1474
- 在線: 231.4小時
- 蟲號: 0
- 注冊: 2012-08-11
- 性別: GG
- 專業(yè): 流體力學(xué)
4樓2017-07-21 13:59:03
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 070300化學(xué)319求調(diào)劑 +5 | 錦鯉0909 2026-03-17 | 5/250 |
|
|
|
[考研] 312求調(diào)劑 +7 | 陌宸希 2026-03-16 | 8/400 |
|
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +13 | 大米飯！ 2026-03-15 | 16/800 |
|
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +4 | 手機(jī)用戶 2026-03-17 | 4/200 |
|
|
|
[基金申請] 被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
|
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +4 | share_joy 2026-03-16 | 4/200 |
|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
|
[考研] 285化工學(xué)碩求調(diào)劑（081700） +9 | 柴郡貓_ 2026-03-12 | 9/450 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
|
[考研] 085600調(diào)劑 +5 | 漾漾123sun 2026-03-12 | 6/300 |
|
|
|
[考研]
 機(jī)械專碩調(diào)劑
+3
|
笨笨兔子 2026-03-12 | 3/150 |
|
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
|
[考研] 265求調(diào)劑 +4 | 威化餅07 2026-03-12 | 4/200 |
|
|
|
[考研] 學(xué)碩285求調(diào)劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
|
[考研] 化工學(xué)碩306求調(diào)劑 +9 | 42838695 2026-03-12 | 9/450 |
|
|
|
[考研] 333求調(diào)劑 +3 | 152697 2026-03-12 | 4/200 |
|