Numerical solution of singular ODE eigenvalue problems in electronic structure computations

doi:10.1016/j.cpc.2010.05.006

2010-09-01

DOI：10.1016/j.cpc.2010.05.006

原文链接

查看原文

摘要

En 中文

We put forward a new method for the solution of eigenvalue problems for (systems of) ordinary differential equations, where our main focus is on eigenvalue problems for singular Schrodinger equations arising for example in electronic structure computations. In most established standard methods, the generation of the starting values for the computation of eigenvalues of higher index is a critical issue. Our approach comprises two stages: First we generate rough approximations by a matrix method, which yields several eigenvalues and associated eigenfunctions simultaneously, albeit with moderate accuracy. In a second stage, these approximations are used as starting values for a collocation method which yields approximations of high accuracy efficiently due to an adaptive mesh selection strategy, and additionally provides reliable error estimates. We successfully apply our method to the solution of the quantum mechanical Kepler, Yukawa and the coupled ODE Stark problems. (C) 2010 Elsevier B.V. All rights reserved.

Keyword:

Electronic structure computation

Polynomial collocation

Fullpotential core solver

Singular eigenvalue problems