arxiv: v2 [math.nt] 16 Mar 2018

Transcription

1 EXPLICIT BOUNDS FOR PRIMES IN ARITHMETIC PROGRESSIONS MICHAEL A BENNETT, GREG MARTIN, KEVIN O BRYANT, AND ANDREW RECHNITZER arxiv: v [mathnt] 6 Mar 08 ABSTRACT We derive eplicit upper bounds for various functions counting primes in arithmetic progressions By way of eample, if q and a are integers with gcda, q = and 3 q 00, and θ; q, a denotes the sum of the logarithms of the primes p a mod q with p, we show that θ; q, a φq < 000 log, for all 3,460,547,899 with sharper constants obtained for individual such moduli q We establish inequalities of the same shape for the other standard prime-counting functions π; q, a and ψ; q, a, as well as inequalities for the nth prime congruent to a mod q when q 500 For moduli q > 0 5, we find even stronger eplicit inequalities, but only for much larger values of Along the way, we also derive an improved eplicit lower bound for L, χ for quadratic characters χ, and an improved eplicit upper bound for eceptional zeros CONTENTS Introduction and statement of results Preparation of the upper bound for ψ; q, a /φq, for q Previous work based on the eplicit formula 6 Upper bounds for V q,m ; δ, H, eploiting verification of GRH up to bounded height 8 3 Further estimates related to vertical distribution of zeros of Dirichlet L-functions 4 Preliminary statement of the upper bound for ψ; q, a /φq 4 3 Optimization of the upper bound for ψ; q, a /φq, for q Elimination of eplicit dependence on zeros of Dirichlet L-functions 9 3 Estimation of integrals using incomplete modified Bessel functions 33 Elementary estimation of the complementary error function erfcu 4 34 Identification of maimum values of bounding functions via calculus 6 35 Assembly of the final upper bound for ψ; q, a /φq 3 4 Deduction of the upper bounds upon θ; q, a /φq and π; q, a li/φq, for q Conversion of bounds for ψ; q, a /φq to bounds for θ; q, a /φq 35 4 Conversion of estimates for θ; q, a to estimates for π; q, a and p n q, a 38 5 Estimation of ψ; q, a /φq, θ; q, a /φq, and π; q, a li/φq, for q Eplicit upper bound for eceptional zeros of quadratic Dirichlet L-functions 4 5 An upper bound for ψ; q, a /φq, including the contribution from a possible eceptional zero Eplicit upper bounds for ψ; q, a /φq and θ; q, a /φq Conversion of estimates for θ; q, a to estimates for π; q, a 53 A Appendi: Computational details 55 A Verification of bound on NT, χ 0 for principal characters χ 0 and the computation of ν 55 A Computations of worst-case error bounds for q 0 5 and for q 55 A3 Computations of the leading constants c ψ, c θ, and c π for q A4 Dominant contributions to c ψ q, c θ q, and c π q for q A5 Computations of ψ q, θ q, θ# q, π q, and 0 q for q A6 Computations of inequalities for π; q, a and p n q, a, for q 500 and very small 63 A7 Computations of error terms for ψ; q, a, θ; q, a, and π; q, a, for q 00 and very small 63 A8 Computations of uniform range of validity for error terms for ψ; q, a, θ; q, a, and π; q, a Mathematics Subject Classification Primary N3, N37, M0, M6; secondary Y35, Y40

2 A9 Computations of lower bounds for L, χ for medium-sized moduli q for Lemma 53 and Proposition 9 64 B Notation reference 65 Acknowledgments 65 References 65 INTRODUCTION AND STATEMENT OF RESULTS The Prime Number Theorem, proved independently by Hadamard [0] and de la Vallée Poussin [37] in 896, states that π = log, p p prime or, equivalently, that θ = p p prime log p and ψ = p n p prime log p, where by f g we mean that lim f/g = Quantifying these statements by deriving eplicit bounds upon the error terms π Li, θ, and ψ 3 is a central problem in multiplicative number theory see for eample Ingham [] for classical work along these lines Here, by Li we mean the function defined by Li = dt log t log Our interest in this paper is the consideration of similar questions for primes in arithmetic progressions Let us define, given positive integers a and q with gcda, q =, θ; q, a = log p and ψ; q, a = p p a mod q log p, 4 p n p n a mod q where the sums are over primes p and prime powers p n, respectively We further let π; q, a = p p a mod q 5 denote the number of primes up to that are congruent to a modulo q We are interested in upper bounds, with eplicit constants, for the analogues to equation 3, namely the error terms Li π; q, a φq, θ; q, a φq, and ψ; q, a φq 6 Such eplicit error bounds can take two shapes The first, which we will term bounds of Chebyshev-type, are upper bounds upon the error terms that are small multiples of the main term in size, for eample inequalities of the form ψ; q, a φq < δ q,a φq, 7 for small positive δ q,a and all suitably large values of The second, which we call bounds of de la Vallée Poussintype, have the feature that the upper bounds upon the error are of genuinely smaller order than the size of the main term and hence, in particular, imply the Prime Number Theorem for the corresponding arithmetic progression, something that is not true of inequality 7 Currently, there are a number of eplicit inequalities of Chebyshev-type in the literature In McCurley [6], we find such bounds for non-eceptional moduli q which is to say, those q for which the associated Dirichlet L-functions have no real zeros near s =, valid for large values of McCurley [7] contains analogous bounds in the case q = 3 Ramaré and Rumely [6] refined these arguments to obtain reasonably sharp bounds of Chebyshev-type for all q 7

3 and various larger composite q 486; the first author [3] subsequently etended these results to primes 73 q 347 Very recently, these results have been sharpened further for all moduli q 0 5 by Kadiri and Lumley [4] Bounds of de la Vallée Poussin-type are rather less common, however, other than the classical case where one considers all primes that is, when q = or, where such inequalities may be found in famous work of Rosser and Schoenfeld [8] see also [9, 3] for subsequent refinements When q 3, however, the only such result currently in the literature in eplicit form may be found in a 00 paper of Dusart [5], which treats the case q = 3 Our goal in the paper at hand is to deduce eplicit error bounds of de la Vallée Poussin-type for all moduli q 3 We prove: Theorem Let q 3 be an integer and let a be an integer that is coprime to q There eist eplicit positive constants c ψ q and ψ q such that ψ; q, a φq < c ψq log for all ψq 8 Moreover, c ψ q and ψ q satisfy c ψ q c 0 q and ψ q 0 q, where 000, if 3 q 00, c 0 q = 340, if 0 q 04, 80, if q > 04, and 4 0 9, if 3 q 6, 0 q = /q, if 7 q 0 5, ep003 q log 3 q, if q > 0 5 Similarly, for θ; q, a and π; q, a we have: 9 0 Theorem Let q 3 be an integer and let a be an integer that is coprime to q There eist eplicit positive constants c θ q and θ q such that θ; q, a φq < c θq log for all θq Moreover, c θ q c 0 q and θ q 0 q, where c 0 q and 0 q are as defined in equations 9 and 0, respectively Theorem 3 Let q 3 be an integer and let a be an integer that is coprime to q There eist eplicit positive constants c π q and π q such that Li π; q, a φq < c πq log for all π q Moreover, c π q c 0 q and π q 0 q, where c 0 q and 0 q are as defined in equations 9 and 0, respectively In each case with 3 q 0 5, eact values of c ψ q, c θ q, c π q, ψ q, θ q, and π q can be found in data files accessible at: See Appendices A3 and A5 for more details We note here that many of our results required considerable computations; the relevant computational details are available at the above URL and are discussed in Appendi A By way of eample, if we restrict our attention to the case q = 3 and a {, }, Theorem yields the inequality θ; 3, a < for all 3,460,547,899 3 log Here the constant sharpens the corresponding value 06 in Dusart [5] by a factor of roughly 500 We remark that 3,460,547,899 is the best-possible range of validity for the error bound 3; indeed this is true for each ψ q, θ q, and π q, for 3 q 0 5 For larger moduli q > 0 5, the inequalities c ψ q c 0 q, c θ q c 0 q, and c π q c 0 q 3

4 are actual equalities by our definitions of the left-hand sides, and similarly ψ q = θ q = π q = 0 q = ep003 q log 3 q, for these large moduli We note that one can obtain a significantly smaller value for 0 q if one assumes that Dirichlet L-functions modulo q have no eceptional zeros see Proposition 58, which sharpens the results of McCurley [6] mentioned above Theorems and, even if one appeals only to the inequalities c ψ q c 0 q and c θ q c 0 q, sharpen Theorem of Ramaré and Rumely [6] for q 3 and every other choice of parameter considered therein An almost immediate consequence of Theorem 3, just from applying the result for q = 3 and performing some routine computations see Appendi A7 for details, is that π Li < log for all 909,08,399 4 While, asymptotically, this result is inferior to the state of the art for this problem, it does provide some modest improvements on results in the recent literature for certain ranges of By way of eample, it provides a stronger error bound than Theorem of Trudgian [36] for all 909,08,399 < and sharpens corresponding results in [4] and [6] in much smaller ranges Eploiting the fact that Li is predictably close to / log, we can readily deduce from Theorem 3 the following two results, which are proved in Section 4 We define p n q, a to be the nth smallest prime that is congruent to a modulo q Theorem 4 Let q 3 be an integer, and let a be an integer that is coprime to q Suppose that c π qφq < Then for > π q, < π; q, a < φq log φq log log We remark that Dusart [5] proved the lower bound in Theorem 4 in the case q = 3 Theorem 5 Let q 3 be an integer, and let a be an integer that is coprime to q Suppose that c π qφq < Then either p n q, a π q or nφq lognφq < p n q, a < nφq lognφq log lognφq 6 Thanks to our computations of the constants c π q, we can produce a very eplicit version of the above two results for certain moduli q see Appendi A6 for details Corollary 6 Let q 500 be an integer, and let a be an integer that is coprime to q For all 50q, we have < π; q, a < + 5 φq log φq log log For all positive integers n such that p n q, a q, we have nφq lognφq < p n q, a < nφq lognφq log lognφq Bounds like those provided by Theorems,, and 3 are of a reasonable size for most purposes, when combined with tractable auiliary computations for the range up to 0 q For very small moduli, however, we may weaken the error bounds to produce analogous results that are easier still to use, in that they apply for smaller values of see Section A7 for the details of the computations involved Corollary 7 Let a and q be integers with q 00 and gcda, q = If 0 3, then ψ; q, a φq < 09 log θ; q, a φq < 040 log Li π; q, a φq < 053 log 4

5 Moreover, if 0 6, then ψ; q, a φq < 00 log θ; q, a φq < 004 log Li π; q, a φq < 007 log In the opposite direction, if we want bounds that are valid for all q and uniform in their statements, but are willing to permit the parameter to be very large, we have the following corollary see Appendi A8 for details of the computation We remark that for q 7 we can weaken the restriction on to ep003 q log 3 q Corollary 8 Let a and q be integers with q 3 and gcda, q = Suppose that ep0 q log 3 q Then ψ; q, a φq < 80 log, θ; q, a φq < 80 log Li π; q, a φq < 80 log The outline of the paper is as follows In Section, we derive an eplicit upper bound for ψ; q, a /φq, valid for the small moduli 3 q 0 5 In Section 3, this bound is carefully refined into a form which is suitable for eplicit calculation; we establish Theorem for these small moduli at the end of Section 35 In Section 4, we move from bounds for approimating ψ; q, a to analogous bounds for θ; q, a and π; q, a In particular, we establish Theorem for these moduli at the end of Section 4, and Theorems 3 5 for small moduli as well as Corollary 6 in Section 4 Section 5 contains upper bounds for ψ; q, a /φq, θ; q, a /φq, and π; q, a li/φq for larger moduli q > 0 5 We establish Theorems and for these large moduli in Section 53 see the remark before Corollary 57, and Theorem 3 for these moduli in Section 54 Indeed, in those sections, we also deduce a number of eplicit results with stronger error terms saving greater powers of log, as well as analogous results for an improved range of that hold under the assumption that there are no eceptional zeros for the relevant Dirichlet L-functions Finally, in Appendi A, we provide details for our eplicit computations, with links to files containing all our data We provide a summary of the notation defined throughout the paper in Appendi B Before we proceed, a few remarks on our methods are in order The error terms 3 depend fundamentally upon the distribution of the zeros of the Riemann zeta function, as evidenced by von Mangoldt s formula: ψ ε + ψ + ε lim = ε 0 ρ ρ ρ log π + log, where the sum is over the zeros ρ of the Riemann zeta function in the critical strip, in order of increasing Im ρ Deriving good approimations for ψ; q, a, θ; q, a, and π; q, a depends in a similar fashion upon understanding the distribution of the zeros of Dirichlet L-functions Note that, as is traditional in this subject, our approach takes as a starting point von Mangoldt s formula, and hence we are led to initially derive bounds for ψ; q, a, from which our estimates for θ; q, a and π; q, a follow The fundamental arguments providing the connection between zeros of Dirichlet L-functions and eplicit bounds for error terms in prime counting functions derive from classic work of Rosser and Schoenfeld [8], as etended by McCurley [6], and subsequently by Ramaré and Rumely [6] and Dusart [5] The main ingredients involved include eplicit zero-free regions for Dirichlet L-functions by Kadiri [] and McCurley [8], eplicit estimates for the zero-counting function for Dirichlet L-functions by Trudgian [35], and the results of large-scale computations of Platt [4], all of which we cite from the literature Other necessary results include lower bounds for L, χ for quadratic characters χ, upper bounds for eceptional zeros of L-functions with associated character χ, and eplicit inequalities for bχ, the constant term in the Laurent epansion of L L s, χ at s = 0 see Definition 56 below In each of these cases, our results sharpen eisting eplicit inequalities and thus might be of independent interest: 5

6 Proposition 9 If χ is a primitive quadratic character with conductor q > 6677, then L, χ > q Proposition 0 Let q 3 be an integer, and let χ be a quadratic character modulo q If β > 0 is a real number for which Lβ, χ = 0, then 40 β q log q Proposition Let q 0 5 be an integer, and let χ be a Dirichlet character mod q Then bχ 0506q log q Proposition 9 is established in Section A9 For larger values of q, we can improve on Proposition 9 by a little more than a factor of 0; see Lemma 9 for a more precise statement Propositions 0 and are established in Sections 5 and 5, respectively Throughout our work, we have made every effort to avoid specifying many of our free parameters, such as a constant R that defines the size of a zero-free region for Dirichlet L-functions and an below which one must verify certain inequalities by eplicit computation even though, at the end of the day, we do make specific choices of these parameters The reason for this is to make it easy to sharpen our bounds in the future when one has available stronger zero-free regions and more computational power The constants present in, for eample, Theorem, decrease roughly as a linear function in R We have chosen to split our small q and large q results at the modulus q = 0 5 even though Platt s calculations etend through the modulus partially due to limitations of computational time, partially because it is a convenient round number, but also partially because it is epedient for us to assume in our auiliary results that H 000, and Platt s verification of GRH is valid for H = 0 8 /q which drops below 000 once q eceeds 0 5 PREPARATION OF THE UPPER BOUND FOR ψ; q, a /φq, FOR q 0 5 In this section, we will derive our initial upper bound upon ψ; q, a /φq for small moduli q, that is, for q 0 5 This bound given as Proposition 0 will turn out to be independent of ecept for a single complicated function F q,m,r ; H, defined in Definition 6, multiplied by various powers of log Our starting point is an eisting eplicit version of the classical eplicit formula for ψ; q, a in terms of zeros of Dirichlet L-functions; by the end of this section, all dependence on the real parts of these zeros will be removed, and the dependence on their imaginary parts will be confined to the single function F q,m,r ; H Along the way, we will use as input eisting eplicit bounds for N χ T see Proposition 6 below, and we will derive an eplicit upper bound, contingent on GRH, for the sum of / β + γ over all zeros β + iγ of a given Dirichlet L-function see Lemma 7 and its variants below We mention also that the eplicit formula we use contains a parameter δ that can be chosen to be constant to obtain bounds of Chebyshev-type However, we must choose δ to be a function of that decreases to 0 in order to obtain our bounds of de la Vallée Poussin-type; we make that choice of δ in the proof of Lemma 9 and motivate our choice in the remarks following that proof Previous work based on the eplicit formula Throughout this paper, q will be a positive integer we will usually assume that q 3, and a will be a positive integer that is relatively prime to q For any Dirichlet character χ mod q, the Dirichlet L-function is defined as usual by Ls, χ = n= χn n s when Rs >, and by analytic continuation for other comple numbers s We adopt the usual convention of letting ρ = β + iγ denote a zero of Ls, χ, so that β = Rρ and γ = Iρ by definition; and we define Zχ = {ρ C: 0 < β <, Lρ, χ = 0} to be the set of zeros of Ls, χ inside the critical strip technically a multiset, since multiple zeros, if any, are included according to their multiplicity Notice in particular that the set Zχ does not include any zeros on the imaginary ais, even when χ is an imprimitive character; consequently, if χ is induced by another character χ, then Zχ = Zχ We recall, by symmetry and the functional equation for Dirichlet L-functions, that if ρ = β + iγ Zχ then also ρ = β +iγ Zχ Finally, we say such an L-function satisfies GRHH, the generalized Riemann hypothesis up to height H, if β + iγ Zχ and γ H = β = 6

7 For eample, Platt [4] has verified computationally that every Dirichlet L-function with modulus q satisfies GRH0 8 /q see [3] for more details of these computations Definition We write N χ T for the standard counting function for zeros of Ls, χ with 0 < β < and γ T In other words, N χ T = #{ρ Zχ: γ T }, counted with multiplicity if there are any multiple zeros In the literature, this function is often written as NT, χ Definition For any positive integer m and any positive real number δ, define A m δ = m m δ m + jδ m+ j We note that A m δ = j=0 m+ k=0 α m,k δ k m, where α m,0 = m and, for k m +, m + m m α m,k = j k k j We quote the following proposition from Ramaré Rumely [6, Theorem 43, p 45] The proposition, which also appears in Dusart s work [5, Theorem, pp 39 40], is a modification of McCurley s arguments [6, Theorem 36] that themselves hearken back to Rosser [7] Proposition 3 Let q be a positive integer, and let a be an integer that is coprime to q Let > and H be real numbers, let m be a positive integer, and let δ be a real number satisfying 0 < δ < m Suppose that every Dirichlet L-function with modulus q satisfies GRH Then φq ψ; q, a φq < U q,m; δ, H + mδ + V q,m; δ, H + W q, 3 where we define U q,m ; δ, H = A m δ V q,m ; δ, H = + mδ j=0 χ mod q ρ Zχ γ >H W q = φq + p q χ mod q ρ Zχ γ H β ρρ + ρ + m p β ρ log + 4 log q To offer some contet, the genesis of this upper bound is the classical eplicit formula for ψ; q, a, smoothed by m-fold integration over an interval near of length δ The term U q,m ; δ, H bounds the contribution of the large zeros to this smoothed eplicit formula in which the factor A m δ arises from some combinatorics of the multiple integration, while the term V q,m ; δ, H bounds the contribution of the small zeros The term mδ arises when recovering the original difference ψ; q, a /φq from its smoothed version Finally, Ramaré Rumely work only with primitive characters, in constrast to McCurley, to avoid the zeros of Ls, χ on the imaginary ais see [6, p 399], although their remark on [6, p 44] is easy to misconstrue This choice, which we follow as evidenced by the definition of Zχ in equation, simplifies the analytic arguments but results in a mild error on the primecounting side, which is bounded by W q In practice, we will be choosing δ so that the first term U q,m ; δ, H is almost eactly δ mδ ; for most moduli q, that term together with the quantity will provide the dominant contribution to our eventual upper bound For very small moduli q, however, it is the term V q,m ; δ, H that provides the dominant contribution We remark that the aforementioned work of Kadiri and Lumely [4] incorporates a different smoothing mechanism that is inherently more fleible than simple repeated integration; such an approach would be a promising avenue for possible sharpening of our results 7

8 In this upper bound, which is a function of for any given modulus q, the parameters m, δ, and H are at our disposal to choose We will quickly make the choice H = 0 8 /q, so that every Dirichlet L-function with modulus q satisfies GRHH by Platt s computations [4]; this choice allows for a strong bound for V q,m ; δ, H In the proof of Lemma 9, we will be motivated to choose δ equal to a specific function of and the other parameters Without some choice of δ that tended to 0 as tends to infinity, it would be impossible to achieve a de la Vallée Poussin-type bound, because of the term mδ in the upper bound; our choice, as it turns out, decays roughly like ep c log for large Finally, after the bulk of the work done to estimate the above upper bound, we will compute the resulting epression for various integer values of m and select the minimal such value It will turn out that we always choose m {8, 9, 0}, for q 0 5, although we have no theoretical eplanation for how we could have predicted these choices to be optimal in practice Upper bounds for V q,m ; δ, H, eploiting verification of GRH up to bounded height The quantities defined in equations 4 and 5 are both sums over zeros of Dirichlet L-functions, and we will require some knowledge of the distribution of those zeros That information is essentially all classical, ecept that of course we require eplicit constants in every estimate, and we can also take advantage of much more etensive modern computations Specifically, we draw information from three sources: Trudgian s work on the zeros of the Riemann ζ-function and Dirichlet L-functions with eplicit constants, Platt s computations of many zeros of Dirichlet L-functions, and direct computation using Rubinstein s lcalc program [30] We begin by a standard partial summation argument relating the inner sum in V q,m ; δ, H to the zero-counting function N χ T ; we state our result in a form that has some fleibility built in Lemma 4 Let q be a positive integer, and let χ be a character with conductor q Let H X be real numbers, and suppose that χ satisfies GRHX Then ρ γ + /4 + N χ H H N χ T + X H X T dt ρ Zχ γ H ρ Zχ γ X Proof Let χ be the primitive character mod q that induces χ, so that Zχ = Zχ For ρ Zχ, if γ X then ρ = γ + / by our assumption of GRHX; on the other hand, if γ > X, then we have the trivial bound ρ > γ Consequently, using partial summation, as required ρ Zχ γ H ρ ρ Zχ γ X = ρ Zχ γ X = ρ Zχ γ X = ρ Zχ γ X = ρ Zχ γ X γ + /4 + ρ Zχ X< γ H H γ + /4 + X γ dn χ T T γ + /4 + N χ T T H X H N χ T d X T γ + /4 + N χ H N χ X H N χ T H X X T dt γ + /4 X + N χ H H N χ T + H X T dt We turn now to eplicit bounds for the zero-counting functions N χ T, beginning with a bound when χ is the principal character 8

9 Proposition 5 If χ 0 is the principal character for any modulus q, then for T > e, T N χ 0 T π log T πe + 7 < 034 log T πT Proof We adopt the standard notation NT for the number of zeros of ζs in the critical strip whose imaginary part lies between 0 and T, as well as ST = π arg ζ + it for the normalized argument of the zeta-function on the critical line Trudgian [33, Theorem ] gives the eplicit estimate ST 07 log T for T > e It is well known that the error term in the asymptotic formula for NT is essentially controlled by ST ; for an eplicit version of this relationship, Trudgian [34, equation 5] gives T NT π log T πe π tan T + T 4π log + 4T + + ST 3πT for T In our notation, N χ0 T is eactly equal to NT since the former counts zeros lying both above and below the imaginary ais Using the inequalities tan y y and log + y y which are valid for y 0, it follows that N χ 0 T T π log T πe which is equivalent to the asserted bound 4π tan 4π T + T 4π T + T 4π log 4T + 3πT + 4T + 3πT + 07 log T log T Proposition 6 Let q be a positive integer, and let χ be a character with conductor q Then for all T, N χ T T π log q T πe < 0399 logq T Proof If χ is nonprincipal, this follows immediately from Trudgian [35, Theorem ] which sharpens McCurley [6, Theorem ] For χ principal, we have q = and the desired inequality is implied by Lemma 5, provided T 04 For T 04, we may verify the bound computationally see Appendi A, completing the proof It is worth mentioning that the main result of [35] contains a number of inequalities like equation 7, with the right-hand side replaced by various terms of the shape a logq T + b for eplicit constants a and b The one we have given here is the best available choice for small values of q T, but could be improved for larger q T ; the end result of such a modification to our proof is rather negligible Before we assemble our bound on V, we need to handle the innermost summation We provide two bounds, one based solely on Proposition 6, and a second bound for 3 q 00 using eplicit computation of zeros produced by Rubinstein s lcalc program Lemma 7 Let q be a positive integer, and let χ be a character with conductor q satisfying GRH For any H, log H ρ < logq H log H 096 log H π π ρ Zχ γ H Proof If χ is the primitive character mod q that induces χ, then Lemma 4 with X = gives ρ γ + /4 + N χ H H N χ T + H T dt ρ Zχ γ H ρ Zχ γ N χ + N χ H H H N χ T + T dt, 9 since / γ + /4 / 0 + /4 = for each summand The three terms on the right-hand side can all be bounded above by Proposition 6: N χ π log q πe logq ; N χ H H 9 π log q H πe H logq H H ;,

10 and ρ Zχ γ H H N χ T T dt H = log H π T π log q T dt πe logq T T logπe π logq H log H π 0399logq H H H The sum of these three epressions yields the upper bound log H ρ + π π logq H log H π π logπe + H H 0399, which implies the statement of the lemma thanks to the inequalities < 4347, π + logπe π log H logπe log H π > 096, and 0676 logπe < 9684 π A considerable source of error is the overestimate of the number of zeros with small height For characters with small conductor, we can estimate those zeros computationally and obtain a substantially superior bound Lemma 8 Let q be a positive integer, and let χ be a character with conductor q satisfying GRH0 3 For any H 0 3, ρ < γ + /4 0 3 ρ Zχ γ H ρ Zχ γ 0 3 ρ Zχ γ H ρ Zχ γ 0 3 logh log q H log H + 95 log H 49 π π Proof If χ is the primitive character mod q that induces χ, then Lemma 4 with X = 0 3 immediately gives ρ γ + / N χ H H N χ T + H 0 T dt 3 We preserve the first term on the right-hand side For the second and third terms, we estimate them just as in the proof of Lemma 7; we omit the straightforward details With these two estimates for the reciprocal sum over zeros of Ls, χ, we will be ready to establish two upper bounds for V q,m ; δ, H a first bound for all moduli under consideration, and a superior bound for moduli up to 00 after one more preliminary lemma Lemma 9 For any positive integer q, φq χ mod q logq H = logqh p q where the sum on the right-hand side is over the distinct prime factors p of q Proof This identity follows immediately from [7, Proposition 33] log p p, We begin by using Lemma 7 to find a bound on V q,m ; δ, H for a convenient value of H 0

11 Lemma 0 Let 3 q 0 8 and m be positive integers, and let > and 0 < δ < m be real numbers If every character with modulus q satisfies GRH0 8 /q, then V q,m ; δ, 0 8 /q < + mδ φq ν q, where V q,m ; δ, H is as defined in equation 5, and where we define ν q = log q log q π + log q π 798 p q log p p 0 Proof Assuming that all characters with modulus q satisfy GRH0 8 /q, we have β = /, and therefore V q,m ; δ, 0 8 /q = + mδ β ρ = + mδ φq χ mod q ρ Zχ γ 0 8 /q φq χ mod q ρ Zχ γ 0 8 /q The summations here can be estimated by applying Lemma 7 with H = 0 8 /q, which is at least thanks to our assumption that q 0 8 : φq ρ χ mod q ρ Zχ γ 0 8 /q < φq χ mod q [ log H = log H log H π = log H log H π π logq H log H π log H + π log H + π ρ ] 096 log H logq H φq χ mod q logqh log p p p q < log q log q log q π + π 798 log p p For small q, we can do better by using Lemma 8 and eplicit computation of zeros in place of Lemma 7 Lemma Let 3 q 00 and m be positive integers, and let > and 0 < δ < m V q,m ; δ, 0 8 /q < + mδ φq ν q, where ν q is given in Table p q be real numbers Then Proof By Platt [4], all characters with modulus at most satisfy GRH0 8 /q whereby, in this range, we have β = / We are now ready to bound V q,m ; δ, 0 8 /q V q,m ; δ, 0 8 /q = + mδ β ρ = + mδ φq χ mod q ρ Zχ γ 0 8 /q φq χ mod q ρ Zχ γ 0 8 /q ρ

12 q ν q q ν q q ν q q ν q q ν q q ν q q ν q TABLE Table of ν q for 3 q 00 Using Lemma 8, with H = 0 8 /q φq χ mod q ρ Zχ γ 0 8 /q ρ < φq = φq χ mod q ρ Zχ γ φq + φq χ mod q χ mod q χ mod q ρ Zχ γ 0 3 γ + /4 0 3 log H log H π log H π π log H π 8800 log q H + 95 log H 49 γ + /4 0 3 logqh p q + 95 log H 49 log p p For each modulus q between 3 and 00, we compute the above bound using Rubinstein s lcalc program to calculate Zχ up to height 0 3 The resulting bounds, rounded up to the nearest thousandth, are given in Table See Appendi A for more details of this computation By way of comparison, we may check that ν 3 > while ν 3 = 40888, and ν 00 > 5034 compared to ν 00 = Further estimates related to vertical distribution of zeros of Dirichlet L-functions We continue by defining certain elementary functions, which we shall use when our analysis calls for upper bounds on the zero-counting functions N χ T from the previous sections, and establishing some simple inequalities for them Definition Let u, H, and q be real numbers, all greater than or equal to Define M q H, u = u qu π log H qh πe π log logq Hu , πe

13 so that u M qh, u = qu π log π u Recall that N χ u N χ H counts the number of zeros of χ with height between H and u The following lemma is the reason we have introduced M q H, u Lemma 3 Let q be a positive integer, let χ be any character mod q, and let u H be real numbers Then N χ u N χ H < M q H, u Proof Let χ mod q be the primitive character that induces χ mod q If q H π, then u q π log u H q πe π log H, πe and we can make use of the triangle inequality: N χ u N χ H = N χ u u π log q u N χ H H πe π log q H u + πe π log q u πe H π log q H πe N χ u u π log q u πe + N χ H H π log q H πe + u π log q u πe H π log q H πe By two applications of Proposition 6, N χ u N χ H < 0399 logq u = 0399 logq u logq H + u π log q u = M q H, u M q H, u, u logq H π log q u πe H π log q H πe πe H π log q H πe since M q H, u is an increasing function of q, for all q If, on the other hand, q H < π, then q 6, and direct calculation with the si primitive characters with conductor at most 6 reveals that N χ H N χ π/q = 0 Thus, by Proposition 6, N χ u N χ H = N χ u < u π log q u πe logq u = M q H, u + H π log q H πe 0399 logq H 5338 Under the hypotheses q H < π and H, we then have N χ u N χ H < M q H, u M q H, u Lemma 4 Let q be an integer and let u H 7 be real numbers Then M q H, u < u π logqu If we assume further that qh 0 8, then Proof For the first inequality, if we set u ɛ = π π logqu M qh, u M q H, H < 3776 logqh = u logπe 0676π 0399π logq u + H log q πe + logh 0399π H, then we need to prove that ɛ > 0 First, assume that u q, so that logq u logu 3 and hence u logπe 0676π 0399π logq u u logπe 0676π 97π logu The right hand side here is increasing in u, for u > 97π/ logπe, whereby, since u H 7 and q, ɛ/h logqh 596π log H H 0676π H 3 > 0

14 If, on the other hand, we have 7 H u < q, then ɛ H 97π log q 0676π + logh H 0399π > 0, which concludes the proof of the first inequality Since M q H, H = 0798 logqh , the inequality is immediate from the assumption that qh Preliminary statement of the upper bound for ψ; q, a /φq Our remaining goal for this section is to establish Proposition 0, which is an upper bound for ψ; q, a /φq in which the dependence on has been confined to functions of a single type to be defined momentarily Building upon the work of the previous two sections, we invoke certain hypotheses on the horizontal distribution of the zeros of Dirichlet L-functions to estimate many of the terms in the upper bound of Proposition 3 We have left these hypotheses in parametric form for much of this paper, in order to facilitate the incorporation of future improvements; for our present purposes, we shall be citing work of Platt and Kadiri see Proposition 34 to confirm the hypotheses for certain values of the parameters Definition 5 Given positive real numbers H, R, and C, we say that a modulus q satisfies Hypothesis ZH, R, C if every nontrivial zero β + iγ of every Dirichlet L-function modulo q satisfies either γ H and β =, or γ > H and β R logq γ /C In other words, zeros with small imaginary part less than H in absolute value lie on the critical line, while zeros with large imaginary part lie outside an eplicit zero-free region We use Hypothesis ZH, R as a shorthand for Hypothesis ZH, R, We say that a modulus q satisfies Hypothesis Z R if every nontrivial zero β + iγ of every Dirichlet L-function modulo q satisfies β R logq ma{, γ }, ecept possibly for a single eceptional zero which, as usual, will necessarily be a real zero of an L-function corresponding to a quadratic character see [0, Sections ] We preserve C in the notation to be consistent with Ramaré Rumely [6], but we will, in fact, always choose the value C = in this paper By way of eample, we note that McCurley [6] establishes Hypothesis Z for q 0 and a closely related statement for smaller moduli, though refined estimates have been given by Kadiri [3] Definition 6 Given positive integers m, q and positive real numbers H and R, define for > and u > the functions Y q,m,r, u = ep log um+ R logqu = u m+ /R logqu F χ,m,r ; H = Y q,m,r, γ ρ Zχ γ >H F q,m,r ; H = φq = φq Note that all of these functions are strictly positive χ mod q χ mod q F χ,m,r ; H ρ Zχ γ >H Y q,m,r, γ G q,m,r ; H = φq F q,m,r ; H + logqh H m = φq logqh 4H m 4 + χ mod q F χ,m,r ; H

15 We begin with a crude upper bound on F χ,m,r ; H and F q,m,r ; H, which will suffice for this section as we dispose of some of the less significant terms in the upper bound in Proposition 3, eventually arriving at Proposition 0 in which F q,m,r ; H represents the only remaining function of in the upper bound In Section 3, we will devote significant energy to developing a more precise upper bound on F q,m,r ; H Lemma 7 If q, H 000, m 3 and >, then F χ,m,r ; H < and the same upper bound holds for F q,m,r ; H ρ Zχ γ >H ρ Zχ γ >H log qh < γ m+ H m, Proof The first inequality is obvious, since the eponential factor in Y from Definition 6 is always less than : F χ,m,r ; H = Y q,m,r, γ < γ m+ Writing as a Riemann Stieltjes integral and integrating by parts, F χ,m,r ; H < H t m+ dn χt N χ H = N χt N χ H t m+ + m + H ρ Zχ γ >H H N χ t N χ H t m+ The boundary term vanishes at both ends since N χ t N χ H t π logqt t log t by Lemmas 3 and 4, and therefore N χ t N χ H F χ,m,r ; H < m + dt < m + π = m + H H t logqt t m+ + m logqh π m H m < log qh H m as claimed, where the last inequality uses the assumptions m 3 and H 000 This upper bound for F χ,m,r ; H is certainly an upper bound for F q,m,r ; H as well, since the latter function is simply an average value of φq functions of the former type t m+ Recall Hypothesis Z from Definition 5 Here and henceforth, we will define { ν q, if 3 q 00, νq = ν q, if q > 00, where ν q is as in equation 0 and ν q is as in Table dt dt 3 Lemma 8 Let 3 q 0 5 be an integer, and let a be an integer that is coprime to q Assume that R is a positive real number such that q satisfies Hypothesis Z0 8 /q, R Let m 3 be an integer, and suppose that > e and δ are positive real numbers satisfying δ < m Then φq ψ; q, a φq A mδg q,m,r ; 08 + mδ q + W q + + mδ φq νq, where A m δ is as in Definition, G q,m,r ; H as in Definition 6, and W q as in equation 6 Proof Set H = 0 8 /q Using Proposition 3 and our bounds for V q,m ; δ, H from Lemmas 0 and, it suffices to show that χ mod q ρ Zχ γ >H β ρρ + ρ + m < 5 χ mod q ρ Zχ γ >H β γ m+ G q,m,r; H,

16 where the first inequality is obvious We begin with the observation that for any fied χ, β γ m+ + β γ m+ = β γ m+ + ρ Zχ γ >H ρ Zχ γ >H ρ Zχ γ >H ρ Zχ γ >H β γ m+ by the functional equation Of the two numbers β and β, one is at most while the other, by Hypothesis ZH, R, is at most /R logq γ Therefore ρ Zχ γ >H β γ m+ + ρ Zχ γ >H β γ m+ < ρ Zχ γ >H log qh H m + / γ m+ + ρ Zχ γ >H ρ Zχ γ >H Y q,m,r, γ, /R logq γ γ m+ by Lemma 7 noting that H 000 since q 0 5 and the definition of Y Summing both sides over all characters χ mod q yields χ mod q which is the desired inequality ρ Zχ γ >H β φq log qh < γ m+ H m + χ mod q ρ Zχ γ >H Y q,m,r, γ, Our net step is to set δ = δ 0, m, q in the previous lemma, where Gq,m,R ; 0 8 /m+ /q δ 0, m, q = 4 + φqνq / The motivation for this choice of δ comes from approimating A m δ by δ m, which is asymptotically correct for small δ, and then minimizing m G q,m,r; H + mδ δ + + mδ φq νq We have checked numerically that, in practice, this value of δ is very close to a local minimum, with any improvement being less than part in 000 Lemma 9 Let 3 q 0 5 be an integer Assume that R is a positive real number such that q satisfies Hypothesis Z0 8 /q, R Let 3 m 475 be an integer and 50 be a real number Then φq ψ; q, a φq A m δ0, m, q G q,m,r ; 08 + mδ 0, m, q +W q + + mδ 0, m, q φq νq, q 5 where δ 0, m, q is as defined in equation 4 We note that we could make our argument work for m =, but only for roughly q < 3000 Since, in practice, we will always choose m {8, 9, 0} and 0, this limitation is not of great concern to us Proof The desired inequality follows immediately from Lemma 8 upon setting δ = δ 0, m, q, as long as we can verify that this choice of δ really is less than m Set H = 08 /q Using Lemma 7, Definition 6 of G q,m,r ; H implies the inequalities φq logqh G q,m,r ; H < 4H m + logqh φq logqh H m = 4H m + 066φq log qh < H m χ mod q for 50 Since + φqνq / >, it suffices to show that G q,m,r ; H /m+ < m, 6

17 or, equivalently, that m+ m+ G q,m,r ; H < m The following quantities are increasing in q, so we can set q = 0 5, and then decreasing in H, so we can set H = 0 3 : 066 φq logqh 066 m+ H m m+ 05 log0 5 H H m 6600 log0 m m m = log0 500 We now successfully compare this epression to m+ m as follows: m m log log0 = m m+ 50 m m+ m+ m+ m m+ =, m where we have used the inequality log m m m+, which is valid for 3 m 475, and the hypothesis 50 We now use the fact that many of the quantities on the right-hand side of the inequality 5 are clearly decreasing functions for the values of we will consider We emphasize that for now, and through the proof of Corollary 340 below, the parameters R and are left unspecified to allow for easy incorporation of any future computational improvements; we apply this work with specific values of R and in the proof of Theorem at the end of Section 35 Proposition 0 Let 3 q 0 5 be an integer, and let a be an integer that is coprime to q Assume that R is a positive real number such that q satisfies Hypothesis Z0 8 /q, R Let 3 m 475 be an integer, and let e m+ be a real number Then, for all, we have φq ψ; q, a φq log m m k=0 α m,k km k+ m+ φq m+ +α m,m+ m+ m+ φq m+ +m +m φq F q,m,r F q,m,r φq ; 08 q ; 08 q F q,m,r F q,m,r ; 08 q ; 08 q log m+ log m+ k+ + log0 8 log m+ k+ m+ k+ 0 8 /q m m+ + log0 8 log m+ m+ m+ m+ 0 8 /q m + φqνq log m+ + log08 log m+ m+ 0 8 /q m + Wq log + φqνq log log /5 + log08 log /5 0 8 /q m m+ + φqνq m k m+ m+ φqνq log 5m+3 5m+5 6 Proof Set H = 0 8 /q Our starting point is Lemma 9, which yields the inequality φq ψ; q, a φq log A m δ0, m, q G q,m,r ; H log + mδ 0, m, q log + W q log + + mδ 0, m, q νq log φq 7 7

18 Note that from Definition and equation 4, A m δ0, m, q G q,m,r ; H = so that by Definition 6, m+ k=0 A m δ0, m, q G q,m,r ; H log = m+ m k=0 α m,k km m+ φq k+ m+ α m,k k m G q,m,r ; H k+/m+ + φqνq / m k/m+, F q,m,r ; H log m+ log qh log m+ k+ k+ m+ k+ + H m + φqνq / m k m+ In the range e m+ e m+/k+, the log m+/k+ / factors are decreasing functions of Additionally, the term + φqνq / m k/m+ is a nonincreasing function of ecept when k = m + Hence, using α m,k > 0 and e m+, A m Gq,m,R ; H m+ Gq,m,R ; H log m m k=0 α m,k km k+ m+ φq m+ F q,m,r ; H log m+ log qh log k+ m+ k+ k+ m+ + H m + φqνq / + α m,m+ m+ m+ φq m+ F q,m,r ; H log m+ log qh log m+ m+ m+ m+ m+ + H m Notice that in the last term we have simply bounded the factor + φqνq / m+ by For the other three terms in the sum on the right hand side of inequality 7, we argue similarly We find that mδ 0, m, q log since qh = 0 8 and, by inspection of equation 6, To bound observe that, by Definition 6 and equation 4, δ 0, m, q log φq = φq = φq m F q,m,r ; H log m+ + log08 log m+ H m F q,m,r ; H + It is a straightforward calculus eercise to check that W q log W q log + mδ 0, m, q νq log φq, log qh H m m+ + φqνq F q,m,r ; H log /5 log qh log /5 + H m log, m+ log m+ + φqνq + φqνq m+ log 5m+3 5m+5 log /5, and are monotone decreasing for > e ; therefore we can bound these functions above by log, + φqνq m+ log 5m+3 5m+5, and log /5, m+ m k m+ m+ log 5m+3 5m+5 8 respectively Together these inequalities establish the proposition 8

19 Observe that the quantity F q,m,r ; H, multiplied by a power of log, appears in four places in the bound of Proposition 0, which is otherwise independent of We have rewritten F q,m,r ; H k+ m+ log as the k+ m+th root of F q,m,r ; H log m+ k+, motivated by the fact that Fq,m,R ; H is a complicated function which we want to analyze in pieces, as it is easier to optimize each such piece separately against the power of log That optimization requires the bulk of Section 3 to carry out, culminating in Proposition 336 below In that proposition the power of log is restricted to be greater than /4, not for deep mathematical reasons but rather because a transition in the behaviour of our upper bound for the function F q,m,r ; H log η at η = /4 would require an inordinate amount of additional notation It is for this reason that we have packaged a small power of log with F q,m,r ; H in equation 8; the specific eponent /5 was chosen empirically after some numerical eploration 3 OPTIMIZATION OF THE UPPER BOUND FOR ψ; q, a /φq, FOR q 0 5 To turn Proposition 0 into something amenable to computation, we are left with the problem of deriving an absolute upper bound for F q,m,r ; H multiplied by various powers of log Indeed, we saw at the end of the last section that we need to bound F q,m,r ; H log r for r of the form m + /k +, where 0 k m + so that, in particular, if m, then 3/4 r m + Our bounds for this function consist of several pieces, each of which can be optimized using calculus; we simply add the individual maima together to deduce a uniform upper bound for F q,m,r ; H That optimization, however, can only take place once we have provided bounds for these pieces of a simpler form Our operating paradigm is that any function that can be easily programmed, and whose values can be calculated to arbitrary precision in a negligible amount of time, is suitable for our purposes, even when there remains a layer of notational compleity that we would find difficult to work with analytically Of course, our choices in this optimization section are guided by our heuristics and hindsight about which pieces of the upper bound are most significant in the end 3 Elimination of eplicit dependence on zeros of Dirichlet L-functions As we see in Definition 6, the function F q,m,r ; H still depends on the vertical distribution of zeros of Dirichlet L-functions mod q A standard partial summation argument, combined with the bounds on N χ T we established in Section 3, allows us to remove that dependence on zeros of L-functions in favor of more elementary functions Definition 3 Let q and m be positive integers, and suppose that and H are real numbers Define H q,m,r = q log ep Rm + and H q,m,r ; H = ma{h, H q,m,r } = { H, if ep m + R log qh, H q,m,r, if ep m + R log qh Lemma 3 Let m, q,, H and R be fied positive real numbers with > The function Y q,m,r, u from Definition 6 is, as a function of u, increasing for /q < u < H q,m,r and decreasing for u > H q,m,r, where H q,m,r is given in Definition 3 Proof This is a straightforward calculus eercise Proposition 33 Let m, q and H 7 Let F q,m,r ; H and Y q,m,r, u be as in Definition 6, and M q H, u as in Definition Then F q,m,r ; H M q H, H q,m,r ; HY q,m,r, H q,m,r ; H 9 + H u M qh, u Y q,m,r, u du 3

20 ; H Then, from Definition 6 and inte- Proof For this proof, write Y u for Y q,m,r, u and H for H gration by parts, F χ,m,r ; H = H Y u d N χ u N χ H q,m,r = lim N χu N χ H Y u N χ H N χ H Y H u = H where the limit equals 0 because N χ u N χ H Y u du, N χ u N χ H < M q H, u ε u +ε, H N χ u N χ H Y u du by Lemmas 3 and 4, while Y u < u m u From Lemma 3, the Y factor is negative when u < H q,m,r and positive when u > H q,m,r Therefore, using N χu N χ H < M q H, u again, and Definition 3, F χ,m,r ; H < < H H N χ u N χ H Y u du M q H, u Y u du = lim M qh, uy u M q H, H Y H u via integration by parts The limit in the last line again equals 0, yielding F χ,m,r ; H M q H, H Y H + u M qh, u Y u du H H u M qh, u Y u du Since this last integrand is positive, we may etend the lower limit of integration from H down to H and still have a valid upper bound The resulting bound for F χ,m,r ; H is again a bound for its average F q,m,r ; H The remainder of this section is devoted to finding an upper bound for the boundary term in equation 3 Other than dealing with two cases depending on the size of relative to H, this optimization is simply a matter of calculus and notation Definition 34 Given integers m and q 3 and positive constants r, H and R, define for > the functions B q,m,r ; r, H = M qh, H Y q,m,r, H log r B q,m,r ; r = qm π = 0798 logqh ep log Hm+ R logqh log r+/ ep m + log Rm + Rm + { r 3776 logqh rr logqh B q,m,r r, H = ma H m+, e log r } πq m+ H m+ Rr m + r log r+ qh Proposition 35 Let q and H 7 be real numbers satisfying qh 0 8, and assume m and 0 < r m+ If 0 < log Rm + log qh, then M q H, H q,m,r ; HY q,m,r, H q,m,r ; H log r = B q,m,r ; r, H, 3 while if log > Rm + log qh, then M q H, H q,m,r ; HY q,m,r, H q,m,r ; H log r < B q,m,r ; r 33 0

21 Proof When 0 < log Rm + log qh, we have H q,m,r ; H = H and so equation 3 follows On the other hand, when log Rm + log qh, we have H q,m,r ; H = H q,m,r H and so by Lemma 4 and M q H, H q,m,r Y q,m,r, H q,m,r ; H = H q,m,r ; H < π = πq H q,m,r = q m+ ep = q m+ ep m+ Therefore, as R Rm+ = m+, we have Rm+ log Rm + log qh q,m,r ep log ep m+ R logqh m + log ep R m + log R log Rm + q,m,r m + log M q H, H q,m,r ; HY q,m,r, H q,m,r ; H log r < qm log ep m + log log r = B q,m,r ; r π Rm + Rm + Lemma 36 Let c, c, λ, and µ be positive real numbers, and define Φu; c, c, λ, µ = c ep c log λ u log µ u Then Φu; c, c, λ, µ, as a function of u, is increasing for < u < u 0 and decreasing for u > u 0, where u 0 = ep µ λc /λ R In particular, Φu; c, c, λ, µ Φu 0 ; c, c, λ, µ = c µ eλc µ/λ for all u Proof This is a straightforward calculus eercise Lemma 37 Let u > and H 7 be real numbers, and let q 3 and m be integers Let 0 < µ m + be a real number, and suppose that qh 0 8 Then with B, B, and B as in Definition 34, we have the following inequalities: i B q,m,r u; µ, H B q,m,rµ, H; ii B q,m,r u; µ B q,m,rµ, H when log u Rm + log qh Proof Using Lemma 36, we find that B q,m,r u; µ, H = Φ u; M qh, H H m+, R logqh,, µ M µ qh, H µr logqh H m+ e Applying the upper bound for M q H, H in Lemma 4 yields µ B 3776 logqh µr logqh q,m,r u; µ, H H m+ B q,m,rµ, H e as claimed

22 Net, observe that and so the u 0 from Lemma 36 is B q,m,r u; µ = Φ u; log u 0 = µ + / m+ m+r q m π m + R, m +, m + R, µ +, // µ + = m + R m + As logqh > µ+ m+, we have Rm + log qh u 0 By the hypothesis of claim ii, we know that Φ is a decreasing function of u Thus, B q,m,r u; µ Φ = ep m + R log qh ; q m π m + R, m +, m + R, µ + πq m+ H m+ Rµ m + µ log µ+ qh B q,m,r µ, H, as claimed We have thus bounded the first term on the right-hand side of equation 3; it remains to treat the second term u M qh, u Y q,m,r, u du, 34 which is the subject of the rest of Section 3 H 3 Estimation of integrals using incomplete modified Bessel functions We follow the strategy of previous work on eplicit error bounds for prime counting functions, going back to Rosser and Schoenfeld [9], of bounding integrals with the form given in equation 34 After some well-chosen changes of variables, we use two Taylor approimations of algebraic functions to construct a bounding integral whose antiderivative we can write down eplicitly Definition 38 Given positive constants n, m, α, β and l, define an incomplete modified Bessel function of the first kind log βt n I n,m α, β; l = t m+ ep α dt log βt l Proposition 39 Let q 3 and m be integers, and let, H and R be positive real numbers If qh 0 8, then u M qh, u Y q,m,r, u du log π I,m R, q; H + π log π log I,m H R, q; H H Proof For this proof, write Y u for Y q,m,r, u If we put α = log /R and β = q, we see from Definition 6 that Y u = ep α um+ log βu Using equation, and writing logqt/π = log βt + log/π, H u M qu qh, u Y u du = log H π π u u log βu ep α du π H um+ log βu + π log π H u since u H, as required m+ ep H α du log βu m+ ep α du, log βu