<div dir="ltr"><div>Dear all,</div><div><br></div><div>I did some calculations for induced magnetism in graphene due to single vacancy (sometimes reffered as V1(5-9))</div><div>using USPP (C.pbe-van_ak.UPF; ecutwfc=30, ecutrho=160) and NCPP (C.pbe-mt_fhi.UPF, ecutwfc=35, ecutrho=140).</div>
<div>Calculations of relaxed structure give magnetic moment in the supercell of about 0.65 \mu_B in both the PP cases.</div><div>However, I did calculations of the same structure using full-potential LAPW code FLEUR where the magnetic moment </div>
<div>(spin_up - spin_down) drops to value of 0.12256. There are also studies where similar structures have been studied </div><div>with SIESTA code:</div><div>depending of supercell size from 1.1 to 1.5 \mu_B [Yazyev et al, Phys Rev B 75, 125408]</div>
<div>Normconserving PPs give about 1.2 \mu_B [<a href="http://arxiv.org/abs/1006.0589">http://arxiv.org/abs/1006.0589</a>]</div><div>GIPAW gives 1.53 \mu_B [<a href="http://www.researchgate.net/publication/235566299_Non-collinear_magnetic_moments_in_graphene_with_vacancies">http://www.researchgate.net/publication/235566299_Non-collinear_magnetic_moments_in_graphene_with_vacancies</a>]</div>
<div><br></div><div>It seems that USPP and NCPP could not be the reason for the discrepancy, or?</div><div>Could you comment please on the origin of the discrepancies between QE and LAPW?</div><div><br></div><div>Best regards,</div>
<div><br></div><div>Martin Gmitra</div><div>Uni Regensburg, Germany</div></div>