RE: energy jump between minimization and simulation

From: Bennion, Brian (Bennion1_at_llnl.gov)
Date: Wed Feb 16 2011 - 10:58:33 CST

Are you using binary restart files created from the minimization when you start your dynamics?

brian

________________________________________
From: owner-namd-l_at_ks.uiuc.edu [owner-namd-l_at_ks.uiuc.edu] On Behalf Of Eva Gonzalez Noya [evanoya_at_uchicago.edu]
Sent: Wednesday, February 16, 2011 8:14 AM
To: Ajasja Ljubetič
Cc: namd-l_at_ks.uiuc.edu
Subject: Re: namd-l: energy jump between minimization and simulation

Thank you very much for your answer.
Actually I tried many things already. I was doing first an NpT (using lagevin pistons) setting the initial velocity close to zero, and also tried NVE but the problem persists
Regads,
Eva

On 2/16/11 3:27 AM, Ajasja Ljubetič wrote:
Hi,
perhaps the energy jumps due to temperature and/or pressure control?
How are you heating the system?

Regards,
Ajasja

On Wed, Feb 16, 2011 at 01:06, Eva Gonzalez Noya <evanoya_at_uchicago.edu<mailto:evanoya_at_uchicago.edu>> wrote:
Hi:

I found that when running a MD simulation after a minimization the total energy suffers a jump when switching from the minimization to the simulation run (the conf file I am using is copied at the end of this message). In addition there two lines labeled with the same step number (the last step of minimization and the first step of the simulation) and the energy does not much for those neither (I am writing energies every step).
I am not rescaling velocities after the minimization. I am using a constraining harmonic potential on the alpha carbons of the protein both in the minimization and simulation runs (which I perform using a single conf file).
Anyone can help with this issue?
Thanks

This archive was generated by hypermail 2.1.6 : Mon Dec 31 2012 - 23:19:50 CST