LXXXIX.2 (1999)
Hankel determinants for the Fibonacci word and Pad´ e approximation
by
Teturo Kamae (Osaka), Jun-ichi Tamura (Tokyo) and Zhi-Ying Wen (Beijing)
1. Introduction. The aim of this paper is to give a concrete and inter- esting example of the Pad´e approximation theory as well as to develop the general theory so as to find a quantitative relation between the Hankel de- terminant and the Pad´e pair. Our example is the formal power series related to the Fibonacci word.
The Fibonacci word ε(a, b) on an alphabet {a, b} is the infinite sequence ε(a, b) = b ε
0ε b
1. . . b ε
n. . .
(1)
:= abaababaabaab . . . (b ε
n∈ {a, b}), which is the fixed point of the substitution
(2) σ : a → ab, b → a.
The Hankel determinants for an infinite word (or sequence) ϕ = ϕ
0ϕ
1. . . (ϕ
n∈ K) over a field K are
(3) H
n,m(ϕ) := det(ϕ
n+i+j)
0≤i,j≤m−1(n = 0, 1, . . . ; m = 1, 2, . . .).
It is known [2] that the Hankel determinants play an important role in the theory of Pad´e approximation for the formal Laurent series
(4) ϕ(z) =
X
∞ k=0ϕ
kz
−k+h.
Let K((z
−1)) be the set of formal Laurent series ϕ as above of z with coeffi- cients in K and h ∈ Z providing a nonarchimedean norm kϕk := exp(−k
0+h) with k
0= inf{k : ϕ
k6= 0}. Let ϕ be as above with h = −1. We say that a
1991 Mathematics Subject Classification: 41A21, 11B39.
This research was partially supported by the Hatori Project of the Mathematical Society of Japan.
[123]
pair (P, Q) ∈ K[z]
2of polynomials of z over K is a Pad´e pair of order m for ϕ if
(5) kQϕ − P k < exp(−m), Q 6= 0, deg Q ≤ m.
A Pad´e pair (P, Q) of order m for ϕ always exists and the rational func- tion P/Q ∈ K(z) is uniquely determined for each m = 0, 1, . . . The element P/Q ∈ K(z) with P, Q satisfying (5) is called the mth diagonal Pad´e ap- proximation for ϕ. A number m is called a normal index if (5) implies deg Q = m. Note that P/Q is irreducible if m is a normal index, although it can be reducible for a general m. A normal Pad´e pair (P, Q), i.e., deg Q is a normal index, is said to be normalized if the leading coefficient of Q is equal to 1. It is a classical result that m is a normal index for ϕ if and only if the Hankel determinant det(ϕ
i+j)
0≤i,j≤m−1is nonzero. Note that 0 is always a normal index and the determinant for the empty matrix is considered to be 1, so that the above statement remains valid for m = 0.
We succeed in obtaining a quantitative relation between the Hankel de- terminant and the normalized Pad´e pair. Namely,
(6) det(ϕ
i+j)
0≤i,j≤m−1= (−1)
[m/2]Y
z; Q(z)=0
P (z)
for any normal index m with the normalized Pad´e pair (P, Q), where Q
z; Q(z)=0
indicates a product taken over all zeros z of Q with their multi- plicity (Theorem 6).
We are specially interested in the Pad´e approximation theory applied to the Fibonacci words ε := ε(1, 0) and ε := ε(0, 1), where 0, 1 are considered as elements in the field Q, since we have the following remark.
Remark 1. Let M be a matrix of size m×m with entries consisting of two independent variables a and b. Then det M = (a − b)
m−1(pa + (−1)
m−1qb), where p and q are integers defined by
p = det M |
a=1, b=0, q = det M |
a=0, b=1.
P r o o f. Subtracting the first column vector from all the other column vectors of M , we see that det M is divisible by (a − b)
m−1as a polynomial in Z[a, b]. Hence, det M = (a − b)
m−1(xa + yb) for integers x, y. Setting (a, b) = (1, 0), (0, 1), we get the assertion.
In Section 2, we study the structure of the Fibonacci word, in particular, its repetition property. The notion of singular words introduced in Z.-X. Wen and Z.-Y. Wen [5] plays an important role.
In Section 3, we give the value of the Hankel determinants H
n,m(ε) and
H
n,m(ε) for the Fibonacci words in some closed forms. It is a rare case where
the Hankel determinants are determined completely. Another such case is
for the Thue–Morse sequence ϕ consisting of 0 and 1, where the Hankel
determinants H
m,n(ϕ) modulo 2 are obtained, and the function H
m,n(ϕ) of (m, n) is proved to be 2-dimensionally automatic (see [1]).
In Section 4, we consider the self-similar property of the values H
n,m(ε) and H
n,m(ε) for the Fibonacci words. The quarter plane {(n, m) : n ≥ 0, m ≥ 1} is tiled by 3 kinds of tiles with the values H
n,m(ε) and H
n,m(ε) on it with various scales.
In Section 5, we develop a general theory of Pad´e approximation. We also obtain the admissible continued fraction expansion of ϕ
εand ϕ
ε, the formal Laurent series (4) with h = −1 for the sequences ε and ε, and determine all the convergents p
k/q
kof the continued fractions. It is known in general that the set of the convergents p
k/q
kfor ϕ is the set of diagonal Pad´e approximations and the set of degrees of q
k’s in z coincides with the set of normal indices for ϕ.
2. Structure of the Fibonacci word. In what follows, σ denotes the substitution defined by (2), and
b
ε = b ε
0ε b
1. . . b ε
n. . . (b ε
n∈ {a, b})
is the (infinite) Fibonacci word (1). A finite word over {a, b} is sometimes considered to be an element of the free group generated by a and b with inverses a
−1and b
−1. For n = 0, 1, . . . , we define the nth Fibonacci word F
nand the nth singular word W
nas follows:
(7) F
n:= σ
n(a) = σ
n+1(b), W
n:= β
nF
nα
−1n, where we put
(8) α
n= β
m=
a (n even, m odd), b (n odd, m even),
and we define W
−2to be the empty word and W
−1:= a for convenience.
Let (f
n; n ∈ Z) be the Fibonacci sequence:
(9) f
n+2= f
n+1+ f
n(n ∈ Z), f
−1= f
0= 1.
Then |F
n| = |W
n| = f
n(n ≥ 0), where |ξ| denotes the length of a finite word ξ.
For a finite word ξ = ξ
0ξ
1. . . ξ
n−1and a finite or infinite word η = η
0η
1. . . over an alphabet, we denote
(10) ξ ≺
kη
if ξ = η
kη
k+1. . . η
k+n−1. We simply write
(11) ξ ≺ η
and say that ξ is a subword of η if ξ ≺
kη for some k. For a finite word
ξ = ξ
0ξ
1. . . ξ
n−1and i with 0 ≤ i < n, we denote the ith cyclic permutation
of ξ by C
i(ξ) := ξ
iξ
i+1. . . ξ
n−1ξ
0ξ
1. . . ξ
i−1. We also define C
i(ξ) := C
i0(ξ) with i
0:= i − n[i/n] for any i ∈ Z.
In this section, we study the structure of the Fibonacci word b ε and discuss the repetition property. The following two lemmas were obtained by Z.-X.
Wen and Z.-Y. Wen [5] and we omit the proofs.
Lemma 1. We have the following statements:
(1) b ε = F
nF
n−1F
nF
n+1F
n+2. . . (n ≥ 1),
(2) F
n= F
n−1F
n−2= F
n−2F
n−1β
n−1α
−1nβ
nα
n(n ≥ 2), (3) F
nF
n≺ b ε (n ≥ 3),
(4) b ε = W
−1W
0W
1W
2W
3. . . , (5) W
n= W
n−2W
n−3W
n−2(n ≥ 1),
(6) W
nis a palindrome, that is, W
nstays invariant under reading the letters from the end (n ≥ −2),
(7) C
i(F
n) ≺ b ε (n ≥ 0, 0 ≤ i < f
n),
(8) C
i(F
n) 6= C
j(F
n) for any i 6= j, moreover , they are different already before their last places (n ≥ 1, 0 ≤ i < f
n),
(9) W
n6= C
i(F
n) (n ≥ 0, 0 ≤ i < f
n),
(10) ξ ≺ b ε and |ξ| = f
nimply that either ξ = C
i(F
n) for some i with 0 ≤ i < f
nor ξ = W
n(n ≥ 0).
Lemma 2. For any k ≥ −1, we have the decomposition of b ε as follows:
b
ε = (W
−1W
0. . . W
k−1)W
kγ
0W
kγ
1. . . W
kγ
n. . . ,
where all the occurrences of W
kin b ε are picked up and γ
nis either W
k+1or W
k−1corresponding to b ε
nis a or b, respectively. That is, any two different occurrences of W
kdo not overlap and are separated by W
k+1or W
k−1.
We introduce another method to discuss the repetition property of b ε.
Let N be the set of nonnegative integers. For n ∈ N, let
(12) n =
X
∞ i=0τ
i(n)f
i,
τ
i(n) ∈ {0, 1} and τ
i(n)τ
i+1(n) = 0 (i ∈ N)
be the regular expression of n in the Fibonacci base due to Zeckendorf. For m, n ∈ N and a positive integer k, we define
(13) m ≡
kn
if τ
i(m) = τ
i(n) for all i < k.
Lemma 3. We have b ε
n= a if and only if τ
0(n) = 0.
P r o o f. We use induction on n. The lemma holds for n = 0, 1, 2. Assume
that it holds for any n ∈ N with n < f
kfor some k ≥ 2. Take any n ∈ N
with f
k≤ n < f
k+1. Then, since 0 ≤ n − f
k< f
k−1, we have
n =
k−1
X
i=0
τ
i(n − f
k)f
i+ f
k,
which gives the regular expression if τ
k−1(n − f
k) = 0. If τ
k−1(n − f
k) = 1, then we have the regular expression n = P
k−2i=0
τ
i(n − f
k)f
i+ f
k+1. In any case, we have τ
0(n) = τ
0(n − f
k). On the other hand, since b ε starts with F
kF
k−1by Lemma 1, we have b ε
n= b ε
n−fk. Hence, b ε
n= a if and only if τ
0(n) = 0 by the induction hypothesis. Thus, we have the assertion for any n < f
k+1, and by induction, we complete the proof.
Lemma 4. Let n = P
∞i=0
n
if
iwith n
i∈ {0, 1} (i ∈ N). Assume that n
in
i+1= 0 for 0 ≤ i < k. Then n
i= τ
i(n) for 0 ≤ i < k.
P r o o f. If there exists i ∈ N such that n
in
i+1= 1, let i
0be the maximum such i. Take the maximum j such that n
i0+1= n
i0+3= n
i0+5= . . . = n
j= 1. Then, replacing f
i0+ f
i0+1+ f
i0+3+ f
i0+5+ . . . + f
jby f
j+1, we have a new expression of n:
n = X
∞ i=0n
0if
i:=
i
X
0−1 i=0n
if
i+ f
j+1+ X
∞ i=j+3n
if
i.
This new expression is unchanged at the indices less than k, and is either regular or has a smaller maximum index i with n
0in
0i+1= 1. By continuing this procedure, we finally get the regular expression of n, which does not differ from the original expression at the indices less than k. Thus, n
i= τ
i(n) for any 0 ≤ i < k.
Lemma 5. For any n ∈ N and k ≥ 0, τ
0(n + f
k) 6= τ
0(n) if and only if either n ≡
k+2f
k+1− 2 or n ≡
k+2f
k+1− 1. Moreover ,
b
ε
n+fk− b ε
n=
(−1)
k−1(a − b) (n ≡
k+2f
k+1− 2), (−1)
k(a − b) (n ≡
k+2f
k+1− 1), where a and b are considered as independent variables.
P r o o f. If k = 0, we can verify the statement by a direct calculation.
Assume that k ≥ 1 and τ
k(n) = 0. Then n + f
k=
k−1
X
i=0
τ
i(n)f
i+ f
k+ X
∞ i=k+1τ
i(n)f
i.
By Lemma 4, we have τ
0(n + f
k) = τ
0(n) if k ≥ 2 or if k = 1 and τ
0(n) = 0.
In the case where k = 1, τ
0(n) = 1 and τ
2(n) = 0, since n + f
k= 1 + 2 +
X
∞ i=3τ
i(n)f
i= f
2+ X
∞ i=3τ
i(n)f
i,
we have τ
0(n + f
k) = 0 by Lemma 4. On the other hand, in the case where
k = 1, τ
0(n) = 1 and τ
2(n) = 1, since
n + f
k= 1 + 2 + 3 + X
∞ i=4τ
i(n)f
i= f
0+ f
3+ X
∞ i=4τ
i(n)f
i, we have τ
0(n + f
k) = 1 by Lemma 4.
Thus, in the case where k ≥ 1 and τ
k(n) = 0, τ
0(n + f
k) 6= τ
0(n) if and only if k = 1, τ
0(n) = 1 and τ
2(n) = 0, or equivalently, if and only if n ≡
k+2f
k+1− 2 with k = 1. Note that n ≡
k+1f
k+1− 1 with k = 1 contradicts τ
k(n) = 0.
Now assume that k ≥ 1 and τ
k(n) = 1. Take the minimum j ≥ 0 such that τ
k(n) = τ
k−2(n) = τ
k−4(n) = . . . = τ
j(n) = 1. Then since 2f
i= f
i+1+ f
i−2for any i ∈ N, we have
n + f
k=
j−3
X
i=0
τ
i(n)f
i+ f
j−2(14)
+ f
j+1+ f
j+3+ f
j+5+ . . . + f
k+1+ X
∞ i=k+2τ
i(n)f
i,
where the first term on the right-hand side vanishes if j = 0, 1, 2. Hence by Lemma 4, τ
0(n + f
k) = τ
0(n) if j ≥ 4.
In the case where j = 3, τ
0(n + f
k) = τ
0(n) holds if τ
0(n) = 0 by (14) and Lemma 4. If τ
0(n) = 1, then by (14) and Lemma 4, τ
0(n + f
k) = 0.
Thus, for j = 3, τ
0(n + f
k) 6= τ
0(n) if and only if τ
0(n) = 1.
If j = 2, then by the assumption on j, we have τ
0(n) = 0. On the other hand, since f
0= 1, by (14) and Lemma 4, we have τ
0(n + f
k) = 1. Thus, τ
0(n + f
k) 6= τ
0(n).
If j = 1, then τ
0(n) = 0 since τ
1(n) = 1 by the assumption on j. On the other hand, since f
−1= 1, we have τ
0(n + f
k) = 1 by (14) and Lemma 4.
Thus, τ
0(n + f
k) 6= τ
0(n).
If j = 0, then by the assumption on j, τ
0(n) = 1. On the other hand, since f
−2= 0, we have τ
0(n + f
k) = 0 by (14) and Lemma 4. Thus, τ
0(n + f
k) 6=
τ
0(n).
By combining all the results as above, we get the first part.
The second part follows from Lemma 3 and the fact that for any k ≥ 0, f
k+1− 1 = f
k+ f
k−2+ . . . + f
iwith i = 0 if k is even and i = 1 if k is odd. Hence, τ
0(f
k+1− 1) = τ
0(f
h+1− 2) =
a (k odd, h even), b (k even, h odd).
Lemma 6. For any k ≥ 0, W
k≺
nε if and only if n ≡ b
k+2f
k+1− 1.
P r o o f. By Lemma 2, the smallest n ∈ N such that W
k≺
nb ε is
f
−1+ f
0+ f
1+ . . . + f
k−1= f
k+1− 1,
which is the smallest n ∈ N such that n ≡
k+2f
k+1− 1. Let n
0:= f
k+1− 1.
Then the regular expression of n
0is
n
0= f
k+ f
k−2+ f
k−4+ . . . + f
d,
where d = (1 − (−1)
k)/2. The next n with n ≡
k+2n
0is clearly n = f
k+2+ f
k+ f
k−2+ . . . + f
d,
which is, by Lemma 2, the next n such that W
k≺
nε since f b
k+f
k+1= f
k+2. For i = 1, 2, . . . , let
n
i= n
0+ X
∞ j=0τ
j(i)f
k+2+j.
Then it is easy to see that n
iis the ith n after n
0such that n ≡
k+2f
k+1− 1.
We prove by induction on i that n
iis the ith n after n
0such that W
k≺
nε. b Assume that it is so for i. Then by Lemma 4, W
kγ
iW
k≺
niε. Hence, the b next n after n
isuch that W
k≺
nε is n b
i+ f
k+ |γ
i|. Thus, we have
n
i+ f
k+ |γ
i| = n
i+ f
k+ f
k+11
εbi=a+ f
k−11
εbi=b= n
i+ f
k+21
τ0(i)=0+ f
k+11
τ0(i)=1= n
i+1, which completes the proof.
Lemma 7. Let k ≥ 0 and n, i ∈ N satisfy n ≡
k+1i.
(1) If 0 ≤ i < f
k, then τ
0(n + j) = τ
0(i + j) for any j = 0, 1, . . . , f
k+2− i − 3.
(2) If f
k≤ i < f
k+1, then τ
0(n + j) = τ
0(i + j) for any j = 0, 1, . . . , f
k+3− i − 3.
P r o o f. (1) We prove the lemma by induction on k. The assertion holds for k = 0. Let k ≥ 1 and assume that the assertion is valid for k − 1. For j = 0, 1, . . . , f
k− i, we have n + j ≡
ki + j and hence, τ
0(n + j) = τ
0(i + j).
Let j
0= f
k− i. Then, since n + j
0≡
ki + j
0≡
k0, we have τ
0(n + j
0+ j) = τ
0(i + j
0+ j) = τ
0(j) for any j = 0, 1, . . . , f
k+1− 3 by the induction hypothesis. Thus, τ
0(n + j) = τ
0(i + j) for any j = 0, 1, . . . , f
k+2− i − 3.
This proves (1).
(2) In this case, τ
k+1(n) = 0. Hence, n ≡
k+2i. Therefore, we can apply (1) with k + 1 for k. Thus, we get (2).
Let n, m, i ∈ N with m ≥ 2 and 0 < i < m. We call n an (m, i)-shift invariant place in b ε if
b
ε
nb ε
n+1. . . b ε
n+m−1= b ε
n+iε b
n+i+1. . . b ε
n+i+m−1.
We call n an m-repetitive place in b ε if there exist i, j ∈ N with i > 0 and
i + j < m such that n + j is an (m, i)-shift invariant place in b ε. Let R
mbe
the set of m-repetitive places in b ε.
Lemma 8. (1) Let n ≡
k+10 for some k ≥ 1. Then n is an (f
k+1− 2, f
k)- shift invariant place in b ε.
(2) Let n ≡
k+1f
kfor some k ≥ 2. Then n is an (f
k+1− 2, f
k−1)-shift invariant place in b ε.
P r o o f. (1) Since the least i ≥ n such that either i ≡
k+2f
k+1− 1 or i ≡
k+2f
k+1− 2 is not less than n + f
k+1− 2, by Lemma 5, we have
b
ε
nb ε
n+1. . . b ε
n+fk+1−3= b ε
n+fkε b
n+fk+1. . . b ε
n+fk+fk+1−3.
(2) Since the minimum i ≥ n such that either i ≡
k+1f
k− 1 or i ≡
k+1f
k− 2 is n + f
k+1− 2, by Lemma 5, we have b
ε
nb ε
n+1. . . b ε
n+fk+1−3= b ε
n+fk−1ε b
n+fk−1+1. . . b ε
n+fk−1+fk+1−3.
Theorem 1. The pair (n, m) of nonnegative integers satisfies n ∈ R
mif one of the following two conditions holds:
(1) f
k+ 1 ≤ m ≤ f
k+1− 2, n − i ≡
k+10 and i ≤ n for some k ≥ 1 and i ∈ Z with f
k+ 1 ≤ m + i ≤ f
k+1− 2.
(2) f
k−1+ 1 ≤ m ≤ f
k+1− 2, i ≤ n and n − i ≡
k+1f
kfor some k ≥ 2 and i ∈ Z with f
k−1+ 1 ≤ m + i ≤ f
k+1− 2.
Remark 2. The “if and only if” statement actually holds in Theorem 1 in place of “if” since we will prove later that H
n,m6= 0 if none of the conditions (1) and (2) hold.
P r o o f (of Theorem 1). Assume (1) and i ≥ 0. By Lemma 8(1), n − i is an (f
k+1− 2, f
k)-shift invariant place. Then n is an (m, f
k)-shift invariant place since i + m ≤ f
k+1− 2. Thus, n ∈ R
mas f
k< m.
Assume (1) and i < 0. Then, since n−i is an (f
k+1−2, f
k)-shift invariant place and m ≤ f
k+2−2, it is an (m, f
k)-shift invariant place. Moreover, since f
k− i < m, n is an m-repetitive place.
Assume (2) and i ≥ 0. Then, n − i is an (f
k+1− 2, f
k−1)-shift invariant place by Lemma 8(2). Then, n is an (m, f
k−1)-shift invariant place since i + m ≤ f
k+1− 2. Thus, n is an m-repetitive place as f
k−1< m.
Assume (2) and i < 0. Then, since n − i is an (f
k+1− 2, f
k−1)-shift invariant place and m ≤ f
k+1− 2, it is an (m, f
k−1)-shift invariant place.
Then n is an m-repetitive place, since f
k−1− i < m. Thus, n ∈ R
m. Corollary 1. The place 0 is m-repetitive for an m ≥ 2 if m 6∈
S
∞k=1
{f
k− 1, f
k}.
Remark 3. The “if and only if” statement actually holds in Corollary 1 in place of “if” since we prove later that H
0,m6= 0 if m ∈ S
∞k=1