Half and full collision matrix methods for scattering: Application to multichannel curve crossing

Yehuda B. Band, Frederick H. Mies

Research output: Contribution to journalArticlepeer-review

14 Scopus citations


We present methods of numerically solving the multichannel Schrödinger equation by propagating exact first-order coupled equations for a specially designed half collision matrix X(r). The method requires choosing a convenient set of reference potentials with which to generate a pair of reference radial functions for each channel. It is easy to tailor the choice of basis to the nature of the exact multichannel interaction matrix in a given region of space, either to enhance the numerical efficiency of propagating the exact X(r), or, to facilitate the introduction of useful approximations. In particular, we define a classical half collision matrix Z(r) by neglecting rapidly oscillating terms in the propagation of X(r). This random phase approximation can only be justified in classically accessible regions of r, and results in a set of first order equations which imposes the unitarity of Z(r)Z†(r) = 1 throughout the propagation. Fortunately, such regions often contain the dominant couplings, particularly for heavy collision partners, and the unitarity constraint generally can be utilized out to r = ∞. We only begin to explore the consequences of this useful result in this paper. Because of the divergence of the reference functions, both the exact X(r) and the approximate Z(r) half collision amplitudes can only be propagated with confidence in classical regions of r. However, using the same basis, we derive a stable algorithm for propagating the full collision matrix Σ(r) = X(X*)-1 which asymptotically defines the scattering matrix. Since Σ(r) is always unitary, the first-order equations for Σ(r) can be propagated in all regions of r. The various methods are demonstrated by applying them to a multichannel curve crossing in atom-diatom vibrational energy transfer collisions. Exact results are obtained from both the complete propagation of Σ(r) and from the partial propagation of X(r) in classical allowed regions. Excellent results are obtained from the unitarized classical half collision matrix Z(r) when reasonable reference potentials are employed.

Original languageEnglish
Pages (from-to)2309-2319
Number of pages11
JournalJournal of Chemical Physics
Issue number4
StatePublished - 1 Jan 1988

ASJC Scopus subject areas

  • General Physics and Astronomy
  • Physical and Theoretical Chemistry


Dive into the research topics of 'Half and full collision matrix methods for scattering: Application to multichannel curve crossing'. Together they form a unique fingerprint.

Cite this