<div dir="ltr">Let me add a little to what Fabien said.<div><br></div><div>In a +U calculation, and also -eece (and perhaps some others) it is possible:</div><div>a) To have different local minima with different spin/orbital states</div><div>b) To get stuck in a trap.</div><div><br></div><div>Case a) is "real" for the functional used. While local minimum are not present in all systems, sometimes they are, for instance low versus high spin for Cr, and parity of the density matrix eigenvectors for DyScO3. It is possible that you have different spin states in your runs.</div><div><br></div><div>Case b) is when :DIS and :PLA appear to have converged, but the density matrix or orbital potential (:MV) has not. Changes in the orbital terms increase :DIS, so appear to be a minimum if that is all you look at, but in reality it is not because the orbital tems are not self-consistent. With 17.1 I added a better scaling of the orbital potential which reduces the possibility of getting stuck in a trap.</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jul 31, 2018 at 6:05 AM, <a href="mailto:tran@theochem.tuwien.ac.at">tran@theochem.tuwien.ac.at</a> <span dir="ltr"><<a href="mailto:tran@theochem.tuwien.ac.at" target="_blank">tran@theochem.tuwien.ac.at</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">So, now you should check what results you obtain with 17.1 when<br>
starting the calculation with the files (.clm/sum/up/dn, .dmat/up/dn,<br>
.vorb/up/dn,) obtained with 16.1.<br>
<br>
<br>
On Tuesday 2018-07-31 12:31, Chouaib AHMANI FERDI wrote:<br>
<br>
>Date: Tue, 31 Jul 2018 12:31:11<br>
<span class="">>From: Chouaib AHMANI FERDI <<a href="mailto:ahmaniferdichouaib@gmail.com">ahmaniferdichouaib@gmail.com</a>><br>
>Reply-To: A Mailing list for WIEN2k users <<a href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.tuwien.ac.<wbr>at</a>><br>
>To: A Mailing list for WIEN2k users <<a href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.tuwien.ac.<wbr>at</a>><br>
</span>>Subject: Re: [Wien] Different GGA+U results for different wien2k versions<br>
<div><div class="h5">><br>
>Dear Mr Tran,<br>
><br>
>Thank you for your response.<br>
><br>
>I did not start the calculations on Version_17.1 using old files from version_16.1<br>
>The only file in common is the .struct file<br>
><br>
>the same parameters are chosen for the calculations in both versions, namely : RKmax, Energy seperation, K-points, approximation : GGA-PBE and SAME U VALUE (J =<br>
>0)<br>
><br>
>Sincerely,<br>
><br>
>On Tue, Jul 31, 2018 at 11:26 AM Chouaib AHMANI FERDI <<a href="mailto:ahmaniferdichouaib@gmail.com">ahmaniferdichouaib@gmail.com</a>> wrote:<br>
>      Dear M. Marks<br>
><br>
>Thank you for your response.<br>
><br>
>VERSION_16.1 :<br>
>the bottom of the file .in1 :<br>
>---<br>
>K-VECTORS FROM UNIT:4  -12.2       1.5  1125   emin / de (emax=Ef+de) / nband<br>
>---<br>
>i didn't modify the lower value<br>
>The energy to separate core and valence states is set to -9.2 Ry in all of the calculations because of P atoms<br>
><br>
>the files .output1up_1 .output1up_2 .output1up_3 and .output1dn_1 ...  gives :<br>
>---<br>
>  0 EIGENVALUES BELOW THE ENERGY  -12.20000<br>
>---<br>
>multiple times as you go through the file<br>
><br>
>VERSION_17.1<br>
>the bottom of the file .in1 :<br>
>---<br>
>K-VECTORS FROM UNIT:4  -12.2       1.5  1101   emin / de (emax=Ef+de) / nband<br>
>---<br>
>the files .output1up_1 .output1up_2 .output1up_3 and .output1dn_1 ...  gives the same as in VERSION_16.1<br>
>---<br>
>  0 EIGENVALUES BELOW THE ENERGY  -12.20000<br>
>---<br>
><br>
>I noticed that nband is a little bit different in the two versions<br>
>I guess there's no ghostbands<br>
><br>
>what so you think ?<br>
><br>
>Sincerely,<br>
><br>
><br>
><br>
><br>
>On Sun, Jul 29, 2018 at 3:58 PM <<a href="mailto:tran@theochem.tuwien.ac.at">tran@theochem.tuwien.ac.at</a>> wrote:<br>
>      Hi,<br>
><br>
>      Did you start the calculation with VERSION_17.1 using the<br>
>      files form the old calculation?<br>
><br>
>      F. Tran<br>
><br>
>      On Sunday 2018-07-29 16:21, Chouaib AHMANI FERDI wrote:<br>
><br>
>      >Date: Sun, 29 Jul 2018 16:21:12<br>
>      >From: Chouaib AHMANI FERDI <<a href="mailto:ahmaniferdichouaib@gmail.com">ahmaniferdichouaib@gmail.com</a>><br>
>      >Reply-To: A Mailing list for WIEN2k users <<a href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.tuwien.ac.<wbr>at</a>><br>
>      >To: A Mailing list for WIEN2k users <<a href="mailto:wien@zeus.theochem.tuwien.ac.at">wien@zeus.theochem.tuwien.ac.<wbr>at</a>><br>
>      >Subject: [Wien] Different GGA+U results for different wien2k versions<br>
>      ><br>
>      >Greetings Wien2k users,<br>
>      ><br>
>      >I have some results of a compound from Wien2k_16.1 a long time ago with GGA+U <br>
>      >When I tried the same calculations on the sam compound with  VERSION_17.1 I found different results.<br>
>      >for example : I had a Band GAP of 1.4 eV with GGA+U U=7eV in Version_16.1<br>
>      >while i have 3.8 eV with GGA+U U=7eV in Version_17.1, I noticed that this value is close to VASP results on the same structure same U value<br>
>      and slightly<br>
>      >different parameters<br>
>      >This remark is valid for another compound where the Version_17.1 gives band gap closer to VASP results for the same value of U<br>
>      ><br>
>      >Another issue is that when I try to calculate the Density Of States, some low energy bands (below 8 eV) do not show, (maybe they are not<br>
>      calculated) while<br>
>      >Version_16.1 they were shown<br>
>      ><br>
>      >Could someone explain ?<br>
>      ><br>
>      >Sincerely,<br>
>      ><br>
>      ><br>
>      >C.Ferdi<br>
>      ><br>
>      ><br>
>      >--<br>
>      >-----------------------------<wbr>------------------------<br>
>      >AHMANI FERDI Chouaïb<br>
>      >Laboratoire Matériaux Nanomatériaux <br>
>      >Nanomagnétisme et Enseignement des Sciences<br>
>      >École Normale Supérieure<br>
>      >Université Mohammed V, Rabat. <br>
>      >Tel : +212 6 94 59 57 60<br>
>      >[signature.gif?u=139409&e=<wbr>29384200&v=<wbr>75b891d71c0edac28160e11bdd3163<wbr>23c4ca280fbfd541add118cf56dd87<wbr>78ff]<br>
>      >_____________________________<wbr>__________________<br>
>      Wien mailing list<br>
>      <a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.<wbr>at</a><br>
</div></div>>      <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.theochem.tuwien.ac.at_mailman_listinfo_wien&d=DwIDaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=ir72_UZvL551lvhmmDxgf_qwJnyRG3w7Dvz38pVu6kQ&s=u8IhsIgnXLrIBBRg9Noo_dtoWNpM0_VGaO3vqwD7Qao&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.<wbr>com/v2/url?u=http-3A__zeus.<wbr>theochem.tuwien.ac.at_mailman_<wbr>listinfo_wien&d=DwIDaQ&c=<wbr>yHlS04HhBraes5BQ9ueu5zKhE7rtNX<wbr>t_d012z2PA6ws&r=U_<wbr>T4PL6jwANfAy4rnxTj8IUxm818jnvq<wbr>KFdqWLwmqg0&m=ir72_<wbr>UZvL551lvhmmDxgf_<wbr>qwJnyRG3w7Dvz38pVu6kQ&s=<wbr>u8IhsIgnXLrIBBRg9Noo_dtoWNpM0_<wbr>VGaO3vqwD7Qao&e=</a><br>
>      SEARCH the MAILING-LIST at:  <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.mail-2Darchive.com_wien-40zeus.theochem.tuwien.ac.at_index.html&d=DwIDaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=ir72_UZvL551lvhmmDxgf_qwJnyRG3w7Dvz38pVu6kQ&s=iwzSm_t804NQzBh7hBfVNMNRsdsjW3h8SRN5rDg6_EA&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.<wbr>com/v2/url?u=http-3A__www.<wbr>mail-2Darchive.com_wien-<wbr>40zeus.theochem.tuwien.ac.at_<wbr>index.html&d=DwIDaQ&c=<wbr>yHlS04HhBraes5BQ9ueu5zKhE7rtNX<wbr>t_d012z2PA6ws&r=U_<wbr>T4PL6jwANfAy4rnxTj8IUxm818jnvq<wbr>KFdqWLwmqg0&m=ir72_<wbr>UZvL551lvhmmDxgf_<wbr>qwJnyRG3w7Dvz38pVu6kQ&s=iwzSm_<wbr>t804NQzBh7hBfVNMNRsdsjW3h8SRN5<wbr>rDg6_EA&e=</a><br>
<div class="HOEnZb"><div class="h5">><br>
><br>
><br>
>--<br>
>-----------------------------<wbr>------------------------<br>
>AHMANI FERDI Chouaïb<br>
>Laboratoire Matériaux Nanomatériaux <br>
>Nanomagnétisme et Enseignement des Sciences<br>
>École Normale Supérieure<br>
>Université Mohammed V, Rabat. <br>
>Tel : +212 6 94 59 57 60<br>
><br>
><br>
><br>
>--<br>
>-----------------------------<wbr>------------------------<br>
>AHMANI FERDI Chouaïb<br>
>Laboratoire Matériaux Nanomatériaux <br>
>Nanomagnétisme et Enseignement des Sciences<br>
>École Normale Supérieure<br>
>Université Mohammed V, Rabat. <br>
>Tel : +212 6 94 59 57 60<br>
><br>
></div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><span style="font-size:12.8px">Professor Laurence Marks</span><br></div><div dir="ltr"><span style="font-size:12.8px">"Research is to see what everybody else has seen, and to think what nobody else has thought", </span><span style="font-size:12.8px">Albert Szent-Gyorgi</span><br><a href="http://www.numis.northwestern.edu" target="_blank">www.numis.northwestern.edu</a> ; <span style="font-size:12.8px">Corrosion in 4D: </span><a href="http://MURI4D.numis.northwestern.edu" style="font-size:12.8px" target="_blank">MURI4D.numis.northwestern.edu</a><div><span style="font-size:12.8px">Partner of the CFW 100% program for gender equity, </span><a href="http://www.cfw.org/100-percent" style="font-size:12.8px" target="_blank">www.cfw.org/100-percent</a></div><div>Co-Editor, Acta Cryst A</div></div></div></div></div></div></div></div></div></div></div>
</div>