<div dir="ltr"><div dir="ltr">On Wed, Dec 23, 2020 at 1:46 PM mkondrin <<a href="mailto:mkondrin@hppi.troitsk.ru">mkondrin@hppi.troitsk.ru</a>> wrote:</div><div dir="ltr"><br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
&SYSTEM<br>
   ibrav = 0<br>
   A =   5.658<br></blockquote><div><br></div><div>A is the experimental lattice parameter and would be fine for ibrav=2 (fcc), but then you set your lattice vectors:<br></div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
CELL_PARAMETERS (alat= 10.69207046)<br>
  -0.510237966   0.000000000   0.510237966<br>
   -0.000000000   0.510237966   0.510237966<br>
   -0.510237966   0.510237966   0.000000000<br></blockquote><div><br></div><div>to a larger cell (they are in "alat" unit, that is, the value specified as "A" or "celldm(1)", NOT the value in the card). With ibrav=2 you get something like 290 cm^-1. Note that the phonons are typically computed at the theoretical equilibrium lattice parameter, that in LDA is typically smaller than the experimental, leading to higher frequencies.<br></div><div><br></div><div> Also note that the convergence threshold conv_thr=1.e-5 in the scf calculation is too large: reduce it to at least 1.0e-8 (it costs close to nothing with respect to the phonon calculation)<br></div><div><br></div><div>Paolo<br></div></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div></div>