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 =  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>&inputPP<br>filplot='inter(1)'<br> 
                        ...........then 'inter(2) then 'inter(3)'
       (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>&plot<br>
nfile = 3             <br>filepp(3)=inter,    <br>weight(3) =<span style="background-color: rgb(255, 102, 102);">weight(1)*rho(1)</span>  -  <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'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>