Orientational order in dense molecular hydrogen: A first principles path integral Monte Carlo calculation