<div style="line-height:1.7;color:#000000;font-size:14px;font-family:arial">Dear Ponce:<br>I am sorry that I may not catch the right meaning of your first question. If you mean how to make extra or less charge for a model, you can use "total_charge" in &system namelist.<br>If you wanna get information of wavefunction, you'd better read Doc/INPUT_PP. For the charge density at specified k-point and band, you can calculate |psi|^2 through pp.x. Moreover, if you want to obtain the charge density for all k-points on specified bands, you can summary them by pp.x, too. P.S.: the original code of pp.x can just allow 7 filepp, but you can change the variable, nfilemax, in PP/chdens.f90, then make pp.x.<br>Good Luck.<br><div>--<br>GAO Zhe<br>CMC Lab, MSE, SNU, Seoul, S.Korea<br>        
</div><br>At 2011-07-29 17:27:27,"Samuel Ponc¨¦" <samuel.pon@gmail.com> wrote:<br> <blockquote id="isReplyContent" style="padding-left: 1ex; margin: 0px 0px 0px 0.8ex; border-left: 1px solid rgb(204, 204, 204);">Dear Quantum espresso user/developers, <br><br>I have two questions:<br><br>1) Is it possible in pwscf to change "by hand" the occupation number ( and impose a charge compensation accordingly) and if so how to do it?<br>
I would like to do so because the sum over the different mode of the partial derivative of the phonon frequency with respect of the occupation number gives the zero point motion renormalization using Brooks theorem. I would like to do this by finite difference (hence doing by hand some calculation at different occupation number).<br>
<br> 2) Is it possible to extract informations in text format form the wavefunctions? What is the right post processing tool to do so? Like having the value of the wavefunction for a defined k points and band index?<br><br>
Thanks you!<br><br>Samuel Ponc¨¦.<br>(phD student, Universit¨¦ Catholique de Louvain, Belgium).<br>
</blockquote></div><br><br><span title="neteasefooter"><span id="netease_mail_footer"></span></span>