From: Ritu Arora (ritu.arora_at_mail.concordia.ca)
Date: Thu Sep 17 2020 - 10:54:33 CDT
Hi Giacomo,
I second Hemanth's advice and this was the first thing I did- checking the trajectory and the restart files. As I mentioned before, right after the first 2ns of simulation time something went wrong. Since the system becomes unstable in the very beginning, I am assuming that the boost is quite high.
These numbers were probably appropriate for the system shown in the tutorial, but are they for your system?
Yes, of course I calculated these parameters for my system as explained here https://www.ks.uiuc.edu/Research/namd/mailing_list/namd-l.2012-2013/1460.html.
I will do trial and error method to choose better aMD parameters and will get back to you all.
Thanks,
Ritu
________________________________
From: Giacomo Fiorin <giacomo.fiorin_at_gmail.com>
Sent: Wednesday, September 16, 2020 5:25 PM
To: NAMD list <namd-l_at_ks.uiuc.edu>; Ritu Arora <ritu.arora_at_mail.concordia.ca>
Cc: HEMANTH H <hemanth.h_at_iitgn.ac.in>
Subject: Re: namd-l: NAN in output of aMD
Hello Ritu, Hemanth's advice is a common and quite sensible one. Check the trajectory or the latest checkpoint (i.e. restart) file prior to the instability occurring to see what is going on.
You also mention using parameters from the AMD tutorial, but have you checked what do they mean in the user's guide?
https://www.ks.uiuc.edu/Research/namd/2.14/ug/node69.html
These numbers were probably appropriate for the system shown in the tutorial, but are they for your system?
Giacomo
On Tue, Sep 15, 2020 at 2:28 PM Ritu Arora <ritu.arora_at_mail.concordia.ca<mailto:ritu.arora_at_mail.concordia.ca>> wrote:
Hi Hemanth,
Thanks for sharing your experience. The given output was just a portion of the log file, I start getting 'nan' energy values right after 2 ns of simulation time. There is no issue with the xsc file; I believe that something is wrong with my input in the configuration file. I am not sure if I am applying the proper boost or if something is wrong with the way I am using dual boost. I calculated accelMDE and accelMDalpha as reported in the tutorials, and considering the first 50 ns conventional MD simulations.
#aMD
accelMD on;
accelMDdihe on;
accelMDE 8581; #followed NAMD tutorial to calculate accelMDE and accelMDalpha
accelMDalpha 533;
accelMDOutFreq 5000;
Any further help from NAMD users or developers will be appreciated.
Thanks and regards,
Ritu
________________________________
From: HEMANTH H <hemanth.h_at_iitgn.ac.in<mailto:hemanth.h_at_iitgn.ac.in>>
Sent: Tuesday, September 15, 2020 1:55 AM
To: namd-l_at_ks.uiuc.edu<mailto:namd-l_at_ks.uiuc.edu> <namd-l_at_ks.uiuc.edu<mailto:namd-l_at_ks.uiuc.edu>>; Ritu Arora <ritu.arora_at_mail.concordia.ca<mailto:ritu.arora_at_mail.concordia.ca>>
Subject: Re: namd-l: NAN in output of aMD
Ms. Arora,
1. I assume that your simulation had worked fine till the step 18050000. In that case, the best way would be to visually check your trajectory to see what is happening there. I had run into something similar before, but my issue was related to the xsc file being not parsed properly. Check your last checkpoint file and you should be able to figure it out.
Regards
Hemanth
On Mon, Sep 14, 2020 at 8:06 PM Ritu Arora <ritu.arora_at_mail.concordia.ca<mailto:ritu.arora_at_mail.concordia.ca>> wrote:
Hello All,
I am trying accelerated MD simulations on my protein system. To do so I added the following commands in the configuration file of production run:
#aMD
accelMD on;
accelMDdihe on;
accelMDE 8581; #followed NAMD tutorial to calculate accelMDE and accelMDalpha
accelMDalpha 533;
accelMDOutFreq 5000;
Unfortunately, my output file contains nan:
TIMING: 18050000 CPU: 253081, 0.0139727/step Wall: 259046, 0.0141721/step, 27.36 hours remaining, 1189.199219 MB of memory in use.
ETITLE: TS BOND ANGLE DIHED IMPRP ELECT VDW BOUNDARY MISC KINETIC TOTAL TEMP POTENTIAL TOTAL3 TEMPAVG PRESSURE GPRESSURE VOLUME PRESSAVG GPRESSAVG
ENERGY: 18050000 nan nan -nan -nan nan nan 0.0000 0.0000 nan nan nan nan nan nan nan nan nan nan nan
I am not sure how-to solve this problem; can anyone help me with this regard? I am sorry if I have re-posted this message.
Thank you.
Best,
Ritu
-- Hemanth H 18310019 Research Scholar IIT Gandhinagar -- Giacomo Fiorin Associate Professor of Research, Temple University, Philadelphia, PA Research collaborator, National Institutes of Health, Bethesda, MD http://goo.gl/Q3TBQU https://github.com/giacomofiorin
This archive was generated by hypermail 2.1.6 : Fri Dec 31 2021 - 23:17:09 CST