<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style></head>
<body class='hmmessage'><div dir='ltr'>
Dear all,<div><br></div><div>Sorry to ask so many questions about computing phonon frequencies but recently i have got some strange results for one of the systems I am working on. My work involves computing phonon frequencies for a 42 carbon at the Gamma point so in total 126 irreducible representations. in a previous calculations, my input was:</div><div><br></div><div><u>SCF FILE</u></div><div><span lang="PT-BR" style="font-size: 10pt; font-family: 'Courier New'; "> </span><span style="font-size: 10pt; font-family: 'Courier New'; "> </span></div><div>
<p class="MsoPlainText"><span lang="PT-BR" style="font-family:"Courier New""> /<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="PT-BR" style="font-family:"Courier New""> &system
<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="PT-BR" style="font-family:"Courier New""> ibrav= 0, celldm(1) =1.88972599, nat=42,
ntyp= 1, ecutwfc =32, ecutrho=128, occupations='smearing', smearing='mp',degauss=0.01<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="PT-BR" style="font-family:"Courier New""> </span></p>
<p class="MsoPlainText"><span lang="PT-BR" style="font-family:"Courier New"">/<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="PT-BR" style="font-family:"Courier New""> &electrons<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="PT-BR" style="font-family:"Courier New""> conv_thr=1.D-6, <o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="PT-BR" style="font-family:"Courier New""> mixing_beta=0.7D0,<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="PT-BR" style="font-family:"Courier New""> diago_david_ndim=2,<o:p></o:p></span></p><p class="MsoPlainText"><span lang="PT-BR" style="font-family:"Courier New""><br></span></p><p class="MsoPlainText"><span lang="PT-BR" style="font-family:"Courier New""><p class="MsoPlainText"><span lang="PT-BR"> <o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="PT-BR">ATOMIC_SPECIES<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="PT-BR"> C
12.0107 C.pbe-rrkjus.UPF<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="PT-BR"><o:p> ...</o:p></span></p><p class="MsoPlainText"><span lang="PT-BR"><o:p><br></o:p></span></p><p class="MsoPlainText"><span lang="PT-BR"><o:p><u>PH FILE:</u></o:p></span></p><p class="MsoPlainText"><br></p><p class="MsoPlainText"><p class="MsoPlainText"><span lang="PT-BR">Phonon
calculation for GB60 at gamma point<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="PT-BR"> &inputph<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="PT-BR"> tr2_ph=1.0d-15,<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="PT-BR"> prefix='GrBo42ph',<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="PT-BR"> fildvscf='GB60dvscf',<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="PT-BR"> elph=.true.,<o:p></o:p></span></p><p class="MsoPlainText"><span lang="PT-BR">.....</span></p><p class="MsoPlainText"><span lang="PT-BR"><br></span></p><p class="MsoPlainText"><span lang="PT-BR">The results obtained were:</span></p><p class="MsoPlainText"><span lang="PT-BR"><br></span></p><p class="MsoPlainText"><span lang="PT-BR"><p class="MsoPlainText"><span lang="PT-BR">omega(
1) = -4.450324 [THz] = -148.446823 [cm-1]<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="PT-BR"> omega( 2) = -3.144068 [THz] = -104.874835 [cm-1]<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="PT-BR"> omega( 3) = 2.963314 [THz] = 98.845504 [cm-1]<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="PT-BR"> omega( 4) = 3.454797 [THz] = 115.239637 [cm-1]<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="PT-BR"> omega( 5) = 3.908880 [THz] = 130.386191 [cm-1]<o:p></o:p></span></p><p class="MsoPlainText">.</p><p class="MsoPlainText">.</p><p class="MsoPlainText">.</p><p class="MsoPlainText"><p class="MsoPlainText"><span lang="PT-BR"> omega(**) = 47.380063 [THz] = 1580.428796 [cm-1]<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="PT-BR"> omega(**) = 47.573318 [THz] = 1586.875087 [cm-1]<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="PT-BR"> omega(**) = 47.821300 [THz] = 1595.146853 [cm-1]<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="PT-BR"> omega(**) = 48.119632 [THz] = 1605.098137 [cm-1]<o:p></o:p></span></p><p class="MsoPlainText"><span lang="PT-BR"><br></span></p><p class="MsoPlainText"><span lang="PT-BR">Where the 2 negative frequencies,as professor Gionazzi repeated several times, could be due to the violation of the ASR rule.</span></p><p class="MsoPlainText"><span lang="PT-BR"><br></span></p><p class="MsoPlainText">For my own piece of mind I repeated the calculation with a better accuracy where ecutrho=256 (8 times ecutwfc which should be the case if an ultrasoft potential is used)and tr2_ph =10^(-17) although mixing_beta was set equal to 0.3 this time. I expected better results but what I got was the following:</p><p class="MsoPlainText"><br></p><p class="MsoPlainText"><p class="MsoNormal"><span lang="PT-BR">omega(
1 - 1) = -77614121.2
[cm-1] --> A' I+R</span></p>
<p class="MsoNormal"><span lang="PT-BR">
omega( 2 - 2) =
-13157144.1 [cm-1] --> A' I+R</span></p>
<p class="MsoNormal"><span lang="PT-BR">
omega( 3 - 3) =
-5577066.3 [cm-1] --> A' I+R</span></p>
<p class="MsoNormal"><span lang="PT-BR">
omega( 4 - 4) =
-4332492.8 [cm-1] --> A' I+R</span></p>
<p class="MsoNormal"><span lang="PT-BR">
omega( 5 - 5) =
-4094962.7 [cm-1] --> A' I+R</span></p>
<p class="MsoNormal"><span lang="PT-BR">
omega( 6 - 6) =
-3368158.6 [cm-1] --> A' I+R</span></p>
<p class="MsoNormal"><span lang="PT-BR">
omega( 7 - 7) =
-2781643.0 [cm-1] --> A' I+R</span></p>
<p class="MsoNormal"><span lang="PT-BR">
omega( 8 - 8) =
-2403295.2 [cm-1] --> A' I+R</span></p>
<p class="MsoNormal"><span lang="PT-BR">
omega( 9 - 9) =
-1777809.3 [cm-1] --> A' I+R</span></p>
<p class="MsoNormal"><span lang="PT-BR">
omega( 10 - 10) =
82.5 [cm-1] --> A'' I+R</span></p>
<p class="MsoNormal"><span lang="PT-BR">
omega( 11 - 11) =
85.8 [cm-1] --> A'' I+R</span></p><p class="MsoNormal"><span lang="PT-BR">.</span></p><p class="MsoNormal"><span lang="PT-BR">.</span></p><p class="MsoNormal"><span lang="PT-BR">.</span></p><p class="MsoNormal"><span lang="PT-BR"><p class="MsoNormal"><span lang="PT-BR"> omega(120 -120) =
2301243.8 [cm-1] --> A' I+R</span></p>
<p class="MsoNormal"><span lang="PT-BR">
omega(121 -121) =
2494732.0 [cm-1] --> A' I+R</span></p>
<p class="MsoNormal"><span lang="PT-BR">
omega(122 -122) =
2788571.0 [cm-1] --> A' I+R</span></p>
<p class="MsoNormal"><span lang="PT-BR">
omega(123 -123) =
3135326.3 [cm-1] --> A' I+R</span></p>
<p class="MsoNormal"><span lang="PT-BR">
omega(124 -124) =
3630869.3 [cm-1] --> A' I+R</span></p>
<p class="MsoNormal"><span lang="PT-BR">
omega(125 -125) =
8808450.4 [cm-1] --> A' I+R</span></p><p class="MsoNormal"><span lang="PT-BR"><br></span></p><p class="MsoNormal"><span lang="PT-BR">What is the reason for those LARGE negative values at the beginning and the LARGE ones at the end? Is it because mixing_beta was changed?</span></p><p class="MsoNormal"><span lang="PT-BR"><br></span></p><p class="MsoNormal"><span lang="PT-BR">Thanks for the help</span></p><p class="MsoNormal"><span lang="PT-BR"><br></span></p><p class="MsoNormal"><span lang="PT-BR">Elie Koukaes</span></p><p class="MsoNormal"><span lang="PT-BR">University of Nottingham</span></p><p class="MsoNormal">NG7 2RD</p><p class="MsoNormal">UK</p></span></p></p></p><p class="MsoPlainText"><span lang="PT-BR"><br></span></p></span></p><p class="MsoPlainText"><span lang="PT-BR"><br></span></p></p></span></p></div> </div></body>
</html>