Solution of a Three-Body Problem in Quantum Mechanics Using Sparse Linear Algebra on Parallel Computers

A complete description of two outgoing electrons following an ionizing collision between a single electron and an atom or molecule has long stood as one of the unsolved fundamental problems in quantum collision theory. In this paper we describe our use of distributed memory parallel computers to calculate a fully converged wave function describing the electron-impact ionization of hydrogen. Our approach hinges on a transformation of the Schrödinger equation that simplifies the boundary conditions but requires solving very ill-conditioned systems of a few million complex, sparse linear equations. We developed a two-level iterative algorithm that requires repeated solution of sets of a few hundred thousand linear equations. These are solved directly by LU-factorization using a specially tuned, distributed memory parallel version of the sparse LU-factorization library Super-LU. In smaller cases, where direct solution is technically possible, our iterative algorithm still gives significant savings in time and memory despite lower megaflop rates.

[1]  M. Seaton,et al.  Ionization of atomic hydrogen by electron impact , 1965, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[2]  Barry Simon,et al.  The definition of molecular resonance curves by the method of exterior complex scaling , 1979 .

[3]  Joseph W. H. Liu,et al.  Modification of the minimum-degree algorithm by multiple elimination , 1985, TOMS.

[4]  R. Dreizler,et al.  High energy electron impact ionisation of H91s in coplanar asymmetric geometry , 1991 .

[5]  Pan,et al.  Angular distributions for near-threshold (e,2e) processes for H, He, and other rare-gas targets. , 1992, Physical review. A, Atomic, molecular, and optical physics.

[6]  Richard Barrett,et al.  Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods , 1994, Other Titles in Applied Mathematics.

[7]  Hudson,et al.  Convergent R matrix with pseudostates calculation for e--He collisions. , 1996, Physical review. A, Atomic, molecular, and optical physics.

[8]  Jung,et al.  Coulomb three-body effects in low-energy impact ionization of H(1s). , 1996, Physical review. A, Atomic, molecular, and optical physics.

[9]  I. Bray.  Close-Coupling Theory of Ionization: Successes and Failures , 1997 .

[10]  C. W. McCurdy,et al.  Making complex scaling work for long-range potentials , 1997 .

[11]  Xiaoye S. Li,et al.  SuperLU Users'' Guide , 1997 .

[12]  T. N. Rescigno,et al.  Approach to electron-impact ionization that avoids the three-body Coulomb asymptotic form , 1997 .

[13]  James Demmel,et al.  Making Sparse Gaussian Elimination Scalable by Static Pivoting , 1998, Proceedings of the IEEE/ACM SC98 Conference.

[14]  Cleve Ashcraft,et al.  SPOOLES: An Object-Oriented Sparse Matrix Library , 1999, PPSC.

[15]  Electron-impact ionization of atomic hydrogen from near threshold to high energies , 1999, physics/9906008.

[16]  Isaacs,et al.  Collisional breakup in a quantum system of three charged particles , 1999, Science.

[17]  Abhishek Kumar Gupta,et al.  Wsmp: watson sparse matrix package , 2000 .

[18]  F. Robicheaux,et al.  Time-dependent close-coupling calculations for the electron-impact ionization of helium , 2000 .

[19]  Patrick Amestoy,et al.  A Fully Asynchronous Multifrontal Solver Using Distributed Dynamic Scheduling , 2001, SIAM J. Matrix Anal. Appl..

[20]  Patrick R. Amestoy,et al.  Analysis and comparison of two general sparse solvers for distributed memory computers , 2001, TOMS.

[21]  J. Colgan,et al.  Ejected-energy differential cross sections for the near-threshold electron-impact ionization of hydrogen , 2001 .

[22]  W. Isaacs,et al.  Doubly differential cross sections for the electron impact ionization of hydrogen , 2001 .

[23]  C. W. McCurdy,et al.  Accurate amplitudes for electron-impact ionization , 2001 .

[24]  Mark Baertschy,et al.  Electron-impact ionization of atomic hydrogen , 2001 .