Thanks, Jia and Gao. Your suggestions are helpful.<br><br><div class="gmail_quote">On Fri, Sep 7, 2012 at 2:49 PM, GAO Zhe <span dir="ltr"><<a href="mailto:flux_ray12@163.com" target="_blank">flux_ray12@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">Dear Peng:<br>To show the difference between DOSes, the best way is put the Fermi level at same place. Normally, (E-Ef) is preferred by lots of researches as the variable of x-axes. Manually shifting DOS for matching is not a good idea, since HOMO should be the standard and adjusted to the same position (although Fermi level given by PWscf is LUMO for semiconductor or insulator).<br>
BTW, in PWscf, nscf calculation also provides Fermi energy, which usually is slightly different with that from scf calculation and displayed by setting verbosity='high' in &control namelist.<br><br><div>--<br>
GAO Zhe<br>CMC Lab, Materials Science & Engineering Department,<br>Seoul National University, South Korea<br>        
 </div><div><div class="h5"><div></div><br>At 2012-09-08 01:16:39,"Peng Chen" <<a href="mailto:pchen@ion.chem.utk.edu" target="_blank">pchen@ion.chem.utk.edu</a>> wrote:<br> <blockquote style="PADDING-LEFT:1ex;MARGIN:0px 0px 0px 0.8ex;BORDER-LEFT:#ccc 1px solid">
Dear All,<div><br></div><div>I calculated the total density of states (in the attachment) for both afm and fm states of an insulator. I shifted the curves so that the band gap starts at 0 eV. In the top figure, the Fermi energy obtained from the scf calculation falls between the gap. So if I plot Fermi surface, can I get reliable result? In the bottom figure, I shifted the spin down DOS of fm state so that it can match with the spin down DOS of afm state. I am not sure if that is the right way to show the DOS difference of these two states.</div>



<div><br></div><div>Another question is about the accuracy of calculation. In the scf calculation, I made the total energy converged within 0.001 Ry with related to the ecutwfc, degauss, kpoints...   Then I did band structure calculation. Can I say the error of calculated band energy is within 0.001 Ry? </div>



<div><br clear="all">
<div><br></div><br>  Best Regards.<br>        Peng  <br>
</div>
</blockquote></div></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>  Best Regards.<br>        Peng  <br>