| 3 | 1/1 | 返回列表 |
| 查看: 2231 | 回復(fù): 2 | ||
[求助]
計算W-Mo合金的熔點,程序lost atoms,請教能人解答 已有1人參與
|
|
參照上海交大那個求Cu或Al熔點的例子,把lattice部分換成了我Mo-0.25W成分合金的原子分布,in文件如下: ++++++++++++++++++++++++++++++++++++++++++++++++++++++ in.melt_temp ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ # LAMMPS Melt_Mo-W_temp units metal boundary p p p atom_style atomic variable x index 10 500 1000 1800 2000 variable y equal 3.1466/15.733 variable z universe 15.733 lattice custom $z a1 0.4 0.0 0.0 a2 0.0 1.0 0.0 a3 0.0 0.0 $y basis 0.0 0.0 0.0 basis 0.25 0.1 0.5 basis 0.0 0.2 0.0 basis 0.0 0.4 0.0 basis 0.0 0.6 0.0 basis 0.0 0.8 0.0 basis 0.25 0.9 0.5 basis 0.5 0.0 0.0 basis 0.75 0.1 0.5 basis 0.75 0.5 0.5 basis 0.5 0.8 0.0 basis 0.75 0.9 0.5 basis 0.5 0.4 0.0 basis 0.5 0.6 0.0 basis 0.5 0.2 0.0 basis 0.25 0.5 0.5 basis 0.25 0.3 0.5 basis 0.25 0.7 0.5 basis 0.75 0.7 0.5 basis 0.75 0.3 0.5 lattice custom 15.733 a1 0.4 0.0 0.0 a2 0.0 1.0 0.0 a3 0.0 0.0 $y basis 0.0 0.0 0.0 basis 0.25 0.1 0.5 basis 0.0 0.2 0.0 basis 0.0 0.4 0.0 basis 0.0 0.6 0.0 basis 0.0 0.8 0.0 basis 0.25 0.9 0.5 basis 0.5 0.0 0.0 basis 0.75 0.1 0.5 basis 0.75 0.5 0.5 basis 0.5 0.8 0.0 basis 0.75 0.9 0.5 basis 0.5 0.4 0.0 basis 0.5 0.6 0.0 basis 0.5 0.2 0.0 basis 0.25 0.5 0.5 basis 0.25 0.3 0.5 basis 0.25 0.7 0.5 basis 0.75 0.7 0.5 basis 0.75 0.3 0.5 lattice custom 15.733 a1 0.4 0.0 0.0 a2 0.0 1.0 0.0 a3 0.0 0.0 0.2 basis 0.0 0.0 0.0 basis 0.25 0.1 0.5 basis 0.0 0.2 0.0 basis 0.0 0.4 0.0 basis 0.0 0.6 0.0 basis 0.0 0.8 0.0 basis 0.25 0.9 0.5 basis 0.5 0.0 0.0 basis 0.75 0.1 0.5 basis 0.75 0.5 0.5 basis 0.5 0.8 0.0 basis 0.75 0.9 0.5 basis 0.5 0.4 0.0 basis 0.5 0.6 0.0 basis 0.5 0.2 0.0 basis 0.25 0.5 0.5 basis 0.25 0.3 0.5 basis 0.25 0.7 0.5 basis 0.75 0.7 0.5 basis 0.75 0.3 0.5 region box block -4.0 4.0 -3.0 3.0 0.1 45.0 create_box 2 box region Mo-0.25W block -4.0 4.0 -3.0 3.0 0.1 45.0 create_atoms 2 region Mo-0.25W basis 1 1 basis 2 1 basis 3 1 basis 4 1 basis 5 1 basis 6 1 basis 7 1 basis 8 1 basis 9 1 basis 10 1 basis 11 1 basis 12 1 basis 13 1 basis 14 1 basis 15 1 basis 16 2 basis 17 2 basis 18 2 basis 19 2 basis 20 2 group Mo-0.25W region Mo-0.25W region Mo-0.25W11 block -4.0 4.0 -3.0 3.0 43.0 45.0 region Mo-0.25W22 block -4.0 4.0 -3.0 3.0 0.1 3.0 group Mo-0.25W11 region Mo-0.25W11 group Mo-0.25W22 region Mo-0.25W22 fix region Mo-0.25W11 setforce 0 0 0 fix region Mo-0.25W22 setforce 0 0 0 mass 1 95.941 mass 2 183.841 #fix wall all wall/reflect zlo zhi timestep 0.01 thermo 100 pair_style hybrid morse 10.0 pair_coeff 1 1 morse 0.77019 1.434 3.012 pair_coeff 1 2 morse 13.68932 1.82699 2.30154 pair_coeff 2 2 morse 0.969483 1.385 3.053 neighbor 0.6 bin neigh_modify every 5 delay 0 check yes velocity all create $x 825577 dist gaussian fix 1 all npt $x $x 2.0 xyz 0.0 0.0 45.0 drag 0.2 #在x溫度下保持外壓為零弛豫 compute 3 all pe/atom compute 4 all ke/atom compute 5 all coord/atom 3.0 run 10000 unfix 1 fix 1 all nvt $x $x 2.0 drag 0.2 #NVT控溫 dump 1 all custom 50 dump_$x.atom id xs ys zs c_3 c_4 c_5 fix 2 all msd 1 msd_Mo-W_$x.dat #輸出MSD文件 run 10000 clear next x jump in.melt_Mo-W_temp +++++++++++++++++++++++++++++++++++++++++++++++++++= 以下是log文件 +++++++++++++++++++++++++++++++++++++++++++++++++++++ LAMMPS (7 Jul 2009) jump in.melt_temp # LAMMPS Melt_Mo-W_temp units metal boundary p p p atom_style atomic variable x index 10 500 1000 1800 2000 variable y equal 3.1466/15.733 variable z universe 15.733 Initial ${z} setting: value 1 on partition 0 lattice custom $z a1 0.4 0.0 0.0 a2 0.0 1.0 0.0 a3 0.0 0.0 $y basis 0.0 0.0 0.0 basis 0.25 0.1 0.5 basis 0.0 0.2 0.0 basis 0.0 0.4 0.0 basis 0.0 0.6 0.0 basis 0.0 0.8 0.0 basis 0.25 0.9 0.5 basis 0.5 0.0 0.0 basis 0.75 0.1 0.5 basis 0.75 0.5 0.5 basis 0.5 0.8 0.0 basis 0.75 0.9 0.5 basis 0.5 0.4 0.0 basis 0.5 0.6 0.0 basis 0.5 0.2 0.0 basis 0.25 0.5 0.5 basis 0.25 0.3 0.5 basis 0.25 0.7 0.5 basis 0.75 0.7 0.5 basis 0.75 0.3 0.5 lattice custom 15.733 a1 0.4 0.0 0.0 a2 0.0 1.0 0.0 a3 0.0 0.0 $y basis 0.0 0.0 0.0 basis 0.25 0.1 0.5 basis 0.0 0.2 0.0 basis 0.0 0.4 0.0 basis 0.0 0.6 0.0 basis 0.0 0.8 0.0 basis 0.25 0.9 0.5 basis 0.5 0.0 0.0 basis 0.75 0.1 0.5 basis 0.75 0.5 0.5 basis 0.5 0.8 0.0 basis 0.75 0.9 0.5 basis 0.5 0.4 0.0 basis 0.5 0.6 0.0 basis 0.5 0.2 0.0 basis 0.25 0.5 0.5 basis 0.25 0.3 0.5 basis 0.25 0.7 0.5 basis 0.75 0.7 0.5 basis 0.75 0.3 0.5 lattice custom 15.733 a1 0.4 0.0 0.0 a2 0.0 1.0 0.0 a3 0.0 0.0 0.2 basis 0.0 0.0 0.0 basis 0.25 0.1 0.5 basis 0.0 0.2 0.0 basis 0.0 0.4 0.0 basis 0.0 0.6 0.0 basis 0.0 0.8 0.0 basis 0.25 0.9 0.5 basis 0.5 0.0 0.0 basis 0.75 0.1 0.5 basis 0.75 0.5 0.5 basis 0.5 0.8 0.0 basis 0.75 0.9 0.5 basis 0.5 0.4 0.0 basis 0.5 0.6 0.0 basis 0.5 0.2 0.0 basis 0.25 0.5 0.5 basis 0.25 0.3 0.5 basis 0.25 0.7 0.5 basis 0.75 0.7 0.5 basis 0.75 0.3 0.5 Lattice spacing in x,y,z = 6.2932 15.733 3.1466 lattice custom 15.733 a1 0.4 0.0 0.0 a2 0.0 1.0 0.0 a3 0.0 0.0 $y basis 0.0 0.0 0.0 basis 0.25 0.1 0.5 basis 0.0 0.2 0.0 basis 0.0 0.4 0.0 basis 0.0 0.6 0.0 basis 0.0 0.8 0.0 basis 0.25 0.9 0.5 basis 0.5 0.0 0.0 basis 0.75 0.1 0.5 basis 0.75 0.5 0.5 basis 0.5 0.8 0.0 basis 0.75 0.9 0.5 basis 0.5 0.4 0.0 basis 0.5 0.6 0.0 basis 0.5 0.2 0.0 basis 0.25 0.5 0.5 basis 0.25 0.3 0.5 basis 0.25 0.7 0.5 basis 0.75 0.7 0.5 basis 0.75 0.3 0.5 lattice custom 15.733 a1 0.4 0.0 0.0 a2 0.0 1.0 0.0 a3 0.0 0.0 0.2 basis 0.0 0.0 0.0 basis 0.25 0.1 0.5 basis 0.0 0.2 0.0 basis 0.0 0.4 0.0 basis 0.0 0.6 0.0 basis 0.0 0.8 0.0 basis 0.25 0.9 0.5 basis 0.5 0.0 0.0 basis 0.75 0.1 0.5 basis 0.75 0.5 0.5 basis 0.5 0.8 0.0 basis 0.75 0.9 0.5 basis 0.5 0.4 0.0 basis 0.5 0.6 0.0 basis 0.5 0.2 0.0 basis 0.25 0.5 0.5 basis 0.25 0.3 0.5 basis 0.25 0.7 0.5 basis 0.75 0.7 0.5 basis 0.75 0.3 0.5 Lattice spacing in x,y,z = 6.2932 15.733 3.1466 lattice custom 15.733 a1 0.4 0.0 0.0 a2 0.0 1.0 0.0 a3 0.0 0.0 0.2 basis 0.0 0.0 0.0 basis 0.25 0.1 0.5 basis 0.0 0.2 0.0 basis 0.0 0.4 0.0 basis 0.0 0.6 0.0 basis 0.0 0.8 0.0 basis 0.25 0.9 0.5 basis 0.5 0.0 0.0 basis 0.75 0.1 0.5 basis 0.75 0.5 0.5 basis 0.5 0.8 0.0 basis 0.75 0.9 0.5 basis 0.5 0.4 0.0 basis 0.5 0.6 0.0 basis 0.5 0.2 0.0 basis 0.25 0.5 0.5 basis 0.25 0.3 0.5 basis 0.25 0.7 0.5 basis 0.75 0.7 0.5 basis 0.75 0.3 0.5 Lattice spacing in x,y,z = 6.2932 15.733 3.1466 region box block -4.0 4.0 -3.0 3.0 0.1 45.0 create_box 2 box Created orthogonal box = (-25.1728 -47.199 0.31466) to (25.1728 47.199 141.597) 1 by 1 by 1 processor grid region Mo-0.25W block -4.0 4.0 -3.0 3.0 0.1 45.0 create_atoms 2 region Mo-0.25W basis 1 1 basis 2 1 basis 3 1 basis 4 1 basis 5 1 basis 6 1 basis 7 1 basis 8 1 basis 9 1 basis 10 1 basis 11 1 basis 12 1 basis 13 1 basis 14 1 basis 15 1 basis 16 2 basis 17 2 basis 18 2 basis 19 2 basis 20 2 Created 42720 atoms group Mo-0.25W region Mo-0.25W 42720 atoms in group Mo-0.25W region Mo-0.25W11 block -4.0 4.0 -3.0 3.0 43.0 45.0 region Mo-0.25W22 block -4.0 4.0 -3.0 3.0 0.1 3.0 group Mo-0.25W11 region Mo-0.25W11 1920 atoms in group Mo-0.25W11 group Mo-0.25W22 region Mo-0.25W22 2880 atoms in group Mo-0.25W22 fix region Mo-0.25W11 setforce 0 0 0 fix region Mo-0.25W22 setforce 0 0 0 mass 1 95.941 mass 2 183.841 #fix wall all wall/reflect zlo zhi timestep 0.01 thermo 100 pair_style hybrid morse 10.0 pair_coeff 1 1 morse 0.77019 1.434 3.012 pair_coeff 1 2 morse 13.68932 1.82699 2.30154 pair_coeff 2 2 morse 0.969483 1.385 3.053 neighbor 0.6 bin neigh_modify every 5 delay 0 check yes velocity all create $x 825577 dist gaussian velocity all create 10 825577 dist gaussian fix 1 all npt $x $x 2.0 xyz 0.0 0.0 45.0 drag 0.2 #在x溫度下保持外壓為零弛豫 fix 1 all npt 10 $x 2.0 xyz 0.0 0.0 45.0 drag 0.2 fix 1 all npt 10 10 2.0 xyz 0.0 0.0 45.0 drag 0.2 compute 3 all pe/atom compute 4 all ke/atom compute 5 all coord/atom 3.0 run 10000 Memory usage per processor = 40.9131 Mbytes Step Temp E_pair E_mol TotEng Press Volume 0 10 -1187595.7 0 -1187540.5 -3305841.1 671447.7 ERROR: Lost atoms: original 42720 current 41050 +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 請教一下大神們,這個問題要怎么解決?我試圖加一段能量最小化的代碼,可是運行后是各種數(shù)據(jù)都輸出為0……同樣不知道怎么解決~ 請教大神,望各位能積極援助,不勝感激~~~ @老虎大王 @月只藍 |

木蟲 (著名寫手)
| 丟失原子了。主要檢查勢函數(shù)和初始構(gòu)型。先做一個簡單的合金構(gòu)型,比如B2,L12之類的,看看能不能正常運算,能量是否正確。 |
新蟲 (小有名氣)
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 化學(xué)工程321分求調(diào)劑 +17 | 大米飯! 2026-03-15 | 21/1050 |
|
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +12 | 冬十三 2026-03-15 | 13/650 |
|
|
[考研] 0703化學(xué) 305求調(diào)劑 +5 | FY_yy 2026-03-14 | 5/250 |
|
|
[考研] 求調(diào)劑 +5 | 十三加油 2026-03-21 | 5/250 |
|
|
[考研] 工科0856求調(diào)劑 +3 | 沐析汀汀 2026-03-21 | 3/150 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
[考研] 324分 085600材料化工求調(diào)劑 +4 | llllkkkhh 2026-03-18 | 4/200 |
|
|
[考研] 一志愿華南師大 070300(化學(xué))304分求調(diào)劑 +3 | 0703武芊慧雪304 2026-03-18 | 3/150 |
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-18 | 6/300 |
|
|
[考研] 一志愿西南交大,求調(diào)劑 +5 | 材化逐夢人 2026-03-18 | 5/250 |
|
|
[考研] 一志愿南昌大學(xué),327分,材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑,總分315(英一) +5 | sbdksD 2026-03-19 | 5/250 |
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
[論文投稿]
申請回稿延期一個月,編輯同意了。但系統(tǒng)上的時間沒變,給編輯又寫郵件了,沒回復(fù)
10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢微 2026-03-15 | 3/150 |
|