Ineffective FIXED_ATOMS constraints in QS MD

Marco Caccin marco.... at gmail.com
Tue Jan 13 16:20:28 UTC 2015


Dear all,

I am trying to run Quickstep MD with some fixed atoms positions (set in 
MOTION / CONSTRAINT / FIXED_ATOMS), but the positions do change over time.
Curiously enough, constrained atoms have zero velocity at frame 0 and zero 
force at every timestep, which would suggest that the constraint is somehow 
applied but ineffective.

This behaviour seems to be machine and version independent:
tested cp2k versions: 2.7 (Development Version) git: 0826abf, 2.6 
(Development Version) svn:14226, 2.5.1 svn:13632 
tested architectures/compilers: gfortran 4.8 on Darwin for cp2k 2.7, xlf 
14.1 on BlueGene/Q for cp2k 2.5 and 2.6

The attachments are a minimal example of the problem, showing 5 timesteps 
of dynamics of 3 atoms of Si of which atom 1 is (or rather, should be) 
fixed.

Is this problem known, and is there any workaround? I know there is for 
GEO_OPT runs as 
in https://groups.google.com/forum/#!topic/cp2k/KAMhME81IbM, but I haven't 
found anything applicable to MD

Kindest regards,

Marco

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20150113/d62ff36c/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test_constrained_MD_minimal.inp
Type: chemical/x-gamess-input
Size: 2809 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20150113/d62ff36c/attachment.inp>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: md1-pos-1.xyz
Type: chemical/x-xyz
Size: 1620 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20150113/d62ff36c/attachment.xyz>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: md1-frc-1.xyz
Type: chemical/x-xyz
Size: 1620 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20150113/d62ff36c/attachment-0001.xyz>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: md1-vel-1.xyz
Type: chemical/x-xyz
Size: 1620 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20150113/d62ff36c/attachment-0002.xyz>


More information about the CP2K-user mailing list