Milne's differential equation and numerical solutions of the Schrodinger equation. II. Complex energy resonance states

Guaranteed resonance enclosures and exclosures for atoms and molecules

In this paper, we confirm, with absolute certainty, a conjecture on a certain oscillatory behaviour of higher auto-ionizing resonances of atoms and molecules beyond a threshold. These results not only definitely settle a more than 30 year old controversy in Rittby et al. (1981 Phys. Rev. A24, 1636–1639 (doi:10.1103/PhysRevA.24.1636)) and Korsch et al. (1982 Phys. Rev. A26, 1802–1803 (doi:10.1103/PhysRevA.26.1802)), but also provide new and reliable information on the threshold. Our interval-arithmetic-based method allows one, for the first time, to enclose and to exclude resonances with guaranteed certainty. The efficiency of our approach is demonstrated by the fact that we are able to show that the approximations in Rittby et al. (1981 Phys. Rev. A24, 1636–1639 (doi:10.1103/PhysRevA.24.1636)) do lie near true resonances, whereas the approximations of higher resonances in Korsch et al. (1982 Phys. Rev. A26, 1802–1803 (doi:10.1103/PhysRevA.26.1802)) do not, and further that there exist two new pairs of resonances as suggested in Abramov et al. (2001 J. Phys. A34, 57–72 (doi:10.1088/0305-4470/34/1/304)).

Accurate complex scaling of three dimensional numerical potentials

The complex scaling method, which consists in continuing spatial coordinates into the complex plane, is a well-established method that allows to compute resonant eigenfunctions of the time-independent Schrödinger operator. Whenever it is desirable to apply the complex scaling to investigate resonances in physical systems defined on numerical discrete grids, the most direct approach relies on the application of a similarity transformation to the original, unscaled Hamiltonian. We show that such an approach can be conveniently implemented in the Daubechies wavelet basis set, featuring a very promising level of generality, high accuracy, and no need for artificial convergence parameters. Complex scaling of three dimensional numerical potentials can be efficiently and accurately performed. By carrying out an illustrative resonant state computation in the case of a one-dimensional model potential, we then show that our wavelet-based approach may disclose new exciting opportunities in the field of computational non-Hermitian quantum mechanics.