AMBER Archive (2008)

Subject: RE: AMBER: Problems with the charge (membrane simulation)

From: Ross Walker (ross_at_rosswalker.co.uk)
Date: Thu Aug 07 2008 - 10:18:12 CDT


Hi Rebeca,

> I am trying to do a membrane+protein simulation using the Amber force
> field. I have ionized with tleap, and the charge of my system was -22.95,
> so Leap added 22 Na+ ions to neutralize:
>
> addions solvated Na+ 0
> 22 Na+ ions required to neutralize.
>
> When I use this force field in NAMD, I can see that the total charge is
> not
> zero, as is logical.
>
> Info: TOTAL CHARGE = -0.94832 e
>
> So, some warnings appear:
>
> Shoud I add 23 Na+ instead of the default Leap applies? In such way, the
> charge won´t be exactly zero, either. Do you think will be there problems,
too?

You should probably try and find out where the missing charge is. 0.052 is
too much to just be rounding errors. This suggests to me that you have an
atom missing, probably a hydrogen, or one too many atoms somewhere. Check
that each of you pieces is neutral. I.e. is the protein neutral? Then what
about each of the lipid molecules? What parameters are you using for the
lipids?

You could simply tell leap to add 23 sodium atoms instead and just ignore
the 0.052 but in this case it is probably telling you about a deeper
problem.

All the best
Ross

/\
\/
|\oss Walker

| Assistant Research Professor |
| San Diego Supercomputer Center |
| Tel: +1 858 822 0854 | EMail:- ross_at_rosswalker.co.uk |
| http://www.rosswalker.co.uk | PGP Key available on request |

Note: Electronic Mail is not secure, has no guarantee of delivery, may not
be read every day, and should not be used for urgent or sensitive issues.

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
      to majordomo_at_scripps.edu