Fourth order real space solver for the time-dependent Schrödinger equation with singular Coulomb potential

Abstract

We present a novel numerical method and algorithm for the solution of the 3D axially symmetric time-dependent Schrödinger equation in cylindrical coordinates, involving singular Coulomb potential terms besides a smooth time-dependent potential. We use fourth order finite difference real space discretization, with special formulae for the arising Neumann and Robin boundary conditions along the symmetry axis. Our propagation algorithm is based on merging the method of the split-operator approximation of the exponential operator with the implicit equations of second order cylindrical 2D Crank-Nicolson scheme. We call this method hybrid splitting scheme because it inherits both the speed of the split step finite difference schemes and the robustness of the full Crank-Nicolson scheme. Based on a thorough error analysis, we verified both the fourth order accuracy of the spatial discretization in the optimal spatial step size range, and the fourth order scaling with the time step in the case of proper high order expressions of the split-operator. We demonstrate the performance and high accuracy of our hybrid splitting scheme by simulating optical tunneling from a hydrogen atom due to a few-cycle laser pulse with linear polarization.

References in corpus (9)

• The Magnus expansion and some of its applications
• QPROP: A Schroedinger-solver for intense laser-atom interaction
• A brief compendium of time-dependent density-functional theory
• Accurate numerical solutions of the time-dependent Schrödinger equation
• Photoionization of helium by attosecond pulses: extraction of spectra from correlated wave functions
• Fourth Order Gradient Symplectic Integrator Methods for Solving the Time-Dependent Schrödinger Equation
• Any order imaginary time propagation method for solving the Schrodinger equation
• Quantum and semiclassical phase-space dynamics of a wave packet in strong fields using initial-value representations
• Oscillations in Quantum Entanglement During Rescattering

Cited by in corpus (8)

• Improved one-dimensional model potentials for strong-field simulations
• Quantum entanglement in strong-field ionization
• Quantum interference in strong-field ionization by a linearly polarized laser pulse, and its relevance to tunnel exit time and momentum
• Density-based one-dimensional model potentials for strong-field simulations in $\text{He}$, ${\text{H}}_2^{+}$ and ${\text{H}}_2$
• One-dimensional model potentials optimized for the calculation of the HHG spectrum
• Diatomic molecule in a strong infrared laser field: level-shifts and bond-length change due to laser-dressed Morse potential
• Coulomb-Corrected Strong Field Approximation without Singularities and Branch Points
• Nonclassical trajectories in head-on collisions