<div dir="ltr"><div><div><div>A bit more spice to Teo's main course,<br></div>One of the reasons why your calculation will have better chances of succeeding when a TS candidate is provided <br></div>is because the band should display minimum lateral motion at the TS (assuming not very stiff springs). The tension at the TS should help stabilizing the band motion. As to the interpolation, it is a bit of a guessing game that cannot be won without using your physical/chemical intuition. Gradients (forces) in Cartesian coordinates tend to vary abruptly from one point to another point in space. Thus, using collective variables (when known) can help a lot to achieve convergence.<br> </div>Carlos<br><div> <br><div><div><div><div><br></div></div></div></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Aug 21, 2013 at 12:19 AM, Teodoro Laino <span dir="ltr"><<a href="mailto:teodor...@gmail.com" target="_blank">teodor...@gmail.com</a>></span> wrote:<br> <blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Jörg,<br> <div class="im"><br> On Aug 21, 2013, at 12:30 AM, Jörg Saßmannshausen <<a href="mailto:j.sassma...@ucl.ac.uk">j.sassma...@ucl.ac.uk</a>> wrote:<br> <br> > Dear all,<br> ><br> > I am in the process of setting up a BAND calculation with 16 frames.<br> > Each frame got 82 atoms and if I recall it correctly I need to have the<br> > cartesian coordinates in all frames in the same order. I assume that is still<br> > correct? As this is a rather large molecule, doing it manually is a bit<br> <br> </div>correct - still, this should not be a limitation at all : during a geometry optimisation or during an MD, the coordinates are never swapped..<br> of course if you build your frames with different atoms order.. well.. then you are asking for trouble… (or at least for more additional work :) )<br> <div class="im"><br> > cumbersome. I have tried to use GaussView for that but but it appears I need a<br> > larger screen to get all 16 molecules in one window. :-(<br> > Is there a better program to do that task?<br> <br> ?<br> <br> ><br> > Also, I want to investigate a dissociation reaction, i.e.:<br> ><br> > A-B ---> A B<br> ><br> > Due to the flexibility of the molecules A and B, the forward reaction might be<br> > different from the reverse reaction. Is there a way to influence that when I am<br> > doing the BAND?<br> <br> </div>A BAND is just a geometry optimisation is a wider space. What you put in, is what you will get out optimised… Rarely there is only one path connecting two points. This means that the success (in finding the minimum energy path) is entirely due to quality of the interpolating frames that you are feeding in.<br> <br> Having said that, if you provide different frames for the opposite reaction you will get probably a different path - unless the landscape is relatively simple.<br> <div class="im"><br> > I got the frames generated by pulling (pushing) A and B apart (together) so I<br> > know the geometry is optimised here. I am doing steps of 0.3 A so that should<br> > be ok as well (I hope). So, if I am using the optimised geometries for the<br> > forward and reverse reaction I should be able to get what I am expecting or is<br> > there something else which I need to bear in mind?<br> <br> </div>Just pulling apart the molecule and optimising them is quite risky:<br> In the rare event that the reaction is only driven by this distance you will have automatically an optimised MEP. But, if the reaction is driven by much more complex reaction coordinates, then you are severely biasing the process.<br> The best approach is to provide A, B and something similar to your TS (or any other crucial point, that you know is characterising your reaction). Let all other frames being interpolated and let the BAND optimisation relax the MEP.<br> <br> Again, the final quality will be very much dependent on how good you'll be able to provide your TS (or these crucial intermediate points).<br> If you do not have any idea about the TS, you can still try with a linear interpolation (or something similar to your procedure) and see what happens, analysing critically the results: most of the time it is easy to understand if you are just throwing ropes between minima, or you really got a MEP.<br> <br> <br> Teo<br> <div class="HOEnZb"><div class="h5"><br> ><br> > I hope that makes sense to you.<br> ><br> > All the best from London!<br> ><br> > Jörg<br> > --<br> > *************************************************************<br> > Jörg Saßmannshausen<br> > University College London<br> > Department of Chemistry<br> > Gordon Street<br> > London<br> > WC1H 0AJ<br> ><br> > email: <a href="mailto:j.sassma...@ucl.ac.uk">j.sassma...@ucl.ac.uk</a><br> > web: <a href="http://sassy.formativ.net" target="_blank">http://sassy.formativ.net</a><br> ><br> > Please avoid sending me Word or PowerPoint attachments.<br> > See <a href="http://www.gnu.org/philosophy/no-word-attachments.html" target="_blank">http://www.gnu.org/philosophy/no-word-attachments.html</a><br> ><br> > --<br> > You received this message because you are subscribed to the Google Groups "cp2k" group.<br> > To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k%2Bun...@googlegroups.com">cp2k+uns...@googlegroups.com</a>.<br> > To post to this group, send email to <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a>.<br> > Visit this group at <a href="http://groups.google.com/group/cp2k" target="_blank">http://groups.google.com/group/cp2k</a>.<br> > For more options, visit <a href="https://groups.google.com/groups/opt_out" target="_blank">https://groups.google.com/groups/opt_out</a>.<br> <br> --<br> You received this message because you are subscribed to the Google Groups "cp2k" group.<br> To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k%2Bun...@googlegroups.com">cp2k+uns...@googlegroups.com</a>.<br> To post to this group, send email to <a href="mailto:cp...@googlegroups.com">cp...@googlegroups.com</a>.<br> Visit this group at <a href="http://groups.google.com/group/cp2k" target="_blank">http://groups.google.com/group/cp2k</a>.<br> For more options, visit <a href="https://groups.google.com/groups/opt_out" target="_blank">https://groups.google.com/groups/opt_out</a>.<br> </div></div></blockquote></div><br></div>