<div dir="ltr">Dear all, I am simulating compounds with rare earths within a collinear calculation. I tried to suggest some occupations of the f orbitals for the Tb atom, using the starting_ns_eigenvalue(m,ispin,I). Unfortunately also with a small value of electronic mixing and big values of U, the imposed occupations was lost during the iterative cycle. So I used also mixing_fixed_ns > electron_maxstep. Also in this case after I obtain the total energy convergence the density matrix of the last iteration is different to the one imposed in the input_file. <div>Can someone of you kindly explain me how the mixing_fixed_ns works ? Do it uses a sort of Lagrange multipliers to force the occupations or is more like a penalty function ? Finally, is it possible to work around the problem, converging to a determined occupation ?<br></div><div><br></div><div>Thanks in advance for every suggestion,</div><div><br></div><div><div><div><div class="gmail_signature"><div dir="ltr"><font color="#888888">Simone Marocchi</font><br style="color:rgb(136,136,136)"><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">S3 Center, Istituto Nanoscienze, CNR</span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">via Campi 213/A, 41125, Modena, Italy</span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">Tel: </span><a value="+390592055300" style="color:rgb(17,85,204)">+39 </a><u>0592055585</u><span style="color:rgb(136,136,136)">;  Skype: jacobi84</span><br style="color:rgb(136,136,136)"><span style="color:rgb(136,136,136)">URL: </span><a href="http://www.nano.cnr.it/" style="color:rgb(17,85,204)" target="_blank">http://www.nano.cnr.it</a><br style="color:rgb(136,136,136)"></div></div>
</div></div></div></div>