arxiv: v1 [math.co] 23 May 2017

Transcription

1 RANK-METRIC CODES AND ZETA FUNCTIONS I. BLANCO-CHACÓN, E. BYRNE, I. DUURSMA, AND J. SHEEKEY arxiv: v [math.co] 23 May 207 Abstract. We define the rank-metric zeta function of a code as a generating function of its normalized q-binomial moments. We show that, as in the Hamming case, the zeta function gives a generating function for the weight enumerators of rank-metric codes. We further prove a functional equation and derive an upper bound for the minimum distance in terms of the reciprocal roots of the zeta function. Finally, we show invariance under suitable puncturing and shortening operators and study the distribution of zeroes of the zeta function for a family of codes.. Introduction In [4, 5], the zeta polynomial of a q-ary linear code of length n and minimum distance d is defined to be the unique polynomial P C (T ) of degree at most n d such that the generating function has expansion around T given by (.) P C (T ) ( T )( qt ) (xt + y( T ))n = + A C(x, y) x n T n d +, q where A C (x, y) is the Hamming weight enumerator of C. For an MDS code C, the weight enumerator A C (x, y) is obtained by setting P C (T ) =. In general the degree of P C (T ) is given by n + 2 d d 0 and is equal to zero if and only if the code is MDS. For codes with high minimum distance and high dual minimum distance, the weight enumerator is completely determined by the polynomial P C (T ) of small degree. In such cases there is an advantage in encoding the weight enumerator A C (x, y) by the smaller polynomial P C (T ). Among the properties of the zeta polynomial P C (T ) is an easy transform between the polynomials P C (T ) and P C (T ) for the code and its dual, and invariance under puncturing and shortening for codes that have a transitive automorphism group. So that for example a cyclic code and its extended code use (.) with the same zeta polynomial P C (T ) but with different n and d. The special form of the generating function is motivated by its interpretation for Reed-Solomon codes. A Reed-Solomon codeword represents the values of a polynomial in the elements of a finite field. The generating function counts the number of zeroes of polynomials and therefore describes the weight distribution of Reed-Solomon codes.

2 2 I. BLANCO-CHACÓN, E. BYRNE, I. DUURSMA, AND J. SHEEKEY In this paper we establish that the zeta polynomial and its properties have q analogues for rank-metric codes. Rank-metric codes have a Singleton type upper bound for the parameters. Codes that attain the bound are called maximum rank distance codes (MRD codes). They include Gabidulin codes that are the q-analogues of Reed-Solomon codes. The rank distance distribution of a MRD code is uniquely determined by the parameters of the code. The rank-metric zeta function has the property that P C (T ) = for MRD codes, and thus in particular for Gabidulin codes. In general the zeta polynomial is of low degree if both the minimum distance and the dual minimum distance of a code are close to the Singleton bound. In Section 2 we summarize some preliminary material and definitions. In Section 3 we introduce the zeta function Z C (T ) of a rank metric code C in terms of its normalized q-binomial moments, and we define the q analogue of the zeta polynomial for a rank-metric code with the generating function Z C (T ). In analogy with the Hamming metric case P C (T ) = if and only if C is an MRD code. In Section 4 we relate the weight enumerator of a rank-metric code C with its zeta polynomial and obtain the rank-metric analogue of (.). Moreover, we show that coefficients of the zeta polynomial coincide with those in the representation of a weight enumerator as a Q-linear combination of MRD weight enumerators. In Section 5 we establish relations between the polynomials P C (T ) and P C (T ) for a code and its dual. We then demonstrate that the zeroes of P C (T ) can be related to the minimum distance of a code by giving an upper bound for the minimum distance in terms of these zeroes. In Section 6 we define operations of puncturing and shortening of rank-metric codes and show that the action of these on the average weight enumerator of a code can be realized in terms of q-derivatives. As in the Hamming metric case the zeta polynomial is invariant under puncturing and shortening. We introduce the normalized weight enumerator and express the action of puncturing and shortening on it in terms of q-commuting operators. In Section 7 we discuss the zeroes of the zeta polynomial. For a self-dual code, the reciprocal zeroes occur in pairs {α, q m /α} and occur as conjugate pairs if and only if both have absolute value q m/2. We consider some classes of codes with zeta polynomials for which all complex zeroes have the same absolute value. 2. Preliminaries We will assume that m, n are positive integers with n m and that q is a prime power. We write F m n q to denote the m n matrices with entries in F q. For any X F m n q, we write ker X or X to denote the (right) nullspace of X in F n q. That is, X := {y F n q : Xy T = 0}.

3 RANK-METRIC CODES AND ZETA FUNCTIONS 3 Unless explicitly stated otherwise, we assume that C is an F q -linear subspace of F m n q. Definition 2.. The dual code of C is the F q -linear code C := {Y F m n q : Tr(XY t ) = 0 for all X C}. The map (X, Y ) Tr(XY t ) defines an inner product on the space F m n q, so we have dim(c ) = mn dim(c) and C = C. Definition 2.2. The rank distance between matrices X, Y F m n q is d(x, Y ) := rk(x Y ). For C 2, the minimum rank distance of C is the integer defined by d(c) := min{d(x, Y ) : X, Y C, X Y }. The weight distribution of C is the integer vector W (C) = (W t (C) : 0 t n), where, for all t {0,..., n}, W t (C) := {X C : rk(x) = t}. The rank metric weight enumerator of C is the bivariate polynomial n (2.) W C (x, y) = W t (C)x n t y t. t=0 Recall that the Gaussian binomial or q-binomial coefficient is defined by [ (q n := r] n )(q n q) (q n q r ) (q r )(q r q) (q r q r if r {,..., n}, ) if r = 0, q 0 otherwise. This quantity counts the number of r-dimensional subspaces of an n-dimensional subspace [ over [ F q. Since q is fixed throughout this paper, for brevity we write n n to mean. The rank-metric analogue of the Singleton bound says that r] r] q C q m(n d+) for any C of minimum rank-distance d [2]. Rank metric codes that meet this bound are called maximum-rank-distance (MRD) codes. It has been known for some decades that such codes exist for all choices of m, n, d [2, 7, 3]. In [2] it was shown that the weight enumerator of an MRD code in F m n q of minimum distance d is uniquely determined and given by n t M n,d (x, y) := x n + ( ) t i q (t i 2 ) n t (q m(i d+) )x n t y t, t i t=d where the coefficient of x n t y t counts the number of matrices in the MRD code of rank t. It is not hard to see that for fixed n the M n,d are linearly independent over Q. Any code C whose dual code has minimum distance at least 2 has weight

4 4 I. BLANCO-CHACÓN, E. BYRNE, I. DUURSMA, AND J. SHEEKEY enumerator that can be expressed as a Q-linear combination of the MRD weight enumerators. That is, there exist p 0,..., p n d Q satisfying W C (x, y) = p 0 M n,d (x, y) + + p n d M n,n (x, y). These coefficients will turn out to define the zeta polynomial of C, that we introduce in the next section. Let P be a partially ordered set (poset, from now on). The Möbius function for P is defined via the recursive formula (2.2) µ(x, x) =, µ(x, z) = x y<z µ(x, y). Lemma 2.3 (Möbius Inversion Formula). Let f, g : P Z be any two functions on P. Then. f(x) = g(y) if and only if g(x) = µ(x, y)f(y). x y x y 2. f(x) = g(y) if and only if g(x) = µ(y, x)f(y). x y x y In particular, for the subspace lattice of F m n q (with partial order defined by set inclusion), and for two subspaces U and V of dimensions u and v, we have that (2.3) µ (U, V ) = ( ) v u q (v u 2 ) if U V 0 otherwise. 3. The Zeta Function We introduce the zeta function of a rank metric code C in terms of its normalized q-binomial moments, following the approach for the Hamming metric case as given in [6]. In order to do so we use the notion of a shortened code (cf. [2]). Definition 3.. Let U F q F n q be a subspace of dimension u. The shortened subcode of C with respect to U is { C U := X C : U X }. The strict shortening of C by U is { Ĉ U := X C : U = X }.

5 RANK-METRIC CODES AND ZETA FUNCTIONS 5 Notice that C U is a subspace (indeed a subcode of C), but ĈU in general is not. Clearly every element of ĈU has rank exactly n u. Definition 3.2. For u 0, the u-th binomial moment of C is defined by B u (C) := ( C U ). dim(u)=n u Lemma 3.4 is a straightforward consequence of the following duality result. The result is used in [2] to obtain a short proof of the MacWilliams Identity for rank-metric codes. Lemma 3.3 ([2], Lemma 28). Let U be a subspace of F n q of dimension n u. Then C U = C C U q m(n u). Proof. Let R U be the subspace of m n matrices over F q with row vectors in U, and let b : C R U be the bilinear form with b(x, Y ) = Tr(XY t ). The left null space is C U and the right null space is R U C = C U. So that C / C U = R U / C U. Lemma 3.4. Let C have dimension k and minimum rank distance d, and let C have minimum distance d. Then 0 if u < d B u (C) = (q k m(n u) n ) if u > n d u Proof. Any X C satisfies rk(x ) n d, so if U has dimension n u > n d then C U = {0}. Similarly, C U = {0} if dim(u ) = u > n d since every element X C has rank at least d. Then from Lemma 3.3 we get that { C U = if u < d q k m(n u) if u > n d. Since n u counts number of subspaces of F n q of dimension n u, the result follows. We remark that in the instance that C is an MRD code, the values B u (C) of Lemma 3.4 were given in [3]. Moreover, we may invoke Lemma 3.3 along with

6 6 I. BLANCO-CHACÓN, E. BYRNE, I. DUURSMA, AND J. SHEEKEY the Singleton bound to deduce that if U has dimension u then C U C /q mu q m(n u d+). Definition 3.5. For u 0, the u-th normalized binomial moment of C is defined by b u (C) := B u+d(c) ]. [ n u+d Then b u (C) counts the average number of non-zero elements of the shortened code C U where U has dimension n u d. We extend the definition of b u to all u Z, by setting { 0 if u < 0. b u (C) = q k m(n u d), if u > n d d. Note that if C is an MRD code of minimum distance d then k = m(n d + ) and d = n + 2 d. Then b u (C) = q m(u+) for all u 0. In particular, the values b u (C) are independent of the minimum distance of C. For this reason we write b u instead of b u (C) in the cases that C is MRD. Definition 3.6. The zeta function of C (c.f. [6]) is defined by Z C (T ) := (q m ) u 0 b u (C)T u. It is straightforward to check that for u / {0,..., n d d + 2} we have the following relation among the b u, namely, (3.) b u (C) (q m + )b u (C) + q m b u 2 (C) = 0. Definition 3.7. For u 0, define p u (C) := (q m ) (b u (C) (q m + )b u (C) + q m b u 2 (C)). The zeta polynomial of C is defined by P C (T ) := n d d+2 u=0 p u (C)T u. Clearly, for any C, the zeta polynomial P C (T ) has degree at most n d +, and p u (C) = 0 for any u > n d d + 2. We will assume that the minimum distance of the dual code is at least 2, and write n d P C = p u (C)T u. u=0

7 RANK-METRIC CODES AND ZETA FUNCTIONS 7 The recursion relates the zeta function and zeta polynomial via Z C (T ) = (q m + )T Z C (T ) q m T 2 Z C (T ) + P C (T ), and hence we obtain that the generating function Z C (T ) satisfies the equation Z C (T ) = P C (T ) ( T )( q m T ). This immediately yields the following for an MRD code. Lemma 3.8. Let C be an MRD code. Then P C (T ) = and Z C (T ) = ( T )( q m T ). Proof. C is MRD if and only if d + d = n + 2, in which case (3.) holds for all u > 0 and P C (T ) =. 4. Weight Enumerators and Zeta Functions We will establish a relation between the weight enumerator of a code and its zeta function, giving the rank metric analogue of [5, Theorem 9.5]. We first approach this using Möbius inversion. Given a subspace U of F n q, we have the relation: (4.) C U = U V ĈV, consequently, by the Möbius inversion formula we get (4.2) ĈU = U V µ(u, V ) C V, for µ(u, V ) defined as in (2.3). Next, we derive a more explicit description of W C (x, y) in terms of the normalized binomial moments. Lemma 4.. Let C have minimum distance d and let be U a subspace of F n q dimension u. Then ĈU = n d u w=0 ( ) w q (w 2) dim(v ) = u + w U V ( C V ) of

8 8 I. BLANCO-CHACÓN, E. BYRNE, I. DUURSMA, AND J. SHEEKEY Proof. For each subspace V, denote its dimension by v. By the Möbius inversion formula, we have ĈU = U V = = ( ) v u q (v u 2 ) CV, U V,v n d U V,v n d n d = using Lemma 3.4. ( ) v u q (v u 2 ) CV + ( ) v u q (v u 2 ) CV + ( ) v u q (v u 2 ) C V + v=u U V U V,v n d+ U V,v n d+ n v=n d+ ( ) v u q (v u 2 ) CV ( ) v u q (v u 2 ), n u ( ) v u q (v u 2 ). v u Adding and subtracting to the inner sum of the left hand side, this becomes n d (4.3) ( ) v u q (v u 2 ) n n u ( C V ) + ( ) v u q (v u 2 ). v u U V v=u Setting v = u + w, (4.3) equals (4.4) n d u w=0 ( ) w q (w 2) U V ( C V ) + v=u Since the right sum vanishes the result follows. n u w=0 [ n u Proposition 4.2. Let C have minimum distance d. Then t W t (C) = ( ) t i q (t i 2 ) n t b i d (C). t i Proof. Let t = n u, and write W n t (C) = and setting u + w = v, this equals n d u w=0 For a fixed V, there are [ u+w u (4.5) W t (C) = n d u w=0 dim(u)=u w ( ) w q (w 2) U V ( C V ). ] ( ) w q (w 2). ĈU. Applying Lemma 4. ] U s contained in V of dimension u. Hence u + w ( ) w q (w 2) B u n (u+w) (C),

9 RANK-METRIC CODES AND ZETA FUNCTIONS 9 which equals t d ( ) w n n t + w t d n t q (w 2) b t w d (C) = ( ) w q (w 2) b t w d (C). n t + w n t t w w=0 w=0 Setting t w = i and observing that [ t w = t t w], the result follows. Equivalently, the following statement holds. Corollary 4.3. Let C have minimum distance d. Then ( n n [ ) W C (x, y) = x n + ( ) t i q (t i 2 ) n t ]x n t y t b i d (C). t i t=i In particular, in the case of an MRD code with weight enumerator M n,d (x, y) we get (4.6) M n,d (x, y) x n = ( n n ( ) t i q (t i 2 ) t=i Definition 4.4. For each r {0,..., n} define [ n t ][ t i ]x n t y t ) (q m(i d+) ). φ n,n r (x, y) := (q m ) (M n,r (x, y) (q m + )M n,r+ (x, y) + q m M n,r+2 (x, y)), and n φ n (T ) := φ n,r (x, y)t r. r=0 (4.7) Therefore, φ n (T ) ( T )( q m T ) (qm ) n (M n,r (x, y) x n ) T n r (mod T n+ ) r=0 Lemma 4.5. For each r {0,..., d}, the coefficient of T n r in the expression, is given by Z C (T )φ n (T ) = P C (T )φ n (T ) ( T )( q m T ) ( n r ) (q m ) p i (C)M n,r+i (x, y) x n.

10 0 I. BLANCO-CHACÓN, E. BYRNE, I. DUURSMA, AND J. SHEEKEY Proof. From (4.7), the coefficient of T n r in the left-hand-side of the above expression is ( n r ) n r (q m ) p i (C)M n,r+i (x, y) x n p i (C). Now from Definition 3.6 we have that n r p i (C) = b n r(c) q m b n r (C) q m. Then since r d, we have n r > n d d, and so b n r (C) = q k m(r d) and b n r (C) = q k m(r+ d). Therefore n r p i (C) = (qk m(r d) ) q m (q k m(r+ d) ) q m proving the claim. An explicit expression for φ n (T ) is given by: =, Lemma 4.6. φ n (T ) = ( n n ( ) t i q (t i 2 ) t=i [ n t ][ t i ]x t y n t ) T n i. Proof. By definition, we have, (q m )φ n,n r (x, y) = M n,r (x, y) (q m + )M n,r+ (x, y) + q m M n,r+2 (x, y). Then by (4.6), we get n = ( ) t r q (t r 2 ) n t x t y n t (q m ) t r t=r n + ( ) t r q (t r 2 ) n t x t y n t (q 2m (q m + )(q m )) t r + t=r+ n n + ( ) t i q (t i 2 ) n t x t y n t (b i r+ (q m + )b i r + q m b i r ) t i i=r+2 t=i n = ( ) t r q (t r 2 ) n t x t y n t (q m ), t r t=r as the third sum vanishes due to the relation (3.).

11 RANK-METRIC CODES AND ZETA FUNCTIONS Therefore, the weight enumerators M n,d (x, y) and W C (x, y) can be expressed as n (4.8) M n,d (x, y) = x n + φ n,n i (x, y)b i d, (4.9) W C (x, y) = x n + n φ n,n i (x, y)b i d (C). In fact, the polynomials φ n,n r (x, y) are related to a well known class of q- polynomials (see, e.g. [9]), and are given by [ n φ n,n r (x, y) = p r] n r (x, y)y r, where p k (x, y) := k (x qj y) = [ k ( )k j q (k j 2 ) k x j] j y k j. Therefore, (4.0) W C (x, y) = x n + n n b i d (C) p i n i (x, y)y i. We now establish a connection between the weight enumerator of a code and its zeta function. Theorem 4.7. The coefficient of T n d in the expression, is given by Z C (T )φ n (T ) = P C (T )φ n (T ) ( T )( q m T ) W C (x, y) x n q m. In particular, the zeta polynomial P C (T ) is the unique polynomial of degree at n d most n d such that W C (x, y) = p i M n,d+i (x, y). Proof. It is immediate from (4.0) that the coefficient of T n d in Z C (T )φ n (T ) is (q m ) (W C (x, y) x n ). From Lemma 4.5, we have that this coefficient is equal to (q m ) ( x n + n d n d p im n,d+i (x, y)), and thus W C (x, y) = p i (C)M n,d+i (x, y). Since the M n,d+i (x, y) are linearly independent over Q, the result follows.

12 2 I. BLANCO-CHACÓN, E. BYRNE, I. DUURSMA, AND J. SHEEKEY Example 4.8. Let us consider m = n = 3. Then we have the following weight enumerators for MRD codes M 3,3 = x 3 + (q 3 )y 3 M 3,2 = x 3 + (q 3 )(q 2 + q + )xy 2 + (q 3 )(q 3 q 2 q)y 3 M 3, = x 3 + (q 3 )(q 2 + q + )x 2 y + (q 3 ) 2 (q 2 + q)xy 2 + (q 3 )(q 3 q)(q 3 q 2 )y 3 Taking q = 2 we get M 3,3 = x 3 + 7y 3 M 3,2 = x xy 2 + 4y 3 M 3, = x x 2 y + 294xy y 3 We formally define M 3,4 = x n, the weight enumerator of the trivial code. Consider a code C of constant rank 2; that is, every nonzero codeword in C has rank-weight two. If C has dimension k, then W C (x, y) = x 3 + (q k )xy 2. Such codes always exist for k = 3, and exist with k = 4, q = 2. For k = 3, we have a code C withe weight enumerator and hence W C (x, y) = ( M3,2 q 2 (q 3 q 2 q)m 3,3 + q 3 ) M 3,4, + q + P C (T ) = (q3 q 2 q)t + q 3 T 2 q 2. + q + For k = 4 and q = 2, we get a code C 2 with zeta polynomial P C2 (T ) = 49 ( 5 30T + 64T 2 ). This reflects the fact that x 3 + 5xy 2 = ( 5(x xy 2 + 4y 3 ) 30(x 3 + 7y 3 ) + 64x 4). 49 From considering an MRD code of 3 3 matrices with minimum distance 2 over F q 2 as a subspace of 6 6 matrices over F q, we produce a code C 3 with weight enumerator W C3 (x, y) = x 6 + (q 6 )(q 4 + q 2 + )x 2 y 4 + (q 6 )(q 6 q 4 q 2 )y 6,

13 RANK-METRIC CODES AND ZETA FUNCTIONS 3 and zeta polynomial P C3 (T ) = + ( q6 + q 4 + q 3 + q 2 + q)t + q 6 T 2 q 4 + q 3 + q 2. + q + Taking q = 2 we get 34T + 64T 2 P C3 (T ) = Duality and an Upper Bound A number of authors have described the duality theory of rank metric codes [2, 8, 2]. In [8], the authors give a direct analogue of MacWilliams duality theorem relating the weight enumerator of a code with that of its dual, or more precisely, with the q-transform of its dual. Define an operation τ on homogeneous polynomials in two variables by (5.) f τ (x, y) := f(x + (q m )y, x y), where f is the q-transform of f, which is a Q-linear map. The following version of MacWilliams duality theorem for rank metric codes can be read in [8, Theorem ]. Theorem 5.. W C (x, y) = C W τ C (x, y). If W C = W C, C is called formally self-dual. Delsarte [2] showed that the dual code of an MRD code is also MRD; the dual of an MRD code with minimum distance d is an MRD code with minimum distance n d+2. Since an MRD code exists for all choices of m, n, d, the family of MRD weight enumerators is closed under τ and hence the MacWilliams identity. More precisely, an immediate corollary of Theorem 5. and the previous observation yields that: Corollary 5.2. M τ n,d (x, y) = qm(n d+) M n,n d+2 (x, y). Corollary 5.3. Let C F m n q be a matrix code with minimal rank distance d and dimension k. Let d be the minimal rank distance of the dual code C. Setting r t W C (x, y) = p i M n,d+i and W C (x, y) = s j M n,d +j, we have that

14 4 I. BLANCO-CHACÓN, E. BYRNE, I. DUURSMA, AND J. SHEEKEY a) r = t = n d d + 2; b) s j = p r j q m(d +j ) k ; Proof. Applying MacWilliams identity to W C (x, y) and using Corollary 5.2 gives r W C (x, y) = q m(n d i+) k p i M n,n d i+2. Setting d + j = n d i + 2, its minimal value (j = 0) happens when i = n d d + 2 and its maximal value (i = 0) when j = n d d + 2. Hence a) follows. Claim b) follows from the fact that the MRD family is a Q-basis, hence the decomposition of W C as a linear combination of MRD polynomials of different minimal distances must be unique. The MacWilliams identity can be translated into a functional relation between the zeta function of a code and that of its dual (or just between the zeta function of the code at different arguments, if the code is formally self dual), as we now show. Theorem 5.4. ( Z C (T ) = q m(n d+)+k T r 2 Z C q m T In particular, if C is formally self-dual, then Z C (T ) = q m(n d+)+k T r 2 Z C ( q m T ). ). r Proof. Let P C (T ) = p i T i and P C (T ) = Corollary 4.3, for 0 j r: r s i T i, for some p i, s i Q. By s j = p r j q m(d +j ) k, hence multiplying by T j r, summing in j and then changing j to r j: r ( ) T r P C (T ) = p j q m(d +r j ) k T j = q m(d +r ) k P C q m. T Thus ( ) P C (T ) = T r q m(n d+) k P C q m. T Taking the quotient with ( T )( q m T ), the result follows.

15 RANK-METRIC CODES AND ZETA FUNCTIONS 5 Example 5.5. We consider the zeta polynomials of the duals of the codes from Example 4.8. The code C 2 has m = n = 3, d = 2, k = 4, q = 2, and zeta polynomial P C2 (T ) = ( 5 30T + 64T 2 ). 49 Thus ( ) P C (T ) = T P C T = ( 4 5T + 60T 2 ) 49 The code C 3 has m = n = 6, d = 4, k = 2, and zeta polynomial Thus 34T + 64T 2 P C3 (T ) =. 3 ( ) P C (T ) = T P C T 34T + 64T 2 = = P C3 (T ). 3 Note that although C 3 and C3 have the same zeta polynomial (and hence the same zeta function), they do not have the same weight enumerator, as they have different minimum distances. In fact, if a code has quadratic zeta polynomial, which occurs if and only if d + d = n, then its dual will have the same zeta polynomial if and only if k = m(n d). We close this section by showing how the zeta function can be used to derive an upper bound on the minimum distance of code, which is the rank-metric analogue of [4, Section 2]. Theorem 5.6. Let P C (T ) = p 0 ( + at +...) be the zeta polynomial of a rank metric code C F m n q, of degree r n d and let a be the negative of the sum of its reciprocal roots. Then d log q [(a + q m + )(q ) + ]. Proof. First, observe that according to Theorem 4.7 the coefficient of T n d in Z C (T )φ n (T ) is (5.2) W C (x, y) x n q m.

16 6 I. BLANCO-CHACÓN, E. BYRNE, I. DUURSMA, AND J. SHEEKEY Second, since ( T )( q m T ) is the common zeta-function of all the MRD weight enumerators (in particular, of those with length n), again due to Theorem 4.7, we have n d φ n (T ) ( T )( q m T ) = q m (M n,d+i (x, y) x n ) T n d i +... Likewise, the coefficient of T n d in Z C (T )φ n (T ) is (5.3) q m r p i (M n,d+i (x, y) x n ). Within the coefficient of T n d in these two coincident expressions, we need to compare the coefficients of x n d y d and x n d y d+. The coefficient of x n d y d in (5.2) is W n d(c) q m, W n d(c) being the number of codewords of rank d in C, which equals the coefficient of x n d y d in (5.3), i.e., in p 0 q m (M n,d(x, y) x n ). Using Theorem 4.7, this last expression can be expanded as p 0 q m ( n ) φ n,n i (x, y)b i d = = p 0 q m n [ n n i ]b i d y i ( ) n i j q (n i j i 2 ) n i x j y n i j, j the coefficient of x n d y d p in which expression is precisely 0 n q [ d] m b0. Comparing with the corresponding coefficient in (5.2) we obtain p 0 0. Observe that the numbers b u are the normalized q-binomial coefficients of an MRD code of length n, and they are independent of the minimal distance d. Now, we compare the coefficients of x n d y d+. First, this coefficient is W n d (C) q m 0 in (4.3). In (5.3), this equals the coefficient of x n d y d+ in the sub-expression q m (p 0(M n,d (x, y) x n ) + p (M n,d+ (x, y) x n )). The coefficient in p q m (M n,d+(x, y) x n ) is, as before p n (5.4) q m b 0. d +

17 RANK-METRIC CODES AND ZETA FUNCTIONS 7 As for the coefficient in p 0 q m (M n,d(x, y) x n ), we expand this expression again as p 0 n q m φ n,n i (x, y)b i d, from which we isolate the sought for coefficient, which is ( [ ) p 0 n n n d (5.5) q m b ]b 0. d + d n d Adding (5.4) and (5.5), dividing by p 0 and taking into account that [ n d+] = q n d [ n q d+ d], yields q n d (5.6) ab 0 q d+ + b q n d n d q d+ b 0 0. n d Since n d n d = q n d q, it holds that, ab 0 + b b 0 q d+ q 0. Since b u = q m(u+) for each u 0, dividing by b 0 yields By taking logarithm the result follows. a + q m + qd+ q. 6. Puncturing and Shortening Puncturing and shortening are fundamental operations of coding theory. In [4], the zeta function was related to the normalized weight enumerator of a code, by successive puncturing and shortening of an MDS weight enumerator. We consider the rank-metric case in what follows. We have already considered shortened subcodes in order to define the zeta function of a code. We will now consider puncturing and shortening as operations that yield codes in F m (n ) q, that is, as projections of codes in F m n q. Such operations have been considered in [], with a slightly different definition. We will establish the invariance of the normalized weight enumerator of a code under puncturing and shortening. Definition 6.. Let H be a hyperplane in F n q. Fix a basis of H and let P H F n (n ) q be the matrix whose columns are the elements this basis of H, in some

18 8 I. BLANCO-CHACÓN, E. BYRNE, I. DUURSMA, AND J. SHEEKEY order. Let h F n q \H. We define the punctured and shortened codes of C with respect to H, respectively by: Π H (C) := {XP H : X C} F m (n ) q (punctured code), Σ h,h (C) := {XP H : X C, Xh T = 0} F m (n ) q (shortened code). Clearly Π H (C) and Σ h,h (C) are only well-defined up to a choice of basis of H. We assert that this is sufficient for our purposes, which is to define operations of puncturing and shortening on normalized [ weight ] [ enumerators. With these definitions, C can by punctured in any of = ways and can be shortened n n n ] n in any of q n = {( h, H) : dim H = n, h / H}. Lemma 6.2. Let H be a hyperplane in F n q and let P H Fq n (n ) that form a basis of H. Let X be non-zero in F m n q. Then { rk(x) if X H rk(xp H ) = rk(x) if X H have columns Proof. Let h F n q \H. Then rk(x) = rk(x[p H, h T ]) = rk([xp H, Xh T ]) = rk(xp H ) if and only if Xh T is contained in the column-space of XP H, which holds if and only if h X + H. X H if and only if F n q = X + H, in which case h X + H and rk(x) = rk(xp H ). If X H then h / X + H = H, so rk(xp H ) = rk(x). Corollary 6.3. Let H be a hyperplane in F n q and let h F n q \H. Let C have parameters [m n, k, d 2] over F q. Then the punctured code Π H (C) has parameters [m (n ), k, d ]. The shortened code Σ h,h (C) has parameters [m (n ), k m, d]. Proof. The lower bounds on the minimum distances of the codes follows directly from Lemma 6.2. Let P H F n (n ) q have columns that form a basis of H. Then XP H = 0 if and only if H X, in which case X = 0, since otherwise dim X = n rk(x) n d n 2. It follows that C and Π H (C) have the same cardinality. Let h F n q \H. Consider the map θ h : F m n q F m q : X Xh T. Then C ker θ h = C h := {X C : Xh T = 0}. Therefore dim C h = dim(c ker θ h )

19 RANK-METRIC CODES AND ZETA FUNCTIONS 9 k m. The result now follows since Σ h,h (C) is obtained as the punctured code of C h with respect to H, which both have the same dimension. Given an arbitrary code, C and hyperplanes H, H 2, the punctured codes Π H (C) and Π H2 (C) may have different weight enumerators. Neither is the weight enumerator of a shortened code Σ h,h (C) necessarily uniquely determined. However, in the case of an MRD code C we observe: Corollary 6.4. Let d 2. The shortened and punctured codes of an MRD code with weight enumerator M n,d (x, y) are also MRD and have weight enumerators M n,d (x, y) and M n,d (x, y), respectively. Proof. Let H be a hyperplane in F n q. Applying Corollary 6.4 together with the rank-metric Singleton bound, we see that the punctured code Π H (C) has parameters [m (n ), m(n d+), d ] and is thus MRD. Similarly the shortened code Σ h,h (C) has parameters [m (n ), m(n d), d] and is also an MRD code. We wish to define algebraic operations on the weight enumerator of a code corresponding to the expected weight enumerator of its punctured or shortened codes. In the Hamming metric case, such operations are defined in terms of partial derivatives [4]. In the rank-metric case, we may implement q-analogues of partial derivatives (c.f. [9], []). We define such operations as follows. Definition 6.5. Let f(x, y) = n f ix n i y i Q[x, y]. The partial q-derivatives of f(x, y) with respect to x and y are defined by: D x (f(x, y)) := D y (f(x, y)) := D q,x (f(x, y)) := f(qx, y) f(x, y) (q )x f(x, qy) f(x, y) (q )y f(qx, qy) f(x, qy) (q )x n n i = f i x n i y i, n i = f i x n i y i, n = f i q i n i x n i y i. It is straightforward to check that these are linear operators. In analogy with [4, Section 3] we now define the operations of puncturing and shortening on the homogeneous polynomials of degree n in Q[x, y] by n n P := (D q,x + D y ) and S := D x.

20 20 I. BLANCO-CHACÓN, E. BYRNE, I. DUURSMA, AND J. SHEEKEY Let us consider the average weight enumerator arrived at over all possible shortenings or puncturings (cf. [4]). Theorem 6.6. The average weight enumerators over all shortened and punctured codes of C are given by n () P(W C (x, y)) = W ΠH (C)(x, y), dim H=n n (2) S(W C (x, y)) = q n W Σh,H (C)(x, y). dim H=n, h H Proof. Let H be a hyperplane in F n q and let P H Fq n (n ) have columns that form a basis of H. From Lemma 6.2, rk(xp h ) = rk(x) if and only if X H. Let X C have rank w > 0. Then rk(xp H ) = w if and only if X H [ ] Since there are n (n w) w = hyperplanes in F n q containing X, n (n w) w X corresponds to a word of rank w in punctured codes of C and to a n w word of rank w in = q w n w punctured codes of C. Therefore, each term x n w y w of the weight enumerator of C yields the contribution q w n w x n w y w w + x n w y w, in the sum of the n different weight enumerators of all possible punctured codes. For any h / H, XP H Σ h,h (C) if and only if h X. There are X \X H such vectors h and so q n w one dimensional [ subspaces ] h such that XP H Σ h,h (C). As outlined above, there are q w n w hyperplanes not containing X. In particular each term x n w y w of the weight enumerator of C yields the contribution q n w q w n w x n w y w = q n n w,

21 RANK-METRIC CODES AND ZETA FUNCTIONS 2 in the sum of the q n n different weight enumerators of all possible shortened codes. Example 6.7. Consider the F 2 -[3 3, 4, 2] code C = 0 0 0, 0 0, 0 0, 0 0, which has weight enumerator W C (x, y) = x 3 + 3xy 2 + 2y 3. Then 6 hyperplanes H yield shortened codes with distribution of weight enumerators similar to those for H = 00, and the hyperplane H = 00 yields 4 shortened codes Σ h,h of order 2. H h W Σh,H (x, y) x 2 + 3y 2 x 2 + y 2 0 x 2 + y 2 0 x 2 + 3y 2 H h W Σh,H (x, y) 00 0 x 2 + y 2 00 x 2 + y 2 x 2 + y 2 0 x 2 + y 2 Then the average weight enumerator over all shortened codes is given by 3 (/4) W Σh,H (C)(x, y) dim H=2, h H = (/28)(28x y 2 ) = x 2 + 3/7y 2 ( ) = S(x 3 + 3xy 2 + 2y ) = x 2 + 3y 2. Lemma 6.8. The zeta polynomial P C (T ) is invariant under the shortening and puncturing operations S and P. Proof. From Corollary 6.4, we have ( n d ) S(W C (x, y)) = S p i (C)M n,d+i (x, y) = and P(W C (x, y)) = P ( n d n d p i (C)M n,d+i (x, y), ) n d p i (C)M n,d+i (x, y) = p i (C)M n,d+i (x, y).

22 22 I. BLANCO-CHACÓN, E. BYRNE, I. DUURSMA, AND J. SHEEKEY Definition 6.9. The normalized weight enumerator of C is defined to be the polynomial, n W C (T ) := (q m ) n W i (C)T i d. i We write WC P(T ) and WS C (T ) to denote the normalized weight enumerators corresponding to P(W C (x, y)) and S(W C (x, y)), respectively. Theorem 6.0. Let C have minimum distance d. Then, () W P C (T ) = qd T W C (qt ) + W C (T ) (q m ) W n (C)q n T n d+, (2) W S C (T ) = W C(T ) (q m ) W n (C)T n d. n Proof. To show (), we apply the operation P = (D q,x + D y ) to W C (x, y). ( ) n n n P(W C (x, y)) = x n + W i (C)q i n i n x n i y i i + W i (C) x n i y i ( n ) = x n n + W i (C)q i n i n x n i y i i + W i (C) x n i y i n ) = x n n + (W i (C)q i n i i + + W i+ (C) x n i y i. Therefore, the associated normalized polynomial WC P (T ) is given by: n ) = (q m ) n (W i (C)q i n i i + n + W i+ (C) T i d+ i n = (q m ) W i (C)q i n n + W i+(c) T i d+ i i + n = (q m ) W i (C)q i n T i d+ (q m ) q n W n (C)T n d+ + i n + (q m ) n W i (C) T i d, i = q d T W C (qt ) + W C (T ) (q m ) q n W n (C)T n d+. The proof of (2) is similar.

23 RANK-METRIC CODES AND ZETA FUNCTIONS 23 Consider the following operations α, ɛ on rational functions f(t ) in indeterminate T. (6.) αf(t ) := T f(t ) = and ɛf(t ) := f(qt ). (6.2) Now α, ɛ form a q-commuting pair [], and obey the relations ɛα = qαɛ, qα = αq, qɛ = ɛq, with respect to composition. Then α, ɛ are non-commuting operators with respect to the q-product determined by the relations in (6.2) and generate the Q-algebra α, ɛ, which acts on the space of rational functions. We thus express Theorem 6.0 () as (6.3) W P C (T ) ( + q d αɛ)w C (T ) = ( + q d ɛα)w C (T ) mod T n d+, An immediate corollary of Theorem 6.0 is as follows. Corollary 6.. Let C have minimum distance d. For 0 i d we have that, () ( + q d αɛ) i W S C (T ) ( + qd αɛ) i W C (T ) mod T n d, (2) ( + q d αɛ) i W P C (T ) ( + qd αɛ) i+ W C (T ) mod T n d+. In particular, for d S = d, d P = d and n P = n S = n we have () ( + q d αɛ) d S W S C (T ) ( + qd αɛ) d W C (T ) mod T n S d S +, (2) ( + q d αɛ) d P W P C (T ) ( + qd αɛ) d W C (T ) mod T n P d P +. We write WC Pr (T ) to denote the normalized weight enumerator that results from applying the puncturing operation r times to W C (T ). Note that the operators ( + q r αɛ) and ( + q s αɛ) commute, so there is no ambiguity in the following statement. Corollary 6.2. Let C have minimum distance d. Then W Pr C r (T ) = ( + q d j αɛ)w C (T ) mod T n d+ r q j(d r+j) r T j W C (q j T ). j Proof. The first equation is derived by repeated applications of (6.3). The well known identity r (+qj y) = r r 2) q(j y j, along with the fact that (αɛ) j = j

24 24 I. BLANCO-CHACÓN, E. BYRNE, I. DUURSMA, AND J. SHEEKEY q (j 2) α j ɛ j, yields the equation r ( + q j (q d r+ αɛ)) = Applying this to W C (T ) results in the 2nd equation. r q j(d r+j) r α j ɛ j. j Example 6.3. Let M m n,d (T ) be the normalized weight enumerator of an F 2 - [m n, m(n d + ), d] code. Then successively puncturing M 7 7,4 (T ), we arrive at the normalized weight enumerator of the whole space W F 7 4(T ) = M 7 4, (T ). 2 M 7 7,4 (T ) = + 98T T T 3 M 7 6,3 (T ) = + 4T T T 3 ( αɛ)m 7 7,4 (T ) mod T 4 = + 4T T T T 4 M 7 5,2 (T ) = + 22T T T 3 ( αɛ)m 7 6,3 (T ) mod T 4 = + 22T T T T 4 ( αɛ)( αɛ)m 7 7,4 (T ) mod T 4 = + 22T T T T T 5 W F 7 4(T ) = + 26T T T 3 2 ( αɛ)m 7 5,2 (T ) mod T 4 = + 26T T T T 4 ( αɛ)( αɛ)m 7 6,3 (T ) mod T 4 ( αɛ)( αɛ)( αɛ)m 7 7,4 (T ) mod T 4 = + 26T T T T T T 6. If we consider the action of α and ɛ on Q[T ]/ T s for some positive integer s, then they may be represented as s s rational matrices: α = and ɛ = diag(, q, q 2,..., q s ),

25 RANK-METRIC CODES AND ZETA FUNCTIONS 25 so that α, ɛ form a sub-algebra of Q s s. Then as an element of Q s s, α is nilpotent so the inverse of the operator ( + q t ɛα) exists and is given by ( + q t ɛα) = q t ɛα + (q t ɛα) ( ) s (q t ɛα) s. For example, with s = 3, we have + q t ɛα = 0 0 q t+ 0 and ( + q t ɛα) = 0 0 q t+ 0 0 q t+2 q 2t+3 q t+2 Recall the following formulae (see [9]): l (T ; q) l := ( q j T ) and (T ; q) l = l + j T j j. Then r ( + qj (q d r+ αɛ)) = ( q d r+ αɛ; q) r and so ( q d r+ αɛ; q) r = ( ) j [ r + j j Lemma 6.4. Let C have minimum distance d. Then ] q j(d r+)+(j 2) α j ɛ j. W C (T ) ( q d r+ αɛ; q) r WC Pr (T ) mod T n d+ = ( ) j r + j q j(d r+)+(j 2) T j WC Pr (q j T ). j In particular, puncturing is an invertible operation on normalized weight enumerators, modulo T n d+. Example 6.5. As i (qm q j l ) is the number of m l matrices of rank i i, we see that W F m l(t ) = q l i (q m q j )T i. i= j= As in Example 6.3, the full space F m l q is an MRD code with parameters [m l, ml, ] and can be obtained by successive puncturings of an MRD code. Let M(T ) be the normalized weight enumerator of an MRD [m (n + r), m(n d + ), d + r] code for some non-negative integer r. Then the

26 26 I. BLANCO-CHACÓN, E. BYRNE, I. DUURSMA, AND J. SHEEKEY previous observations along with Corollary 6.2 implies that M(T ) n d+ i= i (q m q j ) j= n d+ = (q m ) i= n d+ = (q m ) i= d+r ( + q d+r j αɛ) T i mod T n d+ q im (q m ; q) i ( q d r+ αɛ; q) r T i q im (q m ; q) i ( ) j r + j q j(d r+i)+(j 2) j T i+j. Lemma 6.6. Let C have minimum distance d. Then the normalized weight enumerator of the code C satisfies n d W C (T ) (q m ) n T i W d+i mod T n d+. i + d (T ; q) i+ Proof. We have T i (T ; q) i+ = [ i + j j ] T j+i = j T j, i which yields n d n T i n d n j W d+i = W d+i i + d (T ; q) i+ i + d i j=i n d [ = T j n j W d+i = T j n W d+j i + d i j + d j=i T j ] 7. Zeroes of the Zeta Polynomial For a self-dual code C, Theorem 5.4 shows that the (in general complex) zeroes of the zeta polynomial P C (T ) occur in pairs (α, /(q m α)). The two zeroes in a pair are of the same absolute value if and only if α = (q m ) /2. Writing ζ(s) = Z(T = (q m ) s ), this is the case if and only if as zeroes of ζ(s) they satisfy Re s = /2. The latter is the critical line for the classical Riemann zeta function as well as for the zeta function of a curve over a finite field, which we recall here for convenience of the reader.

27 RANK-METRIC CODES AND ZETA FUNCTIONS 27 If C is a non-singular projective curve defined over F q, for k, denote by N k the number of F q k-rational points of C, i.e., the cardinality of C(F q k). The zeta-function of C is the formal power series Z(C, T ) = exp N k k T k. k This expression is well defined as a formal power series (cf. [0, Appendix C]) and satisfies the following properties: Theorem 7. (Weil, Dwork). The zeta function of any non-singular projective curve of genus g can be expressed as Z(C, T ) = P (T ) ( T )( qt ), with P (T ) Z[T ] a polynomial of degree 2g, called the zeta-polynomial of C. Moreover, for each root ω i of P (T ) ( i 2g), we have ω = q /2. As we see, Z(C, T ) is primarily used as a generating function for the number of points on a curve. There is no immediate analogue of this property for linear codes. In another interpretation, the zeta function of a curve describes the growth rate for the dimensions of linear systems on the curve. This is analogous to describing the growth rate of the binomial moments of a linear code. For codes for which the growth rate of the binomial moments is close to the growth rate of dimensions of linear systems the zeroes of the zeta polynomial will lie on the critical line. For the Hamming metric, this occurs for remarkably many codes (more so when the field size is small), including for infinite families of extremal weight enumerators. The zeta polynomial defined in this paper measures the growth rate of binomial moments for rank-metric codes, and we may ask whether this growth rate is such that we can expect to find codes with zeroes of the zeta polynomial on the critical line. We give an example where the zeta polynomial is an exact match to the zeta polynomial of a (Hasse-Weil maximal) elliptic curve. Example 7.2. To construct a formally self-dual 4 4 rank-metric code we take the extended binary QR-code of length 8 in double-circulant form. We puncture at the last coordinate in one circulant block and shorten at the last coordinate in the other block. We then read out the remaining 6 coordinates in order and from each group of four coordinates form a row for a codeword in 4 4 format. The rank distribution for the 256 codewords is 0, 2 2, 3 62, The

28 28 I. BLANCO-CHACÓN, E. BYRNE, I. DUURSMA, AND J. SHEEKEY normalized binomial moments are,, 8/5, 6, 256 and the zeta polynomial is ( + 8T + 6T 2 )/25 = ( + 4T ) 2 /25. The zeroes T = /4 (twice) have absolute value (q m ) /2 = / 6 and thus lie on the critical line. The zeta polynomial is that of a maximal elliptic curve over the field of 6 elements. Remark 7.3. Note that constructing a self-dual rank-metric code from a self-dual linear code in this way is not unique; in other words, equivalent codes of length 6 may lead to inequivalent 4 4 rank-metric codes. The rank-weight distributions obtained need not bear any relations to each other. It requires further research as to what distributions, and hence what zeta polynomials, may occur from codes constructed in this way. Heuristics for the Hamming metric are that a sufficient condition for a matching growth rate, and thus for zeroes on the critical line, are that a (formally self-dual) code has large enough minimum distance and weight distribution close to that of a random code. Precise formulations of these observations and their verification are an open problem. We applied the same construction as in the example to double-circulant codes of length 38 to obtain formally self-dual 6 6 rank-metric codes. Among those, 2 codes had minimum rank distance 3 and thus a quadratic zeta polynomial. The rank distributions for these codes are not far apart and in all 2 cases the zeroes of the zeta polynomial are on the critical line. Among codes with lower minimum rank distance a wider range of rank distributions occurs. The distributions close to the average have zeroes on the critical line and those farther away have pairs of real zeroes. Typically one of the real zeroes among (α, /(q m α)) is close to and the other close to /q m. This agrees with observations for the Hamming metric. These preliminary observations suggest that the zeroes of zeta polynomials for the rank-metric have properties similar to those for the Hamming metric. It remains to make these properties precise and to establish to what extent or under which additional assumptions the properties are similar to those for curves. 7.. Self-dual divisible codes. Here we consider a family of codes which are divisible (that is, the rank-weight of every codeword is divisible by some integer c > ) and also formally self-dual. It is known that for codes in the Hamming metric, such codes can only exist for (q, c) {(2, 2), (2, 4), (3, 3), (4, 2)}. In [4] the zeta polynomials of self-dual divisible codes that are also random were considered. Here we show the existence of formally self-dual divisible codes for c = 2 and all values of q. However we note that this construction never leads to random codes.

29 RANK-METRIC CODES AND ZETA FUNCTIONS 29 The matrix space M m n (F q c) can be naturally embedded in the space M cm cn (F q ). The rank of every matrix in the image of this map is a multiple of c; in fact, we have rk Fq (A) = c rk Fq c (A). Thus from codes over an extension field, we can produce divisible codes easily. This is in contrast to the Hamming metric case, where the Hamming weight of a vector does not behave well when reducing the field of interest. A code C of dimension k over F q c produces a code C of dimension kc over F q. If we wish C to be formally self-dual, we must have k = mnc/2. As k mn, we need c = 2, and thus C is in fact the full space M m n (F q 2). For the image of the space M 2 2 (F q 2) in M 4 4 (F q ), we get the zeta polynomial P (T ) = + ( q4 + q 2 + q)t + q 4 T 2 q 2. + q + For example taking q = 2 we get the zeta polynomial P (T ) = ( 0T +6T 2 )/7. The discriminant of this polynomial is positive for all values of q 2, and so we get two real roots, whose product is q 4. For the image of the space M 3 3 (F q 2) in M 6 6 (F q ), we get the zeta polynomial P (T ) = + ( q6 + q 2 + q)t + ( q 8 q 7 + q 6 + q 4 + q 3 )T 2 + ( q 2 + q 8 + q 7 )T 3 + q 2 T 4 q 4 + q 3 + q 2 + q + For example taking q = 2 we get P (T ) = ( 58T 296T 2 372T T 4 )/5, which has two real roots and one pair of conjugate complex roots. In general, the image C of the space M m m (F q 2) in M 2m 2m (F q ) is a formally self-dual code with minimum rank distance d = 2 and all rank-weights divisible by c = 2. Let ρ = q m. The zeta polynomial P (T ) for C has two real roots. As m, the roots converge to and ρ 2. The rescaled polynomial p(t ) = P (ρt )/P (0) = + + T 2m 2 is self-reciprocal. It has two real roots that converge to ρ and ρ as m and 2m 4 complex roots that lie on the unit circle. As m, the factor p (T ) of p(t ) with two real roots converges to (ρ + ρ)t + T 2. The cofactor p 2 (T ) with complex roots converges to + T 2m 4. The coefficients of p(t ) approach 0 as m except for the coefficients at, T, T 2 and their reciprocals. The coefficients of p(t ) at, T and T 2 can be obtained with (4.6) by expressing the rank-weight enumerator of C as a linear combination of the F q rank-weight enumerators M 2m,d, for d 2. We find lim p(t ) = m (qm ρ(q 2 + q))t (q 2 + q ρ 2 (q 6 + q 4 + q 2 )T (q 2 + q ρ 2 (q 6 + q 4 + q 2 )T 2m 4 (q m ρ(q 2 + q))t 2m 3 + T 2m 2.

30 30 I. BLANCO-CHACÓN, E. BYRNE, I. DUURSMA, AND J. SHEEKEY 400 y x Figure. Complex zeroes for the image of M 9 (F 4 ) in M 8 (F 2 ). As m, and thus ρ = q m 0, p(ρ) 0. Thus p (T ) = (ρ + ρ)t + T 2 divides the limit and the cofactor p 2 (T ) satisfies lim p 2(T ) = + (q 2 + q + )ρt + (q 6 + q 4 + 2q 2 + q + )ρ 2 T m + (q 6 + q 4 + 2q 2 + q + )ρ 2 T 2m 6 + (q 2 + q + )ρt 2m 5 + T 2m 4. The limits converge fast. We find that the complex zeros of the zeta polynomial for the image of M m m (F q 2) in M 2m 2m (F q ) all have the same absolute value ρ and, already for small m, are close to the zeros of + (q m T ) 2m 4 = 0. References [] E. Byrne and A. Ravagnani, Covering radius of matrix codes endowed with the rank metric, SIAM Journal on Discrete Mathematics, to appear, 207. [2] P. Delsarte, Bilinear forms over a finite field, with applications to coding theory, J. Combin. Theory Ser. A, 25(3):226 24, 978. [3] J. G. Dumas, R. Gow, G. McGuire, and J. Sheekey, Subspaces of matrices with special rank properties, Linear Algebra Appl., 433():9 202, 200. [4] I. Duursma, From weight enumerators to zeta functions, Discrete Appl. Math., (- 2):55 73, 200. [5] I. Duursma, Weight distributions of geometric Goppa codes, Trans. Amer. Math. Soc., 35(9): , 999.

31 RANK-METRIC CODES AND ZETA FUNCTIONS 3 [6] I. Duursma, Combinatorics of the two-variable zeta function, in Finite fields and applications, volume 2948 of Lecture Notes in Comput. Sci., pages Springer, Berlin, [7] E. M. Gabidulin, Theory of codes with maximum rank distance, Problemy Peredachi Informatsii, 2():3 6, 985. [8] M. Gadouleau and Z. Yan, Macwilliams identity for the rank metric, CoRR,abs/cs/070097, [9] G. Gasper and R. Mizan, Basic Hypergeometric Series, Volume 96 of Encyclopedia of Mathematics and its Applications. Cambridge University Press, Cambridge, second edition, [0] R. Hartshorne, Algebraic geometry, Springer-Verlag, New York-Heidelberg, 977. Graduate Texts in Mathematics, No. 52. [] T. H. Koornwinder, Special functions and q-commuting variables, In Special Functions, q-series and Related Topics, Fields Institute Communications, pages 09 36, American Mathematical Society, Rhode Island, 997. [2] A. Ravagnani, Rank-metric codes and their duality theory, Des. Codes Cryptogr., 80():97 26, 206. [3] Ron M. Roth, Maximum-rank array codes and their application to crisscross error correction, IEEE Trans. Inform. Theory, 37(2): , 99. I. Blanco-Chacón, E. Byrne, and J. Sheekey School of Mathematics and Statistics, University College Dublin I. Duursma Deptartment of Mathematics, University of Illinois Urbana-Champaign