| 查看: 2197 | 回复: 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 |
回复此楼» 猜你喜欢
面上合作单位盖章
已经有8人回复
-
面上函评意见出来了，像什么等级？
已经有15人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
基金未中，这种答复是模板吗？
已经有7人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有7人回复
3楼2016-04-30 00:03:35
13楼2016-05-03 14:14:18
14楼2016-05-03 15:54:02
20
投票: