Let $(X\_{t})\_{t\geq 0}$ be a Bessel Process starting at $x>0$ of dimension $\delta>0$. Namely 
\begin{align\*}
X\_{t}=x+W\_{t}+\frac{\delta-1}{2}\int\_{0}^{t}\frac{1}{X\_{s}}\, ds.
\end{align\*}
where $(W\_{t})\_{t\geq 0}$ is a Brownian Motion. 
I am interested in how to find the distribution of the Hitting Time $\tau:=\inf\{t>0\,|\, X\_{t}=0\}$. 


Given that the Bessel Process can be expressed as a time changed Brownian Motion, up to the first hitting of the boundary, is it possible to obtain the distribution of $\tau$ by utilising the Reflection Principle for Brownian Motion?



Obtaining the distribution of the First Hitting time of the Bessel Process
Since the Bessel process has Brownian scaling, we have 
$$
\mathbb{P}^x(\tau<t)=F(xt^{-\frac12})
$$
for some unknown function $F$, where $\mathbb{P}^x$ denotes the probability for the Bessel process started from $x$. Now, either from Fokker-Plank equation, or by Ito calculus, or conditioning on the exit time and position from an $\epsilon$-neighborhood of the starting point and sending $\epsilon$ to zero, you can see that $F$ satisfies the following differential equation:
$$
F'(y)\left(y+\frac{\delta-1}{y}\right)+F''(y)=0,
$$
which can be easily solved:
$$
F(y)=c\_1\int^ys^{1-\delta}e^{-s^2/2} ds +c\_2. 
$$
It remains to choose the constants so that $F(0)=1$ and $F(+\infty)=0$. This is of course not possible for $\delta\geq 2$, due to the fact that in this range the process never hits zero.


