Gravitoelectromagnetism, Solar System Test and Weak-Field Solutions in $f(T,B)$ Gravity with Observational Constraints

Gravitomagnetism characterize phenomena in the weak field limit within the context of rotating systems. These are mainly manifested in the geodetic and Lense-Thirring effects. The geodetic effect describes the precession of the spin of a gyroscope in orbit about a massive static central object, while the Lense-Thirring effect expresses the analogous effect for the precession of the orbit about a rotating source. In this work, we explore these effects in the framework of Teleparallel Gravity and investigate how these effects may impact recent and future missions. We find that teleparallel theories of gravity may have an important impact on these effects which may constrain potential models within these theories.


I. INTRODUCTION
General relativity (GR) has passed numerous observational tests since its inception just over a century ago, confirming its predictive power. The detection of gravitational waves in 2015 [1] agreed with the strong field predictions of GR, as does its solar system behaviour [2]. However, GR requires a large portion of dark matter to explain the dynamics of galaxies [3,4] and even greater contributions from dark energy to produce current observations of cosmology [5]. Given the lack of a concrete theoretical explanation of these phenomena we are motivated to explore the possibility of modifying gravity within the observational context.
There are a myriad of ways in which to consider modified theories of gravity [6][7][8], and to constrain them [9][10][11]. These range from extensions to the standard gravity of GR to more exotic directions. One interesting framework that has gained attention in recent years is that of Teleparallel Gravity (TG). TG is formed by first considering a connection that is not curvature-full, i.e. we consider a connection that is distinct from the regular Levi-Civita connection (which forms the Christoffel symbols). In this way, the gravitational contributions to the metric tensor become a source of torsion rather than curvature. This is achieved by replacing the Levi-Civita connection with its Weitzenböck analog. The Weitzenböck connection is torsion-full while being curvature-less and satisfying the metricity condition [12]. Thus, we can construct theories of gravity which express gravitation through torsion rather than curvature. One such theory is the teleparallel equivalent of general relativity (TEGR) which produces the same dynamical equations as GR while being sourced by a different gravitational action, i.e. one that is based on torsion rather than curvature.
TEGR and GR differ in their Lagrangians by a bound- * gabriel.farrugia.11@um.edu.mt † jackson.said@um.edu.mt ‡ andrew.finch.12@um.edu.mt ary term that plays an important role in the extensions of these theories [13][14][15]. The boundary term naturally appears in GR due to the appearance of second-order derivatives in the Lagrangian [16,17], which is the core difference between GR and TEGR at the level of the Action.
In fact, this boundary term is the source of the fourthorder contributions to f (R) theories of gravity. For this reason, TG features a weakened Lovelock theorem [18][19][20], which as a direct result means that many more theories of gravity can be constructed that are generally second-order in their field equation derivatives. It is for this reason that TG is very interesting because it organically avoids Gauss-Ostrogradski ghosts in so many contexts. The TEGR Lagrangian can be immediately generalize to produce f (T ) gravity [21][22][23][24][25] in the same way that the Einstein-Hilbert action leads to f (R) gravity.
A number of f (T ) gravity models have shown promising results in the solar system regime [26][27][28][29], as well as in the galactic [30] and cosmological regimes [13,31,32]. Of particular interest is its effect on weak lensing in galaxygalaxy surveys [33]. However, to fully incorporate f (R) gravity, we must consider f (T, B) gravity where B represents a boundary term that appears as the difference between the Ricci scalar and the torsion scalar (and will be discussed in more detail in §.II).
Gravitomagnetic tests offer offer an ideal vehicle to probe the rotational behaviour of theories of gravity in their weak field limits [34][35][36][37][38][39]. In fact, gravitomagnetic effects are the result of mass currents appearing in the weak field limit of GR where the Einstein field equations take on a form reminiscent of Maxwell's equations [40] (and do not involve actual electromagnetic effects). These effects emerge as a result of a rotating source or observer in a system, which both give independent contributions to the overall observational effect. For the case where an orbiting observer is moving about a stationary source, Geodetic effects emerge [40] where a vector will exhibit precession due to the background spacetime being curved. This is the general relativistic analog of the wellknown Thomas precession exhibited in special relativity [41]. Another closely related relativistic precession phenomenon is that of the Lense-Thirring effect (or framedragging effect) [42] where the neighbourhood of a large rotating source causes precession in nearby gyroscopes. While independent, these effects are often observed as a combined observable phenomenon such as in the Earth-Moon system about the Sun [43], where the precession of the Moon's perigee is caused by this phenomenon [44][45][46].
Motivated by the Gravity Probe B experiment [47], there have been a number of investigations into the behaviour and predictions of modified theories of gravity [48][49][50][51]. However, the accuracy of this experiment is not enough to adequately differentiate between competing models of gravity. There have also been other experimental efforts such as LAGEOS [52,53] which aimed to perform lasers test while in orbit about the Earth. The MGS spacecraft [54,55] tested gravitomagnetism effects about Mars, while there have also been tests about the Sun [56]. For this reason, there have been a number of ambitious proposals put forward in recent years to further test this relativistic effect and to increase the experimental precision of the observations [57][58][59][60].
In this work, we explore the gravitomagentic effects of TG through f (T, B) gravity, as well as the classical solar system tests within this context. We do this by first expanding into the weak field limit of the theory and explore both the Geodetic and Lense-Thirring effects separately. We then compare their combined results against the recent observations. The manuscript is divided as follows, in §.II we briefly review and introduce TG and its f (T, B) gravity extension. In §.III, we explore the weak field regime of f (T, B) gravity and discuss some important properties of the theory in this limit. Perturbations about a statis spherically symmetruc metric are considered in §.IV. The core results associated with gravitomagentism and the classical solar system tests are then determined in §.V, while a comparison with observational values is presented in §.VI. Finally we conclude in §.VII with some remarks and a discussion. Throughout the manuscript, the speed of light is not set to unity for comparison purposes in the electrodynamics analysis in §.III.

II. TELEPARALLEL GRAVITY AND ITS EXTENSION TO f (T, B) GRAVITY
Teleparallel Gravity represents a paradigm shift in the way that gravity is expressed where curvature is replaced by torsion through an exchange of the Levi-Civita connection,Γ σ µν , with its Weitzenböck analog, Γ σ µν , (we use over-dots to represent quantities determined using the Levi-Civita connection) [17]. GR expresses curvature through the Levi-Civita which is torsion-less, while the Weitzenböck connection is curvature-less and also satisfies the metricity condition [61]. In theories based on the Levi-Civita connection, curvature is given a meaningful measure by means of the Riemann tensor on Riemannian manifolds [40]. This formulation of gravity is retained in most popular modified theories of gravity where gravitation continues to be expressed in terms of curvature of a background geometry. However in TG, irrespective of the form of the metric tensor, the Riemann tensor must vanish since the Weitzenböck connection is curvature-less [62]. It is for this reason that TG necessitates a fundamental reformulation of gravitation in order to construct realistic models of gravity.
GR and its variants utilize the metric, g µν , as their fundamental dynamical object, but TG treats this as a derived quantity which emerges from the tetrad, e a µ . The tetrad acts as a soldering agent between the general manifold (Greek indices) and its tangent space (Latin indices) [63]. Through this action, the tetrad (and its inverses e µ a ) can be used to transform indices between these manifolds Moreover, these tetrads observe orthogonality conditions for internal consistency. The Weitzenböck connection is then defined using the tetrad as [12][13][14]62] where ω a bµ is the inertial spin connection. The Weitzenböck connection is the most general linear affine connection that is both curvature-less and satisfies the metricity condition [63]. The appearance of the spin connection is there to retain the covariance of the resulting field equations [64]. This is an issue in TG due to the freedom in the choice of the components of the tetrads, that is, there is an infinite number of tetrads that produce the same metric tensor in Eq. (2). These different tetrads are related by local Lorentz transformations (LLTs). As a result, the spin connection components take on values to account for the LLT invariance of the underlying theory. Thus, there is a particular choice of frames in which the spin connection components are allowed to be zero [13].
In GR, this issue is hidden in the internal structure of the theory [40]. Considering the full breadth of LLTs (boosts and rotations), Λ a b , the spin connection can be represented as ω a bµ = Λ a c ∂ µ Λ c b [14]. Thus, it is the combination of tetrad and an associated spin connection that forms the covariance of TG.
Given a Riemann tensor that measures curvature, we must define a so-called torsion tensor that gives a meaningful measure of torsion, defined as [13] T σ µν : where the square brackets represent the anti-symmetry operator. The torsion tensor represents the field strength of TG, and transforms covariantly under both diffeomorphisms and LLTs [63]. To formulate a gravitational action, we must define two other quantities. Firstly, consider the contorsion tensor which effectively is the difference between the Levi-Civita and Weitzenböck connections, defined as [17,65] which plays a crucial role in relating TG results with Levi-Civita connection based theories. Secondly, we also need the superpotential which is defined as [63,66] S µν a := K µν a − e ν a T αµ α + e µ a T αν α .
This has been shown to potentially relate TG to a gauge current representation of the energy-momentum tensor for gravitation [67,68]. Then, by contracting the torsion and superpotential tensors, the torsion scalar can be defined as which is entirely determine by the Weitzenböck connection along the same vain as the Ricci scalar being determined completely by the Levi-Civita connection. Naturally, the Ricci scalar calculated with the Weitzenböck connection will vanish since it is a measure of curvature. This property in conjunction with the use of the contorsion tensors allows for a relation between the regular Ricci scalar and the torsion scalar defined in Eq. (7) through [13,14,63] where R is the Ricci scalar calculated using the Weitzenböck connection,R is the standard gravity Ricci scalar determined using the regular Levi-Civita connection, and e = det e a µ = √ −g is the determinant of the tetrad. Thus, the standard Ricci and torsion scalars turn out to be equivalent up to a total divergence term where B = 2∇ µ (T σ µ σ ) is a total divergence term. This relation guarantees that the ensuing equations of motion will be equivalent. Thus, the TEGR action can be written as [13,63] where κ 2 = 8πG/c 4 , and L m is the matter Lagrangian. This action leads to the equivalent dynamical equations as the Einstein-Hilbert action, but the difference in their Lagrangians means that the fourth-order boundary terms are not necessary to form a covariant theory within the TG context. While this does not effect the TEGR limit, it will influence the possible theories that can be formed in the modified gravity scenario.
Considering the same reasoning that led to f (R) gravity [6,7], the Lagrangian of TEGR can be immediately generalized to f (T ) gravity [21][22][23][24][25]. The f (T ) gravity setting produces generally second-order field equations in terms of derivatives of the tetrads [13]. This feature is only possible due to a weakening of Lovelock's theorem in the TG setting [18][19][20]. This fact alone guarantees that f (T ) gravity will not exhibit Gauss-Ostrogradsky ghosts since it remains second-order. f (T ) gravity also shares other properties with TEGR such as its GW polarization signature [69,70].
However, to fully encompass the breadth of f (R) gravity, we must consider the generalization to f (T, B) gravity which contains as a subset the limit f (R) = f (−T + B). Thus, f (T, B) gravity is further generalization of f (R) in which the second-and fourth-order contributions to the theory are decoupled [71].
In this work, we investigate the gravitomagnetic effects of f (T, B) gravity and its effect on observational constraints of the theory for particular models of this setting [69,[71][72][73][74][75][76]. To do this, we need the field equations of the theory, which are determined by a variation of the f (T, B) gravitational Lagrangian density, ef (T, B) to give [13,14,77] e λ where subscripts denote derivatives, and Θ ν ρ is the regular energy-momentum tensor. The spin connection is taken to be zero [69,[71][72][73] since this will be a demand in the work that follows. We will revisit this statement at various stages of the analysis to confirm the consistency of the work. Using the contorsion tensor relations, the f (T, B) gravity field can also be represented as whereG µν is the regular Einstein tensor calculated with the Levi-Civita connection. In this setting, the spin connection depends on the choice of tetrad components and so does not produce independent field equations. However, works exist in the literature that consider this scenario such as Refs. [14,78] where a Palatini approach is considered so that a second set of field equations are produced for the spin connection.

A. The Field Equations
Linearised gravity offers a relatively simple procedure to examine the weak-field metric for a given source. As the gravitational field is assumed to be weak, the metric can be expressed as a Minkowski background plus a small (first order) correction, h µν . In other words, the metric tensor can be expanded as with |h µν | ≪ 1. By extension, a similar consideration can be applied for the linearised expansion for the tetrad: a background value γ (0)a µ which yields the Minkowski metric plus some small correction γ with γ Following the methodology considered in Ref. [69], the resulting perturbed torsional quantities and field equations can be derived. Through the relation between the metric and the tetrad given in Eq.(1), the perturbed quantities are interlinked as Given the equations are constructed in the Weitzenböck gauge (ω abµ = 0), this imposes a constraint on the behaviour of γ which when compared to its LLT form reveals that the background tetrad corresponds to the Lorentz matrices. This is expected as this background tetrad represent a trivial frame, one which constructs the Minkowski metric [63]. As the spin connection is zero here, the background tetrad reduces to a constant, i.e. to the class of constant Lorentz matrices. For simplicity, the background tetrad can be chosen to be γ (0)a µ = δ a µ [79] Under these considerations, the torsion tensor Eq.(4) turns out to be a first order quantity in the perturbations Consequently, as both the contorsion Eq.(5) and superpotential Eq.(6) tensors are linearly dependent on the torsion tensor, then these are also of at least first order.
Ultimately, this implies that the torsion scalar is of at least second order. Observe that this result holds true even if the Weitzenböck gauge is not imposed [69].
On the other hand, the boundary term is first order. This is also consistent with the relationR = −T + B as the Ricci scalar is of at least first order. Indeed, the Ricci tensor and Ricci scalar are given to be where h := h µ µ represents the trace. It is remarked that from here onwards, indices are raised and lowered with respect to the Minkowski (background) metric. Moreover, the d'Alembert operator reduces to = ∂ µ ∂ µ .
The next step would be to extract the perturbed field equations. For simplicity, as both T and B are null at a background level, the gravitational Lagrangian f (T, B) is assumed to be Taylor expandable about these latter values, namely Observe that the coefficient f T (0, 0) = 0 as this corresponds to the effective Newtonian gravitational constant as evident from the field equations Eq.(12) (see for instance Refs. [80,81] for detailed discussions in the case of f (T ) gravity). Under this assumption, the zeroth and first order field equations of f (T, B) gravity Eq. (12) are where the resultR = B (which is valid up to this order) has been used, a property which shall be useful in order to simplify the forthcoming equations. The zeroth order equation confirms the absence of a cosmological constant 2Λ ≡ f (0, 0), maintaining consistency with the linearisation regime as the background geometry is Minkowski spacetime.
As mentioned previously, f (R) gravity is a sub-case of f (T, B) gravity. In fact, the resulting perturbed equations Eq.(23) are practically identical in form to those found in f (R) gravity with the only difference being in the form of the coefficients [7,[82][83][84][85][86][87][88][89][90]. Motivated by this, the same procedure as presented in Refs. [70] shall be followed.
First, the quantityh µν defined as is introduced, withh :=h µ µ . As shown in Refs. [88,91], the Lorenz gauge ∂ µh µν = 0 can be imposed. In this way, the field equations Eq.(23) take a relatively simple form The next step is to obtain the form of the perturbed Ricci scalar. Taking the trace of Eq.(23) yields the relation which is of the same form as the Klein-Gordon equation having an effective mass Depending on the form of the source (and hence of the stres-energy tensor), Eqs. (25) and (26) allow for a full determination of the weak-field metric Eq. (24). Observe that in vacuum, these equations give rise to gravitational waves which polarisation states have already been investigated in detail [69].

B. Solving the Field Equations
In general, the solutions for h µν andR can be obtained by making use of a Green's function G(x, x ′ ), which result intō where r = |x − x ′ | and the Green's function G R defined as [85,92] Within the practical application of the weak-field approximation, it is sufficient to consider a slowly rotating source while keeping all terms up to the order of c −3 . Thus, the stress-energy tensorial components would be negligible within this context. In other words, the stressenergy tensor takes the form [59,93] where ρ is the density of the source and v i is the velocity vector. Alternatively, the off-diagonal components can be simply expressed in terms of the mass current vector j i := ρv i . In this way, we therefore find that where Φ and A are the scalar and vector potentials respectively. This yields the weak-field metric where dΣ 2 = dx 2 + dy 2 + dz 2 .
As the main aim of this work is to match with Gravity Probe B and Solar System observations, it is imperative to treat the source as a slowly rotating spherically symmetric static source having a constant mass M , radius R S and angular momentum J with a constant density profile ρ expressed, for simplicity, as Under these assumptions, for distances sufficiently far away from the source (as the field is weak), the integrals can be solved through the Legendre polynomial expansion where L = |x|,L = |x ′ | and Θ is the angle between the two position vectors x and x ′ . This yields the solutions where J represents the angular momentum vector. Therefore, the weak-field metric takes the simple form where we have defined the function

C. Analogy with GEM
From the resulting weak-field metric, we can make a direct analogy with gravitoelectromagnetism (GEM) to generate the corresponding gravitoelectric and gravitomagnetic fields. Whilst these fields remain effectively unchanged in form, the Lorentz force is affected by the scalarR mode similar to what is encountered in f (R) gravity. Following the steps dictated in Ref. [92], the GEM equations and the Lorentz force equation are obtained as follows. Starting from the Lorenz gauge condition ∂ µh µν = 0, we obtain that with the remaining equations ∂ µh µi = 0 are of order O c −4 and therefore neglected. The gravitomagnetic field B and gravitoelectric field E are then defined as It can then be easily shown that using Eqs. (43) and (25), the GEM equations result: On the other hand, the Lorentz force for a test particle of mass m can be obtained starting from its Lagrangian L = −mc ds dt , using the weak-field metric solution Eq.(36) and expanding up to first-order of the potentials. This yields where γ is the Lorentz factor and v = dx dt is the velocity vector. From the equations of motion d dt ∂L ∂v = ∂L ∂x , assuming that the vector potential A is stationary, it can be shown that up to first order in v 2 /c 2 , the force F ≡ dp dt where p = mγv is the relativistic momentum vector obeys Similar to Ref. [92], one obtains the first two terms which are found in GR (except for a gravitational constant rescaling from f T (0, 0)) with a new contribution arising from the scalar mode. However, if the scalar mode is absent (i.e. µ 2 → ∞), the Lorentz force reduces to its GR form.

D. Comparison with a Spherically Symmetric Metric: The Schwarzchild Solution
In the absence of rotation, the resulting weak-field metric Eq.(41) cannot be directly correlated with the Schwarzchild solution due to the preferred choice of coordinates set by the Lorenz gauge. However, the metric can be transformed into a spherically symmetric form which can then be associated to such known solutions, and shall be notably important when discussing the geodetic effect. Here, we follow the procedure shown in Ref. [85]. The aim is to express the weak-field metric into the spherically symmetric form with A(r) and B(r) representing some scalar functions and dΩ 2 represents the polar symmetry. The necessary coordinate transformation is dictated by the conditioñ where the last equality only holds for a weak-field. In particular, for a spherically symmetric static source, we haver In this way, we obtain that up to first order in M/r Observe that the exponential, similar to f (R) gravity, retains the r dependence. On the other hand, B(r) is found to be Evidently, when µ → ∞ (i.e. in the limit of GR or when f (T, B) → f (T )), the metric reduces to its Schwarzchild form.

IV. PERTURBATIONS ON A STATIC SPHERICALLY SYMMETRIC METRIC: f (T ) GRAVITY
In the previous section, we have initially assumed that the gravitational field is weak, for which the relevant weak-field metric for an arbitrary f (T, B) function was obtained. In what follows, a different approach is considered, particularly in the context of f (T ) gravity. Originally considered in Ref. [94] which was further pursued in Ref. [95], the idea is to assume a static spherically symmetric geometry arising due to a spherically symmetric static source of mass M . Then, one solves the field equations Eq.(11) to obtain the corresponding metric. Since no exact solutions have been obtained following this approach (although exact solutions can be found assuming, for instance, that the Lagrangian exhibits a Noether symmetry [96,97]), a perturbative approach is employed where it is assumed that the f (T ) Lagrangian takes the form of where ǫ ≪ 1 represents a small, fiducial, order parameter, which will be omitted once the perturbations are solved. The role of the latter is to represent the small correction to the TEGR Lagrangian. In this way, the scalar functions A(r) and B(r) of the metric Eq.(49) are expected to be in the form of a background solution plus a small correction. This will allow for corrections which were previously omitted in the weak-field regime. For simplicity, the source shall be assumed to be non-rotating as no perturbed solutions to the Kerr metric have been yet obtained.
In this formulation, the TEGR term gives rise to the exact Schwarzchild solution while the small correction sourced by F (T ) yields the first-order correction to the solution. Since the main interest lies in the Gravity Probe B results, an alternative approach is to assume the gravitational field to be weak, meaning the metric can be approximated by a Minkowski spacetime background plus a small correction. Both approaches shall be presented and show that the same results are ultimately recovered, whilst offering a more detailed analysis on the affect of the f (T ) Lagrangian on the geodetic effect.

A. Perturbations on the Schwarzchild Solution
The Schwarzchild correction can be obtained by taking the scalar potentials to be expressed as for some functions A and B. To solve for the corrections, the field equations are perturbed up to first order in ǫ.
For simplicity, the power-law ansatz F (T ) = αT p for some constant α and p is considered. Furthermore, unless otherwise stated, GM c 2 → M .
The solutions for the scalar functions can be obtained from the differential equations Although a general solution is not recovered in general, some special cases are considered. For p = 2, the scalar functions are given to be [94,95] The solutions are well behaved in the sense that in the absence of a source, the solutions reduce to Minkowski space as expected, i.e. when M → 0, A(r), B(r) → 1. If the gravitational field is weak, the solutions follow the order expansion which agrees with the weak-field metric Eq. (41) in the limit of f (T, B) → f (T ) = T + αT p up to first order in M/r. Observe that the α contributions do not appear in the latter as it is a higher-order contribution. Solutions for other values of p are considered in Ref. [95]. For the purpose of the analysis which follows, the solution for p = 3 is given and is listed in Appendix A.

B. An Alternative Approach for a Weak-Field Limit
If the field is assumed to be weak, the scalar functions can be expanded around a Minkowski background according to A(r) = 1 + ǫA 1 (r) + ǫ 2 A 2 (r) + ǫ 3 A 3 (r) + . . . , (63) B(r) = 1 + ǫB 1 (r) + ǫ 2 B 2 (r) + ǫ 3 B 3 (r) + . . . . (64) Once again, assuming that the Lagrangian f (T ) is Taylor expandable about T = 0, and solving the field equations order by order yields where c 1,. . . ,6 are integration constants. To determine these constants, we impose the following constraints. As r → ∞ (i.e. far away from the source), the metric must reduce to Minkowski spacetime and thus sets c 2,4,6 = 0. On the other hand, according to the solution obtained in §.III, namely Eq.(41), we find that c 1 = 2M fT (0) (alternatively, it can be reasoned that in the limit of TEGR, the metric must reduce to the Schwarzchild metric). Finally, the constants c 3,5 have to be zero otherwise the solution does not reduce to its TEGR limit for f (T ) = T . Therefore, the final solution is Taking f (T ) = T + αT 2 recovers the previously obtained weak-field limit solution as expected.
Observe that for f (T ) = T + αT n , n > 2 (n integer) does not reveal any contributions at this order meaning their effects are even smaller. On the other hand, this approach is not applicable for functions which are not expandable about T = 0, for instance f (T ) = T + αT n , n < 0 and even some cosmologically viable ones such as the Linder model f (T ) = T +αT 0 1 − e −p √ T /T0 for some constant p. However, there exist cosmological model Lagrangians which may be further investigated for such weak-field observational tests, such as f (T ) = T + αT 0 (1 − e −pT /T0 ) and f (T ) = T + αT n tanh(T /T 0 ) for appropriate values of p and n [31].
Observe that the result is in agreement with the parameterised post-Newtonian (PPN) approximation investigated in Ref. [98] since there is no deviation up to second order expansion. The first modification appears at third order when the f T T term contributes to the behaviour.

V. OBSERVATIONAL CONSTRAINTS
As we have now closely discussed the theoretical foundations to obtain the necessary metrics, in what follows, we apply those results to observations obtained by Gravity Probe B and from classical Solar System test observations. In particular, we shall focus on the geodetic effect (de Sitter precession), the Lense-Thirring effect, Shapiro time delay, light bending and perihelion precession.

A. Geodetic Effect
The geodetic effect describes the effect of a precessing gyroscope due to its orbit around a massive central body.
Here, we obtain the precession rate following Rindler's approach [99]. Starting from a spherically symmetric metric, we consider the system to be rotating at an angular frequency ω, i.e.
By assuming the gyroscope to lie in a circular polar orbit (at an angle θ = π 2 ) allows us to rewrite the metric in a canonical form where e 2Ψ ≡ A − r 2 ω 2 . As shown in Ref. [99], the angular frequency of the gyroscope is given to be where k ij is the spatial 3-metric and ω i ≡ e −2Ψ ωr 2 δ 3 i , which simplifies to The angle after one full revolution is then given to be α ′ = Ω∆τ , where ∆τ represents the proper time of the gyroscope, which can be obtained directly from the metric Thus, the precession over one orbit is α = 2π − α ′ , which implies that the precession rate per year is given to be

B. Lense-Thirring Precession
It is well known that the Lense-Thirring precession in GR can be simply derived by assuming a freely falling gyroscope initially at rest with an angular spin vector S µ . Taking u µ to represent the gyroscope's rest frame velocity, we have that S µ u µ = 0. Then, the Lense-Thirring precession rate would be then obtained using the geodesic equations In the context of teleparallel gravity, the gyroscope moves according to force-like equations Despite this apparent difference, the above is mathematically equivalent to the geodesic equation due to the fact that K σ µν =Γ σ µν − Γ σ µν . Nonetheless, the force-like equations offer a different interpretation as discussed, for instance, in Refs. [100,101], as the teleparallel force equations allow for a separation between gravitation and inertia which has important implications on the weak equivalence principle (WEP), which lies beyond the scope of this manuscript (see, for instance, Ref. [67] for further discussions on the topic). Within the assumption that the WEP holds, one can follow the same steps encountered in GR. Alternatively, one can work out directly using the torsion and contorsion tensors to obtain the same result.
If the field is weak, the field equations reduce to where Ω k ≡ − 1 2 ǫ kmn ∂ m h 0n defines the angular velocity precession vector of the gyroscope. Following the results obtained in the f (T, B) weak-field solution Eq.(41), we find that the Lense-Thirring precession rate Ω LT remains unaffected except for a Newtonian rescaling, which is expected as the electromagnetic field is identical to that found in GR. However, this result is only valid within the context of weak-fields and thus remains to be investigated in the case of strong gravitational fields.

C. Shapiro Time Delay
The effect of Shapiro time delay [102] can be derived following the steps listed in Ref. [103]. Here, we focus on deriving the α dependent correction for the f (T ) powerlaw model. For the given spherically symmetric metric Eq.(49), the time delay of a radio signal as it travels from the Earth to Mercury and back, as the signal passes through the closest point of approach R ≃ R ⊙ to the Sun is where r ⊕ and represent the Earth and Mercury orbital radii respectively, and t(r, R) is defined as Using the fact that generally, the orbital radii satisfy the conditionr ≫ R, together with the weak-field metric solutions Eqs. (59), (60),(A3),(A3), we find that the α contribution takes the following forms

D. Light Bending
The total deflection angle of light is derived following the method used in Ref. [26]. A photon is assumed to be emitted from some far away source at an angle φ = −π/2. It travels and reaches a point of closest approach r =r ⋆ with respect to some spherical massive source at φ = 0, and then continues to travel away from the source approaching an angle of φ = π/2. To account for the deflection angle due to the gravitational attraction of the source, we start from the spherically symmetric metric Eq.(49) within the equatorial plane θ = π/2 to obtain that the path of the photon obeys the second order differential equation where u = 1/r and u R = 1/R is the inverse impact parameter, with boundary conditions u(0) = u ⋆ = 1/r ⋆ and u(±π/2) = 0. Since the differential equation cannot be solved in general, even for weak-field sources, the perturbative iterative method considered in Ref. [104] is applied.
The approach aims to obtain a perturbative solution by taking the mass M as the perturbation parameter, i.e. we let u = u 0 + u 1 + u 2 + . . . (82) where u i represents the solution up to O M i . As an illustrative example, we derive the perturbative solution for the power-law model with p = 2. In this case, the differential equation up to O(M 3 ) is where we have denoted primes to represent derivatives with respect to φ. This yields the following ordered system of differential equations Once the solution for u is obtained, following Rindler and Ishak's approach [105], it can be shown that the total deflection angle is given to be which, for the quadratic f (T ) Lagrangian yields the following solution A similar analysis for the cubic f (T ) Lagrangian reveals that the total deflection angle is . (87) Observe that in both cases, the GR second order mass correction found in Refs. [104,106,107] is recovered. In general, for a Taylor expandable f (T ) model within the regime of weak gravitational fields, the first deviation from GR appears at O(M 3 ), having the form Naturally, the quadratic weak-field result is recovered while for the cubic case requires the higher order contributions.

E. Perihelion Precession
The effect of α for the power-law ansatz Lagrangian on perihelion precession has been investigated in great detail in Refs. [94,95]. Here, we shall only quote the results: [108] p = 2 ∆φ = 16παM 2 r c 4 (89) where r c represents the circular radius of the orbit. For the n = 2 case, the detailed analysis in Ref. [94] leads to a bound of α 10 20 km 2 .

VI. NUMERICAL RESULTS
In this section, we make use of the weak-field solutions listed in §.III and IV against observations in order to constrain the Lagrangian free model parameters depending on the model considered. It is important, however, to comment on the results for an arbitrary f (T, B) model for the case when µ → ∞.
Although the weak-field metric has been obtained in its spherically symmetric form, the scalar functions A(r) and B(r) are not truly expressed in terms ofr since the relation between r andr is not invertible. This leaves two unknown parameters, the isotropic radial coordinate r and µ. However, r is not measured and hence one must instead opt to impose specific values of µ to determine whether the results would then be consistent. Since the goal is to constrain the Lagrangian parameters (and hence constrain µ through observations), this option is not investigated in detail. Nonetheless, if µ is sufficiently large, the contributions would be small enough that deviations from observations (and hence form GR) are expected to be effectively negligible.
On the other hand, a more thorough investigation can be inferred in the case of f (T ) gravity using the results obtained in §.IV. In particular, we shall make use of the results for the two power-law ansatz values considered, namely p = 2 and 3, which eventually lead to observation constraints on the constant α.

A. Geodetic Effect
In April 2004, Gravity Probe B was launched starting its year and a half flight mission, with the purpose of accurately measuring the geodetic and the frame dragging precession rates while in orbit about the earth. A geodetic precession rate of −6601.8 ± 18.3 mas/yr was measured while in a polar orbit at around 642 km [47].
Through the use of Eq.(74), the α constraints are obtained as listed in Table I. The table also illustrates the α constraint which has to be obeyed for the weak-field approximation to hold (which is a direct consequence of the assumption that the perturbation F (T ) ≪ T ).
Based on the expressions listed in Table I, the corresponding numerical constraints are then obtained as shown in Table II. Evidently, the constraints obtained from observations are well within the expected bounds of the weak-field condition which supports the consistency of the weak-field approach.

B. Classical Solar System Constraints
For Shapiro time delay and light deflection, the PPN formulation together with observations from the Cassini spacecraft pose a viable opportunity to obtain constraints. As illustrated, for instance in Ref. [103,109], the γ PPN parameter appears in the former tests as follows.

VII. CONCLUSION
The main result of this work is that both the classical solar systems and the gravitomagnetic constrains from Gravity Probe B result in a constant on the coupling parameter to |α| 10 16 km 2 for p = 2 and |α| 10 43 km 3 for p = 3. This forms one of the strongest constraints on this parameter (to the best of our knowledge).
Gravitomagnetic effects are imperative for understanding the weak field limit of modified gravity in the context of rotation. In this work, we have explored these related   effects in the TG framework. TG offers a novel possibility of constructing gravitational theories in which the background manifold is torsionful rather than curvatureful. While this is dynamically equivalent to GR in the TEGR limit, modifications of the TEGR action produce theories which may be distinct from what can be constructed in regular curvature-based theories of gravity. This allows for the possibility of totally new models of gravity that may have important consequences for meeting the observational challenges of the coming years. The main crux of the weak field analysis stems from the analysis in §.III where we take an order by order expansion of a general f (T, B) gravity Lagrangian. In Eq.(26), this is found to potentially behave as a massive theory with a mass that is mainly dependent on whether a boundary term contribution is present or not. This approximation is then set into the field equations with a slowly rotating source to find metric solutions in Eqs. (32)(33)(34)(35). In §.III C we go into the details of how this analogy tallies with the well-known GEM effects to produce a Lorentz force-like effect in Eq. (48). Finally, we compare this with the Schwarzschild solution to determine the relation to the effective mass of the general f (T, B) model.
Limiting ourselves to f (T ) gravity, we explore the possibility of perturbative solutions in §.IV where exact solutions are found up to perturbative order in the spherically symmetric setting. These were also investigated in the literature [94,95] and remain an interesting avenue of research in the TG context. In this part of the work, we investigate two possible routes to the perturbative analysis which both agree in their PPN limit.
The traditional gravitomagnetic effects of the geodetic and Lense-Thirring phenomena are determined in §.V.
The Geodetic effect naturally emerges for a static system with a rotating observer. This is achieved by a coordinate transformation, as prescribed in Eq. (69). This eventually produces Eq.(74) which is our result for the geodetic precession rate and the main result of that subsection. The Lense-Thirring effect is then determined for this TG case where the main result result is shown in Eq. (77) which is comparable to the Gravity Probe B mission result. In fact, in §.VI we use the results of this mission to constrain our parameters for the various potential models under investigation.
Gravitomagnetic effects have the potential to have an important impact on understanding which modified theories of gravity are viable and may play an important role in the coming years for developing realistic modified theories of gravity.