> I obtain a paracetemol molecule whos net charge is -0.001 e. This seems
> large to me. I intend to simulate upwards of ~2500 paracetemol molecules,
Round off error / finite precision happens. Try figuring a way to shift
charge by hand slightly to make it net neutral. Although not absolutely
necessary from a code/Ewald point of view since the methods can handle
net-charge, from a reviewer's perspective, best to be neutral. I tend to
look for a way to shift charges to hydrogens to get rid of the net charge.
Tylenol? I would look to smear the charge on the methyl or the other
terminal OH.
> which would yield a net charge -2.5 e. I was wondering if someone could
> comment on what would cause antechamber to not return a neutral molecule.
Finite precision / round off.
> My initial guess was perhaps the initial configuration was poor, but given
No.
--tec3
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Tue May 03 2016 - 22:30:03 PDT
