title               =  slow heating up of a simulation system including protein solvent and ions
cpp                 =  /lib/cpp	  				


; run parameters
integrator          =  md		    			
tinit               =  0		    			
dt                  =  0.001					
nsteps              = 1000000   				


; allocation of the total simulation time to the temperature intervals
annealing           = single 		single			
annealing-npoints   = 4 		4 			 
annealing-temp      = 1 100 293 1	1 100 293 1  	                
annealing-time      = 0 100 1000 1100 	0 100 1000 1100 		


; physical parameters

; pressure coupling
pcoupl              = no       

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


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


; 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