An analytic technique using Fourier integrals over the unit circle to count solutions to additive equations.
IntuitionIntuition: counting by walking around a circle
Suppose you want to count the number of ways to write an integer n as a sum of k elements from a set A of integers (for example, primes, or perfect k-th powers). Encode A by a "wave" running once around the unit circle: for each real number α in [0,1), form the exponential sum FA(α)=∑a∈Ae(aα), where e(x):=e2πix. As α sweeps from 0 to 1, FA(α)k oscillates; near a handful of "resonant" points — rationals a/q with small denominator q — the terms of the sum line up and add constructively (a major arc), while almost everywhere else the phases point in essentially random directions and cancel out (a minor arc). The circle method turns this geometric picture into an exact formula and then estimates it arc by arc.
Interactive unit circle showing a point at angle theta used to build the exponential sum of the circle method.
Dragging θ moves the point e2πiθ/360 around the unit circle; the circle method integrates such points against FA(α)k as θ/360=α runs over [0,1).
UndergraduateDefinition: the generating function on the circle
Definition: Exponential sum and representation count
For a finite set A⊂Z≥0 write FA(α)=∑a∈Ae(aα) with e(x):=e2πix. For n≥0 and k≥1, let rk(n)=#{(a1,…,ak)∈Ak:a1+⋯+ak=n} be the number of ordered representations of n as a sum of k elements of A.
FA(α)=a∈A∑e(aα),e(x):=e2πix
Raising FA to the k-th power and expanding, the coefficient of e(nα) "should" be exactly rk(n); the extraction identity below (Theorem 1) makes this precise: rk(n)=∫01FA(α)ke(−nα)dα. Here FA(α)k is the wave from the picture above, and the integral picks out the one frequency, n, that we care about — this is the entire content of the method: replace a counting problem by an integral estimation problem.
rk(n)=∫01FA(α)ke(−nα)dα
Major arcs versus minor arcs
Feature
Major arc (near a/q, q small)
Minor arc (elsewhere)
Where it sits
A short interval around each rational a/q with q≤Q
What is left of [0,1) after removing all major arcs
Size of FA(α)
Close to its trivial maximum ∣A∣: terms reinforce
Expected to be much smaller than ∣A∣: terms cancel
Role in the estimate
Produces the main term (a "singular series")
Must be bounded above and absorbed into the error term
For every finite A⊂Z≥0, integer k≥1 and n≥0, rk(n)=∫01FA(α)ke(−nα)dα, where rk(n) is the number of ordered k-tuples from A summing to n.
Why is it true?
This is what lets us replace a purely combinatorial counting problem by a question about the size of an analytic integral: if we can show the integral is positive, a representation must exist.
Proof
First note the orthogonality relation ∫01e(mα)dα=1 if m=0, and =0 if m=0: writing e(mα)=cos(2πmα)+isin(2πmα), for m=0 this is a full number of periods of a sine/cosine wave over [0,1), which integrates to 0; for m=0 the integrand is the constant 1.
Now expand the k-th power of the generating function: FA(α)k=(∑a∈Ae(aα))k=∑(a1,…,ak)∈Ake((a1+⋯+ak)α), a finite sum over all ordered k-tuples from A, grouped by the exponent m=a1+⋯+ak.
Multiply both sides by e(−nα) and integrate term by term over [0,1) (legitimate because the sum is finite, so integration and summation commute): ∫01FA(α)ke(−nα)dα=∑(a1,…,ak)∈Ak∫01e((a1+⋯+ak−n)α)dα.
By the orthogonality relation, each summand on the right is 1 exactly when a1+⋯+ak=n and 0 otherwise. So the whole sum collapses to exactly the count of tuples with a1+⋯+ak=n, i.e. rk(n). This proves the identity exactly, with no approximation involved: it is an identity, not an asymptotic.
For every real α and every integer Q≥1, there is a rational a/q with 1≤q≤Q, gcd(a,q)=1, such that α−qa≤q(Q+1)1.
Why is it true?
This guarantees that every point on the circle sits close to some low-denominator rational, which is exactly what lets us partition [0,1) into major arcs — short intervals around the good approximations a/q with q small, where FA is large — and minor arcs, the leftover set.
Proof
Consider the Q+2 numbers 0,{α},{2α},…,{Qα},1, where {x} denotes the fractional part of x; all lie in [0,1]. Partition [0,1] into Q+1 equal subintervals of length 1/(Q+1): [0,Q+11),[Q+11,Q+12),….
We have Q+2 numbers and only Q+1 subintervals, so by the pigeonhole principle two of the numbers {jα} and {iα} (with 0≤i<j≤Q, allowing i=0 so {iα}=0) land in the same subinterval, hence differ by less than 1/(Q+1): ∣{jα}−{iα}∣<Q+11.
Set q=j−i, so 1≤q≤Q. Since {jα}−{iα}=(jα−iα)−(⌊jα⌋−⌊iα⌋)=qα−a for the integer a=⌊jα⌋−⌊iα⌋, we get ∣qα−a∣<Q+11, i.e. α−qa<q(Q+1)1.
Finally, if gcd(a,q)=d>1, divide both a and q by d: the resulting fraction has an even smaller denominator and the same (or a smaller) distance to α, so we may take gcd(a,q)=1 without loss of generality. This completes the proof; it is a finite, constructive pigeonhole argument, with no appeal to any unproved estimate.
UndergraduateReal-World Applications and Worked Examples
Beyond number theory, the same "add up waves and look for resonance" idea is the working principle of the discrete Fourier transform used throughout signal processing and electrical engineering, and the asymptotic p(n)∼4n31exp(π32n) that Hardy and Ramanujan extracted with an early version of the circle method (applied to the partition function p(n)) is used in statistical mechanics to estimate the number of microstates — hence the entropy — of a system of indistinguishable bosonic excitations at fixed total energy n.
Example: A miniature circle method computation
Let A={1,2,…,8}. Using the extraction identity, compute r2(9), the number of ordered pairs (a,b)∈A2 with a+b=9.
Solution
By Theorem 1, r2(9)=∫01FA(α)2e(−9α)dα with FA(α)=∑a=18e(aα); there is no need to actually evaluate the integral analytically — the identity guarantees it equals the direct count, so we may compute the count combinatorially instead and trust the identity.
List every ordered pair (a,b) with a,b∈{1,…,8} and a+b=9: (1,8),(2,7),(3,6),(4,5),(5,4),(6,3),(7,2),(8,1). Each first coordinate from 1 to 8 determines b=9−a uniquely, and b always lands back in {1,…,8} (since 1≤a≤8⇒1≤9−a≤8), so all 8 values of a work.
Hence r2(9) equals 8. This tiny example is exactly the identity of Theorem 1 in action: the analytic integral and the combinatorial count are, by construction, one and the same number — the real content of the method appears only once A becomes an infinite or growing set (like the primes) and one must estimate rather than enumerate.
Example: Best rational approximation via Dirichlet's theorem
Apply Dirichlet's approximation theorem with α=2 and Q=5 to exhibit a fraction a/q, 1≤q≤5, satisfying the bound of Theorem 2, and check it numerically.
Solution
We want a/q with 1≤q≤5 and ∣2−a/q∣<1/(q⋅6) (taking Q+1=6 in the theorem's bound). Try q=5: the nearest integer to 52≈7.0711 is a=7, giving 7/5=1.4.
Check the bound: ∣2−7/5∣=∣1.41421…−1.4∣≈0.01421, while the theorem promises 1/(5⋅6)=1/30≈0.0333; indeed 0.01421<0.0333, so the bound holds comfortably, exactly as guaranteed.
This 7/5 is in fact the convergent of the continued fraction of 2=[1;2,2,2,…] obtained after two steps, which is why it approximates so well — Dirichlet's pigeonhole proof does not need continued fractions to work, but it always lands on comparably good approximants. So 57 is a valid witness.
In the extraction identity rk(n)=∫01FA(α)ke(−nα)dα, why does the integral equal exactly rk(n)?
For A={1,2,…,8}, what is r2(9), the number of ordered pairs (a,b)∈A2 with a+b=9?
According to Dirichlet's approximation theorem, for every real α and integer Q≥1, what is guaranteed to exist?
In statistical mechanics, what does the Hardy–Ramanujan asymptotic p(n)∼4n31exp(π32n) for the partition function help estimate?