A signature invariant geometric algebra framework for spacetime physics and its applications in relativistic dynamics of a massive particle and gyroscopic precession

A signature invariant geometric algebra framework for spacetime physics is formulated. By following the original idea of David Hestenes in the spacetime algebra of signature \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(+,-,-,-)$$\end{document}(+,-,-,-), the techniques related to relative vector and spacetime split are built up in the spacetime algebra of signature \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(-,+,+,+)$$\end{document}(-,+,+,+). The even subalgebras of the spacetime algebras of signatures \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(\pm ,\mp ,\mp ,\mp )$$\end{document}(±,∓,∓,∓) share the same operation rules, so that they could be treated as one algebraic formalism, in which spacetime physics is described in a signature invariant form. Based on the two spacetime algebras and their “common” even subalgebra, rotor techniques on Lorentz transformation and relativistic dynamics of a massive particle in curved spacetime are constructed. A signature invariant treatment of the general Lorentz boost with velocity in an arbitrary direction and the general spatial rotation in an arbitrary plane is presented. For a massive particle, the spacetime splits of the velocity, acceleration, momentum, and force four-vectors with the normalized four-velocity of the fiducial observer, at rest in the coordinate system of the spacetime metric, are given, where the proper time of the fiducial observer is identified, and the contribution of the bivector connection is considered, and with these results, a three-dimensional analogue of Newton’s second law for this particle in curved spacetime is achieved. Finally, as a comprehensive application of the techniques constructed in this paper, a geometric algebra approach to gyroscopic precession is provided, where for a gyroscope moving in the Lense-Thirring spacetime, the precessional angular velocity of its spin is derived in a signature invariant manner.

• For two multivectors A and B in spacetime, their geometric product, inner product, outer product, and commutator product are represented by AB, A · B, A ∧ B , and A × B , respectively; • For a multivector M in spacetime, M and �M� p (p = 0, 1, 2, 3, 4) denote its reverse and p-vector part, respectively, where M 0 is abbreviated as M ; • The Greek letters, denoting the spacetime indices, range from 0 to 3, whereas the Latin letters, denoting the space indices, range from 1 to 3; • The sum should be taken over, when repeated indices appear within a term; • The international system of units is used.
Since the even subalgebras of the two STAs share the same operation rules, we will no longer distinguish them strictly and treat them as one algebraic formalism hereafter. In Appendix B of this paper, a detailed presentation of this algebraic formalism is given. It will be shown that the "common" even subalgebra of the STAs of signatures (±, ∓, ∓, ∓) actually provides a signature invariant GA framework for spacetime physics. In order to give an application paradigm of the two STAs and their "common" even subalgebra, we need to make use of them to study some specific problems in spacetime physics, and gyroscopic precession is such a typical topic.
According to the prediction of General Relativity, the spin of a gyroscope precesses relative to the asymptotic inertial frames as it moves around a rotating spherical source 22 . The conventional method to describe gyroscopic precession under the weak-field and slow-motion (WFSM) approximation in tensor language is presented in Refs. 21,22 . For a uniformly rotating spherical source, the external gravitational field is stationary, and only the leading pole moments need to be considered, so that the spacetime geometry is described by the Lense-Thirring metric 30 . As a result, the corresponding spacetime is known as the Lense-Thirring spacetime. When a torquefree gyroscope is moving in this spacetime, there exist three types of precession for its spin, namely, the de Sitter precession, the Lense-Thirring precession, and the Thomas precession, where these phenomena are, respectively, resulted from gyroscopic motion through the spacetime curved by the mass of the source, rotation of the source, and gyroscopic non-geodesic motion 31 .
In the traditional description for gyroscopic precession based on tensor language, one always needs to work with the components of some tensor in a chosen coordinate frame, which often leads to many equations with a low degree of clarity. The language of STA could provide a physically clear approach to dealing with this topic, since one just involves geometric objects during calculation 32 . As a preliminary attempt, another purpose of the present paper is to handle gyroscopic precession by applying the STAs of signatures (±, ∓, ∓, ∓) and their "common" even subalgebra, so that for a gyroscope moving in the Lense-Thirring spacetime, a signature invariant derivation of the precessional angular velocity of its spin could be achieved. For brevity, in later applications, the signs "±" associated with multivectors and operators will be suppressed, and for equalities like A = F(±B) and C = G(∓D) , the signs " + " and "−" in the former equation correspond to the cases in the signatures (+, −, −, −) and (−, +, +, +) , respectively, and the situation in the latter equation is reverse.
Before analyzing gyroscopic precession, rotor techniques on Lorentz transformation and relativistic dynamics of a massive particle in curved spacetime need to be addressed in the two STAs. Rotor techniques are available in the STA of signature (+, −, −, −) [17][18][19][20] , and however, since the STA of signature (−, +, +, +) is rarely employed, these techniques have not been fully developed in this algebraic formalism, where in particular the expressions of the rotors inducing Lorentz boost and spatial rotation should be clearly established. Being the third purpose of this paper, by virtue of the rotors constructed in the "common" even subalgebra of the two STAs, the general Lorentz boost with velocity in an arbitrary direction and the general spatial rotation in an arbitrary plane are handled in a signature invariant manner. How to study physics in curved spacetime based on STA is a fundamental problem. By following GA techniques for General Relativity formulated in Ref. 33 , the treatment of gyroscopic precession in this paper is able to be put on a solid theoretical footing. To generate the STAs of signatures (±, ∓, ∓, ∓) in a curved spacetime, one just needs to define a local orthonormal tetrad {γ α } by the orthonormalization of a coordinate frame (in either signature), and then, by applying these two STAs and their "common" even subalgebra, the relevant topics in spacetime physics can be dealt with.
Relativistic dynamics of a massive particle in curved spacetime should be studied so as to describe the motion of a gyroscope moving around a gravitating source 34 . We assume that a collection of fiducial observers Scientific Reports | (2022) 12:3981 | https://doi.org/10.1038/s41598-022-06895-0 www.nature.com/scientificreports/ is distributed over space, and each fiducial observer is at rest in the coordinate system of the spacetime metric. For a massive particle, the spacetime splits of the velocity, acceleration, momentum, and force four-vectors with the normalized four-velocity γ 0 of the fiducial observer need to be derived, which is easy when spacetime is flat. However, in curved spacetime, some subtleties appear and ought to be seriously analyzed. For instance, the proper time of fiducial observers should be identified, and the contribution of the bivector connection ω(u) associated with {γ α } (cf. Ref. 33 ) should also be considered. In this paper, after overcoming these difficulties, the results are given, and with them, a three-dimensional analogue of Newton's second law for the particle in curved spacetime is achieved, which is the fourth purpose of the present paper. Besides, the Fermi-Walker derivatives presented in tensor language are recast in the STAs of signatures (±, ∓, ∓, ∓) so that the motion of the spin of a gyroscope can be depicted in these two STAs 21 .
With the aid of the GA techniques constructed before, an efficient treatment of gyroscopic precession could be provided in the two STAs. Considering a gyroscope moving in the Lense-Thirring spacetime, some significant results like the three-dimensional generalized equation of motion for the gyroscope are first given on the basis of relativistic dynamics of a massive particle. Then, the rotor techniques are employed to handle the spin of the gyroscope, and the direct result shows that a bivector field Ω(τ ) along its worldline completely determines the motion of its spin, where τ is the proper time. The bivector field Ω(τ ) is dependent on the rotor L generating the pure Lorentz boost from the gyroscope's four-velocity u to the fiducial observer's four-velocity cγ 0 and the bivector connection ω(u) associated with {γ α } , where c is the velocity of light in vacuum. Just like the Faraday bivector, namely the electromagnetic field strength, the bivector field Ω(τ ) can also be decomposed into the electric part Ω (E) (τ ) and the magnetic part Ω (B) (τ ) . Let {γ β } be the reciprocal tetrad of {γ α } , and technically, if the condition L aLγ 0 = cΩ (E) (τ ) is fulfilled, the spin of the gyroscope always precesses relative to its comoving frame, determined by the pure Lorentz boost generated by the rotor L , with Ω (B) (τ ) as the precessional angular velocity.
The key point is to write down signature invariant expression of the bivector field Ω(τ ) and the spacetime split of the gyroscope's four-acceleration a with the normalized four-velocity γ 0 of the fiducial observer based on the "common" even subalgebra of the two STAs. According to Refs. 33,35 , the bivector connection ω(u) associated with {γ α } can be directly derived, and then, by recasting it in terms of the relative vectors {σ k } , its signature invariant expression and those of its electric part ω (E) (u) and magnetic part ω (B) (u) are obtained. Moreover, by applying the rotor techniques, the pure Lorentz boost L from u to cγ 0 can also be derived. Thus, as noted before, the signature invariant expression of Ω(τ ) and those of Ω (E) (τ ) and Ω (B) (τ ) are completely determined. As to a, its spacetime split with γ 0 could be directly obtained from the relevant conclusion in relativistic dynamics of a massive particle. Thus, with a, L , and Ω (E) (τ ) , one is capable of verifying that the condition L aLγ 0 = cΩ (E) (τ ) holds by means of various operations in the "common" even subalgebra of the two STAs, and hence, the spin of the gyroscope indeed precesses in the comoving frame with Ω (B) (τ ) as the precessional angular velocity. After expanding Ω (B) (τ ) up to 1/c 3 order with 1/c as the WFSM parameter 36 , the gyroscope spin's angular velocities of the de Sitter precession, the Lense-Thirring precession, and the Thomas precession are able to be read out, and their expressions, in the form of geometric objects, are equivalent to their conventional ones in component form, respectively.
The whole derivation implies that the "common" even subalgebra of the STAs of signatures (±, ∓, ∓, ∓) does provide a signature invariant GA framework for spacetime physics, and the rotors, presented in a signature invariant form, can be used to generate Lorentz transformations in these two STAs. The treatment of relativistic dynamics of a massive particle and gyroscopic precession intuitively displays the basic method of dealing with specific topics in curved spacetime within the signature invariant GA framework, which suggests that the GA techniques established in this paper are efficient and reliable. No doubt, if these techniques are directly applied to gyroscopic precession in alternate theories of gravity, such as f(R) gravity 30,37-39 , f (R, G) gravity 40,41 , and f(X, Y, Z) gravity 42 , they will definitely facilitate the relevant studies, where G is the Gauss-Bonnet invariant, X := R is the Ricci scalar, Y := R µν R µν is the quadratic contraction of two Ricci tensors, and Z := R µνσρ R µνσρ is the quadratic contraction of two Riemann tensors. Furthermore, by developing other types of techniques, the method in this paper could also be applied to more fields, and in fact, some topics in classical mechanics and electrodynamics have been described in such a manner. The applications of this method will be expected to be extended to a wider range in the future, so that the study of spacetime physics in the language of GA could be greatly promoted.
This paper is organized as follows. In "STAs of signatures (±, ∓, ∓, ∓) and their "common" even subalgebra", the STAs of signatures (±, ∓, ∓, ∓) and their "common" even subalgebra are formulated. In "Rotor techniques on Lorentz transformation and relativistic dynamics of a massive particle in curved spacetime", rotor techniques on Lorentz transformation and relativistic dynamics of a massive particle in curved spacetime are constructed. In "A GA approach to gyroscopic precession in the Lense-Thirring spacetime", a GA approach to gyroscopic precession in the Lense-Thirring spacetime is given. In "Summary and discussions", some concluding remarks will be made. In Appendix A, operation rules of blades in the STAs of signatures (±, ∓, ∓, ∓) are summarized. In Appendix B, the "common" even subalgebra of the STAs of signatures (±, ∓, ∓, ∓) is introduced in detail. In Appendix C, a local orthonormal tetrad {γ α } and the bivector connection ω(u) associated with it in the Lense-Thirring spacetime are derived.
STAs of signatures (±, ∓, ∓, ∓) and their "common" even subalgebra STA, introduced in the classical literature Space-Time Algebra by David Hestenes (1966), can provide a synthetic framework for relativistic physics 17 , so it has attracted widespread attention in the physical community. Since the establishment of STA, the signature (+, −, −, −) has been widely used, and however, in relativistic physics, one of the main application fields of STA, the opposite signature (−, +, +, +) is often adopted 17,23 . Thus, when one intends to apply STA to relativistic physics, the change of signature from one to another will cause inconvenience even though these two signatures differ only by a minus sign. In fact, the STA of signature (−, +, +, +) www.nature.com/scientificreports/ was also used [24][25][26][27][28][29] , but a lack of long-term attention to it results in that the techniques related to relative vector and spacetime split have not been developed in this algebraic formalism so that its applications are quite limited. In this section, by following the original idea of David Hestenes, we will build up these techniques in the STA of signature (−, +, +, +) so that a more convenient approach to relativistic physics could be given in the language of GA. For the ease of writing, we will directly formulate the STAs of signatures (±, ∓, ∓, ∓) , and analyze the operation rules of multivectors. In spacetime, the STAs of signatures (±, ∓, ∓, ∓) can be generated by corresponding orthogonal vectors {γ ± α } satisfying respectively, where η ± αβ are the Minkowski metrics in the two signatures. With these vector generators {γ ± α } , explicit bases for both the STAs are defined, namely where, in either signature, one scalar, four vectors, six bivectors, four trivectors, and one pseudoscalar are contained. One can perform operations between any two multivectors in spacetime by expanding them in a basis, once operation rules of blades of different grades are given, where the term "blade" here denotes a multivector written as the outer product of a set of vectors (cf. Ref. 17 ). In Appendix A of this paper, a detail list of operation rules of blades in the two STAs is presented, and based on these rules, the "common" even subalgebra of these two STAs will be constructed in the following.
According to Eqs. (A1) and (A7), the orthogonality between the vector generators {γ ± α } implies that the bases (2) can be rewritten as where the geometric products of {γ ± α } are obviously anticommutative, By making use of the anticommutation of {γ ± α } , the pseudoscalars I ± also have the expressions, with ǫ ijk as the three-dimensional Levi-Cività symbol. Among the basis blades, those of even grade, form bases for the even subalgebras of the two STAs. Now, we will first discuss some properties of the bivectors {γ ± 0 γ ± k } . With Eqs. (1), (4), and (A14), one can directly derive the following equalities, where δ ij is the Kronecker symbol, and in the second step of (8), Eqs. (5), (A5), and (A10) have been used. These equalities show that relative vectors, spanning the relative spaces orthogonal to the timelike vectors γ ± 0 , could be defined as {σ ± k = ∓γ ± 0 γ ± k = γ ± k γ 0 ± } with {γ α ± } as the reciprocal frames of {γ ± α } , so that they have the similar algebraic properties to the Pauli matrices,  www.nature.com/scientificreports/ and then, by inserting Eqs. (13) and (14) into Eqs. (7) and (8), respectively, we get which prove once again that the algebraic properties of {σ ± k } are similar to those of the Pauli matrices. In fact, as mentioned in Ref. 32 , {σ + k } or {σ − k } provide a representation-free version of the Pauli matrices. Equations (10) and (12) show that the relative spaces orthogonal to γ ± 0 are both the Euclidean spaces of dimension 3 with {σ ± k } and I ± as orthonormal bases and pseudoscalars, respectively. In relative space, a relative vector, although being a bivector in STA, is actually treated as a multivector of grade 1, and thus, in this sense, the inner product and the cross product between two relative vectors can be defined. Let a ± = a ± i σ ± i and b ± = b ± j σ ± j be relative vectors, and then, with the help of Eqs. (10) and (11), the inner products and the cross products between a ± and b ± are defined as where the commutator products between a ± and b ± , and have been used. Obviously, the above definitions of inner product and cross product are identical to their conventional ones, respectively. The cross products defined in Eqs. (19) determine the handedness of {σ ± k } , and by applying them, one easily gets which clearly suggest that {σ ± k } are both right-handed bases. Next, we will employ relative vectors to reconstruct bases of the even subalgebras of the STAs of signatures (±, ∓, ∓, ∓) . The definitions of {σ ± k } provide and then, by further using Eqs. (1) and (4), there are After inserting Eqs. (23), (24), and (12) into (6), we know that bases of the even subalgebras of the two STAs can be reconstructed as which indicates that {σ ± k } are actually the vector generators of the two subalgebras. Eqs. (11) and (17) imply that equalities hold, and thus, the anticommutation of {σ ± k }, is explicitly obtained. As a consequence, there exist three types of basic homogeneous multivectors (cf. Ref. 17 ) in the even subalgebras of the two STAs, namely, www.nature.com/scientificreports/ In view of (12), a ± × b ± ∧ c ± in (30) are able to be written in the form of multiplications of the pseudoscalars I ± by real numbers, and in fact, from the bases (25), all multivectors of grade 4 could be expressed in such a form. Thus, Eq. (A5) states that the geometric product between any multivector and a pseudoscalar is equivalent to their inner product. Keep this conclusion in mind, and then, with the help of the following formulas, one gets a convenient way to carry out operations involving multivectors of grade 4, where in the derivation of (31), Eqs. (20) and (21) have been used. Eqs. (23) and (8) show that both σ ± k and σ ± i × σ ± j (i � = j) are bivectors in the two STAs, where the former contain timelike components, whereas the latter do not. The geometric products of them also need to be derived, where according to Eq. (A2), we have By further using Eqs. (1), (4), (A8), and (A15), the terms on the right-hand sides of Eqs. (33), (34), and (35) are achieved, With the aid of the above operation rules of the basic homogeneous multivectors, namely Eqs. (31)- (39), one can carry out operations of any two multivectors in the even subalgebras of the STAs of signatures (±, ∓, ∓, ∓) . Evidently, as shown in these formulas, the two even subalgebras share the same operation rules, and thus, when dealing with specific problems, such as relativistic dynamics of a massive particle and gyroscopic precession in the next two sections, we will no longer distinguish them strictly and treat them as one algebraic formalism. In Appendix B of the present paper, a detailed presentation of this "common" even subalgebra of the two STAs is given. It will be shown that this algebraic formalism provides a signature invariant GA framework for spacetime physics.
When STA is used to describe relativistic physics, the techniques on spacetime split are also of significance, where in the STA of signature (+, −, −, −) , these techniques provide an extremely efficient tool for comparing physical effects in different frames 2,17 . Of course, these techniques can also be constructed in the STA of signature (−, +, +, +) . Let b ± = b α ± γ ± α be vectors in spacetime, and the spacetime splits of b ± with γ ± 0 are defined as where b ± = b i ± σ ± i are called the relative vectors of b ± . Besides, as for operators ∂ ± := γ α ± ∂ α , their spacetime splits with γ ± 0 are given by www.nature.com/scientificreports/ where ∂ µ := ∂/∂x µ and ∇ ± := σ k ± ∂ k with x µ and {σ k ± := γ ± 0 γ k ± } as coordinates in spacetime and the reciprocal frames of {σ ± k } in the relative spaces, respectively. As clearly shown, the spacetime splits of b + and ∂ + are indeed the same as those introduced in the STA of signature (+, −, −, −) 2,17 , and the spacetime splits of b − and ∂ − are those defined in the STA of signature (−, +, +, +).
The timelike vectors cγ ± 0 could be recognized as the four-velocities of some observer, so the spacetime split introduced above is observer dependent, and consequently, one of the most powerful applications of the techniques on spacetime split is that they can greatly simplify the study of effects involving different observers 2,17 . Technically, spacetime split actually encodes the crucial geometric relationship between STA and its even subalgebra 2 , where with these techniques, many calculations between vectors in spacetime are able to be transformed into those in the even subalgebra of STA. As a result, based on various operations in this algebraic formalism, a large number of specific problems could be solved efficiently. Moreover, since the even subalgebras of the STAs of signatures (±, ∓, ∓, ∓) share the same operation rules, by resorting to the techniques on spacetime split, one is capable of managing to acquire a signature invariant approach to these problems. We will see that the above advantages of spacetime split play a key role in the following treatment of relevant topics.

Rotor techniques on Lorentz transformation and relativistic dynamics of a massive particle in curved spacetime
It is well known that one of the remarkable advantages of STA is that Lorentz boost and spatial rotation can be handled with rotor techniques in an elegant and highly condensed manner [17][18][19][20] . As shown in classical literatures 21,22 , a knowledge of Lorentz boost and spatial rotation is heavily involved in the description of gyroscopic precession, and hence, it could be expected that a more efficient approach to dealing with this topic will be found in the language of STA. Besides, in "STAs of signatures (±, ∓, ∓, ∓) and their "common" even subalgebra", it is claimed that the "common" even subalgebra of the STAs of signatures (±, ∓, ∓, ∓) provides a signature invariant GA framework for spacetime physics, and thus, when this framework is applied to gyroscopic precession, a signature invariant GA derivation of the precessional angular velocity of the gyroscope spin could be achieved. Therefore, as a preliminary attempt, making use of the two STAs and their "common" even subalgebra to study gyroscopic precession is one objective of the present paper, which, if successful, will definitely become an application paradigm of STA. In view that many relevant techniques need to be constructed in this section, the detailed treatment of gyroscopic precession will be left to the next section.
In the analyse of gyroscopic precession, rotor techniques on Lorentz boost and spatial rotation are widely used, and therefore, these techniques need to be specifically addressed in the two STAs. Rotor techniques are available in the STA of signature (+, −, −, −) [17][18][19][20] , and however, since the STA of signature (−, +, +, +) is rarely employed, these techniques have not been fully developed in this algebraic formalism, where in particular the expressions of the rotors inducing Lorentz boost and spatial rotation should be clearly established. In this section, by constructing the rotors on the basis of the exponential function defined on the "common" even subalgebra of the two STAs, the general Lorentz boost with velocity in an arbitrary direction and the general spatial rotation in an arbitrary plane are handled in a signature invariant manner. In addition, relativistic dynamics of a massive particle in curved spacetime ought be studied so as to describe the motion of a gyroscope moving around a gravitating source 34 . To this end, for a massive particle, the spacetime splits of the velocity, acceleration, momentum, and force four-vectors with the normalized four-velocity of the fiducial observer, at rest in the coordinate system of the spacetime metric, are first derived, and then with these results, a three-dimensional analogue of Newton's second law for this particle in curved spacetime is achieved. Furthermore, in order to describe the motion of the spin of a gyroscope, the Fermi-Walker derivative in the STA of signature (−, +, +, +) is also constructed by following the way in the (+, −, −, −) signature.
In Appendix B of this paper, the signs "±" associated with multivectors have been omitted in the "common" even subalgebra of the STAs of signatures (±, ∓, ∓, ∓) , so that all the formulas in this algebraic formalism are presented in a neat form. Inspired by this, when formulas in the two STAs are involved hereafter, the following convention will be adopted for brevity: The signs "±" associated with multivectors and operators are suppressed, and for equalities like A = F(±B) and C = G(∓D) , the signs " + " and "−" in the former equation correspond to the cases in the signatures (+, −, −, −) and (−, +, +, +) , respectively, and the situation in the latter equation is reverse.
Rotor techniques on Lorentz boost and spatial rotation . In GA, a rotor R is defined as an even multivector satisfying RR = 1 and the property that the map defined by b → RbR transforms any vector into another one 17 . Rotors encode an important geometric object and can provide a more elegant scheme for performing orthogonal transformations in spaces of arbitrary signature, where mathematically, rotor group, formed by the set of rotors, provides a double-cover representation of the connected subgroup of the special orthogonal group. In the present paper, we are only interested in rotors in spacetime, and in such a case, the rotor group in spacetime is a representation of the group of proper orthochronous Lorentz transformations 17 .
In the STA of signature (+, −, −, −) , rotor techniques on Lorentz boost and spatial rotation have been established [17][18][19][20] , which greatly promotes the application of STA in spacetime physics. Of course, in order to complete the necessary discussion on gyroscopic precession in a signature invariant manner, these techniques also need to be explicitly constructed in the STA of signature (−, +, +, +) . To facilitate the writing, as in "STAs of signatures (±, ∓, ∓, ∓) and their "common" even subalgebra", we will directly build up rotor techniques in the STAs of signatures (±, ∓, ∓, ∓).
In Appendix B of this paper, a simple method to construct rotor is presented, and it has been shown that for a real number α and a unit 2-blade B, e αB is a rotor. Here, we will make use of e αB to handle Lorentz boost and spatial rotation in the two STAs. From Eqs. (8) and (23) Clearly, both σ k and σ i × σ j (i � = j) are unit 2-blades, and the signs of their squares are different, which suggests that there are two types of unit 2-blades in spacetime. It is based on the exponential functions of these two types of unit 2-blades that the rotors inducing Lorentz boost and spatial rotation can be constructed. Let v = v k σ k , m = m i σ i , and n = n j σ j be three arbitrary relative vectors. Consider the bivectors v and m × n , and the following results can be easily given by means of Eqs. (42a)-(43b): and The former two equations indicate that both v and m × n are 2-blades, and thus, with the latter two equations, two unit 2-blades are derived, where a direct calculation verifies that According to Ref. 17 , a proper orthochronous Lorentz transformation can be generated by a rotor R in spactime, and under this transformation, a general multivector M will be transformed double-sidedly as M → R −1 MR . Let θ and ϕ be two real numbers, and the corresponding rotors associated with e v and I 2 are constructed as e θ 2 e v and e ϕ 2 I 2 , respectively. When they act on vectors x and y, two new vectors x ′ and y ′ are obtained, In order to analyze the generated Lorentz transformations in the "common" even subalgebra of the two STAs, the techniques on spacetime split need to be applied. From Eqs. (44a)-(46b) and (A1), the orthogonality and anticommutation of {γ α } imply and then, with the help of Eq. (B39), one gets The spacetime splits of x, y, x ′ , and y ′ with γ 0 are provided by applying Eq. (40),  www.nature.com/scientificreports/ As stated in "STAs of signatures (±, ∓, ∓, ∓) and their "common" even subalgebra", the relative space is an Euclidean space of dimension 3, and a relative vector, although being a spacetime bivector, could be treated as a multivector of grade 1, which implies that in terms of the three-dimensional geometric meaning, a relative vector is just a vector 17 . Similarly, the commutator product of two relative vectors, referred to as the relative bivector, also has three-dimensional geometric meaning. After comparing Eq. (B8) with Eq. (A1), one is able to find that the commutator product of two relative vectors serves as the role of the wedge product of two vectors in general finite dimensional GA, and thus, in the three-dimensional relative space, it encodes an oriented plane 3,17 . In such a sense, Eqs. (57a) and (57b) indicate that x � (x ′ � ) and x ⊥ (x ′ ⊥ ) are, respectively, the components of x (x ′ ) parallel and perpendicular to e v .
Of course, I 2 also defines an oriented plane in the relative space. Let and then, together with Eqs. which explicitly show that y � (y ′ � ) and y ⊥ (y ′ ⊥ ) are, respectively, the components of y (y ′ ) parallel and perpendicular to the plane defined by I 2 .
When the relative vectors x, x ′ , y , and y ′ in Eqs. (52a) and (52b) are replaced by their decompositions, namely Eqs. (56a)-(56d), it will be seen that a clear physical explanation of the Lorentz transformations induced by e v and I 2 in Eqs. (48a) and (48b) is able to be achieved. To this end, the following properties of the components of x, x ′ , y , and y ′ need to be first derived by means of the combination of Eqs. In order to handle these two equations, e −θ e v and e −ϕI 2 should be rewritten as  (64a) and (64b), respectively, and then, by using the grade operator �· · · � and the orthogonal projection operator successively, we finally arrive at In the above derivation, Eqs. (57a) and (B8) have been employed, and besides, one also needs to note that in view of Eqs. Here, in order to reasonably interpret relevant equations obtained in this subsection, the active view for Lorentz transformation needs to be adopted [5][6][7] . Moreover, it also needs to be stressed that for the spatial rotation, Eq. (71b) shows that if ϕ > 0 , the relative bivector y � × y ′ � has the same orientation as I 2 in the threedimensional geometry. Let us recall that the relative vectors v, m , and n were chosen arbitrarily in the beginning, and therefore, with the rotors e θ 2 e v and e ϕ 2 I 2 , the general Lorentz boost with velocity in an arbitrary direction and the general spatial rotation in an arbitrary plane can be handled. Furthermore, considering that Eqs. (67a), (67b), (70), (71a), and (71b) are derived in the "common" even subalgebra of the STAs of signatures (±, ∓, ∓, ∓) , all of these equations are presented in a signature invariant form.
According to the previous discussion, the Lorentz boost and the spatial rotation have been first generated in Eqs. (48a) and (48b), and however, until these two equations were transformed into those in the "common" even subalgebra of the two STAs, their physical explanations were achieved in the three-dimensional geometry. In this process, the techniques on spacetime split have been employed, which implies that the intuitive pictures formed in the relative space are observer dependent. In addition, one may also have found that it is since the (66a) β : = tanh θ , www.nature.com/scientificreports/ "common" even subalgebra of the two STAs are independent of the signatures that the original equation (48a) or (48b) has the same three-dimensional meaning in the two signatures, and thus, a signature invariant method for handling Lorentz boost and spatial rotation is gained. In fact, many topics in spacetime physics can be dealt with in such a manner, and inspired by this, we will apply this method to studying gyroscopic precession in the next section, so that a signature invariant GA derivation of the precessional angular velocity of the gyroscope spin could be found.
As the final task of this subsection, the pure Lorentz boost (cf. Ref. 17 ) between two vectors of the same magnitude will be discussed based on the previous results. Assuming that x ′ = cγ 0 , Eqs. x is mapped to x ′ by According to Ref. 17 , the above L in the (+, −, −, −) signature is exactly the rotor that determines the pure Lorentz boost between x and x ′ , and motivated by this, we claim that the above L in the (−, +, +, +) signature also plays the same role. It should be noted that the validity of Eq. (76) is able to be directly verified only by Eqs. (73a) and (75), which does not depend on the selection of the frame {γ α } . In the treatment of gyroscopic precession in the next section, Eqs. (75) and (76) will be used to generate the pure Lorentz boost between a comoving orthonormal frame of the gyroscope and a local orthonormal tetrad at rest in the coordinate system of the spacetime metric, which greatly improves the computational efficiency.
Relativistic dynamics of a massive particle in curved spacetime. As mentioned previously, the description of the motion of a gyroscope requires that relativistic dynamics of a massive particle in curved spacetime should be studied 34 , and to this end, a brief introduction to relevant GA techniques for General Relativity formulated in Ref. 33 needs to be given, so that the treatment of gyroscopic precession in the following can be put on a solid theoretical footing. In order to develop a GA description of curved spacetime, one should define a local orthonormal tetrad {γ α } by the orthonormalization of a coordinate frame and then generate the corresponding STA. Let x µ and {g µ } be local coordinates in a curved spacetime and the associated coordinate frame, respectively. Assume that a collection of fiducial observers is distributed over space, and each fiducial observer is at rest in the coordinate system. Then, the components of the metric with respect to the coordinate frame {g µ }, satisfy the conditions 44 (77) g µν : = g µ · g ν , (78a) ±g 0 · g 0 = ± g 00 > 0, (78b) −(g 1 ∧ g 0 ) · (g 0 ∧ g 1 ) = − det g 00 , g 01 g 10 , g 11 > 0, (78c) ±(g 2 ∧ g 1 ∧ g 0 ) · (g 0 ∧ g 1 ∧ g 2 ) = ± det g 00 , g 01 , g 02 g 10 , g 11 , g 12 g 20 , g 21 45 , where one of its important properties is that it will reduce to ∂ when acting on scalar functions. Suppose that ∇ is the unique torsion-free and metric-compatible derivative operator. Then, according to Ref. 33 , the covariant derivative of a multivector A along a vector b is evaluated by the formula Here, the operator b · ∂ satisfies with {γ β } and φ as the reciprocal tetrad of {γ α } and a scalar field in spacetime, respectively. ω(b) , being the bivector connection associated with {γ α } , is defined by where if b = b µ g µ , the expression of ω(b) is given by 33,35 with With the aid of the corresponding GA technique 3 , {g ν } , as the reciprocal frame of {g µ } , is constructed as where from Eq. (79), the coordinate frame {g µ } can be expanded in the local orthonormal tetrad {γ α }, g 00 , g 01 , g 02 , g 03 g 10 , g 11 , g 12 , g 13 g 20 , g 21 , g 22 , g 23 g 30 , g 31 , g 32 , g 33 . www.nature.com/scientificreports/ Because only the knowledge of covariant derivative and bivector connection will be involved in the discussion of gyroscopic precession, other GA techniques for General Relativity will not be covered here, and the reader wishing to go into more details may consult Ref. 33 . Next, for a massive particle, the spacetime splits of the velocity, acceleration, momentum, and force fourvectors with the normalized four-velocity of the fiducial observer will be discussed, so that relativistic dynamics of this particle in curved spacetime can be studied. Let us first identify the proper time of fiducial observers. As indicated earlier, fiducial observers are at rest in the coordinate system x µ , which means that their worldlines are the coordinate curves with x i = const. (i = 1, 2, 3) , namely, t := x 0 /c coordinate curves. As a consequence, if we let t 0 denote the proper time of each fiducial observer, ±c 2 (dt 0 ) 2 = g 00 c 2 (dt) 2 hold along his worldline, and then, Assuming that x µ (τ ) is the worldline of a massive particle with τ as the proper time, the four-velocity of the particle can be rewritten as 22,30 We will prove that Consider an event P on the particle's worldline. The t coordinate curve with x i = x i (P) (i = 1, 2, 3) passes through P and is the worldline of a fiducial observer. Based on the orthonormal tetrad {γ µ | x i =x i (P) } carried by this fiducial observer, his proper reference frame can be defined, and thus, a local coordinate system y 0 =: ct 0 , y 1 , y 2 , y 3 covering a finite domain near his worldline can also be defined. In this coordinate system, if the worldline of the particle is y µ (τ ) , its four-velocity at the event P is Comparing Eq. (92) with Eq. (90), we get P is an arbitrary event on the particle's worldline, and due to Eq. (89), does not depend on the selection of the coordinate system y µ , so Eq. (91) holds. By applying Eq. (40), the spacetime split of the four-velocity of the particle with γ 0 yields where because of ( uγ 0 ) · (uγ 0 ) = ±u 2 = c 2 , one is able to achieve Since cγ 0 could be identified as the four-velocity of some fiducial observer, u is actually the relative velocity measured in his orthonormal tetrad, which is also able to be inferred from Eq. (92).
(89) dt 0 dt = ±g 00 . www.nature.com/scientificreports/ After clarifying many concepts, we are in a position to derive the spacetime split of the four-acceleration of the particle with γ 0 , which is an essential ingredient in the formalism of relativistic dynamics. The four-acceleration of the particle, a = Du/dτ = u · ∇u , is immediately gained from Eq. (82), and then, by employing Eq. (40), its spacetime split with γ 0 is provided, The first term is in which, Eqs.  Let m be the rest mass of the particle. The spacetime splits of its four-momentum p = mu and the four-force f = Dp/dτ = u · ∇p acting on it also need to be evaluated so that a three-dimensional analogue of Newton's second law in curved spacetime will be achieved. Starting from Eq. (94), the spacetime split of the particle's four-momentum p with γ 0 is where are the energy and the relative momentum of theparticle measured by the fiducial observer (cf. Ref. 17 ), respectively. The relationship between E and p can be directly obtained from ( pγ 0 ) · (pγ 0 ) = m 2 c 2 , which is exactly the same as that in Special Relativity. Assuming that the particle's rest mass remains unchanged as it moves, namely dm/dτ = 0 , the four-force f acting on it is able be expressed as When the spacetime is flat and x µ are coordinates in an inertial frame of reference with g µν = η µν , by definition, fiducial observers reduce to inertial observers. In such a case, Eq. (89) suggests that dt 0 = dt , and the relative force f = f i σ i acting on the particle should be given by f i = dp i /dt 46 . Thus, using σ i = γ i γ 0 , one is capable of recasting f as In curved spacetime, we claim that the corresponding relative force f measured by the fiducial observer is related to Dp/dt 0 in the same way,  (108) is a three-dimensional analogue of Newton's second law in curved spacetime, which constitutes the core content of relativistic dynamics of a massive particle. In the above discussion, the key point is that the relative velocity, relative acceleration, relative momentum, and relative force for the particle could be reasonably defined in the orthonormal tetrad carried by the fiducial observer. Evidently, in terms of the (102c) (110) www.nature.com/scientificreports/ three-dimensional geometric meaning in the relative space, these relative vectors ought to be interpreted as their corresponding three-vectors in tensor language. When the spacetime is flat, the bivector connection ω(u) and its electric part ω (E) (u) and magnetic part ω (B) (u) vanish. In this case, via considering the components of these relative vectors in the rest frame of the fiducial observer, namely {σ k } , one is able to verify that all the above results reduce to those in Special Relativity. Therefore, the formalism of relativistic dynamics of a massive particle constructed in this subsection is an elegant generalization of the classical one in flat spacetime.
In the tetrad formalism of General Relativity 47 , the covariant derivative of a vector b = b α γ α along the coordinate frame vector g µ is given by where are the spin connection coefficients, and due to the metric compatibility condition, they satisfy 48 Using Eqs. (84) and (A14), one obtains which means that the bivector connection ω(g µ ) can be expressed as The above discussion suggests that it could be expected that when the relative vectors in Eq. (108) are expanded in the frame {σ k } , the corresponding generalization of Newton's second law in the tetrad formalism will also be acquired. Compared with those results in the tetrad formalism, the results in this paper are presented in the form of geometric objects, so they are endowed with a higher degree of clarity. Besides, as highlighted before, since the operations in the "common" even subalgebra of the STAs of signatures (±, ∓, ∓, ∓) are independent of the signatures, the relevant results like Eqs. (95), (106), (108), and (109) are able to be handled in a signature invariant manner. As a primary application of the signature invariant GA framework provided by the "common" even subalgebra of the two STAs, the treatment of relativistic dynamics of a massive particle in this subsection provides a paradigm on how to achieve a signature invariant approach to spacetime physics in curved spacetime.
In order to depict the motion of the spin of a gyroscope, the behaviors of vector fields along the worldline of the particle also need to be studied, and here, we only focus our attention on the Fermi-Walker derivatives in the (±, ∓, ∓, ∓) signatures. In fact, their classical forms written in tensor language have been available in Refs. 49,50 , and recasting them in the STAs of the two signatures is a straightforward task. Hence, the results are directly provided as follows: The Fermi-Walker derivatives of a vector field p(τ ) along the particle's worldline in the STAs of signatures (±, ∓, ∓, ∓) are where if D F p(τ )/dτ = 0 , the vector field p(τ ) is said to be Fermi-Walker transported along the particle's worldline. For a torque-free gyroscope moving in spacetime, any nongravitational forces acting on it are applied at its center of mass, and in this case, the spin of the gyroscope experiences the Fermi-Walker transport along its worldline 21 . In the next section, we will regard the transport equation satisfied by the gyroscope spin as the starting point for the discussion of gyroscopic precession. Interestingly, by means of the Leibniz rule and the formula 3 with B as a bivector in spacetime, the above forms of Fermi-Walker derivative can readily be extended to a multivector field A(τ ) along the worldline of the particle, namely, and readers who are interested in this conclusion could attempt to prove it.

A GA approach to gyroscopic precession in the Lense-Thirring spacetime
According to the prediction of General Relativity, the spin of a gyroscope precesses relative to the asymptotic inertial frames as it moves around a rotating spherical source 22 . Conventionally, by following the standard method in tensor language 21,22 , the precessional angular velocity of the gyroscope spin is able to be evaluated under the WFSM approximation. In General Relativity, the time-dependent metric, presented in the form of multipole expansion, for the external gravitational field of a spatially compact supported source is derived under the WFSM approximation in Ref. 30 . Since we are only interested in uniformly rotating spherical sources like the Earth in this paper, the spacetime is stationary, and only the leading pole moments of the source need to be considered. Consequently, in such a case, the metric reduces to the Lense-Thirring metric 30 , and the spacetime is (114) ω µαβ = g µ · ∇γ β · γ α = ω(g µ ) · γ β · γ α = ω(g µ ) · γ β ∧ γ α ), www.nature.com/scientificreports/ accordingly known as the Lense-Thirring spacetime. When a torque-free gyroscope is moving in this spacetime, there exist three types of precession for its spin, namely, the de Sitter precession, the Lense-Thirring precession, and the Thomas precession, where these phenomena are, respectively, resulted from gyroscopic motion through the spacetime curved by the mass of the source, rotation of the source, and gyroscopic non-geodesic motion 31 . Today, the type of experiments designed according to these effects of gyroscopic precession have become an important method to test gravitational theories.
In the traditional description for gyroscopic precession based on tensor language, since one always needs to work with the components of some tensor in a chosen coordinate frame, many equations are given a low degree of clarity. In the language of STA, it could be expected that a physically clear approach to handling this topic will be found, since one just involves geometric objects during calculation 32 . In this section, as a comprehensive application of the STAs of signatures (±, ∓, ∓, ∓) formulated in "STAs of signatures (±, ∓, ∓, ∓) and their "common" even subalgebra" and the GA techniques constructed in "Rotor techniques on Lorentz transformation and relativistic dynamics of a massive particle in curved spacetime", a GA approach to gyroscopic precession will be provided, where for a gyroscope moving in the Lense-Thirring spacetime, the precessional angular velocity of its spin will be derived in a signature invariant manner. The GA description of curved spacetime and the relevant GA techniques for General Relativity introduced at the beginning of "Relativistic dynamics of a massive particle in curved spacetime" will still be adopted, and here, we let x µ and {g µ } be local coordinates in the Lense-Thirring spacetime and the associated coordinate frame, respectively. In addition, it should be pointed out that some physical quantities in this section and Appendix C are presented in the form of the 1/c expansion, where 1/c is used as the WFSM parameter 36 . Since the Lense-Thirring metric is only expanded up to 1/c 3 order, the framework of the linearized General Relativity is sufficient to analyze gyroscopic precession 30,[37][38][39] , and in such a case, the coordinates (x µ ) =: (ct, x i ) are treated as though they were the Minkowski coordinates in flat space 51,52 .
Consider a torque-free gyroscope moving in the Lense-Thirring spacetime, and denote x µ (τ ) as its worldline with τ as the proper time. Assuming that the four-force acting on the gyroscope is f, from Eq. (107), its fouracceleration a is determined by with m as its rest mass. In fact, Eq. (119) should be derived from the Mathisson-Papapetrou-Tulczyjew-Dixon (MPTD) equations, where the term related to the curvature tensor has been omitted because the gyroscope scale is very much smaller than the characteristic dimensions of the gravitational field 34 . In accordance with Refs. 21,22 , the spin s of the gyroscope (i.e., its angular momentum vector) is always orthogonal to its four-velocity u and experiences Fermi-Walker transport along its worldline, It will be seen that starting from the above three equations, the precessional angular velocity of the gyroscope spin can be derived. Besides, gyroscopic precession can also be discussed based on MPTD equations, and interested readers may consult Refs. 53,54 . Since the four-velocity of the gyroscope satisfies by use of Eqs. (A6) and (A7), Eq. (121) is equivalent to Thus, Eqs. (120) and (123) directly result in which means that s 2 remains fixed along the worldline of the gyroscope.
As shown in Appendix C, the Lense-Thirring metric satisfies Eqs. (78a)-(78d), which implies that we are capable of assuming that there exists a collection of fiducial observers who are distributed over space and at rest in the coordinate system x µ , and as a consequence, a local orthonormal tetrad {γ α } in the Lense-Thirring spacetime could be directly defined by means of the corresponding formulas in "Relativistic dynamics of a massive particle in curved spacetime". Based on the detailed calculation in Appendix C, the tetrad {γ α } determined up to 1/c 3 order is given by where the potentials U and U i are, respectively, (122) u 2 = ±c 2 ⇒ u · a = 0, (123) u · ∇s = ∓ 1 c 2 (a · s)u.
(124) ds 2 dτ = u · ∇s 2 = 2s · (u · ∇s) = 0, www.nature.com/scientificreports/ Here, G is the gravitational constant, M and J are the mass and the conserved angular momentum of the gravitating source, respectively, and r := √ x i x i . Before analyzing the motion of the spin s of the gyroscope, its relativistic dynamics needs to be discussed. Let t 0 be the proper time of the fiducial observer, which is related to the coordinate time t by Eq. (89), and from Eq. (C1), the expression of dt 0 /dt up to 1/c 3 order is As in Eqs. (90) and (91), the four-velocity u of the gyroscope can be expanded in the tetrad {γ α }, and then, Eq. (94) indicates that its spacetime split with γ 0 yields where u := u i σ i is the relative velocity measured in the orthonormal tetrad of the fiducial observer. Due to u 2 = ±c 2 , the Lorentz factor γ u has the expression (95), and thus, by expanding it up to 1/c 3 order, one gets Furthermore, based on Eqs. (103) and (108)-(110), the spacetime splits of the four-acceleration of the gyroscope and the four-force acting on it are able to be given, respectively, and in view of Eq. (119), we only give the result of the four-force, In the Lense-Thirring spacetime, after inserting Eqs. (C14) and (C15) into Eqs. (108) and (109), the expressions of the relative force f exerted on the gyroscope and the corresponding power f · u delivered by it up to 1/c 3 order are derived, and with ∇ := σ k ∂ k and V := V i σ i . It could be verified that these two equations are compatible. By plugging the potential U into Eq. (132), one will find that −m∇U is the Newtonian gravitational force acting on the gyroscope, and hence, at the leading order, Eqs. (132) and (133) reduce to the corresponding results in Newtonian gravity, which means that Eq. (132) is a three-dimensional analogue of Newton's second law for the gyroscope in the Lense-Thirring spacetime. Evidently, the terms at the next-leading order fall into three classes that depend on U, V , and f , respectively, and as implied from Eq. (126), they should be resulted from gyroscopic motion through the spacetime curved by the mass of the source, rotation of the source, and gyroscopic non-geodesic motion. It will be seen that due to the same reasons, the spin of the gyroscope also experiences three types of precession. In Eqs. (132) and (133), the corrections to the results in Newtonian gravity are presented in a very elegant way, which intuitively displays the powerful potential of the signature invariant GA framework formulated in "STAs of signatures (±, ∓, ∓, ∓) and their "common" even subalgebra" for application in spacetime physics.
Next, we begin to review the basis process of evaluating the precessional angular velocity of the gyroscope spin in the language of STA. Let {γ (α) } be a local orthonormal frame comoving with the gyroscope, and by definition, the timelike vector γ (0) is given by γ (0) = u/c . In order to determine the other three spacelike vectors γ (i) of {γ (α) } , the pure Lorentz boost between the gyroscope's four-velocity u and the fiducial observer's four-velocity cγ 0 needs to be presented. According to Eqs. (75) and (76), under the pure Lorentz boost generated by the rotor the vector u is mapped to cγ 0 by (131) (138) LL =LL = 1.  This result indicates that the motion of the spin of the gyroscope relative to the comoving frame {γ (α) } is completely determined by the bivector field Ω(τ ) along its worldline, where Ω(τ ) is dependent on the rotor L generating the pure Lorentz boost from the gyroscope's four-velocity u to the fiducial observer's four-velocity cγ 0 and the bivector connection ω(u) associated with the tetrad {γ α } . Like the bivector connection ω(u) in Eqs. (100a)-(101c), the bivector Ω(τ ) is also able to be decomposed into the electric part Ω (E) (τ ) and the magnetic part Ω (B) (τ ) , which clearly suggests that −Ω (B) (τ )I , as a relative vector, is the precessional angular velocity of s ′ in the conventional sense, and because the cross product (denoted by × 3 ) is rarely employed in GA, the relative bivector Ω (B) (τ ) could be regarded as the precessional angular velocity of s ′ . That is to say, in the comoving frame {γ (α) } of the gyroscope, its spin always precesses with Ω (B) (τ ) as the precessional angular velocity. In addition, one should also note that since (164) or (166) has been represented in the "common" even subalgebra of the STAs of signatures (±, ∓, ∓, ∓) , a signature invariant GA derivation of the precessional angular velocity of the gyroscope spin could be found. In order to make a further analysis on Eq. (164), we need to derive the expressions of Ω (E) (τ ) , Ω (B) (τ ) , and a ′ . Let us first evaluate the corresponding results of Ω (E) (τ ) and Ω (B) (τ ) , and as shown in Eqs. (161a)-(161c), their expressions can directly be read out from that of the bivector Ω(τ ) . By plugging Eqs. (136) and (100c) into Eq. (160), the bivector Ω(τ ) is able to be expressed as (172a)  LT (τ ) are resulted from gyroscopic motion through the spacetime curved by the mass of the source and rotation of the source, respectively, and hence, they should be the de Sitter precession and the Lense-Thirring precession 22 . Besides, Ω (B) T (τ ) , associated with the relative force f acting on the gyroscope, explicitly represents the Thomas precession of its spin, which is caused by gyroscopic non-geodesic motion. In the fine structure of atomic spectra, Thomas precession plays a significant role 21 .
Recall that the three-dimensional operator ∇ appearing in Eq. (182) is defined by ∇ = σ k ∂ k (cf. (41)), and on the basis of it, the expression of Ω where r := x i σ i is the relative position vector of the gyroscope, and due to r = √ x i x i , there is r = √ r 2 . In order to deduce the expression of Ω (B) LT (τ ) , some tricks need to be applied. In the language of GA, the relative angular momentum bivector J , is more convenient to describe the rotation of the source. Equations (126) and (C1) indicate that the source is rotating around the x 3 axis, so its relative angular momentum vector is and then, from Eqs. (B14), (B20), (B2), and (B6), its relative angular momentum bivector should be Thus, via Eqs. (126) and (B15), one is able to express V as Keeping in mind that the relative angular momentum bivector J of the source is conserved, the following identity holds,  T (τ ) are capable of being directly transformed into their corresponding expressions in the conventional sense by multiplying −I , namely, Although these expressions presented here seem to be identical to those in Refs. 30,55 , one still needs to note that since the relative velocity u in −Ω T (τ )I is measured in the orthonormal tetrad {γ α } of the fiducial observer instead of in the coordinate frame {g µ } , their above expressions are slightly different from those obtained in tensor language. In despite of this, a straightforward calculation 22,30 shows that the difference between the gyroscope's velocities measured in {γ α } and in {g µ } is at least at 1/c 2 order, so the above −Ω T (τ )I are essentially equivalent to their conventional expressions. These computations in the final part of this section display in detail how to give a signature invariant GA derivation of the precessional angular velocity of the gyroscope spin within the framework provided by the "common" even subalgebra of the STAs of signatures (±, ∓, ∓, ∓) , which could stand as a successful paradigm of the application of this framework in spacetime physics.
In this section, based on the STAs of signatures (±, ∓, ∓, ∓) formulated in "STAs of signatures (±, ∓, ∓, ∓) and their "common" even subalgebra" and the GA techniques constructed in "Rotor techniques on Lorentz transformation and relativistic dynamics of a massive particle in curved spacetime", an efficient treatment of gyroscopic precession is achieved. One significant advantage of GA approach is that only geometric objects are involved during calculation, and thus, many equations are given a degree of clarity which is lost in tensor language. A typical example is that the relationship between the gyroscope spin s and its components s (i) in the comoving frame {γ (α) } is clearly shown by the equation s (i) = s · γ (i) , which could help readers understand that instead of s, it is the spin s (1) , s (2) , s (3) in the frame {γ (α) } that experiences a spatial rotation. However, in the classical derivation with tensor, since one always needs to work with the components of some tensor, the role of s is usually played by its components in the coordinate frame {g µ } , and thus, the above equation is replaced by the corresponding component equations 22,30 , from which, the relationship between s and s (i) can not be explicitly reflected.
It should be noted that the application of the rotor techniques is also very crucial in simplifying the derivation. In the beginning, Eqs. (142)-(144) imply that in order to obtain the precessional angular velocity of the gyroscope spin s (1) , s (2) , s (3) in the frame {γ (α) } , the expression of ds (i) /dτ needs to be given. Then as in Eq. (146), by employing the rotor techniques, the effect of the pure Lorentz boost generated by the rotor L is transformed from γ (i) to s ′ , and as a result, one can deal with the geometric object ds ′ /dτ = ds (i) /dτ γ i rather than ds (i) /dτ . Being a common trick in STA, such an approach is extremely useful for computations. The STAs of signatures (±, ∓, ∓, ∓) and the GA techniques for General Relativity formulated in Ref. 33 are organically integrated in "Relativistic dynamics of a massive particle in curved spacetime", so that physics in curved spacetime is able to be discussed within the signature invariant framework provided in "STAs of signatures (±, ∓, ∓, ∓) and their "common" even subalgebra", which is perhaps the most easily overlooked contribution of the present paper. It is based on the results presented in "Relativistic dynamics of a massive particle in curved spacetime" that relativistic dynamics of the gyroscope and the precession of its spin can be studied in the two STAs. In particular, within the framework provided by the "common" even subalgebra of the two STAs, the three-dimensional generalized equation of motion for the gyroscope and the precessional angular velocity of its spin are able to be derived in a signature invariant manner. The treatment of gyroscopic precession in this section intuitively displays the basic method of dealing with specific problems in curved spacetime within the signature invariant framework. In the future, if the applications of this method could be extended to a wider range, the study of spacetime physics in the language of GA will be greatly promoted.

Summary and discussions
Since the establishment of STA by David Hestenes, the signature (+, −, −, −) has been widely used 2,17 , which may cause inconvenience to the application of STA in relativistic physics because plenty of literatures on relativity adopt the opposite signature (−, +, +, +) . Although the STA of signature (−, +, +, +) was also used 24-29 , a lack of long-term attention to it results in that its applications are quite limited. In this paper, by following the original idea of Hestenes, the techniques related to relative vector and spacetime split are built up in the STA of signature (−, +, +, +) , so that a more convenient approach to relativistic physics could be given in the language of GA. The further research suggests that the two even subalgebras of the STAs of signatures (±, ∓, ∓, ∓) share the same operation rules, so that they could be treated as one algebraic formalism. Consequently, many calculations between vectors involved in a large number of specific problems can be transformed into those in this "common" even subalgebra of the two STAs through the techniques on spacetime split, and then be solved efficiently in a signature invariant manner with the help of various operations provided in Appendix B. Thus, the "common" even subalgebra of the two STAs provides a signature invariant GA framework for spacetime physics. When orthogonal transformations in spaces of arbitrary signature are performed, calculations with rotors are demonstrably more efficient than calculations with matrices, which is a remarkable advantage of GA. Therefore, the topic of rotor techniques on Lorentz transformation should be specifically addressed in the STAs of signatures (±, ∓, ∓, ∓) , and what needs to be pointed out is that since rotor techniques have not been fully developed in the STA of signature (−, +, +, +) , it is significant to explicitly elaborate how to construct the rotors inducing Lorentz boost and spatial rotation in this algebraic formalism. In the present paper, by constructing the rotors on the basis of the exponential function defined on the "common" even subalgebra of the two STAs, the general Lorentz boost with velocity in an arbitrary direction and the general spatial rotation in an arbitrary plane are handled in a signature invariant manner.
Relativistic dynamics of a massive particle in curved spacetime is also studied so as to describe the motion of a gyroscope moving around a gravitating source 34 . To this end, the two STAs and their "common" even subalgebra are first generated by a local orthonormal tetrad, and thus, the corresponding signature invariant GA framework can be set up. Then, after organically integrating the STAs of signatures (±, ∓, ∓, ∓) and the GA techniques for General Relativity formulated in Ref. 33 , physics in curved spacetime is able to be discussed within the signature invariant framework provided in "STAs of signatures (±, ∓, ∓, ∓) and their "common" even subalgebra", which lays the foundation for dealing with gyroscope precession hereafter. With these preparations, for a massive particle, the spacetime splits of the velocity, acceleration, momentum, and force four-vectors with the normalized four-velocity of the fiducial observer are derived, and as a consequence, a three-dimensional analogue of Newton's second law for this particle in curved spacetime is achieved. Since the result is derived in a comoving orthonormal tetrad of the fiducial observer and is presented in the form of geometric objects, it is an elegant generalization of the classical one in flat spacetime.
As a comprehensive application of the GA techniques constructed before, the last task of this paper is to provide an efficient treatment of gyroscopic precession in the STAs of signatures (±, ∓, ∓, ∓) . For a gyroscope moving in the Lense-Thirring spacetime, its relativistic dynamics is first discussed, and some significant results like the three-dimensional generalized equation of motion for the gyroscope are given. Then, by applying the rotor techniques, the geometric object ds ′ /dτ = ds (i) /dτ γ i is able to be directly dealt with instead of ds (i) /dτ , which greatly simplifies the following derivation. The result suggests that if Eq. (165) holds, the spin of the gyroscope always precesses relative to its comoving frame {γ (α) } with Ω (B) (τ ) as the precessional angular velocity. Within the framework provided by the "common" even subalgebra of the two STAs, signature invariant expressions of the relevant physical quantities involved in Eq. (164) are deduced, which clearly indicates that Eq. (165) holds, and therefore, the gyroscope spin indeed precesses in the frame {γ (α) } . After expanding Ω (B) (τ ) up to 1/c 3 order, the gyroscope spin's angular velocities of the de Sitter precession, the Lense-Thirring precession, and the Thomas precession are all directly read out, and their expressions, in the form of geometric objects, are equivalent to their conventional ones in component form, respectively.
All physical laws should be independent of the choice of signature, which implies that many significant techniques constructed in the STA of signature (+, −, −, −) can also be introduced to the STA of signature (−, +, +, +) , and starting from this motivation, we find that the "common" even subalgebra of the two STAs provides a signature invariant GA framework for spacetime physics. In order to pave the way for the applications of these two STAs and their "common" even subalgebra, we elaborate in detail the rotor techniques on Lorentz transformation and the method of handling physics in curved spacetime within the signature invariant framework, and they are of theoretical significance and of practical worth. As two successful paradigms, the treatment of relativistic dynamics of a massive particle and gyroscopic precession clearly shows that the GA techniques constructed in this paper are efficient and reliable. Being straightforward generalizations, these techniques could also be applied to gyroscopic precession in alternative theories of gravity, such as f(R) gravity 30,37-39 , f (R, G) gravity 40,41 , and f(X, Y, Z) gravity 42 . However, since these topics are usually explored by making use of some complicated mathematical tools (e.g., the symmetric and trace-free formalism in terms of the irreducible Cartesian tensors 30 ), it is crucial to develop new techniques to apply these tools in STA. In fact, by generalizing various GA techniques in STA of signature (+, −, −, −) [2][3][4]17 , the approach in this paper could also be applied to other fields, and it has been verified that some topics in classical mechanics and electrodynamics can be described in such a manner. We expect that the applications of this approach will be extended to a wider range in the future, so that the study of spacetime physics in the language of GA could be greatly promoted.