For a geometrically finite, Zariski-dense subgroup Γ\Gamma of the integer points of the Pythagorean group SO(x2+y2−z2)\mathrm{SO}(x^2+y^2-z^2) with no nontrivial unipotent elements and critical exponent δ\delta sufficiently close to 1, the set SS of hypotenuses of the Pythagorean triples in a thin orbit satisfies #(S∩[−N,N])∼CNδ\#(S\cap[-N,N])\sim C N^\delta. Bourgain and Kontorovich proved, by a dispersion method that bypasses the spectral-gap barrier of the classical affine sieve, that SS admits an exponent of distribution α=7/24−ε\alpha=7/24-\varepsilon for any ε>0\varepsilon>0 whenever δ>δ0(ε)\delta>\delta_0(\varepsilon), yielding infinitely many hypotenuses with at most R=4R=4 prime factors. The paper leaves open whether the exponent can be pushed past 7/247/24 — concretely to 7/24+ε07/24+\varepsilon_0 (Remarks 1.9 and 5.18) — how far toward an Elliott–Halberstam-type range the level of distribution of SS can go (Remark 1.19), and whether the requirement that δ\delta be extremely close to 1 can be relaxed. No improvement of the exponent 7/247/24 is known in the literature to date.

Background

A Pythagorean triple is an integral point on the cone F(x,y,z)=x2+y2−z2=0F(x,y,z)=x^2+y^2-z^2=0, and its hypotenuse is the coordinate zz. The connected component of the identity, G=SO∘(2,1)G=\mathrm{SO}^\circ(2,1), of the special orthogonal group of FF preserves integrality on this cone, so for a fixed triple x0x_0 (say (1,0,1)(1,0,1)) and a subgroup Γ<G(Z)\Gamma<G(\mathbb{Z}) one can study the orbit O=x0⋅Γ\mathcal{O}=x_0\cdot\Gamma and the set SS of hypotenuses f(x,y,z)=zf(x,y,z)=z of triples in O\mathcal{O}. The orbit is called thin when Γ\Gamma has infinite index in the Z\mathbb{Z}-points of its Zariski closure — the smallest algebraic variety (zero set of polynomial equations in the matrix entries) containing it; Γ\Gamma is Zariski-dense when that closure is the whole ambient group. Of interest here are finitely generated, geometrically finite (admitting a fundamental region in the hyperbolic plane with finitely many sides), Zariski-dense such Γ\Gamma with no unipotent elements besides the identity, a unipotent element being a matrix all of whose eigenvalues equal 1 — in this two-dimensional setting, a parabolic element fixing a single point on the circle at infinity.

The group carries a critical exponent δ=δΓ∈(0,1]\delta=\delta_\Gamma\in(0,1]: the abscissa of convergence of its Poincaré series, ∑γ∈Γe−s d(i,γi)\sum_{\gamma\in\Gamma}e^{-s\,d(i,\gamma i)} with dd the hyperbolic distance on the upper half-plane and ii the base point; equivalently, it is the Hausdorff dimension of the limit set, the set of accumulation points of any orbit of Γ\Gamma on the boundary circle. Such a group is a lattice — a discrete subgroup of finite hyperbolic covolume — exactly when δ=1\delta=1, and the orbit is thinner the smaller δ\delta is. For δ>1/2\delta>1/2 one has the counting law #(S∩[−N,N])∼CNδ\#(S\cap[-N,N])\sim C N^\delta (Kontorovich; Kontorovich–Oh), so a thin orbit (δ<1\delta<1) produces a genuinely sparse set of hypotenuses. The affine sieve of Bourgain, Gamburd, and Sarnak extends the classical sieve to such orbits, studying the prime factorization of polynomial values on them. One weights the points of O\mathcal{O} by nonnegative coefficients supported on S∩[−N,N]S\cap[-N,N], computes for each modulus qq a main term governed by local densities β(q)\beta(q) (here of linear-sieve type, β(p)≈1/p\beta(p)\approx 1/p on average) and a remainder r(q)r(q), and declares a level of distribution QQ when the total error ∑q<Q∣r(q)∣\sum_{q<Q}|r(q)| is smaller than the total mass by a power saving; writing Q=NαQ=N^\alpha, the number α\alpha is the exponent of distribution. Larger α\alpha feeds the weighted sieve — the classical combinatorial sieve augmented by multiplicative weights attached to the prime factorization of the surviving elements, which sharpens the almost-prime output — and produces points whose ff-value has fewer prime factors (an RR-almost prime is an integer with at most RR prime factors; the least such RR for a sequence is its saturation number).

The standard route to a level of distribution is via the spectral gap of Γ\Gamma acting on congruence covers, the quotients of the hyperbolic plane by the subgroups obtained by reducing the matrix entries of Γ\Gamma modulo qq; the gap is conventionally encoded by a parameter θ\theta, with θ\theta closer to 1/21/2 meaning a stronger gap. Gamburd's gap θ=5/6\theta=5/6 for infinite-index subgroups of SL2(Z)\mathrm{SL}_2(\mathbb{Z}) gives α<1/12\alpha<1/12 and hence R=13R=13 for thin hypotenuses (Kontorovich; Kontorovich–Oh), and even a best-possible gap θ=1/2\theta=1/2 would cap this route at α<1/4\alpha<1/4, giving R=5R=5. Because the spin double cover — the two-to-one covering homomorphism SL2(R)→SO∘(2,1)\mathrm{SL}_2(\mathbb{R})\to\mathrm{SO}^\circ(2,1) — is quadratic, a ball of size NN in the orbit corresponds to matrices of size TT with N∼T2N\sim T^2, so α=1/4\alpha=1/4 corresponds to level Q∼T1/2Q\sim T^{1/2}. Bourgain and Kontorovich call 1/41/4 the natural analogue for this problem of the Bombieri–Vinogradov range (for primes in arithmetic progressions, the level x1/2x^{1/2} of uniform distribution over progressions), and any exponent exceeding 1/41/4 an instance of the Elliott–Halberstam range (the conjectural level x1−εx^{1-\varepsilon} for primes in progressions). Bourgain and Kontorovich broke the spectral-gap barrier unconditionally for thin hypotenuses. Developing bilinear forms and the classical dispersion method (Cauchy–Schwarz bilinearization of incomplete exponential sums) in this thin-group setting, they proved (Theorem 1.6) that for every ε>0\varepsilon>0 there is δ0(ε)<1\delta_0(\varepsilon)<1 such that, whenever δ>δ0(ε)\delta>\delta_0(\varepsilon), the hypotenuse sequence admits exponent of distribution α=7/24−ε\alpha=7/24-\varepsilon; the weighted sieve then yields infinitely many hypotenuses in SS with at most R=4R=4 prime factors (Corollary 1.8). The value 7/247/24 results from a delicate balance of interlocking ranges in the dispersion argument: the modulus sum is split at a height Q0Q_0, and five bilinear sums (denoted S1,…,S5S_1,\dots,S_5) are estimated; the S5S_5 sum is disposed of trivially via ∣S5∣≤∣S2∣|S_5|\le|S_2|, and the authors note (Remarks 1.9 and 5.18) that exploiting the further cancellation in S5S_5 (extra savings of q~η\tilde{q}^{\eta} and (q1q1′)η(q_1q_1')^{\eta}) 'may get the level just a hair higher', i.e. to α=7/24+ε0\alpha=7/24+\varepsilon_0 for a tiny ε0>0\varepsilon_0>0, an improvement for which they state no application is known. They also state that the threshold values δ0(ε)=1−10−15ε\delta_0(\varepsilon)=1-10^{-15}\varepsilon and δ0=1−10−17\delta_0=1-10^{-17} produced by the proof 'can surely be improved' (Remark 1.11), and that the method requires δ\delta extremely close to 1, within the natural baseline δ>1/2\delta>1/2 of the counting and spectral machinery.

Subsequent work has clarified the landscape but not moved the exponent. Hong and Kontorovich discovered that sieving along a larger stabilizer subgroup doubles the exponent of distribution obtainable from expansion alone, improving the almost-prime records for areas and products of coordinates of Pythagorean triples; Ehrman adapted the Bourgain–Kontorovich dispersion machine to those functions and improved the admissible threshold for hypotenuses from 1−10−171-10^{-17} to δ0=0.984\delta_0=0.984, while explicitly not improving the hypotenuse exponent or the value R=4R=4; the companion papers Beyond Expansion II–IV applied the same machinery to low-lying and reciprocal geodesics on the modular surface (the quotient of the hyperbolic plane by SL2(Z)\mathrm{SL}_2(\mathbb{Z})) and to traces of thin semigroups. Recent explicit spectral-gap results for Schottky-type subgroups — finitely generated free groups of hyperbolic isometries (Calderón–Magee; Soares) — strengthen only the expansion route, which, as the source notes, cannot raise the beyond-expansion exponent at all. As of the 2025 literature the exponent 7/247/24 for thin hypotenuses still stands (Li–Liu still quote R=4R=4 as the record), and the ultimate goal of producing primes in thin hypotenuse sets remains out of reach.

Problem Statement

Let F(x,y,z)=x2+y2−z2F(x,y,z)=x^2+y^2-z^2, and let Γ<SOF∘(Z)\Gamma<\mathrm{SO}^\circ_F(\mathbb{Z}) — the group of integer points in the identity component of the special orthogonal group of FF — be finitely generated, geometrically finite (admitting a fundamental region with finitely many sides), and Zariski-dense, with no unipotent elements other than the identity (a unipotent being a matrix all of whose eigenvalues equal 1) and with critical exponent δ∈(1/2,1)\delta\in(1/2,1) (the abscissa of convergence of the Poincaré series ∑γ∈Γe−s d(i,γi)\sum_{\gamma\in\Gamma}e^{-s\,d(i,\gamma i)}, equivalently the Hausdorff dimension of the limit set of Γ\Gamma). Let SS be the set of hypotenuses zz of the Pythagorean triples in the orbit O=x0⋅Γ\mathcal{O}=x_0\cdot\Gamma — via the spin double cover (the two-to-one covering homomorphism SL2(R)→SOF∘(R)\mathrm{SL}_2(\mathbb{R})\to\mathrm{SO}^\circ_F(\mathbb{R})), the values c2+d2c^2+d^2 attained by the bottom rows (c,d)(c,d) of matrices in the lift of Γ\Gamma to SL2(Z)\mathrm{SL}_2(\mathbb{Z}) — so that #(S∩[−N,N])∼CNδ\#(S\cap[-N,N])\sim C N^\delta. Say that SS admits exponent of distribution α\alpha if there exist nonnegative sieve weights supported on S∩[−N,N]S\cap[-N,N], of total mass of order NδN^\delta, whose distribution among residue classes modulo qq has the expected main term with local densities β(q)\beta(q) satisfying ∏w≤p<z(1−β(p))−1≤Clog⁡z/log⁡w\prod_{w\le p<z}(1-\beta(p))^{-1}\le C\log z/\log w and linear-sieve type β(p)≈1/p\beta(p)\approx 1/p, and whose total remainder ∑q<Q∣r(q)∣\sum_{q<Q}|r(q)| is smaller than the total mass by a fixed power saving at level Q=NαQ=N^\alpha — exactly the framework of Theorem 1.6 of Bourgain–Kontorovich. Bourgain and Kontorovich proved that for every ε>0\varepsilon>0 there is δ0(ε)<1\delta_0(\varepsilon)<1 such that SS admits exponent of distribution α=7/24−ε\alpha=7/24-\varepsilon whenever δ>δ0(ε)\delta>\delta_0(\varepsilon), and consequently that SS contains infinitely many R=4R=4 almost primes (integers with at most four prime factors). The problem is to improve the unconditional exponent of distribution of this same sequence; resolving any one of the following linked facets of that single objective, at the strength stated for that facet, resolves the problem. (a) First milestone, from Remarks 1.9 and 5.18 of the source: prove that for some explicit ε0>0\varepsilon_0>0 there is δ0(ε0)<1\delta_0(\varepsilon_0)<1 such that SS admits the exponent α=7/24+ε0\alpha=7/24+\varepsilon_0 whenever δ>δ0(ε0)\delta>\delta_0(\varepsilon_0). (b) Elliott–Halberstam-type range, from Remark 1.19: determine how far past the 'Bombieri–Vinogradov' threshold α=1/4\alpha=1/4 the level of distribution of SS can be pushed and what the true (conjectural) limit for this thin sequence is; a rigorous exponent α>7/24\alpha>7/24, or an established limitation on the achievable level that is strictly stronger than the source's own bounds and that identifies precisely which estimate fails and where, resolves this facet. (c) Range of thinness: the known results require δ\delta extremely close to 1 (the source's δ0=1−10−17\delta_0=1-10^{-17}; the best published value δ0=0.984\delta_0=0.984, due to Ehrman); extend beyond-expansion exponents (any α>1/4\alpha>1/4) to an explicit threshold δ0≤0.9\delta_0\le 0.9, strictly below the best published value, ideally the full natural range δ>1/2\delta>1/2 of the counting framework. An answer must concern the same objects — the hypotenuse set SS of a thin Pythagorean orbit and its affine-sieve level of distribution — and must not replace them by a proxy (a different orbit or function, the area or product functions, or the lattice case δ=1\delta=1). Producing primes in SS, or an asymptotic local–global principle for thin hypotenuses (the statement that SS captures, up to the expected density, the integers that are locally represented, i.e. represented modulo every modulus), would subsume this problem but is not required for it.

References

  1. Jean Bourgain and Alex Kontorovich, The Affine Sieve Beyond Expansion I: Thin Hypotenuses, International Mathematics Research Notices 2015, no. 19 (2015), 9175–9205. DOI: 10.1093/imrn/rnu222, https://doi.org/10.1093/imrn/rnu222 (arXiv:1307.3535, https://arxiv.org/abs/1307.3535)
  2. Jean Bourgain, Alex Gamburd, and Peter Sarnak, Affine linear sieve, expanders, and sum-product, Inventiones Mathematicae 179(3) (2010), 559–644. DOI: 10.1007/s00222-009-0225-3, https://doi.org/10.1007/s00222-009-0225-3
  3. Alex Gamburd, On the spectral gap for infinite index 'congruence' subgroups of SL2(Z), Israel Journal of Mathematics 127 (2002), 157–200. DOI: 10.1007/BF02784530, https://doi.org/10.1007/BF02784530
  4. Alex V. Kontorovich, The hyperbolic lattice point count in infinite volume with applications to sieves, Duke Mathematical Journal 149(1) (2009), 1–36. DOI: 10.1215/00127094-2009-035, https://doi.org/10.1215/00127094-2009-035 (arXiv:0712.1391, https://arxiv.org/abs/0712.1391)
  5. Alex Kontorovich and Hee Oh, Almost prime Pythagorean triples in thin orbits, Journal für die reine und angewandte Mathematik 667 (2012), 89–131. DOI: 10.1515/crelle.2011.128, https://doi.org/10.1515/crelle.2011.128 (arXiv:1001.0370, https://arxiv.org/abs/1001.0370)
  6. Alex Kontorovich, Levels of Distribution and the Affine Sieve, Annales de la Faculté des sciences de Toulouse: Mathématiques 23(5) (2014), 933–966. DOI: 10.5802/afst.1432, https://doi.org/10.5802/afst.1432 (arXiv:1406.1375, https://arxiv.org/abs/1406.1375)
  7. Jiuzu Hong and Alex Kontorovich, Almost prime coordinates for anisotropic and thin Pythagorean orbits, Israel Journal of Mathematics 209(1) (2015), 397–420. DOI: 10.1007/s11856-015-1223-3, https://doi.org/10.1007/s11856-015-1223-3 (arXiv:1401.4701, https://arxiv.org/abs/1401.4701)
  8. Jean Bourgain and Alex Kontorovich, Beyond Expansion II: Low-Lying Fundamental Geodesics, Journal of the European Mathematical Society 19(5) (2017), 1331–1359. DOI: 10.4171/JEMS/694, https://doi.org/10.4171/JEMS/694 (arXiv:1406.1366, https://arxiv.org/abs/1406.1366)
  9. Max Ehrman, Almost Primes in Thin Orbits of Pythagorean Triangles, International Mathematics Research Notices 2019, no. 11 (2019), 3498–3526. DOI: 10.1093/imrn/rnx191, https://doi.org/10.1093/imrn/rnx191 (arXiv:1605.05265, https://arxiv.org/abs/1605.05265)
  10. Jean Bourgain and Alex Kontorovich, Beyond Expansion IV: Traces of Thin Semigroups, Discrete Analysis 2018, no. 6 (2018). DOI: 10.19086/da.3471, https://doi.org/10.19086/da.3471 (arXiv:1310.7190, https://arxiv.org/abs/1310.7190)
  11. Jean Bourgain and Alex Kontorovich, Beyond Expansion III: Reciprocal Geodesics, Duke Mathematical Journal 168(18) (2019), 3413–3435. DOI: 10.1215/00127094-2019-0056, https://doi.org/10.1215/00127094-2019-0056 (arXiv:1610.07260, https://arxiv.org/abs/1610.07260)
  12. Irving Calderón and Michael Magee, Explicit spectral gap for Schottky subgroups of SL(2,Z), Journal of the European Mathematical Society (2025). DOI: 10.4171/jems/1610, https://doi.org/10.4171/jems/1610 (arXiv:2303.17950, https://arxiv.org/abs/2303.17950; published online 13 March 2025)
  13. Louis Soares, Explicit spectral gap for Hecke congruence covers of arithmetic Schottky surfaces, Glasgow Mathematical Journal 68(2) (2026), 345–385. DOI: 10.1017/S0017089525100700, https://doi.org/10.1017/S0017089525100700 (arXiv:2305.02228, https://arxiv.org/abs/2305.02228; conditional on GRH for quadratic L-functions)
  14. Jiamin Li and Jianya Liu, Anisotropic quadratic equations in three variables, Pure and Applied Mathematics Quarterly (2025). DOI: 10.4310/pamq.250421122540, https://doi.org/10.4310/pamq.250421122540 (arXiv:2501.15033, https://arxiv.org/abs/2501.15033)
  15. Jiuzu Hong and Alex Kontorovich, Erratum to 'Almost prime coordinates for anisotropic and thin Pythagorean orbits', Israel Journal of Mathematics 266(1) (2025), 499–500. DOI: 10.1007/s11856-025-2742-1, https://doi.org/10.1007/s11856-025-2742-1