The row identity.
We choose row operations whose first surviving coefficients can be
computed explicitly. The auxiliary lemma applies to more general
products than the moments above; its use here is to identify the
coefficient that remains after each row operation.
Work over 𝐴 =ℤ[[𝑞]]. Let 𝐻(𝑋) =1 +∑𝑠≥1𝑎𝑠(𝑞)𝑋𝑠. For
nonnegative integers 𝑚,𝑡, put
𝑊𝑚(𝑡)=𝑞(𝑚+1)𝑡𝑚∏𝑟=1𝐻(𝑞𝑟),𝐷𝑗=𝑗−1∏𝑟=0(𝐼−𝑞𝑟N),𝐸(𝑚,𝑗)=𝑚𝑗−𝑗(𝑗−1)2,
where N is the backward shift (N𝑓)𝑚 =𝑓𝑚−1 in the index
𝑚. Thus 𝐷𝑗 is Zudilin’s backward-difference
operator in [9].
Write ¯𝐻 =𝐻mod𝑞 and ℎ𝑟 =[𝑋𝑟]¯𝐻(𝑋)−1, with ℎ𝑟 =0 for 𝑟 <0 and ℎ0 =1.
The lemma allows arbitrary coefficients 𝑎𝑠(𝑞) ∈ℤ[[𝑞]]; the only normalisation
imposed on 𝐻 is its constant term
1. It therefore covers the product
series used below, but not an unnormalised series with a different
constant term. All statements here are coefficientwise identities of
formal series, with no analytic convergence assumption. For a general
𝐻, the coefficient in the lemma can
vanish. The application below computes it for the chosen products and
proves the required nonvanishing.
Lemma 3.5 (a coefficient of each transformed row).
For 𝑚 ≥𝑗 ≥0 one has 𝐷𝑗𝑊𝑚(𝑡) ∈𝑞𝐸(𝑚,𝑗)𝐴 and
[𝑞𝐸(𝑚,𝑗)]𝐷𝑗𝑊𝑚(𝑡)=(−1)𝑗ℎ𝑗−𝑡.
Proof. Let 𝑒𝑢 denote a
formal basis vector indexed by 𝑢 ≥0, and suppress the argument 𝑞 in 𝑎𝑠(𝑞). We use the transition rule
Φ𝑒0=∑𝑠≥1𝑎𝑠𝑒𝑠−1,Φ𝑒𝑢=(𝑞𝑢−1)𝑒𝑢−1+𝑞𝑢∑𝑠≥1𝑎𝑠𝑒𝑢+𝑠−1(𝑢>0).
In Φ𝑗𝑒𝑡, a coefficient means
the sum of the weights of length-𝑗
paths from 𝑡 to the specified
endpoint. Each step lowers its index by at most one. A path ending at
𝑢 therefore visits only indices at
most max(𝑡,𝑢 +𝑗), so each
coefficient is a finite sum. No action on the algebraic direct sum is
assumed.
For 𝑚 ≥1, a direct expansion
gives
𝑊𝑚(𝑢)−𝑊𝑚−1(𝑢)=𝑞𝑚∑𝑣(Φ𝑒𝑢)𝑣𝑊𝑚−1(𝑣),
since both sides equal
𝑞𝑚𝑢𝑚−1∏𝑟=1𝐻(𝑞𝑟)(𝑞𝑢𝐻(𝑞𝑚)−1).
For the induction step take 𝑚 ≥𝑗 +1 and note that 𝐸(𝑚,𝑗) −𝐸(𝑚 −1,𝑗) =𝑗. Applying 𝐼 −𝑞𝑗N to 𝑞𝐸(𝑚,𝑗)𝑊𝑚−𝑗(𝑢) therefore gives
𝑞𝐸(𝑚,𝑗)(𝑊𝑚−𝑗(𝑢) −𝑊𝑚−𝑗−1(𝑢)). The
preceding identity contributes one further factor 𝑞𝑚−𝑗, and 𝐸(𝑚,𝑗) +𝑚 −𝑗 =𝐸(𝑚,𝑗 +1). Induction, starting
with 𝐷0 =𝐼, gives
𝐷𝑗𝑊𝑚(𝑡)=𝑞𝐸(𝑚,𝑗)∑𝑢(Φ𝑗𝑒𝑡)𝑢𝑊𝑚−𝑗(𝑢).(13)
The endpoint sum is also locally finite: ord𝑊𝑚−𝑗(𝑢) =(𝑚 −𝑗 +1)𝑢, so
only finitely many endpoints contribute to any fixed power of 𝑞.
Modulo 𝑞 the operator
simplifies: ¯Φ𝑒𝑢 = −𝑒𝑢−1
for 𝑢 >0, and ¯Φ𝑒0 =∑𝑠≥1𝑎𝑠(0)𝑒𝑠−1.
Put 𝑐𝑗,𝑡 =(¯Φ𝑗𝑒𝑡)0.
Then 𝑐𝑗+1,𝑡 = −𝑐𝑗,𝑡−1 for 𝑡 >0 and 𝑐𝑗+1,0 =∑𝑗+1𝑠=1𝑎𝑠(0)𝑐𝑗,𝑠−1,
the sum terminating because a state above 𝑗 cannot reach 0 in 𝑗 steps. Now 𝑐0,𝑡 =𝛿𝑡,0 =ℎ−𝑡, and if 𝑐𝑗,𝑡 =( −1)𝑗ℎ𝑗−𝑡 for all 𝑡 then 𝑐𝑗+1,𝑡 =( −1)𝑗+1ℎ𝑗+1−𝑡 for 𝑡 >0 at once, while for 𝑡 =0 the identity ¯𝐻¯𝐻−1 =1 gives ℎ𝑗+1 = −∑𝑠≥1𝑎𝑠(0)ℎ𝑗+1−𝑠 and
hence 𝑐𝑗+1,0 =( −1)𝑗∑𝑠≥1𝑎𝑠(0)ℎ𝑗+1−𝑠 =( −1)𝑗+1ℎ𝑗+1.
So 𝑐𝑗,𝑡 =( −1)𝑗ℎ𝑗−𝑡 for all
𝑗,𝑡. Since ord𝑊𝑚−𝑗(𝑢) =(𝑚 −𝑗 +1)𝑢 and
𝑚 ≥𝑗, in (13) only
the state 𝑢 =0 contributes at degree
𝐸(𝑚,𝑗), which gives both
assertions. ◻
Apply the lemma to a summand of the normalised remainder:
𝑇𝑚,𝑡=𝑞(𝑚+1)𝑡(𝑞;𝑞)3𝑚(𝑞𝑡+1;𝑞)𝑚(𝑞𝑚+𝑡+1;𝑞)𝑚+1.
Set
𝐻𝑡(𝑋)=(1−𝑋)3(1−𝑞𝑡𝑋)2(1−𝑞𝑡𝑋2)(1−𝑞𝑡+1𝑋2).
Then 𝑇𝑚,𝑡 =(1 −𝑞𝑡+1)−1𝑊𝐻𝑡𝑚(𝑡), by
the telescoping identity
𝑚∏𝑟=1𝐻𝑡(𝑞𝑟)=(𝑞;𝑞)3𝑚(1−𝑞𝑡+1)(𝑞𝑡+1;𝑞)𝑚(𝑞𝑚+𝑡+1;𝑞)𝑚+1.
The two denominator products contribute the consecutive factors 1 −𝑞𝑡+2,…,1 −𝑞𝑡+2𝑚+1. Reducing
modulo 𝑞 gives ¯𝐻𝑡 =(1 −𝑋)3 for 𝑡 >0, so ℎ(𝑡)𝑟 =(𝑟+22), and ¯𝐻0 =(1 −𝑋)4/(1 +𝑋), so ℎ(0)𝑟 =(𝑟 +1)(𝑟 +2)(2𝑟 +3)/6. The scalar
factor (1 −𝑞𝑡+1)−1 has
constant term 1. The rows are 𝑣∗𝑚 =∑𝑡≥0𝑇𝑚,𝑡. In 𝐷𝑗 only the shifts 𝑚,𝑚 −1,…,𝑚 −𝑗 occur. Since 𝑚 ≥𝑗 and ord𝑇𝑚−𝑟,𝑡 =(𝑚 −𝑟 +1)𝑡 ≥(𝑚 −𝑗 +1)𝑡,
only finitely many 𝑡 can affect any
fixed coefficient after any of these shifts. This justifies
interchanging the sum and 𝐷𝑗, and
then extracting its coefficient at degree 𝐸(𝑚,𝑗). Terms with 𝑡 >𝑗 contribute ℎ(𝑡)𝑗−𝑡 =0, so Lemma 3.5 at 𝑚 =𝑗 +ℓ, where 𝐸(𝑗 +ℓ,𝑗) =𝑗(𝑗 +1)/2 +𝑗ℓ, gives
𝐷𝑗𝑣∗𝑗+ℓ=(−1)𝑗(𝑗+1)2(𝑗+2)2𝑞𝑗(𝑗+1)/2+𝑗ℓ+𝑂(𝑞𝑗(𝑗+1)/2+𝑗ℓ+1)(𝑗,ℓ≥0).(14)
The coefficient is obtained by summing
ℎ(0)𝑗+𝑗∑𝑡=1ℎ(𝑡)𝑗−𝑡=(𝑗+1)(𝑗+2)(2𝑗+3)6+(𝑗+23)=(𝑗+1)2(𝑗+2)2.
Proof of Theorem 3.4.
The operators 𝐷𝑗 act by lower
unitriangular row operations, so they leave det(𝑣∗𝑖+𝑗)0≤𝑖,𝑗<𝑁
unchanged. By (14) the entry
in row 𝑗 and column ℓ has order 𝑗(𝑗 +1)/2 +𝑗ℓ. For 𝑁 =2, the matrix of entry orders is
(0012).
The off-diagonal product is the unique term of order one. Its
permutation sign cancels the negative leading coefficient of the second
row, giving 6𝑞 +𝑂(𝑞2). At arbitrary
rank the same argument selects the reversed permutation. For a
permutation 𝜍, the
corresponding Leibniz term has weight ∑𝑗(𝑗(𝑗 +1)/2 +𝑗𝜍(𝑗)).
By the rearrangement inequality, ∑𝑗𝑗𝜍(𝑗) is uniquely minimised
by the reversal 𝜍(𝑗) =𝑁 −1 −𝑗,
the values 𝑗 being distinct. The
minimum weight is ∑𝑗<𝑁𝑗2 =𝑁(𝑁 −1)(2𝑁 −1)/6, so
exactly one Leibniz term attains it and no cancellation is possible
there. The sign of the reversal is ( −1)𝑁(𝑁−1)/2, which cancels ∏𝑗<𝑁( −1)𝑗, and the surviving
coefficient is
𝑁−1∏𝑗=0(𝑗+1)2(𝑗+2)2=(𝑁!)2(𝑁+1)!2𝑁.
◻