Dear Gerhard,<br> You are right in the point that my RMT's are very small. But if you increase the RMT, you will observe error in x nn. I have used the value of Rmt for Bi as 1.75, Se 1 as 1.58 and Se (2) as 1.69 as referred in the paper of Mishra et al ; JPCM 9 (1997) 461 before using the data which I sent you all. So I could not clarify where the data was wrong. Hope you all could solve it out.<br>
Regards<br>Ghimire<br><br><div class="gmail_quote">On Fri, Feb 17, 2012 at 4:18 PM, Fecher, Gerhard <span dir="ltr"><<a href="mailto:fecher@uni-mainz.de">fecher@uni-mainz.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
your RMT's are too small, or in other words: youre structure is wrong<br>
<br>
Ciao<br>
Gerhard<br>
<br>
====================================<br>
Dr. Gerhard H. Fecher<br>
Institut of Inorganic and Analytical Chemistry<br>
Johannes Gutenberg - University<br>
55099 Mainz<br>
________________________________________<br>
Von: <a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a> [<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a>]&quot; im Auftrag von &quot;Madhav Ghimire [<a href="mailto:ghimire.mpg@gmail.com">ghimire.mpg@gmail.com</a>]<br>
Gesendet: Freitag, 17. Februar 2012 08:02<br>
An: A Mailing list for WIEN2k users<br>
Betreff: Re: [Wien] Problems in generating case.struct using shell script and overlapping of spheres with error in lapw1<br>
<div class="im"><br>
Dear Ghosh and wien users,<br>
If you leave Z blank in the structure file, then also it will be taken automatically once you save the structure file. The matter may be something else I suppose. I have tried this also as you mentioned but of no use.<br>
Anyway, one of the wien friend from Russia had asked for the structure file. I am attaching herewith its structure file. Could you please check. Please have a trial.<br>
Thanks and best regards<br>
Ghimire<br>
<br>
</div>2012/2/17 Ghosh SUDDHASATTWA <<a href="mailto:ssghosh@igcar.gov.in">ssghosh@igcar.gov.in</a><mailto:<a href="mailto:ssghosh@igcar.gov.in">ssghosh@igcar.gov.in</a>>><br>
<div class="im">I guess this problem is with R0 which has been discussed before in the mailing list. When you enter the atomic positions through w2web, do not mention Z for the atom. Leave it blank.<br>
<br>
<br>
________________________________<br>
</div>From: <a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a>> [mailto:<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a>>] On Behalf Of Madhav Ghimire<br>
<div class="im">Sent: Friday, February 17, 2012 12:14 PM<br>
<br>
To: A Mailing list for WIEN2k users<br>
Subject: Re: [Wien] Problems in generating case.struct using shell script and overlapping of spheres with error in lapw1<br>
<br>
Dear Ghosh,<br>
Yes, for the space group 166, I have the converted lattice parameters to hexagonal setting. If you are giving Rmt value approx. between 1.4 to 3.0 for Be, then you observe an error in NN distances. So, for this, you must allow for Set automatically Rmt and continue editing. This will take care of your Rmt but the question is on initialization part and error in lapw1. Please have a try on this sample (of course if you find time).<br>
Now regarding the use of w2web, yes I want to use the shell script for generating the case. struct file. Thanks for your information. But could I know where the cif file use to be (or do we need to download separately). If so, please guide.<br>
Thanks in advance<br>
M. P. Ghimire<br>
</div>2012/2/17 Ghosh SUDDHASATTWA <<a href="mailto:ssghosh@igcar.gov.in">ssghosh@igcar.gov.in</a><mailto:<a href="mailto:ssghosh@igcar.gov.in">ssghosh@igcar.gov.in</a>>><br>
<div class="im">Since the space group is 166, have you converted the lattice parameters to hexagonal setting<br>
My limited experience with w2web is that it automatically sets the RMT. It should be around 2.50 for Bi.<br>
<br>
If you do not want to use the graphical interface, you must have the cif file and use the cif2struct command.<br>
<br>
<br>
<br>
________________________________<br>
</div>From: <a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a>> [mailto:<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:wien-bounces@zeus.theochem.tuwien.ac.at">wien-bounces@zeus.theochem.tuwien.ac.at</a>>] On Behalf Of Madhav Ghimire<br>
<div class="im">Sent: Friday, February 17, 2012 11:33 AM<br>
To: A Mailing list for WIEN2k users<br>
Subject: Re: [Wien] Problems in generating case.struct using shell script and overlapping of spheres with error in lapw1<br>
<br>
Dear Jose,<br>
For the generation of structure file, Can I generate a case. struct file without using the graphical interface. In wien2k userguide, it ask us to use w2web to generate the structure file. I want to generate a structure file without using w2web (graphical interface).<br>
In the meantime, as you mentioned in your last mail, I have checked the RmtKmax ranging between 3 to 8. But then also I observe the error as:<br>
Error in LAPW1<br>
'INILPW' - can't open unit: 18<br>
'INILPW' - filename: bise.vsp<br>
'INILPW' - status: old form: formatted<br>
'LAPW1' - INILPW aborted unsuccessfully.<br>
Hope someone might have resolved this. If still there is any possibility, please kindly inform. Thanks for your kind suggestion.<br>
Regards<br>
M. P. Ghimire<br>
</div>2012/2/17 Jose Alfredo Camargo Martinez <<a href="mailto:sork8686@yahoo.es">sork8686@yahoo.es</a><mailto:<a href="mailto:sork8686@yahoo.es">sork8686@yahoo.es</a>>><br>
<div class="im">Us also appeared this error (CORE electrons leak out of MT-sphere) in earlier versions of the code (wien2k_09), was very rare, as its solution.<br>
The solution: In the construction of the case.struct using the graphical environment (w2web) do NOT enter the atomic number (Z) corresponding, simply enter only the atomic symbol. This removed the error in lstart!!!<br>
<br>
For what used old versions, if you already have new version wien2k_11.<br>
<br>
Now, the error in LAPW1, could be solved by adjusting the value of * RMT Kmax, considering (RMT * Kmax) / (RMTmin) = approximately 3.5 ( check mailing-list)<br>
<br>
(excuse my bad english)<br>
José A. Camargo Martínez<br>
Estudiante Doctorado en Ciencias - Física<br>
DEPARTAMENTO DE FÍSICA<br>
CINVESTAV - MÉXICO D.F.<br>
________________________________<br>
</div>De: Madhav Ghimire <<a href="mailto:ghimire.mpg@gmail.com">ghimire.mpg@gmail.com</a><mailto:<a href="mailto:ghimire.mpg@gmail.com">ghimire.mpg@gmail.com</a>>><br>
Para: A Mailing list for WIEN2k users <<a href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.tuwien.ac.at</a>>><br>
<div class="im">Enviado: Jueves 16 de febrero de 2012 22:54<br>
Asunto: Re: [Wien] Problems in generating case.struct using shell script and overlapping of spheres with error in lapw1<br>
<br>
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>2012/2/17 Jose Alfredo Camargo Martinez <<a href="mailto:sork8686@yahoo.es">sork8686@yahoo.es</a><mailto:<a href="mailto:sork8686@yahoo.es">sork8686@yahoo.es</a>>><br>
<div class="im">One solution would be to increase the radii of MT<br>
<br>
José A. Camargo Martínez<br>
Estudiante Doctorado en Ciencias - Física<br>
DEPARTAMENTO DE FÍSICA<br>
CINVESTAV - MÉXICO D.F.<br>
________________________________<br>
</div>De: Madhav Ghimire <<a href="mailto:ghimire.mpg@gmail.com">ghimire.mpg@gmail.com</a><mailto:<a href="mailto:ghimire.mpg@gmail.com">ghimire.mpg@gmail.com</a>>><br>
Para: wien <<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>>><br>
<div><div class="h5">Enviado: Jueves 16 de febrero de 2012 22:24<br>
Asunto: [Wien] Problems in generating case.struct using shell script and overlapping of spheres with error in lapw1<br>
<br>
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: x lstart<br>
Program input is: "13 -6.0 "<br>
<br>
SELECT XCPOT:<br>
<br>
recommended: 13: GGA (Perdew-Burke-Ernzerhof 96)<br>
<br>
5: LSDA<br>
<br>
11: GGA (Wu-Cohen 2006)<br>
<br>
SELECT ENERGY to separate core and valence states:<br>
<br>
recommended: -6.0 Ry (check how much core charge leaks out of MT-sphere)<br>
<br>
WARNING: 0.328380476966373 CORE electrons leak out of MT-sphere !!!!<br>
<br>
WARNING: 1.30946564796119 CORE electrons leak out of MT-sphere !!!!<br>
<br>
WARNING: 1.30946564796119 CORE electrons leak out of MT-sphere !!!!<br>
<br>
LSTART ENDS<br>
<br>
0.537u 0.023s 0:00.59 93.2% 0+0k 8+0io 0pf+0w<br>
<br>
<br>
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>
<br>
start (Fri Feb 17 12:32:41 JST 2012) with lapw0 (40/99 to go)<br>
<br>
<br>
<br>
cycle 1 (Fri Feb 17 12:32:41 JST 2012) (40/99 to go)<br>
<br>
<br>
<br>
> lapw0 (12:32:41) 11.846u 0.116s 0:12.02 99.4% 0+0k 0+4712io 0pf+0w<br>
<br>
> lapw1 -up (12:32:53) 150.227u 0.014s 2:30.80 99.6% 0+0k 0+1088io 0pf+0w<br>
<br>
<br>
<br>
> stop error<br>
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 "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>
Thank you in advance<br>
Best regards<br>
--<br>
M. P. Ghimire<br>
<br>
_______________________________________________<br>
Wien mailing list<br>
</div></div><a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>><br>
<div class="im"><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>
Wien mailing list<br>
</div><a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>><br>
<div class="im"><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>
<br>
--<br>
M. P. Ghimire<br>
<br>
_______________________________________________<br>
Wien mailing list<br>
</div><a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>><br>
<div class="im"><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>
Wien mailing list<br>
</div><a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>><br>
<div class="im"><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>
<br>
--<br>
M. P. Ghimire<br>
<br>
_______________________________________________<br>
Wien mailing list<br>
</div><a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>><br>
<div class="im"><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>
<br>
--<br>
M. P. Ghimire<br>
<br>
_______________________________________________<br>
Wien mailing list<br>
</div><a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><mailto:<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>><br>
<div class="HOEnZb"><div class="h5"><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>
<br>
<br>
--<br>
M. P. Ghimire<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>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>M. P. Ghimire<br><br>