clear
units		lj
dimension	2
boundary        p p p
atom_style	  sphere 
                                                                               

region		mybox block 0 85 0 85 -0.01 0.01
create_box	1 mybox
create_atoms   1 random 5 12345 NULL overlap 0.4 maxtry 50 units box
group		mobile type 1
set		group mobile diameter 1.0
set		group mobile mass  1.0


pair_style	lj/cut 2.5
pair_coeff	1 1 1.0 0.35 0.3928
pair_coeff	2 2 1.0 0.35 0.3928
pair_coeff	1 2 1.0 0.35 0.3928


pair_modify shift yes 
neighbor	0.2 bin
neigh_modify	delay 5 every 1 check yes

fix		    1 all nve/sphere 
fix             4 mobile langevin 1.0 1.0 1.0 437624 omega yes scale 1 3.3 
fix             active1 mobile propel/self velocity 0.5
fix             fric all viscous 1.0
fix		 2 all enforce2d


timestep	0.002

dump DUMP all custom 50 dump.lammpstrj id type x y z 


           
run		1000000