<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Dear Zhen,<div><br></div><div>all of your very legitimate questions have the same answer: aliasing. Real-space force constants are obtained by discrete Fourier transform (FT) from reciprocal-space dynamical matrices. Discrete FTs are periodic, with period {2\pi\over \Delta q}, with hopefully obvious notation. Periodicity tell you what to do with negative indeces in the force constants. If for some physical reason the force constants are even by parity, this should also explain the symmetry you have observed.</div><div><br></div><div>EXERCISE: think about the above and try to answer yourself your questions in as clear a way as you can. I will be glad to read your own explanation on this forum.</div><div><br></div><div>Regards - SB</div><div><br></div><div>PS: please, always indicate the real subject of your post in the subject field (do not limit yourself to hit "reply" to a digest daily dispatch of the forum's posts)</div><div><br><div><div>On Jan 5, 2010, at 4:02 AM, Zhen Huang wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite">Dear Prof. Eyvaz Isaev and users-<br><br>This is Zhen. Thanks for you information. However, I have couple of uncleared questions to follow up. Below is part of the results a friend of mine provided which are the force constants of Au who has FCC structure. Two things made me confused. First, according to your information R=(i-1)*a + (j-1)*b + (k-1)*c, I found that 2 1 1's force constant and 4 1 1's force constants are the same (marked by a line of @ after below). Why is this? On the other hand, 2 2 1 should represent the force constants of x, x between an atom and its third nearest neighbor (distance is (3/2)^0.5a where a is the lattice constant) atom and 3 1 1 represent the force constant between an atom and its fourth nearest neighbor (2^0.5a) (both are marked by $ below). Why the force constant between the atoms and its third nearest neighbor is smaller than the one between it and its fourth nearest's. Is there a particular reason for it? <br>
<br>On the other hand, how do I find the R= -1*a+x*b+y*c. Seems based on R=(i-1)*a + (j-1)*b + (k-1)*c there will never be a negative index. <br><br><br><br>*******************************<br>  1    1  2  7.8423629  0.0000000  0.0000000  0.0000000  0.0000000  0.0000000<br>
           1  'Au '    179527.171937097     <br>    1    1      0.0000000      0.0000000      0.0000000<br> F<br>   4   4   4<br>   1   1   1   1<br>   1   1   1   1.26759193594E-01<br>   2   1   1  -1.87445789062E-02                @@@@@@@@<br>
   3   1   1   1.66957109375E-03                $$$$$$$$$$$$$$<br>   4   1   1  -1.87445789062E-02                @@@@@@@@<br>   1   2   1   6.90391046875E-03<br>   2   2   1  -4.25726562500E-05                $$$$$$$$$$$$$$<br>
   3   2   1  -4.96455781250E-04<br>   4   2   1  -1.87445789063E-02<br>   1   3   1  -2.07408906250E-04<br>   2   3   1  -4.25726562500E-05<br>   3   3   1   1.66957109375E-03<br>   4   3   1  -4.25726562500E-05<br>   1   4   1   6.90391046875E-03<br>
   2   4   1  -1.87445789063E-02<br>   3   4   1  -4.96455781250E-04<br>   4   4   1  -4.25726562500E-05<br>   1   1   2  -1.87445789062E-02<br>******************************<br><br clear="all">Best Regards<br>-------<br>
Zhen (Alex) Huang<br>Ph.D. Student<br>Nanoscale Transport Research Group<br>Laboratory for Computational Methods in Emerging Technologies <br>Cooling Technologies Research Center<br>School of Mechanical Engineering<br>Purdue University<br>
Tel: 765 237 9733<br>
<br><br><div class="gmail_quote">On Sun, Jan 3, 2010 at 2:58 AM,  <span dir="ltr"><<a href="mailto:pw_forum-request@pwscf.org">pw_forum-request@pwscf.org</a>></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;">
Send Pw_forum mailing list submissions to<br>
        <a href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
        <a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
or, via email, send a message with subject or body 'help' to<br>
        <a href="mailto:pw_forum-request@pwscf.org">pw_forum-request@pwscf.org</a><br>
<br>
You can reach the person managing the list at<br>
        <a href="mailto:pw_forum-owner@pwscf.org">pw_forum-owner@pwscf.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of Pw_forum digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
   1. Re: A question about force constants calculation output<br>
      (Eyvaz Isaev)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Sat, 2 Jan 2010 09:04:41 -0800 (PST)<br>
From: Eyvaz Isaev <<a href="mailto:eyvaz_isaev@yahoo.com">eyvaz_isaev@yahoo.com</a>><br>
Subject: Re: [Pw_forum] A question about force constants calculation<br>
        output<br>
To: PWSCF Forum <<a href="mailto:pw_forum@pwscf.org">pw_forum@pwscf.org</a>><br>
Message-ID: <<a href="mailto:980073.51698.qm@web65703.mail.ac4.yahoo.com">980073.51698.qm@web65703.mail.ac4.yahoo.com</a>><br>
Content-Type: text/plain; charset=iso-8859-1<br>
<br>
Dear Zhen,<br>
<br>
In order to understand correctly the content of these files  you should read a textbook on the lattice dynamics theory of semiconductors.   Textbook "Fundamentals of Semiconductors" by Peter Yu and Manuel  Cardona is a very good starting point.<br>

<br>
--- On Wed, 12/30/09, Zhen Huang <<a href="mailto:huang87@purdue.edu">huang87@purdue.edu</a>> wrote:<br>
<br>
> ph.x calculation generated files alas.dynX where X goes from<br>
> 0 to 8. In the file alas.dyn0 it gives the k points.. Seems<br>
> in the other dynX files for example .dyn1 it gives force<br>
> constants of q=(0.0, 0.0, 0.0)? and .dyn2 gives force<br>
> constants of q=(-0.25?? 0.25? -0.25)<br>
<br>
This is correct. 4x4x4 set of q-points gives 8 q-points for which dynamical matrices have to be calculated. These q-points are written in dyn0 file. Files dyn1 - dyn8 are dynamical matrices.<br>
<br>
> What is the T matrix? and Matrix 1 and 2 labeled below.<br>
<br>
There is no T-matrix at all, "T" means .true., i.e. AlAs, in fact, is semi-conducting.<br>
<br>
Next 3x3 matrix is the macroscopic dielectric constants matrix.<br>
<br>
"Matrix 1" and "Matrix 2" contain effective Born charges for Al (atom 1) and As (atom 2).<br>
<br>
For metals you find "F" (.false.) and no dielectric matrix and Born effective charges.<br>
<br>
> Also, the questions are what is 4 4 4 representing the the force<br>
> constants information part  labeled below,? what does the 1 1 1 1 mean > on the second line and also what is the meaning of the index?<br>
> 2?? 1?? 1? -6.26636620736E-03 (for example), I assume -6.266e-03 is<br>
> the force constants but what is 2 1 1.<br>
<br>
You asked this question before. Please have a look at<br>
<br>
<a href="http://www.democritos.it/pipermail/pw_forum/2008-September/010099.html" target="_blank">http://www.democritos.it/pipermail/pw_forum/2008-September/010099.html</a><br>
<br>
The only correction is that R=(i-1)*a + (j-1)*b + (k-1)*c (thanks Stefano d.G).<br>
<br>
> If the forth column is the force constants why it is different from<br>
> Ph.x  calcualtion's output.<br>
<br>
Please do not mix dynamical matrix with IFC. The first one is in q-space, and the second one is in R-space.<br>
<br>
> One more question, if I want to look at the force constants<br>
> between a unit cell and its closest nearest neibhor unit<br>
> cell do I look at the dynX output with k=(0, -1, 0), (-1, 0,<br>
> 0) and such and second nearest unit cells will be in dynX of<br>
> k=(0, -.5, 0), (-.5, 0, 0)?<br>
<br>
This is unclear what do you mean as IFC is for atoms, but not for unit cells. Just keep in mind in the Linear Response theory all calculations are performed inside the unit cell. To calculate the IFC range see a formula for R given above.<br>

<br>
Bests,<br>
Eyvaz.<br>
<br>
 -------------------------------------------------------------------<br>
Prof. Eyvaz Isaev,<br>
Theoretical Physics Department, Moscow State Institute of Steel & Alloys, Russia,<br>
Department of Physics, Chemistry, and Biology (IFM), Linkoping University, Sweden<br>
Condensed Matter Theory Group, Uppsala University, Sweden<br>
<a href="mailto:Eyvaz.Isaev@fysik.uu.se">Eyvaz.Isaev@fysik.uu.se</a>, <a href="mailto:isaev@ifm.liu.se">isaev@ifm.liu.se</a>, <a href="mailto:eyvaz_isaev@yahoo.com">eyvaz_isaev@yahoo.com</a><br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
------------------------------<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>
End of Pw_forum Digest, Vol 31, Issue 2<br>
***************************************<br>
</blockquote></div><br>
_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>http://www.democritos.it/mailman/listinfo/pw_forum<br></blockquote></div><br><div>
<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">---</span></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">Stefano Baroni - SISSA</span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">  </span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">&</span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">  </span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">DEMOCRITOS National Simulation Center - Trieste</span></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; "><a href="http://stefano.baroni.me">http://stefano.baroni.me</a> [+39] 040 3787 406 (tel) -528 (fax) / stefanobaroni (skype)</span></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 12px/normal Helvetica; min-height: 14px; font-size: 12px; "><br style="font-size: 12px; "></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 12px/normal Helvetica; min-height: 14px; font-size: 12px; "><span class="Apple-style-span" style="font-size: 14px; "><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">La morale est une logique de l'action comme la logique est une morale de la pensée - Jean Piaget</span></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; "><br></span></font></div></span></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">Please, if possible, don't</span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">  </span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">send me MS Word or PowerPoint attachments</span></font></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">Why? See:</span></font><font class="Apple-style-span" size="3"><span class="Apple-style-span" style="font-size: 12px; ">  </span></font><font class="Apple-style-span" color="#0023E9" size="3"><span class="Apple-style-span" style="font-size: 12px; "><a href="http://www.gnu.org/philosophy/no-word-attachments.html">http://www.gnu.org/philosophy/no-word-attachments.html</a></span></font></div><div><font class="Apple-style-span" color="#0023E9" size="3"><span class="Apple-style-span" style="font-size: 12px; "><br class="webkit-block-placeholder"></span></font></div></div></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"></span><br class="Apple-interchange-newline">
</div>
<br></div></body></html>