lammps运行 温度无法控制
在模拟氧化铁烧结时 温度无法控制,in文件如下units metal
atom_style full
read_data fe2o316.data
velocity all create 300 87287
group Fe type 1
group O type 2
set group Fe charge 1.2
set group O charge -0.8
pair_style hybrid/overlay buck 5.0 coul/long 5.0
pair_coeff 1 2 buck 1414.6 0.3128 0
pair_coeff 2 2 buck 9547.96 0.2192 32
pair_coeff * * coul/long
kspace_style ewald 1.0e-4
timestep 0.0001
thermo 100
neighbor 0.5 bin
neigh_modify every 5 delay 0 check yes
fix 1 all npt temp 300 300 2.0 iso 0 0 5.0
run 200
unfix 1
fix 1 all nvt temp 300 1800 2.0 drag 2.0
dump 1 all xyz 100 sinterFe2O3.xyz
run 2000
运行后出现了 温度无法控制的情况 程序结果如下图 往前辈们指教
返回小木虫查看更多
fix 1 all npt temp 300 300 0.01试试。2.0有点大
梢园巡绞┐蠹甘叮瑃imestep也可以扩大十倍看看(要小心可能这一步改了会运行不了,那样的话只能还原),你的系统还远没有平衡,所以说不准,1W多K而已,其实结构也许并没有那么大的变化,再跑几万步说不定又好了,所以是你的步数太少,加上步长过短,目前还看不出来。如果这样操作以后还是不能降到指定温度,就得考虑你的势函数单位和数值是否正确了
,