<div style="line-height:1.7;color:#000000;font-size:14px;font-family:Arial"><div>Dear all,</div><div>I have calculated the phonon of GaN. Firstly, I have performed a relaxation  and attained a good optimized result as shown below.</div><div> </div><div>Forces acting on atoms (Ry/au):</div><div>     atom    1 type  2   force =     0.00000000    0.00000000   -0.00000539<br>     atom    2 type  2   force =     0.00000000    0.00000000   -0.00000197<br>     atom    3 type  1   force =     0.00000000   -0.00000000    0.00000498<br>     atom    4 type  1   force =     0.00000000   -0.00000000    0.00000237</div><div>     Total force =     0.000008     Total SCF correction =     0.000001</div><div><br>     entering subroutine stress ...</div><div>          total   stress  (Ry/bohr**3)                   (kbar)     P=    0.04<br>   0.00000001  -0.00000000   0.00000000          0.00     -0.00      0.00<br>  -0.00000000   0.00000001  -0.00000000         -0.00      0.00     -0.00<br>   0.00000000   0.00000000   0.00000071          0.00      0.00      0.10</div><div><br>     bfgs converged in   3 scf cycles and   2 bfgs steps<br>     (criteria: energy <  1.0E-06, force <  1.0E-05)</div><div>     End of BFGS Geometry Optimization<br></div><div>Then, I have performed a scf calculation. The stress convergence was also good as shown below.</div><div> </div><div>   entering subroutine stress ...</div><div>          total   stress  (Ry/bohr**3)                   (kbar)     P=    0.05<br>   0.00000009  -0.00000000   0.00000000          0.01     -0.00      0.00<br>  -0.00000000   0.00000009   0.00000000         -0.00      0.01      0.00<br>   0.00000000  -0.00000000   0.00000078          0.00     -0.00      0.11<br></div><div>Lastly, I have performed the phonon calculation at G point. The results are shown below.</div><div> </div><div> **************************************************************************<br>     freq (    1) =      -0.861081 [THz] =     -28.722558 [cm-1]<br>     freq (    2) =      -0.861081 [THz] =     -28.722558 [cm-1]<br>     freq (    3) =       0.290776 [THz] =       9.699259 [cm-1]<br>     freq (    4) =       4.004845 [THz] =     133.587259 [cm-1]<br>     freq (    5) =       4.004845 [THz] =     133.587259 [cm-1]<br>     freq (    6) =       9.811711 [THz] =     327.283458 [cm-1]<br>     freq (    7) =      15.564963 [THz] =     519.191266 [cm-1]<br>     freq (    8) =      16.231642 [THz] =     541.429307 [cm-1]<br>     freq (    9) =      16.231642 [THz] =     541.429307 [cm-1]<br>     freq (   10) =      16.510116 [THz] =     550.718195 [cm-1]<br>     freq (   11) =      16.510116 [THz] =     550.718195 [cm-1]<br>     freq (   12) =      20.314461 [THz] =     677.617466 [cm-1]<br> **************************************************************************</div><div>     Mode symmetry, C_3v (3m)   point group:</div><div>     freq (  1 -  2) =        -28.7  [cm-1]   --> E    L_3        I+R<br>     freq (  3 -  3) =          9.7  [cm-1]   --> A_1  L_1        I+R<br>     freq (  4 -  5) =        133.6  [cm-1]   --> E    L_3        I+R<br>     freq (  6 -  6) =        327.3  [cm-1]   --> A_1  L_1        I+R<br>     freq (  7 -  7) =        519.2  [cm-1]   --> A_1  L_1        I+R<br>     freq (  8 -  9) =        541.4  [cm-1]   --> E    L_3        I+R<br>     freq ( 10 - 11) =        550.7  [cm-1]   --> E    L_3        I+R<br>     freq ( 12 - 12) =        677.6  [cm-1]   --> A_1  L_1        I+R<br></div><div>However, I have some questions about these results.</div><div>1. The frequences of freq (    1) ~freq (    3)  should be zero, why there are not?</div><div>2. The hexagonal GaN is belong to C_6v (6mm)  point group, why our result is C_3v (3m)   point group as above shown?</div><div>3. For N and Ga elements, there are too many pesedopotential, how to chose the pesedopotentials?</div><div> </div><div>The input files and output files in my calculation are attached in attachments.</div><div> </div><div><div>Thank you in advance!</div><div><br></div><div>Best Rgards!</div><div><br></div><div>Yu Hailin</div><div><br></div><div>College of Physics & Electronic Enginerring,</div><div>Changshu Institute of Technology</div><div>China<br></div></div><div> </div><div> </div><div> </div><div> </div></div><br><br><span title="neteasefooter"><span id="netease_mail_footer"></span></span>