Hi,all
The below is my input file,I want to set the system’s Temperature as 300K,but it displays 158K,I do not know why.
dimension 3
boundary p p p
units metal
neighbor 2.0 bin
neigh_modify every 1 delay 0 check yes
atom_style atomic
lattice fcc 3.615
read_data data.min
group grain1 type 1
group grain2 type 2
pair_style eam/alloy
pair_coeff * * cu.eam Cu Cu
compute p all stress/atom
compute pe all pe/atom
dump 1 all custom 2000 dump.min tag type xs ys zs vx vy vz c_pe c_p[1] c_p[2] c_p[3]
thermo_style custom step pe ke temp enthalpy pxx pyy pzz pxy pyz pxz vol lx ly lz
thermo 100
log log.min
velocity all create 300.0 3456781 units box
fix 1 all npt 300 300 100 xyz 0 0 1000
######fix 2 all temp/rescale 10 300 300 10 1.0
timestep 0.002
run 5000
unfix 1
fix 2 all npt 0.1 0.1 100 xyz 0 0 1000 drag 0.2
run 5000