The time‐dependent quantum wave packet and the quasi‐classical trajectory (QCT) calculations for the title reactions are carried out using three recent‐developed accurate potential energy surfaces of the 11A′, 13A′, and 13A″ states. The two commonly used polarization‐dependent differential cross sections, dσ00/dωt, dσ20/dωt, with ωt being the polar coordinates of the product velocity ω′, and the three angular distributions, P(θr), P(Φr), and P(θr,Φr), with θr, Φr being the polar angles of the product angular momentum, are generated in the center‐of‐mass frame using the QCT method to gain insight into the alignment and the orientation of the product molecules. Influences of the potential energy surface, the collision energy, and the isotope mass on the stereodynamics are shown and discussed. Validity of the QCT calculation has been examined and proved in the comparison with the quantum wave packet calculation. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010