<div dir="ltr"><div><div><div><div><div><div><div><div><div><div>Hello,<br><br></div>I have been working on this myself, so I hope someone with more experience than myself will fill in the gaps and/or correct me. What I have been doing is comparing the DOS with the PDOS and this is how...<br><br></div>Step 1 (after SCF and BANDS calculations) : calculate the DOS with dos.x and an input file using namelist &dos and setting up the relevant parameters (Emin, etc)<br><br></div>Step 2: calculate the projections with projwfc.x and an input file with namelist &projwfc ( I run twice with kresolveddos both true and false )<br><br></div>I plot the DOS and the output files from the projwfc.x calculation which have names like<br><br>PREFIX.pdos_atm#(ATOM)_wfc#(L)<br><br></div>for each L/M combination. These files give the DOS in a variety of ways (still working on the details for each L/M myself, so I cannot say more yet).<br><br></div>The output from the calculation itself ( I redirect STDOUT to a file) provides the following listing of states<br><br>     state #   1: atom   1 (P  ), wfc  1 (l=0 m= 1)<br>     state #   2: atom   1 (P  ), wfc  2 (l=1 m= 1)<br>     state #   3: atom   1 (P  ), wfc  2 (l=1 m= 2)<br>     state #   4: atom   1 (P  ), wfc  2 (l=1 m= 3)<br>    . . . <br><br></div>where this is a snippet from phosphorene. Then following that listing, each state at each k is explicitly listed as in the following snippet<br><br> k =   0.0000000000  0.0000000000  0.0000000000<br>==== e(   1) =   -17.93363 eV ====<br>     psi = 0.223*[#   5]+0.223*[#   9]+0.223*[#  13]+0.222*[#   1]+0.021*[#  10]+<br>          +0.021*[#   2]+0.020*[#   6]+0.020*[#  14]+0.006*[#   3]+0.006*[#   7]+<br>          +0.006*[#  11]+0.006*[#  15]+<br>    |psi|^2 = 0.997<br>==== e(   2) =   -15.46040 eV ====<br>     psi = 0.197*[#   1]+0.197*[#   9]+0.195*[#  13]+0.195*[#   5]+0.037*[#   7]+<br>          +0.037*[#  15]+0.036*[#   3]+0.036*[#  11]+0.016*[#   6]+0.016*[#  14]+<br>          +0.015*[#   2]+0.015*[#  10]+<br>    |psi|^2 = 0.992<br><br></div>which shows you the wavefunction in the basis listed in the snippet above. Where you see #1 or #15, etc, you can map back to the L/M state listed above.<br><br></div>Notice that M does not range from -L to L. Instead you have to read the documentation to see what is actually being represented. This is something I am parsing myself at the moment, so I cannot help much there either.<br><br><br></div>I hope that can get you started.<br><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div>Jeff Mullen, PhD<br><div>North Carolina State University</div></div></div></div></div>
<br><div class="gmail_quote">On Mon, Sep 14, 2015 at 5:58 AM, Sharma SRK Chaitanya Yamijala <span dir="ltr"><<a href="mailto:sharmajncasr@gmail.com" target="_blank">sharmajncasr@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear all,<div><br></div><div>Can some one guide me on how to achieve orbital projected band structure. I have seen projwfc.x and I only k-resolved dos.</div><div><br></div><div>Thanking you for your help,</div><div>Sincerely,</div><div>Sharma.<br clear="all"><div><div><div dir="ltr"><span style="color:rgb(204,51,204)"><font color="#000000"><br><br><br><br><br><br><br></font>********************************************************</span><br><font size="2"><b style="color:rgb(51,51,255)">Chaitanya Sharma,</b></font><br><b><span style="font-family:times new roman,serif">Prof. Pati'</span></b>s group,<br><span style="color:rgb(56,118,29)">Chemistry and Physics Materials unit</span>,<br><span style="color:rgb(255,0,0)">JNCASR</span>, BANGLORE,<br>Lab:: (080-2208) <span style="color:rgb(0,0,255)">2581, 2809</span><br><a href="https://sites.google.com/site/sharmasrkcyamijala/" target="_blank">https://sites.google.com/site/sharmasrkcyamijala/</a><br><span style="color:rgb(204,51,204)">*********************************************************</span><br></div></div></div>
</div></div>
<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" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br></div></div>