Analysis of a modified Schrödinger operator in 2D: Regularity, index, and FEM

Let r=(x"1^2+x"2^2)^1^/^2 be the distance function to the origin [email protected]?R^2, and let us fix @d>0. We consider the ''Schrodinger-type mixed boundary value problem'' [email protected][email protected]^-^[email protected]?H^m^-^1(@W) on a bounded polygonal domain @[email protected]?R^2. The singularity in the potential @dr^-^2 severely limits the regularity of the solution u. This affects the rate of convergence to u of the finite element approximations u"[email protected]?S obtained using a quasi-uniform sequence of meshes. We show that a suitable graded sequence of meshes recovers the quasi-optimal convergence rate @?u-u"[email protected]?"H"^"1"("@W")@?Cdim(S"n)^-^m^/^[email protected][email protected]?"H"^"m"^"-"^"1"("@W"), where S"n are the FE spaces of continuous, piecewise polynomial functions of degree m>=1 associated to our sequence of meshes and u"n=u"S"""[email protected]?S"n are the FE approximate solutions. This is in spite of the fact that [email protected][email protected]?H^m^+^1(@W) in general. One of the main results of our paper is to show that the singularities due to the potential and the singularities due to the singularities of the domain or to the change in boundary conditions can be treated in the same way. Our proof is based on regularity and well-posedness results in weighted Sobolev spaces, with the weight taking into account all singularities (including the ones coming from the potential). Our regularity results apply also to operators with weaker singularities, like the Schrodinger operator [email protected][email protected]^-^1, for which we also obtain Fredholm conditions and a formula for the index. Our a priori estimates also extend to piecewise smooth domains (i.e., curvilinear polygonal domains).