3D rotation group
In mechanics and geometry, the 3D rotation group, often denoted SO(3), is the group of all rotations about the origin of threedimensional Euclidean space [math]\displaystyle{ \R^3 }[/math] under the operation of composition.^{[1]}
By definition, a rotation about the origin is a transformation that preserves the origin, Euclidean distance (so it is an isometry), and orientation (i.e., handedness of space). Composing two rotations results in another rotation, every rotation has a unique inverse rotation, and the identity map satisfies the definition of a rotation. Owing to the above properties (along composite rotations' associative property), the set of all rotations is a group under composition.
Every nontrivial rotation is determined by its axis of rotation (a line through the origin) and its angle of rotation. Rotations are not commutative (for example, rotating R 90° in the xy plane followed by S 90° in the yz plane is not the same as S followed by R), making the 3D rotation group a nonabelian group. Moreover, the rotation group has a natural structure as a manifold for which the group operations are smoothly differentiable, so it is in fact a Lie group. It is compact and has dimension 3.
Rotations are linear transformations of [math]\displaystyle{ \R^3 }[/math] and can therefore be represented by matrices once a basis of [math]\displaystyle{ \R^3 }[/math] has been chosen. Specifically, if we choose an orthonormal basis of [math]\displaystyle{ \R^3 }[/math], every rotation is described by an orthogonal 3 × 3 matrix (i.e., a 3 × 3 matrix with real entries which, when multiplied by its transpose, results in the identity matrix) with determinant 1. The group SO(3) can therefore be identified with the group of these matrices under matrix multiplication. These matrices are known as "special orthogonal matrices", explaining the notation SO(3).
The group SO(3) is used to describe the possible rotational symmetries of an object, as well as the possible orientations of an object in space. Its representations are important in physics, where they give rise to the elementary particles of integer spin.
Length and angle
Besides just preserving length, rotations also preserve the angles between vectors. This follows from the fact that the standard dot product between two vectors u and v can be written purely in terms of length: [math]\displaystyle{ \mathbf{u} \cdot \mathbf{v} = \frac{1}{2} \left(\\mathbf{u} + \mathbf{v}\^2  \\mathbf{u}\^2  \\mathbf{v}\^2\right). }[/math]
It follows that every lengthpreserving linear transformation in [math]\displaystyle{ \R^3 }[/math] preserves the dot product, and thus the angle between vectors. Rotations are often defined as linear transformations that preserve the inner product on [math]\displaystyle{ \R^3 }[/math], which is equivalent to requiring them to preserve length. See classical group for a treatment of this more general approach, where SO(3) appears as a special case.
Orthogonal and rotation matrices
Every rotation maps an orthonormal basis of [math]\displaystyle{ \R^3 }[/math] to another orthonormal basis. Like any linear transformation of finitedimensional vector spaces, a rotation can always be represented by a matrix. Let R be a given rotation. With respect to the standard basis e_{1}, e_{2}, e_{3} of [math]\displaystyle{ \R^3 }[/math] the columns of R are given by (Re_{1}, Re_{2}, Re_{3}). Since the standard basis is orthonormal, and since R preserves angles and length, the columns of R form another orthonormal basis. This orthonormality condition can be expressed in the form
 [math]\displaystyle{ R^\mathsf{T}R = RR^\mathsf{T} = I, }[/math]
where R^{T} denotes the transpose of R and I is the 3 × 3 identity matrix. Matrices for which this property holds are called orthogonal matrices. The group of all 3 × 3 orthogonal matrices is denoted O(3), and consists of all proper and improper rotations.
In addition to preserving length, proper rotations must also preserve orientation. A matrix will preserve or reverse orientation according to whether the determinant of the matrix is positive or negative. For an orthogonal matrix R, note that det R^{T} = det R implies (det R)^{2} = 1, so that det R = ±1. The subgroup of orthogonal matrices with determinant +1 is called the special orthogonal group, denoted SO(3).
Thus every rotation can be represented uniquely by an orthogonal matrix with unit determinant. Moreover, since composition of rotations corresponds to matrix multiplication, the rotation group is isomorphic to the special orthogonal group SO(3).
Improper rotations correspond to orthogonal matrices with determinant −1, and they do not form a group because the product of two improper rotations is a proper rotation.
Group structure
The rotation group is a group under function composition (or equivalently the product of linear transformations). It is a subgroup of the general linear group consisting of all invertible linear transformations of the real 3space [math]\displaystyle{ \R^3 }[/math].^{[2]}
Furthermore, the rotation group is nonabelian. That is, the order in which rotations are composed makes a difference. For example, a quarter turn around the positive xaxis followed by a quarter turn around the positive yaxis is a different rotation than the one obtained by first rotating around y and then x.
The orthogonal group, consisting of all proper and improper rotations, is generated by reflections. Every proper rotation is the composition of two reflections, a special case of the Cartan–Dieudonné theorem.
Complete classification of finite subgroups
The finite subgroups of [math]\displaystyle{ \mathrm{SO}(3) }[/math] are completely classified.^{[3]}
Every finite subgroup is isomorphic to either an element of one of two countably infinite families of planar isometries: the cyclic groups [math]\displaystyle{ C_n }[/math] or the dihedral groups [math]\displaystyle{ D_{2n} }[/math], or to one of three other groups: the tetrahedral group [math]\displaystyle{ \cong A_4 }[/math], the octahedral group [math]\displaystyle{ \cong S_4 }[/math], or the icosahedral group [math]\displaystyle{ \cong A_5 }[/math].
Axis of rotation
Every nontrivial proper rotation in 3 dimensions fixes a unique 1dimensional linear subspace of [math]\displaystyle{ \R^3 }[/math] which is called the axis of rotation (this is Euler's rotation theorem). Each such rotation acts as an ordinary 2dimensional rotation in the plane orthogonal to this axis. Since every 2dimensional rotation can be represented by an angle φ, an arbitrary 3dimensional rotation can be specified by an axis of rotation together with an angle of rotation about this axis. (Technically, one needs to specify an orientation for the axis and whether the rotation is taken to be clockwise or counterclockwise with respect to this orientation).
For example, counterclockwise rotation about the positive zaxis by angle φ is given by
 [math]\displaystyle{ R_z(\phi) = \begin{bmatrix}\cos\phi & \sin\phi & 0 \\ \sin\phi & \cos\phi & 0 \\ 0 & 0 & 1\end{bmatrix}. }[/math]
Given a unit vector n in [math]\displaystyle{ \R^3 }[/math] and an angle φ, let R(φ, n) represent a counterclockwise rotation about the axis through n (with orientation determined by n). Then
 R(0, n) is the identity transformation for any n
 R(φ, n) = R(−φ, −n)
 R(π + φ, n) = R(π − φ, −n).
Using these properties one can show that any rotation can be represented by a unique angle φ in the range 0 ≤ φ ≤ π and a unit vector n such that
 n is arbitrary if φ = 0
 n is unique if 0 < φ < π
 n is unique up to a sign if φ = π (that is, the rotations R(π, ±n) are identical).
In the next section, this representation of rotations is used to identify SO(3) topologically with threedimensional real projective space.
Topology
The Lie group SO(3) is diffeomorphic to the real projective space [math]\displaystyle{ \mathbb{P}^3(\R). }[/math]^{[4]}
Consider the solid ball in [math]\displaystyle{ \R^3 }[/math] of radius π (that is, all points of [math]\displaystyle{ \R^3 }[/math] of distance π or less from the origin). Given the above, for every point in this ball there is a rotation, with axis through the point and the origin, and rotation angle equal to the distance of the point from the origin. The identity rotation corresponds to the point at the center of the ball. Rotation through angles between 0 and −π correspond to the point on the same axis and distance from the origin but on the opposite side of the origin. The one remaining issue is that the two rotations through π and through −π are the same. So we identify (or "glue together") antipodal points on the surface of the ball. After this identification, we arrive at a topological space homeomorphic to the rotation group.
Indeed, the ball with antipodal surface points identified is a smooth manifold, and this manifold is diffeomorphic to the rotation group. It is also diffeomorphic to the real 3dimensional projective space [math]\displaystyle{ \mathbb{P}^3(\R), }[/math] so the latter can also serve as a topological model for the rotation group.
These identifications illustrate that SO(3) is connected but not simply connected. As to the latter, in the ball with antipodal surface points identified, consider the path running from the "north pole" straight through the interior down to the south pole. This is a closed loop, since the north pole and the south pole are identified. This loop cannot be shrunk to a point, since no matter how it is deformed, the start and end point have to remain antipodal, or else the loop will "break open". In terms of rotations, this loop represents a continuous sequence of rotations about the zaxis starting (by example) at the identity (center of the ball), through the south pole, jumping to the north pole and ending again at the identity rotation (i.e., a series of rotation through an angle φ where φ runs from 0 to 2π).
Surprisingly, running through the path twice, i.e., running from the north pole down to the south pole, jumping back to the north pole (using the fact that north and south poles are identified), and then again running from the north pole down to the south pole, so that φ runs from 0 to 4π, gives a closed loop which can be shrunk to a single point: first move the paths continuously to the ball's surface, still connecting north pole to south pole twice. The second path can then be mirrored over to the antipodal side without changing the path at all. Now we have an ordinary closed loop on the surface of the ball, connecting the north pole to itself along a great circle. This circle can be shrunk to the north pole without problems. The plate trick and similar tricks demonstrate this practically.
The same argument can be performed in general, and it shows that the fundamental group of SO(3) is the cyclic group of order 2 (a fundamental group with two elements). In physics applications, the nontriviality (more than one element) of the fundamental group allows for the existence of objects known as spinors, and is an important tool in the development of the spin–statistics theorem.
The universal cover of SO(3) is a Lie group called Spin(3). The group Spin(3) is isomorphic to the special unitary group SU(2); it is also diffeomorphic to the unit 3sphere S^{3} and can be understood as the group of versors (quaternions with absolute value 1). The connection between quaternions and rotations, commonly exploited in computer graphics, is explained in quaternions and spatial rotations. The map from S^{3} onto SO(3) that identifies antipodal points of S^{3} is a surjective homomorphism of Lie groups, with kernel {±1}. Topologically, this map is a twotoone covering map. (See the plate trick.)
Connection between SO(3) and SU(2)
In this section, we give two different constructions of a twotoone and surjective homomorphism of SU(2) onto SO(3).
Using quaternions of unit norm
The group SU(2) is isomorphic to the quaternions of unit norm via a map given by^{[5]} [math]\displaystyle{ q = a\mathbf{1} + b\mathbf{i} + c\mathbf{j} + d\mathbf{k} = \alpha + \beta \mathbf{j} \leftrightarrow \begin{bmatrix}\alpha & \overline \beta \\ \beta & \overline \alpha\end{bmatrix} = U }[/math] restricted to [math]\displaystyle{ a^2+ b^2 + c^2 + d^2 = \alpha^2 +\beta^2 = 1 }[/math] where [math]\displaystyle{ q \in \mathbb{H} }[/math], [math]\displaystyle{ a, b, c, d \in \R }[/math], [math]\displaystyle{ U \in \operatorname{SU}(2) }[/math], and [math]\displaystyle{ \alpha = a+bi \in\mathbb{C} }[/math], [math]\displaystyle{ \beta = c+di \in \mathbb{C} }[/math].
Let us now identify [math]\displaystyle{ \R^3 }[/math] with the span of [math]\displaystyle{ \mathbf{i},\mathbf{j},\mathbf{k} }[/math]. One can then verify that if [math]\displaystyle{ v }[/math] is in [math]\displaystyle{ \R^3 }[/math] and [math]\displaystyle{ q }[/math] is a unit quaternion, then [math]\displaystyle{ qvq^{1}\in \R^3. }[/math]
Furthermore, the map [math]\displaystyle{ v\mapsto qvq^{1} }[/math] is a rotation of [math]\displaystyle{ \R^3. }[/math] Moreover, [math]\displaystyle{ (q)v(q)^{1} }[/math] is the same as [math]\displaystyle{ qvq^{1} }[/math]. This means that there is a 2:1 homomorphism from quaternions of unit norm to the 3D rotation group SO(3).
One can work this homomorphism out explicitly: the unit quaternion, q, with [math]\displaystyle{ \begin{align} q &= w + x\mathbf{i} + y\mathbf{j} + z\mathbf{k} , \\ 1 &= w^2 + x^2 + y^2 + z^2 , \end{align} }[/math] is mapped to the rotation matrix [math]\displaystyle{ Q = \begin{bmatrix} 1  2 y^2  2 z^2 & 2 x y  2 z w & 2 x z + 2 y w \\ 2 x y + 2 z w & 1  2 x^2  2 z^2 & 2 y z  2 x w \\ 2 x z  2 y w & 2 y z + 2 x w & 1  2 x^2  2 y^2 \end{bmatrix}. }[/math]
This is a rotation around the vector (x, y, z) by an angle 2θ, where cos θ = w and sin θ = (x, y, z). The proper sign for sin θ is implied, once the signs of the axis components are fixed. The 2:1nature is apparent since both q and −q map to the same Q.
Using Möbius transformations
The general reference for this section is (Gelfand Minlos). The points P on the sphere
 [math]\displaystyle{ \mathbf{S} = \left \{(x,y,z)\in\R^3: x^2 +y^2 +z^2 = \frac{1}{4} \right \} }[/math]
can, barring the north pole N, be put into onetoone bijection with points S(P) = P' on the plane M defined by z = −1/2, see figure. The map S is called stereographic projection.
Let the coordinates on M be (ξ, η). The line L passing through N and P can be parametrized as
 [math]\displaystyle{ L(t) = N + t(N  P) = \left(0,0,\frac{1}{2}\right) + t \left ( \left(0,0,\frac{1}{2}\right)  (x, y, z) \right ), \quad t\in \R. }[/math]
Demanding that the zcoordinate of [math]\displaystyle{ L(t_0) }[/math] equals −1/2, one finds
 [math]\displaystyle{ t_0 = \frac1{z\frac12}. }[/math]
We have [math]\displaystyle{ L(t_0)=(\xi,\eta,1/2). }[/math] Hence the map
 [math]\displaystyle{ \begin{cases} S:\mathbf{S} \to M \\ P = (x,y,z) \longmapsto P'= (\xi, \eta) = \left(\frac{x}{\frac{1}{2}  z}, \frac{y}{\frac{1}{2}  z}\right) \equiv \zeta = \xi + i\eta \end{cases} }[/math]
where, for later convenience, the plane M is identified with the complex plane [math]\displaystyle{ \Complex. }[/math]
For the inverse, write L as
 [math]\displaystyle{ L = N + s(P'N) = \left(0,0,\frac{1}{2}\right) + s\left( \left(\xi, \eta, \frac{1}{2}\right)  \left(0,0,\frac{1}{2}\right)\right), }[/math]
and demand x^{2} + y^{2} + z^{2} = 1/4 to find s = 1/1 + ξ^{2} + η^{2} and thus
 [math]\displaystyle{ \begin{cases} S^{1}:M \to \mathbf{S} \\ P'= (\xi, \eta) \longmapsto P = (x,y,z) = \left(\frac{\xi}{1 + \xi^2 + \eta^2}, \frac{\eta}{1 + \xi^2 + \eta^2}, \frac{1 + \xi^2 + \eta^2}{2 + 2\xi^2 + 2\eta^2}\right) \end{cases} }[/math]
If g ∈ SO(3) is a rotation, then it will take points on S to points on S by its standard action Π_{s}(g) on the embedding space [math]\displaystyle{ \R^3. }[/math] By composing this action with S one obtains a transformation S ∘ Π_{s}(g) ∘ S^{−1} of M,
 [math]\displaystyle{ \zeta=P' \longmapsto P \longmapsto \Pi_s(g)P = gP \longmapsto S(gP) \equiv \Pi_u(g)\zeta = \zeta'. }[/math]
Thus Π_{u}(g) is a transformation of [math]\displaystyle{ \Complex }[/math] associated to the transformation Π_{s}(g) of [math]\displaystyle{ \R^3 }[/math].
It turns out that g ∈ SO(3) represented in this way by Π_{u}(g) can be expressed as a matrix Π_{u}(g) ∈ SU(2) (where the notation is recycled to use the same name for the matrix as for the transformation of [math]\displaystyle{ \Complex }[/math] it represents). To identify this matrix, consider first a rotation g_{φ} about the zaxis through an angle φ,
 [math]\displaystyle{ \begin{align} x' &= x\cos \phi  y \sin \phi,\\ y' &= x\sin \phi + y \cos \phi,\\ z' &= z. \end{align} }[/math]
Hence
 [math]\displaystyle{ \zeta' = \frac{x' + iy'}{\frac{1}{2}  z'} = \frac{e^{i\phi}(x + iy)}{\frac{1}{2}  z} = e^{i\phi}\zeta = \frac{e^{\frac{i\phi}{2}} \zeta + 0 }{0 \zeta + e^{\frac{i\phi}{2}}}, }[/math]
which, unsurprisingly, is a rotation in the complex plane. In an analogous way, if g_{θ} is a rotation about the xaxis through an angle θ, then
 [math]\displaystyle{ w' = e^{i\theta}w, \quad w = \frac{y + iz}{\frac{1}{2}  x}, }[/math]
which, after a little algebra, becomes
 [math]\displaystyle{ \zeta' = \frac{\cos \frac{\theta}{2}\zeta +i\sin \frac{\theta}{2} }{i \sin\frac{\theta}{2}\zeta + \cos\frac{\theta}{2}}. }[/math]
These two rotations, [math]\displaystyle{ g_{\phi}, g_{\theta}, }[/math] thus correspond to bilinear transforms of R^{2} ≃ C ≃ M, namely, they are examples of Möbius transformations.
A general Möbius transformation is given by
 [math]\displaystyle{ \zeta' = \frac{\alpha \zeta + \beta}{\gamma \zeta + \delta}, \quad \alpha\delta  \beta\gamma \ne 0. }[/math]
The rotations, [math]\displaystyle{ g_{\phi}, g_{\theta} }[/math] generate all of SO(3) and the composition rules of the Möbius transformations show that any composition of [math]\displaystyle{ g_{\phi}, g_{\theta} }[/math] translates to the corresponding composition of Möbius transformations. The Möbius transformations can be represented by matrices
 [math]\displaystyle{ \begin{pmatrix}\alpha & \beta\\ \gamma & \delta\end{pmatrix}, \qquad \alpha\delta  \beta\gamma = 1, }[/math]
since a common factor of α, β, γ, δ cancels.
For the same reason, the matrix is not uniquely defined since multiplication by −I has no effect on either the determinant or the Möbius transformation. The composition law of Möbius transformations follow that of the corresponding matrices. The conclusion is that each Möbius transformation corresponds to two matrices g, −g ∈ SL(2, C).
Using this correspondence one may write
 [math]\displaystyle{ \begin{align} \Pi_u(g_\phi) &= \Pi_u\left[\begin{pmatrix} \cos \phi & \sin \phi & 0\\ \sin \phi & \cos \phi & 0\\ 0 & 0 & 1 \end{pmatrix}\right] = \pm \begin{pmatrix} e^{i\frac{\phi}{2}} & 0\\ 0 & e^{i\frac{\phi}{2}} \end{pmatrix},\\ \Pi_u(g_\theta) &= \Pi_u\left[\begin{pmatrix} 1 & 0 & 0\\ 0 & \cos \theta & \sin \theta\\ 0 & \sin \theta & \cos \theta \end{pmatrix}\right] = \pm \begin{pmatrix} \cos\frac{\theta}{2} & i\sin\frac{\theta}{2}\\ i\sin\frac{\theta}{2} & \cos\frac{\theta}{2} \end{pmatrix}. \end{align} }[/math]
These matrices are unitary and thus Π_{u}(SO(3)) ⊂ SU(2) ⊂ SL(2, C). In terms of Euler angles^{[nb 1]} one finds for a general rotation

[math]\displaystyle{ \begin{align} g(\phi, \theta, \psi) = g_\phi g_\theta g_\psi &= \begin{pmatrix} \cos \phi & \sin \phi & 0\\ \sin \phi & \cos \phi & 0\\ 0 & 0 & 1 \end{pmatrix}\begin{pmatrix} 1 & 0 & 0\\ 0 & \cos \theta & \sin \theta\\ 0 & \sin \theta & \cos \theta \end{pmatrix}\begin{pmatrix} \cos \psi & \sin \psi & 0\\ \sin \psi & \cos \psi & 0\\ 0 & 0 & 1 \end{pmatrix}\\ &= \begin{pmatrix} \cos\phi\cos\psi  \cos\theta\sin\phi\sin\psi & \cos\phi\sin\psi  \cos\theta\sin\phi\cos\psi & \sin\phi\sin\theta \\ \sin\phi\cos\psi + \cos\theta\cos\phi\sin\psi & \sin\phi\sin\psi + \cos\theta\cos\phi\cos\psi & \cos\phi\sin\theta\\ \sin\psi\sin\theta & \cos\psi\sin\theta & \cos\theta \end{pmatrix}, \end{align} }[/math]
(
)
one has^{[6]}

[math]\displaystyle{ \begin{align}\Pi_u(g(\phi, \theta, \psi)) &= \pm \begin{pmatrix} e^{i\frac{\phi}{2}} & 0\\ 0 & e^{i\frac{\phi}{2}} \end{pmatrix} \begin{pmatrix} \cos\frac{\theta}{2} & i\sin\frac{\theta}{2}\\ i\sin\frac{\theta}{2} & \cos\frac{\theta}{2} \end{pmatrix} \begin{pmatrix} e^{i\frac{\psi}{2}} & 0\\ 0 & e^{i\frac{\psi}{2}} \end{pmatrix}\\ &= \pm \begin{pmatrix} \cos\frac{\theta}{2}e^{i\frac{\phi + \psi}{2}} & i\sin\frac{\theta}{2}e^{i\frac{\phi  \psi}{2}}\\ i\sin\frac{\theta}{2}e^{i\frac{\phi  \psi}{2}} & \cos\frac{\theta}{2}e^{i\frac{\phi + \psi}{2}} \end{pmatrix}. \end{align} }[/math]
(
)
For the converse, consider a general matrix
 [math]\displaystyle{ \pm\Pi_u(g_{\alpha,\beta}) = \pm\begin{pmatrix} \alpha & \beta\\ \overline{\beta} & \overline{\alpha} \end{pmatrix} \in \operatorname{SU}(2). }[/math]
Make the substitutions
 [math]\displaystyle{ \begin{align} \cos\frac{\theta}{2} &= \alpha, & \sin\frac{\theta}{2} &= \beta, & (0 \le \theta \le \pi),\\ \frac{\phi + \psi}{2} &= \arg \alpha, & \frac{\psi  \phi}{2} &= \arg \beta. & \end{align} }[/math]
With the substitutions, Π(g_{α, β}) assumes the form of the right hand side (RHS) of (2), which corresponds under Π_{u} to a matrix on the form of the RHS of (1) with the same φ, θ, ψ. In terms of the complex parameters α, β,
 [math]\displaystyle{ g_{\alpha,\beta} = \begin{pmatrix} \frac{1}{2}\left( \alpha^2  \beta^2 + \overline{\alpha^2}  \overline{\beta^2}\right) & \frac{i}{2}\left(\alpha^2  \beta^2 + \overline{\alpha^2} + \overline{\beta^2}\right) & \alpha\beta  \overline{\alpha}\overline{\beta}\\ \frac{i}{2}\left(\alpha^2  \beta^2  \overline{\alpha^2} + \overline{\beta^2}\right) & \frac{1}{2}\left(\alpha^2 + \beta^2 + \overline{\alpha^2} + \overline{\beta^2}\right) & i\left(+\alpha\beta  \overline{\alpha}\overline{\beta}\right)\\ \alpha\overline{\beta} + \overline{\alpha}\beta & i\left(\alpha\overline{\beta} + \overline{\alpha}\beta\right) & \alpha\overline{\alpha}  \beta\overline{\beta} \end{pmatrix}. }[/math]
To verify this, substitute for α. β the elements of the matrix on the RHS of (2). After some manipulation, the matrix assumes the form of the RHS of (1).
It is clear from the explicit form in terms of Euler angles that the map
 [math]\displaystyle{ \begin{cases} p:\operatorname{SU}(2) \to \operatorname{SO}(3)\\ \Pi_u(\pm g_{\alpha \beta}) \mapsto g_{\alpha \beta} \end{cases} }[/math]
just described is a smooth, 2:1 and surjective group homomorphism. It is hence an explicit description of the universal covering space of SO(3) from the universal covering group SU(2).
Lie algebra
Associated with every Lie group is its Lie algebra, a linear space of the same dimension as the Lie group, closed under a bilinear alternating product called the Lie bracket. The Lie algebra of SO(3) is denoted by [math]\displaystyle{ \mathfrak{so}(3) }[/math] and consists of all skewsymmetric 3 × 3 matrices.^{[7]} This may be seen by differentiating the orthogonality condition, A^{T}A = I, A ∈ SO(3).^{[nb 2]} The Lie bracket of two elements of [math]\displaystyle{ \mathfrak{so}(3) }[/math] is, as for the Lie algebra of every matrix group, given by the matrix commutator, [A_{1}, A_{2}] = A_{1}A_{2} − A_{2}A_{1}, which is again a skewsymmetric matrix. The Lie algebra bracket captures the essence of the Lie group product in a sense made precise by the Baker–Campbell–Hausdorff formula.
The elements of [math]\displaystyle{ \mathfrak{so}(3) }[/math] are the "infinitesimal generators" of rotations, i.e., they are the elements of the tangent space of the manifold SO(3) at the identity element. If [math]\displaystyle{ R(\phi, \boldsymbol{n}) }[/math] denotes a counterclockwise rotation with angle φ about the axis specified by the unit vector [math]\displaystyle{ \boldsymbol{n}, }[/math] then
 [math]\displaystyle{ \forall \boldsymbol{u} \in \R^3: \qquad \left. \frac{\operatorname{d}}{\operatorname{d}\phi} \right_{\phi=0} R(\phi,\boldsymbol{n}) \boldsymbol{u} = \boldsymbol{n} \times \boldsymbol{u}. }[/math]
This can be used to show that the Lie algebra [math]\displaystyle{ \mathfrak{so}(3) }[/math] (with commutator) is isomorphic to the Lie algebra [math]\displaystyle{ \R^3 }[/math] (with cross product). Under this isomorphism, an Euler vector [math]\displaystyle{ \boldsymbol{\omega}\in\R^3 }[/math] corresponds to the linear map [math]\displaystyle{ \widetilde{\boldsymbol{\omega}} }[/math] defined by [math]\displaystyle{ \widetilde{\boldsymbol{\omega}}(\boldsymbol{u}) = \boldsymbol{\omega}\times\boldsymbol{u}. }[/math]
In more detail, most often a suitable basis for [math]\displaystyle{ \mathfrak{so}(3) }[/math] as a 3dimensional vector space is
 [math]\displaystyle{ \boldsymbol{L}_x = \begin{bmatrix}0&0&0\\0&0&1\\0&1&0\end{bmatrix}, \quad \boldsymbol{L}_y = \begin{bmatrix}0&0&1\\0&0&0\\1&0&0\end{bmatrix}, \quad \boldsymbol{L}_z = \begin{bmatrix}0&1&0\\1&0&0\\0&0&0\end{bmatrix}. }[/math]
The commutation relations of these basis elements are,
 [math]\displaystyle{ [\boldsymbol{L}_x, \boldsymbol{L}_y] = \boldsymbol{L}_z, \quad [\boldsymbol{L}_z, \boldsymbol{L}_x] = \boldsymbol{L}_y, \quad [\boldsymbol{L}_y, \boldsymbol{L}_z] = \boldsymbol{L}_x }[/math]
which agree with the relations of the three standard unit vectors of [math]\displaystyle{ \R^3 }[/math] under the cross product.
As announced above, one can identify any matrix in this Lie algebra with an Euler vector [math]\displaystyle{ \boldsymbol{\omega} = (x,y,z) \in \R^3, }[/math]^{[8]}
 [math]\displaystyle{ \widehat{\boldsymbol{\omega}} =\boldsymbol{\omega}\cdot \boldsymbol{L} = x \boldsymbol{L}_x + y \boldsymbol{L}_y + z \boldsymbol{L}_z =\begin{bmatrix}0&z&y\\z&0&x\\y&x&0\end{bmatrix} \in \mathfrak{so}(3). }[/math]
This identification is sometimes called the hatmap.^{[9]} Under this identification, the [math]\displaystyle{ \mathfrak{so}(3) }[/math] bracket corresponds in [math]\displaystyle{ \R^3 }[/math] to the cross product,
 [math]\displaystyle{ \left [\widehat{\boldsymbol{u}},\widehat{\boldsymbol{v}} \right ] = \widehat{\boldsymbol{u} \times \boldsymbol{v}}. }[/math]
The matrix identified with a vector [math]\displaystyle{ \boldsymbol{u} }[/math] has the property that
 [math]\displaystyle{ \widehat{\boldsymbol{u}}\boldsymbol{v} = \boldsymbol{u} \times \boldsymbol{v}, }[/math]
where the lefthand side we have ordinary matrix multiplication. This implies [math]\displaystyle{ \boldsymbol{u} }[/math] is in the null space of the skewsymmetric matrix with which it is identified, because [math]\displaystyle{ \boldsymbol{u} \times \boldsymbol{u} = \boldsymbol{0}. }[/math]
A note on Lie algebras
In Lie algebra representations, the group SO(3) is compact and simple of rank 1, and so it has a single independent Casimir element, a quadratic invariant function of the three generators which commutes with all of them. The Killing form for the rotation group is just the Kronecker delta, and so this Casimir invariant is simply the sum of the squares of the generators, [math]\displaystyle{ \boldsymbol{J}_x, \boldsymbol{J}_y, \boldsymbol{J}_z, }[/math] of the algebra
 [math]\displaystyle{ [\boldsymbol{J}_x, \boldsymbol{J}_y] = \boldsymbol{J}_z, \quad [\boldsymbol{J}_z, \boldsymbol{J}_x] = \boldsymbol{J}_y, \quad [\boldsymbol{J}_y, \boldsymbol{J}_z] = \boldsymbol{J}_x. }[/math]
That is, the Casimir invariant is given by
 [math]\displaystyle{ \boldsymbol{J}^2\equiv \boldsymbol{J}\cdot \boldsymbol{J} =\boldsymbol{J}_x^2+\boldsymbol{J}_y^2+\boldsymbol{J}_z^2 \propto \boldsymbol{I}. }[/math]
For unitary irreducible representations D^{j}, the eigenvalues of this invariant are real and discrete, and characterize each representation, which is finite dimensional, of dimensionality [math]\displaystyle{ 2j+1 }[/math]. That is, the eigenvalues of this Casimir operator are
 [math]\displaystyle{ \boldsymbol{J}^2= j(j+1) \boldsymbol{I}_{2j+1}, }[/math]
where j is integer or halfinteger, and referred to as the spin or angular momentum.
So, the 3 × 3 generators L displayed above act on the triplet (spin 1) representation, while the 2 × 2 generators below, t, act on the doublet (spin1/2) representation. By taking Kronecker products of D^{1/2} with itself repeatedly, one may construct all higher irreducible representations D^{j}. That is, the resulting generators for higher spin systems in three spatial dimensions, for arbitrarily large j, can be calculated using these spin operators and ladder operators.
For every unitary irreducible representations D^{j} there is an equivalent one, D^{−j−1}. All infinitedimensional irreducible representations must be nonunitary, since the group is compact.
In quantum mechanics, the Casimir invariant is the "angularmomentumsquared" operator; integer values of spin j characterize bosonic representations, while halfinteger values fermionic representations. The antihermitian matrices used above are utilized as spin operators, after they are multiplied by i, so they are now hermitian (like the Pauli matrices). Thus, in this language,
 [math]\displaystyle{ [\boldsymbol{J}_x, \boldsymbol{J}_y] = i\boldsymbol{J}_z, \quad [\boldsymbol{J}_z, \boldsymbol{J}_x] = i\boldsymbol{J}_y, \quad [\boldsymbol{J}_y, \boldsymbol{J}_z] = i\boldsymbol{J}_x. }[/math]
and hence
 [math]\displaystyle{ \boldsymbol{J}^2= j(j+1) \boldsymbol{I}_{2j+1}. }[/math]
Explicit expressions for these D^{j} are,
 [math]\displaystyle{ \begin{align} \left (\boldsymbol{J}_z^{(j)}\right )_{ba} &= (j+1a)\delta_{b,a}\\ \left (\boldsymbol{J}_x^{(j)}\right )_{ba} &=\frac{1}{2} \left (\delta_{b,a+1}+\delta_{b+1,a} \right ) \sqrt{(j+1)(a+b1)ab}\\ \left (\boldsymbol{J}_y^{(j)}\right )_{ba} &=\frac{1}{2i} \left (\delta_{b,a+1}\delta_{b+1,a} \right ) \sqrt{(j+1)(a+b1)ab}\\ \end{align} }[/math]
where j is arbitrary and [math]\displaystyle{ 1 \le a, b \le 2j+1 }[/math].
For example, the resulting spin matrices for spin 1 ([math]\displaystyle{ j = 1 }[/math]) are
 [math]\displaystyle{ \begin{align} \boldsymbol{J}_x &= \frac{1}{\sqrt{2}} \begin{pmatrix} 0 &1 &0\\ 1 &0 &1\\ 0 &1 &0 \end{pmatrix} \\ \boldsymbol{J}_y &= \frac{1}{\sqrt{2}} \begin{pmatrix} 0 &i &0\\ i &0 &i\\ 0 &i &0 \end{pmatrix} \\ \boldsymbol{J}_z &= \begin{pmatrix} 1 &0 &0\\ 0 &0 &0\\ 0 &0 &1 \end{pmatrix} \end{align} }[/math]
Note, however, how these are in an equivalent, but different basis, the spherical basis, than the above iL in the Cartesian basis.^{[nb 3]}
For higher spins, such as spin 3/2 ([math]\displaystyle{ j=\tfrac{3}{2} }[/math]):
 [math]\displaystyle{ \begin{align} \boldsymbol{J}_x &= \frac{1}{2} \begin{pmatrix} 0 &\sqrt{3} &0 &0\\ \sqrt{3} &0 &2 &0\\ 0 &2 &0 &\sqrt{3}\\ 0 &0 &\sqrt{3} &0 \end{pmatrix} \\ \boldsymbol{J}_y &= \frac{1}{2} \begin{pmatrix} 0 &i\sqrt{3} &0 &0\\ i\sqrt{3} &0 &2i &0\\ 0 &2i &0 &i\sqrt{3}\\ 0 &0 &i\sqrt{3} &0 \end{pmatrix} \\ \boldsymbol{J}_z &=\frac{1}{2} \begin{pmatrix} 3 &0 &0 &0\\ 0 &1 &0 &0\\ 0 &0 &1 &0\\ 0 &0 &0 &3 \end{pmatrix}. \end{align} }[/math]
For spin 5/2 ([math]\displaystyle{ j = \tfrac{5}{2} }[/math]),
 [math]\displaystyle{ \begin{align} \boldsymbol{J}_x &= \frac{1}{2} \begin{pmatrix} 0 &\sqrt{5} &0 &0 &0 &0 \\ \sqrt{5} &0 &2\sqrt{2} &0 &0 &0 \\ 0 &2\sqrt{2} &0 &3 &0 &0 \\ 0 &0 &3 &0 &2\sqrt{2} &0 \\ 0 &0 &0 &2\sqrt{2} &0 &\sqrt{5} \\ 0 &0 &0 &0 &\sqrt{5} &0 \end{pmatrix} \\ \boldsymbol{J}_y &= \frac{1}{2} \begin{pmatrix} 0 &i\sqrt{5} &0 &0 &0 &0 \\ i\sqrt{5} &0 &2i\sqrt{2} &0 &0 &0 \\ 0 &2i\sqrt{2} &0 &3i &0 &0 \\ 0 &0 &3i &0 &2i\sqrt{2} &0 \\ 0 &0 &0 &2i\sqrt{2} &0 &i\sqrt{5} \\ 0 &0 &0 &0 &i\sqrt{5} &0 \end{pmatrix} \\ \boldsymbol{J}_z &= \frac{1}{2} \begin{pmatrix} 5 &0 &0 &0 &0 &0 \\ 0 &3 &0 &0 &0 &0 \\ 0 &0 &1 &0 &0 &0 \\ 0 &0 &0 &1 &0 &0 \\ 0 &0 &0 &0 &3 &0 \\ 0 &0 &0 &0 &0 &5 \end{pmatrix}. \end{align} }[/math]
Isomorphism with 𝖘𝖚(2)
The Lie algebras [math]\displaystyle{ \mathfrak{so}(3) }[/math] and [math]\displaystyle{ \mathfrak{su}(2) }[/math] are isomorphic. One basis for [math]\displaystyle{ \mathfrak{su}(2) }[/math] is given by^{[10]}
 [math]\displaystyle{ \boldsymbol{t}_1 = \frac{1}{2}\begin{bmatrix}0 & i\\ i & 0\end{bmatrix}, \quad \boldsymbol{t}_2 = \frac{1}{2} \begin{bmatrix}0 & 1\\ 1 & 0\end{bmatrix}, \quad \boldsymbol{t}_3 = \frac{1}{2}\begin{bmatrix}i & 0\\ 0 & i\end{bmatrix}. }[/math]
These are related to the Pauli matrices by
 [math]\displaystyle{ \boldsymbol{t}_i \longleftrightarrow \frac{1}{2i} \sigma_i. }[/math]
The Pauli matrices abide by the physicists' convention for Lie algebras. In that convention, Lie algebra elements are multiplied by i, the exponential map (below) is defined with an extra factor of i in the exponent and the structure constants remain the same, but the definition of them acquires a factor of i. Likewise, commutation relations acquire a factor of i. The commutation relations for the [math]\displaystyle{ \boldsymbol{t}_i }[/math] are
 [math]\displaystyle{ [\boldsymbol{t}_i, \boldsymbol{t}_j] = \varepsilon_{ijk}\boldsymbol{t}_k, }[/math]
where ε_{ijk} is the totally antisymmetric symbol with ε_{123} = 1. The isomorphism between [math]\displaystyle{ \mathfrak{so}(3) }[/math] and [math]\displaystyle{ \mathfrak{su}(2) }[/math] can be set up in several ways. For later convenience, [math]\displaystyle{ \mathfrak{so}(3) }[/math] and [math]\displaystyle{ \mathfrak{su}(2) }[/math] are identified by mapping
 [math]\displaystyle{ \boldsymbol{L}_x \longleftrightarrow \boldsymbol{t}_1, \quad \boldsymbol{L}_y \longleftrightarrow \boldsymbol{t}_2, \quad \boldsymbol{L}_z \longleftrightarrow \boldsymbol{t}_3, }[/math]
and extending by linearity.
Exponential map
The exponential map for SO(3), is, since SO(3) is a matrix Lie group, defined using the standard matrix exponential series,
 [math]\displaystyle{ \begin{cases} \exp : \mathfrak{so}(3) \to \operatorname{SO}(3) \\ A \mapsto e^A = \sum_{k=0}^\infty \frac{1}{k!} A^k = I + A + \tfrac{1}{2} A^2 + \cdots. \end{cases} }[/math]
For any skewsymmetric matrix A ∈ 𝖘𝖔(3), e^{A} is always in SO(3). The proof uses the elementary properties of the matrix exponential
 [math]\displaystyle{ \left(e^A\right)^\textsf{T} e^A = e^{A^\textsf{T}} e^A = e^{A^\textsf{T} + A} = e^{A + A} = e^{A  A} = e^A \left(e^A\right)^\textsf{T} = e^0 = I. }[/math]
since the matrices A and A^{T} commute, this can be easily proven with the skewsymmetric matrix condition. This is not enough to show that 𝖘𝖔(3) is the corresponding Lie algebra for SO(3), and shall be proven separately.
The level of difficulty of proof depends on how a matrix group Lie algebra is defined. (Hall 2003) defines the Lie algebra as the set of matrices
 [math]\displaystyle{ \left\{A \in \operatorname{M}(n, \R) \left e^{tA} \in \operatorname{SO}(3) \forall t\right.\right\}, }[/math]
in which case it is trivial. (Rossmann 2002) uses for a definition derivatives of smooth curve segments in SO(3) through the identity taken at the identity, in which case it is harder.^{[11]}
For a fixed A ≠ 0, e^{tA}, −∞ < t < ∞ is a oneparameter subgroup along a geodesic in SO(3). That this gives a oneparameter subgroup follows directly from properties of the exponential map.^{[12]}
The exponential map provides a diffeomorphism between a neighborhood of the origin in the 𝖘𝖔(3) and a neighborhood of the identity in the SO(3).^{[13]} For a proof, see Closed subgroup theorem.
The exponential map is surjective. This follows from the fact that every R ∈ SO(3), since every rotation leaves an axis fixed (Euler's rotation theorem), and is conjugate to a block diagonal matrix of the form
 [math]\displaystyle{ D = \begin{pmatrix}\cos \theta & \sin \theta & 0\\ \sin \theta & \cos \theta & 0\\ 0 & 0 & 1\end{pmatrix} = e^{\theta L_z}, }[/math]
such that A = BDB^{−1}, and that
 [math]\displaystyle{ Be^{\theta L_z}B^{1} = e^{B\theta L_zB^{1}}, }[/math]
together with the fact that 𝖘𝖔(3) is closed under the adjoint action of SO(3), meaning that BθL_{z}B^{−1} ∈ 𝖘𝖔(3).
Thus, e.g., it is easy to check the popular identity
 [math]\displaystyle{ e^{\pi L_x/2} e^{\theta L_z} e^{\pi L_x/2} = e^{\theta L_y}. }[/math]
As shown above, every element A ∈ 𝖘𝖔(3) is associated with a vector ω = θ u, where u = (x,y,z) is a unit magnitude vector. Since u is in the null space of A, if one now rotates to a new basis, through some other orthogonal matrix O, with u as the z axis, the final column and row of the rotation matrix in the new basis will be zero.
Thus, we know in advance from the formula for the exponential that exp(OAO^{T}) must leave u fixed. It is mathematically impossible to supply a straightforward formula for such a basis as a function of u, because its existence would violate the hairy ball theorem; but direct exponentiation is possible, and yields
 [math]\displaystyle{ \begin{align} \exp(\tilde{\boldsymbol{\omega}}) &= \exp(\theta(\boldsymbol{u\cdot L})) = \exp\left(\theta \begin{bmatrix} 0 & z & y \\ z & 0 & x \\ y & x & 0 \end{bmatrix}\right)\\[4pt] &= \boldsymbol{I} + 2cs(\boldsymbol{u\cdot L}) + 2s^2 (\boldsymbol{u \cdot L})^2 \\[4pt] &= \begin{bmatrix} 2 \left(x^2  1\right) s^2 + 1 & 2 x y s^2  2 z c s & 2 x z s^2 + 2 y c s \\ 2 x y s^2 + 2 z c s & 2 \left(y^2  1\right) s^2 + 1 & 2 y z s^2  2 x c s \\ 2 x z s^2  2 y c s & 2 y z s^2 + 2 x c s & 2 \left(z^2  1\right) s^2 + 1 \end{bmatrix}, \end{align} }[/math]
where [math]\displaystyle{ c = \cos\frac{\theta}{2} }[/math] and [math]\displaystyle{ s = \sin\frac{\theta}{2} }[/math]. This is recognized as a matrix for a rotation around axis u by the angle θ: cf. Rodrigues' rotation formula.
Logarithm map
Given R ∈ SO(3), let [math]\displaystyle{ A = \tfrac{1}{2} \left(R  R^\mathrm{T}\right) }[/math] denote the antisymmetric part and let [math]\displaystyle{ \A\ = \sqrt{\frac{1}{2}\operatorname{Tr}\left(A^2\right)}. }[/math] Then, the logarithm of R is given by^{[9]}
 [math]\displaystyle{ \log R = \frac{\sin^{1}\A\}{\A\}A. }[/math]
This is manifest by inspection of the mixed symmetry form of Rodrigues' formula,
 [math]\displaystyle{ e^X = I + \frac{\sin \theta}{\theta}X + 2\frac{\sin^2\frac{\theta}{2}}{\theta^2}X^2, \quad \theta = \X\, }[/math]
where the first and last term on the righthand side are symmetric.
Uniform random sampling
[math]\displaystyle{ SO(3) }[/math] is doubly covered by the group of unit quaternions, which is isomorphic to the 3sphere. Since the Haar measure on the unit quaternions is just the 3area measure in 4 dimensions, the Haar measure on [math]\displaystyle{ SO(3) }[/math] is just the pushforward of the 3area measure.
Consequently, generating a uniformly random rotation in [math]\displaystyle{ \R^3 }[/math] is equivalent to generating a uniformly random point on the 3sphere. This can be accomplished by the following[math]\displaystyle{ (\sqrt{1u_1}\sin(2\pi u_2), \sqrt{1u_1}\cos(2\pi u_2), \sqrt{u_1}\sin(2\pi u_3), \sqrt{u_1}\cos(2\pi u_3)) }[/math]
where [math]\displaystyle{ u_1, u_2, u_3 }[/math] are uniformly random samples of [math]\displaystyle{ [0, 1] }[/math].^{[14]}
Baker–Campbell–Hausdorff formula
Suppose X and Y in the Lie algebra are given. Their exponentials, exp(X) and exp(Y), are rotation matrices, which can be multiplied. Since the exponential map is a surjection, for some Z in the Lie algebra, exp(Z) = exp(X) exp(Y), and one may tentatively write
 [math]\displaystyle{ Z = C(X, Y), }[/math]
for C some expression in X and Y. When exp(X) and exp(Y) commute, then Z = X + Y, mimicking the behavior of complex exponentiation.
The general case is given by the more elaborate BCH formula, a series expansion of nested Lie brackets.^{[15]} For matrices, the Lie bracket is the same operation as the commutator, which monitors lack of commutativity in multiplication. This general expansion unfolds as follows,^{[nb 4]}
 [math]\displaystyle{ Z = C(X, Y) = X + Y + \frac{1}{2} [X, Y] + \tfrac{1}{12} [X, [X, Y]]  \frac{1}{12} [Y, [X, Y]] + \cdots. }[/math]
The infinite expansion in the BCH formula for SO(3) reduces to a compact form,
 [math]\displaystyle{ Z = \alpha X + \beta Y + \gamma[X, Y], }[/math]
for suitable trigonometric function coefficients (α, β, γ).
The (α, β, γ) are given by
 [math]\displaystyle{ \alpha = \phi \cot\left(\frac{\phi}{2}\right) \gamma, \qquad \beta = \theta \cot\left(\frac{\theta}{2}\right)\gamma, \qquad \gamma = \frac{\sin^{1}d}{d}\frac{c}{\theta \phi}, }[/math]
where
 [math]\displaystyle{ \begin{align} c &= \frac{1}{2}\sin\theta\sin\phi  2\sin^2\frac{\theta}{2}\sin^2\frac{\phi}{2}\cos(\angle(u, v)),\quad a = c \cot\left(\frac{\phi}{2}\right), \quad b = c \cot\left(\frac{\theta}{2}\right), \\ d &= \sqrt{a^2 + b^2 + 2ab\cos(\angle(u, v)) + c^2 \sin^2(\angle(u, v))}, \end{align} }[/math]
for
 [math]\displaystyle{ \theta = \frac{1}{\sqrt{2}}\X\,\quad \phi = \frac{1}{\sqrt{2}}\Y\,\quad \angle(u, v) = \cos^{1}\frac{\langle X, Y\rangle}{\X\\Y\}. }[/math]
The inner product is the Hilbert–Schmidt inner product and the norm is the associated norm. Under the hatisomorphism,
 [math]\displaystyle{ \langle u, v\rangle = \frac{1}{2}\operatorname{Tr}X^\mathrm{T}Y, }[/math]
which explains the factors for θ and φ. This drops out in the expression for the angle.
It is worthwhile to write this composite rotation generator as
 [math]\displaystyle{ \alpha X + \beta Y + \gamma[X, Y]\underset{\mathfrak{so}(3)}{=} X + Y + \frac{1}{2} [X, Y] + \frac{1}{12} [X, [X, Y]]  \frac{1}{12} [Y, [X, Y]] + \cdots, }[/math]
to emphasize that this is a Lie algebra identity.
The above identity holds for all faithful representations of 𝖘𝖔(3). The kernel of a Lie algebra homomorphism is an ideal, but 𝖘𝖔(3), being simple, has no nontrivial ideals and all nontrivial representations are hence faithful. It holds in particular in the doublet or spinor representation. The same explicit formula thus follows in a simpler way through Pauli matrices, cf. the 2×2 derivation for SU(2).
The Pauli vector version of the same BCH formula is the somewhat simpler group composition law of SU(2),
 [math]\displaystyle{ e^{i a'\left(\hat{u} \cdot \vec{\sigma}\right)}e^{i b'\left(\hat{v} \cdot \vec{\sigma}\right)} = \exp\left( \frac{c'}{\sin c'} \sin a' \sin b' \left(\left(i\cot b'\hat{u} + i \cot a' \hat{v}\right)\cdot\vec{\sigma} + \frac{1}{2} \left[i \hat{u} \cdot \vec{\sigma}, i \hat{v} \cdot \vec{\sigma}\right]\right) \right), }[/math]
where
 [math]\displaystyle{ \cos c' = \cos a' \cos b'  \hat{u} \cdot\hat{v} \sin a' \sin b', }[/math]
the spherical law of cosines. (Note a', b', c' are angles, not the a, b, c above.)
This is manifestly of the same format as above,
 [math]\displaystyle{ Z = \alpha' X + \beta' Y + \gamma' [X, Y], }[/math]
with
 [math]\displaystyle{ X = i a'\hat{u} \cdot \mathbf{\sigma}, \quad Y = ib'\hat{v} \cdot \mathbf{\sigma} \in \mathfrak{su}(2), }[/math]
so that
 [math]\displaystyle{ \begin{align} \alpha' &= \frac{c'}{\sin c'}\frac{\sin a'}{a'}\cos b' \\ \beta' &= \frac{c'}{\sin c'}\frac{\sin b'}{b'}\cos a' \\ \gamma' &= \frac{1}{2}\frac{c'}{\sin c'}\frac{\sin a'}{a'}\frac{\sin b'}{b'}. \end{align} }[/math]
For uniform normalization of the generators in the Lie algebra involved, express the Pauli matrices in terms of tmatrices, σ → 2i t, so that
 [math]\displaystyle{ a' \mapsto \frac{\theta}{2}, \quad b' \mapsto  \frac{\phi}{2}. }[/math]
To verify then these are the same coefficients as above, compute the ratios of the coefficients,
 [math]\displaystyle{ \begin{align} \frac{\alpha'}{\gamma'} &= \theta\cot\frac{\theta}{2} &= \frac{\alpha}{\gamma}\\ \frac{\beta'}{\gamma'} &= \phi\cot\frac{\phi}{2} &= \frac{\beta}{\gamma}. \end{align} }[/math]
Finally, γ = γ' given the identity d = sin 2c'.
For the general n × n case, one might use Ref.^{[16]}
The quaternion formulation of the composition of two rotations R_{B} and R_{A} also yields directly the rotation axis and angle of the composite rotation R_{C} = R_{B}R_{A}.
Let the quaternion associated with a spatial rotation R is constructed from its rotation axis S and the rotation angle φ this axis. The associated quaternion is given by,
 [math]\displaystyle{ S = \cos\frac{\phi}{2} + \sin\frac{\phi}{2} \mathbf{S}. }[/math]
Then the composition of the rotation R_{R} with R_{A} is the rotation R_{C} = R_{B}R_{A} with rotation axis and angle defined by the product of the quaternions
 [math]\displaystyle{ A = \cos\frac{\alpha}{2} + \sin\frac{\alpha}{2}\mathbf{A}\quad\text{ and }\quad B = \cos\frac{\beta}{2} + \sin\frac{\beta}{2}\mathbf{B}, }[/math]
that is
 [math]\displaystyle{ C = \cos\frac{\gamma}{2} + \sin\frac{\gamma}{2}\mathbf{C} = \left(\cos\frac{\beta}{2} + \sin\frac{\beta}{2}\mathbf{B}\right)\left(\cos\frac{\alpha}{2} + \sin\frac{\alpha}{2}\mathbf{A}\right). }[/math]
Expand this product to obtain
 [math]\displaystyle{ \cos\frac{\gamma}{2} + \sin\frac{\gamma}{2} \mathbf{C} = \left( \cos\frac{\beta}{2}\cos\frac{\alpha}{2}  \sin\frac{\beta}{2}\sin\frac{\alpha}{2} \mathbf{B}\cdot \mathbf{A} \right) + \left( \sin\frac{\beta}{2}\cos\frac{\alpha}{2} \mathbf{B} + \sin\frac{\alpha}{2}\cos\frac{\beta}{2} \mathbf{A} + \sin\frac{\beta}{2}\sin\frac{\alpha}{2} \mathbf{B} \times \mathbf{A} \right). }[/math]
Divide both sides of this equation by the identity, which is the law of cosines on a sphere,
 [math]\displaystyle{ \cos\frac{\gamma}{2} = \cos\frac{\beta}{2}\cos\frac{\alpha}{2}  \sin\frac{\beta}{2}\sin\frac{\alpha}{2} \mathbf{B}\cdot \mathbf{A}, }[/math]
and compute
 [math]\displaystyle{ \tan\frac{\gamma}{2} \mathbf{C} = \frac{\tan\frac{\beta}{2}\mathbf{B} + \tan\frac{\alpha}{2} \mathbf{A} + \tan\frac{\beta}{2}\tan\frac{\alpha}{2} \mathbf{B} \times \mathbf{A}}{1  \tan\frac{\beta}{2}\tan\frac{\alpha}{2} \mathbf{B} \cdot \mathbf{A}}. }[/math]
This is Rodrigues' formula for the axis of a composite rotation defined in terms of the axes of the two rotations. He derived this formula in 1840 (see page 408).^{[17]}
The three rotation axes A, B, and C form a spherical triangle and the dihedral angles between the planes formed by the sides of this triangle are defined by the rotation angles.
Infinitesimal rotations
An infinitesimal rotation matrix or differential rotation matrix is a matrix representing an infinitely small rotation. While a rotation matrix is an orthogonal matrix [math]\displaystyle{ R^\mathsf{T} = R^{1} }[/math] representing an element of [math]\displaystyle{ SO(n) }[/math] (the special orthogonal group), the differential of a rotation is a skewsymmetric matrix [math]\displaystyle{ A^\mathsf{T} = A }[/math] in the tangent space [math]\displaystyle{ \mathfrak{so}(n) }[/math] (the special orthogonal Lie algebra), which is not itself a rotation matrix.
An infinitesimal rotation matrix has the form
 [math]\displaystyle{ I + d\theta \, A, }[/math]
where [math]\displaystyle{ I }[/math] is the identity matrix, [math]\displaystyle{ d\theta }[/math] is vanishingly small, and [math]\displaystyle{ A \in \mathfrak{so}(n). }[/math]
For example, if [math]\displaystyle{ A = L_x, }[/math] representing an infinitesimal threedimensional rotation about the xaxis, a basis element of [math]\displaystyle{ \mathfrak{so}(3), }[/math]
 [math]\displaystyle{ dL_{x} = \begin{bmatrix} 1 & 0 & 0 \\ 0 & 1 & d\theta \\ 0 & d\theta & 1 \end{bmatrix}. }[/math]
Realizations of rotations
We have seen that there are a variety of ways to represent rotations:
 as orthogonal matrices with determinant 1,
 by axis and rotation angle
 in quaternion algebra with versors and the map 3sphere S^{3} → SO(3) (see quaternions and spatial rotations)
 in geometric algebra as a rotor
 as a sequence of three rotations about three fixed axes; see Euler angles.
Spherical harmonics
The group SO(3) of threedimensional Euclidean rotations has an infinitedimensional representation on the Hilbert space
 [math]\displaystyle{ L^2\left(\mathbf{S}^2\right) = \operatorname{span} \left\{ Y^\ell_m, \ell \in \N^+, \ell \leq m \leq \ell \right\}, }[/math]
where [math]\displaystyle{ Y^\ell_m }[/math] are spherical harmonics. Its elements are square integrable complexvalued functions^{[nb 5]} on the sphere. The inner product on this space is given by

[math]\displaystyle{ \langle f,g\rangle = \int_{\mathbf{S}^2}\overline{f}g\,d\Omega = \int_0^{2\pi} \int_0^\pi \overline{f}g \sin\theta \, d\theta \, d\phi. }[/math]
(
)
If f is an arbitrary square integrable function defined on the unit sphere S^{2}, then it can be expressed as^{[19]}

[math]\displaystyle{ f\rangle = \sum_{\ell = 1}^\infty\sum_{m = \ell}^{m = \ell} \leftY_m^\ell\right\rangle\left\langle Y_m^\ellf\right\rangle, \qquad f(\theta, \phi) = \sum_{\ell = 1}^\infty\sum_{m = \ell}^{m = \ell} f_{\ell m} Y^\ell_m(\theta, \phi), }[/math]
(
)
where the expansion coefficients are given by

[math]\displaystyle{ f_{\ell m} = \left\langle Y_m^\ell, f \right\rangle = \int_{\mathbf{S}^2}\overline{{Y^\ell_m}}f \, d\Omega = \int_0^{2\pi} \int_0^\pi \overline{{Y_m^\ell}}(\theta, \phi)f(\theta, \phi)\sin \theta \, d\theta \, d\phi. }[/math]
(
)
The Lorentz group action restricts to that of SO(3) and is expressed as

[math]\displaystyle{ (\Pi(R)f)(\theta(x), \phi(x)) = \sum_{\ell = 1}^\infty\sum_{m = \ell}^{m = \ell}\sum_{m' = \ell}^{m' = \ell} D^{(\ell)}_{mm'} (R) f_{\ell m'}Y^\ell_m \left(\theta\left(R^{1}x\right), \phi\left(R^{1}x\right)\right), \qquad R \in \operatorname{SO}(3), \quad x \in \mathbf{S}^2. }[/math]
(
)
This action is unitary, meaning that

[math]\displaystyle{ \langle \Pi(R)f,\Pi(R)g\rangle = \langle f,g \rangle \qquad \forall f,g \in \mathbf{S}^2, \quad \forall R \in \operatorname{SO}(3). }[/math]
(
)
The D^{(ℓ)} can be obtained from the D^{(m, n)} of above using Clebsch–Gordan decomposition, but they are more easily directly expressed as an exponential of an odddimensional su(2)representation (the 3dimensional one is exactly 𝖘𝖔(3)).^{[20]}^{[21]} In this case the space L^{2}(S^{2}) decomposes neatly into an infinite direct sum of irreducible odd finitedimensional representations V_{2i + 1}, i = 0, 1, ... according to^{[22]}

[math]\displaystyle{ L^2\left(\mathbf{S}^2\right) = \sum_{i = 0}^\infty V_{2i + 1} \equiv \bigoplus_{i=0}^\infty \operatorname{span}\left\{Y_m^{2i+1}\right\}. }[/math]
(
)
This is characteristic of infinitedimensional unitary representations of SO(3). If Π is an infinitedimensional unitary representation on a separable^{[nb 6]} Hilbert space, then it decomposes as a direct sum of finitedimensional unitary representations.^{[19]} Such a representation is thus never irreducible. All irreducible finitedimensional representations (Π, V) can be made unitary by an appropriate choice of inner product,^{[19]}
 [math]\displaystyle{ \langle f, g\rangle_U \equiv \int_{\operatorname{SO}(3)} \langle\Pi(R)f, \Pi(R)g\rangle \, dg = \frac{1}{8\pi^2} \int_0^{2\pi} \int_0^\pi \int_0^{2\pi} \langle \Pi(R)f, \Pi(R)g\rangle \sin \theta \, d\phi \, d\theta \, d\psi, \quad f,g \in V, }[/math]
where the integral is the unique invariant integral over SO(3) normalized to 1, here expressed using the Euler angles parametrization. The inner product inside the integral is any inner product on V.
Generalizations
The rotation group generalizes quite naturally to ndimensional Euclidean space, [math]\displaystyle{ \R^n }[/math] with its standard Euclidean structure. The group of all proper and improper rotations in n dimensions is called the orthogonal group O(n), and the subgroup of proper rotations is called the special orthogonal group SO(n), which is a Lie group of dimension n(n − 1)/2.
In special relativity, one works in a 4dimensional vector space, known as Minkowski space rather than 3dimensional Euclidean space. Unlike Euclidean space, Minkowski space has an inner product with an indefinite signature. However, one can still define generalized rotations which preserve this inner product. Such generalized rotations are known as Lorentz transformations and the group of all such transformations is called the Lorentz group.
The rotation group SO(3) can be described as a subgroup of E^{+}(3), the Euclidean group of direct isometries of Euclidean [math]\displaystyle{ \R^3. }[/math] This larger group is the group of all motions of a rigid body: each of these is a combination of a rotation about an arbitrary axis and a translation, or put differently, a combination of an element of SO(3) and an arbitrary translation.
In general, the rotation group of an object is the symmetry group within the group of direct isometries; in other words, the intersection of the full symmetry group and the group of direct isometries. For chiral objects it is the same as the full symmetry group.
See also
 Orthogonal group
 Angular momentum
 Coordinate rotations
 Charts on SO(3)
 Representations of SO(3)
 Euler angles
 Rodrigues' rotation formula
 Infinitesimal rotation
 Pin group
 Quaternions and spatial rotations
 Rigid body
 Spherical harmonics
 Plane of rotation
 Lie group
 Pauli matrix
 Plate trick
 Threedimensional rotation operator
Footnotes
 ↑ This is effected by first applying a rotation [math]\displaystyle{ g_{\theta} }[/math] through φ about the zaxis to take the xaxis to the line L, the intersection between the planes xy and x'y', the latter being the rotated xyplane. Then rotate with [math]\displaystyle{ g_{\theta} }[/math] through θ about L to obtain the new zaxis from the old one, and finally rotate by [math]\displaystyle{ g_{\psi} }[/math] through an angle ψ about the new zaxis, where ψ is the angle between L and the new xaxis. In the equation, [math]\displaystyle{ g_{\theta} }[/math] and [math]\displaystyle{ g_{\psi} }[/math] are expressed in a temporary rotated basis at each step, which is seen from their simple form. To transform these back to the original basis, observe that [math]\displaystyle{ \mathbf{g}_{\theta} = g_{\phi}g_{\theta}g_{\phi}^{1}. }[/math] Here boldface means that the rotation is expressed in the original basis. Likewise,
 [math]\displaystyle{ \mathbf{g}_{\psi} = g_{\phi}g_{\theta}g_{\phi}^{1} g_{\phi} g_{\psi} \left [ g_{\phi}g_{\theta}g_{\phi}^{1} g_{\phi} \right ]^{1}. }[/math]
 [math]\displaystyle{ \mathbf{g}_{\psi}\mathbf{g}_{\theta}\mathbf{g}_{\phi} = g_{\phi}g_{\theta}g_{\phi}^{1} g_{\phi}g_{\psi} \left [g_{\phi} g_{\theta} g_{\phi}^{1} g_{\phi} \right ]^{1} * g_{\phi}g_{\theta}g_{\phi}^{1}* g_{\phi} = g_{\phi}g_{\theta}g_{\psi}. }[/math]
 ↑ For an alternative derivation of [math]\displaystyle{ \mathfrak{so}(3) }[/math], see Classical group.
 ↑ Specifically, [math]\displaystyle{ \boldsymbol{U} \boldsymbol{J}_{\alpha}\boldsymbol{U}^\dagger=i\boldsymbol{L}_\alpha }[/math] for
 [math]\displaystyle{ \boldsymbol{U}= \left( \begin{array}{ccc} \frac{i}{\sqrt{2}} & 0 & \frac{i}{\sqrt{2}} \\ \frac{1}{\sqrt{2}} & 0 & \frac{1}{\sqrt{2}} \\ 0 & i & 0 \\ \end{array} \right). }[/math]
 ↑ For a full proof, see Derivative of the exponential map. Issues of convergence of this series to the correct element of the Lie algebra are here swept under the carpet. Convergence is guaranteed when [math]\displaystyle{ \X\ + \Y\ \lt \log 2 }[/math] and [math]\displaystyle{ \Z\ \lt \log 2. }[/math] The series may still converge even if these conditions are not fulfilled. A solution always exists since exp is onto in the cases under consideration.
 ↑ The elements of L^{2}(S^{2}) are actually equivalence classes of functions. two functions are declared equivalent if they differ merely on a set of measure zero. The integral is the Lebesgue integral in order to obtain a complete inner product space.
 ↑ A Hilbert space is separable if and only if it has a countable basis. All separable Hilbert spaces are isomorphic.
References
 ↑ Jacobson (2009), p. 34, Ex. 14.
 ↑ n × n real matrices are identical to linear transformations of [math]\displaystyle{ \R^n }[/math] expressed in its standard basis.
 ↑ Coxeter, H. S. M. (1973). Regular polytopes (Third ed.). New York. p. 53. ISBN 0486614808.
 ↑ Hall 2015 Proposition 1.17
 ↑ Rossmann 2002 p. 95.
 ↑ These expressions were, in fact, seminal in the development of quantum mechanics in the 1930s, cf. Ch III, § 16, B.L. van der Waerden, 1932/1932
 ↑ Hall 2015 Proposition 3.24
 ↑ Rossmann 2002
 ↑ ^{9.0} ^{9.1} Engø 2001
 ↑ Hall 2015 Example 3.27
 ↑ See Rossmann 2002, theorem 3, section 2.2.
 ↑ Rossmann 2002 Section 1.1.
 ↑ Hall 2003 Theorem 2.27.
 ↑ Shoemake, Ken (19920101), Kirk, DAVID, ed., "III.6  Uniform Random Rotations" (in en), Graphics Gems III (IBM Version) (San Francisco: Morgan Kaufmann): pp. 124–132, ISBN 9780124096738, https://www.sciencedirect.com/science/article/pii/B9780080507552500361, retrieved 20220729
 ↑ Hall 2003, Ch. 3; Varadarajan 1984, §2.15
 ↑ Curtright, Fairlie & Zachos 2014 Group elements of SU(2) are expressed in closed form as finite polynomials of the Lie algebra generators, for all definite spin representations of the rotation group.
 ↑ Rodrigues, O. (1840), Des lois géométriques qui régissent les déplacements d'un système solide dans l'espace, et la variation des coordonnées provenant de ses déplacements con sidérés indépendamment des causes qui peuvent les produire, Journal de Mathématiques Pures et Appliquées de Liouville 5, 380–440.
 ↑ (Goldstein Poole)
 ↑ ^{19.0} ^{19.1} ^{19.2} Gelfand, Minlos & Shapiro 1963
 ↑ In Quantum Mechanics – nonrelativistic theory by Landau and Lifshitz the lowest order D are calculated analytically.
 ↑ Curtright, Fairlie & Zachos 2014 A formula for D^{(ℓ)} valid for all ℓ is given.
 ↑ Hall 2003 Section 4.3.5.
Bibliography
 Boas, Mary L. (2006), Mathematical Methods in the Physical Sciences (3rd ed.), John Wiley & sons, pp. 120, 127, 129,155ff and 535, ISBN 9780471198260
 Curtright, T. L.; Fairlie, D. B.; Zachos, C. K. (2014), "A compact formula for rotations as spin matrix polynomials", SIGMA 10: 084, doi:10.3842/SIGMA.2014.084, Bibcode: 2014SIGMA..10..084C
 Engø, Kenth (2001), "On the BCHformula in 𝖘𝖔(3)", BIT Numerical Mathematics 41 (3): 629–632, doi:10.1023/A:1021979515229, ISSN 00063835 [1]
 Gelfand, I.M.; Minlos, R.A.; Shapiro, Z.Ya. (1963), Representations of the Rotation and Lorentz Groups and their Applications, New York: Pergamon Press
 Goldstein, Herbert; Poole, Charles P.; Safko, John L. (2002), Classical Mechanics (third ed.), Addison Wesley, ISBN 9780201657029
 Hall, Brian C. (2015), Lie Groups, Lie Algebras, and Representations: An Elementary Introduction, Graduate Texts in Mathematics, 222 (2nd ed.), Springer, ISBN 9783319134666
 Hall, Brian C. (2003). Lie groups, Lie algebras, and representations : an elementary introduction. Graduate Texts in Mathematics. 222. New York: Springer. ISBN 0387401229.
 Jacobson, Nathan (2009), Basic algebra, 1 (2nd ed.), Dover Publications, ISBN 9780486471891
 Joshi, A. W. (2007), Elements of Group Theory for Physicists, New Age International, pp. 111ff, ISBN 9788122409758
 Rossmann, Wulf (2002), Lie Groups – An Introduction Through Linear Groups, Oxford Graduate Texts in Mathematics, Oxford Science Publications, ISBN 0198596839
 van der Waerden, B. L. (1952), Group Theory and Quantum Mechanics, Springer Publishing, ISBN 9783642658624 (translation of the original 1932 edition, Die Gruppentheoretische Methode in Der Quantenmechanik).
 Varadarajan, V. S. (1984). Lie groups, Lie algebras, and their representations. New York: SpringerVerlag. ISBN 9780387909691.
 Veltman, M.; 't Hooft, G.; de Wit, B. (2007). "Lie Groups in Physics (online lecture)". http://www.staff.science.uu.nl/~hooft101/lectures/lieg07.pdf..
Original source: https://en.wikipedia.org/wiki/3D rotation group.
Read more 