Proposed is an efficient quantum jump method that is suitable for the study of subrecoil cooling. The method is an improvement over using the numerical integration of Schrodinger equation to find delay time in which the excited state decay term is incorporated into the Hamiltonian making it a non-Hermitian matrix. It can also be applied to the calculation of three-dimensional velocity selective coherent population trapping proposed by Mauri and Arimondo in F = 3/2 to F = 3/2 transition.
