Question about dynamical matrix calculation
Posted: Tue Jan 26, 2010 3:01 pm
Dear guys,
To learn how to get the prefactor of a reaction by a dynamical matrix calculation, I have built a model of H difusion on Pd (111) surface. After performing two separate dynamical calculations for both initial state and transition state, I issued commands dymmatrix.pl DISPLACECAR OUTCAR to get the two freq.dat files for dymprefactor.pl.
The freq.dat file from initial state looks like:
44.264242 cm^{-1} ... 1
799.984823 cm^{-1} ... 0
1012.129608 cm^{-1} ... 0
The freq.dat file from transition state looks like:
288.652000 cm^{-1} ... 1
0.000013 cm^{-1} ... 0
977.836375 cm^{-1} ... 0
Then I issued the command dymprefactor.pl [freq.dat of initial state] [freq.dat of transition state]. And I received a warning:
Initial state freq.dat has imaginary frequencies.
What should I do to solve this problem and get the value of prefactor?
Thanks a lot.
To learn how to get the prefactor of a reaction by a dynamical matrix calculation, I have built a model of H difusion on Pd (111) surface. After performing two separate dynamical calculations for both initial state and transition state, I issued commands dymmatrix.pl DISPLACECAR OUTCAR to get the two freq.dat files for dymprefactor.pl.
The freq.dat file from initial state looks like:
44.264242 cm^{-1} ... 1
799.984823 cm^{-1} ... 0
1012.129608 cm^{-1} ... 0
The freq.dat file from transition state looks like:
288.652000 cm^{-1} ... 1
0.000013 cm^{-1} ... 0
977.836375 cm^{-1} ... 0
Then I issued the command dymprefactor.pl [freq.dat of initial state] [freq.dat of transition state]. And I received a warning:
Initial state freq.dat has imaginary frequencies.
What should I do to solve this problem and get the value of prefactor?
Thanks a lot.