Masteries · Jason GlowneyArithmoí Book 1

Masteries

1.1.21 Polygonal Rotations in U(3)U(3)

It is es­sen­tial to gain an in­tu­itive grasp of how in­fin­i­tes­i­mal ro­ta­tions of the equi­lat­er­al tri­an­gle might occur. Using the ex­am­ple for our ex­po­nen­ti­at­ed el­e­ment X1X_{1}, and re­call­ing that this per­mu­ta­tion re­flects the points BB and CC of the equi­lat­er­al tri­an­gle about an axis pass­ing through the point AA, which splits the tri­an­gle into equal halves. The U(3)U(3) gen­er­a­tor G=eX1G = e^{X_{1}} acts on the iden­ti­ty el­e­ment ϵ=(A,B,C)′\epsilon = (A,B,C)', and has the matrix-vector prod­uct rep­re­sen­ta­tion

G(X1)∙ϵ=(eiθ000cos⁡(θ)isin(θ)0isin(θ)cos⁡(θ))(ABC)=(Acos⁡(θ)+iAsin(θ)Bcos(θ)+iCsin(θ)Ccos(θ)+iBsin(θ)),\begin{array}{r} G\left( X_{1} \right) \bullet \epsilon = \begin{pmatrix} e^{i\theta} & 0 & 0 \\ 0 & \cos(\theta) & isin(\theta) \\ 0 & isin(\theta) & \cos(\theta) \end{pmatrix}\begin{pmatrix} A \\ B \\ C \end{pmatrix} = \begin{pmatrix} A\cos{(\theta) + iAsin(\theta)} \\ Bcos(\theta) + iCsin(\theta) \\ Ccos(\theta) + iBsin(\theta) \end{pmatrix}, \end{array}
(1.1.124)

where for clar­i­ty, the angle of ro­ta­tion is θ\theta. The re­sult­ing matrix iden­ti­ty-vector prod­uct ro­tates the vertex AA around the unit circle in the X1X_{1}-plane by ∠θ\angle\theta, while the ver­tices B,CB,C are ro­tat­ing in the X2X3X_{2}X_{3}-planes by ∠θ\angle\theta. The ro­ta­tion for the vertex AA around its fixed ge­o­des­ic is the same ro­ta­tion of a point around the com­plex unit circle like the U(1)U(1) case, and, as such, the group U(1)U(1) is em­bed­ded in the G(X1)G\left( X_{1} \right) gen­er­a­tor by de­fault. Since the ro­ta­tion of the B,CB,C ver­tices in the X2/X3X_{2}/X_{3}-plane may be less in­tu­itive; we can con­sid­er it fur­ther in two dif­fer­ent ways. The first is to follow the real and imag­i­nary vector com­po­nents of the re­spec­tive B,CB,C ver­tices in­di­vid­u­al­ly by taking each X2X_{2} and X3X_{3}-plane sep­a­rate­ly. To con­sid­er this geo­met­ri­cal­ly, we set ∠θ=π8\angle\theta = \frac{\pi}{8}, to then cal­cu­late our matrix-vector prod­uct as

(eiπ8000cos⁡(π8)isin(π8)0isin(π8)cos⁡(π8))(ABC)=(Acos⁡(π8)+iAsin(π8)Bcos(π8)+iCsin(π8)Ccos(π8)+iBsin(π8)),\begin{array}{r} \begin{pmatrix} e^{\frac{i\pi}{8}} & 0 & 0 \\ 0 & \cos\left( \frac{\pi}{8} \right) & isin\left( \frac{\pi}{8} \right) \\ 0 & isin\left( \frac{\pi}{8} \right) & \cos\left( \frac{\pi}{8} \right) \end{pmatrix}\begin{pmatrix} A \\ B \\ C \end{pmatrix} = \begin{pmatrix} A\cos{\left( \frac{\pi}{8} \right) + iAsin\left( \frac{\pi}{8} \right)} \\ Bcos\left( \frac{\pi}{8} \right) + iCsin\left( \frac{\pi}{8} \right) \\ Ccos\left( \frac{\pi}{8} \right) + iBsin\left( \frac{\pi}{8} \right) \end{pmatrix}, \end{array}
(1.1.125)

and gath­er­ing like vertex com­po­nents by re­spec­tive planes, we have

xˆ1=(x1,ix1)=A(cos⁡(π8)+isin(π8)),\begin{array}{r} {\widehat{x}}_{1} = \left( x_{1},ix_{1} \right) = A\left( \cos{\left( \frac{\pi}{8} \right) + isin\left( \frac{\pi}{8} \right)} \right), \end{array}
(1.1.126)
xˆ2=(x2,ix2)=B(cos⁡(π8)+isin(π8)),\begin{array}{r} {\widehat{x}}_{2} = \left( x_{2},ix_{2} \right) = B\left( \cos{\left( \frac{\pi}{8} \right) + isin\left( \frac{\pi}{8} \right)} \right), \end{array}
(1.1.127)
xˆ3=(x3,ix3)=C(cos⁡(π8)+isin(π8)).\begin{array}{r} {\widehat{x}}_{3} = \left( x_{3},ix_{3} \right) = C\left( \cos{\left( \frac{\pi}{8} \right) + isin\left( \frac{\pi}{8} \right)} \right). \end{array}
(1.1.128)
Figure
Figure 1.1.12 – G(X1(iπ8))∙ϵG\left( X_{1}\left( \frac{i\pi}{8} \right) \right) \bullet \epsilon for xˆj{\widehat{x}}_{j} vec­tors.

In Figure 1.1.12, we graph these xˆj{\widehat{x}}_{j} vec­tors, seeing that they main­tain the polyg­o­nal form of the equi­lat­er­al tri­an­gle, pic­tured in blue, with ro­tat­ing ver­tices on the sphere's sur­face. The re­la­tion­ship in Eq. 1.1.25 is the result of taking the iden­ti­ty ϵ\epsilon to be equal to the fol­low­ing 3×33 \times 3 matrix

ϵ=(A000B000C),\begin{array}{r} \epsilon = \begin{pmatrix} A & 0 & 0 \\ 0 & B & 0 \\ 0 & 0 & C \end{pmatrix}, \end{array}
(1.1.129)

and G(X1(iπ8))∙ϵG\left( X_{1}\left( \frac{i\pi}{8} \right) \right) \bullet \epsilon be­comes

G(X1(iπ8))∙ϵ=(eiπ8000cos⁡(π8)isin(π8)0isin(π8)cos⁡(π8))(A000B000C)=(Aeiπ8000Bcos(π8)iCsin(π8)0iBsin(π8)Ccos⁡(π8)),\begin{array}{r} G\left( X_{1}\left( \frac{i\pi}{8} \right) \right) \bullet \epsilon = \begin{pmatrix} e^{\frac{i\pi}{8}} & 0 & 0 \\ 0 & \cos\left( \frac{\pi}{8} \right) & isin\left( \frac{\pi}{8} \right) \\ 0 & isin\left( \frac{\pi}{8} \right) & \cos\left( \frac{\pi}{8} \right) \end{pmatrix}\begin{pmatrix} A & 0 & 0 \\ 0 & B & 0 \\ 0 & 0 & C \end{pmatrix} = \begin{pmatrix} Ae^{\frac{i\pi}{8}} & 0 & 0 \\ 0 & Bcos\left( \frac{\pi}{8} \right) & iCsin\left( \frac{\pi}{8} \right) \\ 0 & iBsin\left( \frac{\pi}{8} \right) & C\cos\left( \frac{\pi}{8} \right) \end{pmatrix}, \end{array}
(1.1.130)

where each column rep­re­sents the re­spec­tive xˆj{\widehat{x}}_{j} axis vector. When con­sid­er­ing the matrix gen­er­a­tor G(X1(iπ8))G\left( X_{1}\left( \frac{i\pi}{8} \right) \right) acting on a vector by left mul­ti­pli­ca­tion, per­form­ing a sin­gu­lar value de­com­po­si­tion of the matrix can give us a better idea of what steps happen to make this ro­ta­tion occur. The reader can refer to Sec­tion 1.1.8 to review the SVD of a matrix if needed. As G(X1(iπ8))G\left( X_{1}\left( \frac{i\pi}{8} \right) \right) is di­ag­o­nal­iz­able, we take our Σ\Sigma to be a 3×33 \times 3 matrix with the eigen­val­ues along the di­ag­o­nal, and our uni­tary ma­tri­ces, V,V†V,V^{\dagger}, are just the eigen­vec­tors arranged as columns in the VV matrix, co­in­cid­ing with the re­spec­tive σj\sigma_{j} en­tries in Σ\Sigma. The SVD can then be ex­pressed as

G(X1(θ))=VΣV†,\begin{array}{r} G\left( X_{1}(\theta) \right) = V\Sigma V^{\dagger}, \end{array}
(1.1.131)

where Σ,V,V†\Sigma,V,V^{\dagger} are

Σ=(eiθ000e−iθ000eiθ),   V=(1000−121201212),   V†=(1000−121201212).\begin{array}{r} \Sigma = \begin{pmatrix} e^{i\theta} & 0 & 0 \\ 0 & e^{- i\theta} & 0 \\ 0 & 0 & e^{i\theta} \end{pmatrix},\ \ \ V = \begin{pmatrix} 1 & 0 & 0 \\ 0 & - \frac{1}{\sqrt{2}} & \frac{1}{\sqrt{2}} \\ 0 & \frac{1}{\sqrt{2}} & \frac{1}{\sqrt{2}} \end{pmatrix},\ \ \ V^{\dagger} = \begin{pmatrix} 1 & 0 & 0 \\ 0 & - \frac{1}{\sqrt{2}} & \frac{1}{\sqrt{2}} \\ 0 & \frac{1}{\sqrt{2}} & \frac{1}{\sqrt{2}} \end{pmatrix}. \end{array}
(1.1.132)

Above the eigen­vec­tors in VV are nor­mal­ized, and VV is equal to its con­ju­gate trans­pose V†V^{\dagger}. Set­ting our iden­ti­ty vector to ϵˆ=(A,B,C)T\widehat{\epsilon} = (A,B,C)^{T}, for A,B,C=1A,B,C = 1, choos­ing an angle value θ=π4\theta = \frac{\pi}{4}, and then cal­cu­lat­ing the prod­uct V†ϵˆV^{\dagger}\widehat{\epsilon}, we find

V†ϵˆ=(AC2−B2C2+B2)=(102)=ϵˆ′=(A′B′C′).\begin{array}{r} V^{\dagger}\widehat{\epsilon} = \begin{pmatrix} A \\ \frac{C}{\sqrt{2}} - \frac{B}{\sqrt{2}} \\ \frac{C}{\sqrt{2}} + \frac{B}{\sqrt{2}} \end{pmatrix} = \begin{pmatrix} 1 \\ 0 \\ \sqrt{2} \end{pmatrix} = {\widehat{\epsilon}}^{'} = \begin{pmatrix} A' \\ B' \\ C' \end{pmatrix}. \end{array}
(1.1.133)

This shows that BB is taken to B′=(C2 −B2 )=0B^{'} = \left( \frac{C}{\sqrt{2}\ } - \frac{B}{\sqrt{2}\ } \right) = 0, the origin CC is taken to be C′=(C2 +B2 )=2C^{'} = \left( \frac{C}{\sqrt{2}\ } + \frac{B}{\sqrt{2}\ } \right) = \sqrt{2}, and that the AA vertex re­mains un­changed A′=AA' = A. Since A,B,CA,B,C are uni­tary in mag­ni­tude, the two-norm of the re­sult­ing vector |ϵˆ′|\left| {\widehat{\epsilon}}^{'} \right|is the same as the iden­ti­ty vector's two-norm |ϵˆ|\left| \widehat{\epsilon} \right|, where both are equal to 3\sqrt{3}

|ϵˆ|=|ϵˆ′|→A2+B2+C2=A2+(C2−B2)2+(C2+B2)2=3=3.\begin{array}{r} \left| \widehat{\epsilon} \right| = \left| {\widehat{\epsilon}}^{'} \right| \rightarrow \sqrt{A^{2} + B^{2} + C^{2}} = \sqrt{A^{2} + \left( \frac{C}{\sqrt{2}} - \frac{B}{\sqrt{2}} \right)^{2} + \left( \frac{C}{\sqrt{2}} + \frac{B}{\sqrt{2}} \right)^{2}} = \sqrt{3} = \sqrt{3}. \end{array}
(1.1.134)

For our case, the two-norm is the mag­ni­tude of a vector ∈ℂ3\in \mathbb{C}^{3} with com­po­nents A,B,CA,B,C, rep­re­sent­ing the co­or­di­nate values of the tri­an­gle's ver­tices in the X1,X2,X3X_{1},X_{2},X_{3}-planes re­spec­tive­ly. So, the result of V†ϵˆV^{\dagger}\widehat{\epsilon} is to rotate this tri­an­gle, whose ver­tices rep­re­sent a vector ∈ℂ3\in \mathbb{C}^{3}, with an ini­tial vector start­ing po­si­tion ϵˆ=(1,1,1)\widehat{\epsilon} = (1,1,1), to the co­or­di­nates ϵˆ′=(1,0,2){\widehat{\epsilon}}^{'} = \left( 1,0,\sqrt{2} \right). The ro­ta­tion of the tri­an­gle and its vector rep­re­sen­ta­tion re­sult­ing from V†ϵˆV^{\dagger}\widehat{\epsilon}, is de­pict­ed in Figure 1.1.13, in which the start­ing tri­an­gle-11 with vertex co­or­di­nates (1,1,1)(1,1,1) has the vector rep­re­sen­ta­tion de­pict­ed as (A,B,C)(A,B,C), which both rotate as the result of the V†ϵˆV^{\dagger}\widehat{\epsilon} prod­uct to the tri­an­gle-22, de­pict­ed at co­or­di­nates (1,0,2)(1,0,\sqrt{2}), along with the re­sult­ing vector rep­re­sen­ta­tion (A′,B′,C′).(A',B',C').This is con­sis­tent with what is the ex­pect­ed out­come from the V†V^{\dagger} por­tion of the SVD, which acts to rotate the vector/tri­an­gle by ro­tat­ing the sphere.

Figure
Figure 1.1.13 – The V†ϵˆV^{\dagger}\widehat{\epsilon} ro­ta­tion of the tri­an­gle and its vector rep­re­sen­ta­tion ∈ℂ3\in \mathbb{C}^{3}.

The next step of the SVD is left mul­ti­pli­ca­tion by the Σ\Sigma matrix on the ϵˆ′\widehat{\epsilon}'vector, which yields the ϵˆ′′{\ \widehat{\epsilon}}^{''} vector

Σϵˆ′=(Aeiπ4 (C2−B2)e−iπ4(C2+B2)eiπ4)=(12+i201+i)=ϵˆ′′=(A′′B′′C′′).\begin{array}{r} \Sigma{\widehat{\epsilon}}^{'} = \begin{pmatrix} Ae^{\frac{i\pi}{4}\ } \\ \left( \frac{C}{\sqrt{2}} - \frac{B}{\sqrt{2}} \right)e^{- \frac{i\pi}{4}} \\ \left( \frac{C}{\sqrt{2}} + \frac{B}{\sqrt{2}} \right)e^{\frac{i\pi}{4}} \end{pmatrix} = \begin{pmatrix} \frac{1}{\sqrt{2}} + \frac{i}{\sqrt{2}} \\ 0 \\ 1 + i \end{pmatrix} = {\widehat{\epsilon}}^{''} = \begin{pmatrix} A'' \\ B'' \\ C'' \end{pmatrix}. \end{array}
(1.1.135)

Given that Σ\Sigma is a di­ag­o­nal­ized matrix of the eigen­val­ues of G(X1(θ=π4))G\left( X_{1}\left( \theta = \frac{\pi}{4} \right) \right), the Σϵˆ′\Sigma{\widehat{\epsilon}}^{'} matrix-vector prod­uct scales ϵˆ′{\widehat{\epsilon}}^{'}, and ac­cord­ing­ly, its geo­met­ric tri­an­gle/vector rep­re­sen­ta­tions. As the eigen­val­ues are com­plex and uni­tary, the vector ϵˆ′\widehat{\epsilon}'is ro­tat­ed and scaled, as op­posed to solely the scal­ing seen for a Σ\Sigma-matrix com­posed of real eigen­val­ues.

Figure
Figure 1.1.14 – The SVD de­com­po­si­tion of G(X1(θ))ϵˆG\left( X_{1}(\theta) \right)\widehat{\epsilon}, (A,B,C)=ϵˆ(A,B,C) = \widehat{\epsilon}: tri­an­gle-11, V†ϵˆ=(A′,B′,C′)V^{\dagger}\widehat{\epsilon} = (A',B',C'):

tri­an­gle-2, Σϵˆ′=(A′′,B′′,C′′)\Sigma{\widehat{\epsilon}}^{'} = (A'',B'',C''): tri­an­gle-3, Σϵˆ′′(4)=(A′′′,B′′′,C′′′)T=G(X1(θ))ϵˆ\Sigma{\widehat{\epsilon}}^{''}(4) = (A''',B''',C''')^{T} = G\left( X_{1}(\theta) \right)\widehat{\epsilon}: tri­an­gle-4.

In Figure 1.1.14, we track the path of the ini­tial vector points A,B,CA,B,C, rep­re­sent­ing the tri­an­gle's ver­tices through the SVD de­com­po­si­tion G(X1(θ=π4))=VΣV†.G\left( X_{1}\left( \theta = \frac{\pi}{4} \right) \right) = V\Sigma V^{\dagger}.The matrix-vector prod­uct Vϵˆ′′V\widehat{\epsilon}''brings us to the ul­ti­mate vector result, ϵˆ′′′{\widehat{\epsilon}}^{'''}.

Vϵˆ′′=(Aeiπ4 Bcos(π4)+iCsin(π4)Ccos(π4)+iBsin(π4))=(12+i212+i212+i2)=(A′′′B′′′C′′′)=ϵˆ′′′,\begin{array}{r} V{\widehat{\epsilon}}^{''} = \begin{pmatrix} Ae^{\frac{i\pi}{4}\ } \\ Bcos\left( \frac{\pi}{4} \right) + iCsin\left( \frac{\pi}{4} \right) \\ Ccos\left( \frac{\pi}{4} \right) + iBsin\left( \frac{\pi}{4} \right) \end{pmatrix} = \begin{pmatrix} \frac{1}{\sqrt{2}} + \frac{i}{\sqrt{2}} \\ \frac{1}{\sqrt{2}} + \frac{i}{\sqrt{2}} \\ \frac{1}{\sqrt{2}} + \frac{i}{\sqrt{2}} \end{pmatrix} = \begin{pmatrix} A''' \\ B''' \\ C''' \end{pmatrix} = {\widehat{\epsilon}}^{'''}, \end{array}
(1.1.136)

rep­re­sent­ing the linear trans­for­ma­tion of the ϵˆ\widehat{\epsilon} vector, the result of the matrix-vector prod­uct G(X1(θ=π4))ϵˆG\left( X_{1}\left( \theta = \frac{\pi}{4} \right) \right)\widehat{\epsilon}.