title				= nVT equilibration of a simulation system including protein solvent and ions
cpp                 =  /lib/cpp	  				


; run parameters
integrator  =  md		
tinit       =  0		
dt		    = 0.002		
nsteps		= 500000	


; initial velocity generation
gen_vel		= yes		
gen_temp	= 298		
gen_seed	= -1		


; physical parameters

; pressure coupling
pcoupl              = no       

; temperature coupling
tcoupl              = V-rescale	                    
ref_t               = 298.0  	298.0    	    
tau_t               = 0.1    	0.1       	    
tc_grps             = Protein_Het_vedani_wat  SOL_ION	    
nsttcouple	        = 1               		    


; output parameter   	
nstxout-compressed  	=  2500		
compressed-x-precision 	=  1000		 	
nstenergy           	=  2500		  	
nstxout             	=  500000 		
nstvout             	=  500000		
nstfout             	=  500000		
nstlog              	=  500000		


; neighbor search
nstlist             	=  10			
cutoff-scheme       	=  verlet	 	
verlet-buffer-tolerance =  0.005		
ns_type             	=  grid			
pbc                 	=  xyz	        	
rcoulomb            	=  1.1          	
rvdw                	=  1.1          	
vdw-modifier            = Potential-switch	
rvdw-switch		        = 1.05			


; electrostatics 
coulombtype         =  PME     
pme_order           =  4       
fourierspacing      =  0.12    
ewald_rtol          =  1e-5    
vdwtype             =  cut-off 


; bond parameters
continuation			= no			
constraints             = h-bonds  		
constraint_algorithm    = lincs 		
lincs_iter              = 1          	
lincs_order             = 4          	


; COM motion removal
nstcomm   			= 1000     	
comm_mode 			= Linear    
comm_grps 			= System  	
refcoord_scaling 	= com