Masteries · Jason GlowneyArithmoí Book 1

Masteries

1.5.4 The Baker-Campbell-Hausdorff Formula

When work­ing with ma­tri­ces, we cannot assume matrix prod­ucts com­mute, as is the case for many of the el­e­ments in the Lie al­ge­bra 𝔰3ℂΔ\mathfrak{s}_{3\mathbb{C}}^{\Delta}. Re­gard­ing the matrix ex­po­nen­tial of el­e­ments in the set­ting of non-com­mu­ta­tive Lie al­ge­bras, we can state that the fol­low­ing,

eXeY=eX+Y,   X,Y∈𝔰3ℂΔ\begin{array}{r} e^{X}e^{Y} = e^{X + Y},\ \ \ X,Y \in \mathfrak{s}_{3\mathbb{C}}^{\Delta} \end{array}
(1.5.32)

rep­re­sent­ing com­mu­ta­tive ma­tri­ces, does not hold for non-com­mut­ing matrix groups

eXeY≠eX+Y,\begin{array}{r} e^{X}e^{Y} \neq e^{X + Y}, \end{array}
(1.5.33)

and in­stead, ex­press the group prod­uct, as

eXeY=eZ,   Z∈𝔰3ℂΔ.\begin{array}{r} e^{X}e^{Y} = e^{Z},\ \ \ Z \in \mathfrak{s}_{3\mathbb{C}}^{\Delta}. \end{array}
(1.5.34)

The Baker-Camp­bell-Haus­dorff the­o­rem states just this, for ZZ, as the sum of the Lie al­ge­bra el­e­ments X,YX,Y, in ad­di­tion to a series of Lie brack­et terms de­fined by X,YX,Y. So, for our Lie groups, the BCH for­mu­la proves that the Lie brack­ets of 𝔰3ℂΔ\mathfrak{s}_{3\mathbb{C}}^{\Delta} de­ter­mine the prod­uct op­er­a­tions for the Lie group. We define the ad­joint op­er­a­tor adX(Y)ad_{X}(Y), as the Lie brack­et op­er­a­tion

adX(Y)=[X,Y]≡XY−YX,\begin{array}{r} ad_{X}(Y) = \lbrack X,Y\rbrack \equiv XY - YX, \end{array}
(1.5.35)

and for a pos­i­tive in­te­ger nn, nested com­mu­ta­tors can be de­scribed, as

(adX)n(Y)=[X,…[X,[X,Y]]…]⏟n.\begin{array}{r} \left( ad_{X} \right)^{n}(Y) = \underset{n}{\underbrace{\left\lbrack X,\ldots\left\lbrack X,\lbrack X,Y\rbrack \right\rbrack\ldots \right\rbrack}}. \end{array}
(1.5.36)

The­o­rem (1.5.37) – We can now state and prove that1

eXYe−X=exp⁡(adX)(Y)≡∑n=0∞1n!(adX)n(Y)=Y+[X,Y]+12[X,[X,Y]]+….\begin{array}{r} e^{X}Ye^{- X} = \exp\left( ad_{X} \right)(Y) \equiv \sum_{n = 0}^{\infty}\frac{1}{n!}\left( ad_{X} \right)^{n}(Y) = Y + \lbrack X,Y\rbrack + \frac{1}{2}\left\lbrack X,\lbrack X,Y\rbrack \right\rbrack + \ldots. \end{array}

Proof – We define

Y(t)=etXYe−tX,\begin{array}{r} Y(t) = e^{tX}Ye^{- tX}, \end{array}

and com­put­ing the Taylor series of Y(t)Y(t), around the point, t=0t = 0, shows that Y(0)=YY(0) = Y, and

dY(t)dt=XetXYe−tX−etXYe−tXX=[X,Y(t)]=adX(Y(t)).\begin{array}{r} \frac{dY(t)}{dt} = Xe^{tX}Ye^{- tX} - e^{tX}Ye^{- tX}X = \left\lbrack X,Y(t) \right\rbrack = ad_{X}\left( Y(t) \right). \end{array}

Higher de­riv­a­tives are simple to cal­cu­late and relate to nested ad­joint op­er­a­tors, by

dnY(t)dtn=(adX)n(Y(t)).\begin{array}{r} \frac{d^{n}Y(t)}{dt^{n}} = \left( ad_{X} \right)^{n}\left( Y(t) \right). \end{array}

By sub­sti­tut­ing t=1t = 1 into the Taylor series, we have sat­is­fied the ini­tial state­ment eXYe−X=exp⁡(adX)(Y)e^{X}Ye^{- X} = \exp\left( ad_{X} \right)(Y).

 □\begin{array}{r} \ \square \end{array}

We define the fol­low­ing aux­il­iary func­tions and their power series

f(z)=ez−1z=∑n=0∞zn(n+1)!,  |z|<∞,\begin{array}{r} f(z) = \frac{e^{z} - 1}{z} = \sum_{n = 0}^{\infty}{\frac{z^{n}}{(n + 1)!},\ \ |z| < \infty}, \end{array}
(1.5.38)
g(z)=ln⁡(z)z−1=∑n=0∞(1−z)nn+1,  |1−z|<1,\begin{array}{r} g(z) = \frac{\ln(z)}{z - 1} = \sum_{n = 0}^{\infty}\frac{(1 - z)^{n}}{n + 1},\ \ |1 - z| < 1, \end{array}
(1.5.39)

that sat­is­fy

f(ln⁡(z))g(z)=1,  for |1−z|<1,\begin{array}{r} f\left( \ln(z) \right)g(z) = 1,\ \ for\ |1 - z| < 1, \end{array}
(1.5.40)
f(z)g(ez)=1,  for |z|<∞,\begin{array}{r} f(z)g\left( e^{z} \right) = 1,\ \ for\ |z| < \infty, \end{array}
(1.5.41)

to for­mal­ly state the Baker-Camp­bell-Haus­dorff for­mu­la

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.5.42)

for g(z)g(z), the aux­il­iary func­tion we just de­fined, and the BCH for­mu­la for ln⁡(eXeY)\ln\left( e^{X}e^{Y} \right) will con­verge pro­vid­ed that ‖eXeY−ϵ‖<1\left\| e^{X}e^{Y} - \epsilon \right\| < 1, for ϵ\epsilon, the iden­ti­ty matrix and ‖…‖\left\| \ldots \right\|, the matrix norm. We can now expand the BCH for­mu­la using the Taylor series de­f­i­n­i­tion of g(z)g(z), to show that

eXeY=exp⁡(X+Y+12[X,Y]+112[X,[X,Y]+112[Y,[Y,X]]+…]),\begin{array}{r} e^{X}e^{Y} = \exp{\left( X + Y + \frac{1}{2}\lbrack X,Y\rbrack + \frac{1}{12}\left\lbrack X,\lbrack X,Y\rbrack + \frac{1}{12}\left\lbrack Y,\lbrack Y,X\rbrack \right\rbrack + \ldots \right\rbrack \right),} \end{array}
(1.5.43)

holds for con­ver­gent series. To demon­strate this con­cep­tu­al­ly, we use the fol­low­ing two simple el­e­ments X,YX,Y

X=(iπ20000iπ20iπ20),  Y=(00iπ20iπ20iπ200),\begin{array}{r} X = \begin{pmatrix} \frac{i\pi}{2} & 0 & 0 \\ 0 & 0 & \frac{i\pi}{2} \\ 0 & \frac{i\pi}{2} & 0 \end{pmatrix},\ \ Y = \begin{pmatrix} 0 & 0 & \frac{i\pi}{2} \\ 0 & \frac{i\pi}{2} & 0 \\ \frac{i\pi}{2} & 0 & 0 \end{pmatrix}, \end{array}
(1.5.44)

and com­pute the fol­low­ing matrix ex­po­nen­tials ∈G\in G

eX+Y=((−13+i3)(−13−2i3)(−13+i3)(−13−2i3)(−13+i3)(−13+i3)(−13+i3)(−13+i3)(−13−2i3)),\begin{array}{r} e^{X + Y} = \begin{pmatrix} \left( - \frac{1}{3} + \frac{i}{3} \right) & \left( - \frac{1}{3} - \frac{2i}{3} \right) & \left( - \frac{1}{3} + \frac{i}{3} \right) \\ \left( - \frac{1}{3} - \frac{2i}{3} \right) & \left( - \frac{1}{3} + \frac{i}{3} \right) & \left( - \frac{1}{3} + \frac{i}{3} \right) \\ \left( - \frac{1}{3} + \frac{i}{3} \right) & \left( - \frac{1}{3} + \frac{i}{3} \right) & \left( - \frac{1}{3} - \frac{2i}{3} \right) \end{pmatrix}, \end{array}
(1.5.45)
eZ=eXeY=(00−1−1000−10).\begin{array}{r} e^{Z} = e^{X}e^{Y} = \begin{pmatrix} 0 & 0 & - 1 \\ - 1 & 0 & 0 \\ 0 & - 1 & 0 \end{pmatrix}. \end{array}
(1.5.46)

For eX+Ye^{X + Y}, the matrix log­a­rithm simply equals the sum of the orig­i­nal el­e­ments X,YX,Y

ln⁡(eX+Y)=X+Y=(iπ20iπ20iπ2iπ2iπ2iπ20),\begin{array}{r} \ln{\left( e^{X + Y} \right) =}X + Y = \begin{pmatrix} \frac{i\pi}{2} & 0 & \frac{i\pi}{2} \\ 0 & \frac{i\pi}{2} & \frac{i\pi}{2} \\ \frac{i\pi}{2} & \frac{i\pi}{2} & 0 \end{pmatrix}, \end{array}
(1.5.47)

and we in­ci­den­tal­ly see that eZe^{Z} equals the Lie al­ge­bra −Y3∈𝔰3ℂΔ.- Y_{3} \in \mathfrak{s}_{3\mathbb{C}}^{\Delta}.To find the tan­gent space for eZe^{Z}, we take ad­van­tage of the fact that these ma­tri­ces are di­ag­o­nal­iz­able, ex­press­ing ln⁡(eZ)\ln\left( e^{Z} \right), as

ln⁡(eZ)=VDV−1,\begin{array}{r} \ln\left( e^{Z} \right) = VDV^{- 1}, \end{array}
(1.5.48)

for VV, a uni­tary square matrix of the eigen­vec­tors of eZe^{Z}, and DD, the nat­ur­al log of the di­ag­o­nal matrix of eigen­val­ues of eZe^{Z}. Cal­cu­lat­ing D,VD,V, and V−1V^{- 1}, for eZe^{Z}, gives us the fol­low­ing ma­tri­ces

D=(12−i3200012+i32000−1),V=((−12−i32)3(−12+i32)313(−12+i32)3(−12−i32)313131313),   V−1=((−12+i32)3(−12−i32)313(−12−i32)3(−12+i32)313131313),\begin{array}{r} D = \begin{pmatrix} \frac{1}{2} - \frac{i\sqrt{3}}{2} & 0 & 0 \\ 0 & \frac{1}{2} + \frac{i\sqrt{3}}{2} & 0 \\ 0 & 0 & - 1 \end{pmatrix},V = \begin{pmatrix} \frac{\left( - \frac{1}{2} - \frac{i\sqrt{3}}{2} \right)}{\sqrt{3}} & \frac{\left( - \frac{1}{2} + \frac{i\sqrt{3}}{2} \right)}{\sqrt{3}} & \frac{1}{\sqrt{3}} \\ \frac{\left( - \frac{1}{2} + \frac{i\sqrt{3}}{2} \right)}{\sqrt{3}} & \frac{\left( - \frac{1}{2} - \frac{i\sqrt{3}}{2} \right)}{\sqrt{3}} & \frac{1}{\sqrt{3}} \\ \frac{1}{\sqrt{3}} & \frac{1}{\sqrt{3}} & \frac{1}{\sqrt{3}} \end{pmatrix},\ \ \ V^{- 1} = \begin{pmatrix} \frac{\left( - \frac{1}{2} + \frac{i\sqrt{3}}{2} \right)}{\sqrt{3}} & \frac{\left( - \frac{1}{2} - \frac{i\sqrt{3}}{2} \right)}{\sqrt{3}} & \frac{1}{\sqrt{3}} \\ \frac{\left( - \frac{1}{2} - \frac{i\sqrt{3}}{2} \right)}{\sqrt{3}} & \frac{\left( - \frac{1}{2} + \frac{i\sqrt{3}}{2} \right)}{\sqrt{3}} & \frac{1}{\sqrt{3}} \\ \frac{1}{\sqrt{3}} & \frac{1}{\sqrt{3}} & \frac{1}{\sqrt{3}} \end{pmatrix}, \end{array}
(1.5.49)

and the matrix log of DD, equals

D′=(−iπ3000iπ3000−iπ).\begin{array}{r} D^{'} = \begin{pmatrix} - \frac{i\pi}{3} & 0 & 0 \\ 0 & \frac{i\pi}{3} & 0 \\ 0 & 0 & - i\pi \end{pmatrix}. \end{array}
(1.5.50)

The matrix log­a­rithm cal­cu­la­tion be­comes

ln⁡(eZ)=Z=VD′V−1=(−iπ313(π3−iπ)−13(π3+iπ)−13(π3+iπ)−iπ313(π3−iπ)13(π3−iπ)−13(π3+iπ)−iπ3),\begin{array}{r} \ln{\left( e^{Z} \right) = Z = VD^{'}V^{- 1} = \begin{pmatrix} - \frac{i\pi}{3} & \frac{1}{3}\left( \frac{\pi}{\sqrt{3}} - i\pi \right) & - \frac{1}{3}\left( \frac{\pi}{\sqrt{3}} + i\pi \right) \\ - \frac{1}{3}\left( \frac{\pi}{\sqrt{3}} + i\pi \right) & - \frac{i\pi}{3} & \frac{1}{3}\left( \frac{\pi}{\sqrt{3}} - i\pi \right) \\ \frac{1}{3}\left( \frac{\pi}{\sqrt{3}} - i\pi \right) & - \frac{1}{3}\left( \frac{\pi}{\sqrt{3}} + i\pi \right) & - \frac{i\pi}{3} \end{pmatrix}}, \end{array}
(1.5.51)

and to verify that eln⁡(eZ)e^{\ln\left( e^{Z} \right)} does indeed equal our ini­tial result, we find

eln⁡(eZ)=exp(−iπ313(π3−iπ)−13(π3+iπ)−13(π3+iπ)−iπ313(π3−iπ)13(π3−iπ)−13(π3+iπ)−iπ3)=(00−1−1000−10).\begin{array}{r} e^{\ln\left( e^{Z} \right)} = exp\begin{pmatrix} - \frac{i\pi}{3} & \frac{1}{3}\left( \frac{\pi}{\sqrt{3}} - i\pi \right) & - \frac{1}{3}\left( \frac{\pi}{\sqrt{3}} + i\pi \right) \\ - \frac{1}{3}\left( \frac{\pi}{\sqrt{3}} + i\pi \right) & - \frac{i\pi}{3} & \frac{1}{3}\left( \frac{\pi}{\sqrt{3}} - i\pi \right) \\ \frac{1}{3}\left( \frac{\pi}{\sqrt{3}} - i\pi \right) & - \frac{1}{3}\left( \frac{\pi}{\sqrt{3}} + i\pi \right) & - \frac{i\pi}{3} \end{pmatrix} = \begin{pmatrix} 0 & 0 & - 1 \\ - 1 & 0 & 0 \\ 0 & - 1 & 0 \end{pmatrix}. \end{array}
(1.5.52)

As we have shown, the Baker-Camp­bell-Haus­dorff for­mu­la is the series of Lie brack­et op­er­a­tions of the form

eZ=eXeY=exp⁡(X+Y+12[X,Y]+112[X,[X,Y]+112[Y,[Y,X]]+…]),\begin{array}{r} {e^{Z} = e}^{X}e^{Y} = \exp{\left( X + Y + \frac{1}{2}\lbrack X,Y\rbrack + \frac{1}{12}\left\lbrack X,\lbrack X,Y\rbrack + \frac{1}{12}\left\lbrack Y,\lbrack Y,X\rbrack \right\rbrack + \ldots \right\rbrack \right),} \end{array}
(1.5.53)

in which the por­tion of the series X+YX + Y, is the com­mu­ni­tive al­ge­bra part of ZZ, while the re­main­der is its non-com­mu­ni­tive part. From what we cal­cu­lat­ed in our ex­am­ple, we can see that the non-com­mu­ni­tive por­tion of ZZ is

Z−(X+Y)=(12[X,Y]+112[X,[X,Y]]+112[[Y,[Y,X]]]+…),\begin{array}{r} Z - (X + Y) = \left( \frac{1}{2}\lbrack X,Y\rbrack + \frac{1}{12}\left\lbrack X,\lbrack X,Y\rbrack \right\rbrack + \frac{1}{12}\left\lbrack \left\lbrack Y,\lbrack Y,X\rbrack \right\rbrack \right\rbrack + \ldots \right), \end{array}
(1.5.54)

and this can be ex­pressed, as

Z−(X+Y)=(−5iπ613(π3−iπ)−13(π3+5iπ2)−13(π3+iπ)−5iπ613(π3−5iπ2)13(π3−i5π2)−13(π3+5iπ2)−iπ3).\begin{array}{r} Z - (X + Y) = \begin{pmatrix} - \frac{5i\pi}{6} & \frac{1}{3}\left( \frac{\pi}{\sqrt{3}} - i\pi \right) & - \frac{1}{3}\left( \frac{\pi}{\sqrt{3}} + \frac{5i\pi}{2} \right) \\ - \frac{1}{3}\left( \frac{\pi}{\sqrt{3}} + i\pi \right) & - \frac{5i\pi}{6} & \frac{1}{3}\left( \frac{\pi}{\sqrt{3}} - \frac{5i\pi}{2} \right) \\ \frac{1}{3}\left( \frac{\pi}{\sqrt{3}} - \frac{i5\pi}{2} \right) & - \frac{1}{3}\left( \frac{\pi}{\sqrt{3}} + \frac{5i\pi}{2} \right) & - \frac{i\pi}{3} \end{pmatrix}. \end{array}
(1.5.55)

With a basic un­der­stand­ing of the BCH for­mu­la, we can apply it to 2 el­e­ments in the 𝔰3ℂΔ\mathfrak{s}_{3\mathbb{C}}^{\Delta} Lie al­ge­bra for X1,X2∈U(3)X_{1}{,X}_{2} \in U(3), de­fined

X1=(i0000i0i0),  X2=(00i0i0i00).\begin{array}{r} X_{1} = \begin{pmatrix} i & 0 & 0 \\ 0 & 0 & i \\ 0 & i & 0 \end{pmatrix},{\ \ X}_{2} = \begin{pmatrix} 0 & 0 & i \\ 0 & i & 0 \\ i & 0 & 0 \end{pmatrix}. \end{array}
(1.5.56)

We can de­ter­mine the matrix ex­po­nen­tial eX1+X2e^{X_{1} + X_{2}} in a sim­i­lar way to what we did for the U(3)U(3) group, by using two vari­ables, x1,x2x_{1},x_{2}, in­stead of the three vari­ables. The group form can be trigono­met­ri­cal­ly rep­re­sent­ed as

(eiγ3+2cosh⁡(r)3+∂r∂x1∙2sinh(r)3eiγ3−cosh⁡(r)3−γr∙2sinh(r)3eiγ3−cosh⁡(r)3+∂r∂x2∙2sinh(r)3eiγ3−cosh⁡(r)3−γr∙2sinh(r)3eiγ3+2cosh⁡(r)3+∂r∂x2∙2sinh(r)3eiγ3−cosh⁡(r)3+∂r∂x1∙2sinh(r)3eiγ3−cosh⁡(r)3+∂r∂x2∙2sinh(r)3eiγ3−cosh⁡(r)3+∂r∂x1∙2sinh(r)3eiγ3+2cosh⁡(r)3−γr∙2sinh(r)3),\begin{array}{r} \begin{pmatrix} \frac{e^{i\gamma}}{3} + \frac{{2\cosh}(r)}{3} + \frac{\partial r}{\partial x_{1}} \bullet \frac{2sinh(r)}{3} & \frac{e^{i\gamma}}{3} - \frac{\cosh(r)}{3} - \frac{\gamma}{r} \bullet \frac{2sinh(r)}{3} & \frac{e^{i\gamma}}{3} - \frac{\cosh(r)}{3} + \frac{\partial r}{\partial x_{2}} \bullet \frac{2sinh(r)}{3} \\ \frac{e^{i\gamma}}{3} - \frac{\cosh(r)}{3} - \frac{\gamma}{r} \bullet \frac{2sinh(r)}{3} & \frac{e^{i\gamma}}{3} + \frac{{2\cosh}(r)}{3} + \frac{\partial r}{\partial x_{2}} \bullet \frac{2sinh(r)}{3} & \frac{e^{i\gamma}}{3} - \frac{\cosh(r)}{3} + \frac{\partial r}{\partial x_{1}} \bullet \frac{2sinh(r)}{3} \\ \frac{e^{i\gamma}}{3} - \frac{\cosh(r)}{3} + \frac{\partial r}{\partial x_{2}} \bullet \frac{2sinh(r)}{3} & \frac{e^{i\gamma}}{3} - \frac{\cosh(r)}{3} + \frac{\partial r}{\partial x_{1}} \bullet \frac{2sinh(r)}{3} & \frac{e^{i\gamma}}{3} + \frac{{2\cosh}(r)}{3} - \frac{\gamma}{r} \bullet \frac{2sinh(r)}{3} \end{pmatrix}, \end{array}
(1.5.57)

for r=(x12−x1x2+x22),  γ=(x1+x2),  ∂r∂x1=(x1−x22r),  ∂r∂x2=(x2−x12r).r = \left( \sqrt{x_{1}^{2} - x_{1}x_{2} + x_{2}^{2}} \right),\ \ \gamma = \left( x_{1} + x_{2} \right),\ \ \frac{\partial r}{\partial x_{1}} = \left( \frac{x_{1} - \frac{x_{2}}{2}}{r} \right),\ \ \frac{\partial r}{\partial x_{2}} = \left( \frac{x_{2} - \frac{x_{1}}{2}}{r} \right).

For the X1,X2X_{1},X_{2}, el­e­ments, we would have that x1=ix_{1} = i, x2=ix_{2} = i, and the ex­po­nen­tial map­ping eX1+X2e^{X_{1} + X_{2}} is

eX1+X2=(13(e2i+2cos⁡(1)+isin(1))13(e2i−cos⁡(1)−2isin(1))13(e2i−cos⁡(1)+isin(1))13(e2i−cos⁡(1)−2isin(1))13(e2i+2cos(1)+isin(1))13(e2i−cos⁡(1)+isin(1))13(e2i−cos⁡(1)+isin(1))13(e2i−cos⁡(1)+isin(1))13(e2i+2cos(1)−2isin(1))).\begin{array}{r} e^{X_{1} + X_{2}} = \begin{pmatrix} \frac{1}{3}\left( e^{2i} + {2\cos}(1) + isin(1) \right) & \frac{1}{3}\left( e^{2i} - \cos(1) - 2isin(1) \right) & \frac{1}{3}\left( e^{2i} - \cos(1) + isin(1) \right) \\ \frac{1}{3}\left( e^{2i} - \cos(1) - 2isin(1) \right) & \frac{1}{3}\left( e^{2i} + {2cos}(1) + isin(1) \right) & \frac{1}{3}\left( e^{2i} - \cos(1) + isin(1) \right) \\ \frac{1}{3}\left( e^{2i} - \cos(1) + isin(1) \right) & \frac{1}{3}\left( e^{2i} - \cos(1) + isin(1) \right) & \frac{1}{3}\left( e^{2i} + {2cos}(1) - 2isin(1) \right) \end{pmatrix}. \end{array}
(1.5.58)

For the result of the prod­uct eX1eX2=eZe^{X_{1}}e^{X_{2}} = e^{Z}, having the trigono­met­ric form

exp⁡(ix10000ix10ix10)×exp(,00ix20ix20ix200)=(cos⁡(x2)eix10isin(x2)eix1−sin⁡(x1)sin⁡(x2)cos⁡(x1)eix2isin(x1)cos⁡(x2)isin(x2)cos⁡(x1)isin(x1)eix2cos⁡(x1)cos⁡(x2))\begin{array}{r} \exp\begin{pmatrix} ix_{1} & 0 & 0 \\ 0 & 0 & ix_{1} \\ 0 & ix_{1} & 0 \end{pmatrix} \times exp\left( ,\begin{matrix} 0 & 0 & ix_{2} \\ 0 & ix_{2} & 0 \\ ix_{2} & 0 & 0 \end{matrix} \right) = \begin{pmatrix} {\cos\left( x_{2} \right)e}^{ix_{1}} & 0 & isin\left( x_{2} \right)e^{ix_{1}} \\ - \sin{\left( x_{1} \right)\sin\left( x_{2} \right)} & {\cos\left( x_{1} \right)e}^{ix_{2}} & isin\left( x_{1} \right)\cos\left( x_{2} \right) \\ isin\left( x_{2} \right)\cos\left( x_{1} \right) & isin\left( x_{1} \right)e^{ix_{2}} & \cos\left( x_{1} \right)\cos\left( x_{2} \right) \end{pmatrix} \end{array}
(1.5.59)

where once again, x1=1x_{1} = 1 and x2=1x_{2} = 1. As eZe^{Z} is a di­ag­o­nal­iz­able matrix, eZe^{Z} can be ex­pressed, as

eZ=VDV−1,\begin{array}{r} e^{Z} = VDV^{- 1}, \end{array}
(1.5.60)

where once again for the 3×33 \times 3 square ma­tri­ces, the VV matrix rep­re­sents the eZe^{Z} eigen­vec­tors, and the DD matrix, the di­ag­o­nal­ized eZe^{Z} eigen­val­ues. For the sake of com­plete­ness, any form eX1eX2=eZ,e^{X_{1}}e^{X_{2}} = e^{Z},for X1,X2,∈𝔲(3),X_{1},X_{2}, \in \mathfrak{u}(3),where 𝔲(3)\mathfrak{u}(3) rep­re­sents Lie al­ge­bra el­e­ments of U(3)U(3), the eigen­vec­tor and eigen­val­ues ma­tri­ces take the form

D=(λ1000λ2000λ3),\begin{array}{r} D = \begin{pmatrix} \lambda_{1} & 0 & 0 \\ 0 & \lambda_{2} & 0 \\ 0 & 0 & \lambda_{3} \end{pmatrix}, \end{array}
(1.5.61)

for eigen­val­ues λk\lambda_{k}:

λ1=sin⁡(x1)sin⁡(x2)2+cos⁡(x1)cos⁡(x2)−(3cos⁡(2x1)+3cos⁡(2x2)+5cos⁡(2x1)cos⁡(2x2)+4sin⁡(2x1)sin⁡(2x2)−11))4,\lambda_{1} = \frac{\sin{\left( x_{1} \right)\sin\left( x_{2} \right)}}{2}{+ \cos\left( x_{1} \right)\cos{\left( x_{2} \right) - \frac{\left( \sqrt{3\cos{\left( 2x_{1} \right) + 3\cos{\left( 2x_{2} \right) + 5\cos{\left( 2x_{1} \right)\cos{\left( 2x_{2} \right) + 4\sin{\left( 2x_{1} \right)\sin{\left( 2x_{2} \right) - 11)}}}}}}} \right)}{4}}},
λ2=sin⁡(x1)sin⁡(x2)2+cos⁡(x1)cos⁡(x2)+(3cos⁡(2x1)+3cos⁡(2x2)+5cos⁡(2x1)cos⁡(2x2)+4sin⁡(2x1)sin⁡(2x2)−11))4,\lambda_{2} = \frac{\sin{\left( x_{1} \right)\sin\left( x_{2} \right)}}{2}{+ \cos\left( x_{1} \right)\cos{\left( x_{2} \right) + \frac{\left( \sqrt{3\cos{\left( 2x_{1} \right) + 3\cos{\left( 2x_{2} \right) + 5\cos{\left( 2x_{1} \right)\cos{\left( 2x_{2} \right) + 4\sin{\left( 2x_{1} \right)\sin{\left( 2x_{2} \right) - 11)}}}}}}} \right)}{4}}},
λ3=eix1+ix2,\begin{array}{r} \lambda_{3} = e^{{ix}_{1} + {ix}_{2}}, \end{array}
(1.5.62)

and VV, a matrix com­posed of the eigen­vec­tor columns:

V1=((sin⁡(x1)sin⁡(x2)−2isin(x1)cos⁡(x2)+123cos⁡(2x1)+3cos⁡(2x2)+5cos⁡(2x1)cos⁡(2x2)+4sin⁡(2x1)sin⁡(2x2)−11)2(isin(x1)cos⁡(x2)−sin⁡(x2)(sin⁡(x1)+icos(x1))) (sin⁡(x1)sin⁡(x2)+2isin(x2)cos⁡(x1)−123cos⁡(2x1)+3cos⁡(2x2)+5cos⁡(2x1)cos⁡(2x2)+4sin⁡(2x1)sin⁡(2x2)−11)2(isin(x1)cos⁡(x2)−sin⁡(x2)(sin⁡(x1)+icos(x1)))1),V_{1} = \begin{pmatrix} \frac{(\sin{\left( x_{1} \right)\sin{\left( x_{2} \right) - {2isin\left( x_{1} \right)\cos}{\left( x_{2} \right) + \frac{1}{2}\sqrt{3\cos{\left( 2x_{1} \right) + 3\cos{\left( 2x_{2} \right) + 5\cos{\left( 2x_{1} \right)\cos{\left( 2x_{2} \right) + 4\sin{\left( 2x_{1} \right)\sin{\left( 2x_{2} \right) - 11)}}}}}}}}}}}{2\left( isin\left( x_{1} \right)\cos{\left( x_{2} \right) - \sin{\left( x_{2} \right)\left( \sin{\left( x_{1} \right) + icos\left( x_{1} \right)} \right)}} \right)}\ \\ \frac{(\sin{\left( x_{1} \right)\sin{\left( x_{2} \right) + {2isin\left( x_{2} \right)\cos}{\left( x_{1} \right) - \frac{1}{2}\sqrt{3\cos{\left( 2x_{1} \right) + 3\cos{\left( 2x_{2} \right) + 5\cos{\left( 2x_{1} \right)\cos{\left( 2x_{2} \right) + 4\sin{\left( 2x_{1} \right)\sin{\left( 2x_{2} \right) - 11)}}}}}}}}}}}{2\left( isin\left( x_{1} \right)\cos{\left( x_{2} \right) - \sin{\left( x_{2} \right)\left( \sin{\left( x_{1} \right) + icos\left( x_{1} \right)} \right)}} \right)} \\ 1 \end{pmatrix},
V2=((sin⁡(x1)sin⁡(x2)−2isin(x1)cos⁡(x2)−123cos⁡(2x1)+3cos⁡(2x2)+5cos⁡(2x1)cos⁡(2x2)+4sin⁡(2x1)sin⁡(2x2)−11)2(isin(x1)cos⁡(x2)−sin⁡(x2)(sin⁡(x1)+icos(x1))) (sin⁡(x1)sin⁡(x2)+2isin(x2)cos⁡(x1)+123cos⁡(2x1)+3cos⁡(2x2)+5cos⁡(2x1)cos⁡(2x2)+4sin⁡(2x1)sin⁡(2x2)−11)2(isin(x1)cos⁡(x2)−sin⁡(x2)(sin⁡(x1)+icos(x1)))1),V_{2} = \begin{pmatrix} \frac{(\sin{\left( x_{1} \right)\sin{\left( x_{2} \right) - {2isin\left( x_{1} \right)\cos}{\left( x_{2} \right) - \frac{1}{2}\sqrt{3\cos{\left( 2x_{1} \right) + 3\cos{\left( 2x_{2} \right) + 5\cos{\left( 2x_{1} \right)\cos{\left( 2x_{2} \right) + 4\sin{\left( 2x_{1} \right)\sin{\left( 2x_{2} \right) - 11)}}}}}}}}}}}{2\left( isin\left( x_{1} \right)\cos{\left( x_{2} \right) - \sin{\left( x_{2} \right)\left( \sin{\left( x_{1} \right) + icos\left( x_{1} \right)} \right)}} \right)}\ \\ \frac{(\sin{\left( x_{1} \right)\sin{\left( x_{2} \right) + {2isin\left( x_{2} \right)\cos}{\left( x_{1} \right) + \frac{1}{2}\sqrt{3\cos{\left( 2x_{1} \right) + 3\cos{\left( 2x_{2} \right) + 5\cos{\left( 2x_{1} \right)\cos{\left( 2x_{2} \right) + 4\sin{\left( 2x_{1} \right)\sin{\left( 2x_{2} \right) - 11)}}}}}}}}}}}{2\left( isin\left( x_{1} \right)\cos{\left( x_{2} \right) - \sin{\left( x_{2} \right)\left( \sin{\left( x_{1} \right) + icos\left( x_{1} \right)} \right)}} \right)} \\ 1 \end{pmatrix},
V3=(1 11).\begin{array}{r} V_{3} = \begin{pmatrix} 1\ \\ 1 \\ 1 \end{pmatrix}. \end{array}
(1.5.63)

Our ex­am­ple sim­pli­fies since x1=x2=1,x_{1} = x_{2} = 1,and the ma­tri­ces DD and VV become

D=(12+cos2⁡(1)2−i(6+sin2⁡(2)−6cos⁡(2))400012+cos2⁡(1)2+i(6+sin2⁡(2)−6cos⁡(2))4000e2i),\begin{array}{r} D = \begin{pmatrix} \frac{1}{2} + \frac{\cos^{2}(1)}{2}{- \frac{i\left( \sqrt{6 + \sin^{2}(2) - 6\cos(2)} \right)}{4}} & 0 & 0 \\ 0 & \frac{1}{2} + \frac{\cos^{2}(1)}{2}{+ \frac{i\left( \sqrt{6 + \sin^{2}(2) - 6\cos(2)} \right)}{4}} & 0 \\ 0 & 0 & e^{2i} \end{pmatrix}, \end{array}
(1.5.64)
V=(sin2⁡(1)−isin(2)+i2(6+sin2⁡(2)−6cos⁡(2))cos⁡(2)−1sin2⁡(1)−isin(2)−i2(6+sin2⁡(2)−6cos⁡(2))cos⁡(2)−11sin2⁡(1)+isin(2)−i2(6+sin2⁡(2)−6cos⁡(2))cos⁡(2)−1sin2⁡(1)+isin(2)+i2(6+sin2⁡(2)−6cos⁡(2))cos⁡(2)−11111).\begin{array}{r} V = \begin{pmatrix} \frac{\sin^{2}{(1) - isin(2) + \frac{i}{2}\left( \sqrt{6 + \sin^{2}(2) - 6\cos(2)} \right)}}{\cos(2) - 1} & \frac{\sin^{2}{(1) - isin(2) - \frac{i}{2}\left( \sqrt{6 + \sin^{2}(2) - 6\cos(2)} \right)}}{\cos(2) - 1} & 1 \\ \frac{\sin^{2}{(1) + isin(2) - \frac{i}{2}\left( \sqrt{6 + \sin^{2}(2) - 6\cos(2)} \right)}}{\cos(2) - 1} & \frac{\sin^{2}{(1) + isin(2) + \frac{i}{2}\left( \sqrt{6 + \sin^{2}(2) - 6\cos(2)} \right)}}{\cos(2) - 1} & 1 \\ 1 & 1 & 1 \end{pmatrix}. \end{array}
(1.5.65)

Nu­mer­i­cal­ly the ma­tri­ces above ap­prox­i­mate to

D≈(0.6460−0.7634i0000.6460+0.7634i000−0.4161+0.9093i),\begin{array}{r} D \approx \begin{pmatrix} 0.6460 - 0.7634i & 0 & 0 \\ 0 & 0.6460 + 0.7634i & 0 \\ 0 & 0 & - 0.4161 + 0.9093i \end{pmatrix}, \end{array}
(1.5.66)
V≈(−0.5000−0.4360i−0.5000+1.7202i1−0.5000+0.4360−0.5000−1.7202i1111),\begin{array}{r} V \approx \begin{pmatrix} - 0.5000 - 0.4360i & - 0.5000 + 1.7202i & 1 \\ - 0.5000 + 0.4360 & - 0.5000 - 1.7202i & 1 \\ 1 & 1 & 1 \end{pmatrix}, \end{array}
(1.5.67)

and by cal­cu­lat­ing matrix log­a­rithm of the di­ag­o­nal­ized matrix DD, as simply the log­a­rithms of the in­di­vid­ual en­tries on the main di­ag­o­nal, the result for ln⁡(D)\ln(D), is

ln⁡(D)≈(−0.8685i0000.8685i0002i).\begin{array}{r} \ln(D) \approx \begin{pmatrix} - 0.8685i & 0 & 0 \\ 0 & 0.8685i & 0 \\ 0 & 0 & 2i \end{pmatrix}. \end{array}
(1.5.68)

Solv­ing for ln⁡(eX1eX2)=Z\ {\ln(e}^{X_{1}}e^{X_{2}}) = Z

ln⁡(eX1eX2)=Z= VDV−1≈(0.8391i0.4028+0.3218i−0.4028+0.8391i−0.4028+0.3218i0.8391i0.4028+0.8391i0.4028+0.8391i−0.4028+0.8391i0.3218i),\begin{array}{r} {{\ln(e}^{X_{1}}e^{X_{2}}) = Z = \ }{VDV^{- 1} \approx \begin{pmatrix} 0.8391i & 0.4028 + 0.3218i & - 0.4028 + 0.8391i \\ - 0.4028 + 0.3218i & 0.8391i & 0.4028 + 0.8391i \\ 0.4028 + 0.8391i & - 0.4028 + 0.8391i & 0.3218i \end{pmatrix}}, \end{array}
(1.5.69)

gives us the ZZ matrix for eZe^{Z} ∈U(3)⊂S3ℂΔ\in U(3) \subset S_{3\mathbb{C}}^{\Delta}, rep­re­sent­ing the same result the Baker-Camp­bell-Haus­dorff for­mu­la would yield.

To sum­ma­rize, the Baker-Camp­bell-Haus­dorff for­mu­la can be con­sid­ered as fol­lows2: the el­e­ments X1,X2∈𝔲(3)⊂ 𝔰3ℂΔX_{1},X_{2} \in {\mathfrak{u}(3)\subset \ \mathfrak{s}}_{3\mathbb{C}}^{\Delta} are prim­i­tive; their matrix ex­po­nen­tials eX1,eX2e^{X_{1}},e^{X_{2}} are grou­p­like, as is their prod­uct eX1eX2e^{X_{1}}e^{X_{2}}; so, the matrix log­a­rithm of eX1eX2e^{X_{1}}e^{X_{2}} is prim­i­tive and can be writ­ten as the in­fi­nite sum of Lie brack­et el­e­ments gen­er­at­ed by X1,X2.X_{1},X_{2}.

  1. Haber, Howard – "Notes on the Matrix Ex­po­nen­tial and Log­a­rithm", 2019 http://scipp.ucsc.edu/~haber/web­page/Ma­trix­Ex­pLog.pdf
  2. Nathan, Ja­cob­son,* "Lie Al­ge­bras", John Wiley & Sons, 1966*