Solution (source code)

= Solution

<Uniqueness in law> means that, for a fixed initial value, every two <weak stochastic solutions> of the <stochastic differential equation> have the same distribution as entire solution paths, even if their filtered probability spaces and driving <Brownian motions> differ. For an equation on $(0,\infty)$, initially consider the stopped path and lifetime at the boundary.

We first justify that the boundary cannot intervene when $\gamma\ge2$. For any solution of the displayed Bessel equation, let $\tau_{\varepsilon,L}$ be its first exit from $(\varepsilon,L)$, where $0<\varepsilon<r<L$. Applying Itô to $R^2$ gives
$$
\mathbb E R_{t\wedge\tau_{\varepsilon,L}}^2
=r^2+\gamma\mathbb E(t\wedge\tau_{\varepsilon,L})\le L^2.
$$
Thus the annulus exit has finite expectation and is almost surely finite. The one-dimensional generator is $\frac12\partial_{rr}+\frac{\gamma-1}{2r}\partial_r$. Its harmonic scale functions are $\log r$ for $\gamma=2$ and $r^{2-\gamma}$ for $\gamma>2$, as direct differentiation verifies. They are bounded on the stopped annulus, so optional sampling gives
$$
\mathbb P_r(\tau_\varepsilon<\tau_L)=
\begin{cases}
\dfrac{\log(L/r)}{\log(L/\varepsilon)},&\gamma=2,\\
\dfrac{r^{2-\gamma}-L^{2-\gamma}}{\varepsilon^{2-\gamma}-L^{2-\gamma}},&\gamma>2.
\end{cases}
$$
Both tend to zero as $\varepsilon\downarrow0$. A continuous path hitting zero before $\tau_L$ would hit every inner level first, so that event has probability zero. Exhausting the upper levels proves that zero is not hit at any finite time. Nor can an upper explosion occur: localization of $R^2$ gives $\mathbb P(\tau_L\le t)\le(r^2+\gamma t)/L^2\to0$. These are the <Hitting-zero classification for a Bessel process> arguments needed for the dimension-at-least-two case.

Now suppose $\gamma$ is an integer and take a $\gamma$-dimensional <Brownian motion> $U$ from any deterministic vector $u$ with $|u|=r$. For $\rho=|U|>0$, differentiation of the <Euclidean norm> gives
$$
\partial_i|u|=u_i/|u|,\qquad \Delta|u|=(\gamma-1)/|u|.
$$
Define the radial noise $\beta=\sum_i\int U_i/\rho\,dW_i$, filling the integrand by a fixed unit vector on the zero set, as in part (b). It has bracket $t$ and is a <Brownian motion>. the <Itô formula>, initially before a possible zero, gives
$$
d\rho_t=d\beta_t+\frac{\gamma-1}{2\rho_t}dt,\qquad \rho_0=r.
$$
The boundary argument just given makes this a global positive <weak stochastic solution>. The assumed <uniqueness in law> therefore identifies every solution's path law with that of $|U|$. Orthogonal invariance of <Brownian motion> makes this radial law independent of the selected point $u$ on the sphere. Consequently
$$
\boxed{(R_t)_{t\ge0}\ \stackrel{\mathrm{law}}=\ (|u+W_t|)_{t\ge0},\quad |u|=r.}
$$
This proves equality of process laws, not merely equality of one-time marginal distributions.