<div dir="ltr">Hi Matthias, <div><br></div><div>Thanks very much! I miss understand the point you mention "You should rather check if fixing atoms is really needed as a boundary for a realistic modelling of the reaction you are interested in." Sorry about this. </div><div><br></div><div>Here I would like to reported the performance for my current calculation (~600 atoms). Fix atoms does not improve the speed of BAND calculation. The observation is based on two ionic steps.  </div><div><br></div><div>Best regards, </div><div>J. Ye</div><div><br></div><div><br>On Tuesday, November 8, 2016 at 10:40:24 AM UTC-6, Matthias Krack wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hi<br><br>the default values might be tight, but these defaults should be considered like all other default values only as a first guess or proposal which has to be adapted to the actual application.<br><br>Best regards<br><br>Matthias<br><br>On Tuesday, 8 November 2016 17:12:41 UTC+1, J. Ye  wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>I have another question. The default value for <a href="https://manual.cp2k.org/cp2k-4_1-branch/CP2K_INPUT/MOTION/BAND/CONVERGENCE_CONTROL.html" style="font-family:monospace;text-transform:uppercase" rel="nofollow" target="_blank" onmousedown="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fmanual.cp2k.org%2Fcp2k-4_1-branch%2FCP2K_INPUT%2FMOTION%2FBAND%2FCONVERGENCE_CONTROL.html\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEppGjz2UO84WkxF8os5-UiNw-EAg';return true;" onclick="this.href='https://www.google.com/url?q\x3dhttps%3A%2F%2Fmanual.cp2k.org%2Fcp2k-4_1-branch%2FCP2K_INPUT%2FMOTION%2FBAND%2FCONVERGENCE_CONTROL.html\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNEppGjz2UO84WkxF8os5-UiNw-EAg';return true;">CONVERGENCE_CONTROL</a> for BAND optimization is very difficult to converge. I was curious is why the criterias for BAND optimization is so high? What is the reasonable criteria?   </div><div><div>    &CONVERGENCE_CONTROL</div><div>      MAX_DR 0.0002</div><div>      MAX_FORCE 0.00045</div><div>      RMS_DR 0.0001</div><div>      RMS_FORCE 0.0003</div><div>    &END CONVERGENCE_CONTROL</div></div><div><br></div><div><br></div><div><br>On Tuesday, November 8, 2016 at 9:32:49 AM UTC-6, J. Ye wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Matthias, <div><br></div><div>I am working on a catalyst material that has metal nanoparticles supported on a porous material.  This porous material which is a big system that has about 600 atoms, however the catalytic reaction process only on the metal clusters. We want to study the effect of the pore size on the catalytic reaction. So I want to perform a periodic calculations and compare with cluster model. If the barriers are similar, then I could perform for cluster model instead. </div><div><br></div><div>Your suggestion is peroform QM/MM? </div><div><br></div><div>What's the difference between DIIS and MD? Why DIIS is more efficient? Do you have any books or papers recommend for us to read? </div><div><br></div><div>Thanks very much!</div><div><br></div><div>Best regards,</div><div> </div><div>J. Ye<br><br>On Tuesday, November 8, 2016 at 2:12:42 AM UTC-6, Matthias Krack wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi<br><br>do you have any special reason to use the OPT_TYPE MD? Otherwise I suggest to try the default OPT_TYPE DIIS which should be much more efficient. Moreover, if you want to fix some atoms only because of the computational cost caused by the large size of your system then there will be most likely not much benefit doing so, because the forces for all atoms are calculated anyway and then set to zero for the fixed atoms afterwards. You should rather check if fixing atoms is really needed as a boundary for a realistic modelling of the reaction you are interested in.<br><br>Best regards<br><br>Matthias<br><br>On Tuesday, 8 November 2016 07:01:18 UTC+1, J. Ye  wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"> Hi Fabio, <div><br></div><div>I indeed included both "ALIGN_FRAMES FALSE" and "ROTATE_FRAMES FALSE" in my input(Please see the attachment). But I still found the atoms moves. The methods we use are different. I don't know how does this affect the fixed atoms. Do you have any idea? </div><div><br></div><div>I have tested your input. Atom 1 and 3 indeed do not move. I also tried to switch TI-NEB to CI-NEB, the atoms are not move. So this should not be affected by BAND_TYPE. </div><div><br></div><div>Can any expert of CP2K give us some hints on this problem? </div><div><br></div><div>Thanks very much!</div><div>Best,</div><div><br></div><div>J. Ye</div><div><br>On Monday, November 7, 2016 at 2:02:05 PM UTC-6, <a>elb...@gmail.com</a> wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hello J. Ye!<br><br>Yes, I have checked. If you use only the constraint flag, indeed the fixed atoms move for the other images. But if you include both "ALIGN_FRAMES FALSE" and "ROTATE_FRAMES FALSE" flags in the BAND section, that does not happen anymore, at least in the calculations I have performed so far. Check for example the NEB input attached to this e-mail. It is very simple and it runs very fast (it has no meaning, it is just an example), but you can check that atoms 1 and 3, that are asked to be fixed, do not move for all NEB images. Let me know if you get the same behavior as me!<br><br>All the best!<br></div></blockquote></div></div></blockquote></div></blockquote></div></div></blockquote></div></div></blockquote></div></blockquote></div></div>