vef.py and vef.pl

Vasp transition state theory tools

Moderator: moderators

Post Reply
ibrahim
Posts: 26
Joined: Fri Jun 01, 2018 1:46 am

vef.py and vef.pl

Post by ibrahim »

Dear All,

Recently I found something strange related to vef.py and .pl, the two scripts provide different results for forces, attached is the OUTCAR(1). What happened in that case and which one should I use? With the same vasp version but different installation everything is OK (OUTCAR(2)).

Moreover, when I use newer vasp version (5.4.4), I didn't find "FORCES: max atom, RMS" written in OUTCAR, why?

Many thanks in advance,
Ibrahim
Attachments
files.rar
(1.64 MiB) Downloaded 664 times
graeme
Site Admin
Posts: 2256
Joined: Tue Apr 26, 2005 4:25 am
Contact:

Re: vef.py and vef.pl

Post by graeme »

The newer version of vasp was not built with our vtstcode, which is why that force information is not printed. You can grep the OUTCAR for VTST to check if that code was compiled with your version of vasp.

The difference between vef.pl and vef.py is subtle. The vef.pl script is fine. The vef.py script uses methods in ASE. The issue is that ASE reads any constraints (e.g. frozen atoms) from the POSCAR or CONTCAR file and then applies those constraints to the forces (zeros those components) in the OUTCAR file. If you include a POSCAR or CONTCAR file in the directory with your OUTCAR file, the forces from the two scripts should be the same.
ibrahim
Posts: 26
Joined: Fri Jun 01, 2018 1:46 am

Re: vef.py and vef.pl

Post by ibrahim »

Thanks Prof. Graeme.
Post Reply