<div dir="ltr">Dear Wien2k users,<div><br></div><div>I am doing calculations of a hexaboride films invoking: runsp_c_lapw -so -eece</div><div>My problem is two fold:<br><br></div><div>(i) SCF calculations successfully converged for thin films (9 atomic units) and bulk systems. I do have a problem to converge a thicker film (11 atomic units) even TEMP broadening of 0.02 is used. The observations after 800 iterations are following:</div>
<div><br></div><div>(1) fermi energy is jumping (roughly) between two values</div><div><div>:FER  : F E R M I - ENERGY(FERMI-SM.)=   0.3682855520</div><div>:FER  : F E R M I - ENERGY(FERMI-SM.)=   0.3682855520</div><div>:FER  : F E R M I - ENERGY(FERMI-SM.)=  -0.7779431003</div>
<div>:FER  : F E R M I - ENERGY(FERMI-SM.)=  -0.7779431003</div><div>:FER  : F E R M I - ENERGY(FERMI-SM.)=   0.3683943150</div><div>:FER  : F E R M I - ENERGY(FERMI-SM.)=   0.3683943150</div><div>:FER  : F E R M I - ENERGY(FERMI-SM.)=  -0.7767634505</div>
<div>:FER  : F E R M I - ENERGY(FERMI-SM.)=  -0.7767634505</div></div><div><br></div><div>(2) charge distance non uniformly oscillates between 0.1 and 0.001</div><div><div>:DIS  :  CHARGE DISTANCE       ( 1.1067771 for atom    4 spin 1)      0.1064425</div>
<div>:DIS  :  CHARGE DISTANCE       ( 0.0059078 for atom    6 spin 1)      0.0013328</div><div>:DIS  :  CHARGE DISTANCE       ( 0.0092322 for atom    6 spin 1)      0.0015769</div><div>:DIS  :  CHARGE DISTANCE       ( 0.0092664 for atom    6 spin 1)      0.0015106</div>
<div>:DIS  :  CHARGE DISTANCE       ( 0.5478509 for atom    4 spin 1)      0.0466613</div><div>:DIS  :  CHARGE DISTANCE       ( 0.0078766 for atom    6 spin 1)      0.0014700</div><div>:DIS  :  CHARGE DISTANCE       ( 0.9814094 for atom    4 spin 1)      0.0939173</div>
<div>:DIS  :  CHARGE DISTANCE       ( 0.0070829 for atom    6 spin 1)      0.0015672</div><div>:DIS  :  CHARGE DISTANCE       ( 0.0411319 for atom    4 spin 2)      0.0053433</div><div>:DIS  :  CHARGE DISTANCE       ( 0.0072009 for atom    6 spin 1)      0.0014308</div>
<div>:DIS  :  CHARGE DISTANCE       ( 0.1357701 for atom    4 spin 2)      0.0150618</div></div><div><br></div><div>(3) magnetic moments are zero within spheres whereas there is </div><div>:MMTOT: SPIN MAGNETIC MOMENT IN CELL   =    -0.00010<br>
</div><div>:MMINT: MAGNETIC MOMENT IN INTERSTITIAL =   -0.00003<br></div><div><br></div><div>(4) my case.inm</div><div><div>BROYD  0.0   YES  (BROYD/PRATT, extra charge (+1 for additional e), norm)</div><div>0.005           mixing FACTOR for BROYD/PRATT scheme</div>
</div><div><div>1.00  1.00      PW and CLM-scaling factors</div><div>9999  8         idum, HISTORY</div></div><div><br></div><div>:MIX  :   MSEC3  REGULARIZATION: 2.50E-04  GREED: 0.005  Newton 1.000<br></div><div>with PRATT after 100 iterations</div>
<div>MIX  :   PRATT  REGULARIZATION: 2.50E-04  GREED: 0.005<br></div><div><br></div><div>Do you have any hints how to converge the calculations?</div><div><br></div><div>(ii) I am running the calculations with version 12.1. The runsp_c_lapw script has been updated (debugged) in version 13.1 but SCF diverges in my cases (thin and thicker films).</div>
<div><br></div><div>Does anybody experience similar problem?</div><div><br></div><div>Best regards,</div><div>Martin Gmitra</div><div>Uni Regensburg</div><div><br></div><div><br></div><div><br></div></div>