Dear users,<br>I am doing the phonon calculation of Graphene sheet.<br>In the phonon input file "<span style="color: rgb(51, 51, 255);"><a href="http://graph.ph.in">graph.ph.in</a></span>" i am using <span style="color: rgb(51, 102, 255);">nq(i)= 8x8x1</span>, which gives the following q-points<br>
<br clear="all"> <b><span style="color: rgb(255, 0, 0);">Dynamical matrices for ( 8, 8, 1,) uniform grid of q-points</span><br style="color: rgb(255, 0, 0);"><span style="color: rgb(255, 0, 0);"> ( 10q-points):</span><br style="color: rgb(255, 0, 0);">
<span style="color: rgb(255, 0, 0);"> N xq(1) xq(2) xq(3) </span><br style="color: rgb(255, 0, 0);"><span style="color: rgb(255, 0, 0);"> 1 0.000000000 0.000000000 0.000000000</span><br style="color: rgb(255, 0, 0);">
<span style="color: rgb(255, 0, 0);"> 2 0.000000000 0.144337567 0.000000000</span><br style="color: rgb(255, 0, 0);"><span style="color: rgb(255, 0, 0);"> 3 0.000000000 0.288675135 0.000000000</span><br style="color: rgb(255, 0, 0);">
<span style="color: rgb(255, 0, 0);"> 4 0.000000000 0.433012702 0.000000000</span><br style="color: rgb(255, 0, 0);"><span style="color: rgb(255, 0, 0);"> 5 0.000000000 -0.577350269 0.000000000</span><br style="color: rgb(255, 0, 0);">
<span style="color: rgb(255, 0, 0);"> 6 0.125000000 0.216506351 0.000000000</span><br style="color: rgb(255, 0, 0);"><span style="color: rgb(255, 0, 0);"> 7 0.125000000 0.360843918 0.000000000</span><br style="color: rgb(255, 0, 0);">
<span style="color: rgb(255, 0, 0);"> 8 0.125000000 0.505181486 0.000000000</span><br style="color: rgb(255, 0, 0);"><span style="color: rgb(255, 0, 0);"> 9 0.250000000 0.433012702 0.000000000</span><br style="color: rgb(255, 0, 0);">
<span style="color: rgb(255, 0, 0);"> 10 0.250000000 0.577350269 0.000000000</span></b><br><br>and in the <a href="http://matdyn1.in">matdyn1.in</a> file (input file for <span style="color: rgb(51, 51, 255);">matdyn.x</span>) i used <br>
<b><span style="color: rgb(255, 0, 0);">G 0.0 0.0 0.0</span><br style="color: rgb(255, 0, 0);"><span style="color: rgb(255, 0, 0);">K 0.3333 0.57735 0.0000</span><br style="color: rgb(255, 0, 0);"><span style="color: rgb(255, 0, 0);">M 0.0000 0.57735 0.0000</span><br style="color: rgb(255, 0, 0);">
<span style="color: rgb(255, 0, 0);">G 0.0000 0.00000 0.0000</span></b><br><br>high symmetry points, along with other intermediate point.<br><br>So my question is, when i plot the frequencies (<span style="color: rgb(51, 51, 255);">*.freq</span> file), given after running matdyn.x file, using the above high symmetry points, the dispersion is all positive (the frequencies are all positive)<br>
<br>But when I am doing single <b><span style="color: rgb(255, 0, 0);">q=0.3333 0.57735 0.0000</span></b> point calculation in my <span style="color: rgb(51, 51, 255);"><a href="http://graph.ph.in">graph.ph.in</a></span> file. For this perticular q-point the frequencies are -ve.<br>
<br>Can anybody tell me the reason for this, or I am doing some error?<br><br>-- <br>Bipul Rakshit<br>Research Fellow<br>S N Bose Centre for Basic Sciences,<br>Salt Lake,<br>Kolkata 700 098<br>India<br>