<div dir="ltr">Hi <span style="font-size:12.8px">Rolly!</span><div><span style="font-size:12.8px">I calculated the same object with SIESTA, <a href="https://www.researchgate.net/publication/285619763_Doping_graphene_with_a_monovacancy_bonding_and_magnetism">https://www.researchgate.net/publication/285619763_Doping_graphene_with_a_monovacancy_bonding_and_magnetism</a> and other 49 elements of periodic table</span></div><div><span style="font-size:12.8px">interesting to compare (note spins are not correct).</span></div></div><div class="gmail_extra"><br><div class="gmail_quote">2016-07-19 19:00 GMT+04:00 Rolly Ng <span dir="ltr"><<a href="mailto:rollyng@gmail.com" target="_blank">rollyng@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Prof. Paolo,<br>
<br>
I can confirm that pw.x with CPU works okay for nspin=2, so I would like<br>
to report the bug for the pw-gpu.x (GPU version).<br>
<br>
Thanks for your attention.<br>
<br>
Regards,<br>
<br>
Rolly<br>
<div class="HOEnZb"><div class="h5"><br>
<br>
On 07/18/2016 02:16 PM, Paolo Giannozzi wrote:<br>
> You are using the GPU-enabled version of QE. Try the normal version<br>
> first. You do not need more than a few bands over the Fermi energy.<br>
> Setting mixing_mode='TF-local- improves the convergency.<br>
><br>
> Paolo<br>
><br>
> On Sat, Jul 16, 2016 at 6:34 PM, Rolly Ng <<a href="mailto:rollyng@gmail.com">rollyng@gmail.com</a>> wrote:<br>
>> Dear QE experts,<br>
>><br>
>> I would like to compute the DOS of Nickel doped graphene supercell with QE<br>
>> v5.3.0, and our group has already done that with DMol3.<br>
>><br>
>> Our previous publication can be found here, Ding, N., Lu, X., Wu, C.-M.L.,<br>
>> 2012. Comput. Mater. Sci. 51, 141–145.<br>
>><br>
>> Now, I like to refine the computation with QE, so at first I used "relax" to<br>
>> get the optimized structure of a graphene supercell with of 72 Carbon atoms<br>
>> and a Nickel atom was added at the center to replace one of the carbons.<br>
>><br>
>> Then I run scf with nspin=2 and starting magnetization (1) = 0.5 where (1)<br>
>> refers to the Nickel atom, the scf converged. However, at the very end pw.x<br>
>> stopped with the following error.<br>
>><br>
>>       WARNING: integrated charge=   294.76903202, expected= 302.00000000<br>
>><br>
>><br>
>> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br>
>>       Error in routine electrons (1):<br>
>>       charge is wrong<br>
>><br>
>> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br>
>><br>
>>       stopping ...<br>
>><br>
>> I googled and found that it may be related to smearing and nbnd, I have<br>
>> already set smearing to 0.002 Ry (312K) and I have increased nbnd to 400 but<br>
>> it is still running.<br>
>><br>
>> Can anyone help me to resolve the issue? Please find my input and output<br>
>> attached.<br>
>><br>
>> Thank you very much,<br>
>> Rolly<br>
>><br>
>> --<br>
>> PhD. Research Fellow,<br>
>> Dept. of Physics & Materials Science,<br>
>> City University of Hong Kong<br>
>> Tel: +852 3442 4000<br>
>> Fax: +852 3442 0538<br>
>><br>
>><br>
>> _______________________________________________<br>
>> Pw_forum mailing list<br>
>> <a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
>> <a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
><br>
><br>
<br>
--<br>
PhD. Research Fellow,<br>
Dept. of Physics & Materials Science,<br>
City University of Hong Kong<br>
Tel: +852 3442 4000<br>
Fax: +852 3442 0538<br>
<br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a></div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Best wishes,</div><div>Maxim Arsent'ev, Ph.D. (Chemistry)</div><div>Laboratory of research of nanostructures</div><div>Institute of Silicate Chemistry of RAS</div></div></div></div></div>
</div>