Precise Numerical Solutions of Potential Problems Using Crank-Nicholson Method

Abstract

A new numerical treatment in the Crank-Nicholson method with the imaginary time evolution operator is presented in order to solve the Schrödinger equation. The original time evolution technique is extended to a new operator that provides a systematic way to calculate not only eigenvalues of ground state but also of excited states. This new method systematically produces eigenvalues with accuracies of eleven digits with the Cornell potential that covers non-perturbative regime. An absolute error estimation technique based on a power counting rule is implemented. This method is examined with exactly solvable problems and produces the numerical accuracy down to 10${}^{-11}$.

10 pages, 1 figure. Submitted to Phys. Rev. Lett

References in corpus (3)

• Potential-model calculation of an order-v^2 NRQCD matrix element
• Fast and stable method for simulating quantum electron dynamics
• The Nystrom plus Correction Method for Solving Bound State Equations in Momentum Space

Cited by in corpus (7)

• Schrödinger spectrum generated by the Cornell potential
• $P$-wave quarkonium wavefunctions at the origin in the $\overline{\mathrm{MS}}$ scheme
• Probing QGP droplets with charmonium in high-multiplicity proton-proton collisions
• Simple high-accuracy method for solving bound-state equations with the Cornell potential in momentum space
• Interacting particles in two dimensions: numerical solution of the four-dimensional Schrödinger equation in a hypercube
• Mass spectrum of Mesons via the WKB approximation method
• A method for computing lowest eigenvalues of symmetric polynomial differential operators by semidefinite programming