follower force

Submitted by zamir on Mon, 01/09/2012 - 21:04

Hello! Does anyone have any good ideas if I'm trying to create a follower force. In other words, I want a force that changes direction as a particle rotates. I cannot find any commands to create mobile coordinate systems on particles, is there another way?

Thanks,
Zamir

ckloss's picture

ckloss | Wed, 01/11/2012 - 21:08

You mean a force depending on particle orientation? Starting from fix_addforce.h/cpp, this is pretty straight-forward to do.
You can use fix nve/sphere with the dipole keyword so the orientation is tracked (orientation is stored in atom->mu). If questions arrise, feel free to use the developers forum.

Cheers, Christoph