tcf forces and ccordinate wrapping

From: Arturas (
Date: Tue Apr 09 2013 - 13:23:21 CDT


I need help on tclFoce or tclBC usage.

I want apply force on atoms (best, on COM of their group) with force normal { 0 0 1}, only within {Z_min < z < Z_max} - kind of a pump and non-eq. MD. Using tclForce, I have increasing coordinates - not wrapped, which breaks down my Z-based atom tagging. If I use tclBS implementation, it is much slower than tclForces (large system) and have 'weird' wrapping, where z coordinates of tagged atoms are different from the input coordinates (e.g. -30 instead of 30; serial IDs are correct).

Could anyone hint or suggest how to effectively use tclForce with atom wrapping, or tclBC with 'right' wrapping?


This archive was generated by hypermail 2.1.6 : Tue Dec 31 2013 - 23:23:07 CST