I have an unbounded symmetric operator, and I would like to find its self-adjoint extension if possible. First off, what properties does such an operator need in order to have a self-adjoint extension, and what are the steps I need to take to get this extension?
For concreteness, the operator I am looking at is an infinite dimensional matrix $M_{nm}$ acting on $l^2$ which reads
$$M_{nm}=\frac{(2n+1)!!}{(2n)!!}\Big((4n+3)\delta_{nm}+2n\delta_{n-1m}+(2n+3)\delta_{n+1m}\Big)$$
This is symmetric with respect to the inner product $(u,v)=\sum_{n=0}^{\infty}u^*_nv_n$.