A simple recursion polynomial expansion of the Green’s function with absorbing boundary conditions. Application to the reactive scattering

Abstract

The new recently introduced [J. Chem. Phys 102, 7390 (1995)] empirical recursion formula for the scattering solution is here proved to yield an exact polynomial expansion of the operator [E−(Ĥ+Γ̂)]−1, Γ̂ being a simple complex optical potential. The expansion is energy separable and converges uniformly in the real energy domain. The scaling of the Hamiltonian is trivial and does not involve complex analysis. Formal use of the energy-to-time Fourier transform of the ABC (absorbing boundary conditions) Green’s function leads to a recursion polynomial expansion of the ABC time evolution operator that is global in time. Results at any energy and any time can be accumulated simultaneously from a single iterative procedure; no actual Fourier transform is needed since the expansion coefficients are known analytically. The approach can be also used to obtain a perturbation series for the Green’s function. The new iterative methods should be of a great use in the area of the reactive scattering calculations and other related fields.

Vibrational de-excitation of v=1 H2 during collisions with a Cu(100) surface

“…16 The time propagation of the wave packet is evaluated using the absorbing boundary condition ͑ABC͒ evolution operator. 25 The ABC operator contains an optical potential that absorbs the wave packet at the edges of the grid corresponding to large values of Z and r. The initial wave packet is written as…”

Low-energy fusion dynamics of weakly bound nuclei: A time dependent perspective

“…The time evolution operator is represented as a convergent series of modified Chebyshev polynomials [12,13]. In the absence of the imaginary potentials, the norm of the wave function is preserved with an accuracy of ∼ 10 −14 .…”

Unimolecular rovibrational bound and resonance states for large angular momentum: J=20 calculations for HO2

“…In their modified version of Chebyshev propagation, Mandelshtam and Taylor 15 proposed a real damped Chebyshev polynomial recursion to impose the outgoing boundary conditions:…”

“…Kouri and co-workers 14 derived a new time-independent ͑TI͒ wavepacket Lippmann-Schwinger equation and presented Chebyshev expansion expressions for both Green operator and Dirac delta function. Mandelshtam and Taylor 15 introduced a real damping scheme into the Chebyshev recursion which made the real wave-packet method possible for dissipative systems. The real Chebychev propagation method can be viewed in an alternative way as a modification of the timedependent Schrödinger equation.…”