[AMBER] accelerated molecular dynamics when iamd=1 or iamd=3

From: Sonia Ziada <sonia.ziada.gmail.com>
Date: Sun, 4 Jun 2017 17:39:21 +0200

Dear Amber users,



Concerning accelerated MD, I am not sure I get the point about the
conditions that must be verified to add the boost potential (deltaV(r)) and
the expression of the boost potential in the case where iamd=1 or 3.



Can you please tell me if the following is correct:



E_dih_boost = E_total_boost = E_boost = 0



If ( totdih_ene < EthreshD )

            # expression of the boost potential when iamd=2

            E_dih_boost = (EthreshD - totdih_ene )2 / (alphaD + EthreshD -
totdih_ene)



If ( iamd==1 or iamd==3 )

If ( tot_potenergy *+ E_dih_boost* < EthreshP ) * # correct ?*

totalenergy = tot_potenergy + *E_dih_boost*

E_total_boost = (EthreshP - totalenergy)2 / (alphaP + EthreshP -
totalenergy ) *# expression of the boost potential when iamd=1 ?*



If (iamd==3)

E_boost = E_total_boost + E_dih_boost *# expression of the boost potential
when iamd=3** ?*





Thanks in advance.


Best regards,
Sonia ZIADA
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sun Jun 04 2017 - 09:00:03 PDT
Custom Search