sella icon indicating copy to clipboard operation
sella copied to clipboard

dihedral range bug?

Open wenyan4work opened this issue 1 year ago • 0 comments

ase returns dihedral values in [0, 2pi], but sella internally uses [-pi,pi]: https://github.com/zadorlab/sella/blob/723756ffa263e7dccf154d01e5e3201e5ab9cc8e/sella/internal.py#L718

but when converting from ase constraints, the range of ase dihedrals seem remain [0, 2pi]: https://github.com/zadorlab/sella/blob/723756ffa263e7dccf154d01e5e3201e5ab9cc8e/sella/internal.py#L1095

As a result, if the ase constraint is added to a dihedral with 300 degrees in ase, this constraint may crash sella, because in sella dihedral range does not go beyond 180 and this constraint can never be satisfied.

I found this seems to be true in my code, and can be worked around if I manually added constraints instead of calling the 'merge_ase_constraints' method. Could you confirm?

wenyan4work avatar Sep 05 '23 21:09 wenyan4work