El teorema general es XIII.96 en Reed/Simon "Methods of Modern Mathematical Physics, Vol. IV: Analysis of Operators" (Academic Press, 1978). Un operador de multiplicación $V$ es $-\Delta$ -con límite relativo cero si es uniformemente local $L^p$ donde $p = 2$ si dimensión $n \leq 3$ y $p>n/2$ si $n \geq 4$ . Entonces se puede aplicar Kato-Rellich (Th. X.12 en Reed/Simon Vol. II) para obtener un operador autoadjunto $H = -\Delta + V$ con dominio $H^2(\mathbb{R}^n)$ como bien dices.
Para su potencial de Coulomb singular el $L^p$ -significa la siguiente integral (sobre una pequeña región acotada alrededor del origen, por lo que elegimos una bola $B$ con radio 1) tiene que converger. $$ \int_B |V(x)|^p d^n x $$
Para $V(x)=-|x|^{-1}$ equivale a $$ \int_0^1 r^{-p+n-1} dr < \infty. $$
Obtenemos la condición $-p+n-1 > -1$ así $n > p$ . Así que su afirmación es cierta en dimensiones $n \geq 3$ . Para los potenciales limitados desde abajo, tiene aún más opciones (es decir, localmente $L^2$ para todas las dimensiones), véase http://www.math.caltech.edu/papers/bsimon/xliv.pdf p.3 para una introducción y referencias.
La idea de utilizar formas cuadráticas podría utilizarse para incluir potenciales aún más singulares. Pero yo tendría cuidado con el dominio (se necesita algo así como $\sqrt{|V|} \psi \in L^2$ para $\psi \in H^1$ ) de la forma cuadrática, tal vez puedas encontrar algunas pistas en otro artículo de B. Simon: http://www.math.caltech.edu/SimonPapers/13.pdf
¡Que aproveche! ;-)