Hello!<div><br></div><div>I am sorry, I did not answer before!</div><div>well I saw the "pp" file (the *.charge in my case) and how it is a large file I have to use a big command line for see what are the minimum an maximum values for the elf calculation, and the range was from 0 to 1, so they are right! :) . How said Giovanni, xcrysden interpolates the data, for this reason, the values were out of range. Probably <meta http-equiv="content-type" content="text/html; charset=utf-8"><span class="Apple-style-span" style="border-collapse: collapse; font-family: arial, sans-serif; font-size: 13px; ">contour and </span><span class="Apple-style-span" style="border-collapse: collapse; font-family: arial, sans-serif; font-size: 13px; ">plotrho </span><span class="Apple-style-span" style="border-collapse: collapse; font-family: arial, sans-serif; font-size: 13px; ">output_format do the same. </span></div>
<div><span class="Apple-style-span" style="border-collapse: collapse; font-family: arial, sans-serif; font-size: 13px; "><br></span></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">Thanks for help me! :)</span></font></div>
<div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">best regards,</span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">Claudia Loyola<br>
</span></font><div><br><div class="gmail_quote">On Tue, Feb 15, 2011 at 10:16 AM, Stefano de Gironcoli <span dir="ltr"><<a href="mailto:degironc@sissa.it">degironc@sissa.it</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div text="#000000" bgcolor="#ffffff">
I see...<br>
how large is your calculation? if it is not too big, would you mind
send me (at my email) one set of input and output to have a look ? <br><font color="#888888">
stefano</font><div><div></div><div class="h5"><br>
<br>
On 02/15/2011 04:42 PM, Claudia Loyola wrote:
<blockquote type="cite">
<pre>Hello!
yes, you are right if I define the vectors in function of "a" and I define
alat=a , but let me show you how I define the vectors in the scf input file:
&SYSTEM
ibrav = 0 ,
celldm(1) =1.0 ,
.
.
.
CELL_PARAMETERS (alat=1.0)
19.109476287 -0.014917030 0.000000953
-9.542147268 16.554932043 0.000000170
0.000000706 0.000000344 14.136501037
so in this case the vector basis are the real large in bohr, and the alat is
1.0 like celldm(1). Indeed, the value of the vectors for pp.x are in
function of alat, but my alat is 1.0, so I can write directly in bohr the
values of the planes that I take, in fact I take "two replicates" in each
direction. I think this is correct too.
Thank you very much for your comments :)
regards,
Claudia Loyola
On Tue, Feb 15, 2011 at 2:01 AM, Stefano de Gironcoli <a href="mailto:degironc@sissa.it" target="_blank"><degironc@sissa.it></a>wrote:
</pre>
<blockquote type="cite">
<pre>Dear Claudia,
I may not remember correctly what the format adopted by pp.x to
define planes is... but i think it is Cartesian coordinates in unit of
alat (aka celldm(1)).
so when you write
e1(1) =38.218, e1(2)=0.0, e1(3) = 0.0,
e2(1) =0.0, e2(2)=33.108, e2(3) = 0.0,
this means presumably 38 x 33 UNIT CELLS, sampled with 100x100 points
(that is about 2-3 points in each unit cell). No surprise the FFT
interpolation used gets confused.
I think you should use more sensible values for e1(1) and e2(2)... what
is your value of celldm(1)?
At least this is my guess.
stefano
On 02/14/2011 06:57 PM, Claudia Loyola wrote:
</pre>
<blockquote type="cite">
<pre>Well, this is one of the input file:
&inputpp
&plot
nfile = 1 ,
filepp(1) = 'Ba38_33_10.charge' ,
weight(1) = 1.0 ,
iflag = 2 ,
output_format = 3 ,
fileout = 'ba.rho38_33_10.xsf' ,
x0(1) =0.0, x0(2)=0.0, x0(3) = 10.602,
e1(1) =38.218, e1(2)=0.0, e1(3) = 0.0,
e2(1) =0.0, e2(2)=33.108, e2(3) = 0.0,
nx=100, ny=100
/
</pre>
</blockquote>
<pre>_______________________________________________
Pw_forum mailing list
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a>
</pre>
</blockquote>
<pre></pre>
<pre><fieldset></fieldset>
_______________________________________________
Pw_forum mailing list
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a>
</pre>
</blockquote>
<br>
</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"><br>-- <br>Claudia Loyola Canales<br>Postdoc in Computational Physics<br>Iowa State University<br>Iowa, USA<br>http://<a href="http://www.lpmd.cl/claudial" target="_blank">www.lpmd.cl/claudial</a><div>
<a href="http://www.lpmd.cl/claudial" target="_blank"></a>http://<a href="http://www.lpmd.cl" target="_blank">www.lpmd.cl</a><br><a href="http://cosmic.mse.iastate.edu/" target="_blank">http://cosmic.mse.iastate.edu/</a><br>
http://<a href="http://www.gnm.cl" target="_blank">www.gnm.cl</a></div><br>
</div></div>