We derive an algorithm to calculate quantum transmission and reflection amplitudes (S‐matrix elements) for rearrangement reactions using the invariant imbedding method. Our approach can incorporate closed channels, treat complex potentials (to simulate decay), and can be used within an adiabatic representation (in which derivative coupling terms are present) or within a diabatic representation. The method can also be used to determine bound‐state energies when only closed channels are present. The method is demonstrated on a rearrangement problem with two coupled potential surfaces (four channels) as a function of energy in the energy region from the lowest channel threshold through the highest channel threshold.  

1.
R. Kalaba and D. H. Kagiwada, Integral Equations via Imbedding Methods (Addison-Wesely, Reading, 1974).
2.
R. Scott, Invariant Imbedding and Its Application to Ordinary Differential Equations (Addison-Wesely, Reading, 1973).
3.
D.
Secrest
,
Methods Comput. Phys.
10
,
243
(
1971
).
4.
S. J.
Singer
,
K. F.
Freed
, and
Y. B.
Band
,
J. Chem. Phys.
77
,
1942
(
1982
).
5.
S. J.
Singer
,
S.
Lee
,
K. F.
Freed
, and
Y. B.
Band
,
J. Chem. Phys.
87
,
4762
(
1987
).
6.
I.
Tuvi
and
Y. B.
Band
,
J. Chem. Phys.
99
,
9697
(
1993
).
7.
Y. B.
Band
and
F. H.
Mies
,
J. Chem. Phys.
88
,
2309
(
1988
).
8.
P. S. Julienne and F. H. Mies, in Electron and Atomic Collisions, edited by F. Rostas (de Gruyter, Berlin, 1985), Vol. 3, p. 393;
P. S. Julienne and F. H. Mies, in Electron and Atomic Collisions, edited by D. C. Lorents, W. E. Meyerhof, and J. R. Peterson (Elsevier, Amsterdam, 1986).
9.
R. L.
Dubs
,
P. S.
Julienne
, and
F. H.
Mies
,
J. Chem. Phys.
93
,
8784
(
1990
).
10.
R. L.
Dubs
and
P. S.
Julienne
,
J. Chem. Phys.
95
,
4177
(
1991
).
11.
Y. B.
Band
and
I.
Tuvi
,
J. Chem. Phys.
99
,
9704
(
1993
).
12.
B. L.
Christensen-Dalsgaard
,
Phys. Rev. A
29
,
470
(
1984
).
13.
F. S. Action, Numerical Methods that Usually Work (Harper & Row, New York, 1970), pp. 148–151;
W. H. Press, B. P. Flannery, S. A. Teukolsky, and W. T. Vetterling, Numerical Recipes (Cambridge University, Cam bridge, 1986).
14.
For example, in one of the problems treated below, the integral over the barrier, ∫ k(x)dx, is very large, and therefore the linearly independent solutions of the Schrödinger equation, which can be written as linear combinations of exp[±∫ k(x)dx], become numerically linearly dependent in the numerical propagation of the wave function.
15.
J. C.
Light
and
R. B.
Walker
,
J. Chem. Phys.
65
,
4272
(
1976
);
B. R.
Johnson
,
J. Chem. Phys.
67
,
4086
(
1977
).,
J. Chem. Phys.
16.
F.
Mrugala
and
D.
Secrest
,
J. Chem. Phys.
78
,
5954
(
1983
).
17.
R. Bellman and G. Wing, An Introduction to Invariant Imbedding (Wiley, New York, 1975).
18.
Note the different convention used here for the first derivative matrix Aγ,γ′(r) as compared to Ref. 11, in which Aγ,γ′(r) included a factor of−μ−1.
19.
P. S. Julienne, K.-A. Suominen, and Y. B. Band, Phys. Rev. A (in press).
20.
E. F.
van Dishoeck
,
M. C.
van Hemert
,
A. C.
Allison
, and
A.
Dalgarno
,
J. Chem. Phys.
81
,
5709
(
1984
).
This content is only available via PDF.
You do not currently have access to this content.