Hi Huazhong,<div><br></div><div>Only problem I can spot is that the smearing is a little bit too large. Room temperature is only 23meV, degauss = 0.002Ry would be more reasonable. I guess the situation is that the single electron should be a d state of Ti, but PBE functional makes it a delocalized state and the system becomes metallic, then big smearing completely eliminates magnetization... Just my guess...</div>
<div><br><div class="gmail_quote">On Thu, May 17, 2012 at 9:58 PM, Huazhong Liu <span dir="ltr"><<a href="mailto:huazhorg@163.com" target="_blank">huazhorg@163.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="line-height:1.7;font-size:14px;font-family:arial"><div>Dear all,</div><div>     I have got a problem of spin-polarization calculation. I put one H atom onto the stoichiometric TiO2 surface. I think that is must have magnetic properties after I optimized the adsorption system. But I didn't find the magnetism in the system after the calculation. Here are the parameters which set for spin-polarization calculation:</div>
<div>...................</div><div> /</div><div> &system</div><div>    ibrav=..........</div><div>    nat=49  , ntyp=3,</div><div>    nspin=2,</div><div>    starting_magnetization(1)= 0.5</div><div>    starting_magnetization(2)= 0.5</div>
<div>    starting_magnetization(3)= 0.5    </div><div>    ecutwfc=30.0,</div><div>    ecutrho=300.0,</div><div>    occupations='smearing', smearing='gauss', degauss=0.02,</div><div> /</div><div>ATOMIC_SPECIES</div>
<div> Ti   47.867       Ti.pbe-sp-van_ak.UPF</div><div> O    15.9994      O.pbe-van_ak.UPF</div><div> H    1.00794      H.pbe-van_ak.UPF</div><div>...........................</div><div><br></div><div>I got the result as below:</div>
<div>      The total energy is the sum of the following terms:</div><div><br></div><div>     one-electron contribution =  -12795.07240972 Ry</div><div>     hartree contribution      =    6472.21004247 Ry</div><div>     xc contribution           =    -456.66736582 Ry</div>
<div>     ewald contribution        =    3879.97911645 Ry</div><div>     smearing contrib. (-TS)   =      -0.01329486 Ry</div><div><br></div><div>     total magnetization       =     0.00 Bohr mag/cell</div><div>     absolute magnetization    =     0.00 Bohr mag/cell</div>
<div>    </div><div> convergence has been achieved in  10 iterations</div><div>..................................</div><div><br></div><div>Will you please find any problem in my parameters? Why I can not get the magnetism of the system?  Thank you!</div>
<div><br></div><div>regards,</div><div>Huazhong Liu</div></div><br><br><span title="neteasefooter"><span></span></span><br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br>Jia Chen<br><br>Dept of Chemistry<br>Princeton University<br>Princeton, NJ 08544<br>
</div>