[gmx-users] Continuing inturepted REMD simulation

Nizar Masbukhin nizar.fkub08 at gmail.com
Tue Dec 2 06:59:37 CET 2014


Dear gromacs users,
I have finished an REMD simulation, tough it had interupted (by electrical
instabilty that lead computer shutdown) several times.  I continue the
interupted simulation using this command:
mpirun -np 8 mdrun_mpi -v -multidir md[01234567] -replex 2000 -cpt 3 -cpi
state.cpt

however, the remd statistics show unexpected exchange probability.
Replica exchange statistics
Repl  3365 attempts, 1683 odd, 1682 even
Repl  average probabilities:
Repl     0    1    2    3    4    5    6    7
Repl
Repl  number of exchanges:
Repl     0    1    2    3    4    5    6    7
Repl        0    0    0    0    0    0    0
Repl  average number of exchanges:
Repl     0    1    2    3    4    5    6    7
Repl      .00  .00  .00  .00  .00  .00  .00


Repl                                Empirical Transition Matrix
Repl       1       2       3       4       5       6       7       8
Repl  1.0000  0.0000  0.0000  0.0000  0.0000  0.0000  0.0000  0.0000  0
Repl  0.0000  1.0000  0.0000  0.0000  0.0000  0.0000  0.0000  0.0000  1
Repl  0.0000  0.0000  1.0000  0.0000  0.0000  0.0000  0.0000  0.0000  2
Repl  0.0000  0.0000  0.0000  1.0000  0.0000  0.0000  0.0000  0.0000  3
Repl  0.0000  0.0000  0.0000  0.0000  1.0000  0.0000  0.0000  0.0000  4
Repl  0.0000  0.0000  0.0000  0.0000  0.0000  1.0000  0.0000  0.0000  5
Repl  0.0000  0.0000  0.0000  0.0000  0.0000  0.0000  1.0000  0.0000  6
Repl  0.0000  0.0000  0.0000  0.0000  0.0000  0.0000  0.0000  1.0000  7

That statistics was very different when i try remd simulation with short
timesteps, and uninterupted one. the exchange probabilty was 0.3 as i
expected.
Could anyone tell me whats wrong?
--


More information about the gromacs.org_gmx-users mailing list