 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 2942 | 回復(fù): 5 | ||
HitYingph新蟲 (初入文壇)
|
[求助]
Pair style hybrid設(shè)置的問題 已有2人參與
|
|
我剛接觸分子動力學(xué)模擬不久，用lammps做碳納米管充氦得仿真，C-C采用airebo力場，C-He 和He-He采用LJ 6 12勢，在進行初步的能量最小化弛豫輸出的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ù)此樓» 猜你喜歡
268求調(diào)劑
已經(jīng)有3人回復(fù)
-
302求調(diào)劑
已經(jīng)有5人回復(fù)
-
311求調(diào)劑
已經(jīng)有10人回復(fù)
-
被我言中：新模板不強調(diào)格式了，假專家開始管格式了
已經(jīng)有4人回復(fù)
-
303求調(diào)劑
已經(jīng)有5人回復(fù)
-
考研求調(diào)劑
已經(jīng)有4人回復(fù)
-
085601專碩，總分342求調(diào)劑，地區(qū)不限
已經(jīng)有4人回復(fù)
-
085601材料工程專碩求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有6人回復(fù)
-
能源材料化學(xué)課題組招收碩士研究生8-10名
已經(jīng)有6人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 跪求各位大神，幫忙看看小弟的多體勢命令編寫是否合適，求指點��！
已經(jīng)有0人回復(fù)
- Bad Fene bond
已經(jīng)有2人回復(fù)
- 【討論】lammps中利用velocity命令無法控制住速度問題
已經(jīng)有11人回復(fù)
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
青丘白淺
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1.5
- 帖子: 22
- 在線: 32.2小時
- 蟲號: 3943407
- 注冊: 2015-06-27
- 專業(yè): 高分子材料的加工與成型
5樓2018-05-16 10:43:22
6樓2018-06-08 10:44:20
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 HitYingph 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 工科材料085601 279求調(diào)劑 +3 | 困于星晨 2026-03-17 | 3/150 |
|
|
|
[考研] 08工科 320總分 求調(diào)劑 +4 | 梨花珞晚風(fēng) 2026-03-17 | 4/200 |
|
|
[考研] 286求調(diào)劑 +3 | lemonzzn 2026-03-16 | 5/250 |
|
|
|
[考研] 一志愿華中師范071000，325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
|
[考研] 308 085701 四六級已過求調(diào)劑 +7 | 溫喬喬喬喬 2026-03-12 | 14/700 |
|
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
|
[考研] 一志愿中科院，化學(xué)方向，295求調(diào)劑 +4 | 一氧二氮 2026-03-11 | 4/200 |
|
|
|
[考研] 求調(diào)劑（材料與化工327） +4 | 愛吃香菜啦 2026-03-11 | 4/200 |
|
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
|
[考研] （081700）化學(xué)工程與技術(shù)-298分求調(diào)劑 +12 | 11啦啦啦 2026-03-11 | 35/1750 |
|
|
|
[考研] 315求調(diào)劑 +9 | 小羊小羊_ 2026-03-11 | 10/500 |
|
|
|
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9 | 2026調(diào)劑俠 2026-03-12 | 9/450 |
|
|
|
[考研] 材料與化工085600調(diào)劑求老師收留 +9 | jiaanl 2026-03-11 | 9/450 |
|
|
|
[考研] 材料專碩350 求調(diào)劑 +4 | 王金科 2026-03-12 | 4/200 |
|
|
|
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|