From: Boyang Wang (pkuwangboyang_at_gmail.com)
Date: Thu Jan 19 2012 - 17:45:41 CST
I have a comment -- I would care more about how to calculate the force at
the next time step. It is like a wavepacket, and the historical information
of forces are stored in the eigenmodes of the system.
On Fri, Jan 20, 2012 at 1:44 AM, Gungor Ozer <gungor_at_bu.edu> wrote:
> Dear all,
> I have been trying to adapt a TCL script that would use TCLforces to
> rescale force on specified atoms. The idea is simple:
> 1) Calculate current force on a particular atom, f ---
> 2) Rescale it using some factor, x.
> 3) Add force, (x-1)*f --- (addforce)
> I was planning to do this using the built-in loadtotalforces and addforce
> commands. However, there is one issue: loadtotalforces reads the total
> forces on each requested atom in the "previous timestep" but addforce adds
> the force at the "current time step".
> Is there anyway to get around this? Any input will be helpful and greatly
> Thanks in advance,
-- Boyang Wang, Ph.D.
This archive was generated by hypermail 2.1.6 : Mon Dec 31 2012 - 23:21:08 CST