A Split-Step Numerical Method for the Time-Dependent Dirac Equation in 3-D Axisymmetric Geometry

Abstract

A numerical method is developed to solve the time-dependent Dirac equation in cylindrical coordinates for 3-D axisymmetric systems. The time evolution is treated by a splitting scheme in coordinate space using alternate direction iteration, while the wave function is discretized spatially on a uniform grid. The longitudinal coordinate evolution is performed exactly by the method of characteristics while the radial coordinates evolution uses Poisson's integral solution, which allows to implement the radial symmetry of the wave function. The latter is evaluated on a time staggered mesh by using Hermite polynomial interpolation and by performing the integration analytically. The cylindrical coordinate singularity problem at $r=0$ is circumvented by this method as the integral is well-defined at the origin. The resulting scheme is reminiscent of non-standard finite differences. In the last step of the splitting, the remaining equation has a solution in terms of a time-ordered exponential, which is approximated to a higher order than the time evolution scheme. We study the time evolution of Gaussian wave packets, and we evaluate the eigenstates of hydrogen-like systems by using a spectral method. We compare the numerical results to analytical solutions to validate the method. In addition, we present three-dimensional simulations of relativistic laser-matter interactions, using the Dirac equation.

39 pages, 12 figures

References in corpus (8)

• Chiral tunneling and the Klein paradox in graphene
• Pair production in laser fields oscillating in space and time
• Numerical methods for nonlinear Dirac equation
• Quark-antiquark production from classical fields in heavy ion collisions: 1+1 dimensions
• Calculation of electron-positron production in supercritical uranium-uranium collisions near the Coulomb barrier
• Staggered grid leap-frog scheme for the (2+1)D Dirac equation
• Solution of two-center time-dependent Dirac equation in spherical coordinates: Application of the multipole expansion of the electron-nuclei interaction
• Visualizing the kinematics of relativistic wave packets

Cited by in corpus (15)

• Error estimates of numerical methods for the nonlinear Dirac equation in the nonrelativistic limit regime
• Numerical methods and comparison for the Dirac equation in the nonrelativistic limit regime
• The split-operator technique for the study of spinorial wavepacket dynamics
• Identifying the Stern-Gerlach force of classical electron dynamics
• Relativistic ionization dynamics for a hydrogen atom exposed to super-intense XUV laser pulses
• Krylov subspace methods for the Dirac equation
• Pseudospectral computational methods for the time-dependent Dirac equation in static curved spaces
• Galerkin method for unsplit 3-D Dirac equation using atomically/kinetically balanced B-spline basis
• A spectral method for half-integer spin fields based on spin-weighted spherical harmonics
• Plunging in the Dirac sea using graphene quantum dots
• Efficient state initialization by a quantum spectral filtering algorithm
• Gauge invariant canonical symplectic algorithms for real-time lattice strong-field quantum electrodynamics
• Landau levels for discrete-time quantum walks in artificial magnetic fields
• Error estimates of finite difference methods for the Dirac equation in the massless and nonrelativistic regime
• Short-term evolution of electron wave packet in a constant crossed field with radiative corrections