Differences

This shows you the differences between two versions of the page.

--- exercises:2018_uzh_cmest:path_optimization_neb [2018/10/23 15:52] – [Vibrational analysis] abussy
+++ exercises:2018_uzh_cmest:path_optimization_neb [2020/08/21 10:15] (current) – external edit 127.0.0.1
@@ Line 5: / Line 5: @@
 In this exercise, we will therefore perform Nudged Elastic Band (NEB) calculations using the same molecule as before and investigate the energy path between the two geometries.
-Following are three geometry files you should put/create in a new exercise directory:
+Following are four geometry files you should put/create in a new exercise directory:
 <code - ethane_1_opt.xyz>
@@ Line 300: / Line 300: @@
 </code>
-Once this run completes, you should find a file ''ethane_TS_va-VIBRATIONS-1.mol''.
+Once this run completes, you can grep for the modes frequency with:
-Now we are going to use the application //molden// (which you can load using ''module load molden'') to visualize the vibrational modes:
 <code>
-$ molden ethane_TS_va-VIBRATIONS-1.mol
+$ grep "VIB|Frequency" ethane_TS_va.out
 </code>
-<note important>Molden is currently not available on tcopt2, it should be fixed soon. Sorry for the inconvenience.</note>
+The presence of a  negative (imaginary) mode means that it is actually a transition state (and not stable). Are there any of that kind for this structure ? If so, what are the atomic displacements associated with it ?
-Click the //Norm. Mode// checkbox in the //Molden Control// window to list all the modes. What is the lowest frequency you get? By clicking on it you can visualize it.
-The presence of a  negative (imaginary) mode means that it is actually a transition state (and not stable).
-Now repeat the same steps presented here for the bead with the lowest energy. What is now the first frequency you get in the list? Is this geometry stable?
-Please note: while you should get only 18 different frequencies you get 21 instead. That means that 3 frequencies are global rotations instead of modes in the molecule and should be ignored when looking for negative frequencies to identify whether a conformer is stable or not.
+Note that cp2k also produces a ''ethane_TS_va-VIBRATIONS-1.mol'' file. It is intended to be read by the ''Molden'' program, which makes the visualization of the modes possible.