<div>Dear users,</div>
<div> </div>
<div>I wanted to compute the TOTAL DOS for a non-spin-orbit spin-polarized calculation (FM and AFM) </div>
<div>so I did the following after scf convergence:</div>
<div> </div>
<div>x lapw2 -up -qtl</div>
<div>x lapw2 -dn -qtl </div>
<div> </div>
<div>Then I used the following <a href="http://case.int">case.int</a> file </div>
<div> </div>
<div>++++++++++++++++++++++++++++++++++++++++++++++++++++++</div>
<div>+   Title                                                                                               +<br>+  -0.50 0.002 1.500 0.02     EMIN, DE, EMAX, Gauss-broadening(&gt;de)   +<br>+    1                      NUMBER OF DOS-CASES specified below             +<br>
+    0   1   total         atom, case=column in qtl-header, label                   +</div>
<div>++++++++++++++++++++++++++++++++++++++++++++++++++++++</div>
<div> </div>
<div>to compute the DOS using tetra: </div>
<div> </div>
<div>x tetra -up</div>
<div>x tetra -dn</div>
<div> </div>
<div> </div>
<div>I expected case.dos1evup and case.dos1evdn plots to look exactly the same </div>
<div>but they do not. The same is true for the spin-up and spin-down band structure </div>
<div>(with jatom=0, jtype=1 in case.insp). I find this puzzling. Can someone kindly </div>
<div>explain?</div>
<div> </div>
<div>Thank you.</div>
<div> </div>
<div>F. Garcia</div>
<div> </div>
<div> </div>