<div dir="ltr">Hi Qiang,<div><br></div><div>I think that dissociation into partially charged ions is the 'correct' incorrect behaviour using standard functionals - see</div><div><br></div><div>'@article{ruzsinszky2006spurious,<div>  title={Spurious fractional charge on dissociated atoms: Pervasive and resilient self-interaction error of common density functionals},</div><div>  author={Ruzsinszky, Adrienn and Perdew, John P and Csonka, G{\'a}bor I and Vydrov, Oleg A and Scuseria, Gustavo E},</div><div>  journal={The Journal of chemical physics},</div><div>  volume={125},</div><div>  number={19},</div><div>  pages={194112},</div><div>  year={2006},</div><div>  publisher={AIP Publishing}</div><div>}'</div><div><br></div><div>for instance. </div><div><br></div><div>Also, it is quite hard to converge wavefunctions for the stretched molecule - I'm not sure your set up really converges. I suspect that Gaussian converges to a different solution - though you don't say what setup you are using there.</div><div><br></div><div>So, actually what you are trying to do is quite difficult, rather than an easy task. A better approach is to model the NaCl molecule near equilibrium then calculate the Na and Cl atoms separately for most thermodynamic applications.</div><div><br></div><div>Matt</div><br>On Sunday, November 13, 2016 at 7:47:13 AM UTC, Qiang Li wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Dear Vladimir,<div><br></div><div>The UKS does not help to separate the charges, as I expected. I did a calculation with same system by Gausssian 09 for a test, and the Mulliken charge analysis shows that the Na(+1)Cl(-1) are well separated. I do expect similar result from CP2k, since it also describes the valence electrons with gaussian basis sets.  I am uploading the input file for my cp2k run, which you can have a quick look whether I did something wrong or not.</div><div><br></div><div>Thank you for being replying.</div><div><br></div><div>Best regards,</div><div>Qiang<br><br>在 2016年11月10日星期四 UTC+8下午7:51:26，Vladimir Rybkin写道：<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Li,<br><br>this means you have a "covalent" dissociation into atoms, not ions. Try to do a spin-unrestricted calculation, i.e. UKS (alternative name - LSD) calculation in the DFT section.<br><br>Yours,<br><br>Vladimir<br><br>четверг, 10 ноября 2016 г., 12:01:16 UTC+1 пользователь Qiang Li написал:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><br>Hi, Vladimir,<div>Thanks for your reply, that was quite helpful.</div><div>I did some tests before I reply back. Simply, I calculated Na+Cl- system (distance is set to 10 angstrom, using a 20x10x10 box), and did the Mulliken charge analysis. I found that using  GPW method, i.e. using the GTH potential and DZVP basis set (same method in this paper: <a href="http://dx.doi.org/10.1088/2053-1583/3/2/025026" rel="nofollow" target="_blank" onmousedown="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fdx.doi.org%2F10.1088%2F2053-1583%2F3%2F2%2F025026\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHH8gN7c8Wb30YuIidD86CBzxggcg';return true;" onclick="this.href='http://www.google.com/url?q\x3dhttp%3A%2F%2Fdx.doi.org%2F10.1088%2F2053-1583%2F3%2F2%2F025026\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHH8gN7c8Wb30YuIidD86CBzxggcg';return true;">http://dx.doi.org/10.1088/<wbr>2053-1583/3/2/025026</a>), the charges was Na(0.5)Cl(-0.5), which is very strange to me. Alternatively, the charges are well separated when I apply GAPW method (Na+Cl-).</div><div>What's the physics behind this?</div><div>Best regards,</div><div>Li</div><div><br></div><div><br></div>在 2016年11月8日星期二 UTC+8下午11:03:55，Vladimir Rybkin写道：<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear Li,<br><br>there are Mulliken and Hirshfeld population analyses in CP2K, which provide partial charges on atoms. <br>To keep charges separated you need a large simulation cell. Alternatively, you may consider a charged cell (no counter ions). This is simpler to do, although, the large cell is even more a must.<br><br>Yours,<br><br>Vladimir <br><br>вторник, 8 ноября 2016 г., 15:31:07 UTC+1 пользователь Qiang Li написал:<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>I found several publications dealing with metal ions in the AIMD calculations, Zn(2+) for example. The way I see it is that the authors use <span>Cl</span><span>- </span><span>counter ions to neutralize the system. However, how one can make sure that the charge of the two ions are well separated, i.e. Zn (2+) and Cl (-1) ? is there any charge analysis available  in CP2K 4.1 ? Any suggestions are appreciated.</span></div><div>All the bests,</div><div>Li</div> </div></blockquote></div></blockquote></div></blockquote></div></blockquote></div></div></blockquote></div></div>