			# REAX potential for Fe/O/H system
			# .....

			units real

			atom_style charge
			read_data data.FeOH3

			pair_style reax/c lmp_control
			pair_coeff * * ffield.reax.Fe_O_C_H 2 3 4

			neighbor 2 bin
			neigh_modify every 10 delay 0 check no

			fix 1 all nve
			fix 2 all qeq/reax 1 0.0 10.0 1e-6 param.qeq
			fix 3 all temp/berendsen 500.0 500.0 100.0

			timestep 0.25

			dump 1 all atom 30 dump.reax.feoh

			run 3000

lammps/examples/reax/Fe_O_H/in.FeOH3
d8fe7a91f045lammm-master

in.FeOH3
View Options