The size of πβπ at a fixed
base
The preceding bounds leave an error of order π in the logarithm. We now determine its
linear, logarithmic and constant terms. There are two calculations:
compare every node configuration with the first π nodes, then multiply the resulting
moment weights. The first is independent of the particular weights under
the hypotheses below; the second produces the Lambert-value correction
for πβπ. Throughout this
subsection 0 <π <1 is fixed,
and we write
π=(π;π)β,M(π)=βπβ₯1(1βππ)βπ,πΏ=πΉ(1/π)=βπβ₯1ππ1βππ,
with π΅π =βπ<ππ2 and πΆπ =βπ<πππ =(π!)2(π +1)!/2π as
above. Both products converge to positive finite values, and the
abbreviation πΏ is used only in this
subsection; it is the value at π of
the divisor series L of
SectionΒ 10.
The estimates in this subsection concern real determinants at fixed
π. Denominator valuations for the
2004 forms were treated separately in SectionΒ 2.
For ππ =1, the moments are (1 βππ+1)β1, so Cauchyβs determinant
formula gives the model calculation
π·(0)π:=det(11βππ+π+1)0β€π,π<π=ππ΅πΞπβ0β€π,π<π(1βππ+π+1),Ξπ=πβ1βπ=1(1βππ)2(πβπ).(15)
The denominator has limiting multiplicity π at 1 βππ, while
Ξππ2π=βπ<π(1βππ)β2πβπβ₯π(1βππ)β2πβΆM(π)2.
The convergence follows from
βπ|logβ‘(1 βππ)| <β and
the bound ππ(πππ) for the
logarithm of the second product. Hence
π·(0)πβΌM(π)3ππ΅ππ2π.(16)
Two copies of M(π) come from the squared
Vandermonde and one from the Cauchy denominator. We shall compare
general weights with this exact model. Fixed displacements of the last
few node indices have weight ratios tending to 1; the polynomial bound in the theorem
controls the sum over all displacements.
Theorem 3.6 (determinants of geometric moments).
Let ππ >0 satisfy, for fixed
constants πΆ and π
,
ππ+βππβΆ1(πββ)Β for each fixedΒ ββ₯0,ππ+βππβ€πΆ(1+β)π
(π,ββ₯0).
Put ππ =βπβ₯0πππ(π+1)π and π·π =det(ππ+π)0β€π,π<π. Then
π·πβΌM(π)3ππ΅ππ2ππβ1βπ=0ππ(πββ).
Proof of TheoremΒ 3.6.
We may increase πΆ to at least 1 and π
to at least 0. The bound πβ β€πΆπ0(1 +β)π
makes βπβ₯0πππππΏππ a finite
positive measure with moments ππ.
Apply CauchyβBinet to finite truncations. Entrywise convergence gives
convergence of their determinants, and monotone convergence applies to
the nonnegative tuple sums. Hence
π·π=β0β€π0<β―<ππβ1βπ<πππππππβπ<π<π(πππβπππ)2.(17)
The tuple ππ =π contributes ππ΅πΞπβπ<πππ. It does
not account for the whole asymptotic. For example, replace (0,1,β¦,π β2,π β1) by (0,1,β¦,π β2,π). The ratio of its
contribution to that of the original tuple is
πππππβ1πβ1βπ=1(1βππ+11βππ)2=πππππβ1(1βππ1βπ)2βΆπ(1βπ)2.(18)
Only pairs involving the moved node change, so
the product telescopes. At π =1/2
the limiting ratio is 2: the
displaced tuple contributes twice as much asymptotically as (0,1,β¦,π β1).
Subtracting (0,1,β¦,π β1)
gives the shifts ππ =ππ βπ.
Since ππ+1 β₯ππ +1, these are
nonnegative and weakly increasing; reversing their order gives a
partition. For π =4, for example,
(0,1,3,5)ξ½β΅β΅ξ
β΅β΅ξΎ(ππ)βΌ(0,0,1,2)ξ½β΅β΅ξ
β΅β΅ξΎ(ππ)βΌ(2,1,0,0)ξ½β΅β΅ξ
β΅β΅ξΎ(ππ).
In general, put
ππ =ππβπ for 1 β€π β€π, so that π1 β₯β― β₯ππ β₯0. Division by
the contribution of ππ =π turns the
summand for this tuple into
ππ(π)=πβπβ€π(2πβ1)ππβπβ€πππβπ+ππππβπβ
β1β€π<πβ€π(1βππβπ+ππβππ1βππβπ)2.(19)
We order the indices increasingly, so their
geometric nodes decrease towards zero. The πth index from the right contributes once
through πππ and twice through
each of its π β1 pairs with later
indices in the squared Vandermonde. Its displacement therefore has
exponent 1 +2(π β1) =2π β1.
Set ππ(π) =0 when π has more than π positive parts, so that every sum has
the same index set. Fix a partition with β positive parts. It affects only the
last β indices of (0,1,β¦,π β1): for instance, (2,1) affects the last two for every
π β₯2. Those indices tend to
infinity, so the weight ratios in (19)
tend to 1 by the fixed-shift
hypothesis.
The Vandermonde quotient still contains a growing number of factors,
but its dependence on π telescopes.
The factors with π <π β€β are
fixed, and those with π >β
equal 1. For π β₯β and each π β€β, the remaining factors satisfy
πβπ=β+11βππβπ+ππ1βππβπ=(ππ+1βπ;π)ππ(πβ+1βπ;π)ππβΆ1(πβ+1βπ;π)ππ.(20)
Thus each ππ(π) has the same limit as for
constant weights.
To pass this limit through the sum over partitions, we need a bound
independent of π that remains
summable as both the number and size of the parts grow. Only factors
with π β€β can differ from 1. For each such π, the denominator product in the
Vandermonde quotient is at least π,
and its numerator product is at most 1. Together with the polynomial shift
bound, this gives
ππ(π)β€(πΆπβ2)βββπ=1(1+ππ)π
π(2πβ1)ππ.
Let π΄ =βπ’β₯1(1 +π’)π
ππ’β1 <β. Dropping the ordering of the positive parts,
we obtain
ββ(π)=βsupπππ(π)β€(πΆπβ2)βββπ=1βπ’β₯1(1+π’)π
π(2πβ1)π’β€(πΆπβ2π΄)βπβ2.
Here π(2πβ1)π’ β€π2πβ1ππ’β1 and 1 +3 +β― +(2β β1) =β2. Since π <1, the bound (πΆπβ2π΄)βπβ2 is summable
over β, independently of π.
Dominated convergence gives the same limiting partition sum as for
ππ =1. By (15), this
sum tends to M(π), the
reciprocal of the limiting Cauchy denominator. Together with Ξπ/π2π βM(π)2, it
accounts for all three copies of M(π). Restoring the principal tuple gives
π·ππ·(0)πβπ<πππβΆ1.
The model asymptotic (16)
now proves the theorem. The comparison explains why the normalised limit
retains the geometric spacing of the nodes but loses the individual
weights.Β β»
We now return to ππ =π4πΎπ. A factorisation into two
finite Gaussian multinomial sums supplies both the shift hypotheses and
the first correction to ππ/ππ.
The shift hypotheses allow application of the determinant theorem; the
correction determines the product βπ<πππ.
Proposition 3.7 (the moment weights as a product of
two finite sums). For π β₯1
and π β₯0 set
π
(π)π(π)=βπ1+β―+ππ=π(π;π)πβπβ€π(π;π)ππ,
the sum being over
compositions into π nonnegative
parts, so that π
(π)π is the sum
of all π-multinomial coefficients
of degree π and length π, equivalently the multivariate
RogersβSzegΕ polynomial at unit argumentsΒ [42]. These are polynomial identities in the
indeterminate π; at a prime power,
π =2 gives the Galois number
counting the subspaces of a π-dimensional space over the field of
that order, and π =3 counts the
flags π1 βπ2 in that
space. The source is cited for those identifications; the factorisation
below is proved here. Then
πΎπ(π)=π
(2)π(π)π
(3)π(π)(π;π)π.
The product π
(2)ππ
(3)π lies
in β€β₯0[π], has degree βπ2/4β +βπ2/3β, and has coefficient sum 6π. Consequently (π;π)ππΎπ ββ€[π], and the weights
ππ =π4πΎπ satisfy
π5ππβ€ππβ€πβ1ππ(πβ₯0),
together with ππ+β/ππ β€πβ6(1 +β)3, so the
hypotheses of TheoremΒ 3.6
hold with π
=3.
Proof. Write E(π§) =(π§;π)β1β. Eulerβs expansion gives E(π§)π =βπβ₯0π
(π)ππ§π/(π;π)π. We seek the product
π
(2)ππ
(3)π at a common
index π. Splitting π =π +π keeps π
(3)π fixed while the sum over ((π;π)π(π;π)π)β1 supplies π
(2)π/(π;π)π. The operator πππ(π§) =(π(π§) βπ(ππ§))/π§ produces
the needed index shift: πππ§π =(1 βππ)π§πβ1 cancels the
last factor of (π;π)π. Thus
ππE=E,πππE(π§)3=βπβ₯0π
(3)π+π(π;π)ππ§π.
The π-Leibniz rule, obtained by induction
from the product rule, is
πππ(ππ)(π§)=πβπ=0[ππ]π(ππππ)(ππβππ§)(ππβπππ)(π§).
Write βπ(π§) =βππ=0[ππ]π(π§;π)π.
Applying the rule to E β
E and using E(πππ§) =(π§;π)πE(π§) gives πππE(π§)2 =E(π§)2βπ(π§); at π =1 both
sides are E(π§)2(2 βπ§).
Applying it once more, to E β
E2, gives
πππE(π§)3=E(π§)3πβπ=0[ππ]π(π§;π)πβπ(π§),
Hence
βπβ₯0π€π(π;π)ππππE(π§)3=E(π€)E(π§)3βπβ₯0π€π(π§;π)π(π;π)πβπ(π§).
The left-hand
side is βπ,πβ₯0π
(3)π+ππ€ππ§π/((π;π)π(π;π)π).
On setting π§ =π€, all pairs with the
same sum contribute to a single coefficient. Thus the coefficient of
π€π is π
(3)πβπβ€π((π;π)π(π;π)πβπ)β1 =π
(2)ππ
(3)π/(π;π)π. On the right-hand side, put π§ =π€, expand βπ, and write π =π +π‘. The identity (π€;π)π+π‘ =(π€;π)π(π€ππ;π)π‘ gives
E(π€)4βπβ₯0π€π(π€;π)2π(π;π)πβπ‘β₯0(π€ππ;π)π‘(π;π)π‘π€π‘.
The inner sum
is evaluated by the π-binomial
theorem:
βπ‘β₯0(π€ππ;π)π‘(π;π)π‘π€π‘=(π€2ππ;π)β(π€;π)β.
Finally, (π€;π)π =(π€;π)β/(π€ππ;π)β
cancels two of the five factors E(π€), leaving
E(π€)3βπβ₯0π€π(π;π)π(π€2ππ;π)β(π€ππ;π)2β=πΊπ(π€).
Every step holds coefficientwise, and absolutely for
|π€| small.
A π-multinomial coefficient has
nonnegative integral coefficients and degree βπ<πππππ. At π =1 it equals the corresponding ordinary
multinomial coefficient. Summing over compositions gives coefficient sum
ππ for π
(π)π and hence 6π for the product. The degree is
maximised by balanced parts, which gives βπ2/4β for π =2 and βπ2/3β for π =3; no leading cancellation occurs
because all coefficients are nonnegative.
For the bounds, write π(π)π =[π§π](π§;π)βπβ, so that
π
(π)π =(π;π)ππ(π)π and
πΎπ =(π;π)ππ(2)ππ(3)π.
Eulerβs expansion exhibits π(π)π as an π-fold convolution of the sequence 1/(π;π)π, whose terms lie in [1,πβ1]; the number of compositions of
π into π nonnegative parts is (π+πβ1πβ1), so
(π+πβ1πβ1)β€π(π)πβ€πβπ(π+πβ1πβ1).
Combining these at π =2,3 with π β€(π;π)π β€1 and ππ =(π +1)(π+22) gives π5ππ β€ππ β€πβ1ππ. The shift
bound follows from ππ+β/ππ β€(1 +β)3.
It remains to show that ππ+β/ππ β1 for each fixed β. The two convolutions grow at least
linearly and quadratically; it suffices to bound their increments by a
constant and a linear function, respectively. Put ππ =1/(π;π)π. Since (π;π)π+1 =(1 βππ+1)(π;π)π, these
terms satisfy 1 =π0 β€π1 β€β― β€πβ1, and the convolutions read π(2)π =βπβ€πππππβπ and
π(3)π =βπβ€ππππ(2)πβπ. Hence
π(2)π+1βπ(2)π=βπβ€πππ(ππ+1βπβππβπ)+ππ+1,π(3)π+1βπ(3)π=βπβ€πππ(π(2)π+1βπβπ(2)πβπ)+ππ+1.
Every bracket is nonnegative. In the first line
the brackets add up to ππ+1 β1 <πβ1 and each ππ β€πβ1, so 0 β€π(2)π+1 βπ(2)π β€πβ2 +πβ1 β€2πβ2. Inserting this bound into the π +1 terms of the second line gives 0 β€π(3)π+1 βπ(3)π β€2(π +1)πβ3 +πβ1 β€3(π +1)πβ3.
Against the lower bounds π(2)π β₯π +1 and π(3)π β₯(π +1)(π +2)/2 displayed above,
the ratio π(π)π+1/π(π)π
lies in [1,1 +2πβ2/(π +1)] for
π =2 and in [1,1 +6πβ3/(π +2)] for π =3. Both ratios tend to 1, and therefore so does π(π)π+β/π(π)π, a product of
β consecutive ratios. Finally (π;π)π+β/(π;π)π =ββπ=1(1 βππ+π) β1,
and ππ+β/ππ =πΎπ+β/πΎπ is the
product of these three ratios. The Lean proof obtains the fixed-shift
limit by the same steps.Β β»
The factorisation proves positivity at each fixed real π, not coefficientwise positivity in
π: for example, π0 =π4 has π-coefficient β4. The cubic factors ππ multiply to πΆπ, but ππ/ππ β1 alone would not determine the
product of π weights. The first
relative correction has order 1/π
and accumulates as a harmonic sum; the remaining ππ(πβ2) terms determine a finite
constant. The two convolutions supply this first correction
explicitly.
Theorem 3.8 (the size of πβπ at a fixed base). Let πΎπΈ be Eulerβs constant. The
product
A(π)=πβ8πΎπΈπΏβπβ₯0(πππππ8πΏ/(π+1))
converges to a positive finite value. With πΎ(π) =A(π)M(π)3,
πβπ(π)βΌπΎ(π)πΆπππ΅ππ2ππβ8πΏ,
equivalently
logβ‘πβπ(π)=π΅πlogβ‘π+logβ‘πΆπ+2πlogβ‘πβ8πΏlogβ‘π+logβ‘πΎ(π)+π(1).
Proof of TheoremΒ 3.8.
Put πΏ =πΉ(1/π), ππ =(π;π)β1π and πΈ =πβ1. Replacing every ππ by πΈ in the two convolution sums gives πΈ2(π +1) and πΈ3(π+22). The factors π +1 and (π+22) count the nonnegative
pairs and triples with sum π,
respectively. For a fixed index π,
there is one pair (π,π βπ) and there
are π βπ +1 triples whose first
coordinate is π. Thus a summable
deficit in the sequence ππ
produces a constant correction in the two-fold sum and a linear
correction in the three-fold sum. We calculate that deficit and bound
terms with two or more deficits.
Writing ππ =πΈ βππ, the identity
ππ+1 βππ =ππ+1ππ+1 gives
ππ=βπ>πππππ,0β€ππβ€πΈππ+11βπ.
Thus βππ and βπππ converge. To evaluate the
former, use Eulerβs identity E(π§) =βπβ₯0πππ§π =(π§;π)β1β for |π§| <1. Tonelliβs theorem and
logarithmic differentiation at π§ =π
give
π:=βπβ₯0ππ=βπβ₯1πππππ=πEβ²(π)=πΈβπβ₯1ππ1βππ=πΈπΏ.
In the double sum, ππππ occurs once for each π =0,β¦,π β1, which explains the factor
π. Logarithmic differentiation is
justified by locally uniform convergence on |π§| <1.
Let π(π)π =[π§π]E(π§)π. The finite-sum formula in PropositionΒ 3.7
becomes
ππ=π4(π;π)ππ(2)ππ(3)π.
We expand each finite convolution using
ππ =πΈ βππ. For two factors,
π(2)π=πΈ2(π+1)β2πΈπβπ=0ππ+πβπ=0ππππβπ=πΈ2(π+1β2πΏ)+ππ((π+1)ππ).
For three factors, the terms with
one ππ contribute β3πΈ2βππ=0(π βπ +1)ππ. Here the sum
equals (π +1)π +ππ(1) because βπππ <β. Terms with two
ππ are bounded by 3πΈπ2, and the term with three is ππ((π +1)2ππ). It follows that
π(3)π=πΈ3((π+22)β3πΏ(π+1))+ππ(1).
Dividing the correction β3πΏ(π +1) by
the leading term (π+22)
gives the relative correction β6πΏ/(π +2). The two-fold convolution
contributes β2πΏ/(π +1) in the same
way. Since (π;π)π =π(1 +ππ(ππ)),
the leading constants in ππ =π4(π;π)ππ(2)ππ(3)π multiply to π5πΈ5 =1. Replacing (π +2)β1 by (π +1)β1 changes only the ππ((π +1)β2) error. Thus the
coefficients 2πΏ and 6πΏ add, and with ππ =(π +1)2(π +2)/2 we obtain
ππππ=(1β2πΏπ+1+ππ(ππ))(1β6πΏπ+2+ππ((π+1)β2))(1+ππ(ππ))=1β8πΏπ+1+ππ((π+1)β2).
In particular ππ+1/ππ β1, and multiplication of
finitely many adjacent ratios gives ππ+β/ππ β1 for each fixed β. The uniform shift bound follows from
PropositionΒ 3.7,
so TheoremΒ 3.6
applies.
For large π, taking logarithms
preserves the first-order term and leaves the summable error logβ‘(ππ/ππ) +8πΏ/(π +1) =ππ((π +1)β2).
The finitely many earlier factors are positive, so the regularised
product also has a positive finite limit: 0 <A(π) <β. Since βπ<πππ =πΆπ, the harmonic
contribution separates as
logβ‘βπ<ππππΆπ=β8πΏβπ<π1π+1+βπ<π(logβ‘ππππ+8πΏπ+1).
The second sum tends to logβ‘A(π) +8πΎπΈπΏ. Using βπ<π(π +1)β1 =logβ‘π +πΎπΈ +π(1), the Euler-constant terms cancel, giving
βπ<πππ βΌA(π)πΆππβ8πΏ. TheoremΒ 3.6
now gives the stated determinant asymptotic. Taking logarithms yields
logβ‘πβπ=π΅πlogβ‘π+logβ‘πΆπ+2πlogβ‘πβ8πΏlogβ‘π+logβ‘πΎ(π)+π(1).
β»
The convolution proof explains 8πΏ as 2πΏ +6πΏ, one contribution from each finite
sum. An alternative reads the same correction from the generating
function near π€ =1. Its pole of
order four accounts for the cubic growth of ππ; the next principal-part coefficient
accounts for the relative term of order 1/π.
An alternative coefficient calculation. Write Ξ (π€) =(ππ€;π)β, so that (π€;π)β =(1 βπ€)Ξ (π€), and for π‘ β₯1 put
ππ‘(π€)=(ππ‘π€2;π)β(π;π)π‘(ππ‘π€;π)2β,S(π€)=βπ‘β₯1π€π‘(ππ‘(π€)βπβ1).
On each
compact subset of |π€| <πβ1 the
bracket is ππ(ππ‘) uniformly, so
S is analytic there.
Separating the summand π‘ =0, whose
denominator carries the factor (1 βπ€)2, and subtracting the constant
part of the remaining summands gives
D(π€):=π4(1βπ€)4πΊπ(π€)=π4Ξ (π€)3((1+π€)(ππ€2;π)βΞ (π€)2+π€π+(1βπ€)S(π€)).
Hence D is analytic on a disc of
radius greater than 1.
At π€ =1 one has Ξ (1) =π and Ξ β²(1)/Ξ (1) = βπΏ, since Ξ β²/Ξ = ββπβ₯1ππ/(1 βπππ€).
Also (ππ€2;π)β equals π at π€ =1, and the logarithmic derivative of
(ππ€2;π)β/Ξ (π€)2 at π€ =1 is β2πΏ +2πΏ =0. Therefore the bracket has value
3/π and derivative (2 βπΏ)/π at π€ =1; here S(1) =πΏ/π, because ππ‘(1) =(π(1 βππ‘))β1 and
βπ‘β₯1((1 βππ‘)β1 β1) =πΏ.
Combining with the factor π4Ξ β3, whose logarithmic derivative
at 1 is 3πΏ, gives
D(1)=3,Dβ²(1)=3(3πΏ+2βπΏ3)=2+8πΏ.
The first
two principal-part terms are determined by D(1) and Dβ²(1). The terms from its
second and third derivatives have coefficients of order at most π; after the full principal part is
subtracted, the remainder is analytic on a larger disc and has
exponentially decreasing coefficients. Keeping only the first two terms
therefore gives
ππ=3(π+33)β(2+8πΏ)(π+22)+ππ(π),ππππ=1β8πΏπ+1+ππ((π+1)β2),
the
second line because 3(π+33) =ππ(π +3)/(π +1) and (π+22) =ππ/(π +1).
The logarithmic error is again summable. The product argument in the
first proof therefore gives 0 <A(π) <β and, using βπ<π(π +1)β1 =logβ‘π +πΎπΈ +π(1),
βπ<πππβΌA(π)πΆππβ8πΏ.
PropositionΒ 3.7
supplies the uniform shift bound and the fixed-shift limit, so
TheoremΒ 3.6
applies to (ππ) and gives πβπ βΌM(π)3ππ΅ππ2πβπ<πππ. Substituting the product
asymptotic proves both displays.Β β»
Comparison with the cubic bound.
An extension with limiting weight ratio