<div dir="ltr">Dear <span style="font-family:arial,sans-serif;font-size:13px"><font color="#000000">Bonfa',</font></span><div><span style="font-family:arial,sans-serif;font-size:13px"><font color="#000000"><br></font></span></div>
<div><span style="font-family:arial,sans-serif;font-size:13px"><font color="#000000">Thank you for providing the code, it is really helpful. I think the codes </font></span>"  if nspin > 1 : ..." need to be changed for the noncollinear case. In this case, nspin=4, but there is no spin component part in the xml file. </div>
<div>
<p style="margin:0px"><br></p></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Apr 1, 2014 at 3:22 PM, Pietro Bonfa' <span dir="ltr"><<a href="mailto:pietro.bonfa@fis.unipr.it" target="_blank">pietro.bonfa@fis.unipr.it</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Dr. Peng Chen,<br>
<br>
I made a very simple python script that parses projwfc.x xml output (a<br>
file that is called atomic_proj.xml in the outdir directory) and writes<br>
atomic orbital contributions to the eigenfunction at each kpoint (and<br>
eigenvalues too).<br>
I uploaded it here: <a href="https://gitorious.org/qe-stuff/aopbs/" target="_blank">https://gitorious.org/qe-stuff/aopbs/</a><br>
<br>
You can use it to produce the kind of figures attached (Fe orbitals in<br>
simple cubic Fe).<br>
I used this script just for debugging purposes and never for data<br>
analysis. Please check if it gives reliable outputs before using it.<br>
<br>
To get the figure attached you need to run a band structure calculation<br>
and then use the projwfc.x postprocessing tool.<br>
My script is very badly documented but it should easy to understand.<br>
<br>
If you find errors please let me know.<br>
<br>
Best regards,<br>
Pietro Bonfa'<br>
<div><div class="h5"><br>
<br>
<br>
On 03/31/2014 09:43 PM, Peng Chen wrote:<br>
> Dear QE users,<br>
><br>
> Is it possible to project band structure onto atoms like projected DOS?<br>
> Eg. for the slab calculation, I am interested in how much the first<br>
> layer atoms contributes to the surface band. For the heretostructure 1<br>
> monolayer A+substrate B, I am wondering how much A atoms contribute to<br>
> the interface band structure. I tried k-resolved DOS, but the result<br>
> image is too fuzzy to identify the bands. Any suggestions are greatly<br>
> appreciated!<br>
><br>
> Another problem is related to Fermi surface calculation. Can I calculate<br>
> Fermi surface for a slab? I followed example08, but always get error:<br>
>      Error in routine irrek_nc (1):<br>
>      Internal problem with k points<br>
> I think the only difference is that bulk is used in example08.<br>
><br>
><br>
</div></div>> _______________________________________________<br>
> Pw_forum mailing list<br>
> <a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
> <a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
><br>
<span class="HOEnZb"><font color="#888888"><br>
<br>
--<br>
Pietro Bonfa' - PhD student<br>
Dipartimento di Fisica e Scienze della Terra "Macedonio Melloni"<br>
Viale delle Scienze 7A<br>
43124 Parma - Italy<br>
</font></span><br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br></div>