arxiv: v2 [math.co] 2 Dec 2008

Transcription

1 FOURIER ANALYSIS, LINEAR PROGRAMMING, AND DENSITIES OF DISTANCE AVOIDING SETS IN R n arxiv: v2 [math.co] 2 Dec 2008 FERNANDO MÁRIO DE OLIVEIRA FILHO AND FRANK VALLENTIN Abstract. In this paper we derive new upper bounds for the densities of measurable sets in R n which avoid a finite set of prescribed distances. The new bounds come from the solution of a linear programming problem. We apply this method to obtain new upper bounds for measurable sets which avoid the unit distance in dimensions 2,...,24. This gives new lower bounds for the measurable chromatic number in dimensions 3,..., 24. We apply it to get a short proof of a variant of a recent result of Bukh which in turn generalizes theorems of Furstenberg, Katznelson, and Weiss and Bourgain and Falconer about sets avoiding many distances.. Introduction Let d,..., d N be positive real numbers. We say that a subset A of the n- dimensional Euclidean space R n avoids the distances d,..., d N if the distance between any two points in A is never d,..., d N. We define the upper density of a Lebesgue measurable set A R n as δ(a) = limsup T vol(a [ T,T] n ) vol[ T,T] n. In this expression [ T,T] n denotes the regular cube in R n with side 2T centered at the origin. We denote the extreme density which a measurable set in R n that avoids the distances d,..., d N can have by m d,...,d N (R n ) = sup{δ(a) : A R n is measurable and avoids distances d,..., d N }. In this paper we derive upper bounds for this extreme density from the solution of a linear programming problem. To formulate our main theorem we consider the function Ω n given by ( n )( 2 2 () Ω n (t) = Γ 2 t) (n 2) J (n 2)(t), for t > 0, Ω n(0) =, 2 where J 2 (n 2) is the Bessel function of the first kind with parameter (n 2)/2. To fix ideas we plotted the graph of the function Ω 4 in Figure.. Date: November 9, Mathematics Subject Classification. 42B05, 52C0, 52C7, 90C05. Key words and phrases. Nelson-Hadwiger problem, measurable chromatic number, linear programming, almost periodic functions, autocorrelation function. The first author was partially supported by CAPES/Brazil under grant BEX 242/04-6. The second author was partially supported by the Deutsche Forschungsgemeinschaft (DFG) under grant SCHU 503/4. The second author thanks the Hausdorff Research Institute for Mathematics (Bonn) for its hospitality and support.

2 2 F.M. de Oliveira Filho, F. Vallentin Ω 4 (t) t Figure.. Graph of the function Ω 4 (t) = 2 t J (t) Theorem.. Let d,..., d N be positive real numbers. Let A R n be a measurable set which avoids the distances d,..., d N. Suppose there are real numbers z 0, z,..., z N which sum up to at least one and which satisfy z 0 +z Ω n (td )+z 2 Ω n (td 2 )+ +z N Ω n (td N ) 0 for all t 0. Then, the upper density of A is at most z 0. In Section 2 we provide a proof where we make essential use of basic harmonic analysis, which we briefly recall. In the sections that follow we apply the main theorem in a variety of situations: sets avoiding one distance, sets avoiding two distances, and sets avoiding many distances. For the history of these Euclidean distance problems we refer to the surveys by Székely [9] and Székely and Wormald [8] and the references therein. Sets avoiding one distance have been studied by combinatorialists because of their relation to the measurable chromatic number of the Euclidean space. This is the minimum number of colors one needs to color all points in R n so that any two points at distance receive different colors and so that points receiving the same color form Lebesgue measurable sets; it will be denoted by χ m (R n ). Since every color class of a coloring provides a measurable set which avoids the distance, we have the inequality (2) m (R n ) χ m (R n ). For the plane it is only known that 5 χ m (R 2 ) 7, where the lower bound is due to Falconer [0] and the upper bound comes e.g. from a coloring one constructs using a tiling by regular hexagons with circumradius slightly less than. Erdős conjectured that m (R 2 ) < /4 so that (2) would yield an alternative proof of Falconer s result. So far the best known results on m (R 2 ) are the lower bound m (R 2 ) by Croft [7] and the upper bound m (R 2 ) 2/ by

3 Fourier analysis, linear programming, and densities of distance avoiding sets 3 Székely[8]. In Section 3 we compute new upper bounds for m (R n ) for dimensions n = 2,..., 24 based on a strengthening of our main theorem by extra inequalities. These new upper bounds for m (R n ) imply by (2) new lower bounds for χ m (R n ) in dimensions 3,..., 24. If one considers sets which avoid more than one distance one can ask how N distances can be chosen so that the extreme density becomes as small as possible: What is the value of inf{m d,...,d N (R n ) : d,..., d N > 0} for fixed N? For planar sets avoidingtwo distances Székely[8] showed that inf{m d,d 2 (R 2 ) : d,d 2 > 0} m, 3 (R 2 ) 2/ In Section 4 we improve his result and show that inf{m d,d 2 (R 2 ) : d,d 2 > 0} Recently, Bukh [5], using harmonic analysis and ideas resembling Szémeredi s regularity lemma, showed that inf{m d,...,d N (R n ) : d,..., d N > 0} drops to zero exponentially in N: He shows that there is a number r, strictly greater than, which depends only on N and n so that if d 2 /d > r, d 3 /d 2 > r,..., d N /d N > r, then m d,...,d N (R n ) (m (R n )) N. This implies a theorem of Furstenberg, Katznelson, and Weiss [3] that for every subset A in the plane which has positive upper density there is a constant d so that A does not avoid distances larger than d. Their original proof used tools from ergodic theory and measure theory. Alternative proofs have been proposed by Bourgain [4] using elementary harmonic analysis and by Falconer and Mastrand [2] using geometric measure theory. Bukh s result also implies that m d,...,d N (R n ) becomes arbitrarily small if the distances d, d 2,..., d N become arbitrarily small. This is originally due to Bourgain [4] and Falconer []. In Section 5 we give a short proof of a variant of Bukh s result using our main theorem, where we replace (m (R n )) N by the weaker estimate 2 N. We could improve this considerably, but we cannot get (m (R n )) N. Still our estimate is strong enough to give all the implications mentioned. Furthermore, our proof has the additional advantage that it easily provides quantitative estimates about the spacing r between the distances. The idea of linear programming bounds for packing problems of discrete point sets in compact metric spaces goes back to Delsarte [9] and it has been successfully applied to a variety of situations. Cohn and Elkies [6] were the first who were able to set up a linear programming bound for packing problems in non-compact spaces; by then no less than 30 years since Delsarte s fundamental contribution had gone by. Our main theorem can be viewed as a continuous analogue to their linear programming bound. 2. Proof of the main theorem For the proof of our main theorem elementary notions from harmonic analysis will be important. We recall these here. For details we refer to, e.g., the book by Katznelson [4]. A measurable, complex valued function f: R n C is called periodic if it is invariantunderann-dimensionaldiscretesubgroupofr n or,inotherwords,ifthere is a basis b,..., b n of R n so that for all α,..., α n Z we have f(x+ n i= α ib i ) = f(x). The set L = { n i= α ib i : α i Z} is called the period lattice of f and L = {u R n : x u Z for all x L} is called the dual lattice of L.

4 4 F.M. de Oliveira Filho, F. Vallentin The mean value of a periodic function f is given by M(f) = lim T vol[ T,T] n f(x)dx. [ T,T] n For two periodic functions f and g we write f,g = M(fg). We say that f is square-integrable if f,f <. By f = f,f we denote its norm. If f and g are both square-integrable, then f,g exists. For u R n we define the Fourier coefficient f(u) = f,e iu x. Here, x y denotes the standard inner product on R n. Notice that the support of f is a discrete set, namely it lies in the dual lattice of the period lattice of f, scaled by 2π. If we let f y (x) = f(y + x) for a vector y R n, then f y (u) = f(u)e iu y. For square-integrable, periodic functions f and g Parseval s formula f,g = u R n f(u)ĝ(u) holds. By writing the latter sum we mean that we sum over the intersection of the supports of f and ĝ. Proof of Theorem.. Let A be a measurable subset of R n that avoids distances d,..., d N. By A we denote its characteristic function A : R n {0,} whose support is precisely A. Without loss of generality we can assume that A is a periodic function; in this case we say that A is periodic. Indeed, from any measurable set A which avoids distances d,..., d N we can construct a periodic set which avoids distances d,..., d N and with upper density arbitrarily close to the one of A. To do this we intersect A with a regular cube of side 2T so that vol(a [ T,T] n )/vol[ T,T] n is close to the upper density δ(a) and so that vol([ T +d,t d] n )/vol[ T,T] n, with d = max{d,...,d N }, differs from only negligibly. Then we construct the new periodic set by tiling R n with copies of A [ T +d,t d] n centered at the points of the lattice 2TZ n. Notice that, for a periodic set A, one may replace the limsup in the definition of δ(a) by a simple limit. By A y we denote the translation of the set A by the vector y R n so that its characteristic function satisfies A y (x) = A (x+y) = ( A ) y (x). The following two properties are crucial: (3) (4) A, = δ(a), A y, A = δ(a (A y)), for all y R n. In particular, we have A, A = δ(a) and A y, A = 0 for all vectors y of Euclideannormd,...,d N. Notice A, = A (0). ByapplyingParseval sformula to (4), we can express it in terms of the Fourier coefficients of A, thus obtaining u R n A (u) 2 e iu y Now we consider the function (5) ϕ(y) = A (0) = δ(a), = δ(a (A y)) for all y R n. u R n A (u) 2 e iu y = δ(a (A y)),

5 Fourier analysis, linear programming, and densities of distance avoiding sets 5 which is called the autocorrelation function (or two-point correlation function) of A. By taking spherical averages we construct from it a radial function f whose values only depend on the norm of the vectors. In other words, we set f(y) = ω n S n ϕ( y ξ)dω(ξ). Here ω denotes the standard surface measure on the unit sphere S n = {ξ R n : ξ ξ = } and ω n = ω(s n ) = (2π n/2 )/Γ(n/2). Clearly, f(0) = δ(a), and f(y) = 0 whenever y {d,...,d N }. Because of the formula (cf. Schoenberg [7, (.6)], see () for an explicit expression for Ω n ) e iu ξ dω(ξ) = Ω n ( u ) ω n S n we can represent f in the form f(y) = t 0α(t)Ω n (t y ), where α(t) is the sum of A (u) 2 for vectors u having norm t, so the α(t) s are real and nonnegative. Furthermore, α(0) = A (0) 2 = δ(a) 2 and t 0α(t) = f(0) = δ(a). So the following linear program in the variables α(t) gives an upper bound for the upper density of any measurable set which avoids the distances d,..., d N : (6) sup { α(0) : α(t) 0 for all t 0, α(t) =, t 0 α(t)ω n (td k ) = 0 for k =,..., N }. t 0 Above, all but a countable subset of the α(t) s are zero. Note moreover that we used the normalization t 0α(t) =. This linear program has infinitely many variables α(t) but only N + equality constraints. A dual of it is (7) inf { z 0 : z 0 +z +z 2 + +z N, z 0 +z Ω n (td )+z 2 Ω n (td 2 )+ +z N Ω n (td N ) 0 for all t > 0 }, which has N + variables z 0,z,z 2,...,z N and infinitely many constraints. As usual, weak duality holds between the pair of linear programs (6) and (7): If α(t) satisfies the conditions in (6) and if (z 0,z,...,z N ) satisfies the conditions in (7), then α(0) t 0α(t)(z 0 +z Ω n (td )+z 2 Ω n (td 2 )+ +z N Ω n (td N )) = z 0, which finishes the proof of our main theorem.

6 6 F.M. de Oliveira Filho, F. Vallentin 3. Sets avoiding one distance It is notable that the linear programming bounds for the extreme density of sets avoiding exactly one distance allow for an analytic optimal solution. Since this problem is scaling invariant we can assume that we consider sets avoiding the unit distance d =. Let j α,k be the k-th positive zero of the Bessel function J α. It is known that the absolute minimum of the function Ω n is attained at j n/2, (see Askey, Andrews, Roy [, (4.6.2)], and Watson [2, Chapter 5, 3]). So, the point (z 0,z ) which is determined by the equations z 0 +z = z 0 +z Ω n (j n/2, ) = 0 provides the optimal solution for the linear program in Theorem.. Hence, (8) z 0 = Ω n (j n/2, )/(Ω n (j n/2, ) ) m (R n ), and this gives by (2) a lower bound for the measurable chromatic number, namely χ m (R n ) /Ω n (j n/2, ). It is interesting to notice that this lower bound coincides with the one provided by Bachoc, Nebe, Oliveira, and Vallentin [2, Corollary 8.2], albeit by a shift of one dimension. This shift of one dimension is due to the fact that Bachoc, Nebe, Oliveira, and Vallentin [2] study the problem of sets avoiding one distance on the (n )-dimensional unit sphere S n R n and the lower bound for the measurable chromatic number χ m (R n ) was obtained by upper bounding the density of sets in the unit sphere which avoid the distance d where d goes to zero. So, we see now that this limit process gives a lower bound for the measurable chromatic number of R n and not only for the one of R n. 3.. Adding extra inequalities. It is possible to strengthen the main theorem and the resulting bound (8) by introducing extra inequalities. Consider a regular simplex in R n with edge length having vertices v,...,v n+. A set A R n which avoids the unit distance can only contain one vertex of this regular simplex. So we have for the autocorrelation function ϕ of the characteristic function A defined in (5) that (9) ϕ(v )+ +ϕ(v n+ ) = δ(a (A v ))+ +δ(a (A v n+ )) δ(a) = ϕ(0). Let O(R n ) be the n-dimensional orthogonal group, that is, the set of all n n real matrices Z such that Z t Z = I. Let µ denote the Haar measure over O(R n ) normalized by µ(o(r n )) =. Taking spherical averages of ϕ is the same as symmetrizing ϕ with respect to the orthogonal group, i.e., for all y R n, f(y) = ω n ϕ( y ξ)dω(ξ) = S n O(R n ) ϕ(zy) dµ(z) (this follows, e.g., from Theorem 3.7 in the book by Mattila [5]). Let f be, as above, the radial function obtained by the symmetrization of ϕ. For a nonnegative real number t, we write f(t) for the common value of f for vectors of length t. Then, by symmetrizing both sides of (9) with respect to the orthogonal group, and since distances are preserved by the action of O(R n ), we conclude that the inequality (0) f( v )+ +f( v n+ )

7 Fourier analysis, linear programming, and densities of distance avoiding sets 7 can be used to strengthen our original linear program. Here, observe that we already took into account the normalization f(0) =, introduced in (6). If we center a regular simplex at the origin, the above inequality specializes to (n+)f( /2 /(2n+2)), which gives the following strengthening of the dual formulation (7) { inf z 0 +z c : z c 0, z 0 +z +z c (n+), z 0 +z Ω n (t)+z c (n+)ω n (t /2 /(2n+2)) 0 } for all t 0. In Table 3. we give the new upper bounds on m (R n ) we get for n = 4,...,24 by solving the linear program on a computer (we discuss numerical issues at the end of this section) which are improvements over the values which Székely and Wormald give in [20]. This in turn gives new lower bounds for the measurable chromatic number for n = 4,...,24. n best upper bound new upper best lower bound new lower for m (R n ) bound for for χ m (R n ) bound for previously known m (R n ) previously known χ m (R n ) [8] [0] [20] [0] [20] [20] [20] [20] [20] [20] [20] [20] [20] [20] [20] [20] [20] [2] [20] [2] [20] [2] [20] [2] [20] [2] [20] [2] [20] [2] [20] [2] [20] [2] [20] [2] [20] [2] [20] [2] [20] [2] [20] [2] [20] [2] 206 Table 3.. Upper bounds for m (R n ) and lower bounds for χ m (R n ).

8 8 F.M. de Oliveira Filho, F. Vallentin However, in dimension 2 we only get an upper bound of To improve Székely s bound of 2/ in the plane, we replace the regular triangle centered at the origin by more triangles. We use the following three triples of squarednorms( v 2, v 2 2, v 3 2 ) for(0): (2.4,2.4, ),(3.,3., ) (3.7,3.7,7.474),where the last coordinate of (a,b,c) is a root of 3(a 2 +b 2 +c 2 + ) (a+b+c+) 2. This condition assures that the determinant of the positive semidefinite Gram matrix a 2 (a+b ) 2 (a+c ) 2 (a+b ) b 2 (b+c ) 2 (a+c ) 2 (b+c ) c of the points v,v 2,v 3 of a corresponding regular simplex vanishes. Solving the corresponding linear program yields the new upper bound of We found thethreetriplesbyconsideringalltriples(a,b,c)with a,b = 0.j with j = 0,...,40. In dimension 3 we use three quadruples ( v 2, v 2 2, v 3 2, v 4 2 ) of squared norms for (0): (0.3, 0.4, 0.4, ), (.9,.9,.9, ), (2, 2, 2, ), where the last coordinate of (a,b,c,d) is a root of 3(a 2 +b 2 +c 2 +d 2 +) 2(ab+ ac+ad+bc+bd+cd) 2(a+b+c+d). Solving the corresponding linear programming problem yields the new upper bound of We found the three quadruples by considering all triples (a,b,c,d) with a,b,c = 0.j with j = 0,..., Numerical calculations. A few technical remarks concerning the numerical calculations are in order. For solving the linear programs we use the software lpsolve [3] and we generate the input using the program GP/PARI [6]. We discretize the conditions of the form z 0 +z Ω n (t)+z c (n+)ω n (t /2 /(2n+2)) 0 for all t 0 by discretizing the interval [0, 20] into steps of size Now we demonstrate in the case n = 4 how we turn the numerical calculations into a rigorous mathematical proof: The linear program has the optimal numerical solution z 0 = , z = , z c = A lower bound of the minimum of the function z(t) = z 0 +z Ω 4 (t)+5z c Ω 4 ( 2/5t) in t [0,20] is The function z(t) is positive for t 20 because there Ω 4 (t) 0.02 and Ω 4 ( 2/5t) 0.04 holds. Thus by adding to z 0 we make sure that the new function z(t) is nonnegative. This only slightly effects the value of the bound. 4. Planar sets avoiding two distances In this section we quickly report on the problem of finding the smallest extreme density a measurable set in the plane can have which avoids exactly two distances, i.e., inf{m d,d 2 (R 2 ) : d,d 2 > 0}. Székely [8] showed that this number is at most 2/ by giving an upper bound for m, 3 (R 2 ). By solving the corresponding linear program on the computer we improve his bound to m, 3 (R 2 ) By adjusting the distances we can improve this further: m,j,2/j, (R 2 ) where j, and j,2 are the first two positive zeros of the Bessel function J.

9 Fourier analysis, linear programming, and densities of distance avoiding sets 9 By combining Bukh s result, which we explained in the introduction, with our new bound on m (R 2 ) from the previous section we can improve on this even further: inf{m d,d 2 (R 2 ) : d,d 2 > 0} ( m (R 2 ) ) Sets avoiding many distances In this section we give a proof of a variant of Bukh s result [5, Theorem ] about densities of sets avoiding many distances. His proof is based on a so-called zooming out lemma which resembles Szemerédi s regularity lemma for dense graphs, whereas our proof is an easy consequence of Theorem. and simple properties of the function Ω n. Theorem 5.. For every positive integer N there is a number r = r(n) strictly greater than such that for distances d,..., d N with () d 2 /d > r, d 3 /d 2 > r,..., d N /d N > r we have m d,...,d N (R n ) 2 N. In the proof of Theorem 5. some facts about the function Ω n will be useful. First, we have (2) J 0 (t), and J α (t) / 2 for all α > 0 and t 0 (cf. (4.9.3) in Andrews, Askey, and Roy []). From this, it follows at once that lim t Ω n (t) = 0 for n > 2. For n = 2 the same follows, e.g., from the asymptotic expansion for J α (cf. (4.8.5) in Andrews, Askey, and Roy []). Moreover, (3) Ω n (t) /2, for all n 2 and t 0. To see this, set α = (n 2)/2. It is known that for α > /2 we have J α (t)+j α+ (t) = 2α t J α(t) (cf. Andrews, Askey, and Roy [, (4.6.5)]). Combining this identity with () we obtain ( 2 ) α t Ω n (t) = Γ(α+) t 2α (J α (t)+j α+ (t)) (4) ( 2 ) α Jα+ = Ω n 2 (t)+γ(α) (t). t Now, recall that the globalminimum ofω n is attained at j α+,, the first positive zero of J α+ (cf. Section 3). From (4), we have that Ω n (j α+, ) = Ω n 2 (j α+, ). It follows that the minimum of Ω n is at least the minimum of Ω n 2. To finish the proof of (3) we have to check Ω 2 and Ω 3, what can be easily accomplished. Proof of Theorem 5.. Given N > 0, set ε = /(N2 N+ ). Since Ω n (0) = and since Ω n is continuous, there is a number t 0 > 0 such that Ω n (t) > ε for t t 0. Likewise, since lim t Ω n (t) = 0, there is a number t > t 0 such that Ω n (t) < ε for t t.

10 0 F.M. de Oliveira Filho, F. Vallentin Set r = r(n) = t /t 0 and let distances d,..., d N be given such that () is satisfied. With this we claim that, for j N, N 2 N i+ Ω n(td i ) (N j)ε. 2N j+2 i=j Before we prove the claim, we show how to apply it. By taking j = in the claim, and since by our choice of ε we have (N )ε /2 N+, it follows that N i= 2 N i+ Ω n(td i ) 2 N. Now we may set z 0 = /2 N and z i = /2 N i+ for i =,..., N and apply Theorem., proving our result. To finish, we prove the claim by induction. For j = N, the statement follows immediately from (3). Now, suppose the statement is true for some < j N. We show that the statement is also true for j by distinguishing two cases. First, for t t 0 /d j, we have from the choice of t 0 that Ω n (td j ) > ε. Using this and the induction hypothesis, we then have that N i=j 2 N i+ Ω n(td i ) = 2 N j+2 Ω n(td j )+ N i=j 2 N i+ Ω n(td i ) ε 2 N j+2 (N j)ε 2N j+2 (N j +)ε. 2N j+3 Now suppose t t 0 /d j. Observe that, for j i N, we have td i t 0 d i /d j t 0 r = t, hence Ω n (td i ) < ε. So, by using (3), we have N i=j 2 N i+ Ω n(td i ) = 2 N j+2 Ω n(td j )+ (N j +)ε, 2N j+3 finishing the proof of the claim. Acknowledgements N i=j 2 N i+ Ω n(td i ) We thank Christine Bachoc for helpful discussions, Noam D. Elkies for enlightening comments, and Boris Bukh and Lex Schrijver for their valuable suggestions on our manuscript. References [] G.E. Andrews, R. Askey, R. Roy, Special functions, Cambridge University Press, 999. [2] C. Bachoc, G. Nebe, F.M. de Oliveira Filho, F. Vallentin, Lower bounds for measurable chromatic numbers, arxiv: v2 [math.co], 8 pages, to appear in Geom. Funct. Anal. [3] M. Berkelaar, K. Eikland, P. Notebaert, lpsolve version 5.5, available from the web site [4] J. Bourgain, A Szemerédi type theorem for sets of positive density in R k, Israel J. Math. 54 (986),

11 Fourier analysis, linear programming, and densities of distance avoiding sets [5] B. Bukh, Measurable sets with excluded distances, arxiv:math/ v2 [math.co], 23 pages, to appear in Geom. Funct. Anal. [6] H. Cohn, N.D. Elkies, New upper bounds on sphere packings I, Ann. of Math. 57 (2003), [7] H.T. Croft, Incidence incidents, Eureka (Cambridge) 30 (967), [8] H.T. Croft, K.J. Falconer, R.K. Guy, Unsolved problems in geometry, Springer, 99. [9] P. Delsarte, An algebraic approach to the association schemes of coding theory, Philips Res. Rep. Suppl. (973), vi+97. [0] K.J. Falconer, The realization of distances in measurable subsets covering R n, J. Combin. Theory Ser. A 3 (98), [] K.J. Falconer, The realization of small distances in plane sets of positive measure, Bull. London Math. Soc. 8 (986), [2] K.J. Falconer, J.M. Marstrand, Plane sets with positive density at infinity contain all large distances, Bull. London Math. Soc. 8 (986), [3] H. Furstenberg, Y. Katznelson, B. Weiss, Ergodic theory and configurations in sets of positive density, p in Mathematics of Ramsey theory (J. Nešetřil, V. Rödl ed.), Springer, 989. [4] Y. Katznelson, An introduction to harmonic analysis, Dover Publications, 976. [5] P. Mattila, Geometry of Sets and Measures in Euclidean Spaces: Fractals and rectifiability, Cambridge University Press, 995. [6] The PARI-group, PARI/GP version 2.3.3, Bordeaux 2006, available from the web site [7] I.J. Schoenberg, Metric spaces and completely monotone functions, Ann. of Math. 39(938), [8] L.A. Székely, Measurable chromatic number of geometric graphs and sets without some distances in Euclidean space, Combinatorica 4 (984), [9] L.A. Székely, Erdős on unit distances and the Szemerédi-Trotter theorems, p in Paul Erdős and his mathematics (G. Halász, L. Lovász, M. Simonovits, V.T. Sós ed.), Springer, [20] L.A. Székely, N.C. Wormald, Bounds on the measurable chromatic number of R n, Discrete Math. 75 (989), [2] G.N. Watson, A treatise on the theory of Bessel function, Cambridge Mathematical Library, 995. F.M. de Oliveira Filho, Centrum voor Wiskunde en Informatica (CWI), Kruislaan 43, 098 SJ Amsterdam, The Netherlands address: f.m.de.oliveira.filho@cwi.nl F. Vallentin, Centrum voor Wiskunde en Informatica (CWI), Kruislaan 43, 098 SJ Amsterdam, The Netherlands address: f.vallentin@cwi.nl