CCL: Charge



 Sent to CCL by: Peeter Burk [peeter.burk-.-ut.ee]
 Hi,
 
Can we even discuss which charge partitioning scheme is better (aside of convergency with increasing basis set size) as AFAIK atomic charges are not a measurable property of atoms in molecules or materials?
 
Chemists (including myself) like to use them to explain properties like reactivity, but as long as you use them within series most schemes give the similar answers (if you compare trends and changes, not absolute values)...
 Best regards
 Peeter
 On 28/06/2020 19:56, Thomas Manz thomasamanz(a)gmail.com wrote:
 
 Hi Partha,
 Thank you for your question.
 
Atom-in-material properties should not explicitly depend on the basis set used to perform the quantum chemistry calculation. The reason for this is that nature corresponds to the limit of an arbitrarily large basis set (aka 'complete basis set limit'). For reasons of computational efficiency, in practical calculations a finite basis set is usually used to get an answer quickly. So, we want to use population analysis methods that are not too sensitive to the choice of basis set, so that reasonably sized basis sets will produce answers that are not much different than what would be computed if an arbitrarily large basis set (aka 'complete basis set limit') were used.
 
The Mulliken method explicitly depends on the basis set choice and has no complete basis set limit. In other words, the Mulliken populations are not defined for arbitrarily large basis set. This is very bad, because it means that as the accuracy/precision of the quantum chemistry calculation is improved by using larger and larger basis set, the Mulliken charges get worse and worse with no mathematical limits or meaning. So, you should not use Mulliken populations. It also means that if you run a similar calculation twice, but using different basis sets for each calculation, the Mulliken charges that you get from the two similar calculations may not be similar. For example, Mulliken population analysis may predict the cations have turned into anions, and vice versa, even though it is exactly the same material. In other words, Mulliken population analysis often gets cations and anions confused; it does not know the distinction between them.
 
Between the two methods you mentioned, natural population analysis (NPA) is the better choice. An even better choice would be to use a method like DDEC6 for which the atom-in-material properties are a functional of the electron and spin density distributions with no explicit basis set dependence, and for which the atom-in-material properties are chemically meaningful across a broad range of material types.
 Sincerely,
 Tom
 
On Sun, Jun 28, 2020 at 8:43 AM Partha Sengupta anapspsmo%x%gmail.com <http://gmail.com>; <owner-chemistry*_*ccl.net <mailto:owner-chemistry*_*ccl.net>> wrote:
     Sir, For a metal complexes involving Cu[N,O donor], [Chlorine and
     fluorine atoms are in the benzene ring]. I found that cu has
     0.96091 charges( Natural Population analysis) while Mulliken
     atomic charge is 0.457333. On the same process N has -0.53717 and
     -0.29941 respectively. What is the reason behind this? What will
     be the better choice to represent?
     Partha Sengupta
 
-- */Dr. Partha Sarathi Sengupta
     Associate Professor
     Vivekananda Mahavidyalaya, Burdwan/*