############          simulation model and settings        ############
units real 
dimension 3
boundary p p p
atom_style full
box tilt   large
timestep 0.25
############          force field setting        ############
bond_style      harmonic
angle_style     cosine/periodic 
dihedral_style  harmonic
improper_style  fourier
pair_style      lj/cut/coul/long 12.5   
pair_modify     tail yes mix geometric 

special_bonds   lj/coul 0.0 0.0 0.5

kspace_style    ewald 1.0e-6
kspace_modify	gewald 0.5

read_data data.data

neighbor        2.0 bin
neigh_modify    every 1 delay 0 check yes

############          energy minimization        ############
reset_timestep  0

velocity	all create 0.0 1
thermo		50

min_style       cg
min_modify      line quadratic dmax 0.001
minimize        1e-10 1e-10 10000 1000000

write_data opt.data