[lammps-users] How to do rotational oscillation


I am trying to apply a rotational oscillation to a group of atoms. That is
when a group of atoms rotates by angle \phi as a function of time --around
a predefined rotation axis and origin.

\phi = \phi_0 * Sin (\omega*t)

I am aware of *fix move* but it seems having an oscillation is not
implemented on the *rotate* style.
Is there other any command to perform such a simulation, other than
intervening on *fix_move.cpp* source file?


fix move could be used for this purpose with the variable keyword.
However, constructing atom style variables for the velocity of the affected atoms could be rather complex.
You need to see how complex a set of variables you need to realize that.
Creating a custom fix may be a simpler solution.