The Prime Number Theorem and Dirichlet's Theorem

  • Prime Number Theorem: π(x) ~ x/log x asymptotically; probability a random integer near x is prime ≈ 1/log x
  • Proved 1896 by Hadamard and de la Vallée Poussin using complex analysis; conjectured by Gauss and Legendre ~1800
  • Dirichlet's Theorem (1837): primes equidistribute among invertible residue classes modulo any integer in arithmetic progressions
  • Launched analytic number theory by resolving discrete arithmetic problems via continuous analysis and zeta functions
  • Riemann zeta function and Dirichlet L-functions connect prime distribution to complex analysis; Riemann Hypothesis remains open

The Prime Number Theorem and Dirichlet's Theorem

Overview

The distribution of prime numbers among the integers represents one of the deepest mysteries in mathematics, transforming from empirical observation to rigorous theory through the nineteenth century. The Prime Number Theorem (PNT) and Dirichlet's Theorem on Primes in Arithmetic Progressions together establish the fundamental asymptotic laws governing how primes are distributed, both globally and within arithmetic structures. While individual primes appear chaotic, these theorems reveal profound regularity in their collective behavior. The Prime Number Theorem asserts that the number of primes up to xx, denoted π(x)\pi(x), is asymptotic to x/log⁡xx/\log x, meaning the "probability" that a random integer near xx is prime is approximately 1/log⁡x1/\log x. This result, conjectured by Gauss and Legendre around 1800 based on numerical evidence, remained unproved until 1896 when Hadamard and de la Vallée Poussin independently established it using complex analysis. Dirichlet's earlier theorem (1837) showed that primes are equidistributed among the invertible residue classes modulo any integer, establishing arithmetic progressions as the natural setting for understanding prime distribution. These theorems inaugurated analytic number theory, demonstrating that questions about discrete arithmetic objects could be resolved through continuous analysis. The interplay between the Riemann zeta function, Dirichlet L-functions, and prime distribution remains one of the most fertile areas in mathematics, with the Riemann Hypothesis representing the most famous open problem concerning the finer structure of prime distribution.

💡 The Core Insight

Primes, though individually unpredictable, obey precise statistical laws: globally, their density near xx is 1/log⁡x1/\log x, and locally, they distribute equally among all allowable arithmetic progressions. The bridge between primes and analysis is built through generating functions that encode arithmetic data in their analytic properties.

Definition

The Prime Number Theorem and Dirichlet's Theorem: 1. (Prime Number Theorem) Let π(x)=#{p≤x:p prime}\pi(x) = \#\{p \leq x : p \text{ prime}\} denote the prime counting function. Then

π(x)∼xlog⁡xas x→∞,\pi(x) \sim \frac{x}{\log x} \quad \text{as } x \to \infty,

meaning lim⁡x→∞π(x)x/log⁡x=1\lim_{x \to \infty} \frac{\pi(x)}{x/\log x} = 1. 2. (Dirichlet's Theorem) Let a,q∈Za, q \in \mathbb{Z} with gcd⁡(a,q)=1\gcd(a,q) = 1. Then the arithmetic progression {a+nq:n∈N}\{a + nq : n \in \mathbb{N}\} contains infinitely many primes. More precisely, let π(x;q,a)=#{p≤x:p prime,p≡a(modq)}\pi(x; q, a) = \#\{p \leq x : p \text{ prime}, p \equiv a \pmod{q}\}. Then

π(x;q,a)∼1ϕ(q)⋅xlog⁡xas x→∞,\pi(x; q, a) \sim \frac{1}{\phi(q)} \cdot \frac{x}{\log x} \quad \text{as } x \to \infty,

where ϕ(q)\phi(q) is the <span class="wiki-link-missing">Euler totient function</span>.

Core Theory

TheoremMain Result

Theorem [1] (Prime Number Theorem — Analytic Form): The Riemann zeta function ζ(s)=∑n=1∞n−s\zeta(s) = \sum_{n=1}^{\infty} n^{-s} for ℜ(s)>1\Re(s) > 1 extends to a meromorphic function on C\mathbb{C} with a simple pole at s=1s=1, and satisfies ζ(s)≠0\zeta(s) \neq 0 for ℜ(s)=1\Re(s) = 1. This implies π(x)∼x/log⁡x\pi(x) \sim x/\log x.

Theorem [2] (Dirichlet's Theorem): For gcd⁡(a,q)=1\gcd(a,q) = 1, let χ\chi be a Dirichlet character modulo qq, and define the Dirichlet L-function

L(s,χ)=∑n=1∞χ(n)ns,ℜ(s)>1.L(s, \chi) = \sum_{n=1}^{\infty} \frac{\chi(n)}{n^s}, \quad \Re(s) > 1.

Then L(1,χ)≠0L(1, \chi) \neq 0 for all non-principal characters χ\chi, which implies infinitely many primes in each residue class coprime to qq.

Proof

Proof of Theorem [1] (Sketch of classical approach): Step 1: Euler product and logarithmic connection. For ℜ(s)>1\Re(s) > 1, the Euler product formula gives

ζ(s)=∏p prime11−p−s.\zeta(s) = \prod_{p \text{ prime}} \frac{1}{1 - p^{-s}}.

Taking logarithms,

log⁡ζ(s)=−∑plog⁡(1−p−s)=∑p(p−s+O(p−2s))=∑pp−s+O(1).\log \zeta(s) = -\sum_p \log(1 - p^{-s}) = \sum_p \left( p^{-s} + O(p^{-2s}) \right) = \sum_p p^{-s} + O(1).

Step 2: Connection to prime counting. Define the Chebyshev function ψ(x)=∑pk≤xlog⁡p\psi(x) = \sum_{p^k \leq x} \log p. By Perron's formula and Contour Integration,

ψ(x)=12πi∫c−i∞c+i∞(−ζ′(s)ζ(s))xssds\psi(x) = \frac{1}{2\pi i} \int_{c - i\infty}^{c + i\infty} \left(-\frac{\zeta'(s)}{\zeta(s)}\right) \frac{x^s}{s} ds

for suitable c>1c > 1. The PNT is equivalent to ψ(x)∼x\psi(x) \sim x. Step 3: Analytic continuation. The functional equation

π−s/2Γ(s/2)ζ(s)=π−(1−s)/2Γ((1−s)/2)ζ(1−s)\pi^{-s/2} \Gamma(s/2) \zeta(s) = \pi^{-(1-s)/2} \Gamma((1-s)/2) \zeta(1-s)

shows ζ(s)\zeta(s) extends meromorphically to C\mathbb{C} with only a simple pole at s=1s=1 (residue 1). Step 4: Zero-free region. The critical step: showing ζ(s)≠0\zeta(s) \neq 0 on ℜ(s)=1\Re(s) = 1. For σ>1\sigma > 1, the inequality

ζ(σ)3∣ζ(σ+it)∣4∣ζ(σ+2it)∣≥1\zeta(\sigma)^3 |\zeta(\sigma + it)|^4 |\zeta(\sigma + 2it)| \geq 1

follows from 3+4cos⁡θ+cos⁡2θ≥03 + 4\cos\theta + \cos 2\theta \geq 0. As σ→1+\sigma \to 1^+, if ζ(1+it)=0\zeta(1+it) = 0, the left side would tend to 0, contradiction. Step 5: Shifting the contour. The zero-free region on ℜ(s)=1\Re(s) = 1 allows shifting the integration contour in Step 2 to ℜ(s)=1\Re(s) = 1, picking up the residue at s=1s=1 (which is xx) and showing error terms are o(x)o(x). Thus ψ(x)∼x\psi(x) \sim x, which implies π(x)∼x/log⁡x\pi(x) \sim x/\log x by partial summation.

Proof of Theorem [2] (Sketch): Step 1: Orthogonality of characters. The Dirichlet characters modulo qq form an orthogonal basis for class functions on (Z/qZ)∗(\mathbb{Z}/q\mathbb{Z})^*. By orthogonality,

∑χ  qχ(a)‾L(s,χ)=∑n≡a(modq)1ns.\sum_{\chi \bmod q} \overline{\chi(a)} L(s, \chi) = \sum_{n \equiv a \pmod{q}} \frac{1}{n^s}.

Step 2: Showing L(1,χ)≠0L(1, \chi) \neq 0. For the principal character χ0\chi_0, we have L(s,χ0)=ζ(s)∏p∣q(1−p−s)L(s, \chi_0) = \zeta(s) \prod_{p \mid q}(1 - p^{-s}), which has a pole at s=1s=1. For non-principal χ\chi:

  • If χ2=χ0\chi^2 = \chi_0 (quadratic), use class number formula connections
  • In general, if L(1,χ)=0L(1, \chi) = 0, analyzing ∏χL(s,χ)\prod_{\chi} L(s, \chi) near s=1s=1 via the product formula leads to contradictions with growth estimates

Step 3: Conclusion. Since L(1,χ0)L(1, \chi_0) has a pole and L(1,χ)≠0L(1, \chi) \neq 0 for non-principal χ\chi, the orthogonality sum shows infinitely many primes p≡a(modq)p \equiv a \pmod{q}, with the asymptotic following from refined analysis. □

∎
⚠️ Common mistake

Common error: The PNT does NOT say π(x)=x/log⁡x\pi(x) = x/\log x, only that the ratio tends to 1. The better approximation is the <span class="wiki-link-missing">logarithmic integral</span> li(x)=∫2xdtlog⁡t\text{li}(x) = \int_2^x \frac{dt}{\log t}, which satisfies li(x)∼x/log⁡x\text{li}(x) \sim x/\log x but li(x)−π(x)=O(xe−clog⁡x)\text{li}(x) - \pi(x) = O(x e^{-c\sqrt{\log x}}). Also, Dirichlet's theorem requires gcd⁡(a,q)=1\gcd(a,q) = 1; there are obviously no primes p≡0(modq)p \equiv 0 \pmod{q} for q>1q > 1 except possibly p=qp = q itself.

mathematician

Peter Gustav Lejeune Dirichlet (1837) — Proved the first deep result on primes in arithmetic progressions using the analytic properties of L-functions he introduced, creating the field of analytic number theory.

visualization

x π(x) 0 50 100 150 π(x) x/log x li(x)

Examples and Applications

Example1

Computing π(100) and comparing with PNT approximation The primes up to 100 are: 2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97. Thus π(100)=25\pi(100) = 25. The PNT approximation gives:

100log⁡100=100ln⁡100=1004.605...≈21.71\frac{100}{\log 100} = \frac{100}{\ln 100} = \frac{100}{4.605...} \approx 21.71

The relative error is ∣25−21.71∣25≈13.2%\frac{|25 - 21.71|}{25} \approx 13.2\%. The better approximation via the logarithmic integral:

li(100)=∫2100dtlog⁡t≈30.12\text{li}(100) = \int_2^{100} \frac{dt}{\log t} \approx 30.12

This overshoots but is known to be systematically higher than π(x)\pi(x) for small xx (though this pattern eventually reverses, conjecturally around x≈10316x \approx 10^{316} by Skewes' number). For larger values: π(106)=78498\pi(10^6) = 78498 while 106/log⁡(106)≈7238210^6/\log(10^6) \approx 72382 (error ~8%), showing convergence of the asymptotic.

Example2

Primes in arithmetic progressions modulo 4 Consider q=4q = 4. The invertible residue classes are a∈{1,3}a \in \{1, 3\} with ϕ(4)=2\phi(4) = 2. Dirichlet's theorem predicts equal distribution:

π(x;4,1)∼π(x;4,3)∼12⋅xlog⁡x\pi(x; 4, 1) \sim \pi(x; 4, 3) \sim \frac{1}{2} \cdot \frac{x}{\log x}

For x=100x = 100:

  • Primes ≡1(mod4)\equiv 1 \pmod{4}: 5, 13, 17, 29, 37, 41, 53, 61, 73, 89, 97 — total: 11
  • Primes ≡3(mod4)\equiv 3 \pmod{4}: 3, 7, 11, 19, 23, 31, 43, 47, 59, 67, 71, 79, 83 — total: 13

Prediction: 12⋅21.71≈10.85\frac{1}{2} \cdot 21.71 \approx 10.85 for each class. The counts (11 and 13) are close and roughly equal, as predicted. The equidistribution becomes more pronounced as xx increases. This connects to the fact that the non-principal character modulo 4 is χ(n)={1n≡1(mod4)−1n≡3(mod4)\chi(n) = \begin{cases} 1 & n \equiv 1 \pmod{4} \\ -1 & n \equiv 3 \pmod{4} \end{cases}, and L(1,χ)=1−1/3+1/5−1/7+...=π/4≠0L(1, \chi) = 1 - 1/3 + 1/5 - 1/7 + ... = \pi/4 \neq 0.

Counterexample

Why gcd condition is necessary in Dirichlet's theorem Consider the arithmetic progression 2,6,10,14,18,22,...2, 6, 10, 14, 18, 22, ... (i.e., a=2,q=4a = 2, q = 4). Here gcd⁡(2,4)=2≠1\gcd(2, 4) = 2 \neq 1. This progression contains only one prime: p=2p = 2 itself. All other terms are composite (even numbers >2> 2). This shows the coprimality condition gcd⁡(a,q)=1\gcd(a,q) = 1 cannot be dropped. The underlying reason is group-theoretic: characters are defined on the unit group (Z/qZ)∗(\mathbb{Z}/q\mathbb{Z})^*, and orthogonality relations only apply to invertible classes. Elements sharing a common factor with qq don't participate in the multiplicative structure that L-functions detect. More generally, if d=gcd⁡(a,q)>1d = \gcd(a,q) > 1, then every element of a+qZa + q\mathbb{Z} is divisible by dd, so at most one prime (namely p=dp = d if it divides aa) can occur.

🧠 Intuition

The sieve perspective: Imagine "sieving" the integers, removing multiples of each prime in turn (Sieve of Eratosthenes). After removing multiples of all primes up to x\sqrt{x}, approximately x∏p≤x(1−1/p)x \prod_{p \leq \sqrt{x}}(1 - 1/p) numbers remain. By Mertens' theorems, this product behaves like C/log⁡xC/\log x for some constant CC, explaining the x/log⁡xx/\log x density. For arithmetic progressions, imagine the integers arranged in columns by residue class modulo qq. Sieving acts independently on each column (except for primes dividing qq), so primes distribute equally among the ϕ(q)\phi(q) allowable columns. The L-function machinery makes this heuristic precise through harmonic analysis on the finite abelian group (Z/qZ)∗(\mathbb{Z}/q\mathbb{Z})^*.

Probabilistic model: If we randomly assign each integer nn to be "prime" with probability 1/log⁡n1/\log n, independently, we'd expect about ∫2xdtlog⁡t∼x/log⁡x\int_2^x \frac{dt}{\log t} \sim x/\log x primes up to xx. The PNT says the actual deterministic distribution matches this random model asymptotically—a deep instance of pseudorandomness in number theory.

Historical Development

The empirical investigation of prime distribution began with Gauss and Legendre around 1792-1798. Gauss, through extensive hand computation, conjectured that π(x)∼li(x)\pi(x) \sim \text{li}(x), recorded in his unpublished notes. Legendre proposed the form x/(Alog⁡x+B)x/(A\log x + B) based on his tables.

Dirichlet proved his theorem on primes in arithmetic progressions in 1837, published in "Beweis des Satzes, dass jede unbegrenzte arithmetische Progression...". This work introduced Dirichlet characters, Dirichlet L-functions, and established the paradigm of using analytic properties (specifically non-vanishing at s=1s=1) to deduce arithmetic results. His techniques relied heavily on Fourier analysis on finite abelian groups.

Riemann's revolutionary 1859 paper "Über die Anzahl der Primzahlen unter einer gegebenen Größe" connected π(x)\pi(x) explicitly to zeros of ζ(s)\zeta(s) via the explicit formula

π(x)=li(x)−∑ρli(xρ)+lower order terms\pi(x) = \text{li}(x) - \sum_{\rho} \text{li}(x^\rho) + \text{lower order terms}

where ρ\rho ranges over non-trivial zeros of ζ(s)\zeta(s). He conjectured these zeros all have real part 1/21/2 (Riemann Hypothesis). The first complete proof of PNT came in 1896 independently from Jacques Hadamard ("Sur la distribution des zéros de la fonction ζ(s)\zeta(s)...") and Charles-Jean de la Vallée Poussin ("Recherches analytiques sur la théorie des nombres premiers"), both using complex analysis to establish the zero-free region ℜ(s)=1\Re(s) = 1. Their methods were refined by Landau in his systematic treatise.

Erdős and Selberg shocked the mathematical world in 1949 by discovering elementary proofs (avoiding complex analysis), though these are technically intricate. The analytic approach remains standard for its power and generalizations. Modern developments include effective versions (explicit error bounds), the explicit formula, connections to random matrix theory (zeros behaving like eigenvalues of random Hermitian matrices), and vast generalizations to automorphic L-functions and Langlands program.

Connections to Other Areas

connection

Random Matrix Theory and the Montgomery-Odlyzko Law: In the 1970s, Hugh Montgomery discovered that the pair correlation of zeros of ζ(s)\zeta(s) (when rescaled appropriately on the critical line) matches the pair correlation of eigenvalues of random matrices from the Gaussian Unitary Ensemble (GUE). Andrew Odlyzko's extensive numerical computations in the 1980s confirmed this with stunning precision. This connection suggests deep relationships between prime number distribution and Quantum Chaos, since GUE eigenvalues model energy levels of quantum systems with chaotic classical limits. The correspondence remains mysterious but has led to using random matrix theory to predict statistics of zeros and primes, with applications to cryptography and models in quantum physics.

What You Need First

What This Opens Up

See Also

Riemann Zeta Function · Dirichlet L-functions · Prime Number Distribution · Analytic Number Theory · Riemann Hypothesis · Sieve Theory · Explicit Formulas · Class Number Formula · Mertens Theorems · Chebyshev Functions