On the computation of matrix elements between numerical wave functions: The canonical functions method

✍ Scribed by H. Kobeissi; M. Dagher; A. El-Hajj; M. Kobeissi

• Publisher: John Wiley and Sons
• Year: 1989
• Tongue: English
• Weight: 727 KB
• Volume: 10
• Category: Article
• ISSN: 0192-8651
• DOI: 10.1002/jcc.540100309

✦ Synopsis

The problem of the computation of the matrix elements Z ( v , v ' ; k ) = [ q U ( r ) ( r -rJkqUr(r)dr, is considered when q u ( r ) and q J r ) are eigenfunctions related to a diatomic potential of the RKR type (defined by the coordinates of its turning points Pi with polynomial interpolations). The eigenfunction Y(r) is computed by the canonical functions method making use of the abscissas rj of Pi uniquely. This limited number of points allows the storage of Tu(ri) for all the required levels v, and reduces greatly the computational effort when v, v', and k are varying. The present method maintains all the advantages of a highly accurate numerical method (even for levels near the dissociation), and reduces greatly the computing time. Furthermore, it is shown that it may be extended to analytical potentials like Morse and Lennard-Jones functions, to vibrational-rotational eigenfunctions and to matrix elements between eigenfunctions related to two different potentials. Numerical applications are presented and discussed.