<div dir="ltr"><div><div><div><div><div>Dear all,<br><br></div><font size="2"><span style="font-family:arial,helvetica,sans-serif">To show interface dipoles, I am trying to calculate charge density for a 2D combined system containing Ag adsorbed on graphene sheet, like what is done in a sample attached (rho_tot). <span>The
vacuum space is provided along the z axis. </span> Is it right to do the following steps?</span></font><br><br></div>1- Running a SCF<br><br></div>2- Performing a PP calculation:<br></div>    &inputpp<br></div>          plot_num=0,<br><div><div>    &plot<br>          nfile=1<br>          filepp(1)='sys.charge'<br>          iflag=2,<br>          output_format=3,<br>          e1(1)=1.0, e1(2)=0.0, e1(3)=0.0,<br>          e2(1)=0.0, e2(2)=1.0, e2(3)=0.0,<br>          x0(1)=0.0, x0(2)=0.0, x0(3)=0.0,<br>          nx=40, ny=40<br>          fileout='sys.charge001.dat'<br>     /<br><br>========================<br><br></div><div>What if I set   iflag=3,  and finally run average.x ?<br><br></div><div>Any help will be greatly appreciated.<br><br></div><div>Best,<br><br></div><div>Mohammad,<br><br></div><div>Shahid Chamran University of Ahvaz<br></div></div></div>