I have finished this extension of my essay paper in a rough form. This is still a bit of an outline of things. As I push into this business I find it encompasses considerable depth. This extended version is attached here.
In reference to the three color problem, I attach a file with diagrams lifted from Lisi's paper, since you made mention of it, which illustrates how three colors interact. This is a QCD matter, and even in Lisi's paper he works out how this fits in E_8 as QCD with g_2 = su(3) + 3 + 3-bar = gluons plus quarks. The G_2 group sits inside E_8 and the F_4 is stable or invariant under G_2 action. So for F_4 is the octahedral net of the 24-cell, which is dual to the D_4 subgroup, or 3 sets of 8 rectified tetrahedral cells. This is a source of the F_4 triality. The D_4 group is of course SO(8), which in a non-Euclidean form SO(7,1) decomposes into SO(3,1) and SO(4). Lisi at this point identifies SO(4) with the Pati-Salam weak interactions, but I will not make such an identification at this time. We could just as well identify this with the loop variable gauge-like curvature F^i_{ab} where the upper index is SO(4). The 24-cell is the minimal sphere packing arrangement for a four dimensional spacetime.
What is of course interesting is that F_4 and G_2 are the central groups (centralizers) in E_8 so that the action of one group leaves elements of the other group alone. The G_2 group is also interesting in its own right. For one thing if we have the Jordan exceptional algebra J^3(O) as I outline in my paper, the diagonal elements are scalars, which are the "ones" from each of the octonions. This gives the 3x8 + 3x1 = 27 dimensional exceptional algebra. Now suppose I decide to construct a similar matrix but with the "imaginary" or components of i, j, k, ii, jj, kk, l basis. This will be a 3x8 + 3x7 = 45 dimensional anti-Hermitian matrix. We now do our restrictions on the Hermitian J^3(O) to get 26 dimensions, and we do the same for the anti-Hermitian *J^3(O) with 38 dimensions. The tensor product of this element with the G_2 in 14 dimensions ( *J^3(O)xG_2 ) gives 38 + 14 = 52 dimensions, which is the dimension of the F_4 group. A further tensor product with the J^3(O) extends this to 52 + 26 = 78 dimensions, which is the size of the E_6 group. A similar trick once again with G_2 gives us 7x26 + 52 + 14 = 248. In this way we have that J^3(O) are the real octonions, *J^3(O)xG_2 the complex octonions, *J^3(O)xG_2xJ^3(O), the quaternionic octonions and eventually up to the octonionic octonions.
The exceptional group G2 is the automorphism on O, or equivalently that F_4 x G_2 defines a centralizer on E8. The fibration G_2 - -> S^7 is completed with SO(8), where the three O's satisfy the triality condition in SO(8). So the three octonions in the exceptional J^3(O) exhibit a triality condition inherited from SO(8) by the fibration of G_2. BTW, the triality condition of SO(8) can be seen with the Dynkin diagram with a central o surrounded by three other oriented 2π/3 from each other. So we might consider the role of G_2, which has as it maximal subgroup SU(3) and thus has QCD-like physics, as a source for the three colors. Indeed G_2 was considered as a possible candidate for hadronic gauge theory.
G_2 is the isotropy group of three forms in 7 dimensions. This definition and early structure was advanced by Fredrich Engel and his student Walter Reichel in the early 20 th century. G_2 acts on seven dimensions, the S^7 sphere with 3^2 - 1 generator for a total of 14 dimensions. How this comes about is the following. Consider the complex space C^5 and then for every point a ∈ C5, let the 2-plane π_a exist in the tangent space T_aC^5 t as the zero set of the Pfaffian system
Dx_3 = x_1 dx_2 − x_2 dx_1,
Dx_4 = x_2 dx_3 − x_3 dx_2,
Dx_5 = x_3 dx_1 − x_1 dx_3.
There then 14 vector fields on C^5 whose local flows map the planes π_a to each other satisfy commutator relations of the Lie algebra g_2. Engel derived it from the first by a contact transformation, while Cartan identified g_2 as the symmetries of the solution space of the system of second order partial differential equations5 (f = f (x, y))
f_{xx} = 4/3(f_{yy})^3, f_{xy} = (f_{yy})^2.
The complex Lie group G2 has two non-conjugate 9-dimensional parabolic subgroups P_1 and P_2, and G_2 acts on the two compact homogeneous spaces M5i := G2/Pi , i = 1, 2. The lattice inside the maximal linearly independent, pairwise commuting elements of G_2 are spanned by the 12 roots of g_2, the root lattice, with 2 additional roots at the origin. The roots define a hexagonal planar lattice. The root diagram is with the attached gif file. The 9-dimensional parabolic groups P_1 and P_2 subgroups of G_2 with Lie algebras p_1 and p_2. By general results, the space G_2/P_i is a compact homogeneous variety, and a projectivization of the representation space V_i with highest weight ω_i as the G_2 orbit of some distinguished vector v_i ; but ω_i generates the 7-dimensional representation, here spanned by the six short roots and zero with multiplicity one. ω_2 is the highest weight of the adjoint representation that is spanned by all roots and zero with multiplicity two.. This means that
M^5_1 = G2/P1 = G2 • [v1] ⊂ P(C^7) = CP^6,
M^5_2 = G2/P2 = G2 • [v2] ⊂ P(g_2) = CP^13.
The first space M^5_1 is thus a quadric in CP^6.
There are two real 9-dimensional subgroups P_ i inside the noncompact real form G_2 corresponding to the complex parabolic groups Pi ⊂ G2; but they have no counterparts in the compact Lie group G_2 (roughly speaking, g_2 ⊂ so(7) consists of skew symmetric matrices, while parabolics are always upper triangular): a maximal subgroup of G_2 is isomorphic to SU(3) and thus 8-dimensional. M^5_1 is a projective variety over generic three forms in 7 dimensions, where the three forms have dimension
dimG = dimGL(7,C) − dimΛ^3(C^7) = 14.
The second algebraic varieity M^5_2 = G2/P2 ⊂ CP^13 is a more complicated projective algebraic variety: it has degree 18 and its complete intersection with three hyperplanes is a K3 surface of genus 10. Geometrically G2/P2 in terms of ω is seen in the 21-dimensional representation Λ^2C^7 splits G_2 into g_2 ⊕ C_7 Consequently G2/P2 is a subvariety of P(Λ^2C^7) as well. By the Plücker embedding, the 14-dimensional Grassmann variety G(2, 7) of 2-planes in C^7 lies in P(Λ^2C^7). This then leads to a GL(7,C) group structure on 2-planes, which is a gauge-like fibration.
So G_2 can then act as an internal gauge structure which has a triality on the octonions. BTW in Conway and Sloane on page 559 there is discussion of the Mathieu group automorphism on the Fischer Greiss group with a high order triality condition. As time is getting a bit late I will have to continue this later. Yet I think that this could how the gauge structure of this might look, and it would be comparatively elementary. This three color system would then be an internal gauge structure on the octonions which results in three-way products with nonassociative structure.
Cheers LCAttachment #1: jordan.pdfAttachment #2: 3color.GIF