<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Dear Dr. Blaha,</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Thank you so much for your great help!</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
I was unaware of the new feature HELOs and HDLO. I just downloaded WIen19 version. I will try to use them in my calculation. I am sorry that I forgot to mention that I need states both a few eV above the Fermi level up to several Rydberg above the Fermi level,
which covers 130 bands from 3s,3p above. I want to describe all of them accurately. My system is FeNi3. I created my own weight file (which is read into my modified wien code, see an example below) and nearly all the band states are partially occupied.
You are correct in general APW is quite flexible, but not in my case with so many states . I tested several compounds and found that LAPW basis leads to fewer FERMI errors, but I am not entirely whether this is due to the difference between APW and LAPW
basis or not. Also my case.in1 does not lead to the same Fermi error message when I work on the ground state calculation or fewer states occupied in conduction bands. 3d-lo is not needed in the beginning. It is when I have more states occupied in the conduction
band range, where this error message appears. </div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Here is my </div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
case.weightdn</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
1st k point</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span> ENERGY OCC wien weight </span></div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span> -7.028525570298 0.99999655769771D+00 0.13717373905318D-02<br>
</span>
<div> -7.026380887859 0.99999774915802D+00 0.13717390249081D-02<br>
</div>
<div> -7.026364863146 0.99999641420436D+00 0.13717371936960D-02<br>
</div>
<div> -6.963105821346 0.99999684097144D+00 0.13717377791103D-02<br>
</div>
<div> -6.960884433752 0.99999241965826D+00 0.13717317142089D-02<br>
</div>
<div> -6.960863985385 0.99999681157351D+00 0.13717377387840D-02<br>
</div>
<div> -5.955580145632 0.99999727476037D+00 0.13717383741569D-02<br>
</div>
<span></span>...</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span> 0.551102177130 0.99654806344721D+00 0.13670069457438D-02<br>
</span>
<div> 0.556450855740 0.99979910735537D+00 0.13714665395821D-02<br>
</div>
<div> 0.557215914130 0.99891542054497D+00 0.13702543491701D-02<br>
</div>
<div> 0.598361667748 0.86195838931551D+00 0.11823846218320D-02 <---Fermi energy<br>
</div>
<div> 0.614426820474 0.21955408776520D-01 0.30117158815528D-04<br>
</div>
<div> 0.616258443618 0.61299145119382D-02 0.84086618819453D-05<br>
</div>
<div> 0.620265725818 0.63525870932193D-04 0.87141112389839D-07<br>
</div>
<div> 0.671961904610 0.20208405018240D-03 0.27720720189629D-06<br>
</div>
<div>....</div>
<div><span> 2.523661999828 0.60120208311812D-04 0.82469421552554D-07<br>
</span>
<div> 2.535036138282 0.77696873303722D-04 0.10658007311896D-06<br>
</div>
<div> 2.557885989252 0.19321304918263D-03 0.26503847624503D-06<br>
</div>
<div> 2.561145153263 0.16466297198602D-03 0.22587513303980D-06<br>
</div>
<div> 2.570553631703 0.47163224554491D-04 0.64695781281880D-07<br>
</div>
<div> 2.573893684513 0.96953048385584D-04 0.13299457940409D-06<br>
</div>
<div> 2.764459110536 0.14599878318928D-04 0.20027267927199D-07<br>
</div>
<div> 2.765906180448 0.17919946692148D-04 0.24581545531067D-07<br>
</div>
<div> 2.786646061603 0.12925884883145D-04 0.17730980635316D-07<br>
</div>
<div> 2.787413941376 0.40356598962982D-05 0.55358846314104D-08<br>
</div>
<div> 2.790360218136 0.35673794370040D-04 0.48935246049437D-07<br>
</div>
<div> 2.808193555479 0.17264394477140D-04 0.23682296950809D-07</div>
</div>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
One side mark on WIEN update: Is it possible that you keep those old comments and remarks inside the modified codes? I found the new version has fewer comment lines. </div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
I am grateful if you give me further comments and suggestions.</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Thanks a lot again!</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Best regards,</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Guoping</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div>
<div id="appendonsend"></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> Wien <wien-bounces@zeus.theochem.tuwien.ac.at> on behalf of Peter Blaha <pblaha@theochem.tuwien.ac.at><br>
<b>Sent:</b> Wednesday, July 24, 2019 4:40 PM<br>
<b>To:</b> wien@zeus.theochem.tuwien.ac.at <wien@zeus.theochem.tuwien.ac.at><br>
<b>Subject:</b> Re: [Wien] Larger basis at higher energy scales</font>
<div> </div>
</div>
<div class="BodyFragment"><font size="2"><span style="font-size:11pt">
<div class="PlainText">Hi,<br>
<br>
I don't know which system it is, but definitely the case.in1 file cannot <br>
work.<br>
Also a couple of other statements are definitely wrong, see below.<br>
<br>
> I am interested in the constraint excited states calculation a few eV <br>
> about the Fermi level, so I have to add a few local orbitals at those <br>
> higher energies (see my case.in1 below). I also use LAPW basis instead <br>
> of APW (which is not flexible for delocalized states). But since wien2k <br>
> is restricted to one local basis for each atom, I can not increase <br>
> anymore. <br>
<br>
Why do you think that APW is not flexible for delocalized states ????<br>
I don't think this is true.<br>
<br>
If you just want states "a few eV" above EF, probably nothing is <br>
necessary, but eventually you could add a HDLO.<br>
<br>
Why do you think WIEN2k is restricted to one local basis function/atom <br>
??? Since a couple of years one can add multiple LOs (see eg. the NMR <br>
code), but they have to be sufficiently separated in energy .<br>
<br>
As a result, WIEN2k gave me<br>
> FERMI - Error. The program stops at LAPW2, with uplapw2.error as<br>
> Error in LAPW2<br>
> 'FERMI' - EFERMI OUT OF ENERGY RANGE<br>
> 'FERMI' - STOP IN EFI<br>
> 'FERMI' - ENERGY OF LOWER BOUND : 0.59917<br>
> 'FERMI' - NUMBER OF STATES AT THE LOWER BOUND : 70.04360<br>
> 'FERMI' - ENERGY OF UPPER BOUND : 0.59917<br>
> 'FERMI' - NUMBER OF STATES AT THE UPPER BOUND : 70.04377<br>
> 'FERMI' - ADD 69.72319<br>
> 'FERMI' - SOS 0.0000.0000.0000.0000.0000.0000.0000.0000.0000.000<br>
> 'FERMI' - NOS **************************************************<br>
> ** testerror: Error in Parallel LAPW2<br>
<br>
This definitely comes because of a wrong in1 file.<br>
> <br>
> The following is my case.in1<br>
> <br>
> WFFIL EF=.599146891400 (WFFIL, WFPRI, ENFIL, SUPWF)<br>
> 9.00 10 4 (R-MT*K-MAX; MAX L IN WF, V-NMT<br>
> 0.30 6 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global <br>
> APW/LAPW)<br>
> 0 0.30 0.000 CONT 0<br>
> 0 -6.50 0.001 STOP 0<br>
> 1 0.30 0.000 CONT 0<br>
> 1 -3.97 0.001 STOP 0<br>
> 2 0.30 0.005 CONT 0<br>
> 2 1.00 0.000 CONT 0 --> this cannot work. Put a HDLO<br>
> 0.30 6 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global <br>
> APW/LAPW)<br>
> 0 0.30 0.000 CONT 0<br>
> 0 -7.93 0.001 STOP 0<br>
> 1 0.30 0.000 CONT 0<br>
> 1 -4.96 0.001 STOP 0<br>
> 2 0.30 0.005 CONT 0<br>
> 2 1.20 0.000 CONT 0 --> HDLO<br>
> 0.30 6 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global <br>
> APW/LAPW)<br>
> 0 0.30 0.000 CONT 0<br>
> 0 -7.93 0.001 STOP 0<br>
> 1 0.30 0.000 CONT 0<br>
> 1 -4.96 0.001 STOP 0<br>
> 2 0.30 0.005 CONT 0<br>
> 2 1.20 0.000 CONT 0 --> HDLO<br>
> K-VECTORS FROM UNIT:4 -12.0 5.0 145 emin / de (emax=Ef+de) / <br>
> nband<br>
<br>
You cannot put a d-LO so close to the normal APW d-energy parameter, but <br>
if you think that the linearization of the d-states may be insufficient, <br>
you can add a HDLO.<br>
<br>
-- <br>
--------------------------------------------------------------------------<br>
Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna<br>
Phone: +43-1-58801-165300 FAX: +43-1-58801-165982<br>
Email: blaha@theochem.tuwien.ac.at WIEN2k: <a href="http://www.wien2k.at">http://www.wien2k.at</a><br>
WWW: <br>
<a href="http://www.imc.tuwien.ac.at/tc_blaha-------------------------------------------------------------------------">http://www.imc.tuwien.ac.at/tc_blaha-------------------------------------------------------------------------</a>
<br>
<br>
_______________________________________________<br>
Wien mailing list<br>
Wien@zeus.theochem.tuwien.ac.at<br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at: <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html">
http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</div>
</span></font></div>
</div>
</body>
</html>