Statistics
| Revision:

root / examples / Test / Zmat / README @ 5

History | View | Annotate | Download (1.8 kB)

1 1 equemene
As HCN can almost be considered as a 2D case because the CN bond does not
2 1 equemene
 change much during the HCN -> CNH path, one can visualize the path
3 1 equemene
by projecting it on a 2D PES.
4 1 equemene
This is done by analyzing a path.
5 1 equemene
So, first, you have to launch run_Paht.
6 1 equemene
Then you can analyse it either using:
7 1 equemene
**)  AnaPath (as usual)
8 1 equemene
../../../utils/AnaPath HCN_zmat.out 50 Path 12
9 1 equemene
**) Dep_tscan_p
10 1 equemene
./Dep_tscan_p HCN_zmat.out 50 Path PES
11 1 equemene
This will create a Path_p.gplot file that can be ran to visualizd the
12 1 equemene
 evolution of the Path projected on the 2D PES (in green). The actual
13 1 equemene
 IRC is also plotted in red.
14 1 equemene
To go to the next iteration, you have to press a key.
15 1 equemene
16 1 equemene
**) Dep_tscan_m
17 1 equemene
./Dep_tscan_m HCN_zmat.out 50 Path PES
18 1 equemene
This will create a Path_m.gplot file that can be ran to visualizd the
19 1 equemene
 evolution of the Path projected on the 2D PES (in green). The actual
20 1 equemene
 IRC is also plotted in red, as with Dep_tscan_p.
21 1 equemene
What is new is that you have 2 windows: the second one shows the energy
22 1 equemene
 along the path.
23 1 equemene
To go to the next iteration, you have to press a key.
24 1 equemene
25 1 equemene
**) Dep_tscan_movie
26 1 equemene
./Dep_tscan_movie HCN_zmat.out 50 Path PES
27 1 equemene
This will create a Path_movie.gplot file that can be ran to visualizd the
28 1 equemene
 evolution of the Path projected on the 2D PES (in green). The actual
29 1 equemene
 IRC is also plotted in red, as with Dep_tscan_p.
30 1 equemene
Here the evolution is shown as a movie, there is no need to press a key to go
31 1 equemene
to the next iteration.
32 1 equemene
33 1 equemene
** Remarks:
34 1 equemene
For all Dep_tscan_XX scripts,
35 1 equemene
The "PES" argument can be omitted if you do not want to plot the 2D PES contours in green.
36 1 equemene
If you have asked (using Path.valid, see manual) to print the tangents, then you can add
37 1 equemene
 "T" as a last  argument, this will print the path and its tangents.
38 1 equemene
All three scripts use the files HCN_irc.dat to plot the actual IRC, as well as HCN_PES_2D.dat
39 1 equemene
 to plot the 2D PES.