<div dir="ltr"><div>The frequencies you report are 0 by definition, almost 0 when computed: <a href="http://www.quantum-espresso.org/faq/phonons/#7.2">http://www.quantum-espresso.org/faq/phonons/#7.2</a> . No available solution other than imposing the Acoustic Sum Rule (ASR) afterwards.<br>Basically, the ASR violation is numerical noise, and as such, rather unpredictable.<br><br></div><div>Paolo<br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Sep 11, 2016 at 8:10 PM, XIAOMING ZHANG <span dir="ltr"><<a href="mailto:xiaom.zhang@utah.edu" target="_blank">xiaom.zhang@utah.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
<div style="direction:ltr;font-family:Tahoma;color:#000000;font-size:10pt">
<pre>Dear all,</pre>
<pre>I have several questions when calculating phonons by ph.x:</pre>
<pre>(1) Using the same input files, I got different output files of dynamical matrices when using different computing cores; </pre>
<pre>take q = ( 0.000000000 0.000000000 0.000000000 ) as an example:</pre>
<pre>Using 96 computing cores, the results a<span style="font-size:10pt">fter the line "Dynamical Matrix in cartesian axes"</span></pre>
<pre>and the q-value are:</pre>
<pre> ******************************<wbr>******************************<wbr>**************
freq ( 1) = -1.940214 [THz] = -64.718584 [cm-1]
( 0.000000 0.000000 -0.000000 0.000000 -0.237593 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 -0.237593 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 -0.237593 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 -0.237593 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 -0.237593 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 -0.237593 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 -0.209687 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 -0.209687 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 -0.209687 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 -0.209687 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 -0.209687 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 -0.209687 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 -0.363998 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 -0.363998 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 -0.363998 0.000000 )
freq ( 2) = 0.946381 [THz] = 31.567870 [cm-1]
( 0.000000 0.000000 -0.000000 0.000000 -0.101981 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 -0.101981 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 0.101979 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 0.101979 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 0.000001 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 0.000001 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 -0.384498 0.000000 )
( -0.000000 0.000000 -0.000000 0.000000 -0.384498 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 0.384493 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 0.384493 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 0.000005 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 0.000005 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 0.428402 0.000000 )
( -0.000000 0.000000 -0.000000 0.000000 -0.428397 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 -0.000006 0.000000 )
freq ( 3) = 0.946381 [THz] = 31.567870 [cm-1]
( 0.000000 0.000000 -0.000000 0.000000 -0.058877 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 -0.058877 0.000000 )
( 0.000000 0.000000 -0.000000 0.000000 -0.058879 0.000000 )</pre>
<pre>However, when I using 128 computing cores, the corresponding results turn to be:</pre>
<pre>******************************<wbr>******************************<wbr>**************
freq ( 1) = -1.941478 [THz] = -64.760743 [cm-1]
( -0.000000 0.000000 -0.000000 0.000000 0.237504 0.000000 )
( -0.000000 0.000000 -0.000000 0.000000 0.237504 0.000000 )
( -0.000000 0.000000 -0.000000 0.000000 0.237504 0.000000 )
( -0.000000 0.000000 -0.000000 0.000000 0.237504 0.000000 )
( -0.000000 0.000000 -0.000000 0.000000 0.237504 0.000000 )
( -0.000000 0.000000 -0.000000 0.000000 0.237504 0.000000 )
( -0.000000 0.000000 -0.000000 0.000000 0.209612 0.000000 )
( -0.000000 0.000000 -0.000000 0.000000 0.209612 0.000000 )
( -0.000000 0.000000 -0.000000 0.000000 0.209612 0.000000 )
( -0.000000 0.000000 -0.000000 0.000000 0.209612 0.000000 )
( -0.000000 0.000000 -0.000000 0.000000 0.209612 0.000000 )
( -0.000000 0.000000 -0.000000 0.000000 0.209612 0.000000 )
( -0.000000 0.000000 -0.000000 0.000000 0.364201 0.000000 )
( -0.000000 0.000000 -0.000000 0.000000 0.364201 0.000000 )
( -0.000000 0.000000 -0.000000 0.000000 0.364201 0.000000 )
freq ( 2) = -0.797967 [THz] = -26.617308 [cm-1]
( -0.217955 0.000000 -0.132999 0.000000 -0.000000 0.000000 )
( -0.217955 0.000000 -0.132999 0.000000 -0.000000 0.000000 )
( -0.212134 0.000000 -0.133209 0.000000 0.000000 0.000000 )
( -0.212134 0.000000 -0.133209 0.000000 0.000000 0.000000 )
( -0.215227 0.000000 -0.138145 0.000000 0.000000 0.000000 )
( -0.215227 0.000000 -0.138145 0.000000 0.000000 0.000000 )
( -0.230258 0.000000 -0.129133 0.000000 -0.000000 0.000000 )
( -0.230258 0.000000 -0.129133 0.000000 -0.000000 0.000000 )
( -0.205571 0.000000 -0.130026 0.000000 0.000000 0.000000 )
( -0.205571 0.000000 -0.130026 0.000000 0.000000 0.000000 )
( -0.218687 0.000000 -0.150959 0.000000 0.000000 0.000000 )
( -0.218687 0.000000 -0.150959 0.000000 0.000000 0.000000 )
( -0.228822 0.000000 -0.140474 0.000000 0.000000 0.000000 )
( -0.224086 0.000000 -0.140645 0.000000 -0.000000 0.000000 )
( -0.226602 0.000000 -0.144661 0.000000 -0.000000 0.000000 )
freq ( 3) = -0.797967 [THz] = -26.617308 [cm-1]
( -0.136570 0.000000 0.212256 0.000000 -0.000000 0.000000 )
( -0.136570 0.000000 0.212256 0.000000 -0.000000 0.000000 )
( -0.136360 0.000000 0.218077 0.000000 0.000000 0.000000 )</pre>
<pre>So, why do different computing cores lead to different results? Does someone encounter the same problems?</pre>
<pre>By the way, is it reasonable for the negative frequency (such as <span style="font-size:13.3333px">freq ( 1) = -1.940214 [THz] = -64.718584 [cm-1]</span><span style="font-size:10pt">) ?</span></pre>
<pre><span style="font-size:10pt"><br></span></pre>
<pre><span style="font-size:10pt">(2) With different computing cores, I always get the p</span>honon spectrum with small imaginary frequency no matter how I <span style="font-size:10pt">adjust the parameters suggested by other guys</span><span style="font-size:10pt">. </span></pre>
<pre><span style="font-size:10pt">S</span><span style="font-size:10pt">o can someone give me some useful suggestions to </span>eliminate<span style="font-size:10pt"> the </span><span style="font-size:13.3333px">imaginary frequency</span><span style="font-size:10pt">?</span></pre>
<pre><br></pre>
<pre>Thanks for <font size="3">your time,</font></pre>
<pre><font size="3">Xiaoming</font></pre>
<pre><span style="white-space:normal"><font size="3">Department of Materials Science and Engineering, University of Utah</font></span></pre>
<pre><span style="font-size:10pt"><br></span></pre>
</div>
</div>
<br>______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a><br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>
</div>