How to evaluate the limit of multifactorial $\lim_{n\to 0} \sqrt[n]{n!!!!\cdots !}$

Question

It is well known that $\displaystyle \lim_{n\to \infty}\sqrt[n]{n!}=\infty$, however, if we let $n\to 0$ we have a different result with a beautiful combination of $e$ and $\gamma$, that is

$$\lim_{n\to 0}\sqrt[n]{n!}= e^{-\gamma}\tag{1}\label{1}$$

To prove result \eqref{1}, we observe that the limit attains the form of $1^{\infty}$ so we can write it as $$\lim_{n\to 0} \exp\left(\frac{\ln\Gamma(n+1)}{n}\right)\underbrace{=}_{\text{L'Hopital's rule}}\lim_{n\to 0} e^{\Gamma'(n+1)}=e^{\psi_0(1)}= e^{-\gamma}$$

Now I wish to know the limit of the following multifactorial form for $k\in\mathbb {Z^+}$

$$\lim_{n\to 0}\sqrt[n]{n\smash[b]{\underbrace{!! !!\cdots !}_{k}}}={?}\tag{2}\label{2}\\$$

For $k=1$ we are done above and for $k=2$ we get the limit $\sqrt{2} e^{-\frac{\gamma}{2}}$. To prove this we use the double factorial argument (see equation (5)) \begin{align}\lim_{n\to 0}\sqrt[n]{n!!}&=\lim_{n\to 0} \left(2^{\frac{n}{2}+\frac{1-\cos(\pi n)}{4}}\pi^{\frac{\cos(\pi n)-1}{4}}\Gamma\left(1+\frac{n}{2}\right)\right)^{\frac{1}{n}}\\&=\sqrt{2}\lim_{n\to 0} \sqrt[n]{\Gamma\left(1+\frac{n}{2}\right)}\\&=\sqrt 2\exp\lim_{n\to 0}\left(2^{-1} \Gamma\left(1+\frac{n}{2}\right)\psi_0\left(1+\frac{n}{2}\right)\right)\tag{L'Hopital's rule}\\&=\sqrt{2}e^{-\frac{\gamma}{2}}\end{align}

since for $k=1,2$ we have evaluated the limit. How to evaluate the limit of equation \eqref{2} for all $k>2$?

Answer 1

The main task in this problem is to construct an analytic continuation of the multifactorial function $$n!^{(k)}=n(n-k)(n-2k)\cdots$$ over the real numbers. On the set of integers, we can take advantage of modulo arithmetic to arrive at $$n!^{(k)}_i=i\cdot k^{(n-i)/k}\frac{\Gamma(1+n/k)}{\Gamma(1+i/k)}$$ which is defined whenever $n\equiv i\pmod k$. Therefore, we need to find a function which interpolates the values of $$1,\frac{k^{-1/k}}{\Gamma(1+1/k)},\frac{2k^{-2/k}}{\Gamma(1+2/k)},\cdots,\frac{(k-1)k^{-(k-1)/k}}{\Gamma(1+(k-1)/k)}$$ so we consider the ansatz $$f(x)=\prod_{i=1}^{k-1}\left(\frac{ik^{-i/k}}{\Gamma(1+i/k)}\right)^{g(x)}$$ where $g(x)\equiv1$ when $x-i\equiv0\pmod k$ and is zero otherwise. This suggests the analytic function $$\frac{\sin(\pi(x-i))}{k\sin(\pi(x-i)/k)}$$ which is unfortunately negative when $(x-i)/k$ is an odd integer. We can circumvent this by introducing a factor of $\cos(\pi(x-i)/k)$ as its sign agrees with that of the above function when $x-i\equiv0\pmod k$. It follows that we can extend the multifactorial function to the reals through $$x!^{(k)}=k^{x/k}\Gamma\left(1+\frac xk\right)\prod_{i=1}^{k-1}\left(\frac{ik^{-i/k}}{\Gamma(1+i/k)}\right)^{\sin(\pi(x-i))\cot(\pi(x-i)/k)/k}\tag1$$ It is worth noting that this is not unique, as we can multiply by a factor of $\cos(j\pi(x-i)/k)$ for some integer $j$ — here we just took $j=1$. We can now determine the limit by considering each term separately \begin{align}\lim_{x\to0}x!^{(k)/x}&=k^{1/k}\lim_{x\to0}\exp\left(\frac{\log\Gamma(1+x/k)}x\right)\lim_{x\to0}\prod_{i=1}^{k-1}\left(\frac{ik^{-i/k}}{\Gamma(1+i/k)}\right)^{\sin(\pi(x-i))\cot(\pi(x-i)/k)/kx}\\&=k^{1/k}e^{-\gamma/k}\prod_{i=1}^{k-1}\left(\frac{ik^{-i/k}}{\Gamma(1+i/k)}\right)^{h(k)}\end{align} where \begin{align}h(k)&=\lim\limits_{x\to0}\frac{\sin(\pi(x-i))\cot(\pi(x-i)/k)}{kx}=-\cot\frac{\pi i}k\lim\limits_{x\to0}\frac{\sin\pi x\cos\pi i}{kx}\\&=-(-1)^i\frac\pi k\cot\frac{\pi i}k.\end{align} Using this definition of $g(x)$, the limit evaluates to $$\lim_{x\to0}x!^{(k)/x}=\left[\frac k{e^\gamma}\prod_{i=1}^{k-1}\left(\frac{\Gamma(1+i/k)}{ik^{-i/k}}\right)^{\pi(-1)^i\cot\frac{\pi i}k}\right]^{1/k}.$$ Note that the multifactorial function in $(1)$ can be extended onto the complex plane. It is holomorphic everywhere except at negative integer multiples of $k$, similar to the gamma function.

Answer 2

I tried with this definition of multifactorial

The calculation comes quite quickly, I'll just give you the final result

$$\lim_{z\to 0}z!_{(\alpha)}^{\frac{1}{z}}=\alpha^{\frac{1}{\alpha}}e^{-\frac{\gamma}{\alpha}}\prod_{j=1}^{\alpha-1}\left(\frac{\alpha^{\frac{\alpha-j}{\alpha}}}{\Gamma\left(\frac{j}{\alpha}\right)}\right)^{w\left(j\right)}$$

Where:

$$s(j)=\frac{4\pi}{a^{2}}\sum_{k=1}^{\left\lfloor\frac{\alpha-1}{2}\right\rfloor}k\sin\left(\frac{2k\pi}{\alpha}j\right)$$

For example:

if $\alpha=3$ the limit is: $\displaystyle 3^{\frac{1}{3}}e^{-\frac{\gamma}{3}}\left(\frac{3^{\frac{1}{3}}\Gamma\left(\frac{2}{3}\right)}{\Gamma\left(\frac{1}{3}\right)}\right)^{\frac{2\sqrt{3}\pi}{9}}$