# Beurling integers with RH and large oscillation

Beurling integers with RH and large oscillation FREDERIK BROUCKE, GREGORY DEBRUYNE, AND JASSON VINDAS Abstract. We construct a Beurling generalized number system satisfying the Riemann hypothesis and whose integer counting function displays extremal oscillation in the follow- ing sense. The prime counting function of this number system satisﬁes π(x) = Li(x) + O( x), while its integer counting function satisﬁes the oscillation estimate N(x) = ρx + Ω x exp(−c log x log log x) for some c > 0, where ρ > 0 is its asymptotic density. The construction is inspired by a classical example of H. Bohr for optimality of the con- vexity bound for Dirichlet series, and combines saddle-point analysis with the Diamond- Montgomery-Vorhauer probabilistic method via random prime number system approxima- tions. 1. Introduction In [10], H. G. Diamond, H. L. Montgomery, and U. M. A. Vorhauer constructed a Beurling generalized number system with very regularly distributed integers, but whose distribution of prime numbers has large oscillation. In fact, given 1/2 < θ < 1, they showed the existence of Beurling numbers with integer counting function N satisfying the asymptotic estimate (1.1) N(x) = ρx + O(x ), for some ρ > 0, and whose primes realize the de la Vall´ee Poussin error term in the Prime Number Theorem (PNT), namely, with prime counting function satisfying the oscillation estimate π(x) = Li(x) + Ω x exp −c log x for some c > 0, where Li(x) stands for the logarithmic integral and the notation f(x) = ′ ′ Ω (g(x)) means that there is c > 0 such that the inequalities f(x) > c g(x) and f(x) < −c g(x) hold inﬁnitely often for arbitrary large values of x. Their Beurling number system has the additional feature that its associated zeta function, deﬁned as Z Z ∞ ∞ −s −s (1.2) ζ(s) = x dN(x) = exp x dΠ(x) 1 1 −1 1/j with Π(x) = j π(x ), also realizes the classical de la Vall´ee Poussin zero-free region; j=1 in particular, the Riemann hypothesis (RH) fails for it. 2020 Mathematics Subject Classiﬁcation. Primary 11M41, 11N80; Secondary 11M26, 11N05. Key words and phrases. Generalized integers with Riemann hypothesis; generalized integers with large oscillation; Beurling generalized prime numbers; saddle-point method; Bohr’s extremal example for convexity bound; Diamond-Montgomery-Vorhauer probabilistic method. F. Broucke was supported by the Ghent University BOF-grant 01J04017. G. Debruyne acknowledges support by Postdoctoral Research Fellowships of the Research Foundation– Flanders (grant number 3E006818) and the Belgian American Educational Foundation. The latter one allowed him to do part of this research at the University of Illinois at Urbana-Champaign. J. Vindas was partly supported by Ghent University through the BOF-grant 01J04017 and by the Research Foundation–Flanders through the fwo-grant 1510119N. arXiv:2004.11501v2 [math.NT] 26 May 2020 2 F. BROUCKE, G. DEBRUYNE, AND J. VINDAS The existence of such a number system proves the optimality of E. Landau’s classical PNT [15], which, recasted in the language of Beurling numbers, states that (1.1) always implies the validity of the PNT in the form (1.3) π(x) = Li(x) + O x exp −c log x , for some c > 0. That Landau’s theorem was optimal is part of a long-standing open problem, posed by P. T. Bateman and Diamond in [4, 13B, p. 199]. The number system constructed by Diamond, Montgomery, and Vorhauer also provides the valuable information that a zeta function might not have a wider zero-free region than that of de la Vall´ee Poussin if we only require that Beurling’s property that the integers have multiplicative structure and (1.1) hold. Several arguments from [10] have been sharpened by W.-B. Zhang in [19]. Interestingly, Zhang complemented these results by showing that there are also Beurling number systems for which, in contrast, the RH and the asymptotic estimate (1.1) both hold. In this work we shall establish the existence of a counterpart of the result from [10] for number systems whose integers display large oscillation in the sense of the following theorem, which is in fact our main result. The present article is fully devoted to its proof. Let us ﬁrst recall that, following A. Beurling [5] (see also [4, 11]), a set of generalized primes is simply an unbounded sequence of real numbers P : p ≤ p ≤ . . . subject to the only 1 2 requirement p > 1. Its associated set of generalized integers is the multiplicative semigroup generated by P and 1 (where multiplicities according to diﬀerent representations as products of generalized primes are taken into account). The symbols N(x) and π(x) already used above are the functions that count the number of generalized integers and primes, respectively, not exceeding a given number x. Theorem 1.1. There exists a Beurling generalized number system such that (1.4) π(x) = Li(x) + O( x) and, for any constant c > 2 2, (1.5) N(x) = ρx + Ω x exp(−c log x log log x) , where ρ > 0 is the asymptotic density of N. Theorem 1.1 proves, apart from the numerical value of the constant c, the optimality of the ensuing result due to T. W. Hilberdink and M. L. Lapidus [14] (see also [13]), which, in turn, might be regarded as the analog of Landau’s result in the reverse direction, that is, for number systems with very regular prime counting function (namely, satisfying (1.6)). Theorem 1.2 ([14]). Suppose the generalized prime counting function satisﬁes (1.6) π(x) = Li(x) + O(x ) for some 0 < θ < 1. Then, there are constants ρ > 0 and c > 0 such that its associated generalized integer counting function satisﬁes (1.7) N(x) = ρx + O x exp −c log x log log x . Theorem 1.2 is an improvement to a result of P. Malliavin [17], who obtained the weaker error term O x exp −c log x in (1.7) also under the hypothesis (1.6). Hilberdink and √ √ Lapidus have actually proved that (1.7) holds with the constant c = min{ 1 − θ/4, 2/8}. We mention that, using a variant of M. Balazard’s technique from [3], one might improve BEURLING INTEGERS WITH RH AND LARGE OSCILLATION 3 this value to any c < min{ 2(1 − θ), 1}. We include details about how to obtain this improvement in Appendix A. Note that (1.4) implies the Beurling number system satisﬁes the RH, that is, its zeta function analytically extends to Re s > 1/2, except for a simple pole located at s = 1, and has no zeros in this half-plane. In this regard, it is worthwhile to compare our generalized number system from Theorem 1.1 with earlier examples by E. Balanzario [2] and Zhang [19]. On the one hand, in Balanzario’s example the oscillation estimate (1.5) holds for N, but π only satisﬁes the weaker asymptotic relation (1.3). On the other hand, Zhang’s example quoted above has generalized prime counting function satisfying (1.4) (hence the RH holds here), but its generalized integer counting function is too regular for our purposes, i.e., (1.1) holds for it. Our Beurling number system enjoys the most extremal features of Balanzario’s and Zhang’s instances and, in turn, neither of their constructions is able to simultaneously deliver (1.4) and (1.5). Our method for showing Theorem 1.1 is ﬁrst to construct a continuous analog of a number system having the desired properties and then to ﬁnd a suitable discrete approximation to it, yielding the sought discrete Beurling number system. As in classical number theory, the key property linking N and Π for Beurling numbers is the zeta function identity (1.2), or equivalently, the relation dN = exp (dΠ), where the exponential is taken with respect to the multiplicative convolution of measures [9, 11]. The latter exponential identity then makes sense for not necessarily atomic measures dΠ and dN (hereafter supported in [1,∞) and non-negative), giving rise to the possibility to consider ‘continuous number systems’. The idea of using such continuous analogs to show optimality of results of course goes back to Beurling’s seminal paper [5], and has extensively been exploited by several authors since then . We shall show the next result. Theorem 1.3. There exists an absolutely continuous prime distribution function Π with associated integer distribution function N (determined by dN = exp (dΠ )) such that C C C x −1 1 − u (1.8) Π (x) = du + O(1) log u and, for any constant c > 2 2, (1.9) N (x) = ρx + Ω x exp −c log x log log x , C ± where ρ > 0 is the asymptotic density of N . It is worth pointing out that Theorem 1.3 explains the fundamental diﬀerence between Zhang’s work [19] and ours. Zhang’s example of a Beurling number system satisfying the RH can basically be considered as a discretization of the absolutely continuous prime distribution −1 P (x) = (1 − u )/ log u du, which is also the starting point of numerous constructions in the theory of Beurling primes due to fact that it leads to a highly regular continuous number system [11]. Indeed, the exponential of dP is exp (dP ) = δ + dx, with δ the 1 1 Dirac delta concentrated at 1 and dx the Lebesgue measure, so that its associated integer distribution function is simply the function x on [1,∞), while its associated zeta function is the meromorphic function s/(s − 1). In order to demonstrate Theorem 1.3, one needs to Balanzario constructs a ‘continuous’ example in [2], but F. A. Al-Maamori has recently shown [1] via probabilistic arguments that the example of Balanzario can be discretized. A large number of instances has been collected in the monograph [11]; see also the recent work [8]. 4 F. BROUCKE, G. DEBRUYNE, AND J. VINDAS show the existence of a perturbation of dP whose exponential displays the extremal behavior (1.9). It is quite remarkable that such a small perturbation in the primes can reinforce itself in such a way, analogous to the butterﬂy eﬀect, to produce such a big discrepancy in the integers. Nevertheless, it cannot fully destroy the integer law. The construction and analysis of a continuous example possessing the properties stated in Theorem 1.3 that we give is quite involved, and is the subject of Sections 2 through 5. We refer to Section 2 for a sketch of the proof of Theorem 1.3, and some insights into the motivation for our considerations. We conﬁne ourselves here to mention that the starting template for our example is an old construction from H. Bohr’s thesis [6], which we shall manipulate to achieve the desired properties. The estimate (1.8) will automatically be satisﬁed by construction; the challeng- ing part is to match it with the oscillation estimate (1.9), which we shall actually deduce from a certain extremal behavior of the associated zeta function that we will generate in our construction. As a matter of fact, most of our work in the subsequent sections is a detailed saddle-point analysis of this zeta function. After establishing that our continuous example satisﬁes all requirements from Theorem 1.3, we proceed to carry out a discretization proce- dure for it in Section 6. This procedure will be accomplished by adapting to our problem the Diamond-Montgomery-Vorhauer probabilistic scheme [10] based upon approximations by random Beurling primes. Our adaptation of this scheme delivers suﬃciently strong bounds for the modulus of the relevant zeta functions; however, for our application, we also need to keep good control on the argument of the randomly found zeta function, for which the direct bounds from the Diamond-Montgomery-Vorhauer method appear to be insuﬃcient. We will resolve this issue with a new idea of adding ﬁnitely many well-chosen primes to our number system. We end this introduction by placing Theorem 1.1 in the context of a long-standing open problem. In fact, the Beurling number system we exhibit in this article and the Diamond- Montgomery-Vorhauer example from [10] are just two pieces of a fascinating unsolved puzzle essentially raised by Bateman and Diamond in [4, 13B, p. 199] and having its roots in the work of Malliavin. In [17], Malliavin discovered that the two asymptotic relations (P ) π(x) = Li(x) + O(x exp(−c log x)) and ′ β (N ) N(x) = ρx + O(x exp(−c log x)) (ρ > 0), for some c > 0 and c > 0, are closely related to each other in the sense that if (N ) holds for a given 0 < β ≤ 1, then (P ∗) is satisﬁed for a certain α , and vice versa the relation (P ) α α ∗ ∗ ∗ for a given 0 < α ≤ 1 implies that (N ∗) holds for a certain β . Writing α (β) and β (α) for the best possible exponents in these implications, we have: ∗ ∗ Problem 1.4. Given any α, β ∈ (0, 1], ﬁnd the best exponents α (β) and β (α). Theorem 1.1 and Theorem 1.2 together then yield β (1) = 1/2, while the work of Diamond, Montgomery, and Vorhauer in combination with Landau’s PNT gives α (1) = 1/2. These are the only two cases where a solution to Problem 1.4 is known, and for the remaining values 0 < α < 1 and 0 < β < 1 the question remains wide open. It has been conjectured by Bateman and Diamond that α (β) = β/(1 + β). However, the best known admissible value 3 ∗ ∗ To be precise, the suprema over all admissible values α and β in these implications, respectively. BEURLING INTEGERS WITH RH AND LARGE OSCILLATION 5 [11, Theorem 16.8, p. 187] when 0 < β < 1 is α ≈ β/(6.91 + β); this falls far short of the conjectural exponent. It is also believed that β (α) = α/(α + 1) for each 0 < α < 1, which is suggested by the work of Diamond, who showed in [9] that the hypothesis (P ) actually ensures a slightly better asymptotic estimate than (N ), namely, α/(α+1) α+1 N(x) = ρx + O(x exp(−c (log x log log x) ) for some ρ, c > 0. Al-Maamori [1] has recently found an upper bound, providing α/(α+1) ≤ β (α) ≤ α when 0 < α < 1, but we strongly believe that there is still room for improvement here. In fact, it is worth noting that Theorem 1.1 leads to β (α) ≤ 1/2, which is better than Al-Maamori’s upper bound in the range 1/2 < α < 1. The authors thank Harold G. Diamond for his useful comments and remarks. 2. Setup and overview of construction for the continuous example Before introducing our continuous analog of a Beurling number system that will satisfy the properties stated in Theorem 1.3, let us ﬁrst sketch the motivation for its deﬁnition. In the complex analysis proof of Theorem 1.2 from [14], the error term in (1.7) comes from an integral of the zeta function over the contour given by log log|t| s(t) = 1 − + it, |t| ≥ 3. log|t| In order to generate an example for which this error term is reached, one might thus attempt to ﬁnd a zeta function having certain extremal growth properties along this contour. Taking a closer look at the proof of Theorem 1.2, one sees that the bound obtained in [14] for log ζ is essentially the convexity bound (cf. [18, Theorem 1.19 and Theorem 1.20, pp. 201– 202]) for the Dirichlet series of log ζ. The core of our idea is to construct a prime counting function whose Mellin-Stieltjes transform attains this convexity bound. The inspiration for our construction goes back to Bohr, who showed in his thesis [6] via an ingenious example that the convexity bound for Dirichlet series is basically optimal. The example we shall now study is in fact a subtle variant of Bohr’s example, modiﬁed in such a way that it indeed gives rise to an absolutely continuous prime distribution function having the desired properties to deliver a proof of Theorem 1.3. Let us set up our construction. We consider a positive sequence (τ ) with τ ≥ 3 and k k 0 rapidly increasing to ∞, a positive sequence (δ ) tending to 0, and a sequence (ν ) that k k k k takes values between 2 and 3. For x > 1, we then deﬁne x 1+δ ν k k 1 − 1/u sin(τ log x) for τ < x ≤ τ , k k Π (x) := du+ R (x), with R (x) := C k k log u 0 otherwise. k=0 We make a choice for the sequences (τ ) , (δ ) , (ν ) such that the following (technical) k k k k k k properties (whose relevance will become clear in later stages of our analysis) hold. First, we set log log τ + a k k δ = , log τ 4 ′ In [14], they actually prove the convexity bound for −ζ /ζ, the Mellin-Stieltjes transform of ψ(x) = log u dΠ(u), and then derive a bound for log ζ on the given contour via integration. One can however start from Π instead of ψ and directly prove the convexity bound for log ζ. 6 F. BROUCKE, G. DEBRUYNE, AND J. VINDAS for a sequence (a ) taking values between log 6 and log 6 + 1 say, and a sequence (τ ) to k k k k be deﬁned below, and we deﬁne a sequence (x ) via k k log x log log x k k (2.1) log τ = . The deﬁnition of the sequence (x ) is of course reminiscent of the error term in (1.7), and k k it is in fact on this sequence that the average N (u) du will display a desired deviation from the main term ρx /2 as explained below (see (2.3)). We also mention that τ will be assumed to be suﬃciently large as needed in some of our future arguments. Then, we require the following properties: (a) τ > (2τ ) ; k+1 k (b) (1 + δ )τ log τ ∈ 2πZ and ν τ log τ ∈ 2πZ; k k k k k k (c) τ log x ∈ 2πZ if k is even while τ log x ∈ π + 2πZ when k is odd; k k k k (d) √ p log x 1 + 2 log log x / log x 1/32 k k k d 1 − ,Z < . 3/4 (1 + δ ) log τ (1 + δ ) log τ (log τ ) k k k k k Here d( · ,Z) denotes the distance to the nearest integer. The existence of such sequences is stated in the ensuing lemma, whose proof will be postponed to Section 5. Lemma 2.1. There exist sequences (τ ) , (a ) , (ν ) such that, with the above deﬁnitions of k k k k k k (δ ) and (x ) , the properties (a)-(d) are satisﬁed. k k k k It is obvious that Π satisﬁes (1.8). The function Π is indeed an absolutely continuous C C prime distribution function: Lemma 2.2. The function Π is absolutely continuous and non-decreasing. Proof. That Π is absolutely continuous is a simple consequence of its deﬁnition and Property 1+δ ν k k (b). For τ < x < τ , k k 1 − 1/x τ cos(τ log x) 1 k k −δ Π (x) = + ≥ − τ ≥ 0, C k log x x 2ν log τ k k by the deﬁnition of δ and since a ≥ log(2ν ). Hence Π is non-decreasing. k k k C As in the statement of Theorem 1.3, we deﬁne the associated continuous integer distribu- tion function N via dN = exp (dΠ ) and set C C C Z Z ∞ ∞ −s −s ζ (s) := x dN (x) = exp x dΠ (x) . C C C 1 1 Then, we have using (b), 1−(1+δ )s 1−ν s 1−(1+δ )s 1−ν s ∞ k k k k s 1 τ − τ τ − τ −s k k k k log ζ (s) = x dΠ (x) = log + + . C C s − 1 2 s − iτ s + iτ k k k=0 This Mellin-Stieltjes transform is absolutely convergent for σ > 1, and from the above formula, we immediately see that log ζ has an analytic continuation to any simply connected region contained in σ > 0 which does not contain 1, and furthermore ζ has a meromorphic continuation to σ > 0 with a single simple pole at s = 1. We remark that log ζ indeed C BEURLING INTEGERS WITH RH AND LARGE OSCILLATION 7 reaches the convexity bound (cf. [18, Theorem 1.19, p. 201]): if t = τ , then log ζ (σ + it) ≫ k C 1−(1+δ )σ t . To go from the zeta function back to N , one uses Perron’s inversion formula. For the sake of technical simplicity, we will employ Perron’s formula for the primitive of N , so that the integral converges absolutely. For κ > 1, we have Z Z Z x κ+i∞ κ+i∞ s+1 s+1 1 x 1 x N (u) du = ζ (s) ds = exp . . . ds. C C 2πi s(s + 1) 2πi (s − 1)(s + 1) 1 κ−i∞ κ−i∞ We shift the contour of integration to the left, more speciﬁcally to the contour considered by Hilberdink and Lapidus in their proof of Theorem 1.2. Set e e Γ := { 1 − 1/e + it : 0 < t < e } ∪ { 1 − log log t/ log t + it : t > e }. H-L Set ρ := Res ζ (s). By the residue theorem, we get a contribution from the pole at s = 1: s=1 C Z Z s+1 ρ 1 x (2.2) N (u) du = x + Im exp . . . ds, 2 π (s − 1)(s + 1) 1 Γ H-L where we have used ζ (s) = ζ (s) to restrict the path of integration to the upper half plane. C C 1−(1+δ )s In order to estimate the remaining integral, we will exploit the fact that τ /(s − iτ ) becomes relatively small when t is far from τ . Speciﬁcally, for each k we will choose a suitable x (namely x = x deﬁned above in (2.1)) so that the integral near t = τ will give k k a contribution of order x exp(−c log x log log x), and so that the rest of the integral is of lower order. Making this explicit is however a technically challenging problem. One of the diﬃculties that arises comes from taking the exponential: lower bounds on |log ζ | do not necessarily imply lower bounds on |ζ |, and furthermore, exponentiation also introduces a lot of oscillation. In order to extract the contribution of the integral near t = τ , we will use the saddle-point method. Let us brieﬂy review some ideas connected with the saddle-point method. Given a region Ω ⊆ C, a contour Γ ⊆ Ω, and analytic functions f and g on Ω, one might proceed as follows f(z) to estimate the integral g(z)e dz. First, one computes the saddle points of f; these are the points s ∈ Ω for which f (s ) = 0. Near a saddle point, the graph of Re f looks like j j a saddle surface. The idea is to shift the contour Γ while ﬁxing the endpoints to a contour Γ which passes through the saddle points in such a way that on the new contour Ref(s) reaches a maximum at the saddle points – whether this is possible of course depends on the speciﬁc situation. Approximating f by its second order Taylor polynomial near saddle points, one gets the following approximation: Z Z X X 2π ′′ 2 f(z) f(s ) f (s )(s−s ) /2 f(s ) j j j j g(z)e dz ≈ g(s )e e ds ≈ g(s )e . j j ′′ −f (s ) Γ j j j It might still be of interest to try to carry out the computation with the Perron formula for N instead of the one for its primitive, since it appears that on the sequence x˜ deﬁned via log τ = log x˜ log log x˜ , k k k k the contribution of the saddle points (see Section 3) is at least x˜ exp(−2 log x˜ log log x˜ + lower order ). k k k If the remainder of the Perron integral could be adequately estimated, this would improve the range for the constant c in Theorem 1.3 to any c > 2. 8 F. BROUCKE, G. DEBRUYNE, AND J. VINDAS λf(z) Often, one wants to estimate g(z)e dz for a parameter λ tending to ∞, and under certain assumptions, one can deduce that 2π λf(z) λf(s ) g(z)e dz ∼ g(s )e , as λ → ∞. ′′ −λf (s ) We refer to [7, Chapters 5 and 6] for a classical account of the method, and to [12, Section 3.6] for a distributional approach to this technique. In our case, we will apply the saddle-point method to a portion of the integral near t = τ with exp · · · 1−(1+δ )s 1 τ f(s) = f (s) := (s + 1) log x + , g(s) = g (s) := , k k k 2 s − iτ (s − 1)(s + 1) 1−(1+δ )s where we use the notation to indicate that we exclude the term τ /(2(s − iτ )) from the summation. In the absence of a parameter λ as in the standard situation mentioned above, we are led to make a rather explicit and detailed analysis of the integral term appearing in the right-hand side of (2.2). Summarizing, we will get N (u) du = x + contribution from saddle points + remainder . In Section 3, we will deal with the saddle points and show that their contribution is ≫ x exp(−c log x log log x). In Section 4, we will show that the remaining part of the Perron integral is of strictly lower order. Summing up all results, we will prove that 2 2 (2.3) N (u) du = x + Ω x exp −c log x log log x C ± for any c > 2 2. The relation (2.3) readily implies (1.9), so this will ﬁnally establish Theorem 1.3. From now on, we ﬁx a speciﬁc k, and investigate the Perron integral for x = x given by (2.1). For ease of notation, we will drop the index k everywhere, unless we need to make the distinction between the speciﬁc τ and the other τ , j 6= k. k j 3. The contribution from the saddle points 3.1. The saddle points. Recall that we have set 1−(1+δ)s 1 τ log x log log x f(s) = (s + 1) log x + , log τ = . 2 s − iτ 2 Also, log log τ + a δ = . log τ We have 1−(1+δ)s 1 τ 1 (3.1) f (s) = log x − (1 + δ) log τ + . 2 s − iτ s − iτ BEURLING INTEGERS WITH RH AND LARGE OSCILLATION 9 We will show that f has a saddle point on the line t = τ. Due to the periodicity of 1−(1+δ)s τ , there will also be saddle points near t = τ + 2πm/((1 + δ) log τ), when m ∈ Z is not too large. Making this explicit, set 2πm + π/2 2πm − π/2 + − : : (3.2) t = τ + , t = τ + , m m (1 + δ) log τ (1 + δ) log τ ± ± and set V to be the rectangle with vertices 1/2 + it , 1 + it , m ∈ Z. m m Lemma 3.1. Suppose |m| < ε log τ for some suﬃciently small (ﬁxed) ε > 0. Then f has a unique saddle point (of multiplicity 1) in the interior of V . Proof. Starting with the vertical edge on the right, and continuing in a counter clockwise fashion, we name the edges of ∂V as E , j = 1, . . ., 4. By the assumption on m, arg(s−iτ) < m j ′ ′ ε for some small ε when s ∈ ∂V . On the ﬁrst segment E , the ﬁrst term in (3.1) is dominant m 1 and so the argument of f is close to 0. On E , the second term becomes more and more signiﬁcant, and the argument increases from about 0 to about π/2. On E the argument increases further from about π/2 to about 3π/2. Finally, on E , the argument increases further from about 3π/2 to about 2π, as the ﬁrst term becomes once again dominant. Hence the winding number of the curve f (∂V ) around the origin is 1, and the lemma follows from the argument principle. 6 3/4 Suppose now that|m| < (log τ) . By the above lemma, for every such m, f has a unique saddle point s in the rectangle V , which is of multiplicity 1. We have f (s ) = 0, which m m m is equivalent to 1−(1+δ)s 1 τ log x (3.3) = . 2 s − iτ (1 + δ) log τ + s −iτ By taking logarithms, one sees that for each m, there is an integer n such that the following implicit equations for the real and imaginary part of s = σ + it hold: m m m (3.4a) 1 1 1 σ = 1 − log log x + log 2 + log|s − iτ| − log (1 + δ) log τ + , m m 1 + δ log τ s − iτ 1 1 (3.4b) t = arg (1 + δ) log τ + − arg(s − iτ) + 2πn . m m (1 + δ) log τ s − iτ 3/4 Let us ﬁrst look at the equation (3.4a) for σ . Since s ∈ V and |m| < (log τ) , m m m −1 log|s − iτ| ≪ 1 and log (1 + δ) log τ + (s − iτ) = log log τ + O(1). Inserting this in m m (3.4a) and using a Taylor approximation for 1/(1 + δ), we see that log log τ + a log log τ 1 σ = 1 − + O 1 − log log x − log log τ + O(1) log τ log τ log τ log log x 1 = 1 − 2 + O √ . log x log x log log x The choice for the exponent 3/4 is a bit arbitrary. It allows us to achieve error terms of decent quality with relatively modest eﬀort (compared to the bound ε log τ). Later on, we will impose a stronger restriction on m, which will appear naturally. We have chosen not to impose this restriction here already and to begin with this rather arbitrary one, because we think the intrinsic nature of the additional restriction might get clouded otherwise. 10 F. BROUCKE, G. DEBRUYNE, AND J. VINDAS 3/4 Using this approximation for σ (and again |m| < (log τ) ), one sees that log|s − iτ| ≪ m m −1/8 −1 1/2 (log x log log x) and log (1 + δ) log τ + (s − iτ) = log log τ + O((log log x/ log x) ), so that log log x 2(a + log 2) 1 (3.5) σ = 1 − 2 − √ . + O 5/8 log x log x log log x (log x log log x) Also, repeating the argument for σ we obtain log log x 2(a + log 2) log log x (3.6) σ = 1 − 2 √ − + O . log x log x log log x log x Let us now look at equation (3.4b) for t . By assumption (b), (1 + δ)τ log τ = 2πM for some integer M. Hence, we see that t = τ satisﬁes the equation with n = M. To see 0 0 that indeed Im s = τ, one can check using a continuity argument that equation (3.4a) for m = 0 and τ = t has a solution σ between 1/2 and 1. The point σ + iτ then satisﬁes both 0 0 0 (3.4a) and (3.4b), and by uniqueness (Lemma 3.1), we must have s = σ + iτ. For general 0 0 m, diﬀerent from 0, we will again approximate the solutions. The arguments appearing in (3.4b) can be written as follows: t − τ α := arg(σ + i(t − τ)) = arctan , m m m t −τ 1 2 2 σ +(t −τ) β = arg (1 + δ) log τ + = arctan . σ + i(t − τ) (1 + δ) log τ + m m 2 2 σ +(t −τ) 3/4 Using the bound arctan x ≪|x|, and the fact that t ∈ V (|m| < (log τ) ), we ﬁnd that m m n = M + m, and we get a ﬁrst approximation for t : m m 2πm 1 t = τ + + O . 5/4 (1 + δ) log τ (log τ) Employing the above approximation and (3.5), we can get ﬁner estimates for α and β : m m 1 2πm 1 |m| α = + O + O 5/4 3 σ (1 + δ) log τ (log τ) (log τ) 2πm |m| log log x 1 |m| = + O + + , 5/4 3 (1 + δ) log τ log τ log x (log τ) (log τ) β = O . 5/4 (log τ) They, in turn, yield a better asymptotic estimate for t : 2πm 1 1 |m| log log τ |m| t = τ + 1 − + O + + . 9/4 3 4 (1 + δ) log τ (1 + δ) log τ (log τ) (log τ) (log τ) BEURLING INTEGERS WITH RH AND LARGE OSCILLATION 11 Repeating the procedure one ﬁnal time, we obtain 1 2πm 1 1 |m| log log τ |m| α = 1 − + O + + 9/4 3 3 σ (1 + δ) log τ (1 + δ) log τ (log τ) (log τ) (log τ) 2πm log log x |m| |m| = 1 + 2 + O + , 2 3 (1 + δ) log τ log x (log τ) (log τ) 2πm |m| log log τ |m| β = − + O + , 2 2 3 4 (1 + δ) (log τ) (log τ) (log τ) so that log log x 1 + 2 2πm log x |m| |m| (3.7) t = τ + 1 − + O + . 3 4 (1 + δ) log τ (1 + δ) log τ (log τ) (log τ) For future computations, it is useful to have an approximation for f and f near the saddle points. Lemma 3.2. There are continuous functions λ (s) and λ (s) such that m m ′′ f (s ) f(s) = f(s ) + (s − s ) (1 + λ (s)), m m m ′ ′′ f (s) = f (s )(s − s )(1 + λ (s)), m m m and with the property that for each ε > 0 there exists an η > 0 independent of τ and m such that |s − s | < =⇒ λ (s) + λ (s) < ε. m m m log τ Proof. We will show the assertion for f, the proof of the statement concerning f is similar. As f is analytic in a neighborhood of s and f (s ) = 0, we have, for s suﬃciently close to m m s , ′′ 3 f (s ) M |s − s | m f,s,m m f(s) − f(s ) − (s − s ) ≤ , m m 2 6 (3) where M is the maximum of |f | on the line between s and s. The derivatives of f f,s,m m for n ≥ 2 are 1 n (−1) l! (n) 1−(1+δ)s n−l n−l f (s) = τ (−1) ((1 + δ) log τ) . l+1 2 l (s − iτ) l=0 Therefore, since s − iτ ∼ 1, M |s − s | f,s,m m |λ (s)| ≤ ≪ |s − s | log τ, m m ′′ 3|f (s )| where the implicit constant is independent of m and τ. The continuity of λ is obvious. m 12 F. BROUCKE, G. DEBRUYNE, AND J. VINDAS 3.2. The paths of steepest descent. In order to estimate the contribution of the saddle points, we shall shift, near each saddle point, the contour to the so-called path of steepest descent. This is a contour through the saddle point which, when starting at the saddle point, displays the biggest decrease in Re f(s) among all possible paths. Intuitively, this path connects the two “valleys” on both sides of the saddle point in the most economical way. Starting in one of the “valleys”, the tangent vector along this path is at ﬁrst a positive multiple of ∇ Ref(s), as Re f(s) increases to a maximum at the saddle point. After passing the saddle point, the tangent vector along this path is a positive multiple of −∇ Re f(s), as Re f(s) decreases. Using the Cauchy-Riemann equations, one sees that Im f(s) is constant along this path. It is worth mentioning that there is another path through the saddle point on which Im f(s) is constant, namely the path of steepest ascent (which displays the opposite behavior of the descent path). We will show that for each m, there is a path of steepest descent which goes from the bottom horizontal edge of V to its upper horizontal edge. The situations for m = 0 and m 6= 0 are a bit diﬀerent; let us ﬁrst describe it for m = 0. Write θ = (t − τ)(1 + δ) log τ, so that θ varies between −π/2 and π/2 as t varies between − + t and t ; see (3.2). The equation Im f(s) = Im f(s ) is equivalent to 0 0 1−(1+δ)σ 1 τ (3.8) (t − τ) cos θ + σ sin θ = (t − τ) log x. 2 2 2 σ + (t − τ) Trivial solutions are given by the line s = σ+iτ. It is however readily seen that this line is the path of steepest ascent, by examining the behavior of Re f. Consider now t 6= τ ﬁxed (hence also θ 6= 0 ﬁxed). Then, the equation (3.8) has a unique solution for σ in the range 1/2 < σ < 1: ﬁrst of all sgn(LHS) = sgn(RHS), and second we have that |LHS| is monotonically (1−δ)/2 decreasing in that range, with LHS(σ = 1) ≪ |t − τ|, and LHS(σ = 1/2) ≫ τ |t − τ| (1−δ)/2 when θ ≤ π/4 and LHS(σ = 1/2) ≫ τ when θ ≥ π/4. This gives the existence of another path of constant imaginary part through s , which crosses the horizontal edges of V , and which is necessarily the path of steepest descent. The situation in the case m 6= 0 is less straightforward. Consider the case m > 0 (the case m < 0 is analogous). It is convenient to write θ = (t − τ)(1 + δ) log τ − 2πm − + so that, as t varies between t and t , the quantity θ varies between −π/2 and π/2; see m m again (3.2). Using (3.3), the equation Im f(s) = Im f(s ) is equivalent to 1−(1+δ)σ 1 τ (3.9) (t − τ) cos θ + σ sin θ = (t − t ) log x − v , m m 2 2 2 σ + (t − τ) where t −τ 2 2 log x σ +(t −τ) v = Im = log x . 1 2 2 σ t −τ m m (1 + δ) log τ + (1 + δ) log τ + + s −iτ 2 2 2 2 σ +(t −τ) σ +(t −τ) m m m m Note that log x (3.10) v ≪ m . (log τ) Denote the left hand side of (3.9) by l (σ), and the right hand side by r . t t BEURLING INTEGERS WITH RH AND LARGE OSCILLATION 13 it (2) it (1) it it it 0 1 1/2 Figure 1. Qualitative behavior of steepest paths for m > 0 (green for descent, red for ascent). The blue line is given by σ = (t − τ) cot|θ|. In Figure 1, the qualitative behavior of the steepest decent/ascent paths are depicted, as well as the relative positions of these paths with respect to some other points and curves. The (1) (2) (1) points t and t are deﬁned as the solution of r = 0, i.e., t = t +v / log x, respectively m m t m m (2) t = τ + 2πm/((1 + δ) log τ), corresponding to θ = 0. (Note that (3.7) and (3.10) imply (1) (2) − + t < t < t < t < t .) The blue line is the set of points for which l (σ) = 0. m m m t m m One can show the existence of the path of steepest descent by a careful analysis of (3.9). It is however more convenient to use the approximations for f and f given in Lemma 3.2 to show this existence, at least in a suﬃciently small neighborhood of the saddle point. We have 1 1 ′′ f (s ) = log x (1 + δ) log τ + + s − iτ (s − iτ) (1 + δ) log τ + m m s −iτ (3.11) = log x (1 + δ) log τ + O(1) , ′′ iφ so arg f (s ) ≪ 1/ log τ. Now ﬁx an ε > 0 suﬃciently small and write s − s = re . By m m Lemma 3.2, there are a function λ and η > 0 such that for r < η/ log τ, we have ′′ f (s ) 2 2iφ f(s) = f(s ) + r e (1 + λ (s)), λ (s) < ε. m m m How small we need ε to be will be determined later, but it is important to note that ε and later also η can be chosen independently of τ and m. 14 F. BROUCKE, G. DEBRUYNE, AND J. VINDAS ′′ 2 2iφ Set g (s) := f (s )r e (1 + λ (s))/2. The path of steepest descent is given by Im g = 0 m m m m under the constraint Re g < 0. Since ′′ f (s ) Re g (s) = r (1 + Re λ (s)) cos 2φ − (Im λ (s)) sin 2φ + O(1/ log τ) , m m m we see that by choosing ε suﬃciently small and τ large enough we must necessarily have φ ∈ (−7π/8,−π/8) ∪ (π/8, 7π/8) in order to satisfy the condition Re g (s) < 0. On the other hand, ′′ f (s ) Im g (s) = r (1 + Re λ (s)) sin 2φ + (Im λ (s)) cos 2φ + O(1/ log τ) . m m m For each r < η/ log τ, there is at least one solution φ near π/2 and at least one solution near iφ −π/2 of the equation Im g (re ) = 0. For example, by selecting ε suﬃciently small and τ large enough, there is a solution for φ in (2π/5, 3π/5) and in (−3π/5,−2π/5). This guarantees the existence of the path of steepest descent in the range |θ| ≤ η/2 say (2) 5/4 (since t = t + O(1/(log τ) )). Note that this part of the path lies inside the rectangle V , since B(s , η/ log τ) ⊆ V for suﬃciently large τ. When θ > η/2, r is positive; hence, m m m t by monotonicity of l (σ) and calculating the values at σ = 1/2 and σ = 1, we deduce that for every θ ∈ (η/2, π/2] there is a unique solution σ ∈ (1/2, 1) of (3.9). We conclude that we can extend the path of steepest descent from θ = η/2 upwards to θ = π/2, corresponding to t = t . It remains to treat the case θ < −η/2. Then r < 0, so that (3.9) can only have solutions when σ > (t − τ) cot|θ|, and we remark that for θ < −η/2 and suﬃciently large τ, (t − τ) cot|θ| < 1/2. Since l ((t − τ) cot|θ|) = 0 and l (σ) → 0 as σ → ∞, l has at least one t t t (local) minimum for σ > (t − τ) cot|θ|. In fact, l has precisely one minimum: 1−(1+δ)σ 1 τ l (σ) = (sin θ − ψ (σ)), 2 2 2 σ + (t − τ) with 2σ ψ (σ) = (t − τ) cos θ + σ sin θ (1 + δ) log τ + . 2 2 σ + (t − τ) Now 2 2 2 2σ 2(σ + (t − τ) ) − 4σ ψ (σ) = sin θ (1 + δ) log τ + + (t− τ) cos θ + σ sin θ , 2 2 2 2 2 σ + (t − τ) (σ + (t − τ) ) and the second term is bounded in absolute value by 2 2 2σ − 2(t − τ) 2σ|sin θ| , if σ ≥ (t − τ) cot|θ| . 2 2 2 (σ + (t − τ) ) This implies that ψ is monotonically decreasing for σ > (t−τ) cot|θ|, so that l has a unique t t local extremum in this range, which must be a minimum. Now |r | ≍ log x/ log τ, while t η (1−δ)/2 l (1/2) ≫ τ , so that for each ﬁxed θ < −η/2, (3.9) has two distinct solutions for σ, t η one on either side of 1/2. In particular , we can extend the path of steepest descent from θ = −η/2 (where it was on the right of 1/2) downwards to θ = −π/2, corresponding to t = t . This path cannot cross the line σ = 1, since l (1) ≪ 1/ log τ, which is of strictly lower order than r . The solution to the left of 1/2 corresponds to the path of steepest ascent. BEURLING INTEGERS WITH RH AND LARGE OSCILLATION 15 Now that we have shown the existence of the paths of steepest descent, we will also deduce some information about the argument of the tangent vector along these paths. Denote the − + steepest path by Γ , and let γ : [y , y ] → Γ be a unit speed parametrization of Γ m m m m m m ′ ± ± ′ (i.e. |γ | = 1) with Im γ (y ) = t and γ (0) = s . We have that γ (y) is a positive m m m m m m m ′ ′ multiple of ∇ Re f(γ (y)) = f (γ (y)) for y < 0, while γ (y) is a negative multiple of m m ∇ Ref(γ (y)) = f (γ (y)) for y > 0. m m iφ Write again s − s = re . From the above discussion, we know that for |θ| < η/2 this path lies in the cone φ ∈ (−3π/5,−2π/5) ∪ (2π/5, 3π/5). By appealing to Lemma 3.2 once more, we see that there exists a function λ such that ′ ′′ iφ f (s) = f (s )re (1 + λ (s)) m m −iπ/2 ′ with |λ (s)| < ε for r < η/ log τ. From this we deduce that arg(e γ (y)) < π/5 in the range|θ| < η/2, provided we choose ε suﬃciently small and τ suﬃciently large. For the range η/2 <|θ| ≤ π/2, we estimate the argument of f directly. We have 1−(1+δ)σ 1 τ f (s) = log x + (A − iB), 2 2 2 σ + (t − τ) with (t − τ) cos θ + σ sin θ (t − τ) A := (t − τ) sin θ − σ cos θ (1 + δ) log τ + + , 2 2 2 2 σ + (t − τ) σ + (t − τ) σ cos θ − (t − τ) sin θ (t − τ) B := (t − τ) cos θ + σ sin θ (1 + δ) log τ + + . 2 2 2 2 σ + (t − τ) σ + (t − τ) Whence we see that (with γ (y) = σ + it) 2 2 σ +(t−τ) log x · 2 + A 1−(1+δ)σ −iπ/2 ′ τ arg(e γ (y)) = arctan . The relation (3.9) holds on the path of steepest descent, so (t−τ) cos θ+σ sin θ + A (t−t )−v / log x m m −iπ/2 ′ arg(e γ (y)) = arctan . −9/4 −2 Now by (3.10), v / log x ≪ (log τ) and t − t = θ/((1 + δ) log τ) + O(m(log τ) ), so m m that, for η/2 <|θ| ≤ π/2, (t − τ) cos θ + σ sin θ (1 + δ) log τ 3/4 = (t − τ) cos θ + σ sin θ + O ((log τ) ) (t − t ) − v / log x θ m m σ(1 + δ) log τ sin θ 3/4 = + O ((log τ) ). Furthermore, we have 3/4 3/4 A = −σ(1 + δ) log τ cos θ + O((log τ) ) and B = σ(1 + δ) log τ sin θ + O((log τ) ), so sin θ −1/4 − cos θ + O ((log τ) ) −iπ/2 ′ θ arg(e γ (y)) = arctan −1/4 sin θ + O((log τ) ) −1/4 = arctan 1/θ − cot θ + O ((log τ) ) . η 16 F. BROUCKE, G. DEBRUYNE, AND J. VINDAS Now, 1/θ − cot θ is bounded on [−π/2, π/2] with 1/θ − cot θ ≤ 2/π, so we see that −iπ/2 ′ | arg(e γ (y))| ≤ π/5 on the range η/2 <|θ| ≤ π/2, and combining this with the estimate above, we see that this inequality holds in the entire range|θ| ≤ π/2. 3.3. The contribution from the saddle points. We will now estimate the contribution to the Perron integral in (2.2) coming from the integrals over the paths of steepest descent Γ . Since we take the imaginary part of this Perron integral, we need to control the argument of , and see that this is close to π/2, or at least suﬃciently far from 0 and π. For this we use the following lemma. Lemma 3.3. Let a < b and suppose that g : [a, b] → C is integrable. If there exist θ and η −iθ with 0 ≤ η < π/2 such that arg(ge ) ≤ η, then i(θ +ϕ) g(y) dy = ρe for some real numbers ρ and ϕ satisfying ρ ≥ (cos η) g(y) dy and |ϕ| ≤ η. iθ(y) Proof. Assume that g is not identically zero (that case is trivial) and write g(y) = R(y)e with θ(y) − θ ≤ η. Then, Z Z Z b b b iθ g(y) dy = e R(y) cos θ(y) − θ dy + i R(y) sin θ(y) − θ dy . 0 0 a a a The modulus of this expression is larger than R(y) cosη dy, while R(y) sin θ(y) − θ dy sin η |ϕ| = arctan ≤ arctan = η. cos η R(y) cos θ(y) − θ dy We have that X X ′ ′ Z Z exp . . . y exp . . . f(s) f(s ) f(γ (y))−f(s ) ′ m m m e ds = e e γ (y) dy. (s − 1)(s + 1) (γ (y) − 1)(γ (y) + 1) m m Γ y m m Let us ﬁrst focus on the argument of this integral. Using (3.3), log x f(s ) = (1 + σ ) log x + it log x + , m m m (1 + δ) log τ + s −iτ whence, in view of (3.7) and the bound on v (3.10), we obtain log log x 1 + 2 3 2πm log x |m| log x |m| log x Im f(s ) = log x τ + 1 − + O + . 3 4 (1 + δ) log τ (1 + δ) log τ (log τ) (log τ) arctan(2/π) ≈ 0.18π. BEURLING INTEGERS WITH RH AND LARGE OSCILLATION 17 By the technical assumptions imposed on the sequence (τ ) , more speciﬁcally properties (c) j j and (d), we see that the main term of the above expression has distance at most π/16 from an even multiple of π when k is even or from an odd multiple of π when k is odd. Furthermore, 1/3 2/3 by restricting the range for m to |m| ≤ c(log x) (log log x) for some suﬃciently small absolute constant c > 0, we force the error term to be at most π/16. Hence, for |m| ≤ 1/3 2/3 c(log x) (log log x) , π π d Im f(s ), 2Zπ < for k even and d Im f(s ), π + 2Zπ < for k odd. m m 8 8 Second, in the range|t − τ | ≤ 1, 1/2 ≤ σ ≤ 1 say, we have by property (a) that −(1+δ )s −ν s ∞ −(1+δ )s −ν s ∞ j j j j X X X X ′ τ − τ τ − τ 1 1 j j j j −1/2 = τ + τ ≪ τ . j j 2 s − iτ 2 s + iτ j j j6=k j=0 j=0 We can make this smaller than π/16 by choosing τ suﬃciently large. Third, in the same range we have 1 1 1 1 1 = + O , arg = π + O . 2 3 (s − 1)(s + 1) τ τ (s − 1)(s + 1) τ k k −iπ/2 ′ Finally, by the results of the previous subsection, |arg(e γ (y))| < π/5. Combining all of this, and using the fact that on Γ , f(s) − f(s ) is real, we get by Lemma 3.3 that m m exp . . . f(s) k+1 i(π/2+ϕ ) e ds = (−1) R e , (s − 1)(s + 1) where|ϕ | < 2π/5 and R is a positive number satisfying m m Z + Re f(s ) y m m R ≫ exp Re(f(γ (y)) − f(s )) dy. m m m k+1 We see that the imaginary part of the integral over Γ always has sign (−1) , so the k+1 sum over m of the imaginary parts of all these integrals has also sign (−1) and is in absolute value larger than the contribution of the integral over Γ . Therefore, the saddle point contribution of the integrals Γ , m 6= 0, cannot destroy the contribution of Γ , and as m 0 we have better control over σ compared to general σ ((3.6) versus (3.5)), we will proceed 0 m to get a lower bound for the contribution of the saddle points only using the integral Γ . To estimate R , we restrict the integration interval to the interval which corresponds to the part of Γ which lies inside B(s , η/ log τ) for some η > 0. By Lemma 3.2 and (3.11), we see 0 0 that exp Re(f(γ (y)) − f(s )) dy ≫ √ . 0 0 log x log τ Now by (3.3), we have 1−(1+δ)σ 1 τ log x (3.12) = 2 + O(1), 2 σ log log x With some more calculations, one can show that the size of the contributions of the saddle points s to the Perron integral decays relatively slowly with respect to m. Therefore, if the phases of the contributions were arbitrary, the contribution of the saddle point s might get cancelled by the others. As a consequence we were led to this intricate analysis of the phases of the contributions. 18 F. BROUCKE, G. DEBRUYNE, AND J. VINDAS and by using (3.6), we see that Re f(s ) √ p √ e log x = exp 2 log x− 2 log x log log x− 2(a+log 2−1) −2 log τ+O(log log x) . τ log log x 1/3 2/3 Hence, we can conclude that the contribution of saddle points with |m| ≤ c(log x) (log log x) k+1 has sign (−1) and is bounded in absolute value from below as follows: exp . . . f(s) Im e ds π (s − 1)(s + 1) √ √ log x (3.13) ≫ exp (1 + σ ) log x − 2 log x log log x + 2 + O(log log x) log log x √ p log x (3.14) = x exp −2 2 log x log log x − 2(a + log 2 − 1) + O(log log x) . log log x 4. The remainder Recall that the Perron integral we are considering is given by Z exp . . . f(s) (4.1) e ds. (s − 1)(s + 1) In this section we will show that one can integrate over a contour Γ which incorporates the paths of steepest descent Γ , in such a way that the integral over Γ\∪ Γ is of strictly lower m m m order than the contribution of the saddle points. For brevity, we will omit the integrand when writing integrals. The omitted integrand is always meant to be the integrand of (4.1). Let us begin by bounding the zeta function on the contour of Hilberdink and Lapidus. Lemma 4.1. For t ≥ e , the function −(1+δ )s −ν s −(1+δ )s −ν s j j j j τ − τ τ − τ j j j j j + , resp. 2 s − iτ s + iτ j j j=0 −(1+δ )s −ν s −(1+δ )s −ν s j j ∞ j j X X τ − τ τ − τ τ τ j j j j j j + , 2 s − iτ 2 s + iτ j j j6=k j=0 is bounded in the region log log t 1/5 s = σ + it : σ ≥ 1 − , resp. { s = σ + it : σ ≥ 1/2, τ ≤ t ≤ τ }. k k log t Proof. For σ ≥ 1 − log log t/ log t, we have log log τ + a log log t log log τ + a log log t 1−(1+δ )σ j j j j j τ ≤ exp log τ − + + log τ log t log τ log t j j log log t = exp − log log τ − a + log τ + log log τ + a . j j j j j log t BEURLING INTEGERS WITH RH AND LARGE OSCILLATION 19 Since log log t/ log t ≤ 1/2, this is ≪ τ / log τ . Furthermore, if τ /2 ≤ t, the above is j j j log(log τ − log 2) ≪ exp − log log τ + log τ + O(1) ≪ 1. j j log τ − log 2 1−(1+δ )σ j 2/3 Also, for σ ≥ 1/2, τ ≤ τ say. For the ﬁrst case, we bound s − iτ from below by j j 1/2 if τ /2 ≤ t ≤ 2τ and by τ /2 otherwise, while for the second case, we bound s − iτ , j j j j j 6= k, from below by τ /2, in view of property (a). Hence, we see that the functions are s ! ∞ ∞ X X 1 1 ≪ + 1, resp. , 1/3 τ log τ j j j=0 j=0 j which are bounded (by property (a)). The lemma implies that we may indeed shift the contour from the line σ = κ to the Hilberdink-Lapidus contour Γ as described in Section 2. Indeed, the integral over the H-L segments log log|T| 1 − ± iT, κ ± iT log|T| κ+1 2 is O (x /T ), which tends to 0 as T → ∞. 1/3 2/3 4.1. Connecting the steepest paths. Set m := ⌊c(log x) (log log x) ⌋, and max 2πm + π/2 max ± ± T := τ ± = t . 1 ±m max (1 + δ) log τ In this subsection we will show that one can connect the diﬀerent paths of steepest descent − + Γ to form one contour whose imaginary part ranges from T to T . 1 1 1/3 2/3 ± ± For m in the range|m| ≤ c(log x) (log log x) , set σ = Re γ (y ). By (3.9) (the value m m of θ is here either π/2 or −π/2 correspondingly, cf. (3.2)), we see that these numbers satisfy 1 1 (t − t ) − v / log x m m ± ± 2 ± 2 σ = 1− log log x+log 2+log (σ ) +(t −τ) +log . m m m 1 + δ log τ ±σ Using that 1/3 π 1 (log log τ) v 1 ± ± 1/2 < σ < 1, t − t = ± + O , ≪ , m m 4/3 9/4 2 (1 + δ) log τ (log τ) log x (log τ) we get log log x 2(a + log 2 + log π/2) (4.2) σ = 1 − 2 − + O . 2/3 1/3 log x log x log log x (log x) (log log x) + + − − Consider now the contour Υ which connects σ + it with σ + it via a vertical and m m m+1 m+1 − − − − + + + + horizontal line: Υ := [σ + it , σ + it ] ∪ [σ + it , σ + it ]. We have m m m m+1 m m+1 m+1 m+1 1−(1+δ)s 1−(1+δ)σ 1 τ 1 τ Re = σ cos((t − τ)(1 + δ) log τ) − (t − τ) sin((t − τ)(1 + δ) log τ) . 2 2 2 s − iτ 2 σ + (t − τ) On this contour, 1/3 log log τ σ cos (t − τ)(1 + δ) log τ − (t − τ) sin (t − τ)(1 + δ) log τ ≤ t − τ ≪ , |m|+1 log τ 20 F. BROUCKE, G. DEBRUYNE, AND J. VINDAS since the argument of cos belongs to [π/2, 3π/2] + 2πZ. Therefore, by (4.2), 1−(1+δ)s 1−(1+δ)σ 1 τ 1 τ Re ≤ (t − τ) 2 2 |m|+1 2 s − iτ 2 σ + (t − τ) 1/6 1/3 log log x log x 1−(1+δ)σ ≪ τ ≪ . log x log log x Hence, in view of (4.2) and using Lemma 4.1 to bound the series , we see that the integrand of (4.1) over Υ is 1/3 √ p √ log x π log x ≪ x exp −2 2 log x log log x− 2 a + log 2 + log +O . 2 log log x log log x 1/3 2/3 Summing over |m| ≤ c(log x) (log log x) , we get that is of lower order than the ∪ Υ m m main contribution given by (3.14). 4.2. Returning to the original contour. Finally we show that we can connect the contour ∪ Γ ∪ Υ to Γ (see Figure 2). First we go from the endpoints of Γ to the line m m m H-L ±m max σ = σ : set − − − − + + + + : : Δ = [σ + iT , σ + iT ], Δ = [σ + iT , σ + iT ]. 0 0 0 −mmax 1 1 0 mmax 1 1 On this contour, we have 1−(1+δ)s 1 τ Re < 0, 2 s − iτ R √ which together with Lemma 4.1 implies that ± ≪ exp (1+σ ) log x− 2 log x log log x , which is admissible in view of (3.13). Suppose that 1/6 2πm log log x max |t − τ| ≥ ≥ cπ 2 . (1 + δ) log τ log x Then 1−(1+δ)σ 1−(1+δ)σ 0 0 1 τ 1 τ 1 = q 2 2 σ 2 σ + i(t − τ) 1 + (t − τ)/σ 1/3 1−(1+δ)σ 1 τ log log x ≤ 1 − c , 2 σ log x for some constant c > 0. By (3.12), we thus have 1/6 1−(1+δ)σ 1 τ log x log x ′′ ≤ 2 − c + O(1), 2 log log x log log x σ + i(t − τ) ′′ for some constant c > 0. We now set 1/6 ′′ c log x ± ± ± ± (4.3) T := τ ± exp , Δ := [σ + iT , σ + iT ]. 0 0 2 1 1 2 2 log log x BEURLING INTEGERS WITH RH AND LARGE OSCILLATION 21 iT iT iτ iT iT σ σ 0 0 1 Figure 2. Returning to the original contour of Hilberdink and Lapidus. On Δ , t ≍ τ and the series is bounded by Lemma 4.1. Therefore, we have 1/6 ′′ c log x 1+σ ≪ exp x × 2 log log x 1/6 √ p √ log x log x ′′ exp − 2 log x log log x + 2 − c + O(1) , log log x log log x which is of lower order than the contribution of the saddle points (3.13). Next, set ′′ 1 2 c σ := 1 − 1/3 2/3 1 + δ 2 (log x) (log log x) ′′ 2 c log log x = 1 − + O . 1/3 2/3 2 (log x) (log log x) log x 22 F. BROUCKE, G. DEBRUYNE, AND J. VINDAS ′ + − We point out that σ < σ (for x suﬃciently large). Then, for t ≥ T or t ≤ T , 2 2 1/6 1/6 1−(1+δ)σ ′′ ′′ 1 τ c log x c log x (4.4) ≪ exp − + = 1. 2 σ + i(t − τ) 2 log log x 2 log log x R √ ± ± ± ′ 1+σ Set Δ := [σ +iT , σ +iT ]. We obtain that ≪ x exp − 2 log x log log x , which 2 2 2 is again admissible, by comparing with (3.13). + − ± ± ± 5 1/5 ′ ′ Consider now T := τ , T := τ , and Δ := [σ + iT , σ + iT ]. We get (again 3 3 3 2 3 bounding using Lemma 4.1 and using (4.4)) 2/3 2 log x 2 ′′ ≪ x exp − c + O( log x log log x) , 2 log log x which is negligible. ± ′ ± ± ± ± Next we move to the contour Γ : set Δ := [σ +iT , 1−log log T / log T +iT ]. Now H-L 4 3 3 3 3 1/5 log log τ 9 log log x ≤ 1 − 2 1 − 1/5 log τ 4 log x R √ say, so again by (4.4) and Lemma 4.1 ≪ x exp − 2 log x log log x , which is admis- Δ 4 sible with respect to (3.14). Also, 1 5 2 = exp − log x log log x , τ 2 so + is admissible as well. − 1/5 + 5 : : : : Finally, set Γ = { s ∈ Γ t ≤ τ }, Γ = { s ∈ Γ t ≥ τ }. By Lemma 4.1, H-L H-L H-L H-L the series is bounded on these contours. We get Z Z (log τ)/5 1/5 log u log log τ 2 2 ≪ x exp − log x − u du ≪ x exp − log x 1/5 u log τ Γ e H-L √ p ≪ x exp − 2 log x log log x , Z Z log u 2 2 ≪ x exp − log x − u du ≪ x exp(−5 log τ) Γ 5 log τ H-L ≪ x exp − 2 log x log log x , and these integrals are therefore also negligible. 5. Conclusion of analysis of continuous example: proof of Lemma 2.1 The results of the previous two sections now yield the relation (2.3), which as we have already remarked at the end of Section 2 suﬃces to conclude the proof of Theorem 1.3. All of this only remains true provided that Lemma 2.1 holds. Let us now prove Lemma 2.1. 11 e 2−1/e The part of the integral for t < e is ≪ x . BEURLING INTEGERS WITH RH AND LARGE OSCILLATION 23 Proof of Lemma 2.1. We will show that for any T one can ﬁnd τ , a , and ν which obey k k k τ > T and the properties (b)-(d). This will allow us to inductively deﬁne the sequences for which also (a) holds, by selecting T := (2τ ) in step k + 1. Set α := log 6 + 1/2, and pick a number ξ > e such that log ξ √ × 1 1 1 log 2 log ξ log log ξ + log log ξ + log log log ξ + α − k k k k 2 2 2 log log ξ 1 + 2 log ξ 1 − ∈ Z 1 1 1 log 2 log ξ log log ξ + log log ξ + log log log ξ + α − k k k k 2 2 2 (corresponding to property (d)). This is possible since this function is continuous and tends to ∞ as ξ → ∞. Next, deﬁne x := exp(log ξ + ε ), with ε the smallest positive number k k k k k such that 1 2πZ if k is even log x exp √ log x log log x ∈ k k k π + 2πZ if k is odd. Applying the mean value theorem to the function y exp y log y , we see that ε ≪ √ √ . log ξ log log ξ exp log ξ log log ξ k k k k Hence, by replacing log ξ by log x = log ξ + ε , we introduce an error in the condition for k k k k (d) of order y 1 ≪ ε ≪ √ log y log ξ log log ξ exp log ξ log log ξ k k k k y=log ξ k 2 ≍ , log x log log x exp log x log log x k k k k which is admissible. Next, set a = α + η , with η the smallest positive number such that k k k r r log x log log x log x log log x 1 1 log 2 k k k k exp + log log x + log log log x − +α+η k k k 2 2 2 2 2 belongs to 2πZ, corresponding to the ﬁrst requirement of property (b). Then η ≪ exp −√ log x log log x , k k k and replacing α by a = α + η , the newly introduced error in the condition for (d) is k k ≪ √ , log log x exp log x log log x k k k log x log log x , and we can choose a which is admissible. Finally, we set τ = exp k k k value for ν between 2 and 3 to satisfy the second part of property (b). All the properties are now fulﬁlled.  24 F. BROUCKE, G. DEBRUYNE, AND J. VINDAS 6. Discretization: proof of Theorem 1.1 This last section is devoted to completing the proof of Theorem 1.1. We will apply the probabilistic approach of Diamond, Montgomery, and Vorhauer [10] and complement it with a new procedure of adding ﬁnitely many well-chosen primes. This allows us to obtain a suitable random approximation to the continuous prime measure dΠ that we have been studying in the previous sections and will enable us to select a random discrete Beurling prime number system having the desired properties (1.4) and (1.5). Let 1 = v < v < v < . . . be a ﬁxed sequence of real numbers tending slowly to ∞ (how 0 1 2 slowly will be speciﬁed later), and set q = dΠ (u). j C j−1 We will include the number v as a prime in our discrete prime number system with prob- ability q , where as our ﬁrst requirement on the v we ask that they increase suﬃciently j j slowly such that all q < 1. To make this precise, let (X ) be a sequence of independent j j j>0 Bernoulli variables with parameters q of success on a ﬁxed probability space. Given a point ω in the probability space, denote by P(ω) the set of those v for which X (ω) = 1. The j j idea is to show that the probability that a prime number system P(ω) satisﬁes the bounds we need is nonzero. Denote Z + Z y y −it −it S(y; t) = S (y; t) := u dπ(u), S (y; t) := u dΠ (u), ω C C 1 1 where π = π is the prime counting function of the discrete system P(ω). We require the following bounds: (A) π(y) = Π (y) + O( y) for y suﬃciently large; (B) S(y; t) = S (y; t) + O( y log|t|) uniformly for y and|t| suﬃciently large; 1/4 (C) S(y; t) = S (y; t) + O( y(log τ ) ) uniformly in the range C k 1/6 ′′ c log x |t − τ | ≤ exp , 2 log log x ′′ for y and k suﬃciently large. The constant c is the same constant as the one appearing in (4.3). Lemma 6.1. For every ε > 0, there exist Y , T , and K such that the probability that a ε ε ε prime number system P(ω) satisﬁes the bounds (A)-(C) with y ≥ Y , |t| ≥ T , and k ≥ K ε ε ε is at least 1 − ε. We will only prove the lemma for the bounds (C). The proof of the validity of the bounds (A) and (B) is identical to that of [10, Lemma 9]. In fact, the only assumption that is needed in the proof of [10, Lemma 9] is that the measure dΠ satisﬁes du 2 du ≪ dΠ (u) ≤ , log(2u) log 2u which of course holds in our case as well. We will employ the following inequality, which follows from an equality of Kolmogorov for sums of independent random variables (see [16, Chapter V] or [10, Lemma 8, p. 17]). Let X be independent Bernoulli variables with j BEURLING INTEGERS WITH RH AND LARGE OSCILLATION 25 parameter q , and suppose that r are real numbers with r ≤ 1 (j = 1, . . . , J). Set j j j X = r X . If j j j=1 (6.1) 0 ≤ v ≤ 2 q (1 − q ), j j j=1 then −v (6.2) P (X ≥ E(X) + v) ≤ exp . 4 q (1 − q ) j j j=1 Proof of Lemma 6.1. As previously indicated, we only show the bound (C), and we will assume that the bound (A) holds . By the trivial estimates S(y; t), S (y; t) ≪ y/ log y, we may assume that y ≥ C log τ (log log τ ) for some ﬁxed but arbitrarily large constant k k 1/4 C > 0. We apply the inequality (6.2) with r = cos(t log v ) and v = y(log τ ) . Let J k k k be such that v ≤ y < v . Using that (provided that q ≤ 1/2) J J+1 j Π (v ) ≤ q (1 − q ) ≤ Π (v ), C J j j C J j=1 we see that (6.1) holds, so that by (6.2), y log τ 1 1/4 P (ReS(y; t) ≥ E(ReS(y; t))+ y(log τ ) ) ≤ exp − ≤ exp − log y log τ . k k 4Π (y) 8 The fact that in the range 1/6 ′′ c log x |t − τ | ≤ exp and y ≥ C log τ (log log τ ) k k k 2 log log x E(S(y; t)) is close to S (y; t) can be proven in exactly the same way as in [10, pp. 21–22], although we need a diﬀerent choice for the sequence v . The reader may check that the 1/4 choice v = (log j) (j ≥ j ) is adequate for obtaining E(S(y; t)) − S (y; t) ≤ y in the j 0 C given ranges (and also works ﬁne for the proof of the bounds (A) and (B)). Hence, √ 1 1/4 P ReS(y; t) ≥ Re S (y; t) + 2 y(log τ ) ≤ exp − log y log τ . C k k Applying the same argument for r = − cos(t log v ) and r = ± sin(t log v ) gives that j j j j √ 1 1/4 P S(y; t) − S (y; t) ≥ 4 y(log τ ) ≤ 4 exp − log y log τ . C k k 1/4 Let C denote the event S(m; n) − S (m; n) ≥ 4 m(log τ ) for some integer n in the mk C range 1/6 ′′ c log x |n − τ | ≤ exp . 2 log log x 12 c c c If P (A ) ≤ ε, we can bound the probability of an event D as P (D) = P (D|A)P (A) + P (D|A )P (A ) ≤ P (D|A)P (A) + ε. 26 F. BROUCKE, G. DEBRUYNE, AND J. VINDAS ′′′ Using the relation (2.1) between τ and x , we get for some constant c k k 1/3 X p p 1 1 log τ ′′′ P (C ) ≪ exp − log m log τ ≪ exp − log m log τ + c mk k k 8 8 log log τ ≤ exp − log m log τ , provided that k is suﬃciently large, but otherwise independent of m. Thus, by the rapid growth of the sequence τ (namely, property (a) from Section 2), X X X X 1 k P (C ) ≪ exp − log m · 5 log τ mk 0 k≥K m≥M k≥K m≥M ≤ exp − log τ − 1 log(M − 1) < ∞. k≥K We conclude that [ [ ∀ε > 0∃M , K ∈ N: P C < ε. ε ε mk m≥M k≥K ε ε T T Consider now an event in the complement ω ∈ C . Then the bound (C) m≥M k≥K mk ε ε holds for integral y and t. One then uses that|y − y | ≤ 1 =⇒ S(y ; t) = S(y ; t)+O( y ) 1 2 1 2 1 (since π(y) = Π (y) + O( y)) and that i(t −t ) i(t −t )−1 2 1 2 1 S(y; t ) = y S(y; t ) − i(t − t ) S(u; t )u du 1 2 2 1 2 (and similarly for S ) to see that the bound (C) also holds for non-integral y and t. We now ﬁx an event ω with corresponding prime number system P = P(ω ) for which 0 0 0 the bounds (A)-(C) hold. Denote the zeta function of this prime number system by ζ . The bounds (A)-(C) imply the following bounds for ζ . Lemma 6.2. The function log ζ (s)− log ζ (s) admits an analytic continuation to σ > 1/2. 0 C Uniformly for σ ≥ 1/2 + ε we have log ζ (s) = log ζ (s) + O ( log(|t| + 2)), 0 C ε while for k suﬃciently large, 1/4 log ζ (s) = log ζ (s) + O ((log τ ) ), 0 C ε k ′ ′ 1/4 (log ζ (s)) = (log ζ (s)) + O ((log τ ) ) 0 C ε k uniformly in the range 1/6 ′′ c log x σ ≥ 1/2 + ε, |t − τ | ≤ exp . 2 log log x One can check that the above bounds are strong enough so that one can repeat the proof of Theorem 1.3 with ζ instead of ζ along the same contour, except for the estimation of 0 C the contribution from the saddle points. For the argument in Subsection 3.3 to go through, we would also require that on the paths of steepest descent (6.3) d(Im(log ζ(s) − log ζ (s)), 2πZ) < π/20, C BEURLING INTEGERS WITH RH AND LARGE OSCILLATION 27 say, but it is unclear whether this holds true for ζ = ζ . We will thus modify the prime number system P by adding a ﬁnite number of primes such that this bound holds, at least inﬁnitely often for subsequences of (τ ) and (τ ) . 2k k 2k+1 k Set π π π π S = m − , m + + 2πZ, for m = 0, 1, . . ., 159. 80 160 80 160 By the pigeonhole principle, there is a number m (resp. l) such that for inﬁnitely many even k (resp. odd k), Im(log ζ (1 + iτ ) − log ζ (1 + iτ )) ∈ S (resp. S ). 0 k C k m l Suppose without loss of generality that l ≤ m. We will add m times the prime p to P , −s where p is a well chosen number around 80/π. This changes log ζ (s) by −m log(1 − p ). This additional term and its derivative are O(1), so Lemma 6.2 still holds for the new zeta function. In 1 + iτ , the imaginary part of log ζ changes by k 0 −1 p sin(τ log p) −1−iτ −m arg(1 − p ) = −m arctan . −1 1 − p cos(τ log p) Let α be a solution of sin α l = , 0 ≤ α ≤ π/2. 1 − cos α m We set 80 1 p := e , where ε = ε , ε ≪ . k k π τ k=0 We deﬁne the numbers ε inductively: suppose ε , ε , . . . , ε are deﬁned. Set ε = λ /τ , k 0 1 k−1 k k k with λ ∈ [0, 2π) the unique number such that 80 π τ log + ε + ε + · · · + ε ∈ + 2πZ ( resp. ∈ α + 2πZ), k 0 1 k π 2 for k even (resp. odd). Suppose now that k is even (the reasoning for odd k is completely analogous). Using the rapid growth of the sequence τ , that is, property (a) from Section 2, we get 80 τ τ log p = τ log + ε + ε + · · · + ε + O k k 0 1 k π τ k+n n=1 −4 = + 2πM + O(τ ), for some integer M . Then, −8 −4 −1 sin(τ log p) = 1 + O(τ ), cos(τ log p) = O(τ ), p = + O(τ ), k k k k 0 so that −1 p sin(τ log p) π −1 = + O(τ ). −1 1 − p cos(τ log p) 80 Since|arctan u − u| < 3|u| for|u| < 1, we have (for τ suﬃciently large) for every even k π π π −1−iτ Im(−m log(1 − p )) + m < 6m < , 80 80 40 28 F. BROUCKE, G. DEBRUYNE, AND J. VINDAS and similarly for every odd k π π −1−iτ Im(−m log(1 − p )) + l < . 80 40 −s Set F (s) := log ζ (s) − log ζ (s) − m log(1 − p ). Then, for an inﬁnite number of even 0 C k and an inﬁnite number of odd k, d(Im F (1 + iτ ), 2πZ) < 5π/160. To see that such a (slightly weaker) bound also holds on the corresponding paths of steepest descent, write F (s) = F (1 + iτ ) + F (z) dz. 1+iτ For s on such a path we have 1/3 2/3 log log τ (log x ) (log log x ) k k k |s − 1 − iτ | ≪ + , log τ log τ k k and using the bound on the derivatives from Lemma 6.2, we get 2/3 1/3 1/3 (log τ ) (log log τ ) (log log τ ) k k k ′ 1/4 F (z) dz ≪ (log τ ) ≪ , 1/12 log τ (log τ ) k k 1+iτ so that for an inﬁnite number of even k (resp. odd k), we have d(ImF (s), 2πZ) < π/20 for s on the corresponding paths of steepest descent. −s −m Therefore, the bound (6.3) holds for ζ(s) = ζ (s)(1−p ) , the zeta function correspond- ing to the number system P obtained by adding m times the prime p to P . This allows one to estimate the contribution from the saddle points as in Subsection 3.3, and hence to deduce that (1.5) holds for the counting function of the set of generalized integers associated to P. The prime counting function of this number system obviously satisﬁes (1.4) as well. This concludes the proof of Theorem 1.1. Appendix A. Improving the constant in Theorem 1.2 We indicate here how Balazard’s method [3] yields an improvement for the value of c in Theorem 1.2 over the value given by Hilberdink and Lapidus [14]. Theorem A.1. Suppose that the generalized Riemann prime counting function satisﬁes −1 1 − u Π(x) = du + O(x ), log u for some 0 ≤ θ < 1. Then, there is ρ > 0 such that the generalized integer counting function satisﬁes, for each c < 2(1 − θ), N(x) = ρx + O x exp −c log x log log x . Naturally, when θ ≥ 1/2, Theorem A.1 implies Theorem 1.2 and improves the values of the constant c. In its proof, we shall make extensive use of the operational calculus for multiplicative convolution of measures for which we can refer the reader to [11, Chapter 2]. Proof. We deﬁne dP as the “natural” continuous prime measure and dE as the deviation between dΠ and dP , that is, −1 1 − u dP = du, dE = dΠ − dP. log u BEURLING INTEGERS WITH RH AND LARGE OSCILLATION 29 By hypothesis, E(x) := dE = O(x ). We recuperate N as Z Z Z x x x ∗ ∗ ∗ −1 ∗ N(x) = exp (dΠ) = exp (dP ) ∗ exp (dE) = x u exp (dE(u)) − − − 1 1 1 Z Z ∗n x x dE(u) 1 dE(u) = x exp = x , u n! u − − 1 1 n=0 where we have used the identity exp (dP ) = δ + du, where δ denotes the Dirac measure 1 1 concentrated at 1. We are therefore led to the integrals ∗n dE(u) I = . −(1−θ) −2 The ﬁrst two cases are trivial, I = 1 and I = b + O(x ), where b = u E(u)du. 0 1 For larger n, we claim that there is an absolute constant A for which n n−1 A (log log x) (A.1) I = b + O , (1−θ)/n and where the O-constant is independent of n. We proceed to show the claim via induction. Employing the Dirichlet hyperbola method we decompose I as S +S −S , where y shall be taken as the free variable to be optimized n+1 1 2 3 at our will. The ﬁrst term is   Z Z   Z ∗n y x/u y n n−1 (1−θ)/n dE(u) dE(v) A (log log x) u dE(u)   S := = b + O (1−θ)/n u v x u 1 1 1 n n n−1 b A (log log x) n+1 −1+(1−θ)/n = b + O + O u |dE(u)| . 1−θ (1−θ)/n y x Subsequently, as |dE| ≤ dΠ + dP = 2dP + dE, the integral in the ﬁnal error term can, for large enough x, be estimated by Z Z y x −1+(1−θ)/n (1−θ)/n −1 (1−θ)/n u |dE(u)| ≤ y u (2dP (u) + dE(u)) ≤ 3y log log x. 1 1 The second summand is estimated using the representation   Z Z Z ∗n ∗n x/y x/v x/y 1−θ dE(v) dE(u) v dE(v)   S = = b + O 1−θ v u x v 1 1 1 ! ! ∗n n n (1−θ)/n x/y A (log log x) y 1 |dE(v)| n+1 = b + O + O (1−θ)/n 1−θ x y v Considering the integral involving the n-folded convolution of |dE|, we obtain, for x large enough, ! ! n n Z Z Z ∗n x/y x/y x/y |dE(v)| |dE(v)| 2dP (v) + dE(v) ≤ ≤ v v v 1 1 1 ≤ (3 log log x) . 30 F. BROUCKE, G. DEBRUYNE, AND J. VINDAS Finally, S is bounded using Z Z ∗n y x/y n n n−1 (1−θ)/n dE(u) dE(v) b A (log log x) y n+1 S := = b + O + O . 1−θ (1−θ)/n u v y x 1 1 1/(n+1) Upon choosing y = x and A suﬃciently large such that the extra factor A can absorb all the absolute O-constants, the claim (A.1) follows. Now it remains to analyze ∞ ∞ n n−1 X X I A (log log x) (A.2) = e + O . (1−θ)/n n! n!x n=0 n=1 We now select the n, say n , for which 2 times the n-th term in the above series reaches max its maximum. A few standard calculations allows one to ﬁnd an approximation for n , max log x log log log x n = 2(1 − θ) 1 + O , max log log x log log x and after inserting this into (A.2), one obtains   1/2 p p (log x) log log log x   N(x) = e x + O x exp − 2(1 − θ) log x log log x + O . 1/2 (log log x) References [1] F. A. Al-Maamori, Examples of Beurling prime systems, Math. Slovaca 67 (2017), 321–344. [2] E. P. Balanzario, An example in Beurling’s theory of primes, Acta Arith. 87 (1998), 121–139. [3] M. Balazard, La version de Diamond de la m´ethode de l’hyperbole de Dirichlet, Enseign. Math. 45 (1999), 253–270. [4] P. T. Bateman, H. G. Diamond, Asymptotic distribution of Beurling’s generalized prime numbers, in: Studies in number theory, W. J. LeVeque (ed.), pp. 152–210, Mathematical Association of America, [5] A. Beurling, Analyse de la loi asymptotique de la distribution des nombres premiers g´en´eralis´es, Acta Math. 68 (1937), 255–291. [6] H. Bohr, Bidrag til de Dirichlet’ske Raekkers Theori, Thesis, Copenhagen, 1910. [7] N. G. de Bruijn, Asymptotic methods in analysis, Third edition, Dover Publications, Inc. New York, [8] G. Debruyne, J.-C. Schlage-Puchta, J. Vindas, Some examples in the theory of Beurling’s generalized prime numbers, Acta Arith. 176 (2016), 101–129. [9] H. G. Diamond, Asymptotic distribution of Beurling’s generalized integers, Illinois J. Math. 14 (1970), 12–28. [10] H. G. Diamond, H. L. Montgomery, U. M. A. Vorhauer, Beurling primes with large oscillation, Math. Ann. 334 (2006), 1–36. [11] H. G. Diamond, W.-B. Zhang, Beurling generalized numbers, Mathematical Surveys and Monographs series, American Mathematical Society, Providence, RI, 2016. [12] R. Estrada, R. P. Kanwal, A distributional approach to asymptotics. Theory and applications, Second edition, Birkh¨auser, Boston, 2002. [13] R. Garunkˇstis, L. Kaziulyte˙, Asymptotic distribution of Beurling integers, Int. J. Number Theory 14 (2018), 2555–2569. [14] T. W. Hilberdink, M. L. Lapidus, Beurling zeta functions, generalised primes, and fractal membranes, Acta Appl. Math 94 (2006), 21–48. [15] E. Landau, Neuer Beweis des Primzahlsatzes und Beweis des Primidealsatzes, Math. Ann. 56 (1903), 645–670. BEURLING INTEGERS WITH RH AND LARGE OSCILLATION 31 [16] M. Loeve, Probability theory. I, Graduate Texts in Mathematics 45, Fourth edition, Springer-Verlag, New York-Heidelberg, 1977. [17] P. Malliavin, Sur le reste de la loi asymptotique de r´epartition des nombres premiers g´en´eralis´es de Beurling, Acta Math. 106 (1961), 281–298. [18] G. Tenenbaum, Introduction to analytic and probabilistic number theory, Third edition, Graduate Studies in Mathematics, 163, American Mathematical Society, Providence, RI, 2015. [19] W.-B. Zhang, Beurling primes with RH and Beurling primes with large oscillation, Math. Ann. 337 (2007), 671–704. Department of Mathematics: Analysis, Logic and Discrete Mathematics, Ghent Univer- sity, Krijgslaan 281, 9000 Gent, Belgium E-mail address: fabrouck.broucke@ugent.be E-mail address: gregory.debruyne@ugent.be E-mail address: jasson.vindas@ugent.be http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Mathematics arXiv (Cornell University)

# Beurling integers with RH and large oscillation

, Volume 2020 (2004) – Apr 24, 2020
31 pages

/lp/arxiv-cornell-university/beurling-integers-with-rh-and-large-oscillation-75pA2ti6jA
ISSN
0001-8708
eISSN
ARCH-3343
DOI
10.1016/j.aim.2020.107240
Publisher site
See Article on Publisher Site