| 查看: 2264 | 回复: 17 | |||
[交流]
用LAMMPS模拟Ni3AL降温，体积和能量没有突变
|
|||
|
我想在2500K到300K降温，每100K降温10W步，弛豫10W步，一共450W步。 但是第一次计算的时候发现体系在2500到2400K之间体积突变，明显不对，又换了个势函数，这下直接没有突变了。 附上in文件 # LAMMPS Ni3Al units metal boundary p p p atom_style atomic read_data Ni3Al1000.dat pair_style eam/alloy pair_coeff * * Al-Ni2009.eam.alloy Ni Al timestep 0.001 thermo 50 neighbor 0.6 bin neigh_modify every 20 delay 0 check yes variable N equal step compute 3 all pe/atom compute 4 all ke/atom compute 5 all coord/atom 4.0 velocity all create 2500 825577 dist gaussian dump 1 all cfg 20000 a.*.cfg id type xs ys zs c_3 c_4 c_5 dump 2 all atom 20000 dump.2500 fix 1 all npt temp 2500 2500 0.02 iso 0 0 0.2 drag 1 run 100000 unfix 1 write_restart restart.$N undump 1 undump 2 label loop variable x index 2500 2400 2300 2200 2100 2000 1900 1800 1700 1600 1500 1400 1300 1200 1100 1000 900 800 700 600 500 400 variable y index 2400 2300 2200 2100 2000 1900 1800 1700 1600 1500 1400 1300 1200 1100 1000 900 800 700 600 500 400 300 dump 1 all cfg 5000 a$y.*.cfg id type xs ys zs c_3 c_4 c_5 dump 2 all atom 5000 dump.$y # Cooling steps: fix 1 all npt temp $x $y 0.01 iso 0 0 0.2 drag 1 run 100000 unfix 1 write_restart restart.$y.$N # Equilibrium steps: fix 2 all npt temp $y $y 0.02 iso 0 0 0.2 drag 1 run 100000 unfix 2 undump 1 undump 2 write_restart restart.$N.$y next x next y jump Ni3Al.in loop |
回复此楼» 猜你喜欢
【长期招聘】复旦大学 AI 制药实验室 博士后招聘
已经有11人回复
-
2027博士申请-临床检验诊断学、病原生物学、微生物方向
已经有5人回复
-
各位基础老师们，看完中央《从严治党》的决定后，对学术界更加有信心了吧
已经有4人回复
3楼2016-04-30 00:03:35
13楼2016-05-03 14:14:18
14楼2016-05-03 15:54:02
15楼2016-05-04 10:01:34
16楼2016-05-04 10:02:13
17楼2016-05-04 10:04:38
18楼2017-05-15 19:44:41
简单回复
2016-04-29 23:44
回复
潇潇雨歇99(金币+1): 谢谢参与
2016-04-30 00:09
回复
潇潇雨歇99(金币+1): 谢谢参与
xhmaohan5楼
2016-04-30 02:04
回复
潇潇雨歇99(金币+1): 谢谢参与
2016-04-30 05:44
回复
潇潇雨歇99(金币+1): 谢谢参与
2016-04-30 07:23
回复
潇潇雨歇99(金币+1): 谢谢参与
2016-04-30 07:28
回复
潇潇雨歇99(金币+1): 谢谢参与
2016-04-30 07:33
回复
潇潇雨歇99(金币+1): 谢谢参与
yangdw12345710楼
2016-04-30 07:40
回复
潇潇雨歇99(金币+1): 谢谢参与
yangdw12345811楼
2016-04-30 07:45
回复
潇潇雨歇99(金币+1): 谢谢参与
yangdw12345912楼
2016-04-30 07:50
回复