# Goldston–Pintz–Yıldırım sieve

The Goldston–Pintz–Yıldırım sieve (also called GPY sieve or GPY method) is a sieve method and variant of the Selberg sieve with generalized, multidimensional sieve weights. The sieve led to a series of important breakthroughs in analytic number theory.
It is named after the mathematicians Dan Goldston, János Pintz and Cem Yıldırım. They used it in 2005 to show that there are infinitely many prime tuples whose distances are arbitrarily smaller than the average distance that follows from the prime number theorem.
The sieve was then modified by Yitang Zhang in order to prove a finite bound on the smallest gap between two consecutive primes that is attained infinitely often. Later the sieve was again modified by James Maynard (who lowered the bound to $600$) and by Terence Tao.

## Goldston–Pintz–Yıldırım sieve
### Notation

Fix a $k\in \mathbb {N}$ and the following notation:
* $\mathbb {P}$ is the set of prime numbers and $1_{\mathbb {P} }(n)$ the characteristic function of that set,
* $\Lambda (n)$ is the von Mangoldt function,
* $\omega (n)$ is the small prime omega function (which counts the distinct prime factors of $n$)
* ${\mathcal {H}}=\{h_{1},\dots ,h_{k}\}$ is a set of distinct nonnegative integers $h_{i}\in \mathbb {Z} _{+}\cup \{0\}$.
* $\theta (n)$ is another characteristic function of the primes defined as
       $\theta (n)={\begin{cases}\log(n)&{\text{if }}n\in \mathbb {P} \\0&{\text{else.}}\end{cases}}$
    Notice that $\theta (n)=\log((n-1)1_{\mathbb {P} }(n)+1)$.
For an ${\mathcal {H}}$ we also define
* ${\mathcal {H}}(n):=(n+h_{1},\dots ,n+h_{k})$,
* $P_{\mathcal {H}}(n):=(n+h_{1})(n+h_{2})\cdots (n+h_{k})$
* $\nu _{p}({\mathcal {H}})$ is the amount of distinct residue classes of ${\mathcal {H}}$ modulo $p$. For example $\nu _{3}(\{0,2,4\})=3$ and $\nu _{3}(\{0,2\})=2$ because $\{0,2,4\}{\stackrel {\pmod {3}}{=}}\{0,1,2\}$ and $\{0,2\}{\stackrel {\pmod {3}}{=}}\{0,2\}$.
If $\nu _{p}({\mathcal {H}})<k$ for all $p\in \mathbb {P}$, then we call ${\mathcal {H}}$ admissible.

### Construction

Let ${\mathcal {H}}=\{h_{1},\dots ,h_{k}\}$ be admissible and consider the following sifting function
    ${\mathcal {S}}(N,c;{\mathcal {H}}):=\sum \limits _{n=N+1}^{2N}\left(\sum \limits _{h_{i}\in {\mathcal {H}}}1_{\mathbb {P} }(n+h_{i})-c\right)w(n)^{2},\quad w(n)\in \mathbb {R} ,\quad c>0,$
where $w(n)$ is a weight function we derive later.
For each $n\in [N+1,2N]$ this sifting function counts the primes of the form $n+h_{i}$ minus some threshold $c$, so if ${\mathcal {S}}>0$ then there exist some $n$ such that at least $\lfloor c\rfloor +1$ are prime numbers in ${\mathcal {H}}(n)$.
Since $1_{\mathbb {P} }(n)$ has not so nice analytic properties one chooses rather the following sifting function
    ${\mathcal {S}}(N;{\mathcal {H}}):=\sum \limits _{n=N+1}^{2N}\left(\sum \limits _{h_{i}\in {\mathcal {H}}}\theta (n+h_{i})-\log(3N)\right)w(n)^{2}.$
Since $\log(N)<\theta (n+h_{i})<\log(2N)$ and $c=\log(3n)$, we have ${\mathcal {S}}>0$ only if there are at least two prime numbers $n+h_{i}$ and $n+h_{j}$. Next we have to choose the weight function $w(n)$ so that we can detect prime k-tuples.

#### Derivation of the weights

A candidate for the weight function is the generalized von Mangoldt function
    $\Lambda _{k}(n)=\sum \limits _{d\mid n}\mu (d)\left(\log \left({\frac {n}{d}}\right)\right)^{k},$
which has the following property: if $\omega (n)>k$, then $\Lambda _{k}(n)=0$. This functions also detects factors which are proper prime powers, but this can be removed in applications with a negligible error.$^{:826}$
So if ${\mathcal {H}}(n)$ is a prime k-tuple, then the function
    $\Lambda _{k}(n;{\mathcal {H}})={\frac {1}{k!}}\Lambda _{k}(P_{\mathcal {H}}(n))$
will not vanish. The factor $1/k!$ is just for computational purposes. The (classical) von Mangoldt function can be approximated with the truncated von Mangoldt function
    $\Lambda (n)\approx \Lambda _{R}(n):=\sum \limits _{\begin{array}{c}d\mid n\\d\leq R\end{array}}\mu (d)\log \left({\frac {R}{d}}\right),$
where $R$ now no longer stands for the length of ${\mathcal {H}}$ but for the truncation position. Analogously we approximate $\Lambda _{k}(n;{\mathcal {H}})$ with
    $\Lambda _{R}(n;{\mathcal {H}})={\frac {1}{k!}}\sum \limits _{\begin{array}{c}d\mid P_{\mathcal {H}}(n)\\d\leq R\end{array}}\mu (d)\left(\log \left({\frac {R}{d}}\right)\right)^{k}$
For technical purposes we rather want to approximate tuples with primes in multiple components than solely prime tuples and introduce another parameter $0\leq \ell \leq k$ so we can choose to have $k+\ell$ or less distinct prime factors. This leads to the final form
    $\Lambda _{R}(n;{\mathcal {H}},\ell )={\frac {1}{(k+\ell )!}}\sum \limits _{\begin{array}{c}d\mid P_{\mathcal {H}}(n)\\d\leq R\end{array}}\mu (d)\left(\log \left({\frac {R}{d}}\right)\right)^{k+\ell }$
Without this additional parameter $\ell$ one has for a distinct $d=d_{1}d_{2}\cdots d_{k}$ the restriction $d_{1}\leq R,d_{2}\leq R,\dots ,d_{k}\leq R$ but by introducing this parameter one gets the more looser restriction $d_{1}d_{2}\dots d_{k}\leq R$.$^{:827}$ So one has a $k+\ell$-dimensional sieve for a $k$-dimensional sieve problem.

### Goldston–Pintz–Yıldırım sieve

The GPY sieve has the following form
    ${\mathcal {S}}(N;{\mathcal {H}},\ell ):=\sum \limits _{n=N+1}^{2N}\left(\sum \limits _{h_{i}\in {\mathcal {H}}}\theta (n+h_{i})-\log(3N)\right)\Lambda _{R}(n;{\mathcal {H}},\ell )^{2},\qquad |{\mathcal {H}}|=k$
with
    $\Lambda _{R}(n;{\mathcal {H}},\ell )={\frac {1}{(k+\ell )!}}\sum \limits _{\begin{array}{c}d\mid P_{\mathcal {H}}(n)\\d\leq R\end{array}}\mu (d)\left(\log \left({\frac {R}{d}}\right)\right)^{k+\ell },\quad 0\leq \ell \leq k$.$^{:827–829}$

### Proof of the main theorem by Goldston, Pintz and Yıldırım

Consider $({\mathcal {H}}_{1},\ell _{1},k_{1})$ and $({\mathcal {H}}_{2},\ell _{2},k_{2})$ and $1\leq h_{0}\leq R$ and define $M:=k_{1}+k_{2}+\ell _{1}+\ell _{2}$. In their paper, Goldston, Pintz and Yıldırım proved in two propositions that under suitable conditions two asymptotic formulas of the form
    $\sum \limits _{n\leq N}\Lambda _{R}(n;{\mathcal {H}}_{1},\ell _{1})\Lambda _{R}(n;{\mathcal {H}}_{2},\ell _{2})=C_{1}\left({\mathcal {S}}({\mathcal {H}}^{i})+o_{M}(1)\right)N$
and
    $\sum \limits _{n\leq N}\Lambda _{R}(n;{\mathcal {H}}_{1},\ell _{1})\Lambda _{R}(n;{\mathcal {H}}_{2},\ell _{2})\theta (n+h_{0})=C_{2}\left({\mathcal {S}}({\mathcal {H}}^{j})+o_{M}(1)\right)N$
hold, where $C_{1},C_{2}$ are two constants, ${\mathcal {S}}({\mathcal {H}}^{i})$ and ${\mathcal {S}}({\mathcal {H}}^{j})$ are two singular series whose description we omit here.
Finally one can apply these results to ${\mathcal {S}}$ to derive the theorem by Goldston, Pintz and Yıldırım on infinitely many prime tuples whose distances are arbitrarily smaller than the average distance.$^{:827–829}$