<div dir="ltr">Dear Juerg,<br><br><br>thanks for your answer. Whilst I am using an SCF convergence of 10e-6 as you suggest, I am not using the conjugate gradient OT minimizer but DIIS instead. Perhaps something to try out, I will give it a go. What I was wondering about is, if there is a possibility to control these sort of parameters for individual replicas.<br>The electronic structure will change quite a bit as one traverses along the profile between starting and last point, hence I wouldn't expect too much of a problem with the first point, but more towards the transition state. Anyways, I guess it is more reasonable to control the SCF globally for all points.<br><br>Another question I have in this context, is there a rule of thumb for how many intermediate structures for NEB should be provided? At the moment I am using a starting structure, a guess for the TS plus a final structure. With 8 replicas defined, the algorithm should interpolate between these. I was successful to localise<br>one transition state with this already, however, for other transition states things are more complicated. Also, in cases I observe rather strange energy profiles, in which the first (sometimes also the last) point lies higher in energy then the remaining points (see attachment). This is somewhat surprising, since the first point corresponds to an optimised structure (I am assuming that the first replica correlates with my first starting structure). Inspection of the movie.xyz generated from the individual BAND files reassures me that I am looking at a reasonable trajectory. Is is worth waiting for this calculation to converge to some stationary point, or should I come up with a better setup in terms of structures?    <br><br>Thanks for the advice<br><br><br>Tobias<br>    <br>   <br><br>  <br><br><br><br>On Wednesday, July 29, 2015 at 8:33:36 AM UTC+1, jgh wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;">Hi Tobias
<br>
<br>I would start using a very conservative setting for the SCF
<br>optimizer. For example
<br>
<br>    &SCF
<br>      MAX_SCF 20
<br>      SCF_GUESS ATOMIC
<br>      EPS_SCF  1.e-06
<br>      &OT ON
<br>        MINIMIZER CG
<br>        ENERGY_GAP 0.002
<br>        LINESEARCH 3PNT
<br>        PRECONDITIONER FULL_ALL
<br>      &END OT
<br>      &OUTER_SCF
<br>        MAX_SCF 20
<br>        EPS_SCF  1.e-06
<br>      &END OUTER_SCF
<br>    &END SCF
<br>
<br>But maybe you already do this?
<br>
<br>Juerg
<br>
<br>--------------------------------------------------------------
<br>Juerg Hutter                         Phone : ++41 44 635 4491
<br>Institut für Chemie C                FAX   : ++41 44 635 6838
<br>Universität Zürich                   E-mail: <a href="javascript:" target="_blank" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">hut...@chem.uzh.ch</a>
<br>Winterthurerstrasse 190
<br>CH-8057 Zürich, Switzerland
<br>---------------------------------------------------------------
<br>
<br>-----<a href="javascript:" target="_blank" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a> wrote: -----To: cp2k <<a href="javascript:" target="_blank" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>>
<br>From: Tobias Kraemer <<a href="javascript:" target="_blank" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">161br...@gmail.com</a>>
<br>Sent by: <a href="javascript:" target="_blank" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>
<br>Date: 07/28/2015 12:44PM
<br>Subject: [CP2K:6828] SCF convergence in NEB calculation
<br>
<br>Hi all,
<br>
<br>
<br>I am running CI-NEB calculations in order to locate transition states along a reaction profile involving organometallic species in the solid state. More precisely, this is a hydrogenation reaction of double bonds, i.e. addition of H2 to an alkene. This should be a straightforward task. Anyway, for some cases it seems the SCF encounters problems to converge, or converges to the wrong electronic state. In particular, those replicas are affected which lead to the  transition
<br> state. I can see that the electronic structure may be a bit tricky, 
<br>since in a TS notable rearrangement of the electrons takes place. It goes without saying, that having a wrong electronic structure in the TS leads to unreasonably high barriers. I should add that these type of calculations haven't imposed any problems so far using other codes. What causes some additional puzzlement is the fact that for the same TS search slightly different input  geometries lead to very different behaviour of the SCF. In one case there is reasonable convergence to a decent electronic structure, in the other case the SCF goes completely havoc.  
<br>
<br>So my question really is, how can I control the SCF convergence for a particular replica? Can this be done in the main input file?
<br>
<br>your comments are appreciated
<br>
<br>
<br>Tobi
<br>
<br>
<br>
<br>
<br>
<br>-- 
<br>
<br>You received this message because you are subscribed to the Google Groups "cp2k" group.
<br>
<br>To unsubscribe from this group and stop receiving emails from it, send an email to <a href="javascript:" target="_blank" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp2k+...@googlegroups.com</a>.
<br>
<br>To post to this group, send email to <a href="javascript:" target="_blank" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>.
<br>
<br>Visit this group at <a href="http://groups.google.com/group/cp2k" target="_blank" rel="nofollow" onmousedown="this.href='http://groups.google.com/group/cp2k';return true;" onclick="this.href='http://groups.google.com/group/cp2k';return true;">http://groups.google.com/group/cp2k</a>.
<br>
<br>For more options, visit <a href="https://groups.google.com/d/optout" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;">https://groups.google.com/d/optout</a>.
<br>
<br></blockquote></div>