Six-dimensional quantum dynamical calculations are presented for the reaction of (v, j) H 2 on Cu͑100͒, at normal incidence, for vϭ0 -1 and jϭ0 -5. The dynamical calculations employed a potential energy surface computed with density functional theory, using the generalized gradient approximation and a slab representation for the adsorbate-substrate system. The aim of the calculations was to establish signatures from which experiments could determine the dominant reaction site of H 2 on the surface and the dependence of the reaction site on the initial rovibrational state of H 2 . Two types of signatures were found. First, we predict that, at energies near threshold, the reaction of (vϭ1) H 2 is rotationally enhanced, because it takes place at the top site, which has an especially late barrier and a reaction path with a high curvature. On the other hand, we predict the reaction to be almost independent of j for (vϭ0) H 2 , which reacts at the bridge site. Second, we predict that, at collision energies slightly above threshold for which the reaction probabilities of the (vϭ0) and (vϭ1) states are comparable, the rotational quadrupole alignment of (vϭ1) reacting molecules should be larger than that of (vϭ0) reacting molecules, for jϭ1, 4, and 5. For ( j ϭ2) H 2 , the opposite should be true, and for ( jϭ3) H 2 , the rotational quadrupole alignment should be approximately equal for (vϭ1) and (vϭ0) H 2 . These differences can all be explained by the difference in the predicted reaction site for (vϭ1) and (vϭ0) H 2 ͑top and bridge͒ and by the differences in the anisotropy of the potential at the reaction barrier geometries associated with these sites. Our predictions can be tested in associative desorption experiments, using currently available experimental techniques.