Dear All,<br><br>I am trying to calculate a spin polarized state using USPP with GGA+U.<br>But the calculation was aborted with an message as the following:<br><br>&nbsp;&nbsp;&nbsp;&nbsp; iteration #&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp; ecut=&nbsp;&nbsp;&nbsp; 37.00 Ry&nbsp;&nbsp;&nbsp;&nbsp; beta=0.20<br>&nbsp;&nbsp;&nbsp;&nbsp; CG style diagonalization<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 4 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 3 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 5 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 2 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 1 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 4 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 4 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 1 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 5 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 3 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 3 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 4 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 1 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 4 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 4 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 3 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 4 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 5 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 5 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 3 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 5 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 5 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 2 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 3 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 4 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 3 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 3 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 4 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 3 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 5 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 1 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 5 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 1 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 3 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 5 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 5 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 5 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 5 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 4 eigenvalues not converged<br>&nbsp;&nbsp;&nbsp;&nbsp; c_bands:&nbsp; 5 eigenvalues not converged<br><br>&nbsp;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br>&nbsp;&nbsp;&nbsp;&nbsp; from davcio : error #&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 10<br>&nbsp;&nbsp;&nbsp;&nbsp; error while reading from file<br>&nbsp;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br><br>&nbsp;&nbsp;&nbsp;&nbsp; stopping ... <br><br><br>Please let me know how to avoid this problem.<br>
<br>
<br>Just in case, I also pasted a part of my input file (see below).<br>
<br>&amp;control<br>&nbsp;&nbsp;&nbsp; calculation='scf'<br>&nbsp;&nbsp;&nbsp; restart_mode='restart',<br>&nbsp;&nbsp;&nbsp; prefix='vso'<br>&nbsp;&nbsp;&nbsp; pseudo_dir = '$PSEUDO_DIR/',<br>&nbsp;&nbsp;&nbsp; outdir='$TMP_DIR/'<br>&nbsp;/<br>&nbsp;&amp;system<br>&nbsp;&nbsp;&nbsp; ibrav= 6, A= 11.2338, B= 11.2338, C= 7.8711, cosAB= 0, cosAC= 0, cosBC= 0, nat= 72, ntyp= 3,<br>&nbsp;&nbsp;&nbsp; ecutwfc =37.0, ecutrho= 296, nbnd= 342, occupations='smearing', smearing='gauss', degauss=0.01,<br>&nbsp;&nbsp;&nbsp; nspin= 2, starting_magnetization(1)= 0.5, report=1,<br>&nbsp;&nbsp;&nbsp; lda_plus_u=.true.&nbsp; Hubbard_U(1)=1<br>&nbsp;/<br>&nbsp;&amp;electrons<br>&nbsp;&nbsp;&nbsp; electron_maxstep= 300<br>&nbsp;&nbsp;&nbsp; mixing_mode = 'plain'<br>&nbsp;&nbsp;&nbsp; mixing_beta = 0.2<br>&nbsp;&nbsp;&nbsp; conv_thr =&nbsp; 1.0d-6<br>&nbsp;&nbsp;&nbsp; diagonalization= 'cg'<br>&nbsp;&nbsp;&nbsp; mixing_fixed_ns= 0<br>&nbsp;/<br><br><br>Best regards,<br>Sang-Hwan Kim<br><br><BR><BR>Department of Chemistry<br>University of Houston<br>136 Fleming Building<br>Houston, TX 77204-5003</BR></BR>