Explicit constructions of RIP matrices and related problems

Jean Bourgain, S. J. Dilworth, Kevin Ford, Sergei Konyagin, Denka Kutzarova

Introduction

While most authors work with real signals and matrices, in this paper we work with complex matrices for convenience. Given a complex matrix Φ\Phi satisfying (1.1), the 2n×2N2n\times 2N real matrix Φ′\Phi^{\prime}, formed by replacing each element a+iba+ib of Φ\Phi by the 2×22\times 2 matrix (ab−ba)(\begin{smallmatrix}a&b\\ -b&a\end{smallmatrix}), also satisfies (1.1) with the same parameters k,δk,\delta.

Given n,N,δn,N,\delta, we wish to find n×Nn\times N RIP matrices of order kk with constant δ\delta, and with kk as large as possible. If the entries of Φ\Phi are independent Bernoulli random variables with values ±1/n\pm 1/\sqrt{n}, then with high probability, Φ\Phi will have the required properties forFor convenience, we utilize the Vinogradov notation a≪ba\ll b, which means a=O(b)a=O(b), and the Hardy notation a≍ba\asymp b, which means b≪a≪bb\ll a\ll b.

See ; also for a proof based on the Johnson-Lindenstrauss lemma . The first result of similar type for these matrices is due to Kashin . See also for RIP matrices with rows randomly selected from the rows of a discrete Fourier transform matrix and for other random constructions of RIP matrices. The parameter kk cannot be taken larger; in fact

It is an open problem to find good explicit constructions of RIP matrices; see T. Tao’s Weblog for a discussion of the problem. We mention here that all known explicit examples of RIP matrices are based on constructions of systems of unit vectors (the columns of the matrix) with small coherence.

Matrices whose columns are unit vectors with small coherence are connected to a number of well-known problems, a few of which we describe below. Systems of vectors with small coherence are also known as spherical codes. Some other applications of matrices with small coherence may be found in .

Suppose that u1,…,uN{\mathbf{u}}_{1},\dots,{\mathbf{u}}_{N} are the columns of a matrix Φ\Phi and have coherence μ\mu. Then Φ\Phi satisfies RIP of order kk with constant δ=(k−1)μ\delta=(k-1)\mu.

For any kk-sparse vector x{\mathbf{x}},

All explicit constructions of matrices with small coherence are based on number theory. There are many constructions producing matrices with

In particular, such examples have been constructed by Kashin , Alon, Goldreich, Håstad and Peralta , DeVore , and Nelson and Temlyakov . By Proposition 1, these matrices satisfy RIP with constant δ\delta and order

It follows from random constructions of Erdős and Rényi for Turán’s problem (see Proposition 2 and (1.15) below) that for any n,Nn,N there are vectors with coherence

By contrast, there is a universal lower bound

valid for 2log⁡N≤n≤N/22\log N\leq n\leq N/2 and all Φ\Phi, due to Levenshtein (see also and ). Therefore, by estimating RIP parameters in terms of the coherence parameter we cannot construct n×Nn\times N RIP matrices of order larger than n\sqrt{n} and constant δ<1\delta<1.

Using methods of additive combinatorics, we construct RIP matrices of order kk with n=o(k2)n=o(k^{2}).

There is an effective constant ε0>0\varepsilon_{0}>0 and an explicit number n0n_{0} such that for any positive integers n≥n0n\geq n_{0} and n≤N≤n1+ε0n\leq N\leq n^{1+\varepsilon_{0}}, there is an explicit n×Nn\times N RIP matrix of order ⌊n12+ε0⌋\lfloor n^{\frac{1}{2}+\varepsilon_{0}}\rfloor with constant n−ε0n^{-\varepsilon_{0}}.

For application to sparse signal recovery, it is sufficient to take fixed δ<2−1\delta<\sqrt{2}-1 , and one needs an upper bound on nn in terms of k,Nk,N. By Theorem 1, for some ε0′>0\varepsilon_{0}^{\prime}>0, large NN and N1/2−ε0′≤k≤N1/2+ε0′N^{1/2-\varepsilon_{0}^{\prime}}\leq k\leq N^{1/2+\varepsilon_{0}^{\prime}}, we construct explicit RIP matrices with n≤k2−ε0′n\leq k^{2-\varepsilon_{0}^{\prime}}.

The proof of Theorem 1 uses a result on additive energy of sets (Corollary 2, Theorem 4), estimates for sizes of sumsets in product sets (Theorem 5), and bounds for exponential sums over products of sets possessing special additive structure (Lemma 10).

We now return to the problem of constructing matrices with small coherence. By (1.6), the bound (1.4) cannot be improved if log⁡n≫log⁡N\log n\gg\log N, but there is a gap between bounds (1.6) and (1.4) when log⁡n=o(log⁡N)\log n=o(\log N). For example, (1.4) is nontrivial only for n≫(log⁡N/log⁡log⁡N)2n\gg(\log N/\log\log N)^{2}. Of particular interest in coding theory is the range n=O(log⁡CN)n=O(\log^{C}N) for fixed CC, where there have been some improvements made to (1.4). A construction obtained by concatenating algebraic-geometric codes with Hadamard codes (see e.g. [23, Corollary 3] and Section 3 of ) produces matrices with coherence

which is nontrivial for n≫log⁡Nn\gg\log N, and is better than (1.4) when log⁡N≪n≪(log⁡Nlog⁡log⁡N)4\log N\ll n\ll(\frac{\log N}{\log\log N})^{4}. In the range (log⁡Nlog⁡log⁡N)5/2≪n≪(log⁡Nlog⁡log⁡N)5(\frac{\log N}{\log\log N})^{5/2}\ll n\ll(\frac{\log N}{\log\log N})^{5}, Ben-Aroya and Ta-Shma improved both (1.4) and (1.7) by constructing binary codes (vectors with entries ±1/n\pm 1/\sqrt{n}) with coherence

In this paper, we introduce very elementary constructions of matrices with coherence which matches (up to a log⁡log⁡N\log\log N factor) the bound (1.7). Our constructions, which are based on a method of Ajtai, Iwaniec, Komlós, Pintz and Szemerédi , have the added utility of applying to Turán’s power-sum problem and to the problem of finding thin sets with small Fourier coefficients. For the last two problems, our construction gives better estimates than existing explicit constructions in certain ranges of the parameters.

Roughly speaking, a set with small Fourier coefficients can be used to construct a set of numbers for Turán’s problem, and a set of numbers in Turán’s problem can be used to produce a matrix with small coherence. This is made precise below.

We next describe the problem of explicitly constructing thin sets with small Fourier coefficients. If NN is a positive integer and SS is a set (or multiset) of residues modulo NN, we let

Given NN, we wish to find a small set SS with ∣fS∣|f_{S}| also small.

Turán’s problem concerns the estimation of the function

where n,Nn,N are positive integers. There is a vast literature related to Turán’s problem; see, e.g., , , (chapter 5), , .

If S={t1,…,tn}S=\{t_{1},\dots,t_{n}\} is a multiset of integers modulo NN and zj=e2πitj/Nz_{j}=e^{2\pi it_{j}/N} for 1≤j≤n1\leq j\leq n, we see that

We also have the following easy connection between Turán’s problem and coherence.

Given any vector z=(z1,…,zn){\mathbf{z}}=(z_{1},\ldots,z_{n}) with ∣zj∣=1|z_{j}|=1 for all jj, the coherence μ\mu of the n×Nn\times N matrix with the columns

satisfies μ=n−1MN−1(z)\mu=n^{-1}M_{N-1}({\mathbf{z}}).

Combining (1.9) and Proposition 2, for any multiset SS of residues modulo NN, the vectors (1.10) satisfy

An application of Dirichlet’s approximation theorem shows that a set SS with ∣S∣<log⁡N|S|<\log N must have ∣fS∣≫1|f_{S}|\gg 1. In , sets which are not much larger are explicitly constructed so that ∣fS∣|f_{S}| is small. Specifically, by [1, (1),(2)], for each primeA corresponding result when NN is composite is given in . NN there is a set SS with ∣S∣=O(log⁡N(log⁡∗N)13log⁡∗N)|S|=O(\log N(\log^{*}N)^{13\log^{*}N}) and

where log⁡∗N\log^{*}N is the integer kk so that the kk-th iterate of the logarithm of NN lies in [1,e)[1,e). The proof uses an iterative procedure. By modifying this procedure, and truncating after two steps, we prove the following. To state our results, for brevity write

For sufficiently large prime NN and μ\mu such that

a set SS of residues modulo NN can be explicitly constructed so that

The method from , if applied without modification (with two iterations of the basic lemma), produces a conclusion in Theorem 2 with

The bound on ∣S∣|S| in Theorem 2 is better than (1.12) for μ≫L1−1/2L2\mu\gg L_{1}^{-1/2}L_{2}.

Together, the construction for Theorem 2 and (1.9) give explicit sets z{\mathbf{z}} for Turán’s problem. By further modifying the construction, we can do better.

For sufficiently large positive integer NN and μ\mu such that

a multiset z={z1,…,zn}{\mathbf{z}}=\{z_{1},\dots,z_{n}\} such that ∣z1∣=⋯=∣zn∣=1|z_{1}|=\dots=|z_{n}|=1, can be explicitly constructed so that

To put Theorem 3 in context, we briefly review what is known about T(n,N)T(n,N). P. Erdős and A. Rényi used probabilistic methods to prove an upper estimate

Using the character sum bound of Katz , J. Andersson gave explicit examples of sets z{\mathbf{z}} which give

One can see that (1.16) supersedes (1.15) for log⁡N≪log⁡2n\log N\ll\log^{2}n. Also, combining (1.16) with Proposition 2 provides yet another construction of matrices with coherence satisfying (1.4). On the other hand, by (1.6) and Proposition 2, we have the lower estimate

By comparison, the constructions in Theorem 3 are better than (1.16) in the range n≪L14/L28n\ll L_{1}^{4}/L_{2}^{8}, that is, throughout the range (1.14) (our constructions require nn to be prime, however).

The constructions in Theorem 3 also produce, by Proposition 2, explicit examples of matrices with coherence

which is close to the bound (1.7). By Proposition 1, these matrices satisfy RIP with constant δ\delta and order

We prove Theorem 1 in Sections 2–6, Theorem 2 in Section 7 and Theorem 3 in Section 8.

Construction of the matrix in Theorem 1

We fix a large even number mm. A value of mm can be specified; it depends on the constant c0c_{0} in an estimate from additive combinatorics (Proposition 3, Section 4). Also, the value mm can be reduced if one proves a better version of the Balog–Szemerédi–Gowers lemma (Lemma 6 below).

and the sets \curlyA,\curlyB\curly A,\curly B will be defined below. Notice that the matrix Φp\Phi_{p} can be extended to a n×Nn\times N matrix Φ\Phi by adding n−pn-p zero rows. Clearly, the matrices Φp\Phi_{p} and Φ\Phi have the same RIP parameters.

We notice that all elements of \curlyB\curly B are at most p/2p/2, and

For n≤N≤n1+β/2n\leq N\leq n^{1+\beta/2}, take Φ\Phi to be the matrix formed by the first NN columns of Φp\Phi_{p}, padded with n−pn-p rows of zeros.

In the next four sections, we show that Φ\Phi has the required properties for Theorem 1. First, in Section 3, we show that in (1.1) we need only consider vectors x{\mathbf{x}} whose components are 0 or 1 (emphflat vectors). We prove the following.

Let k≥210k\geq 2^{10} and ss be a positive integer. Assume that the coherence parameter of the matrix Φ\Phi is μ≤1/k\mu\leq 1/k. Also, assume that for some δ≥0\delta\geq 0 and any disjoint J1,J2⊂{1,…,N}J_{1},J_{2}\subset\{1,\dots,N\} with ∣J1∣≤k,∣J2∣≤k|J_{1}|\leq k,|J_{2}|\leq k we have

Then Φ\Phi satisfies the RIP of order 2sk2sk with constant 44sδlog⁡k44s\sqrt{\delta}\log k.

and, for a∈\curlyAa\in\curly A and a1,…,a2m∈\curlyA∖{a}a_{1},\dots,a_{2m}\in\curly A\setminus\{a\},

with some γ>0\gamma>0, where E(S,S)E(S,S) is the number of solutions of s1+s2=s3+s4s_{1}+s_{2}=s_{3}+s_{4} with each si∈Ss_{i}\in S.

holds where ε1=c0γ/20−43α/m\varepsilon_{1}=c_{0}\gamma/20-43\alpha/m.

The proof of Lemma 2 is quite involved, and will be handled in three subsequent sections. We next demonstate how Theorem 1 may be deduced from it.

We first prove (2.4) for the specific set \curlyA\curly A defined in (2.1), provided that p>(2m)8m2p>(2m)^{8m^{2}} (and thus L≥2mL\geq 2m). We have to show that for any distinct x,x1…,xn∈{1,…,L}x,x_{1}\dots,x_{n}\in\{1,\dots,L\} and any nonzero integers λ1,…,λn\lambda_{1},\dots,\lambda_{n} such that n≥2mn\geq 2m and ∣λ1∣+⋯+∣λn∣≤2m,|\lambda_{1}|+\cdots+|\lambda_{n}|\leq 2m, the sum

All summands in the right-hand side of (2.6) but the first one are divisible by x+x1+Ux+x_{1}+U. For the first summand we have

This shows that V1≠0  ( mod   x0+x1+U)V_{1}\neq 0\;(\bmod\;x_{0}+x_{1}+U). Therefore, V≠0V\neq 0. By assumption, p∤D1p\nmid D_{1}, and

Condition (2.5) is satisfied due to Corollary 4 of Section 5 with γ=β/50\gamma=\beta/50. If m>86000c0−1m>86000c_{0}^{-1} then Lemma 2 gives a nontrivial estimate with ε1>0\varepsilon_{1}>0. Thus, Φp\Phi_{p} satisfies the conditions of Corollary 1 with k=⌊p⌋≥n/2k=\lfloor\sqrt{p}\rfloor\geq\sqrt{n/2} and δ=p−ε1≤(n/2)−ε1\delta=p^{-\varepsilon_{1}}\leq(n/2)^{-\varepsilon_{1}} (using p≥0.9np\geq 0.9n for large nn, which follows from the prime number theorem). Let ε0=ε1/5\varepsilon_{0}=\varepsilon_{1}/5. Let n≤N≤n1+ε0n\leq N\leq n^{1+\varepsilon_{0}}, and let Φ\Phi be the n×Nn\times N matrix formed by taking the first NN columns of Φp\Phi_{p}, then adding n−pn-p rows of zeros. Clearly, Φ\Phi satisfies the conditions of Corollary 1 with the same parameters as Φp\Phi_{p}. By Lemma 1 with s=⌊pε1/4⌋s=\lfloor p^{\varepsilon_{1}/4}\rfloor, Theorem 1 follows.

In Section 4 we introduce some notation and recall standard estimates in additive combinatorics, which will be applied to subsets of \curlyB\curly B. Section 5 is devoted to the sumset theory of \curlyB\curly B, from which we deduce (2.5). The completion of the proof of Lemma 2 is in Section 6. We give some preliminaries here.

where the summands with a1=a2a_{1}=a_{2} are excluded from the summation. We next break Ω1,Ω2\Omega_{1},\Omega_{2} into balanced sets. For a∈\curlyAa\in\curly A and i=1,2i=1,2, let

whenever M1,M2M_{1},M_{2} are powers of two and, for i=1,2i=1,2 and for any ai∈Aia_{i}\in A_{i},

Indeed, there are O(log⁡2p)O(\log^{2}p) choices for M1,M2M_{1},M_{2}. To prove the cancellation in (2.8), we basically split into two cases: (i) some B′=Ωi(aj)B^{\prime}=\Omega_{i}(a_{j}) has additive structure (that is, E(B′,B′)E(B^{\prime},B^{\prime}) is large), where the cancellation comes from the sum over b1,b2b_{1},b_{2} (with a1,a2a_{1},a_{2} fixed), and (ii) when B′B^{\prime} does not have additive structure, in which case one gets dispersion of the phases from the dilation weights 1/(a1−a2)1/(a_{1}-a_{2}) (taking a large moment and using (2.4)). Incidentally, oscillations of the factor (a1−a2p)(\frac{a_{1}-a_{2}}{p}) play no role in the argument.

The Flat-RIP property

Let u1,…,uN{\mathbf{u}}_{1},\ldots,{\mathbf{u}}_{N} be the columns of an n×Nn\times N matrix Φ\Phi. Suppose that for every jj, ∥uj∥2=1\|{\mathbf{u}}_{j}\|_{2}=1. We say that Φ\Phi satisfies the flat RIP of order kk with constant δ\delta if for any disjoint J1,J2⊂{1,…,N}J_{1},J_{2}\subset\{1,\dots,N\} with ∣J1∣≤k,∣J2∣≤k|J_{1}|\leq k,|J_{2}|\leq k we have

For technical reasons, it is more convenient to work with the flat-RIP than with the RIP. However, flat-RIP implies RIP with an increase in δ\delta. The flat-RIP property is closely related to the property that (1.1) holds for any x{\mathbf{x}} with entries which are zero or one and at most kk ones (see the calculation at the end of this section).

Let k≥210k\geq 2^{10} and ss be a positive integer. Suppose that Φ\Phi satisfies flat-RIP of order kk with constant δ\delta. Then Φ\Phi satisfies RIP of order 2sk2sk with constant 44sδlog⁡k44s\delta\log k.

First, by a convexity-type argument and our assumption,

provided that ∣J1∣≤k,∣J2∣≤k|J_{1}|\leq k,|J_{2}|\leq k, 0≤xj,yj≤10\leq x_{j},y_{j}\leq 1 for all jj. Next, suppose ∣J1∣≤k,∣J2∣≤k|J_{1}|\leq k,|J_{2}|\leq k, and 0≤xj,yj0\leq x_{j},y_{j} for all jj. Without loss of generality assume that ∥x∥2=∥y∥2=1\|{\mathbf{x}}\|_{2}=\|{\mathbf{y}}\|_{2}=1, where ∥⋅∥2\|\cdot\|_{2} denotes the l2l_{2} norm. For a positive integer ν\nu let

Applying (3.2) to sets J1,ν,J2,νJ_{1,\nu},J_{2,\nu}, we get

Let t=⌊3+log⁡k/(2log⁡2)⌋t=\lfloor 3+\log k/(2\log 2)\rfloor. By the Cauchy–Schwarz inequality we infer that

For the next step, suppose xj,yjx_{j},y_{j} take arbitrary complex values, ∣J1∣≤sk|J_{1}|\leq sk and ∣J2∣≤sk|J_{2}|\leq sk. We partition J1J_{1} and J2J_{2} into ss subsets of cardinality at most kk each: J1=∪μ=1sJ1,μJ_{1}=\cup_{\mu=1}^{s}J_{1,\mu}, J2=∪μ=1sJ2,μ.J_{2}=\cup_{\mu=1}^{s}J_{2,\mu}. Next, for any jj we have

where xj,ν,yj,νx_{j,\nu},y_{j,\nu} are non-negative. By (3.4) and the Cauchy–Schwarz inequality,

For any disjoint J1,J2⊂{1,…,N}J_{1},J_{2}\subset\{1,\dots,N\} with ∣J1∣≤k,∣J2∣≤k|J_{1}|\leq k,|J_{2}|\leq k we have

Using the assumptions of the Lemma 1 directly rather than reducing it to Lemma 3, one can get a better constant for RIP; However, we do not need a stronger version of the corollary for our purposes.

Some definitions and results from additive combinatorics

For an (additive) abelian group GG we define the sum and the difference of subsets A,B⊂GA,B\subset G:

We will use the following lemma which is a particular case of Plünecke – Ruzsa estimates (, Exercise 6.5.15).

For any nonempty set A⊂GA\subset G we have ∣A+A∣≤∣A−A∣2/∣A∣|A+A|\leq|A-A|^{2}/|A|.

If A,B⊂GA,B\subset G, we define the (additive) energy E(A,B)E(A,B) of the sets AA and BB as the number of solutions of the equation

Next, let F⊂A×BF\subset A\times B. The FF-restricted sum of AA and BB is defined as

Trivially E(A,A)≤∣A∣3.E(A,A)\leq|A|^{3}. If E(A,A)E(A,A) is close to ∣A∣3|A|^{3} then AA must have a special additive structure.

(, Lemma 2.30) If E(A,A)≥∣A∣3/KE(A,A)\geq|A|^{3}/K then there exists F⊂A×AF\subset A\times A such that ∣F∣≥∣A∣2/(2K)|F|\geq|A|^{2}/(2K) and ∣A+FA∣≤2K∣A∣|A+_{F}A|\leq 2K|A|.

The following lemma is a version of the Balog–Szemerédi–Gowers lemma which plays a very important role in additive combinatorics.

If F⊂A×AF\subset A\times A, ∣F∣≥∣A∣2/L|F|\geq|A|^{2}/L and ∣A+FA∣≤L∣A∣|A+_{F}A|\leq L|A|. Then there exists a set A′⊂AA^{\prime}\subset A such that ∣A′∣≥∣A∣/(10L)|A^{\prime}|\geq|A|/(10L) and ∣A′−A′∣≤104L9∣A∣|A^{\prime}-A^{\prime}|\leq 10^{4}L^{9}|A|.

Combining Lemma 5 and Lemma 6 gives the following.

If E(A,A)≥∣A∣3/KE(A,A)\geq|A|^{3}/K then there exists a set A′⊂AA^{\prime}\subset A such that ∣A′∣≥∣A∣/(20K)|A^{\prime}|\geq|A|/(20K) and ∣A′−A′∣≤107K9∣A∣|A^{\prime}-A^{\prime}|\leq 10^{7}K^{9}|A|.

By 1A1_{A} we denote the indicator function of the set AA. With this notation, we have

An explicit version of Proposition 3, with c0=1/10430c_{0}=1/10430, is given in .

Note that if ∣A∣<∣B∣|A|<|B|, we may decompose BB as a disjoint union of at most 2∣B∣/∣A∣2|B|/|A| sets BjB_{j} with ∣A∣/2<∣Bj∣≤∣A∣|A|/2<|B_{j}|\leq|A| and apply (4.2) for each BjB_{j}. Hence

Applying the Cauchy–Schwarz inequality we get

It would be interesting to find best possible value for c0c_{0} in Proposition 3. The example A=B={1,…,[p]}A=B=\{1,\dots,[\sqrt{p}]\} shows that c0<1c_{0}<1.

Put λ(p)=0\lambda(p)=0, and let bb be a permutation of {1,…,p}\{1,\ldots,p\} such that λ(b1)≥⋯≥λ(bp)=0\lambda(b_{1})\geq\dots\geq\lambda(b_{p})=0. By (4.3), for 1≤j≤p−11\leq j\leq p-1 we have Sj≪GjS_{j}\ll G_{j}, where

Denote u0=∥λ∥2−2u_{0}=\|\lambda\|_{2}^{-2}. Notice that 1≤u0≤p1\leq u_{0}\leq p since ∥λ∥1=1\|\lambda\|_{1}=1. Separately considering j≤u0j\leq u_{0} and j>u0j>u_{0} and using the Cauchy–Schwarz inequality, we get

Using a parameter Δ≥1\Delta\geq 1 which will be specified later we define the sets

Decompose μ=μ−+μ0+μ+\mu=\mu_{-}+\mu_{0}+\mu_{+} where

The contribution to the sum in the theorem from μ−\mu_{-} and μ+\mu_{+} is negligible. First,

Using Young’s inequality (cf , Theorem 4.8), we find that

So, it suffices to estimate the contribution of μ0\mu_{0}. We have

Hence, ∣A∣≤∥μ∥2−2Δ2|A|\leq\|\mu\|_{2}^{-2}\Delta^{2}. Now we can use Corollary 2:

Combining the last inequality with (4.6) – (4.8) we get

Taking Δ=max⁡(1,S1/7)\Delta=\max(1,S^{1/7}) completes the proof of the theorem. ∎

A sumset estimate in product sets

The main result of this section is the following.

Then for any subsets A,B⊂\curlyCA,B\subset\curly C we have

Observe that for A=B=\curlyCA=B=\curly C we have ∣A+B∣=∣A∣τ′∣B∣τ′|A+B|=|A|^{\tau^{\prime}}|B|^{\tau^{\prime}} where

By Theorem 5, τ≤τ′\tau\leq\tau^{\prime}. On the other hand, τ>1/2\tau>1/2. If M→∞M\to\infty then

So, the asymptotic behavior of 2τM−12\tau_{M}-1 as M→∞M\to\infty is sharp. Likely, inequality (5.1) holds with τ=τ′\tau=\tau^{\prime}. This was proved in the case M=2M=2 by Woodall .

For positive integers K,LK,L we define an UR−UR-path as a sequence of pairs of integers \curlyP=((i1,j1)=(0,0),…,(iK+L−1,jK+L−1)=(K−1,L−1))\curly P=((i_{1},j_{1})=(0,0),\ldots,(i_{K+L-1},j_{K+L-1})=(K-1,L-1)) such that for any nn either in+1=in+1,jn+1=jni_{n+1}=i_{n}+1,j_{n+1}=j_{n}, or in+1=in,jn+1=jn+1i_{n+1}=i_{n},j_{n+1}=j_{n}+1.

Let KL≤M2KL\leq M^{2}, u0≥⋯≥uK−1≥0u_{0}\geq\dots\geq u_{K-1}\geq 0, v0≥⋯≥vL−1≥0v_{0}\geq\dots\geq v_{L-1}\geq 0, τ=τM\tau=\tau_{M}. Then there exists an UR−UR-path \curlyP\curly P such that

We proceed by induction on K+LK+L. For K=1K=1 or L=1L=1 the assertion is obvious. We prove it for K,LK,L with min⁡(K,L)≥2\min(K,L)\geq 2, KL≤M2KL\leq M^{2} supposing that it holds for (K,L)(K,L) replaced by (K−1,L)(K-1,L) and (K,L−1)(K,L-1). Without loss of generality we assume that

By the induction supposition, there exists an UR−UR-path \curlyP\curly P such that i1=1,j1=0i_{1}=1,j_{1}=0 and

Similarly, S≥(u0v0)τ+(1−v0)τ.S\geq(u_{0}v_{0})^{\tau}+(1-v_{0})^{\tau}. Thus, S≥w2τ+(1−w)τS\geq w^{2\tau}+(1-w)^{\tau} where

The function f(x)=x2τ+(1−x)τ−1f(x)=x^{2\tau}+(1-x)^{\tau}-1 has negative third derivative on $andandf(0)=f(1/M)=f(1)=0.ByRolle’stheorem,. By Rolle’s theorem,fhasnootherzerosonhas no other zeros on,andsince, and sincef(u)>0forforucloseto1,close to 1,f(x)\geq 0forfor1/M\leq x\leq 1.Therefore,. Therefore,f(w)\geq 0$ as desired. ∎

We will need Lemma 7 only for K=L=MK=L=M (although for the proof it was convenient to have varying K,LK,L).

Let U0,…,UM−1U_{0},\dots,U_{M-1}, V0,…,VM−1V_{0},\dots,V_{M-1} be non-negative numbers, and τ=τM\tau=\tau_{M}. Then

Lemma 8 has some similarity with inequality (2.1) from .

We order U0,…,UM−1U_{0},\dots,U_{M-1} and V0,…,VM−1V_{0},\dots,V_{M-1} in the descending order u0≥⋯≥uM−1u_{0}\geq\dots\geq u_{M-1} and v0≥⋯≥vM−1v_{0}\geq\dots\geq v_{M-1}, respectively, where for some permutations π\pi and σ\sigma of the set {0,…,M−1}\{0,\dots,M-1\} we have ui=Uπi,vj=Vσju_{i}=U_{\pi_{i}},v_{j}=V_{\sigma_{j}}. We consider an arbitrary UR−UR-path \curlyP\curly P with K=L=MK=L=M. Since ∣{πi1,…,πin}∣=in+1|\{\pi_{i_{1}},\dots,\pi_{i_{n}}\}|=i_{n}+1 and ∣{σj1,…,σjn}∣=jn+1|\{\sigma_{j_{1}},\dots,\sigma_{j_{n}}\}|=j_{n}+1,

Consequently, there is a permutation ψ\psi of {0,…,2M−2}\{0,\dots,2M-2\} so that

Thus, for some κ0∈{πi1,…,πin}\kappa_{0}\in\{\pi_{i_{1}},\dots,\pi_{i_{n}}\} and λ0∈{σj1,…,σjn}\lambda_{0}\in\{\sigma_{j_{1}},\dots,\sigma_{j_{n}}\} we have

But Uκ0=uiU_{\kappa_{0}}=u_{i} for some i∈{i1,…,in}i\in\{i_{1},\dots,i_{n}\}. Recalling that i1≤i2≤…i_{1}\leq i_{2}\leq\dots and u1≥u2≥…u_{1}\geq u_{2}\geq\dots we obtain Uκ0≥uinU_{\kappa_{0}}\geq u_{i_{n}}. Similarly, Vλ0≥vjnV_{\lambda_{0}}\geq v_{j_{n}}. Therefore,

Now we are ready to prove Theorem 5. We proceed by induction on rr. For r=0r=0 the set \curlyCM,r\curly C_{M,r} is a singleton, and there is nothing to prove. Now suppose that the assertion holds for rr replaced by r−1≥0r-1\geq 0. We consider arbitrary subsets A,B⊂\curlyC=\curlyCM,rA,B\subset\curly C=\curly C_{M,r}. For i=0,…,M−1i=0,\dots,M-1 we denote

Let D=A+BD=A+B. For n=0,…,2M−2n=0,\dots,2M-2 we denote

By the induction supposition, ∣Ai+Bj∣≥(∣Ai∣∣Bj∣)τ|A_{i}+B_{j}|\geq(|A_{i}||B_{j}|)^{\tau}. Hence,

The set −B-B is a translate of some set B′⊂\curlyBB^{\prime}\subset\curly B, and \curlyB\curly B is Freiman isomorphic to \curlyCM,r\curly C_{M,r}. Hence, for any B⊂\curlyBB\subset\curly B we have ∣B−B∣=∣B+B′∣≥∣B∣2τM|B-B|=|B+B^{\prime}|\geq|B|^{2\tau_{M}}. If ∣B∣>p1/4|B|>p^{1/4} then ∣B−B∣≥∣p∣(2τM−1)/4∣B∣|B-B|\geq|p|^{(2\tau_{M}-1)/4}|B|. By (5.2) and a short calculation using M≥215M\geq 2^{15}, p(2τM−1)/4≥pβ/5p^{(2\tau_{M}-1)/4}\geq p^{\beta/5}. ∎

Let E(S,S)=∣S∣3/KE(S,S)=|S|^{3}/K. By Corollary 1, there is a set B⊂SB\subset S such that ∣B∣≥∣S∣/(20K)|B|\geq|S|/(20K) and ∣B−B∣≤107K9∣S∣|B-B|\leq 10^{7}K^{9}|S|. If K≤pβ/50<p1/24K\leq p^{\beta/50}<p^{1/24} and pp is so large that 107≤pβ/5010^{7}\leq p^{\beta/50} then we get contradiction with Corollary 3. ∎

The proof of Lemma 2

We may assume ε1>0\varepsilon_{1}>0, otherwise there is nothing to prove. Adopt the notation (Ai,Mi,Ωi(a)A_{i},M_{i},\Omega_{i}(a)) from Section 2. If ∣A1∣M1<p1/2−γ/10|A_{1}|M_{1}<p^{1/2-\gamma/10}, then by (2.9), ∣S(A1,A2)∣≤2p1−γ/10|S(A_{1},A_{2})|\leq 2p^{1-\gamma/10} and (2.8) holds (recall that c0<1c_{0}<1, hence ε1<γ/20\varepsilon_{1}<\gamma/20). Thus, we can assume that ∣A1∣M1≥p1/2−γ/10|A_{1}|M_{1}\geq p^{1/2-\gamma/10}, which implies, by (2.3), that

Let WW denote the double sum over b1,b2b_{1},b_{2}. By the Cauchy–Schwarz inequality,

Another application of the Cauchy–Schwarz inequality gives

A third application of the Cauchy–Schwarz inequality, followed by Parseval’s identity yields a well-known inequality (cf. , Problem 14(a) for Chapter 6)

By (6.1), ∣Ωi(ai)∣≥p1/3|\Omega_{i}(a_{i})|\geq p^{1/3}, and by Lemma 9 and (2.5),

Thus, if ∣A1∣<pγ/2|A_{1}|<p^{\gamma/2} and ∣A2∣<pγ/2|A_{2}|<p^{\gamma/2}, then ∣S(A1,A2)∣≤4p1−γ/8|S(A_{1},A_{2})|\leq 4p^{1-\gamma/8} and (2.8) follows. Otherwise, without loss of generality we may assume that

The following lemma gives the necessary estimates to complete the proof of Lemma 2. For a1∈A1a_{1}\in A_{1}, set

The proof of Lemma 10 applies to more general sums, e.g. in T(A,B)T(A,B) one may replace the Legendre symbol (a1−a2p)(\frac{a_{1}-a_{2}}{p}) with arbitrary complex numbers ψ(a1,a2)\psi(a_{1},a_{2}) with modulus ≤1\leq 1, and one may replace 1a1−a2\frac{1}{a_{1}-a_{2}} with different quantities g(a1,a2)g(a_{1},a_{2}) having the dissociative property (the analog of (2.4) holds).

Postponing the proof of Lemma 10, we show first how to deduce Lemma 2.

We take a maximal subset B0⊂Ω1(a1)B_{0}\subset\Omega_{1}(a_{1}) so that (6.5) holds for B=B0B=B_{0}. Denote B1=Ω1(a1)∖B0B_{1}=\Omega_{1}(a_{1})\setminus B_{0}. By Lemma 9, (2.9), and (2.3) we have

Now assume that (6.6) does not hold. By (2.9), we get

Applying now Corollary 1 and (2.9) we obtain the existence of a set B1′⊂B1B_{1}^{\prime}\subset B_{1} such that

and ∣B1′−B1′∣≤107p27α∣B1∣≤p28α∣B1∣|B_{1}^{\prime}-B_{1}^{\prime}|\leq 10^{7}p^{27\alpha}|B_{1}|\leq p^{28\alpha}|B_{1}|. Using Lemma 10 we get inequality (6.5) for B=B1′B=B_{1}^{\prime}. Therefore, (6.5) is also satisfied for B=B0∪B1′B=B_{0}\cup B_{1}^{\prime}, contradicting the choice of B0B_{0}.

Thus, we have shown that (6.6) must hold. Using (6.5) for B=B0B=B_{0} and (6.7) we get

Summing on a1∈A1a_{1}\in A_{1} and using (2.3) and (2.9), we obtain

Hence, for some complex numbers εy,ξ\varepsilon_{y,\xi} of modulus ≤1\leq 1,

Then ∣ζ′(z)∣≤ζ(z)|\zeta^{\prime}(z)|\leq\zeta(z). By Hölder’s inequality,

As ζ(z)=∑ξ1B−B(z/ξ)\zeta(z)=\sum_{\xi}1_{B-B}(z/\xi), we have by the triangle inequality,

Define the probability measure λ1\lambda_{1} by

with a(1),…,a(2m)∈A2a^{(1)},\dots,a^{(2m)}\in A_{2}. By (2.4), this has only trivial solutions and thus

On the other hand, it follows from (6.3) and (6.4) that

Subsequent application of (6.9), (6.10) and (6.11) gives

By (6.3), p1/4≤∣B∣1/2p3αp^{1/4}\leq|B|^{1/2}p^{3\alpha}. Recalling γ≤α\gamma\leq\alpha, (2.9), (6.2) and (6.4), we conclude that

Plugging the last estimate into (6.8), we get

Thin sets with small Fourier coefficients

Denote by (a−1)m(a^{-1})_{m} the inverse of aa modulo mm. It is easy to see for relatively prime integers a,ba,b that

Let P≥4P\geq 4, S≥2S\geq 2, and RR be a positive integer. Suppose that for every prime p≤Pp\leq P, SpS_{p} is a set of integers in (−p/2,p/2)(-p/2,p/2). Suppose qq is a prime satisfying q≥RP2q\geq RP^{2}. Then the numbers r+s(p)(p−1)qr+s^{(p)}(p^{-1})_{q}, where 1≤r≤R,P/2<p≤P,s(p)∈Sp1\leq r\leq R,P/2<p\leq P,s^{(p)}\in S_{p}, are distinct modulo qq.

Multiplying both sides by p1p2p_{1}p_{2} gives

The right side is divisible by p1p2p_{1}p_{2} and the absolute value of the right side is <p1p2<p_{1}p_{2}, hence both sides are zero, r1=r2r_{1}=r_{2}, p1=p2p_{1}=p_{2} and s1(p1)=s2(p2)s_{1}^{(p_{1})}=s_{2}^{(p_{2})}. ∎

For brevity, we write e(z)e(z) for e2πize^{2\pi iz} is what follows.

Let P≥4P\geq 4, S≥2S\geq 2, and RR be a positive integer. Suppose that for every prime p∈(P/2,P]p\in(P/2,P], SpS_{p} is a multiset of integers in (−p/2,p/2)(-p/2,p/2), ∣Sp∣=S|S_{p}|=S and ∣fSp∣≤ε|f_{S_{p}}|\leq\varepsilon. Suppose qq is a prime satisfying q>Pq>P. Then the multiset

where VV is the number of primes in (P/2,P](P/2,P].

Since ∣fT(k)∣=∣fT(q−k)∣|f_{T}(k)|=|f_{T}(q-k)|, we may assume without loss of generality that 1≤k<q/21\leq k<q/2. We have

If k≥q/3k\geq q/3, we use the trivial bound ∣B(p,k)∣≤S|B(p,k)|\leq S and conclude

Now assume k≤q/3k\leq q/3. If p∣kp|k, then ∣B(p,k)∣≤S|B(p,k)|\leq S. When p∤kp\nmid k, by (7.1),

Since there are ≤log⁡klog⁡(P/2)\leq\frac{\log k}{\log(P/2)} primes p∣kp|k with p>P/2p>P/2, we have

Combining our estimates for ∣A(k)∣|A(k)| and ∣B(p,k)∣|B(p,k)|, we arrive at

For a specific choice of SpS_{p}, the inequality (7.2) can be strengthened.

Let P≥4P\geq 4 and RR be a positive integer. For every prime p∈(P/2,P]p\in(P/2,P] denote by SpS_{p} the set of all integers in (−p/2,p/2)(-p/2,p/2). Suppose qq is a prime satisfying q>Pq>P. Then the multiset

where VV is the number of primes in (P/2,P](P/2,P] and W=4log⁡(q/2)log⁡(P/2)W=4\frac{\log(q/2)}{\log(P/2)}.

Again, we may assume without loss of generality that 1≤k<q/21\leq k<q/2. We use notation from the proof of Lemma 12. If p∣kp|k, we use the trivial estimate ∣B(p,k)∣≤∣Sp∣≤P|B(p,k)|\leq|S_{p}|\leq P. Now there are ≤log⁡(q/2)log⁡(P/2)\leq\frac{\log(q/2)}{\log(P/2)} primes p∣kp|k with p>P/2p>P/2. When p∤kp\nmid k, by (7.1),

where it is assumed that k(q−1)p∈(−p/2,p/2)k(q^{-1})_{p}\in(-p/2,p/2). For a=1,…,[(P−1)/2]a=1,\dots,[(P-1)/2] we denote

Taking into account that ∣e(u)−1∣−1≤1/(4u)|e(u)-1|^{-1}\leq 1/(4u) for u∈(0,1/2]u\in(0,1/2] we get

If k(q−1)p=±ak(q^{-1})_{p}=\pm a then k±aqk\pm aq is divisible by pp. But ∣k±aq∣≤Pq/2|k\pm aq|\leq Pq/2. Therefore, the number of prime divisors p>P/2p>P/2 of any number k±aqk\pm aq is at most log⁡qlog⁡P/2+1\frac{\log q}{\log P/2}+1 and for any aa we get

Combining our estimates for ∣A(k)∣|A(k)| and ∣B(p,k)∣|B(p,k)| ((7.3) and (7.5)), we arrive at

Applying Lemma 12 for all primes qq in a dyadic interval, we can then feed these multisets T=TqT=T_{q} back into the lemma and iterate.

Using explicit estimates for counts of prime numbers , we have

For P≥250P\geq 250, there are more than 2P5log⁡(P/2)\frac{2P}{5\log(P/2)} primes in (P/2,P](P/2,P]. For any P>2P>2, there are at most 0.76P/log⁡P0.76P/\log P primes in (P/2,P](P/2,P].

Using Proposition 4 we obtain a more convenient version of Lemma 13.

Let P≥250P\geq 250. For every prime p∈(P/2,P]p\in(P/2,P] denote by SpS_{p} the set of all nonzero integers in (−p/2,p/2)(-p/2,p/2). Suppose qq is a prime satisfying q>Pq>P and suppose R≥1+log⁡(1+0.26P/log⁡(2q))/2R\geq 1+\log(1+0.26P/\log(2q))/2 is a positive integer. Then the multiset

We use the notation of Lemma 13. By Proposition 4 we have

On the other hand, using Proposition 4 again we get

Now the inequality (7.6) follows from (7.7) and (7.4). ∎

Using just one iteration one can get the following effective result on thin sets with small Fourier coefficients, of nearly the same strength as (1.12).

For sufficiently large prime NN and μ\mu such that N−1/2log⁡2N≤μ<1N^{-1/2}\log^{2}N\leq\mu<1 there is a set TT of residues modulo NN so that

Clearly, R≪1+log⁡(1/μ)R\ll 1+\log(1/\mu). Let TT be the multiset constructed in Lemma 14. We have ∣fT∣≤μ|f_{T}|\leq\mu. By Lemma 11, TT is a set. Moreover,

We choose real parameters P0P_{0}, P1P_{1} and positive integers R0R_{0}, R1R_{1} so that

For P0/2<p≤P0P_{0}/2<p\leq P_{0}, let SpS_{p} be the set of integers in (−p/2,p/2)(-p/2,p/2). By Lemmas 11, 14 and (7.8), for each prime q∈(P1/2,P1]q\in(P_{1}/2,P_{1}], there is a set T=SqT=S_{q} of residues modulo qq such that

By an application of Lemmas 11 and 12 with P=P1P=P_{1}, ε=ε1\varepsilon=\varepsilon_{1}, q=Nq=N, and S=R0∑P0/2<p≤P0pS=R_{0}\sum_{P_{0}/2<p\leq P_{0}}p, together with (7.9), there is a set TT of residues modulo NN so that

so that (7.9) follows immediately. The condition (1.13) implies (7.8) for large enough NN. ∎

Theorem 2 supersedes Corollary 5 for μ≫L1−1/2L21/2\mu\gg L_{1}^{-1/2}L_{2}^{1/2}.

An explicit construction for Turán’s problem

We follow the proof of Theorem 2 and Lemma 12. We choose real parameters P0P_{0}, P1P_{1} and a positive integer R0R_{0}, so that

For P0/2<p≤P0P_{0}/2<p\leq P_{0}, let SpS_{p} be the set of integers in (−p/2,p/2)(-p/2,p/2). By Lemma 14 and (8.1), for each prime q∈(P1/2,P1]q\in(P_{1}/2,P_{1}], there is a multiset T=SqT=S_{q} of residues modulo qq such that

We have ∣Sq∣=S|S_{q}|=S for all qq, where S=R0∑P0/2<p≤P0pS=R_{0}\sum_{P_{0}/2<p\leq P_{0}}p. Now define a multiset {z1,…,zn}\{z_{1},\dots,z_{n}\} as a union of multisets {e(s/q):s∈Sq,q∈(P1/2,P1]}\{e(s/q):s\in S_{q},q\in(P_{1}/2,P_{1}]\}. We have, for 1≤k≤N1\leq k\leq N,

If q∣kq|k, then B(q,k)=SB(q,k)=S. When q∤kq\nmid k, by (8.3), ∣B(q,k)∣≤ε1S|B(q,k)|\leq\varepsilon_{1}S. Therefore,

The sum over q∣kq|k is estimated at the same way as in Lemma 12:

Combining (8.4), (8.5) and using Proposition 4 we arrive at

as required. Moreover, by Proposition 4 we have

Now we take R0,P0,P1R_{0},P_{0},P_{1} the same as in the proof of Theorem 2 so that (8.2) follows immediately. The condition (1.14) implies (8.1) for large enough NN. ∎

As in , one can construct thin sets TT modulo NN with ∣T∣=o(L1L2)|T|=o(L_{1}L_{2}) and ∣fT∣|f_{T}| small, by iterating Lemma 12. Roughly speaking, applying Lemma 14 followed by rr iterations of Lemma 12 produces sets TT, with small ∣fT∣|f_{T}|, as small as ∣T∣=O(L1Lr+1)|T|=O(L_{1}L_{r+1}), where LjL_{j} is the jj-th iterate of the logarithm of NN. We omit the details.

Acknowledgments. The authors thank Ronald DeVore, Zeev Dvir, Venkatesan Guruswami, Piotr Indyk, Sina Jafarpour, Boris Kashin, Howard Karloff, Imre Leader, Igor Shparlinski and Avi Wigderson for helpful conversations.

References