The large sieve inequality is a fundamental tool in analytic number theory, bounding exponential sums over well-spaced frequencies. It underlies the Bombieri-Vinogradov theorem and many results on the distribution of primes.
Statement
Theorem6.3Large Sieve Inequality
Let aM+1β,β¦,aM+Nβ be complex numbers and let Ξ±1β,β¦,Ξ±RββR satisfy β₯Ξ±rββΞ±sββ₯β₯Ξ΄>0 for rξ =s (where β₯β β₯ denotes distance to the nearest integer). Then:
βr=1Rβββn=M+1M+Nβanβe2ΟinΞ±rββ2β€(Nβ1+Ξ΄β1)βn=M+1M+Nββ£anββ£2
We seek a trigonometric polynomial T(Ξ±)=ββ£kβ£β€Kβckβe2ΟikΞ± with K=βΞ΄β1β such that T(Ξ±rβ)=1 for all r and β«01ββ£T(Ξ±)β£2dΞ± is minimized.
By duality, the large sieve is equivalent to: there exist b1β,β¦,bRβ with ββ£brββ£2β€1 such that ββrβbrβe2ΟinΞ±rβββ22ββ€N+Ξ΄β1β1.
Direct proof via Beurling-Selberg functions. Consider the Hermitian form βr,sβBrsβzrβzΛsβ where Brsβ=βn=M+1M+Nβe2Οin(Ξ±rββΞ±sβ). The large sieve is equivalent to showing the operator norm β₯Bβ₯β€N+Ξ΄β1β1.
For the diagonal: Brrβ=N. For off-diagonal terms, β£Brsββ£=β£sin(ΟN(Ξ±rββΞ±sβ))/sin(Ο(Ξ±rββΞ±sβ))β£.
Hilbert-type inequality. The key estimate: βrξ =sββ£sin(Ο(Ξ±rββΞ±sβ))β£β£zrββ£β£zsββ£ββ€(Ξ΄β1β1)ββ£zrββ£2 when β₯Ξ±rββΞ±sββ₯β₯Ξ΄. This follows because the quadratic form βrξ =sβzrβzΛsβ/sin(Ο(Ξ±rββΞ±sβ)) is bounded by (Ξ΄β1β1) times the identity, which can be proved using the Beurling-Selberg majorant/minorant for the characteristic function of an interval.
Combining: βrββ£S(Ξ±rβ)β£2=βr,sβBrsβzrβzΛsββ€(N+Ξ΄β1β1)ββ£anββ£2 where we identified zrβ=βanβe2ΟinΞ±rβ... more precisely, expanding βrββ£S(Ξ±rβ)β£2 and using the bound on the Hermitian form.
Alternatively (Montgomery-Vaughan): βrββ£S(Ξ±rβ)β£2β€βnββ£anββ£2βrββββ£kβ£β€Kβck(r)βe2Οiknβ where ck(r)β are chosen optimally, leading to the stated bound with constant Nβ1+Ξ΄β1. β‘
β
ExampleLarge Sieve for Characters
Specializing to Dirichlet characters: Ξ±rβ=a/q for Farey fractions with qβ€Q. The spacing Ξ΄=1/Q2 (Farey fractions are spaced at least 1/Q2 apart). The character form of the large sieve: βqβ€QβΟ(q)qββΟmodqββββn=M+1M+NβanβΟ(n)β2β€(N+Q2)ββ£anββ£2.
RemarkSharpness
The constant Nβ1+Ξ΄β1 is the best possible (it cannot be replaced by max(N,Ξ΄β1) in general). The extremal case is achieved by an arithmetic progression of frequencies Ξ±rβ=rΞ΄.