<p dir="ltr">Please look in the UG at clmextrapol, and also what is done in the optimize.job script created by x optimize. The clmextrapol utility is designed for this type of issue.</p>
<p dir="ltr">___________________________<br>
Professor Laurence Marks<br>
Department of Materials Science and Engineering<br>
Northwestern University<br>
<a href="http://www.numis.northwestern.edu">www.numis.northwestern.edu</a><br>
<a href="http://MURI4D.numis.northwestern.edu">MURI4D.numis.northwestern.edu</a><br>
Co-Editor, Acta Cryst A<br>
"Research is to see what everybody else has seen, and to think what nobody else has thought"<br>
Albert Szent-Gyorgi</p>
<div class="gmail_quote">On Apr 23, 2015 4:26 AM, "Madesis Ioannis(John)" <<a href="mailto:imadesis@physics.uoc.gr">imadesis@physics.uoc.gr</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">First of all, I want to apologize for sending yesterdays answer to a new<br>
topic.As it probably needs to be erased, you will find it as<br>
"msg12195.html". Hopefully, I am not doing the same mistake now.<br>
Continuing, I will copy that text here for proper continuity of the<br>
discussion.<br>
<br>
This is why I specified that I work with phonopy. Phonopy, (if you don't<br>
already know) requires all generated structures to be of P symmetry. To<br>
be honest, this is not always possible since I cannot override Wien2K<br>
algorith, which detects possible symmetries (initialization runs in<br>
batch mode, inside a home-made script). For phonopy, all atoms must be<br>
inequivalent, for proper results. Since I run FIRST all the<br>
initializations, and keep a log file, I am in position to tell if the<br>
spacegroup of all structures is P1, which again, is a primary target for<br>
a good phonopy calculation.<br>
<br>
So, your suggestion to decrease the duration of the convergence for the<br>
n+1 cases is to<br>
<br>
cp case_n/case_n.clmsum case_n+1/case_n+1.clmsum for 1<n<=Number of<br>
displacements<br>
<br>
after the initialization of the cases and the run of case_1.<br>
Are there any other files that could be transfered, therefor, enhance<br>
this idea even more?<br>
<br>
Today, I came to the realization, that the idea worked! I am currently<br>
running a calculation with 78 displacements.<br>
At #15 and greater, I replaced the recommended file, with the one<br>
proposed. The results so far, are the following:<br>
The first 14 structures required required 20-21 cycles for -cc 0.001 -ec<br>
0.001 -fc 1 convergence. So, at this point I regard them as the "rule"<br>
for the duration.<br>
Of following 15 structures that have run since the replacement, 11<br>
required only 12 cycles, with the same convergence criteria. Checking<br>
the time of the last modification, I can see that the 12 cycles require<br>
54 minutes, whereas the 20 cycles require 86. If one combines the<br>
"success rate" (so far) with the improvement in time for convergence<br>
(For the particular sample with 52 atoms), the result is 25%<br>
improvement, that could probably go up to 37%, depending on the<br>
following "success rate".<br>
<br>
--<br>
Ioannis Madesis (Μαδέσης Ιωάννης)<br>
PhD Student<br>
Atomic and Molecular Physics<br>
Department of Physics<br>
University of Crete<br>
(0030)-210-6503598<br>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" target="_blank">http://zeus.theochem.tuwien.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" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</blockquote></div>