Efficient solver for time-dependent Schrödinger equation with interaction between atoms and strong laser field*

Abstract

We present a parallel numerical method of simulating the interaction of atoms with a strong laser field by solving the time-depending Schrödinger equation (TDSE) in spherical coordinates. This method is realized by combining constructing block diagonal matrices through using the real space product formula (RSPF) with splitting out diagonal sub-matrices for short iterative Lanczos (SIL) propagator. The numerical implementation of the solver guarantees efficient parallel computing for the simulation of real physical problems such as high harmonic generation (HHG) in these interaction systems.

Intensity Effect of Laser-assisted Double Ionization of Helium Atoms by Ultrashort XUV Pulses

“…The initial singlet-spin state is obtained by replacing the real time t in TDSE with the imaginary time τ = it [33]. The numerical solution of the final wave function can be accurately obtained by application of FEDVR [28,32,33,[36][37][38].…”

Laser-assisted XUV double ionization of helium atoms: Intensity dependence of joint angular distributions*

“…[36] The numerical solutions of the final wave functions can be accurately obtained by application of FE-DVR. [31,35,36,[39][40][41] Then the DI probability for double ionization corresponding to the final state with momenta 𝑘 1 and 𝑘 2 can be given by projecting the final wave function to the asymptotic twoelectron wave functions for a long time after the termination of the pulses. For the purpose of removing spurious contribution caused by nonorthogonality of the approximate asymptotic wave function and the initial state, we rewrite the final wave function with exclusion of the overlap between the initial state and the final state.…”