Harnack inequality for Bessel operators

We prove uniqueness results and Harnack inequality for Bessel operators \begin{align*} %\label{def L transf alpha} D_t-Δ_{x} -2a\cdot\nabla_xD_y- D_{yy}- \frac cy D_y % \nonumber \\[1ex]&=y^α\sum_{i,j=1}^{N+1}a_{ij}D_{ij}+y^{α-1}\left(v,\nabla\right)-by^{α-2}. \end{align*} in the strip $[0,T]\times \mathbb{R}^{N+1}_+=\{0 \leq t \leq T, x \in \mathbb{R}^N, y>0\}$ under Neumann boundary conditions at $y=0$.