The N( 4 S)1H 2 reaction and its isotopic variants have been investigated by means of time-dependent quantum wave packet with split operator method on the ground state 1 4 A 00 potential energy surface (Zhai and Han, J. Chem. Phys. 2011, 135, 104314). The reaction probabilities, integral cross sections, branching ratio of the integral cross sections, and effect of vibrational excitation of H 2 , HD, and D 2 diatomic molecules are presented and discussed. The results reveal that the intramolecular isotopic effect is greater than the intermolecular one, and that the vibrational excitation of the diatomic molecules can promote the progress of this reaction. In addition, a limited number of rigorous Coriolis coupling calculations of the integral cross sections of the N( 4 S)1H 2 reaction have been carried out. Also shown is that since the Coriolis coupling plays a small part in this accurate quantum calculation, the cheaper centrifugal sudden calculations here reported are effective for this reactive system.