# Factorial of infinity

$$\zeta'(0) = \log\sqrt\frac 1 {2\pi}$$

And that: $$e^{-\zeta'(0)} = 1\cdot2\cdot3\cdots \infty = \infty! = \sqrt{2\pi}$$

I found this result very strange and tried to simplify if but i couldn't. Can somebody help me?

Thanks!

• I suggest to start with simpler problems: $$\zeta(0)=1+1+1+\ldots=-\frac12,$$ and also $$\zeta(-1)=1+2+3+\ldots+\infty=\frac{\infty(\infty+1)}{2}=-\frac{1}{12}.$$ – Start wearing purple May 13 '13 at 14:50
• In the linked page there is (at least) one wrong value. It is stated that $\eta(1)=\infty$. However, it is $\eta(1)=\log(2)$. I didn't check the other entries around this... – Gottfried Helms May 13 '13 at 15:06
• A related question. – J. M. isn't a mathematician May 13 '13 at 17:26

Let me try to explain what is going on with a simpler example. Consider geometric series $$F(s)=\sum_{n=0}^{\infty}s^n,$$ which converges inside the unit circle $|s|<1$ to the function $F(s)=\frac{1}{1-s}$. Note, however, that $F(s)$ can be analytically continued to the whole complex plane with excluded point $s=1$, which is a simple pole of $F(s)$. This point was in fact the origin of radius of convergence equal to $1$. So we are tempted to write identities like $$F(2)=1+2+4+8+\ldots=\frac{1}{1-2}=-1.$$ It is important to understand that they do not "really" hold. The logic is the following: the initial series defines some function in a "small" domain, but then sometimes this function can be (uniquely!) continued to a larger domain where the original series no longer makes sense. In general, this continuation will have singularities.

Sometimes - actually, almost always - one cannot construct a continuation to the whole complex plane with isolated singular points because of specific singularities called natural boundaries. An example of this is encountered if we consider instead of $F(s)$ the series $$G(s)=\sum_{n=0}^{\infty}s^{n^2},$$ which converges inside the unit circle even better than $F(s)$ but cannot be continued elsewhere: $|s|=1$ is a dense set of singularities for $G(s)$.

Now let us go back to initial question. If we consider the series $$\zeta(s)=\sum_{n=1}^{\infty}n^{-s},$$ it converges to some function $\zeta(s)$ in the halfplane $\mathrm{Re}\,s>1$. It turns out that this function can be (again, uniquely!) continued to the whole complex plane except the point $s=1$ where it has simple pole, similarly to our previous example. This continued function, as well as its derivatives, can be evaluated for any $s\neq1$ (for example, at $s=-1$ and $0$). One of the ways to do it is to systematically use integral representation $$\Gamma(s)\zeta(s)=\int_0^{\infty}\frac{x^{s-1}dx}{e^x-1},$$ which for $\mathrm{Re}\,s>1$ reproduces the original series definition.

• I suppose one could be charitable and say that the result in the OP is the Dirichlet regularization of the factorial. In related news, the "product" of the prime numbers is $4\pi^2$. – J. M. isn't a mathematician May 13 '13 at 17:25

One way to crosscheck (not prove/disprove (!)) that results approximative and numerically is to use the alternating zeta at the same parameter, which can be summed using simple summation procedures for divergent summation. If we write $$\eta(s) = \frac 1{1^s} - \frac 1{2^s} + \frac 1{3^s} - \cdots$$ and furtherly $$\eta(s) = \exp(-s\ln 1) - \exp(- s\ln2) + \exp(-s\ln 3) - \cdots$$

then we can write the derivative as $$\eta(s)' = (-\ln1)\cdot \exp(-s\ln 1) - (-\ln2)\cdot \exp(- s\ln2) + (-\ln3)\cdot \exp(-s\ln 3) - \cdots$$ and letting $\lim s \to 0$ this gives $$\begin{eqnarray} \eta(0)' &=& (-\ln1) - (-\ln2) + (-\ln3) - \cdots \\ \eta(0)' &=& - (\ln1 - \ln2 + \ln3 - \cdots) \tag 1 \end{eqnarray}$$ The latter can be handled by a simple procedure for divergent summation; for instance Euler-summation $\mathfrak E$ which gives then $$\eta(0)' \underset{\mathfrak E}{=} 0.225791352645\ldots$$ Now in (1) we see the alternating signs but for the equivalent value of the zeta-function this series is non-alternating. So $$\zeta(0)' = - (\ln1 + \ln2 + \ln3 + \cdots) \tag 2$$ We assume now, that the value $\zeta(0)' = s$ exists and that we can handle it into the divergent sum. Here we allude to Euler's trick of the transformation-formula between zeta() and eta() and exploit that the logs at even arguments are negative in the eta-formula. But $\ln(2m) = \ln2 + \ln m$ and thus we might write $$\begin{eqnarray} s &=& - ( &(\ln 1 + \ln 2 + \ln 3 + \ln 4 + \ln 5 + \ln 6 + \ln 7\cdots) \\ s &=& - ( &(\ln 1 - \ln 2 + \ln 3 - \ln 4 + \ln 5 - \ln 6 + \ln 7\cdots) \\ & & &+2 \cdot ( \ln 2 + \ln 4 + \ln 6 + \ln 8 + \cdots) ) \\ s&=& &\eta(0)' - 2 \cdot ( (\ln2 + \ln1) + (\ln2 + \ln 2) + (\ln 2 + \ln 3) + \cdots)) \\ s&=& &\eta(0)' - 2 \cdot ( \zeta(0) \cdot \ln 2 - s ) \tag 3\\ \end{eqnarray}$$ where we have accumulated the constant $\ln 2$ which occurs as summand at each consecutive logarithm in the last equation (3) as $\zeta(0) \cdot \ln 2$. Then we have a determination for $s$ by $$\begin{eqnarray} &s&=& \eta(0)' - 2 \cdot ( \zeta(0) \cdot \ln 2 - s ) \\ &-s&=& \eta(0)' - 2 \cdot \zeta(0) \cdot \ln 2 \\ \zeta(0)'&=s&=& -\eta(0)' - \ln 2 \tag 4 \end{eqnarray}$$ $$\begin{eqnarray} &\zeta(0)' &\sim - 0.225791352645\ldots - 0.693147180560 \\ && \sim -0.918938533205\ldots \tag 5 \\ \end{eqnarray}$$ And indeed, numerically we get using Pari/GP

log(1/sqrt(2*Pi))
%393 = -0.918938533205


and also

h=1e-12;(zeta(0+h/2)-zeta(0-h/2))/h
%394 = -0.918938533205


This all required assumptions about the general summability of the divergent terms. But there is a fairly well developed theory of divergent summation, from which one may take the confirmation about the acceptability of the above transformations. They are intended to show one path to obtain the result; perhaps/likely there are others more elegant or even rigid ones.

Appendix: The Pari/GP code for the numerical evaluation of $\eta(0)'$ in (3) is

-sumalt(k=0,(-1)^k*log(1+k))
%387 = 0.225791352645


I am riemannium, author of the above quoted entry and the blog The Spectrum Of Riemannium (TSOR). I have just corrected the typos in my article quoted above. I mean, I wrote $\eta (1)=\log 2$. It was a typo while I was writing the article. I think there are no additional mistakes but if you find more "errors", let me know. I am taking my blog as a serious specialized "job", so I don't like to input errors. I am also a perfectionist. I also thank you for the links so some pretty references about regularized infinite products I did not know till now! :).

• You're welcome. I couldn't post a comment on your site (don't have wordpress-id or such and didn't find a email-adress) – Gottfried Helms May 14 '13 at 15:15
• I think you should be able to post as anonymous...At current time I don't like to release my e-mail there. I don't want more spam or face the same trolls I had to erase by only admitting posts that I can read previously. A blog like mine is young and not as "big" yet...But I think I have done good posts about some topics...The only "but" is that the inclusion of Mathematics sometimes can scare people, but it is a payment I am decided to do. The Universe/Polyverse is written in mathematical tongue, move away from my blog if you dislike it :P. – riemannium May 14 '13 at 15:24
• What is interesting about the factorial of infinity is that it is the same number that the gaussian integral $$\int_{-\infty}^\infty e^{-x^2/2}dx=\sqrt{2\pi}=\infty !$$ ;) – riemannium May 15 '13 at 21:49

My method is the following: $$\prod_{n=1}^\infty(n)= (2\pi)^{0.5}$$ $$\prod_{n=1}^\infty(n)= e^{\sum_{n=1}^\infty \ln(n)}$$ First we write it down in it's alternating form. This is done in order to make it Cesàro sumable. Since the periode doesn't matter i take a period of $2$. $$\sum_{n=1}^\infty \ln(n)= \sum_{n=1}^\infty 2\ln(2n/4)-\ln(n/4) =\sum_{n=1}^\infty (-1)^n\ln(n/4)$$ $$\sum_{n=1}^\infty -(-1)^n\ln(4)=1/2 (\ln(4))$$

$$\sum_{n=1}^\infty (-1)^n\ln(n)=$$ $$\frac 1 2 \left(\sum_{n=1}^{1/2(m)} \ln(2n)-\sum_{n=1}^{1/2(m)} \ln(2n-1) + \sum_{n=1}^{1/2(m-1)} \ln(2n) -\sum_{n=1}^{1/2(m+1)} \ln(2n-1) \right)=$$ $$(2m-1)/2\ln(2)+\ln((1/2\cdot m)!)+\ln((1/2(m-1))!)-\ln(m!)=$$ $$-\frac 1 2 \ln(2)+\ln\left(\frac{((1/2)m)!\cdot((1/2\cdot(m-1))!)\cdot2^m}{m!}\right) = \frac 1 2 \ln\frac \pi 2$$

So here you go: $$\prod_{n=1}^\infty(n)= (2\pi)^{0.5}$$