Joseph,<br><br>I&#39;ve never tried computing quite that many bands (only gone up to about 1200 bands) or with that many atoms per cell. However, I&#39;ve never not been able to get the bands to converge using both &#39;cg&#39; diagonalization and also with increasing the cutoff. You mentioned that you have tried increasing the cutoff, but perhaps you haven&#39;t gone high enough? <br>
<br>Here is a quote from Paolo that might be relevant in your situation:<br><br><br><blockquote style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;" class="gmail_quote"><pre>the algorithm used in PWscf is based on the assumption that<br>
the number of Kohn-Sham states (aka bands) is much smaller<br>than the number of basis functions (i.e. plane waves). If this<br>assumption doesn&#39;t hold (and it doesn&#39;t if you calculate 500<br>bands in fcc Si unless you use a very large cutoff) you may<br>
run into trouble. Moreover iterative diagonalization becomes<br>slower and more memory-consuming than conventional<br>diagonalization. If you really need to do that, you need to<br>modify the code<br></pre></blockquote><div>
<br>Best,<br>Brad Malone<br>UC Berkeley <br></div>