1 Introduction
In combinatorics, a composition of a positive integer
$n$
is a representation of
$n$
as a sum of positive integers in which the order of the summands matters. For example,
$4=4=3+1=1+3=2+2=2+1+1=1+2+1=1+1+2=1+1+1+1$
. Thus there are 8 compositions of 4. It is easy to see that in general there are exactly
$2^{n-1}$
compositions of
$n$
. Analogously, we let
$c(n)$
denote the number of representations of
$n$
as a product of primes, in which the order of the factors matters. Thus
$c(p)=1$
,
$c(p_{1}p_{2})=2$
, and
$c(p_{1}^{2}p_{2})=3$
, and so
$c(n)$
is not a multiplicative function. In general, if
$n=p_{1}^{\unicode[STIX]{x1D708}_{1}}p_{2}^{\unicode[STIX]{x1D708}_{2}}\cdots p_{r}^{\unicode[STIX]{x1D708}_{r}}$
is the canonical factorization of
$n$
and
$\unicode[STIX]{x1D6FA}(n)=\sum _{j}\unicode[STIX]{x1D708}_{j}$
, then
$c(n)$
is the multinomial coefficient
John G. Thompson, in the course of some of his (as yet unpublished) research on the Riemann zeta function, defined the Dirichlet series
and asked for its abscissa of convergence, while Greg Simay in a private communication noted the combinatorial interpretation of the coefficients
$c(n)$
and inquired about their asymptotic mean value.
Theorem 1.1. For
$\unicode[STIX]{x1D70E}:=\Re s>1$
, and non-negative integers
$k$
, let
Let
$\unicode[STIX]{x1D70E}_{0}\approx 1.399433329$
denote the unique real number
${>}1$
such that
With
$c(n)$
and
$D(s)$
defined as in (1) and (2), we have
for
$\unicode[STIX]{x1D70E}>\unicode[STIX]{x1D70E}_{0}$
, and
$\unicode[STIX]{x1D70E}_{0}$
is the abscissa of convergence of this Dirichlet series.
Proof. Suppose that
$\unicode[STIX]{x1D70E}>1$
. By the multinomial theorem,
$$\begin{eqnarray}G_{0}(s)^{m}=\mathop{\sum }_{\substack{ n\geqslant 1 \\ \unicode[STIX]{x1D6FA}(n)=m}}\frac{c(n)}{n^{s}}.\end{eqnarray}$$
We sum this over
$m=0,1,2,\ldots$
to obtain (2). Since
$\lim _{\unicode[STIX]{x1D70E}\rightarrow \unicode[STIX]{x1D70E}_{0}+}D(\unicode[STIX]{x1D70E})=+\infty$
, and the coefficients are non-negative, it follows that
$D(s)$
does not converge at
$s=\unicode[STIX]{x1D70E}_{0}$
. Hence
$\unicode[STIX]{x1D70E}_{0}$
is the abscissa of convergence of
$D(s)$
.◻
By the Möbius inversion formula it is evident that, for
$\unicode[STIX]{x1D70E}:=\Re s>1$
,
and these series actually converge for
$\unicode[STIX]{x1D70E}>0$
if the pole and zeros of the zeta-function are suitably avoided. Moreover
$G_{1}(s)$
may be continued to the half-plane
$\unicode[STIX]{x1D70E}>0$
as a meromorphic function with poles of order 1. Since
$\unicode[STIX]{x1D701}(s)$
can be evaluated by the Euler–MacLaurin formula (or by the Riemann–Siegel formula), one can calculate
$G_{0}(s)$
without a detailed knowledge of the distribution of prime numbers. The function
$D(s)$
has poles at points where
$G_{0}(s)=1$
; the abscissæ of such points are dense in
$(1,\unicode[STIX]{x1D70E}_{0})$
. In addition,
$D(s)$
has a transcendental singularity at points of the form
$1/d$
or
$\unicode[STIX]{x1D71A}/d$
for squarefree
$d$
, where
$\unicode[STIX]{x1D71A}$
runs through the non-trivial zeros of the zeta-function. By slitting the complex plane suitably one may continue
$D(s)$
to
$\unicode[STIX]{x1D70E}>0$
, but the imaginary axis is a natural boundary for
$D(s)$
, since Landau and Walfisz [Reference Landau and Walfisz4] showed that
$\unicode[STIX]{x1D70E}=0$
is a natural boundary of
$G_{0}(s)=1-1/D(s)$
.
2 The summatory function
With
$c(n)$
as in (1), let
By a basic theorem concerning Dirichlet series (as in Montgomery and Vaughan [Reference Montgomery and Vaughan5, Theorem 1.3] or Tenenbaum [Reference Tenenbaum7, Theorem II.1.13]) it follows from Theorem 1.1 that
As
$c(n)\geqslant 0$
for all
$n$
, and
is continuous in the closed half-plane
$\unicode[STIX]{x1D70E}\geqslant \unicode[STIX]{x1D70E}_{0}$
, it follows by the Wiener–Ikehara theorem (as discussed, for example, in Montgomery and Vaughan [Reference Montgomery and Vaughan5, §8.3]) that
as
$x\rightarrow \infty$
. Here we see that the average size of
$c(n)$
is
$cn^{\unicode[STIX]{x1D70E}_{0}-1}$
, but
for all
$n$
, so
for almost all
$n$
(in the sense of natural density), where, here and in the sequel,
$\log _{k}$
denotes the
$k$
-fold iterated logarithm.
The estimate (8) can be made quantitative by using further information about
$D(s)$
. For example, suppose we put
$\unicode[STIX]{x1D6FC}:=(\log 2)/\log \,3$
. Since
$\unicode[STIX]{x1D6FC}$
is irrational, it follows that there is no
$\unicode[STIX]{x1D70F}\neq 0$
such that
$2^{\text{i}\unicode[STIX]{x1D70F}}=3^{\text{i}\unicode[STIX]{x1D70F}}=1$
. We note that Salykhov [Reference Salikhov6] has given a measure of the irrationality of
$\unicode[STIX]{x1D6FC}$
. In the improved from of Wu and Wang [Reference Wu and Wang8], it asserts that if
$\unicode[STIX]{x1D6FF}=4.117$
, then
Here
$\Vert \unicode[STIX]{x1D717}\Vert$
denotes the distance from the real number
$\unicode[STIX]{x1D717}$
to the set of integers. We note that if
$s:=\unicode[STIX]{x1D70E}+\text{i}\unicode[STIX]{x1D70F}$
with
$\unicode[STIX]{x1D70E}\geqslant 0$
, then
$$\begin{eqnarray}\displaystyle |1-G(\unicode[STIX]{x1D70E}_{0}+s)| & = & \displaystyle \biggl|\mathop{\sum }_{p}p^{-\unicode[STIX]{x1D70E}_{0}}(1-p^{-s})\biggr|\geqslant \Re \mathop{\sum }_{p}p^{-\unicode[STIX]{x1D70E}_{0}}(1-p^{-s})\nonumber\\ \displaystyle & \gg & \displaystyle \mathop{\sum }_{p}p^{-\unicode[STIX]{x1D70E}_{0}}\sin ^{2}\biggl(\frac{1}{2}\unicode[STIX]{x1D70F}\log p\biggr)\nonumber\\ \displaystyle & {\geqslant} & \displaystyle \sin ^{2}\biggl(\frac{1}{2}\unicode[STIX]{x1D70F}\log 2\biggr)+\sin ^{2}\biggl(\frac{1}{2}\unicode[STIX]{x1D70F}\log 3\biggr).\nonumber\end{eqnarray}$$
On putting
$\unicode[STIX]{x1D717}:=\unicode[STIX]{x1D70F}(\log 3)/(2\unicode[STIX]{x1D70B})$
, we find that this is
If
$\unicode[STIX]{x1D70F}_{0}<|\unicode[STIX]{x1D70F}|\ll 1$
, the above expression is
$\gg 1$
. Otherwise there is some
$q\in \mathbb{Z}^{+}$
such that
$\unicode[STIX]{x1D700}:=\min (|\unicode[STIX]{x1D717}-q|,|\unicode[STIX]{x1D717}+q|)\leqslant {\textstyle \frac{1}{2}}$
. If
$\unicode[STIX]{x1D700}\leqslant \unicode[STIX]{x1D705}/q^{\unicode[STIX]{x1D6FF}}$
with sufficiently small, absolute
$\unicode[STIX]{x1D705}>0$
then, by (9), we have
If
$\unicode[STIX]{x1D700}>\unicode[STIX]{x1D705}/q^{\unicode[STIX]{x1D6FF}}$
, then
$\Vert \unicode[STIX]{x1D717}\Vert \gg 1/q^{\unicode[STIX]{x1D6FF}}$
. Thus the right-hand side of (10) is
$\gg 1/q^{2\unicode[STIX]{x1D6FF}}$
for all
$\unicode[STIX]{x1D717}\geqslant 1$
. Since
$q\asymp \unicode[STIX]{x1D717}\asymp \unicode[STIX]{x1D70F}$
, we deduce that
Inserting this information into the effective Ikehara–Ingham–Delange Tauberian theorem of Tenenbaum [Reference Tenenbaum7, Theorem II.7.13], we find that, if
$K:=1/(4\unicode[STIX]{x1D6FF}+1)\approx 0.05724$
, then
as
$x\rightarrow \infty$
. Although the error term here is not as small as we would like, it raises the question as to whether lower bounds for linear forms in logarithms of small primes might be used to obtain a better error term, possibly even a very good error term.
Due to the many singularities of
$D(s)$
, the asymptotic estimate (12) cannot be sharpened so as to save a power of
$x$
. Indeed, since (i) the
$c(n)$
are non-negative, (ii)
$D(s)$
is regular on the real axis in the interval
$(1,\unicode[STIX]{x1D70E}_{0})$
, and (iii)
$D(s)$
has poles with abscissæ
${>}\unicode[STIX]{x1D70E}_{0}-\unicode[STIX]{x1D700}$
, it follows by the Phragmén–Landau theorem (cf. Montgomery and Vaughan [Reference Montgomery and Vaughan5, §15.1] or Tenenbaum [Reference Tenenbaum7, Theorems 1.11 and 1.12]) that
for all
$\unicode[STIX]{x1D700}>0$
.
3 A pole-free region for
$D(s)$
Tauberian theorems are easy to employ because they do not require any knowledge of an analytic continuation for the Dirichlet series involved. As pointed out by Karamata [Reference Karamata3], the corresponding drawback is that, when available, their remainder terms are necessarily weak. We now determine a pole-free region for
$D(s)$
in order to obtain, via contour integration, a better error term in the asymptotic estimate for
$C(x)$
.
For every
$\unicode[STIX]{x1D70F}\geqslant 10$
, there exists an
$X_{0}(\unicode[STIX]{x1D70F})$
such that if
$x\geqslant X_{0}(\unicode[STIX]{x1D70F})$
, then
Of course, by the prime number theorem the left-hand side is
${\sim}2x$
as
$x$
tends to infinity. The issue is how large must we make
$x$
in order to be assured of obtaining a positive fraction of this amount.
Lemma 3.1. Let
$X_{0}(\unicode[STIX]{x1D70F})$
be defined as in (14). There exists an absolute constant
$c>0$
such that
If the Riemann Hypothesis (RH) is true, then
Proof. We may clearly assume
$\unicode[STIX]{x1D70F}$
is arbitrarily large. Write
$$\begin{eqnarray}\begin{array}{@{}ll@{}}\displaystyle \mathscr{L}(z):=\exp \biggl\{\frac{(\log z)^{3/5}}{(\log _{2}z)^{1/5}}\biggr\}\quad & (z\geqslant 3),\\ \unicode[STIX]{x1D6FD}(T):=(\log T)^{-2/3}(\log _{2}T)^{-1/3}\quad & (T\geqslant 3).\end{array}\end{eqnarray}$$
Using the Vinogradov–Korobov zero-free region
$\unicode[STIX]{x1D70E}\geqslant 1-c_{1}\unicode[STIX]{x1D6FD}(T)$
for the zeta function with
$s=\unicode[STIX]{x1D70E}+i\unicode[STIX]{x1D70F}$
,
$T:=3+|\unicode[STIX]{x1D70F}|$
, where
$c_{1}$
is an absolute constant (see e.g. [Reference Ivić2, Ch. 6]), and the validity of the bound
$\unicode[STIX]{x1D701}^{\prime }(s)/\unicode[STIX]{x1D701}(s)\ll \log T$
in the same region, we may employ standard complex integration to get, for
$x\geqslant 3$
,
$\unicode[STIX]{x1D70F}\in \mathbb{R}$
,
$\log T\leqslant (\log x)^{3/2}/(\log _{2}x)^{2}$
,
where
$c_{2}$
is absolute. We omit the details which are very similar to those in the proof of [Reference Tenenbaum7, Lemma III.5.16]. This immediately implies (15).
Assuming RH, we note that
$$\begin{eqnarray}\displaystyle & & \displaystyle \mathop{\sum }_{x\leqslant n\leqslant 3x}\unicode[STIX]{x1D6EC}(n)\{1-\cos (\unicode[STIX]{x1D70F}\log n)\}\nonumber\\ \displaystyle & & \displaystyle \quad \geqslant \frac{1}{x}\mathop{\sum }_{x\leqslant n\leqslant 3x}(x-|n-2x|)\unicode[STIX]{x1D6EC}(n)\{1-\cos (\unicode[STIX]{x1D70F}\log n)\}\nonumber\\ \displaystyle & & \displaystyle \quad =x+o(x)-\frac{1}{x}\mathop{\sum }_{x\leqslant n\leqslant 3x}(x-|n-2x|)\unicode[STIX]{x1D6EC}(n)\cos (\unicode[STIX]{x1D70F}\log n).\nonumber\end{eqnarray}$$
Thus it suffices to show that
The above sum is
$$\begin{eqnarray}\displaystyle & & \displaystyle =\frac{-1}{2\unicode[STIX]{x1D70B}\text{i}}\int _{2-\text{i}\infty }^{2+\text{i}\infty }\frac{\unicode[STIX]{x1D701}^{\prime }}{\unicode[STIX]{x1D701}}(s+\text{i}\unicode[STIX]{x1D70F})\frac{(3x)^{s+1}-2(2x)^{s+1}+x^{s+1}}{s(s+1)}\,ds\nonumber\\ \displaystyle & & \displaystyle \ll \frac{x^{2}}{\unicode[STIX]{x1D70F}^{2}}+\mathop{\sum }_{\unicode[STIX]{x1D6FE}}\frac{x^{\unicode[STIX]{x1D6FD}+1}}{(\unicode[STIX]{x1D6FE}-\unicode[STIX]{x1D70F})^{2}+1}.\nonumber\end{eqnarray}$$
Assuming RH, the sum over
$\unicode[STIX]{x1D6FE}$
is
$\ll x^{3/2}\log \unicode[STIX]{x1D70F}$
since the number of
$\unicode[STIX]{x1D6FE}$
in
$[u,u+1]$
is
$\ll \log (3+|u|)$
. The desired bound follows for all
$x\geqslant X_{0}(\unicode[STIX]{x1D70F})$
by selecting
$X_{0}(\unicode[STIX]{x1D70F})=C(\log \unicode[STIX]{x1D70F})^{2}$
with
$C$
a sufficiently large absolute constant.◻
Theorem 3.2. There exists an absolute constant
$c>0$
such that if
$\unicode[STIX]{x1D70F}\geqslant 10$
, then
uniformly for
Assuming RH, if
$\unicode[STIX]{x1D70F}\geqslant 10$
, then
uniformly for
Proof. Since
$\sum _{p}p^{-\unicode[STIX]{x1D70E}_{0}}=1$
, it follows that
$$\begin{eqnarray}\displaystyle 1-\Re \mathop{\sum }_{p}p^{-s} & = & \displaystyle \mathop{\sum }_{p}p^{-\unicode[STIX]{x1D70E}_{0}}-\mathop{\sum }_{p}p^{-\unicode[STIX]{x1D70E}}\cos (t\log p)\nonumber\\ \displaystyle & = & \displaystyle \mathop{\sum }_{p}p^{-\unicode[STIX]{x1D70E}}\{1-\cos (\unicode[STIX]{x1D70F}\log p)\}-\mathop{\sum }_{p}(p^{-\unicode[STIX]{x1D70E}}-p^{-\unicode[STIX]{x1D70E}_{0}}).\end{eqnarray}$$
Here the last term is
Now
$G_{1}$
is a decreasing function, and
$G_{1}(1.3)=2.475639<2.5$
, so if
$\unicode[STIX]{x1D70E}\geqslant 1.3$
, then the above is
${\leqslant}\frac{5}{2}(\unicode[STIX]{x1D70E}_{0}-\unicode[STIX]{x1D70E})$
. Let
$X_{0}(\unicode[STIX]{x1D70F})$
be defined as in (14). Then the expression (22) is
We now suppose that
$\unicode[STIX]{x1D70E}$
is close enough to
$\unicode[STIX]{x1D70E}_{0}$
to ensure that the second term above is at most one half the first term. For such
$\unicode[STIX]{x1D70E}$
we have a lower bound for the real part of
$1-G_{0}(s)=1/D(s)$
, and hence an upper bound for
$|D(s)|$
. To complete the argument it suffices to substitute the upper bounds for
$X_{0}(\unicode[STIX]{x1D70F})$
derived in Lemma 3.1.◻
4 A stronger quantitative mean value theorem
Theorem 4.1. There is a constant
$a>0$
such that
uniformly for
$x\geqslant 3$
. Assuming RH,
Proof. Let
We move the path of integration to form a rectilinear contour
$\mathscr{C}$
from
$\unicode[STIX]{x1D70E}_{0}-\text{i}\infty$
to
$\unicode[STIX]{x1D70E}_{0}-\text{i}T$
to
$\unicode[STIX]{x1D70E}_{0}-\unicode[STIX]{x1D6FF}-\text{i}T$
to
$\unicode[STIX]{x1D70E}_{0}-\unicode[STIX]{x1D6FF}+\text{i}T$
to
$\unicode[STIX]{x1D70E}_{0}+\text{i}T$
to
$\unicode[STIX]{x1D70E}_{0}+\text{i}\infty$
where
Thus the above is
We note that
We take
$a$
to be small enough to ensure that the right-hand side above is
${<}\frac{1}{3}\log _{2}x$
. Thus
Consequently,
$$\begin{eqnarray}\displaystyle & \displaystyle \int _{\unicode[STIX]{x1D70E}_{0}-\unicode[STIX]{x1D6FF}-\text{i}T}^{\unicode[STIX]{x1D70E}_{0}-\unicode[STIX]{x1D6FF}+\text{i}T}\ll x^{\unicode[STIX]{x1D70E}_{0}+1}\exp \{-(\log x)^{2/3}\}(\log x)^{1/3}, & \displaystyle \nonumber\\ \displaystyle & \displaystyle \int _{\unicode[STIX]{x1D70E}_{0}-\unicode[STIX]{x1D6FF}+\text{i}T}^{\unicode[STIX]{x1D70E}_{0}+\text{i}T}\ll \frac{x^{\unicode[STIX]{x1D70E}_{0}+1}}{T^{2}},\qquad \int _{\unicode[STIX]{x1D70E}_{0}+\text{i}T}^{\unicode[STIX]{x1D70E}_{0}+\text{i}\infty }\ll \frac{x^{\unicode[STIX]{x1D70E}_{0}+1}}{T^{2}}(\log x)^{1/3}. & \displaystyle \nonumber\end{eqnarray}$$
Thus
Since
$c(n)\geqslant 0$
for all
$n$
, it follows that
We note that
Hence
We take
$h=x\exp (-\frac{1}{2}b(\log _{2}x)^{3/2}(\log _{3}x)^{-2})$
to obtain the upper bound part of our estimate. The lower bound is derived similarly.
When we assume RH we argue similarly, with
Thus
$x^{-\unicode[STIX]{x1D6FF}}=\exp ((\log x)^{-5/9})=1/T$
.◻
5 Large values of
$c(n)$
From (12) or Theorem 4.1 it follows that
$c(n)=o(n^{\unicode[STIX]{x1D70E}_{0}})$
as
$n\rightarrow \infty$
. To complete our discussion we show that there are
$n$
for which
$c(n)$
is nearly this large.
Theorem 5.1. Let
$\unicode[STIX]{x1D70E}_{0}$
be defined as in (4). There is an absolute constant
$C>0$
such that
for infinitely many integers
$n$
.
From the above it is evident that
Here the order of magnitude is more precise than in (13), but the above does not guarantee that the error becomes large in both signs.
Numerical studies indicate that large values of
$c(n)$
occur only when
$n$
is quite large. For example, the least
$n$
for which
$c(n)>n$
is
for which
Proof. Let
$y$
be a large real number, and consider
$n=\prod _{p\leqslant y}p^{\unicode[STIX]{x1D708}(p)}$
where the
$\unicode[STIX]{x1D708}(p)$
are of the form
Here
$\unicode[STIX]{x1D70E}>1$
is chosen later, but we assume throughout our calculations that
$\unicode[STIX]{x1D70E}$
is bounded. For brevity we put
First we note that
Similarly,
Hence
We appeal to Stirling’s formula in the crude form
Thus
$$\begin{eqnarray}\displaystyle \log \{\unicode[STIX]{x1D6FA}(n)!\} & = & \displaystyle \biggl\{y^{\unicode[STIX]{x1D70E}}S_{0}(\unicode[STIX]{x1D70E})+O\biggl(\frac{y}{\log y}\biggr)\biggr\}\biggl\{\unicode[STIX]{x1D70E}\log y+\log S_{0}(\unicode[STIX]{x1D70E})+O\biggl(\frac{y^{1-\unicode[STIX]{x1D70E}}}{\log y}\biggr)\biggr\}\nonumber\\ \displaystyle & & \displaystyle -\,y^{\unicode[STIX]{x1D70E}}S_{0}(\unicode[STIX]{x1D70E})+O\biggl(\frac{y}{\log y}\biggr)\nonumber\\ \displaystyle & = & \displaystyle \unicode[STIX]{x1D70E}S_{0}(\unicode[STIX]{x1D70E})y^{\unicode[STIX]{x1D70E}}\log y+S_{0}(\unicode[STIX]{x1D70E})y^{\unicode[STIX]{x1D70E}}\log S_{0}(\unicode[STIX]{x1D70E})-S_{0}(\unicode[STIX]{x1D70E})y^{\unicode[STIX]{x1D70E}}+O(y).\end{eqnarray}$$
From (31) we also find that
$$\begin{eqnarray}\displaystyle \log \{\unicode[STIX]{x1D708}(p)!\} & = & \displaystyle \biggl\{\frac{y^{\unicode[STIX]{x1D70E}}}{p^{\unicode[STIX]{x1D70E}}}+O(1)\biggr\}\biggl\{\unicode[STIX]{x1D70E}\log \biggl(\frac{y}{p}\biggr)+O\biggl(\frac{p^{\unicode[STIX]{x1D70E}}}{y^{\unicode[STIX]{x1D70E}}}\biggr)\biggr\}-\frac{y^{\unicode[STIX]{x1D70E}}}{p^{\unicode[STIX]{x1D70E}}}+O\biggl(\log \frac{2y}{p}\biggr)\nonumber\\ \displaystyle & = & \displaystyle \frac{\unicode[STIX]{x1D70E}y^{\unicode[STIX]{x1D70E}}\log (y/p)}{p^{\unicode[STIX]{x1D70E}}}-\frac{y^{\unicode[STIX]{x1D70E}}}{p^{\unicode[STIX]{x1D70E}}}+O\biggl(\log \frac{2y}{p}\biggr).\nonumber\end{eqnarray}$$
Hence
On combining this with (32) we deduce that
From (28) we deduce that this is
where
Since
$S_{k}^{\prime }(\unicode[STIX]{x1D70E})=-S_{k+1}(\unicode[STIX]{x1D70E})$
, we find that
$$\begin{eqnarray}\displaystyle f^{\prime }(\unicode[STIX]{x1D70E}) & = & \displaystyle 1-\log S_{0}(\unicode[STIX]{x1D70E})-1+\frac{S_{0}(\unicode[STIX]{x1D70E})S_{2}(\unicode[STIX]{x1D70E})}{S_{1}(\unicode[STIX]{x1D70E})^{2}}\log S_{0}(\unicode[STIX]{x1D70E})\nonumber\\ \displaystyle & = & \displaystyle \biggl(\frac{S_{0}(\unicode[STIX]{x1D70E})S_{2}(\unicode[STIX]{x1D70E})}{S_{1}(\unicode[STIX]{x1D70E})^{2}}-1\biggr)\log S_{0}(\unicode[STIX]{x1D70E}).\nonumber\end{eqnarray}$$
Now
$S_{1}(\unicode[STIX]{x1D70E})^{2}\leqslant S_{0}(\unicode[STIX]{x1D70E})S_{2}(\unicode[STIX]{x1D70E})$
by Cauchy’s inequality, and indeed
$S_{1}(\unicode[STIX]{x1D70E})^{2}<S_{0}(\unicode[STIX]{x1D70E})S_{2}(\unicode[STIX]{x1D70E})$
since the summands are not proportional. Thus
$f^{\prime }(\unicode[STIX]{x1D70E})>0$
when
$S_{0}(\unicode[STIX]{x1D70E})>1$
and
$f^{\prime }(\unicode[STIX]{x1D70E})<0$
when
$S_{0}(\unicode[STIX]{x1D70E})<1$
. Hence
$f(\unicode[STIX]{x1D70E})$
is maximized by choosing
$\unicode[STIX]{x1D70E}$
so that
$S_{0}(\unicode[STIX]{x1D70E})=1$
. This
$\unicode[STIX]{x1D70E}$
depends on
$y$
, but is very close to
$\unicode[STIX]{x1D70E}_{0}$
. We obtain the same result—and the estimates are easier—if we take the slightly inferior choice
$\unicode[STIX]{x1D70E}=\unicode[STIX]{x1D70E}_{0}$
. Since
it follows that
We now obtain the stated result by taking
$\unicode[STIX]{x1D70E}=\unicode[STIX]{x1D70E}_{0}$
in (33).◻
Although the choice
$\unicode[STIX]{x1D70E}=\unicode[STIX]{x1D70E}_{0}$
suffices asymptotically, it may still be the case in numerical studies with
$y$
of moderate size that better results are obtained by choosing
$\unicode[STIX]{x1D70E}$
so that
$S_{0}(\unicode[STIX]{x1D70E})=1$
.
The estimate (28) can be refined: with more care it can be shown that
However, this does not seem to lead to a sharper result in the present context.