<p dir="ltr">Your question is based upon a common misunderstanding of minimizing forces. You are thinking that the time to optimize depends upon the number of positions (3*Number of Pt atoms that move) which is wrong. With PORT it scales as the number of clusters of phonon frequencies and weakly with the width of the clusters; with MSR1a it scales as the number of clusters of eigenvalues of the joint electron-phonon problem and their width.</p>
<p dir="ltr">You will gain far more by generating a proper symmetric model with P-3m1 (maybe P-31m) symmetry rather than what you appear to have, a P1 structure. This will reduce the number of uniques phonon etc modes by something like 6-12, and will be much faster and converge better.</p>
<p dir="ltr">---<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>
<a href="http://www.numis.northwestern.edu">http://www.numis.northwestern.edu</a><br>
Corrosion in 4D <a href="http://MURI4D.numis.northwestern.edu">http://MURI4D.numis.northwestern.edu</a><br>
Partner of the CFW 100% gender equity project, <a href="http://www.cfw.org/100-percent">www.cfw.org/100-percent</a><br>
Co-Editor, Acta Cryst A<br>
</p>
<div class="gmail_quote">On Jul 14, 2016 06:25, "Saleem Ayaz" <<a href="mailto:saleemayaz_hu@yahoo.com">saleemayaz_hu@yahoo.com</a>> wrote:<br type="attribution"><blockquote class="quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
<div style="color:#000;background-color:#fff;font-family:bookman old style,new york,times,serif;font-size:16px">
<div>Dear users</div>
<div dir="ltr">I want to minimize the forces for Pt111 surface. I want to restrict the structure by minimizing the forces for first two layers and keeping the other layers bulk-like. Please some one can help me how to change
case.inM file. Should I make some numbers equal to zero?</div>
<div dir="ltr">Looking forward</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">regards</div>
<div dir="ltr">Saleem</div>
<div dir="ltr"><br>
</div>
<div dir="ltr">PORT 2.00 0.35 # PORT/NEWT; tolf, Initial Trust Radius</div>
<div dir="ltr">1.0 1.0 1.0 1.0 #Atom 1 Generated by pairhess</div>
<div dir="ltr">1.0 1.0 1.0 1.0 #Atom 2 Generated by pairhess</div>
<div dir="ltr">1.0 1.0 1.0 1.0 #Atom 3 Generated by pairhess</div>
<div dir="ltr">1.0 1.0 1.0 1.0 #Atom 4 Generated by pairhess</div>
<div dir="ltr">1.0 1.0 1.0 1.0 #Atom 5 Generated by pairhess</div>
<div dir="ltr">1.0 1.0 1.0 1.0 #Atom 6 Generated by pairhess</div>
<div dir="ltr">1.0 1.0 1.0 1.0 #Atom 7 Generated by pairhess</div>
<div dir="ltr">1.0 1.0 1.0 1.0 #Atom 8 Generated by pairhess</div>
<div dir="ltr">1.0 1.0 1.0 1.0 #Atom 9 Generated by pairhess</div>
<div dir="ltr">1.0 1.0 1.0 1.0 #Atom 10 Generated by pairhess</div>
<div dir="ltr">1.0 1.0 1.0 1.0 #Atom 11 Generated by pairhess</div>
<div dir="ltr">1.0 1.0 1.0 1.0 #Atom 12 Generated by pairhess</div>
<div dir="ltr">1.0 1.0 1.0 1.0 #Atom 13 Generated by pairhess</div>
<div dir="ltr">1.0 1.0 1.0 1.0 #Atom 14 Generated by pairhess</div>
<div dir="ltr">1.0 1.0 1.0 1.0 #Atom 15 Generated by pairhess</div>
<div dir="ltr">1.0 1.0 1.0 1.0 #Atom 16 Generated by pairhess</div>
<div dir="ltr">1.0 1.0 1.0 1.0 #Atom 17 Generated by pairhess</div>
<div dir="ltr">1.0 1.0 1.0 1.0 #Atom 18 Generated by pairhess</div>
<div dir="ltr">1.0 1.0 1.0 1.0 #Atom 19 Generated by pairhess</div>
<div dir="ltr">1.0 1.0 1.0 1.0 #Atom 20 Generated by pairhess</div>
<div dir="ltr">1.0 1.0 1.0 1.0 #Atom 21 Generated by pairhess</div>
<div dir="ltr">1.0 1.0 1.0 1.0 #Atom 22 Generated by pairhess</div>
<div dir="ltr">1.0 1.0 1.0 1.0 #Atom 23 Generated by pairhess</div>
<div dir="ltr">1.0 1.0 1.0 1.0 #Atom 24 Generated by pairhess</div>
<div dir="ltr">1.0 1.0 1.0 1.0 #Atom 25 Generated by pairhess</div>
<div dir="ltr">1.0 1.0 1.0 1.0 #Atom 26 Generated by pairhess</div>
<div dir="ltr">1.0 1.0 1.0 1.0 #Atom 27 Generated by pairhess</div>
</div>
</div>
</blockquote></div>