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