**From:** Brian Radak (*brian.radak_at_gmail.com*)

**Date:** Thu Mar 15 2018 - 08:24:07 CDT

**Next message:**Giacomo Fiorin: "Re: statistical error in ABF calculations"**Previous message:**Srijita Paul: "(no subject)"**In reply to:**Srijita Paul: "(no subject)"**Next in thread:**Artur Hermano: "Re:"**Reply:**Artur Hermano: "Re:"**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ] [ attachment ]

I think you are looking for the *mean* acceptance probability, that is, the

average probability with which an exchange is accepted. You can compute

this in two different ways which, to my knowledge, are essentially

equivalent for any reasonably large number of REMD cycles:

1) take the expectation of the Metropolis criterion P = min[1,

exp(-Delta_ij)] where Delta_ij contains the observed energies and

temperatures -- this method is cumbersome and requires sifting through a

lot of data

2) just divide the number of success by the number of attempts = (50 / 450)

= 0.11 = 11%

The reasonableness of the observed value depends on what your exchange

scheme is. I assume that your script implements a nearest neighbor

sampling? If the neighbors are chosen stochastically with equal

probabilities, then 11% is very close to optimal. If the neighbors are

chosen in the deterministic "up/down" strategy then something closer to 20%

is preferred.

If you are unhappy with your acceptance probability you have two options:

1) assume that you have bad statistics and keep running until the

performance numbers change or 2) start over and choose more replicas over

the temperature range.

HTH,

BKR

On Thu, Mar 15, 2018 at 7:44 AM, Srijita Paul <srijitap91_at_gmail.com> wrote:

*> Hi,
*

*>
*

*> Can anybody explain me the output file obtained from a remd simulation
*

*> .job0.restart900.0.tcl.
*

*>
*

*> array set replica {index.b 4 index 3 temperature.a 283.48
*

*> exchanges_attempted 450 loc.a 4 temperature.b 289.07 temperature 286.26
*

*> exchanges_accepted 50 loc.b 18 index.a 2}
*

*>
*

*> exchanges_attempted 450
*

*> exchanges_accepted 50
*

*>
*

*> Is it a good result for remd? How can I find acceptance probability for my
*

*> system?
*

*>
*

**Next message:**Giacomo Fiorin: "Re: statistical error in ABF calculations"**Previous message:**Srijita Paul: "(no subject)"**In reply to:**Srijita Paul: "(no subject)"**Next in thread:**Artur Hermano: "Re:"**Reply:**Artur Hermano: "Re:"**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ] [ attachment ]

*
This archive was generated by hypermail 2.1.6
: Tue Dec 31 2019 - 23:19:45 CST
*