Dynamic matrix calculation gave me all 0 frequencies, can't figure out why

Vasp transition state theory tools

Moderator: moderators

Post Reply
hanc0129
Posts: 4
Joined: Tue Feb 23, 2016 3:23 pm

Dynamic matrix calculation gave me all 0 frequencies, can't figure out why

Post by hanc0129 »

I have tried to carry out a dynamic matrix calculation of a simple carbon monoxide molecule (CO), in an attempt to figure out its zero point energy (ZPE). I have carried out many attempts but the frequencies I got from using dymmatrix.pl were always all 0. It's the first time I'm doing a dynamic matrix calculation so can you tell me where went wrong in my code? I'm attaching the INCAR, POSCAR, DISPLACECAR and the OUTCAR that comes after the calculation. The molecular geometry was already optimized in the POSCAR file.

Thank you!
Attachments
OUTCAR.txt
(185.53 KiB) Downloaded 805 times
POSCAR.txt
(566 Bytes) Downloaded 814 times
INCAR.txt
(3.16 KiB) Downloaded 827 times
DISPLACECAR.txt
(59 Bytes) Downloaded 818 times
graeme
Site Admin
Posts: 2280
Joined: Tue Apr 26, 2005 4:25 am
Contact:

Re: Dynamic matrix calculation gave me all 0 frequencies, can't figure out why

Post by graeme »

You need to set ICHAIN=1 to use our dynamical matrix code.

There are additional instructions here: http://theory.cm.utexas.edu/vtsttools/dymmat.html
hanc0129
Posts: 4
Joined: Tue Feb 23, 2016 3:23 pm

Re: Dynamic matrix calculation gave me all 0 frequencies, can't figure out why

Post by hanc0129 »

Hi. I have changed value of ICHAIN from 0 to 1. But the calculation still gives me all 0 frequencies. Are there additional problems with my code? Should the atoms in POSCAR all set to "T" or "F" for this to work?

Many thanks
graeme
Site Admin
Posts: 2280
Joined: Tue Apr 26, 2005 4:25 am
Contact:

Re: Dynamic matrix calculation gave me all 0 frequencies, can't figure out why

Post by graeme »

Ah, good point, yes, set the atomic degrees of freedom to T and change the sign of your EDIFFG. With everything frozen the force will be precisely zero and the run will stop after one iteration.
Post Reply