  We study equal weight numerical integration, or Quasi Monte Carlo (QMC)
rules, for functions in a Sobolev space $H^s(S^d)$ with smoothness parameter
$s>d/2$ defined over the unit sphere $S^d$ in $R^{d+1}$. Focusing on $N$-point
sets that achieve optimal order QMC error bounds (as is the case for efficient
spherical designs), we are led to introduce the concept of QMC designs: these
are sequences of $N$-point node sets $X_N$ on $S^d$ such that the worst-case
error of the corresponding QMC rules satisfy a bound of order $O(N^{-s/d})$ as
$N\to\infty$ with an implied constant that depends on the $H^s(S^d)$-norm.
  We provide methods for generation and numerical testing of QMC designs. As a
consequence of a recent result of Bondarenko et al. on the existence of
spherical designs with appropriate number of points, we show that minimizers of
the $N$-point energy for the reproducing kernel for $H^s(S^d)$, $s>d/2$, form a
sequence of QMC designs for $H^s(S^d)$. Furthermore, without appealing to the
Bondarenko et al. result, we prove that point sets that maximize the sum of
suitable powers of the Euclidean distance between pairs of points form a
sequence of QMC designs for $H^s(S^d)$ with $s\in(d/2,d/2+1)$.
  Numerical experiments suggest that many familiar sequences of point sets on
the sphere (equal area, spiral, minimal [Coulomb or log.] energy, and Fekete
points) are QMC designs for appropriate values of $s$. For comparison purposes
we show that sets of random points that are independently and uniformly
distributed on the sphere do not constitute QMC designs for any $s>d/2$.
  If $(X_N)$ is a sequence of QMC designs for $H^s(S^d)$, we prove that it is
also a sequence of QMC designs for $\mathbb{H}^{s'}(S^d)$ for all
$s'\in(d/2,s)$. This leads to the question of determining the supremum of such
$s$, for which we provide estimates based on computations for the
aforementioned sequences.
