# Why does the Gamma-function complete the Riemann Zeta function?

Defining $$\xi(s) := \pi^{-s/2}\ \Gamma\left(\frac{s}{2}\right)\ \zeta(s)$$ yields $\xi(s) = \xi(1 - s)$ (where $\zeta$ is the Riemann Zeta function).

Is there any conceptual explanation - or intuition, even if it cannot be made into a proof - for this? Why of all functions does one have to put the Gamma-function there?

Whoever did this first probably had some reason to try out the Gamma-function. What was it?

(Best case scenario) Is there some uniform way of producing a factor out of a norm on the rationals which yields the other factors for the p-adic norms and the Gamma factor for the absolute value?

-
Have you ever read Emil Artin's monograph about the gamma function? – Harry Gindi Dec 3 '09 at 14:20
They are both conceptually related to sums of powers. The $\zeta$ function itself is defined as a non-alternating sum of powers for $\Re(z)>1$, and as an alternating sum of powers (times a certain factor) for $\Re(x)\in(0,1)$ On the other hand, geometric shapes of the form $x^n+y^m=1$, called superellipses or Lame curves, are also bounded sums of powers. But by integrating $y=\sqrt[m]{1-x^n}$ or $x=\sqrt[n]{1-y^m}$ on $(0,1)$ we get the multiplicative inverse of the binomial coefficient ${m+n\choose n}={m+n\choose m}$, which is obviously expressible in terms of the $\Gamma$ function. – Lucian Jun 1 '14 at 16:52

To the best of my understanding, the answer is yes, and this uniform way consists of some integration over the local field. This is explained in John Tate's dissertation. One starts with a certain smooth rapidly decreasing function, for which one takes the characteristic function of the p-adic integers in the nonarchimedean case and the function $e^{-|x|^2}$ for an archimedean field. This is being multiplied with $|x|^s$ (approximately) and integrated over the Haar measure of the additive group of the field. This produces the $\Gamma$-factor for an archimedean field and $(1-p^{-s})^{-1}$ for a p-adic field.

-
That is the content of John Tate's PhD thesis. It is in Cassels & Froehlich - Algebraic Number Theory (Last chapter). In my opinion, this is still the best reference for this matter. – Marc Palm Feb 9 '11 at 10:45
The factor $\pi ^{-s/2}\Gamma (s/2)$ is the Mellin transform of $e^{-\pi x^2}$ which is its own Fourier transform. This yields a functional equation for the theta function whose Mellin transform is the Zeta function (together with the Gamma factor); therefore, this zeta function gets a functional equation. – Venkataramana Oct 25 '13 at 14:38
The expression $\pi^{-s/2} \Gamma(s/2)$ looks similar to the volume of the $s$-dimensional unit ball (say $s$ is an integer). Is this a coincidence? – mlbaker Jan 23 '15 at 7:02

One way to get started is to look at the integral for the gamma function: $$\Gamma(s) = \int_0^\infty t^{s-1} e^{-t}\,dt$$ Subsitute $t=nx$ in the integral to arrive at $$\frac{\Gamma(s)}{n^s} = \int_0^\infty e^{-nx}x^{s-1}\,dx$$ which we then sum up to get $$\Gamma(s)\zeta(s)=\int_0^\infty \frac{x^{s-1}}{e^x-1}\,dx$$ which already shows that there is some connection between the gamma and zeta functions, and it does in fact allow us to extend the definition of the zeta function into the critical strip.

What comes next is far less obvious, but the idea is to introduce a branch cut for $x^{s-1}$ along the positive real axis, and to replace the above integral by one running from $+\infty$ along the bottom of the positive real axis, around the origin, and back to $+\infty$ along the top of the real axis. This introduces an extra factor $1-e^{2\pi i s}$. Now start expanding the circle around the origin, taking account of the poles of the integrand along the imaginary axis as we go, and end up with $$\Gamma(s)\zeta(s)=(2\pi)^{s-1}\Gamma(1-s)\sin(\tfrac12\pi s)\zeta(1-s).$$ From there, some cleanup still remains. As I said, this is not terribly intuitive, so it doesn't answer your question, but the first paragraph should at least give you a notion how the gamma and zeta functions are interrelated.

-
So far, this answer has received four upvotes and two downvotes. I am curious about the reason for the downvotes: I thought they were intended for off topic or wrong answers, especially the sort of answers you want to discourage, and I can't see that this answer is either. Perhaps I should have left the second half out of it, since it does not contribute much to the why question, but to me, that doesn't seem sufficient reason for a downvote. – Harald Hanche-Olsen Dec 3 '09 at 15:24
I guess Riemann's opinions about the Riemann zeta function aren't good enough for some people? – Ryan Budney Dec 3 '09 at 22:18
I was one of the downvotes. I'll say, I don't think your answer is wrong or problematic, I just think there's a much better answer, which hasn't been written properly: "Q has a real prime." I'm not familiar enough with the subject to write a good answer like that, but I'm familiar enough with the subject to say I don't think an answer which leaves it out should be at the top. – Ben Webster Dec 3 '09 at 23:14
Ryan- Actually yes. I think Tate understands the Riemann zeta function a lot better than Riemann ever did, though of course, that involved a lot of standing on the shoulders of giants (specifically, Tate knew about class field theory, and Riemann never had a chance). – Ben Webster Dec 3 '09 at 23:18
@Ben Webster: though it is true that it is the real prime of Q that allows for the appearance of the gamma factor, the question was "Why the Gamma function?", not "Why is there another factor?". The real prime could be considered as a reason to have another factor, not a reason for that factor to be the gamma function. Harald's answer illustrates how the gamma function arises in the proof. – Rob Harron Dec 3 '09 at 23:40

As has been explained above, the zeta function has a factor for each completion of $\mathbb{Q}$. The factor at $\mathbb{R}$ has to do with integrating $e^{- \pi x^2}$ and the factor at $\mathbb{Q}_p$ has to do with integrating the characteristic function of $\mathbb{Z}_p$.

Some people might wonder why these two functions were chosen. The answer is simple: they are both their own Fourier transforms.

Also, I don't think anyone has recommended Terry Tao's expository post on this material yet. It is quite good.

-
Ah, Terry Tao's post is clearly the answer I was looking for. – Ben Webster Dec 4 '09 at 22:32
But there is a lot of functions which are their own Fourier transform, isn't it ? why this one in particular ? what would happen with other functions with the same property ? – Simon Henry Oct 26 '15 at 9:59

There may also be some interest in the point of the "local functional equation", namely, that in fact the Gamma function (with the power of $\pi$) is just one (optimized) possibility, and somehow making a suboptimal choice doesn't really matter:

For a Schwartz function $f$ on $\mathbb R$, let $\Gamma(f,s)=\int_{\mathbb R^\times} |x|^s\,f(x)\;{dx\over |x|}$. The usual Gamma factor is obtained by taking a Gaussian. The local functional equation (proven by changing variables in the defining integrals, in the range $0<{\rm Re}(s)<1$, is $$\Gamma(f,s)\cdot \Gamma(\hat{g},1-s)\;=\; \Gamma(\hat{f},1-s)\cdot \Gamma(g,s)$$ for any two Schwartz functions $f,g$. And Riemann's argument proves $$\Gamma(f,s)\cdot \zeta(s) \;=\; \Gamma(\hat{f},1-s)\cdot \zeta(1-s)$$ for any Schwartz $f$.

-
If I'm thinking correctly, $\zeta(s)$ is formally $\Gamma(g,1-s)$ where $g$ is a sum of $\delta$ functions at the positive integers. Is there some sense in which $g = \hat{g}$ here? – David Speyer Sep 10 '14 at 16:47
I guess that would be the Dirac comb. See mathoverflow.net/a/39187/6818 – Vít Tuček Sep 10 '14 at 16:53
For the record, the Hermite functions yield just a polynomial multiple of the standard Gamma factor. For example $H_{16}$ leads to $2027025 + 16 (-1 + s) s (274455 + 2 (-1 + s) s (27051 + 8 (-1 + s) s (127 + (-1 + s) s)))$ – Vít Tuček Sep 10 '14 at 22:46
Starting from $\Gamma(f,s)\zeta(s) = \Gamma(\widehat{f},1-s)\zeta(1-s)$ and multiplying both sides by $\Gamma\left(\frac{s}{2}\right)\Gamma\left(\frac{1-s}{2}\right)$ I obtain, after cancellation of the classical functional equation for $\xi(s)$, the following equation $$\Gamma\left(\frac{1-s}{2}\right)\Gamma\left(\frac{1+s}{2}\right) = -\imath \frac{\pi}{\sin(\pi s)},$$ which can't be true since for real $s$ the left hand side is real whereas the right hand side is imaginary. Where did I make mistakes? – Vít Tuček Sep 11 '14 at 22:00
@VítTuček, ah, one little point is that the game is to integrate over the whole $\mathbb R^\times$, so using odd Schwartz function $xe^{-\pi x^2}$ simply produces $0$, avoiding the seeming paradox. (I'd wager that using the second Hermite polynomial produces no contradiction!) For quadratic imaginary fields, for example, the range of choices of equivariance under the circle action (thinking of Iwasawa-Tate set-up) gives many chances to cancel-and-be-zero, etc. – paul garrett Sep 11 '14 at 22:07

Whoever did this first probably had some reason to try out the Gamma-function. What was it?

The first one to do this was, precisely, Riemann in his famous (and 150 years old) paper: Ueber die Anzahl der Primzahlen unter einer gegebenen Grösse. There he proved the functional equation as well, with the method that Harald explained above.

-
Related to mathoverflow.net/questions/58004/…. – Tom Copeland May 19 '12 at 15:33
Actually, not quite. Define $\eta(z)=(1-2^{1-z})\zeta(z)=\sum_{n\ge 1}\frac{(-1)^{n+1}}{n^z}$. This sum converges (conditionally) when Re$\,z>0$, thus $\eta$ is defined in the same half-plane (modulo considerations for Re$\,z=1$. The functional equation for $\zeta$ leads to a functional equation for $\eta$. The latter makes sense without complex analysis since $\eta(s)$ and $\eta(1-s)$ are both defined if $0<s<1$. This functional equation was already published by Euler! See: 1) E. Landau: Euler und die Funktionalgleichung der Riemannschen Zetafunktion. 2) A. Weil: Prehistory of the zeta-func. – M Mueger Nov 11 '15 at 22:01

"Why of all functions does one have to put the Gamma-function there?"

zeta(s) has trivial zeroes at -2, -4, -6, etc. zeta(1-s) thus has trivial zeroes at s=3, 5, 7, etc - a completely different set of zeroes.

To make a reflection formula where zeta(s) is somehow equal zeta(1-s), you have to get rid of the two differing sets of trivial zeroes. Multiplying by the gamma is perfect for this since its poles will cancel out those zeroes. For example, gamma(s/2) has poles at 0, 2, 4, 6, etc. and should go with zeta(s). gamma((1-s)/2) has poles at s=1, 3, 5, etc. and should go with zeta(1-s).

It's possible to prove that gamma is the right choice, but Euler no doubt discovered that gamma is the right function through numerical experimentation - when he discovered the zeta reflection formula like 250 years ago.

-
Nice, thank you! I can imagine that this was maybe the first hint, then leading to the connection pointed out by Harald Hanche-Olsen... – Peter Arndt Aug 2 '10 at 10:34
Peter: that was not the first hint. That the zeta-function even makes sense for negative numbers (in a rigorous sense) was first worked out by Riemann through his proof of analytic continuation. Before that there was not a known relation between zeta(s) and zeta(1-s) to motivate using the Gamma-function. Although Euler, long before Riemann, had derived a non-rigorous formula that is equivalent to the functional equation of the zeta-function just at integers, I don't think it was something that influenced Riemann's work which brought in the Gamma-function explicitly. – KConrad Mar 10 '11 at 5:45
I posted an MSE question asking about Riemann's thinking on symmetrizing the functional equation, math.stackexchange.com/questions/143449/…. – Tom Copeland May 10 '12 at 23:03

I'm not sure of the history of the gamma factor, though I would suggest that no one "tried it out", but rather it simply arose in trying to prove of the functional equation. Riemann was the first to prove the functional equation, and his proof essentially follows that in Harald Hanche-Olsen's answer, which makes my explanation plausible. Alternatively, the functional equation of the zeta function comes out of the functional equation of a theta series, and the Mellin transform of a theta series gives rise to a Gamma function. This latter explanation arises more naturally for modular forms: the L-function of a modular form is also completed by a gamma factor to obtain a functional equation; in this case, the completed L-function is simply the Mellin transform of the modular form itself.

Furthermore, as Leonid Positselski answers, it is indeed true that Tate's thesis provides a uniform way of obtain the gamma factors at infinity in the same manner as one obtains the local L-factors at finite places.

More generally, there is a recipe given an arbitrary motive for the expected gamma factors that should give a functional equation for the motivic L-functions. These are due to Deligne and Serre (I believe) and are determined by the Hodge structure of the motive (see Deligne's corvallis article "Valeurs de fonctions L..."). This shows that there's a uniform way of obtaining the gamma factors as one varies the L-function one is studying, an orthogonal question to the one Leonid Positselski answered.

-

Gamma function arises when we consecutively differentiate an Appell sequence. An example of Appell polynomials are Bernoulli polynomials. When we differentiate it, the factors combine with themselves:

$$B_n'(x)=nB_{n-1}(x)$$

$$B_n''(x)=n(n-1)B_{n-2}(x)$$

$$B_n'''(x)=n(n-1)(n-2)B_{n-3}(x)$$

They are just another name for Hurwitz Zeta function:

$$B_n(x) = -n \zeta(1-n,x)$$

Thus, for $f(s,q)=\zeta(s,-q)$

$$\frac\partial{\partial q}f(s,q)= s f(s+1,q)$$

$$\frac{\partial^2}{\partial q^2}f(s,q)= s(s+1) f(s+2,q)$$

$$\frac{\partial^3}{\partial q^3}f(s,q)= s(s+1)(s+2) f(s+3,q)$$

Since Reihmann zeta is Hurwitz zeta evaluated at $q=1$, the expression you give is apparently consecutive derivative of Hurwitz Zeta, with factor $\pi^{-s}$ appearing if we normalize Hurwitz Zeta by stretching it horizontally by factor of pi.

Consecutive derivatives of Hurwitz Zeta in turn are nothing more than just polygamma function.

For instance, here is the function $-1/x$:

If we add infinitely many similar functions with a shift of pi/2 each in both directions, we get $\tan x$. But if we do the same only in one direction, we get "incomplete tangent":

The yellow one is $\operatorname{pg}(x)=\frac 1\pi \psi (\frac x\pi)$, the blue one is $\operatorname{cpg}(x)=-\frac 1\pi \psi (1-\frac x\pi)$. They obey $\operatorname{cpg}(x)+\operatorname{pg}(x)=-\cot(x)$.

Now if we differentiate cpg(x) we get:

$$(\operatorname{cpg}(x))^{(s-1)}=\pi^{-s}\Gamma(s)\zeta(s,1-\frac x\pi)$$

Compare it with yours formula:

$$\xi(2s) = \pi^{-s}\Gamma\left(s\right)\zeta(2s)$$

-