46 lines
1.4 KiB
Plaintext
46 lines
1.4 KiB
Plaintext
newton on
|
|
log cnt.log
|
|
echo both
|
|
units metal
|
|
lattice sc 1.0
|
|
boundary p p fs
|
|
neighbor 1.0 bin
|
|
neigh_modify every 5 delay 0 check yes
|
|
|
|
atom_style mesont
|
|
# cut, path, BendingMode, TPMType
|
|
pair_style mesont/tpm 25.0 ../../../potentials/CNT_10_10 0 0
|
|
read_data system.init
|
|
pair_coeff * *
|
|
|
|
velocity all create 600.0 2019
|
|
timestep 0.010
|
|
fix 1 all nve
|
|
#fix 1 all nvt temp 300.0 300.0 1.0
|
|
thermo_modify flush yes
|
|
thermo 10
|
|
reset_timestep 0
|
|
|
|
compute Es all mesont/Es
|
|
compute Eb all mesont/Eb
|
|
compute Et all mesont/Et
|
|
compute Ek all ke/atom
|
|
compute Es_tot all mesont/Es_tot
|
|
compute Eb_tot all mesont/Eb_tot
|
|
compute Et_tot all mesont/Et_tot
|
|
compute Ep_tot all pe
|
|
compute Ek_tot all ke
|
|
variable time_ equal time
|
|
variable Ep_ equal c_Ep_tot
|
|
variable Ek_ equal c_Ek_tot
|
|
variable Etot_ equal v_Ek_+v_Ep_
|
|
variable Es_ equal c_Es_tot
|
|
variable Eb_ equal c_Eb_tot
|
|
variable Et_ equal c_Et_tot
|
|
|
|
dump out_dump all custom 50 config_E.dump id type x y z c_Es c_Eb c_Et c_Ek ix iy iz
|
|
fix out_info all print 10 "${time_} ${Etot_} ${Ek_} ${Ep_} ${Es_} ${Eb_} ${Et_}" file "E.txt" screen no
|
|
|
|
run 50
|
|
write_data system_E.data
|