Abstract
An efficient way of evolving a solution to an ordinary differential equation is presented. A finite element method is used where we expand in a convenient local basis set of functions that enforce both function and first derivative continuity across the boundaries of each element. We also implement an adaptive step-size choice for each element that is based on a Taylor series expansion. This algorithm is used to solve for the eigenpairs corresponding to the one-dimensional soft Coulomb potential, 1/sqrt[x(2)+β(2)], which becomes numerically intractable (because of extreme stiffness) as the softening parameter (β) approaches zero. We are able to maintain near machine accuracy for β as low as β = 10(-8) using 16-digit precision calculations. Our numerical results provide insight into the controversial one-dimensional hydrogen atom, which is a limiting case of the soft Coulomb problem as β → 0.
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 callMore From: Physical review. E, Statistical, nonlinear, and soft matter physics
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
