<div dir="ltr"><div><div>Dear  Lima,<br><br></div>Thank you very much for your replay.  i will go through it and get back to you.<br><br></div>Thank you.<br><div><br><br><br></div><div class="gmail_extra"><br><br><div class="gmail_quote">

On Tue, Apr 1, 2014 at 7:42 PM, Filipe Camargo Dalmatti Alves Lima <span dir="ltr"><<a href="mailto:flima@if.usp.br" target="_blank">flima@if.usp.br</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

<div dir="ltr">It is because you did not finish all the steps to analyse the data.<div><br></div><div>I got a similar problem before and I remember it is possible to obtain a file that has as the vectors of the modes that can be viewed in the xcrysden. Unfortunately I don't run simulations phonon simulations anymore, I couldn't remember everything I did in the past. <br>



</div><div><br></div><div><div>I am pasting a manual section that explains it in details. (I also change the color to red the point I think might solve your issue.)</div><div><br></div><div><br></div><div><h2 style="font-family:'Times New Roman'">



<a name="1451da359f46b9ed_SECTION00052000000000000000">4.2 Calculation of interatomic force constants in real space</a></h2><p style="font-size:medium;font-family:'Times New Roman'">First, dynamical matrices are calculated and saved for a suitable uniform grid of <b>q</b>-vectors (only those in the Irreducible Brillouin Zone of the crystal are needed). Although this can be done one <b>q</b>-vector at the time, a simpler procedure is to specify variable <tt>ldisp=.true.</tt> and to set variables <tt>nq1</tt>, <tt>nq2</tt>,<tt>nq3</tt> to some suitable Monkhorst-Pack grid, that will be automatically generated, centered at <img src="http://www.quantum-espresso.org/wp-content/uploads/Doc/ph_user_guide/img1.png" alt="$ \bf q$" height="35" width="16" align="MIDDLE" border="0"> = 0 .</p>



<p style="font-size:medium;font-family:'Times New Roman'">Second, code <tt>q2r.x</tt> reads the dynamical matrices produced in the preceding step and Fourier-transform them, writing a file of Interatomic Force Constants in real space, up to a distance that depends on the size of the grid of <b>q</b>-vectors. Input documentation in the header of <tt>PHonon/PH/q2r.f90</tt>.</p>



<p style="font-family:'Times New Roman';font-size:medium"><font color="#ff0000"><u>Program <tt>matdyn.x</tt> may be used to produce phonon modes and frequencies at any <b>q</b> using the Interatomic Force Constants file as input. Input documentation in the header of <tt>PHonon/PH/matdyn.f90</tt>.</u><br>



</font></p><p style="font-size:medium;font-family:'Times New Roman'">See Example 02 for a complete calculation of phonon dispersions in AlAs. </p><p style="font-size:medium;font-family:'Times New Roman'">



source: <a href="http://www.quantum-espresso.org/wp-content/uploads/Doc/ph_user_guide/node1.html" target="_blank">http://www.quantum-espresso.org/wp-content/uploads/Doc/ph_user_guide/node1.html</a><br></p><p style="font-size:medium;font-family:'Times New Roman'">



<br></p><p style="font-size:medium;font-family:'Times New Roman'">I hope it will solve your issue.</p><p style="font-size:medium;font-family:'Times New Roman'"><br></p><p style="font-size:medium;font-family:'Times New Roman'">



Best regards,</p><p style="font-size:medium;font-family:'Times New Roman'"><br></p><p style="font-size:medium;font-family:'Times New Roman'">Filipe</p></div></div><div class="gmail_extra"><div><div class="h5">



<br><br><div class="gmail_quote">On Tue, Apr 1, 2014 at 8:49 AM, Peram sreenivasa reddy <span dir="ltr"><<a href="mailto:peramsreenivas@gmail.com" target="_blank">peramsreenivas@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">



<div dir="ltr"><div><div><div><div>Dear xirainbow,<br><br></div>Thank you very much for your replay.<br><br></div>In that case.dyn1 file also it is giving same only as in the case.ph.out file.<br><br></div>Please give me more clarification.<br>





<br></div>Thanking you.<div><div><br><div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Apr 1, 2014 at 3:15 PM, xirainbow <span dir="ltr"><<a href="mailto:nkxirainbow@gmail.com" target="_blank">nkxirainbow@gmail.com</a>></span> wrote:<br>





<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">You can find the answer at the end of  *.dyn* files.<br>
<div><div><br>
On Mon, Mar 31, 2014 at 11:10 PM, Peram sreenivasa reddy<br>
<<a href="mailto:peramsreenivas@gmail.com" target="_blank">peramsreenivas@gmail.com</a>> wrote:<br>
> Dear Pwscf,<br>
><br>
> In my phonon calculations i got one imaginary mode in gamma to X direction.<br>
> I want to know which mode it is? and due to which atom it is happening?<br>
><br>
> I checked the case.ph.out file. In this file it is given like below.<br>
><br>
> Mode symmetry, O_h (m-3m)  point group:<br>
>   omega(  1 -  3) =        -16.6  [cm-1]   --> T_1u G_15  G_4- I<br>
>   omega(  4 -  6) =        122.3  [cm-1]   --> T_1u G_15  G_4- I<br>
>   omega(  7 -  9) =        145.0  [cm-1]   --> T_2g G_25' G_5+ R<br>
>   omega( 10 - 12) =        258.5  [cm-1]   --> T_1u G_15  G_4- I<br>
><br>
> Here i am attaching the case.agr file<br>
><br>
> Here all acoustic modes are given with same frequency. How can know that<br>
> particular mode name.<br>
><br>
><br>
>  My system is X2YZ type. How can i know this imaginary mode is due to which<br>
> atom. I plotted partial phonon density of states also. But all atoms are<br>
> giving same contribution in this acoustic mode region.<br>
><br>
>  Thank you in advance.<br>
><br>
><br>
><br>
</div></div>> _______________________________________________<br>
> Pw_forum mailing list<br>
> <a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br>
> <a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
<span><font color="#888888"><br>
<br>
<br>
--<br>
____________________________________<br>
Hui Wang<br>
School of physics, Henan University of Science and Technology, Henan, China<br>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
</font></span></blockquote></div><br><br clear="all"><br><br></div></div></div></div>
<br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br><br clear="all"><div><br></div>-- <br></div></div><div dir="ltr">_________________________________________<span class="HOEnZb"><font color="#888888"><br>



Filipe Camargo Dalmatti Alves Lima<br>PhD Student<br>University of São Paulo, Physics Institute, Materials Physics Department, Nanomol Group, Brazil.<br>Phones:    (11) 3091-6881  (USP)<br>                (11) 97408-2755 (Vivo)<br>



</font></span></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://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br><br clear="all"><br><br></div></div>