<div dir="ltr">please consider it as a speculation:<br>I guess you are using GGA for your calculations. It may make sense to check whether such problem appears within LDA or not.<br><br>Javad Hashemifar<br><br><br><div class="gmail_quote">
On Tue, Jul 29, 2008 at 9:19 AM, Stefaan Cottenier <span dir="ltr">&lt;<a href="mailto:Stefaan.Cottenier@fys.kuleuven.be">Stefaan.Cottenier@fys.kuleuven.be</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
Hello again,<br>
<br>
I have been too optimistic, yesterday. The reported problem<br>
(&quot;eigenvalues below...&quot;) did indeed disappear when I modified the<br>
sphere radii and/or put deep semi-core states in the core, as reported<br>
yesterday. But when I applied the same modification to a Cu atom at a<br>
different site in Ge (a case that was working well with the old<br>
settings), exactly the same problem now appeared there...<br>
<br>
This happens, in the very first iteration, in a directory that was<br>
freshly initialized with not old files lingering around:<br>
<br>
an error message in lapw2_3.error: &#39;l2main&#39; - QTL-B.GT.15.,<br>
Ghostbands, check scf files<br>
<br>
end of case.scf2_3 :<br>
<br>
 &nbsp; &nbsp;QTL-B VALUE .EQ. ********** &nbsp; in Band of energy &nbsp; -5.71041 &nbsp; ATOM=<br>
 &nbsp; &nbsp;5 &nbsp; L= &nbsp;2<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>
<br>
But atom 5 (which is Ge) does not have states at that energy:<br>
<br>
 &nbsp; 0.30 &nbsp; &nbsp;4 &nbsp;0 &nbsp; &nbsp; &nbsp;(GLOBAL E-PARAMETER WITH n OTHER CHOICES, global APW/LAPW)<br>
 &nbsp;2 &nbsp; -1.82 &nbsp; &nbsp; &nbsp;0.010 CONT 1<br>
 &nbsp;2 &nbsp; &nbsp;0.30 &nbsp; &nbsp; &nbsp;0.000 CONT 1<br>
 &nbsp;0 &nbsp; &nbsp;0.30 &nbsp; &nbsp; &nbsp;0.000 CONT 1<br>
 &nbsp;1 &nbsp; &nbsp;0.30 &nbsp; &nbsp; &nbsp;0.000 CONT 1<br>
<br>
The only other element in the cell is Cu (ATOM 1), which has these<br>
linearization energies:<br>
<br>
 &nbsp; 0.30 &nbsp; &nbsp;4 &nbsp;0 &nbsp; &nbsp; &nbsp;(GLOBAL E-PARAMETER WITH n OTHER CHOICES, global APW/LAPW)<br>
 &nbsp;1 &nbsp; &nbsp;0.30 &nbsp; &nbsp; &nbsp;0.000 CONT 1<br>
 &nbsp;1 &nbsp; -5.36 &nbsp; &nbsp; &nbsp;0.005 STOP 1<br>
 &nbsp;2 &nbsp; &nbsp;0.30 &nbsp; &nbsp; &nbsp;0.010 CONT 1<br>
 &nbsp;0 &nbsp; &nbsp;0.30 &nbsp; &nbsp; &nbsp;0.000 CONT 1<br>
<br>
It looks like eigenstates for Cu and Ge are somehow interchanged by<br>
the error message...?<br>
<br>
The end of case.scf1_3 :<br>
<br>
:EIG00936: &nbsp; &nbsp; &nbsp; 1.9587652 &nbsp; &nbsp;1.9605022 &nbsp; &nbsp;1.9620300 &nbsp; &nbsp;1.9651681 &nbsp; &nbsp;1.9677904<br>
:EIG00941: &nbsp; &nbsp; &nbsp; 1.9709236 &nbsp; &nbsp;1.9759120 &nbsp; &nbsp;1.9783975 &nbsp; &nbsp;1.9822565 &nbsp; &nbsp;1.9843245<br>
:EIG00946: &nbsp; &nbsp; &nbsp; 1.9858109 &nbsp; &nbsp;1.9880890 &nbsp; &nbsp;1.9903135 &nbsp; &nbsp;1.9927309 &nbsp; &nbsp;1.9944311<br>
<br>
:EIG00951: &nbsp; &nbsp; &nbsp; 1.9963998 &nbsp; &nbsp;1.9968078 &nbsp; &nbsp;1.9992533<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 2 EIGENVALUES BELOW THE ENERGY &nbsp; -9.00000<br>
 &nbsp; &nbsp; &nbsp; &nbsp;********************************************************<br>
<br>
:KPT &nbsp; : &nbsp; &nbsp; &nbsp;NUMBER OF K-POINTS: &nbsp; &nbsp; 1<br>
<br>
All other case.scf1_* files have &quot;0 eigenvalues below&quot;. And these are<br>
the deepest eigenvalues in case.scf1:<br>
<br>
:RKM &nbsp;: MATRIX SIZE 5615LOs: 908 &nbsp;RKM= 6.53 &nbsp;WEIGHT= 8.00 &nbsp;PGR:<br>
 &nbsp; &nbsp; &nbsp; &nbsp;EIGENVALUES ARE:<br>
:EIG00001: &nbsp; &nbsp; &nbsp;-4.9618775 &nbsp; -4.9618775 &nbsp; -4.9618749 &nbsp; -1.4433134 &nbsp; -1.4425104<br>
:EIG00006: &nbsp; &nbsp; &nbsp;-1.4425104 &nbsp; -1.4424955 &nbsp; -1.4415851 &nbsp; -1.4415851 &nbsp; -1.4414788<br>
:EIG00011: &nbsp; &nbsp; &nbsp;-1.4414788 &nbsp; -1.4414291 &nbsp; -1.4414284 &nbsp; -1.4412709 &nbsp; -1.4412709<br>
:EIG00016: &nbsp; &nbsp; &nbsp;-1.4411669 &nbsp; -1.4411669 &nbsp; -1.4410760 &nbsp; -1.4410224 &nbsp; -1.4409746<br>
<br>
Here are the linearization energies that lapw1 determined (no in1new):<br>
<br>
For Cu:<br>
<br>
<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ATOMIC SPHERE DEPENDENT PARAMETERS FOR ATOM &nbsp;Cu1<br>
:e__0001: OVERALL ENERGY PARAMETER IS &nbsp; &nbsp;0.3000<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; OVERALL BASIS SET ON ATOM IS LAPW<br>
:E1_0001: E( 1)= &nbsp; &nbsp;0.3000<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;APW+lo<br>
:E1_0001: E( 1)= &nbsp; -4.9650 &nbsp; E(BOTTOM)= &nbsp; -4.995 &nbsp; E(TOP)= &nbsp; -4.935<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;LOCAL ORBITAL<br>
:E2_0001: E( 2)= &nbsp; -0.0450 &nbsp; E(BOTTOM)= &nbsp; -0.270 &nbsp; E(TOP)= &nbsp; &nbsp;0.180<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;APW+lo<br>
:E0_0001: E( 0)= &nbsp; &nbsp;0.3000<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;APW+lo<br>
<br>
For Ge:<br>
<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ATOMIC SPHERE DEPENDENT PARAMETERS FOR ATOM &nbsp;Ge3<br>
:e__0004: OVERALL ENERGY PARAMETER IS &nbsp; &nbsp;0.3000<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; OVERALL BASIS SET ON ATOM IS LAPW<br>
:E2_0004: E( 2)= &nbsp; -1.4450 &nbsp; E(BOTTOM)= &nbsp; -1.530 &nbsp; E(TOP)= &nbsp; -1.360<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;APW+lo<br>
:E2_0004: E( 2)= &nbsp; &nbsp;0.3000<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;LOCAL ORBITAL<br>
:E0_0004: E( 0)= &nbsp; &nbsp;0.3000<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;APW+lo<br>
:E1_0004: E( 1)= &nbsp; &nbsp;0.3000<br>
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;APW+lo<br>
<br>
I used RKMax=6.53 (Kmax=3.00), and 4 k-points -- values that should be<br>
more than good enough for reasonable results. Changing the initial<br>
linearization energies did not help to the values reported case.scf1<br>
did not help.<br>
<br>
The real reason of this behaviour must lie elsewhere...?<br>
<br>
Stefaan<br>
<br>
Disclaimer: <a href="http://www.kuleuven.be/cwis/email_disclaimer.htm" target="_blank">http://www.kuleuven.be/cwis/email_disclaimer.htm</a><br>
<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>
</blockquote></div><br><br clear="all"><br>-- <br>======================================<br>Seyed Javad Hashemifar, Ph.D.<br>current address:<br>Physics Department, University of Duisburg-Essen, <br>47057 Duisburg, Germany<br>
Tel:+49 203 379 4743 Fax:+49 203 379 4742<br>permanent address:<br>Physics Department, Isfahan University of Technology<br>84154 Isfahan, Iran<br>Tel: +98 311 391 2375 Fax:+98 311 3912376<br>---------------------------------------------------------------------------<br>

</div>