<div dir="ltr">Hi Tobias,<div><br></div><div>One thing you need to keep in mind is that in CP2K, the vibrational modes are calculated numerically, i.e. they will be sensitive to the DX value you specified in your &VIBRATIONAL_ANALYSIS section. You will also need to use extremely tight energy and force convergence criteria to obtain a very good geometry for the VIBRATIONAL_ANALYSIS calculation.</div><div><br></div><div>SL</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On 10 September 2015 at 12:46, Tobias Kraemer <span dir="ltr"><<a href="mailto:161brun...@gmail.com" target="_blank">161brun...@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hello,<br><br><br>I have tried out your suggestions regarding the dimer method, and in principle I have now successfully run this calculation. Unfortunately I didn't get the desired result, i.e. I didn't eliminate the spurious imagianary mode in the ts structure. What surprises me is the fact that there are hardly any noticable changes in the geometry. Changes affect the 3. and 4. decimal place of bond distances and angles (negliglible). What I expected to see are somewhat more noticable movements in those parts of the molecule which correspond to the small imaginary mode and should refine to a minimum (say small rotations of certain groups or other soft modes). This hasn't happened apparently, and hence the new structure is almost identical to the input structure, and hence - unsurprisingly - yields the same undesired imaginary modes again. I want to point out that this spurious mode is around 25 cm-1, a very soft mode, but it must be possible to get rid of it (and in fact, the constrained optimisation which I have used in another but related example did exaclty this). Could you advice me here in terms of parameters I should change. I have played around with the angle tolerance parameter, what about DR? Shoud I use tighter geometry convergence? At the moment I am more worried about the fatc that there is hardly any change in the geometry during the dimer optimisation. What protocols do you normally employ when you optimise to a saddle point of higher order? How do you refine the TS?<br><br>Thanks for your help (inputs plut other files attached)<span class="HOEnZb"><font color="#888888"><br><br>Tobias</font></span><span class=""><br><br> <br><br><br><br> <br><br>On Tuesday, September 1, 2015 at 5:11:13 PM UTC+1, Tobias Kraemer wrote:</span><div><div class="h5"><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi all,<br><br><br>I am trying to re-optimise a transition state structure which has a spurious (small) imaginary mode on top of the large imaginary mode which corresponds to the transition state. I've also tested doing a constrained optimisation fixing selected atoms involved in the transition state, which also worked in this particular case. Howeve, I would also like to try our the dimer method for the purpose of re-optimising a transition state structure. To this end, I have composed an input, however, the job crashes and I keep getting an error message about certain line Search type not being implemented (looks quite obvious from the message). However, I can't seem to find what I need to change to make it work. Your advice would be much appreciated. The motion section can be found below, full input file also uploaded. <br><br>Any ideas on how to correctly set up a calculation using the dimer method are welcome<br><br><br>Thanks for your help<br><br>Tobias<br><br><br> ************************************************************************<br> *** 16:48:48 ERRORL2 in cg_utils:cg_linmin processor 0 :: err=-1000 ***<br> *** UNIMPLEMENTED, Line Search type not yet implemented in CG for TS ***<br> *** search. ***<br> ************************************************************************<br><br><br> ===== Routine Calling Stack =====<br><br> 10 cg_linmin<br> 9 cp_cg_main<br> 8 geoopt_cg<br> 7 cp_rot_opt<br> 6 cp_eval_at_ts<br> 5 cp_eval_at<br> 4 cp_cg_main<br> 3 geoopt_cg<br> 2 cp_geo_opt<br> 1 CP2K<br><br><br><br><br>&MOTION<br> &GEO_OPT<br> TYPE TRANSITION_STATE<br> &TRANSITION_STATE<br> METHOD DIMER<br> &DIMER<br> &ROT_OPT<br> OPTIMIZER CG<br> &END ROT_OPT<br> &END DIMER<br> &END TRANSITION_STATE<br> OPTIMIZER CG<br> MAX_FORCE 8.0E-4<br> MAX_ITER 500<br> &END GEO_OPT<br> &PRINT<br> &TRAJECTORY<br> &EACH<br> GEO_OPT 1<br> &END EACH<br> &END TRAJECTORY<br> &END<br>&END MOTION<br><br></div></blockquote></div></div></div><div class="HOEnZb"><div class="h5">
<p></p>
-- <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+uns...@googlegroups.com" target="_blank">cp2k+uns...@googlegroups.com</a>.<br>
To post to this group, send email to <a href="mailto:cp...@googlegroups.com" target="_blank">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/d/optout" target="_blank">https://groups.google.com/d/optout</a>.<br>
</div></div></blockquote></div><br></div>