From: Francesco Pietra (chiendarret_at_gmail.com)
Date: Thu Jul 12 2018 - 01:58:45 CDT
I was also perplexed at the performance degrading as the lambda changes,
which occurred soon or later, not always at the same lambda value, and to
the same extent when either one (ligand alone) or four nodes
(ligand+protein) are involved.
As I said, the rmsd is good, in particular the structure and pose of the
ligand (a polycyclic diterpenoid with a mobile side chain, rather exotic
structure) is well conserved during FEP. The ligand was parameterized
charmm36 with dih fitting at HF/6-31G* level and MD equilibration was
pretty long (>100ns) with absolutely flat rmsd/frame.
The only I can do (actually I am just doing that) is decreasing the number
of steps per lambda in order to keep the calculation within 70 hours (which
still requires a special permission at the cluster). Hopefully it will not
bring the calculation out of pseudo-convergence. Which occurred, as
expected, when I tried by decreasing the number of windows, while
increasing the number of steps per window.
Unfortunately there is little specific recent literature with namd/FEP for
complicated organic ligands. This is why I asked you about topogromacs to
compare with gromacs running charmm36. However, even the literature of FEP
with gromacs is limited to rather simple organic ligands and, what
surprised me very much, in accordance with experiments while the ligands
had been parameterized with gaff ff at semiempirical level. Probably I'll
see all these affairs with a different eye when my experience is ripe.
On Thu, Jul 12, 2018 at 1:23 AM Vermaas, Joshua <Joshua.Vermaas_at_nrel.gov>
> Colvars are indeed driven by a single CPU. Most of the colvars perform
> well if the number of atoms involved isn't too big, and bond lengths and
> angles are typical examples of that. But if you are asking for colvars that
> involve many atoms in a complicated relationship, performance isn't all
> that good. To me, the weird thing is that the performance degrades only as
> the lambda changes. Are you getting any absurd bonds as the trajectory
> On 2018-07-11 15:32:12-06:00 Francesco Pietra wrote:
> Thanks for your answer.
> 36 core Intel® Xeon® Broadwell/node, memory 115Gb/node, so that the
> problems are to look for elsewhere.
> 50555 atoms, including waters, whereby 4 nodes proved to be the best
> choice for MD
> where the performance was excellent.
> With the ligand alone in water the best choice proved to be one node.
> In retrospect, are colvars driven by a single CPU? Is that the problem? I
> could not set less colvars that I described in order
> to maintain the ligand in place.
> On Wed, Jul 11, 2018 at 7:21 PM Vermaas, Joshua <Joshua.Vermaas_at_nrel.gov>
>> What is the hardware on your cluster? FEP is not accelerated with GPUs.
>> Neither are colvars, which is I think where the problem may actually be.
>> How many atoms are in your colvar definitions?
>> On 2018-07-10 23:36:51-06:00 owner-namd-l_at_ks.uiuc.edu wrote:
>> I am observing a marked decrease in the performance of a NextScale
>> cluster running a FEP for protein-ligand, previously equilibrated for over
>> 100ns. No such problems when running MD equilibration on the same system--000000000000c057450570c7e5ac--
This archive was generated by hypermail 2.1.6 : Thu Sep 19 2019 - 23:19:49 CDT