<div dir="ltr">Inlined is my intuition, which does not have to be completely right.<div class="gmail_extra"><br><div class="gmail_quote">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><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><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><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><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><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><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><div>See my answer to 1) </div><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></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></div>