[GaMD-discuss] GaMD Equilibration Time Query

[Christine Q. <chr...@uc...>]

Dear GaMD developers,

I am running GaMD on several protein (a monomer of 42 residues)-membrane systems with waters and ions. The systems have sizes of about 170000 atoms. This is the parameters that I'm using for my equilibration:
# GaMD Equilibration Params _________________________
accelMD                                         on
accelMDdual                                     on
accelMDdihe                                     on
accelMDG                                        on
accelMDGiE        1
accelMDGRestart  off

accelMDGcMDPrepSteps    200000
accelMDGEquiPrepSteps   200000
accelMDGStatWindow      20000

accelMDGcMDSteps        1000000
accelMDGEquiSteps       25000000
accelMDOutFreq          1000

accelMDGsigma0P 6.0
accelMDGsigma0D 6.0

run 26000000
I've noticed from the log file for all of my systems that the k0d and k0p values reach 1 really early (at about step 1210000) and then stay at 1 throughout. However, I'm afraid that it's a bit too fast. Is this behavior normal? Do my parameters look ok?

If everything looks alright, is it ok that I shorten the equilibration time a bit, or should I keep it at 50 ns?

Thank you so much for your help!

All the best,
Christine Qian