I'm curious if there is a known asymptotic scaling for the return-to-origin (i.e. recurrence) probability for a random on $Z^d$ as a function of $d$? 


Mathworld gives the recurrence probability:
<http://mathworld.wolfram.com/PolyasRandomWalkConstants.html>


$p(d) = 1-\frac{1}{u(d)} = 1- \left(\int\_{t=0}^{\infty}I\_0(\frac{t}{d})^de^{-t} dt\right)^{-1}$


Where $I\_0$ is a modified Bessel function of the first kind: <http://mathworld.wolfram.com/ModifiedBesselFunctionoftheFirstKind.html>


Taking a wild guess, I found that a quartic fit seems to work well for values of $p(d)$ for $d = 1$ to $8$: <http://www.wolframalpha.com/input/?i=quartic+fit+%7B3%2C+0.340537%7D%2C%7B4%2C+0.193206%7D%2C%7B5%2C+0.135178%7D%2C%7B6%2C+0.104715%7D%2C%7B7%2C+0.0858449%7D%2C%7B8%2C+0.0729126%7D>


Also, what's the best way to approximate the integral expression necessary to calculate $p(d)$ if we only need $k$ digits of precision?




---


Sorry folks: [Pólya's Random Walk Constants at infinity](https://mathoverflow.net/questions/83317/polyas-random-walk-constants-at-infinity)


I should have done a more thorough job checking for previous questions about this topic. =(



Is there a known asymptotic scaling for the probability of recurrence for a walk on $Z^d$?
Asymptotic in the sense $d \to \infty$? From Maple I get:
$$
p(d) = \frac{1}{2 d} + \frac{1}{2 d^{2}} + \frac{7}{8 d^{3}} + \frac{35}{16 d^{4}} + O \Bigl(d^{-5}\Bigr)
$$
as $d \to \infty$.


