Hi all, I am using g_potential to calculate the charge distribution of my system which consists of a protein layer covalently bound to some fatty acids. The protein layer by itself would have an overall negative charge but they system as a whole is neutral. My question is this: When I specify different values for the number of slices, and I sum up the charges for the protein only, I always get different results with different number of slices. But shouldn't the overall charge of the protein layer be the same regardless of the number of slices or the size of each bin? The difference is also quite significant so I don't think it is simply a rounding error. For example with a 0.1 nm bin width, I get an overall charge of -3, but with a bin width of 0.05nm, the overall charge is -7. Could someone possibly explain to me what is going on?
Thanks! ________________________________ This email is confidential and may be privileged. If you are not the intended recipient, please delete it and notify us immediately. Please do not copy or use it for any purpose, or disclose its contents to any other person. Thank you.
_______________________________________________ gmx-users mailing list [email protected] http://www.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to [email protected]. Can't post? Read http://www.gromacs.org/mailing_lists/users.php
