more work investigating issue with rotating, appears to be only at higher simulation rates, tied to Quaternion __pow__ function with small delta times perhaps