Dear PWscf,<br><br>I wish to get some help regarding the calculation of the charge difference.<br>For example, if H2O was adsorbed onto TiO2 i want to see the charge difference this is how i would go about it<br><br>
<span style="font-weight: bold;">Charge density of </span><br>a. (H2O + TiO2) <br>b. (TiO2)<br>c. (H2O)<br>Charge density difference =&nbsp; a - (b+c)<br><br><span style="font-weight: bold;">So i would save the the a, b, c, intermediate files all in 
e.g.</span><br>separate calculation<br><br>&amp;inputPP<br>filplot=&#39;inter(1)&#39;<br>&nbsp;
&nbsp; &nbsp; &nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ...........then &#39;inter(2) then &#39;inter(3)&#39;
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (i guess they will all have to be in different names? inter(1),
inter(2), inter(3) as a, b, c, respectively)
<br><br><span style="font-weight: bold;">Now, here is where i am lost.....</span><br>&amp;plot<br>
nfile = 3 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; <br>filepp(3)=inter, &nbsp;&nbsp; <br>weight(3) =<span style="background-color: rgb(255, 102, 102);">weight(1)*rho(1)</span>&nbsp; -&nbsp; <span style="font-weight: bold;">(</span><span style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 51);">

weight(2)*rho(2)</span>+<span style="background-color: rgb(51, 102, 255);">weight(3)*rho(3)</span><span style="font-weight: bold;">)</span>,<br><br>i don&#39;t know what to do in this section, in regards to how to use weight(i), filepp(i)
<br><br><br>Thank you very much,<br><br>regards, Ken.<br>