lammps在轨迹中，不同单位下在轨迹中输出原子温度
・温度
metal单位时
compute KE all ke/atom variable KB equal 8.625e-5 variable TEMP atom c_KE/v_KB fix properties all ave/atom 50 60 3000 v_TEMP dump 1 all custom 3000 dump.lammpstrj id type x y z f_properties[1]
real单位时
compute KE all ke/atom variable KB equal 1.9872041e-3 variable TEMP atom c_KE/v_KB fix properties all ave/atom 50 60 3000 v_TEMP dump 1 all custom 3000 dump.lammpstrj id type x y z f_properties[1]
应力
compute str all stress/atom NULL compute p all voronoi/atom variable stressx atom c_str[1]/10000 variable stressy atom c_str[2]/10000 variable stressz atom c_str[3]/10000 variable stressxy atom c_str[4]/10000 variable stressyz atom c_str[5]/10000 variable stressxz atom c_str[6]/10000 variable stressHY atom -(v_stressx/c_p[1]+v_stressy/c_p[1]+v_stressz/c_p[1])/3 variable stressVon atom sqrt((v_stressx/c_p[1]-v_stressy/c_p[1])^2+(v_stressy/c_p[1]-v_stressz/c_p[1])^2+(v_stressz/c_p[1]-v_stressx/c_p[1])^2+6*((v_stressxy/c_p[1])^2+(v_stressyz/c_p[1])^2+(v_stressxz/c_p[1])^2))/sqrt(2) fix properties all ave/atom 50 60 3000 v_stressHY v_stressVon dump 1 all custom 3000 dump.lammpstrj id type x y z f_properties[1] f_properties[2]
