Hello Wien Users,<br><div class="gmail_quote"> &nbsp;Many thanks first to the forum for providing guidance to the wien package. I am interested in calculating the surface state of Bi, using a 1x1x10 supercell (10 bilayers) with a 5x5x1 k-mesh. &nbsp;After initializing the calculation and accepting the geometry offered by nn, I ran the scf with spin orbit coupling. The calculation ran for 21 iterations, until at the end of the 22nd the calculation stopped:
<br><i> cycle 22         (Fri Dec &nbsp;7 14:45:49 EST 2007)         (19/19 to go)<br><br>&gt; &nbsp; lapw0         (14:45:49) 228.126u 1.693s 3:56.95 96.9%        0+0k 0+0io 24pf+0w<br>&gt; &nbsp; lapw1 &nbsp;-c &nbsp;        (14:49:47) 9935.601u 146.192s 7:40:52.37 36.4%        0+0k 0+0io 640167pf+0w
<br>&gt; &nbsp; lapwso &nbsp;-c        (22:30:44) 2084.313u 31.731s 1:26:14.72 40.8%        0+0k 0+0io 117209pf+0w<br>&gt; &nbsp; lapw2 -c -so &nbsp;        (23:57:03) 1310.015u 166.919s 31:28.94 78.1%        0+0k 0+0io 4578pf+0w<br><br>&gt; &nbsp; stop error<br></i> <br>Tha 
lapw2.error reads:<i> &#39;l2main&#39; - QTL-B.GT.15., Ghostbands, check scf files</i><br><br>and the end of scf2 file reads: &nbsp; &nbsp;<i>QTL-B VALUE .EQ. ********** &nbsp; in Band of energy &nbsp; -4.74977 &nbsp; ATOM= &nbsp; 44 &nbsp; L= &nbsp;0<br> &nbsp; &nbsp;Check for ghostbands or EIGENVALUES BELOW XX messages
<br> &nbsp; &nbsp;Adjust your Energy-parameters or use -in1new switch, check RMTs &nbsp;!!!<br>:WARN : QTL-B value eq.*******&nbsp; in Band of energy&nbsp; -4.74977&nbsp;&nbsp; ATOM=&nbsp;&nbsp; 44&nbsp;&nbsp; L=&nbsp; 0<br>:WARN : You should change the E-parameter in case.in1 or use -in1new switch
<br><br></i>I&nbsp; have used the same RMT values and energy-parameters before to do bulk calculations with the Bi unit cell, and that calculation has ran well without any problems. My in1 file reads:<br><i>WFFIL&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (WFPRI, SUPWF)
<br>&nbsp; 7.00&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 10&nbsp;&nbsp;&nbsp; 4 (R-MT*K-MAX; MAX L IN WF, V-NMT<br>&nbsp; 0.30&nbsp;&nbsp;&nbsp; 5&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global APW/LAPW)<br>&nbsp;2&nbsp;&nbsp; -1.54&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.010 CONT...<br>0.000 CONT 1<br>K-VECTORS FROM UNIT:4&nbsp;&nbsp; -9.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
3.0&nbsp; 519&nbsp;&nbsp; emin/emax/nband #red<br></i>And my inso file reads:<br><i>WFFIL<br>&nbsp;4&nbsp; 1&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; llmax,ipr,kpot <br>&nbsp;-10.0000&nbsp;&nbsp; 1.50000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; emin,emax (output energy window)<br>&nbsp;&nbsp; 1.&nbsp; 1.&nbsp; 1.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; direction of magnetization (lattice vectors)
<br>&nbsp;0&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; number of atoms for which RLO is added<br>&nbsp;0&nbsp;&nbsp; -4.97&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.005&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; atom number,e-lo,de (case.in1), repeat NX times<br>&nbsp;0 0 0 0 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; number of atoms for which SO is switch off; atoms
<br><br></i>The SCF file reads, for the last iteration as:<br>&nbsp;&nbsp; ---------<br><i>:ITE001:&nbsp; 1. ITERATION<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>:NATO :&nbsp;&nbsp; 60 INDEPENDEND AND&nbsp;&nbsp; 60 TOTAL ATOMS IN UNITCELL<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; LATTICE&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = H&nbsp;&nbsp; 
<br>:POT&nbsp; : POTENTIAL OPTION&nbsp;&nbsp; 13<br>:LAT&nbsp; : LATTICE CONSTANTS=&nbsp; 8.56991&nbsp; 8.56991227.25225<br>:VOL&nbsp; : UNIT CELL VOLUME =&nbsp;&nbsp; 14454.10292<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; MODE OF CALCULATION IS&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = RELA<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NON-SPINPOLARIZED CALCULATION
<br>
:IFFT&nbsp; : FFT-parameters:&nbsp;&nbsp; 36&nbsp;&nbsp; 36&nbsp; 900 Factor: 1.00<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ATOMNUMBER=&nbsp; 1 Bi1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; VCOUL-ZERO =&nbsp; 0.21765E+00<br>:EFG001:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; EFG&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp;&nbsp;&nbsp; -0.00400&nbsp;&nbsp; *10**21&nbsp; V / m**2<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>

OVERALL ENERGY PARAMETER IS&nbsp;&nbsp;&nbsp; 0.3000<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; OVERALL BASIS SET ON ATOM IS LAPW<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; E( 2)=&nbsp;&nbsp; -1.2200&nbsp;&nbsp; E(BOTTOM)=&nbsp;&nbsp; -1.400&nbsp;&nbsp; E(TOP)=&nbsp;&nbsp; -1.040<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; APW+lo<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; E( 2)=&nbsp;&nbsp;&nbsp; 0.3000<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; LOCAL ORBITAL
<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; E( 0)=&nbsp;&nbsp; -0.7300&nbsp;&nbsp; E(BOTTOM)=&nbsp;&nbsp; -1.090&nbsp;&nbsp; E(TOP)= -200.000<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; APW+lo<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; E( 0)=&nbsp;&nbsp;&nbsp; 0.3000<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; LOCAL ORBITAL<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; E( 1)=&nbsp;&nbsp;&nbsp; 0.3000<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; APW+lo<br><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; K=&nbsp;&nbsp; 
0.00000&nbsp;&nbsp; 0.00000&nbsp;&nbsp; 0.00000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1<br>:RKM&nbsp; : MATRIX SIZE 6227LOs: 900&nbsp; RKM= 6.98&nbsp; WEIGHT= 1.00&nbsp; PGR:&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; EIGENVALUES ARE:<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -1.3853569&nbsp;&nbsp; -1.3853569&nbsp;&nbsp; -1.3817454&nbsp;&nbsp; -1.3817454&nbsp;&nbsp; -1.3808292<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -
1.3808292&nbsp;&nbsp; -1.3783348&nbsp;&nbsp; -1.3783348&nbsp;&nbsp; -1.3762158&nbsp;&nbsp; -1.3762158<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; -1.3741646&nbsp;&nbsp; -1.3741646&nbsp;&nbsp; -1.3729990&nbsp;&nbsp; -1.3707902&nbsp;&nbsp; ....<br><br></i>Could someone please suggest where the problem might be? Any guidance will be much appreciated.
<br><br>Thank you.<br><font color="#888888"><br>&nbsp;<br><i></i>
</font></div><br>