The Golden Code for 2×22 \times 2 MIMO

From Existence to Construction

Chapter 12 closed with the Zheng-Tse theorem: the DMT curve of an nt×nrn_t \times n_r i.i.d. Rayleigh channel is d(r)=(ntr)(nrr)d^*(r) = (n_t - r)(n_r - r), achievable in principle by a random Gaussian codebook of block length Lnt+nr1L \ge n_t + n_r - 1. The point is that this is a non-constructive existence proof — a capacity-achieving argument that tells us what is achievable but not how to build the code. The practical question is: can we write down an explicit codebook — a structured, algebraically clean family of nt×Ln_t \times L complex matrices — that achieves the full Zheng-Tse curve, decodes in polynomial-in-MM complexity, and does not require the random-coding infinities?

For nt=2n_t = 2, Belfiore, Rekaya, and Viterbo answered this in 2005 with a breathtakingly simple construction: the Golden code. Four complex information symbols a,b,c,da, b, c, d drawn from a QAM constellation are packed into a 2×22 \times 2 matrix built from the golden ratio θ=(1+5)/2\theta = (1 + \sqrt{5})/2 — the same irrational number that governs Fibonacci growth and classical æsthetic proportions, now pressed into service as the algebraic basis for a DMT-optimal space-time code. The code achieves rate 22 complex symbols per channel use (full multiplexing on a 2×22 \times 2 channel), full diversity 44 at r=0r = 0, and — most importantly — a non- vanishing determinant δmin=1/5\delta_{\min} = 1/5 that does not shrink as we grow the input constellation. This last property is the algebraic fingerprint of DMT optimality; it is what separates the Golden code from naive lattice codes whose determinants collapse at large QAM sizes.

This section defines the Golden code, proves the non-vanishing determinant via the algebraic-norm argument, and unpacks what "coding gain" means when the constellation scales. §2 generalises to arbitrary ntn_t via cyclic division algebras; §3 presents the Perfect-codes family; §4 proves approximate universality; §5 makes the NVD property precise.

, ,

Historical Note: The 2005 Golden Code

2005

The Golden code was introduced in a 4-page IEEE Transactions on Information Theory correspondence by Jean-Claude Belfiore, Ghaya Rekaya, and Emanuele Viterbo (April 2005). At the time, Zheng and Tse's DMT theorem (2003) was only two years old and the question of whether any explicit full-rate 2×22 \times 2 space-time code could achieve the full DMT curve was open. Prior full-rate designs (threaded algebraic space-time codes, the Damen-Tewfik- Belfiore code) had determinants that decayed as 1/M1/\sqrt{M} with constellation size — acceptable at small MM, fatal at large MM.

Belfiore-Rekaya-Viterbo's key move was to use the number field Q(j,5)=Q(j,θ)\mathbb{Q}(j, \sqrt{5}) = \mathbb{Q}(j, \theta) — a degree-2 extension of the Gaussian rationals Q(j)\mathbb{Q}(j) — and exploit the fact that the algebraic norm on this field satisfies N(a+θb)=(a+θb)(a+θˉb)=a2+abb2Z[j]N(a + \theta b) = (a + \theta b)(a + \bar\theta b) = a^2 + ab - b^2 \in \mathbb{Z}[j], a non-zero Gaussian integer whenever a+θb0a + \theta b \neq 0. The minimum absolute value of a non-zero Gaussian integer is 11 — hence the determinant of any non-zero codeword difference is bounded below by 1/51/|5|, and the 1/51/5 factor comes from the overall normalisation chosen to make the average codeword energy equal to the transmit power.

The paper christened the construction "the Golden code" in a reference to the golden ratio θ\theta that sits at its algebraic heart. It became the most-cited space-time code in wireless communications, adopted as a reference design in the DVB-NGH mobile TV standard (2012) and used as a benchmark for every subsequent DMT-optimal construction. Its fame is deserved: the Golden code is, in Tse and Viswanath's phrase, "as close to a perfect space-time code as one can hope to find."

Definition:

The Golden Code

Let θ=(1+5)/2\theta = (1 + \sqrt{5})/2 (the golden ratio) with Galois conjugate θˉ=1θ=(15)/2\bar\theta = 1 - \theta = (1 - \sqrt{5})/2. Let α=1+jθˉ\alpha = 1 + j\bar\theta and αˉ=1+jθ\bar\alpha = 1 + j\theta be the two units of the Golden-code construction. Let a,b,c,da, b, c, d be four information symbols drawn from a scaled QAM constellation QMZ[j]\mathcal{Q}_M \subset \mathbb{Z}[j] of size MM.

The Golden codeword matrix is the 2×22 \times 2 complex matrix \ntnXGold=15(α(a+θb)α(c+θd)jαˉ(c+θˉd)αˉ(a+θˉb)).\ntn{X}_{\rm Gold} = \frac{1}{\sqrt{5}} \begin{pmatrix} \alpha (a + \theta b) & \alpha (c + \theta d) \\ j \, \bar\alpha (c + \bar\theta d) & \bar\alpha (a + \bar\theta b) \end{pmatrix}.

The factor 1/51/\sqrt{5} is the energy normalisation; α2=1+θˉ2|\alpha|^2 = 1 + \bar\theta^2 and αˉ2=1+θ2|\bar\alpha|^2 = 1 + \theta^2 satisfy α2+αˉ2=5|\alpha|^2 + |\bar\alpha|^2 = 5. The golden code transmits two complex symbols per channel use over nt=2n_t = 2 antennas during L=2L = 2 channel uses — it is full-rate (log2M4/2=2log2M\log_2 M^4 / 2 = 2 \log_2 M bits per channel use).

Intuitively, what happens is this. The top row is built from the number-field element a+θba + \theta b on the first channel use and c+θdc + \theta d on the second. The bottom row is the Galois conjugate of the top row — a+θba + \theta b becomes a+θˉba + \bar\theta b — which provides the diversity. The factor jj in the bottom- left entry is what breaks the "transposed" symmetry and lets the determinant stay away from zero; the units α,αˉ\alpha, \bar\alpha are what make the columns and rows orthogonal on average so the code has uniform average energy across antennas.

,

Example: An Explicit Golden Codeword

Compute the Golden codeword \ntnXGold\ntn{X}_{\rm Gold} for information symbols (a,b,c,d)=(1,j,1,1+j)(a, b, c, d) = (1, j, -1, 1 + j) and verify that det(\ntnXGold)20|\det(\ntn{X}_{\rm Gold})|^2 \ne 0.

Golden Codeword Entries Viewed as a C\mathbb{C}-Valued Constellation

Each information quadruple (a,b,c,d)QM4(a, b, c, d) \in \mathcal{Q}_M^4 produces a 2×22 \times 2 complex Golden codeword with four complex-valued entries. This plot scatters those four entries in the complex plane for the full constellation QM4\mathcal{Q}_M^4 (M=4,16,64M = 4, 16, 64) and reveals the dense lattice structure of the Golden-code output — the algebraic substrate that makes the non-vanishing-determinant property work.

Parameters

Golden Code: Sphere Decoding of a 2×22 \times 2 Codeword

Animation of the Golden-code receiver: an information quadruple (a,b,c,d)(a, b, c, d) is encoded into \ntnXGold\ntn{X}_{\rm Gold}, corrupted by the 2×22 \times 2 i.i.d. Rayleigh channel H\mathbf{H} and AWGN w\mathbf{w}, and decoded by the Viterbo-Boutros sphere decoder. The sphere decoder searches through the lattice cosets {H\ntnX:\ntnXC}\{\mathbf{H} \ntn{X} : \ntn{X} \in \mathcal{C}\} inside a shrinking radius until the nearest lattice point is found — the ML estimate.
The sphere radius (red circle) shrinks after each candidate examined; the ML codeword (green) is the unique lattice point inside the final sphere. Complexity scales as O(Mnt2/2)=O(M2)O(M^{n_t^2/2}) = O(M^2) for nt=2n_t = 2 on average.

Theorem: Non-Vanishing Determinant of the Golden Code

Let \ntnXGold,\ntnX^Gold\ntn{X}_{\rm Gold}, \hat{\ntn{X}}_{\rm Gold} be two distinct Golden codewords built from information quadruples (a,b,c,d)(a, b, c, d) and (a^,b^,c^,d^)(\hat a, \hat b, \hat c, \hat d) in QM4Z[j]4\mathcal{Q}_M^4 \subset \mathbb{Z}[j]^4. Let Δ=\ntnXGold\ntnX^Gold\boldsymbol{\Delta} = \ntn{X}_{\rm Gold} - \hat{\ntn{X}}_{\rm Gold}. Then det(Δ)2    15,|\det(\boldsymbol{\Delta})|^2 \;\ge\; \frac{1}{5}, with the bound achieved; equivalently δmin=1/5\delta_{\min} = 1/5. The bound is independent of the constellation size MM — it does not shrink as MM grows. Hence the Golden code has the non-vanishing-determinant (NVD) property.

The determinant of a Golden codeword factors through the algebraic norm on the number field Q(θ)\mathbb{Q}(\theta). Norms of non-zero algebraic integers are non-zero rational integers; the minimum absolute value of a non-zero Gaussian integer is 11. The 1/51/5 prefactor is the square of the 1/51/\sqrt{5} energy normalisation. Crucially, scaling to larger QAM only changes the magnitudes within Z[j]\mathbb{Z}[j] — it never produces smaller non-zero integers — so the bound survives at every MM.

,

What the Golden Code Achieves

The Golden code is full-rate (rmax=2r_{\max} = 2, the channel's spatial degrees of freedom), full-diversity (dmax=4=ntnrd_{\max} = 4 = n_t n_r for 2×22 \times 2), and has non-vanishing determinant (δmin=1/5\delta_{\min} = 1/5 at every MM). Combined with the rank and determinant criteria of Chapter 11 and the Zheng-Tse DMT upper bound of Chapter 12, these three properties imply that the Golden code achieves the Zheng-Tse DMT curve d(r)=(2r)(2r)d^*(r) = (2 - r)(2 - r) for all r[0,2]r \in [0, 2]. It is the first explicit 2×22 \times 2 construction with this property, and — for its dimension — the benchmark against which all subsequent space-time codes are measured.

,

Key Takeaway

The Golden Code. \ntnXGold=15(α(a+θb)α(c+θd)jαˉ(c+θˉd)αˉ(a+θˉb))\ntn{X}_{\rm Gold} = \tfrac{1}{\sqrt{5}} \begin{pmatrix} \alpha(a + \theta b) & \alpha(c + \theta d) \\ j\bar\alpha(c + \bar\theta d) & \bar\alpha(a + \bar\theta b) \end{pmatrix} over Q(j,θ)\mathbb{Q}(j, \theta) is full-rate, full- diversity, and has δmin=1/5\delta_{\min} = 1/5 independent of MM. It is the first explicit DMT-optimal code for nt=2n_t = 2. The next section generalises this construction to arbitrary dimensions via cyclic division algebras.