I have a question about harmonic functions and hitting probabilities.
Let $d \ge 3$ be an integer.
Let $D=\mathbb{R}^d \times (-1,1)$ and denote points $z \in D$ by $z=(r,\theta,y),$ where $(r,\theta)$ denotes the polar coordinates in $\mathbb{R}^d$ and $y \in (-1,1)$. Let $D_0=\{(r,\theta,y) \in D \mid r<1\}$ denote the unit cylinder in $D$.
Let $X=(X_t,P_z)$ be the (normally) reflecting Brownian motion on $\bar{D}$. We define $\sigma=\inf \{t \ge 0 \mid X_t \in D_0\}$. Then, $P_{z}(\sigma<\infty)$ is a harmonic function on $D \setminus D_0$ and satisfies the Neumann boundary condition on $\partial D \setminus \partial D_0$.
There are some harmonic functions on $D \setminus D_0$ and satisfies the Neumann boundary condition on $\partial D \setminus \partial D_0$. For example, if we define $h : \bar{D} \mapsto \mathbb{R}$ by \begin{align*} h(z)&=r^{-d+2}, \quad (z=(r,\theta,y) \in \bar{D},\ r\ge 1), \\ h(z)&=1,\quad (z=(r,\theta,y) \in \bar{D},\ r<1), \end{align*} $h$ is a desired one.
The following should hold: \begin{equation*} (1)\quad P_z(\sigma<\infty)=h(z), \, z \in \bar{D}. \end{equation*}
My question
Is there a theorem which justifies (1)?
What are sufficient conditions for uniqueness of harmonic functions?