QUESTION ABOUT COULOMBIC POTENTIAL



SUBJECT: QUESTION ABOUT COULOMBIC POTENTIAL
 Dear all,
 I am trying to write a program to calculate the electrostatic contribution of
 the lattice energy for molecular crystals. I have run plane-wave DFT
 calculations (using the CASTEP package) to calculate the total lattice energy
 and using Mulliken charges I hoped to estimate the electrostatic contribution
 by coding a simple Coulombic potential. But, I get really weird results!! I am
 overestimating my lattice energy by many orders of magnitude.
 I know that Mulliken charges are not the ideal point charge to use for this but
 I had not expected my numbers to be that wrong. Could anybody help or if not,
 could anybody recommend and existing program to use instead (preferably free)?
 Many thanks,
 Patricia Lozano-Casal (p.lozano-casal|at|ed.ac.uk)
 --
 Room 85, School of Chemistry,
 Joseph Black Building,
 Kings Buildings,
 West Mains Road,
 Edinburgh, EH9 3JJ
 United Kingdom
 Tel: +441316504804