# Density of hitting time for a two-sided barrior for Brownian motion with drift

For the following Brownian motion with drift $$X_t = X_0 + \mu t + \sigma B_t$$ where $$\mu \in \mathbb{R}$$, $$\sigma > 0$$ and $$X_0 \in \mathbb{R}$$ which is a solution to the stochastic differential equation: \begin{align*} dX_t = \mu dt + \sigma dB_t \ . \end{align*} I am interested in finding the probability density funtion of the following stopping time: \begin{align*} \tau = \inf[t: X_t \notin (a,b) ] \ . \end{align*} Where we for the initial position have $$X_0 \in (a,b)$$. The corresponding forward Kolmologov equation to the above SDE is: \begin{align*} \frac{\partial u}{\partial t} = \frac{\sigma^2}{2} \frac{\partial^2u}{\partial x^2} - \mu \frac{\partial u}{\partial x} \ . \end{align*} From the answers in the link Density of first hitting time of Brownian motion with drift it should be possible to find the density of $$\tau$$ as a function of $$t$$ by solving the forward Kolmologov equation with the initial condition $$u(x,t=0) = \delta(x-x_0)$$ and with boundary condition $$u(a,t)=u(b,t) = 0$$ for all $$t>0$$. I have already tried the same strategy as mentioned in the link, but it doesn't seem to work. Do anyone have a suggesten to another strategy to solve this problem?

• For two-sided boundary condition the solution is an infinite series, like the function $q(x,y,t)$ in this answer. I think you can find it also in the book Brownian Motion and Stochastic Calculus by Karatzas & Shreve. Commented Sep 18, 2022 at 13:03
• Thanks for your comment Kurt G. I have used your link to find the solution above, but there is a mistake somewhere. I get a wrong sign and normalization constant in the expression of $f_{\tau}(t)$.
– Max
Commented Sep 22, 2022 at 11:19
• I have now found all mistakes and posted the answer below. Thanks again for your link Kurt it helped alot! :)
– Max
Commented Sep 23, 2022 at 10:15

We use separation of variables to find a solution we define $$u(x,t) = f(x)g(t)$$ then we get the two equations: \begin{align*} \begin{cases} D\partial_x^2 f - \mu \partial_x f + \lambda f = 0 \\ \partial_t g = - \lambda g \end{cases} \end{align*} where we have introduced the constant $$-\lambda$$ and where $$D=σ^2/2$$. We start by solving the time dependent equation, this equation has the general solution: \begin{align*} g(t) = g(t_0) \exp\big(-\lambda(t-t_0) \big) \ . \end{align*} For the spatial dependent solution $$f$$, we put $$f(x) = \exp(rx)$$ then we get: \begin{align*} D r^2 \exp(rx) - \mu r \exp(rx) + \lambda \exp(rx) = 0 \ . \end{align*} We divide by $$\exp(rx)$$ and solve the corresponding second order polynomial: \begin{align*} r_1 &= \frac{\mu + \sqrt{\mu^2 - 4 D \lambda}}{2 D} \ , \\ r_2 &= \frac{\mu - \sqrt{\mu^2 - 4 D \lambda}}{2 D} \ . \end{align*} Now if $$\mu^2 - 4 D \lambda>0$$ and $$r_1 \neq r_2$$ is not equal then we have the solution: \begin{align*} f(x) = A \exp(r_1 x) + B \exp(r_2 x) \ . \end{align*} If $$\mu^2 - 4 D \lambda=0$$ then $$r_1 = r_2$$ and we then have the solution: \begin{align*} f(x) = A \exp(r x) + B x\exp(r x) \ . \end{align*} If $$\mu^2 - 4 D \lambda<0$$ then the roots are complex, thus we have $$r_1 = \frac{\mu + i \sqrt{4 D \lambda - \mu^2}}{2D}$$, $$r_2 = \frac{\mu - i \sqrt{4 D \lambda - \mu^2}}{2D}$$ and $$\beta = \frac{\sqrt{4 D \lambda - \mu^2}}{2D}$$. The solution is: \begin{align*} f(x) = \exp\left(\frac{\mu x}{2D}\right) \left(A \sin(\beta x) + B \cos(\beta x) \right) \ . \end{align*} Using the boundary condition only the last expression for $$f$$ will give us a non-trivial solution. Thus we have: \begin{align*} f(a) &= \exp\left(\frac{\mu a}{2D}\right) \left(A \sin(\beta a) + B \cos(\beta a) \right) = 0 \ \Rightarrow \ A = - B \frac{\cos(\beta a)}{ \sin(\beta a)} \ . \end{align*} Using the other boundary point we get: \begin{align*} f(b) &= \exp\left(\frac{\mu b}{2D}\right) \left( - B \frac{\cos(\beta a)}{ \sin(\beta a)} \sin(\beta b) + B \cos(\beta b) \right) \ . \end{align*} Form this we conclude that: \begin{align*} B \left ( \cos(\beta b) - \frac{\cos(\beta a)}{ \sin(\beta a)} \sin(\beta b) \right) = 0 \ . \end{align*} If $$B$$ is zero we get a trivial solution, thus $$B$$ is not zero and then we conclude that $$\cos(\beta b) - \frac{\cos(\beta a)}{ \sin(\beta a)} \sin(\beta b) = 0$$. By multiplication with $$-\sin(\beta a)$$ and using the following formula: \begin{align*} \sin(\beta b)\cos(\beta a) - \cos(\beta b) \sin(\beta a) = \sin(\beta ( b - a) ) \ , \end{align*} we see that: \begin{align*} \sin(\beta ( b - a) ) = 0 \ . \end{align*} This is true for $$\beta ( b - a) = n \pi$$ where $$n \in \mathbb{Z}$$. By the formula $$\beta = \frac{\sqrt{4 D \lambda - \mu^2}}{2D}$$ we see that there is not just one possible constant $$\lambda$$, but one for every $$n \in \mathbb{Z}$$: \begin{align*} \lambda_n = D \left( \left( \frac{n \pi}{b-a} \right)^2 + \left( \frac{\mu}{2D} \right)^2 \right) \ , \ \ n \in \mathbb{Z} \ . \end{align*} Where we have introduced a subscript $$n$$ to the possible $$\lambda$$ constants to separate them. Now by the linearity of the equation we get a Fourier series as a solution: \begin{align*} u(x,t) = \exp\left(\frac{\mu x}{2D}\right) \sum_{n \in \mathbb{Z}} \left(A_n \sin\left(\frac{n \pi x}{b-a} \right) + B_n \cos\left(\frac{n \pi x}{b-a} \right) \right) \exp\big(-\lambda_n(t-t_0) \big) \ . \end{align*} Where $$g(t_0)$$ has been left out since it can be absorbed into the constants $$A_n$$ and $$B_n$$. Now because of the odd and even properties of sine and cosine we see that: \begin{align*} A_n \sin\left(\frac{n \pi x}{b-a} \right) + A_{-n} \sin\left(\frac{-n \pi x}{b-a} \right) &= (A_n - A_{-n}) \sin\left(\frac{n \pi x}{b-a} \right) \ , \ n \in \mathbb{N} \ , \\ B_n \cos\left(\frac{n \pi x}{b-a} \right) + B_{-n} \cos\left(\frac{-n \pi x}{b-a} \right) &= (B_n + B_{-n}) \cos\left(\frac{n \pi x}{b-a} \right) \ , \ n \in \mathbb{N} \ . \end{align*} Because of these formulas and the squared $$n^2$$ in the expression of $$\lambda_n$$ we can limit our sum to $$n \in \mathbb{N}_0$$ thus: \begin{align*} u(x,t) = \exp\left(\frac{\mu x}{2D}\right) \sum_{n \in \mathbb{N}_0} \left(A_n \sin\left(\frac{n \pi x}{b-a} \right) + B_n \cos\left(\frac{n \pi x}{b-a} \right) \right) \exp\big(-\lambda_n(t-t_0) \big) \ . \end{align*} One can show the following for sine and cosine: \begin{align*} \int_{a}^b \cos\left(\frac{n \pi x}{b-a} \right) \sin\left(\frac{m \pi x}{b-a} \right) dx &= 0 \ , \ \forall n,m \ . \\ \int_a^b \cos\left(\frac{n \pi x}{b-a} \right) \cos\left(\frac{n \pi x}{b-a} \right) dx &= 0 \ , \ n \neq m \ . \\ \int_a^b \sin\left(\frac{n \pi x}{b-a} \right) \sin\left(\frac{n \pi x}{b-a} \right) dx &= 0 \ , \ n \neq m \ . \\ \end{align*} For the two last expressions, in the case $$m=n$$, we have: \begin{align*} \int_a^b \cos^2\left(\frac{n \pi x}{b-a} \right) dx &= \frac{b-a}{2} \ , \\ \int_a^b \sin^2\left(\frac{n \pi x}{b-a} \right) dx &= \frac{b-a}{2} \ . \end{align*} From these formulas one can deduce that: \begin{align*} B_n = \frac{2}{b-a} \exp\big(\lambda_n(t-t_0) \big) \int_a^b u(x,t) \cos\left(\frac{n \pi x}{b-a} \right) \exp\left(\frac{-\mu x}{2D}\right) dx \ , \\ A_n = \frac{2}{b-a} \exp\big(\lambda_n(t-t_0) \big) \int_a^b u(x,t) \sin\left(\frac{n \pi x}{b-a} \right) \exp\left(\frac{-\mu x}{2D}\right) dx \ . \end{align*} Since this has to be true for all $$t \geq t_0$$ we can put $$t=t_0$$ then we get: \begin{align*} B_n = \frac{2}{b-a} \int_a^b \delta(x-x_0) \cos\left(\frac{n \pi x}{b-a} \right) \exp\left(\frac{-\mu x}{2D}\right) dx = \frac{2}{b-a} \cos\left(\frac{n \pi x_0}{b-a} \right) \exp\left(\frac{-\mu x_0}{2D}\right) \ , \\ A_n = \frac{2}{b-a} \int_a^b \delta(x-x_0) \sin\left(\frac{n \pi x}{b-a} \right) \exp\left(\frac{-\mu x}{2D}\right) dx = \frac{2}{b-a} \sin\left(\frac{n \pi x_0}{b-a} \right) \exp\left(\frac{-\mu x_0}{2D}\right) \ . \end{align*} Inserting this back into our solution we get: \begin{align*} u(x,t) = \frac{2}{b-a} \exp\left(\frac{\mu(x-x_0)}{2D}\right) \sum_{n \in \mathbb{N}_0} \cos\left(\frac{n \pi (x-x_0)}{b-a} \right) \exp\big(-\lambda_n(t-t_0) \big) \end{align*} We now define the survival probability as: \begin{align*} S(t) = \int_a^b u(x,t) dx \ . \end{align*} We see that: \begin{align*} S(t) &= \frac{2}{b-a} \sum_{n \in \mathbb{N}_0} \exp\big(-\lambda_n(t-t_0) \big) \int_a^b \exp\left( \frac{\mu (x-x_0)}{2D} \right) \cos\left( \frac{n \pi (x-x_0)}{b-a} \right) dx \ . \end{align*} Calculating the integral we have: \begin{align*} \int_a^b \exp\left( \frac{\mu (x-x_0)}{2D} \right) &\cos\left( \frac{n \pi (x-x_0)}{b-a} \right) dx = \\ \frac{D}{\lambda_n} \Bigg( \frac{\mu}{2D} & \left( \exp\left(\frac{\mu(b-x_0)}{2D} \right) \cos\left( \frac{n \pi (b-x_0)}{b-a} \right) - \exp\left(\frac{\mu(a-x_0)}{2D} \right) \cos\left( \frac{n \pi (a-x_0)}{b-a} \right) \right) \\ + \frac{n \pi}{b-a} & \left( \exp\left(\frac{\mu(b-x_0)}{2D} \right) \sin\left( \frac{n \pi (b-x_0)}{b-a} \right) - \exp\left(\frac{\mu(a-x_0)}{2D} \right) \sin\left( \frac{n \pi (a-x_0)}{b-a} \right) \right) \Bigg)\ . \end{align*} By inserting this in the expression for $$S(t)$$ we get: \begin{align*} S(t) &= \frac{2}{b-a} \exp\left(\frac{\mu(b-x_0)}{2D} \right) \sum_{n \in \mathbb{N}_0} \frac{D}{\lambda_n} \frac{\mu}{2D} \cos\left( \frac{n \pi (b-x_0)}{b-a} \right) \exp\big(-\lambda_n(t-t_0) \big) \\ &- \frac{2}{b-a} \exp\left(\frac{\mu(a-x_0)}{2D} \right) \sum_{n \in \mathbb{N}_0} \frac{D}{\lambda_n} \frac{\mu}{2D} \cos\left( \frac{n \pi (a-x_0)}{b-a} \right) \exp\big(-\lambda_n(t-t_0) \big) \\ &+ \frac{2}{b-a} \exp\left(\frac{\mu(b-x_0)}{2D} \right) \sum_{n \in \mathbb{N}_0} \frac{D}{\lambda_n} \frac{n \pi}{b-a} \sin\left( \frac{n \pi (b-x_0)}{b-a} \right) \exp\big(-\lambda_n(t-t_0) \big) \\ &- \frac{2}{b-a} \exp\left(\frac{\mu(a-x_0)}{2D} \right) \sum_{n \in \mathbb{N}_0} \frac{D}{\lambda_n} \frac{n \pi}{b-a} \sin\left( \frac{n \pi (a-x_0)}{b-a} \right) \exp\big(-\lambda_n(t-t_0) \big) \ . \end{align*} Now the first hitting time density $$f_{\tau}(t)$$ can be found by the expression: \begin{align*} f_{\tau}(t) = - \frac{dS}{dt} \ . \end{align*} $$f_\tau$$ is also known as the first passage time density. we see that: \begin{align*} f_{\tau}(t) &= \frac{2}{b-a} \exp\left(\frac{\mu(b-x_0)}{2D} \right) \sum_{n \in \mathbb{N}_0} \frac{\mu}{2} \cos\left( \frac{n \pi (b-x_0)}{b-a} \right) \exp\big(-\lambda_n(t-t_0) \big) \\ &- \frac{2}{b-a} \exp\left(\frac{\mu(a-x_0)}{2D} \right) \sum_{n \in \mathbb{N}_0} \frac{\mu}{2} \cos\left( \frac{n \pi (a-x_0)}{b-a} \right) \exp\big(-\lambda_n(t-t_0) \big) \\ &+ \frac{2}{b-a} \exp\left(\frac{\mu(b-x_0)}{2D} \right) \sum_{n \in \mathbb{N}_0} D \frac{n \pi}{b-a} \sin\left( \frac{n \pi (b-x_0)}{b-a} \right) \exp\big(-\lambda_n(t-t_0) \big) \\ &- \frac{2}{b-a} \exp\left(\frac{\mu(a-x_0)}{2D} \right) \sum_{n \in \mathbb{N}_0} D \frac{n \pi}{b-a} \sin\left( \frac{n \pi (a-x_0)}{b-a} \right) \exp\big(-\lambda_n(t-t_0) \big) \ . \end{align*} For convenience we define: \begin{align*} v(x,t) = \frac{2}{b-a} \exp\left(\frac{\mu(x-x_0)}{2D} \right) \sum_{n \in \mathbb{N}_0} \frac{n \pi}{b-a} \sin\left( \frac{n \pi (x-x_0)}{b-a} \right) \exp\big(-\lambda_n(t-t_0) \big) \ . \end{align*} Then we have: \begin{align*} f_{\tau}(t) &= \frac{\mu}{2} \Big( u(b,t) - u(a,t) \Big) + D \Big( v(b,t) - v(a,t) \Big) \ . \end{align*} By our boundary condition $$u(b,t) = u(a,t) = 0$$ thus: \begin{align*} f_{\tau}(t) &= D \Big( v(b,t) - v(a,t) \Big) \ . \end{align*}
The following plots are for $$a=-1$$, $$b=1$$, $$x_0=0$$ and $$t_0 = 0$$. In the plots we have $$u:=\mu$$.