Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Lie Groups II

Adjoint Representation and Campbell Identity

We define AdX(Y)\mathrm{Ad}_{X}(Y) by:

AdX(Y)=XYX1\mathrm{Ad}_{X}(Y) = X Y X^{-1}

and define a matrix function f(t)f(t) by:

f(t)=AdetX(Y)=etXYetX.f(t) = \mathrm{Ad}_{e^{tX}}(Y) = e^{tX} Y e^{-tX}\,.

Differentiate f(t)f(t):

f(t)=ddtAdetX(Y)=ddtetXYetX=f'(t) = {d\over d t} \mathrm{Ad}_{e^{tX}}(Y) = {d\over d t} e^{tX} Y e^{-tX} =
=XetXYetXetXYXetX== X e^{tX} Y e^{-tX} - e^{tX} Y X e^{-tX} =
=etXXYetXetXYXetX== e^{tX} X Y e^{-tX} - e^{tX} Y X e^{-tX} =
=etX(XYYX)etX== e^{tX} (XY - YX) e^{-tX} =
=etX[X,Y]etX.= e^{tX} [X, Y] e^{-tX}\,.

We apply it recursively for the second derivative:

f(t)=d2dt2AdetX(Y)=d2dt2etXYetX=f''(t) = {d^2\over d t^2} \mathrm{Ad}_{e^{tX}}(Y) = {d^2\over d t^2} e^{tX} Y e^{-tX} =
=ddtetX[X,Y]etX== {d\over d t} e^{tX} [X, Y] e^{-tX} =
=etX[X,[X,Y]]etX.= e^{tX} [X, [X, Y]] e^{-tX}\,.

We evaluate these derivatives at t=0t=0:

f(0)=Y,f(0)=Y\,,

f(0)=[X,Y],f'(0)=[X, Y]\,,

f(0)=[X,[X,Y]].f''(0)=[X, [X, Y]]\,.

We expand f(t)f(t) using a Taylor series around the point t=0t=0:

f(t)=f(0)+11!f(0)t+12!f(0)t2+13!f(0)t3+f(t)= f(0) + {1\over 1!}f'(0) t + {1\over 2!}f''(0) t^2 + {1\over 3!}f'''(0) t^3 + \cdots

and we get:

f(t)=etXYetX=Y+11![X,Y]t+12![X,[X,Y]]t2+13![X,[X,[X,Y]]]t3+f(t) = e^{tX} Y e^{-tX} = Y + {1\over 1!} [X,Y] t + {1\over 2!} [X,[X,Y]] t^2 + {1\over 3!} [X,[X,[X,Y]]] t^3 + \cdots

We can set t=1t=1:

f(1)=eXYeX=Y+11![X,Y]+12![X,[X,Y]]+13![X,[X,[X,Y]]]+f(1) = e^{X} Y e^{-X} = Y + {1\over 1!} [X,Y] + {1\over 2!} [X,[X,Y]] + {1\over 3!} [X,[X,[X,Y]]] + \cdots

This last equation is called the Campbell identity. We can also write it as:

eXYeX=Y+11![X,Y]+12![X,[X,Y]]+13![X,[X,[X,Y]]]+=e^{-X} Y e^{X} = Y + {1\over 1!} [-X,Y] + {1\over 2!} [-X,[-X,Y]] + {1\over 3!} [-X,[-X,[-X,Y]]] + \cdots =
=Y+11![Y,X]+12![[Y,X],X]+13![[[Y,X],X],X]+= Y + {1\over 1!} [Y,X] + {1\over 2!} [[Y,X],X] + {1\over 3!} [[[Y,X],X],X] + \cdots

Lie Algebra

Lie algebra gg for a Lie group GG that is a topologically closed subgroup of GL(n, R) is defined by:

g={X:exp(tX)GtR}.g = \{X: \exp(t X) \in G \quad \forall t \in \mathbb{R}\}\,.

The following can be proven:

  1. gg is a vector space

  2. X,Yg[X,Y]gX, Y \in g \Rightarrow [X, Y] \in g

  3. gg is the tangent space to GG at II.

  4. exp:gG\exp: g \rightarrow G is locally invertible

Invariant subgroup aa of a group gg: for every aa and gg the following holds:

gag1=a,g a g^{-1} = a'\,,

where aa' is a member of aa.

Invariant abelian subgroup: is an invariant subgroup that is abelian: aa=aaa a' = a' a, or [a,a]=0[a, a'] = 0.

Lie group is:

Typical examples of semisimple groups are the direct products of simple groups.

Using the structure constants [Li,Lj]=CijkLk[L_i, L_j] = C_{ij}{}^k L_k we can define a metric tensor (also called Killing form):

gij=gji=CiklCjlk.g_{ij} = g_{ji} = C_{ik}{}^l C_{jl}{}^k\,.

The regular representation is:

(Li)jk=Cijk.(L_i)_{j}{}^k = C_{ij}{}^k\,.

From this it follows:

gij=CiklCjlk=(Li)kl(Lj)lk=Tr(Li,Lj)(Li,Lj).g_{ij} = C_{ik}{}^l C_{jl}{}^k = (L_i)_{k}{}^l (L_j)_{l}{}^k = \mathrm{Tr} (L_i, L_j) \equiv (L_i, L_j)\,.

This fulfills all properties of a metric, e.g.:

(Li,Lj)=(Lj,Li),(L_i,L_j) = (L_j, L_i)\,,

(Li+Lj,Lk)=(Li,Lk)+(Lj,Lk),(L_i+L_j,L_k) = (L_i, L_k) + (L_j, L_k)\,,

but it is not necessarily positive definite (i.e., (Li,Li)>0(L_i, L_i) > 0).

Theorem: a Lie algebra is semisimple if and only if detgij0\det g_{ij} \neq 0.

The center of a group G is a set of elements that commute with every element of G:

Z(G)={zGgG,zg=gz}.Z(G) = \{z ∈ G \quad | \quad ∀g ∈ G, \quad zg = gz\}\,.

The center is a normal subgroup of G, so the quotient G/Z(G)G / Z(G) is a group.

To construct every other connected Lie group with the same Lie algebra su(2) as SU(2), such as SO(3), we take quotients of SU(2) by discrete subgroups of its center Z(SU(2))={I,I}Z(SU(2))=\{I,−I\}. The possible quotients are:

SU(2)/{I}SU(2),SU(2)/\{I\}≅SU(2),
SU(2)/{I,I}SO(3).SU(2)/\{I,−I\}≅SO(3).

Utilities

Computing exp(θJ3)P1exp(θJ3)\exp(\theta J_3) P_1 \exp(-\theta J_3)

Where JiJ_i are generators of the SO(3) group and P1P_1 are generators of the translation group. We do not know their matrices, we only know their commutators. We use the formula above with t=θt = \theta, X=J3X=J_3 and Y=P1Y=P_1:

eθJ3P1eθJ3=P1+11![J3,P1]θ+12![J3,[J3,P1]]θ2+13![J3,[J3,[J3,P1]]]θ3+.e^{\theta J_3} P_1 e^{-\theta J_3} = P_1 + {1\over 1!} [J_3,P_1] \theta + {1\over 2!} [J_3,[J_3,P_1]] \theta^2 + {1\over 3!} [J_3,[J_3,[J_3,P_1]]] \theta^3 + \cdots\,.

We need to compute the commutators. The Lie algebra commutator is:

[Ji,Pj]=ϵijkPk.[J_i, P_j] = \epsilon_{ijk} P_k\,.

So:

[J3,P1]=ϵ312P2=P2,[J_3, P_1] = \epsilon_{312} P_2 = P_2\,,

[J3,[J3,P1]]=[J3,P2]=ϵ321P1=P1,[J_3, [J_3, P_1]] = [J_3, P_2] = \epsilon_{321} P_1 = -P_1\,,

[J3,[J3,[J3,P1]]]=[J3,P1]=ϵ312P2=P2.[J_3, [J_3, [J_3, P_1]]] = [J_3, -P_1] = -\epsilon_{312} P_2 = -P_2\,.

And higher terms repeat. The commutators alternate between P2P_2 and P1P_1 and their signs:

eθJ3P1eθJ3=P1+11!P2θ+12!(P1)θ2+13!(P2)θ3+14!P1θ4+15!P2θ5+=e^{\theta J_3} P_1 e^{-\theta J_3} = P_1 + {1\over 1!} P_2 \theta + {1\over 2!} (-P_1) \theta^2 + {1\over 3!} (-P_2) \theta^3 + {1\over 4!} P_1 \theta^4 + {1\over 5!} P_2 \theta^5 + \cdots =
=P1(112!θ2+14!θ4+)+P2(11!θ13!θ3+15!θ5+)== P_1 (1 - {1\over 2!} \theta^2 + {1\over 4!} \theta^4 + \cdots) + P_2 ({1\over 1!} \theta - {1\over 3!} \theta^3 + {1\over 5!} \theta^5 + \cdots) =
=P1cos(θ)+P2sin(θ).= P_1 \cos(\theta) + P_2 \sin(\theta)\,.

Any vector operator has the same Lie bracket, for example for JiJ_i we have [Ji,Jj]=ϵijkJk[J_i, J_j] = \epsilon_{ijk} J_k and obtain:

eθJ3J1eθJ3=J1cos(θ)+J2sin(θ).e^{\theta J_3} J_1 e^{-\theta J_3} = J_1 \cos(\theta) + J_2 \sin(\theta)\,.

Computing exp(θJ3)P2exp(θJ3)\exp(\theta J_3) P_2 \exp(-\theta J_3)

Similarly, let’s compute:

eθJ3P2eθJ3=P2+11![J3,P2]θ+12![J3,[J3,P2]]θ2+13![J3,[J3,[J3,P2]]]θ3+=e^{\theta J_3} P_2 e^{-\theta J_3} = P_2 + {1\over 1!} [J_3,P_2] \theta + {1\over 2!} [J_3,[J_3,P_2]] \theta^2 + {1\over 3!} [J_3,[J_3,[J_3,P_2]]] \theta^3 + \cdots =
=P2+11!(P1)θ+12!(P2)θ2+13!P1θ3+14!P2θ4+15!(P1)θ5+== P_2 + {1\over 1!} (-P_1) \theta + {1\over 2!} (-P_2) \theta^2 + {1\over 3!} P_1 \theta^3 + {1\over 4!} P_2 \theta^4 + {1\over 5!} (-P_1) \theta^5 + \cdots =
=P1(θ1!+θ33!θ55!+)+P2(1θ22!+θ44!+)== P_1 \left(-{\theta\over 1!} + {\theta^3\over 3!} - {\theta^5\over 5!} + \cdots\right) + P_2 \left(1 - {\theta^2\over 2!} + {\theta^4\over 4!} + \cdots \right) =
=P1sin(θ)+P2cos(θ).= -P_1 \sin(\theta) + P_2 \cos(\theta)\,.

Computing exp(θJ3)P3exp(θJ3)\exp(\theta J_3) P_3 \exp(-\theta J_3)

Finally:

eθJ3P3eθJ3=P3+11![J3,P3]θ+12![J3,[J3,P3]]θ2+=P3,e^{\theta J_3} P_3 e^{-\theta J_3} = P_3 + {1\over 1!} [J_3,P_3] \theta + {1\over 2!} [J_3,[J_3,P_3]] \theta^2 + \cdots = P_3\,,

since all the commutators are zero.

Summary

We have thus shown that for any operator satisfying [Ji,Pj]=ϵijkPk[J_i, P_j] = \epsilon_{ijk} P_k we get:

eθJ3P1eθJ3=P1cos(θ)+P2sin(θ),eθJ3P2eθJ3=P1sin(θ)+P2cos(θ),eθJ3P3eθJ3=P3.\begin{aligned} e^{\theta J_3} P_1 e^{-\theta J_3} &= P_1 \cos(\theta) + P_2 \sin(\theta)\,, \\ e^{\theta J_3} P_2 e^{-\theta J_3} &= -P_1 \sin(\theta) + P_2 \cos(\theta)\,, \\ e^{\theta J_3} P_3 e^{-\theta J_3} &= P_3\,. \end{aligned}

We can thus write these as a single equation in index notation:

eθJ3PjeθJ3=(R3)ijPi,e^{\theta J_3} P_j e^{-\theta J_3} =(R_3)^i{}_j P_i\,,

where the matrix (R3)ij(R_3)^i{}_j is a rotation around the zz-axis:

(R3(θ))ij=(cosθsinθ0sinθcosθ0001).(R_3(\theta))^i{}_j =\begin{pmatrix} \cos \theta &-\sin \theta & 0 \\ \sin \theta & \cos \theta & 0 \\ 0 & 0 & 1 \\ \end{pmatrix}\,.

We are interested in computing aP=ajPj{\bf a} \cdot {\bf P} = a^j P_j:

eθJ3ajPjeθJ3=(R3)ijajPie^{\theta J_3} a^j P_j e^{-\theta J_3} =(R_3)^i{}_j a^j P_i

and using the below indentity:

eθJ3eajPjeθJ3=eeθJ3ajPjeθJ3=e(R3)ijajPi,e^{\theta J_3} e^{a^j P_j} e^{-\theta J_3} =e^{e^{\theta J_3} a^j P_j e^{-\theta J_3}} =e^{(R_3)^i{}_j a^j P_i}\,,

or:

eθJ3T(a)eθJ3=T(R3(θ)a).e^{\theta J_3} T({\bf a}) e^{-\theta J_3} =T(R_3(\theta) {\bf a})\,.

In an analogous way we can obtain:

eθJ1T(a)eθJ1=T(R1(θ)a),e^{\theta J_1} T({\bf a}) e^{-\theta J_1} =T(R_1(\theta) {\bf a})\,,
eθJ2T(a)eθJ2=T(R2(θ)a).e^{\theta J_2} T({\bf a}) e^{-\theta J_2} =T(R_2(\theta) {\bf a})\,.

Putting it all together we get:

eθ1J1eθ2J2eθ3J3T(a)eθ3J3eθ2J2eθ1J1=e^{\theta_1 J_1} e^{\theta_2 J_2} e^{\theta_3 J_3} T({\bf a}) e^{-\theta_3 J_3} e^{-\theta_2 J_2} e^{-\theta_1 J_1}=
=eθ1J1eθ2J2T(R3(θ3)a)eθ2J2eθ1J1== e^{\theta_1 J_1} e^{\theta_2 J_2} T(R_3(\theta_3) {\bf a}) e^{-\theta_2 J_2} e^{-\theta_1 J_1}=
=eθ1J1T(R2(θ2)R3(θ3)a)eθ1J1= e^{\theta_1 J_1} T(R_2(\theta_2) R_3(\theta_3) {\bf a}) e^{-\theta_1 J_1}
=T(R1(θ1)R2(θ2)R3(θ3)a).= T(R_1(\theta_1) R_2(\theta_2) R_3(\theta_3) {\bf a}) \,.

In other words, we can see how a general element of the rotation group eθ1J1eθ2J2eθ3J3e^{\theta_1 J_1} e^{\theta_2 J_2} e^{\theta_3 J_3} acts on a general element of the translation group T(a)T({\bf a}). The result is another translation T(R1(θ1)R2(θ2)R3(θ3)a)T(R_1(\theta_1) R_2(\theta_2) R_3(\theta_3) {\bf a}) and the translation vector a{\bf a} gets rotated with the three rotation matrices aR1(θ1)R2(θ2)R3(θ3)a{\bf a} \rightarrow R_1(\theta_1) R_2(\theta_2) R_3(\theta_3) {\bf a}.

For T=J

As a special case, we’ll use Pi=JiP_i = J_i, and we’ll use n\bf n instead of a\bf a:

eθJ3njJjeθJ3=(R3)ijnjJie^{\theta J_3} n^j J_j e^{-\theta J_3} =(R_3)^i{}_j n^j J_i

In an analogous way we can obtain:

eθJ1njJjeθJ1=(R1)ijnjJie^{\theta J_1} n^j J_j e^{-\theta J_1} =(R_1)^i{}_j n^j J_i
eθJ2njJjeθJ2=(R2)ijnjJie^{\theta J_2} n^j J_j e^{-\theta J_2} =(R_2)^i{}_j n^j J_i

Putting it all together we get:

eθ1J1eθ2J2eθ3J3njJjeθ3J3eθ2J2eθ1J1=e^{\theta_1 J_1} e^{\theta_2 J_2} e^{\theta_3 J_3} n^j J_j e^{-\theta_3 J_3} e^{-\theta_2 J_2} e^{-\theta_1 J_1}=
=eθ1J1eθ2J2(R3)ijnjJieθ2J2eθ1J1== e^{\theta_1 J_1} e^{\theta_2 J_2} (R_3)^i{}_j n^j J_i e^{-\theta_2 J_2} e^{-\theta_1 J_1}=
=eθ1J1(R2)ij(R3)jknkJieθ1J1== e^{\theta_1 J_1} (R_2)^i{}_j (R_3)^j{}_k n^k J_i e^{-\theta_1 J_1} =
=(R1)ij(R2)jk(R3)klnlJi.= (R_1)^i{}_j (R_2)^j{}_k (R_3)^k{}_l n^l J_i \,.

This shows how to “encode” the 3D vector n\bf n into any jj representation. The RHS is encoding the rotated vector, the LHS is encoding the original vector and then transforming it using the “sandwich” formula. Thus showing that the “sandwich” transformation is equivalent to a rotation. It is thus showing how any SU(2) transformation corresponds to SO(3) transformation, the “minus” cancels out in the sandwich, so 360 degree rotation transforms n\bf n into itself. To obtain the 360 rotation equal to -1 we need to act with just the left part of the sandwich on some “spinor”, and it would transform it to -1.

Thus we have proven:

U(θ1,θ2,θ3) aJ U(θ1,θ2,θ3)1=(R(θ1,θ2,θ3)a)J,U(\theta_1, \theta_2, \theta_3) \ {\bf a} \cdot {\bf J}\ U(\theta_1, \theta_2, \theta_3)^{-1} = (R(\theta_1, \theta_2, \theta_3) {\bf a}) \cdot {\bf J}\,,

where:

U(θ1,θ2,θ3)=eθ1J1eθ2J2eθ3J3,U(\theta_1, \theta_2, \theta_3) =e^{\theta_1 J_1} e^{\theta_2 J_2} e^{\theta_3 J_3}\,,
R(θ1,θ2,θ3)=R1(θ1)R2(θ2)R3(θ3).R(\theta_1, \theta_2, \theta_3) = R_1(\theta_1) R_2(\theta_2) R_3(\theta_3) \,.

Spin 1/2

For spin j=1/2 we have:

Ji=i12σiJ_i = -i {1\over 2} \sigma_i

Now we would like to show that:

U(θ1,θ2,θ3)=eiθ12σ1eiθ22σ2eiθ32σ3=U(n,θ)=eiθ2(nσ),U(\theta_1, \theta_2, \theta_3) =e^{-i{\theta_1\over2} \sigma_1} e^{-i{\theta_2\over2} \sigma_2} e^{-i{\theta_3\over2} \sigma_3} =U({\bf n}, \theta) = e^{-i {\theta\over2}\left({\bf n}\cdot {\bf\sigma} \right)}\,,

where n{\bf n} is a unit vector representing the axis of rotation and θ\theta is the rotation around this axis. This follows from the fact that the LHS is a general element of the SU(2) group, and so is the RHS. Constructively, this is proven in a separate file.

TODO: We also want to prove that the three rotation matrices can be written as just one matrix:

R(θ1,θ2,θ3)=R1(θ1)R2(θ2)R3(θ3)=R(n,θ).R(\theta_1, \theta_2, \theta_3) = R_1(\theta_1) R_2(\theta_2) R_3(\theta_3) = R({\bf n}, \theta)\,.

Then we get:

U(n,θ) aσ U(n,θ)1=(R(n,θ)a)σ,U({\bf n}, \theta) \ {\bf a} \cdot {\bf \sigma}\ U({\bf n}, \theta)^{-1} = (R({\bf n}, \theta) {\bf a}) \cdot {\bf \sigma}\,,

Here UU can be expanded as:

U(n,θ)=eiθ2(nσ)=cos(θ2)isin(θ2)(nσ).U({\bf n}, \theta) = e^{-i {\theta\over2}\left({\bf n}\cdot {\bf\sigma} \right)} = \cos\left({\theta\over 2}\right) - i \sin\left({\theta\over2}\right) \left({\bf n}\cdot {\bf\sigma} \right)\,.

We can now identify this relation with quaternions as follows:

U(n,θ)=cos(θ2)+sin(θ2)(n1(iσ1)+n2(iσ2)+n3(iσ3)),U({\bf n}, \theta) = \cos\left({\theta\over 2}\right) + \sin\left({\theta\over2}\right) \left(n^1 (-i\sigma_1) + n^2 (-i\sigma_2) + n^3 (-i\sigma_3) \right)\,,
q=cos(θ2)+sin(θ2)(n1i+n2j+n3k).q = \cos\left({\theta\over 2}\right) + \sin\left({\theta\over2}\right) \left(n^1 i + n^2 j + n^3 k \right)\,.

And the rotation equation is:

q(a1i+a2j+a3k)q1=(R(q)a)(i,j,k)q (a^1 i + a^2 j + a^3 k) q^{-1} = (R(q){\bf a}) \cdot (i, j, k)

In other words, the LHS is encoding the vector a{\bf a} into a quaternion as a1i+a2j+a3ka^1 i + a^2 j + a^3 k, and the RHS gives us the rotated quaternion vector that we can decode back.

Identity geXg1=egXg1g e^X g^{-1} = e^{g X g^{-1}}

The identity

geXg1=egXg1g e^X g^{-1} = e^{g X g^{-1}}

can be proven using:

geXg1=g(n=0Xnn!)g1=n=0gXng1n!=n=0(gXg1)nn!=egXg1.g e^X g^{-1}=g \left( \sum_{n=0}^\infty {X^n\over n!} \right) g^{-1} =\sum_{n=0}^\infty {g X^n g^{-1} \over n!} =\sum_{n=0}^\infty {(g X g^{-1})^n \over n!} =e^{g X g^{-1}}\,.

Computing exp(θ3J3)exp(θ1P1)exp(θ3J3)\exp(\theta_3 J_3) \exp(\theta_1 P_1) \exp(-\theta_3 J_3)

Using the identity geXg1=egXg1g e^X g^{-1} = e^{g X g^{-1}} we get:

eθ3J3eθ1P1eθ3J3=exp(eθ3J3(θ1P1)eθ3J3)=eθ1(P1cosθ3+P2sinθ3).e^{\theta_3 J_3} e^{\theta_1 P_1} e^{-\theta_3 J_3} =\exp\left(e^{\theta_3 J_3} (\theta_1 P_1) e^{-\theta_3 J_3}\right) =e^{\theta_1 \left(P_1 \cos\theta_3 + P_2 \sin\theta_3\right)} \,.

Rodrigues’ Rotation Formula

If v\mathbf{v} is a vector in R^3 and k\mathbf{k} is a unit vector (axis of rotation) and θ\theta the angle of rotation, then the vector v\mathbf{v} gets rotated to vrot\mathbf{v}_\mathrm{rot}.

We define a matrix K\mathbf{K} using:

Kij=ϵijkkk=(0k3k2k30k1k2k10)\mathbf{K}_{ij} = -\epsilon_{ijk} k^k = \begin{pmatrix} 0 & -k^3 & k^2 \\ k^3 & 0 & -k^1 \\ -k^2 & k^1 & 0 \end{pmatrix}

Note that k2=k12+k22+k33=1|\mathbf{k}|^2 = k_1^2 + k_2^2 + k_3^3 = 1. We get:

(Kv)i=Kijvj=ϵijkkkvj=ϵikjkkvj=(k×v)i.(\mathbf{K}\mathbf{v})_i =\mathbf{K}_{ij} v^j = -\epsilon_{ijk} k^k v^j = \epsilon_{ikj} k^k v^j = (\mathbf{k}\times\mathbf{v})_i\,.

We will also need:

(K2)ij=KikKkj=(ϵiklkl)(ϵkjmkm)=(δjlδmiδjiδml)klkm=kikjδij,\left(\mathbf{K}^2\right)_{ij} =K_i{}^k K_{kj} =(-\epsilon_i{}^k{}_l\, k^l)(-\epsilon_{kjm} k^m) =(\delta_{jl} \delta_{mi}-\delta_{ji}\delta_{ml}) k^l k^m =k_i k_j - \delta_{ij}\,,

where we used δmlklkm=1\delta_{ml}k^l k^m=1.

Now we compute the characteristic polynomial:

P(t)=det(KtI)=t3t(k12+k22+k32)=t3t,P(t) = \det(\mathbf{K}-t\mathbf{I})=-t^3-t(k_1^2+k_2^2+k_3^2) = -t^3-t\,,
P(K)=0,P(\mathbf{K}) = 0\,,
K3K=0,-\mathbf{K}^3-\mathbf{K} = 0\,,
K3=K.\mathbf{K}^3=-\mathbf{K}\,.

From the last equation we get K4=K2\mathbf{K}^4=-\mathbf{K}^2 and K5=K3=K\mathbf{K}^5=-\mathbf{K}^3=\mathbf{K} and so on.

Then the rotation matrix R\mathbf{R} is:

R=eθK=k=0(θK)kk!=I+θK+(θK)22!+(θK)33!+=\mathbf{R} = e^{\theta \mathbf{K}} = \sum_{k=0}^\infty {(\theta\mathbf{K})^k \over k!} = \mathbf{I} + \theta\mathbf{K} + {(\theta\mathbf{K})^2 \over 2!} + {(\theta\mathbf{K})^3 \over 3!} + \cdots =
=I+K(θθ33!+θ55!)+K2(θ22!θ44!+θ66!)== \mathbf{I} + \mathbf{K}\left(\theta - {\theta^3\over3!} + {\theta^5\over5!} - \cdots\right) + \mathbf{K}^2\left({\theta^2\over2!} - {\theta^4\over4!} + {\theta^6\over6!} - \cdots\right) =
=I+Ksinθ+K2(1cosθ).= \mathbf{I} + \mathbf{K}\sin\theta + \mathbf{K}^2(1-\cos\theta)\,.

So we got:

R=I+Ksinθ+K2(1cosθ).\mathbf{R} = \mathbf{I} +\mathbf{K} \sin\theta +\mathbf{K}^2 (1-\cos\theta)\,.

Now we can write:

vrot=Rv\mathbf{v}_\mathrm{rot} = \mathbf{R} \mathbf{v}

and get:

vrot=v+Kvsinθ+K2v(1cosθ).\mathbf{v}_\mathrm{rot} = \mathbf{v} +\mathbf{K}\mathbf{v} \sin\theta +\mathbf{K}^2\mathbf{v} (1-\cos\theta)\,.

We use:

Kv=k×v,\mathbf{K}\mathbf{v} = \mathbf{k}\times\mathbf{v}\,,
K2v=K(Kv)=k×(k×v)=k(kv)(kk)v=k(kv)v\mathbf{K}^2\mathbf{v} = \mathbf{K}(\mathbf{K}\mathbf{v}) = \mathbf{k}\times(\mathbf{k}\times\mathbf{v}) = \mathbf{k}(\mathbf{k}\cdot\mathbf{v}) - (\mathbf{k}\cdot\mathbf{k})\mathbf{v} = \mathbf{k}(\mathbf{k}\cdot\mathbf{v}) - \mathbf{v}

and get:

vrot=vcosθ+(k×v)sinθ+k(kv)(1cosθ).\mathbf{v}_\mathrm{rot} = \mathbf{v} \cos\theta +(\mathbf{k}\times\mathbf{v}) \sin\theta +\mathbf{k}(\mathbf{k}\cdot\mathbf{v}) (1-\cos\theta)\,.

This is called the Rodrigues’ formula. In components:

Rij=Iij+Kijsinθ+(K2)ij(1cosθ).\mathbf{R}_{ij} = \mathbf{I}_{ij} +\mathbf{K}_{ij} \sin\theta +\left(\mathbf{K}^2\right)_{ij} (1-\cos\theta)\,.

We can simplify:

Rij=δijϵijknksinθ+(ninjδij)(1cosθ),R_{ij} = \delta_{ij} -\epsilon_{ij}{}^k n_k \sin\theta +(n_i n_j - \delta_{ij}) (1-\cos\theta)\,,
Rij=cosθδij+(1cosθ)ninjsinθϵijknk.R_{ij}=\cos\theta\delta_{ij}+(1-\cos\theta)n_i n_j -\sin\theta\,\epsilon_{ij}{}^k n_k\,.