An Efficient and Accurate Grid Method for Solving the Time-Dependent Schroedinger Equation: Application of Coulomb Wave Function DVR to Atomic Systems in Strong Laser Fields

Abstract

We present an efficient and accurate grid method for solving the time-dependent Schrödinger equation of atomic systems interacting with intense laser pulses. As usual, the angular part of the wave function is expanded in terms of spherical harmonics. Instead of the usual finite difference (FD) scheme, the radial coordinate is discretized using the discrete variable representation which is constructed from the Coulomb wave function. For an accurate description of the ionization dynamics of atomic systems, the Coulomb wave function discrete variable representation (CWDVR) method needs 3-10 times less grid points than the FD method. The resultant grid points of CWDVR distribute unevenly so that one has finer grid near the origin and coarser one at larger distances. The other important advantage of the CWDVR method is that it treats the Coulomb singularity accurately and gives a good representation of continuum wave functions. The time propagation of the wave function is implemented using the well-known Arnoldi method. As examples, the present method is applied to the multiphoton ionization of both H and H${}^{-}$ in intense laser fields. Short-time excitation and ionization dynamics of H by static electric fields is also investigated. For a wide range of photon energies and laser intensities, ionization rates calculated using this method are in excellent agreement with those from other theoretical calculations.

Submitted to Journal of Chemical Physics

References in corpus (10)

• QPROP: A Schroedinger-solver for intense laser-atom interaction
• Exact field ionization rates in the barrier suppression-regime from numerical TDSE calculations
• Strong-field ionization of diatomic molecules and companion atoms: strong-field approximation and tunneling theory including nuclear motion
• Lithium Ionization by a Strong Laser Field
• Dynamic tunnelling ionization of $H_2^{+}$ in intense fields
• Optimal reflection-free complex absorbing potentials for quantum propagation of wave packets
• Decay in a uniform field: an exactly solvable model
• A discrete time-dependent method for metastable atoms in intense fields
• Simple Soluble Molecular Ionization Model
• Comments on the Discrete Variable Representation

Cited by in corpus (9)

• Gaussian continuum basis functions for calculating high-harmonic generation spectra
• Microscopic model for the higher-order nonlinearity in optical filaments
• Extracting photoelectron spectra from the time-dependent wave function: Comparison of the projection onto continuum states and window-operator methods
• Photoionization of atoms and molecules studied by Crank-Nicolson method
• Nondipole effects in atomic dynamic interference
• A systematic construction of Gaussian basis sets for the description of laser field ionization and high-harmonic generation
• A GPGPU based program to solve the TDSE in intense laser fields through the finite difference approach
• Simulation of a Hydrogen Atom in Laser Field Using the Time-Dependent Variational Principle
• Coulomb wave discrete variable method for calculation of the fully differential cross sections of antiproton impact ionization of hydrogen atom