Construction of New DelayTolerant
SpaceTime Codes
Abstract
Perfect SpaceTime Codes (STC) are optimal codes in their original construction for Multiple Input Multiple Output (MIMO) systems. Based on Cyclic Division Algebras (CDA), they are fullrate, fulldiversity codes, have NonVanishing Determinants (NVD) and hence achieve DiversityMultiplexing Tradeoff (DMT). In addition, these codes have led to optimal distributed spacetime codes when applied in cooperative networks under the assumption of perfect synchronization between relays. However, they loose their diversity when delays are introduced and thus are not delaytolerant. In this paper, using the cyclic division algebras of perfect codes, we construct new codes that maintain the same properties as perfect codes in the synchronous case. Moreover, these codes preserve their fulldiversity in asynchronous transmission.
I Introduction and Problem Statement
During the past decade, MIMO techniques have experienced a great interest in wireless communication systems. Using multiple antennas at the transmitter and the receiver provides high data rates and exploits the spatial diversity in order to fight channel fadings and hence improve the link reliability. Lately, cooperative diversity has emerged as a new form of spatial diversity via cooperation of multiple users in the wireless system [1]. While preserving the same MIMO benefits, it counteracts the need of incorporating many antennas into a single terminal, especially in cellular systems and adhoc sensor networks, where it can be impractical for a mobile unit to carry multiple antennas due to its size, power and cost limitations.
In cooperative networks, users communicate cooperatively to transmit their information by using distributed antennas belonging to other independent terminals. This way, a virtual MIMO scheme is created, where a transmitter is also acting as a relay terminal, with or without some processing, assisting another transmitter to convey its messages to a destination. The cooperative schemes have been widely investigated by analyzing their performance through different cooperative protocols [1, 2, 3]. These protocols fall essentially into two families: AmplifyandForward (AF) and DecodeandForward (DF). In order to achieve the cooperative diversity, spacetime coding techniques of MIMO systems have also been applied yielding many designs of distributed spacetime codes under the assumption of synchronized relay terminals [2, 3, 4].
However, this a priori condition on synchronization can be quite costly in terms of signaling and even hard to handle in relay networks [5, 6]. Unlike conventional MIMO transmitter, equipped with one antenna array using one local oscillator, distributed antennas are dispersed on different terminals, each one with its local oscillator. Thus, they are not sharing the same timing reference, resulting in an asynchronous cooperative transmission.
On the other hand, in a synchronous transmission, the distributed STCs are constructed basically according to the rank and determinant criteria [7] and hence aim at achieving full diversity. Note that the rows of the codeword matrix represent the different relay terminals (antennas). So, when asynchronicity is evoked, delays are introduced between transmitted symbols from different distributed antennas shifting the matrix rows. This matrix misalignment can cause rank deficiency of the spacetime code, and thus performance degradation.
Therefore, the codes previously designed are no more effective unless they tolerate asynchronicity. Furthermore, an efficient code design should satisfy the fulldiversity order for any delay profile. This intends to guarantee fullrank codewords distance matrix i.e., its rank equal to the number of involved relays, hence leading to the socalled delaytolerant distributed spacetime codes [6].
Ii DelayTolerant Distributed SpaceTime Codes
The first designs of such codes were presented by Li and Xia [6] as fulldiversity binary SpaceTime Trellis Codes (STTC) based on the HammonsEl Gamal stacking construction, its generalization to LuKumar multilevel spacetime codes, and the extension of the latter codes for more diverse AMPSK constellations [8, 9]. Systematic construction including the shortest STTC with minimum constraint length was also proposed in [10], as well as some delaytolerant short binary SpaceTime Block Codes (STBC) [11]. Recently, Damen and Hammons extended the Threaded Algebraic SpaceTime (TAST) codes to asynchronous transmission [12]. The delaytolerant TAST codes are based on three different thread structures where the threads are separated by using different algebraic or transcendental numbers that guarantee a nonzero determinant of the codewords distance matrix. An extension of this TAST framework to minimum delay length codes was considered in [13].
Meanwhile, perfect spacetime block codes that are optimal codes originally constructed for MIMO systems [14, 15, 16, 17], were also investigated for wireless relay networks. In [18, 19], the authors provided optimal coding schemes in the sense of DMT tradeoff [20], based cyclic division algebras for any number of users and for different cooperative strategies. Nevertheless, all these schemes assumed perfect synchronization between users. Then, it was in [21] that Petros and Kumar discussed the delaytolerant version of the optimal perfect code variants for asynchronous transmission. They stated that delaytolerant diagonallyrestricted CDA codes and delaytolerant fullrate CDA codes can be obtained from previous designs by multiplying the codeword matrix by a random unitary matrix. This matrix can be taken specifically from an infinite set of unitary matrices that do not have elements in the code field.
In this paper, we construct delaytolerant distributed codes based on the perfect codes algebras from a different point of view. The new construction is obtained from the tensor product of two number fields, one of them being the field used for the perfect code. The codes are designed in such a way to maintain the same properties of their corresponding perfect codes in the synchronous transmission, namely fullrate, fulldiversity and nonvanishing minimum determinant. In addition, unlike the perfect codes, the new codes preserve the full diversity in the asynchronous transmission.
Iii Background
Before addressing the STC construction, we dedicate this section to briefly review the remarkable properties of the perfect codes as analyzed in [14, 15, 16, 17]. Then, following the framework of [6], we present the cooperative communication model of interest.
Iiia Perfect SpaceTime Block Codes
The concept of Perfect Code was originally proposed in [14, 15] for transmit antennas to describe a square linear dispersion STC . The perfect codes are constructed from cyclic division algebras of degree defined by

and are number fields and the corresponding ring of integers. is called the base field and taken as or since the ST code transmits QAM or HEX information symbols for or , respectively. Thus, the constellations can be seen as finite subsets of the ring of Gaussian integers or Eisenstein integers , respectively.

is a cyclic Galois extension of of degree with or a field extension appropriately chosen in order to get an existing lattice and a division algebra, and an algebraic number.

is the generator of the Galois group , . For an element , the conjugates of are . So, the norm and the trace are defined respectively as
(1) 
the set of nonzero elements of . It is a nonnorm element suitable for the cyclic extension [15].
The cyclic division algebra is then expressed as a right space
(2) (3)
The Perfect Codes satisfy the criteria:

Fullrate: The code transmits symbols drawn from QAM or HEX constellation and thus has a rate of symbols per channel use (spcu).

Fulldiversity: According to the rank criterion [7], the determinant of the codeword distance matrix for any two distinct codewords is nonzero. By code linearity, it can be reduced to
(4) 
Nonvanishing minimum determinant: The minimum determinant of any codeword distance matrix, prior to SNR normalization, is lower bounded by a constant that is independent of the constellation size
(5) 
Cubic shaping: The QAM or HEX constellations are normalized according to the power at the transmitter so that the real vectorized codeword vectors are isomorphic to cubic lattices or . In other words, the rotation matrix encoding the information symbols into each layer is required to be unitary to guarantee the energy efficiency of the codes. The shaping constraint leads thus to two other properties. The first one is the Uniform average transmitted energy per antenna. The second one is the Information losslessness as the unitary linear dispersion matrix allows to preserve the mutual information of the MIMO channel.
Thanks to prominent results on diversitymultiplexing tradeoff [20], the perfect codes also verify two other equivalent properties:

DMT optimality: In [16], Elia et al. proved that the fullrate STCs from cyclic division algebra having NVD property achieve the optimal DMT in Rayleigh fading channel.

Approximate universality: Being CDAbased codes with NVD property, the perfect codes are approximately universal and achieve DMT for arbitrary channel fading distribution.
Satisfying all these criteria, the perfect codes showed to improve the performance in terms of error probability upon the best known codes.
IiiB Cooperative System Model
In the sequel, we consider a cooperative system with a source communicating to a destination via relays in two phases as in Figure 1, and without direct links between the source and the destination. In the first phase, the source broadcasts its message to the potential relays. In the second phase, the relays use the DF protocol to detect the source message then if successfully detected transmit it to the destination. We assume that all the relays are able to achieve error free decoding which could be possible by selecting the sourcerelays links, and consider only the links that are not in outage. Note that it could also be possible that not all the relays may successfully decode the original message, so the number of transmitting relays is usually assumed as a random variable. Since the relays transmission overlap in time and frequency, they can cooperatively implement a distributed spacetime code.
Considering only the second phase of transmission, the system is equivalent to a MIMO scheme where the distributed perfect spacetime code is used by the relays, with transmit antennas one by relay, and receive antennas at the destination. Every time slot , the relays send the column vector of the codeword and the destination receives
(6) 
where is the additive white Gaussian noise with i.i.d complex Gaussian variables with zeromean and variance , , , being the noise variance per real dimension. represents the complex channel matrix modeled as i.i.d Gaussian random variables with zero mean and unit variance . The channel is assumed quasistatic with constant fadings during a transmitted codeword and independent fadings between subsequent codewords. Dealing with square STCs , the codeword matrix contains information symbols carved from twodimensional QAM or HEX finite constellations denoted by .
IiiC Asynchronous Cooperative Diversity
The above expression (6) is valid only when relays are synchronized. In the presence of asynchronicity, the codeword transmission is spanned on more than symbol intervals due to delays. Although the symbol synchronization is not required, we assume that the relays are synchronized at the framecodeword level, which can be provided by means of network feedback signaling from the destination. Therefore, the start and the end of each codeword are aligned for different relays by transmitting zero symbols, and hence there is no interference between codewords transmission. We further assume that the timing errors between different relays are integer multiples of the symbol duration and the fractional timing errors are absorbed in the channel dispersion. In the codeword matrix, these delays are also filled with zeros; they are known at the receiver but not at the transmitting relays [6].
Denoting a delay profile by , a delay corresponds to the relative delay of the received signal from the relay as referenced to the earliest received relay signal. Let denotes the maximum of the relative delays, then from the receiver perspective, the codeword matrix was sent instead of the spacetime code.
IiiD Motivation of the Code Construction
The diversity order of any spacetime code is defined by the minimum rank of the distance codeword matrix over all pairs of distinct codewords [7]. The distributed perfect codes are fullrate fulldiversity for the synchronous transmission between the relays and the destination. Note that in general, a transmission between source, halfduplex relays and destination will result in rate loss. When asynchronicity is introduced, the code is no more fullrate since it is spanned on time instants. Moreover, certain delay profiles can result in linearly dependent rows, thus the code will loose its fulldiversity property. Let us illustrate this by the following example.
Example of Golden Code
We consider the distributed Golden code transmitting information QAM symbols from two synchronized relays with the codeword matrix.
(7) 
The Golden code is designed on a cyclic field extension of degree over the base field . Using the generator matrix of the corresponding complex dimensional lattice, the codeword elements are lattice points obtained by linear combination of pairs of symbols.
Now, let the first relay be delayed by one symbol period with respect to the second , such that the new asynchronous codeword matrix be
(8) 
Suppose we have two distinct codewords and with and the other symbols equal i.e., . The difference between matrix codewords is defined in both synchronous and asynchronous cases as
(9) 
It can be seen that is a fullrank matrix whereas has rank one, so the Golden code is not a delaytolerant code.
In fact, it can be seen from the asynchronous codeword matrix that some symbols are aligned at the same instant due to delays loosing thus diversity. In order to resolve this problem of rank deficiency, our solution consists in transmitting from each antenna (relay) at each transmission time a different combination of all the information symbols. This way, in the presence of delays, we ensure that any combined symbol sent from the relays arrives at the destination in at least different instants, hence guaranteeing the fulldiversity order of the spacetime code.
A new STC will have then the shifted codeword matrix
(10) 
Now, to get these linear combinations of the symbols, we need a higher dimensional lattice compared to the dimensional lattice used for the Golden code. So, we propose to obtain the corresponding lattice generator matrix by the tensor product of two field extensions of , one of them being the field extension of the Golden code.
Following this idea, we aim at constructing, in general, new codes that are based on CDA of the perfect codes such that they maintain the same optimal properties as perfect codes in the synchronous case. But also, these codes preserve their fulldiversity in asynchronous transmission and thus are delaytolerant for arbitrary delay profile.
Iv Construction of DelayTolerant Distributed Codes
based Perfect Codes Algebras
Iva General Construction
The approach consists in constructing a division algebra isomorphic to the tensor product (also called Kronecker product or crossproduct) of two number fields of lower degrees. Other constructions based on the crossedproduct algebras have been investigated in [22, 23] either for prime or coprime degrees of the composite algebras. In these constructions, the spacetime code was built on the cyclic product algebra. However, in the present construction, the higher degree algebra is only used to derive appropriately the spacetime code.
Since we intend to construct a fullrate spacetime code that is based on the CDA of the fullrate perfect code, then the first algebra to be considered is the cyclic division algebra of the perfect code of degree over the base field . For sake of simplicity, we analyze in the sequel the case of Gaussian Field to explain the construction. Indeed, we consider the cyclic field extension of degree over , being an algebraic number. The principal ideal is generated by an element and its integral basis is (or if unitary, it is given by ). The basis of the complex algebraic lattice is obtained by applying the canonical embedding to . Consequently, the generator matrix corresponds to the rotation matrix in
(11) 
where is a normalization factor used to guarantee the matrix unitarity.
Now, we consider another Galois extension over of the same degree such that its discriminant is coprime to the one of i.e., . Let with an algebraic number. The Galois group is generated by as . The principal ideal of the algebra is such that and thus its integral basis is given by . The canonical embedding of gives another complex rotated lattice of that is generated by the unitary matrix with the normalization factor,
(12) 
The tensor product of both field extensions allows to build a rotated lattice in higher dimension corresponding to the complex unitary matrix based on the previous constructions. According to [24],
Proposition
: Let be the compositum of the above Galois extensions, of order over as presented in Figure 2.
Since and have coprime discriminants, the corresponding lattice generator matrix can be obtained as the tensor product of the previous unitary generator matrices.
(13) 
Consequently,
Proposition
: Let the order of the extensions, then the discriminant of is . The minimum product distance of the lattice is derived from the discriminant of as
(14) 
Using the matrix , the spacetime coded components are given by the linear combination where is the information symbol vector carved from a QAM constellation . Then, the spacetime codeword matrix is defined by distributing the components with appropriate constant factors . It can be represented as a Hadamard product
(15) 
The key idea in the code construction is to determine the coefficients that allow one to preserve the same properties of the corresponding perfect codes in synchronous transmission (Section IIIA).

On one side, it can be seen that the new code transmits information symbols and thus is fullrate with spcu for a relaysdestination transmission phase.

On the other side, we need to find the factors that satisfy the rank criterion (4) in order to have fulldiversity codes.

Moreover, the perfect codes have nonvanishing minimum determinants. Then, we are interested in deriving ST codes that have not only nonzero determinants, but also these determinants do not vanish when constellation size increases.

In order to guarantee uniform energy distribution in the codeword, we ask that verify . Choosing further the coefficients yields better determinants as obtained for the nonnorm elements of the perfect codes [15]. This restricts the values of to .

It can also be noticed that the new code satisfies the cubic shaping property since the generator matrix of the dimensional lattice is unitary, and hence the code is information lossless.
In addition, when asynchronicity between relays is involved, the rank criterion should be also verified for the shifted matrix and another criterion will be analyzed that is the nonzero product distance of the codeword matrix in order to prove that the new codes are delaytolerant, and thus keep their fulldiversity in asynchronous transmission.
V New DelayTolerant Codes from dimensional Perfect Codes
Based on the previous approach, we consider the perfect codes proposed in [14, 15] for dimensions to construct the new delaytolerant codes. Then, in the next section, we apply this construction for the perfect codes presented for any number of antennas in [Elia2:2005].
Va Code based on Golden Code
The Golden Code was constructed in [14] using the cyclic division algebra of degree over . is a Galois extension of degree . It is a dimensional vector space of with basis , being the Golden number. Its Galois group is generated by . In order to get a rotated lattice of , the principal ideal generated by was found. Its basis is and its unitary generator matrix is given by
(16) 
with and the respective conjugates of and .
Let the cyclotomic extension of degree over with the primitive root of unity. Its discriminant and it is coprime to the one of since . The Galois group is generated by and the integral basis of is . The corresponding unitary generator matrix is
(17) 
Therefore, is the compositum of Galois extensions of degree each, with coprime discriminants. The unitary matrix is obtained by the tensor product of previous matrices as
(18) 
and the codeword matrix is defined by
(19) 
where are the components of the vector with are QAM symbols. We propose now to determine the coefficients that satisfy the nonvanishing determinant criterion.
VA1 Nonvanishing minimum determinant
The determinant of this codeword matrix equals
(20) 
By developing and , we obtain
(21)  
(22) 
with
(23) 
It is interesting to note that the Golden codeword given by matrix (7) has a determinant of
(24) 
Therefore, by choosing and , the determinant of the new code is equal to the Golden code determinant, and does not vanish when increasing the size of the QAM constellation carved from . Hence, the new code achieves the diversitymultiplexing tradeoff [20, 16].
It can also be noticed that the coefficients can be changed equivalently to the coefficients of the Fourier matrix where is the primitive root of unity. For dimension , we have
(25) 
Furthermore, we have find fixed unitary matrices and such that for all values of with
(26) 
VA2 Delaytolerance
In the distributed setup, each row of the code matrix is transmitted by a different relay (Section IIIB). In practical scenarios, the two relays do not share a common timing reference, and therefore, the arrival of packets is not synchronous. As we assume synchronization at the symbol level, the distributed code can still achieve full diversity if the differences between matrix codewords are full rank even when the different rows are arbitrarily shifted. In what follows, we prove that the new code satisfy this condition.
Consider the shifted codeword matrix of
(27) 
we need to guarantee that it is full rank when i.e., for any from the constellation . This restricts to show that the submatrix
is full rank i.e., its determinant when .
More generally, having delay profiles or , the problem turns to prove that the product distance in the rotated constellation associated with the matrix of is nonzero over , so that any component product is nonzero. This product distance is evaluated as
(28) 
with for .
As a direct consequence from the tensor product construction, Equation (14) gives
Thus, the minimum product distance is nonzero. It can also be verified in by setting . So, is nonzero unless , and consequently the submatrix is full rank since unless .
Therefore, the new code unlike the Golden code keeps its fulldiversity in the case of asynchronous relays. However, we cannot guarantee the nonvanishing determinant property in the asynchronous case because the expression of can be interpreted as a Diophantine approximation of by rational numbers which can be made tighter by using larger constellation size.
VB Code based on Perfect Code
In order to construct the delaytolerant code, we consider the base field and we use HEX symbols. Let , with the root of unity. The perfect code was constructed using the cyclic division algebra of order [15], where the relative extension and the generator of the cyclic extension with . The integral basis is given by and the complex lattice is a rotated version of . It is generated by
(29) 
The relative discriminant of is . Another extension of of degree that has coprime discriminant with is the cyclotomic extension with the primitive root of unity and . Its Galois group is generated by . The integral basis of is and the lattice generator matrix is
(30) 
The compositum of both extensions is of order over . Then, the corresponding dimensional complex lattice is generated by the unitary matrix
(31) 
and the spacetime code is defined by the matrix
(32) 
where are the components of vector , being the information symbol vector carved from HEX constellation.
VB1 Nonvanishing minimum determinant
By proceeding as previously, we need to determine the coefficients that guarantee the nonvanishing minimum determinant. In order to get so that a uniform average energy is transmitted per antenna, and to obtain better values of the determinant, we limit the choice of to .
By developing the code determinant using symbolic computation under Mathematica, we find that it has the same expression as the perfect code determinant by choosing as the Fourier matrix coefficients in
(33) 
Therefore, the infinite code has nonvanishing minimum determinant equal to
(34) 
VB2 Delaytolerance
On the other hand, to prove the delaytolerance of this code, we should guarantee that the corresponding shifted codeword matrices are full rank. Therefore, it suffices to verify that for each asynchronous matrix there exists a square matrix that is full rank i.e., its determinant is nonzero. In fact, if we enumerate all the delay profiles, it can be noticed that the problem of guaranteeing fullrank shifted matrices turns to guarantee that

All component products are nonzero. This condition is always verified since the product distance over as .

All minors of are nonzero that is equivalent to verify that the entries of the cofactor matrix of are nonzero.
In order to prove the second condition, we find two unitary matrices and such that the codeword matrix can be written as for all , with is the perfect code matrix and and are defined by
(35) 
Let define the cofactor matrix of the perfect code by . Since is a finite subset of the cyclic division algebra , is also a subset of taken from the lattice with and is the ring of integers of . Hence, the cofactor matrix can be represented as a codeword matrix. For simplicity, we denote by and , the conjugates of an entry of the codeword matrix. The cofactor codeword matrix is then defined by
(36) 
where each diagonal .
Since , we denote its cofactor matrix. It is given by and satisfies
(37) 
with
(38) 
Developing the cofactor matrix , we get