WATERWAVES: wave particles dynamics on a complex triatomic potential

The WATERWAVES program suite performs complex scattering calculations by propagating a wave packet in a complex, full-dimensional potential for non-rotating (J = 0) but vibrating triatomic molecules. Potential energy and decay probability surfaces must be provided. Expectation values of geometric quantities can be calculated, which are useful for following the wave packet motion. The programs use a local complex potential approximation (LCP) for the Hamiltonian and Jacobi coordinates. The bottleneck of the calculation is the application of each term of the Hamiltonian to the wave packet. To solve this problem the programs use a different representation for each term: normalized associated Legendre polynomials PK j (x) as a functional basis for the angular kinetic term and an evenly spaced grid for the radial kinetic term yielding a fully point-wise representation of the wave functions. The potential term is treated using an efficient Discrete Variable Representation (DVR) being diagonal in the coordinate representation. The radial kinetic term uses a fast Fourier transform (FFT) to obtain an operator which is diagonal in the momentum space. To avoid artificial reflection at the boundaries of the grid a complex absorbing potential is included for calculating continuum quantities. Asymptotic analysis is performed to obtain scattering observables such as cross sections and other dynamical properties.