[Pw_forum] Projected density of states

Pingo Mutombo mutombo at fzu.cz
Thu Mar 31 20:06:44 CEST 2005


Hello Gerardo,
Many thanks for the explanation. What I need is to project the 
wavefunctions onto atomic states of particular atoms. In my case, I want 
that "fraction" of the wavefunctions which are localised on adatoms and 
rest atoms. For some reasons , the program" projwfc" stops without 
giving any error message. I thought it was related to the size of the 
unit cell.
Best regards,
                                             Pingo


Gerardo Ballabio wrote:

> On 03/31/05 19:14:40, Pingo Mutombo wrote:
>
>> How can I compute the projected density of states for selected  atoms 
>> of a given supercell? I can calculate the total DOS for a   
>> Si(111)-7x7 surface. However  I am unable to calculate the   
>> projected density of states  probably because of the big size of  the 
>> unit cell(249 atoms). At least I would like to calculate  the  
>> projected density  of states for selected Si atoms (adatom and rest  
>> atoms).
>
>
> Hi Pingo,
> "density of states for selected atoms" isn't a well defined concept  
> in density-functional theory. Indeed, wavefunctions computed by PWscf  
> (or any other DFT code) are "collective" and can't be associated with  
> any particular atom.
>
> You may project the wavefunctions onto atomic states:
>
>   |psi> = (sum over atomic states) <atomic state|psi> |atomic state>
>
> where "atomic state" is s, px, py, pz... for a particular atom; this  
> gives you the "fraction" of wavefunction localized on that atom. Keep  
> in mind however that this is an approximation, because atomic states  
> of different atoms aren't orthogonal in general, thus they don't form  
> an orthonormal set. As a result, the sum of such "fractions" over all  
> atoms may not be exactly 1, and the sum of all "atomic" DOS may not  
> equal the total DOS.
>
> Gerardo
>
> _______________________________________________
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://www.democritos.it/mailman/listinfo/pw_forum
>


-- 
Quantum Potes Tantum Aude
Odvaz se,kolik jen muzes.
As much as you are able,that much dare to do.
Aie de l'audace,autant que tu le peux.

                    St.Thomas Aquinas
-------------------------------------------------------------------------------------------
Dr.Pingo Mutombo
Institute of Physics of the Academy of Sciences of the Czech republic
Cukrovarnicka 10
162 53 Praha 6
Czech republic
Tel:+420-220-318-528		    
		    
		    




More information about the Pw_forum mailing list