Asymptotic bounds for the sizes of constant dimension codes and an improved lower bound
Daniel Heinlein, Sascha Kurz

TL;DR
This paper investigates asymptotic bounds for constant dimension codes used in network coding, reviews existing bounds, and introduces an improved linkage construction that surpasses previous bounds for certain code parameters.
Contribution
It provides a refined linkage construction that yields larger constant dimension codes, improving lower bounds beyond the MRD bound for specific parameters.
Findings
Improved lower bounds for constant dimension codes.
Relations between known upper bounds are clarified.
Construction exceeds MRD bounds for certain code parameters.
Abstract
We study asymptotic lower and upper bounds for the sizes of constant dimension codes with respect to the subspace or injection distance, which is used in random linear network coding. In this context we review known upper bounds and show relations between them. A slightly improved version of the so-called linkage construction is presented which is e.g. used to construct constant dimension codes with subspace distance , dimension of the codewords for all field sizes , and sufficiently large dimensions of the ambient space, that exceed the MRD bound, for codes containing a lifted MRD code, by Etzion and Silberstein.
| bklb | mrdb | bkub | lold | lnew | ea | |
|---|---|---|---|---|---|---|
| 6 | 77 | 71 | 77 | 65 | 65 | |
| 7 | 333 | 291 | 381 | 257 | 265 | 301 |
| 8 | 1326 | 1179 | 1493 | 1033 | 1101 | 1117 |
| 9 | 5986 | 4747 | 6205 | 4929 | 4929 | 4852 |
| 10 | 23870 | 19051 | 24698 | 21313 | 21313 | 18924 |
| 11 | 97526 | 76331 | 99718 | 85249 | 85257 | 79306 |
| 12 | 385515 | 305579 | 398385 | 383105 | 383105 | 309667 |
| 13 | 1597245 | 1222827 | 1597245 | 1532417 | 1532425 | 1287958 |
| 14 | 6241665 | 4892331 | 6387029 | 6241665 | 6241665 | 4970117 |
| 15 | 24966665 | 19571371 | 25562941 | 24966657 | 24966665 | 20560924 |
| 16 | 102223681 | 78289579 | 102243962 | 102223681 | 102223681 | 79608330 |
| 17 | 408894729 | 313166507 | 409035142 | 408894721 | 408894729 | |
| 18 | 1635578957 | 1252682411 | 1636109361 | 1635578889 | 1635578957 | |
| 19 | 6542315853 | 5010762411 | 6544674621 | 6542315597 | 6542315853 | 5200895489 |
| bklb | mrdb | bkub | lold | lnew | ea | |
|---|---|---|---|---|---|---|
| 6 | 1.203125 | 1.109375 | 1.203125 | 1.015625 | 1.015625 | |
| 7 | 1.300781 | 1.136719 | 1.488281 | 1.003906 | 1.035156 | 1.175781 |
| 8 | 1.294922 | 1.151367 | 1.458008 | 1.008789 | 1.075195 | 1.090820 |
| 9 | 1.461426 | 1.158936 | 1.514893 | 1.203369 | 1.203369 | 1.184570 |
| 10 | 1.456909 | 1.162781 | 1.507446 | 1.300842 | 1.300842 | 1.155029 |
| 11 | 1.488129 | 1.164719 | 1.521576 | 1.300797 | 1.300919 | 1.210114 |
| 12 | 1.470623 | 1.165691 | 1.519718 | 1.461430 | 1.461430 | 1.181286 |
| 13 | 1.523252 | 1.166179 | 1.523252 | 1.461427 | 1.461434 | 1.228292 |
| 14 | 1.488129 | 1.166423 | 1.522786 | 1.488129 | 1.488129 | 1.184968 |
| 15 | 1.488129 | 1.166545 | 1.52367 | 1.488129 | 1.488129 | 1.225527 |
| 16 | 1.523252 | 1.166606 | 1.523554 | 1.523252 | 1.523252 | 1.186257 |
| 17 | 1.523252 | 1.166636 | 1.523775 | 1.523252 | 1.523252 | |
| 18 | 1.523252 | 1.166651 | 1.523746 | 1.523252 | 1.523252 | |
| 19 | 1.523252 | 1.166659 | 1.523801 | 1.523252 | 1.523252 | 1.210928 |
| bklb | mrdb | bkub | lold | lnew | ea | |
|---|---|---|---|---|---|---|
| 6 | 1.084507 | 1.0 | 1.084507 | 0.915493 | 0.915493 | |
| 7 | 1.144330 | 1.0 | 1.309278 | 0.883162 | 0.910653 | 1.034364 |
| 8 | 1.124682 | 1.0 | 1.266327 | 0.876166 | 0.933842 | 0.947413 |
| 9 | 1.261007 | 1.0 | 1.307141 | 1.038340 | 1.038340 | 1.022119 |
| 10 | 1.252953 | 1.0 | 1.296415 | 1.118734 | 1.118734 | 0.993334 |
| 11 | 1.277672 | 1.0 | 1.306389 | 1.116833 | 1.116938 | 1.038975 |
| 12 | 1.261589 | 1.0 | 1.303705 | 1.253702 | 1.253702 | 1.013378 |
| 13 | 1.306190 | 1.0 | 1.306190 | 1.253176 | 1.253182 | 1.053263 |
| 14 | 1.275806 | 1.0 | 1.305519 | 1.275806 | 1.275806 | 1.015900 |
| 15 | 1.275673 | 1.0 | 1.306140 | 1.275672 | 1.275673 | 1.050561 |
| 16 | 1.305712 | 1.0 | 1.305972 | 1.305712 | 1.305712 | 1.016845 |
| 17 | 1.305678 | 1.0 | 1.306127 | 1.305678 | 1.305678 | |
| 18 | 1.305661 | 1.0 | 1.306085 | 1.305661 | 1.305661 | |
| 19 | 1.305653 | 1.0 | 1.306124 | 1.305653 | 1.305653 | 1.037945 |
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Asymptotic bounds for the sizes of constant dimension codes and an improved lower bound
Daniel Heinlein and Sascha Kurz
University of Bayreuth
[email protected] The work was supported by the ICT COST Action IC1104 and grants KU 2430/3-1, WA 1666/9-1 – “Integer Linear Programming Models for Subspace Codes and Finite Geometry” – from the German Research Foundation.
Abstract
We study asymptotic lower and upper bounds for the sizes of constant dimension codes with respect to the subspace or injection distance, which is used in random linear network coding. In this context we review known upper bounds and show relations between them. A slightly improved version of the so-called linkage construction is presented which is e.g. used to construct constant dimension codes with subspace distance , dimension of the codewords for all field sizes , and sufficiently large dimensions of the ambient space, that exceed the MRD bound, for codes containing a lifted MRD code, by Etzion and Silberstein.
Keywords: constant dimension codes, subspace distance, injection distance, random network coding
1 Introduction
Let be a -dimensional vector space over the finite field with elements. By we denote the set of all -dimensional subspaces in , where , which has size . More general, the set of all subspaces of forms a metric space with respect to the subspace distance defined by , see [32], and the injection distance defined by , see [40]. Coding Theory on is motivated by Kötter and Kschischang [32] via error correcting random network coding, see [4]. In this context it is natural to consider codes where each codeword, i.e., each element of , has the same dimension , called constant dimension code (cdc), since this knowledge can be exploited by decoders. For constant dimension codes we have , so that we will only consider the subspace distance in this paper. By we denote a cdc in with minimum (subspace) distance and size , where the dimensions of each codeword is . As usual, a cdc has the minimum distance , if for all and equality is attained at least once. If , we set the minimum distance to . The corresponding maximum size is denoted by , where we allow the minimum distance to be larger than . In [32] the authors provided lower and upper bounds for which are less than a factor of apart. For increasing field size this factor tends to . Here, we tighten the corresponding analysis and end up in a factor of less than for the binary field and a strictly better factor for larger values of . With respect to lower bounds, we slightly generalize the so-called linkage construction by Gluesing-Luerssen, Troha / Morrison [22, 21] and Silberstein, (Horlemann-)Trautmann [38]. This improvement then gives the best known lower bounds for for many parameters, cf. the online tables http://subspacecodes.uni-bayreuth.de associated with [23]. For codes containing a lifted maximum rank distance (LMRD) code as a subcode an upper bound on the size has been presented in [16] for some infinite series of parameters. Codes larger than this MRD bound are very rare. Based on the improved linkage construction we give an infinite series of such examples.
In this context we mention the following asymptotic result based on the non-constructive probabilistic method. If the subspace distance and the dimension of the codewords is fixed, then the ratio of the lower and upper bound tends to as the dimension of the ambient space approaches infinity, see [18, Theorem 4.1], which is implied by a more general result of Frankl and Rödl on hypergraphs. The same result, with an explicit error term, was also obtained in [8, Theorem 1]. If and is fixed we have the same result due to the orthogonal code. If the parameter can vary with the dimension , then our asymptotic analysis implies there is still a gap of almost between the lower and the upper bound of the code sizes for and , which is the worst case.
The remaining part of the paper is organized as follows. In Section 2 we collect the basic facts and definitions for constant dimension codes. Upper bounds on the achievable codes sizes are reviewed in Section 3. Here, we partially extend the current knowledge on the relation between these bounds. While most of them are known around 2008 there are some recent improvements for the subclass of partial spreads, where , which we summarize in Subsection 3.1. In Section 4 we present the mentioned improvement of the linkage construction. Asymptotic bounds for the ratio between lower and upper bounds for code sizes are studied in Section 5. We continue with the upper bound for constant dimension codes containing a lifted MRD code in Section 6, including some numerical results, before we draw a short conclusion in Section 7.
2 Preliminaries
For the remainder of the paper we set , where is a prime power. By we denote the dimension of . Using the language of projective geometry, we will call the -dimensional subspaces of points and the -dimensional subspaces lines. First, we observe that the -binomial coefficient indeed gives the cardinality of . To this end, we associate with a subspace a unique matrix in row reduced echelon form (rref) having the property that and denote the corresponding bijection
[TABLE]
by . An example is given by , where is a line that contains the three points , , and . Counting those matrices gives
[TABLE]
for all integers . Especially, we have . Given a non-degenerate bilinear form, we denote by the orthogonal subspace of a subspace , which then has dimension . Then, we have , so that . The recurrence relation for the usual binomial coefficients generalize to . In order to remove the restriction , we set for and , whenever or . This extension goes in line with the interpretation of the number of -dimensional subspaces of and respects the orthogonality relation. In order to write for positive integers in later formulas, we apply the definition of also in cases where is not a prime power and set for .
Using the bijection we can express the subspace distance between two -dimensional subspaces via the rank of a matrix:
[TABLE]
Using as vertex set, we obtain the so-called Grassmann graph, where two vertices are adjacent iff the corresponding subspaces intersect in a space of dimension . It is well-known that the Grassmann graph is distance regular. The injection distance corresponds to the graph distance in the Grassmann graph. Considered as an association scheme one speaks of the -Johnson scheme.
If is a cdc with minimum subspace distance , we speak of a constant dimension code. In the special case of one speaks of so-called partial spreads, i.e., collections of -dimensional subspaces with pairwise trivial intersection.
Besides the injection and the subspace distance we will also consider the Hamming distance , for two vectors , and the rank distance , for two matrices . The latter is indeed a metric, as observed in [20]. A subset is called a rank metric code. If the minimum rank-distance of is given by , we will also speak of an rank metric code in order to specify its parameters. A rank metric code is called linear if forms a subspace of , which implies that has to be a power of the field size .
Theorem 1**.**
*(see [20])
Let be positive integers, a prime power, and be a rank metric code with minimum rank distance . Then, .*
Codes attaining this upper bound are called maximum rank distance (MRD) codes. They exist for all (suitable) choices of parameters, which remains true if we restrict to linear rank metric codes, see [20]. If or , then only is possible, which can be achieved by a zero matrix and may be summarized to the single upper bound . Using an identity matrix as a prefix one obtains the so-called lifted MRD codes.
Theorem 2**.**
[39, Proposition 4]** For positive integers with , , and even, the size of a lifted MRD code in with subspace distance is given by
[TABLE]
If , then we have .
The Hamming distance can be used to lower bound the subspace distance between two codewords (of the same dimension). To this end let denote the pivot positions of the matrix in rref. For our example we we have . Slightly abusing notation we also write for subspaces instead of .
Lemma 1**.**
[15*, Lemma 2]**
For two subspaces , we have .*
3 Upper Bounds
In this section we review and compare known upper bounds for the sizes of constant dimension codes. Here we assume that , , and are integers with , , and even in all subsequent results. The bound just ensures that is non-empty. Note that and is even for all . Restricting to the set case, we trivially have for or , so that we assume and , which then implies and . We remark that some of the latter bounds are also valid for parameters outside the ranges of non-trivial parameters considered by us. Since the maximum size of a code with certain parameters is always an integer and some of the latter upper bounds can produce non-integer values, we may always round them down. To ease the notation we will commonly omit the final rounding step.
The list of known bounds has not changed much since [29], see also [17]. Comparisons of those bounds are scattered among the literature and partially hidden in comments, see e.g. [6]. Additionally some results turn out to be wrong or need a reinterpretation at the very least.
Counting -dimensional subspaces having a large intersection with a fixed -dimensional subspace gives:
Lemma 2**.**
For integers and we have
[TABLE]
where and .
Proof.
Let us denote by , where . With this, the number of choices for is given by
[TABLE]
Finally apply the convention for integers with or .
∎
Note that is equivalent to . The fact that the Grassmann graph is distance-regular implies:
Theorem 3**.**
(Sphere-packing bound)[32, Theorem 6]
[TABLE]
We remark, that we can obtain the denominator of the formula of Theorem 3 by setting , in Lemma 2 and applying . The right hand side is symmetric with respect to orthogonal subspaces, i.e., the mapping leaves it invariant.
By defining a puncturing operation one can decrease the dimension of the ambient space and the codewords. Since the minimum distance decreases by at most two, we can iteratively puncture times, so that since . Considering either the code or its orthogonal code gives:
Theorem 4**.**
(Singleton bound)[32, Theorem 9]
[TABLE]
Referring to [32] the authors of [29] state that even a relaxation of the Singleton bound is always stronger than the sphere packing bound for non-trivial codes. However, for , , , and , the sphere-packing bound gives an upper bound of while the Singleton bound gives an upper bound of . For , , , and it is just the other way round, i.e., the Singleton bound gives and the sphere-packing bound gives . Examples for the latter case are easy to find. For both bounds coincide and for the Singleton bound is always stronger than the sphere-packing bound since . The asymptotic bounds [32, Corollaries 7 and 10], using normalized parameters, and [32, Figure 1] suggest that there is only a small range of parameters where the sphere-packing bound can be superior to the Singleton bound.111By a tedious computation one can check that the sphere-packing bound is strictly tighter than the Singleton bound iff , and .
Given an arbitrary metric space , an anticode of diameter is a subset whose elements have pairwise distance at most . Since the -Johnson scheme is an association scheme the Anticode bound of Delsarte [11] can be applied. As a standalone argument we go along the lines of [2] and consider bounds for codes on transitive graphs. By double-counting the number of pairs , where , we obtain:
Lemma 3**.**
[2, Lemma 1]**, cf. [3, Theorem 1’] Let be a graph that admits a transitive group of automorphisms and let be arbitrary subsets of the vertex set . Then, there exists a group element such that
[TABLE]
Corollary 1**.**
[2, Corollary 1]**, cf. [3, Theorem 1] Let be a code with (injection or graph) distances from . Then, for an arbitrary subset there exists a code with distances from such that
[TABLE]
If is a constant dimension code with minimum injection distance , i.e., , and is an anticode with diameter , we have , so that we obtain Delsarte’s Anticode bound
[TABLE]
The set of all elements of which contain a fixed -dimensional subspace is an anticode of diameter with elements. By orthogonality, the set of all elements of which are contained in a fixed -dimensional subspace is also an anticode of diameter with elements. Frankl and Wilson proved in [19, Theorem 1] that these anticodes have the largest possible size, which implies:
Theorem 5**.**
(Anticode bound)
[TABLE]
Using different arguments, Theorem 5 was proved in [42, Theorem 5.2] by Wang, Xing, Safavi-Naini in 2003. Codes that can achieve the (unrounded) value are called Steiner structures. It is a well-known and seemingly very hard problem to decide whether a Steiner structure for , , and exists. For the best known bounds are . Additionally it is known that a code attaining the upper bound can have automorphisms of at most order , see [30]. So far, the only known Steiner structure corresponds to [9]. The reduction to Delsarte’s Anticode bound can be found e.g. in [17, Theorem 1].
Since the sphere underlying the proof of Theorem 3 is also an anticode, Theorem 3 is implied by Theorem 5. For both bounds coincide. In [43, Section 4] Xia and Fu verified that the Anticode bound is always stronger than the Singleton bound for the ranges of parameters considered by us.
Mimicking a classical bound of Johnson on binary error-correcting codes with respect to the Hamming distance, see [28, Theorem 3] and also [41], Xia and Fu proved:
Theorem 6**.**
*(Johnson type bound I) [43, Theorem 2]
If , then*
[TABLE]
However, the required condition of Theorem 6 is rather restrictive and can be simplified considerably.
Proposition 1**.**
For , the bound in Theorem 6 is applicable iff and . Then, it is equivalent to
[TABLE]
Proof.
If we have , so that we assume in the following. If and , then
[TABLE]
If and , then
[TABLE]
If , , and , then it can be easily checked that the condition of Theorem 6 is satisfied and we obtain the proposed formula after simplification.
∎
For Theorem 6 gives which is trivially satisfied with equality. In Subsection 3.1 we will provide tighter upper bounds for the special case where , i.e., partial spreads. Indeed, the bound stated in Proposition 1 corresponds to the most trivial upper bounds for partial spreads that is tight iff divides , as we will see later on. So, due to orthogonality, Theorem 6 is dominated by the partial spread bounds discussed later on.
While the previously mentioned generalization of a classical bound of Johnson on binary error-correcting codes yields the rather weak Theorem 6, generalizing [28, Inequality (5)], see [43] yields a very strong upper bound:
Theorem 7**.**
(Johnson type bound II) [43, Theorem 3], [17, Theorem 4,5]
[TABLE]
Note that for Inequality (3) gives since we have by definition. Similarly, for , Inequality (4) gives .
Some sources like [43, Theorem 3] list just Inequality 3 and omit Inequality 4. This goes in line with the treatment of the classical Johnson type bound II for binary error-correcting codes, see e.g. [35, Theorem 4 on page 527], where the other bound is formulated as Problem (2) on page 528 with the hint that ones should be replaced by zeros. Analogously, we can consider orthogonal codes:
Proposition 2**.**
Inequality (3) and Inequality (4) are equivalent using orthogonality, cf. [17, Section III, esp. Lemma 13].
Proof.
We have
[TABLE]
which is Inequality (4), and
[TABLE]
which is Inequality (3).
∎
Of course, the bounds in Theorem 7 can be applied iteratively. In the classical Johnson space the optimal order of the corresponding inequalities is unclear, see e.g. [35, Research Problem 17.1]. Denoting the maximum size of a binary constant-weight block code of length , Hamming distance and weight by , the two corresponding variants of the inequalities in Theorem 7 are and . Applying the first bound yields
[TABLE]
while applying the second bound yields
[TABLE]
using the numerical bounds from
http://webfiles.portal.chalmers.se/s2/research/kit/bounds/cw.html, cf. [1].
The authors of [17, 29] state that the optimal choice of Inequality (3) or Inequality (4) is unclear, too. However, this question is much easier to answer for constant dimension codes.
Proposition 3**.**
For we have
[TABLE]
where equality holds iff .
Proof.
By considering orthogonal codes we obtain equality for . Now we assume and show
[TABLE]
which implies the proposed statement. Considering the size of the LMRD code we can lower bound the right hand side of Inequality (5) to
[TABLE]
Since
[TABLE]
we can use the Anticode bound to upper bound the left hand side of Inequality (5) to
[TABLE]
where . Thus, it suffices to verify
[TABLE]
where we have divided by
[TABLE]
Since , we have . Since and , we have . Since , we have , which is at most for . Thus, Inequality (6) is valid for all .
If and , then and , so that Inequality (6) is satisfied.
In the remaining part of the proof we assume and . If , then and . If , then and . If , then , , and due to . Thus, Inequality (6) is valid in all cases.
∎
Knowing the optimal choice between Inequality (3) and Inequality (4), we can iteratively apply Theorem 7 in an ideal way initially assuming :
Corollary 2**.**
(Implication of the Johnson type bound II)
[TABLE]
We remark that this upper bound is commonly stated in an explicit version, where is inserted, see e.g. [17, Theorem 6], [29, Theorem 7], and [43, Corollary 3]. However, currently much better bounds for partial spreads are available.
It is shown in [43] that the Johnson bound of Theorem 7 improves on the Anticode bound in Theorem 5, see also [6]. To be more precise, removing the floors in the upper bound of Corollary 2 and replacing by gives
[TABLE]
which is the right hand side of the Anticode bound for . So, all upper bounds mentioned so far are (weakly) dominated by Corollary 2, if we additionally assume . As a possible improvement [2, Theorem 3] was mentioned as [29, Theorem 8]. Here, we correct typos and give a slightly enlarged proof, thanks to a personal communication with Aydinian.
Theorem 8**.**
[2, Theorem 3]** For integers , , and we have
[TABLE]
Proof.
Let be a fixed subspace with and define
[TABLE]
so that is given by Lemma 2. Consider a code and take noting that the latter has a minimum distance of at least . Two arbitrary codewords have distance , where we write and for integers . Replacing each codeword of by an arbitrary -dimensional subspace, we obtain a cdc with a minimum distance of at least . Since we have , so that Corollary 1 gives the proposed upper bound.
∎
As Theorem 8 has quite some degrees of freedom, we partially discuss the optimal choice of parameters. For and , we obtain , which is the -fold iteration of Inequality (4) of the Johnson bound (without rounding). Thus, is the best choice for , yielding a bound that is equivalent to Inequality (4). For and the bound can be rewritten to , see the proof of Proposition 4. For the bound remains valid but is strictly weaker than for . Choosing gives the trivial bound . For the range of parameters , , limited facing nerve-jangling numerical pitfalls, and , where is of course a prime power and is even, the situation is as follows. If , there are no proper improvements with respect to Theorem 7. For the case , i.e., partial spreads treated in the next subsection, we have some improvements compared to which is the most trivial bound for partial spreads. Within our numerical range, most of them are covered by the following proposition, where we apply Theorem 8 with and to . The other cases are due to the fact that Theorem 14 is tighter than Theorem 16 for larger values of . In no case a proper improvement with respect to the tighter bounds from the next subsection emerged.
Proposition 4**.**
For and we have
[TABLE]
Proof.
Note that implies . The left hand side simplifies to
[TABLE]
Then we apply Theorem 16 with , , and , which yields for .
∎
We remark that applying Theorem 14 and Theorem 16 directly is at least as good as the application of Theorem 8 with and for .
The Delsarte linear programming bound for the -Johnson scheme was obtained in [12]. However, numerical computations indicate that it is not better than the Anticode bound, see [6]. For , i.e., the non-partial spread case, besides the stated bound only the following two specific bounds, based on extensive computer calculations, are known:
Theorem 9**.**
[26, Theorem 1]**
Proposition 5**.**
[24]**
As the authors of [24] have observed, the Johnson bound of Theorem 7 does not improve upon Corollary 2 when applied to Theorem 9 or Proposition 5.
If we additionally restrict ourselves to constant dimension codes, that contain a lifted MRD code, another upper bound is known:
Theorem 10**.**
[16, Theorem 10 and 11]** Let be a constant dimension code, with and minimum subspace distance , that contains a lifted MRD code.
- •
If and , then ;
- •
if , where is even, then .
3.1 Upper bounds for partial spreads
The case of constant dimension codes with maximum possible subspace distance is known under the name partial spreads. Counting points, i.e., -dimensional subspaces, in and gives the obvious upper bound . In the case of equality one speaks of spreads, for which a handy existence criterion is known from the work of Segre in 1964.
Theorem 11**.**
[37, §VI]** contains a spread if and only if is a divisor of .
If is not a divisor of , far better bounds are known including some recent improvements, which we will briefly summarize. For a more detailed treatment we refer to e.g. [27]. The best known parametric construction was given by Beutelspacher in 1975:
Theorem 12**.**
[7*]**
For positive integers satisfying , and we have with equality for .*
The determination of for was achieved more than 30 years later in [14] and continued to for and arbitrary in [34]. Besides the parameters of , for , see [14] for an example showing , no partial spreads exceeding the lower bound from Theorem 12 are known.
For a long time the best known upper bound on was the one obtained by Drake and Freeman in 1979:
Theorem 13**.**
[13*, Corollary 8]**
If with , then*
[TABLE]
where .
Quite recently this bound has been generalized to:
Theorem 14**.**
[33*, Theorem 2.10]**
For integers , , , with , , , , and , we have .*
The construction of Theorem 12 is asymptotically optimal for , as recently shown by Năstase and Sissokho:
Theorem 15**.**
[36*, Theorem 5]**
Suppose with and . If then .*
Applying similar techniques, the result was generalized to :
Theorem 16**.**
[33, Theorem 2.9]** For integers , , , and with we have , where and .
Using Theorem 14 the restriction can be removed from Theorem 16, see [27].
Currently, Theorem 11, Theorem 14, and Theorem 16 constitute the tightest parametric bounds for . The only known improvements, by exactly one in every case, are given by the specific bounds stated in [33], which are based on the linear programming method applied to projective -divisible linear error-correcting codes over with respect to the Hamming distance, see [27]. As this connection seemed to be overlooked before, it may not be improbable that more sophisticated methods from classical coding theory can improve further values, which then imply improved upper bounds for constant dimension codes via the Johnson bound of Theorem 7.
4 The linkage construction revisited
A very effective and widely applicable construction of constant dimension codes was stated by Gluesing-Luerssen and Troha:
Theorem 17**.**
[22*, Theorem 2.3]**, cf. [38, Corollary 39]
Let be a constant dimension code for and let be a linear rank metric code. Then*
[TABLE]
is a constant dimension code.
Here denotes the concatenation of two matrices with the same number of rows and denotes the -matrix consisting entirely of zeros. The resulting code depends on the choice of the codes , , and their representatives within isomorphism classes, so that one typically obtains many isomorphism classes of codes with the same parameters.
We remark that [38, Theorem 37] corresponds to the weakened version of Theorem 17 where the codewords from the cdc are not taken into account, cf. [21, Theorem 5.1]. In [38, Corollary 39] Silberstein and (Horlemann-)Trautmann obtain the same lower bound, assuming , which is indeed the optimal choice, and .222It can be verified that for the optimal choice of in [38, Corollary39] is given by . In that case the construction is essentially the union of an LMRD code with an code . Note that for the constructed code is an embedded code .
The main idea behind Theorem 17 is to consider two sets of codewords with disjoint pivot vectors across the two sets and to utilize the interplay between the rank and the subspace distance for a product type construction. Using Lemma 1 the restriction of the disjointness of the pivot vectors can be weakened, which gives the following improvement:
Theorem 18**.**
Let be a constant dimension code for , and let be a linear rank metric code. Then
[TABLE]
is a constant dimension code.
Proof.
The dimension of the ambient space and the codewords of directly follow from the construction. Since the constructed matrices all are in rref and pairwise distinct, is well defined and we have . It remains to lower bound the minimum subspace distance of .
Let and . If , we have
[TABLE]
using Equation (1) in the first step. If but , we have
[TABLE]
For applying Equation (1) gives
[TABLE]
Last, for two codewords and , we can use the shape of the pivot vectors and apply Lemma 1. The pivot vector has its ones in the first positions and the pivot vector has its ones not in the first positions, so that the ones can coincide at most at the positions . Thus, . Lemma 1 then gives .
∎
An example where Theorem 18 yields a larger code than Theorem 17 is e.g. given for the parameters , , , and . In order to apply Theorem 17 we have to choose , , and . For we obtain and . Since the size of the rank metric code is bounded by , the constructed code has a size of at most . For the roles of and interchange. Since the size of the rank metric code is bounded by , the constructed code has a size of at most . In Theorem 18 we can choose , so that we can drop one column of the zero matrix preceding the matrices of the second set of codewords, i.e., . Choosing and we can achieve and . Since the size of the rank metric code can attain we can construct a code of size . While for these parameters sill larger codes are known, the situation significantly changes in general. Considering the range of parameters , , and , where is of course a prime power and is even, Theorem 17 provides the best known lower bound for in % of the cases, while Theorem 18 provides the best known lower bound in % of the cases. Since the sizes of both constructions can coincide, the sum of both fractions gives more than %. In just % of the cases strictly superior constructions are known compared to Theorem 18, where most of them arose from the so-called Echelon-Ferrers construction or one of their variants, see [23] and the corresponding webpage.333Entries of type improved_linkage(m) correspond to Corollary 4 with chosen as parameter.
If one is interested in codes of large size, then one should choose the parameters , , , and , in Theorem 18, as small as possible in order to maximize the sizes , , and , i.e., we can assume . Moreover, the codes , , and should have the maximum possible size with respect to their specified parameters. For the maximum possible size is and for the maximum possible size is , where .
Corollary 3**.**
For positive integers and we have .
Instead of or we may also insert any lower bound of these commonly unknown values. By a variable transformation we obtain:
Corollary 4**.**
For positive integers and we have .
For the parameters of spreads the optimal choice of the parameter in Corollary 4 can be determined analytically:
Lemma 4**.**
If and divides , then Corollary 4 gives for all and smaller sizes otherwise.
Proof.
Using for all integers being divisible by , we obtain
[TABLE]
if divides . Otherwise, gives a lower bound.
∎
We remark that the tightest implications of Corollary 4 can be evaluated by dynamic programming. To this end we consider fixed parameters , , and use the abbreviations and for integers , so that the inequality of Corollary 4 reads
[TABLE]
For a given maximal value we initialize the values for by the best known lower bounds for from other constructions. Then we loop over from to and eventually replace by
[TABLE]
By an arithmetic progression we can use (7) in order to obtain a lower bound for given just two initial -values.
Proposition 6**.**
For positive integers , , and , we have
[TABLE]
If additionally, and , then we have
[TABLE]
Proof.
Using Inequality (7) with and gives
[TABLE]
By induction, we obtain
[TABLE]
for all .
For the second part, applying Inequality (7) with and gives
[TABLE]
By induction, we obtain
[TABLE]
for all .
If and , then
[TABLE]
so that
[TABLE]
∎
Example 1**.**
Using [9] and [23], applying Proposition 6 with gives
[TABLE]
and
[TABLE]
for all .
In the next section we will see that the first lower bound almost meets the Anticode bound.
We remark that Theorem 18 can be easily generalized to a construction based on a union of sets of codewords.
Corollary 5**.**
For positive integers , , and let
- •
* be an constant dimension code,*
- •
, ,
- •
* be a linear rank metric code, where and ,*
- •
, , , and .
Then
[TABLE]
is a constant dimension code with
- •
,
- •
, and
- •
.
Proof.
We prove by inductively applying Theorem 18 times. Denote
[TABLE]
for , i.e., is a padded constant dimension code. Applying Theorem 18 for and with yields a constant dimension code. If the first codes, yield an constant dimension code , then performing Theorem 18 for this code and with yields an constant dimension code.
∎
Since the proof uses multiple applications of Theorem 18 this code can be found by the dynamic programming approach based on Theorem 18, i.e., Corollary 5 is redundant. However, it can be used to prove:
Corollary 6** (cf. [22, Theorem 4.6]).**
Let be an linear MRD code, where , for and let be an constant dimension codes for . Then
[TABLE]
is a constant dimension code. Note that implies .
Proof.
Applying Corollary 5 with
- •
- •
, ,
- •
(i.e., an constant dimension code)
- •
- •
- •
(i.e., an rank metric code)
- •
an MRD code
yields an constant dimension code:
[TABLE]
∎
We remark that is a linear MRD code, since each codeword has . The other direction is not necessarily true, e.g., , where is a non-zero column, cannot be split in two matrices and both having rank distance at least for . Hence, this corollary constructs codes of the same size as Theorem 4.6 in [22] but these codes are not necessarily equal.
5 Asymptotic bounds
Kötter and Kschischang have stated the bounds
[TABLE]
for in [32, Lemma 4] for the -binomial coefficients. They used this result in order to prove that the lifted MRD codes, they spoke about linearized polynomials, have at least a size of a quarter of the Singleton bound if tends to infinity. Actually, they have derived a more refined bound, which can best expressed using the so called -Pochhammer symbol specializing to :
[TABLE]
where denotes , cf. the estimation for the Anticode bound in the proof of Proposition 3. The sequence is monotonically increasing with and approaches for large , see e.g. [32] and [29] for some numerical values.
Lemma 5**.**
For each we have .
Proof.
Plugging in the definition of the -binomial coefficient, we obtain
[TABLE]
∎
Using this asymptotic result we can compare the size of the lifted MRD codes to the Singleton and the Anticode bound.
Proposition 7**.**
For the ratio of the size of an LMRD code divided by the size of the Singleton bound converges for monotonically decreasing to .
Proof.
Setting and the ratio is given by , so that Lemma 5 gives the proposed limit. The sequence is monotonically decreasing, since we have and
[TABLE]
∎
Proposition 8**.**
For the ratio of the size of an LMRD code divided by the size of the Anticode bound converges for monotonically decreasing to .
Proof.
The LMRD code has cardinality and the Anticode bound is . From Lemma 5 we conclude
[TABLE]
so that the limit follows. The subsequent inequalities follow from , the monotonicity of , and .
It remains to show the monotonicity of the sequence
[TABLE]
Using the abbreviation we define
[TABLE]
and observe that is strictly monotonically increasing in , so that . Using routine manipulations of -binomial coefficients we compute
[TABLE]
∎
In other words the ratio between the best known lower bound and the best known upper bound for constant dimension codes is strictly greater than for all parameters and the most challenging parameters are given by , , and .
Replacing the Anticode bound by the Johnson bound of Theorem 2 does not change the limit behavior of Proposition 8 when tends to infinity. As stated above, we obtain the Anticode bound if we remove the floors in Corollary 2 and replace by . First we consider the latter weakening. Applying the lower bound of Theorem 12 for , where with , we consider
[TABLE]
If , then the subtrahend on the right hand side is at most . Otherwise we have , so that . Since the subtrahend on the right hand side is at most . Thus, the ratio between the lower and the upper bound for partial spreads tends to if tends to infinity. Since
[TABLE]
the ratio between Corollary 2 and the Anticode bound tends to as tends to infinity.
Next, we consider the ratio between the lower bound from the first construction of Proposition 6 and the Anticode bound when tends to infinity.
Proposition 9**.**
For integers satisfying the conditions of Proposition 6, and , we have
[TABLE]
Proof.
For and we have , so that
[TABLE]
where we abbreviate and . Thus,
[TABLE]
Plugging in the definition of the -binomial coefficients gives
[TABLE]
so that
[TABLE]
Dividing both derived limits gives the proposed result.
∎
For Example 1 with and , we obtain a ratio of
[TABLE]
for with using , i.e., the Anticode bound is almost met by the underlying improved linkage construction.
6 Codes better than the MRD bound
For constant dimension codes that contain a lifted MRD code, Theorem 10 gives an upper bound which is tighter than the Johnson bound of Theorem 7. In [5] two infinite series of constructions have been given where the code sizes exceed the MRD bound of Theorem 10 for , , and . Given the data available from [23] we mention that, besides , , the only other case where the MRD bound was superseded is , see [10]. Next, we show that for and the MRD bound can be superseded for all field sizes if is large enough. For the limit of the achievable ratio we obtain:
Proposition 10**.**
For we have .
Proof.
For , [25, Theorem 4] gives
[TABLE]
With this, we conclude
[TABLE]
from Corollary 4 choosing . Applying Proposition 6 with gives
[TABLE]
for , so that for all .
From Lemma 5 we conclude
[TABLE]
Since
[TABLE]
the statement follows for .
∎
For the estimations of the proof of Proposition 10 are too crude in order to obtain a factor larger than one. However, for the binary case better codes with moderate dimensions of the ambient space have been found by computer searches – with the prescription of automorphisms as the main ingredient in order to reduce the computational complexity, see e.g. [31].
Proposition 11**.**
For we have .
Proof.
Applying Proposition 6 with and using gives for all and , so that
[TABLE]
Using [23], [10], [10], and [9] we apply Corollary 4 with to obtain lower bounds for with . For these values of the minimum of the right hand side of Inequality (9) is attained at with value .
∎
Note that the application of Proposition 6 was used in a rather crude estimation in the proof of Proposition 11. Actually, we do not use the codewords generated by the codewords of cdc in Theorem 18, so that we might have applied [38, Theorem 37] directly for this part of the proof – similarly for Proposition 10, which then allows to consider just one instead of starters. In the latter part of the proof of Proposition 11 the use of Corollary 4 is essential in order to obtain large codes for medium sized dimensions of the ambient space from and relatively good lower bounds for small dimensions. This is a relative typical behavior of Corollary 4 and Proposition 6, i.e., the first few applications yield a significant improvement which quickly bottoms out – in a certain sense. As column bklb of Table 3 suggests, we may slightly improve upon the value stated in Proposition 11 by some fine-tuning effecting the omitted less significant digits.
In Tables 1, 2, and 3 we compare the sizes of different constructions with the LMRD and the best known upper bound. Here bklb and bkub stand for best known lower and upper bound respectively. The values of Theorem 10 are given in column mrdb. Applying Theorem 17 and Theorem 18 to the best known codes give the columns lold and lnew, respectively. The results obtained in [5] are stated in column ea. The achieved ratio between the mentioned constructions and the MRD bound can be found in Table 3. Since differences partially are beyond the given accuracy, we give absolute numbers in Table 1. Note that the values in column bklb of Table 3 show that Proposition 11 is also valid for , while we have a smaller ratio for . The relative advantage over lifted MRD codes is displayed in Table 2.
To conclude this section, we remark that an application of Corollary 4 with using a lifted MRD in the cdc cannot generate a code that exceeds the MRD bound of Theorem 10.
Lemma 6**.**
Using the notation of Theorem 18, let , a linear MRD code, , and contains a lifted MRD code (in ). Then, the codes constructed in Theorem 18 contain a lifted MRD code (in ).
Proof.
Let be the lifted MRD code in . Since is a MRD code, we have . The first set of the construction contains
[TABLE]
in which forms a rank metric code of size , hence it is a maximum rank metric code.
∎
7 Conclusion
In this paper we have considered the maximal sizes of constant dimension codes. With respect to constructive lower bounds we have improved the so-called linkage construction, which then yields the best known codes for many parameters, see Footnote 2. With respect to upper bounds there is a rather clear picture. The explicit Corollary 2, which refers back to bounds for partial spreads, is the best known parametric bound in the case of , while Theorem 8 or the linear programming method may possibly yield improvements. Since Theorem 8 implies the Johnson bound and so Corollary 2, it would be worthwhile to study whether it can be strictly sharper than Theorem 7 for at all. Compared to Corollary 2, the only two known improvements are given for the specific parameters from Theorem 9 and Proposition 5. In the case of partial spreads we have reported the current state-of-the-art mentioning that further improvements are far from being unlikely.
In general we have shown that the ratio between the best-known lower and upper bound is strictly larger than for all parameters. The bottleneck is formed by the parameters , , and , where no known method can properly improve that factor, see Footnote 2 for the linkage construction. For , and general field sizes we have applied the improved linkage construction in order to show that is by a factor, depending on , larger than the MRD bound for sufficiently large dimensions .
Acknowledgement
The authors would like to thank Harout Aydinian for providing an enlarged proof of Theorem 8, Natalia Silberstein for explaining the restriction in [38, Corollary 39], Heide Gluesing-Luerssen for clarifying the independent origin of the linkage construction, and Alfred Wassermann for discussions about the asymptotic results of Frankl and Rödl.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] E. Agrell, A. Vardy, and K. Zeger. Upper bounds for constant-weight codes. IEEE Transactions on Information Theory , 46(7):2373–2395, 2000.
- 2[2] R. Ahlswede and H. Aydinian. On error control codes for random network coding. In Network Coding, Theory, and Applications, 2009. Net Cod’09. Workshop on , pages 68–73. IEEE, 2009.
- 3[3] R. Ahlswede, H. K. Aydinian, and L. H. Khachatrian. On perfect codes and related concepts. Designs, Codes and Cryptography , 22(3):221–237, 2001.
- 4[4] R. Ahlswede, N. Cai, S.-Y. R. Li, and R. W. Yeung. Network information flow. IEEE Transactions on Information Theory , 46(4):1204–1216, 2000.
- 5[5] J. Ai, T. Honold, and H. Liu. The expurgation-augmentation method for constructing good plane subspace codes. ar Xiv preprint 1601.01502 , 2016.
- 6[6] C. Bachoc, A. Passuello, and F. Vallentin. Bounds for projective codes from semidefinite programming. Advances in Mathematics of Communications , 7(2):127–145, 2013.
- 7[7] A. Beutelspacher. Partial spreads in finite projective spaces and partial designs. Mathematische Zeitschrift , 145(3):211–229, 1975.
- 8[8] S. R. Blackburn and T. Etzion. The asymptotic behavior of grassmannian codes. IEEE Transactions on Information Theory , 58(10):6605–6609, 2012.
