Hi, Paul<br>I think this two choice wonot make difference for FCC crystal.&nbsp; <br>And I do agree with your last statement .<br>Thanks&nbsp; Paolo and Axel, and Happy new year !<br><br>best wishes,<br>hai-ping<br><br><div class="gmail_quote">
On Tue, Dec 30, 2008 at 12:56 PM, Paul M. Grant <span dir="ltr">&lt;<a href="mailto:w2agz@pacbell.net">w2agz@pacbell.net</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
OK, let&#39;s get specific...suppose I want to generate a uniformly distributed<br>
set of k-points for the reciprocal lattice of the primitive cell of an fcc<br>
lattice. &nbsp;The output of an scf run (and common sense matrix algebra) tells<br>
me the reciprocal lattice is defined as<br>
<br>
reciprocal axes: (cart. coord. in units 2 pi/a_0)<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; b(1) = ( &nbsp;1.000000 &nbsp;0.000000 &nbsp;0.000000 )<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; b(2) = ( &nbsp;0.000000 &nbsp;1.000000 &nbsp;0.000000 )<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; b(3) = ( &nbsp;0.000000 &nbsp;0.000000 &nbsp;1.000000 )<br>
<br>
OK (again)...I now generate a(n) uniformly distributed grid over the &quot;b&quot;<br>
vectors, and plot the points via xcrysden with PRIMVEC as above and see that<br>
grid fills the entire basic reciprocal cell (that is, with translation<br>
vector = O). &nbsp;No specification of celldm(1).<br>
<br>
Now, I want to input this grid to an nscf run via K_POINTS...is my grid<br>
&quot;tpiba&quot; or &quot;crystal?&quot;<br>
<br>
BTW, we (this community) owe a huge debt to Paolo and Axel.<br>
<div class="Ih2E3d"><br>
Paul M. Grant, PhD<br>
Principal, W2AGZ Technologies<br>
Visiting Scholar, Applied Physics, Stanford (2005-2008)<br>
EPRI Science Fellow (Retired)<br>
IBM Research Staff Member Emeritus<br>
<a href="mailto:w2agz@pacbell.net">w2agz@pacbell.net</a><br>
<a href="http://www.w2agz.com" target="_blank">http://www.w2agz.com</a><br>
&nbsp;<br>
&nbsp;<br>
<br>
</div><div><div></div><div class="Wj3C7c">-----Original Message-----<br>
From: <a href="mailto:pw_forum-bounces@pwscf.org">pw_forum-bounces@pwscf.org</a> [mailto:<a href="mailto:pw_forum-bounces@pwscf.org">pw_forum-bounces@pwscf.org</a>] On<br>
Behalf Of <a href="mailto:giannozz@democritos.it">giannozz@democritos.it</a><br>
Sent: Monday, December 29, 2008 12:29 PM<br>
To: <a href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</a><br>
Subject: Re: [Pw_forum] K_POINTS, tpiba or crystal<br>
<br>
Quoting &quot;Paul M. Grant&quot; &lt;<a href="mailto:w2agz@pacbell.net">w2agz@pacbell.net</a>&gt;:<br>
<br>
&gt; Frankly, I don&#39;t understand what &quot;crystal&quot; means wrt to reciprocal space<br>
<br>
it means &quot;in units of the reciprocal lattice primitive vectors&quot;<br>
<br>
Paolo<br>
<br>
<br>
----------------------------------------------------------------<br>
This message was sent using IMP, the Internet Messaging Program.<br>
<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>
_______________________________________________<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>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>Hai-Ping Lan <br>Department of Electronics ,<br>Peking University , Bejing, 100871<br><a href="mailto:lanhaiping@gmail.com">lanhaiping@gmail.com</a>, <a href="mailto:hplan@pku.edu.cn">hplan@pku.edu.cn</a><br>