 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 2943 | 回復(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ù)此樓» 猜你喜歡
281求調(diào)劑（0805）
已經(jīng)有8人回復(fù)
-
環(huán)境領(lǐng)域全國重點實驗室招收博士1-2名
已經(jīng)有3人回復(fù)
-
材料專碩306英一數(shù)二
已經(jīng)有10人回復(fù)
-
301求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有7人回復(fù)
-
302求調(diào)劑
已經(jīng)有6人回復(fù)
-
26博士申請
已經(jīng)有3人回復(fù)
-
268求調(diào)劑
已經(jīng)有3人回復(fù)
-
311求調(diào)劑
已經(jīng)有10人回復(fù)
-
被我言中：新模板不強調(diào)格式了，假專家開始管格式了
已經(jīng)有4人回復(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ā)表 | |
|---|---|---|---|---|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +3 | 脫穎而出 2026-03-16 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +6 | 世界首富 2026-03-11 | 6/300 |
|
|
|
[考研] 274求調(diào)劑0856材料化工 +13 | z2839474511 2026-03-11 | 14/700 |
|
|
|
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
|
[考研] 化學(xué)調(diào)劑0703 +8 | 啊我我的 2026-03-11 | 8/400 |
|
|
|
[基金申請]
今年的國基金是打分制嗎？
 50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +8 | 秋有木北 2026-03-14 | 8/400 |
|
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +5 | d如愿上岸 2026-03-12 | 8/400 |
|
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對?/a> 2026-03-13 | 5/250 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油！? 2026-03-15 | 5/250 |
|
|
|
[考研] 材料080500調(diào)劑求收留 +3 | 一顆meteor 2026-03-13 | 3/150 |
|
|
|
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機器學(xué)習(xí) 一志愿哈工大 +4 | 林yaxin 2026-03-12 | 4/200 |
|
|
|
[考研] （081700）化學(xué)工程與技術(shù)-298分求調(diào)劑 +12 | 11啦啦啦 2026-03-11 | 35/1750 |
|
|
|
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9 | 2026調(diào)劑俠 2026-03-12 | 9/450 |
|
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +7 | 18880831720 2026-03-11 | 7/350 |
|
|
|
[論文投稿]
投稿問題
5+4
|
星光燦爛xt 2026-03-12 | 6/300 |
|
|
|
[考研] 0856化學(xué)工程280分求調(diào)劑 +4 | shenzxsn 2026-03-11 | 4/200 |
|