Return
Kepler Predictor-Corrector Algorithm: Scattering Dynamics with One-Over-R Singular Potentials
DOI:10.1021/ct200452h.png)
Abstract
En 中文
An accurate and efficient algorithm for dynamics simulations of particles with attractive 1/r singular potentials is introduced. The method is applied to semiclassical dynamics simulations of electron proton scattering processes in the Wigner-transform time-dependent picture, showing excellent agreement with full quantum dynamics calculations. Rather than avoiding the singularity problem by using a pseudopotential, the algorithm predicts the outcome of close-encounter two-body collisions for the true 1/r potential by solving the Kepler problem analytically and corrects the trajectory for multiscattering with other particles in the system by using standard numerical techniques (e.g., velocity Verlet, or Gear Predictor corrector algorithms). The resulting integration is time-reversal symmetric and can be applied to the general multibody dynamics problem featuring close encounters as occur in electron ion scattering events, in particle antiparticle dynamics, as well as in classical simulations of charged interstellar gas dynamics and gravitational celestial mechanics.
Keywords:
MOLECULAR-DYNAMICS
SCHRODINGER-EQUATION
COMPUTER-SIMULATION
PERTURBATION THEORY
SPLIT INTEGRATION
SYMPLECTIC METHOD
SPECTRAL METHOD
ELECTRON
STATES
FORM
AI Summary
Key information extracted from the uploaded paper, including a brief overview, abstract, background, key highlights, visual analysis, and future outlook.
Journal
IF:
5.5
Papers:
1.1W
Citations:
5.4W

