alphaFeedback

WorksMathematics

On the van der Corput property for cubes

AI contributions
  • A large language model · LLM
View details

Abstract

We prove a power-saving form of the van der Corput property for cubes. As a consequence, we obtain a power-saving bound ∣A∣≪εN59/60+ε|A|\ll_\varepsilon N^{59/60+\varepsilon} for sets A⊆{1,…,N}A\subseteq\{1,\ldots,N\} with no nonzero cube differences.

Introduction

Answering a question of Lovász, Furstenberg [5] and Sárközy [20] independently proved that if A⊆[N]:={1,…,N}A \subseteq[N] := \{1,\ldots,N\} has no nonzero square difference, then ∣A∣=o(N)|A| = o(N). Sárközy’s proof gave

∣A∣≪N(log⁡log⁡N)2/3(log⁡N)−1/3.|A| \ll N(\log\log N)^{2/3}(\log N)^{-1/3}.

Pintz, Steiger, and Szemerédi [16 Theorem 1] used a more elaborate version of Sárközy’s density increment strategy to prove

∣A∣≪N(log⁡N)−112log⁡log⁡log⁡log⁡N.|A| \ll N(\log N)^{-\frac{1}{12}\log\log\log\log N}.

Bloom and Maynard [2] later improved this bound to

∣A∣≪N(log⁡N)−clog⁡log⁡log⁡N|A| \ll N(\log N)^{-c\log\log\log N}

for some absolute constant c>0c > 0. Recently, Green and Sawhney [7] combined Sárközy’s density increment strategy with their arithmetic level-dd inequality to show that, for some absolute c>0c > 0,

∣A∣≪Nexp⁡(−clog⁡N).|A| \ll N\exp(-c\sqrt{\log N}).

Fan and Lott [4] subsequently used the van der Corput property for squares to improve this to

∣A∣≪Nexp⁡(−clog⁡Nlog⁡log⁡N).(1)|A| \ll N\exp(-c\sqrt{\log N\log\log N}). \tag*{(1)}

for some absolute constant c>0c > 0.

A longstanding question is whether there is an absolute constant c>0c > 0 such that every A⊆[N]A \subseteq[N] with no nonzero square difference satisfies the power-saving bound ∣A∣≪N1−c|A| \ll N^{1-c}. Very recently, OpenAI [14] released a preprint claiming such a bound.

The lower bound of Fan and Lott [4 Theorem 1.3] shows that the van der Corput method cannot yield a power saving for sets with no nonzero square difference. For cubes, however, we obtain a power saving by exploiting the oddness of x↦x3x \mapsto x^3, which makes the complete cubic exponential sums real.

Write e(t)=exp⁡(2πit)e(t) = \exp(2\pi i t). Let F⊆ZF \subseteq\mathbb{Z}, let NN be a positive integer, and let δ>0\delta> 0. Consider a real trigonometric polynomial

P(θ)=a0+∑n=1Nancos⁡(2πnθ),a0,a1,…,aN∈R,(2)P(\theta) = a_0 + \sum_{n=1}^{N} a_n \cos(2\pi n\theta), \qquad a_0, a_1, \ldots, a_N \in\mathbb{R}, \tag*{(2)}

such that

  1. P(θ)≥0P(\theta) \ge0 for every θ∈R/Z\theta\in\mathbb{R}/\mathbb{Z},

  2. P(0)=1P(0) = 1,

  3. an=0a_n = 0 if 1≤n≤N1 \le n \le N and n∉F∪(−F)n \notin F \cup(-F),

  4. 0≤a0≤δ0 \le a_0 \le\delta.

We call a polynomial of the form (2) satisfying (i)–(iii) a van der Corput witness for FF of degree at most NN. Its nonconstant coefficients may have either sign. In the trigonometric characterization of Kamae and Mendès France [9] and Ruzsa [18], FF has the van der Corput property if for every δ>0\delta> 0 there exist an integer NN and a witness for FF satisfying (iv).

Suppose that A⊆[N]A \subseteq[N] and (A−A)∩F⊆{0}(A-A) \cap F \subseteq\{0\}. For any PP satisfying (i)–(iv), Fourier orthogonality gives

∣A∣≤(2N+1)a0≤(2N+1)δ.(3)|A| \le(2N+1)a_0 \le(2N+1)\delta. \tag*{(3)}

We give the argument in the proof of Corollary 1.2.

For the set of shifted primes F={p−1:p prime}F = \{p-1 : p\text{ prime}\}, Green [6, Theorems 1.1 and 1.2] constructed a witness of degree at most NN with a0≪N−ca_0 \ll N^{-c} for some absolute c>0c > 0. It follows that ∣A∣≪N1−c|A| \ll N^{1-c} whenever A⊆[N]A \subseteq[N] contains no two elements differing by one less than a prime.

Several other sets of forbidden differences admit power savings by the van der Corput method. Fan and Lott [3] proved that if A−AA-A contains no nonzero sum of two squares, then ∣A∣≪εN7/8+ε|A| \ll_{\varepsilon} N^{7/8+\varepsilon} for every ε>0\varepsilon> 0. Younis [22] proved that ∣A∣≪N1−c|A| \ll N^{1-c} for some absolute c>0c > 0 if A−AA-A contains no integer of the form n⌊21/3n⌋n\lfloor2^{1/3}n\rfloor with n∈N={1,2,…}n \in\mathbb{N} = \{1,2,\ldots\}. Kalogirou, Lott, Magyar, and Singha Roy [8] proved that for each odd k≥3k \ge3 and even s≥2s \ge2 there is a ck,s>0c_{k,s} > 0 such that ∣A∣≪k,sN1−ck,s|A| \ll_{k,s} N^{1-c_{k,s}} if A−AA-A contains no sum of ss positive kkth powers.

Matolcsi and Ruzsa [12, Theorem 2.5 and Corollary 2.6] proved that if B⊆Z/qZB \subseteq\mathbb{Z}/q\mathbb{Z} contains no two elements whose difference is a nonzero cube modulo qq, then ∣B∣≤q1−κ|B| \le q^{1-\kappa}, where κ=0.1195357…\kappa= 0.1195357\ldots. Our construction, like theirs, uses compatible witnesses for cubic residues modulo powers of each prime.

Ninčević and Slijepčević [13, Theorems 1 and 2] studied the van der Corput property for the values of odd polynomials h∈Z[x]h \in\mathbb{Z}[x] of degree k≥3k \ge3 with positive leading coefficient. They constructed witnesses of degree at most NN with nonnegative coefficients and a0≪h(log⁡N)−1/ka_0 \ll_h (\log N)^{-1/k}. For h(n)=nkh(n) = n^k they also proved that every such witness satisfies a0≫k1/log⁡Na_0 \gg_k 1/\log N.

By allowing the cosine coefficients to have either sign, we obtain a power saving for cubes.

Theorem 1.1. Let X≥1X \ge1 and let ε>0\varepsilon> 0. There exists a trigonometric polynomial with real coefficients of the form

PX(θ)=a0+∑1≤n≤Xancos⁡(2πn3θ)(4)P_X(\theta) = a_0 + \sum_{1 \le n \le X} a_n \cos(2\pi n^3\theta) \tag*{(4)}

such that

PX(θ)≥0(θ∈R/Z),PX(0)=1,0≤a0≪εX−1/20+ε.(5)P_X(\theta) \ge0 \quad(\theta\in\mathbb{R}/\mathbb{Z}), \qquad P_X(0) = 1, \qquad0 \le a_0 \ll_{\varepsilon} X^{-1/20+\varepsilon}. \tag*{(5)}

Taking X=⌊N1/3⌋X = \lfloor N^{1/3} \rfloor and applying Theorem 1.1 with 3ε3\varepsilon in place of ε\varepsilon gives a witness of degree at most NN with a0≪εN−1/60+εa_0 \ll_{\varepsilon} N^{-1/60+\varepsilon}. Applying (3) gives the following corollary.

Corollary 1.2. Let N≥1N \ge1 be an integer and let A⊆[N]A \subseteq[N] have no nonzero cube difference. Then for every ε>0\varepsilon> 0,

∣A∣≪εN59/60+ε.|A| \ll_{\varepsilon} N^{59/60+\varepsilon}.

Previously, the best upper bound was due to Adajar et al. [1 Theorem 1.1]. They combined the method of Green and Sawhney with exponential sum estimates in the style of Rice [17] to show that for every 0<μ<1/20 < \mu< 1/2 there is a cμ>0c_{\mu} > 0 such that

∣A∣≤Nexp⁡(−cμ(log⁡N)μ)|A| \le N \exp\left(-c_{\mu}(\log N)^{\mu}\right)

for all sufficiently large NN.

In the first version of their paper on squares [4 Section 1], Fan and Lott announced the van der Corput power saving for cube differences that we prove here. After they submitted that version to arXiv, OpenAI released a preprint claiming a power saving for intersective polynomial differences by a different method [15 Theorem 1.1]. In the cubic case, Corollary 1.2 gives ∣A∣≪cN1−c|A| \ll_c N^{1-c} for every 0<c<1/600 < c < 1/60, whereas the parameters in that preprint give c=c3≤[72006(1201!)48001]−1c = c_3 \le[72006(1201!)^{48001}]^{-1}.1

Let D3(N)D_3(N) denote the largest size of a subset of [N][N] with no nonzero cube difference. Using a construction of Ruzsa [19 Theorem 2] and a computer-assisted search, Lewko [11 Theorem 4] proved that

D3(N)≫Nβ,β=23+log⁡143log⁡91=0.861681….D_3(N) \gg N^{\beta}, \qquad\beta= \frac{2}{3} + \frac{\log14}{3\log91} = 0.861681\ldots.

There is still a considerable gap between this exponent β\beta and the exponent 59/6059/60 in Corollary 1.2. One may ask whether D3(N)∼CNγD_3(N) \sim CN^{\gamma} for some C>0C > 0 and γ∈[β,59/60]\gamma\in[\beta, 59/60]. A weaker conjecture is that the limit

lim⁡N→∞log⁡D3(N)log⁡N\lim_{N \to\infty} \frac{\log D_3(N)}{\log N}

exists. Ruzsa [19 Section 4] asked the corresponding question for squares.

Acknowledgments

Andrew Lott and Ákos Magyar thank the Erdős Center of the HUN-REN Alfréd Rényi Institute of Mathematics for its support and hospitality during the spring 2026 Simons Semester. Máté Matolcsi and Ákos Magyar were supported by grant NKFIH-154121. Imre Z. Ruzsa was supported by grant NKFIH-146387. This work was funded by a Simons Foundation International grant administered by the Simons Foundation [SFI-MPS-TSM-00026118, Á.M.].

AI disclosure

A large language model (LLM) was used as a sounding board during the developmental phase of the project. The LLM was also used for copyediting. However, the final argument is the authors’ own, and the manuscript was written and verified by the authors, who take full responsibility for the contents of the proof.

Sketch of the proof.

Our starting point is to define local weight functions Wp,KW_{p,K} on the groups Z/pKZ\mathbb{Z}/p^K\mathbb{Z}, where pp is prime and K∈NK \in\mathbb{N}, attached to a parameter 0<γp<10 < \gamma_p < 1. These functions are densities, relative to uniform measure, of geometric probability distributions defined by Pr⁡(pj∣n)=γpj\Pr(p^j \mid n) = \gamma_p^j. Extending them as pKp^K-periodic functions, we show that the corresponding exponential sums

Fp,m(a)=1pm∑n mod pmWp,K(n)e(an3/pm)(K=⌈m/3⌉)F_{p,m}(a) = \frac{1}{p^m}\sum_{n \bmod p^m} W_{p,K}(n)e(an^3/p^m) \qquad(K = \lceil m/3\rceil)

are nonnegative. With suitable optimized choices of γp\gamma_p, these sums recover the self-compatible modular witness functions constructed by Matolcsi–Ruzsa; see [12 Lemma 2.4 and Theorem 2.5] and the remark below. We will use simpler parameter choices. The local weights Wp,KW_{p,K} are constant on the “valuation shells” {n∈Z/pKZ:vp(n)=k}\{n \in\mathbb{Z}/p^K\mathbb{Z}: v_p(n) = k\}, where vp(n)=jv_p(n) = j if pj∥np^j \mathrel{\|} n. Thus they can be expressed as linear combinations of the Ramanujan sums cpk(n)c_{p^k}(n):

Wp,K(n)=∑k=0Kbpkcpk(n).W_{p,K}(n) = \sum_{k=0}^{K} b_{p^k}c_{p^k}(n).

with coefficients bpkb_{p^k} independent of the underlying group Z/pKZ\mathbb{Z}/p^K\mathbb{Z}. This allows one to define the coefficients bd=∏pk∥dbpkb_d = \prod_{p^k\mathrel{\|}d} b_{p^k} and the global truncated weight function

wR(n):=∑d≥1bd≥R−1bdcd(n).w_R(n) := \sum_{\substack{d\ge1\\ b_d\ge R^{-1}}} b_d c_d(n).

By exploiting the multiplicative properties of the coefficients bdb_d and the functions cd(n)c_d(n) we obtain the crucial estimate

SR(a,q):=EnwR(n)e(an3/q)≥−CϵR−1qϵ,(a,q)=1S_R(a,q) := \mathbb{E}_n w_R(n)e(an^3/q) \ge-C_{\epsilon}R^{-1}q^{\epsilon}, \qquad(a,q)=1

where by Enf(n)\mathbb{E}_n f(n) we denote the mean value of a periodic function over its period. The reason behind this estimate is that SR(a,q)S_R(a,q) can be approximated via the complete exponential sum

EnWD(n)e(an3/q)=∏pm∥qFp,m(ap)≥0(D=∏pm∥qp⌈m/3⌉).\mathbb{E}_n W_D(n)e(an^3/q) = \prod_{p^m\mathrel{\|}q} F_{p,m}(a_p) \ge0 \qquad\left(D = \prod_{p^m\mathrel{\|}q} p^{\lceil m/3\rceil}\right).

where ap≡a(q/pm)−1(modpm)a_p \equiv a(q/p^m)^{-1} \pmod{p^m} for pm∥qp^m \mathrel{\|} q, and WD(n)=∑d∣Dbdcd(n)=∏pK∥DWp,K(n)W_D(n) = \sum_{d\mid D} b_d c_d(n) = \prod_{p^K\mathrel{\|}D} W_{p,K}(n).

For X≥1X \ge1 and R≥1R \ge1, define the weighted exponential sum

TX(θ)=∑1≤n≤XhX(n)wR(n)e2πin3θ,hX(t)=6t2X3(1−t3X3)(0≤t≤X).T_X(\theta) = \sum_{1\le n\le X} h_X(n)w_R(n)e^{2\pi i n^3\theta}, \qquad h_X(t) = \frac{6t^2}{X^3}\left(1-\frac{t^3}{X^3}\right) \qquad(0\le t\le X).

We will take R=X1/20R = X^{1/20}. A major/minor arc argument then yields the uniform estimate

TX(0)=1+o(1),Re⁡TX(θ)≥−CεX−1/20+ε.T_X(0) = 1 + o(1), \qquad\operatorname{Re} T_X(\theta) \ge-C_{\varepsilon}X^{-1/20+\varepsilon}.

Taking η=CεX−1/20+ε\eta= C_{\varepsilon}X^{-1/20+\varepsilon} with a sufficiently large constant, we add η\eta and normalize:

PX(θ):=Re⁡TX(θ)+ηTX(0)+η.P_X(\theta) := \frac{\operatorname{Re} T_X(\theta) + \eta}{T_X(0) + \eta}.

Then PX≥0P_X \ge0, PX(0)=1P_X(0) = 1, and its constant Fourier coefficient is

a0=ηTX(0)+η≪εX−1/20+ε.a_0 = \frac{\eta}{T_X(0) + \eta} \ll_{\varepsilon} X^{-1/20+\varepsilon}.

providing the witness polynomial in Theorem 1.1. A short, well-known argument then gives Corollary 1.2.

The arithmetic factors

We define the finite exponential sums that appear in the major-arc approximation to TX(θ)T_X(\theta) and prove their required properties. We begin with nonnegative local modular witnesses attached to primes, following the modular construction of Matolcsi–Ruzsa; see [12 Section 2].

Local weight functions attached to a prime.

Fix a prime pp and a parameter 0<γ=γp<10 < \gamma= \gamma_p < 1, whose value will be determined later. For K≥1K \ge1, consider the probability distribution on Z/pKZ\mathbb{Z}/p^K\mathbb{Z} that is uniform on each valuation shell {vp(n)=j}\{v_p(n) = j\}, defined by

Pr⁡(pj∣n)=γj(0≤j≤K).(6)\Pr(p^j \mid n) = \gamma^j \qquad(0 \le j \le K). \tag*{(6)}

The uniform distribution corresponds to γ=p−1\gamma= p^{-1}. On Z/pKZ\mathbb{Z}/p^K\mathbb{Z}, we use the valuation: vp(n)=jv_p(n) = j when pj∥np^j \parallel n and j<Kj < K, while vp(0)=Kv_p(0) = K. Then

Pr⁡(vp(n)=j)=(1−γ)γj(0≤j<K),Pr⁡(vp(n)=K)=γK.(7)\Pr(v_p(n) = j) = (1-\gamma)\gamma^j \qquad(0 \le j < K), \qquad\Pr(v_p(n) = K) = \gamma^K. \tag*{(7)}

We define the local weight function Wp,KW_{p,K} on this group as the density of this probability distribution relative to uniform distribution, namely as

Wp,K(n)={p(1−γ)p−1(pγ)j,vp(n)=j<K,(pγ)K,pK∣n.(8)W_{p,K}(n) = \begin{cases} \frac{p(1-\gamma)}{p-1}(p\gamma)^j, & v_p(n) = j < K,\\ (p\gamma)^K, & p^K \mid n. \end{cases} \tag*{(8)}

We consider this function as a pKp^K-periodic function on the integers and set Wp,0=1W_{p,0} = 1. It is clear from the construction that it has mean value one over any of its periods, i.e. for every m≥Km \ge K, the numbers p−mWp,K(n)p^{-m}W_{p,K}(n), n mod pmn \bmod p^m also form a probability distribution. The probability distributions modulo pKp^K are compatible as KK varies, namely one has the compatibility property

1p∑t=0p−1Wp,K+1(n+tpK)=Wp,K(n)(n mod pK).\frac{1}{p}\sum_{t=0}^{p-1} W_{p,K+1}(n + tp^K) = W_{p,K}(n) \qquad(n \bmod p^K).

If vp(n)=j<Kv_p(n)=j<K, all lifts n+tpKn+tp^K of nn have valuation jj, and their weights equal Wp,K(n)W_{p,K}(n). For the zero residue modulo pKp^K, one lift has weight (pγ)K+1(p\gamma)^{K+1} and the other p−1p-1 lifts have weight p(1−γ)(pγ)K/(p−1)p(1-\gamma)(p\gamma)^K/(p-1). Hence their average is

1p((pγ)K+1+p(1−γ)(pγ)K)=(pγ)K=Wp,K(0).\frac{1}{p}\left((p\gamma)^{K+1}+p(1-\gamma)(p\gamma)^K\right)=(p\gamma)^K=W_{p,K}(0).

This compatibility underlies the construction of the global weight on Z\mathbb{Z}.

Positivity of local modular witness functions.

For p∤ap\nmid a and s≥1s\geq1, define

Ups(a)=1φ(ps)∑umodpsp∤ue(au3/ps).U_{p^s}(a)=\frac{1}{\varphi(p^s)}\sum_{\substack{u\bmod p^s\\p\nmid u}}e(au^3/p^s).

where φ\varphi is the Euler totient function. The change of variable u↦−uu\mapsto-u shows that these are real and lie in [−1,1][-1,1]. We include a proof of the following standard estimate for completeness.

Lemma 3.1. If p≠3p\neq3, then Ups(a)=0U_{p^s}(a)=0 for s≥2s\geq2. If p=3p=3, then U3s(a)=0U_{3^s}(a)=0 for s≥3s\geq3. Moreover, for p≠3p\neq3,

Up(a)≥−1+2pp−1.U_p(a)\geq-\frac{1+2\sqrt{p}}{p-1}.

Proof. For p≠3p\neq3 and s≥2s\geq2, first partition the units uu into residue classes modulo ps−1p^{s-1}. On the residue class represented by uu, replacing uu by u+tps−1u+tp^{s-1} multiplies the phase by e(3au2t/p)e(3au^2t/p). Summing over t mod pt\bmod p gives zero.

For p=3p=3 and s≥3s\geq3, the permutation u↦u+3s−2u\mapsto u+3^{s-2} of the units changes the phase by e(au2/3)=e(a/3)≠1e(au^2/3)=e(a/3)\neq1. The other terms in the binomial expansion are divisible by 3s3^s, so the sum again vanishes.

Finally, the complete sum ∑x mod pe(ax3/p)\sum_{x\bmod p}e(ax^3/p) has absolute value at most 2p2\sqrt{p}. Indeed, when p≡2(mod3)p\equiv2\pmod3, the map u↦u3u\mapsto u^3 is a permutation and the sum is zero. When p≡1(mod3)p\equiv1\pmod3, the cubic exponential sum is the sum of the two Gauss sums associated with the nontrivial cubic characters, each of absolute value p\sqrt{p}.

Subtracting the term x=0x=0 and dividing by p−1p-1 proves the asserted lower bound.

For m≥1m\geq1, put K=⌈m/3⌉K=\lceil m/3\rceil and define the local modular witness functions

Fp,m(a)=1pm∑n mod pmWp,K(n)e(an3/pm),Fp,0(a)=1.(9)F_{p,m}(a)=\frac{1}{p^m}\sum_{n\bmod p^m}W_{p,K}(n)e(an^3/p^m),\qquad F_{p,0}(a)=1. \tag*{(9)}

Lemma 3.2 (Local positivity). For any 0<γ<10<\gamma<1 and p∤ap\nmid a, we have

Fp,m+3(a)=γFp,m(a)(m≥0),(10)F_{p,m+3}(a)=\gamma F_{p,m}(a)\qquad(m\geq0), \tag*{(10)}

with initial values

Fp,1(a)=γ+(1−γ)Up(a),Fp,2(a)=γ+(1−γ)Up2(a).F_{p,1}(a)=\gamma+(1-\gamma)U_p(a),\qquad F_{p,2}(a)=\gamma+(1-\gamma)U_{p^2}(a).

In particular, with the choices

γ3=12,γp=1+2pp+2p(p≠3),(11)\gamma_{3}=\frac{1}{2},\qquad\gamma_{p}=\frac{1+2\sqrt{p}}{p+2\sqrt{p}}\quad(p\ne3), \tag*{(11)}

one has Fp,m(a)≥0F_{p,m}(a)\ge0 for every m≥0m\ge0 and every aa coprime to pp.

Proof. Note that if pK∣np^{K}\mid n then e(an3/pm)=1\mathrm{e}(an^{3}/p^{m})=1. If vp(n)=j<Kv_{p}(n)=j<K, write n=pjun=p^{j}u. The conditional average of the phase on this shell is Upm−3j(a)U_{p^{m-3j}}(a). Thus, summing over 0≤j<K0\le j<K gives

Fp,m(a)=γK+(1−γ)∑j=0K−1γjUpm−3j(a).(12)F_{p,m}(a)=\gamma^{K}+(1-\gamma)\sum_{j=0}^{K-1}\gamma^{j}U_{p^{m-3j}}(a). \tag*{(12)}

The first term is the contribution of the residues satisfying pK∣np^{K}\mid n.

For m≥1m\ge1, separating the j=0j=0 term for Fp,m+3F_{p,m+3} and shifting the remaining indices gives

Fp,m+3(a)=(1−γ)Upm+3(a)+γFp,m(a).F_{p,m+3}(a)=(1-\gamma)U_{p^{m+3}}(a)+\gamma F_{p,m}(a).

The term Upm+3(a)U_{p^{m+3}}(a) vanishes by Lemma 3.1. The lemma gives Fp,3(a)=γF_{p,3}(a)=\gamma, proving (10) also for m=0m=0. Formula (12) gives the stated initial values for m=1,2m=1,2.

For p≠3p\ne3, our choice in (11) satisfies

γp1−γp=1+2pp−1.\frac{\gamma_{p}}{1-\gamma_{p}}=\frac{1+2\sqrt{p}}{p-1}.

Hence Lemma 3.1 implies Fp,1(a)≥0F_{p,1}(a)\ge0, and it is easy to see that Fp,2(a)=γp>0F_{p,2}(a)=\gamma_{p}>0.

For p=3p=3, the bound U3,U9≥−1U_{3},U_{9}\ge-1 and the choice γ3=1/2\gamma_{3}=1/2 suffice for both initial values to be nonnegative. Finally, the recurrence (10) implies positivity at every higher modulus pmp^{m}.

Remark. The sums Fp,m(a)F_{p,m}(a) can be evaluated explicitly. If m=3j+rm=3j+r with r∈{0,1,2}r\in\{0,1,2\} and p∤ap\nmid a, then

Fp,m(a)={γpj,r=0,γpj(γp+(1−γp)Up(a)),r=1,γpj(γp+(1−γp)Up2(a)),r=2.(13)F_{p,m}(a)= \begin{cases} \gamma_{p}^{j}, & r=0,\\ \gamma_{p}^{j}\bigl(\gamma_{p}+(1-\gamma_{p})U_{p}(a)\bigr), & r=1,\\ \gamma_{p}^{j}\bigl(\gamma_{p}+(1-\gamma_{p})U_{p}^{2}(a)\bigr), & r=2. \end{cases} \tag*{(13)}

We have chosen the parameters γ=γp\gamma=\gamma_{p} above for their simplicity, and have not optimized them. If

tp=max⁡(0,max⁡p∤a(−Up(a)),max⁡p∤a(−Up2(a))),t_{p}=\max\left(0,\max_{p\nmid a}(-U_{p}(a)),\max_{p\nmid a}(-U_{p}^{2}(a))\right),

then Fp,m(a)≥0F_{p,m}(a)\ge0 for all m≥0m\ge0 and p∤ap\nmid a if and only if γ≥tp/(1+tp)\gamma\ge t_{p}/(1+t_{p}). Extend (9) to all numerators aa by the same formula. One may replace Wp,⌈m/3⌉W_{p,\lceil m/3\rceil} by Wp,mW_{p,m} in that formula: the weights agree below valuation ⌈m/3⌉\lceil m/3\rceil, and the remaining tail has the same mass and phase 1. Compatibility then gives Fp,m(pa)=Fp,m−1(a)F_{p,m}(pa)=F_{p,m-1}(a) for m≥1m\ge1. For p≠3p\ne3, choosing γ=tp/(1+tp)\gamma=t_{p}/(1+t_{p}) recovers the self-compatible family in [12 Lemma 2.4], with the base witness averaged under multiplication by nonzero cubes. For p=3p=3, the optimal choice in this family is γ=cos⁡(π/9)/(1+cos⁡(π/9))\gamma=\cos(\pi/9)/(1+\cos(\pi/9)); at modulus 27 it gives the invariant optimal base witness used in [12 Theorem 2.5], and the higher powers agree with its recursive extension. Thus the geometric probability laws give a realization of those modular witnesses. We keep the simpler choices in (11); optimizing them does not improve the exponent obtained from the estimates below. The properties of the parameters γp\gamma_{p} that we need are summarized in the following lemma.

Lemma 3.3. For all primes pp, we have

(i)1/p<γp≤p−1/3,and(ii)γp≪p−1/2.(14)\text{(i)}\quad1/p < \gamma_p \le p^{-1/3},\qquad\text{and}\qquad\text{(ii)}\quad\gamma_p \ll p^{-1/2}. \tag*{(14)}

Proof. The second assertion and the inequality 1/p<γp1/p < \gamma_p follow immediately from (11). The upper bound in (i) is also immediate for p=3p=3. For p≠3p\ne3, that upper bound is equivalent to

2p(p1/3−1)≤p−p1/3=(p1/3−1)(p2/3+p1/3),2\sqrt{p}\left(p^{1/3}-1\right)\le p-p^{1/3}=\left(p^{1/3}-1\right)\left(p^{2/3}+p^{1/3}\right),

which follows from the arithmetic–geometric mean inequality 2p≤p2/3+p1/32\sqrt{p}\le p^{2/3}+p^{1/3}.

Ramanujan sums and the global weight function wRw_R.

Recall the definition of the Ramanujan sums

cd(n)=∑rmodd(r,d)=1e(rn/d)=∑f∣dμ(d/f)f1f∣n,c1(n)=1.(15)c_d(n)=\sum_{\substack{r\bmod d\\(r,d)=1}}e(rn/d)=\sum_{f\mid d}\mu(d/f)f1_{f\mid n},\qquad c_1(n)=1. \tag*{(15)}

In particular for d=pkd=p^k

cpk(n)=pk1pk∣n−pk−11pk−1∣n(k≥1).c_{p^k}(n)=p^k1_{p^k\mid n}-p^{k-1}1_{p^{k-1}\mid n}\qquad(k\ge1).

We now express Wp,KW_{p,K} as a linear combination of cpkc_{p^k}, 0≤k≤K0\le k\le K. These K+1K+1 functions form an orthogonal basis for the space of functions constant on valuation shells, with respect to the inner product

⟨f,g⟩K=1pK∑n mod pKf(n)g(n)‾.\langle f,g\rangle_K=\frac{1}{p^K}\sum_{n\bmod p^K}f(n)\overline{g(n)}.

To be more precise, one has

⟨cpk,cpℓ⟩K={φ(pk),k=ℓ,0,k≠ℓ,0≤k,ℓ≤K.(16)\langle c_{p^k},c_{p^\ell}\rangle_K= \begin{cases} \varphi(p^k),& k=\ell,\\ 0,& k\ne\ell, \end{cases} \qquad0\le k,\ell\le K. \tag*{(16)}

Indeed, after expanding the two Ramanujan sums, a pair of characters has nonzero contribution only if their frequencies agree modulo one. Reduced fractions with denominators pkp^k and pℓp^\ell cannot agree mod 1 unless k=ℓk=\ell; in that case there are exactly φ(pk)\varphi(p^k) matching pairs.

Thus one may write

Wp,K(n)=∑k=0Kbpkcpk(n),bpk=⟨Wp,K,cpk⟩Kφ(pk).(17)W_{p,K}(n)=\sum_{k=0}^{K}b_{p^k}c_{p^k}(n),\qquad b_{p^k}=\frac{\langle W_{p,K},c_{p^k}\rangle_K}{\varphi(p^k)}. \tag*{(17)}

The constant coefficient is b1=1b_1=1, since Wp,KW_{p,K} has mean one. For 1≤k≤K1\le k\le K, direct calculation (or the probability interpretation) gives

⟨Wp,K,cpk⟩K=pkPr⁡(pk∣n)−pk−1Pr⁡(pk−1∣n)=pkγpk−pk−1γpk−1.\langle W_{p,K},c_{p^k}\rangle_K=p^k\Pr(p^k\mid n)-p^{k-1}\Pr(p^{k-1}\mid n) =p^k\gamma_p^k-p^{k-1}\gamma_p^{k-1}.

Dividing by φ(pk)=pk−1(p−1)\varphi(p^k)=p^{k-1}(p-1) yields

bpk=pγp−1p−1γpk−1(k≥1).(18)b_{p^k}=\frac{p\gamma_p-1}{p-1}\gamma_p^{k-1}\qquad(k\ge1). \tag*{(18)}

A crucial fact is that the coefficients bpkb_{p^k} do not depend on the truncation level K≥kK \ge k. This allows us to extend them multiplicatively, defining bd:=∏pk∥dbpkb_d := \prod_{p^k\parallel d} b_p^k with b1=1b_1=1. Note that for p≠3p \ne3, our choice of the γp\gamma_p parameters gives bp=2/(p+2)b_p = 2/(\sqrt{p}+2), while b3=1/4b_3=1/4. Since 0<bp≤γp0 < b_p \le\gamma_p, we obtain

0<bd≤∏pk∥dγpk≤d−1/3.(19)0 < b_d \le\prod_{p^k\parallel d}\gamma_p^k \le d^{-1/3}. \tag*{(19)}

For R≥1R \ge1, define the global weight function

wR(n):=∑d∈ERbdcd(n),ER:={d≥1:bd≥R−1}.(20)w_R(n) := \sum_{d\in\mathcal{E}_R} b_d c_d(n), \qquad\mathcal{E}_R := \{d \ge1 : b_d \ge R^{-1}\}. \tag*{(20)}

By (19), ER⊆[1,R3]\mathcal{E}_R \subseteq[1,R^3], so the defining sum is finite. Although every local weight Wp,KW_{p,K} is nonnegative, the truncated global weight wRw_R may have either sign.

Lower bounds for the global witness function at rational points a/qa/q.

Our aim is to show that the periodic mean of wR(n)e(an3/q)w_R(n)e(an^3/q) is bounded below by a small negative quantity. We write Enf(n)\mathbb{E}_n f(n) for the mean of a periodic function f:Z→Cf:\mathbb{Z}\to\mathbb{C} over any of its periods; the value is independent of the period chosen.

Proposition 3.4. Let (a,q)=1(a,q)=1. For every ε>0\varepsilon>0, uniformly for all reduced fractions a/qa/q and R≥1R\ge1, we have

SR(a,q):=EnwR(n)e(an3/q)≥−CεR−1qε.(21)S_R(a,q):=\mathbb{E}_n w_R(n)e(an^3/q)\ge-C_\varepsilon R^{-1}q^\varepsilon. \tag*{(21)}

Moreover SR(0,1)=1S_R(0,1)=1.

Let

Cd(a,q):=Encd(n)e(an3/q)(22)C_d(a,q):=\mathbb{E}_n c_d(n)e(an^3/q) \tag*{(22)}

denote the mean value of the periodic phase function cd(n)e(an3/q)c_d(n)e(an^3/q) over its period. By (20),

SR(a,q)=∑d∈ERbdCd(a,q).S_R(a,q)=\sum_{d\in\mathcal{E}_R}b_dC_d(a,q).

The means Cd(a,q)C_d(a,q) and SR(a,q)S_R(a,q) are real: replacing nn by −n-n conjugates the phase and leaves the weights unchanged. The key is to show that SR(a,q)S_R(a,q) is a truncated divisor sum.

Lemma 3.5. For (a,q)=1(a,q)=1, let

D(q):=∏pm∥qp⌈m/3⌉,D(1)=1.(23)D(q):=\prod_{p^m\parallel q}p^{\lceil m/3\rceil}, \qquad D(1)=1. \tag*{(23)}

Then

Cd(a,q)=0unless d∣D(q),C_d(a,q)=0 \quad\text{unless } d\mid D(q),

and ∣Cd(a,q)∣≤2ω(d)|C_d(a,q)|\le2^{\omega(d)}, where ω(d)\omega(d) is the number of distinct prime factors of dd.

Proof. A common period of the functions cd(n)c_d(n) and e(an3/q)e(an^3/q) is Δ=[d,q]\Delta=[d,q]. For every prime p∣Δp\mid\Delta, write kp=vp(d)k_p=v_p(d), mp=vp(q)m_p=v_p(q), and ℓp=max⁡(kp,mp)\ell_p=\max(k_p,m_p). If p∣qp\mid q, set

qp=q/pmp,ap≡aqp−1(modpmp).q_p=q/p^{m_p}, \qquad a_p\equiv aq_p^{-1}\pmod{p^{m_p}}.

The inverse qp−1q_p^{-1} (mod pmpp^{m_p}) exists, and apa_p is coprime to pp. The additive character and the Ramanujan sum factor as

e(an3/q)=∏p∣qe(apn3/pmp),cd(n)=∏p∣dcpkp(n).\mathrm{e}(an^3/q)=\prod_{p\mid q}\mathrm{e}(a_pn^3/p^{m_p}),\qquad c_d(n)=\prod_{p\mid d}c_{p^{k_p}}(n).

Indeed, the first identity follows from the fact that 1/q≡∑p∣qqp−1/pmp(mod1)1/q\equiv\sum_{p\mid q}q_p^{-1}/p^{m_p}\pmod{1}. The Chinese remainder theorem applied to Δ=∏p∣Δpβp\Delta=\prod_{p\mid\Delta}p^{\beta_p} yields the factorization

Cd(a,q)=∏p∣ΔLp(kp,mp;ap),(24)C_d(a,q)=\prod_{p\mid\Delta}\mathcal{L}_p(k_p,m_p;a_p), \tag*{(24)}

where, for m≥1m\ge1 and ℓ=max⁡(k,m)\ell=\max(k,m),

Lp(k,m;ap)=1pℓ∑x mod pℓcpk(x)e(apx3/pm).(25)\mathcal{L}_p(k,m;a_p)=\frac{1}{p^\ell}\sum_{x\bmod p^\ell}c_{p^k}(x)\mathrm{e}(a_px^3/p^m). \tag*{(25)}

with cp0=1c_{p^0}=1. For m=0m=0, the exponential is interpreted as 1 and no numerator is needed.

If m=0m=0 and k≥1k\ge1, the local factor is the mean of cpkc_{p^k}, hence zero. Suppose now that m≥1m\ge1 and k≥1k\ge1. For 0≤r≤k0\le r\le k, put

Ar=prpℓ∑xmodpℓpr∣xe(apx3/pm).A_r=\frac{p^r}{p^\ell}\sum_{\substack{x\bmod p^\ell\\p^r\mid x}}\mathrm{e}(a_px^3/p^m).

This is a normalized exponential sum, so ∣Ar∣≤1|A_r|\le1. The identity cpk(n)=pk1pk∣n−pk−11pk−1∣nc_{p^k}(n)=p^k1_{p^k\mid n}-p^{k-1}1_{p^{k-1}\mid n} implies Lp(k,m;ap)=Ak−Ak−1\mathcal{L}_p(k,m;a_p)=A_k-A_{k-1}. If k>⌈m/3⌉k>\lceil m/3\rceil, then 3(k−1)≥m3(k-1)\ge m, so the phase is one on both supports and Ak=Ak−1=1A_k=A_{k-1}=1. This proves the asserted support property. It also gives ∣Lp∣≤2|\mathcal{L}_p|\le2 when k≥1k\ge1; when k=0k=0, the local factor is an average of numbers of modulus one and has absolute value at most one. The stated bound follows from (24).

Proof of Proposition 3.4. Let D=D(q)D=D(q). Multiplicativity and (17) give

∑d∣Dbdcd(n)=∏pm∥qWp,⌈m/3⌉(n).\sum_{d\mid D}b_dc_d(n)=\prod_{p^m\parallel q}W_{p,\lceil m/3\rceil}(n).

Thus, by the Chinese remainder theorem,

∑d∣DbdCd(a,q)=1q∑n mod q(∏pm∥qWp,⌈m/3⌉(n))e(an3/q)=∏pm∥qFp,m(ap)≥0.(26)\begin{aligned} \sum_{d\mid D}b_dC_d(a,q)=\frac{1}{q}\sum_{n\bmod q}\left(\prod_{p^m\parallel q}W_{p,\lceil m/3\rceil}(n)\right)\mathrm{e}(an^3/q) \\ &=\prod_{p^m\parallel q}F_{p,m}(a_p)\ge0. \tag*{(26)} \end{aligned}

Indeed, each local weight function Wp,⌈m/3⌉W_{p,\lceil m/3\rceil} is pmp^m-periodic and Lemma 3.2 applies with a=apa=a_p. By Lemma 3.5,

SR(a,q)=∑d∣Dd∈ERbdCd(a,q).S_R(a,q)=\sum_{\substack{d\mid D\\d\in\mathcal{E}_R}}b_dC_d(a,q).

Every omitted coefficient satisfies bd<R−1b_d < R^{-1}, so the absolute value of the omitted contribution is at most

∑d∣Dd∉ERbd∣Cd(a,q)∣≤R−1∑d∣D2ω(d)=R−1∏pK∥D(1+2K)≤R−1τ(D)2≤R−1τ(q)2≪εR−1+ε.\sum_{\substack{d\mid D\\d\notin\mathcal E_R}} b_d\lvert C_d(a,q)\rvert\le R^{-1}\sum_{d\mid D}2^{\omega(d)} =R^{-1}\prod_{p^K\parallel D}(1+2K)\le R^{-1}\tau(D)^2\le R^{-1}\tau(q)^2\ll_{\varepsilon}R^{-1+\varepsilon}.

Together with the nonnegativity in (26), this proves (21). Finally, for q=1q=1 only d=1d=1 contributes, so SR(0,1)=1S_R(0,1)=1.

Lifting to integer cubes via the circle method.

Using a major/minor arc decomposition, we approximate Re⁡TX(θ)\operatorname{Re} T_X(\theta) on the major arcs by the product of the real arithmetic factor SR(a,q)S_R(a,q) and a nonnegative Archimedean term. The two elementary formulas for cdc_d in (15) serve different purposes: the divisor formula is used in the major-arc approximation, and the character formula allows us to apply Weyl’s inequality on the minor arcs.

A preliminary bound.

The quantity

AR=∑d∈ERbd∑f∣d∣μ(d/f)∣f(27)A_R=\sum_{d\in\mathcal E_R}b_d\sum_{f\mid d}\lvert\mu(d/f)\rvert f \tag*{(27)}

appears in both the major-arc approximation and the minor-arc estimate. We will use the following bound.

Lemma 4.1. For every ε>0\varepsilon>0 and R≥1R\ge1,

AR≪εR3+ε.(28)A_R\ll_{\varepsilon}R^{3+\varepsilon}. \tag*{(28)}

Proof. Note that on ER\mathcal E_R one has 1≤(Rbd)3+ε1\le(R b_d)^{3+\varepsilon}. Extending a sum of non-negative terms to all integers and using multiplicativity, we obtain

AR≤R3+ε∑d≥1bd4+ε∑f∣d∣μ(d/f)∣f=R3+ε∏p(1+∑k≥1bpk4+ε(pk+pk−1)).A_R\le R^{3+\varepsilon}\sum_{d\ge1}b_d^{4+\varepsilon}\sum_{f\mid d}\lvert\mu(d/f)\rvert f =R^{3+\varepsilon}\prod_p\left(1+\sum_{k\ge1}b_{p^k}^{4+\varepsilon}(p^k+p^{k-1})\right).

This Euler product converges. Indeed, put rp=pγp4+εr_p=p\gamma_p^{4+\varepsilon}. By Lemma 3.3,

rp≤p−(1+ε)/3<1,r_p\le p^{-(1+\varepsilon)/3}<1,

and, since bpk≤γpkb_{p^k}\le\gamma_p^k,

∑k≥1bpk4+ε(pk+pk−1)≤(1+p−1)∑k≥1rpk=(1+p−1)rp1−rp.\sum_{k\ge1}b_{p^k}^{4+\varepsilon}(p^k+p^{k-1})\le(1+p^{-1})\sum_{k\ge1}r_p^k=(1+p^{-1})\frac{r_p}{1-r_p}.

For sufficiently large pp, the estimate γp≪p−1/2\gamma_p\ll p^{-1/2} gives rp=Oε(p−1−ε/2)r_p=O_{\varepsilon}(p^{-1-\varepsilon/2}). Thus the local factor is 1+Oε(p−1−ε/2)1+O_{\varepsilon}(p^{-1-\varepsilon/2}), and the finitely many remaining factors are finite. This proves (28).

Approximation on the major arcs.

We derive the major-arc approximation by partial summation.

Lemma 4.2. Let (a,q)=1(a,q)=1 and θ=a/q+β\theta=a/q+\beta. Then

Re⁡TX(θ)=SR(a,q)J(βX3)+O(qARX(1+∣β∣X3)).\operatorname{Re} T_X(\theta)=S_R(a,q)J(\beta X^3)+O\left(\frac{qA_R}{X}(1+|\beta|X^3)\right).

where

J(t)=2∫01(1−u)cos⁡(2πtu) du=(sin⁡πtπt)2,J(0)=1.J(t)=2\int_0^1(1-u)\cos(2\pi tu)\,\mathrm{d}u=\left(\frac{\sin\pi t}{\pi t}\right)^2,\qquad J(0)=1.

In particular, 0≤J(t)≤10\leq J(t)\leq1.

Proof. We use the divisor formula

cd(n)=∑f∣dμ(d/f)f1f∣n.c_d(n)=\sum_{f\mid d}\mu(d/f)f\mathbf{1}_{f\mid n}.

For a fixed divisor ff, the function f1f∣ne(an3/q)f\mathbf{1}_{f\mid n}e(an^3/q) has period fqfq and its mean value over its period is

Af(a,q)=1q∑x mod qe(af3x3/q).A_f(a,q)=\frac{1}{q}\sum_{x\bmod q}e(af^3x^3/q).

Its partial sums, after subtracting the mean times the number of terms, are O(fq)O(fq). Indeed, the sum of the absolute values over a period is fqfq, and ∣Af(a,q)∣≤1|A_f(a,q)|\leq1. After subtracting the mean, complete periods cancel and the remaining terms contribute at most O(fq)O(fq).

The function H(t)=hX(t)e(βt3)H(t)=h_X(t)e(\beta t^3) vanishes at both endpoints of [0,X][0,X] and satisfies

∫0X∣H′(t)∣ dt≪X−1+∣β∣X2=X−1(1+∣β∣X3).\int_0^X|H'(t)|\,\mathrm{d}t\ll X^{-1}+|\beta|X^2=X^{-1}(1+|\beta|X^3).

Thus summation by parts gives

∑n≤XhX(n)f1f∣ne(θn3)=Af(a,q)∫0XhX(t)e(βt3) dt+O(fqX(1+∣β∣X3)).\sum_{n\leq X}h_X(n)f\mathbf{1}_{f\mid n}e(\theta n^3)=A_f(a,q)\int_0^Xh_X(t)e(\beta t^3)\,\mathrm{d}t+O\left(\frac{fq}{X}(1+|\beta|X^3)\right).

Multiplying by bdμ(d/f)b_d\mu(d/f) and summing over f∣df\mid d and d∈ERd\in\mathcal{E}_R, we obtain the arithmetic factor

∑d∈ERbd∑f∣dμ(d/f)Af(a,q)=∑d∈ERbdCd(a,q)=SR(a,q).\sum_{d\in\mathcal{E}_R}b_d\sum_{f\mid d}\mu(d/f)A_f(a,q)=\sum_{d\in\mathcal{E}_R}b_dC_d(a,q)=S_R(a,q).

The total error is bounded using (27). Finally, SR(a,q)S_R(a,q) is real. Taking real parts and substituting u=t3/X3u=t^3/X^3 gives

Re⁡∫0XhX(t)e(βt3) dt=2∫01(1−u)cos⁡(2πβX3u) du=J(βX3),\operatorname{Re}\int_0^Xh_X(t)e(\beta t^3)\,\mathrm{d}t=2\int_0^1(1-u)\cos(2\pi\beta X^3u)\,\mathrm{d}u=J(\beta X^3),

as required. □\square

Minor-arc estimates.

We use the standard cubic Weyl inequality. If (a,q)=1(a,q)=1 and ∣θ−a/q∣≤q−2|\theta-a/q|\le q^{-2}, then

∣∑n≤Ye(θn3+ξn)∣≪εY1+ε(q−1+Y−1+qY−3)1/4(Y≥1),(29)\left|\sum_{n\le Y}e(\theta n^3+\xi n)\right|\ll_{\varepsilon}Y^{1+\varepsilon}\left(q^{-1}+Y^{-1}+qY^{-3}\right)^{1/4}\qquad(Y\ge1), \tag*{(29)}

uniformly in the real linear coefficient ξ\xi, see [21 Lemma 2.4] or [10 equation (1)].

Lemma 4.3. Suppose 1≤Q≤X1\le Q\le X, (a,q)=1(a,q)=1, ∣θ−a/q∣≤q−2|\theta-a/q|\le q^{-2}, and Q<q≤X3/QQ<q\le X^3/Q. Then

∣TX(θ)∣≪εARXεQ−1/4.|T_X(\theta)|\ll_{\varepsilon}A_RX^{\varepsilon}Q^{-1/4}.

Proof. For every 1≤Y≤X1\le Y\le X, the right side of (29) is at most

CεX1+ε(q−1+X−1+qX−3)1/4≪εX1+εQ−1/4.C_{\varepsilon}X^{1+\varepsilon}\left(q^{-1}+X^{-1}+qX^{-3}\right)^{1/4}\ll_{\varepsilon}X^{1+\varepsilon}Q^{-1/4}.

Indeed, the expression Y4/q+Y3+qYY^4/q+Y^3+qY increases with YY and the second inequality follows from q>Qq>Q, X≥QX\ge Q, and q≤X3/Qq\le X^3/Q.

Summation by parts against hXh_X, whose total variation is O(X−1)O(X^{-1}), gives

∣∑n≤XhX(n)e(θn3+ξn)∣≪εXεQ−1/4,\left|\sum_{n\le X}h_X(n)e(\theta n^3+\xi n)\right|\ll_{\varepsilon}X^{\varepsilon}Q^{-1/4},

uniformly in ξ\xi. Using the character formula for the Ramanujan sum, we may write

TX(θ)=∑d∈ERbd∑rmodd(r,d)=1∑n≤XhX(n)e(θn3+rn/d).T_X(\theta)=\sum_{d\in\mathcal{E}_R}b_d\sum_{\substack{r\bmod d\\(r,d)=1}}\sum_{n\le X}h_X(n)e(\theta n^3+rn/d).

Since

φ(d)=∑f∣dμ(d/f)f≤∑f∣d∣μ(d/f)∣f,\varphi(d)=\sum_{f\mid d}\mu(d/f)f\le\sum_{f\mid d}|\mu(d/f)|f,

we have ∑d∈ERbdφ(d)≤AR\sum_{d\in\mathcal{E}_R}b_d\varphi(d)\le A_R. Therefore

∣TX(θ)∣≪εARXεQ−1/4.|T_X(\theta)|\ll_{\varepsilon}A_RX^{\varepsilon}Q^{-1/4}.

Proof of the main result.

We are now in a position to prove the main result. First we show the existence of the witness function PX(θ)P_X(\theta).

Proof of Theorem 1.1. It suffices to treat 0<ε<1/200<\varepsilon<1/20 and sufficiently large XX. Apply the preceding estimates with ε/2\varepsilon/2 in place of ε\varepsilon, and put R=X1/20R=X^{1/20}, Q=X4/5Q=X^{4/5}. Given θ∈T=R/Z\theta\in\mathbb{T}=\mathbb{R}/\mathbb{Z}, Dirichlet approximation with H=⌊X3/Q⌋H=\lfloor X^3/Q\rfloor gives a reduced fraction a/qa/q such that

1≤q≤H≤X3/Q,∣β∣:=∣θ−a/q∣≤1qH≤2QqX3.1\le q\le H\le X^3/Q,\qquad|\beta|:=|\theta-a/q|\le\frac{1}{qH}\le\frac{2Q}{qX^3}.

Also ∣β∣≤q−2|\beta| \le q^{-2} because q≤Hq \le H.

If q≤Qq \le Q, Proposition 3.4 and Lemma 4.2, together with 0≤J≤10 \le J \le1, imply

Re⁡TX(θ)≥−CεR−1Qε/2−O(ARQ/X).\operatorname{Re} T_X(\theta) \ge-C_{\varepsilon}R^{-1}Q^{\varepsilon/2} - O(A_RQ/X).

If q>Qq > Q, Lemma 4.3 yields ∣TX(θ)∣≪εARXε/2Q−1/4|T_X(\theta)| \ll_{\varepsilon} A_RX^{\varepsilon/2}Q^{-1/4}. Lemma 4.1 therefore reduces the three errors, up to factors bounded by XεX^{\varepsilon}, to

R−1=X−1/20,R3Q/X=X−1/20,R3Q−1/4=X−1/20.R^{-1} = X^{-1/20}, \qquad R^3Q/X = X^{-1/20}, \qquad R^3Q^{-1/4} = X^{-1/20}.

This implies the lower bound Re⁡TX(θ)≥−CεX−1/20+ε\operatorname{Re} T_X(\theta) \ge-C_{\varepsilon}X^{-1/20+\varepsilon} uniformly in θ\theta.

If θ=0\theta= 0, we have a=0a = 0, q=1q = 1, and β=0\beta= 0 in Lemma 4.2. Since SR(0,1)=J(0)=1S_R(0,1) = J(0) = 1, it follows

TX(0)=1+O(AR/X)=1+Oε(X−17/20+ε).T_X(0) = 1 + O(A_R/X) = 1 + O_{\varepsilon}(X^{-17/20+\varepsilon}).

Take η=CεX−1/20+ε\eta= C_{\varepsilon}X^{-1/20+\varepsilon}, with CεC_{\varepsilon} large enough that Re⁡TX(θ)+η≥0\operatorname{Re} T_X(\theta) + \eta\ge0 for all θ\theta. Since TX(0)=1+o(1)T_X(0) = 1 + o(1), the denominator TX(0)+ηT_X(0) + \eta is positive for sufficiently large XX. Define

PX(θ)=Re⁡TX(θ)+ηTX(0)+η.P_X(\theta) = \frac{\operatorname{Re} T_X(\theta) + \eta}{T_X(0) + \eta}.

Then PX≥0P_X \ge0, PX(0)=1P_X(0) = 1, and its constant Fourier coefficient is

a0=ηTX(0)+η≪εX−1/20+ε.a_0 = \frac{\eta}{T_X(0) + \eta} \ll_{\varepsilon} X^{-1/20+\varepsilon}.

All nonzero frequencies are ±n3\pm n^3 with 1≤n≤X1 \le n \le X, and the cosine coefficients are real because hX(n)wR(n)h_X(n)w_R(n) is real. For XX in a bounded range, use PX≡1P_X \equiv1 and enlarge the implied constant. This completes the proof for all X≥1X \ge1. □

Corollary 1.2 now follows from the following short, well-known argument.

Proof of Corollary 1.2. Put X=⌊N1/3⌋X = \lfloor N^{1/3}\rfloor and apply Theorem 1.1 with 3ε3\varepsilon in place of ε\varepsilon. Embed AA in the cyclic group G=Z/LZG = \mathbb{Z}/L\mathbb{Z}, where L=2N+1L = 2N + 1, and set 1A^(j)=∑a∈Ae(aj/L)\widehat{1_A}(j) = \sum_{a \in A} e(aj/L). By orthogonality we have

1L∑j=0L−1PX(j/L)∣1A^(j)∣2=a0∣A∣.(30)\frac{1}{L}\sum_{j=0}^{L-1} P_X(j/L)|\widehat{1_A}(j)|^2 = a_0|A|. \tag*{(30)}

Indeed, a nonconstant term could contribute only if a−b≡±n3(modL)a-b \equiv\pm n^3 \pmod L for some a,b∈Aa,b \in A. Since ∣a−b∣≤N−1|a-b| \le N-1 and n3≤Nn^3 \le N, the congruence would be an equality in Z\mathbb{Z}, contradicting our assumption on the set AA.

The terms on the left of (30) are nonnegative; hence, using the fact that PX(0)=1P_X(0) = 1, we have

∣A∣2L≤a0∣A∣.\frac{|A|^2}{L} \le a_0|A|.

For A≠∅A \ne\varnothing, this implies

∣A∣≤La0≪εNX−1/20+3ε≪εN1−1/60+ε.|A| \le La_0 \ll_{\varepsilon} NX^{-1/20+3\varepsilon} \ll_{\varepsilon} N^{1-1/60+\varepsilon}.
□\square

Steve Fan, Andrew Lott, and Ákos Magyar

Department of Mathematics, University of Georgia, Athens, GA 30602, USA

[email protected], [email protected], [email protected]

Máté Matolcsi

HUN-REN Alfréd Rényi Institute of Mathematics, Budapest, Hungary

Department of Analysis and Operations Research, Institute of Mathematics, Budapest University of Technology and Economics, Budapest, Hungary

[email protected]

Imre Z. Ruzsa

HUN-REN Alfréd Rényi Institute of Mathematics, Budapest, Hungary

[email protected]

References

  1. [1]C. F. E. Adajar, R. Agrawal, M. R. Choudhuri, C. Y. Chuah, S. Fan, S. Hegde, A. Lott, K. Nandakumar, and N. R. Pongandla, Extensions of the Furstenberg–Sárközy theorem via the arithmetic level-d inequality, arXiv:2605.16216v2 (2026).arxiv.org/abs/2605.16216
  2. [2]T. F. Bloom and J. Maynard, A new upper bound for sets with no square differences, Compos. Math. 158 (2022), 1777–1798.arxiv.org/abs/2011.13266
  3. [3]S. Fan and A. Lott, The van der Corput property for sums of two squares, arXiv:2606.29185v2 (2026).arxiv.org/abs/2606.29185
  4. [4]S. Fan and A. Lott, The optimal van der Corput property for squares, arXiv:2610.07765v2 (2026).arxiv.org/abs/2610.07765
  5. [5]H. Furstenberg, Ergodic behavior of diagonal measures and a theorem of Szemerédi on arithmetic progressions, J. Analyse Math. 31 (1977), 204–256.DOI
  6. [7]B. Green and M. Sawhney, New bounds for the Furstenberg–Sárközy theorem, arXiv:2411.17448v2 (2025).arxiv.org/abs/2411.17448
  7. [8]A. Kalogirou, A. Lott, Á. Magyar, and A. Singha Roy, The Furstenberg–Sárközy theorem for sums of an even number of odd powers, arXiv:2609.16595 (2026).arxiv.org/abs/2609.16595
  8. [9]T. Kamae and M. Mendès France, Van der Corput’s difference theorem, Israel J. Math. 31 (1978), 335–342.DOI
  9. [10]B. Kerr, On the cubic Weyl sum, arXiv:2101.08058 (2021).arxiv.org/abs/2101.08058
  10. [11]M. Lewko, An improved lower bound related to the Furstenberg–Sárközy theorem, Electron. J. Combin. 22 (2015), Paper P1.32.DOI
  11. [12]M. Matolcsi and I. Z. Ruzsa, Difference sets and positive exponential sums. II: Cubic residues in cyclic groups, Proc. Steklov Inst. Math. 314 (2021), 138–143.DOI
  12. [14]OpenAI, A power saving for square-difference-free sets, preprint (2026), with accompanying Lean formalization.hub.valency.io/works/W_a9yjw7uw
  13. [15]OpenAI, A power saving for inters­ective polynomial differences with an exponent depending only on the degree, preprint (2026).hub.valency.io/works/W_86qmxgsn
  14. [16]J. Pintz, W. L. Steiger, and E. Szemerédi, On sets of natural numbers whose difference set contains no squares, J. London Math. Soc. (2) 37 (1988), 219–231.DOI
  15. [17]A. Rice, A maximal extension of the best-known bounds for the Furstenberg–Sárközy theorem, Acta Arith. 187 (2019), 1–41.arxiv.org/abs/1612.01760
  16. [18]I. Z. Ruzsa, Connections between the uniform distribution of a sequence and its differences, in Topics in Classical Number Theory, vol. II, Colloq. Math. Soc. János Bolyai, vol. 34, North-Holland, Amsterdam, 1984, pp. 1419–1443.
  17. [19]I. Z. Ruzsa, Difference sets without squares, Period. Math. Hungar. 15 (1984), 205–209.DOI
  18. [20]A. Sárközy, On difference sets of sequences of integers. I, Acta Math. Acad. Sci. Hungar. 31 (1978), 125–149.DOI
  19. [21]R. C. Vaughan, The Hardy–Littlewood Method, 2nd ed., Cambridge Tracts in Mathematics, vol. 125, Cambridge University Press, Cambridge, 1997.DOI
  20. [22]K. Younis, Sets whose differences avoid a bracket quadratic, arXiv:2608.30078 (2026).arxiv.org/abs/2608.30078

Paper details

vlnc:W_dh87yzgy
v1 · latestHistory (1)

published · Fresh — added to the Hub within the last 30 days

added ·

Usage

views · 7

downloads · 0

A view counts once per browser per day, a download each time the PDF is requested.

CC BY 4.0

Original source

  • Mathematics

A large language model (LLM) · View statement →

Looking for related papers…
Contents