On the classification of unitals on 28 points of low rank

The classification of unitals with parameters 2-(28, 4, 1) according to the 2-rank of their incidence matrices was initiated by McGuire, Tonchev and Ward, who proved that the 2-rank of any unital on 28 points is greater than or equal to 19, and up to isomorphism, there is a unique unital with 2-rank equal to 19. Jaffe and Tonchev investigated the next two 2-ranks, 20 and 21, and showed that there are no unitals on 28 points with 2-rank equal to 20, and there are exactly 4 isomorphism classes of unitals of rank 21. The subject of this paper is the classification of unitals having 2-rank 22, 23 and 24.

Let V be a set of v points and B be a collection of k-subsets of V called blocks.
The number r is called the replication number of D . An isomorphism between two designs D ′ , D ′′ is a bijection between their point sets that maps every block of D ′ to a block of D ′′ . An automorphism of a design D is any isomorphism of D to itself.
On the classification of unitals on 28 points of low rank introduced a new technique called paramodification to construct new designs from known ones and found 173 new 2- (28,4,1) designs. For the rest of the paper, we will focus on the case q = 3 and we will use the term unital synonymously for a 2- (28,4,1) design. Since r = 9, = 1 and r − = 8 = 2 3 , the p-rank of a 2-(28, 4, 1) design can be smaller than 27 only if p = 2 , and we will write rank instead of 2-rank.
Another way to classify 2- (28,4,1) designs is according to their rank. In 1981, Brouwer [7] computed the ranks and the binary codes of the classical Hermitian unital that has rank 21, the unital found by Lüneburg which Brouwer called the Ree unital and has 2-rank 19, and 136 other unitals found by computer, with ranks ranging from 21 to 27.
In [18] McGuire, Tonchev, and Ward proved that the rank of a unital on 28 points is greater than or equal to 19, with equality achieved by the Ree unital only.
Although the unique 2-(28, 4, 1) design of rank 19 is not embeddable in any of the projective planes of order 9, this design has a nice geometric interpretation, namely, it is isomorphic to a design associated with a maximal arc of degree 4 in the projective plane of order 8. Recall that a maximal arc of degree d in a projective plane Π of order q = sd is a set A of (s(d − 1) + 1)d points, such that every line is either disjoint from A or intersects A in exactly d points. If d > 1 , the non-empty intersections of lines of Π with A are the blocks of a 2-((s(d − 1) + 1)d, d, 1) design. The set of lines of P which have no points in common with A determine a maximal arc of degree s in the dual plane Π ⟂ , called the dual arc A ⟂ . A hyperoval is a maximal arc of degree 2. In particular, if A is a hyperoval in the projective plane of order 8, then the dual arc A ⟂ is a maximal arc of degree 4 that determines a 2-(28, 4, 1) design having 2-rank 19. This design is a member of an infinite class of 2-(2 2t−1 − 2 t−1 , 2 t−1 , 1) designs associated with dual arcs of hyperovals in the Desarguesian plane of order 2 t , PG(2, 2 t ) . It was proved by Mackenzie [17] that the 2-rank of any such design is greater than or equal to 3 t − 2 t , and Carpenter [8] proved that the 2-rank is equal to 3 t − 2 t if the related hyperoval is regular. It is conjectured that the 2-rank of every maximal arc of degree 2 t−1 in PG(2, 2 t ) is equal to 3 t − 2 t [26].
In [11] Jaffe and Tonchev proved that there are no unitals of rank 20, and classified the unitals of rank 21. There are exactly four nonisomorphic unitals of rank 21. One of these unitals is the Hermitian unital H(3) with a full automorphism group of order 12,096, while the remaining three unitals have automorphism groups of order 192, 24, and 6, respectively. The four binary codes of dimension 21 of the four unitals of rank 21 are pairwise inequivalent.
In this paper, we report the classification of unitals of ranks 22, 23, and 24. There are 12 isomorphism classes of unitals of 2-rank 22, 78 isomorphism classes of unitals of 2-rank 23, four of which have a trivial automorphism group, and 298 isomorphism classes of unitals of 2-rank 24, 38 of which have a trivial automorphism group.

Classification method
Our classification method extends the work of Jaffe and Tonchev [11], but we used different software tools. An abstract description of the classification is given in Algorithm 1, where C denotes the binary [28, rk ] code of a 2-(28, 4, 1) design of rank rk , and d is the minimum weight of the dual code C ⟂ . Algorithm 1 Input: A fixed rank rk ≥ 19.
1. Classify the binary linear [28, 28 − rk , d ≥ 10] codes C ⟂ that meet all known restrictions. 2. For each code C ⟂ from 1) do: a. Compute the set of all codewords of C of weight 4. b. Enumerate all unitals having incidence matrices with rows being codewords of C of weight 4. c. Classify the unitals obtained from C up to isomorphism.

Classify the unitals obtained from all codes obtained in 1).
Remark We applied Algorithm 1 successively for rk = 22, 23, 24 and compared the outcome of each rank with the previous results on smaller ranks. The check if the determined numbers of unitals of rank < rk are confirmed by known results is a necessary condition for the computer programs to be correct.

Code restrictions
A binary code is even if all codewords have even weight, and doubly-even if all codewords are of weight divisible by 4. An [n, k, d] code C is a linear code of length n, dimension k, and minimum weight d. The dual code C ⟂ of an [n, k, d] code is the linear code of length n and dimension n − k consisting of all vectors that are orthogonal to every codeword of C. A code C is self-orthogonal if C ⊆ C ⟂ . A linear binary code is self-complementary if it contains the all-one vector. It was pointed out by Brouwer [7] that every codeword of weight 12 in the dual code of the binary code of any unital U on 28 points is the sum of the incidence vectors of three pairwise disjoint blocks of U.
The binary codes of 2-(28, 4, 1) designs have been studied in [7,11,18]. Here, we summarize the properties of such codes that we used in the computation. 1. The binary code C of a unital of 2-rank rk is an even self-complementary [28, rk ] code. 2. Its dual code C ⟂ is an even self-complementary [28, 28 − rk , d ≥ 10] code.

3
On the classification of unitals on 28 points of low rank [28, t, 12] code for some t ≤ 28 − rk.

The codewords of weight divisible by 4 form a doubly-even self-complementary
In order to determine the possible values of t the following result can be used. Corollary 1 If C is the binary code of a 2-(28, 4, 1) design then D is a doubly-even subcode of C ⟂ of dimension smaller than or equal to t. This implies that 2 t is one of the numbers in (1).

Unitals of rank 22
The code C of a 2-(28, 4, 1) design of rank 22 is a linear binary code of dimension 22 and its dual code C ⟂ has dimension 6. It follows from Corollary 1 that the dimension of the self-complementary doubly-even subcode of C ⟂ is greater than or equal to 4. We classified the self-complementary doubly even [28,4,12] codes with the algorithm described in [6], using the software tools solvediophant by the second author [28,29] and codecan by Feulner [9]. There are four equivalence classes of such codes.
These four codes can be extended to 128 self-complementary even [28, 5, d ≥ 10] codes. This was done with an algorithm similar to the one described in [5,15] but with a less sophisticated approach to suppress the generation of isomorphic copies. Again, solvediophant and codecan were used. The 128 self-complementary even [28, 5, d ≥ 10] codes were extended to 15 007 self-complementary even [28, 6, d ≥ 10] codes.
For each of these codes, the set of all codewords of weight 4 was constructed and tested if it contains unitals. Clearly, the incidence matrix of any unital consists of 63 codewords of weight 4 that form a constant weight code with minimum distance 6. Using the exact cover solver dlx by Knuth [13], it turns out that there are 661 504 unitals that arise from 124 codes, and these unitals fall into 17 isomorphism classes. For the latter computation, the GAP design package by Soicher [23] was used. The 17 isomorphism classes consist of one unital of rank 19, four unitals of rank 21, and 12 unitals of rank 22, confirming the results of [11,18]. The 12 unitals of rank 22 all have nontrivial automorphisms and therefore appear already in [14].

3
The 12 unitals of rank 22 span ten inequivalent codes of dimension 22: nine codes contain a single isomorphism class of unitals of 2-rank 22, and one code contains three isomorphism classes of unitals of 2-rank 22. The distribution of unitals among the ten codes according to their rank and automorphism group order is given in Table 1, where the second column contains the order of the automorphism group and 2-rank of the unitals contained in the code. The entry (48, 22) 2 in the last row of Table 1 indicates that the code contains two isomorphism classes of unitals of rank 22 having full automorphism group of order 48.
As a check of correctness, we computed the codes of all 17 unitals and searched for all unitals contained in these codes. The GAP design package reduced the list of unitals to the same set of 17 unitals.

Unitals of rank 23
The code C of a 2-(28, 4, 1) design of rank 23 is a linear binary code of dimension 23 and its dual code C ⟂ has dimension 5.
It follows from Corollary 1 that the self-complementary doubly-even subcode of C ⟂ has dimension 3.
In contrast with the rank 22 case, a straightforward application of dlx for the search of unitals in the sets of codewords of weight 4 was not feasible for the codes of dimension 23, because the resulting exact cover problems were too time consuming. Instead, we first classified all derived 1-(27, 3, 1) designs with respect to point 28 under the stabilizer of the automorphism group of the code. In other words, all feasible, non-isomorphic combinations of 9 blocks through point 28 were enumerated, and each of these 9-block-configurations was extended to unitals using dlx.
Among the 880 [28,23] codes, 126 codes contain together 16 911 unitals which fall into 95 isomorphism classes. The set of 95 nonisomorphic unitals consists of one unital of rank 19, four unitals of rank 21, 12 unitals of rank 22, and 78 unitals of rank 23. Again, the results of [11,18] and of the rank 22 case are confirmed. Among the 78 unitals of rank 23, there are exactly four unitals with trivial automorphism group. The 78 unitals of rank 23 span 30 inequivalent [28,23] codes, see Table 2, where an entry (a, b) denotes a unital of rank b with a full automorphism group of order a, and (a, b) s denotes s nonisomorphic unitals of rank b with a full automorphism group of order a.

Unitals of rank 24
The code C of a 2-(28, 4, 1) design of rank 24 is a linear binary code of dimension 24 and its dual code C ⟂ has dimension 4. It follows from Corollary 1 that the self-complementary doubly-even subcode of C ⟂ has dimension 2. It is easy to see that, up to equivalence, there is only one such code, and this code can be extended to nine self-complementary even [28, 3, d ≥ 10] codes and 82 self-complementary even [28, 4, d ≥ 10] codes.
The straightforward enumeration of unitals in the sets of codewords of weight 4 as in the rank 22 case was not possible here, and even the approach of the rank 23 case was too time consuming. In order to reduce the enumeration time even further, for each code we determined a transversal of all block orbits containing point 28 under the stabilizer of 28 of the automorphism group of the code, i.e. all possibilities were determined to fix one block. For each of these blocks, all 1-(27, 3, 1) designs that contain that block were enumerated, and all possibilities to fix nine blocks were classified. Finally, each of these sets of nine blocks is extended to unitals.
Again, the results of [11,18] and the rank 22 and rank 23 cases were confirmed, leaving 298 unitals of rank 24. It turns out that 38 unitals of rank 24 have a trivial automorphism group. The 298 unitals of rank 24 span 31 inequivalent codes, see Table 3.
In all cases, computations involving groups were done with GAP [25]. The complete list of unitals up to rank 24 can be downloaded from [27].
The distribution of the unitals according to their rank and order of the automorphism group is given Table 4. The numbers of unitals with non-trivial automorphism groups coincide with the results of Krčadinac [14], and that gives further confidence in the correctness of the computations.
The enumeration of unitals of rank 25 and 26, as well as a full classification seems to be still out of reach with the approach described in this paper.
Funding Open Access funding enabled and organized by Projekt DEAL.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/ licenses/by/4.0/.