This paper describes a parallel algorithm for quantum dynamics simulat
ions of three body Coulomb problems. The algorithm has been implemente
d on a SIMD (single instruction multiple data) computer (MasPar MPI).
The algorithm employs the finite-difference and split operator (FD-SPO
) method and variable time-step technique. A transformation of wavefun
ction to polar coordinates and a modified finite-difference method are
introduced to reduce the number of the dimension and to handle singul
arities. The variable-time-step scheme improves the efficiency four ti
mes.