Statistics
| Revision:

root / examples / Test / Zmat / README @ 5

History | View | Annotate | Download (1.8 kB)

1
As HCN can almost be considered as a 2D case because the CN bond does not
2
 change much during the HCN -> CNH path, one can visualize the path 
3
by projecting it on a 2D PES.
4
This is done by analyzing a path.
5
So, first, you have to launch run_Paht.
6
Then you can analyse it either using:
7
**)  AnaPath (as usual)
8
../../../utils/AnaPath HCN_zmat.out 50 Path 12 
9
**) Dep_tscan_p 
10
./Dep_tscan_p HCN_zmat.out 50 Path PES
11
This will create a Path_p.gplot file that can be ran to visualizd the
12
 evolution of the Path projected on the 2D PES (in green). The actual
13
 IRC is also plotted in red.
14
To go to the next iteration, you have to press a key.
15

    
16
**) Dep_tscan_m
17
./Dep_tscan_m HCN_zmat.out 50 Path PES
18
This will create a Path_m.gplot file that can be ran to visualizd the
19
 evolution of the Path projected on the 2D PES (in green). The actual
20
 IRC is also plotted in red, as with Dep_tscan_p.
21
What is new is that you have 2 windows: the second one shows the energy
22
 along the path.
23
To go to the next iteration, you have to press a key.
24

    
25
**) Dep_tscan_movie
26
./Dep_tscan_movie HCN_zmat.out 50 Path PES
27
This will create a Path_movie.gplot file that can be ran to visualizd the
28
 evolution of the Path projected on the 2D PES (in green). The actual
29
 IRC is also plotted in red, as with Dep_tscan_p.
30
Here the evolution is shown as a movie, there is no need to press a key to go
31
to the next iteration.
32

    
33
** Remarks:
34
For all Dep_tscan_XX scripts, 
35
The "PES" argument can be omitted if you do not want to plot the 2D PES contours in green.
36
If you have asked (using Path.valid, see manual) to print the tangents, then you can add
37
 "T" as a last  argument, this will print the path and its tangents.
38
All three scripts use the files HCN_irc.dat to plot the actual IRC, as well as HCN_PES_2D.dat
39
 to plot the 2D PES.