Convert atomic velocities from DL_POLY to cp2k format

iskarmou dcya... at hotmail.com
Wed Jun 17 16:49:17 UTC 2015


Dear all,

I am trying to convert the velocities from a pre-equilibrated configuration
obtained with DL_POLY file to the cp2k format.
The velocities in the DLPOLY REVCON file are expressed in A/ps.
According to the cp2k manual the velocity unit is:


*bohr/au_time*




*The au_time is derived from the hbar value (1.054e-34 J*sec) and the value 
of the hartree unit of energy (27.21 eV or 4.359e-18 J) as hbar/Ehartree = 
2.42e-17 sec = 0.02042 fs*With a simple calculation (taking int account 
that 1 Bohr=0.52918 A) the cp2k velocity unit should be
0.52918/0.0242 A/fs = 25.914789 * 10**3 A/ps

However, when I am doing this convertion the initial temperature that I get 
in my cp2k BO-MD run is much lower
than the desired one.

Am I doing a wrong conversion?

Thanks in advance.


Yannis Skarmoutsos
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20150617/eebe923d/attachment.htm>


More information about the CP2K-user mailing list