计算原子的温度显示到ovito(lammps)、如何将三斜盒子变为正交盒子
###############计算单个原子的温度#################
compute KE pp ke/atom
variable KB equal 8.625e-5 # 玻尔兹曼常数,eV/K, 这个数要根据计算unit来确定,(1)如果是real单位,能量单位为kcal/mol,则KB=1.987252e-3 Kcal/K;(2)如果是metal单位,能量单位为ev,则KB=8.6173713 ev/K;(3)百度百科的单位为KB=1.380658j/K
variable TEMP atom c_KE/1.5/${KB} # 单原子温度
dump 4 all custom ${interval} ${filename1}.lammpstrj id type x y z c_coord c_ke_atom c_pe_atom c_displace[1] c_displace[2] c_displace[3] c_displace[4] v_TEMP
##############lammps如何将三斜盒子变为正交盒子######################
在data文件中添加下面一行数据即可
0.0 0.0 0.0 xy xz yz
##############MS如何将斜晶胞改为正晶胞######################
Build——symmetry——Redefine Lattice——A:2 1 0