[gmx-users] Hamiltonian replica exchange not working in 4.6

Sanku M msanku65 at yahoo.com
Wed Aug 7 00:52:59 CEST 2013


Dear Gromacs user,
  I am not sure whether this is any potential bug. But I found that in gromacs4.6.3, hamiltonian replica exchange simulation is not
working. The SAME simulation works perfectly fine in gromacs 4.5.4
But, when I try to run the simulation using gromacs4.6.3, I get following error:
"The properties of the 10 systems are all the same, there is nothing to exchange
For more information and tips for troubleshooting, please check the GROMACS"


Here is the details on the error:
   I have first defined A and B state by scaling the relevant parameter.. Then I have  generated different tpr file by using mdp file which only differ by 'init-lambda' parameter:
For example, I have 10 different mdp files: Free energy part of two of them are shown below:
for hremd0.mdp, the relevant part for free energy part is 
  ; Free energy control stuff
free_energy  = yes
init-lambda              = 0
delta_lambda = 0

for hremd1.mdp
; Free energy control stuff
free_energy  = yes
init-lambda              = .080
delta_lambda = 0
..
and so on..
The rest of mdp portions are same for all the mdp files...



However, with gromacs4.5.4, the same simulation works without problem. I do not get the error that all the systems are identical.

In both cases, I am using same system with OPLS forcefield.. 
The way I run the simulations is:
mpirun -np 10 mdrun_463 -multi 10 -v -s hremd -replex 250 -c &>log_repl


Any help will be appreciated.

Sanku



More information about the gromacs.org_gmx-users mailing list