Edge dislocations are crucial in understanding both mechanical and electrical transport in solid and are modeled as line distributions of dipole moments. The calculation of the electronic spectrum for the two dimensional dipole, represented by the potential energy V(r,θ)=pcosθ/r, has been the topic of several studies that show significant difficulties in obtaining accurate results. In this work, we demonstrate that the source of these difficulties is a logarithmic contribution to the behavior of the wave function at the origin that was neglected by previous authors. By taking into account this non-analytic deviation of the solution of Schrödinger's equation, superior results, with the expected rate of convergence, are obtained. This goal is accomplished by "adapting" general algorithms for solving partial derivative differential equations to include the desired asymptotic behavior. We illustrate this principle for the variational principle and finite difference methods. Accurate energies and wave functions are obtained not only for the ground state but also for the first eleven excited states and are useful for designing nanoelectronic devices. This paper demonstrates that augmentary knowledge about analytic properties of the solutions leads to the improved convergence and stability of numerical methods.
Keywords: edge dislocations; electronic states; numerical methods.