#####################  LAMMPS INPUT ##########################
#####################  AUTHOR -> MANASH  ###################################
atom_style	full
units	real
dimension	3
boundary	p p p
neighbor	2.0 bin
neigh_modify	every 1 delay 0 check yes

####################   ATOM TYPES   ########################################
#atom type numbering
variable Artype equal	1

##############################  LJ PARAMETERS  #############################

# Define interaction parameters
pair_style lj/cut 12.0
pair_modify mix arithmetic

#############################  READ GEOMETRY ##############################

read_data	geometry.data


###################### PAIR Coefficients  ############################################
pair_coeff ${Artype} ${Artype} 0.238 3.405 # O-O

############################  BASIC ATOM  GROUPS ###################################

#groups
group ar        type ${Artype}

########################## THERMO STYLES  #####################

thermo 100
thermo_style custom step temp press density pe etotal

fix nvtall all nvt temp 120 120 100

compute pair_atoms ar property/local patom1 patom2
compute force_dist ar pair/local dist force
dump pair ar local 1000 dump.pair.* c_pair_atoms[1] c_pair_atoms[2] c_force_dist[1] c_force_dist[2]
dump atoms ar custom 1000 dump.full.* id x y z

run 50000