[CP2K:7431] Nonsense VELOCITY DCD files with the default unit bohr/au_time

Krack Matthias (PSI) matthia... at psi.ch
Mon Feb 15 17:44:54 CET 2016

Hi Norman

The DCD format provides the information (coordinates, velocities, and forces) as 4 byte floating point numbers in binary format and thus the accuracy should be fine. You should just dump the 32 bit floating point numbers from the binary DCD in a proper format. Indeed, the standard floating point output using just 6 digits might not be appropriate for some cases. In principle, the DCD format is here even superior, i.e. much more compact/efficient than the usual ASCII XMOL format. I have added a –eformat flag to dumpdcd.f90 code in cp2k/tools which optionally allows the dump in scientific format (ES14.6 instead of F14.6). I hope this option will provide enough significant digits for further processing incl. velocity data.



From: cp... at googlegroups.com [mailto:cp... at googlegroups.com] On Behalf Of NORMAN GEIST
Sent: 15 February 2016 15:07
To: cp2k
Subject: [CP2K:7431] Nonsense VELOCITY DCD files with the default unit bohr/au_time

Just to warn people and to have this searchable I want to inform that using DCD VELOCITY files using the default unit bohr/au_time
does not represent a reasonable floating point range when saved to DCD format. Cause of the fixed point format obviously used. most values are zero, highest values are 0.001


Norman Geist
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+uns... at googlegroups.com<mailto:cp2k+uns... at googlegroups.com>.
To post to this group, send email to cp... at googlegroups.com<mailto:cp... at googlegroups.com>.
Visit this group at https://groups.google.com/group/cp2k.
For more options, visit https://groups.google.com/d/optout.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20160215/9c6b6ae2/attachment.html>

More information about the CP2K-user mailing list