The Golden Ratio Is Hiding in Projective Planes


collinear1 collinear out of 5 others = 1/5 011111001101010110000100𝔽₂³
-2-10121/Ο†βˆ’Ο†0.618-1.618xΒ² + 1.00x βˆ’ 1 = 0Ξ± = 1.000
ΞΆ Ο†

Click points on either view. Drag to interpolate between ΞΆ and Ο†.

Left: The Fano plane PG(2,F2)\mathrm{PG}(2, \mathbb{F}_2) β€” 7 points, 7 lines, 3 points per line. Click any two points to highlight the unique line through them: there’s always exactly 1 collinear point out of 5 remaining, giving the ratio 1/5=1/(p2+pβˆ’1)1/5 = 1/(p^2 + p - 1) at p=2p = 2.

Right: The vector space F23\mathbb{F}_2^3 that the Fano plane comes from. Each nonzero vertex of the cube is a projective point; the origin 000000 is the one that gets quotiented out. The gold triangles are the Fano lines β€” three nonzero vectors that XOR to zero.

Below: The polynomial x2+Ξ±xβˆ’1x^2 + \alpha x - 1. At Ξ±=0\alpha = 0 its roots are Β±1\pm 1 (Riemann zeta territory). At Ξ±=1\alpha = 1 they land at 1/Ο†1/\varphi and βˆ’Ο†-\varphi β€” the golden ratio, forced by the combinatorics of the projective plane. Drag the slider to watch the roots move; press Play to animate.

I’ve been playing around with projective planes over finite fields and stumbled into something I wasn’t expecting: the golden ratio shows up in a completely natural way, and it drags Lucas numbers and Euler products along with it. Here’s the story.

The setup

Take a prime pp and build the projective plane PG(2,Fp)\mathrm{PG}(2, \mathbb{F}_p) β€” the set of one-dimensional subspaces of Fp3\mathbb{F}_p^3. It has p2+p+1p^2 + p + 1 points, the same number of lines, and every line passes through p+1p + 1 points. Standard stuff.

Now ask a simple question: pick two distinct points PP and QQ. They determine a unique line. If you pick a third point RR at random from the remaining p2+pβˆ’1p^2 + p - 1 points, what’s the probability that RR lands on the same line?

There are pβˆ’1p - 1 other points on the line through PP and QQ (the p+1p + 1 on the line, minus PP and QQ), so the collinearity probability is

collp=pβˆ’1p2+pβˆ’1.\mathrm{coll}_p = \frac{p - 1}{p^2 + p - 1}.

That denominator, p2+pβˆ’1p^2 + p - 1, is where things get interesting.

The golden ratio appears

The polynomial x2+xβˆ’1x^2 + x - 1 has roots

x=βˆ’1Β±52,x = \frac{-1 \pm \sqrt{5}}{2},

which are 1/Ο†1/\varphi and βˆ’Ο†-\varphi, where Ο†=(1+5)/2\varphi = (1 + \sqrt{5})/2 is the golden ratio. So we get an exact factorization:

p2+pβˆ’1=(pβˆ’1Ο†)(p+Ο†).p^2 + p - 1 = \left(p - \frac{1}{\varphi}\right)(p + \varphi).

This isn’t an approximation. It’s algebraically exact for every prime pp. The golden ratio just… shows up, uninvited, from the combinatorics of projective planes.

Building an Euler product

With an interesting polynomial in hand, the natural move is to build an Euler product β€” take the product over all primes:

PPG(2)=∏pp2+pβˆ’1p2=∏p(1+1pβˆ’1p2).\mathcal{P}_{\mathrm{PG}(2)} = \prod_p \frac{p^2 + p - 1}{p^2} = \prod_p \left(1 + \frac{1}{p} - \frac{1}{p^2}\right).

This product diverges, but in a controlled way. Expanding the logarithm:

log⁑(1+1pβˆ’1p2)=1pβˆ’32p2+O(pβˆ’3),\log\left(1 + \frac{1}{p} - \frac{1}{p^2}\right) = \frac{1}{p} - \frac{3}{2p^2} + O(p^{-3}),

so the divergence comes entirely from βˆ‘p1/p\sum_p 1/p (which diverges by Mertens’ theorem), and everything else converges to a constant Cβ‰ˆβˆ’0.5323C \approx -0.5323.

primes 6
$p$$1 + 1/p - 1/p^2$running $\prod$
21.2500001.250000
31.2222221.527778
51.1600001.772222
71.1224491.989229
111.0826452.153628
131.0710062.306548
0.51.01.5Ξ£ 1/p (diverges)remainder (converges)23571113

The running product grows without bound β€” but the log splits into a divergent piece (Ξ£ 1/p) and a convergent remainder.

Factoring into convergent pieces

The golden ratio factorization of p2+pβˆ’1p^2 + p - 1 splits the Euler product into two pieces:

∏pp2+pβˆ’1p2=∏p(1βˆ’1Ο†p)Γ—βˆp(1+Ο†p).\prod_p \frac{p^2 + p - 1}{p^2} = \prod_p \left(1 - \frac{1}{\varphi p}\right) \times \prod_p \left(1 + \frac{\varphi}{p}\right).

Both individual products diverge (the first to 00, the second to +∞+\infty), but you can tame them. The trick is a fractional-power Mertens normalization β€” you compare each factor against the corresponding power of the Riemann zeta function’s Euler product. This gives two convergent products:

C1=∏p1βˆ’1/(Ο†p)(1βˆ’1/p)1/Ο†β‰ˆ1.0956,\mathcal{C}_1 = \prod_p \frac{1 - 1/(\varphi p)}{(1 - 1/p)^{1/\varphi}} \approx 1.0956,

C2=∏p1+Ο†/p(1+1/p)Ο†β‰ˆ0.8745.\mathcal{C}_2 = \prod_p \frac{1 + \varphi/p}{(1 + 1/p)^{\varphi}} \approx 0.8745.

That second one is suspiciously close to 7/87/8. Within 0.06%0.06\%. And 7=L47 = L_4, the fourth Lucas number. Coincidence? Maybe.

Lucas numbers everywhere

The real punchline is the constant CC. It has an exact series representation:

C=βˆ‘n=2∞(βˆ’1)n+1LnnP(n),C = \sum_{n=2}^{\infty} (-1)^{n+1} \frac{L_n}{n} P(n),

where LnL_n is the nn-th Lucas number and P(n)=βˆ‘ppβˆ’nP(n) = \sum_p p^{-n} is the prime zeta function. Written out:

C=βˆ’32P(2)+43P(3)βˆ’74P(4)+115P(5)βˆ’186P(6)+β‹―C = -\frac{3}{2}P(2) + \frac{4}{3}P(3) - \frac{7}{4}P(4) + \frac{11}{5}P(5) - \frac{18}{6}P(6) + \cdots

The Lucas numbers 3,4,7,11,18,29,…3, 4, 7, 11, 18, 29, \ldots show up as numerators. This isn’t a coincidence either β€” it comes from the power series expansion of log⁑(1+tβˆ’t2)\log(1 + t - t^2), and the fact that 1+tβˆ’t21 + t - t^2 has roots at Ο†\varphi and βˆ’1/Ο†-1/\varphi means the Binet formula naturally produces Lucas numbers as the coefficients.

terms 4
C-0.50.00.5n=2n=3n=4n=5n=6n=7n=8n=9n=10n=11n=12n=13n=14n=15partial sum = -0.4987
-3/2
L2 = 3
+4/3
L3 = 4
-7/4
L4 = 7
+11/5
L5 = 11
-18/6
L6 = 18
+29/7
L7 = 29
-47/8
L8 = 47
+76/9
L9 = 76
-123/10
L10 = 123
+199/11
L11 = 199
-322/12
L12 = 322
+521/13
L13 = 521
-843/14
L14 = 843
+1364/15
L15 = 1364

Each term's numerator is a Lucas number. The partial sums oscillate and converge to C β‰ˆ -0.5284.

Where does regular zeta fit?

The Riemann zeta function ΞΆ(s)=∏p(1βˆ’pβˆ’s)βˆ’1\zeta(s) = \prod_p (1 - p^{-s})^{-1} comes from the polynomial p2βˆ’1=(pβˆ’1)(p+1)p^2 - 1 = (p-1)(p+1), which counts the multiplicative group Fp2βˆ—\mathbb{F}_{p^2}^*. The roots are Β±1\pm 1 β€” rational.

Our projective plane polynomial is p2+pβˆ’1p^2 + p - 1, with irrational roots involving Ο†\varphi. The difference between the two polynomials is just pp, which represents the contribution from the line at infinity in projective geometry.

You can actually interpolate between them. Consider fΞ±(x)=x2+Ξ±xβˆ’1f_\alpha(x) = x^2 + \alpha x - 1 for α∈[0,1]\alpha \in [0, 1]. At Ξ±=0\alpha = 0 the roots are Β±1\pm 1 (zeta territory), and at Ξ±=1\alpha = 1 they’re 1/Ο†1/\varphi and βˆ’Ο†-\varphi (golden territory). The roots slide continuously from one to the other.

And there’s a uniqueness result: Ξ±=1\alpha = 1 is the only value for which the associated constant C(Ξ±)C(\alpha) converges. For any other Ξ±\alpha, you get a leftover (Ξ±βˆ’1)βˆ‘p1/p(\alpha - 1)\sum_p 1/p term that diverges. The golden ratio isn’t just a coincidence from projective geometry β€” it’s analytically forced.

The Fano plane

The smallest case is p=2p = 2: the Fano plane, with 7 points and 7 lines. Here p2+pβˆ’1=5p^2 + p - 1 = 5, a Fibonacci number, and the point count 77 is a Lucas number. Both the point count and the collinearity denominator land in the Fibonacci-Lucas family. For larger primes this stops happening pretty quickly, which makes the Fano plane a nice minimal example where the golden ratio structure is maximally visible.

Open questions

A few things I haven’t resolved:

  • Is C2=7/8\mathcal{C}_2 = 7/8 exactly, or is there a tiny correction?
  • Do these geometric zeta functions have functional equations like the Riemann zeta function?
  • What happens for higher-dimensional projective spaces PG(n,Fp)\mathrm{PG}(n, \mathbb{F}_p) with n>2n > 2? What algebraic numbers show up?
  • Is there an analogue of the Riemann hypothesis here?

Whether any of this leads somewhere deep or is just a nice curiosity, I’m not sure yet. But the golden ratio emerging from projective plane combinatorics and then pulling in Lucas numbers and Euler products feels like it’s pointing at something.