<div dir="ltr">Hi Josip,<div><br></div><div>Mulliken charges are not real observables - they are strongly dependent on the basis set, and they do not converge as you make the basis set complete. Check a quantum chemistry textbook for a good discussion ('Computational Chemistry', Errol Lewars for instance).</div><div><br></div><div>Something more quantitative can come from something like Bader charge analysis (which needs external processing of a total_charge_density cube file).</div><div><br></div><div>Matt</div><div><br>On Monday, December 19, 2016 at 3:56:55 PM UTC, Josip Lovrić wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr"><div><span style="font-family:arial,sans-serif;font-size:small">Dear CP2K community,</span></div><div><br></div><div>Recently I was working on NO2-dimer and now I need help from you. Thing is that I am running gas phase calculation of NO2 cis-dimer. After program successfully  finish GEOMETRIC OPTIMIZATION output configuration is consistent with already published papers. Problem is in the MULLIKEN charges. Charges are too small and not consistent with published values(I am using same basis set).</div><div>!-----------------------------<wbr>------------------------------<wbr>------------------!</div><div>                     Mulliken Population Analysis</div><div><br></div><div> #  Atom  Element  Kind  Atomic population                           Net charge</div><div>       1     N        1          4.968657                              0.031343</div><div>       2     O        2          5.852787                              0.147213</div><div>       3     O        2          6.191534                             -0.191534</div><div>       4     N        1          4.928901                              0.071099</div><div>       5     O        2          6.056358                             -0.056358</div><div>       6     O        2          6.001764                             -0.001764</div><div> # Total charge                             34.000000                  0.000000</div><div><br></div><div> !----------------------------<wbr>------------------------------<wbr>-------------------!</div><div>I am using BLYP functional with dispersion corrections and TZV2P basis set with pseudo-potentials. Charges are becoming more physical when performing all electron calculation, they increase for factor ~10. </div><div>For DZVP basis set charges increase little, but still to low:</div><div> !----------------------------<wbr>------------------------------<wbr>-------------------!</div><div>                     Mulliken Population Analysis</div><div><br></div><div> #  Atom  Element  Kind  Atomic population                           Net charge</div><div>       1     N        1          4.951639                              0.048361</div><div>       2     O        2          5.822911                              0.177089</div><div>       3     O        2          6.202144                             -0.202144</div><div>       4     N        1          4.858609                              0.141391</div><div>       5     O        2          6.103330                             -0.103330</div><div>       6     O        2          6.061367                             -0.061367</div><div> # Total charge                             34.000000                 -0.000000</div><div><br></div><div> !----------------------------<wbr>------------------------------<wbr>-------------------!</div><div>  Am I doing something wrong? Do you have any explanation for this results?</div><div>I would really appreciate. </div><div><br></div><div>I am putting in attachment my input and output files. </div><div><br></div><div>Best wishes</div><div>Josip</div></div></blockquote></div></div>