The number of primes in short intervals and numerical calculations for Harman's sieve

Basic Information

• Paper ID: 2308.04458
• Title: The number of primes in short intervals and numerical calculations for Harman's sieve
• Author: Runbo Li
• Classification: math.NT (Number Theory)
• Publication Date: August 2023 (Latest version October 2025)
• Paper Link: https://arxiv.org/abs/2308.04458

Abstract

The author establishes nontrivial upper and lower bounds for the number of primes in the interval $[x - x^θ, x]$ where $0.52 \leqslant θ \leqslant 0.525$, and proves that for all sufficiently large $x$, the interval $[x - x^{0.52}, x]$ contains a prime. This improves upon the results of Baker, Harman, and Pintz (2001) and provides an affirmative answer to the Harman-Pintz conjecture. The paper effectively employs new arithmetic information, refined sieve decompositions, various techniques from Harman's sieve, and precise integral estimates.

Research Background and Motivation

Core Problem

This research addresses the classical problem of prime distribution in short intervals: determining the asymptotic behavior of the number of primes in intervals of the form $[x - x^θ, x]$, particularly finding the minimal exponent $θ$ such that such intervals necessarily contain a prime.

Significance of the Problem

1. Theoretical Significance: Prime distribution in short intervals is a central problem in analytic number theory, directly related to refinements of the Prime Number Theorem
2. Historical Importance: Since Cramér's 1937 conjecture, this problem has attracted the attention of numerous mathematicians
3. Technical Challenges: Requires combining sieve methods, exponential sum estimates, and deep results on L-function zero distribution

Historical Development

• 1930 Hoheisel: First proved the existence of $θ < 1$ for which asymptotic formulas hold ($θ ≥ 1 - 1/33000$)
• 1972 Huxley: Reduced $θ$ to $7/12$ through zero density estimates
• 2024 Guth-Maynard: Latest zero density results reduced $θ$ to $17/30$
• Sieve Method Development: Iwaniec-Jutila introduced sieve methods in 1979; subsequently BHP (2001) achieved $θ = 0.525$

Existing Limitations

Although the groundbreaking work of Baker-Harman-Pintz (BHP) achieved $θ = 0.525$, it faced the following issues:

1. Computational details were almost entirely omitted, making verification and improvement difficult
2. Kumchev pointed out that $θ = 0.52$ exceeded the scope of available techniques at that time
3. Lack of a systematic numerical computation framework

Core Contributions

1. Breaking the 0.52 barrier: First rigorously proved the feasibility of $θ = 0.52$, resolving the technical challenge raised by Kumchev
2. Providing complete computational details: Filled in the numerous omitted calculations in BHP's work, making results verifiable and reproducible
3. Developing new arithmetic information: Proved new arithmetic lemmas (Lemmas 4.4-4.5), extending the applicability of sieve methods
4. Optimizing sieve decomposition: Significantly improved loss estimates through role reversal, reverse Buchstab identities, and other techniques
5. Establishing a numerical computation framework: Used C++ and Mathematica for cross-verification to ensure computational accuracy

Detailed Methodology

Task Definition

Input: Positive integer $x$ and parameter $θ \in [0.52, 0.525]$Output: Upper and lower bounds for the number of primes $π(x) - π(x - x^θ)$ in the interval $[x - x^θ, x]$Objective: Prove the existence of constants $LB(θ), UB(θ)$ such that $LB(θ) \frac{x^θ}{\log x} ≤ π(x) - π(x - x^θ) ≤ UB(θ) \frac{x^θ}{\log x}$

Core Architecture: Harman's Sieve

1. Basic Setup

Define key parameters:

• $A = \{a : a ∈ \mathbb{Z}, x - y ≤ a < x\}$, where $y = x^{θ+ε}$
• $B = \{b : b ∈ \mathbb{Z}, x - y_1 ≤ b < x\}$, where $y_1 = x \exp(-3(\log x)^{1/3})$
• $S(C, z) = \sum_{a∈C, (a,P(z))=1} 1$

The objective is transformed into estimating $S(A, x^{1/2})$.

2. Sieve Decomposition Strategy

Recursive decomposition via Buchstab identity: $S(C, z) = S(C, w) - \sum_{w≤p<z} S(C_p, p)$

Key idea: Establish asymptotic relations between sparse set $A$ and dense set $B$: $S(A, z) = \frac{y}{y_1}(1 + o(1))S(B, z)$

3. Region Partition

Divide the integration domain into three subregions:

• Region A: $\{α_2 : 1/4 ≤ α_1 ≤ 2/5, 1/3(1-α_1) ≤ α_2 ≤ \min(α_1, 1/2(1-α_1))\}$
• Region B: $\{α_2 : 1/3 ≤ α_1 ≤ 1/2, \max(1/2 α_1, 1-2α_1) ≤ α_2 ≤ 1/2(1-α_1)\}$
• Region C: Remaining regions

Technical Innovations

1. New Arithmetic Information (Lemmas 4.4-4.5)

Extended traditional Type-II information to Type-II5 and Type-II6 through refined Hölder inequalities and mean value theorems:

Lemma 4.4: For five-dimensional cases, asymptotic formulas can be obtained under specific conditions: $\sum_{l_1l_2l_3l_4l_5∈A} a_{l_1}b_{l_2}c_{l_3}d_{l_4}e_{l_5}$

2. Role Reversal Technique

In four-dimensional and six-dimensional summations, exchange variable roles to replace larger 4D losses with smaller 6D losses.

3. Reverse Buchstab Identity

$\sum_{α_4} S(Ap_1p_2p_3p_4, (x/(p_1p_2p_3p_4))^{1/2}) = \sum_{α_4} S(Ap_1p_2p_3p_4, p_4) - \sum_{α_4} \sum_{α_4<α_5<1/2(1-α_1-α_2-α_3-α_4)} S(Ap_1p_2p_3p_4p_5, p_5)$

Through this "splitting" technique, almost primes become visible, thereby reducing losses.

Experimental Setup

Computing Environment

• Hardware: Intel(R) Xeon(R) Platinum 8383C CPU, 160 threads
• Software: C++ (primary computation) + Mathematica 14 (cross-verification)
• Parallelization: 80 Wolfram kernels

Numerical Integration Strategy

1. Hierarchical Computation: Calculate integrals in order of increasing dimensionality
2. Error Control: Use upper bound function $ω_1(u)$ and lower bound function $ω_0(u)$
3. Cross-Verification: Independent computation by C++ and Mathematica followed by comparison

Region Definition

Defined 13 different integration regions ($UC_{01}$ to $UC_{13}$, $VC_1$ to $VC_7$, etc.), each corresponding to specific sieve decomposition paths.

Experimental Results

Main Results

Theorem 1: For all sufficiently large $x$, the interval $[x - x^{0.52}, x]$ contains a prime.

Theorem 2: For $0.52 ≤ θ ≤ 0.525$ and $ε > 0$, $LB(θ) \frac{x^{θ+ε}}{\log x} ≤ π(x) - π(x - x^{θ+ε}) ≤ UB(θ) \frac{x^{θ+ε}}{\log x}$

Specific numerical results:

Technical Improvement Effects

Lower Bound Improvement

Through new decomposition strategies, total loss is controlled at:

• Region A: < 0.241 (2-dimensional loss)
• Region C: < 0.514 (including high-dimensional compensation)
• Total Loss: < 0.996, ensuring lower bound > 0.004

Upper Bound Calculation

Upper bound calculation is more complex, with total loss < 1.874, yielding upper bound < 2.874.

Numerical Verification

Comparison of C++ and Mathematica computational results demonstrates good consistency with errors within acceptable ranges.

Related Work

Historical Development

1. Classical Methods: Analytic methods based on Riemann zeta function zero distribution
2. Sieve Revolution: Iwaniec-Jutila (1979) introduced sieve methods
3. BHP Breakthrough: Combined Harman's sieve with Watt's power mean value theorem
4. Subsequent Developments: Various technical improvements and numerical optimizations

Position of This Work

• Deep optimization within the BHP framework
• First provision of complete verifiable computations
• Establishes foundation for further improvements

Applications and Extensions

The paper demonstrates 15 important applications:

1. Primes in Arithmetic Progressions (Theorem 3)

$π(x; q, a) - π(x - x^{0.52}; q, a) ≥ 0.004 \frac{x^{0.52}}{φ(q) \log x}$

2. Prime Gap Problem (Theorem 4)

Short intervals contain ≫ $x^{0.52}(\log x)^{-k}$ pairs of consecutive primes with bounded gaps.

3. Goldbach Problem

• Theorem 6: Almost all even numbers in intervals $[x, x + x^{13/225}]$ are Goldbach numbers
• Theorem 7: Intervals $[x, x + x^{26/1075}]$ contain Goldbach numbers

4. Other Applications

Including prime ideal distribution in imaginary quadratic fields, Carmichael number estimates, Linnik constants, etc.

Conclusions and Discussion

Main Conclusions

1. Technical Breakthrough: Successfully broke through the θ = 0.52 barrier
2. Methodological Innovation: Developed multiple sieve optimization techniques
3. Computational Verification: Established a complete numerical verification framework
4. Broad Applications: Achieved improvements on multiple number-theoretic problems

Limitations

1. Method Limitations: Still within the Harman sieve framework, difficult to achieve qualitative leaps
2. Computational Complexity: High-dimensional integral calculations are extremely complex, limiting further improvements
3. Constant Issues: While improving exponents, constants still have room for optimization

Future Directions

1. New Sieve Techniques: Explore new methods beyond existing frameworks
2. Computational Optimization: Develop more efficient numerical computation techniques
3. Theoretical Breakthroughs: Seek new connections with L-function zero distribution

In-Depth Evaluation

Strengths

1. Technical Depth: Achieved substantial breakthrough within an extremely complex technical framework
2. Computational Rigor: Provided unprecedented computational details and verification
3. Rich Applications: Demonstrated extensive number-theoretic applications
4. Methodological Innovation: Developed multiple valuable technical improvements

Weaknesses

1. Reading Difficulty: Technical details are extremely complex with high barriers to understanding
2. Limited Improvement: While achieving breakthrough, improvements are relatively modest
3. Limited Generalizability: Methods are highly specialized, difficult to extend to other problems

Impact

1. Academic Value: Resolved an important technical problem in the field
2. Methodological Contribution: Provided new technical tools for sieve theory
3. Practical Value: Generated direct applications to multiple related problems

Applicable Scenarios

This method is particularly suitable for number-theoretic problems requiring precise control of sieve losses, especially those involving prime distribution in short intervals.

References

The paper includes rich references spanning the complete development from classical Prime Number Theorem to latest sieve techniques, providing readers with guidance for in-depth study.

This paper represents an important advance in sieve theory within analytic number theory. Through refined technical improvements and extensive numerical computations, it successfully overcame long-standing technical barriers, laying a solid foundation for further developments in related fields.