A theoretical time-dependent wave-packet dynamics method is applied to calculate the distribution of vibrational states of hydrogen molecular ions produced in high-energy ionization processes of hydrogen molecules. The isotope effect is elucidated in agreement with the available experimental measurements. Our proposed method should be readily applied in other atomic and molecular processes considering great advances in electronic computation science and technology.