On Mon, 2014-11-10 at 13:51 +0000, CHAMI F. wrote:
> Hi Jason,
>
> I run sqm with the option "printdipole=1"
> the dipole moment that SQM output is for the AM1 before the bond charge correction ..AM1-BCC is
> final partial charges in prepi file
>
> in sqm.out for example the AM1 charge on nitrogen in cyano group of biphenylcyano molecule is
> so different from AM1-BCC
>
>
> sqm.out 14 N -0.038
> ANTECHAMBER_AM1BCC_PRE.AC
> ATOM 14 N 2 -8.532 1.150 0.912 -0.038000 n1
> ANTECHAMBER_AM1BCC.AC
> ATOM 14 N 2 -8.532 1.150 0.912 -0.363800 25
>
> this mean i have to work out the dipole moment using the partial charges in prepi file
There are two dipole moments here -- one from the QM calculation based
directly on the density matrix and another from the partial atomic
charges you derive _from_ the calculation. In an ideal world, these
would be equal. However, point charges in general are insufficient to
reproduce exact QM dipole moments (except, of course, when the dipole
moment is exactly 0 by symmetry).
So if you want the QM dipole moment, use the value you get from sqm. If
you want the MM dipole moment, as I mentioned before, you have to
compute it yourself from the charges and coordinates in the prepin file.
HTH,
Jason
--
Jason M. Swails
BioMaPS,
Rutgers University
Postdoctoral Researcher
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Nov 10 2014 - 06:30:02 PST
