##########################################################
#  Units, BC, Forcefield Parameters and Initialization   #
##########################################################
units           real
boundary	    p p p

atom_style      full #sphere
bond_style      harmonic
angle_style     charmm
dihedral_style  charmm
improper_style  harmonic
pair_style      hybrid lj/charmm/coul/long 10 12 colloid 100.0
pair_modify     mix arithmetic

kspace_style    pppm 1e-4
read_data       ./data_read/01.data

special_bonds   charmm
neighbor	2.0 multi
neigh_modify	delay 0 every 1 check yes one 100000 page 1000000
comm_modify	mode multi

pair_coeff	1  11  colloid 57.86  3.47 100.0    0.0 52.3
pair_coeff	2  11  colloid 105.21 3.47 100.0    0.0 63.3
pair_coeff	3  11  colloid 84.45  3.47 100.0    0.0 62.8		
pair_coeff	4  11  colloid 48.26  3.47 100.0    0.0 64.96
pair_coeff	5  11  colloid 63.84  3.47 100.0    0.0 65.0
pair_coeff	6  11  colloid 75.34  3.47 100.0    0.0 65.3
pair_coeff	7  11  colloid 50.47  3.47 100.0    0.0 60.3
pair_coeff	8  11  colloid 41.79  3.47 100.0    0.0 60.3
pair_coeff	9  11  colloid 93.45  3.47 100.0    0.0 62.8
pair_coeff	10 11  colloid 85.31  3.47 100.0    0.0 62.5
#########################################################
#                group materials                        #
#########################################################
group		pmma type  3 4 5 6 7 8 9 10
group		water type 1 2
group		colloid type 11
###################################################################
#		Wall group definition for anchring treatment	          #
###################################################################
region          anchored block INF INF INF INF -2.5 2.5 units box
group           anchoredH_gross region anchored
group 			anchoredH_net subtract anchoredH_gross water

group			pgroup subtract all anchoredH_net			
group			dilation subtract all anchoredH_net			

group			mobile_pmma subtract pmma anchoredH_net
########################################################
#        Make center slab of PMMA immobile              #
#########################################################
velocity	anchoredH_net zero angular
velocity	anchoredH_net zero linear
neigh_modify	exclude group anchoredH_net anchoredH_net
delete_bonds    anchoredH_net multi
########################################################
fix		1 pgroup shake 1e-6 500 0 m 1.0 a 15
#########################################################
compute         t0 pgroup temp

fix             a1 water nvt temp 300.0 300.0 200.0
fix             a2 mobile_pmma nvt temp 300.0 300.0 200.0
fix             b4 colloid nve

compute 	t1 all pressure t0
#########################################################
#                      Outputs                          #        
#########################################################
thermo          100
thermo_style    custom step c_a1_temp c_a2_temp c_t1[3] press
########################################################
timestep	2.0
run		    10000