<div dir="ltr"><div><div><br></div>Dear Simone<br><br></div>I try to answer your questions below.<br><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Nov 7, 2014 at 3:32 PM, simone marocchi <span dir="ltr"><<a href="mailto:simone.roz@gmail.com" target="_blank">simone.roz@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><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). </div></blockquote><div><br><br></div><div>why do you need to do that? if you want to force different values of occupation on states that are equivalent by symmetry you will not be successful. If this is the case you have to do something that makes the system loose that symmetry operation connecting the two states. <br></div><div><br> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Unfortunately also with a small value of electronic mixing and big values of U, the imposed occupations was lost during the iterative cycle. </div></blockquote><div><br><br></div><div>this can happen. a dirty trick I learned recently is to suggest a starting_ns_eigenvalue(m,<div dir="ltr">ispin,I) bigger than one in input (1.2 or 1.3 maybe). This does not make any physical sense of course. However the Hubbard potential becomes more attractive for the specific eigenvector of the occupation matrix you want to fill completely and the code takes more time to "come back" to a physical value. If that is a state it likes (at least a local minimum of the energy) it might fall into it. Of course you have to check at the end that the occupation has gone back to a value <= 1.<br></div></div><div><br> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">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></blockquote><div><br><br></div><div>this is (possibly) strange. How different is it? How well are you converging?<br></div><div>The fact that it is different is not surprising: the routine that prints the ns, always prints the ones that are computed from the KS states just obtained from the new diagonalization. These can be different from the ones e.g. used in contructing the Hubbard potential.<br></div><div>However if your calculation is well converged this difference should not be big: both KS wfcs and their occupations should be converged reasonably well.<br></div><div>If this does not happen and you still see a significant difference maybe it means that the values you are trying to impose is not consistent with what the system wants to do (e.g., you have less symmetry than the code finds for your crystal)<br></div><div><br><br> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><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></blockquote><div><br><br></div><div>no with the current version of the code. you could implement some (e.g., quadratic) constraint and try. I think I once tried (long time ago) and seem to remember problems in convergence.<br><br></div><div>Best,<br><br>Matteo<br><br></div><div><br> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div></div><div><br></div><div>Thanks in advance for every suggestion,</div><div><br></div><div><div><div><div><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>
<br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br></div></div>