Hi Annika,
Well that doesn't look good... It's hard to be able to tell what is wrong
though from what you list. I.e. which force field are you using, how did
you build your prmtop and inpcrd file? Are the PBradii set correctly for
IGB=8 etc.
I notice you don't specify a cutoff - this means, I think, that a default
cutoff of 8 angstroms will be used. This is too short for GB runs. You
should set this to something more like 16 or 18 or even 99999 to give you
no cut off. This could be your problem. This will also slow down the
simulations massively. Unless you have a specific reason to be running in
implicit solvent I would suggest trying this in explicit solvent.
All the best
Ross
On 8/7/14, 10:16 PM, "Annika Wittenbecher" <annika.wittenbecher.gmail.com>
wrote:
>Dear amber users,
>
>I have a problem with my molecule. It's unfolding a lot after 6 stages of
>heating to 300K.
>
>Here the input file of Stage 6:
>
>Stage 6 heating of ArsM 250 to
>300K
> &cntrl
>
> imin=0, irest=1,
>ntx=5,
> nstlim=10000,
>dt=0.0005,
> ntc=2,
>ntf=2,
> ntt=1,
>tautp=1.0,
> tempi=250.0,
>temp0=300.0,
> ntpr=50, ntwx=50,
>ntxo=2,iwrap=0
> ntb=0,
>igb=8,
>
>cut=999.0,rgbmax=999.
>
> /
>
>Below 3 pictures of 3 different stages:
>http://speedy.sh/QDuvp/hAs3MT-beforeheating.tga
>http://speedy.sh/yvQnZ/hAs3MT-heat1.tga
>http://speedy.sh/wVWd6/hAs3MT-afterheating.tga
>
>Is there anything you can suggest?
>Thank you for your help.
>
>
>Regards,
>Annika
>_______________________________________________
>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 Thu Aug 07 2014 - 23:00:02 PDT
