Numerical integration of the time-dependent Schrödinger equation for laser-driven helium

The full time-dependent Schrodinger equation for 2-electron atoms in intense laser fields is solved using a mixed finite-difference/basis set approach. We discuss Krylov subspace techniques for the propagation of the equation in time, and numerical methods for optimizing the description of the system on a finite-difference lattice. We describe the implementation of a parallelized code based on these numerical methods, and review performance and scaling results of the code on the Cray T3D/E.

[1]  H. Tal-Ezer,et al.  An accurate and efficient scheme for propagating the time dependent Schrödinger equation , 1984 .

[2]  Henry C. Kapteyn,et al.  GENERATION OF COHERENT SOFT X RAYS AT 2.7 NM USING HIGH HARMONICS , 1997 .

[3]  R. Heather An asymptotic wavefunction splitting procedure for propagating spatially extended wavefunctions: application to intense field photodissociation of H + 2 , 1991 .

[4]  Andr'e Nauts,et al.  New Approach to Many-State Quantum Dynamics: The Recursive-Residue-Generation Method , 1983 .

[5]  W. Arnoldi The principle of minimized iterations in the solution of the matrix eigenvalue problem , 1951 .

[6]  N. H. Burnett,et al.  Coherent XUV Generation from Gases Ionized by Several Cycle Optical Pulses , 1997 .

[7]  H. V. D. Vorst,et al.  An iterative solution method for solving f ( A ) x = b , using Krylov subspace information obtained for the symmetric positive definite matrix A , 1987 .

[8]  C. Clark,et al.  Study of a plane-wave final-state theory of above-threshold ionization and harmonic generation , 1996 .

[9]  T. Park,et al.  Unitary quantum time evolution by iterative Lanczos reduction , 1986 .

[10]  G. Stewart,et al.  Reorthogonalization and stable algorithms for updating the Gram-Schmidt QR factorization , 1976 .

[11]  Yousef Saad,et al.  Efficient Solution of Parabolic Equations by Krylov Approximation Methods , 1992, SIAM J. Sci. Comput..

[12]  Jennifer A. Scott,et al.  An Arnoldi code for computing selected eigenvalues of sparse, real, unsymmetric matrices , 1995, TOMS.