arxiv:math/ v1 [math.nt] 19 Oct 2001

Transcription

1 arxiv:math/007v [math.nt] 9 Oct 00 THE STATISTICS OF THE TRAJECTORY IN A CERTAIN BILLIARD IN A FLAT TWO-TORUS FLORIN P. BOCA RADU N. GOLOGAN AND ALEXANDRU ZAHARESCU Abstract. We study a billiard in a modified suare with pockets of small size ε. The length respectively the number of reflections in the side cushions of the trajectory that originates at one corner under angle θ is denoted by l εθ respectively by R εθ. We prove that the probability measures associated with the random variables εl ε and εr ε are convergent as ε ց 0 and provide explicit formulas for the limits.. Introduction and main results For each ε > 0 we consider D ε = 0 [ε ε] {0} {0} [ε ε] that is the unit suare [0] with horizontal and vertical segments of length ε removed from each corner. The object of study of this paper is the billiard in D ε where the trajectory starts at the origin O is straigthforward and rebounds elastically in the side cushions until it reaches one of the removed segments cf. Figure. ε ε O ε θ ε ε ε Figure. The billiard in D ε For each θ 0 π lε θ denotes the length of the trajectory that starts under angle θ. It is not hard to see that l ε θ coincides actually with the length of the trajectory in a right isosceles triangle associated by Sinai [6 pp ] to a system of two points with eual masses m = m that move in the interval [0] starting from opposite end points with velocities v and v and rebounding elastically at the end points of the interval and when they collide. Received by the editors October Research partially supported by ANSTI grant C469/000.

2 F.P. BOCA R.N. GOLOGAN A. ZAHARESCU The moments of l ε as ε ց 0 were evaluated in [3] where it was proved that the r th power of the length of the trajectory is proportional to cos r x when averaging on every interval [β] [ 0 π ]. More precisely these authors have proved in [3] that for every r > 0 the estimate. ε r holds whenever δ > 0 with and D r is given by r + ln rr + l ε θ r dθ = c r β + O rδ ε /6 δ c r β = D r π r r + r+ r + + r k= dx cos r x k r k k k + r + k= k r+ k k k. In this paper we show that the existence of the main term in. arises as a result of the convergence of the probability measures µ βε on [0 defined by µ βε f = f εl ε θ dθ f C c [0 β as ε ց 0. The value of the repartition function F βε of µ βε at > 0 is β µ βε ] = χ β [0] εlε θ dθ = {θ [β]; εl εθ }. β Using techniues from [] [] and [3] which rely ultimately on estimates [5][4] of Weil type [7] for Kloosterman sums with non-prime modulus we prove Theorem.. Given 0 < β π 4 and 0 the limit F β = lim F βε exists. εց0 Moreover if we set ψt = t t + log t t

3 TRAJECTORY STATISTICS IN A CERTAIN BILLIARD 3 then β π F β euals β cos β + cos cos β / β cos β + cos β arccos t dt = sin β sin if [ 0 cos ] arccos t dt + cos / = sin β + arccos 4 + cos β β + ψtdt + β + β + / cos β / / ψtdt + ψtdt cos cos β cos β ψt arccos t dt cos / ψt arccos t dt if [ cos cos β] ψt arccos t dt if [ cos β cos ] ψt arccos t dt if [ cos cos β] if [ cos β. Lemma 3. in [3] implies that l ε θ [ ] ε ε. Hence all functions F βε and thus their limit F β are constant on the interval [ and this constant is eual to. As a result we gather t t t + log dt = π t. t With the change of variable t = u we arrive at Corollary.. log u π du = u + u log. The continuous function F β represents the repartition function of a probability measure f β d supported on [ 0 cos β]. That is Fβ = f β tdt. Besides we have 0

4 4 F.P. BOCA R.N. GOLOGAN A. ZAHARESCU Corollary.3. µ βε µ β weakly as ε ց 0 where dµ β d = f β with f β given by sinβ sin if [ 0 cos ] arccos 4 sinβ + ψcos xcos xdx if [ cos cos β] π β arccos ψcos xcos xdx if [ cos β cos ] ψcos xcos xdx if [ cos cos β] 0 if [ cos β Figure. The repartition function F 0π/4 and the density function f 0π/ Figure 3. The repartition function F 0π/6 and the density function f 0π/6 By the previous corollary µ βε t r µ β t r as ε ց 0 that is ε r. lim l ε θ r dθ = εց0 β 0 t r f β tdt = t r f β tdt. 0

5 TRAJECTORY STATISTICS IN A CERTAIN BILLIARD Figure 4. The repartition function F π/6π/4 and the density function f π/6π/4 Now we infer from. and. that.3 c r β = β 0 t r f β tdt. The number of reflections R ε θ has also been estimated in [3] where it was found that for all δ > 0 as ε ց 0 where ε R ε θdθ = cβ + O δ ε /6 δ cβ = 6 β + log cos cos β log π. As in the case of the length function we shall consider the probability measures ν βε on [0 defined by with repartition function ν βε f = f εr ε θ dθ β G βε = ν βε ] = {θ [β]; εr ε θ } β In this case we prove [0.

6 6 F.P. BOCA R.N. GOLOGAN A. ZAHARESCU Theorem.4. Given any 0 < β π 4 and 0 the limit G β = lim εց0 G βε exists. Moreover π β G β is eual to β + log sinβ + cosβ sin + cos = dx + tanx + arctan + β + sin log β + cosβ +tan + ψt arctan t dt β β / +tan β + β / +tan β + β / β + +tan ψtdt ψtdt / ψtdt / / +tan ψtdt + ψtdt + +tan β +tan β if [ 0 +tan ] + if [ +tan +tan β ] ψt arctan t dt if [ + tan + tan β] ψt arctan t dt if [ + tan β. if [ + tan + tan β] Corollary.5. ν βε ν β weakly as ε ց 0 where dν β d by dx + tanx = g β with β g β given if [ 0 +tan ] arctan dx + tanx + ψ +tan x if [ +tan + tanx +tan β ] arctan ψ +tan x + tanx dx if [ +tan β + tan ] ψ +tan x dx if [ + tan + tan β] + tanx arctan 0 if [ + tan β.

7 TRAJECTORY STATISTICS IN A CERTAIN BILLIARD Figure 5. The repartition function G 0π/4 and the density function g 0π/4. Formulas for sectors ending at Farey points For each integer let F denote the set of Farey fractions of order. It is well-known that if a < a are two succesive elements in F then a a = and + +. Conversely if {... } + + and a a = for some a {... } and a {... } then a a and are successive elements in F. We keep 0 < β π 4 and > 0 fixed throughout and set = [ ε]. The analysis in the proof of Lemma 3. in [3] and the triangle ineuality show that. S β ε ε β F βε S β + ε ε where we set. S β ε = a/ F [tan tan β] +a < with θ a denoting the angle determined by the trajectories which end near the point a. Let a < a < a be three consecutive fractions in F. The relation a a = gives = ā mod where ā denotes the multiplicative inverse of a mod. Since ] ā is the uniue multiple of in the interval ā ā]. Hence ā = [ ā] and therefore [ ] ā.3 = + ā. θ a Employing a a = and ] we arrive in a similar way at [ ] + ā.4 = ā. and Taking into account.3 and.4 we see that whence > ā > ā

8 8 F.P. BOCA R.N. GOLOGAN A. ZAHARESCU Lemma.. min > ā. < < ā max +. < < ā min. > max > ā >. To evaluate the sum S β ε we first need to estimate θ a in each of the four cases from Lemma. using the formula arctanx + h arctan x = for x [0] and small h > 0. We shall use h + x + Oh = ε = + O = + Oε and the ineualities + ε + ε > and + ε >. Case I. min >. h + x + h + Oh a+ ε " a"- ε a a - ε a- ε O Figure 6. Case I In this case a ε < a +ε < a < a ε < a+ε thus arctan a ε ε ε = + a ε + O ε ε = + a ε + O = ε + a + O ε = ε ε + a + O ε = + a + O. θ a = arctan a + ε Case II. < <. εa + a + ε

9 TRAJECTORY STATISTICS IN A CERTAIN BILLIARD 9 a+ ε " a"- ε a a + ε a- ε O In this case a ε < a +ε < a < a ε θ a = arctan a + ε arctan a + ε = = ε ε + a + ε + O Figure 7. Case II < a+ε a+ε. With h = ε + a+ε = ε + a we gather ε + O εa + O + a + ε. a +ε = ε Since > and + > we have >. Hence the error in the formula for θ a above is ε + ε ε3 + ε ε and we find θ a = ε ε + a + O = + a + O + ε = + ε + a + O. Since = ā mod and < we have = ā and so + ā ε θ a = ā + a + O. Case III. > >. " a"- ε a+ ε a a- ε a + ε " a"- ε a+ε a a- ε a +ε O O Figure 8. Case III

10 0 F.P. BOCA R.N. GOLOGAN A. ZAHARESCU In this case = ā. Moreover we have a ε < a +ε we may take h = a a ε = ε and gather θ a = arctan a ε = Case IV. > max. arctan a ε = ε. ā ā + a + O ε + a ε < a < a ε ε + O < a+ε. As a result a+ ε a+ε O " a"- ε a +ε a a- ε O " a"- ε a + ε a a- ε with In this case a ε < a +ε < a < a ε < a+ε θ a = arctan a ε Figure 9. Case IV θ a = arctan a ε arctan a = = + a + O and we gather arctan a + ε = θ a + θ a ε θ a = arctan a arctan a + ε = = + a + O ε. Since in this case = ā and = ā we arrive at θ a = = ε + a + O ε + a + O ā + a + ā + a + O + a ā ā + O ε +. ε ε ε +

11 TRAJECTORY STATISTICS IN A CERTAIN BILLIARD 3. The existence of the limiting repartition function To estimate the sum S β ε we shall employ Lemma 3.. [3 Lemma.] Let I J R be intervals of length lesser than let f : I J R be a C function and let T > and δ > 0. Then a I b J ab=mod fab = ϕ I J fxydxdy + O δ T /+δ f + T 3/+δ Df + I J Df T where = I J. If we put I = [ tan tan β] [0 ] then [ tan tan β] if [cos β] I = [ tan ] if [cos βcos ] if [cos. We also set I = [ ] min + I = [ ] [min + max ] I 3 = [ ] max = I f j xy = f j xy = F jy + x F y = y I F y = + y y where y I 3 F 3 y = y y = F y y I F 4 y = y y y I. Since ϕ = ε the analysis of cases I-IV and Lemma. provide S β ε = Sβ I ε + SII β ε + SIII β ε + SIV β ε + Oε

12 F.P. BOCA R.N. GOLOGAN A. ZAHARESCU where S I β ε = S II β ε = S III β ε = a/ F tan a/ tan β /3 +a < ā I a/ F tan a/ tan β +a < ā I 3 a/ F tan a/ tan β +a < ā I S IV β ε = a/ F tan a/ tan β >/3 +a < ā I + a = + ā ā + a = ā + a ā = + a ā ā = min/3 cos a I b I ab=mod min cos a I b I 3 ab=mod f ab f ab min cos a I b I ab=mod f 3 ab /3< min cos a I b I ab=mod f 4 ab. Lemma 3.. For all j { 3 4} i f j. ii Df j +. Proof. i In case I F y =. In case II y > thus 0 < +y y = F y <. In case III y <. Hence y < y and 0 < F 3 y = y y < so f 3. In case IV > y >. Hence y > and so < miny y. If y > then 0 F 4y = y y < y <. If y < then 0 F 4y < y < so F 4 <. Hence f j =. y = y ii In case II F < < + +. In case III y + and F 3 y = = y + since y >. In case IV F 4 y = y + y < + since y + and y +. Summarizing we get f j x = x + x F jy

13 TRAJECTORY STATISTICS IN A CERTAIN BILLIARD 3 and Hence f j y = + x F jy +. Df j + + = +. ED We shall next apply Lemma 3. to the functions f f f 3 f 4 I = I J one of the intervals I I I 3 and T = /6. For we first employ Lemma 3. to estimate the error term and note that I J T /+δ f j T T T = 3/+δ Df j T = = So if we set /+δ T /+δ = /6+δ = = = T 3/+δ + / /+δ + + T /< /+δ + T /+δ + T /+δ log /3+δ log + = T log /6+δ. β if [cos β] L = + x dx = arccos if [cos βcos ] I 0 if [cos + y H = I = dy H = dy y I I 3 y H 3 = y dy = H I H 4 = I dy = y y I dy y

14 4 F.P. BOCA R.N. GOLOGAN A. ZAHARESCU then S I β ε = S II β ε = S III β ε = S IV β ε = min/3 cos min cos min cos /3< min cos ϕ LH + O δ /6+δ ϕ LH + O δ /6+δ ϕ LH + O δ /6+δ ϕ LH 4 + O δ /6+δ. Since the contribution of = to each of the previous four sums is 0 we may assume next that <. If + then I = I = [ ] H = and I 3 =. If + < < + 3 then I = [ I = [ ] and I 3 = ]. Moreover H = + log and H = 3. If + 3 < then I = [ +] I = [ + ] and I 3 = ]. Moreover H = + log and H 4 = log +. Summarizing we find 3. S β ε = = ϕ V + O δ /6+δ whenever δ > 0 where V = LW and W is given by if min cos + + log if + < min cos log if + < min cos if > min cos. It is clear that L and W. Hence 3.3 V. Moreover W = log. Hence W and we gather 3.4 V d L W + W L d.

15 TRAJECTORY STATISTICS IN A CERTAIN BILLIARD 5 So we may apply [ Lemma.3] and employ and 3.4 to get 3.5 S β ε = ζ and Using now the ineualities 0 < 0 < V d + O δ ε /6 δ. log log d = log + d d together with and log log d d d S β ε = ζ d = log we arrive at 0 LW d + O δ ε /6 δ where W euals if 0 min/cos log if < mincos 0 if > mincos. The change of variable = t together witj 3.6 and 3.7 provide 3.8 S β ε = LtW tdt + O δ ε /6 δ ζ 0 where β if t [0cos β] 3.9 Lt = arccos t if t [cos βcos ] 0 if t [cos

16 6 F.P. BOCA R.N. GOLOGAN A. ZAHARESCU and if t [0mincos /] 3.0 W t = ψt if t [/mincos ] 0 if t [mincos with ψt as in Theorem.. Then and. prove Theorem.. To estimate G βε we follow literally the proof of Theorem. replacing only S β ε by S β ε = a/ F [tan tan β] +a< θ a thus I being replaced by Ĩ = [ tan tan β] [0 ]. Clearly We arrive at [ tan tan β] if [ +tan β] Ĩ = [ tan ] if [ +tan β +tan ] if > +tan. 3. G βε = Lt π W tdt + O δ ε /6 δ 0 where 3. β if t [ 0min +tan β ] Lt = arctan t if t [ min +tan ] 0 if t > min +tan and if t [ 0min +tan ] 3.3 W t = ψt if t [ min +tan ] 0 if t > +tan.

17 TRAJECTORY STATISTICS IN A CERTAIN BILLIARD 7 By and 3.3 we find that π β G βε is eual to β + β + tanβ + β + tanβ + +tan +tan β / +tan β β + β + +tan β / +tan β / / arctan arctan ψtdt + ψtdt + ψtdt t dt t dt + +tan +tan β +tan β ψt +tan / if [ 0 +tan ] ψt arctan t dt if [ +tan +tan β ] arctan t dt if [ +tan β + tan ] ψt arctan t dt if [ + tan β. if [ + tan + tan β] Theorem.4 follows now from the previous euality and some straightforward computations. References. F. P. Boca C. Cobeli and A. Zaharescu Distribution of lattice points visible from the origin Commun. Math. Phys F.P. Boca C. Cobeli and A. Zaharescu A conjecture of R.R. Hall on Farey points J. Reine Angew. Mathematik F.P. Boca R.N. Gologan and A. Zaharescu The average length of a trajectory in a certain billiard in a flat two-torus preprint arxiv math.nt/ T. Estermann On Kloosterman s sums Mathematika C. Hooley An asymptotic formula in the theory of numbers Proc. London Math. Soc Ya. G. Sinai Introduction to Ergodic Theory translated by V. Scheffer Mathematical Notes Princeton University Press Princeton NJ A. Weil On some exponential sums Proc. Nat. Acad. U.S.A FPB RNG and AZ: Institute of Mathematics of the Romanian Academy P.O. Box -764 Bucharest RO Romania FPB and AZ: Department of Mathematics University of Illinois at Urbana-Champaign Urbana IL 680 USA fboca@math.uiuc.edu rgologan@theta.ro zaharesc@math.uiuc.edu