<div dir="ltr">Dear Wien2k team,<div><br></div><div>I am trying to apply the on-site hybrid functional method to treat the Tb 4f electrons in TbMn6Sn6, together with the SO coupling. I use</div><div><br></div><div>runsp_lapw -p -eece -so -ec 0.00001 -cc 0.0001</div><div><br></div><div>My case.ineece is this:</div><div><span style="background-color:transparent"><br></span></div><div><span style="background-color:transparent">-12.0  1       Emin natom</span></div>1  1 3      index of atom, number of l, l<br>HYBR          HYBR / EECE<br>0.25          amount of exact exchange<div><br></div><div>I would like to have a hybrid solution with the following f-orbital occupation matrices:</div><div><br></div><div>Density matrix UPUP block, real part.  L= 3<br>         0.98992  0.00000  0.00000  0.00000  0.00000  0.00000 -0.00047<br>         0.00000  0.99177  0.00000  0.00000  0.00000  0.00000  0.00000<br>         0.00000  0.00000  0.99246  0.00000  0.00000  0.00000  0.00000<br>         0.00000  0.00000  0.00000  0.99251  0.00000  0.00000  0.00000<br>         0.00000  0.00000  0.00000  0.00000  0.99261  0.00000  0.00000<br>         0.00000  0.00000  0.00000  0.00000  0.00000  0.99181  0.00000<br>        -0.00047  0.00000  0.00000  0.00000  0.00000  0.00000  0.99071</div><div> Density matrix DNDN block, real part.  L= 3<br>         0.98393  0.00000  0.00000  0.00000  0.00000  0.00000  0.01313<br>         0.00000  0.00237  0.00000  0.00000  0.00000  0.00000  0.00000<br>         0.00000  0.00000  0.00251  0.00000  0.00000  0.00000  0.00000<br>         0.00000  0.00000  0.00000  0.00303  0.00000  0.00000  0.00000<br>         0.00000  0.00000  0.00000  0.00000  0.00240  0.00000  0.00000<br>         0.00000  0.00000  0.00000  0.00000  0.00000  0.00212  0.00000<br>         0.01313  0.00000  0.00000  0.00000  0.00000  0.00000  0.00201</div><div><br></div><div>I stabilize this solution in a separate calculation with fixed occupation matrices run as</div><div><br></div><div>runsp_lapw -p -orb -so -ec 0.00001 -cc 0.0001</div><div><br></div><div>I start my hybrid (-eece) calculation from a converged fixed matrices (-orbc) calculation.  The problem is that DNDN occupation matrix changes dramatically in the very first iteration to this:</div><div>         0.49295  0.00000  0.00000  0.00000  0.00000  0.00000  0.01313<br>         0.00000  0.00223  0.00000  0.00000  0.00000  0.00000  0.00000<br>         0.00000  0.00000  0.00247  0.00000  0.00000  0.00000  0.00000<br>         0.00000  0.00000  0.00000  0.00303  0.00000  0.00000  0.00000<br>         0.00000  0.00000  0.00000  0.00000  0.00247  0.00000  0.00000<br>         0.00000  0.00000  0.00000  0.00000  0.00000  0.00223  0.00000<br>         0.01313  0.00000  0.00000  0.00000  0.00000  0.00000  0.49295</div><div><br></div><div>My mixing factor is small: 0.05 for the Pratt scheme.</div><div><br></div><div>I would appreciate your insights and suggestions.</div><div><br></div><div>Thank you,</div><div>Kateryna</div><div class="gmail-yj6qo"></div><div class="gmail-adL"><br></div></div>