# input file for SiO2

units real
dimension 3
processors  * * *
boundary p p p


#read data
atom_style charge
read_data Sio2.txt

#potential
include pot.mod

# outputs
thermo 100
Thermo_style custom step temp pe press vol lx density
dump	2 all custom md.lammpstrj id type x y x 

# Initial minimization
minimize 1.0e-10 1.0e-10 100000 10000

# Initial mixing
fix 1 all nvt temp 5000 5000 100
run 200000
unfix 1

fix 1 all npt temp 5000 5000 100 iso 0 0 1000
run 100000
unfix 1

# cooling 
fix 1 all npt temp 5000 300 100 iso 0 0 1000
run 4700000
unfix 1

# final relaxation
fix 1 all npt temp 300 3000 100 iso 0 0 1000
run 100000
unfix 1

Write_data 5300.dat
Write_restart 5300K.rest