From: Boyang Wang (pkuwangboyang_at_gmail.com)
Date: Thu Jan 19 2012 - 17:45:41 CST
Hi Gungor,
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.
Best,
Boyang
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 ---
> (loadtotalforces)
> 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
> appreciated.
>
> Thanks in advance,
> Gungor
>
-- Boyang Wang, Ph.D.
This archive was generated by hypermail 2.1.6 : Mon Dec 31 2012 - 23:21:08 CST