Dear community,   
   
I am currently running REMD with reaction field. I started with the berendsen 
thermostat, and my md average temperatures were 10 to 15K warmer than the ones 
specified in the input file. I had thermostat time coupling constants of 0.1 
for the protein and 0.4 for the solvent, that I did not want to change. I tried 
playing with the cutoff (originally set to 1.4nm) and frequence of neighbour 
list update, and finally, with the former set to 2.5nm, had satisfying 
results.. Changing for the Nose-Hoover thermostat, I obtained satisfying 
results without having to increase the cutoff.. Is there any explanation for 
that?  
   
Thank you very much for your help!! And many thanks to Berk and David, for all 
their help on that coul-LR term !! 
   
Pascal   
   
   
******************************************************************************* 
  
Pascal Baillod (PhD student)    
******************************************************************************* 
  
Swiss Federal Institute of Technology EPFL              Tel: +41-(0)21-693-0322 
  
Institute of Chemical Sciences and Engineering ,        Fax: +41-(0)21-693-0320 
  
Laboratory of Computational Chemistry and Biochemistry  [EMAIL PROTECTED]   
Room BCH 4121, Avenue Forel,                            http://lcbcpc21.epfl.ch 
  
CH-1015 Lausanne           
******************************************************************************* 
  
   
_______________________________________________
gmx-users mailing list    gmx-users@gromacs.org
http://www.gromacs.org/mailman/listinfo/gmx-users
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [EMAIL PROTECTED]
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

Reply via email to