Almost crystallographic groups are finitely generated nilpotent-by-finite groups satisfying certain properties. There are several group theoretic characterizations known for almost crystallographic groups; for example, a group G is almost crystallographic if and only if G is a finitely generated nilpotent-by-finite group such that the only finite normal subgroup in G is the trivial subgroup. Further, the group G is called an almost Bieberbach group if it is almost crystallographic and torsion free.
Thus an almost crystallographic group G has a finitely generated nilpotent normal subgroup N of finite index. Clearly, N is polycyclic and thus has a polycyclic series. The number of infinite cyclic factors in such a series is an invariant of G: the Hirsch length of G.
Further, for each of the almost crystallographic groups of Hirsch length 3 and 4 there exists a representation as a rational matrix group (in dimension 4 or 5). These representations can be seen as representations into the affine group of dimension 3 or 4. Via these representations, the almost crystallographic groups act (properly discontinuously) on R3 or R4. That is one reason to define the dimension of an almost crystallographic group as its Hirsch length.
The 3-dimensional and a part of the 4-dimensional almost crystallographic groups have been classified by K. Dekimpe, see KD. This classification includes all almost Bieberbach groups in dimension 3 and 4.
This package gives access to the resulting library of groups. The groups are available in two different representations. First, they can be accessed as rational matrix groups. This is also the representation used in KD. Secondly, all 3- and 4-dimensional almost crystallographic groups are polycyclic and they can also be obtained as polycyclically presented groups using the methods from the Polycyclic package of GAP. This second description is probably more useful for computations with the almost crystallographic groups. The package Polycyclic must be installed to use this feature.
Almost crystallographic groups can be seen as natural generalizations of crystallographic groups. A library of crystallographic groups is also available in the GAP package CrystCat. The GAP packages Cryst and CARAT can be used to compute with crystallographic groups.
This documentation is organized as follows. We first recall basic definitions and background on almost crystallographic groups in Section Almost crystallographic groups. Then we outline briefly the main access functions to this library in Sections Rational matrix groups and Polycyclically presented groups. A more extended outline of these access functions is given in Section More about almost crystallographic groups. This library is an electronic version of the classification given in KD and we relate the two version of the library in Section The electronic versus the printed library. Finally, we outline various interesting example computations in Sections Example computations I - Example computations III.
1.1 Almost crystallographic groups
Almost crystallographic groups were first discussed in the theory of actions on Lie groups. We recall the original definition here briefly (see AUS, KD, LEE for more details).
Let L be a connected and simply connected nilpotent Lie group. For
example, the 3-dimensional Heisenberg group, consisting of all upper
unitriangular 3×3--matrices with real entries is of this type.
Then L\rtimes Aut(L) acts affinely (on the left) on L via
|
Let C be a maximal compact subgroup of Aut(L). Then a subgroup G of L \rtimes C is said to be an almost crystallographic group if and only if the action of G on L, induced by the action of L\rtimes Aut(L), is properly discontinuous and the quotient space G \L is compact. One recovers the situation of the ordinary crystallographic groups by taking L=\Bbb Rn, for some n, and C=O(n), the orthogonal group.
More generally, we say that an abstract group is an almost crystallographic group if it can be realized as a genuine almost crystallographic subgroup of some L \rtimes C. In the following theorem we outline some algebraic characterizations of almost crystallographic groups; see Theorem 3.1.3 of KD. Recall that the Fitting subgroup Fitt(G) of a polycyclic-by-finite group G is its unique maximal normal nilpotent subgroup.
Theorem. The following are equivalent for a polycyclic-by-finite group G:
In particular, if G is almost crystallographic, then G / Fitt(G) is finite. This factor is called the holonomy group of G.
The dimension of an almost crystallographic group equals the dimension of the Lie group L above which coincides also with the Hirsch length of the polycyclic-by-finite group. This library therefore contains families of virtually nilpotent groups of Hirsch length 3 and 4.
The following three main functions are available to access the library of almost crystallographic groups as rational matrix groups.
AlmostCrystallographicGroup(
dim,
type,
parameters )
AlmostCrystallographicDim3(
type,
parameters )
AlmostCrystallographicDim4(
type,
parameters )
dim is the dimension of the required group. Thus dim must be either 3 or 4. The inputs type and parameters are used to define the desired group as described in KD. We outline the possible choices for type and parameters here briefly. A more extended description is given later in Section More about almost crystallographic groups or can be obtained from KD.
type specifies the type of the required group. There are 17 types
in dimension 3 and 95 types in dimension 4. The input type can either
be an integer defining the position of the desired type among all types;
that is, in this case type is a number in [1..17]
in dimension 3 or a
number in [1..95]
in dimension 4. Alternatively, type can be a string
defining the desired type. In dimension 3 the possible strings are
"01"
, "02"
, ¼, "17"
. In dimension 4 the possible strings
are listed in the list ACDim4Types
and thus can be accessed from GAP.
ACDim3Param V
ACDim4Param V
parameters is a list of integers. Its length depends on the type of
the chosen group. The lists ACDim3Param
and ACDim4Param
contain
at position i the length of the parameter list for the type number i.
Every list of integers of this length is a valid parameter input.
Alternatively, one can input false
instead of a parameter list. Then
GAP will chose a random parameter list of suitable length.
gap> G := AlmostCrystallographicGroup( 4, 50, [ 1, -4, 1, 2 ] ); <matrix group of size infinity with 5 generators> gap> DimensionOfMatrixGroup( G ); 5 gap> FieldOfMatrixGroup( G ); Rationals gap> GeneratorsOfGroup( G ); [ [ [ 1, 0, -1/2, 0, 0 ], [ 0, 1, 0, 0, 1 ], [ 0, 0, 1, 0, 0 ], [ 0, 0, 0, 1, 0 ], [ 0, 0, 0, 0, 1 ] ], [ [ 1, 1/2, 0, 0, 0 ], [ 0, 1, 0, 0, 0 ], [ 0, 0, 1, 0, 1 ], [ 0, 0, 0, 1, 0 ], [ 0, 0, 0, 0, 1 ] ], [ [ 1, 0, 0, 0, 0 ], [ 0, 1, 0, 0, 0 ], [ 0, 0, 1, 0, 0 ], [ 0, 0, 0, 1, 1 ], [ 0, 0, 0, 0, 1 ] ], [ [ 1, 0, 0, 0, 1 ], [ 0, 1, 0, 0, 0 ], [ 0, 0, 1, 0, 0 ], [ 0, 0, 0, 1, 0 ], [ 0, 0, 0, 0, 1 ] ], [ [ 1, -4, 1, 0, 1/2 ], [ 0, 0, -1, 0, 0 ], [ 0, 1, 0, 0, 0 ], [ 0, 0, 0, 1, 1/4 ], [ 0, 0, 0, 0, 1 ] ] ] gap> G.1; [ [ 1, 0, -1/2, 0, 0 ], [ 0, 1, 0, 0, 1 ], [ 0, 0, 1, 0, 0 ], [ 0, 0, 0, 1, 0 ], [ 0, 0, 0, 0, 1 ] ] gap> ACDim4Types[50]; "076" gap> ACDim4Param[50]; 4
1.3 Polycyclically presented groups
All the almost crystallographic groups considered in this package are polycyclic. Hence they have a polycyclic presentation and this can be used to facilitate efficient computations with the groups. To obtain the polycyclic presentation of an almost crystallographic group we supply the following functions. Note that the package Polycyclic must be installed to use these functions.
AlmostCrystallographicPcpGroup(
dim,
type,
parameters )
AlmostCrystallographicPcpDim3(
type,
parameters )
AlmostCrystallographicPcpDim4(
type,
parameters )
The input is the same as for the corresponding matrix group functions. The output is a pcp group isomorphic to the corresponding matrix group. An explicit isomorphism from an almost crystallographic matrix group to the corresponding pcp group can be obtained by the following function.
IsomorphismPcpGroup(
G )
We can use the polycyclic presentations of almost crystallographic groups to exhibit structure information on these groups. For example, we can determine their Fitting subgroup and ask group-theoretic questions about this nilpotent group. The factor G / Fitt(G) of an almost crystallographic group G is called its holonomy group. We provide access to this factor of a pcp group via the following functions. Let G be an almost crystallographic pcp group.
HolonomyGroup(
G )
NaturalHomomorphismOnHolonomyGroup(
G )
The following example shows applications of these functions.
gap> G := AlmostCrystallographicPcpGroup( 4, 50, [ 1, -4, 1, 2 ] ); Pcp-group with orders [ 4, 0, 0, 0, 0 ] gap> Cgs(G); [ g1, g2, g3, g4, g5 ] gap> F := FittingSubgroup( G ); Pcp-group with orders [ 0, 0, 0, 0 ] gap> Centre(F); Pcp-group with orders [ 0, 0 ] gap> LowerCentralSeries(F); [ Pcp-group with orders [ 0, 0, 0, 0 ], Pcp-group with orders [ 0 ], Pcp-group with orders [ ] ] gap> UpperCentralSeries(F); [ Pcp-group with orders [ 0, 0, 0, 0 ], Pcp-group with orders [ 0, 0 ], Pcp-group with orders [ ] ] gap> MinimalGeneratingSet(F); [ g2, g3, g4 ] gap> H := HolonomyGroup( G ); Pcp-group with orders [ 4 ] gap> hom := NaturalHomomorphismOnHolonomyGroup( G ); [ g1, g2, g3, g4, g5 ] -> [ g1, identity, identity, identity, identity ] gap> U := Subgroup( H, [Pcp(H)[1]^2] ); Pcp-group with orders [ 2 ] gap> PreImage( hom, U ); Pcp-group with orders [ 2, 0, 0, 0, 0 ]
1.4 More about almost crystallographic groups
Each group from this library knows that it is almost crystallographic and it knows its type and defining parameters.
IsAlmostCrystallographic(
G )
AlmostCrystallographicInfo(
G )
We note here that these two properties are defined in GAP for groups from the AC library only. In particular, it is not possible so far, to test if an arbitrary group is almost crystallographic or to determine its type. Methods for this purpose might be added in a later version of this library.
gap> G := AlmostCrystallographicGroup( 4, 70, false ); <matrix group of size infinity with 5 generators> gap> IsAlmostCrystallographic(G); true gap> AlmostCrystallographicInfo(G); rec( dim := 4, type := 70, param := [ 1, -4, 1, 2, -3 ] )
gap> G := AlmostCrystallographicPcpGroup( 4, 70, false ); Pcp-group with orders [ 6, 0, 0, 0, 0 ] gap> IsAlmostCrystallographic(G); true gap> AlmostCrystallographicInfo(G); rec( dim := 4, type := 70, param := [ -3, 2, 5, 1, 0 ] )
Next we consider the types of almost crystallographic groups in more detail. The almost crystallographic groups in dimensions 3 and 4 fall into three families
These families are split up further into subfamilies in KD and to
each subfamily is assigned a type; that is, a string which is used to
identify the subfamily. As mentioned above, for the 3-dimensional almost
crystallographic groups the type is a string representing the numbers from
1 to 17, i.e. types are "01"
, "02"
, ¼, "17"
.
For the 4-dimensional almost
crystallographic groups with a Fitting subgroup of class 2 the type is a
string of 3 or 4 characters. In general, a string of 3 characters representing
the number of the table entry in KD is used. So possible types are
"001"
, "002"
, ¼. The reader is warned however that not all
possible numbers are used, e.g. there are no groups of type "016"
. Also,
the types do not appear in their natural order in KD. Moreover, for
certain numbers there is more than one family of groups listed in KD.
For example, the 3 families of groups corresponding to number 19 on pages
179-180 of KD have types "019"
, "019b"
and "019c"
(the order is
the one given in KD).
For the last category of groups, the 4-dimensional almost crystallographic
groups with a Fitting subgroup of class 3, the type is a string of 2 or 3
characters, where the first character is always the letter "B"
. This "B"
is followed by the number of the table entry as found in KD,
eventually followed by a "b"
or "c"
as in the previous case.
ACDim3Funcs V
ACDim4Funcs V
ACPcpDim3Funcs V
For each type of almost crystallographic group contained in the library
there exists a function taking a parameter list as input and returning
the desired matrix or pcp group. These functions can be accessed
from GAP using the lists ACDim3Funcs
, ACDim4Funcs
, ACPcpDim3Funcs
and ACPcpDim4Funcs
which consist of the corresponding functions.
Although we include these direct access functions here for completeness, we note that the user should in general use the higher-level functions introduced above to obtain almost crystallographic groups from the library. In particular, these low-level access functions return matrix or pcp groups, but the almost crystallographic info flags will not be attached to them.
gap> ACDim3Funcs[15]; function( k1, k2, k3, k4 ) ... end gap> ACDim3Funcs[15](1,1,1,1); <matrix group with 5 generators> gap> ACPcpDim3Funcs[1](1); Pcp-group with orders [ 0, 0, 0 ]
1.5 The electronic versus the printed library
The package ACLib can be considered to be the electronic version of Chapter 7 of KD. In this section we outline the relationship between the library presented in this manual and the printed version in KD. First we consider an example. On page 175 of KD, we find the following groups in the table starting with entry ``13''.
13. Q=P2/c
|
|
|
|
|
The number ``13'' at the beginning of this entry is the type of the almost crystallographic group in this library. This family of groups with type 13 depends on 6 parameters k1, k2, ¼, k6 and these are the parameters list in this library. The rational matrix representation in GAP corresponds exactly to the printed version in KD where it is named l. In the example below, we consider the group with parameters (k1,k2,k3,k4,k5,k6)=(8,0,1,0,1,0).
gap> G:=AlmostCrystallographicDim4("013",[8,0,1,0,1,0]); <matrix group with 6 generators> gap> G.5; [ [ 1, 4, 0, 0, 1/2 ], [ 0, -1, 0, 0, 0 ], [ 0, 0, 1, 0, 0 ], [ 0, 0, 0, -1, 1/2 ], [ 0, 0, 0, 0, 1 ] ] gap> G.6; [ [ 1, 8, 0, 0, 1/2 ], [ 0, -1, 0, 0, 0 ], [ 0, 0, -1, 0, 0 ], [ 0, 0, 0, -1, 0 ], [ 0, 0, 0, 0, 1 ] ]
For a 4-dimensional almost crystallographic group the matrix group is built up such that { a, b, c, d, a, b, g} as described in KD forms the defining generating set of G. For certain types the elements a, b or g may not be present. Similarly, for a 3-dimensional group we have the generating set { a, b, c, a, b} and a and b may be absent.
To obtain a polycyclic generating sequence from the defining generators of the matrix group we have to order the elements in the generating set suitably. For this purpose we take the subsequence of (g, b, a, a, b, c, d) of those generators which are present in the defining generating set of the matrix group. This new ordering of the generators is then used to define a polycyclic presentation of the given almost crystallographic group.
Using the functions available for pcp groups in the package Polycyclic it is now easy to redo some of the calculations of KD. As a first example we check whether the groups indicated as torsion free in KD are also determined as torsion free ones by GAP. In KD these almost Bieberbach groups are listed as ``AB-groups''. So for type ``013'' these are the groups with parameters (k,0,1,0,1,0) where k is an even integer. Let's look at some examples in GAP:
gap> G:=AlmostCrystallographicPcpDim4("013",[8,0,1,0,1,0]); Pcp-group with orders [ 2, 2, 0, 0, 0, 0 ] gap> IsTorsionFree(G); true gap> G:=AlmostCrystallographicPcpDim4("013",[9,0,1,0,1,0]); Pcp-group with orders [ 2, 2, 0, 0, 0, 0 ] gap> IsTorsionFree(G); false
Further, there is also some cohomology information in the tables of KD. In fact, the groups in this library were obtained as extensions E of the form
|
where, in the 4-dimensional case Q = E/ád ñ. The cohomology information for the particular example above shows that the groups determined by a parameter set (k1,k2,k3,k4,k4,k6) are equivalent as extensions to the groups determined by the parameters (k1, k2 mod 2, k3 mod 2, k4 mod 2, k5 mod 2, 0). This is also visible in finding torsion:
gap> G:=AlmostCrystallographicPcpDim4("013",[10,0,2,0,1,0]); Pcp-group with orders [ 2, 2, 0, 0, 0, 0 ] gap> IsTorsionFree(G); false gap> G:=AlmostCrystallographicPcpDim4("013",[10,0,3,0,1,9]); Pcp-group with orders [ 2, 2, 0, 0, 0, 0 ] gap> IsTorsionFree(G); true
The computation of cohomology groups played an important role in the classification of the almost Bieberbach groups in KD. Using GAP, it is now possible to check these computations. As an example we consider the 4-dimensional almost crystallographic groups of type 85 on page 202 of KD. This group E has 6 generators. In the table, one also finds the information
|
for Q=E/ád ñ as above. Moreover, the Q--module Z is in fact the group ád ñ, where the Q-action comes from conjugation inside E. In the case of groups of type 85, Z is a trivial Q-module. The following example demonstrates how to (re)compute this two-cohomology group H2(Q,Z).
gap> G:=AlmostCrystallographicPcpGroup(4, "085", false); Pcp group with orders [ 2, 4, 0, 0, 0, 0 ] gap> GroupGeneratedByd:=Subgroup(G, [G.6] ); Pcp group with orders [ 0 ] gap> Q:=G/GroupGeneratedByd; Pcp group with orders [ 2, 4, 0, 0, 0 ] gap> action:=List( Pcp(Q), x -> [[1]] ); [ [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ] ] gap> C:=CRRecordByMats( Q, action);; gap> TwoCohomologyCR( C ).factor.rels; [ 2, 2, 4, 0 ]
This last line gives us the abelian invariants of the second cohomology group H2(Q,Z). So we should read this line as
|
which indeed coincides with the information in KD.
As another application of the capabilities of the combination of ACLib and Polycyclic we check some computations of DM.
Section 5 of the paper DM is completely devoted to an example of the computation of the P-localization of a virtually nilpotent group, where P is a set of primes. Although it is not our intention to develop the theory of P-localization of groups at this place, let us summarize some of the main results concerning this topic here.
For a set of primes P, we say that n Î P if and only if n is a product of primes in P. A group G is said to be P-local if and only if the map mn:G® G: g ® gn is bijective for all n Î P¢, where P¢ is the set of all primes not in P. The P-localization of a group G, is a P-local group GP together with a morphism a:G ® GP which satisfy the following universal property: For each P-local group L and any morphism j: G ® L, there exists a unique morphism y:GP ® L, such that y°a = j.
This concept of localization is well developed for finite groups and for nilpotent groups. For a finite group G, the P-localization is the largest quotient of G, having no elements with an order belonging to P¢ (the morphism a, mentioned above is the natural projection).
In DM a contribution is made towards the localization of virtually nilpotent groups. The theory developed in the paper is then illustrated in the last section of the paper by means of the computation of the P-localization of an almost crystallographic group. For their example the authors have chosen an almost crystallographic group G of dimension 3 and type 17. For the set of parameters (k1,k2,k3,k4) they have considered all cases of the form (k1,k2,k3,k4)=(2,0,0,k4).
Here we will check their computations in two cases k4=0 and k4=1
using the set of primes P={2}. The holonomy group of these almost
crystallographic groups G is the dihedral group D6 of order
12. Thus there is a short exact sequence of the form
|
First case: The parameters are (k1,k2,k3,k4)=(2,0,0,0)
gap> G := AlmostCrystallographicPcpGroup(3, 17, [2,0,0,0] ); Pcp group with orders [ 2, 6, 0, 0, 0 ] gap> projection := NaturalHomomorphismOnHolonomyGroup( G ); [ g1, g2, g3, g4, g5 ] -> [ g1, g2, identity, identity, identity ] gap> F := HolonomyGroup( G ); Pcp group with orders [ 2, 6 ] gap> IPprimeD6 := Subgroup( F , [F.2^2] ); Pcp group with orders [ 3 ] gap> K := PreImage( projection, IPprimeD6 ); Pcp group with orders [ 3, 0, 0, 0 ] gap> PrintPcpPresentation( K ); pcp presentation on generators [ g2^2, g3, g4, g5 ] g2^2 ^ 3 = identity g3 ^ g2^2 = g3^-1*g4^-1 g3 ^ g2^2^-1 = g4*g5^-2 g4 ^ g2^2 = g3*g5^2 g4 ^ g2^2^-1 = g3^-1*g4^-1*g5^2 g4 ^ g3 = g4*g5^2 g4 ^ g3^-1 = g4*g5^-2 gap> Gamma3K := CommutatorSubgroup( K, CommutatorSubgroup( K, K )); Pcp group with orders [ 0, 0, 0 ] gap> quotient := G/Gamma3K; Pcp group with orders [ 2, 6, 3, 3, 2 ] gap> S := SylowSubgroup( quotient, 3); Pcp group with orders [ 3, 3, 3 ] gap> N := NormalClosure( quotient, S); Pcp group with orders [ 3, 3, 3 ] gap> localization := quotient/N; Pcp group with orders [ 2, 2, 2 ] gap> PrintPcpPresentation( localization ); pcp presentation on generators [ g1, g2, g3 ] g1 ^ 2 = identity g2 ^ 2 = identity g3 ^ 2 = identityThis shows that GP @ Z23.
Second case: The parameters are (k1,k2,k3,k4)=(2,0,0,1)
gap> G := AlmostCrystallographicPcpGroup(3, 17, [2,0,0,1]);; gap> projection := NaturalHomomorphismOnHolonomyGroup( G );; gap> F := HolonomyGroup( G );; gap> IPprimeD6 := Subgroup( F , [F.2^2] );; gap> K := PreImage( projection, IPprimeD6 );; gap> Gamma3K := CommutatorSubgroup( K, CommutatorSubgroup( K, K ));; gap> quotient := G/Gamma3K;; gap> S := SylowSubgroup( quotient, 3);; gap> N := NormalClosure( quotient, S);; gap> localization := quotient/N; Pcp group with orders [ 2, 2, 2 ] gap> PrintPcpPresentation( localization ); pcp presentation on generators [ g1, g2, g3 ] g1 ^ 2 = identity g2 ^ 2 = g3 g3 ^ 2 = identity g2 ^ g1 = g2*g3 g2 ^ g1^-1 = g2*g3
In this case, we see that GP=D4.
The reader can check that these results coincide with those obtained in DM. Note also that we used a somewhat different scheme to compute this localization than the one used in DM. We invite the reader to check the same computations, tracing exactly the steps made in DM.
ACLib manual