# Name of the script - Example

units		metal
atom_style	atomic
boundary	p p p


#--------Reading Configuration----------

read_data     alpha8108Silica5.txt

#--Masses--##

mass 1 28.0855
mass 2 15.9994

#------interaction coefficients e,sigma,cutouff----------


#pair_style tersoff
#pair_coeff * * SiO.tersoff  Si O 

pair_style buck/coul/cut 10.0  
pair_style buck/coul/cut 10.0 8.0 
pair_coeff 1 1 1282.9 0.321 10.61 10.0
pair_coeff 1 1 1282.9 0.321 10.61 10.0 8.0

#----initial velocity-----

velocity all create 5000 14 dist gaussian

#-------neighbor list rebuild parameters----------

timestep        0.002    
neighbor        0.6 bin                             					
neigh_modify every 1 delay 0 check yes

#-------integration method(temp Tstart Tstop Tdamp)---------------

fix 1 all langevin 5000 300 0.125 14
#fix 1 all langevin 5000 270 0.02 10
fix 2 all nve
#fix 3 all viscous 0.1

#minimize 1.0e-7 1.0e-5 100 100
#min_style cg

#------output files-------------------------------

##dump        Saketh all custom 20 dump.movie id type x y z vx vy vz
dump    SilicaDump all custom 20 dump.starting id type x y z vx vy vz
dump SilicaDump3 all xyz 100 dump.starting1.xyz



thermo_style custom step temp press pe etotal lx ly lz xlo xhi ylo yhi zlo zhi
thermo  100

#run 100000
run 200000