Hello All,<br>I am trying to calculate the phonon dispersion relations and phonon density of states of UO2.<br>Please find the attached files <a href="http://UO2.scf.in">UO2.scf.in</a> and <a href="http://UO2.ph.in">UO2.ph.in</a><br>
<br>The calculations assume GGA and include spin-polarized effects.<br clear="all"><br>There are 6 representations and nine modes. The problem I have is: When I run the input file <a href="http://UO2.ph.in">UO2.ph.in</a>, the second  representation of the second mode for q (-0.1666667   0.1666667  -0.1666667), the<span style="color: rgb(255, 0, 0);"> <b style="color: rgb(255, 0, 0);">scf </b></span><b style="color: rgb(255, 0, 0);">calculations do not converge.</b><br>
<br>Notes:<br>1- My runs assume UO2  is metallic, in fact UO2 is insulator with band gap of 2eV.( However, when I do the calculations using the VASP code I got very good result only by including the polarization effects) <br>
2- In order to account for the band gap, one should implement the Hubbard term (LDA/GGA+U)<br>3- I tried different occupations  (smearing, tetrahedra)<br>4- I tried different starting_magnetization= 0.5, 1.0<br>5- the lattice parameter a= 10.285, is based on optimizing the structure and including spin-polarization effects.<br>
<br><br>Any adivce to overcome this problem is highly welcomed.<br><br>Kindest Regards,<br>______________________________<br>aIYAD I. AL-QASIR, PhD<br>Research Associate <br><br>Department of Nuclear Engineering<br>North Carolina State University<br>
Campus Box 7909<br>2500 Stinson Dr.<br>Raleigh, NC 27695-7909<br>