Dear all<br><br>I have some questions about the post processing the results from pw.x<br><br>first, i want to know the nature of bond (ionic or covalent). i think i should do overlap population analysis to identify the bond type.<br>
As Segall, M. D. reported in <em>Phys. Rev. B</em>, <strong>54</strong>, 16317, there is a criterion for bonding between atoms.  A high value of the bond <br>population indicates a covalent bond, while a low 
value indicates an ionic interaction. I know projwfc.x do Lowdin population analysis<br>instead of Mulliken. I only could find Lowdin charges directly in output file. and i wonder if projwfc.x could give us the information <br>
about overlap population. how to extract that information?<br><br>second, I want to plot the bonding charge density of water adsorbed on TiO2 surfaces as the Exercise 4  of postproc.pdf in <br>'hands on tutorial of QE'. I think i should  relax the water/TiO2 system and then run 3 scf calculations based on the relaxed systems, <br>
and extract the charge density for the water/TiO2  system, the water alone and the TiO2 surface alone. But when water dissociate <br>on the surface, the scf calculation of dissociated water is not converged. how could we plot the bonding charge density for this case.<br>
<br>At last i still have one more question, if I want to calculation the noble metal on TiO2 surfaces, such as Pt, Au or Pd etc., how could I<br>deal with the 'nbnd' options? should we add 20% more states to the system?  and do you think whether it's necessary to add empty<br>
bands for bare TiO2  suraface or water adsorbed on TiO2 surfaces calculation? If no empty band added, do you think the results are<br>creditable?<br><br>Finally, I want to report something abnormal with neb restart process. When I restart the calculation of water dissociation barrier <br>
using NEB. I found all the coordinate of fixed atoms change to zero in the .path files under /tmp/ and output folder when I triggered the<br>restart command, and the restart calculation crashed after a couple of minute. So I added the all the coordinate of fixed atoms manually <br>
before the calculation was down. Then the calculation could restart normally. My QE is 4.0.3 version. Does any one encounter the same<br>problem? <br><br>thank you for reading. any hints will be deeply appreciated.<br><br>
vega<br clear="all"><br>-- <br>==================================================================================<br>Vega Lew ( weijia liu)<br>PH.D Candidate in Chemical Engineering<br>State Key Laboratory of Materials-oriented Chemical Engineering<br>
College of Chemistry and Chemical Engineering<br>Nanjing University of Technology, 210009, Nanjing, Jiangsu, China<br>******************************************************************************************************************<br>
Email: <a href="mailto:vegalew@gmail.com">vegalew@gmail.com</a><br>Office: Room A705, Technical Innovation Building, Xinmofan Road 5#, Nanjing, Jiangsu, China  <br>****************************************************************************************************************** <br>