Plectis

Long paper

Geometric Moments and Rational Lambert Values: Proofs and Further Results

Summary. The record proves the fixed-base Hankel asymptotic for Zudilin's normalised moments and the rational-base irrationality criterion log b/log a<0.4056830213840605…. The determinant proof separates geometric node configurations from positive weights; the irrationality proof instead specialises Zudilin's 2004 polynomial forms after cancellation. Coefficient moments, evaluated lattices and congruence conditions distinguish integer divisibility from nonzero remainder decay. A broader ratio-limit argument and a rational-value model remain unformalised. None of these arguments settles F(3/2).

This paper covers the sharp fixed-base Hankel asymptotic and rational-base irrationality proofs, coefficient moments and arithmetic obstructions, with ordinary ratio-limit and rational-value models distinguished from checked results.

It does not establish the validity of claims tagged Lean, which belongs to the cited kernel-checked source, or a solution to ErdΕ‘s #1049, which remains open.

In this paper

The size of VNβˆ— at a fixed base

The preceding bounds leave an error of order N in the logarithm. We now determine its linear, logarithmic and constant terms. There are two calculations: compare every node configuration with the first N 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 VNβˆ—. Throughout this subsection 0<q<1 is fixed, and we write

P=(q;q)∞,M(q)=∏dβ‰₯1(1βˆ’qd)βˆ’d,L=F(1/q)=βˆ‘rβ‰₯1qr1βˆ’qr,

with BN=βˆ‘j<Nj2 and CN=∏k<Nck=(N!)2(N+1)!/2N as above. Both products converge to positive finite values, and the abbreviation L is used only in this subsection; it is the value at q of the divisor series L of SectionΒ 10. The estimates in this subsection concern real determinants at fixed q. Denominator valuations for the 2004 forms were treated separately in SectionΒ 2.

For ak=1, the moments are (1βˆ’qm+1)βˆ’1, so Cauchy’s determinant formula gives the model calculation

(15)DN(0):=det(11βˆ’qi+j+1)0≀i,j<N=qBNΞ”N∏0≀i,j<N(1βˆ’qi+j+1),Ξ”N=∏d=1Nβˆ’1(1βˆ’qd)2(Nβˆ’d).

The denominator has limiting multiplicity d at 1βˆ’qd, while

Ξ”NP2N=∏d<N(1βˆ’qd)βˆ’2d∏dβ‰₯N(1βˆ’qd)βˆ’2N⟢M(q)2.

The convergence follows from βˆ‘d|log⁑(1βˆ’qd)|<∞ and the bound Oq(NqN) for the logarithm of the second product. Hence

(16)DN(0)∼M(q)3qBNP2N.

Two copies of M(q) 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 ak>0 satisfy, for fixed constants C and ΞΊ,

ak+hak⟢1(kβ†’βˆž)Β for each fixedΒ hβ‰₯0,ak+hak≀C(1+h)ΞΊ(k,hβ‰₯0).

Put Mm=βˆ‘kβ‰₯0akq(m+1)k and DN=det(Mi+j)0≀i,j<N. Then

DN∼M(q)3qBNP2N∏k=0Nβˆ’1ak(Nβ†’βˆž).

Proof of TheoremΒ 3.6. We may increase C to at least 1 and ΞΊ to at least 0. The bound ah≀Ca0(1+h)ΞΊ makes βˆ‘kβ‰₯0akqkΞ΄qk a finite positive measure with moments Mm. Apply Cauchy–Binet to finite truncations. Entrywise convergence gives convergence of their determinants, and monotone convergence applies to the nonnegative tuple sums. Hence

(17)DN=βˆ‘0≀k0<β‹―<kNβˆ’1∏i<Nakiqki∏i<j<N(qkiβˆ’qkj)2.

The tuple ki=i contributes qBNΞ”N∏i<Nai. It does not account for the whole asymptotic. For example, replace (0,1,…,Nβˆ’2,Nβˆ’1) by (0,1,…,Nβˆ’2,N). The ratio of its contribution to that of the original tuple is

(18)qaNaNβˆ’1∏d=1Nβˆ’1(1βˆ’qd+11βˆ’qd)2=qaNaNβˆ’1(1βˆ’qN1βˆ’q)2⟢q(1βˆ’q)2.

Only pairs involving the moved node change, so the product telescopes. At q=1/2 the limiting ratio is 2: the displaced tuple contributes twice as much asymptotically as (0,1,…,Nβˆ’1).

Subtracting (0,1,…,Nβˆ’1) gives the shifts Ξ»i=kiβˆ’i. Since ki+1β‰₯ki+1, these are nonnegative and weakly increasing; reversing their order gives a partition. For N=4, for example,

(0,1,3,5)⏟(ki)⟼(0,0,1,2)⏟(λi)⟼(2,1,0,0)⏟(μj).

In general, put ΞΌj=Ξ»Nβˆ’j for 1≀j≀N, so that ΞΌ1β‰₯β‹―β‰₯ΞΌNβ‰₯0. Division by the contribution of ki=i turns the summand for this tuple into

(19)WN(ΞΌ)=qβˆ‘j≀N(2jβˆ’1)ΞΌj∏j≀NaNβˆ’j+ΞΌjaNβˆ’jβ‹…βˆ1≀j<k≀N(1βˆ’qkβˆ’j+ΞΌjβˆ’ΞΌk1βˆ’qkβˆ’j)2.

We order the indices increasingly, so their geometric nodes decrease towards zero. The jth index from the right contributes once through qki and twice through each of its jβˆ’1 pairs with later indices in the squared Vandermonde. Its displacement therefore has exponent 1+2(jβˆ’1)=2jβˆ’1.

Set WN(ΞΌ)=0 when ΞΌ has more than N 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,…,Nβˆ’1): for instance, (2,1) affects the last two for every Nβ‰₯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 N telescopes. The factors with j<k≀ℓ are fixed, and those with j>β„“ equal 1. For Nβ‰₯β„“ and each j≀ℓ, the remaining factors satisfy

(20)∏k=β„“+1N1βˆ’qkβˆ’j+ΞΌj1βˆ’qkβˆ’j=(qN+1βˆ’j;q)ΞΌj(qβ„“+1βˆ’j;q)ΞΌj⟢1(qβ„“+1βˆ’j;q)ΞΌj.

Thus each WN(ΞΌ) has the same limit as for constant weights.

To pass this limit through the sum over partitions, we need a bound independent of N that remains summable as both the number and size of the parts grow. Only factors with j≀ℓ can differ from 1. For each such j, the denominator product in the Vandermonde quotient is at least P, and its numerator product is at most 1. Together with the polynomial shift bound, this gives

WN(ΞΌ)≀(CPβˆ’2)β„“βˆj=1β„“(1+ΞΌj)ΞΊq(2jβˆ’1)ΞΌj.

Let A=βˆ‘uβ‰₯1(1+u)ΞΊquβˆ’1<∞. Dropping the ordering of the positive parts, we obtain

βˆ‘β„“(ΞΌ)=β„“supNWN(ΞΌ)≀(CPβˆ’2)β„“βˆj=1β„“βˆ‘uβ‰₯1(1+u)ΞΊq(2jβˆ’1)u≀(CPβˆ’2A)β„“qβ„“2.

Here q(2jβˆ’1)u≀q2jβˆ’1quβˆ’1 and 1+3+β‹―+(2β„“βˆ’1)=β„“2. Since q<1, the bound (CPβˆ’2A)β„“qβ„“2 is summable over β„“, independently of N.

Dominated convergence gives the same limiting partition sum as for ak=1. By (15), this sum tends to M(q), the reciprocal of the limiting Cauchy denominator. Together with ΔN/P2N→M(q)2, it accounts for all three copies of M(q). Restoring the principal tuple gives

DNDN(0)∏k<Nak⟢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 ak=P4γk. A factorisation into two finite Gaussian multinomial sums supplies both the shift hypotheses and the first correction to ak/ck. The shift hypotheses allow application of the determinant theorem; the correction determines the product ∏k<Nak.

Proposition 3.7 (the moment weights as a product of two finite sums). For rβ‰₯1 and kβ‰₯0 set

Rk(r)(q)=βˆ‘n1+β‹―+nr=k(q;q)k∏j≀r(q;q)nj,

the sum being over compositions into r nonnegative parts, so that Rk(r) is the sum of all q-multinomial coefficients of degree k and length r, equivalently the multivariate Rogers–SzegΕ‘ polynomial at unit argumentsΒ [42]. These are polynomial identities in the indeterminate q; at a prime power, r=2 gives the Galois number counting the subspaces of a k-dimensional space over the field of that order, and r=3 counts the flags V1βŠ†V2 in that space. The source is cited for those identifications; the factorisation below is proved here. Then

Ξ³k(q)=Rk(2)(q)Rk(3)(q)(q;q)k.

The product Rk(2)Rk(3) lies in Zβ‰₯0[q], has degree ⌊k2/4βŒ‹+⌊k2/3βŒ‹, and has coefficient sum 6k. Consequently (q;q)kΞ³k∈Z[q], and the weights ak=P4Ξ³k satisfy

P5ck≀ak≀Pβˆ’1ck(kβ‰₯0),

together with ak+h/ak≀Pβˆ’6(1+h)3, so the hypotheses of TheoremΒ 3.6 hold with ΞΊ=3.

Proof. Write E(z)=(z;q)βˆžβˆ’1. Euler’s expansion gives E(z)r=βˆ‘nβ‰₯0Rn(r)zn/(q;q)n. We seek the product Rk(2)Rk(3) at a common index k. Splitting k=n+j keeps Rk(3) fixed while the sum over ((q;q)n(q;q)j)βˆ’1 supplies Rk(2)/(q;q)k. The operator βˆ‚qf(z)=(f(z)βˆ’f(qz))/z produces the needed index shift: βˆ‚qzj=(1βˆ’qj)zjβˆ’1 cancels the last factor of (q;q)j. Thus

βˆ‚qE=E,βˆ‚qnE(z)3=βˆ‘jβ‰₯0Rn+j(3)(q;q)jzj.

The q-Leibniz rule, obtained by induction from the product rule, is

βˆ‚qn(fg)(z)=βˆ‘j=0n[nj]q(βˆ‚qjf)(qnβˆ’jz)(βˆ‚qnβˆ’jg)(z).

Write hk(z)=βˆ‘j=0k[kj]q(z;q)j. Applying the rule to Eβ‹…E and using E(qmz)=(z;q)mE(z) gives βˆ‚qkE(z)2=E(z)2hk(z); at k=1 both sides are E(z)2(2βˆ’z). Applying it once more, to Eβ‹…E2, gives

βˆ‚qnE(z)3=E(z)3βˆ‘k=0n[nk]q(z;q)khk(z),

Hence

βˆ‘nβ‰₯0wn(q;q)nβˆ‚qnE(z)3=E(w)E(z)3βˆ‘kβ‰₯0wk(z;q)k(q;q)khk(z).

The left-hand side is βˆ‘n,jβ‰₯0Rn+j(3)wnzj/((q;q)n(q;q)j). On setting z=w, all pairs with the same sum contribute to a single coefficient. Thus the coefficient of wk is Rk(3)βˆ‘j≀k((q;q)j(q;q)kβˆ’j)βˆ’1=Rk(2)Rk(3)/(q;q)k. On the right-hand side, put z=w, expand hk, and write k=j+t. The identity (w;q)j+t=(w;q)j(wqj;q)t gives

E(w)4βˆ‘jβ‰₯0wj(w;q)j2(q;q)jβˆ‘tβ‰₯0(wqj;q)t(q;q)twt.

The inner sum is evaluated by the q-binomial theorem:

βˆ‘tβ‰₯0(wqj;q)t(q;q)twt=(w2qj;q)∞(w;q)∞.

Finally, (w;q)j=(w;q)∞/(wqj;q)∞ cancels two of the five factors E(w), leaving

E(w)3βˆ‘jβ‰₯0wj(q;q)j(w2qj;q)∞(wqj;q)∞2=Gq(w).

Every step holds coefficientwise, and absolutely for |w| small.

A q-multinomial coefficient has nonnegative integral coefficients and degree βˆ‘i<jninj. At q=1 it equals the corresponding ordinary multinomial coefficient. Summing over compositions gives coefficient sum rk for Rk(r) and hence 6k for the product. The degree is maximised by balanced parts, which gives ⌊k2/4βŒ‹ for r=2 and ⌊k2/3βŒ‹ for r=3; no leading cancellation occurs because all coefficients are nonnegative.

For the bounds, write bk(r)=[zk](z;q)βˆžβˆ’r, so that Rk(r)=(q;q)kbk(r) and Ξ³k=(q;q)kbk(2)bk(3). Euler’s expansion exhibits bk(r) as an r-fold convolution of the sequence 1/(q;q)n, whose terms lie in [1,Pβˆ’1]; the number of compositions of k into r nonnegative parts is (k+rβˆ’1rβˆ’1), so

(k+rβˆ’1rβˆ’1)≀bk(r)≀Pβˆ’r(k+rβˆ’1rβˆ’1).

Combining these at r=2,3 with P≀(q;q)k≀1 and ck=(k+1)(k+22) gives P5ck≀ak≀Pβˆ’1ck. The shift bound follows from ck+h/ck≀(1+h)3.

It remains to show that ak+h/akβ†’1 for each fixed h. The two convolutions grow at least linearly and quadratically; it suffices to bound their increments by a constant and a linear function, respectively. Put en=1/(q;q)n. Since (q;q)n+1=(1βˆ’qn+1)(q;q)n, these terms satisfy 1=e0≀e1≀⋯≀Pβˆ’1, and the convolutions read bk(2)=βˆ‘i≀keiekβˆ’i and bk(3)=βˆ‘i≀keibkβˆ’i(2). Hence

bk+1(2)βˆ’bk(2)=βˆ‘i≀kei(ek+1βˆ’iβˆ’ekβˆ’i)+ek+1,bk+1(3)βˆ’bk(3)=βˆ‘i≀kei(bk+1βˆ’i(2)βˆ’bkβˆ’i(2))+ek+1.

Every bracket is nonnegative. In the first line the brackets add up to ek+1βˆ’1<Pβˆ’1 and each ei≀Pβˆ’1, so 0≀bk+1(2)βˆ’bk(2)≀Pβˆ’2+Pβˆ’1≀2Pβˆ’2. Inserting this bound into the k+1 terms of the second line gives 0≀bk+1(3)βˆ’bk(3)≀2(k+1)Pβˆ’3+Pβˆ’1≀3(k+1)Pβˆ’3. Against the lower bounds bk(2)β‰₯k+1 and bk(3)β‰₯(k+1)(k+2)/2 displayed above, the ratio bk+1(r)/bk(r) lies in [1,1+2Pβˆ’2/(k+1)] for r=2 and in [1,1+6Pβˆ’3/(k+2)] for r=3. Both ratios tend to 1, and therefore so does bk+h(r)/bk(r), a product of h consecutive ratios. Finally (q;q)k+h/(q;q)k=∏j=1h(1βˆ’qk+j)β†’1, and ak+h/ak=Ξ³k+h/Ξ³k 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 q, not coefficientwise positivity in q: for example, a0=P4 has q-coefficient βˆ’4. The cubic factors ck multiply to CN, but ak/ckβ†’1 alone would not determine the product of N weights. The first relative correction has order 1/k and accumulates as a harmonic sum; the remaining Oq(kβˆ’2) terms determine a finite constant. The two convolutions supply this first correction explicitly.

Theorem 3.8 (the size of VNβˆ— at a fixed base). Let Ξ³E be Euler’s constant. The product

A(q)=eβˆ’8Ξ³EL∏kβ‰₯0(akcke8L/(k+1))

converges to a positive finite value. With K(q)=A(q)M(q)3,

VNβˆ—(q)∼K(q)CNqBNP2NNβˆ’8L,

equivalently

log⁑VNβˆ—(q)=BNlog⁑q+log⁑CN+2Nlog⁑Pβˆ’8Llog⁑N+log⁑K(q)+o(1).

Proof of TheoremΒ 3.8. Put L=F(1/q), ek=(q;q)kβˆ’1 and E=Pβˆ’1. Replacing every ek by E in the two convolution sums gives E2(k+1) and E3(k+22). The factors k+1 and (k+22) count the nonnegative pairs and triples with sum k, respectively. For a fixed index i, there is one pair (i,kβˆ’i) and there are kβˆ’i+1 triples whose first coordinate is i. Thus a summable deficit in the sequence ei 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 dk=Eβˆ’ek, the identity ek+1βˆ’ek=qk+1ek+1 gives

dk=βˆ‘j>kqjej,0≀dk≀Eqk+11βˆ’q.

Thus βˆ‘dk and βˆ‘kdk converge. To evaluate the former, use Euler’s identity E(z)=βˆ‘kβ‰₯0ekzk=(z;q)βˆžβˆ’1 for |z|<1. Tonelli’s theorem and logarithmic differentiation at z=q give

S:=βˆ‘kβ‰₯0dk=βˆ‘jβ‰₯1jqjej=qEβ€²(q)=Eβˆ‘rβ‰₯1qr1βˆ’qr=EL.

In the double sum, qjej occurs once for each k=0,…,jβˆ’1, which explains the factor j. Logarithmic differentiation is justified by locally uniform convergence on |z|<1.

Let bk(r)=[zk]E(z)r. The finite-sum formula in PropositionΒ 3.7 becomes

ak=P4(q;q)kbk(2)bk(3).

We expand each finite convolution using ek=Eβˆ’dk. For two factors,

bk(2)=E2(k+1)βˆ’2Eβˆ‘i=0kdi+βˆ‘i=0kdidkβˆ’i=E2(k+1βˆ’2L)+Oq((k+1)qk).

For three factors, the terms with one di contribute βˆ’3E2βˆ‘i=0k(kβˆ’i+1)di. Here the sum equals (k+1)S+Oq(1) because βˆ‘idi<∞. Terms with two di are bounded by 3ES2, and the term with three is Oq((k+1)2qk). It follows that

bk(3)=E3((k+22)βˆ’3L(k+1))+Oq(1).

Dividing the correction βˆ’3L(k+1) by the leading term (k+22) gives the relative correction βˆ’6L/(k+2). The two-fold convolution contributes βˆ’2L/(k+1) in the same way. Since (q;q)k=P(1+Oq(qk)), the leading constants in ak=P4(q;q)kbk(2)bk(3) multiply to P5E5=1. Replacing (k+2)βˆ’1 by (k+1)βˆ’1 changes only the Oq((k+1)βˆ’2) error. Thus the coefficients 2L and 6L add, and with ck=(k+1)2(k+2)/2 we obtain

akck=(1βˆ’2Lk+1+Oq(qk))(1βˆ’6Lk+2+Oq((k+1)βˆ’2))(1+Oq(qk))=1βˆ’8Lk+1+Oq((k+1)βˆ’2).

In particular ak+1/ak→1, and multiplication of finitely many adjacent ratios gives ak+h/ak→1 for each fixed h. The uniform shift bound follows from Proposition 3.7, so Theorem 3.6 applies.

For large k, taking logarithms preserves the first-order term and leaves the summable error log⁑(ak/ck)+8L/(k+1)=Oq((k+1)βˆ’2). The finitely many earlier factors are positive, so the regularised product also has a positive finite limit: 0<A(q)<∞. Since ∏k<Nck=CN, the harmonic contribution separates as

log⁑∏k<NakCN=βˆ’8Lβˆ‘k<N1k+1+βˆ‘k<N(log⁑akck+8Lk+1).

The second sum tends to log⁑A(q)+8Ξ³EL. Using βˆ‘k<N(k+1)βˆ’1=log⁑N+Ξ³E+o(1), the Euler-constant terms cancel, giving ∏k<Nak∼A(q)CNNβˆ’8L. TheoremΒ 3.6 now gives the stated determinant asymptotic. Taking logarithms yields

log⁑VNβˆ—=BNlog⁑q+log⁑CN+2Nlog⁑Pβˆ’8Llog⁑N+log⁑K(q)+o(1).

β—»

The convolution proof explains 8L as 2L+6L, one contribution from each finite sum. An alternative reads the same correction from the generating function near w=1. Its pole of order four accounts for the cubic growth of ak; the next principal-part coefficient accounts for the relative term of order 1/k.

An alternative coefficient calculation. Write Ξ (w)=(qw;q)∞, so that (w;q)∞=(1βˆ’w)Ξ (w), and for tβ‰₯1 put

ft(w)=(qtw2;q)∞(q;q)t(qtw;q)∞2,S(w)=βˆ‘tβ‰₯1wt(ft(w)βˆ’Pβˆ’1).

On each compact subset of |w|<qβˆ’1 the bracket is Oq(qt) uniformly, so S is analytic there. Separating the summand t=0, whose denominator carries the factor (1βˆ’w)2, and subtracting the constant part of the remaining summands gives

D(w):=P4(1βˆ’w)4Gq(w)=P4Ξ (w)3((1+w)(qw2;q)∞Π(w)2+wP+(1βˆ’w)S(w)).

Hence D is analytic on a disc of radius greater than 1.

At w=1 one has Ξ (1)=P and Ξ β€²(1)/Ξ (1)=βˆ’L, since Ξ β€²/Ξ =βˆ’βˆ‘jβ‰₯1qj/(1βˆ’qjw). Also (qw2;q)∞ equals P at w=1, and the logarithmic derivative of (qw2;q)∞/Ξ (w)2 at w=1 is βˆ’2L+2L=0. Therefore the bracket has value 3/P and derivative (2βˆ’L)/P at w=1; here S(1)=L/P, because ft(1)=(P(1βˆ’qt))βˆ’1 and βˆ‘tβ‰₯1((1βˆ’qt)βˆ’1βˆ’1)=L. Combining with the factor P4Ξ βˆ’3, whose logarithmic derivative at 1 is 3L, gives

D(1)=3,Dβ€²(1)=3(3L+2βˆ’L3)=2+8L.

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 k; 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

ak=3(k+33)βˆ’(2+8L)(k+22)+Oq(k),akck=1βˆ’8Lk+1+Oq((k+1)βˆ’2),

the second line because 3(k+33)=ck(k+3)/(k+1) and (k+22)=ck/(k+1).

The logarithmic error is again summable. The product argument in the first proof therefore gives 0<A(q)<∞ and, using βˆ‘k<N(k+1)βˆ’1=log⁑N+Ξ³E+o(1),

∏k<Nak∼A(q)CNNβˆ’8L.

PropositionΒ 3.7 supplies the uniform shift bound and the fixed-shift limit, so TheoremΒ 3.6 applies to (ak) and gives VNβˆ—βˆΌM(q)3qBNP2N∏k<Nak. Substituting the product asymptotic proves both displays.Β β—»

About this paper

Authorship and AI use. Will Cook built and directed the research infrastructure and maintains the public release. He reviewed claims when he could. AI agents did most of the research and drafting. Cook did not independently verify every claim.

Cite and contact. Cite this paper by its title, author and date above, with its PDF; cite the earlier sources it uses for a mathematical result. For software, use the release citation and give the commit used. Contact Will with questions or corrections.

Prefer the manuscript?

Equations typeset from the exact TeX. Rendered from the LaTeX manuscript source with file fingerprint sha256:1f04580702ed8f8c. Its file fingerprint matched the published source record, so the PDF above, the LaTeX, and this page are one manuscript. Redeploying the site regenerates this page from whatever the public repository holds at that moment.