In dimensions 1,3 we have that for $T=\int_{[0,\infty)}1_{B_{t}\in B(0,1)}dt$ has the Laplace transform $$E[e^{-\lambda T}]=sech(\sqrt{2\lambda}).$$

Inverting this in Mathematica didn't give a clean answer, but I found in "An Atlas of Functions: with Equator, the Atlas Function Calculator" that

$$\int_{a-i\infty}^{a+i\infty}\frac{sech(\nu \sqrt{\lambda})}{\sqrt{\lambda}}e^{\lambda s}\frac{ds}{2\pi i}=\frac{1}{\nu}\hat{\theta}_{2}(\frac{1}{2},\frac{t}{\nu^{2}}).$$

Has there been any further work on relating this to behaviour for Brownian motion?