<div dir="ltr"><div><div><div><div><div>Dear Prof. Marks,<br><br></div>   Thank you very much for the explanations. They will be useful now and in the future.<br></div>   My last query (I promise !!): Could you, please, clarify the meaning of the last two columns in a :MIX line ? For example :<br><br>:MIX  :   MSE1   REGULARIZATION: 4.93E-06  GREED: 0.437  Newton 1.00  1.14<br><br></div>   Thank you again !<br></div>   All the best,<br></div>                   Luis<br><br></div><div class="gmail_extra"><br><div class="gmail_quote">2017-01-26 11:01 GMT-02:00 Laurence Marks <span dir="ltr"><<a href="mailto:L-marks@northwestern.edu" target="_blank">L-marks@northwestern.edu</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="auto"><div>A clarification on some terms/types. For certain -eece and probably also +U can stop at what should be called "traps". This is when :DIS and the plane waves (:PLA) are self-consistent, but :MV is not. You can also have a trap where :DIS is small but :PLA is not, and if you start from a converged density but the forces are large. The latest mixer does a better job of escaping these traps. I would not be surprised to find that there are published results which are really traps. (Not just with Wien2k, with other DFT codes as well.)<div dir="auto"><br></div><div dir="auto">Local minima are "real" and "right". Wien2k (and most other codes) will stop at local minima, only a few search for the global minimum. For certain one can have local minima for spin states; it is not impossible to have local minima for atomic positions. MSR1(a) can escape from small traps, i.e. ones which are only present for a small range of densities/positions. Sometimes small traps are due to noise in the algorithms.</div><div dir="auto"><br></div><div dir="auto">Your InP case could find a local minimum, or a trap. Traps you can check by looking to see if terms in addition to :DIS are converged.</div><div dir="auto"><br></div><div dir="auto">mBJ is a special case that I don't understand -- I have not worked with it enough to know if it has traps. I also don't know if the forces with mBJ can be trusted. From some simple tests I did I suspect that they can't, but I am not certain.</div><div dir="auto"><br></div><div dir="auto">To escape a trap or a small local minimum you don't want to add SLOW, decrease the GREED etc. -- use the larger default parameters.</div><div dir="auto"><br></div><div dir="auto">To escape true local minima (spin or positions) you have to change to a new starting point that is outside the radius of convergence (in spin/position/density space) of your earlier calculation. There are some emails on the list about changing spin states by hand.<br><br><div data-smartmail="gmail_signature" dir="auto">---<span class=""><br>Professor Laurence Marks<br>"Research is to see what everybody else has seen, and to think what nobody else has thought", Albert Szent-Gyorgi<br></span><a href="http://www.numis.northwestern.edu" target="_blank">http://www.numis.northwestern.<wbr>edu</a><br>Corrosion in 4D <a href="http://MURI4D.numis.northwestern.edu" target="_blank">http://MURI4D.numis.<wbr>northwestern.edu</a><br>Partner of the CFW 100% gender equity project, <a href="http://www.cfw.org/100-percent" target="_blank">www.cfw.org/100-percent</a><br>Co-Editor, Acta Cryst A<br>    </div></div><br><div class="gmail_extra"><br><div class="gmail_quote"><span class="">On Jan 26, 2017 04:18, "Luis Ogando" <<a href="mailto:lcodacal@gmail.com" target="_blank">lcodacal@gmail.com</a>> wrote:<br type="attribution"></span><blockquote class="m_-1646327723223284828quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">



<div><span class="">
<div dir="ltr">
<div>
<div>
<div>
<div>
<div>
<div>
<div>Dear Prof. Marks,<br>
<br>
</div>
   Thank you very much for your answers ! I am pretty sure that your "intuition" will save a lot of computation time !<br>
</div>
   Just one last question: I have another system , a supercell formed by 15 InP zinc blend cells along [111] (hexagonal representation of the cubic lattice) and 3 InP wurtzite cells. They are aligned along the hexagonal "c" axis.<br>
</div>
                              <wbr>       Comparing the gaps of the respective bulks and this supercell (same calculation parameters) , I believe that the SCF cycle (here, it is not a lattice optimization) stopped at a local minimum.<br>
</div>
                              <wbr>        My question is : in this case (regular SCF cycle and InP cells), would you change any of your previous answers ?<br>
</div>
   Thank you again !!<br>
</div>
   All the best,<br>
</div>
                 Luis<br>
<br>
</div>
</span><div class="gmail_extra"><br>
<div class="gmail_quote"><span class=""><div class="m_-1646327723223284828elided-text">2017-01-25 17:56 GMT-02:00 Laurence Marks <span dir="ltr">
<<a href="mailto:L-marks@northwestern.edu" target="_blank">L-marks@northwestern.edu</a>></span>:<br>
</div></span><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div class="m_-1646327723223284828elided-text">
<div dir="ltr">Inlined is my intuition, which does not have to be completely right.
<div class="gmail_extra"><br>
<div class="gmail_quote"><span>On Wed, Jan 25, 2017 at 11:32 AM, Luis Ogando
<span dir="ltr"><<a href="mailto:lcodacal@gmail.com" target="_blank">lcodacal@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
<div dir="ltr">
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>Dear Prof. Marks (and Wien2k community),<br>
<br>
</div>
   After a recent discussion about "difficult" optimizations in this mailing list (subject: "Mixer surprise when using PBE0 hybrid on-site functional"), I would like ask you for an advice.<br>
</div>
   I have a system with:<br>
* 5 rings with C, H and N atoms<br>
* 100 atoms with P1 symmetry <br>
</div>
* the rings are out of a plane<br>
</div>
* vacuum along y and z<br>
</div>
   I know that this is a very hard optimization problem, so I would like to kindly ask:<br>
<br>
</div>
1) Do you believe that MSEC3a will work better than MSR1a in such a case ?<br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
<div><br>
</div>
</span>
<div>I would switch to MSEC3a, use SLOW or reduce the GREED to 0.1 only if you see indications of problems. If the system is a decent insulator and the experimental positions are quite good you may have no problems. At the end I would switch back to MSR1a certainly
 for a system with OH as the positions of the hydrogens can be quite soft. It can take quite some time to get the rotations of the OH bond distance right as Wien2k uses cartesian coordinates not polars. Which of the 3 is best -- I am not sure.</div>
<span>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
<div dir="ltr">
<div>
<div>
<div>
<div>
<div>
<div><br>
</div>
2) Do you recommend using -it, -vec2pratt and -noHinv options with run_lapw ?<br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
<div><br>
</div>
</span>
<div>I use -it -noHinv -vec2pratt. Sometimes you need to do an occasional full diagonalization (touch .fulldiag) as the iterative method is less stable (it adds a little noise). In my personal version I have added back the old -itn option so this is done automatically
 every few steps. </div>
<span>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
<div dir="ltr">
<div>
<div>
<div>
<div>
<div><br>
</div>
3) Should I reduce TRUST to 0.5 (I am using LDA and experimental values for the initial atomic positions) ?<br>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
<div><br>
</div>
</span>
<div>Probably not. TRUST 0.5 would be if it is taking much too large steps which tends (in my experience) to occur more with soft electronic modes such as one has with d and f electrons. </div>
<span>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
<div dir="ltr">
<div>
<div>
<div>
<div><br>
</div>
4) Should I use SLOW in case.inm ?   <br>
</div>
</div>
</div>
</div>
</div>
</blockquote>
<div><br>
</div>
</span>
<div>See my answer to 1) </div>
<span>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
<div dir="ltr">
<div>
<div>
<div><br>
</div>
   Many thanks in advance.<br>
</div>
   All the best,<br>
</div>
                      Luis<br>
<br>
</div>
</div>
</blockquote>
</span></div>
<span class="m_-1646327723223284828m_-412137898862597576HOEnZb"><font color="#888888"><br>
<br clear="all">
<div><br>
</div>
-- <br>
<div class="m_-1646327723223284828m_-412137898862597576m_2547980773262265899gmail_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">C<wbr>orrosion 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="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.cfw.org_100-2Dpercent&d=CwMFaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=MNscdaGYE15HljeOPwtzDZH4S_JeFNBLtGtQWhVb-FY&s=kL1N7z6QnaXOzYazrQU6rH7eGamLhLbu0ncgZ2Hfz0E&e=" style="font-size:12.8px" target="_blank">www.cfw.org/100-percen<wbr>t</a></div>
<div>Co-Editor, Acta Cryst A</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</font></span></div>
</div>
<br></div></div></div><span class="">
______________________________<wbr>_________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.a<wbr>t</a><br>
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.theochem.tuwien.ac.at_mailman_listinfo_wien&d=CwMFaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=MNscdaGYE15HljeOPwtzDZH4S_JeFNBLtGtQWhVb-FY&s=UDSAKya4XEowVr79jb00mjlqKhzg84E5l-IbK3bng8c&e=" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac<wbr>.at/mailman/listinfo/wien</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=CwMFaQ&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=MNscdaGYE15HljeOPwtzDZH4S_JeFNBLtGtQWhVb-FY&s=iutaZBD8d3Ad8lUmlUTPnDkvBCtab7ay4LhoDtLRwqY&e=" rel="noreferrer" target="_blank">
http://www.mail-archive.com/wi<wbr>en@zeus.theochem.tuwien.ac.at/<wbr>index.html</a><br>
<br>
</span></blockquote>
</div>
<br>
</div>
</div>

</blockquote></div><br></div></div></div>
<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>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.<wbr>ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at:  <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/<wbr>wien@zeus.theochem.tuwien.ac.<wbr>at/index.html</a><br>
<br></blockquote></div><br></div>