Using Rotations to Control Observable Relativistic Effects

: This paper examines the possibility of controlling the outcome of measured (flat space-time) relativistic effects, such as time dilation or length contractions, using pure rotations and their nontrivial interactions with Lorentz boosts in the isometry group SO + ( 3,1 ) . In particular, boost contributions may annihilate leaving only a geometric phase (Wigner rotation), which we see in the complex solutions of the generalized Euler decomposition problem in R 3 . We consider numerical examples involving specific matrix factorizations, along with possible applications in special relativity, electrodynamics and quantum scattering. For clearer interpretation and simplified calculations we use a convenient projective biquaternion parametrization which emphasizes the geometric phases and for a large class of problems allows for closed-form solutions in terms of only rational functions.


Introduction
The peculiar role of rotations in special relativity (SR) has been puzzling theoretical physicists for more than a century, starting with the famous Erenfest's paradox: a seeming contradiction between classical rigid body mechanics ('prohibiting' shape deformations) and Lorentz length contraction of a spinning disk (cf.[1]).On top of that, we have the infamous Sagnac effect of light circulating around a rotating closed path with different time periods, depending on the direction.This for a while led people to believe there must be something wrong with the postulates of relativity, which was Sagnac's original intention (see [2]).Nowadays the phenomenon is not only well integrated into SR but also plays a major role in laser interferometers and precise Doppler gyroscopes [3][4][5].Another interesting discovery from the early days of relativity was the Thomas precession [6] explaining how electron spins are affected by acceleration in non-inertial frames.This mechanism of boostrotation entanglement, known more generally as Wigner rotation, is quite similar to the Coriolis effect and may be attributed to a classical geometric phase [7,8].The same pattern is observed in many areas of theoretical physics and engineering: from classical mechanics to geometric optics and even quantum computation [9][10][11][12][13].A particularly curious context with plenty of practical implementations lies within the realm of electrodynamics [14][15][16], where the magnetic force itself may be regarded as a consequence of such inertial effect and thus linked to a geometric phase following a mechanical analogy [16][17][18].We refer to [19][20][21][22] for a deeper investigation of the problem of rotating frames in relativity, both from a purely mechanical perspective and with an account of electromagnetic interactions.Instead of diving into the physics and various applications of this topic, the present study discusses a unified geometric approach based on the structure of proper Lorentz transformations SO + (3, 1) ∼ = SO(3, C) and the corresponding Wigner little groups, starting with the compact subgroup of euclidean rotations in R 3 .A convenient representation proposed by Rodrigues back in 1840, based on Clifford's geometric algebras but preceeding their invention (cf.[23,24]) was rediscovered by Fedorov [25] almost a century and a half later in the context of SR and high energy physics.This paper exploits it for the description of interactions between rotations and Lorentz boosts.We begin by briefly revising this useful tool and then explore its applications to the matter of interest: from matrix decompositions to relativistic kinematics, electrodynamics and quantum scattering, providing also examples.

Preliminaries
Rotations in 3D Euclidean space provide a first glimpse of the weirdness of motion due to the Coriolis effect which we may attribute to a classical geometric phase, associated with parallel transport (see Figure 1).The complex projective representation of the proper Lorentz group SO + (3, 1) ∼ = SO(3, C) allows us to extend this weirdness quite smoothly to the realm of special relativity, where the interactions of pure rotations (the compact subgroup) and boosts (the quotient rapidity space B 3 ) create a whole spectrum of effects due to Wigner, Sagnac, Hall, etc. [26,27].The projective biquaternion (complex Rodrigues' vector) c ∈ CP 3 is a natural extension of the rational parametrization of SO(3) (see [25,26] for details) and inherits the group composition law compatible with matrix multiplication where the SO(3) or SO(3, C) representation is given, respectively, via the Cayley transform In the compact (real) setting this coincides with the well-known Rodrigues' rotation formula with the Euler substitution c = tan φ 2 n, i.e., the angle of rotation and the unit normal to the rotation plane is expressed as φ = 2 arctan |c| and n = |c| −1 c, while c = q −1 0 q is the projection of the unit quaternion q ∈ S 3 describing the lift of our rotation up to the spin covering group SU(2) ∼ = Spin(3).Hence, q does not need to be a unit in the first place: we can start with any element of H × .This simplifies our construction significantly, apart from having to deal with only rational expressions, and it remains true also in the complex setting, as the pseudoscalar I in the even part of the space-time algebra Cliff belongs to the center, and therefore, we identify it with the imaginary unit i = √ −1.Unlike in the compact real case, however, the group does not cover the whole projective space: one has to remove the quadric c 2 + 1 = 0, which corresponds to the boundary of the unit ball B 3 .This representation is quite convenient both for practical purposes (due to the reduced number of computations) and theoretical studies, e.g., for the classification of transformations as elliptic (c 2 > 0), hyperbolic (c 2 < 0), parabolic (c 2 = 0) or loxodromic ) seen here as the Möbius group of the 2-sphere S 2 ∼ = CP 2 .The realization in Minkowski space-time R 3,1 is also straightforward (see [25] for details): and denoting Λ = Λ − η Λ t η with η = diag(1, 1, 1, −1) we easily obtain the inverse map However, one needs to keep in mind that there is a non-trivial Plücker condition (namely Im c 2 = 0) to ensure the existence of an invariant plane Σ, be it space-like, time-like or null (isotropic).It is determined by the real and imaginary parts of the projective blade [28] we use this to realize the embedding into higher-dimensional groups via the twistor correspondence).Note also that Formula (4) is ill-defined for isotropic directions of infinite scale.

Generalized Euler Decompositions
The classical factorization result obtained by Euler for SO(3) is shown to generalize for any set of rotation axes as long as the one in the middle is perpendicular to the other two (Davenport's condition).Otherwise, such decomposition is also possible in many cases but does not cover the entire group.In [29] we use the projective bivector representation described above to derive a closed form solution with the necessary and sufficient condition where {a i } represent the invariant unit vectors and the angles are φ i = 2 arctan τ i .Denoting as well as r ij = a i • Ra j and g ij = a i • a j , we use the eigenvector property to express the τ i 's (unless in the so-called 'gimbal lock' setting a 3 = ±R a 1 , which has a singular solution, also studied in [29]) where ε ijk stands for the Levi-Civita symbol.Finally, τ i ∈ R ∪ ∞ if and only if (see Figure 1) while if the above is not satisfied, the solutions (7) have imaginary counterparts.Before discussing this in detail, let us point out that half-turns correspond to infinite elements in RP 3 , so (1) and (2) should be applied with the aid of l'Hôpital's rule.Consider, for example, and the set of invariant axes, given by the unit vectors Using our algorithm, we obtain for R = R 3 R 2 R 1 only one solution (since we have ∆ = 0) and the corresponding rotation matrices in the decomposition (13) take the following form: As mentioned above, there is a singular 'gimbal lock' setting, determined by the condition a 3 = ±R a 1 , to which (7) does not apply, but instead we have a degenerate solution, namely where a [i b j] = a i b j − a j b i is the alternator and ( 8) is reduced to g 12 = ±g 23 .This expression is valid also in the case of two axes R = R 2 R 1 (with φ 3 = 0) as long as the necessary and sufficient condition for that r 21 = g 21 is satisfied.For closed rotation sequence R 3 R 2 R 1 = 1 we have the trivial solution φ i = 0 for any axis configuration, a singular one with φ 1 ± φ 3 = φ 2 = 0 if g 13 = ±1 (a 3 = ±a 1 ) and for ω 2 ̸ = 0 we also have a regular solution in the form , Most of the above construction is easily transferred to the hyperbolic case of SO + (2, 1) where we have an isotropic singularity, apart from the gimbal lock, which is still present.

Complex Solutions
Whenever the condition (8) is not satisfied, one ends up with complex solutions and according to Formula (7), they come in conjugate pairs which clearly correspond to Lorentz transformations of the loxodromic type.Since there is no invariant plane for these factors, their physical interpretations are somewhat ambiguous (compared to the so-called Wigner little groups).However, we may still factorize using (1): and thus obtain the parameter values allowing for such pairs of transformations in the form Note, that on the unit circle, S 1 : τ = e iϑ one has γ = 1 (quarter-turn) and δ = tan ϑ 2 , so that (1) yields the Cayley transform.In that case τ cannot be purely imaginary, as it leads to an ill-defined Lorentz boost with τ 2 = −1.Apart from these two singular cases, the above decomposition yields a pair of well-defined mutually commuting rotations and boosts This means on the one hand, that a sequence of finite Lorentz transformations may result in a pure rotation, while on the other, it seems possible to eliminate the resultant relativistic effect even for a sequence of pure boosts if we compose them with appropriate rotations.Certainly, not all such sequences would allow it, e.g., if the Davenport condition g 12 = g 23 = 0 in the SO(3, C) representation of the boosts holds, then the procedure described above cannot result in a pure rotation.
Next, we consider an example of such decomposition, using the same system of axes {a i }, but this time the compound rotation will be a (positively oriented) quarter-turn about the x-axis.Since r 31 = g 31 = −1, both the gimbal lock condition and the one ensuring the decomposition hold, but (8) does not, so we end up with a pair of complex solutions: where we have used (12) to factorize the loxodromic transformation with invariant axis a 1 into a rotation and a pure Lorentz boost.Note, however, that the boosts generated by ±ia 3 lead to ill-defined transformations with c 2 + 1 = 0, so one may attempt to avoid this by composing c = ⟨c ± , ±3 ± 1 2 a 1 ⟩ where c± = ⟨±ia 3 , i3 ±1 a 1 ⟩, but c+ hits infinity in an isotropic direction while c2 − = −1, i.e., they are both unfit in this case, yet we may factorize: using the fact that c± in our example satisfy the Plücker condition.However, both solutions yield a 'forbidden' boost in the form of −ia 3 so in this case the decomposition is impossible even in SO(3, C) due to the isotropic singularity.That is a somewhat exotic configuration and thus worth taking into account.We shall provide more numerical examples at the end.

Cancellation of Boosts
So far we discussed only a quite specific mechanism for boost cancellation associated with the complex solutions in the generalized Euler decomposition problem.Our logic was somewhat backward: constructing the sequence of commuting boost-rotation pairs that produce the desired result as a combination of rotational counterpart and a geometric phase.More generally, one may pose the question 'What conditions should a set of Lorentz transformations satisfy so that their composition is a pure rotation?',but even for just two factors the algebraic system is quite cumbersome.However, if we impose once more the additional restriction c k = τ k n (n ∈ S 2 ), since c 2 × c 1 = 0, the condition is greatly simplified On the other hand, one may also demand from the individual factors to be pure boosts, i.e., c k = iµ k a k with a k ∈ S 2 and µ k ∈ R/{±1}.In that setting decompositions into pairs are no longer an option: for collinear axes, one has µ 1 + µ 2 = 0, which yields a trivial result, while for a 1 × a 2 ̸ = 0 the system is inconsistent.In the case of three factors we use (1) to express and consider only odd power terms as they contribute to the imaginary part.It is simplified greatly if we demand the vectors a k to be coplanar, as in that case we have c 3 ∧ c 2 ∧ c 1 = 0, which yields the system of quadratic equations for the unknown parameters µ i in the form Note, that the latter is unsolvable in the Davenport setting, but other symmetric configurations allow it, e.g., for g 12 = g 23 = 1 2 and g 13 = 1 we easily obtain a dependent solution with µ 2 = −2µ 1 and µ 3 = µ −1 2 .On the other hand, if the condition g 12 = g 23 = 0 that guarantees (8) holds, e.g., in the classical Euler setting a 1 = a 3 ⊥ a 2 , we typically end up with either the trivial solution, some contradiction, or the 'forbidden quadric' c 2 + 1 = 0.
One may come to the conclusion that the Davenport condition prevents us from having a complex solution in the decomposition problem, but it would be wrong, for example, if the axes {a i } form an orthonormal basis in R 3 (the Bryan setting), there is an interesting solution τ 1 τ 3 = −τ 2 = 1, such that all three parameters lie on the unit circle, so we can write τ 1 = τ3 = e iϑ and align {a k } with the standard basis, which yields the factorization: But how is this even possible when real solutions are guaranteed for any set of for orthogonal axes?In short, the result is a half-turn about a 1 + a 3 , so we are once more in a 'gimbal lock' mode, which allows for complex continuation as long as φ 1 + φ 3 ∈ R in Formula ( 9).Now, let us further factorize the above loxodromic transformations into pairs of mutually commuting boosts and rotations using Formulas (11) and (12).As pointed out earlier, for τ = e iϑ we have γ = 1, δ = tan ϑ 2 , which amounts in this example to a multiplication of the quarter-turn product R = R(a 3 )R 2 R(a 1 ) = R with a pair of boosts Λ 1,3 generated, respectively, by iδa 1 and −iδa 3 .To be explicit, we also provide the R 3,1 matrix representation Alternatively, one may apply the boosts first, since they commute with the corresponding quarter-turns and the transformation P → Λ 3 P Λ 1 preserves both R and R 2 , as it should.

Clifford's Perspective
Let us point out that the matrix representations of our solutions are provided here only because it is a common practice: all calculations are made in the Clifford's geometric algebra Cliff • 3,1 ∼ = Cliff 3 ∼ = H C or rather, its projective version (projective quaternions were discovered by Rodrigues years before Hamilton gave his version, but he had no clue of Rodrigues' invention which solves the 'half-angle problem' and so would help in the scientific debate) introduced in [26].For example, the matrix transformations in the above decomposition may be expressed in Cliff 3 as follows: (1 + cos ϑ e 12 − sin ϑ e 3 )(1 − e 13 )(1 + cos ϑ e 23 + sin ϑ e 1 ) = 2 cos ϑ(e 12 + e 23 ) and as the scalar part of the product is zero, projection onto the hyperplane ξ 0 = 1 in H C ∼ = C 4 yields a half-turn (geometrically, an infinite point in RP 3 ) in the plane determined by e 12 + e 23 .Here, we identify {a i } with the standard orthonormal basis {e i } and denote: A pure rotation in this framework is represented by an even, invertible (i.e., non-vanishing) element ξ of Cliff 3 expressed in homogeneous coordinates and after projection onto RP 3 as where ⟨•⟩ k denotes the grade projection, e.g., pure bivectors ⟨ξ⟩ 0 = 0 are mapped to the plane at infinity and correspond to half-turns in SO(3), c 2 = 1 yields a quarter-turn and c = 0 represents the trivial element, while inversion is just multiplication with −1 since ⟨c, −c⟩ = 0.If we allow for odd grades as well (vectors and pseudo-scalars), this results in where ξ stands for the bivector part, as before.Projections in this setting involve the center of the Lipschitz group Cliff × 3 , which includes the pseudo-scalar I and is isomorphic to C × : Here, ι denotes the embedding into CP 3 and by abuse of notation we let ξ, η represent both pure quaternions and real vectors.The sign for the pseudo-scalar part is chosen for convenience and I itself is mapped to the imaginary unit i. Note, that in the case of pure rotation (ξ 0 η = η 0 ξ) we end up with the familiar expression ξ −1 0 ξ, while in order to have a pure Lorentz boost, the condition ξ 0 ξ + η 0 η = 0 needs to be met, which leads similarly to a parametrization of the unit ball B 3 with iξ −1 0 η.The Plücker relation is expressed in the form and in the case η 0 = 0 it is reduced to ξ ⊥ η.Note, that the above construction yields more than the proper Lorentz group.For example, if (19) holds and ξ • η = −1, we hit the singularity c 2 + 1 = 0 which can be written in this case also as ξ η + I = 0. Setting ξ 0 = η 0 = 0 leads to a proper half-turn if η = 0 and a 'forbidden' boost in the case ξ 2 = η 2 .

Kinematics and Relativity
Here, we discuss the above ideas in the continuous setting, so a convenient representation of the spatial kinematics is needed: either as a linear quaternion equation or a Ricatti ODE in the projective setting derived as follows: consider the Maurer-Cartan form associated with a smooth spatial trajectory R ∋ t → R(t) ∈ SO(3), which yields the angular velocity respectively, in the stationary (lab) and co-moving (body) frame.Straightforward differentiation of the Cayley transform (2) allows for expressing the above formulas equivalently as and solving for ċ via inversion, we end up with the (vector) Riccati equation Note, that the non-commutative term c × ċ in ( 21) is associated with the geometric phase, e.g., Coriolis effect or Thomas precession [12,27].In particular, when it is real, for a real ω it suffices to have Im( ċ) = 0 even if c is itself complex, as long as the Plücker condition Im(c 2 ) = 0 is satisfied, providing a pure rotation on top of a Lorentz boost.One has ċ ∥ ω ± ∈ R 3 also from (22) in that setting.Now, let us try to find other complex paths c(t) in the Lorentz group that lead to real angular velocities (21) as we did in the discrete case using Euler decomposition.Denoting c = α + iβ with α, β ∈ RP 3 and α ⊥ β, we arrive at where the non-commutative terms in the real part account for the Coriolis and Thomas effects, while those in the imaginary part -for the Sagnac and Hall effects (see [26,27] for more details).Making the latter disappear means finding a subset of trajectories, for which i.e., the two contributions cancel each other.This effectively links the two paths α(t) and β(t) on RP 3 .In particular, for a given α(t), this results in a linear ODE for β(t) in the form which is resolved with a matrix exponent (possibly time-ordered or via Magnus expansion).
If the unit vector m processes uniformly in a circular orbit, then α has fixed direction and linear growth, like a resonance amplitude.This corresponds to a plane rotation evolving with the inverse tangent of the time parameter, and the angular velocity remains real along the complex trajectories.Thus, certain kinematic effects of β should be 'canceled out' by α.
The above solutions have quite straightforward interpretations in a relativistic context, e.g., a processing pure boost may be complemented by a moment-like rotation vector in order to eliminate the imaginary part of the generalized angular velocity vector (23).However, simply adding the two components in the complex vector parameter c = α + iβ is by no means equivalent to a composition, so we use (1) to express the corresponding factorization which leads to a self-contradicting system of equations (implying α Hence, one cannot preserve the boost part in such decompositions, at least in the Plücker setting, while on the other hand, it is possible to write instead, as already pointed out in ( 13) Now, combining the above with (25) we can easily express the decomposition into a pair of mutually orthogonal rotations and boost corresponding to a Lorentz transformation satisfying the Plücker relation, such that the angular velocity vector given by Formula ( 21) is real.But how to interpret this real angular velocity corresponding to a complex projective bivector?It helps to go back to the linear quaternion equations from which (22) are derived.They clearly show it depends on the initial condition whether one has a pure rotation or more general Lorentz kinematics even for ω ± ∈ R 3 , so this merely yields invariance of the compact subgroup SO(3).However, one can use the above solutions to manipulate the observable relativistic effects via boost-rotation interactions.In physical terms we express where v is the velocity with v = ||v|| and ψ the corresponding rapidity, c the speed of light in vacuum and n, φ the axis-angle coordinates of the rotational motion.Thus, we easily evaluate the effect on the Lorentz factor of time dilation ∆t → γ∆t with γ = cosh ψ or the Thomas precession phase adjusting the electron spin in non-inertial relativistic frames as: Note, that the Plücker condition in this context means that the translation motion takes place in the rotation plane.We also discussed the setting of screw motion (11) in which α ∥ β.Both cases are interesting and have relevant applications but only the former is associated with an invariant plane in the R 3,1 representation.In particular, for massless particles, e.g., gluons, the Lorentz transformation ( 2) is parabolic, as c 2 = 0, and reduces to Note, also that Formula (28) yields a straightforward linear relation between the boosts i β± : confirmed by the Cayley representation (2).As illustrated by the examples below, in some cases this technique provides the usual polar decomposition.However, it only changes the orientation of the invariant plane, keeping the rapidity (hence, the Lorentz factor) constant.

Geometric Phases in Electrodynamics
One may consider possible applications in electrodynamics based on the so-called complex representation of the electromagnetic field [14,30] F = B + iE where E, B are pure quaternions corresponding to the electric and magnetic component (some authors put them in reverse order as they work in a split-quaternion basis spanned by the Pauli matrices).We think of the magnetic field as a relativistic effect (Thomas precession) associated with non-inertial reference frames (see [16]) as a consequence of Lorentz covariance.Note, that the Maxwell field (treated also as a tensor, principal bundle curvature, or an element of Clifford's algebra of differential forms) is not a representative of SO(3, C), but rather, of its Lie algebra sl 2 so, in general, it cannot be directly associated to a projective biquaternion (vector-parameter) c = α + iβ.However, in the isotropic case (31) one has equivalence between Cayley and exponential maps, namely: so at least for plane waves our representation remains valid and we can write where Θ is proportional to the Maxwell tensor.Even though (1) may not be relevant in this context, the above kinematic considerations are based solely on the Cayley map and thus still apply.But here ( 23), ( 25) and ( 26) have different interpretations-they describe the aforementioned relativistic concept of magnetism, hence various interactions with the particles' spins.Manipulating the imaginary part of ω in this context is equivalent to controlling the Lorentz force and the measured Hall effect.There are other interesting applications of geometric phases in electrodynamics-some related to electronics and engineering, others to the theory of light propagation (we refer to [11] for an extensive review on the matter).
One particular example is the interaction of electric fields with particle spins which allows for more convenient hardware in quantum computers (see [31]).In [32] one may find the simple mathematical treatment of the Thomas factor in the spin-orbit interaction, which in our approach emerges naturally from the complex extension of the algebra (see also [26]).

The 2 + 1 Setting
From a practical perspective, it would be interesting to see also how the cancellation and phase effects take place in the non-compact real form (Wigner little group [6]) of SO(3, C): related to classical hyperbolic geometry and quantum mechanical scattering.It is a toy model for special relativity with real implementations, for instance in graphenes (see [33]).
In [27] we consider the Wigner rotation in this setting and its various physical applications: from quantum scattering to the falling cat problem [9,10].The latter may be regarded as yet another application of classical geometric phases.The cat's rotation during free fall seemingly violates the conservation principle for angular momenta while it is a result of a Coriolis-type effect emerging in a sequence of shear transformations.In the plane projection (e.g., Arnold's cat) we may think of composing non-parallel SL(2, R) boosts generating a U(1) phase (see [9,27] for details).The monodromy matrix approach in quantum scattering leads to a quite similar mathematical description.Namely, one defines to be the monodromy matrix corresponding to a given scattering potential.Its entries are given by the complex transition and reflection coefficients (denoted, respectively, as t and r) relating the left and right free particle asymptotic solutions in the linear scattering problem: Thus, a superposition of two scattering potentials V 1 and V 2 (in this order) yields a compound monodromy matrix M = M 2 M 1 .In particular, if we choose M 1,2 to be purely hyperbolic, i.e., t 1,2 ∈ R, the Wigner rotation angle is given by the elliptic factor M • in the composition M = MM • , e.g., Note, that as a real form of SO + (3, 1) the three-dimensional Lorentz group inherits the composition law (1) and the Cayley map construction (2) along with the generalized Euler decomposition (both real solutions for time-like axes and imaginary for space-like ones yield rotations) and the kinematic analogy (22).However, the non-euclidean metric η = diag(1, 1, −1) takes part in calculating projectors, dot and cross products, which leads to another type of singularity, apart from the gimbal lock, referred to as 'isotropic' in [29].
One particular advantage of the 3D setting is that we can use a dual extension and the transfer principle to describe the non-homogeneous Poincaré group, like in the Euclidean case [34].

Back to Factorizations
Now, let us come back to matrix factorizations with broadened horizons and more examples.If we consider the generalized Euler decomposition problem for the proper Lorentz group in the SO(3, C) representation, the results given in ( 7) and ( 9) still apply.However, even if we choose normalized vectors a k ∈ C 3 satisfying the Plücker condition a 2 k ∈ R, the complex parameters τ k in the solution may violate it, leaving us with no invariant plane.Thus, one may start with a 2 k ∈ S 2 and then use (11) in order to expand into pairs of mutually commuting pure boosts and rotations (this approach has been demonstrated also in [28]).
Another possibility is to work with fixed a 1 and define a 2 = a 1 × R a 1 , thus satisfying the condition r 21 = g 21 as long as a 1 is not an eigenvector of R (or one has a singular solution): Take for example c = (1 + 2i, 2i − 1, −1) t and a 1 = (1, −1, 0) t .Then τ 1 = 1 and τ 2 a 2 = 1 3 (2i + 1, 2i + 1, 4i − 1) t so the SO(3, 1) factorization (37) yields for Λ(c) the result where we use (28) for the second equality to show the boost-rotation interactions.Our next example is borrowed from [35], where we propose and investigate a generalization of the classical Wigner decomposition.Consider a factorization of proper Lorentz transformations into a boost and rotation, for which a convenient vector-parametric solution yields (cf.[35]) Let us now take a specific example given by the matrix where we use (28), so in the end Λ(c • ) and Λ(α) are conjugated with a (symmetric) boost.Our last example is numerical and concerns a Euler-type decomposition into four factors for which the additional parameter ψ allows for optimization of the path length in SO(3) or avoiding the 'gimbal lock' singularity.If we choose the unit vectors a i , n ∈ S 2 expressed as in spherical coordinates and the angle of compound rotation to be φ = 111 • , a generalized version of Formula (7) yields our solution in the form (rounding here is just for convenience) Obviously complex angles correspond to loxodromic Lorentz transformations with parallel rotation and boost counterpart, so one may easily factorize further using ( 11) and ( 12), but the interesting feature here is that ψ is real, so we effectively control the relativistic effect (represented by the imaginary part of R) using a pure rotation in the sequence changing the value of ψ results in a generic Lorentz transformation instead of the special orthogonal matrix R(n, φ).In other words, inserting an appropriately chosen rotation in the sequence of SO + (3, 1) mappings eliminates the net relativistic effect from the observer's point of view.

Discussion
The present study suggests a consistent geometric approach with applications in various branches of physics and engineering.Specific examples are given but without too many technical details that would distract the main focus.One particular benefit is in the context of matrix factorizations, such as the Iwasawa decomposition for PSL(2, R) or the polar decomposition for the monodromy matrix considered above.These results generalize easily to SO(4, C) and its other real forms like SO (2,2) or SO * (4), and in the Plücker setting allow for higher-dimensional embeddings as well (once again, we shall refer to [28] for details).
The method used here has many advantages we attempt to emphasize with examples.Firstly, it emerges directly from the Clifford algebra and is thus easily interpretable in geometric terms.Secondly, it provides a quite natural parametrization of the rotation and Lorentz groups, preserving the projective space topology, unlike Euler angles for instance (which use the torus for a model).Hence, no coordinate singularities appear in this representation.Thirdly, it is quite efficient both computation-wise and from a theoretical perspective: no redundant parameters are required and numerical algorithms run much faster (compared for instance with the usual matrix multiplication), while at the same time, the simplicity of the analytic expressions (typically in terms of rational functions) allows for obtaining solutions to non-trivial problems in a casual manner.Lie group factorizations, geometric phases in kinematics and the relativistic extension via complexification are just a few examples to support this.In the present paper we focus mainly on the latter, as it yields deep geometric insight for both theoretical speculations based on analogies and practical problem-solving techniques in engineering, but at the same time, has been neglected in other research papers that aim only at obtaining the end result for a given specific task.Our goal here is to fill this gap (at least partially), discussing mostly complex extensions of already existing solutions in the real setting and their interpretation in the more general Lorentzian context.On the engineering front, this provides plenty of useful applications, such as manipulating the outcome of measurable relativistic effects using pure rotations and vice versa, as well as analogous results in optics and micro-electronics, quantum computation and scattering of elementary particles.Since these solutions have been obtained in other studies, we rely extensively on them without providing detailed derivations and proofs, but that is the price one has to pay for keeping the volume within reasonable limits.
use(5) to obtain c = (4, 3, 4i) t and according to (39) arrive at above allows us to express the pure rotation Λ(c • ) in the form