<div dir="ltr">Dear Paolo;<div><br></div><div>Thank you very much for your answer</div><div><br></div><div><br></div><div>However, it is difficult for me to understand your answer. Can you be more specific?</div><div><br></div><div>For example, singularity of acoustic phonon at gamma and  zero coefficient...</div><div><br></div><div><br></div><div><br></div><div>and to make my original question more specific, </div><div><br></div><div>following is my output of ph.x</div><div><br></div><div><br></div><div><br></div><div>at Gamma point (q=0)</div><div>-----------------------------------------------------------------------------------------------</div><div><br></div><div><div>Mode symmetry, D_3h (-62m) point group:</div><div><br></div><div><br></div><div>     Atomic displacements:</div><div>     There are   6 irreducible representations</div><div><br></div><div>     Representation     1      1 modes -A'_1  To be done</div><div><br></div><div>     Representation     2      2 modes -E'  To be done</div><div><br></div><div>     Representation     3      2 modes -E'  To be done</div><div><br></div><div>     Representation     4      1 modes -A''2  To be done</div><div><br></div><div>     Representation     5      1 modes -A''2  To be done</div><div><br></div><div>     Representation     6      2 modes -E''  To be done</div></div><div><br></div><div><br></div><div><br></div><div>...</div><div>...</div><div><br></div><div><br></div><div><div>     Diagonalizing the dynamical matrix</div><div><br></div><div>     q = (    0.000000000   0.000000000   0.000000000 )</div><div><br></div><div> **************************************************************************</div><div>     freq (    1) =      -0.266351 [THz] =      -8.884502 [cm-1]</div><div>     freq (    2) =       0.470018 [THz] =      15.678118 [cm-1]</div><div>     freq (    3) =       0.470018 [THz] =      15.678119 [cm-1]</div><div>     freq (    4) =       8.700163 [THz] =     290.206203 [cm-1]</div><div>     freq (    5) =       8.700164 [THz] =     290.206217 [cm-1]</div><div>     freq (    6) =      11.760562 [THz] =     392.290107 [cm-1]</div><div>     freq (    7) =      11.760562 [THz] =     392.290126 [cm-1]</div><div>     freq (    8) =      12.300285 [THz] =     410.293340 [cm-1]</div><div>     freq (    9) =      14.310518 [THz] =     477.347485 [cm-1]</div><div> **************************************************************************</div><div><br></div><div>     Mode symmetry, D_3h (-62m) point group:</div><div><br></div><div>     freq (  1 -  1) =         -8.9  [cm-1]   --> A''2            I</div><div>     freq (  2 -  3) =         15.7  [cm-1]   --> E'              I+R</div><div>     freq (  4 -  5) =        290.2  [cm-1]   --> E''             R</div><div>     freq (  6 -  7) =        392.3  [cm-1]   --> E'              I+R</div><div>     freq (  8 -  8) =        410.3  [cm-1]   --> A'_1            R</div><div>     freq (  9 -  9) =        477.3  [cm-1]   --> A''2            I</div></div><div><br></div><div><br></div><div>----------------------------------------------------------------------------------------------------</div><div><br></div><div><br></div><div>At the gamma point as you can see from above, two dimensional irrep, E', appears twice as "Representation 2" and "Representation 3".</div><div><br></div><div>and i get <span style="color:rgb(0,0,0)">elph.1.2.xml , </span><span style="color:rgb(0,0,0)">elph.1.3.xml for them, which may be printed in order according to elphon.f90 code ( i am interested in extracting electron-phonon matrix element).</span></div><div><span style="color:rgb(0,0,0)"><br></span></div><div><span style="color:rgb(0,0,0)"><br></span></div><div><span style="color:rgb(0,0,0)">but the two irrep E'  are definitely inequivalent (different freq.) and i do not know whether the frequency of </span>"Representation 2" is "15.7" or "392.3".</div><div><br></div><div><br></div><div>This problem makes me hard to identify  <span style="color:rgb(0,0,0)">elph.1.2.xml and </span><span style="color:rgb(0,0,0)">elph.1.3.xml since i cannot assign the frequency.</span></div><div><span style="color:rgb(0,0,0)"><br></span></div><div><span style="color:rgb(0,0,0)"><br></span></div><div><span style="color:rgb(0,0,0)"><br></span></div><div><span style="color:rgb(0,0,0)">This problem becomes even worse when the site symmetry of q-point is reduced</span></div><div><span style="color:rgb(0,0,0)"><br></span></div><div><span style="color:rgb(0,0,0)"><br></span></div><div><span style="color:rgb(0,0,0)"><br></span></div><div><font color="#000000"> q = (    0.000000000   0.128300065   0.000000000 )</font><br></div><div><span style="color:rgb(0,0,0)">-------------------------------------------------------------------------------------------</span></div><div><span style="color:rgb(0,0,0)"><br></span></div><div><div style><font color="#000000">     Mode symmetry, C_2v (mm2)  point group:</font></div><div style><font color="#000000"><br></font></div><div style><font color="#000000"><br></font></div><div style><font color="#000000">     Atomic displacements:</font></div><div style><font color="#000000">     There are   9 irreducible representations</font></div><div style><font color="#000000"><br></font></div><div style><font color="#000000">     Representation     1      1 modes -A_1  D_1  S_1  To be done</font></div><div style><font color="#000000"><br></font></div><div style><font color="#000000">     Representation     2      1 modes -A_1  D_1  S_1  To be done</font></div><div style><font color="#000000"><br></font></div><div style><font color="#000000">     Representation     3      1 modes -A_1  D_1  S_1  To be done</font></div><div style><font color="#000000"><br></font></div><div style><font color="#000000">     Representation     4      1 modes -A_2  D_2  S_2  To be done</font></div><div style><font color="#000000"><br></font></div><div style><font color="#000000">     Representation     5      1 modes -B_1  D_3  S_3  To be done</font></div><div style><font color="#000000"><br></font></div><div style><font color="#000000">     Representation     6      1 modes -B_1  D_3  S_3  To be done</font></div><div style><font color="#000000"><br></font></div><div style><font color="#000000">     Representation     7      1 modes -B_2  D_4  S_4  To be done</font></div><div style><font color="#000000"><br></font></div><div style><font color="#000000">     Representation     8      1 modes -B_2  D_4  S_4  To be done</font></div><div style><font color="#000000"><br></font></div><div style><font color="#000000">     Representation     9      1 modes -B_2  D_4  S_4  To be done</font></div><div style="color:rgb(0,0,0)"><br></div></div><div style="color:rgb(0,0,0)"><br></div><div style="color:rgb(0,0,0)">...</div><div style="color:rgb(0,0,0)">...</div><div style="color:rgb(0,0,0)">...</div><div style="color:rgb(0,0,0)"><br></div><div style><div style><font color="#000000">     Mode symmetry, C_2v (mm2)  point group:</font></div><div style><font color="#000000"><br></font></div><div style><font color="#000000">     freq (  1 -  1) =         11.2  [cm-1]   --> B_2  D_4  S_4</font></div><div style><font color="#000000">     freq (  2 -  2) =         56.2  [cm-1]   --> B_1  D_3  S_3</font></div><div style><font color="#000000">     freq (  3 -  3) =         90.1  [cm-1]   --> A_1  D_1  S_1</font></div><div style><font color="#000000">     freq (  4 -  4) =        292.5  [cm-1]   --> A_2  D_2  S_2</font></div><div style><font color="#000000">     freq (  5 -  5) =        300.9  [cm-1]   --> B_2  D_4  S_4</font></div><div style><font color="#000000">     freq (  6 -  6) =        385.3  [cm-1]   --> A_1  D_1  S_1</font></div><div style><font color="#000000">     freq (  7 -  7) =        390.9  [cm-1]   --> B_1  D_3  S_3</font></div><div style><font color="#000000">     freq (  8 -  8) =        414.1  [cm-1]   --> A_1  D_1  S_1</font></div><div style><font color="#000000">     freq (  9 -  9) =        466.7  [cm-1]   --> B_2  D_4  S_4</font></div><div style><font color="#000000"><br></font></div><div style="color:rgb(0,0,0)">--------------------------------------------------------------------------------------------</div></div><div style="color:rgb(0,0,0)"><br></div><div style="color:rgb(0,0,0)"><br></div><div style="color:rgb(0,0,0)">As you can see, irrep A_1 appears triply. </div><div style="color:rgb(0,0,0)"><br></div><div style="color:rgb(0,0,0)"><br></div><div style="color:rgb(0,0,0)">How can i resolve this problem ?</div><div style="color:rgb(0,0,0)"><br></div><div style="color:rgb(0,0,0)"><br></div><div style="color:rgb(0,0,0)">(I found similar thread in pw_forum </div><div style="color:rgb(0,0,0)"><br></div><div style><font color="#000000"><a href="http://qe-forge.org/pipermail/pw_forum/2009-March/086475.html">http://qe-forge.org/pipermail/pw_forum/2009-March/086475.html</a></font><br></div><div style><font color="#000000"><br></font></div><div style><font color="#000000">but i couldn't find solution.)</font></div><div style="color:rgb(0,0,0)"><br></div><div style="color:rgb(0,0,0)"><br></div><div style="color:rgb(0,0,0)">Thank you in advance.</div><div style="color:rgb(0,0,0)"><br></div><div style="color:rgb(0,0,0)"><br></div><div style="color:rgb(0,0,0)"><br></div><div style="color:rgb(0,0,0)">Best regards</div><div style="color:rgb(0,0,0)"><br></div><div style="color:rgb(0,0,0)">Jaehong Park</div><div> </div><div><br></div><div><br></div><div><br></div></div>