<br />
Dear Prof. Scott Beckman<br />
<br />
Thank you for the enlightenment !!<br />
<br />
About selection of isovalue:<br />
<br />
I would like to know, if I have to compare two systems (say, Si and Ge), which iso-value should I select for both of them? (so that comparison is valid.)<br />
<br />
Thank you once again!!!<br />
<br />
Thanks & Regards<br />
Amit<br />
<br />
<br />
On Mon, 02 Aug 2010 18:11:00 +0530  wrote<br />
><br />
<br />
<br />
The charge is not calculated at EFermi, but it is the sum of the charge<br />
from the Kohn-Sham states integrated over the Brillouin zone and summed<br />
from the lowest state to the EFermi.  <br />
<br />
<br />
<br />
The integral over the<br />
charge distribution (on your XSF grid) should equal "nelec" from PWSCF. <br />
<br />
<br />
<br />
<br />
The minimum and maximum value tells you the min. and max. value on the<br />
XSF grid.  So if you have a uniform charge concentration min=max.  For<br />
Si it is likely that 0.0014 will be located somewhere in an<br />
interstitial region and 0.87 will be located in a bond-charge.  <br />
<br />
<br />
<br />
I wouldn't say that there is any "ideal" isovalue.  I would fiddle with<br />
this value and look at the charge distribution.  This will allow you to<br />
develop qualitative picture of how your system works.  <br />
<br />
<br />
<br />
<br />
<br />
<br />
<br />
<br />
<br />
<br />
<br />
Yamit wrote:<br />
Hello All<br />
<br />
  <br />
<br />
I am using Xcrysden to plot charge density of Silicon (example 05:<br />
Quantum Esspresso-4.2). <br />
<br />
*) Could anyone please tell the details of "xsf" file?<br />
<br />
( This charge density is calculated at Fermi Energy. Please also<br />
correct me, if I am wrong.)<br />
<br />
  <br />
<br />
*)When for Silicon: Minimum grid value :0.001411 and <br />
<br />
Maximum grid value : 0.87371<br />
<br />
What should be the ideal isovalue? Are there any rules to select the<br />
isovalues, from the values in *.xsf file? <br />
<br />
  <br />
<br />
(For isovalue of 0.001411(min) & 0.87371(max), no charge density<br />
has been found. )<br />
<br />
  <br />
<br />
Thank you in advance!!!<br />
<br />
  <br />
<br />
Regards<br />
<br />
Amit N. Harode<br />
<br />
CRL INDIA<br />
<br />
  <br />
<br />
P.S.<br />
<br />
Input file for QE's pp.x is as follows, &inputpp<br />
<br />
prefix = 'si'<br />
<br />
filplot = 'sicharge'<br />
<br />
plot_num= 0<br />
<br />
/<br />
<br />
&plot<br />
<br />
nfile = 1<br />
<br />
filepp(1) = 'sicharge'<br />
<br />
weight(1) = 1.0<br />
<br />
iflag = 3<br />
<br />
output_format = 5<br />
<br />
fileout = 'si.rho.xsf'<br />
<br />
/<br />
<br />
  <br />
<br />
  <br />
    <br />
      <br />
        <br />
      <br />
    <br />
  <br />
  <br />
<br />
_______________________________________________<br />
XCrySDen mailing list<br />
XCrySDen@democritos.it<br />
http://www.democritos.it/mailman/listinfo/xcrysden<br />
  <br />
<br />
<br />
<br />
_______________________________________________<br />
<br />
XCrySDen mailing list<br />
<br />
XCrySDen@democritos.it<br />
<br />
http://www.democritos.it/mailman/listinfo/xcrysden<br />
<br />
<br><Table border=0 Width=644 Height=57 cellspacing=0 cellpadding=0 style="font-family:Verdana;font-size:11px;line-height:15px;"><TR><td><A HREF="http://sigads.rediff.com/RealMedia/ads/click_nx.ads/www.rediffmail.com/signatureline.htm@Middle?" target="_blank"><IMG SRC="http://sigads.rediff.com/RealMedia/ads/adstream_nx.ads/www.rediffmail.com/signatureline.htm@Middle"></A></td></TR></Table>