 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2945 | 回復(fù): 5 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
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ù)
HitYingph
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 375
- 帖子: 12
- 在線: 13.3小時
- 蟲號: 5894243
- 注冊: 2017-03-09
- 專業(yè): 固體力學(xué)
3樓2017-07-18 10:22:52
2樓2017-07-18 08:53:10
匿名 
用戶注銷 (著名寫手)
- 應(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
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 303求調(diào)劑 +3 | 睿08 2026-03-17 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 296求調(diào)劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|
[考研] 268求調(diào)劑 +7 | 好運連綿不絕 2026-03-12 | 8/400 |
|
|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-16 | 8/400 |
|
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +16 | 不知道叫什么！ 2026-03-15 | 18/900 |
|
|
|
[考研] 290求調(diào)劑 +6 | 孔志浩 2026-03-12 | 11/550 |
|
|
|
[考研] 0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +3 | 小李0854 2026-03-16 | 3/150 |
|
|
|
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
|
[碩博家園] 深圳大學(xué)碩士招生（2026秋，傳感器方向，僅錄取第一志愿） +4 | xujiaoszu 2026-03-11 | 9/450 |
|
|
|
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +6 | SUICHILD 2026-03-12 | 6/300 |
|
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
|
[考研]
 機械專碩調(diào)劑
+3
|
笨笨兔子 2026-03-12 | 3/150 |
|
|
|
[基金申請]
有必要更換申報口嗎
 20+3
|
fannyamoy 2026-03-11 | 3/150 |
|
|
|
[考研] 308求調(diào)劑 +5 | 是Lupa啊 2026-03-11 | 5/250 |
|
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
|
[考研] 考研調(diào)劑 +4 | 芬達46 2026-03-12 | 4/200 |
|
|
|
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-12 | 4/200 |
|