<div dir="ltr"><div>Dear Wien2k community,</div><div><br></div><div>     recently I have been focusing on EFG calculations and I have studied the essay written by Katrin Koch and Stefaan Cottenier (August 10, 2011) on the EFG switch in LAPW2, within it there is an example of bcc In in which the case.struct file and the energy convergence criterion along with Rmt*Kmax and the lattice parameters are given. Now after running the SCF when I look at my case.scf file the values I get for the Fermi energy and also the EFG components do not match with the predicted ones in the paper(up to a factor of 3). Since I have double checked my input files and etc. I do not know where the problem is originating from, this is the original data:</div>
<div><br></div><div>   <font lang="KO" face="CMR12"><p align="LEFT">These are the input files:</p>
<p align="LEFT">case.struct:</p>
</font><font face="CMTT12"><p align="LEFT">In</p>
<p align="LEFT">B LATTICE,NONEQUIV.ATOMS: 1139 I4/mmm</p>
<p align="LEFT">MODE OF CALC=RELA unit=ang</p>
<p align="LEFT">6.145960 6.145960 9.346780 90.000000 90.000000 90.000000</p>
<p align="LEFT">ATOM -1: X=0.00000000 Y=0.00000000 Z=0.00000000</p>
<p align="LEFT">MULT= 1 ISPLIT= 8</p>
<p align="LEFT">In NPT= 781 R0=0.00001000 RMT= 2.5000 Z: 49.0</p>
<p align="LEFT">LOCAL ROT MATRIX: 1.0000000 0.0000000 0.0000000</p>
<p align="LEFT">0.0000000 1.0000000 0.0000000</p>
<p align="LEFT">0.0000000 0.0000000 1.0000000</p>
<p align="LEFT">0 NUMBER OF SYMMETRY OPERATIONS</p>
</font><font lang="KO" face="CMR12"><p align="LEFT">case.inst:</p>
</font><font face="CMTT12"><p align="LEFT">In 1</p>
<p align="LEFT">Kr 4 5</p>
<p align="LEFT">4, 2,2.0 N</p>
<p align="LEFT">4, 2,2.0 N</p>
<p align="LEFT">4,-3,3.0 N</p>
<p align="LEFT">4,-3,3.0 N</p>
<p align="LEFT">5,-1,1.0 N</p>
<p align="LEFT">5,-1,1.0 N</p>
<p align="LEFT">5, 1,1.0 N</p>
<p align="LEFT">5, 1,0.0 N</p>
<p align="LEFT">*** End of Input</p>
<p align="LEFT">*** End of Input</p>
</font><font lang="KO" face="CMR12"></font><p align="LEFT"><font lang="KO" face="CMR12">Run this case to self consistency using GGA (PBE) and with a separation energy of </font><i><font lang="JA" face="CMSY10">−</font></i><font lang="KO" face="CMR12">6.5</font></p>
<font lang="KO" face="CMR12">
</font><p align="LEFT"><font lang="KO" face="CMR12">Ry in lstart. Use e.g. R</font><font face="CMR8" size="1"><font face="CMR8" size="1">mt</font></font><font lang="KO" face="CMR12">*K</font><font face="CMR8" size="1"><font face="CMR8" size="1">max</font></font><font lang="KO" face="CMR12">=7 and 40 k-points in the IBZ (=400 points as input for</font></p>
<font lang="KO" face="CMR12">
<p align="LEFT">kgen), which is not too accurate but very fast for a try-out. With these settings and an</p>
</font><p align="LEFT"><font lang="KO" face="CMR12">energy convergence criterion of 10</font><i><font lang="JA" face="CMSY8" size="1"><font lang="JA" face="CMSY8" size="1">−</font></font></i><font face="CMR8" size="1"><font face="CMR8" size="1">6 </font></font><font lang="KO" face="CMR12">Ry, you will end up with :EFG001=3</font><i><font lang="JA" face="CMMI12">.</font></i><font lang="KO" face="CMR12">74 </font><i><font lang="JA" face="CMSY10">・ </font></i><font lang="KO" face="CMR12">10</font><font face="CMR8" size="1"><font face="CMR8" size="1">21 </font></font><font lang="KO" face="CMR12">V/m</font><font face="CMR8" size="1"><font face="CMR8" size="1">2</font></font><font lang="KO" face="CMR12">.</font></p>
<font lang="KO" face="CMR12">
</font><p align="LEFT"><font lang="KO" face="CMR12">Inspection of the following part of case.scf learns that there are six 4p-electrons at </font><i><font lang="JA" face="CMSY10">−</font></i><font lang="KO" face="CMR12">5.02</font></p>
<font lang="KO" face="CMR12">
</font><p align="LEFT"><font lang="KO" face="CMR12">Ry, 10 4d-electrons at about </font><i><font lang="JA" face="CMSY10">−</font></i><font lang="KO" face="CMR12">0.65 Ry and some tendency to localize s-electrons at -0.24</font></p>
<font lang="KO" face="CMR12">
<p>Ry (the Fermi energy :FER=0.41):</p></font></div><div><br></div><div>so I should end up with an EFG of 3.74*10^21 and FER of 0.41, but this is not the case.</div><div><br></div><div>I truly appreciate your assistance</div>
<div><br></div><div>with regards,</div><div><br></div><div>Farshad</div><div><b><font face="CMSSBX10" size="7"><font face="CMSSBX10" size="7"><p align="LEFT"><br></p>
</font></font></b><font face="CMR17" size="5"><font face="CMR17" size="5"><p align="LEFT"><br></p></font></font></div></div>