<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: http://dx.doi.org/10.1088/2053-1583/3/2/025026), 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>