<div>Dear Lorenzo Paulatto</div>
<div> </div>
<div>Thank you for your quick reply.</div>
<div>I have read INPUT_PP.html. </div>
<div>I don't know which plot_num can show wavefunction or how to set the input file to make it generate the files content info of wavefunction and kpoints</div>
<div>My purpose is to print out the info of wavefunction, k points and elph_matrix indexes after pw.x or ph.x  are completed.</div>
<div>Therefore, I think I have to know which matrixes collectes these info.</div>
<div>Then I can simply modify the codes to achieve my purpose.</div>
<div> </div>
<div>So far, I print out the elph_matrix by modifying elphon.f90</div>
<div> </div>
<div>elphmat(j,i)=<psi_{k+q,j}|dvscf_q*psi_{k, i}></div>
<div> ipert=           1<br> ibnd=           1<br>   ig                     evc(complex, real img)<br>    1            0.6683811960749346 -0.4860617886287616<br>    2            0.2651733872655152 -0.1928400315088636<br>
    3            0.2651733872655152 -0.1928400315088628<br>    4           -0.0005605875763412  0.0004076718519918<br>    5           -0.0005605875763366  0.0004076718519851<br>    6           -0.0005605875763413  0.0004076718519924<br>
                             .</div>
<div>                             .</div>
<div>                             .</div>
<div> </div>
<div>My question is what the ipert, psi and ig mean. Can I list the elphon matrix following (k+q,j) and (k, i) ?</div>
<div><br>Thank you for your help</div>
<div> </div>
<div>IAMS, Academic Sinica</div>
<div>Research Assistant<br></div>
<div class="gmail_quote">  2009/5/21 Lorenzo Paulatto <span dir="ltr"><<a href="mailto:paulatto@sissa.it">paulatto@sissa.it</a>></span><br>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">
<div class="im"><br>On Thu, May 21, 2009 05:08, DdrSmiLe wrote:<br>>  Does anyone have any suggestion or idea?<br><br></div>You can already do it with pp.x, please read the guide Doc/INPUT_PP.html<br>for details. Unless you take the wavefunctions from the Gamma point,<br>
wavefunctions are complex, hence the square module is saved. Nevertheless,<br>it is trivial to modify the code (PP/post_proc.f90) to make it save the<br>real and imaginary part.<br><br>At the moment (and probably forever) you can only extract wavefunctions<br>
for kpoint (what you call reciprocal vectors) that where included in the<br>self-consistent calculation. If you want an arbitrary kpoint you have to<br>pregenerate the wavefuntions doing a non-self-consistent calculation.<br>
<br>best regards<br><br>--<br>Lorenzo Paulatto<br>SISSA  &  DEMOCRITOS (Trieste)<br>phone: +39 040 3787 511<br>skype: paulatz<br>www:   <a href="http://people.sissa.it/~paulatto/" target="_blank">http://people.sissa.it/~paulatto/</a><br>
<br><br><br>----------------------------------------------------------------<br> SISSA Webmail <a href="https://webmail.sissa.it/" target="_blank">https://webmail.sissa.it/</a><br> Powered by SquirrelMail <a href="http://www.squirrelmail.org/" target="_blank">http://www.squirrelmail.org/</a><br>
<br>_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br><a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
</blockquote></div><br>