# Goldberg–Coxeter construction

The Goldberg–Coxeter construction or Goldberg–Coxeter operation (GC construction or GC operation) is a graph operation defined on regular polyhedral graphs with degree 3 or 4. It also applies to the dual graph of these graphs, i.e. graphs with triangular or quadrilateral "faces". The GC construction can be thought of as subdividing the faces of a polyhedron with a lattice of triangular, square, or hexagonal polygons, possibly skewed with regards to the original face: it is an extension of concepts introduced by the Goldberg polyhedra and geodesic polyhedra. The GC construction is primarily studied in organic chemistry for its application to fullerenes, but it has been applied to nanoparticles, computer-aided design, basket weaving, and the general study of graph theory and polyhedra.
The Goldberg–Coxeter construction may be denoted as $GC_{k,\ell }(G_{0})$, where $G_{0}$ is the graph being operated on, $k$ and $l$ are integers, $k>0$, and $\ell \geq 0$.

## History

Michael Goldberg introduced the Goldberg polyhedron in 1937. Buckminster Fuller coined the term "geodesic dome" in the 1940s, although he largely kept the mathematics behind the domes a trade secret. Geodesic domes are the geometric dual of (a section of) a Goldberg polyhedron: a full geodesic dome can be thought of as a geodesic polyhedron, dual to the Goldberg polyhedron. In 1962, Donald Caspar and Aaron Klug published an article on the geometry of viral capsids that applied and expanded upon concepts from Goldberg and Fuller. H.S.M. Coxeter published an article in 1971 covering much of the same information. Caspar and Klug were the first to publish the most general correct construction of a geodesic polyhedron, making the name "Goldberg–Coxeter construction" an instance of Stigler's law of eponymy.
The discovery of Buckminsterfullerene in 1985 motivated research into other molecules with the structure of a Goldberg polyhedron. The terms "Goldberg–Coxeter fullerene" and "Goldberg–Coxeter construction" were introduced by Michel Deza in 2000. This is also the first time the degree 4 case was considered.

## Construction

This section largely follows Deza et al.'s two articles.

### Master polygons

Lattices :
| n-Regular | 3 | 4 |
| Domain | Eisenstein $a+bu$ | Gaussian $a+bi$ |
| Adjoined unit | ${\begin{aligned}u&={\frac {1}{2}}\left(1+i{\sqrt {3}}\right)\\&=e^{{\frac {1}{3}}\pi i}=\omega +1\end{aligned}}$ | $i={\sqrt {-1}}$ |
| Norm $t(a,b)$ | $a^{2}+ab+b^{2}$ | $a^{2}+b^{2}$. |
| Master polygon | $(0,x,xu)$ | $(0,x,x(1+i),xi)$ |
Regular lattices over the complex plane can be used to create "master polygons". In geodesic dome terminology, this is the "breakdown structure" or "principal polyhedral triangle" (PPT). The 4-regular case uses the square lattice over the Gaussian integers, and the 3-regular case uses triangular lattice over the Eisenstein integers. For convenience, an alternate parameterization of the Eisenstein integers is used, based on the sixth root of unity instead of the third. The usual definition of Eisenstein integers uses the element ω=1/2 (−1+i√3)=e$^{2/3πi}$=u−1{\textstyle \omega ={\frac {1}{2}}\left(-1+i{\sqrt {3}}\right)=e^{{\frac {2}{3}}\pi i}=u-1}. A norm, $t(k,\ell )$, is defined as the square of the absolute value of the complex number. For 3-regular graphs this norm is the T-number or triangulation number used in virology.
The master polygon is an equilateral triangle or square laid over the lattice. The table to the right gives formulas for the vertices of the master polygons in the complex plane, and the gallery below shows the (3,2) master triangle and square. So that the polygon can be described by a single complex number, one vertex is fixed at 0. There are multiple numbers that can describe the same polygon: these are associates of each other: if $x$ and $y$ are associates, then $x=u^{n}y$ in the Eisensteins or $x=i^{n}y$ in the Gaussians for some integer $n$. The set of elements that are associates of each other is an equivalence class, and the element of each equivalence class that has $a>0$ and $b\geq 0$ is the normal form.
* (3,2) master triangle over triangular grid
* (3,2) master square over square grid
Master polygons, and the operator $GC_{k,\ell }\left(G_{0}\right)$, can be classified as follows:
* Class I: $\ell =0.$
* Class II: $k=\ell .$
* Class III: all other. Class III operators exist in chiral pairs: $GC_{k,\ell }\left(G_{0}\right)$ is the chiral pair of $GC_{\ell ,k}(G_{0})$.
Below are tables of master triangles and squares. Class I corresponds to the first column, and Class II corresponds to the diagonal with a slightly darker background.

#### Master polygons for triangles

Master triangles through (8,8) :
* 1
* 2
* 3
* 4
* 5
* 6
* 7
* 8

#### Master polygons for squares

Master squares through (8,8) :
* 1
* 2
* 3
* 4
* 5
* 6
* 7
* 8
Composition of Goldberg–Coxeter operations corresponds to multiplication of complex numbers. If and only if $GC_{a,b}\left(GC_{c,d}\left(G_{0}\right)\right)=GC_{e,f}\left(G_{0}\right)$ (i.e. the series of operations on the left produces a graph isomorphic to the one on the right), then for a 3-regular graph $(a+bu)(c+du)$ is in the equivalence class of $(e+fu)$, and for a 4-regular graph $(a+bi)(c+di)$ is in the equivalence class of $(e+fi)$. There are some useful consequences of this:
* The application of repeated Goldberg–Coxeter operations is commutative and associative.
* Complex conjugation of the element $a+bi$ or $a+bu$ corresponds to reflection of the constructed graph.
* Since the Gaussian integers and Euclidean integers are both Euclidean domains, elements of those domains can be uniquely factored into prime elements. Therefore, it also makes sense to decompose a Goldberg–Coxeter operator into a sequence of "prime" Goldberg–Coxeter operators, and this sequence is unique (up to rearrangement).

### Performing the GC construction

The steps of performing the GC construction $GC_{k,\ell }(G_{0})$are as follows:
* Determine the master polygon, based on $k$, $\ell$, and $G_{0}$
* If operating on a 3- or 4-regular graph (instead of a graph with triangular/quadrilateral faces), take its dual graph. This dual graph will have triangular or quadrilateral faces.
* Replace the faces of the triangulated/quadrangulated graph with the master polygon. Be aware that planar graphs have an "external" face that must be replaced as well. In the below example, this is done by attaching it to one side of the graph and connecting other sides as necessary. This temporarily introduces overlapping edges into the graph, but the resulting graph is planar. The vertices can be rearranged so that there are no overlapping edges.
* If the original graph was a 3- or 4-regular graph, take the dual of the result of step 3. Otherwise, the result of step 3 is the GC construction.
Below is an example, where $GC_{1,1}(G_{0})$ is constructed on the skeleton of a cube. In the last two graphs, blue lines are edges of $G_{0}$, while black lines are edges of $GC_{1,1}(G_{0})$. (Dotted lines are normal graph edges, just drawn differently to make overlapping graph edges more visible.) Red vertices are in $G_{0}$ and remain in $GC_{1,1}(G_{0})$, while blue vertices are newly created by the construction and are only in $GC_{1,1}(G_{0})$.
* Master square (1,1)
* Starting polyhedron (Cube)
* $G_{0}$, the skeleton of the cube
* Intermediate step of constructing $GC_{1,1}(G_{0})$.
* The result $GC_{1,1}(G_{0})$, after rearrangement
* Embedding of the result (rhombic dodecahedron)

## Extensions

The Goldberg–Coxeter construction can be easily extended to some non-planar graphs, such as toroidal graphs. Class III operators, because of their chirality, require a graph that can be embedded on an orientable surface, but class I and II operators can be used on non-orientable graphs.