We present an alteative approach for determining the photoionization rate of hydrogen molecules under the interaction of intense light, by calculating the spatial overlap integral between the potential function and the time-dependent wavefunction. The suggested method was applied to varying excitation pulse shapes: square envelope and chirped hyperbolic secant envelope. The computed results confirmed that our method was robust and could be extended to general molecular dynamics calculations.