In my last mail, rotmat should be rotaxis, of course. So: rotation_90_around_X=rotaxis(pi/2.0, Vector(1,0,0)) translation=array((1,0,0), 'f') atom.transform(rotation_90_around_X, translation) -Thomas