<div dir="ltr"><div>Hi,</div><div><br></div>apologies for patronising. Which published values are you comparing to?<div><br></div><div>At a glance your values don't look unreasonable to me.</div><div><br></div><div>Your input looks OK, except the cutoff is rather small - to get away with that small a value you need quite heavy smoothing of the XC - do the values change if you increase the cutoff?<br><div><br></div><div>Matt<br><br>On Monday, December 19, 2016 at 4:09:18 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">Dear Matt,<div><br></div><div>I totally agree with you and I am aware of fact that atomic charge are just model. </div><div>Maybe I was not clear enough but my problem is that charges are not consistent with published values.</div><div><br></div><div>P.S. input files are named "*_all_electrons*". This is mistake, calculation is not all electron, it is a relict from previous calculation.</div><div><br></div><div>Josip<br><br>On Monday, December 19, 2016 at 5:03:50 PM UTC+1, Matt W wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><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></blockquote></div></div></blockquote></div></div></div>