A split-step numerical method for the time-dependent Dirac equation in 3-D axisymmetric geometry
From MaRDI portal
Publication:349374
DOI10.1016/J.JCP.2014.03.068zbMATH Open1349.65292arXiv1303.3781OpenAlexW2110498046MaRDI QIDQ349374FDOQ349374
Authors: François Fillion-Gourdeau, Emmanuel Lorin, André D. Bandrauk 
Publication date: 5 December 2016
Published in: Journal of Computational Physics (Search for Journal in Brave)
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 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.
Full work available at URL: https://arxiv.org/abs/1303.3781
Recommendations
- Avoiding the coordinate singularity problem in the numerical solution of the Dirac equation in cylindrical coordinates
- Numerical solution of the time-dependent Dirac equation in coordinate space without fermion-doubling
- Galerkin method for unsplit 3-D Dirac equation using atomically/kinetically balanced B-spline basis
- High-fidelity numerical solution of the time-dependent Dirac equation
- A simple pseudospectral method for the computation of the time-dependent Dirac equation with perfectly matched layers
Cites Work
- FFT-split-operator code for solving the Dirac equation in \(2+1\) dimensions
- Quantum field theory techniques in graphical enumeration
- Title not available (Why is that?)
- Title not available (Why is that?)
- Handbook of linear partial differential equations for engineers and scientists
- Group Velocity in Finite Difference Schemes
- A time-splitting spectral scheme for the Maxwell-Dirac system
- General theory of fractal path integrals with applications to many-body theories and statistical physics
- Pole condition for singular problems: The pseudospectral approximation
- Lattice Boltzmann equation for quantum mechanics
- Quantum dynamics of relativistic electrons
- An efficient and stable numerical method for the Maxwell--Dirac system
- Solution of the Schrödinger equation by a spectral method
- Accelerating the Fourier split operator method via graphics processing units
- Staggered grid leap-frog scheme for the \((2+1)D\) Dirac equation
- Numerical solution of the time-dependent Dirac equation in coordinate space without fermion-doubling
- Title not available (Why is that?)
- A B-spline Galerkin method for the Dirac equation
- Numerical methods for nonlinear Dirac equation
- Title not available (Why is that?)
- A highly accurate technique for the treatment of flow equations at the polar axis in cylindrical coordinates using series expansions.
- Numerical treatment of polar coordinate singularities
- Stabilized finite element method for the radial Dirac equation
- Applications of nonstandard finite difference schemes. Papers from the minisymposium on nonstandard finite difference schemes: theory and applications, SIAM annual meeting, Atlanta, GA, USA, 1999
- The Dirac equation in orthogonal coordinate systems. I. The local representation
- Atoms in strong laser fields: challenges in relativistic quantum mechanics
- Title not available (Why is that?)
- Physical constraints on the coefficients of Fourier expansions in cylindrical coordinates
- High-accuracy Yee algorithm based on nonstandard finite differences: new developments and verifications
- Title not available (Why is that?)
- A simple and accurate mixed \(P^{0}-Q^{1}\) solver for the Maxwell-Dirac equations
Cited In (19)
- A simple pseudospectral method for the computation of the time-dependent Dirac equation with perfectly matched layers
- Numerical methods and comparison for the Dirac equation in the nonrelativistic limit regime
- Stationary state computation for nonlinear Dirac operators
- Error estimates of finite difference methods for the Dirac equation in the massless and nonrelativistic regime
- Explicit and structure-preserving exponential wave integrator Fourier pseudo-spectral methods for the Dirac equation in the simultaneously massless and nonrelativistic regime
- Numerical solution of the time-dependent Dirac equation in coordinate space without fermion-doubling
- Galerkin method for unsplit 3-D Dirac equation using atomically/kinetically balanced B-spline basis
- Numerical solution of the Dirac equation by a mapped Fourier grid method
- Frozen Gaussian Approximation for the Dirac Equation in Semiclassical Regime
- Krylov subspace methods for the Dirac equation
- A uniformly accurate exponential wave integrator Fourier pseudo-spectral method with structure-preservation for long-time dynamics of the Dirac equation with small potentials
- Spatial resolution of different discretizations over long-time for the Dirac equation with small potentials
- Pseudospectral computational methods for the time-dependent Dirac equation in static curved spaces
- High-fidelity numerical solution of the time-dependent Dirac equation
- Absorbing boundary conditions for relativistic quantum mechanics equations
- FFT-split-operator code for solving the Dirac equation in \(2+1\) dimensions
- Computational performance of simple and efficient sequential and parallel Dirac equation solvers
- Landau levels for discrete-time quantum walks in artificial magnetic fields
- Avoiding the coordinate singularity problem in the numerical solution of the Dirac equation in cylindrical coordinates
Uses Software
This page was built for publication: A split-step numerical method for the time-dependent Dirac equation in 3-D axisymmetric geometry
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q349374)
