GAPS BETWEEN LOGS JENS MARKLOF AND ANDREAS STRÖMBERGSSON

Similar documents
The two-point correlation function of the fractional parts of n is Poisson

arxiv: v2 [math.nt] 27 Mar 2018

arxiv: v2 [math.nt] 21 Jun 2017

Applications of measure rigidity: from number theory to statistical mechanics

Abstract. 2. We construct several transcendental numbers.

Applications of homogeneous dynamics: from number theory to statistical mechanics

Solving a linear equation in a set of integers II

BURGESS INEQUALITY IN F p 2. Mei-Chu Chang

Notes on Equidistribution

ON THE GAPS BETWEEN ZEROS OF TRIGONOMETRIC POLYNOMIALS

ON SPECTRAL CANTOR MEASURES. 1. Introduction

The Berry-Tabor conjecture

arxiv: v1 [math.co] 7 Jul 2014

HOMEWORK ASSIGNMENT 6

A combinatorial problem related to Mahler s measure

REAL AND COMPLEX ANALYSIS

Longest element of a finite Coxeter group

Character sums with Beatty sequences on Burgess-type intervals

A Remark on Sieving in Biased Coin Convolutions

arxiv: v1 [math.nt] 24 Mar 2019

DIGITAL EXPANSION OF EXPONENTIAL SEQUENCES

On Systems of Diagonal Forms II

Chapter 8. P-adic numbers. 8.1 Absolute values

Lecture II: The Lorentz gas

Theorem 5.3. Let E/F, E = F (u), be a simple field extension. Then u is algebraic if and only if E/F is finite. In this case, [E : F ] = deg f u.

NORMAL NUMBERS AND UNIFORM DISTRIBUTION (WEEKS 1-3) OPEN PROBLEMS IN NUMBER THEORY SPRING 2018, TEL AVIV UNIVERSITY

arxiv: v1 [math.nt] 5 Mar 2019

NON-UNIVERSALITY OF THE NAZAROV-SODIN CONSTANT

Reciprocity formulae for general Dedekind Rademacher sums

Quadratic reciprocity (after Weil) 1. Standard set-up and Poisson summation

Decomposition of Riesz frames and wavelets into a finite union of linearly independent sets

Real Analysis Math 131AH Rudin, Chapter #1. Dominique Abdi

A Characterization of (3+1)-Free Posets

CHAPTER 3 Further properties of splines and B-splines

On Some Mean Value Results for the Zeta-Function and a Divisor Problem

The density of rational points on non-singular hypersurfaces, I

MORE ON CONTINUOUS FUNCTIONS AND SETS

Variations on Quantum Ergodic Theorems. Michael Taylor

We denote the space of distributions on Ω by D ( Ω) 2.

ON THE LIMIT POINTS OF THE FRACTIONAL PARTS OF POWERS OF PISOT NUMBERS

INTRODUCTION TO REAL ANALYTIC GEOMETRY

Fourier series, Weyl equidistribution. 1. Dirichlet s pigeon-hole principle, approximation theorem

Arithmetic progressions in sumsets

Logarithmic scaling of planar random walk s local times

UNIFORM DISTRIBUTION MODULO 1 AND THE UNIVERSALITY OF ZETA-FUNCTIONS OF CERTAIN CUSP FORMS. Antanas Laurinčikas

COMBINATORIAL CURVE NEIGHBORHOODS FOR THE AFFINE FLAG MANIFOLD OF TYPE A 1 1

Quantitative Oppenheim theorem. Eigenvalue spacing for rectangular billiards. Pair correlation for higher degree diagonal forms

Ergodic Theorems. Samy Tindel. Purdue University. Probability Theory 2 - MA 539. Taken from Probability: Theory and examples by R.

THE AUTOCORRELATION OF THE MÖBIUS FUNCTION AND CHOWLA S CONJECTURE FOR THE RATIONAL FUNCTION FIELD

arxiv: v5 [math.ca] 23 Jan 2016

The expansion of random regular graphs

Discrete Math, Fourteenth Problem Set (July 18)

Mi-Hwa Ko. t=1 Z t is true. j=0

Course 311: Michaelmas Term 2005 Part III: Topics in Commutative Algebra

Zeros of lacunary random polynomials

The Degree of the Splitting Field of a Random Polynomial over a Finite Field

On the modification of the universality of the Hurwitz zeta-function

Low-discrepancy sequences obtained from algebraic function fields over finite fields

WEIERSTRASS THEOREMS AND RINGS OF HOLOMORPHIC FUNCTIONS

A NOTE ON SIMPLE DOMAINS OF GK DIMENSION TWO

On the law of the iterated logarithm for the discrepancy of lacunary sequences

THE S 1 -EQUIVARIANT COHOMOLOGY RINGS OF (n k, k) SPRINGER VARIETIES

A NOTE ON ZERO SETS OF FRACTIONAL SOBOLEV FUNCTIONS WITH NEGATIVE POWER OF INTEGRABILITY. 1. Introduction

On stable inversion of the attenuated Radon transform with half data Jan Boman. We shall consider weighted Radon transforms of the form

The Hilbert Transform and Fine Continuity

The primitive root theorem

Computing a Lower Bound for the Canonical Height on Elliptic Curves over Q

P -adic root separation for quadratic and cubic polynomials

LIMITS FOR QUEUES AS THE WAITING ROOM GROWS. Bell Communications Research AT&T Bell Laboratories Red Bank, NJ Murray Hill, NJ 07974

arxiv: v1 [math.co] 3 Nov 2014

Plane quartics and. Dedicated to Professor S. Koizumi for his 70th birthday. by Tetsuji Shioda

MEAN CURVATURE FLOW OF ENTIRE GRAPHS EVOLVING AWAY FROM THE HEAT FLOW

ON THE LINEAR INDEPENDENCE OF THE SET OF DIRICHLET EXPONENTS

AN EXPLORATION OF THE METRIZABILITY OF TOPOLOGICAL SPACES

LIMITING CASES OF BOARDMAN S FIVE HALVES THEOREM

Radial balanced metrics on the unit disk

Distinguishing Cartesian powers of graphs

Integral Bases. 1. Resultants

CHAPTER 0 PRELIMINARY MATERIAL. Paul Vojta. University of California, Berkeley. 18 February 1998

ZEROS OF SPARSE POLYNOMIALS OVER LOCAL FIELDS OF CHARACTERISTIC p

MATH 8253 ALGEBRAIC GEOMETRY WEEK 12

1 Definition of the Riemann integral

Propagating terraces and the dynamics of front-like solutions of reaction-diffusion equations on R

Van der Corput sets with respect to compact groups

MASTERS EXAMINATION IN MATHEMATICS SOLUTIONS

the neumann-cheeger constant of the jungle gym

Weighted Sums of Orthogonal Polynomials Related to Birth-Death Processes with Killing

f(x n ) [0,1[ s f(x) dx.

CLASSIFICATIONS OF THE FLOWS OF LINEAR ODE

ESTIMATES FOR ELLIPTIC HOMOGENIZATION PROBLEMS IN NONSMOOTH DOMAINS. Zhongwei Shen

Bounds for the Eventual Positivity of Difference Functions of Partitions into Prime Powers

PAIR CORRELATION FOR QUADRATIC POLYNOMIALS MOD 1

Möbius Randomness and Dynamics

MATH 425, FINAL EXAM SOLUTIONS

u( x) = g( y) ds y ( 1 ) U solves u = 0 in U; u = 0 on U. ( 3)

The Polynomial Numerical Index of L p (µ)

Empirical Processes: General Weak Convergence Theory

MATH 51H Section 4. October 16, Recall what it means for a function between metric spaces to be continuous:

Quadratic reciprocity (after Weil) 1. Standard set-up and Poisson summation

Real Analysis Problems

Transcription:

GAPS BETWEE LOGS JES MARKLOF AD ADREAS STRÖMBERGSSO Abstract. We calculate the limiting gap distribution for the fractional parts of log n, where n runs through all positive integers. By rescaling the sequence, the proof quickly reduces to an argument used by Barra and Gaspard in the context of level spacing statistics for quantum graphs. The key ingredient is Weyl equidistribution of irrational translations on multi-dimensional tori. Our results extend to logarithms with arbitrary base; we deduce explicit formulas when the base is transcendental or the rth root of an integer. If the base is close to one, the gap distribution is close to the exponential distribution.. The gap distribution is a popular measure to quantify the degree of randomness in a given deterministic sequence. Rudnick and Zaharescu have shown [2] that the gap distribution of the fractional parts of 2 n α (n =,..., ), for almost all α, converges in the limit to an exponential distribution the gap distribution of a Poisson point process. (2 n may be replaced by any lacunary sequence.) The same is expected to hold for the fractional parts of n k α, for any integer k 2 and α of bounded type, but this remains unproven [,, 8, 4]. umerical experiments also suggest that the fractional parts of n β may have an exponential gap distribution provided β (, ) (, ). In the 2 2 case β =, Elkies and McMullen [3] have calculated an explicit formula for the limiting 2 gap distribution, which is evidently not exponential. Remarkably, their distribution coincides with the limiting gap distribution of directions of points in the affine lattice Z 2 +a for any fixed vector a / Q 2 [9]. Both of these findings follow from the equidistribution of translates of certain curves on the homogeneous space ASL(2, Z)\ ASL(2, R). In the present paper we show that also the fractional parts of n have a non-exponential limiting gap distribution, which we calculate explicitly. Our derivation reduces quickly to an argument used by Barra and Gaspard [] in the context of spectral statistics of quantum graphs. The key ingredient here is Weyl equidistribution on multi-dimensional tori, similar in spirit to [5, Chap. 3]. 2. To state our main results, let us denote by ξ n the fractional parts of n. We denote by {ξ n, } n=,..., the ordered set of the first elements of ξ n, so that () ξ, ξ 2,... ξ, <. For purely notational reasons it will be convenient to set ξ +, := ξ, +, and also consider the gap between the first and last element mod ; this additional gap has of course no effect on the existence and form of the limiting gap distribution. For given, the gap distribution P (s) is defined as the probability density on R, (2) P (s) = δ(s (ξ n+, ξ n, )). n= Date: 6 January 23. 2 Mathematics Subject Classification. K6, J7. J.M. is supported by a Royal Society Wolfson Research Merit Award and ERC Advanced Grant HFAKT. A.S. is a Royal Swedish Academy of Sciences Research Fellow supported by a grant from the Knut and Alice Wallenberg Foundation.

2 JES MARKLOF AD ADREAS STRÖMBERGSSO..8.6.4.2.5..5 2. 2.5 3. Figure. The distribution of gaps between the fractional parts of log n, where n =,..., 4. The piecewise continuous curve is the limit distribution P (s) of Theorem for b = e...8.6.4.2.5..5 2. 2.5 3. Figure 2. The distribution of gaps between the fractional parts of n with base b = e /5 =.2242758..., where n =,..., 4. The blue curve is the limit distribution P (s) of Theorem, the red curve is the exponential distribution e s. We denote by T the set of transcendental numbers b > ; b = e is a natural example, cf. Figure. The technique we present here also works for algebraic b, but in general leads to more intricate limit distributions. In we discuss the simple case of integer base b, and in the case when the base is the rth root of an integer. Theorem. Let b T. For any bounded continuous f : R R, (3) lim f(s) P (s) ds = f(s) P (s) ds,

GAPS BETWEE LOGS 3 where s 2 (F (s ) F (sb )) if < s < (4) P (s) = s 2 F (sb ) if < s < if s > b (5) F (a) = a a (a; b ) (a; b ) and (6) (a; q) = (7), ( aq n ) (q-pochhammer symbol). n= The derivative of the Pochhammer symbol is explicitly a (a; q) = (a; q) q j a. ote that P (s) has jump discontinuities at s = and b ; however both the left and right limits exist at these points, and are finite. We also note that lim s + P (s) =. b A more elementary formulation, which is equivalent to the above, is that for every s (8) lim j= {n : (ξ n+, ξ n, ) > s} = s P (s ) ds. We will show in the final paragraph 2 of this paper that, in the limit b, P (s) converges to the exponential distribution e s, cf. Figure 2. 3. To gain a better insight into the limiting gap distribution P (s) for fixed b, let us recall that the fractional parts of n are not uniformly distributed mod [6]. However, the sequence of η n (n =,..., ) given by the fractional parts of (n/) = n has a limit density mod, and evidently the same gap distribution (to see this, recall that we have added the gap between the first and last element of the sequence mod, and the ordering in the sequence remains the same). It is an easy exercise to show that for any interval A [, ], (9) lim {n : ηn A} = ρ(x) dx with density () ρ(x) = b bx. 4. A more refined statement of Theorem is to consider the joint distribution of η n and the subsequent gap; define the corresponding probability density on [, ] R by () P (x, s) = δ(x η n, ) δ(s (η n+, η n, )), n= where {η n, } n=,..., is again the ordered set of the first elements of η n. The joint distribution allows us to compare the statistics of gaps between elements near one point x of the unit interval with the gap statistics at a second point x. A b

4 JES MARKLOF AD ADREAS STRÖMBERGSSO The explicit formula for the limit distribution will involve the density 2 (2) R(a, b) = (b ; b ) δ(a ) + a (a; 2 b ) if < a < if a >. The second derivative of the Pochhammer symbol is 2 (3) a (a; q) = (a; q) 2 (q j a)(q k a). j,k= j k Theorem 2. Let b T. For any bounded continuous f : [, ] R R, (4) lim where f(x, s) P (x, s) ds dx = (5) P (x, s) = (b x ) 2 R(sb x, b). f(x, s) P (x, s) ds dx, 5. To unravel the previous statements further, let us consider a strictly increasing function [, ] [, ] which maps the sequence η j to a new sequence η j which now is uniformly distributed mod. This function is given by (cf. [7]) (6) x and thus x ρ(y) dy = bx b, (7) η n = bηn b, η n, = bηn, b. We have, explicitly, for n [b k, b k+ ) and k Z, (8) η n = bk n (b ). We will see in 8 that the knowledge of the joint distribution P (x, s) for the original sequence and the analogue for the rescaled and ordered sequence, (9) P (x, s) = δ(x η n, ) δ(s ( η n+, η n, )) are equivalent: n= Theorem 3. Let b T. For any bounded continuous f : [, ] R R, (2) lim where f(x, s) P (x, s) ds dx = (2) P (s) = ( b ) 2 R ( ( b )s, b ). f(x, s) P (s) ds dx, ote that the limit distribution is independent of x. This is not a consequence of η j being uniformly distributed mod, but reflects the fact that the gap distribution is the same if we restrict the sequence to an arbitrary subinterval of [, ]. We will first prove Theorem 3 in 6 and 7, and then retrace our steps by subsequently showing in 8 that Theorem 3 implies Theorem 2, where (22) P (x, s) = P (ρ(x)s)ρ(x) 2 = (b x ) 2 R(sb x, b).

GAPS BETWEE LOGS 5 We obtain Theorem from Theorem 2 by taking an x-independent test function. 6. We now turn to the proof of Theorem 3 and study a more general setting. Consider a finite sequence ω,..., ω J of positive real numbers (the frequencies ). With each ω j we associate a real number β j. We are interested in the statistical properties of the multiset J (23) S = (β j + ω Z). j= Here denotes multiset sum, i.e. a union where we keep track of the multiplicity of each element (this is only relevant if the sets β j + ω j Z, j =,..., J are not mutually disjoint). The number of elements of S that fall into the interval I = [t L, t + L] is 2 2 (24) (t, L) = J j (t, L), j= j j (t, L) = I (βj + ω j Z). Evidently j (t, L) is periodic in t with period ω j. Furthermore (25) ω j L < j (t, L) ω j L +. Therefore (26) lim L (t, L) L = J ω j, which means that the asymptotic density of S in R equals j ω j. For any x R we write x Z [, ] for the distance between x and the nearest 2 integer. We also set (27) n(x, L) = [x L 2, x + L 2 ] Z. The verification of the following lemma is a simple exercise. Lemma. Given x R, L and k Z, we have n(x, L) = k if and only if x k 2 Z is both (k L) and > (L k). 2 2 In particular, if x is picked at random with respect to the uniform probability measure P on [, ], then { k L if L < k < L + (28) P (n(x, L) = k) = E (k, L) := otherwise. The following observation is due to Barra and Gaspard []. Theorem 4. Assume that the frequencies ω,..., ω J are linearly independent over Q. Then, for < a < b < and k Z, meas({t [at, bt ] : (t, L) = k}) (29) lim T (b a)t where (3) E(k, L) = k +...+k J =k j= j= J ( E kj, ω j L ). = E(k, L)

6 JES MARKLOF AD ADREAS STRÖMBERGSSO In particular, { J ( j= ωj L ) if L < max(ω,..., ω J ) (3) E(, L) = otherwise. The convergence in (29) is uniform over all choices of β,..., β J ω,..., ω J fixed). R (but keeping Proof. Decomposing [at, bt ] into intervals of length ω and using (t, L) = J j= j(t, L), we see that meas({t [at, bt ] : (t, L) = k}) equals (32) k +...+k J =k n b n=n a ω meas({s [, ) : j (ω (s + n), L) = k j for j =,..., J}) + E, where n a = ω at, n b = ω bt, and E is a real number whose absolute value is bounded above by the measure of the symmetric difference between [at, bt ] and [n a ω, n b ω ]; thus E 2ω. Let us define α j,n R/Z through the relation β j + ω j Z = nω + ω j (α j,n + Z). Then j (ω (s + n), L) = k j holds if and only if [ω s L, 2 ω s + L] 2 ω j (α j,n + Z) = k j, and by Lemma this holds if and only if 2 ω j ω s α j,n k 2 j Z is both k j ω j L and > ω j L k j. It follows that meas({t [at, bt ] : (t, L) = k}) equals (33) k +...+k J =k n b n=n a ω f k,...,k J (α,n, α 2,n,..., α J,n ) + E, where f k,...,k J : (R/Z) J [, ] is defined by J ( f k,...,k J (α, α 2,..., α J ) = I 2 ω j s α j ) ω 2 k (34) j > k j ω j L ds. Z j= ote that each function f k,...,k J is continuous on (R/Z) J. It follows from our definition of α j,n that α j,n+ = α j,n ω j ω mod Z. In particular α,n is independent of n; in fact α,n = α := ω β in R/Z for all n. On the other hand our assumption implies that the numbers, ω 2 ω, ω 3 ω,..., ω J ω are linearly independent over Q, and hence by Weyl equidistribution we have nb (35) n=n a f k,...,k J (α,n, α 2,n,..., α J,n ) n b n a f k,...,k J (α, α 2,..., α J ) dα 2... dα J (R/Z) J as n b n a. ote that the convergence here is uniform over all choices of β,..., β J R. The right hand side equals J ( I 2 ω j s α j ) ω 2 k (36) j > k j ω j L ds, Z (R/Z) J j= and changing order of integration and then substituting α j := x j + ω j ω s k 2 j (x j R/Z) we see that the expression factors as J j= E (k j, ω j L). Finally noticing also that n b n a ω (b a)t as T, we obtain (29). We now turn to the gap distribution. We order the elements of S and label them as (37)... λ 2 λ λ λ λ 2....

GAPS BETWEE LOGS 7 Theorem 5. Assume that the frequencies ω j are linearly independent over Q, and that ω is the largest among ω,..., ω J. Then, for < a < b <, and any s, s ω, (38) lim T (b a)t {j Z : λ j [at, bt ] and λ j+ λ j > s} = s P ω (s )ds where J [ J ( J ( (39) P ω (s) = ω ω ) ω h ω i ωj s )] if < s < ω j δ(s ω )+ h,i= j= ω j=2 h i j h,i if s > ω The convergence in (38) is uniform over all choices of β,..., β J ω,..., ω J fixed).. R (but keeping Proof. This is a corollary of Theorem 4, cf. e.g. Theorem 2.2 in [7]. The limiting gap density is given by the formula (4) P ω (s) = d2 E(, s) ds 2. 7. Proof of Theorem 3. It suffices to prove that (2) holds when f is the characteristic function of a box [, a] [, A] with A ( b ), i.e. to prove that for any fixed a [, ] and A, A ( b ), we have (4) G a,a () a A P (s) ds as, where G a,a () := { (42) n {,..., } : η n, [, a], ( η n+, η n, ) [, A]} (cf., e.g., [2, Example 2.3]). Set ω j = b j (b ) and β () j = β () = for j =, 2,.... For any given J, Z b we note that by (8), { η,..., η } equals the multiset sum of (β () +ω j Z) [, ) for j =, 2,..., J and a remaining multiset of cardinality b J consisting of η n for n = and all n < b J. Hence if we order the elements of S := J j=(β () + ω j Z) as... λ 2 λ λ λ λ 2..., then G a,a () = { ( (43) j Z : λ j [, a], λ j+ λ j [, A]} + O + b J). Since b is transcendental, the frequencies ω,..., ω J are linearly independent over Q; thus Theorem 5 applies, and we obtain, for any fixed J Z, (44) lim sup G a,a () a lim inf G a,a() a A A P ω,j (s ) ds + O(b J ); P ω,j (s ) ds O(b J ),

8 JES MARKLOF AD ADREAS STRÖMBERGSSO where P ω,j (s) = if s > ω while for s < ω we have J ( P ω,j (s) = ω ω ) j (45) δ(s ω ) + d2 ω ds 2 We now compute that (46) A j=2 P ω,j (s ) ds = J ( ωj s ). j= if A > ω J ω j J J ω j ( ω i A) if A < ω. j= j= i= (i j) Hence, letting J in (44) and using ω j = b j (b ), we conclude (47) lim G a,a() = a + a if A > b ( ( ( b )s; b )) if A < b s s=a b b = a A P (s) ds, where the last equality follows by a direct computation using (2) and (2). Hence (4) holds, and Theorem 3 is proved. 8. Let us now prove that Theorem 3 implies Theorem 2. (An entirely analogous argument also shows that Theorem 2 implies Theorem 3, so that the statements of these two theorems are in fact equivalent.) Given a bounded continuous function f : [, ] R R, we wish to prove that (4) holds. Since the sequence of probability measures P (x, s) ds dx is tight (cf., e.g., Lemma 2. in [7]), a familiar approximation argument shows that without loss of generality we may assume that f has compact support, i.e. there is some B > such that f(x, s) = whenever s B. Let T : [, ] [, ] be the map in (6) and define T : [, ] R [, ] R through T(x, s) = (T (x), T (x)s). This is a bijection, with inverse given by (48) T (x, s) = ( T (x), s T (T (x)) ow apply Theorem 3 with the test function f T. Then in (2) we have f T (x, s) P (x, s) ds dx = (49) f(η n,, s n, ), where s n, := T (η n+,) T (η n, ). Let F T be the set of those n {,..., } for (η n, ) which η n+, η n, Bb. Then since T (x) = ρ(x) is continuous and bounded away from zero on [, ] we have s n, = (η n+, η n, ) + o() uniformly over all n F as, and hence, since f is uniformly continuous, (5) n= f(η n,, s n, ) = f(η n,, (η n+, η n, )) + o(), uniformly over all n F. On the other hand if n / F then (η n+, η n, ) > Bb > B and also s n, > Bb inf T = B, so that both sides of (5) vanish. Hence (5) in sup T fact holds uniformly over all n {,..., } as, and it follows that the limit in the left hand side of (4) exists, and equals the limit in (2) (with f T ). Finally substituting (x, s) = T(x, s ) in the right hand side of (2) we see (via (22)) that this limit equals the right hand side of (4), and the proof is complete. ).

GAPS BETWEE LOGS 9 2.5 2..5..5.5..5 2. 2.5 3. Figure 3. The distribution of gaps between the fractional parts of n with base b =, where n =,..., 4. The piecewise continuous curve is the limit distribution P (s) in (56). 9. Theorem is now a direct consequence of Theorem 2. As to the explicit formula for the gap distribution, (5) P (s) = = P (x, s) dx (b x ) 2 R(sb x, b) dx. With the variable substitution a = sb x, da = sb x () 2 dx = () a dx we have s (52) P (s) = R(a, b) a da. s 2 sb Recall the definition of R(a, b) in (2). The integral over the first term in (2) yields (b ;b ) if sb < < s and otherwise. For the second term, we use integration by parts, s 2 a [ a 2 a a (a; 2 b ) da = a ] a a a (a; b ) a a a (a; b ) da (53) [ = a ] a a (a; b ) (a; b ) where a = min{, sb } and a = min{, s }. This yields (4).. As we have noted, the technique presented here works also for algebraic b. Let us consider briefly the simplest case of b > being an integer. In this case, the multiset sum S = J j=(β () + ω j Z) considered in the proof of Theorem 3 has all its points lying in β () + ω Z, i.e. the only gap lengths appearing are and ω = ( b ) ; furthermore S is periodic with period ω J = b J ω. One finds that the average multiplicity of β () + ω k in S for varying k Z equals b J. Hence following the b proof of Theorem 3 and taking the limit J, we conclude that for b > an integer, a,

JES MARKLOF AD ADREAS STRÖMBERGSSO.2..8.6.4.2.5..5 2. 2.5 3. Figure 4. The distribution of gaps between the fractional parts of n with base b =, where n =,..., 4. The piecewise continuous curve is the limit distribution P (s) in (62) with r = 2..2..8.6.4.2.5..5 2. 2.5 3. Figure 5. The distribution of gaps between the fractional parts of log n (= n with base b = / =.2589254...), where n =,..., 4. The blue curve is the limit distribution P (s) in (62) with r =. The red curve is the exponential distribution e s. we have a limit result as in (2) in Theorem 3, but with the limit density being given simply by (54) P (s) = b δ(s) + ( b )δ ( s ( b ) ). As to the limit of the unscaled sequence η n, the first equality in (22) yields (55) P (x, s) = b bx δ(s) + ()b x δ (s b x ),

and therefore (56) P (s) = b δ(s) + GAPS BETWEE LOGS { s 2 if < s < otherwise. See Figure 3 for an illustration in the case b =.. The findings of the previous paragraph can be generalised to bases of the form b = m /r, with m and r positive integers such that m /p is not an integer for any prime divisor p of r. In this case, the multiset sum S = J j=(β () + ω j Z) considered in the proof of Theorem 3 has all points in the set r j=(β () + ω j Z). Our assumption implies that the polynomial X r m is irreducible over Q (since any monic polynomial of degree d < r dividing X r m must have zeroth coefficient of absolute value m d/r, which is irrational). Hence ω,..., ω r are linearly independent over Q, and we now see by the same argument as in the proof of Theorem 4 that the limit (29) exists in the case k = and is given by { r ( j= b j (b )L ) if L < ( b ) (57) E(, L) = otherwise. Using the finite q-pochhammer symbol r (58) (a; q) r = ( aq n ), we can write this as (59) E(, L) = We note that for L < ( b ) n= b { (( b )L; b ) r if L < ( b ) otherwise. (6) d dl E(, L) = ( r b )E(, L) b j ( b )L and thus (6) d E(, L) dl = ( b ) L= j= r j= b j = b r. The fact that this value is less than is due to the non-trivial multiplicity of values in our sequence, which results in a positive density of zero gaps. The statements of Theorems 3 therefore hold with the following limit distributions. The limiting gap distribution of Theorem is s 2 (F r (s ) F r (sb )) if < s < (62) P (s) = b r δ(s) + s 2 F r (sb ) if < s < if s > b with (63) F r (a) = a a (a; b ) r (a; b ) r. The joint limiting gap distribution of Theorem 2 is (64) P (x, s) = b bx r δ(s) + (b x ) 2 R r (sb x, b), b

2 JES MARKLOF AD ADREAS STRÖMBERGSSO and the rescaled gap distribution of Theorem 3 is (65) P (s) = b r δ(s) + ( b ) 2 R r ( ( b )s, b ), with 2 (66) R r (a, b) = (b ; b ) r δ(a ) + a (a; 2 b ) r if < a < if a >. See Figures 4 and 5 for illustrations in the cases b = /2 and b = /. 2. We conclude this investigation with a few remarks on the family of limit distributions which we have obtained for different values of b T. It is helpful to introduce the random point process in R given by the sequence (67) ϕ (b) = {( η n + t + m) : n, m Z} where t is uniformly distributed in [, ]. This process is clearly stationary and has intensity one. By using the same strategy as in the proof of Theorem 3, one can deduce from Theorem 4 a limit law for the one-dimensional distribution of ϕ (b) : For every closed interval A R, k Z, (68) lim P( ϕ(b) A = k) = E(b) (k, meas(a)), where (69) E (b) (k, L) = lim J k +...+k J =k j= J ( E kj, (b )b j L ), with E (k, L) as in (28). In the case k =, relation (68) in fact follows directly from Theorem 3 (exploiting again Theorem 2.2 in [7]), where (7) P (b) (s) = d2 E (b) (, s) ds 2. Although (68) requires b T, the family of distributions (69) is well defined for any real b (, ). Let us briefly analyse the limiting cases b and b. For b, E (b) (k, L) converges to the statistics of the point process given by the randomly shifted integer lattice Z + t, with t uniformly distributed in [, ]. That is, E (b) (k, L) E (k, L) as in (28). The rescaled gap distribution satisfies (7) P (b) (s) R(s, ) = δ(s ). The limit of the raw P (b) (s) distribution is slightly more involved since the asymptotic density ρ (b) (x) in () converges to δ(x ). We therefore cannot expect P (b) (s) to converge to a non-trivial limit without further rescaling. Let us extend ρ (b) (x) to a probability density on R by setting ρ (b) (x) = for x >. We have the scaling limit ( ) b x (72) ρ(b) e x, which, however, is not a probability density on R. Working directly from (4), we see that for s fixed, ( ) { s (73) P (b) if < s < s 2 if s >.

GAPS BETWEE LOGS 3 In the limit b, E (b) (k, L) converges to the Poisson distribution. That is, (74) E (b) (k, L) Lk k! e L, P (b) (s) e s. To see this, note that for L fixed and b sufficiently close to, (75) log E (b) (, L) = log( (b )b j L) = L + O(b ), j= which proves the claim for k = and the gap distribution P (s) (via Theorem 2.2 in [7]). Again, for L fixed and b sufficiently close to, we have for the remaining cases k, E (b) (k, L) = E (b) E (, (b (, L) )b j L) E (, (b L) )b jk E j <...<j (, (b )b j L) E (, (b )b j k L) k (76) = (b ) k L k E (b) (, L) b (j+...+jk) ( (b )b j j <...<j L) ( (b )b j kl) k ( ) = L k E (b) (, L) + O(b ). k! As to the limit of the unscaled sequence η n, we have for the density (77) ρ (b) (x), x [, ], i.e., we have uniform distribution mod one. In view of (22), P (b) (x, s) e s and thus, for the raw gap distribution, P (b) (s) e s. The above discussion can be readily adapted to bases of the form b = m /r considered in. In particular, we again observe that both the raw and rescaled limiting gap distributions (62) and (65) converge to the exponential distribution e s when r for fixed m, cf. Fig. 5. Acknowledgments. We thank Jon Keating and Zeév Rudnick for their comments on the first draft of this paper. References [] F. Barra and P. Gaspard, On the level spacing distribution in quantum graphs. J. Statist. Phys. (2) 283-39. [2] P. Billingsley, Convergence of probability measures, John Wiley & Sons Inc., ew York, 999. [3].D. Elkies and C.T. McMullen, Gaps in n mod and ergodic theory. Duke Math. J. 23 (24) 95 39. [4] D.R. Heath-Brown, Pair correlation for fractional parts of αn 2. Math. Proc. Cambridge Philos. Soc. 48 (2) 385-47. [5] M. Kac, Statistical independence in probability, analysis and number theory. The Carus Mathematical Monographs, o. 2. Published by the Mathematical Association of America. Distributed by John Wiley and Sons Inc., ew York 959. [6] L. Kuipers and H. iederreiter, Uniform distribution of sequences. Pure and Applied Mathematics. Wiley-Interscience, ew York-London-Sydney, 974. [7] J. Marklof, Distribution modulo one and Ratner s theorem. Equidistribution in number theory, an introduction, 27-244, ATO Sci. Ser. II Math. Phys. Chem., 237, Springer, Dordrecht, 27. [8] J. Marklof and A. Strömbergsson, Equidistribution of Kronecker sequences along closed horocycles. Geom. Funct. Anal. 3 (23) 239-28. [9] J. Marklof and A. Strömbergsson, The distribution of free path lengths in the periodic Lorentz gas and related lattice point problems. Ann. of Math. (2) 72 (2) 949-233. [] Z. Rudnick and P. Sarnak, The pair correlation function of fractional parts of polynomials, Comm. Math. Phys. 94 (998) 6-7.

4 JES MARKLOF AD ADREAS STRÖMBERGSSO [] Z. Rudnick, P. Sarnak and A. Zaharescu, The distribution of spacings between the fractional parts of n 2 α. Invent. Math. 45 (2) 37-57. [2] Z. Rudnick and A. Zaharescu, The distribution of spacings between fractional parts of lacunary sequences. Forum Math. 4 (22) 69-72. School of Mathematics, University of Bristol, Bristol BS8 TW, U.K. j.marklof@bristol.ac.uk Department of Mathematics, Box 48, Uppsala University, SE-756 Uppsala, Sweden astrombe@math.uu.se