VASP-MD velocities and VACF

Vasp transition state theory tools

Moderator: moderators

santana
Posts: 1
Joined: Fri Jun 17, 2011 6:13 pm

VASP-MD velocities and VACF

Post by santana »

Dear all,
I want to evaluate the VACF and frequencies using MD from VASP. I found a post in the VASP forum where the xdat2vdat.pl was recommended to generate the velocities using the XDATCAR.
Is this script not longer available?
graeme
Site Admin
Posts: 2276
Joined: Tue Apr 26, 2005 4:25 am
Contact:

Re: VASP-MD velocities and VACF

Post by graeme »

This script may not have been updated for vasp5.2. Vasp has changed the format of the xdatcar, and I don't think anyone has used this script since then. I'll attach a copy of it here in case it is of use. But do check the output to make sure that it makes sense. It requires a massfile as well, with a simple list of the atom number and the mass, one per line. It may also be limited to orthogonal boxes, which is another reason we may have removed it from the tar file.
Attachments
xdat2vdat.pl.gz
(1.23 KiB) Downloaded 8438 times
stauff
Posts: 4
Joined: Mon Feb 13, 2012 4:59 pm

Re: VASP-MD velocities and VACF

Post by stauff »

The xdat2vdat.pl script has been updated for vasp5.2-format XDATCAR files. It also requires the associated OUTCAR file. It can be found with the other vtstscripts.
vimal
Posts: 14
Joined: Mon Nov 22, 2010 10:19 pm

Re: VASP-MD velocities and VACF

Post by vimal »

I am using this script to calculate the VACF.
I end up with the VDATCAR file which looks similar to
XDATCAR with the following additional lines for each poscar.

KINETIC ENERGY (eV): 9.591403906717447
TEMP EFF (K): 290.989883855749042

my question is how should i proceed from here to get the VACF.
Could you please help.

best regards,
vimal
graeme
Site Admin
Posts: 2276
Joined: Tue Apr 26, 2005 4:25 am
Contact:

Re: VASP-MD velocities and VACF

Post by graeme »

We do not have a script to do this. Now that you have the velocities, taking a bunch of dot products to get the autocorrelation function is not too difficult. If you write such a script, we would be happy to include it in the set.
vimal
Posts: 14
Joined: Mon Nov 22, 2010 10:19 pm

Re: VASP-MD velocities and VACF

Post by vimal »

Dear Graeme,

I find that the KE and Temperature are listed in
the OUTCAR file already.

kinetic Energy EKIN = 10.100042 (temperature 305.22 K)
kinetic Energy EKIN = 10.307751 (temperature 311.50 K)
kinetic Energy EKIN = 10.466675 (temperature 316.30 K)

But it really differs from the one calculated from xdat2vdat script that
is attached on the VDATCAR.

Why this happens?
May be i am wrong. Please correct.

Best regards,
Vimal
graeme
Site Admin
Posts: 2276
Joined: Tue Apr 26, 2005 4:25 am
Contact:

Re: VASP-MD velocities and VACF

Post by graeme »

Can you post your OUTCAR, POSCAR, XDATCAR, and VDATCAR so that we can look into this?
anoop
Posts: 1
Joined: Thu Feb 06, 2014 2:20 pm

Re: VASP-MD velocities and VACF

Post by anoop »

Dear all,

I had faced a similar problem mentioned in this post.
From OUTCAR
kinetic Energy EKIN = 3.786618 (temperature 305.15 K)
kinetic Energy EKIN = 3.913352 (temperature 315.36 K)

Where as after using the xdat2vdat script shows the lower values.
KINETIC ENERGY (eV): 1.310660628736558
TEMP EFF (K): 105.622127787439126

Could you please sugest, where It went wrong. Many thanks in advance.

Thanking you,

Regards,
Anoop
graeme
Site Admin
Posts: 2276
Joined: Tue Apr 26, 2005 4:25 am
Contact:

Re: VASP-MD velocities and VACF

Post by graeme »

If you can post a .tar.gz file of the calculation, I'll take a look. This is not a widely used script and there may be some problem with it.
tyuttian
Posts: 2
Joined: Sat Apr 26, 2014 9:37 am

Re: VASP-MD velocities and VACF

Post by tyuttian »

VASP version for my compute is 4.6.3, i use the script (see attachment) and have a "outcar"," massfile ",''xdatcar'' in the same file.
but the vdatcr can not be produced completely(see attachment).

The formation of 'massfile': (for the example of CH4):

1 1 12.0
2 1 1.0
3 1 1.0
4 1 1.0
5 1 1.0
Attachments
VDATCAR.zip
vdatcar
(3.1 KiB) Downloaded 9556 times
xdat2vdat.pl.gz
scripts file
(1.23 KiB) Downloaded 8934 times
graeme
Site Admin
Posts: 2276
Joined: Tue Apr 26, 2005 4:25 am
Contact:

Re: VASP-MD velocities and VACF

Post by graeme »

I would need all the input files to reproduce the issue and see what is wrong.
It is possible that the script only works for vasp.5.*
graeme
Site Admin
Posts: 2276
Joined: Tue Apr 26, 2005 4:25 am
Contact:

Re: VASP-MD velocities and VACF

Post by graeme »

anoop, there does clearly seem to be a problem here. Two new students will be taking a look at this and hopefully make an update with a corrected script.
tyuttian
Posts: 2
Joined: Sat Apr 26, 2014 9:37 am

Re: VASP-MD velocities and VACF

Post by tyuttian »

graeme, thank you very much for replying and helping. These days are national holiday of my country, so my replay is late.

The script of "xdat2vdat.pl" have two versions: One for vasp.5.*, and another of vasp.4.*.(Right?). I selected the second one.

The input files are too large to upload. I delete some content, only save two steps of 5000.
Attachments
input.zip
inputfiles
(75.43 KiB) Downloaded 9823 times
wc5879
Posts: 21
Joined: Thu May 22, 2014 7:49 pm

Re: VASP-MD velocities and VACF

Post by wc5879 »

Hi Y'all
It turned out that the script
xdat2vdat.pl
cannot read mass properly from the OUTCAR
this has been corrected

It has been verified that the script calculates velocities and kinetic energies and effective temperatures correctly. The calculation is strictly based on XDATCAR therefore its result may deviate from VASP by a little bit. The deviation is most likely due to the induced velocity changes done by VASP, as a result of certain flags such as SMASS

So, the script should do exactly what it is meant for.
If there is any other problems, let me know!

Ray
keduemml
Posts: 1
Joined: Thu Nov 12, 2020 6:25 pm

Re: VASP-MD velocities and VACF

Post by keduemml »

Hello,
Where can I find the most updated version of the script? And does anyone know if it is compatible with VASP 6?
Post Reply