相同的材料构型,在不同的系综下产生如下不同的结果:
1.NVE系综下:
2.NVT系综如下:
产生的不同结果如上,原脚本如下:
1.
shell mkdir initial-comp
shell mkdir config-comp
shell mkdir restart-comp
#--------------------------DEFINITION-------------------
units metal
dimension 3
boundary s s s
atom_style atomic
read_data Ti.data
region 1 block INF INF INF INF INF 27 units box
region 2 block INF INF INF INF 319 INF units box
group bottom region 1
group upper region 2
group mobile subtract all bottom upper
group boundary union bottom upper
#---------- Define Interatomic Potential ---------------------
pair_style eam/fs
pair_coeff * * Ti.eam.fs Ti
#-----------------------------setting------------------------