#############################################
#		Liquid Precursor Modelling of		#
#			Phenol-Formaldehyde				#
#		No. of monomers of Phenol=450		#
#		 No. of monomers of CH2O=470		#
#############################################

 clear
 echo both

 
#################################################
#	 			INITIALIZATION					#
#################################################
 read_restart restart\restart_1.restart
 kspace_style ewald 1e-6
 

#################################################
#	 			SIMULATION						#
#################################################
 compute 1 all rdf 100 1 1 5 5 1 5 5 1
 
#################################################
#	 			Equilibration					#
#################################################
 log log\nvt298_1.txt
 velocity all scale 298.0
 thermo 50
 thermo_style custom step temp press density vol pe ke etotal enthalpy
 fix 1 all nvt temp 298.0 298.0 100.0
 fix 2 all ave/time 100 1 100 c_1[*] file tmp.txt mode vector
 dump 1 all atom 5000 strj\nvt298_1.lammpstrj
 reset_timestep 0
 timestep 1
 run 200000
 unfix 1
 undump 1
 
 write_restart restart\nvt298_1.restart
 write_data data\nvt298_1.data