Abstract
We describe a rapidly converging algorithm for solving the Schrödinger equation with local potentials in real space. The algorithm is based on evolving the Schrödinger equation in imaginary time by factorizing the evolution operator e -εH to fourth order with purely positive coefficients. The states |ψj> and the associated energies extracted from the normalisation factor e -εEj converge as [Formula: see text]. Our algorithm is at least a factor of 100 more efficient than existing second order split operator methods. We apply the new scheme to a spherical jellium cluster with 20 electrons. We show that the low-lying eigenstates converge very rapidly and that the algorithm does not lose any of its effectiveness for very steep potentials.
Sign-in/Register to access full text options 
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a call
