Masteries · Jason GlowneyArithmoí Book 1

Masteries

1.6.5 Outer Products and Infinities

Math­e­mati­cians and physi­cists tend to reduce math­e­mat­ics equa­tions into the sim­plest and most com­pact forms, but some­times taking an ex­pand­ed ap­proach to these equa­tions can yield valu­able in­sights, hint­ing at an el­e­gant struc­ture in the math­e­mat­ics that might have oth­er­wise been ob­scured by re­duc­tion­ism.

An outer prod­uct is the prod­uct of two co­or­di­nate vec­tors that form a matrix and is re­lat­ed to tensor al­ge­bra. If we have two vec­tors uiˆ, vjˆ\widehat{u_{i}},\ \widehat{v_{j}}, for uˆi{\widehat{u}}_{i}, a column vector and vˆj{\widehat{v}}_{j}, a row vector, their tensor prod­uct, as uˆ⨂vˆ\widehat{u}\bigotimes_{}^{}\widehat{v}, is an i×ji \times j matrix, wijw_{ij}. So, for vec­tors uˆi=(a1, a2, a3, a4)T{\widehat{u}}_{i} = \left( a_{1},\ a_{2},\ a_{3},\ a_{4} \right)^{T}, vˆj=(b1, b2, b3, b4){\widehat{v}}_{j} = \left( b_{1},\ b_{2},\ b_{3},\ b_{4} \right), their tensor prod­uct is

uˆ⊗vˆ=(a1a2a2a4)(b1b2b3b4)=[a1b1a1b2a1b3a1b4a2b1a2b2a2b3a2b4a3b1a3b2a3b3a3b4a4b1a4b2a4b3a4b4]=wij.\begin{array}{r} \widehat{u} \otimes \widehat{v} = \begin{pmatrix} a_{1} \\ a_{2} \\ a_{2} \\ a_{4} \end{pmatrix}\begin{pmatrix} b_{1} & b_{2} & b_{3} & b_{4} \end{pmatrix} = \begin{bmatrix} a_{1}b_{1} & a_{1}b_{2} & a_{1}b_{3} & a_{1}b_{4} \\ a_{2}b_{1} & a_{2}b_{2} & a_{2}b_{3} & a_{2}b_{4} \\ a_{3}b_{1} & a_{3}b_{2} & a_{3}b_{3} & a_{3}b_{4} \\ a_{4}b_{1} & a_{4}b_{2} & a_{4}b_{3} & a_{4}b_{4} \end{bmatrix} = w_{ij}. \end{array}
(1.6.120)

An im­por­tant dis­tinc­tion, re­vers­ing the order of vec­tors in the outer prod­uct will result in a dif­fer­ent matrix result. The matrix re­sults will have de­ter­mi­nants equal to zero, for outer prod­ucts formed by two vec­tors of the same di­men­sions greater than one. Ten­sors are vital to math­e­mat­i­cal physics and are used fre­quent­ly in quan­tum me­chan­ics to de­scribe quan­tum system in­ter­ac­tions very well, where the goal in their use is to de­scribe how basic things in­ter­act with one an­oth­er.

When it comes to the in­fi­nite series that en­com­pass­es the ex­po­nen­tial map­ping in Lie theory, tensor al­ge­bras can be con­sid­ered by taking these in­fi­nite series to be co­or­di­nate vec­tors. By doing so, we can take a more ex­pan­sive look at how these in­fi­nite math­e­mat­i­cal forms in­ter­act with each other. The in­fi­nite series that con­verges to eixe^{ix}, is com­posed of two dis­tinct in­fi­nite series, the in­fi­nite power series and rec­i­p­ro­cal fac­to­r­i­al series, whose dot prod­uct leads to Euler's for­mu­la

eix=∑n=0∞(ix)n∙1n!=1+ix−x22−ix36+x424+…,\begin{array}{r} e^{ix} = \sum_{n = 0}^{\infty}{(ix)^{n} \bullet \frac{1}{n!} = 1 + ix - \frac{x^{2}}{2} - \frac{ix^{3}}{6} + \frac{x^{4}}{24} + \ldots}, \end{array}
(1.6.121)

Euler's for­mu­la is com­posed of a geo­met­ric series por­tion, being

∑n=0∞(ix)n=1+ix−x2−ix3+x4+…=i(i+x),   (for |ix|<1) ,\begin{array}{r} \sum_{n = 0}^{\infty}{(ix)^{n} = 1 + ix - x^{2} - ix^{3} + x^{4} + \ldots = \frac{i}{(i + x)},\ \ \ \left( for\ |ix| < 1 \right)\ ,} \end{array}
(1.6.122)

and the rec­i­p­ro­cal fac­to­r­i­al series por­tion

∑n=0∞1n!=1+1+12+16+124+…=e.\begin{array}{r} \sum_{n = 0}^{\infty}\frac{1}{n!} = 1 + 1 + \frac{1}{2} + \frac{1}{6} + \frac{1}{24} + \ldots = e. \end{array}
(1.6.123)

From the two dis­tinct series, we set umu_{m} equal to the sum of the geo­met­ric power series

um=∑m=0∞(ix)m,\begin{array}{r} u_{m} = \sum_{m = 0}^{\infty}{(ix)^{m},} \end{array}
(1.6.124)

and vnv_{n}, as the sum of the in­fi­nite rec­i­p­ro­cal fac­to­r­i­al series

vn=∑n=0∞1n!.\begin{array}{r} v_{n} = \sum_{n = 0}^{\infty}\frac{1}{n!}. \end{array}
(1.6.125)

We define the column vector uˆm{\widehat{u}}_{m}, and the row vector vˆn{\widehat{v}}_{n}, as in­fi­nite-di­men­sion­al, com­posed of the el­e­ments from um,vnu_{m},v_{n}

 uˆm=(1,ix,−x2,−ix3,x4,ix5,−x6,−ix7,…),\begin{array}{r} \ {\widehat{u}}_{m} = \left( 1,ix, - x^{2}, - ix^{3},x^{4},ix^{5}, - x^{6}, - ix^{7},\ldots \right), \end{array}
(1.6.126)
vˆn=(1,1,12,16,124,1120,1720,15040,…),\begin{array}{r} {\widehat{v}}_{n} = \left( 1,1,\frac{1}{2},\frac{1}{6},\frac{1}{24},\frac{1}{120},\frac{1}{720},\frac{1}{5040},\ldots \right), \end{array}
(1.6.127)

and the outer prod­uct of these vec­tors uˆ⊗vˆ\widehat{u} \otimes \widehat{v}, form the matrix XX

X=uˆ⊗vˆ=(𝟏112161241120172015040⋯ix𝐢𝐱ix2ix6ix24ix120ix720ix5040⋯−x2−x2−𝐱𝟐𝟐−x26−x224−x2120−x2720−x25040⋯−ix3−ix3−ix32−𝐢𝐱𝟑𝟔−ix324−ix3120−ix3720−ix35040⋯x4x4x42x46𝐱𝟒𝟐𝟒x4120x4720x45040⋯ix5ix5ix52ix56ix524𝐢𝐱𝟓𝟏𝟐𝟎ix5720ix55040⋯−x6−x6−x62−x66−x624−x6120−𝐱𝟔𝟕𝟐𝟎−x65040⋯−ix7−ix7−ix72−ix76−ix724−ix7120−ix7720−𝐢𝐱𝟕𝟓𝟎𝟒𝟎⋯⋮⋮⋮⋮⋮⋮⋮⋮⋱).\begin{array}{r} X = \widehat{u} \otimes \widehat{v} = \begin{pmatrix} \mathbf{1} & 1 & \frac{1}{2} & \frac{1}{6} & \frac{1}{24} & \frac{1}{120} & \frac{1}{720} & \frac{1}{5040} & \cdots \\ ix & \mathbf{ix} & \frac{ix}{2} & \frac{ix}{6} & \frac{ix}{24} & \frac{ix}{120} & \frac{ix}{720} & \frac{ix}{5040} & \cdots \\ - x^{2} & - x^{2} & \mathbf{-}\frac{\mathbf{x}^{\mathbf{2}}}{\mathbf{2}} & - \frac{x^{2}}{6} & - \frac{x^{2}}{24} & - \frac{x^{2}}{120} & - \frac{x^{2}}{720} & - \frac{x^{2}}{5040} & \cdots \\ - ix^{3} & - ix^{3} & - \frac{ix^{3}}{2} & \mathbf{-}\frac{\mathbf{i}\mathbf{x}^{\mathbf{3}}}{\mathbf{6}} & - \frac{ix^{3}}{24} & - \frac{ix^{3}}{120} & - \frac{ix^{3}}{720} & - \frac{ix^{3}}{5040} & \cdots \\ x^{4} & x^{4} & \frac{x^{4}}{2} & \frac{x^{4}}{6} & \frac{\mathbf{x}^{\mathbf{4}}}{\mathbf{24}} & \frac{x^{4}}{120} & \frac{x^{4}}{720} & \frac{x^{4}}{5040} & \cdots \\ ix^{5} & ix^{5} & \frac{ix^{5}}{2} & \frac{ix^{5}}{6} & \frac{ix^{5}}{24} & \frac{\mathbf{i}\mathbf{x}^{\mathbf{5}}}{\mathbf{120}} & \frac{ix^{5}}{720} & \frac{ix^{5}}{5040} & \cdots \\ - x^{6} & - x^{6} & - \frac{x^{6}}{2} & - \frac{x^{6}}{6} & - \frac{x^{6}}{24} & - \frac{x^{6}}{120} & \mathbf{-}\frac{\mathbf{x}^{\mathbf{6}}}{\mathbf{720}} & - \frac{x^{6}}{5040} & \cdots \\ - ix^{7} & - ix^{7} & - \frac{ix^{7}}{2} & - \frac{ix^{7}}{6} & - \frac{ix^{7}}{24} & - \frac{ix^{7}}{120} & - \frac{ix^{7}}{720} & \mathbf{-}\frac{\mathbf{i}\mathbf{x}^{\mathbf{7}}}{\mathbf{5040}} & \cdots \\ \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \ddots \end{pmatrix}. \end{array}
(1.6.128)

The first thing that stands out is that the trace of the XX matrix, tr(X)tr(X), shown as the bold el­e­ments along the main di­ag­o­nal, is the in­fi­nite series for eixe^{ix}, tr(X)=eixtr(X) = e^{ix}. In terms of the eigen­val­ues and eigen­vec­tors of the XXmatrix, there is a single eigen­val­ue λ\lambda, equal to eixe^{ix}

λ=1+ix−x22−ix36+x424+ix5120−x6720−ix75040+…=eix,\begin{array}{r} \lambda = 1 + ix - \frac{x^{2}}{2} - \frac{ix^{3}}{6} + \frac{x^{4}}{24} + \frac{ix^{5}}{120} - \frac{x^{6}}{720} - \frac{ix^{7}}{5040} + \ldots = e^{ix}, \end{array}
(1.6.129)

that would be in an eigen­val­ue matrix DD, in the row/column po­si­tion d∞,∞d_{\infty,\infty}. Note­wor­thy, is that the eigen­val­ue λ=eix\lambda = e^{ix}, can take on an in­fi­nite number of values for any ar­gu­ment for x∈ℝx\in \mathbb{R}, on the unit circle, en­com­pass­ing the essence of Lie theory and the in­fin­i­tes­i­mal ro­ta­tions. In con­trast to the lone eigen­val­ue, there are an in­fi­nite number of eigen­vec­tors vˆi{\widehat{v}}_{i}, rep­re­sent­ing the columns in the eigen­vec­tor matrix VV

V=(−1−12−16−124−1120−1720−15040…(ix)−∞10000000(ix)−∞+101000000(ix)−∞+200100000(ix)−∞+300010000(ix)−∞+400001000(ix)−∞+500000100(ix)−∞+600000010(ix)−∞+7⋮⋮⋮⋮⋮⋮⋮⋱⋮00000000ix).\begin{array}{r} V = \begin{pmatrix} \begin{matrix} - 1 & - \frac{1}{2} & - \frac{1}{6} & - \frac{1}{24} & - \frac{1}{120} & - \frac{1}{720} & - \frac{1}{5040} & \ldots & (ix)^{- \infty} \\ 1 & 0 & 0 & 0 & 0 & 0 & 0 & 0 & (ix)^{- \infty + 1} \\ 0 & 1 & 0 & 0 & 0 & 0 & 0 & 0 & (ix)^{- \infty + 2} \\ 0 & 0 & 1 & 0 & 0 & 0 & 0 & 0 & (ix)^{- \infty + 3} \\ 0 & 0 & 0 & 1 & 0 & 0 & 0 & 0 & (ix)^{- \infty + 4} \\ 0 & 0 & 0 & 0 & 1 & 0 & 0 & 0 & (ix)^{- \infty + 5} \\ 0 & 0 & 0 & 0 & 0 & 1 & 0 & 0 & (ix)^{- \infty + 6} \\ 0 & 0 & 0 & 0 & 0 & 0 & 1 & 0 & (ix)^{- \infty + 7} \\ \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 & ix \end{matrix} \end{pmatrix}. \end{array}
(1.6.130)

The eigen­vec­tor columns in VV, except for the last column, arrange in a pat­tern con­sis­tent with the first row equal­ing the fac­to­r­i­al el­e­ment −1k!- \frac{1}{k!}, for kk, the column number, and the number one in the row number j=k+1j = k + 1. The last eigen­vec­tor column as a sum, rep­re­sents a geo­met­ric power series, which, for |x|>1|x| > 1, con­verges to

∑k=1∞(ix)−k=1ix−1,   (for  |x|>1).\begin{array}{r} \sum_{k = 1}^{\infty}(ix)^{- k} = \frac{1}{ix - 1},\ \ \ \left( for\ \ |x| > 1 \right). \end{array}
(1.6.131)

If we return and con­sid­er the in­fi­nite sum of the el­e­ments ∈X\in X, it is easily seen that each row in the matrix will sum to e(ix)j−1e(ix)^{j - 1}, for j,j,the row number. This leads to a con­ver­gent in­fi­nite row sum for XX, when |x|<1|x| < 1

∑k=0∞e(ix)k=iex+i,   (for |x|<1),\begin{array}{r} \sum_{k = 0}^{\infty}{e(ix)^{k}} = \frac{ie}{x + i},\ \ \ \left( for\ |x| < 1 \right), \end{array}
(1.6.132)

and is con­sis­tent with the prod­uct result we would find be­tween the sum­ma­tion of the vector el­e­ments from uiˆ, vjˆ\widehat{u_{i}},\ \widehat{v_{j}}

(∑uiˆ,)(∑ vjˆ,)=iex+i,   (for |x|<1).\begin{array}{r} \left( \sum_{}^{}{\widehat{u_{i}},} \right)\left( \sum_{}^{}{\ \widehat{v_{j}},} \right) = \frac{ie}{x + i},\ \ \ \left( for\ |x| < 1 \right). \end{array}
(1.6.133)

From the stand­point of matrix theory, the el­e­ments on the main di­ag­o­nal rep­re­sent the sym­met­ric por­tion of the prod­uct uˆ⊗vˆ\widehat{u} \otimes \widehat{v}, which sum to eixe^{ix}, while other el­e­ments then rep­re­sent the an­ti­sym­met­ric por­tion of the prod­uct. Taken to­geth­er, the sym­met­ric and an­ti­sym­met­ric el­e­ments are prod­uct rep­re­sen­ta­tions for every el­e­ment in each in­fi­nite series with one an­oth­er, rep­re­sent­ing columns, rows, and off-di­ag­o­nals, which are also in­fi­nite series. To better un­der­stand the sym­met­ric, and asym­met­ric por­tions of XX, it is im­por­tant to re­al­ize that uˆ⊗vˆ≠vˆ⊗uˆ\widehat{u} \otimes \widehat{v} \neq \widehat{v} \otimes \widehat{u}, mean­ing our vec­tors are non-com­mu­ta­tive as outer prod­ucts. This be­comes clear by show­ing an ab­bre­vi­at­ed ver­sion of the outer prod­uct matrix vˆ⊗uˆ\widehat{v} \otimes \widehat{u},

X′=vˆ⊗uˆ=(1ix−x2−ix3x4⋯1ix−x2−ix3x4⋯12ix2−x22−ix32x42⋯16ix6−x26−ix36x46⋯124ix24−x224−ix324x424…⋮⋮⋮⋮⋮⋱),\begin{array}{r} X^{'} = \widehat{v} \otimes \widehat{u} = \begin{pmatrix} 1 & ix & - x^{2} & - ix^{3} & x^{4} & \cdots \\ 1 & ix & - x^{2} & - ix^{3} & x^{4} & \cdots \\ \frac{1}{2} & \frac{ix}{2} & - \frac{x^{2}}{2} & - \frac{ix^{3}}{2} & \frac{x^{4}}{2} & \cdots \\ \frac{1}{6} & \frac{ix}{6} & - \frac{x^{2}}{6} & - \frac{ix^{3}}{6} & \frac{x^{4}}{6} & \cdots \\ \frac{1}{24} & \frac{ix}{24} & - \frac{x^{2}}{24} & - \frac{ix^{3}}{24} & \frac{x^{4}}{24} & \ldots \\ \vdots & \vdots & \vdots & \vdots & \vdots & \ddots \end{pmatrix}, \end{array}
(1.6.134)

for which the el­e­ments on the main di­ag­o­nal are the same as the first outer prod­uct uˆ⊗vˆ\widehat{u} \otimes \widehat{v}, defin­ing the sym­met­ric por­tion.

The an­ti­sym­met­ric por­tion of vˆ⊗uˆ,\widehat{v} \otimes \widehat{u},are the non-di­ag­o­nal el­e­ments ∈X′\in X^{'}, and are trans­posed in com­par­i­son to those in XX. This leads us back to the geo­met­ric prod­uct of vec­tors, to con­sid­er its sym­met­ric dot prod­uct and an­ti­sym­met­ric wedge prod­uct in terms of our in­fi­nite vector func­tion's outer prod­uct. The trace of our outer prod­uct matrix is con­sis­tent with an inner prod­uct, and could also be for the dot prod­uct, if we were to define the vec­tors as in­fi­nite in their co­or­di­nates. To con­sid­er the wedge prod­uct in re­la­tion to the outer prod­uct XX, we con­sid­er the ex­pres­sion that re­lates the wedge to outer prod­uct

uˆ∧vˆ=uˆ⊗vˆ−vˆ⊗uˆ,\begin{array}{r} \widehat{u} \land \widehat{v} = \widehat{u} \otimes \widehat{v} - \widehat{v} \otimes \widehat{u}, \end{array}
(1.6.135)

which, for the dif­fer­ence of outer prod­ucts should remind us of the com­mu­ta­tor op­er­a­tion. So, in terms of a geo­met­ric al­ge­bra analog, we can ex­press the geom­e­try prod­uct of the vec­tors uˆ,vˆ,\widehat{u},\widehat{v},as

uˆvˆ=uˆ∙vˆ+uˆ∧vˆ=tr(vˆ⊗uˆ)+(uˆ⊗vˆ−vˆ⊗uˆ),\widehat{u}\widehat{v} = \widehat{u} \bullet \widehat{v} + \widehat{u} \land \widehat{v} = tr\left( \widehat{v} \otimes \widehat{u} \right) + \left( \widehat{u} \otimes \widehat{v} - \widehat{v} \otimes \widehat{u} \right),

which re­sults in

X=uˆvˆ=eix+(01−ix12+x216+ix3124−x4⋯ix−10ix2+x2ix6+ix3ix24−x4⋯−x2−12−x2−ix20−x26+ix32−x224−x42⋯−ix3−16−ix3−ix6−ix32+x260−ix324−x46⋯x4−124x4−ix24x42+x224x46+ix3240⋯⋮⋮⋮⋮⋮⋱),\begin{array}{r} \mathrm{X} = \widehat{u}\widehat{v} = e^{ix} + \begin{pmatrix} 0 & 1 - ix & \frac{1}{2} + x^{2} & \frac{1}{6} + ix^{3} & \frac{1}{24} - x^{4} & \cdots \\ ix - 1 & 0 & \frac{ix}{2} + x^{2} & \frac{ix}{6} + ix^{3} & \frac{ix}{24} - x^{4} & \cdots \\ - x^{2} - \frac{1}{2} & - x^{2} - \frac{ix}{2} & 0 & - \frac{x^{2}}{6} + \frac{ix^{3}}{2} & - \frac{x^{2}}{24} - \frac{x^{4}}{2} & \cdots \\ - ix^{3} - \frac{1}{6} & - ix^{3} - \frac{ix}{6} & - \frac{ix^{3}}{2} + \frac{x^{2}}{6} & 0 & - \frac{ix^{3}}{24} - \frac{x^{4}}{6} & \cdots \\ x^{4} - \frac{1}{24} & x^{4} - \frac{ix}{24} & \frac{x^{4}}{2} + \frac{x^{2}}{24} & \frac{x^{4}}{6} + \frac{ix^{3}}{24} & 0 & \cdots \\ \vdots & \vdots & \vdots & \vdots & \vdots & \ddots \end{pmatrix}, \end{array}
(1.6.136)

in which, the sym­met­ric por­tion of the outer prod­uct ma­tri­ces X,X′X,X', van­ish­es in their dif­fer­ence, and we now have ad­di­tive in­verse re­la­tion­ships be­tween en­tries and their trans­pos­es. Given this ad­di­tive in­verse re­la­tion­ship, we now find the grand sum of the en­tries ∈\in X\mathrm{X}, equals zero.

Re­turn­ing to the XX matrix, we can apply matrix group theory to our ex­po­nen­tial series tensor, by defin­ing a di­ag­o­nal­ized matrix DD, rep­re­sent­ing the di­ag­o­nal el­e­ments ∈\in X

D=diag(x0, ix1,−x22!,−ix33!,x44!,ix55!,−x66!,−ix77!…).\begin{array}{r} D = diag\left( x^{0},\ {ix}^{1}, - \frac{x^{2}}{2!}, - \frac{{ix}^{3}}{3!},\frac{x^{4}}{4!},\frac{{ix}^{5}}{5!}, - \frac{x^{6}}{6!}, - \frac{ix^{7}}{7!}\ldots \right). \end{array}
(1.6.137)

We en­ter­tain the matrix-prod­uct D−1XD^{- 1}X , to find the re­sult­ing matrix SS, which is a com­po­si­tion of in­vers­es ∈ℂ\in \mathbb{C}, given that each entry's trans­pose is its in­verse, and this is a bi­jec­tive map­ping, with a matrix prod­uct in­verse SD=X.SD = X.

S=D−1X=(1−ix−1x2ix31x4−ix5−1x6ix7⋯ix1−ix−1x2ix31x4−ix5−1x6⋯−x2ix1−ix−1x2ix31x4−ix5⋯−ix3−x2ix1−ix−1x2ix31x4⋯x4−ix3−x2ix1−ix−1x2ix3⋯ix5x4−ix3−x2ix1−ix−1x2⋯−x6ix5x4−ix3−x2ix1−ix⋯−ix7−x6ix5x4−ix3−x2ix1⋯⋮⋮⋮⋮⋮⋮⋮⋮⋱).\begin{array}{r} S = D^{- 1}X = \begin{pmatrix} 1 & - \frac{i}{x} & - \frac{1}{x^{2}} & \frac{i}{x^{3}} & \frac{1}{x^{4}} & - \frac{i}{x^{5}} & - \frac{1}{x^{6}} & \frac{i}{x^{7}} & \cdots \\ ix & 1 & - \frac{i}{x} & - \frac{1}{x^{2}} & \frac{i}{x^{3}} & \frac{1}{x^{4}} & - \frac{i}{x^{5}} & - \frac{1}{x^{6}} & \cdots \\ - x^{2} & ix & 1 & - \frac{i}{x} & - \frac{1}{x^{2}} & \frac{i}{x^{3}} & \frac{1}{x^{4}} & - \frac{i}{x^{5}} & \cdots \\ - ix^{3} & - x^{2} & ix & 1 & - \frac{i}{x} & - \frac{1}{x^{2}} & \frac{i}{x^{3}} & \frac{1}{x^{4}} & \cdots \\ x^{4} & - ix^{3} & - x^{2} & ix & 1 & - \frac{i}{x} & - \frac{1}{x^{2}} & \frac{i}{x^{3}} & \cdots \\ ix^{5} & x^{4} & - ix^{3} & - x^{2} & ix & 1 & - \frac{i}{x} & - \frac{1}{x^{2}} & \cdots \\ - x^{6} & ix^{5} & x^{4} & - ix^{3} & - x^{2} & ix & 1 & - \frac{i}{x} & \cdots \\ - ix^{7} & - x^{6} & ix^{5} & x^{4} & - ix^{3} & - x^{2} & ix & 1 & \cdots \\ \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \ddots \end{pmatrix}. \end{array}
(1.6.135)

A closer look at the SS-matrix re­veals the con­ver­gence de­tails of the power series por­tion ∑n=0∞(ix)n\sum_{n = 0}^{\infty}(ix)^{n}, from the func­tion eixe^{ix}, having a radius of con­ver­gence |x|<1|x| < 1, for which we define ss, as

s=∑n=0∞(ix)n=ii+x,   (for |x|<1),\begin{array}{r} s = \sum_{n = 0}^{\infty}(ix)^{n} = \frac{i}{i + x},\ \ \ \left( for\ |x| < 1 \right), \end{array}
(1.6.136)

and the limit, as x→1x \rightarrow 1, as

limx→1⁡s(x) =1+i2.\begin{array}{r} \lim_{x \rightarrow 1}{s(x)}\ = \frac{1 + i}{2}. \end{array}
(1.6.137)

In the limit limx→1⁡s(x)\lim_{x \rightarrow 1}{s(x)}, we find a result which rep­re­sents the per­fect bal­ance of real and imag­i­nary num­bers, fit­ting the de­f­i­n­i­tion of a cone number from Chap­ter 4. The radius of con­ver­gence for values of |x|<1|x| < 1, holds for the columns formed by en­tries in the lower-left tri­an­gle of the SS-matrix. The main di­ag­o­nal, com­posed of ones, rep­re­sents the crit­i­cal point of the com­plex unit circle's bound­ary, which leads to the con­ver­gence radius tran­si­tion for the func­tion eixe^{ix}, rep­re­sent­ed by the el­e­ments in the upper-right tri­an­gle, hold­ing for |x|>1|x| > 1.

The notion of an un­ex­tend­ed state and the tran­si­tion to an ex­tend­ed one, are math­e­mat­i­cal­ly present in this outer prod­uct matrix ex­am­ple, rem­i­nis­cent of con­cepts re­lat­ed to the in­fi­nite sphere of tan­gents spaces. This notion par­al­lels the log­a­rithm func­tion ln⁡(x)\ln(x), from Figure 1.6.5, and the real number sign tran­si­tion oc­cur­ring at ln⁡(1)=0\ln(1) = 0, the point where the func­tion's curve breach­es the hor­i­zon­tal axis.

Math­e­mat­i­cal­ly speak­ing, a sphere is the bound­ary rep­re­sen­ta­tion of a ball, where the n-ball is de­fined by its nn number of di­men­sions and bound­ed by an (n-1) sphere. The 22-ball, de­pict­ed below in Figure 1.6.9, has for its bound­ary the 11-sphere ∈ℂ\in \mathbb{C}, rep­re­sent­ing the U(1)U(1) group.

Figure
Figure 1.6.9 – The n-ball and its bound­ary the (n-1)-sphere, Sn−1S^{n - 1}.

We can vi­su­al­ize the un­ex­tend­ed state, the ex­tend­ed state, and their phase tran­si­tion by plot­ting the values xj={0.9, 1.0, 1.1}x_{j} = \{ 0.9,\ 1.0,\ 1.1\}, to rep­re­sent the nu­mer­i­cal phase tran­si­tion in the neigh­bor­hood of the crit­i­cal point at ln⁡(x)=0\ln(x) = 0. By taking suc­ces­sive powers of the outer prod­uct matrix X= uˆ⊗vˆX = \ \widehat{u} \otimes \widehat{v}, for n=60n = 60 it­er­a­tions, we are cal­cu­lat­ing the outer prod­uct equiv­a­lent to taking suc­ces­sive powers of (eix)n =eixn\left( e^{ix} \right)^{n}\ = e^{ixn}. This has the effect of ad­vanc­ing the point on the com­plex unit circle by xjx_{j} ra­di­ans for each power it­er­a­tion, for which, the di­ag­o­nal en­tries will remain on the 2-ball's bound­ary.

The un­ex­tend­ed phase plot for xj=0.9x_{j} = 0.9, is seen in Figure 1.6.10, where the di­a­mond points on the unit circle rep­re­sent the el­e­ments of the trace(X)=eixtrace(X) = e^{ix}. Lines con­nect each point to dis­cern the path that el­e­ment takes during each ro­ta­tion. The points are either on the bound­ary of the disc or inside the disc, given that |xj|<1\left| x_{j} \right| < 1, and is topo­log­i­cal­ly de­fined as closed.

The phase tran­si­tion occurs at the crit­i­cal point, when xj=1.0x_{j} = 1.0 (Figure 1.6.11), which rep­re­sents the bound­ary be­tween non-ex­ten­sion and ex­ten­sion, where the tran­si­tion is ev­i­denced by the mor­pho­log­ic change of pat­tern in the closed topol­o­gy, to that of re­peat­ing ring shapes em­a­nat­ing from the center of the 2-ball, ex­tend­ing to its bound­ary.

Figure
Figure 1.6.10 – Outer prod­uct: X= vˆ⊗uˆX = \ \widehat{v} \otimes \widehat{u} , for xj=0.9.x_{j} = 0.9.
Figure
Figure 1.6.11 – Outer prod­uct: X= uˆ⊗vˆX = \ \widehat{u} \otimes \widehat{v} X,X,for x=1.0.x = 1.0.

This math­e­mat­i­cal phase tran­si­tion can be con­sid­ered as the min­i­mal an­gu­lar fre­quen­cy nec­es­sary to extend beyond the bound­ary of the 2-ball. The ex­tend­ed-phase for x=1.1x = 1.1, is plot­ted in Figure 1.6.12, where ray-like spin­dles extend past the ring of di­a­monds rep­re­sent­ing the 2-ball's bound­ary.

Figure
Figure 1.6.12 – Outer prod­uct: X= uˆ⊗vˆX = \ \widehat{u} \otimes \widehat{v} X,X,for x=1.1.x = 1.1.

We can take the outer prod­uct a step fur­ther, by ap­ply­ing it to the matrix groups we have been study­ing, and in par­tic­u­lar, the outer prod­ucts of the in­fi­nite series rep­re­sen­ta­tions for non-com­mu­ta­tive al­ge­bras. We take eXe^{X}, to be the matrix ex­po­nen­tial of the X1X_{1}-type Lie al­ge­bra, and eYe^{Y}, the matrix ex­po­nen­tial of the X2X_{2}-type, for eX=θ1(3,X1)=e(iθ(3,X1)),e^{X} = \theta_{1}^{\left( 3,X_{1} \right)} = e^{\left( i\theta_{\left( 3,X_{1} \right)} \right)},and eY=θ1(3,X2)=e(iθ(3,X2))e^{Y} = \theta_{1}^{\left( 3,X_{2} \right)} = e^{\left( i\theta_{\left( 3,X_{2} \right)} \right)}. We can then present the groups eX,eYe^{X},e^{Y}, as the fol­low­ing in­fi­nite matrix ex­po­nen­tial series

eX=exp⁡(eiθX1(1)0000eiθX1(2)0eiθX1(3)0)=∑n=0∞X00!+X11!+X22!+X33!+…,\begin{array}{r} e^{X} = \exp\begin{pmatrix} e^{i\theta_{X_{1}(1)}} & 0 & 0 \\ 0 & 0 & e^{i\theta_{X_{1}(2)}} \\ 0 & e^{i\theta_{X_{1}(3)}} & 0 \end{pmatrix} = \sum_{n = 0}^{\infty}{\frac{X^{0}}{0!} + \frac{X^{1}}{1!} + \frac{X^{2}}{2!} + \frac{X^{3}}{3!} + \ldots,} \end{array}
(1.6.138)
eY=exp⁡(00eiθX2(1)0eiθX2(2)0eiθX2(3)00)=∑n=0∞Y00!+Y11!+Y22!+Y33!+….\begin{array}{r} e^{Y} = \exp\begin{pmatrix} 0 & 0 & e^{i\theta_{X_{2}(1)}} \\ 0 & e^{i\theta_{X_{2}(2)}} & 0 \\ e^{i\theta_{X_{2}(3)}} & 0 & 0 \end{pmatrix} = \sum_{n = 0}^{\infty}{\frac{Y^{0}}{0!} + \frac{Y^{1}}{1!} + \frac{Y^{2}}{2!} + \frac{Y^{3}}{3!} + \ldots.} \end{array}
(1.6.139)

We know that the prod­ucts XnYk≠YkXnX^{n}Y^{k} \neq Y^{k}X^{n}, for n×k=oddn \times k = odd, will not com­mute, where­as, XnYk=YkXnX^{n}Y^{k} = Y^{k}X^{n} will com­mute for n×k=evenn \times k = even, since even powers of the ma­tri­ces X,YX,Y will only have en­tries on the main di­ag­o­nal, en­sur­ing prod­uct com­mu­ta­tiv­i­ty. We form the ZZ-matrix by the outer prod­uct of the vec­tors Xˆ,Yˆ\widehat{X},\widehat{Y}, which rep­re­sent the in­fi­nite series for eX,eYe^{X},e^{Y}, as

Z=Xˆ⊗Yˆ,\begin{array}{r} Z = \widehat{X} \otimes \widehat{Y}, \end{array}
(1.6.140)

and the matrix group prod­uct eXeY=eZe^{X}e^{Y} = e^{Z}, as the grand sum of the outer prod­uct el­e­ments ∑Xˆ⊗Yˆ\sum_{}^{}{\widehat{X} \otimes \widehat{Y}}.

Main­tain­ing the matrix prod­uct arrange­ment is im­por­tant for XY,XY,as XX and YY do not nec­es­sar­i­ly com­mute, as seen for the ZZ-matrix in Figure 1.6.13. Our at­ten­tion should be drawn to the trace of ZZ, tr(Z)tr(Z), which rep­re­sents the fol­low­ing in­fi­nite series

tr(Z)=∑n=0∞XnYnΓ(n+1)2,\begin{array}{r} tr(Z) = \sum_{n = 0}^{\infty}\frac{X^{n}Y^{n}}{\Gamma(n + 1)^{2}}, \end{array}
(1.6.141)

and if X,YX,Y, com­mute, XY=YXXY = YX, the series is a so­lu­tion of the mod­i­fied Bessel func­tion of the first kind

∑n=0∞XnYnΓ(n+1)2=I0(2XY)=1π∫0πe2XYcosθdθ.\begin{array}{r} \sum_{n = 0}^{\infty}{\frac{X^{n}Y^{n}}{\Gamma(n + 1)^{2}} = I_{0}\left( 2\sqrt{XY} \right) = \frac{1}{\pi}\int_{0}^{\pi}{e^{2\sqrt{XY}cos\theta}d\theta.}} \end{array}
(1.6.142)
(𝐗𝟎𝐘𝟎X0Y1X0Y22!X0Y33!X0Y44!X0Y55!X0Y66!…X0Y∞−1(∞−1)!X0Y∞∞!X1Y0𝐗𝟏𝐘𝟏X1Y22!X1Y33!X1Y44!X1Y55!X1Y66!…X1Y∞−1(∞−1)!X1Y∞∞!X2Y02!X2Y12!𝐗𝟐𝐘𝟐(𝟐!)(𝟐!)X2Y3(3!)(2!)X2Y4(4!)(2!)X2Y5(5!)(2!)X2Y6(6!)(2!)…X2Y∞−1(∞−1)!2!X2Y∞(∞!)(2!)X3Y03!X3Y13!X3Y2(2!)(3!)𝐗𝟑𝐘𝟑(𝟑!)(𝟑!)X3Y4(4!)(3!)X3Y5(5!)(3!)X3Y6(6!)(3!)…X3Y∞−1(∞−1)!3!X3Y∞(∞!)(3!)X4Y04!X4Y14!X4Y2(2!)(4!)X4Y3(3!)(4!)𝐗𝟒𝐘𝟒(𝟒!)(𝟒!)X4Y5(5!)(4!)X4Y6(6!)(4!)…X4Y∞−1(∞−1)!4!X4Y∞(∞!)(4!)X5Y05!X5Y15!X5Y2(2!)(5!)X5Y3(3!)(5!)X5Y4(4!)(5!)𝐗𝟓𝐘𝟓(𝟓!)(𝟓!)X5Y6(6!)(5!)…X5Y∞−1(∞−1)!5!X5Y∞(∞!)(5!)X6Y06!X6Y16!X6Y2(2!)(6!)X6Y3(3!)(6!)X6Y4(4!)(6!)X6Y5(5!)(6!)𝐗𝟔𝐘𝟔(𝟔!)(𝟔!)…X6Y∞−1(∞−1)!6!X6Y∞(∞!)(6!)⋮⋮⋮⋮⋮⋮⋮⋱⋮⋮X∞−1Y0(∞−1)!X∞−1Y1(∞−1)!X∞−1Y22!(∞−1)!X∞−1Y33!(∞−1)!X∞−1Y44!(∞−1)!X∞−1Y55!(∞−1)!X∞−1Y66!(∞−1)!…𝐗∞−𝟏𝐘∞−𝟏(∞−𝟏)!(∞−𝟏)!X∞−1Y∞∞!(∞−1)!X∞Y0∞!X∞Y1∞!X∞Y2(2!)(∞!)X∞Y3(3!)(∞!)X∞Y4(4!)(∞!)X∞Y5(5!)(∞!)X∞Y6(6!)(∞!)…X∞Y∞−1(∞−1)!∞!𝐗∞𝐘∞(∞!)(∞!))\left( \begin{array}{r} \begin{matrix} \mathbf{X}^{\mathbf{0}}\mathbf{Y}^{\mathbf{0}} & X^{0}Y^{1} & \frac{X^{0}Y^{2}}{2!} & \frac{X^{0}Y^{3}}{3!} & \frac{X^{0}Y^{4}}{4!} & \frac{X^{0}Y^{5}}{5!} & \frac{X^{0}Y^{6}}{6!} & \ldots & \frac{X^{0}Y^{\infty - 1}}{(\infty - 1)!} & \frac{X^{0}Y^{\infty}}{\infty!} \\ X^{1}Y^{0} & \mathbf{X}^{\mathbf{1}}\mathbf{Y}^{\mathbf{1}} & \frac{X^{1}Y^{2}}{2!} & \frac{X^{1}Y^{3}}{3!} & \frac{X^{1}Y^{4}}{4!} & \frac{X^{1}Y^{5}}{5!} & \frac{X^{1}Y^{6}}{6!} & \ldots & \frac{X^{1}Y^{\infty - 1}}{(\infty - 1)!} & \frac{X^{1}Y^{\infty}}{\infty!} \\ \frac{X^{2}Y^{0}}{2!} & \frac{{X^{2}Y}^{1}}{2!} & \frac{{\mathbf{X}^{\mathbf{2}}\mathbf{Y}}^{\mathbf{2}}}{(\mathbf{2}!)(\mathbf{2}!)} & \frac{{X^{2}Y}^{3}}{(3!)(2!)} & \frac{{X^{2}Y}^{4}}{(4!)(2!)} & \frac{{X^{2}Y}^{5}}{(5!)(2!)} & \frac{{X^{2}Y}^{6}}{(6!)(2!)} & \ldots & \frac{X^{2}Y^{\infty - 1}}{(\infty - 1)!2!} & \frac{X^{2}Y^{\infty}}{(\infty!)(2!)} \\ \frac{X^{3}Y^{0}}{3!} & \frac{{X^{3}Y}^{1}}{3!} & \frac{{X^{3}Y}^{2}}{(2!)(3!)} & \frac{{\mathbf{X}^{\mathbf{3}}\mathbf{Y}}^{\mathbf{3}}}{(\mathbf{3}!)(\mathbf{3}!)} & \frac{{X^{3}Y}^{4}}{(4!)(3!)} & \frac{{X^{3}Y}^{5}}{(5!)(3!)} & \frac{{X^{3}Y}^{6}}{(6!)(3!)} & \ldots & \frac{X^{3}Y^{\infty - 1}}{(\infty - 1)!3!} & \frac{X^{3}Y^{\infty}}{(\infty!)(3!)} \\ \frac{X^{4}Y^{0}}{4!} & \frac{{X^{4}Y}^{1}}{4!} & \frac{{X^{4}Y}^{2}}{(2!)(4!)} & \frac{X^{4}Y^{3}}{(3!)(4!)} & \frac{{\mathbf{X}^{\mathbf{4}}\mathbf{Y}}^{\mathbf{4}}}{(\mathbf{4}!)(\mathbf{4}!)} & \frac{{X^{4}Y}^{5}}{(5!)(4!)} & \frac{{X^{4}Y}^{6}}{(6!)(4!)} & \ldots & \frac{X^{4}Y^{\infty - 1}}{(\infty - 1)!4!} & \frac{X^{4}Y^{\infty}}{(\infty!)(4!)} \\ \frac{X^{5}Y^{0}}{5!} & \frac{{X^{5}Y}^{1}}{5!} & \frac{{X^{5}Y}^{2}}{(2!)(5!)} & \frac{{X^{5}Y}^{3}}{(3!)(5!)} & \frac{{X^{5}Y}^{4}}{(4!)(5!)} & \frac{{\mathbf{X}^{\mathbf{5}}\mathbf{Y}}^{\mathbf{5}}}{(\mathbf{5}!)(\mathbf{5}!)} & \frac{{X^{5}Y}^{6}}{(6!)(5!)} & \ldots & \frac{X^{5}Y^{\infty - 1}}{(\infty - 1)!5!} & \frac{X^{5}Y^{\infty}}{(\infty!)(5!)} \\ \frac{X^{6}Y^{0}}{6!} & \frac{{X^{6}Y}^{1}}{6!} & \frac{{X^{6}Y}^{2}}{(2!)(6!)} & \frac{{X^{6}Y}^{3}}{(3!)(6!)} & \frac{{X^{6}Y}^{4}}{(4!)(6!)} & \frac{{X^{6}Y}^{5}}{(5!)(6!)} & \frac{{\mathbf{X}^{\mathbf{6}}\mathbf{Y}}^{\mathbf{6}}}{(\mathbf{6}!)(\mathbf{6}!)} & \ldots & \frac{X^{6}Y^{\infty - 1}}{(\infty - 1)!6!} & \frac{X^{6}Y^{\infty}}{(\infty!)(6!)} \\ \mathbf{\vdots} & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \mathbf{\ddots} & \vdots & \vdots \\ \frac{X^{\infty - 1}Y^{0}}{(\infty - 1)!} & \frac{X^{\infty - 1}Y^{1}}{(\infty - 1)!} & \frac{X^{\infty - 1}Y^{2}}{2!(\infty - 1)!} & \frac{X^{\infty - 1}Y^{3}}{3!(\infty - 1)!} & \frac{X^{\infty - 1}Y^{4}}{4!(\infty - 1)!} & \frac{X^{\infty - 1}Y^{5}}{5!(\infty - 1)!} & \frac{X^{\infty - 1}Y^{6}}{6!(\infty - 1)!} & \ldots & \frac{\mathbf{X}^{\mathbf{\infty - 1}}\mathbf{Y}^{\mathbf{\infty - 1}}}{\left( \mathbf{\infty - 1} \right)\mathbf{!}\left( \mathbf{\infty - 1} \right)\mathbf{!}} & \frac{X^{\infty - 1}Y^{\infty}}{\infty!(\infty - 1)!} \\ \frac{X^{\infty}Y^{0}}{\infty!} & \frac{X^{\infty}Y^{1}}{\infty!} & \frac{X^{\infty}Y^{2}}{(2!)(\infty!)} & \frac{X^{\infty}Y^{3}}{(3!)(\infty!)} & \frac{X^{\infty}Y^{4}}{(4!)(\infty!)} & \frac{X^{\infty}Y^{5}}{(5!)(\infty!)} & \frac{X^{\infty}Y^{6}}{(6!)(\infty!)} & \ldots & \frac{X^{\infty}Y^{\infty - 1}}{(\infty - 1)!\infty!} & \frac{\mathbf{X}^{\mathbf{\infty}}\mathbf{Y}^{\mathbf{\infty}}}{\left( \mathbf{\infty!} \right)\left( \mathbf{\infty!} \right)} \end{matrix} \end{array} \right)
Figure 1.6.13 – Outer Prod­uct Matrix, Z=Xˆ⊗YˆZ = \widehat{X} \otimes \widehat{Y}.

The mod­i­fied Bessel func­tion of the first kind In(x)I_{n}(x), for n=0n = 0, is the so­lu­tion of the fol­low­ing mod­i­fied Bessel dif­fer­en­tial equa­tion

x2d2ydx2+xdydx−x2y=0,\begin{array}{r} x^{2}\frac{d^{2}y}{dx^{2}} + x\frac{dy}{dx} - x^{2}y = 0, \end{array}
(1.6.143)

and these func­tions rep­re­sent the radial com­po­nents of the vi­bra­tion modes on the cir­cu­lar drum. To vi­su­al­ize the Bessel func­tion, we take the outer prod­uct Xˆ⊗Yˆ\widehat{X} \otimes \widehat{Y}, for two U(1)U(1)Lie groups with vari­ables X=ix,Y=iy,X = ix,Y = iy,to yield a Bessel func­tion of the first kind, as

∑n=0∞(ix)n(iy)nΓ(n+1)2=J0(2xy).\begin{array}{r} \sum_{n = 0}^{\infty}{\frac{(ix)^{n}(iy)^{n}}{\Gamma(n + 1)^{2}} = J_{0}\left( 2\sqrt{xy} \right).} \end{array}
(1.6.144)

In Figure 1.6.14 we plot this Bessel func­tion for both its real and imag­i­nary radial com­po­nents, where wave forms can be seen prop­a­gat­ing along the sur­face of the drum and dis­pers­ing as the dis­tance from the im­pulse in­creas­es.

Figure
Figure 1.6.14 – Radial com­po­nents of the first-order Bessel func­tion, J0(2xy)J_{0}\left( 2\sqrt{xy} \right).

The ZZ-matrix can also be mul­ti­plied on the left by the in­verse di­ag­o­nal com­po­nents of X,X, (DX×Z){(D}_{X} \times Z), and mul­ti­plied on the right by the in­verse di­ag­o­nal com­po­nents of Y,(Z×DY)Y,{(Z \times D}_{Y}), for DX,DYD_{X},D_{Y} de­fined, as

DX=diag(X0, X1,X22!,X33!,X44!,X55!,X66!,…),D_{X} = diag\left( X^{0},\ X^{1},\frac{X^{2}}{2!},\frac{X^{3}}{3!},\frac{X^{4}}{4!},\frac{X^{5}}{5!},\frac{X^{6}}{6!},\ldots \right),
DY=diag(Y0, Y1,Y22!,Y33!,Y44!,Y55!,Y66!,…),\begin{array}{r} D_{Y} = diag\left( Y^{0},\ Y^{1},\frac{Y^{2}}{2!},\frac{Y^{3}}{3!},\frac{Y^{4}}{4!},\frac{Y^{5}}{5!},\frac{Y^{6}}{6!},\ldots \right), \end{array}
(1.6.145)

to find the fol­low­ing two-vari­able tan­gent space field matrix SXYS_{XY}, com­posed of X,YX,Y

SXY=DY−1ZDX−1=(1X1Y−1X2Y−2X3Y−3X4Y−4X5Y−5X6Y−6…X−1Y11X1Y−1X2Y−2X3Y−3X4Y−4X5Y−5…X−2Y2X−1Y11X1Y−1X2Y−2X3Y−3X4Y−4…X−3Y3X−2Y2X−1Y11X1Y−1X2Y−2X3Y−3…X−4Y4X−3Y3X−2Y2X−1Y11X1Y−1X2Y−2…X−5Y5X−4Y4X−3Y3X−2Y2X−1Y11X1Y−1…X−6Y6X−5Y5X−4Y4X−3Y3X−2Y2X−1Y11…⋮⋮⋮⋮⋮⋮⋮⋱).\begin{array}{r} {S_{XY} = D}_{Y}^{- 1}ZD_{X}^{- 1} = \begin{pmatrix} 1 & X^{1}Y^{- 1} & X^{2}Y^{- 2} & X^{3}Y^{- 3} & X^{4}Y^{- 4} & X^{5}Y^{- 5} & X^{6}Y^{- 6} & \ldots \\ X^{- 1}Y^{1} & 1 & X^{1}Y^{- 1} & X^{2}Y^{- 2} & X^{3}Y^{- 3} & X^{4}Y^{- 4} & X^{5}Y^{- 5} & \ldots \\ X^{- 2}Y^{2} & X^{- 1}Y^{1} & 1 & X^{1}Y^{- 1} & X^{2}Y^{- 2} & X^{3}Y^{- 3} & X^{4}Y^{- 4} & \ldots \\ {X^{- 3}Y}^{3} & X^{- 2}Y^{2} & X^{- 1}Y^{1} & 1 & X^{1}Y^{- 1} & X^{2}Y^{- 2} & X^{3}Y^{- 3} & \ldots \\ {X^{- 4}Y}^{4} & {X^{- 3}Y}^{3} & X^{- 2}Y^{2} & X^{- 1}Y^{1} & 1 & X^{1}Y^{- 1} & X^{2}Y^{- 2} & \ldots \\ {X^{- 5}Y}^{5} & {X^{- 4}Y}^{4} & {X^{- 3}Y}^{3} & X^{- 2}Y^{2} & X^{- 1}Y^{1} & 1 & X^{1}Y^{- 1} & \ldots \\ {X^{- 6}Y}^{6} & {X^{- 5}Y}^{5} & {X^{- 4}Y}^{4} & {X^{- 3}Y}^{3} & X^{- 2}Y^{2} & X^{- 1}Y^{1} & 1 & \ldots \\ \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \ddots \end{pmatrix}. \end{array}
(1.6.146)

and the matrix SXYS_{XY}, is like the outer-prod­uct matrix SS in Eq. 1.6.135, having en­tries with in­vers­es for their trans­pos­es. We define the fol­low­ing non-com­mut­ing ma­tri­ces X,YX,Y

X=(rX1(1)eiθX1(1)0000rX1(2)eiθX1(2)0rX1(3)eiθX1(3)0),   Y=(00rX2(1)eiθX2(1)0rX2(2)eiθX2(2)0rX2(3)eiθX2(3)00),\begin{array}{r} X = \begin{pmatrix} {r_{X_{1}(1)}e}^{i\theta_{X_{1}(1)}} & 0 & 0 \\ 0 & 0 & {r_{X_{1}(2)}e}^{i\theta_{X_{1}(2)}} \\ 0 & r_{X_{1}(3)}e^{i\theta_{X_{1}(3)}} & 0 \end{pmatrix},\ \ \ Y = \begin{pmatrix} 0 & 0 & r_{X_{2}(1)}e^{i\theta_{X_{2}(1)}} \\ 0 & r_{X_{2}(2)}e^{i\theta_{X_{2}(2)}} & 0 \\ r_{X_{2}(3)}e^{i\theta_{X_{2}(3)}} & 0 & 0 \end{pmatrix}, \end{array}
(1.6.147)

for X=θr(3,X1)=rX1eiθ(3,X1)X = \theta_{r}^{\left( 3,X_{1} \right)} = r_{X_{1}}e^{i\theta_{\left( 3,X_{1} \right)}}, and Y=θr(3,X2)=rX2eiθ(3,X2)Y = \theta_{r}^{\left( 3,X_{2} \right)} = {r_{X_{2}}e}^{i\theta_{\left( 3,X_{2} \right)}}, The fol­low­ing prod­ucts would result:

XX∗,YY∗=I3,XX^{*},YY^{*} = I_{3},
X−nYnXnY−n→Y2, (for n, even),X^{- n}Y^{n}X^{n}Y^{- n} \rightarrow Y_{2},\ (for\ n,\ even),
X−nYnXnY−n→Y3,(for n, odd),\begin{array}{r} X^{- n}Y^{n}X^{n}Y^{- n} \rightarrow Y_{3},(for\ n,\ odd), \end{array}
(1.6.148)

for I3I_{3}, the 3×33 \times 3 iden­ti­ty matrix, and Y2,Y3Y_{2},Y_{3}, the YY-type cyclic matrix groups. A prod­uct ex­am­ple in­volv­ing X,YX,Y, is the Lie al­ge­bra prod­uct X−3Y3X3Y−3{X^{- 3}Y}^{3}X^{3}Y^{- 3}, which will gen­er­ate the fol­low­ing Y2Y_{2}-type matrix

X−3Y3X3Y−3=(0rX1(2)(rX2(1))2(rX2(1))2rX2(3)(rX1(1))3(rX2(2))3ei(−3(θX1(1)+θX2(2))+2(θX1(3)+θX2(1))+(θX1(2)+θX2(3)))000rX2(1)(rX2(3))2rX1(2)(rX1(3))2ei(3θX1(1)−2θX1(3)−θX1(2))(rX2(2))3(rX1(2))2rX1(3)ei(3θX2(2)−2θX2(1)−θX2(3))00),\begin{array}{r} {X^{- 3}Y}^{3}X^{3}Y^{- 3} = \begin{pmatrix} 0 & \frac{{r_{X_{1}(2)}\left( r_{X_{2}(1)} \right)^{2}\left( r_{X_{2}(1)} \right)}^{2}r_{X_{2}(3)}}{\left( r_{X_{1}(1)} \right)^{3}\left( r_{X_{2}(2)} \right)^{3}}e^{i\left( \begin{array}{r} - 3\left( {\theta_{X_{1}(1)} + \theta}_{X_{2}(2)} \right) \\ + 2\left( {\theta_{X_{1}(3)} + \theta}_{X_{2}(1)} \right) \\ + \left( {\theta_{X_{1}(2)} + \theta}_{X_{2}(3)} \right) \end{array} \right)} & 0 \\ 0 & 0 & \frac{{r_{X_{2}(1)}\left( r_{X_{2}(3)} \right)}^{2}}{r_{X_{1}(2)}\left( r_{X_{1}(3)} \right)^{2}}e^{i\left( \begin{array}{r} 3\theta_{X_{1}(1)} \\ - 2\theta_{X_{1}(3)} \\ {- \theta}_{X_{1}(2)} \end{array} \right)} \\ \frac{\left( r_{X_{2}(2)} \right)^{3}}{\left( r_{X_{1}(2)} \right)^{2}r_{X_{1}(3)}}e^{i\left( \begin{array}{r} 3\theta_{X_{2}(2)} \\ - 2\theta_{X_{2}(1)} \\ {- \theta}_{X_{2}(3)} \end{array} \right)} & 0 & 0 \end{pmatrix}, \end{array}
(1.6.149)

given nn, odd-valued, and the X−nYnXnY−nX^{- n}Y^{n}X^{n}Y^{- n} prod­uct arrange­ment. The evenly dis­trib­uted rr terms could have also been omit­ted in the matrix result, since they are evenly dis­trib­uted in the nu­mer­a­tor and de­nom­i­na­tors and cancel in each ra­tio­nal. Each entry in the field SXYS_{XY}, is com­prised of tan­gent space el­e­ments for their ar­gu­ments, which can be ap­pre­ci­at­ed by eval­u­at­ing the prod­uct for the SXY(2,1)S_{XY(2,1)} entry, X−1Y1X^{- 1}Y^{1},

X−1Y1=(rX2(3)rX1(3)00rX2(1)rX1(1)(ei(θX2(1)−θX1(1)))(ei(θX2(3)−θX1(3)))000rX2(2)rX1(2)(ei(θX2(2)−θX1(2)))0),\begin{array}{r} X^{- 1}Y^{1} = \left( \frac{r_{X_{2}(3)}}{r_{X_{1}(3)}}\begin{matrix} 0 & 0 & \frac{r_{X_{2}(1)}}{r_{X_{1}(1)}}\left( e^{i\left( {\theta_{X_{2}(1)} - \theta}_{X_{1}(1)} \right)} \right) \\ \left( e^{i\left( {\theta_{X_{2}(3)} - \theta}_{X_{1}(3)} \right)} \right) & 0 & 0 \\ 0 & \frac{r_{X_{2}(2)}}{r_{X_{1}(2)}}\left( e^{i\left( {\theta_{X_{2}(2)} - \theta}_{X_{1}(2)} \right)} \right) & 0 \end{matrix} \right), \end{array}
(1.6.150)

where the in­verse al­ge­bra el­e­ments X−1X^{- 1} rotate in the op­po­site di­rec­tion to the Y1Y^{1} al­ge­bra. So, the matrix X−1Y1X^{- 1}Y^{1}, can be thought of as a su­per­po­si­tion of the Lie al­ge­bras ∈X−1,Y1\in X^{- 1},Y^{1}, in the ex­po­nen­tial, mul­ti­plied by ratios of the al­ge­bra vari­ables (rX2(j)rX1(j))\left( \frac{r_{X_{2}(j)}}{r_{X_{1}(j)}} \right), as lead­ing co­ef­fi­cients. Math­e­mat­i­cal­ly we can con­cep­tu­al­ize the outer prod­uct eX⊗eYe^{X} \otimes e^{Y}, as a field formed, as SXYS_{XY}, in which the Lie groups eX, eYe^{X},\ e^{Y} op­er­ate in.

Con­cep­tu­al­ly, we can see move­ment in the SXYS_{XY} field, by taking it­er­a­tion prod­ucts for XX, on the left, and YY, on the right, to map or trans­form from a start­ing po­si­tion to an end po­si­tion. For ex­am­ple, the fol­low­ing se­quence of prod­ucts takes us from a start­ing point at 𝐗−𝟏𝐘𝟏\mathbf{X}^{- \mathbf{1}}\mathbf{Y}^{\mathbf{1}}, to the end po­si­tion at 𝐗𝟐𝐘−𝟐\mathbf{X}^{\mathbf{2}}\mathbf{Y}^{- \mathbf{2}}:

𝐗−𝟏𝐘𝟏→(X1X−1Y1Y−1=𝟏)→(X1X1X−1Y1Y−1Y−1=𝐗𝟏𝐘−𝟏)→(X1X1X1X−1Y1Y−1Y−1Y−1=𝐗𝟐𝐘−𝟐),\begin{array}{r} \mathbf{X}^{- \mathbf{1}}\mathbf{Y}^{\mathbf{1}} \rightarrow {(X}^{1}X^{- 1}Y^{1}Y^{- 1} = \mathbf{1}) \rightarrow \left( {X^{1}X}^{1}X^{- 1}Y^{1}Y^{- 1}Y^{- 1} = \mathbf{X}^{\mathbf{1}}\mathbf{Y}^{- \mathbf{1}} \right) \\ \rightarrow \left( {X^{1}X^{1}X}^{1}X^{- 1}Y^{1}Y^{- 1}Y^{- 1}Y^{- 1} = \mathbf{X}^{\mathbf{2}}\mathbf{Y}^{- \mathbf{2}} \right), \end{array}
(1.6.151)

This move­ment se­quence is shown in bold with di­rec­tion­al arrows, rep­re­sent­ing the path-con­nec­tion

(1X1Y−1X2Y−2X3Y−3X4Y−4X5Y−5X6Y−6…𝐗−𝟏𝐘𝟏→𝟏→𝐗𝟏𝐘−𝟏→→𝐗𝟐𝐘−𝟐X3Y−3X4Y−4X5Y−5…X−2Y2X−1Y11X1Y−1→X2Y−2X3Y−3X4Y−4…X−3Y3X−2Y2X−1Y11X1Y−1X2Y−2X3Y−3…X−4Y4X−3Y3X−2Y2X−1Y11X1Y−1X2Y−2…X−5Y5X−4Y4X−3Y3X−2Y2X−1Y11X1Y−1…X−6Y6X−5Y5X−4Y4X−3Y3X−2Y2X−1Y11…⋮⋮⋮⋮⋮⋮⋮⋱).\begin{array}{r} \begin{pmatrix} 1 & X^{1}Y^{- 1} & X^{2}Y^{- 2} & X^{3}Y^{- 3} & X^{4}Y^{- 4} & X^{5}Y^{- 5} & X^{6}Y^{- 6} & \ldots \\ \mathbf{X}^{- \mathbf{1}}\mathbf{Y}^{\mathbf{1}}\mathbf{\rightarrow} & \mathbf{1} \rightarrow & \mathbf{X}^{\mathbf{1}}\mathbf{Y}^{- \mathbf{1}} \rightarrow & \rightarrow \mathbf{X}^{\mathbf{2}}\mathbf{Y}^{- \mathbf{2}} & X^{3}Y^{- 3} & X^{4}Y^{- 4} & X^{5}Y^{- 5} & \ldots \\ X^{- 2}Y^{2} & X^{- 1}Y^{1} & 1 & X^{1}Y^{- 1} & \rightarrow X^{2}Y^{- 2} & X^{3}Y^{- 3} & X^{4}Y^{- 4} & \ldots \\ {X^{- 3}Y}^{3} & X^{- 2}Y^{2} & X^{- 1}Y^{1} & 1 & X^{1}Y^{- 1} & X^{2}Y^{- 2} & X^{3}Y^{- 3} & \ldots \\ {X^{- 4}Y}^{4} & {X^{- 3}Y}^{3} & X^{- 2}Y^{2} & X^{- 1}Y^{1} & 1 & X^{1}Y^{- 1} & X^{2}Y^{- 2} & \ldots \\ {X^{- 5}Y}^{5} & {X^{- 4}Y}^{4} & {X^{- 3}Y}^{3} & X^{- 2}Y^{2} & X^{- 1}Y^{1} & 1 & X^{1}Y^{- 1} & \ldots \\ {X^{- 6}Y}^{6} & {X^{- 5}Y}^{5} & {X^{- 4}Y}^{4} & {X^{- 3}Y}^{3} & X^{- 2}Y^{2} & X^{- 1}Y^{1} & 1 & \ldots \\ \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \ddots \end{pmatrix}. \end{array}
(1.6.152)

As a whole, move­ment in the outer prod­uct field Z=Xˆ⊗YˆZ = \widehat{X} \otimes \widehat{Y}, occurs by di­rec­tion­al in­te­gra­tion and dif­fer­en­ti­a­tion, where sim­i­lar to our pre­vi­ous ex­am­ple; from a start­ing point at X−1Y1∈ZX^{- 1}Y^{1} \in Z, moving three field po­si­tions to the right is the result of three in­te­gra­tions of X0Y1X^{0}Y^{1}with re­spect to XX, which, as­sum­ing no con­stants of in­te­gra­tion, is

X0Y1→X3Y13!=∭X0Y1dXdXdX =X3Y13!.\begin{array}{r} X^{0}Y^{1} \rightarrow \frac{X^{3}Y^{1}}{3!} = \iiint_{}^{}{X^{0}Y^{1}}dXdXdX\ = \frac{X^{3}Y^{1}}{3!}. \end{array}
(1.6.153)

Next, by taking the par­tial de­riv­a­tive, ∂Y\partial Y, of this in­te­gral result, the field po­si­tion shifts upward by one

X3Y13!→X3Y03!=ddY(X3Y13!)=X3Y03!,\begin{array}{r} \frac{X^{3}Y^{1}}{3!} \rightarrow \frac{X^{3}Y^{0}}{3!} = \frac{d}{dY}\left( \frac{X^{3}Y^{1}}{3!} \right) = \frac{X^{3}Y^{0}}{3!}, \end{array}
(1.6.154)

and the path-con­nect­ed move­ment se­quence is shown in bold in the ZZ-matrix, as

(X0Y0X1Y0X2Y02!𝐗𝟑𝐘𝟎𝟑!X4Y04!X5Y05!X6Y06!…𝐗𝟎𝐘𝟏→∫dX1𝐗𝟏𝐘𝟏→∫dX2𝐗𝟐𝐘𝟏𝟐!→∫dX3𝐗𝟑𝐘𝟏𝟑! ↑ddYX4Y14!X5Y15!X6Y16!…X0Y22!X1Y22!X2Y2(2!)(2!)X3Y2(3!)(2!)X4Y2(4!)(2!)X5Y2(5!)(2!)X6Y2(6!)(2!)…X0Y33!X1Y33!X2Y3(2!)(3!)X3Y3(3!)(3!)X4Y3(4!)(3!)X5Y3(5!)(3!)X6Y3(6!)(3!)…X0Y44!X1Y44!X2Y4(2!)(4!)X3Y4(3!)(4!)X4Y4(4!)(4!)X5Y4(5!)(4!)X6Y4(6!)(4!)…X0Y55!X1Y55!X2Y5(2!)(5!)X3Y5(3!)(5!)X4Y5(4!)(5!)X5Y5(5!)(5!)X6Y5(6!)(5!)…X0Y66!X1Y66!X2Y6(2!)(6!)X3Y6(3!)(6!)X4Y6(4!)(6!)X5Y6(5!)(6!)X6Y6(6!)(6!)…⋮⋮⋮⋮⋮⋮⋮⋱).\begin{array}{r} \begin{pmatrix} X^{0}Y^{0} & X^{1}Y^{0} & \frac{X^{2}Y^{0}}{2!} & \frac{\mathbf{X}^{\mathbf{3}}\mathbf{Y}^{\mathbf{0}}}{\mathbf{3}!} & \frac{X^{4}Y^{0}}{4!} & \frac{X^{5}Y^{0}}{5!} & \frac{X^{6}Y^{0}}{6!} & \ldots \\ \mathbf{X}^{\mathbf{0}}\mathbf{Y}^{\mathbf{1}}\overset{\int_{}^{}{dX^{1}}}{\rightarrow} & \mathbf{X}^{\mathbf{1}}\mathbf{Y}^{\mathbf{1}}\overset{\int_{}^{}{dX^{2}}}{\rightarrow} & \frac{\mathbf{X}^{\mathbf{2}}\mathbf{Y}^{\mathbf{1}}}{\mathbf{2}!}\overset{\int_{}^{}{dX^{3}}}{\rightarrow} & \frac{\mathbf{X}^{\mathbf{3}}\mathbf{Y}^{\mathbf{1}}}{\mathbf{3}!}\ \mathbf{\uparrow}\frac{d}{dY} & \frac{X^{4}Y^{1}}{4!} & \frac{X^{5}Y^{1}}{5!} & \frac{X^{6}Y^{1}}{6!} & \ldots \\ \frac{{X^{0}Y}^{2}}{2!} & \frac{{X^{1}Y}^{2}}{2!} & \frac{{X^{2}Y}^{2}}{(2!)(2!)} & \frac{{X^{3}Y}^{2}}{(3!)(2!)} & \frac{{X^{4}Y}^{2}}{(4!)(2!)} & \frac{{X^{5}Y}^{2}}{(5!)(2!)} & \frac{{X^{6}Y}^{2}}{(6!)(2!)} & \ldots \\ \frac{{X^{0}Y}^{3}}{3!} & \frac{{X^{1}Y}^{3}}{3!} & \frac{{X^{2}Y}^{3}}{(2!)(3!)} & \frac{{X^{3}Y}^{3}}{(3!)(3!)} & \frac{{X^{4}Y}^{3}}{(4!)(3!)} & \frac{{X^{5}Y}^{3}}{(5!)(3!)} & \frac{{X^{6}Y}^{3}}{(6!)(3!)} & \ldots \\ \frac{X^{0}Y^{4}}{4!} & \frac{{X^{1}Y}^{4}}{4!} & \frac{{X^{2}Y}^{4}}{(2!)(4!)} & \frac{{X^{3}Y}^{4}}{(3!)(4!)} & \frac{{X^{4}Y}^{4}}{(4!)(4!)} & \frac{{X^{5}Y}^{4}}{(5!)(4!)} & \frac{{X^{6}Y}^{4}}{(6!)(4!)} & \ldots \\ \frac{X^{0}Y^{5}}{5!} & \frac{{X^{1}Y}^{5}}{5!} & \frac{{X^{2}Y}^{5}}{(2!)(5!)} & \frac{{X^{3}Y}^{5}}{(3!)(5!)} & \frac{{X^{4}Y}^{5}}{(4!)(5!)} & \frac{{X^{5}Y}^{5}}{(5!)(5!)} & \frac{{X^{6}Y}^{5}}{(6!)(5!)} & \ldots \\ \frac{X^{0}Y^{6}}{6!} & \frac{{X^{1}Y}^{6}}{6!} & \frac{{X^{2}Y}^{6}}{(2!)(6!)} & \frac{{X^{3}Y}^{6}}{(3!)(6!)} & \frac{{X^{4}Y}^{6}}{(4!)(6!)} & \frac{{X^{5}Y}^{6}}{(5!)(6!)} & \frac{{X^{6}Y}^{6}}{(6!)(6!)} & \ldots \\ \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \ddots \end{pmatrix}. \end{array}
(1.6.155)

Math­e­mat­i­cal­ly speak­ing, move­ment in the ZZ-matrix away from the iden­ti­ty at ϵ=X0Y0\epsilon = X^{0}Y^{0} is an­a­lyt­i­cal­ly con­ver­gent, since its basis is rooted in the ab­solute con­ver­gence of the ex­po­nen­tial series. So, the limits as n→∞n \rightarrow \infty, for in­fi­nite in­te­gra­tions to the right (dX)(dX), down­ward (dY)(dY), and di­ag­o­nal­ly (dXdY)(dXdY), all con­verge to zero

limn→∞⁡XnYmn!m!=0,\begin{array}{r} \lim_{n \rightarrow \infty}{\frac{X^{n}Y^{m}}{n!m!} = 0}, \end{array}
(1.6.156)
limm→∞⁡XnYmn!m!=0,\begin{array}{r} \lim_{m \rightarrow \infty}\frac{X^{n}Y^{m}}{n!m!} = 0, \end{array}
(1.6.157)
limn,m→∞⁡XnYmn!m!=0.\begin{array}{r} \lim_{n,m \rightarrow \infty}\frac{X^{n}Y^{m}}{n!m!} = 0. \end{array}
(1.6.158)

This holds, given the sum of el­e­ments for any row, column, or di­ag­o­nal, con­verge to an ex­po­nen­tial map­ping of Lie al­ge­bras X,YX,Y, in varied pro­por­tions. Dif­fer­en­tial move­ment toward the origin is anal­o­gous to its in­te­gra­tion coun­ter­part, where the re­spec­tive limits will all con­verge to the iden­ti­ty ϵ\epsilon, for n,m,(n,m)→0.n,m,(n,m) \rightarrow 0.

We have seen that the Baker-Camp­bell-Haus­dorff for­mu­la ex­press­es the matrix log­a­rithm for the prod­uct of two matrix groups as an in­fi­nite series of nested Lie brack­ets

ln⁡(eXeY)=Y+∫01g[exp⁡(tadX)exp⁡(adY)](X) dt,\begin{array}{r} \ln\left( e^{X}e^{Y} \right) = Y + \int_{0}^{1}{g\left\lbrack \exp{\left( tad_{X} \right)\exp\left( ad_{Y} \right)} \right\rbrack}(X)\ dt, \end{array}
(1.6.159)

and the outer prod­uct of two matrix groups in in­fi­nite series form eX⊗eYe^{X} \otimes e^{Y}, rep­re­sents the result eXeY=eZe^{X}e^{Y} = e^{Z}, as the sum­ma­tion of an in­fi­nite number of com­mu­ta­tive and or non-com­mu­ta­tive Lie al­ge­bras

eZ=∑n=0∞(X+Y)nn!=∑n=0∞(∑k=0nX(k)Y(n−k)k!(n−k)! ).\begin{array}{r} e^{Z} = \sum_{n = 0}^{\infty}\frac{(X + Y)^{n}}{n!} = \sum_{n = 0}^{\infty}\left( \sum_{k = 0}^{n}\frac{X^{(k)}Y^{(n - k)}}{k!(n - k)!\ } \right). \end{array}
(1.6.160)

The sum­ma­tion above can be de­fined by a bi­no­mi­al fac­to­r­i­al series, from the well-known bi­no­mi­al the­o­rem

(X+Y)n=∑k=0nn!(n−k)!k!Xn−kYk=∑k=0n(nk)Xn−kYk,\begin{array}{r} (X + Y)^{n} = \sum_{k = 0}^{n}{\frac{n!}{(n - k)!k!}X^{n - k}Y^{k} = \sum_{k = 0}^{n}{\binom{n}{k}X^{n - k}Y^{k}}}, \end{array}
(1.6.161)

with the ex­cep­tion that the rec­i­p­ro­cal fac­to­r­i­al series is present in the de­nom­i­na­tor, and we must keep in mind that X,YX,Y, will not nec­es­sar­i­ly com­mute as they are typ­i­cal­ly treat­ed in the com­bi­na­torics of com­mut­ing bi­no­mi­al theory. Ap­ply­ing the outer prod­uct ap­proach to the matrix log­a­rithm also en­ables us to study the fields that they form. By choos­ing vec­tors de­fined by the log­a­rithm func­tion's in­fi­nite series

∑n=1∞(−1)n−1(X−1)nn,\begin{array}{r} \sum_{n = 1}^{\infty}\frac{( - 1)^{n - 1}(X - 1)^{n}}{n}, \end{array}
(1.6.162)

we set the vector uˆ\widehat{u} to be the geo­met­ric series in the nu­mer­a­tor, and the vector vˆ\widehat{v} to be the al­ter­nat­ing har­mon­ic series

uˆ=((X−1)1,(X−1)2,(X−1)3,(X−1)4,(X−1)5,(X−1)6,(X−1)7… ),\begin{array}{r} \widehat{u} = \left( (X - 1)^{1},(X - 1)^{2},(X - 1)^{3},(X - 1)^{4},(X - 1)^{5},(X - 1)^{6},(X - 1)^{7}\ldots\ \right), \end{array}
(1.6.163)
vˆ=(1,−12,13,−14,15,−16,17,… ).\begin{array}{r} \widehat{v} = \left( 1, - \frac{1}{2},\frac{1}{3}, - \frac{1}{4},\frac{1}{5}, - \frac{1}{6},\frac{1}{7},\ldots\ \right). \end{array}
(1.6.164)

We take the outer prod­uct of these two vec­tors, as X=uˆ⊗vˆX = \widehat{u} \otimes \widehat{v}, which re­sults in the fol­low­ing XX-matrix

X=((X−1)−(X−1)2(X−1)3−(X−1)4(X−1)5−(X−1)6(X−1)7⋯(X−1)2−(X−1)22(X−1)23−(X−1)24(X−1)25−(X−1)26(X−1)27⋯(X−1)3−(X−1)32(X−1)33−(X−1)34(X−1)35−(X−1)36(X−1)37⋯(X−1)4−(X−1)42(X−1)43−(X−1)44(X−1)45−(X−1)46(X−1)47⋯(X−1)5−(X−1)52(X−1)53−(X−1)54(X−1)55−(X−1)56(X−1)57⋯(X−1)6−(X−1)62(X−1)63−(X−1)64(X−1)65−(X−1)66(X−1)67⋯(X−1)7−(X−1)72(X−1)73−(X−1)74(X−1)75−(X−1)76(X−1)77⋯⋮⋮⋮⋮⋮⋮⋮⋱),\begin{array}{r} X = \begin{pmatrix} (X - 1) & - \frac{(X - 1)}{2} & \frac{(X - 1)}{3} & - \frac{(X - 1)}{4} & \frac{(X - 1)}{5} & - \frac{(X - 1)}{6} & \frac{(X - 1)}{7} & \cdots \\ (X - 1)^{2} & - \frac{(X - 1)^{2}}{2} & \frac{(X - 1)^{2}}{3} & - \frac{(X - 1)^{2}}{4} & \frac{(X - 1)^{2}}{5} & - \frac{(X - 1)^{2}}{6} & \frac{(X - 1)^{2}}{7} & \cdots \\ (X - 1)^{3} & - \frac{(X - 1)^{3}}{2} & \frac{(X - 1)^{3}}{3} & - \frac{(X - 1)^{3}}{4} & \frac{(X - 1)^{3}}{5} & - \frac{(X - 1)^{3}}{6} & \frac{(X - 1)^{3}}{7} & \cdots \\ (X - 1)^{4} & - \frac{(X - 1)^{4}}{2} & \frac{(X - 1)^{4}}{3} & - \frac{(X - 1)^{4}}{4} & \frac{(X - 1)^{4}}{5} & - \frac{(X - 1)^{4}}{6} & \frac{(X - 1)^{4}}{7} & \cdots \\ (X - 1)^{5} & - \frac{(X - 1)^{5}}{2} & \frac{(X - 1)^{5}}{3} & - \frac{(X - 1)^{5}}{4} & \frac{(X - 1)^{5}}{5} & - \frac{(X - 1)^{5}}{6} & \frac{(X - 1)^{5}}{7} & \cdots \\ (X - 1)^{6} & - \frac{(X - 1)^{6}}{2} & \frac{(X - 1)^{6}}{3} & - \frac{(X - 1)^{6}}{4} & \frac{(X - 1)^{6}}{5} & - \frac{(X - 1)^{6}}{6} & \frac{(X - 1)^{6}}{7} & \cdots \\ (X - 1)^{7} & \frac{- (X - 1)^{7}}{2} & \frac{(X - 1)^{7}}{3} & \frac{- (X - 1)^{7}}{4} & \frac{(X - 1)^{7}}{5} & \frac{- (X - 1)^{7}}{6} & \frac{(X - 1)^{7}}{7} & \cdots \\ \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \ddots \end{pmatrix}, \end{array}
(1.6.165)

where anal­o­gous to the ex­po­nen­tial series, the tr(X)tr(X) is the in­fi­nite series that de­fines the log­a­rithm func­tion

tr(X)=∑n=1∞(−1)n−1(X−1)nn=ln⁡(X).\begin{array}{r} tr(X) = \sum_{n = 1}^{\infty}\frac{( - 1)^{n - 1}(X - 1)^{n}}{n} = \ln(X). \end{array}
(1.6.166)

Next, we define DXD_{X}, as a matrix of the el­e­ments along the main di­ag­o­nal of the XX-matrix, as

DX=diag((X−1)1,−(X−1)22,(X−1)33,−(X−1)44,(X−1)55,−(X−1)66,(X−1)77,…),\begin{array}{r} D_{X} = diag\left( (X - 1)^{1}, - \frac{(X - 1)^{2}}{2},\frac{(X - 1)^{3}}{3}, - \frac{(X - 1)^{4}}{4},\frac{(X - 1)^{5}}{5}, - \frac{(X - 1)^{6}}{6},\frac{(X - 1)^{7}}{7},\ldots \right), \end{array}
(1.6.167)

which, when taken as a sum, rep­re­sents the log­a­rithm func­tion series ln⁡(X).\ln(X).Next, we take the prod­uct SX=XDX−1{S_{X} = XD}_{X}^{- 1}, to find the geo­met­ric log­a­rithm field SXS_{X}

SX=XDX−1=(11(X−1)1(X−1)21(X−1)31(X−1)41(X−1)51(X−1)61(X−1)7⋯(X−1)11(X−1)1(X−1)21(X−1)31(X−1)41(X−1)51(X−1)6⋯(X−1)2(X−1)11(X−1)1(X−1)21(X−1)31(X−1)41(X−1)5⋯(X−1)3(X−1)2(X−1)11(X−1)1(X−1)21(X−1)31(X−1)4⋯(X−1)4(X−1)3(X−1)2(X−1)11(X−1)1(X−1)21(X−1)3⋯(X−1)5(X−1)4(X−1)3(X−1)2(X−1)11(X−1)1(X−1)2⋯(X−1)6(X−1)5(X−1)4(X−1)3(X−1)2(X−1)11(X−1)⋯(X−1)7(X−1)6(X−1)5(X−1)4(X−1)3(X−1)2(X−1)1⋯⋮⋮⋮⋮⋮⋮⋮⋮⋱),\begin{array}{r} {S_{X} = XD}_{X}^{- 1} = \begin{pmatrix} 1 & \frac{1}{(X - 1)} & \frac{1}{(X - 1)^{2}} & \frac{1}{(X - 1)^{3}} & \frac{1}{(X - 1)^{4}} & \frac{1}{(X - 1)^{5}} & \frac{1}{(X - 1)^{6}} & \frac{1}{(X - 1)^{7}} & \cdots \\ (X - 1) & 1 & \frac{1}{(X - 1)} & \frac{1}{(X - 1)^{2}} & \frac{1}{(X - 1)^{3}} & \frac{1}{(X - 1)^{4}} & \frac{1}{(X - 1)^{5}} & \frac{1}{(X - 1)^{6}} & \cdots \\ (X - 1)^{2} & (X - 1) & 1 & \frac{1}{(X - 1)} & \frac{1}{(X - 1)^{2}} & \frac{1}{(X - 1)^{3}} & \frac{1}{(X - 1)^{4}} & \frac{1}{(X - 1)^{5}} & \cdots \\ (X - 1)^{3} & (X - 1)^{2} & (X - 1) & 1 & \frac{1}{(X - 1)} & \frac{1}{(X - 1)^{2}} & \frac{1}{(X - 1)^{3}} & \frac{1}{(X - 1)^{4}} & \cdots \\ (X - 1)^{4} & (X - 1)^{3} & (X - 1)^{2} & (X - 1) & 1 & \frac{1}{(X - 1)} & \frac{1}{(X - 1)^{2}} & \frac{1}{(X - 1)^{3}} & \cdots \\ (X - 1)^{5} & (X - 1)^{4} & (X - 1)^{3} & (X - 1)^{2} & (X - 1) & 1 & \frac{1}{(X - 1)} & \frac{1}{(X - 1)^{2}} & \cdots \\ (X - 1)^{6} & (X - 1)^{5} & (X - 1)^{4} & (X - 1)^{3} & (X - 1)^{2} & (X - 1) & 1 & \frac{1}{(X - 1)} & \cdots \\ (X - 1)^{7} & (X - 1)^{6} & (X - 1)^{5} & (X - 1)^{4} & (X - 1)^{3} & (X - 1)^{2} & (X - 1) & 1 & \cdots \\ \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \ddots \end{pmatrix}, \end{array}
(1.6.168)

com­posed of (X−1)±n(X - 1)^{\pm n} el­e­ments, having the in­verse re­la­tion­ship to their re­spec­tive trans­pose. Mul­ti­ply­ing the XX-matrix on the left by DX−1D_{X}^{- 1}, re­sults in the al­ter­nat­ing har­mon­ic log­a­rithm field SH,S_{H},seen below:

SH=DX−1X=(1−1213−1415−1617⋯−21−2312−2513−27⋯3−321−3435−1237⋯−42−431−4523−47⋯5−5253−541−5657⋯−63−232−651−67⋯7−7273−7475−761⋯⋮⋮⋮⋮⋮⋮⋮⋱).\begin{array}{r} {S_{H} = D}_{X}^{- 1}X = \begin{pmatrix} 1 & - \frac{1}{2} & \frac{1}{3} & - \frac{1}{4} & \frac{1}{5} & - \frac{1}{6} & \frac{1}{7} & \cdots \\ - 2 & 1 & - \frac{2}{3} & \frac{1}{2} & - \frac{2}{5} & \frac{1}{3} & - \frac{2}{7} & \cdots \\ 3 & - \frac{3}{2} & 1 & - \frac{3}{4} & \frac{3}{5} & - \frac{1}{2} & \frac{3}{7} & \cdots \\ - 4 & 2 & - \frac{4}{3} & 1 & - \frac{4}{5} & \frac{2}{3} & - \frac{4}{7} & \cdots \\ 5 & - \frac{5}{2} & \frac{5}{3} & - \frac{5}{4} & 1 & - \frac{5}{6} & \frac{5}{7} & \cdots \\ - 6 & 3 & - 2 & \frac{3}{2} & - \frac{6}{5} & 1 & - \frac{6}{7} & \cdots \\ 7 & - \frac{7}{2} & \frac{7}{3} & - \frac{7}{4} & \frac{7}{5} & - \frac{7}{6} & 1 & \cdots \\ \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \ddots \end{pmatrix}. \end{array}
(1.6.169)

We can see that SHS_{H} is com­posed of el­e­ments of the al­ter­nat­ing har­mon­ic series in the first row, and the en­su­ing rows are it­er­a­tive­ly scaled in mag­ni­tude by the se­quence of nat­ur­al num­bers, ℕ\mathbb{N}. For the columns, when moving from the left to right, each row number ∈ℕ\in \mathbb{N}, is it­er­a­tive­ly mul­ti­plied the column number value (−1)k−1k\frac{( - 1)^{k - 1}}{k}, rep­re­sent­ing the en­tries in the se­quence of the al­ter­nat­ing har­mon­ic series. So, the outer prod­uct that forms the SHS_{H} field can be rep­re­sent­ed as fol­lows

(∑k=1∞(−1)k+1k,)⊗(∑k=1∞(−1)k+1k).\begin{array}{r} \left( \sum_{k = 1}^{\infty}{( - 1)^{k + 1}k,} \right) \otimes \left( \sum_{k = 1}^{\infty}\frac{( - 1)^{k + 1}}{k} \right). \end{array}
(1.6.170)

The num­bers that form the SHS_{H} field have some very in­ter­est­ing prop­er­ties, prompt­ing us to ex­plore this in more detail in the next sec­tion. For func­tions with in­fi­nite series rep­re­sen­ta­tions, the outer prod­uct of the con­stituents in its series of el­e­ments can lead to new in­sights into the said func­tions. In a sim­i­lar manner to the ex­po­nen­tial series outer prod­ucts, the outer prod­uct of two log­a­rithm series vec­tors Xˆ,Yˆ\widehat{X},\widehat{Y}, form the outer prod­uct SXY=Xˆ⊗YˆS_{XY} = \widehat{X} \otimes \widehat{Y},

Xˆ=((X−1)11,−(X−1)22,(X−1)33,−(X−1)44,(X−1)55,−(X−1)66,…),\begin{array}{r} \widehat{X} = \left( \frac{(X - 1)^{1}}{1},\frac{{- (X - 1)}^{2}}{2},\frac{(X - 1)^{3}}{3},\frac{{- (X - 1)}^{4}}{4},\frac{(X - 1)^{5}}{5}, - \frac{(X - 1)^{6}}{6},\ldots \right), \end{array}
(1.6.171)
Yˆ=((Y−1)11,−(Y−1)22,(Y−1)33,−(Y−1)44,(Y−1)55,−(Y−1)66,…),\begin{array}{r} \widehat{Y} = \left( \frac{(Y - 1)^{1}}{1},\frac{{- (Y - 1)}^{2}}{2},\frac{(Y - 1)^{3}}{3},\frac{{- (Y - 1)}^{4}}{4},\frac{(Y - 1)^{5}}{5}, - \frac{(Y - 1)^{6}}{6},\ldots \right), \end{array}
(1.6.172)
SX,Y=((X−1)1(Y−1)11×1−(X−1)1(Y−1)21×2(X−1)1(Y−1)31×3−(X−1)1(Y−1)41×4(X−1)1(Y−1)51×5…−(X−1)2(Y−1)12×1(X−1)2(Y−1)22×2−(X−1)2(Y−1)32×3(X−1)2(Y−1)42×4−(X−1)2(Y−1)52×5…(X−1)3(Y−1)13×1−(X−1)3(Y−1)23×2(X−1)3(Y−1)33×3−(X−1)3(Y−1)43×4(X−1)3(Y−1)53×5…−(X−1)4(Y−1)14×1(X−1)4(Y−1)24×2−(X−1)4(Y−1)34×3(X−1)4(Y−1)44×4−(X−1)4(Y−1)54×5…(X−1)5(Y−1)15×1−(X−1)5(Y−1)25×2(X−1)5(Y−1)35×3−(X−1)5(Y−1)45×4(X−1)5(Y−1)55×5…⋮⋮⋮⋮⋮⋱).\begin{array}{r} S_{X,Y} = \begin{pmatrix} \frac{{(X - 1)^{1}(Y - 1)}^{1}}{1 \times 1} & - \frac{{(X - 1)^{1}(Y - 1)}^{2}}{1 \times 2} & \frac{{(X - 1)^{1}(Y - 1)}^{3}}{1 \times 3} & - \frac{{(X - 1)^{1}(Y - 1)}^{4}}{1 \times 4} & \frac{{(X - 1)^{1}(Y - 1)}^{5}}{1 \times 5} & \ldots \\ - \frac{{(X - 1)^{2}(Y - 1)}^{1}}{2 \times 1} & \frac{{(X - 1)^{2}(Y - 1)}^{2}}{2 \times 2} & - \frac{{(X - 1)^{2}(Y - 1)}^{3}}{2 \times 3} & \frac{{(X - 1)^{2}(Y - 1)}^{4}}{2 \times 4} & - \frac{{(X - 1)^{2}(Y - 1)}^{5}}{2 \times 5} & \ldots \\ \frac{{(X - 1)^{3}(Y - 1)}^{1}}{3 \times 1} & - \frac{{(X - 1)^{3}(Y - 1)}^{2}}{3 \times 2} & \frac{{(X - 1)^{3}(Y - 1)}^{3}}{3 \times 3} & - \frac{{(X - 1)^{3}(Y - 1)}^{4}}{3 \times 4} & \frac{{(X - 1)^{3}(Y - 1)}^{5}}{3 \times 5} & \ldots \\ - \frac{{(X - 1)^{4}(Y - 1)}^{1}}{4 \times 1} & \frac{{(X - 1)^{4}(Y - 1)}^{2}}{4 \times 2} & - \frac{{(X - 1)^{4}(Y - 1)}^{3}}{4 \times 3} & \frac{{(X - 1)^{4}(Y - 1)}^{4}}{4 \times 4} & - \frac{{(X - 1)^{4}(Y - 1)}^{5}}{4 \times 5} & \ldots \\ \frac{{(X - 1)^{5}(Y - 1)}^{1}}{5 \times 1} & - \frac{{(X - 1)^{5}(Y - 1)}^{2}}{5 \times 2} & \frac{{(X - 1)^{5}(Y - 1)}^{3}}{5 \times 3} & - \frac{{(X - 1)^{5}(Y - 1)}^{4}}{5 \times 4} & \frac{{(X - 1)^{5}(Y - 1)}^{5}}{5 \times 5} & \ldots \\ \vdots & \vdots & \vdots & \vdots & \vdots & \mathbf{\ddots} \end{pmatrix}. \end{array}
(1.6.173)

The trace of the SXYS_{XY} outer prod­uct matrix is

tr(SXY)=∑k=1∞(X−1)k(Y−1)kk2,\begin{array}{r} tr\left( S_{XY} \right) = \sum_{k = 1}^{\infty}{\frac{(X - 1)^{k}(Y - 1)^{k}}{k^{2}},} \end{array}
(1.6.174)

which is also a well-de­fined series known as Jon­quiere's func­tion. This func­tion is de­fined as the poly­log­a­rithm Lin(z)=∑k=1∞zkknLi_{n}(z) = \sum_{k = 1}^{\infty}\frac{z^{k}}{k^{n}}, and for the above, we find

∑k=1∞(X−1)k(Y−1)kk2=Li2((X−1)(Y−1)),   (for |(X−1)(Y−1)|<1).\begin{array}{r} \sum_{k = 1}^{\infty}\frac{(X - 1)^{k}(Y - 1)^{k}}{k^{2}} = Li_{2}\left( (X - 1)(Y - 1) \right),\ \ \ \left( for\ \left| (X - 1)(Y - 1) \right| < 1 \right). \end{array}
(1.6.175)

In the case of the poly­log­a­rithm above, Li2(z)Li_{2}(z), it is known as Spence's func­tion, a dilog­a­rithm with the in­te­gral rep­re­sen­ta­tion

Li2(z)=−∫0zln⁡(1−u)udu,z∈ℂ,\begin{array}{r} Li_{2}(z) = - \int_{0}^{z}{\frac{\ln(1 - u)}{u}du,z\in \mathbb{C}}, \end{array}
(1.6.176)

which is com­mon­ly en­coun­tered in par­ti­cle physics in the cal­cu­la­tions of ra­dia­tive cor­rec­tions. It can be noted that in the above equa­tion, ln⁡(1−u)u\frac{\ln(1 - u)}{u} is equal to the Li1(z)Li_{1}(z) poly­log­a­rithm, and that for Li2(z)Li_{2}(z), when z=(X−1)(Y−1)=1z = (X - 1)(Y - 1) = 1, it is equal to the Rie­mann zeta func­tion ζ(2),\zeta(2),which is con­ver­gent and sums to π26\ \frac{\pi^{2}}{6},

Li2(1)=ζ(2)=π26.\begin{array}{r} Li_{2}(1) = \zeta(2) = \frac{\pi^{2}}{6}. \end{array}
(1.6.177)

The fact that π2\pi^{2} ap­pears in the result is amaz­ing, and this was first proved by Leonard Euler. The sum of the en­tries in the matrix SXYS_{XY}, we find, is equal to ln⁡(Xˆ)ln⁡(Yˆ)\ln\left( \widehat{X} \right)\ln\left( \widehat{Y} \right)

∑SXY=∑n=1∞(∑k=1∞(−1)n+k (X−1)k(Y−1)nnk)=ln⁡(X)ln⁡(Y),   (for X=2∨|1−X|<1 ∧ Y=2∨ |1−Y|<1).\begin{array}{r} \sum_{}^{}{S_{XY} =}\sum_{n = 1}^{\infty}\left( \sum_{k = 1}^{\infty}\frac{( - 1)^{n + k\ }(X - 1)^{k}(Y - 1)^{n}}{nk} \right) = \ln(X)\ln(Y),\ \ \ \left( \begin{array}{r} for\ X = 2 \vee |1 - X| < 1\ \\ \land \ Y = 2 \vee \ |1 - Y| < 1 \end{array} \right). \end{array}
(1.6.178)

and for the change of vari­ables X′=X+1, Y′=Y+1,X^{'} = X + 1,\ Y^{'} = Y + 1,like we saw with the Mer­ca­tor series, we find

∑SXY′=∑n=1∞(∑k=1∞(−1)n+k (X′)k(Y′)nnk)=ln⁡(X+1)ln⁡(Y+1),  (for |X′|<1∧ |Y′|<1).\begin{array}{r} \sum_{}^{}{S_{XY}^{'} =}\sum_{n = 1}^{\infty}\left( \sum_{k = 1}^{\infty}\frac{( - 1)^{n + k\ }\left( X^{'} \right)^{k}\left( Y^{'} \right)^{n}}{nk} \right) = \ln(X + 1)\ln{(Y + 1),}\ \ \left( \begin{array}{r} for\ |X'| < 1 \land \ |Y'| < 1 \end{array} \right). \end{array}
(1.6.179)