<div dir="ltr">It works...but (sorry Peter) it can be bad science. There are special cases where it is appropriate, but it can also be very inappropriate.<div><br></div><div>There is a long history of people fixing some atomic positions to bulk positions when performing surface calculations. However, the strain field of any surface decays exponentially in the selvedge layer, and removing this can be very wrong. For instance, there are many known cases where the surface induces a ferroelectric distortion which fixed positions will prevent from happening. The rigorous method is to check that the atomic positions in the center of the slab are "close enough" to those in the bulk. If they are not then either the slab needs to be thicker or there has been a surface induced phase transition.</div><div><br></div><div>And, fixing positions for a surface <i style="font-weight:bold">in general does not improve the convergence</i>. If you are lucky it will, but not in general. The reason (as I have mentioned before) is that the convergence scales with the number and width of the eigenvalue clusters, not the number of variables. If the convergence of Wien2k and other DFT codes scaled with the number of density variables, the methods would be useless! DFT is a useful method because it will typically converge in 10-100 scf iterations because there are only a few eigenvector clusters of the dielectric response (in density only calculations) and the force matrix (with PORT). MSR1a is typically efficient because the number and width of the clusters of the general dielectric response with both density and electronic terms included is smaller than the product of the clusters for the electronic and atomic components.</div><div><br></div><div>Sorry for being so pedantic...but rigor matters.</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, May 19, 2018 at 7:56 AM, Peter Blaha <span dir="ltr"><<a href="mailto:pblaha@theochem.tuwien.ac.at" target="_blank">pblaha@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"><span class="">> - How can I lock some atoms in position during structure optimization?<br>
<br>
</span>In case.inM  put the 1.0 values of x,y and/or z of the corresponding <br>
atom (line) to 0.0.<br>
<br>
It works well. We do it all the time....<br>
<br>
> <br>
> Best,<br>
> Lukasz<br>
> <br>
<br>
-- <br>
------------------------------<wbr>------------------------------<wbr>--------------<br>
Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna<br>
Phone: +43-1-58801-165300             FAX: +43-1-58801-165982<br>
Email: <a href="mailto:blaha@theochem.tuwien.ac.at">blaha@theochem.tuwien.ac.at</a>    WIEN2k: <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.wien2k.at&d=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=HGBSlLfU5SbOh_JfodIoRVa6I69Rsj216R2i2idPWPA&s=Z18mJHtLP0OuMpkNbZyk6BppzH7a7iqaKjSmbAxk0Cw&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.<wbr>com/v2/url?u=http-3A__www.<wbr>wien2k.at&d=DwICAg&c=<wbr>yHlS04HhBraes5BQ9ueu5zKhE7rtNX<wbr>t_d012z2PA6ws&r=U_<wbr>T4PL6jwANfAy4rnxTj8IUxm818jnvq<wbr>KFdqWLwmqg0&m=HGBSlLfU5SbOh_<wbr>JfodIoRVa6I69Rsj216R2i2idPWPA&<wbr>s=<wbr>Z18mJHtLP0OuMpkNbZyk6BppzH7a7i<wbr>qaKjSmbAxk0Cw&e=</a><br>
WWW: <br>
<a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.imc.tuwien.ac.at_tc-5Fblaha-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-2D&d=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=HGBSlLfU5SbOh_JfodIoRVa6I69Rsj216R2i2idPWPA&s=GxPbnc8632eZdJgUG2LT9xuLo25j8wmiwf7mi35vykc&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.<wbr>com/v2/url?u=http-3A__www.imc.<wbr>tuwien.ac.at_tc-5Fblaha-2D-2D-<wbr>2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-<wbr>2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-<wbr>2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-<wbr>2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-<wbr>2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-<wbr>2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-<wbr>2D-2D-2D-2D-2D-2D-2D-2D-2D-2D-<wbr>2D&d=DwICAg&c=<wbr>yHlS04HhBraes5BQ9ueu5zKhE7rtNX<wbr>t_d012z2PA6ws&r=U_<wbr>T4PL6jwANfAy4rnxTj8IUxm818jnvq<wbr>KFdqWLwmqg0&m=HGBSlLfU5SbOh_<wbr>JfodIoRVa6I69Rsj216R2i2idPWPA&<wbr>s=<wbr>GxPbnc8632eZdJgUG2LT9xuLo25j8w<wbr>miwf7mi35vykc&e=</a> <br>
<span class=""><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>
</span><a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__zeus.theochem.tuwien.ac.at_mailman_listinfo_wien&d=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=HGBSlLfU5SbOh_JfodIoRVa6I69Rsj216R2i2idPWPA&s=ax2aySwdZVqULGgL-LOIT7EHit4gKJxeBu2ErcMmvEY&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=DwICAg&c=<wbr>yHlS04HhBraes5BQ9ueu5zKhE7rtNX<wbr>t_d012z2PA6ws&r=U_<wbr>T4PL6jwANfAy4rnxTj8IUxm818jnvq<wbr>KFdqWLwmqg0&m=HGBSlLfU5SbOh_<wbr>JfodIoRVa6I69Rsj216R2i2idPWPA&<wbr>s=ax2aySwdZVqULGgL-<wbr>LOIT7EHit4gKJxeBu2ErcMmvEY&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=DwICAg&c=yHlS04HhBraes5BQ9ueu5zKhE7rtNXt_d012z2PA6ws&r=U_T4PL6jwANfAy4rnxTj8IUxm818jnvqKFdqWLwmqg0&m=HGBSlLfU5SbOh_JfodIoRVa6I69Rsj216R2i2idPWPA&s=ke6SIVUCmWlbQaSt9_oxTYR55RBkTN_pIRsUznD79B8&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=DwICAg&c=<wbr>yHlS04HhBraes5BQ9ueu5zKhE7rtNX<wbr>t_d012z2PA6ws&r=U_<wbr>T4PL6jwANfAy4rnxTj8IUxm818jnvq<wbr>KFdqWLwmqg0&m=HGBSlLfU5SbOh_<wbr>JfodIoRVa6I69Rsj216R2i2idPWPA&<wbr>s=ke6SIVUCmWlbQaSt9_<wbr>oxTYR55RBkTN_pIRsUznD79B8&e=</a><br>
</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>