arxiv: v1 [math.ca] 31 Mar 2014

Transcription

1 An asymptotic epansion for the generalised quadratic Gauss sum revisited arxiv: v1 [math.ca] 31 Mar 214 R. B. Paris School of Engineering, Computing and Applied Mathematics, University of Abertay Dundee, Dundee DD1 1HG, UK Abstract An asymptotic epansion for the generalised quadratic Gauss sum S N,θ) = N epij 2 ijθ), where, θ are real and N is a positive integer, is obtained as and N such that N is finite. The form of this epansion holds for all values of N θ and, in particular, in the neighbourhood of integer values of Nθ. A simple bound for the remainder in the epansion is derived. Numerical results are presented to demonstrate the accuracy of the epansion and the sharpness of the bound. Mathematics Subject Classification: 11L5, 41A3, 41A6, 41A8 Keywords: Quadratic Gauss sum, Asymptotics, Curlicues 1. Introduction We consider the asymptotic epansion of the generalised quadratic Gauss sum S N,θ) = N fj), ft) := epit 2 iθt), < < 1, 1 2 θ 1 2, 1.1) where N is a positive integer, as and N, such that the quantity N is finite. Applications of the above eponential sum arise in various number-theoretic contets and in the study of disorder in dynamical systems. The sum S N,θ) has a long history that goes back to Gauss, who evaluated the sum corresponding to = 2/N when θ =. The results of Gauss were generalised for rational = M/N, where M and N are relatively prime, into the well-known Cauchy-Kronecker formula [1] S N,) = ei/4 S M 1 ), = M/N, MN even). When its terms are regarded as unit vectors in the comple plane, the patterns produced by the partial sums of 1.1) for fied as N often result in a superposition of spirals or curlicues ) that can be highly intricate; see [2, 3, 4, 9]. The scalings of this hierarchy of spirals are found to depend delicately on the arithmetic nature of [3]. When = p/q, where p and q are relatively prime, and θ = the trace is relatively simple: when pq is even the spiral pattern is regular and diffuses in the comple plane away from the origin in blocks, whereas when pq 1

2 Asymptotics of the Gauss sum 2 is odd the pattern is periodic and repeats itself indefinitely as N. When is irrational a more complicated pattern emerges that seems to ehibit a random-walk behaviour; see [3, 12]. Estimates for the growth of S N,θ) when N is large and is fied in the range < < 1 are obtained by employing a renormalisation process based on the approimate functional relation [7] /i/4 S N,θ) = e iθ2 S N 1, θ ) ) 1 θ O. 1.2) This transformation shows that the sum S N,θ) over N terms can be approimated by a similar sum taken over N terms with the variable replaced by 1/ and θ by θ/. Repeated application of 1.2), making use of the simple symmetry properties satisfied by 1.1) to maintain in the interval < < 1 at each stage, enables the representation of S N,θ) in terms of a steadily decreasing number of terms. In this way it was shown in [7] that S N,θ) = on) for any irrational, with more precise order estimates depending on the detailed arithmetic structure of. The problem that concerns us here is the asymptotic estimation of S N,θ) for when N suchthatnisfinite. Anearlypaperdealingwith estimatesfors N,θ) when < < 1 is that of Fiedler et al. [6], and more recently that in [8, 2.2], but their error terms are too large for our purposes when. Following on from the gross estimates in [9], the leading terms in the epansion in the case θ = were obtained in [13] when N < 1. An epansion for S N,) valid as and finite N was obtained in [3] and [5, Theorem 4], although the remainder term was left as an order estimate. In this paper, we revisit the epansion of S N,θ) as and N such that N = O1) obtained in [12]. The sum S N,θ) is epressed eactly as a series of complementary error functions with argument proportional to 1/2, so that in the small- limit we may employ the well-known asymptotics of the complementary error function in the form [11, 7.12i)] where e z2 erfcz) = 1 ) rγr 1 2 ) Γ 1 2 ) z 2r 1 ˆT n z) z ), 1.3) r= ˆT n z) Γn 1 2 ) z 2 arg z 1 4) 1.4) and n is a positive integer; see also [1, p. 111]. In [12], the coefficients in the resulting epansion were epressed in terms of even-order derivatives of cotξ, where ξ = Nθ, which presented a complication when ξ passes through integer values. In addition, the remainder in the epansion was not epressed as a convenient bound. Here we remedy these deficiencies and give the epansion in a form with coefficients that do not present any difficulty in computation in the neighbourhood of integer values of ξ. In order to make the paper reasonably self-contained, we repeat in Section 2 the derivation of the representation of S N,θ) in terms of complementary error functions given in [12]. In Section 3, we establish the central result of the paper in the following theorem. Theorem 1 Let S N,θ) be the sum defined in 1.1), where < < 1, 1 2 θ 1 2 and N is a positive integer. Further, let ξ := Nθ, M = [ξ] be the nearest integer part of ξ and ǫ := ξ M, where 1 2 < ǫ 1 2. Then, as and N, such that N is finite, we have the epansion valid for M and n 1 /i/4 S N,θ) e iθ2 S M 1, θ ) 1 2 fn) 1) ei/4 {Eθ) fn)eǫ) = 1 Γr 1 2 ) ) rcr i Γ 1 r= 2 ) R n, 1.5) i

3 Asymptotics of the Gauss sum 3 where Et) := e it2 / erfce i/4 t /). The coefficients C r are given by and the remainder R n satisfies the bound C r = fn) r ǫ) r θ) r ) 1.6) R n 1 2 ) n where the quantities ± r λ) are defined in 3.1) and 3.2). ) n n ǫ) n θ)) n 1), 1.7) In Section 4, we present numerical results to demonstrate the accuracy of the above epansion and also the sharpness of the bound on the remainder term R n. 2. A representation for S N,θ) Let ξ := N θ, M = [ξ], ǫ = ξ [ξ], where [ξ] denotes the nearest integer part of ξ and 1 2 < ǫ 1 2. Define also the function Et) := e it2 / erfcωt /), ω = e i/4, E) = 1, E t) = 2e it2 / Et), 2.1) where erfc is the complementary error function. The reflection formula follows from the wellknown result erfcz) = 2 erfc z). From 1.3), we have the epansion for 1/2 t where from 1.4) Et) = 1 ) rγr 1 2 ) Γ 1 2 ) r= T n t) Γn 1 2 ) An application of Cauchy s theorem shows that N 1 fj) = 1 2i ) r 1 i 2 Tn t 2 t) n = 1,2,...), 2.2) C ) n1 2 t ) cott)ft) dt, where ft) is defined in 1.1) and C is a closed path encircling only the poles of the integrand at t = 1,2,...,N 1. We deform the path C into a parallelogram with two sides inclined at 1 4 to the real ais; see [12]. The vertices are situated at ±Pe i/4, N ±Pe i/4 P > ) and there are semi-circular indentations of radius δ < 1 around the points t = and t = N. Then, denoting the upper and lower halves of the contour by C 1 and C 2 respectively, we find following the discussion given in [1, p. 29] that N 1 fj) = N δ δ ft)dt C 1 ft) 1 e itdt C 2 ft) e it 1 dt. Now let P, so that the contributions from the parts of C 1 and C 2 parallel to the real ais vanish on account of the eponential decay of the factor epit 2 ), and let δ. The integrals around the indentation linking δe i/4 with δ and δ with δe i/4 then tend to 1 8 f) and 3 8 f), respectively; similarly for the indentation at t = N the integrals contribute 1 2 fn). Thus we obtain S N,θ) = N fj) = 1 2 fn) 1)J N e i/4 I N I ), 2.4)

4 Asymptotics of the Gauss sum 4 where the integral and we have defined with J N := N I j := ft)dt = ei/4 2 {Eθ) fn)eξ) 2.5) F j τ) e ωτ 1 dτ j =,N) F j τ) := fj τe i/4 ) fj τe i/4 ) = 2e τ2 fj)sinh{jθ)ωτ. It now remains to evaluate the integrals I and I N. If we epand the factor e ωτ 1) 1 as a finite geometric series together with a remainder we find, for positive integer K, I j = K e kωτ F j τ)dτ e Kωτ F j τ) e ωτ dτ. 2.6) 1 The first term on the left-hand side of this epression becomes upon insertion of the definition of F j τ) where 2fj) K e τ2 kωτ sinh{jθ)ωτdτ The remainder term in 2.6) is given by H K := = fj) e Kωτ F j τ) e ωτ dτ = fj) 1 K {Ek j θ) Ek jθ). 2.7) e ττω) e K j θ)ωτ Gτ)dτ, Gτ) := e jθ)ωτ sinh{jθ)ωτ/sinhωτ). Now G) = 2j θ) and Gτ) e ωτ as τ. It is also easy to see we omit these details) that Gτ) G) for τ [, ). Then, provided K > jθ j =,N) it follows that H K < e K j θ)ωrτ Gτ) dτ 2 1 2G) K j θ) where ω r = 1/ 2, with the result that H K as K. Therefore, from 2.7), we obtain I j = fj) {Ek j θ) Ek jθ) j =,N). 2.8) Substitution of 2.8) with j =,N into 2.4) then yields the desired representation of S N,θ) in terms of complementary error functions. From 2.2) we see that the terms in 2.8) are Ok 2 ) as k. We define the quantities for λ < 1 3. The epansion of S N,θ) as with N finite r λ) := ζ2r 1,1λ)ζ2r 1,1 λ) r 1), 3.1)

5 Asymptotics of the Gauss sum 5 and cotλ λ 1 r = ) r λ) := ζ2r1,1λ) ζ2r1,1 λ) r 1), where ζs,a) = k= k a) s Rs) > 1) is the Hurwitz zeta function. Note that r ) = for r and r ) = 2ζ2r1), where ζs) is the Riemann zeta function. When j =, we have from 2.8) I = 1 {Ek θ) Ek θ). In the limit, the arguments of the complementary error functions contained in Ek ± θ) have large modulus for k 1 and phase equal to 1 4, since 1 2 θ 1 2. Employing the epansion 2.2), we then obtain 3.2) where c r θ) = I = ei/4 k r=1 k θ) 2r 1 = Γr 1 2 ) ) rcr Γ 1 2 ) θ)r n θ) ), 3.3) i k 1θ) 2r 1 k 1 θ) 2r 1 k= k= = r θ). 3.4) In the case r = the sums must be interpreted in the principal value sense lim s s k= s a k to yield the evaluation c θ) = cotθ 1/θ. The remainder term R n θ) is given by R n θ) = 1 {T n k θ) T n k θ) and, from 2.3), therefore satisfies the bound R n θ) 1 2 ) n ) n k k θ 2 = 1 2 ) n where a) r = Γar)/Γa) is the Pochhammer symbol. Proceeding in a similar manner when j = N, we have I N = fn) = fn) {Ek ξ) Ek ξ) { M {2e ik ξ)2 / Eξ k) k=m1 Ek ξ) ) n n θ), 3.5) Ek ξ). Here we have made use of the reflection formula in 2.1) to separate off the error functions in Ek ξ) corresponding to k M when M 1). Upon noting that fn)e ik ξ)2 / = e iθ2 / e ik2 /kiθ/, we obtain when M 1 / I N = e iθ2 S M 1, θ ) { fn) 2 M Ek ξ) Eξ k) k=m1 Ek ξ). 3.6)

6 Asymptotics of the Gauss sum 6 If we now etract from the second sum in curly braces in 3.6) the error function Eξ k) corresponding to k = M that is, Eǫ)) and use the evaluation of the integral J N in 2.5), we can write e i/4 /i/4 I N J N = e iθ2 S M 1, θ ) ei/4 {Eθ) fn)eǫ) ei/4 fn) { k= M 1 Ek ξ) Eξ k) k=m1 Ek ξ), 3.7) where the term involving fn)eξ) from J N has been absorbed into the first sum in curly braces. Then in a similar manner to the determination of the epansion of I in 3.3) we find { 1 2 M 1 Ek ξ) Eξ k) Ek ξ) k= = ei/4 where, recalling that ξ = M ǫ, M = [ξ], c r ǫ) = k M r= k ξ) 2r 1 = and the remainder R n ǫ) satisfies the bound R n ǫ) 1 2 ) n = 1 2 ) n ) n k M k=m1 Γr 1 2 ) ) rcr Γ 1 2 ) ǫ)r n ǫ), 3.8) i k k ξ 2 = 1 2 ) n k ǫ) 2r 1 = r ǫ) 3.9) ) n k k ǫ 2 ) n n ǫ). 3.1) Combination of 3.7) and 3.8) then yields the epansion when M 1 e i/4 /i/4 I N J N = e iθ2 S M 1, θ ) ei/4 {Eθ) fn)eǫ) 1 Γr 1 2 ) ) rcr i Γ 1 r= 2 ) ǫ)r n ǫ). 3.11) i In the case M = when ξ = ǫ), the sum S M and from 3.6) we have { e i/4 I N J N = ei/4 {Eθ) fn)eǫ) ei/4 fn) 2 Ek ǫ) Ek ǫ). It is easily seen that we obtain the same epansion as 3.11). The form of the coefficients in 3.4) and 3.9) with r 1 presents no difficulty in computation in the neighbourhood of integer values of ξ where ǫ, in contrast to those given in [12] which involved even derivatives of cotξ. Although the coefficients c ǫ) and c θ) have a removable singularity at ǫ = and θ = their computation is straightforward. If we now define the coefficients C r and the remainder R n by C r := fn)c r ǫ) c r θ), R n := e i/4 {fn)r n ǫ) R n θ),

7 Asymptotics of the Gauss sum 7 then we see that C r = fn) r ǫ) r θ) r ) 3.12) and R n R n ǫ) R n θ) 1 2 ) n ) n{ n ǫ) nθ) n 1). 3.13) Combination of 2.4), 3.3) and 3.11), together with the above definitions of C r and the bound on R n, then gives the epansion of S N,θ) stated in Theorem 1. We remark that the terms Eθ) and Eǫ) have been left unepanded as in 3.11) and in Theorem 1, since for small values of θ and ǫ = o 1/2 ) these quantities can no longer be approimated by 2.2). 4. Numerical results and discussion In order to demonstrate the accuracy of the epansion in Theorem 1, we define the quantity S by /i/4 S := S N,θ) e iθ2 S M 1, θ ) 12 ei/4 fn) 1) 2 {Eθ) fn)eǫ). 4.1) Then from Theorem 1 we have the epansion as and N such that N is finite S = 1 Γr 1 2 ) ) rcr i Γ 1 r= 2 ) R n, 4.2) i where the coefficients C r are defined in 1.6) and the remainder R n satisfies the bound in 1.7). We remark that the bound in 1.7) is eplicitly independent of N. In Table 1, we show the absolute value of the error in the computation of S using the epansion 4.2) truncated after n terms for two different sets of values of, θ, summation inde N and different levels n. The eact value of S N,θ) was obtained by high-precision summation of 1.1). In Table 2, we compare the absolute values of the remainder R n calculated from 4.2) and its bound to illustrate the sharpness of 1.7). Table 1: Values of the absolute error in the computation of S by 4.2) for different truncation inde n. = 1/25 ) = 1/25 ) = 1/25 3) N = 73, θ =.125 N = 743, θ =.25 N = 6, θ = n ξ. = ξ. = ξ. = In the case of the classical quadratic Gauss sum θ = ), we have ξ = N = M ǫ with 1 2 < ǫ 1 2 when M 1). From 1.5) as, N such that N is finite, we obtain the epansion S N,) = ei/4 S M 1 ), 1 ei/4 fn) 1) 2 {1 fn)eǫ)

8 Asymptotics of the Gauss sum 8 Table 2: The absolute values of R n and the bound in 1.7) for different truncation inde n. = 1/25 ), θ =.125 = 1/25 ), θ =.25 N = 73, ξ. = N = 743, ξ. = n R n Bound R n Bound fn) i r= Γr 1 2 ) ) rcr Γ 1 2 ) ǫ)r i n, 4.3) where fn) = epin 2 ). From 1.7) and the fact that I when θ =, we have R n 1 2 ) n ) n n ǫ). 4.4) We emphasise that the epansion in 4.3) holds for all finite values of N; see Table 1. When M = [ξ] = and < ǫ 1 2 that is, when N < 1 2 the sum S M =. If, in addition, Eǫ) in 4.3) is epanded by means of 2.2) then we obtain an epansion equivalent to that in [5, Theorem 4], albeit with a bound for the remainder rather than an order estimate and coefficients c r epressed in a different form. However, the epansion of Eǫ) by 2.2) is only applicable when ǫ 1/2 ; that is, when N 1/2. Finally, we remark that since C r 1 ǫ ) 2r 1 for r 1 when ǫ is bounded away from zero), the optimal truncation inde r of the sum in 4.2) corresponding to truncation at, or near, the term of least magnitude) is given by r 1 ǫ ) 2 /. This shows that the values of the truncation inde n in Table 1 are highly sub-optimal and also gives an indication of the enormous accuracy that could be obtained from the epansion 4.2). References [1] T.M. APOSTOL, Introduction to Analytic Number Theory, Springer, New York, [2] E.A. COUTSIAS and N.D. KAZARINOFF, Disorder, renormalizability, theta functions and Cornu spirals, Physica, 26D 1987), [3] E.A. COUTSIAS and N.D. KAZARINOFF, The approimate functional formula for the theta function and Diophantine Gauss sums, Trans. Amer. Math. Soc., ), [4] F.M. DEKKING and M. MENDÈS-FRANCE, Uniform distribution modulo one: a geometrical viewpoint, J. Reine Angew. Math., ), [5] R. EVANS, M. MINEI and B. YEE, Incomplete higher-order Gauss sums, J. Math. Anal. Appl ), [6] H. FIEDLER, W. JURKAT and O. KÖRNER, Asymptotic epansions of finite theta series, Acta Arithmetica, ),

9 Asymptotics of the Gauss sum 9 [7] G.H. HARDY and J.E. LITTLEWOOD, Some problems of Diophantine approimation. II: The trigonometrical series associated with the elliptic ϑ-functions, Acta Math., ), Reprinted in Collected Papers of G.H. Hardy, Vol. 1, , Oford University Press, Oford. [8] E. KRÄTZEL, Analytische Funktionen in der Zahlentheorie, Teubner, Stuttgart-Leipzig, 2. [9] D.H. LEHMER, Incomplete Gauss sums, Mathematika, ), [1] F.W.J. OLVER, Asymptotics and Special Functions, Academic Press, New York, Reprinted A.K. Peters, Massachussets, [11] F.W.J. OLVER, D.W. LOZIER, R.F. BOISVERT and C.W. CLARK eds.) NIST Handbook of Mathematical Functions, Cambridge University Press, Cambridge, 21. [12] R.B. PARIS, An asymptotic epansion for the generalised quadratic Gauss sum, Appl. Math. Sci. 2 28), [13] F. SULLIVAN and U. ZANNIER, A note on incomplete quadratic Gauss sums, Boll. Un. Mat. Ital. A6 1992),