Hi,
I'm trying to implement a zero angle potential with minima at 0 or 2*pi respectively.
For this I just added a phase shift into the former angle_harmonic.hpp potential:
cos_phi = scalar(vec21, vec31) / (vec21_magn * vec31_magn) - (3.141592653589793238462 / 2.000000000000000);
However, in my model I would need also flat angle and right angle definitions.
Moreover, Espresso doesn't allow me using two bond angle potentials.
Please, is there anybody who would share an idea how to implement the zero angles?
Thanks in advance
Kind regards
Dusan
--
____________________
Ing. Dusan Racko, PhD