<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=us-ascii"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:DengXian;
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:"\@DengXian";
panose-1:2 1 6 0 3 1 1 1 1 1;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=EN-US link="#0563C1" vlink="#954F72" style='word-wrap:break-word'><div class=WordSection1><p class=MsoNormal>Dear all,<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>I encountered one problem when doing phonon calculations. <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>I tried to do the test calculation provided in the ‘qe-6.8/PHonon/examples/GRID_example’. <o:p></o:p></p><p class=MsoNormal>After the whole calculations finished, the phonon frequencies of Gamma point are shown below:<o:p></o:p></p><p class=MsoNormal>###################################################<o:p></o:p></p><p class=MsoNormal>&plot nbnd= 6, nks= 161 /<o:p></o:p></p><p class=MsoNormal> 0.000000 0.000000 0.000000<o:p></o:p></p><p class=MsoNormal> -0.0000 -0.0000 -0.0000 375.5206 375.5206 410.5615<o:p></o:p></p><p class=MsoNormal>…<o:p></o:p></p><p class=MsoNormal>…<o:p></o:p></p><p class=MsoNormal>###################################################<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>The input file is <o:p></o:p></p><p class=MsoNormal>###################################################<o:p></o:p></p><p class=MsoNormal>&input<o:p></o:p></p><p class=MsoNormal> asr='simple',<o:p></o:p></p><p class=MsoNormal> flfrc='AlAs.fc',<o:p></o:p></p><p class=MsoNormal> flfrq='AlAs.freq',<o:p></o:p></p><p class=MsoNormal> flvec='AlAs.modes',<o:p></o:p></p><p class=MsoNormal> q_in_band_form=.true.,<o:p></o:p></p><p class=MsoNormal> amass(1)=26.98,<o:p></o:p></p><p class=MsoNormal> amass(2)=74.92,<o:p></o:p></p><p class=MsoNormal> /<o:p></o:p></p><p class=MsoNormal> 6<o:p></o:p></p><p class=MsoNormal> 0.000 0.0 0.0 40<o:p></o:p></p><p class=MsoNormal> 1.0 0.0 0.0 20<o:p></o:p></p><p class=MsoNormal> 1.0 0.5 0.0 20<o:p></o:p></p><p class=MsoNormal> 1.0 1.0 0.0 40<o:p></o:p></p><p class=MsoNormal> 0.00 0.0 0.0 40<o:p></o:p></p><p class=MsoNormal> 0.5 0.5 0.5 1<o:p></o:p></p><p class=MsoNormal>###################################################<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>However, if I used the first dynamic matrix, which corresponds to that of Gamma point, and do phonon calculation for Gamma point via dynmat.x with the input file:<o:p></o:p></p><p class=MsoNormal>###################################################<o:p></o:p></p><p class=MsoNormal>&input<o:p></o:p></p><p class=MsoNormal> fildyn='AlAs.dynG', asr='simple',<o:p></o:p></p><p class=MsoNormal> fileig='AlAs.eigfile', filout='AlAs.normal'<o:p></o:p></p><p class=MsoNormal> amass(1)=26.98,<o:p></o:p></p><p class=MsoNormal> amass(2)=74.92,<o:p></o:p></p><p class=MsoNormal> /<o:p></o:p></p><p class=MsoNormal>###################################################<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>The output file is <o:p></o:p></p><p class=MsoNormal>###################################################<o:p></o:p></p><p class=MsoNormal> Program DYNMAT v.6.8 starts on 3Feb2022 at 16:54:12<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal> This program is part of the open-source Quantum ESPRESSO suite<o:p></o:p></p><p class=MsoNormal> for quantum simulation of materials; please cite<o:p></o:p></p><p class=MsoNormal> "P. Giannozzi et al., J. Phys.:Condens. Matter 21 395502 (2009);<o:p></o:p></p><p class=MsoNormal> "P. Giannozzi et al., J. Phys.:Condens. Matter 29 465901 (2017);<o:p></o:p></p><p class=MsoNormal> "P. Giannozzi et al., J. Chem. Phys. 152 154105 (2020);<o:p></o:p></p><p class=MsoNormal> URL <a href="http://www.quantum-espresso.org">http://www.quantum-espresso.org</a>",<o:p></o:p></p><p class=MsoNormal> in publications or presentations arising from this work. More details at<o:p></o:p></p><p class=MsoNormal> <a href="http://www.quantum-espresso.org/quote">http://www.quantum-espresso.org/quote</a><o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal> Parallel version (MPI & OpenMP), running on 4 processor cores<o:p></o:p></p><p class=MsoNormal> Number of MPI processes: 4<o:p></o:p></p><p class=MsoNormal> Threads/MPI process: 1<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal> MPI processes distributed on 1 nodes<o:p></o:p></p><p class=MsoNormal> R & G space division: proc/nbgrp/npool/nimage = 4<o:p></o:p></p><p class=MsoNormal> 121728 MiB available memory on the printing compute node when the environment starts<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal> Reading Dynamical Matrix from file AlAs.dynG<o:p></o:p></p><p class=MsoNormal> ...Force constants read<o:p></o:p></p><p class=MsoNormal> ...epsilon and Z* read<o:p></o:p></p><p class=MsoNormal> A direction for q was not specified:TO-LO splitting will be absent<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal> Polarizability (A^3 units)<o:p></o:p></p><p class=MsoNormal> multiply by 0.190550 for Clausius-Mossotti correction<o:p></o:p></p><p class=MsoNormal> 43.491364 0.000000 -0.000000<o:p></o:p></p><p class=MsoNormal> 0.000000 43.491364 0.000000<o:p></o:p></p><p class=MsoNormal> -0.000000 0.000000 43.491364<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal> IR activities are in (D/A)^2/amu units<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal># mode [cm-1] [THz] IR<o:p></o:p></p><p class=MsoNormal> 1 -0.00 -0.0000 0.0000<o:p></o:p></p><p class=MsoNormal> 2 0.00 0.0000 0.0000<o:p></o:p></p><p class=MsoNormal> 3 0.00 0.0000 0.0000<o:p></o:p></p><p class=MsoNormal> 4 375.52 11.2578 7.6122<o:p></o:p></p><p class=MsoNormal> 5 375.52 11.2578 7.6122<o:p></o:p></p><p class=MsoNormal> 6 375.52 11.2578 7.6122<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal> DYNMAT : 0.01s CPU 0.03s WALL<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal> This run was terminated on: 16:54:12 3Feb2022<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>=------------------------------------------------------------------------------=<o:p></o:p></p><p class=MsoNormal> JOB DONE.<o:p></o:p></p><p class=MsoNormal>=------------------------------------------------------------------------------=<o:p></o:p></p><p class=MsoNormal>###################################################<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>The highest phonon frequencies (375.52 cm $^{-1}$ vs 410.5615 cm $^{-1}$) differ a lot. <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Could anyone please tell me why there is such a hug difference? And how to eliminate the difference?<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>I also tried calculations on other materials, the single point phonon calculation at Gamma point gives no imaginary frequencies while when doing the full phonon dispersion calculations, there exist huge imaginary frequencies at Gamma point. <o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Best regards<o:p></o:p></p><p class=MsoNormal>Zhishuo Huang<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p></div></body></html>