<br>Dear Forum,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; I have been trying to calculate the energy of a CO crystal using the&nbsp; <br>&nbsp;coordinates given in cst-www.nrl.navy.mul/lattice/struk.xmol/b21.pos.&nbsp; <br>This has 14 CO molecules. I used the format shown in example03 with the proper<br>
&nbsp;changes for the number of atoms and crystal structure.&nbsp; Unfortunately the cal-<br>culation does not converge.&nbsp; Increasing the number of iterations and ecutwfc and<br>ecutrho improved things a little but the convergence is still about 10^-4 Ry. I am<br>
particularly interested in the interactions between molecules so this is not satis-<br>factory.&nbsp; Is there anything I can do to improve the calculation to reach convergence<br>to about 10^-7 Ry?<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Incidentally the calculated triple bond energy for a single CO molecule is about <br>
43 Ry while the literature value given in texts is about 54 Ry.&nbsp; Is this to be expected?<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Hopefully someone has ideas about how to improve the convergence.<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Amos Leffler<br>
<br>