Statistiques
| Révision :

root / examples / Test / Zmat / README @ 8

Historique | Voir | Annoter | Télécharger (1,8 ko)

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
**) The files created by Opt'n Path: Path.dat and Path.gplot
8
**)  AnaPath (as usual)
9
../../../utils/AnaPath HCN_zmat.out 
10
**) Dep_tscan_p 
11
./Dep_tscan_p HCN_zmat.out 50 Path PES
12
This will create a Path_p.gplot file that can be ran to visualizd the
13
 evolution of the Path projected on the 2D PES (in green). The actual
14
 IRC is also plotted in red.
15
To go to the next iteration, you have to press a key.
16

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

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

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