Dear Jose,<br> I have tried all those also. when you give that then you will observe error in nn distances. Its not possible.<br><div class="gmail_quote">2012/2/17 Jose Alfredo Camargo Martinez <span dir="ltr"><<a href="mailto:sork8686@yahoo.es">sork8686@yahoo.es</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div style="font-size:10pt;font-family:arial,helvetica,sans-serif"><div><span lang="en"><span>One</span> <span>solution would be</span> <span>to increase</span> <span>the radii of</span> <span>MT</span></span></div>
<div> </div><div><span style="font-style:italic">José A. Camargo Martínez</span><br>Estudiante Doctorado en Ciencias - Física<br><span style="font-weight:bold">DEPARTAMENTO DE FÍSICA</span><br style="font-weight:bold"><span style="font-weight:bold">CINVESTAV - MÉXICO D.F.<br>
</span><br></div> <div style="font-family:arial,helvetica,sans-serif;font-size:10pt"> <div style="font-family:times new roman,new york,times,serif;font-size:12pt"> <div dir="ltr"> <font face="Arial"> <hr size="1"> <b><span style="font-weight:bold">De:</span></b> Madhav Ghimire
<<a href="mailto:ghimire.mpg@gmail.com" target="_blank">ghimire.mpg@gmail.com</a>><br> <b><span style="font-weight:bold">Para:</span></b> wien <<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a>> <br>
<b><span style="font-weight:bold">Enviado:</span></b> Jueves 16 de febrero de 2012 22:24<br> <b><span style="font-weight:bold">Asunto:</span></b> [Wien] Problems in generating case.struct using shell script and overlapping of spheres with error in lapw1<br>
</font> </div><div><div class="h5"> <br><div>Dear Peter Blaha and wien users,<br> I am trying to run a very simple compound Bi2Se3 with the given lattice parameters and atomic positions for space group R3m (166).<br>
The atomic positions along (x, y, z) are as follows:<br>
Bi = 0, 0. 0.40046<br>Se (1) = 0, 0, 0.2097<br>Se (2) = 0, 0, 0<br>With the above given parameters, I initiated the calculations, but when I select l start =-6 or -9 or any value between -1 to -9.9 for energy to separate core and valence states, it shows that core charge leaks out of the spheres. Then the calculations stuck there itself. The error to this is as listed below:<br>
Commandline: <b>x lstart </b><br>
Program input is: <b>"13
-6.0
"</b>
<br>
<br>
<pre> SELECT XCPOT:
recommended: 13: GGA (Perdew-Burke-Ernzerhof 96)
5: LSDA
11: GGA (Wu-Cohen 2006)
SELECT ENERGY to separate core and valence states:
recommended: -6.0 Ry (check how much core charge leaks out of MT-sphere)
WARNING: 0.328380476966373 CORE electrons leak out of MT-sphere !!!!
WARNING: 1.30946564796119 CORE electrons leak out of MT-sphere !!!!
WARNING: 1.30946564796119 CORE electrons leak out of MT-sphere !!!!
LSTART ENDS
0.537u 0.023s 0:00.59 93.2% 0+0k 8+0io 0pf+0w
</pre> When I perform for other compounds complexer than Bi2Se3, I resolve the problem easily by changing its RMT value or lattice parameters or the energy to separate core/valence states. Here I do not understand why this case arises. This case occurs especially for wien2k 05 to 09 version. I did the same calculation using wien2k 11 version, here I got through this problem by selecting 1.0 as the alternative option given as ALTERNATIVELY: specify charge localization (between 0.97 and 1.0) to select core state. Then once I complete initialization part and give a run scf cycles, I observe error in show dayfile as:<br>
<pre>start (Fri Feb 17 12:32:41 JST 2012) with lapw0 (40/99 to go)
cycle 1 (Fri Feb 17 12:32:41 JST 2012) (40/99 to go)
> lapw0 (12:32:41) 11.846u 0.116s 0:12.02 99.4% 0+0k 0+4712io 0pf+0w
> lapw1 -up (12:32:53) 150.227u 0.014s 2:30.80 99.6% 0+0k 0+1088io 0pf+0w
> stop error</pre>To clarify this, I checked the error file and error is observed in lapw1 as listed below:<br>Error in LAPW1<br> 'SELECT' - no energy limits found for L= 2 <br>
'SELECT' - E-bottom -10.46100 E-top -200.00000 <br>I do not understand why this case is observed in this particular compound although there had been a work on the same material using wien2k code. <br>
Apart from this, I have a question regarding <b>"How to generate a case.struct without using the w2web". I hope the experts in wien2k generally use shell script rather than the interface so hopefully you can instruct the step by step commands for generating case. struct.<br>
</b>Thank you in advance<br>Best regards<br>-- <br>M. P. Ghimire<br><br>
</div><br></div></div>_______________________________________________<br>Wien mailing list<br><a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">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>
<br><br> </div> </div> </div></div><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>
<br></blockquote></div><br><br clear="all"><br>-- <br>M. P. Ghimire<br><br>