<div dir="ltr">Hi,<div><br></div><div>just out of curiosity could you try using the new TZV basis sets (for Mo they are in BASIS_MOLOPT_UCL) and see how that looks?<br><br>Mulliken charges don't really mean anything much, so if band gap and geometry are correct then it just indicates the basis set on Mo is better than S.<br><br>Matt<br><br>On Thursday, June 29, 2017 at 5:52:34 PM UTC+1, Chanwoo Noh wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Dear all,<div><br></div><div>Recently, I started to study on MoS2 system and performed DFT calculation of MoS2 monolayer using CP2K.</div><div>At first, geometry optimization calculation was performed. </div><div>The calculated MoS2 monolayer structure seems to be right, but the calculated Mulliken charge of atoms are wrong.</div><div>As I know, Mo charge should be positive, and S charge should have  negative value.</div><div>However, in calculation result, their charges are opposite.</div><div>Since the optimized structure of the MoS2 layer is right, I do not think that my calculation is completely wrong, but I wonder what causes the atomic charge to be opposite.</div><div>In any case, how can I check my calculations?</div><div>Please give me some advice.</div><div>I would appreciate if you give me any comment.</div><div><br></div><div>I used PBE functional, DZVP-MOLOPT-SR basis, and cutoff as 500Ry. <br></div><div>I attach my input and output files.</div><div> </div></div></blockquote></div></div>