The Connection Between the Riemann Zeta Function and Prime Distribution

---

The Eternal Question of Primes

One of the most foundational pillars of mathematics is the realization that prime numbers never end. This was established by the geometer Euclid around 300 BCE. For many centuries, this discovery was seen as the final word on the matter.

"If there are infinitely many primes, is there anything left to discover about them?"

The answer is a resounding yes. To understand why, we must look at the "density" of different sets of numbers.

Quantifying "Commonality"

Consider the difference between even numbers and powers of 2. While both are infinite sets, they occupy the number line very differently.

We can mathematically express this difference by summing their reciprocals:

• Even Numbers: $\frac{1}{2} + \frac{1}{4} + \frac{1}{6} + \frac{1}{8} + \cdots = \infty$ (Divergent)
• Powers of 2: $\frac{1}{2} + \frac{1}{2^2} + \frac{1}{2^3} + \frac{1}{2^4} + \cdots = 1$ (Convergent)

In essence, powers of 2 grow so rapidly that their reciprocals shrink fast enough to hit a finite ceiling. Even numbers, however, stay "close" enough together that their sum climbs forever. This led Leonhard Euler to wonder about the primes:

$\text{Is } \sum_{p\text{ prime}} \frac{1}{p} = \frac{1}{2} + \frac{1}{3} + \frac{1}{5} + \frac{1}{7} + \frac{1}{11} + \cdots \text{ finite or infinite?}$

---

Introducing the Riemann Zeta Function

To tackle this, Euler utilized a specific mathematical object. Although it is now called the Riemann zeta function (named after Bernhard Riemann, who expanded upon it much later), Euler was the one who first uncovered its secrets.

The function is defined as the infinite sum: $\zeta(s) = \sum_{n=1}^{\infty} \frac{1}{n^s} = 1 + \frac{1}{2^s} + \frac{1}{3^s} + \frac{1}{4^s} + \cdots$

The Euler Product Formula

Euler's breakthrough was realizing that $\zeta(s)$ could be rewritten as a product involving only prime numbers. Imagine multiplying these infinite series together:

$\left(1+\frac{1}{2^s}+\frac{1}{4^s}+\cdots\right)\left(1+\frac{1}{3^s}+\frac{1}{9^s}+\cdots\right)\left(1+\frac{1}{5^s}+\frac{1}{25^s}+\cdots\right) \cdots$

When you distribute these terms, you generate every possible combination of prime powers. For example:

• $1 \cdot 1 \cdot 1 = 1$
• $\frac{1}{2^s} \cdot 1 \cdot 1 = \frac{1}{2^s}$
• $\frac{1}{2^s} \cdot \frac{1}{3^s} \cdot 1 = \frac{1}{6^s}$
• $\frac{1}{4^s} \cdot \frac{1}{9^s} \cdot \frac{1}{5^s} = \frac{1}{180^s}$

Because of the Fundamental Theorem of Arithmetic (which states every integer has a unique prime factorization), every single $1/n^s$ appears exactly once.

# Conceptual logic of the Euler Product expansion
def expand_zeta(primes, s):
    # Each prime p contributes a geometric series: 1 + p^-s + p^-2s...
    # The product of these series equals the sum of all n^-s
    return "Sum of all 1/n^s"

This gives us the Euler Product Formula: $\sum_{n=1}^{\infty} \frac{1}{n^s} = \prod_{p\text{ prime}} \left(1 + \frac{1}{p^s} + \frac{1}{p^{2s}} + \frac{1}{p^{3s}} + \cdots\right)$

---

Simplifying via Geometric Series

The terms inside the product are geometric series. We can simplify them using a standard algebraic trick. If we let $X$ be the sum: $X = 1 + \frac{1}{p^s} + \frac{1}{p^{2s}} + \cdots$ Then: $\frac{X}{p^s} = \frac{1}{p^s} + \frac{1}{p^{2s}} + \frac{1}{p^{3s}} + \cdots$

Subtracting the two: $X - \frac{X}{p^s} = 1 \implies X\left(1 - \frac{1}{p^s}\right) = 1 \implies X = \frac{1}{1 - \frac{1}{p^s}}$

Thus, the zeta function becomes: $\zeta(s) = \prod_{p\text{ prime}} \frac{1}{1 - p^{-s}}$

The Bridge to Calculus

This formula is a critical "bridge." It connects Number Theory (the study of primes) with Calculus (the study of infinite series and continuous functions).

Crossing the Bridge

To solve the original problem—the sum of the reciprocals of primes—we need to turn the product into a sum. The standard tool for this is the logarithm.

Proof Roadmap:

• Define $\zeta(s)$ as a sum.
• Convert $\zeta(s)$ to a product over primes.
• Simplify using geometric series.
• Apply logarithms to convert the product back to a sum.
• Use Taylor series to isolate the prime reciprocals.

By taking the log of the Euler product, we get: $\log \zeta(s) = \sum_{p \text{ prime}} \log \left( \frac{1}{1 - \frac{1}{p^s}} \right)$

While these logarithms look daunting, they can be further simplified using the Taylor series expansion for $\log(\frac{1}{1-x}) \dots$