Dan Romik on the Riemann zeta function

This post, about the Riemann zeta function, which is among the most important and mysterious mathematical objects was kindly written by Dan Romik. It is related to his paper Orthogonal polynomial expansions for the Riemann xi function,  that we mentioned in this post.

Dan Romik on the Riemann zeta function

Recently when I was thinking about the Riemann zeta function, I had the double thrill of discovering some new results about it, and then later finding out that my new ideas were closely related to some very classical ideas due to two icons of twentieth-century mathematics, George Pólya and Pál Turán. When you are trying to stand on the shoulders of giants, it’s nice to see other giants right there beside you trying to do the same!

It all goes back to one of the most famous problems in mathematics, the Riemann Hypothesis (RH). Both Pólya and Turán were rather enamored with this problem and published about it extensively; Pólya was said to have been preoccupied with the problem to the very end of his life.(1) And they both recognized that an important first step in trying to prove something about the zeros of the zeta function is having a good representation
for the Riemann zeta function. After all, there are many different formulas that can be used to define or compute the zeta function. If you don’t choose the right one, you probably won’t get very far with your analysis.

Pólya in one of his famous attacks on the problem considered the representation of the zeta function (or more precisely of the Riemann xi function, which is a symmetrized and better-behaved version of the zeta function; see below) as a Fourier transform—a standard representation due (essentially) to Riemann. I’ll have more to say about that later.

Turán also looked at the Riemann xi function, and instead of working with one of the standard “named” representations such as the Fourier transform or Taylor series, looked around a bit more intentionally for a representation of the function that seemed particularly suited to answering the specific question of whether the zeros all lie on a line. In a 1950 address to the Hungarian Academy of Sciences, he put forward his ideas about what he thought was the correct representation to look at: the infinite series expansion of the xi function in the Hermite polynomials. About eighty years after Turán’s discovery, my own investigations led me to discover [5] that the Hermite polynomials are not the only polynomials in which it’s interesting to expand the Riemann xi function. It turns out that there are at least two other families of polynomials for which the respective expansions are no less (and, in some ways, more) well-behaved. My motto for these polynomial families, which are known to experts in special functions but have until now been somewhat esoteric (though I hope that is about to change), is that they are “the coolest polynomials that you never heard of.”

Let’s look at some of the technical details so that I can explain why these new expansions are interesting, and how they relate to Turán’s work and ultimately back to Pólya’s ideas and one of the particular threads that grew out of them. First, define the Riemann xi function as

\displaystyle \xi(s) = \frac12 s(s-1) \pi^{-s/2} \Gamma\left(\frac{s}{2}\right) \zeta(s) \qquad (s\in\mathbb{C}),

where {\Gamma(z)} is the Euler gamma function and {\zeta(s)} is the Riemann zeta function. It’s also common to denote
\displaystyle \Xi(t) = \xi\left(\frac12+it\right) \qquad (t\in\mathbb{C}).

This is Riemann’s “capital xi” function, which is still usually referred to as Riemann’s xi function. (This seems reasonable: the two functions are the same up to a trivial linear change of variables.) The main point of these definitions is that {\Xi(t)} is an entire function of the complex variable {t}, and that RH can now be reformulated as the statement that the zeros of {\Xi(t)} all lie on the real line. Moreover, the famous functional equation satisfied by the Riemann zeta function maps to the statement that {\Xi(t)} is an even function.
Now consider the following four ways of representing the xi function:

\displaystyle \Xi(t) = \sum_{n=0}^\infty (-1)^n a_{2n} t^{2n},~~~~~(1)

\displaystyle \Xi(t) = \sum_{n=0}^\infty (-1)^n b_{2n} H_{2n}(t),~~~~~(2)

\displaystyle\Xi(t) = \sum_{n=0}^\infty (-1)^n c_{2n} f_{2n}\left(\frac{t}{2}\right),~~~~~(3)

\displaystyle \Xi(t) = \sum_{n=0}^\infty (-1)^n d_{2n} g_{2n}\left(\frac{t}{2}\right).~~~~~(4)

Here, the first representation (1) is simply the Taylor expansion of {\Xi(t)}, which contains only even terms since {\Xi(t)} is an even function. The numbers {a_{2n}} are (up to the {(-1)^n} sign factor) the Taylor coefficients. Some attempts have been made to understand them, and one interesting and fairly trivial observation (again going back to facts already known to Riemann) is that they are all positive. Some additional and less trivial things can be said—see for example Section 6.1 of my paper [5], and the recent paper by Griffin, Ono, Rolen and Zagier [2]. But at the end of the day, no one has yet succeeded in using the Taylor expansion to prove anything new about the location of the zeros.

The second representation (2) is the infinite series expansion of {\Xi(t)} in the classical sequence of Hermite polynomials, defined by the well-known formula

\displaystyle H_n(t) = (-1)^n e^{t^2} \frac{d^n}{dt^n} \left( e^{-t^2} \right).

This is the representation whose use was advocated by Turán. His reasoning was that expanding a function of a complex variable (for example, in the simplest case, a polynomial) in monomials {t^n} doesn’t provide useful information to easily decide if the function has only real zeros, because the monomials have, roughly speaking, radial symmetry: their level curves are concentric circles. The Hermite polynomials on the other hand, at least heuristically, have level curves that are closer to being straight lines parallel to the real axis, Turán argued; thus, they are more suited to the geometry of the problem we are trying to solve.

Turán’s case for supporting the Hermite polynomials as the right basis to use is quite detailed—you can read about it in his papers [6,7,8] (and no, he was not able to actually prove anything about the zeros of {\Xi(t)}; this is a common theme in most of the attacks on RH to date…). I’ll simply mention that again one interesting and fairly easy observation is that the coefficients {b_{2n}} in the expansion (2)—adjusted through the introduction of the sign factor {(-1)^n}—end up being positive numbers. Their asymptotic behavior can also be analyzed: I prove a result about this in my paper (though it’s not particularly pretty).

Now comes the part that to me seems the most exciting, involving the expansions (3) and (4). These are the expansions in the more exotic families of polynomials

\displaystyle f_n(x)=(-i)^n \sum_{k=0}^n 2^k\binom{n+\frac12}{n-k}\binom{-\frac34+ix}{k},

\displaystyle g_n(x)= (-i)^n \sum_{k=0}^n \frac{(n+k+1)!}{(n-k)!(3/2)_k^2} \binom{-\frac34+ix}{k}

(where {(3/2)_n} is a Pochhammer symbol), mildly rescaled by replacing {x}  with {t/2}. In the terminology of the theory of orthogonal polynomials, the family {f_n(x)} is a special case of a two-parameter family {P_n^{(\lambda)}(x;\phi)} known as the Meixner-Pollaczek polynomials, with the parameters taking the particular values {\phi=\frac{\pi}{2}, \lambda=\frac{3}{4}}. Similarly, the family {g_n(x)} is a special case of the four-parameter family {p_n(x;a,b,c,d)} known as the continuous Hahn polynomials, with the parameters taking the particular values {a=b=c=d=\frac{3}{4}}. Their main characterizing property is that they are orthogonal sequences of polynomials for two specific weight functions on {\mathbb{R}}: the {f_n(x)} are orthogonal with respect to the weight function {w_1(x)=\left|\Gamma\left(\frac34+ix\right)\right|^2}, and the {g_n(x)} are orthogonal with respect to {w_2(x)=\left|\Gamma\left(\frac34+ix\right)\right|^4}. Again, fairly esoteric. But interesting!

There are several things that make the expansions (3)–(4) well-behaved. First, the coefficients {c_{2n}}, {d_{2n}} are again positive. This actually seems pretty relevant for questions like RH: for example, if we consider “toy” versions of (1)–(3) in which the coefficient sequences {a_n}, {b_n} and {c_n} are replaced by the sequence {\alpha^n} for fixed {0<\alpha<1}, all three expansions sum up to rescaled cosines, which are entire functions that of course have only real zeros. (Without the {(-1)^n} factor, we would get a hyperbolic cosine, which has imaginary zeros.)

Second, one can derive asymptotics for {c_{2n}} and {d_{2n}}, and they are quite a bit nicer than the asymptotic formulas for the Taylor and Hermite expansion coefficients. In my paper, I proved that \displaystyle c_{2n} \sim A \sqrt{n} e^{-B \sqrt{n}}, \qquad d_{2n} \sim C n^{4/3} e^{-D n^{2/3}} \qquad \textrm{as }n\rightarrow\infty, where {A,B,C,D} are the constants \displaystyle A = 16\sqrt{2}\pi^{3/2}, \qquad B = 4\sqrt{\pi}, \qquad C = \frac{128 \times 2^{1/3} \pi^{2/3} e^{-2\pi /3}}{\sqrt{3}}, \qquad D = 3 (4\pi)^{1/3}.

Third, the expansions have some conceptual meaning: (3) turns out to be equivalent to the expansion of the elementary function {\frac{d^2}{du^2} (u \coth(\pi u))}, {u>0}, in an orthogonal basis of functions related to the Laguerre polynomials {L_n^{1/2}(x)}. And analogously, (4) arises out of the expansion of a certain auxiliary function {\tilde{\nu}(u)} (I won’t define it here) in yet another classical family of orthogonal polynomials, the Chebyshev polynomials of the second kind.

Fourth (and fifth, sixth, …): the expansions are just… nice, in the sense that they arise in a way that seems natural when one asks certain questions, that they have excellent convergence properties, and that the coefficients {c_{2n}} and {d_{2n}} have several elegant formulas, each revealing something interesting about them. Read the paper to understand more.

I said I will get back to Pólya’s work on RH. This post is already quite long so I will say only a little bit about this. One of Pólya’s major discoveries was that there are operations on entire functions that (under certain mild assumptions) preserve the property of the function having only real zeros. Specifically this is the case for the operation of multiplying the Fourier transform of the function by the factor {e^{\lambda u^2}} for {\lambda>0}  (where {u} is the frequency variable). This opens the way to defining a family of deformations {\Xi_\lambda(t)} of the Riemann xi function arising out of this operation, and trying to generalize RH by asking for which values of {\lambda} it is the case that {\Xi_\lambda(t)} has only real zeros. Since Pólya’s work, and important later extensions of it by De Bruijn and Newman, this has become a very active topic of research, nowadays referred to under the name of the De Bruijn-Newman constant.
See the recent survey of Newman and Wu [3], a 2018 paper by Rodgers and Tao [4] proving a major conjecture of Newman, and the recent paper [9] by the Polymath15 project (mentioned by Gil in his earlier post), for the latest on this subject.

The connection I found between this topic and the idea of expanding the Riemann xi function in families of orthogonal polynomials is the following: expansions such as (2)–(4) suggest yet another natural way of “deforming” the Riemann xi function by adding a parameter {\alpha}: simply multiply the {n}th term in the expansion by {\alpha^{2n}} (the linear operator that does this is called the Poisson kernel, and generalizes the standard Poisson kernel from complex analysis and the theory of harmonic functions). It turns out—and is actually easy to prove, and really isn’t terribly surprising in the grand scheme of things—that in the case of the Hermite expansion (2), this family of deformations is the same, up to some trivial reparametrization, as the family of deformations {\Xi_\lambda(t)} that was studied in connection with the work of Pólya, De Bruijn, Newman and their successors. A nice connection between two threads of research that were not previously recognized as being related to each other, I think. Furthermore, this suggests that the Poisson kernel and associated deformations may yet have an important role to play in the context of the new expansions in the orthogonal polynomial families {f_n} and {g_n}, where we get genuinely new families of deformations of the Riemann xi function. I explore this idea in my paper and it leads to some interesting things.

So let’s summarize. The key questions you are no doubt wondering about are: where does any of this lead? And do these new ideas say anything really useful or especially relevant for the Riemann hypothesis? The answer is that I don’t know (and I’m wondering about the same things). That being said, these orthogonal polynomial expansions seem quite interesting in their own right. The Riemann zeta function is a mysterious object, and there are other things we wish to understand about it beside where its zeros are, so it’s always good to have additional points of view from which to approach it. Moreover, even on the question of the zeros there are reasons to be cautiously optimistic that this approach may have something useful to offer; see Chapter 7 of my paper for a brief discussion of why that is the case.


[1] D. Albers and G. L. Alexanderson, editors. Mathematical People: Profiles and Interviews. A K Peters, 2008.

[2] M. Griffin, K. Ono, L. Rolen and D. Zagier. Jensen polynomials for the Riemann zeta function and other sequences. Preprint (2019), arXiv:1902.07321.

[3] C. M. Newman. Constants of de Bruijn-Newman type in analytic number theory and statistical physics. To appear in Bull. Amer. Math. Soc.

[4] B. Rodgers and T. Tao. The De Bruijn-Newman constant is nonnegative. Preprint (2018), arXiv:1801.05914.

[5] D. Romik. Orthogonal polynomial expansions for the Riemann xi function. Preprint (2019), arXiv:1902.06330.

[6] P. Turán. Sur l’algèbre fonctionelle. Pages 279–290 in: Comptes Rendus du Premier Congrès des Mathématiciens Hongrois, 27 Août–2 Septembre 1950. Akadémiai Kiadó, 1952. Reprinted in Collected Papers of Paul Turán, Ed. P. Erdos, Vol. 1, pp. 677–688. Akadémiai Kiadó, 1990. An English translation of the paper by Dan Romik On functional algebra.

[7] P. Turán. Hermite-expansion and strips for zeros of polynomials. Arch. Math. 5 (1954), 148–152. Reprinted in Collected Papers of Paul Turán, Ed. P. Erdos, Vol. 1, pp. 738–742. Akadémiai Kiadó, 1990.

[8]  P. Turán. To the analytical theory of algebraic equations. Bulgar. Akad. Nauk. Otd. Mat. Fiz. Nauk. Izv. Mat. Inst. 3 (1959), 123–137. Reprinted in Collected Papers of Paul Turán, Ed. P. Erdos, Vol. 2, pp. 1080–1090. Akadémiai Kiadó, 1990.

[9] D.H.J. Polymath. Effective approximation of heat flow evolution of the Riemann ξ function, and a new upper bound for the de Bruijn-Newman constant. Preprint (2019), arXiv:1904.12438.


(1) Alexanderson writes in [1, p. 259]: “A week or so before he died, Pólya asked me to look on his desk at home for some papers on which he said he had written down some interesting ideas he had for proving RH. Of course I could find no such notes, but until the day he died he was thinking about that famous problem.”


(2) Turán’s Hungarian Academy of Sciences talk was published in a rather obscure French-language paper [6] that seems to have been largely forgotten. It’s an interesting read nonetheless, and to make it more accessible to anyone who may be interested, I recently translated it to English.


(3) Turán mentions in [8] that he discovered the results on the Hermite expansion in 1938–39, but they were not published until much later. Clearly this was not a convenient time in history for publishing such discoveries; Turán, a Hungarian Jew, spent much of World War II interned in labor camps in Hungary.

This entry was posted in Combinatorics, Guest blogger, Number theory and tagged , , , , . Bookmark the permalink.

2 Responses to Dan Romik on the Riemann zeta function

  1. vznvzn says:

    nice, thx for sharing! so which is harder? P vs NP or Riemann?
    suspect the concept of fractals can be helpful in analyzing both problems but very few seem to be pursuing this angle. this interesting 2010 paper by May has 26 refs to fractals in the contents but not in the title or abstract. whats the deal?

    Click to access king2010.pdf

  2. Monti E. ARRIS says:

    Fundamental to the problem are two issues. The first is that the duality already seen in the
    relationship between zeta and the prime products is already the duality transform one is seeking.
    That is the system that 5th decodes the zeta zeros is the distribution of numerical primes itself, so
    seeking an analogue from other mathematical areas cannot necessarily simplify the problem.
    Secondly these GUE systems may show similarities in their statistics to zeta’s zeros, because they
    share overall features combining structured constraints and pseudo-randomness in common with
    the primes and zeta zeros, without necessarily being isomorphic to them. In a sense there is a
    regress occurring, in which attempting to model GUE systems to zeta’s zeros results in more
    elaborate mathematical constructions which share zeta’s characteristics but neither provide a
    breakthrough in proving RH, nor result in a real valued quantum operator.
    The quantum is direct

    y = 0.57721 56649 01532 86060 65120 90082 40243

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out /  Change )

Google photo

You are commenting using your Google account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

Connecting to %s