echo            both
 units           real
 boundary        p p p
 
 timestep        1
 atom_style      full


 read_data       wetting.data

#group atoms to molecules
 group    O type 2
 group    H type 1
 group    carbon type 3
 group    water type 1 2

 group    drop molecule 1 2184
 
# water model: spc/e
 set  type 1 charge 0.4238
 set  type 2 charge -0.8476


 pair_style	lj/cut/coul/long 12.0

 pair_coeff     3  3     0.01544    3.214    #C-C
 pair_coeff     1  1     0.0     0.0     #H-H   
 pair_coeff     2  2    0.1553    3.166    #O-O  
	
	
 pair_coeff     1  3      0.0     1.607     #H-C
 pair_coeff     2  3      0.04896     3.19     #O-C
 pair_coeff     1  2    0.0     0.0     #H-O    
   
   
 bond_style    	harmonic              
 angle_style   	harmonic
 bond_coeff    	1 1000.00 1.000  # H-O
 angle_coeff   	1 100.0 109.47   #H-O-H

 neighbor	      2 bin
 neigh_modify     delay 0 every 1 check yes


velocity carbon set 0 0 0 units box

velocity water set 0.0 0.0 0.0
fix nomove water momentum 1 linear 1 1 0 angular

kspace_style pppm 1.0e-4
#########################################
#Fixing water molecules via a harmonic potential:
region sph1 sphere 50.0 50.0 28.0 27.0 
fix wall1 drop wall/region sph1 harmonic 1.0e8 0.0 1.0 

#Define a soft substrate:
region subwall plane 0.0 0.0 0.0 0.0 0.0 1.0 side in
fix wall2 drop wall/region subwall harmonic 6.66e-7 3.1 6.0 

#Fixing carbon atoms via a harmonic potential:
fix wall3 carbon spring/self 90000.0 xy
fix wall5 carbon spring/self 9.8122e-7 z   



#####################################
minimize 1.0e-4 1.0e-6 1000 10000

fix 1 water shake 1.0e-4 20 10000 b 1 a 1


#pair_modify tail yes

fix NVT1 water nvt temp 300.0 300.0 100.0  #Tdamp 100


compute com water com

 reset_timestep 0
 thermo_style custom step c_NVT1_temp pe ke c_com[*]
 thermo 1000 
 log logequi.txt 
 
 dump 1 all custom 10000 dump1.lammpstrj id element type xu yu zu
 dump_modify 1 element C O H
 
 run 50000
 undump 1
 reset_timestep 0
#######################################################
unfix   nomove
unfix wall1

velocity water set 0.0 0.0 NULL


 thermo_style custom step c_NVT1_temp pe ke c_com[*]
 thermo 1000 
 
 dump 2 all custom 10000 dump2.lammpstrj id type xu yu zu
 
 run 500000
 undump 2
 reset_timestep 0