
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">

<HTML><HEAD>

<META http-equiv=Content-Type content="text/html; charset=iso-8859-1">

<META content="MSHTML 6.00.2900.2873" name=GENERATOR>

<STYLE></STYLE>

</HEAD>

<BODY bgColor=#ffffff>

<DIV><FONT face=Arial size=2>Dear M. Xia&nbsp; <BR>&nbsp;&nbsp;&nbsp; I don't 

have much experience in slab calculations, but in the case of bulk calculations 

I observe similar problems regarding this type of error (ghost bands), just 

because of the optimization of the atomic coordinates.&nbsp;&nbsp; 

<BR>&nbsp;&nbsp;&nbsp; This is why I guess that your problem is a common one, 

since you are moving from a bulk to a slab calculation, and you will have 

changes in the values of the linearization energies.&nbsp; 

<BR>&nbsp;&nbsp;&nbsp; According to the error descriptions in the file 

lapw2.error, you must choose a better energy parameters in the case.in1(c) file. 

Or, If you want to do that in a automatically, you should add the "-in1new x" 

option to your normal executing procedure e.g. "runsp -cc 0.0001 -so&nbsp; 

-in1new 3 ...."&nbsp;&nbsp; <BR>&nbsp;&nbsp; I repeat, I dont' have much 

experience in slabs calculations, but this is what I can tell you.&nbsp; 

<BR>Best regards&nbsp; <BR>Ricardo</FONT></DIV>

<DIV><FONT face=Arial size=2></FONT>&nbsp;</DIV>

<DIV>-------------------------------------------------------------------------<BR>-----&nbsp;&nbsp; 

Dr. Ricardo Faccio<BR>&nbsp;<BR>&nbsp; Mail: Cryssmat-Lab., Cátedra de Física, 

DETEMA<BR>&nbsp; Facultad de Química, Universidad de la 

República<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Av. Gral. Flores 2124, C.C. 

1157<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; C.P. 11800, Montevideo, 

Uruguay.<BR>&nbsp; E-mail: <A 

href="mailto:rfaccio@fq.edu.uy">rfaccio@fq.edu.uy</A><BR>&nbsp; Phone: 598 2 

9241860 Int. 

109<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

598 2 9290705<BR>&nbsp; Fax:&nbsp;&nbsp;&nbsp; 598 2 9241906<BR>&nbsp; 

Web:&nbsp; <A 

href="http://cryssmat.fq.edu.uy/ricardo/ricardo.htm">http://cryssmat.fq.edu.uy/ricardo/ricardo.htm</A></DIV>

<BLOCKQUOTE 

style="PADDING-RIGHT: 0px; PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: #000000 2px solid; MARGIN-RIGHT: 0px">

  <DIV style="FONT: 10pt arial">----- Original Message ----- </DIV>

  <DIV 

  style="BACKGROUND: #e4e4e4; FONT: 10pt arial; font-color: black"><B>From:</B> 

  <A title=xiay@princeton.edu href="mailto:xiay@princeton.edu">Matthew Xia</A> 

  </DIV>

  <DIV style="FONT: 10pt arial"><B>To:</B> <A 

  title=wien@zeus.theochem.tuwien.ac.at 

  href="mailto:wien@zeus.theochem.tuwien.ac.at">A Mailing list for WIEN2k 

  users</A> </DIV>

  <DIV style="FONT: 10pt arial"><B>Sent:</B> Sunday, December 09, 2007 6:01 

  PM</DIV>

  <DIV style="FONT: 10pt arial"><B>Subject:</B> [Wien] Error Calculating Bi 

  Surface Bands</DIV>

  <DIV><BR></DIV>Hello Wien Users,<BR>&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> 'l2main' - 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>MX</FONT> 

  <P>

  <HR>


  <P></P>_______________________________________________<BR>Wien mailing 

  list<BR>Wien@zeus.theochem.tuwien.ac.at<BR>http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien<BR>

  <P>

  <HR>


  <P></P>No virus found in this incoming message.<BR>Checked by AVG Free 

  Edition. <BR>Version: 7.5.503 / Virus Database: 269.16.17/1179 - Release Date: 

  09/12/2007 11:06<BR></BLOCKQUOTE></BODY></HTML>