6

I was reading about Glaisher-Kinkelin Constant and came across the following formulas $$ \sqrt{2\pi} = \lim_{n\rightarrow\infty} \frac{n!}{n^{\frac{2n+1}{2}}e^{-n}} $$ $$ A = \lim_{n\rightarrow\infty} \frac{H(n)}{n^{\frac{6n^2+6n+1}{12}}e^{-\frac{n^2}{4}}} $$

  • I came across the Stirling numbers of the second kind $S(p,k)$ in research, in the form of $a_{p,k} = k! S(p,k)$.
  • Now in the formulas above, we see that there are the polynomials $2n+1$ and $6n^2+6n+1$, whose coefficients are $a_{2,\cdot}$ and $a_{3,\cdot}$ respectively.
  • Is there a known constant that uses the polynomial $24n^3+36n^2+14n+1$, the coefficients for $a_{4,\cdot}\ ?$.
  • My guess is that the numerator in the limit would have something like $H_2(n) = \prod k^{k^k}$.
Felix Marin
  • 89,464
gist076923
  • 1,266

1 Answers1

2

Okay so after a lot of complicated math (Euler-Maclaurin formula, as well as some umbral calculus), I can answer my own question kinda. So suppose you want $C,p,q,r$ s.t. $$ C = \lim_{n\to\infty} \frac{\prod_{k=1}^n k^{p(k)}}{n^{q(n)}e^{r(n)}} $$

Instead of looking at $a_{p,k}$, it turns out Bernoulli polynomials occur pretty naturally as $q$'s

So we want to estimate $w(n) = \prod_{k=1}^n k^{p(k)}$

Taking logs, we get $\ln w(n) = \sum_{k=1}^n p(k) \ln(k)$

Letting $\phi(x) = p(x) \ln x$, Euler Maclaurin then says $$ \ln w(n) = \phi(1) + \int_1^n \phi(x) dx + \sum_{k=1}^p \frac{B_k^+}{k!} \left(\phi^{(k-1)}(n) - \phi^{(k-1)}(1)\right) + R_p $$

Now note we want something of the form $w(n) = c_nn^{q(n)}e^{r(n)}$, and so the $q(n)$ term arises out of the $\ln n$ terms of $\ln w(n)$.

Now by integration by parts, $$\begin{align*} &\int_1^n p(x) \ln(x) dx \\ &= p^{(-1)}(n) \ln n - \int_1^n \frac{p^{(-1)}(x)}{x} dx \end{align*}$$

Note we assume the constant term of $p^{(-1)}$ to be zero, so that the right integral does not produce any $\ln$'s

Also, if $f(x) = h(x) \ln(x)$ for some polynomial $h$, we have $$f'(x) = h'(x) \ln(x) + h(x)/x$$

Noticing that $\phi^{(|p|)}$ is the last derivative with any $\ln n$ terms, we obtain $$q(n) = \sum_{n=0}^{|p|+1} \frac{B_k^+}{k!} p^{(k-1)}(n)$$

Defining $P = p^{(-1)}$, and using the umbral calculus linear operator defined by $L(z^k) = B_k^-$, we obtain

$$ \begin{align*} &q(n) \\ &=\sum_{n=0}^{|P|} \frac{B_k(1)}{k!} P^{(k)}(n) \\ &=\sum \frac{L((z+1)^k)}{k!} P^{(k)}(n) \\ &=L \left( \sum \frac{(z+1)^k}{k!} P^{(k)}(n)\right) \\ &=L(P(z+n+1)) \\ &=\sum_{k=0}^{|P|} P_k L((z+n+1)^k) \\ &=\sum_{k=0}^{|P|} P_k B_k(n+1) \end{align*} $$

This implies that if $q(n) = B_k(n+1)$, we obtain $P(n) = n^k$, thus $p(n) = kn^{k-1}$

gist076923
  • 1,266