in.reaxc.tatb
d1bbe9546c02lammm-devel

in.reaxc.tatb
View Options

			# ReaxFF potential for TATB system
			# this run is equivalent to reax/in.reax.tatb,

			units real

			atom_style charge
			read_data data.tatb

			pair_style reax/c control.reax_c.tatb
			pair_coeff * * ffield.reax C H O N

			compute reax all pair reax/c

			variable eb equal c_reax[1]
			variable ea equal c_reax[2]
			variable elp equal c_reax[3]
			variable emol equal c_reax[4]
			variable ev equal c_reax[5]
			variable epen equal c_reax[6]
			variable ecoa equal c_reax[7]
			variable ehb equal c_reax[8]
			variable et equal c_reax[9]
			variable eco equal c_reax[10]
			variable ew equal c_reax[11]
			variable ep equal c_reax[12]
			variable efi equal c_reax[13]
			variable eqeq equal c_reax[14]

			neighbor 2.5 bin
			neigh_modify delay 0 every 5 check no

			fix 1 all nve
			fix 2 all qeq/reax 1 0.0 10.0 1.0e-6 reax/c
			fix 4 all reax/c/bonds 5 bonds.reaxc

			thermo 5
			thermo_style custom step temp epair etotal press &
			v_eb v_ea v_elp v_emol v_ev v_epen v_ecoa &
			v_ehb v_et v_eco v_ew v_ep v_efi v_eqeq

			timestep 0.0625

			#dump 1 all custom 100 dump.reaxc.tatb id type q x y z

			#dump 2 all image 5 image.*.jpg type type &
			# axes yes 0.8 0.02 view 60 -30
			#dump_modify 2 pad 3

			#dump 3 all movie 5 movie.mpg type type &
			# axes yes 0.8 0.02 view 60 -30
			#dump_modify 3 pad 3

			fix 3 all reax/c/species 1 5 5 species.tatb

			run 25