Hello,<br>&nbsp; Thank you Professor Blaha and Faccio for your messages. About the crystal structure, I am using a supercell comprised of 10 unit cells in the z direction, with 6 atoms per unit cell, so there are 60 atoms. In otherewords, it&#39;s 60 layers instead of 10. The first few entries generated by supercell is given by:
<br><br><i>H&nbsp;&nbsp; LATTICE,NONEQUIV. ATOMS 10&nbsp;&nbsp;&nbsp;&nbsp; <br>MODE OF CALC=RELA unit=bohr<br>&nbsp; 8.569908&nbsp; 8.569908227.252250 90.000000 90.000000120.000000<br>ATOM&nbsp;&nbsp; 1: X=0.00000000 Y=0.00000000 Z=0.07525172<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; MULT= 6&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ISPLIT=15
<br>ATOM&nbsp;&nbsp; 1: X=0.00000000 Y=0.00000000 Z=0.02298812<br>ATOM&nbsp;&nbsp; 1: X=0.66666667 Y=0.33333333 Z=0.00975849<br>ATOM&nbsp;&nbsp; 1: X=0.66666667 Y=0.33333333 Z=0.05573474<br>ATOM&nbsp;&nbsp; 1: X=0.33333333 Y=0.66666667 Z=0.04250511<br>ATOM&nbsp;&nbsp; 1: X=
0.33333333 Y=0.66666667 Z=0.08848135<br>Bi&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NPT=&nbsp; 781&nbsp; R0=0.00000500 RMT=&nbsp;&nbsp;&nbsp; 2.0000&nbsp;&nbsp; Z: 83.0<br>LOCAL ROT MATRIX:&nbsp;&nbsp;&nbsp; 1.0000000 0.0000000 0.0000000<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0000000 1.0000000 0.0000000<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
0.0000000 0.0000000 1.0000000<br><br></i><div style="text-align: left;">After running symmetry, I adopted the struct generated by the file, given by:<br></div><i>H&nbsp;&nbsp; LATTICE,NONEQUIV.ATOMS: 60 156 P3m1&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; 
<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; RELA&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; <br>&nbsp; 8.569908&nbsp; 8.569908227.252250 90.000000 90.000000120.000000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>ATOM&nbsp; -1: X=0.00000000 Y=0.00000000 Z=0.00000000<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; MULT= 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ISPLIT= 4<br>Bi1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NPT=&nbsp; 781&nbsp; R0=0.00000500 RMT=&nbsp;&nbsp;&nbsp; 2.5000&nbsp;&nbsp; Z: 83.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>LOCAL ROT MATRIX:&nbsp;&nbsp;&nbsp; 1.0000000 0.0000000 0.0000000<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0000000 1.0000000 0.0000000
<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0000000 0.0000000 1.0000000<br>ATOM&nbsp; -2: X=0.00000000 Y=0.00000000 Z=0.94773640<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; MULT= 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ISPLIT= 4<br>Bi2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NPT=&nbsp; 781&nbsp; R0=0.00000500 RMT=&nbsp;&nbsp;&nbsp; 2.5000&nbsp;&nbsp; Z: 83.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
<br>LOCAL ROT MATRIX:&nbsp;&nbsp;&nbsp; 1.0000000 0.0000000 0.0000000<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0000000 1.0000000 0.0000000<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0000000 0.0000000 1.0000000<br></i>So the number of inequivalent atoms has been increased by symmetry. Should I not accept the struct file generated by it? Also, as a side question along the same vein, for the surface bands is 1x1 unit cell in the x-y direction sufficient? Ans is the 5x5x1 k-mesh sufficient for that purpose? I used 2x2x10 before but I am trying to minimize the number of unit cells used to save calculation time.
<br><br>The calculation doesn&#39;t seem to be converging,<br>DIS&nbsp; :&nbsp; CHARGE DISTANCE&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ( 0.1641572 for atom&nbsp;&nbsp; 37 spin 1)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.1545264<br>DIS&nbsp; :&nbsp; CHARGE DISTANCE&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ( 0.1913422 for atom&nbsp;&nbsp;&nbsp; 3 spin 1)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.1529912
<br>...<br>DIS&nbsp; :&nbsp; CHARGE DISTANCE&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ( 1.1335284 for atom&nbsp;&nbsp; 51 spin 1)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.5305745<br>DIS&nbsp; :&nbsp; CHARGE DISTANCE&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ( 1.1275002 for atom&nbsp;&nbsp; 18 spin 1)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.6506100<br>...<br>DIS&nbsp; :&nbsp; CHARGE DISTANCE&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ( 0.3828902
 for atom&nbsp;&nbsp; 49 spin 1)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.1680688<br>DIS&nbsp; :&nbsp; CHARGE DISTANCE&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ( 0.4418754 for atom&nbsp;&nbsp; 46 spin 1)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.2222958<br>DIS&nbsp; :&nbsp; CHARGE DISTANCE&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ( 0.4079805 for atom&nbsp;&nbsp; 47 spin 1)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.2201738<br><br>I will try to run a calculation without -scf. I did not do so before because each calculation takes at least 2 weeks, and the system is unphysical without spin-orbit coupling. Could you please explain why the calculation is not converging?
<br><br>Thank you.<br>Sincerely yours,<br>Matthew Xia<br><br><div class="gmail_quote">On Dec 10, 2007 2:03 AM, Peter Blaha &lt;<a href="mailto:pblaha@theochem.tuwien.ac.at">pblaha@theochem.tuwien.ac.at</a>&gt; wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
First of all, please check your convergence (grep :DIS case.scf)<br>Is it converging or diverging ?<br><br>Second: I&#39;m not sure if it matters, but remove the line<br><div class="Ih2E3d">&nbsp;&gt; &nbsp;0 &nbsp; -4.97 &nbsp; &nbsp; &nbsp;0.005 &nbsp; &nbsp; &nbsp; atom number,e-lo,de (
case.in1), repeat NX<br></div>from your inso file (You specified correctly &quot;0&quot; RLOs, so this line<br>should not be there).<br><br>Did you run without -so before ? Usually this should be the first step.<br><br>Your error ocured for atom 44 ? Why do you have all 60 atoms
<br>non-equivalent ? I&#39;d expect some symmetry ?<br>Ad I do not understand either why a 10 layer cell has 60 atoms ?<br><br><br>Matthew Xia schrieb:<br><div><div></div><div class="Wj3C7c">&gt; Hello Wien Users,<br>&gt; &nbsp;Many thanks first to the forum for providing guidance to the wien
<br>&gt; package. I am interested in calculating the surface state of Bi, using a<br>&gt; 1x1x10 supercell (10 bilayers) with a 5x5x1 k-mesh. &nbsp;After initializing<br>&gt; the calculation and accepting the geometry offered by nn, I ran the scf
<br>&gt; with spin orbit coupling. The calculation ran for 21 iterations, until<br>&gt; at the end of the 22nd the calculation stopped:<br>&gt; / cycle 22 (Fri Dec &nbsp;7 14:45:49 EST 2007) (19/19 to go)<br>&gt;<br>&gt; &nbsp;&gt; &nbsp; lapw0 (14:45:49) 
228.126u 1.693s 3:56.95 96.9% 0+0k 0+0io 24pf+0w<br>&gt; &nbsp;&gt; &nbsp; lapw1 &nbsp;-c &nbsp; (14:49:47) 9935.601u 146.192s 7:40:52.37 36.4% 0+0k<br>&gt; 0+0io 640167pf+0w<br>&gt; &nbsp;&gt; &nbsp; lapwso &nbsp;-c (22:30:44) 2084.313u 31.731s 1:26:14.72
 40.8% 0+0k 0+0io<br>&gt; 117209pf+0w<br>&gt; &nbsp;&gt; &nbsp; lapw2 -c -so &nbsp; (23:57:03) 1310.015u 166.919s 31:28.94 78.1% 0+0k<br>&gt; 0+0io 4578pf+0w<br>&gt;<br>&gt; &nbsp;&gt; &nbsp; stop error<br>&gt; /<br>&gt; Tha lapw2.error reads:/ &#39;l2main&#39; - 
QTL-B.GT.15., Ghostbands, check scf<br>&gt; files/<br>&gt;<br>&gt; and the end of scf2 file reads: &nbsp; &nbsp;/QTL-B VALUE .EQ. ********** &nbsp; in<br>&gt; Band of energy &nbsp; -4.74977 &nbsp; ATOM= &nbsp; 44 &nbsp; L= &nbsp;0<br>&gt; &nbsp; &nbsp;Check for ghostbands or EIGENVALUES BELOW XX messages
<br>&gt; &nbsp; &nbsp;Adjust your Energy-parameters or use -in1new switch, check RMTs &nbsp;!!!<br>&gt; :WARN : QTL-B value eq.******* &nbsp;in Band of energy &nbsp;-4.74977 &nbsp; ATOM=<br>&gt; 44 &nbsp; L= &nbsp;0<br>&gt; :WARN : You should change the E-parameter in 
case.in1 or use -in1new switch<br>&gt;<br>&gt; /I &nbsp;have used the same RMT values and energy-parameters before to do<br>&gt; bulk calculations with the Bi unit cell, and that calculation has ran<br>&gt; well without any problems. My in1 file reads:
<br>&gt; /WFFIL &nbsp; &nbsp; &nbsp; &nbsp;(WFPRI, SUPWF)<br>&gt; &nbsp; 7.00 &nbsp; &nbsp; &nbsp; 10 &nbsp; &nbsp;4 (R-MT*K-MAX; MAX L IN WF, V-NMT<br>&gt; &nbsp; 0.30 &nbsp; &nbsp;5 &nbsp;0 &nbsp; &nbsp; &nbsp;(GLOBAL E-PARAMETER WITH n OTHER CHOICES, global<br>&gt; APW/LAPW)<br>&gt; &nbsp;2 &nbsp; -1.54 &nbsp; &nbsp; &nbsp;0.010
 CONT...<br>&gt; 0.000 CONT 1<br>&gt; K-VECTORS FROM UNIT:4 &nbsp; -9.0 &nbsp; &nbsp; &nbsp; 3.0 &nbsp;519 &nbsp; emin/emax/nband #red<br>&gt; /And my inso file reads:<br>&gt; /WFFIL<br>&gt; &nbsp;4 &nbsp;1 &nbsp;0 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;llmax,ipr,kpot<br>&gt; &nbsp;-10.0000
 &nbsp; 1.50000 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; emin,emax (output energy window)<br>&gt; &nbsp; &nbsp;1. &nbsp;1. &nbsp;1. &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; direction of magnetization (lattice vectors)<br>&gt; &nbsp;0 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; number of atoms for which RLO is added<br>&gt; &nbsp;0 &nbsp; -
4.97 &nbsp; &nbsp; &nbsp;0.005 &nbsp; &nbsp; &nbsp; atom number,e-lo,de (case.in1), repeat NX times<br>&gt; &nbsp;0 0 0 0 0 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;number of atoms for which SO is switch<br>&gt; off; atoms<br>&gt;<br>&gt; /The SCF file reads, for the last iteration as:
<br>&gt; &nbsp; &nbsp;---------<br>&gt; /:ITE001: &nbsp;1. ITERATION<br>&gt;<br>&gt; :NATO : &nbsp; 60 INDEPENDEND AND &nbsp; 60 TOTAL ATOMS IN UNITCELL<br>&gt; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; LATTICE &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= H<br>&gt; :POT &nbsp;: POTENTIAL OPTION &nbsp; 13<br>
&gt; :LAT &nbsp;: LATTICE CONSTANTS= &nbsp;8.56991 &nbsp;8.56991227.25225<br>&gt; :VOL &nbsp;: UNIT CELL VOLUME = &nbsp; 14454.10292<br>&gt; &nbsp; &nbsp; &nbsp; &nbsp;MODE OF CALCULATION IS &nbsp; &nbsp; &nbsp; = RELA<br>&gt; &nbsp; &nbsp; &nbsp; &nbsp;NON-SPINPOLARIZED CALCULATION<br>&gt; :IFFT &nbsp;: FFT-parameters: &nbsp; 36 &nbsp; 36 &nbsp;900 Factor: 
1.00<br>&gt; &nbsp; &nbsp; &nbsp; &nbsp;ATOMNUMBER= &nbsp;1 Bi1 &nbsp; &nbsp; &nbsp; &nbsp;VCOUL-ZERO = &nbsp;0.21765E+00<br>&gt; :EFG001: &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;EFG &nbsp; &nbsp; &nbsp; &nbsp; = &nbsp; &nbsp;-0.00400 &nbsp; *10**21 &nbsp;V /<br>&gt; m**2<br>&gt;<br>&gt; OVERALL ENERGY PARAMETER IS &nbsp; &nbsp;0.3000<br>
&gt; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; OVERALL BASIS SET ON ATOM IS LAPW<br>&gt; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; E( 2)= &nbsp; -1.2200 &nbsp; E(BOTTOM)= &nbsp; -1.400 &nbsp; E(TOP)= &nbsp; -1.040<br>&gt; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;APW+lo<br>&gt; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; E( 2)= &nbsp; &nbsp;0.3000<br>&gt; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;LOCAL ORBITAL
<br>&gt; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; E( 0)= &nbsp; -0.7300 &nbsp; E(BOTTOM)= &nbsp; -1.090 &nbsp; E(TOP)= -200.000<br>&gt; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;APW+lo<br>&gt; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; E( 0)= &nbsp; &nbsp;0.3000<br>&gt; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;LOCAL ORBITAL<br>&gt; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; E( 1)= &nbsp; &nbsp;0.3000<br>&gt; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;APW+lo
<br>&gt;<br>&gt; &nbsp; &nbsp; &nbsp; &nbsp;K= &nbsp; 0.00000 &nbsp; 0.00000 &nbsp; 0.00000 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;1<br>&gt; :RKM &nbsp;: MATRIX SIZE 6227LOs: 900 &nbsp;RKM= 6.98 &nbsp;WEIGHT= 1.00 &nbsp;PGR:<br>&gt; &nbsp; &nbsp; &nbsp; &nbsp;EIGENVALUES ARE:<br>&gt; &nbsp; &nbsp; &nbsp; &nbsp; -1.3853569 &nbsp; -1.3853569 &nbsp; -1.3817454
 &nbsp; -1.3817454 &nbsp; -1.3808292<br>&gt; &nbsp; &nbsp; &nbsp; &nbsp; - 1.3808292 &nbsp; -1.3783348 &nbsp; -1.3783348 &nbsp; -1.3762158 &nbsp; -1.3762158<br>&gt; &nbsp; &nbsp; &nbsp; &nbsp; -1.3741646 &nbsp; -1.3741646 &nbsp; -1.3729990 &nbsp; -1.3707902 &nbsp; ....<br>&gt;<br>&gt; /Could someone please suggest where the problem might be? Any guidance
<br>&gt; will be much appreciated.<br>&gt;<br>&gt; Thank you.<br>&gt; MX<br>&gt;<br>&gt;<br></div></div>&gt; ------------------------------------------------------------------------<br><div><div></div><div class="Wj3C7c">
&gt;<br>&gt; _______________________________________________<br>&gt; Wien mailing list<br>&gt; <a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br>&gt; <a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>_______________________________________________<br>Wien mailing list<br><a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br></div></div></blockquote></div><br>