 |
 |
 |
| 24小時(shí)熱門(mén)版塊排行榜 |
- >論壇更新日志 (425)
- >導(dǎo)師招生 (42)
- >考研 (38)
- >基金申請(qǐng) (8)
- >考博 (7)
- >蟲(chóng)友互識(shí) (5)
- >碩博家園 (4)
- >招聘信息布告欄 (2)
- >博后之家 (2)
- >專(zhuān)業(yè)外語(yǔ) (1)
- >材料綜合 (1)
- >有機(jī)交流 (1)
- >找工作 (1)
- >論文投稿 (1)
- >休閑灌水 (1)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3078 | 回復(fù): 6 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
[求助]
【求助】lammps模擬碰撞竟然丟失原子，難道是勢(shì)函數(shù)的問(wèn)題？？ 已有2人參與
|
||
 建了一個(gè)模型碰撞的時(shí)候一直丟失兩個(gè)原子，都要瘋了，為了檢查錯(cuò)誤，直接建立了一個(gè)簡(jiǎn)單盒子模擬碰撞，竟然還是報(bào)錯(cuò)丟失兩個(gè)原子！請(qǐng)各路大神幫忙瞄瞄，很簡(jiǎn)單的小程序，咋會(huì)一直報(bào)錯(cuò)呢？請(qǐng)不要用什么截?cái)喟霃竭^(guò)小，還是飛出盒子這樣的話回復(fù)我。因?yàn)槲矣X(jué)得完全不可能飛出盒子，步長(zhǎng)也調(diào)到過(guò)0.00001ps，但是不解決問(wèn)題，而且是周期邊界，丟失原子覺(jué)得很奇怪��！ velocity 0 0 500的時(shí)候沒(méi)問(wèn)題，但是速度一提高就報(bào)錯(cuò)，比如velocity 0 0 1000。實(shí)在是不知道咋辦了，所以請(qǐng)各路大神幫忙看看，小弟先謝過(guò)了�？赡苡行┩瑢W(xué)打開(kāi)附件有問(wèn)題，所以我直接把代碼粘貼出來(lái)了，勢(shì)能是lammps自帶的勢(shì)能 units metal boundary p p p atom_style atomic lattice bcc 2.855 #orient x 1 1 1 orient y -1 -1 2 orient z 1 -1 0 region box block -16 16 -15 15 -20 20 create_box 1 box mass 1 56.0 create_atoms 1 box pair_style eam/alloy pair_coeff * * Fe_mm.eam.fs Fe timestep 0.001 neighbor 1.5 bin neigh_modify every 30 delay 0 check no #displace_box all z delta -2 2 remap none compute pe all pe/atom compute dis all displace/atom dump 1 all custom 500 2.lammpstrj id type x y z vx vy vz c_pe c_dis[1] fix 3 all nvt temp 300 300 0.1 run 1000 unfix 3 write_restart dislocation.restart ——————————————————下面是restart文件 read_restart dislocation.restart lattice bcc 2.855 #orient x 1 1 1 orient y -1 -1 2 orient z 1 -1 0 pair_style eam/alloy pair_coeff * * Fe_mm.eam.fs Fe timestep 0.0001 neighbor 1.5 bin neigh_modify every 10 delay 0 check no compute pe all pe/atom compute dis all displace/atom dump 2 all custom 100 4.lammpstrj id type x y z vx vy vz c_pe c_dis[1] group kpa id 62410 velocity kpa set 0.0 0.0 1000.0 units box fix 5 all nve run 1000 unfix 5 |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲(chóng)僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : 超簡(jiǎn)單碰撞.zip
2016-04-01 10:16:49, 328.17 K
» 猜你喜歡
299求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料專(zhuān)業(yè)求調(diào)劑
已經(jīng)有6人回復(fù)
-
311求調(diào)劑
已經(jīng)有5人回復(fù)
-
317求調(diào)劑
已經(jīng)有13人回復(fù)
-
296求調(diào)劑
已經(jīng)有10人回復(fù)
-
274求調(diào)劑
已經(jīng)有10人回復(fù)
-
22408 344分 求調(diào)劑 一志愿 華電計(jì)算機(jī)技術(shù)
已經(jīng)有4人回復(fù)
-
材料專(zhuān)碩英一數(shù)二306
已經(jīng)有7人回復(fù)
-
南京大學(xué)化學(xué)376求調(diào)劑
已經(jīng)有6人回復(fù)
-
南昌大學(xué)材料專(zhuān)碩311分求調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 【求助】lammps關(guān)于fix ave/spatial命令應(yīng)用問(wèn)題
已經(jīng)有4人回復(fù)
- 【求助】lammps中是否有關(guān)于Bond-orientational order參數(shù)的
已經(jīng)有3人回復(fù)
- 【求助】Lammps使用？
已經(jīng)有0人回復(fù)
- 【求助】LAMMPS-運(yùn)行報(bào)錯(cuò)
已經(jīng)有9人回復(fù)
- 【求助】lammps勢(shì)的選擇
已經(jīng)有3人回復(fù)
- 【求助】lammps模擬聚合物刷如何建模？
已經(jīng)有1人回復(fù)
- 【求助】LAMMPS構(gòu)造晶界
已經(jīng)有0人回復(fù)
- 【求助】lammps或者dlpoly能否做粗粒化模擬嗎？
已經(jīng)有3人回復(fù)
- 【求助】Lammps能做布朗動(dòng)力學(xué)模擬嗎？
已經(jīng)有3人回復(fù)
- 【求助】lammps怎么算聲子譜？
已經(jīng)有1人回復(fù)
- 【求助】lammps的使用
已經(jīng)有9人回復(fù)
- 【求助】lammps 安裝問(wèn)題
已經(jīng)有1人回復(fù)
- 【求助】Lammps 中加外力的問(wèn)題
已經(jīng)有1人回復(fù)
- 【求助】LAMMPS中添加Lees-Edwards邊界條件
已經(jīng)有2人回復(fù)
- 【求助】lammps這段代碼的意思是？
已經(jīng)有1人回復(fù)
- 【求助】lammps計(jì)算heat/flux出錯(cuò)
已經(jīng)有1人回復(fù)
- 【求助】lammps如何測(cè)量模擬以后材料的晶格參數(shù)？
已經(jīng)有7人回復(fù)
- 【求助】LAMMPS－REAX計(jì)算出錯(cuò)
已經(jīng)有3人回復(fù)
- 【求助】lammps中多元化合物的晶格創(chuàng)建問(wèn)題
已經(jīng)有14人回復(fù)
- 【求助】lammps中消除平動(dòng)對(duì)溫度統(tǒng)計(jì)的影響問(wèn)題
已經(jīng)有2人回復(fù)
- 【求助】Lammps 多機(jī)并行的問(wèn)題
已經(jīng)有11人回復(fù)
5樓2016-04-01 18:13:55
hooge
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 1790.1
- 散金: 236
- 紅花: 1
- 帖子: 975
- 在線: 135.1小時(shí)
- 蟲(chóng)號(hào): 89636
- 注冊(cè): 2005-08-29
- 專(zhuān)業(yè): 金屬材料的力學(xué)行為
2樓2016-04-01 13:09:58
3樓2016-04-01 13:15:09
whyaquarius
新蟲(chóng) (初入文壇)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 253.8
- 散金: 70
- 帖子: 30
- 在線: 34.7小時(shí)
- 蟲(chóng)號(hào): 2175398
- 注冊(cè): 2012-12-09
- 性別: GG
- 專(zhuān)業(yè): 核技術(shù)及其應(yīng)用
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感謝參與，應(yīng)助指數(shù) +1
大懶貓王浩: 金幣+300, ★★★★★最佳答案, 師兄竟然來(lái)賺金幣了，，，，不好意思，只能給師兄了。。。 2016-04-01 17:29:35
月只藍(lán): IP檢查，初步認(rèn)定為馬甲行為，涉嫌轉(zhuǎn)移金幣，違反論壇相關(guān)管理規(guī)定：http://emuch.net/bbs/viewthread.php?tid=3711903 2016-04-02 10:16:56
月只藍(lán): 經(jīng)查，您和大懶貓王浩的IP地址十分相近，請(qǐng)于2天內(nèi)通過(guò)站內(nèi)短消息解釋您和大懶貓王浩的關(guān)系！逾期按馬甲處理，并可能追加額外金幣處罰。 2016-04-02 10:18:04
月只藍(lán): 金幣-50, 認(rèn)定為金幣轉(zhuǎn)移行為，核減金幣300×0.8=240個(gè)，現(xiàn)扣除50個(gè)，還需要扣除190個(gè)。 2016-04-07 13:38:49
月只藍(lán): 金幣-50, 認(rèn)定為金幣轉(zhuǎn)移行為，核減金幣300×0.8=240個(gè)，現(xiàn)扣除50個(gè)，還需要扣除140個(gè)。 2016-04-07 13:39:08
月只藍(lán): 金幣-50, 認(rèn)定為金幣轉(zhuǎn)移行為，核減金幣300×0.8=240個(gè)，現(xiàn)扣除50個(gè)，還需要扣除90個(gè)。 2016-04-07 13:39:34
月只藍(lán): 金幣-50, 認(rèn)定為金幣轉(zhuǎn)移行為，核減金幣300×0.8=240個(gè)，現(xiàn)扣除50個(gè)，還需要扣除40個(gè)。 2016-04-07 13:40:20
月只藍(lán): 金幣-40, 認(rèn)定為金幣轉(zhuǎn)移行為，核減金幣300×0.8=240個(gè)，現(xiàn)已扣除完畢。 2016-04-07 13:40:42
月只藍(lán): 考慮到用戶可能不了解馬甲行為規(guī)定且是初犯，免除額外處罰。 2016-04-07 13:42:00
感謝參與，應(yīng)助指數(shù) +1
大懶貓王浩: 金幣+300, ★★★★★最佳答案, 師兄竟然來(lái)賺金幣了，，，，不好意思，只能給師兄了。。。 2016-04-01 17:29:35
月只藍(lán): IP檢查，初步認(rèn)定為馬甲行為，涉嫌轉(zhuǎn)移金幣，違反論壇相關(guān)管理規(guī)定：http://emuch.net/bbs/viewthread.php?tid=3711903 2016-04-02 10:16:56
月只藍(lán): 經(jīng)查，您和大懶貓王浩的IP地址十分相近，請(qǐng)于2天內(nèi)通過(guò)站內(nèi)短消息解釋您和大懶貓王浩的關(guān)系！逾期按馬甲處理，并可能追加額外金幣處罰。 2016-04-02 10:18:04
月只藍(lán): 金幣-50, 認(rèn)定為金幣轉(zhuǎn)移行為，核減金幣300×0.8=240個(gè)，現(xiàn)扣除50個(gè)，還需要扣除190個(gè)。 2016-04-07 13:38:49
月只藍(lán): 金幣-50, 認(rèn)定為金幣轉(zhuǎn)移行為，核減金幣300×0.8=240個(gè)，現(xiàn)扣除50個(gè)，還需要扣除140個(gè)。 2016-04-07 13:39:08
月只藍(lán): 金幣-50, 認(rèn)定為金幣轉(zhuǎn)移行為，核減金幣300×0.8=240個(gè)，現(xiàn)扣除50個(gè)，還需要扣除90個(gè)。 2016-04-07 13:39:34
月只藍(lán): 金幣-50, 認(rèn)定為金幣轉(zhuǎn)移行為，核減金幣300×0.8=240個(gè)，現(xiàn)扣除50個(gè)，還需要扣除40個(gè)。 2016-04-07 13:40:20
月只藍(lán): 金幣-40, 認(rèn)定為金幣轉(zhuǎn)移行為，核減金幣300×0.8=240個(gè)，現(xiàn)已扣除完畢。 2016-04-07 13:40:42
月只藍(lán): 考慮到用戶可能不了解馬甲行為規(guī)定且是初犯，免除額外處罰。 2016-04-07 13:42:00
|
來(lái)來(lái)來(lái)，我來(lái)賺金幣啦~ 在輻照損傷的模擬中，在離子間距離比較小時(shí)，已經(jīng)不能用平衡態(tài)下的作用勢(shì)來(lái)描述，一般用 ZBL 來(lái)描述這個(gè)碰撞過(guò)程。比如說(shuō)常用到的 tersoff/zbl 作用勢(shì)就勢(shì)就是將 Tersoff 作用勢(shì)與 ZBL 作用勢(shì)通過(guò)一個(gè)費(fèi)米函數(shù)連接起來(lái)，一般 認(rèn)為小于 1Å 時(shí)需要用 ZBL 描述。 P.S. eam勢(shì)木有官方的zbl銜接版本~ |
4樓2016-04-01 17:20:32
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
|
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
|
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分，有科研經(jīng)歷，求調(diào)劑 +4 | 一瓶苯甲酸 2026-03-14 | 4/200 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 ，六級(jí)已過(guò)，有科研經(jīng)歷 +13 | 曦熙兮 2026-03-15 | 13/650 |
|
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
|
[基金申請(qǐng)]
學(xué)校已經(jīng)提交到NSFC，還能修改嗎？
 40+4
|
babangida 2026-03-19 | 8/400 |
|
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
|
[考研] 324分 085600材料化工求調(diào)劑 +3 | llllkkkhh 2026-03-18 | 3/150 |
|
|
|
[考研] 332求調(diào)劑 +3 | ydfyh 2026-03-17 | 3/150 |
|
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
|
[考研] 293求調(diào)劑 +11 | zjl的號(hào) 2026-03-16 | 16/800 |
|
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
|
[考研] 一志愿，福州大學(xué)材料專(zhuān)碩339分求調(diào)劑 +3 | 木子momo青爭(zhēng) 2026-03-15 | 3/150 |
|
|
|
[考研] 機(jī)械專(zhuān)碩325，尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|