Galactic Rotation Dynamics in f(T) gravity

We investigate galactic rotation curves in $f(T)$ gravity, where $T$ represents a torsional quantity. Our study centers on the particular Lagrangian $f(T)=T+\alpha{T^n}$, where $|n|\neq 1$ and $\alpha$ is a small unknown constant. To do this we treat galactic rotation curves as being composed from two distinct features of galaxies, namely the disk and the bulge. This process is carried out for several values of the index $n$. The resulting curve is then compared with Milky Way profile data to constrain the value of the index $n$ while fitting for the parameter $\alpha$. These values are then further tested on three other galaxies with different morphologies. On the galactic scale we find that $f(T)$ gravity departs from standard Newtonian theory in an important way. For a small range of values of $n$ we find good agreement with data without the need for exotic matter components to be introduced.


I. INTRODUCTION
Over the past several decades, the consistent missing mass or dark matter problem has attracted increasing interest in modified and alternative theories of gravity. This is one of the biggest potential contentions between general relativity (GR) and observational astronomy [1]. In fact, in essentially every observed galaxy, it appears that the expected rotational velocities of Newtonian gravity do not conform with the measured values. Since the discrepancy was first noted, observational techniques have drastically improved, making it possible to study and constrain the motions of luminous matter with much greater precision [2] then ever before. The ΛCDM model solves this problem by assuming a much larger dark form of matter present in every galaxy and cluster of galaxies. However, over the decades no observation has confirmed the existence of this material either way [3][4][5].
The unmodified Newtonian picture treats each galaxy as a collection of N individual sources with a typical mass M , and then combines the respective Newtonian potentials. This generates a global outline of the Newtonian potential. Assuming a thin disk shaped galaxy with an exponential radial light distribution Σ(R) = Σ 0 e −R/β d , where β d is the disk scale radius and Σ 0 characterizes the central surface brightness, the resulting circular velocity for a test particle at a radial distance R from the center of the disk * andrew.finch.12@um.edu.mt † jackson.said@um.edu.mt is given by the Freeman formula [6,7] v e T EGRd where G is Newton's constant of gravity, and I n and K n (n = 0, 1) are modified Bessel functions of the first and second kind respectively [8].
The central surface brightness distribution can not be expressed by the disk on its own. As such another component of the galaxy, the bulge component, is defined. These bulges have a higher stellar concentration than the stellar disk. The bulge mass distribution is approximately spherically symmetric, with stellar orbits being roughly circular. Thus the individual sources would follow the velocity relation [10] where the mass, M b (R), is determined through the surface density calculation [11] ρ(R) = 1 π and Σ b (R) is the de Vaucouleurs profile for the surface mass density suit and treat these two observable regions separately.
Teleparallel gravity is one alternative to GR where the mechanism by which gravitation is communicated is torsional rather than curvature based. This theory was first proposed by Einstein himself in order to unify gravitation and electromagnetism [12]. The equivalent reformulation remained dormant until the relatively recent resurgence in alternative and modified theories of gravity. Now, the resulting theory is equivalent at the level of equations and phenomenology. However, the resulting action can be generalized much like the f (R) proposition. This is the origin of the distinction between the two theories in that the equivalence no longer holds. In fact, while the f (R) field equations are fourth order, the f (T ) field equations are second order.
Galactic rotation curves were first tackled in f (T ) gravity in Ref. [13]. While this study provides promising results there are some issues with the procedure and theory. Firstly, in the work, the separate regions of galaxies are not treated individually. Secondly, the mass profile is considered to be spherically symmetric. This may be true for the surrounding dark matter halo but not for the luminous matter segment of the galaxy. Finally, the choice of tetrad, e a µ , for the metric should be treated with somewhat different field equations. We expand on this in section II but the main point is that certain choices of tetrad require an associated quantity (spin connection) to be defined so that the theory continues to respect local Lorentz invariance, as explained in Ref. [14]. However, the work does result in flat rotation disks in general.
The paper is divided as follows. In section II we give a brief overview of f (T ) gravity with an emphasis on the relevant solution to the field equations. In section III the mechanics of galactic rotation curves are worked through for the current setting. The resulting dynamics are then applied to Milky Way data for several instances of the general model under consideration in section IV. The best model is highlighted, and then a three-parameter fitting for the bulge and disk masses as well as the new coupling parameter is described and determined. In section V, we tackle the broader problem of galaxies with other morphologies. Finally, the results are discussed in section VI.

II. f (T ) GRAVITY
GR and its modifications are largely based on the metric tensor, g µν , which solves the field equations and acts as the fundamental dynamical variable. The metric acts as a potential quantity with curvature being represented through the Levi-Civita connection (torsion-free), Γ λ µν . In teleparallelism this connection is replaced by the Weitzenböck connection,Γ λ µν . This new connection is curvature-free and is based on two fundamental dynamical variables, namely the tetrads (or vierbein) and the spin connection. The tetrads, e a µ , are four orthonormal vectors that transform inertial and global frames in that they build the metric up from the Minkowski metric by means of an application of this transformation. They also observe the metricity condition. Physically, they represent the observer and can be related to the metric tensor by means of [15] g µν = η ab e a µ e b ν .
These conditions are not enough to fully constrain the tetrad, and so there is an element of choice in forming the tetrad frames.
The other necessary ingredient to describe f (T ) gravity is the spin connection, ω b aµ [14,16]. This is not a tensor and its particular form depends heavily on the system under consideration, that is, it accounts for the coordinate system such that the theory remains covariant. The choice of tetrad plays a deciding factor in whether the spin connection vanishes or not. This leads to a division in tetrads [17], there are pure tetrads whose associated spin connection vanishes, and impure tetrads who spin connection gives a nonzero contribution. In the current case, we will consider pure tetrads and so will not consider the spin connection any further.
With the introduction of the spin connection, the tetrad remains the fundamental dynamical field on the manifold since every tetrad ansatz produces a well-defined associated spin connection. The question of the inheritability of solutions in GR to teleparallel gravity is then resolved since TEGR is equivalent to GR at the level of equations [14]. As vacuum GR solutions are also solutions to f (R) gravity, TEGR (or GR) solutions are also solutions of f (T ) gravity.
In GR, curvature is communicated between tangent spaces through the Levi-Civita connection.
While the resulting TEGR theory is equivalent to GR at the level of equations, the ingredients leading up to this are not. The difference between the Weitzenböck and the Levi-Civita connections is represented by the contorsion tensor Lastly, the superpotential tensor is introduced This is defined purely for convenience in the resulting equations [15] but plays an important role in the gravitational energy-momentum tensor of teleparallel gravity. These tensors can be contracted to form the torsion scalar, T = T a µν S µν a , which is the Lagrangian for TEGR.
The distinction between GR and TEGR can now be made clearer, i.e. the relationship between the Ricci scalar, R and the torsion scalar, T , can be expressed explicitly. The difference between the two quantities obviously lies in a boundary term since they produce the same theory at the level of equations [15]. The distinction in the difference can be quantified through [19] R(e) = −T + B, where B = 2 e ∂ µ eT λ λµ = 2∇ µ T λ λµ is the boundary term. This relation represents the source of the disparity between the f (R) and f (T ) generalizations since the boundary term no longer remains a total divergence term when the generalization is taken. Thus, GR and TEGR would be indistinguishable in terms of observations. However, while astrophysical observations cannot test distinctions between curvature and torsion, they can compare the predictions of the models available. Once enough tests have been compiled, a multi-test survey would be the best way to compare and contrast the individual models of the two theories.
Thus, taking the Lagrangian −T + B will precisely reproduce the Ricci scalar. As with the generalization of the GR Lagrangian to the f (R) class of theories [20,21], TEGR can also be generalized to f (T ). However, the relation between the resulting theories stops being equivalent since f (R) = f (−T ) + B (f (R) = f (−T + B)). In fact, the ensuing field equations are unique in that, out of the three possible quantities involved, namely R, T , and B, f (T ) is the only Lagrangian that produces second order field equations [19,22]. Therefore the action with arbitrary functional form of the torsion scalar, f (T ), is given by whereκ = 4πG and e = det e a µ . Taking a variable with respect to the tetrad results in the field equations [15] where Θ µ a ≡ 1 e δLm δe a µ , f T and f T T denote the first and second derivatives of f (T ) with respect to T , and L m is the matter Lagrangian.
In Ref. [23], the power-law instance of the Lagrangian is investigated, i.e. f (T ) = T + αT n where α is a coupling constant and |n| = 1 is any other real number. The following weak field solution is found where and For this case the torsion scalar is defined as follows In the limit of vanishing α, GR is again recovered. By determining the galactic rotation curve dynamics for this solution, we will investigate the effect of various n values on the resulting behavior.
In order to determine the rotational curve profile we consider a point particle with energy, E, and angular momentum, L, performing orbits about the galactic core. In the following we investigate this type of orbit with a focus on the effective potential. We then use this result to determine the velocity profile for the disk and the bulge separately.

A. Effective Potential
In order to obtain an effective potential we follow the procedure described in Wald [24]. We start by determining the conservation relations for the energy and angular momentum from the background metric in Eq. (14). These are given by and where ζ µ = (δ/δt) µ and Ψ µ = (δ/δφ) µ are the static and rotational killing vectors respectively. The background metric naturally leads to the radial differential relation below where Eqs. (19)(20) were used.
The effective potential, V e , can be read off from Eq.(21) by comparison with Ref. [24] dr dτ As a gravitational source, the system under consideration is not being treated as having any rotation, and so we can set L = 0. Finally, if we assume roughly circular orbits, i.e. dr dτ = 0, and simplify the resulting expression, then we find an effective potential Here we see that the effective potential includes the GR potential as well as an extra f (T ) contribution which can be divided as follows and Since the velocity contribution of the GR potential is known, we will now continue with the derivation for the f (T ) component.
With the effective potential in hand, it is now possible to obtain the velocity curve profile. To do this, we consider the centripetal and gravitational acceleration equations [25] and Assuming circular paths for the orbiting stars and dust, the effective velocity profile for the central mass turns out to be [6] v where the potential is inherently negative.

B. Disk and Bulge
Given the geometric diversity of the galactic disk and bulge regions, the two sectors are treated separately in the following calculations. In particular, the core contrast in the treatment is related to the difference in their mass density distributions which clearly affects the whole calculation due to the stark change in the effective potential.

Disk
In order to determine the velocity curve profile of the disk component of galaxies we follow the method developed in [26][27][28]. Consider a system of N galactic bodies. The calculation of the velocity profile will necessarily involve the sum of the combined potential of the individual sources within the galaxy. To measure the potential for a particular position with radius R, all other source will be summed together. Consider the n th source with radius R from the galactic center; the distance between the position where the potential is being measured and the n th source will be denoted by r. This is depicted in Fig.(1) where cylindrical coordinates are used.
FIG. 1: For a position with radius R, the relative distances of an n th source are shown, where the radius of the source is denoted by R , and the r represents the distance between the source and the position where the potential is being measured. All distances are in cylindrical coordinates.
With this picture in mind, we can proceed to express the unknown radius r in terms of other radial terms as follows where z is the height of the n th source, and (φ, z) represents an arbitrary reference position which we choose to take as the origin.
In order to obtain the combined potential of the whole disk we integrate over the whole range of source positions available where ρ is the mass density distribution given by [29] ρ where M 0 represents one solar mass, N is the total number of sources in the disk sector, and β d is the scale radius of the galactic disk.
It is at this point that the index n cannot be left arbitrary in value, i.e. we must consider a set of values in order to proceed. The division is as follows: Integer values in the range −∞ < n < 0, n = 0, all values in the ranges 0 < n < 1 and 1 < n < 3 2 . The core of the problem has to do with the expansion of the radial factor in the last integral of Eq. (30). For the instances of integer values in the range −∞ < n < 0 and n = 0, this results in the integral The calculations for the other ranges are presented in the appendices since they are more intricate. The resulting velocity curve profiles for the disk sector are given below for various n values 0 < n < 1 : This range of n results in an integral that does not have a solution. Thus we solve it numerically in the next section.

Bulge
The bulge velocity profile calculation differs from the disk significantly in that it can be completed independently of the value of n. Following Refs. [29,30], we initially treat the bulge as a spherical mass which straightforwardly leads to For this region of the galaxy, the spherical mass distribution can be described through the de Vaucouleurs profile shown in Eq.(4) which directly leads to the modified velocity profile where the mass is calculated through the density distribution in Eq. (3) and which can turned into an integration over the individual spherical shells, giving The contribution from the modified f (T ) gravity terms to the velocity profile then turns out to be where M g(R) is the Meijer G-function defined through [31] M g(R) ≡ For the TEGR case, using the potential in Eq. (24) in the calculation in Eq.(39) leads directly to the velocity profile v e T EGRb 2 = 8388608κG 1 where β b is the bulge scalar radius [30]. Now that both velocity profile sectors are derived, the full velocity profile can be derived for the f (T ) = T + αT n Lagrangian. This results in the velocity profile where both TEGR and modified f (T ) gravity components are added for the disk and bulge segments of the galaxy.

IV. THE MILKY WAY GALAXY
In the first part of this section we compare results for the velocity profile with Milky Way data using the combined velocity equations in Eq.(46) to determine the best range of values of n in the Lagrangian, f (T ) = T + αT n . At this stage, we simply fit for the constant α. In the second part of the section the determined best range of n will be used to fit for the surface mass density of the bulge Σ be , the mass of the disk, M and the coupling constant α to further test the velocity profile against real world data.
A. Determination of best range of index n The data set being utilized is collected from two sources, namely Refs. [30,32]. A representative sample of this data is shown in Table I; this spans the breath of the region under consideration. Here, R is the radial galactic distance shown in Eq. (29), v is the rotational velocity of the sources at that radial distance, and e u and e l are the upper and the lower error bars of the velocity values respectively.
In Table II we present the necessary values of the constants to calculate the rotation curve profile contribution of the Milky Way bulge and disks for the various ranges of n respectively [30].   The velocity profile presented in Eq.(46) can now be used to find suitable values of the coupling parameter, α. We do this by employing a least-squares approach to minimize the difference between the predicted and observed values. This is performed by setting v(R m ) e αb 2 = αu(R m ) e αb 2 and v(R m ) e αd 2 = αu(R m ) e αd 2 and using the relation (47) Here D max represents the size of the data set, i.e. the number of velocity points considered.
In Fig.(2) we present plots for a range of values of n with best fits for α; this ranges from n = −4 to n = 1.4. The corresponding GR plot is also being shown. Each plot consists of data points with associated error bars, the GR prediction for the Milky Way galaxy, and the best fit curve for the f (T ) model being considered where the coupling parameter α is being fitted. On a similar note, every plot has an embedded figure showing the difference between the best fit and the GR result. As both predicted curves plateau so does their associated difference in the embedded figure.
The first row shows negative values of n where the behavior of the best fits are fairly similar in that they quickly diverge for data points in the disk region. Moreover, the predicted curves are completely nonphysical. In the following row, the special case of n = 0 is considered. This corresponds to Einstein's GR with a constant similar to the cosmological constant. However, in this case the constant is designed to account for the anomalous rotation curves of galaxies. Naturally, the constant gives a near linear increase in velocity against the galactic radius. Again, the resulting velocity curve is unrealistic. The linear behavior is best seen in the embedded figure. As with the instances of negative n, the behavior performs worse when compared with GR.
Next we consider the 0 < n < 1 range. The predictions from the f (T ) functional model become much more promising in this region as compared with the preceding values. As with the previous cases, the behavior in the core of the galaxy is relatively well behaved. Although as larger radial values are considered, we find that the predicted curve overshoots the velocity curve data points. The predicted curves get better as the value of n approaches unity. However, this value of n was excluded from the solution presented in Eq.(14) [23]. In the limit of n actually taking on the unity value, the Lagrangian tends to a re-scaling of GR.
Finally the 1 < n < 3 2 range is investigated against the Milky Way velocity curve. These provide the best behavior of graphs from all the ranges considered in this work. As in the other plots, a cut-off value for the maximum radius considered is set to R = 3 × 10 18 . The predicted velocity profile does not have much interest for us beyond that point. The basic result from this range of values of n is that as n gets closer to unity the resulting velocity curve profiles behave better as compared to the observational points in question.
The best behaved value of index of the f (T ) Lagrangian is given by n = 1.00001 which results in the best fit coupling parameter α = −2.428 +1.303 −1.577 × 10 −34 km 2.00002 kg −1 s −2 . This is shown in Fig.(3). The resulting curve agrees with GR in the core part of the galaxy with slight variations in the central region. The difference between the two curves then plateaus, with the f (T ) function behaving better against observation than GR without dark matter contributions.

B. The three parameter fitting
Heaving determined that the best range of values for n is 1 < n < 3 2 , we now move to the next stage of the process which is to produce a three parameter fitting for the f (T ) rotation curve. This method should fit the surface mass density of the bulge Σ be , the Mass of the Disk M and the coupling constant α. In this case, the fitting function proves to be nonlinear and as such an iterative method accompanied by a cross validation fitting is employed to determine these three parameters.
The results of this fit are presented in Table III where the first column consists of the tested values of n, the second column is the minimized sum of the square of the residuals for the fitting function used, S, the third and fourth columns are the fits for the surface mass density of the galactic bulge, Σ be , and the mass of the disk, M , respectively and the final column is the fit for the coupling constant α. The resulting plots from these values can be seen in Fig.(4). From the plots it is possible to narrow down the best fit to two values of n, namely n = 1.00001 and n = 1.000001. The choice between these two instances can then be made using the second column of Table III. The S values for n = 1.00001 and for n = 1.000001 are 2.78 × 10 9 and 1.43 × 10 29 respectively. This clearly indicates that n = 1.00001 provides us with the best f (T ) rotation curve fit. The best fit value for n results in a bulge surface mass density of 5.94 × 10 6 kg km −2 , a disk mass of 4.08 × 10 10 M and a coupling constant of −4.51 × 10 −34 km 2n kg −1 s −2 . The Surface mass density of the bulge Σ be is directly related to the bulge mass through the equation ter. Ref. [36] reports a total combined mass of 5.22 × 10 11 M of the Milky Way which includes the dark matter contribution. In Ref. [35], it is given that the bulge mass is 0.91±0.07×10 10 M and the disk mass is 5.17 ± 1.11 × 10 10 M giving a total mass of 6.08 ± 1.14 × 10 10 M . The total mass of these contributions to the galaxy obtained in this work agrees with the order of the total masses from both Ref. [35] and Ref. [36]. Furthermore, the values of the total mass, and the disk mass obtained here even fall within the error bars of the respective masses given by Ref. [35].
Finally we compare this fit, which for ease of reference is given in Fig.(5), with the best fit while just fitting for α, Fig.(3). With this new three parameter fitting the f (T ) rotation curve successfully interprets the rotational data profile at the inner parts of the galaxy while still retaining the plateau effect obtained in the first fit. Having a value of n that is greater than 1 also prevents the possibility that the velocity curve will tend to infinity at larger radii, which would lead to an unrealistic prediction.

V. WHAT ABOUT OTHER GALAXY MORPHOLOGIES?
In this section we consider three other galaxies apart from the Milky way which adhere to separate galaxy morphologies. Specifically we consider the bright spiral galaxy NGC 3198, the low surface brightness galaxy UGC 128 and the irregular dwarf galaxy DDO 154. These were chosen so as to test the fit on different types of galaxies since these clearly exhibit their own separate morphologies.
The rotation curve data for these galaxies is shown in Table V. The data used for NGC 3198 and UGC 128 was obtained from Ref. [33], and the data for DDO154 was obtained from Refs. [33,34]. In these three cases the surface mass density of the galaxies' bulges, Σ be , is assumed to be zero as they are variants of disk galaxies. The mass of their disks, M , was fit while keeping the coupling constant α fixed at the value obtained by the best fit of the Milky Way in the previous section. The results for the fit can be seen in Table IV. The first column of Table IV represents the name of the galaxy considered, the second column gives the disk scale radii, β d , the third column gives the galactic masses as predicted by GR, M GR and the fourth gives the value of the constant n at which these galaxies were fit. The fifth and sixth columns provide the fitting values for the disk masses and the value of the coupling constant α considered respectively. Finally the last column provides the data sources for β d and M GR . Using these results the three plots shown in Fig.(6) were generated.
The first galaxy considered, NGC 3198, is a bright spiral galaxy [37]. The data available clearly shows that the rotational velocities plateau at around 150 km s −1 . Using the same fitting method used to fit the f (T ) rotation curve for the Milky Way results in a luminous matter disk mass of 2.49 × 10 10 M . The mass obtained is of the expected order and is smaller than the mass predicted by GR. Overall the f (T ) fit is by far more accurate than the General Relativistic one. Up to a galactic radius of around 3×10 17 km the rotation curve fits perfectly with the data, while further out the predicted curve falls below the velocities of the data points. The reason for this drop in the fit curve may be attributed to the fact that no spiral arm equation was included in the fit which could also be the partly the reason why the mass is somewhat less than expected.
The second galaxy considered was UGC 128, a low surface brightness disk galaxy [40]. Here the data points once again produce an evident plateau at around 140 km s −1 . Once again the order of the disk mass is as expected though the magnitude is smaller than that predicted by the GR values. It should be noted that in both this case and the case of NGC 3198 the masses where fit while also fitting for the mass of non-luminous matter. As such, the larger values can, at least in part, be attributed to a missing distribution of mass from non-luminous to luminous. From the plot in Fig.(6.2) we observe that the f (T ) fit is accurate throughout accept for a drop in the data velocity magnitudes between 3×10 17 km and 7 × 10 17 km.
The third and last galaxy considered, DDO 154, is an irregular dwarf galaxy. In this case it proved not possible to fit the f (T ) rotation curve to the data points. This can clearly be seen in the last plot presented in Fig.(6.3). The reason for this is that most of DDO 154's mass is dominated by a distribution of gasses [39] that can neither be classified as being part of the disk nor a type of bulge. As such in order to correctly test this galaxy and others like it one would have to derive the correct velocity contributions of such gas distributions. The red dashed line is the GR disk contribution which also behaves much worse than the other galaxy morphology types.
FIG. 6: In these plots the full green curve represents the f (T ) gravity rotational velocity curves with n = 1.00001, with α = −4.51 × 10 −34 km 2n kg −1 s −2 , and with masses as shown in Table IV. The red dashed curve represent the general relativistic rotation curve and the data points are the rotational velocity data sets provided in Table V.

VI. CONCLUSION
In this paper we have considered galactic rotation curves in f (T ) gravity while taking f (T ) = T + αT n as the working model of the gravitational action. A weak field metric [23] was utilized for this purpose. We compared various regions of values for the index n with observations of the Milky Way, and three other galaxies which represented other galaxy morphologies. In all cases, the galactic velocity profile was dealt with in terms of a bulge and a disk radial region. This way of decomposing a galaxy better explains the various segments that make up the profile. For every fitted parameter being considered, Eq.(46) was used to determine the combined effect on the final velocity profile.
The Milky Way galaxy was then used to determine which regions of n are more realistic than others. For each index value, we fit the unknown coupling parameter α. Besides the effect on the value of this constant, the various values of n also had an effect on the units of the constant. In the grid figures of Fig.(2), we show plots of all the regions under consideration. The index regions of negative n, n = 0, and 0 < n < 1 are discarded due to their velocity profile never vanishing for very large radii. Moreover, the intermediate region behaves very poorly in several instances of those regions.
The most promising region for n is 1 < n < 3 2 . For values of n close to 1, rotation curves were produced which fit very well with large values of R. These curves also harbored the possibility of eventually going to zero if the limitations of the functions used were to be surpassed. Such a property is also important as it means that galactic potentials would not interfere on the cosmological scale which would cause serious problems in the theory. The only disadvantage of this range of values of n with this type of fitting is that there is a slight overshooting of the velocity profile for the region of R between 0.1 × 10 18 km and 0.4 × 10 18 km. That being said, the curve is still within the error bars for a large part of the data points available. One potential solution would be to consider a variation of the teleparallel Lagrangian being investigated in this work, possibly a combination of power laws. Unfortunately, as of yet metrics for such Lagrangians do not exist in teleparallel gravity but are intended to be developed and tested in the near future.
For the Milky Way, the best fit galactic velocity profile in this formulation of f (T ) gravity was found when applying a multi parameter fitting instead of just a fitting for the coupling constant α. Through this fitting we conclude that n = 1.0001 gives the best fit accompanied by a bulge surface mass density of 5.94 × 10 6 kg km −2 , a disk mass of 4.08 × 10 10 M and a coupling constant of −4.51 × 10 −34 km 2n kg −1 s −2 . With these parameters no problems are observed in the small R limit and it fits perfectly with the whole data range given while conserving the possibility that the curve eventually falls to zero. The resulting masses for the galactic components were also consistent with current luminous matter estimates.
The other three galaxies that were considered were NGC 3198 (bright spiral galaxy [37]), UGC 128 (low surface brightness disk galaxy [40]), and DDO 154 (an irregular dwarf galaxy). These were considered to test the dexterity of the general result for large variance in the components being considered. In Figs.(6.1-6.2), a bright spiral galaxy and a low surface brightness disk galaxy are fit consistently for parameter values also found for the Milky Way. This shows the broadness of the result. In Ref. [41] a series of low surface brightness galaxies were considered in the power-law incarnation of the generalized f (R) class of theories. The results in that case were similarly relatively good when compared to GR. These types of galaxies are supposed to be dark matter dominated in comparison to other types of galaxies however using power law models in either f (R) or f (T ) (present case) the rotation curve profiles can be almost entirely accounted for using a modified gravitational action.
In the case of the irregular galaxy shown in Fig.(6.3), the fit is very poor due to the dominance of a third component in the galaxy dynamics, namely the role of gas. It would be very interesting to model this component using the f (T ) Lagrangian being advanced in this work. In Ref. [42], this very case is considered for the power-law model in f (R) theory. The results are very promising and cast further doubt on the need for a dark matter contribution to describe the galactic rotation curve profile. This would be the natural next step for our analysis.
While the work conducted is promising, it is of paramount importance that such results are tested with regards to larger galaxies surveys that include weighted proportions of the various galaxy morphology types. The parameters obtained here are to be compared with those obtained in each case and if the model allows, a best fit for all galaxy types is to be acquired. This proposed analysis would also test the universal-ity of the coupling parameters, namely n and α. A larger multi-galaxy analysis of the results obtained here is thus intended to be conducted in a future work. Moreover, it would be of utmost interest to compare this model with other modified theories of gravity, such as the f (R) gravity models advanced in Refs. [41,42] among others. Beyond the question of the need for dark matter, the galactic rotation curve profile problem may also advance the question of a preferred model of gravity.
For the region 0 < n < 1, the velocity curve can be represented through the integral Noting that we define a new variable x such that Assuming that the value of n is restricted between but not including 0 and 1, we define the following inverse transformation for the function With this new expression for x 1−n in hand, we substitute it into Eq.(50). Assuming a flat disk and noting that we obtain the following form for this equation Substituting for x and integrating with respect to φ we find that Integrating with respect to k we obtain the following form of the potential this last integral has to be worked through numerically since no analytic solution exists. This is shown in the respective velocity profiles shown in Fig.(2).
Appendix II: Calculating the velocity profile for the range 1 < n < 3 2 For the case of Lagrangian index in the range 1 < n < 3 2 , the velocity curve profile takes on the following effective potential As in the first appendix, we set the following transformation and inverse transformation pair s Given that the integral involves a Bessel function, we need the following relation to move forward in the calculation [8] J 0 (kr) = Considering the galaxy as a flat disk and using ∞ −∞ δ(z )e kz dz = 1, we find that We now note that m = 0 is the only value for which the integral gives a non-zero value for the φ integral. After integrating over R and k as in Ref. [8,28] we obtain the final form of the effective potential contribution