Hi,
From my experience, the problem is most certainly in your input files. I
would think there are steric clashes (overlapping spheres) that could
potentially cause this problem.
How are you inspecting your input files prior to minimization?
On Fri, Jun 2, 2023 at 10:27 AM Bhattacharjee, Sinjini via AMBER <
amber.ambermd.org> wrote:
> Hi, I am trying to minimise a protein-membrane system using sander and
> pmemd in Amber20. However, the calculation keep crashing in the first cycle
> with the following error:
>
>
> NSTEP ENERGY RMS GMAX NAME NUMBER
> 1 NaN NaN 1.0669E+24 C19 90066
>
> BOND = 31692.7352 ANGLE = 68711.0955 DIHED =
> 59597.8979
> VDWAALS = NaN EEL = -12389.2862 HBOND =
> 0.0000
> 1-4 VDW = NaN 1-4 EEL = -149672.0786 RESTRAINT =
> 0.0000
>
>
>
> This is my initial input:
>
> Minimize
> System minimization:
> &cntrl
> imin=1, ntmin=1, nmropt=0
> maxcyc=2000, ncyc=1000,
> ntx=1, irest=0,
> ntpr=100, ntwr=100, iwrap=0,
> ntf=1, ntb=1, cut=10.0, nsnb=20,
> igb=0,
> ibelly=0, ntr=1,
> restraintmask=“.CA", restraint_wt=20.0,
> &end
> /
>
>
> I added the membrane using packmol-memgen.
> The system size is huge and I checked and removed clashes in the protein
> and the non-standard residues.
> Can you please tell me what is the problem and how can I solve this? Is it
> an issue with the parameters or the structure?
>
>
> Many thanks,
> —
> Sinjini Bhattacharjee
>
> PhD Student
> Molecular Theory and Spectroscopy
> Max-Planck-Institut für Kohlenforschung
> 45470 Mülheim an der Ruhr, Germany
>
> email: bhattacharjee.kofo.mpg.de<mailto:bhattacharjee.kofo.mpg.de>
> ---
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri Jun 02 2023 - 08:00:04 PDT