<div dir="ltr">Dear Yun,<div><br></div><div>You can try to compute the charges using self-consistent hirshfeld partitioning (use SELF_CONSISTENT T inside the section &HIRSHFELD..&END) and using SHAPE_FUNCTION DENSITY. It has been highlighted before that there is in general a dependence of the charges on the basis set on MoS2 (see https://groups.google.com/forum/#!topic/cp2k/gZgiyqZURs4)</div><div><br></div><div>For the DZVP-MOLOPT-SR basis, setting self consistent partitioning will give the right sign of the charges. We have also noticed that using this scheme there is not a significant dependence of the charges wrt to the size of the basis set. However, it is worth checking the dependence on the basis set size of the geometry and of other properties that may be relevant for your specific case.</div><div><br></div><div>For that, you may want to check the file BASIS_MOLOPT_UCL, which is an extended molopt basis set library.</div><div><br></div><div>https://sourceforge.net/p/cp2k/code/HEAD/tree/trunk/cp2k/data/BASIS_MOLOPT_UCL<br></div><div><br></div><div>Best</div><div>Gabriele</div><div><br><br>On Monday, November 20, 2017 at 2:26:15 PM UTC+1, any...@gmail.com wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr"><div>Dear all,</div><div><br></div><div>I run the cell_opt  for 4x4x1 MoS2 system and I want to analyze the Hirshfeld charges, the result shows the Net charge of Mo and S are -1.243 and 0.622, respectively.  But I thought the charge of Mo should be positive while S should be negative.  </div><div><br></div><div>Can anybody help me to see if there are any problems or I misunderstood it?  Thank you very much in advance. </div><div><br></div><div>Kind regards,</div><div>Yun</div></div></blockquote></div></div>