<div dir="ltr"><p class="MsoNormal">Dear Wien2k authors and users</p>

<p class="MsoNormal">I’m interested to calculate total, Coulomb, and XC potential
using wien2k code. I performed the scf calculations using B3PW91 with alpha=0.2
using following command:</p>

<p class="MsoNormal">runsp_lapw -p -so -eece -in1ef -i 900 -cc 0.00001</p>

<p class="MsoNormal">Then I used lapw5 to calculate the potentials using the
following case.in5:</p>

<p class="MsoNormal">0           0    10000000      10000000</p>

<p class="MsoNormal">    10000000           0   
10000000      10000000</p>

<p class="MsoNormal">           0    10000000   
10000000      10000000</p>

<p class="MsoNormal">3 2 3</p>

<p class="MsoNormal">1000 1000</p>

<p class="MsoNormal">RHO ADD</p>

<p class="MsoNormal">ATU VAL NODEBUG</p>

<p class="MsoNormal">ORTHO</p>

<p class="MsoNormal">Then I used x lapw5 –d command and edit the lapw5.def as
follows and run lapw5 lapw5.def</p>

<p class="MsoNormal">For total potential:</p>

<p class="MsoNormal">5 ,'Case.in5', 'old',   
'formatted',0</p>

<p class="MsoNormal">6 ,'Case.output5',  
'unknown','formatted',0</p>

<p class="MsoNormal">8 ,'Case.struct',   
'old',    'formatted',0</p>

<p class="MsoNormal">9 ,'Case.vtotal',   
'old',    'formatted',0</p>

<p class="MsoNormal">10,'Case.tmp',      
'unknown','unformatted',0</p>

<p class="MsoNormal">11,'Case.vtotaldn', 
'unknown','formatted',0</p>

<p class="MsoNormal">12,'Case.sigma',    
'unknown','formatted',0</p>

<p class="MsoNormal">20,'Case.rho_onedim','unknown','formatted',0</p>

<p class="MsoNormal">21,'Case.rho',      
'unknown','formatted',0</p>

<p class="MsoNormal"><br></p><p class="MsoNormal">For Coulomb potential:</p>

<p class="MsoNormal">5 ,'Case.in5', 'old',   
'formatted',0</p>

<p class="MsoNormal">6 ,'Case.output5',  
'unknown','formatted',0</p>

<p class="MsoNormal">8 ,'Case.struct',   
'old',    'formatted',0</p>

<p class="MsoNormal">9 ,'Case.vcoul',   
'old',    'formatted',0</p>

<p class="MsoNormal">10,'Case.tmp',      
'unknown','unformatted',0</p>

<p class="MsoNormal">11,'Case.vcoul', 
'unknown','formatted',0</p>

<p class="MsoNormal">12,'Case.sigma',    
'unknown','formatted',0</p>

<p class="MsoNormal">20,'Case.rho_onedim','unknown','formatted',0</p>

<p class="MsoNormal">21,'Case.rho',      
'unknown','formatted',0</p>

<p class="MsoNormal"><br></p><p class="MsoNormal">And for XC potential:</p>

<p class="MsoNormal">5 ,'Case.in5', 'old',   
'formatted',0</p>

<p class="MsoNormal">6 ,'Case.output5',  
'unknown','formatted',0</p>

<p class="MsoNormal">8 ,'Case.struct',   
'old',    'formatted',0</p>

<p class="MsoNormal">9 ,'Case.r2v',   
'old',    'formatted',0</p>

<p class="MsoNormal">10,'Case.tmp',      
'unknown','unformatted',0</p>

<p class="MsoNormal">11,'Case.r2vdn', 
'unknown','formatted',0</p>

<p class="MsoNormal">12,'Case.sigma',    
'unknown','formatted',0</p>

<p class="MsoNormal">20,'Case.rho_onedim','unknown','formatted',0</p>

<p class="MsoNormal">21,'Case.rho',      
'unknown','formatted',0</p>

<p class="MsoNormal">I expected that the range of Coulomb potential must be very
higher than XC potential and XC term is very little. However, after potential
calculations I saw that between atoms the range of Coulomb potential is very
close to XC potential.  A small positive Coulomb
potential and small negative XC potential leads to the very small negative
total potential.</p>

<p class="MsoNormal">Could you help me please to know that whether are my
calculations true or not?</p>

<p class="MsoNormal">Is this correct that XC and Coulomb terms are in the same
range?</p>

<p class="MsoNormal"> </p>

<p class="MsoNormal"> </p></div>