Dear Gabriele,<div><br><div><div>The attachment is the charge density given by pp.x. Left pic is no doping and the right pic is N doped. These pictures present differently, so I think the charge density shouldn't be zero.</div>
<div><br></div><div><br></div><div><br><br><div class="gmail_quote">2012/9/19 Gabriele Sclauzero <span dir="ltr"><<a href="mailto:gabriele.sclauzero@epfl.ch" target="_blank">gabriele.sclauzero@epfl.ch</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><div><br></div>It seems to me that you have not completely understood how to use pp.x.<div>
If you take the difference between SRO-up.rho and SROup-doping.rho obtained with the inputs here below you will get zero, because they will contain the same charge density. You see why?</div><div><br></div><div>GS<br><div>
<br></div><div><br><div><div>Il giorno 19/set/2012, alle ore 10.15, Yue-Wen Fang ha scritto:</div><div><div class="h5"><br><blockquote type="cite">The input files are listed here. In fact, there is no difference.<br><br><div class="gmail_quote">
2012/9/19 Yue-Wen Fang <span dir="ltr"><<a href="mailto:yuewen.fang@gmail.com" target="_blank">yuewen.fang@gmail.com</a>></span><br>
<blockquote class="gmail_quote" style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div><br></div>
<div>&inputpp</div><div>prefix = 'SRO'</div><div>outdir = './'</div><div>filplot = 'SROcharge'</div>
<div>plot_num= 0</div><div>spin_component=1</div><div>/</div><div>&plot</div>
<div>nfile = 1</div><div>filepp(1) = 'SROcharge'</div><div>weight(1) = 1.0</div><div>iflag = 2</div><div>output_format = 2</div><div>fileout = 'SRO-up.rho'</div><div><div>e1(1) =1.0, e1(2)=0.0, e1(3) = 0.0,</div>


<div>e2(1) =0.0, e2(2)=0.0, e2(3) = 1.0,</div><div>x0(1) =0.0, x0(2)=0.50, x0(3) = 0.0,</div><div>nx=56, ny=40</div><div>/</div><div><br></div></div><div><div>&inputpp</div><div>prefix = 'SRO'</div><div>outdir = './'</div>


<div>filplot = 'SROcharge'</div><div>plot_num= 0</div><div>spin_component=1</div><div>/</div><div>&plot</div><div>nfile = 1</div><div>filepp(1) = 'SROcharge'</div><div>weight(1) = 1.0</div><div>iflag = 2</div>


<div>output_format = 2</div><div>fileout = 'SROup-doping.rho'</div><div><div>e1(1) =1.0, e1(2)=0.0, e1(3) = 0.0,</div><div>e2(1) =0.0, e2(2)=0.0, e2(3) = 1.0,</div><div>x0(1) =0.0, x0(2)=0.50, x0(3) = 0.0,</div>
<div>nx=56, ny=40</div>
<div>/</div><div><br></div></div></div><div><div><br><div class="gmail_quote">2012/9/19 Gabriele Sclauzero <span dir="ltr"><<a href="mailto:gabriele.sclauzero@epfl.ch" target="_blank">gabriele.sclauzero@epfl.ch</a>></span><br>

<blockquote class="gmail_quote" style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<div style="word-wrap:break-word"><div><br></div><div>I mean the header (the first 4 lines at least) of the density files that you are supplying to pp.x to make the plot ( filepp(1) and filepp(2) ), not the xsf files.</div>


<div>Those lines contain information about the fft grid and are read by pp.x to check consistency between the grids in the two files.</div><div><br></div><div>GS </div></div></blockquote></div></div></div></blockquote></div>
</blockquote></div></div></div><br><div class="im"><div>
<span style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:auto;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:medium;white-space:normal;font-family:Helvetica;word-spacing:0px"><span style="text-indent:0px;letter-spacing:normal;font-variant:normal;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:medium;white-space:normal;font-family:Helvetica;word-spacing:0px"><span style="text-indent:0px;letter-spacing:normal;font-variant:normal;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:medium;white-space:normal;font-family:Helvetica;word-spacing:0px"><div style="word-wrap:break-word">
<span style="text-indent:0px;letter-spacing:normal;font-variant:normal;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:medium;white-space:normal;font-family:Helvetica;word-spacing:0px"><div style="word-wrap:break-word">
<span style="text-indent:0px;letter-spacing:normal;font-variant:normal;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:medium;white-space:normal;font-family:Helvetica;word-spacing:0px"><div style="word-wrap:break-word">
<span style="text-indent:0px;letter-spacing:normal;font-variant:normal;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:medium;white-space:normal;font-family:Helvetica;word-spacing:0px"><div style="word-wrap:break-word">
<span style="text-indent:0px;letter-spacing:normal;font-variant:normal;font-style:normal;font-weight:normal;line-height:normal;border-collapse:separate;text-transform:none;font-size:medium;white-space:normal;font-family:Helvetica;word-spacing:0px"><div style="word-wrap:break-word">
<div><span style="color:rgb(126,126,126);font-size:16px;font-style:italic"><br>§ Gabriele Sclauzero, EPFL SB ITP CSEA</span></div><div><font color="#7E7E7E"><i>   PH H2 462, Station 3, CH-1015 Lausanne</i></font></div></div>
</span><br></div></span><br></div></span><br></div></span><br></div></span><br></span><br></span>
</div>
<br></div></div></div></div><br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
<br></blockquote></div><br><br clear="all"><div><br></div>-- <br>----<div>Yue-Wen Fang</div><div><br></div><div><br></div><br>
</div></div></div>