| 24小時熱門版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 2237 | 回復(fù): 2 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
[求助]
計算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……同樣不知道怎么解決~ 請教大神,望各位能積極援助,不勝感激~~~ @老虎大王 @月只藍(lán) |

新蟲 (小有名氣)
木蟲 (著名寫手)
| 丟失原子了。主要檢查勢函數(shù)和初始構(gòu)型。先做一個簡單的合金構(gòu)型,比如B2,L12之類的,看看能不能正常運算,能量是否正確。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求化學(xué)調(diào)劑 +4 | wulanna 2026-03-28 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 070300求調(diào)劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
[考研] 085701環(huán)境工程,267求調(diào)劑 +16 | minht 2026-03-26 | 16/800 |
|
|
[考研] 292求調(diào)劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
[考研] 291求調(diào)劑 +7 | 孅華 2026-03-22 | 7/350 |
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
[考研] 材料292調(diào)劑 +12 | 橘頌思美人 2026-03-23 | 12/600 |
|
|
[考研] 322求調(diào)劑 +4 | 我真的很想學(xué)習(xí) 2026-03-23 | 4/200 |
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
[考研] 351求調(diào)劑 +4 | 麥克阿磊 2026-03-24 | 4/200 |
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
[考研] 329求調(diào)劑 +5 | 1() 2026-03-22 | 5/250 |
|
|
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +5 | 星沉uu 2026-03-26 | 6/300 |
|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +6 | cccchenso 2026-03-23 | 6/300 |
|
|
[考研] 材料調(diào)劑 +3 | iwinso 2026-03-23 | 3/150 |
|
|
[考研] B區(qū)考研調(diào)劑 +4 | yqdszhdap- 2026-03-22 | 5/250 |
|
|
[考研] 340求調(diào)劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|
|
[考研] 一志愿吉大化學(xué)322求調(diào)劑 +4 | 17501029541 2026-03-23 | 6/300 |
|
|
[考研] 341求調(diào)劑(一志愿湖南大學(xué)070300) +5 | 番茄頭--- 2026-03-22 | 6/300 |
|
|
[考研] 328求調(diào)劑 +4 | LHHL66 2026-03-23 | 4/200 |
|