[Wien] Total density of states versus projected density of states

Jian-Xin Zhu jxzhu at lanl.gov
Thu Jun 18 05:30:07 CEST 2009


Dear Prof. Blaha,

I got the point.
Thank you very much.

Jian-Xin


On Jun 17, 2009, at 2:48 PM, Peter Blaha wrote:

The UG says:

Please note: The total DOS is equal to the sum over the atoms of the  
total-atomic DOS (inside spheres) and the interstitial-DOS. (Thus in  
the total-atomic DOS the “multiplicity” of an atom is considered). On  
the other hand, in the partial (lm-like) DOS the multiplicity is not  
considered and one obtains the total-atomic DOS as a sum over all  
partial DOS times the multiplicity.

So the total Fe DOS is already multiplied by 2.
On the other hand the Fe-s,p,d,eg,... DOS would NOT be multiplied by  
MULT.

Just use "common sense": The total DOS is 6.4; the sum of the DOS  
inside spheres is 5.6, so obviously 0.8 is the interstitial.

You can also check case.outputt, where next to the DOS also the  
integrated DOS is listed. Compare the integrals at EF with the partial
charges in case.scf (:QTL and :CHA0xx ).

> I did read the TETRA section of the UG.  I still need to understand  
> better.
> Here  I send you several files of the density of states and
> and the input parameter files.   The structure file is also attached.
> The calculation is done with the spin-polarized and spin-orbit  
> coupling.
> If I look at the DOS at the Fermi energy, e.g., at -0.00476eV,
> the second column in the file of UFe5As3-FM.dos1evup,
> the total DOS value is 6.41498709,
> and the following columns for the total DOS of U, Fe1  Fe2  Fe3 Fe4  
> Fe5  As1 As2 As3, are respectively
> 1.28828049    0.60202104    0.69573355    0.86723715    0.88412511  
> 0.91823441  0.14181454    0.11641405    0.10355555
> (The last three values are located from UFe5As3-FM.dos2evup  
> according the setting in UFe5As3-FM.int.)
> If add them up, the sum is 5.61741584.
> But I know each of these atoms has MULT=2.
> So if I time the sum by 2, I will get 5.61741584x2 = 11.23483168  
> which is way too bigger than 6.41498709.
> This disagreement makes me wonder:
> Are you implying that the total DOS for a specified atom (e.g., U  
> etc) is also for the whole unit cell, which has included the  
> multiplicity MULT,
> while  lm decomposed DOS, e.g., 5f, 6d,7s, is per single atom?
> In the TETRA section of UG, it say
> jatom=0 means total DOS
> then
> jcol=1 also means total, which is for a specific type atom when  
> jatom is not equal to 0.
> So I can ask the question in another way ---
> When jatom is set to nonzero, does the "total" for the setting  
> jcol=1 mean the total DOS of this type of atom (including already  
> MULT) in the unit cell?
> I very much appreciate for your answer and look forward to hearing  
> from you.
> Best regards,
> Jianxin
> On Jun 16, 2009, at 3:17 PM, Peter Blaha wrote:
> The total DOS is per unit cell.
> The partial, lm decomposed DOS is per single atom (and you have to  
> multiply it by MULT to get it per unit cell). Please read the TETRA  
> section of the UG.
> Jian-Xin Zhu schrieb:
>> Dear Prof. Blaha and wien2k users,
>> Recently I did an electronic structure calculation for UFe5As3 in a  
>> ferromagnetic state.
>> Although the chemical formula is UFe5As3 involving 9 atoms, the  
>> unit cell for the calculation contains
>> 2 U atoms,  10 Fe atoms, and 6 As atoms, i.e., 18 atoms in total.
>> My question:
>> Is the obtained value of the total density of states for the sum of  
>> the contributions from 9 atoms (i.e., per chemical
>> formula) or 18 atoms (per unit cell)?
>> Is the projected density of states for each individual atom or each  
>> individual type of atom?
>> Here each individual atom is really meaning one single atom while  
>> each individual type atom is meaning the
>> nonequivalent type, which can have multiple equivalent atoms. For  
>> example, the atomic type 1, which is U, can
>> have two atoms.
>> If I understand that the projected density of states is for one  
>> single atom, I would understand the total density
>> of states is for one chemical formula.
>> Thanks for the help and sharing.
>> Jian-Xin Zhu
>> _______________________________________________
>> Wien mailing list
>> Wien at zeus.theochem.tuwien.ac.at
>> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
_______________________________________________
Wien mailing list
Wien at zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien








More information about the Wien mailing list