I’m not entirely sure if this has been reported before, but using the set_angle() method on linear triatomic molecules seems return a norm(v)==0 division by zero error. Currently what I do is move the position of one of the molecules by a certain amount along a coordinate and then set the angle as now it is no longer linear. I just wish there was a way this method could check if …