Joseph,<br><br>I've never tried computing quite that many bands (only gone up to about 1200 bands) or with that many atoms per cell. However, I've never not been able to get the bands to converge using both 'cg' diagonalization and also with increasing the cutoff. You mentioned that you have tried increasing the cutoff, but perhaps you haven'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't hold (and it doesn'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>