Galactic rotation dynamics in a new f ( R ) gravity model

We propose to test the viability of the recently introduced f ( R ) gravity model in the galactic scales. For this purpose we consider test particles moving in stable circular orbits around the galactic center. We study the Palatini approach of f ( R ) gravity via Weyl transformation, which is the frame transformation from the Jordan frame to the Einstein frame. We derive the expression of rotational velocities of test particles in the new f ( R ) gravity model. For the observational data of samples of high surface brightness and low surface brightness galaxies, we show that the predicted rotation curves are well ﬁtted with observations, thus implying that this model can explain ﬂat rotation curves of galaxies. We also study an ultra diffuse galaxy, AGC 242019 which has been claimed in literature to be a dark matter dominated galaxy similar to low surface brightness galaxies with a slowly rising rotation curve. The rotation curve of this galaxy also ﬁts well with the model prediction in our study. Furthermore, we studied the Tully-Fisher relation for the entire sample of galaxies and found that the model prediction shows the consistency with the data.


I. INTRODUCTION
A longstanding challenge in astrophysics and cosmology is the mystery of dark matter (DM) [1][2][3][4][5].The missing mass problem was first predicted in the early 1930s by J. H. Oort [6] while studying the motion of stars in the Milky Way.Around the same time similar evidences were reported by Swiss astronomer Fritz Zwicky, who studied the Coma Cluster [7,8].Two prime evidences that insist on the existence of DM are the velocity rotation curve of galaxies [9][10][11] and the gravitational lensing [10,12].American astronomer Vera Rubin, did pioneering work by conducting a study of rotation curves of 60 isolated galaxies [13] and established the idea of "missing mass".Using the Planck data [14] on the Cosmic Microwave Background (CMB) radiation, measurements of the cosmological parameters infer that the Universe is made up of ∼ 4 − 5% baryons, ∼ 25% non-baryonic dark matter, and ∼ 70% dark energy.A few particles claimed as DM candidates [15,16] are, namely, weakly interacting massive particles (WIMPs), standard model (SM) neutrinos, sterile neutrinos, axions, supersymmetric candidates (neutralinos, sneutrinos, gravitinos, axinos), etc.However, after almost eight decades since the development of the concept of DM, the DM particle is still missing from the table of elementary particles of nature, i.e. the fundamental nature of DM remains a mystery, and the problem of DM persists.
Over the past few decades various other issues, especially the flatness and horizon problems [24], and the riddle related to the late time cosmic acceleration [17,18] have come to light, which specify that the standard cosmological model based on Einstein's General Relativity (GR) and the particle physics standard model fails to explain the Universe at large scales.This has generated an increasing interest to explore alternative theories of gravity (ATGs) [19][20][21][22][23], where gravitational interactions other than the ones described by GR were proposed.Within the broad area of ATGs, here we refer the modified theories of gravity (MTGs) [25][26][27][28][29] as the ATGs those were proposed to modify GR to provide solutions to emerging issues.The simplest class of MTGs is the f (R) gravity [30,31].In these gravity theories the modification is made to the geometry part of Einstein's field equations.This is done by replacing the Ricci scalar R of the Einstein-Hilbert action with a function f (R) of R.There are two main variational approaches to derive the field equations in f (R) gravity, the metric formalism and the palatini formalism.In the metric formalism, matter is minimally coupled with the metric, and the energy-momentum tensor is independently conserved.In Palatini formalism, the metric as well as the connection are treated as independent variables.Here, the Riemann tensor as well as the Ricci tensor are constructed with the independent connection.Moreover, there are other formalisms which are found in the literature are the metric-affine formalism [40] and the hybrid metric-Palatini formalism [41].In the metric-affine formalism the matter action is considered as variable with respect to connection in contrast to the case of the Palatini formalism.The hybrid metric-Palatini formalism is the combination of the suitable elements of both metric and Palatini formalisms.
A plethora of research works have been dedicated to explain the effects of DM in MTGs [32][33][34][35][36][37][38].As the concept of DM has mainly been indicated by irregularities in the galactic rotation curves, several MTGs have been proposed to study galactic rotation curves.Harko [39] investigated galactic rotation curves in MTGs with non-minimal coupling between matter and geometry.Capozziello et al. [41] studied rotation curves in hybrid metric-Palatini gravity model.Gergely et al. [42] considered the asymptotic behaviour of galactic rotation curves in brane world models.Several other studies [43][44][45][46][47][48][49][50] have been carried out in different MTGs to explain galactic rotation curves.These studies have motivated us to study and explain the rotation curves in a viable f (R) gravity model.Here we propose to investigate the rotation curves of galaxies in the Palatini f (R) gravity taking into consideration the conformal transformation of the metric.It is basically a frame transformation from the Jordan frame to the Einstein frame [51][52][53][54] via a conformal factor.To this aim, we consider the new model [55,56] of f (R) gravity, which has been recently introduced as a viable dark energy model of the theory.The main motive of the present study is to test the viability of this recent model in the galactic scales.We start by considering test particles around galaxies moving in stable circular orbits.The rotational velocity obtained by employing the new f (R) gravity model has the Newtonian term as well as that coming from the modified geometry.The predicted rotational velocity from the model is fitted with observations of a few samples of high surface brightness (HSB), low surface brightness (LSB) and dwarf galaxies.Another addition to our work is the study of the ultra diffuse galaxies (UDGs).They are a fascinating class of galaxies with unusual properties, such as very high [57], or very low content of DM [58].Infact, UDGs are difficult to observe [59][60][61][62] as well as to analyze [63,64].We have investigated the behaviour of the rotation curve of a particular UDG, AGC 242019 [65,66].The rotational velocity predicted the model fits well for this UDG.Furthermore, we also derive the Tully-Fisher relation [67][68][69] for the new f (R) gravity model.
Our work is organised as follows.In section II, we discuss the simplest type of modified gravity, i.e. the f (R) gravity in the Palatini formalism.Here, we obtain the modified field equations in f (R) gravity.The Weyl transformation from the Jordan frame to the Einstein frame and its implications on the Palatini f (R) gravity is also discussed in this section.In section III, we study the dynamics of a test particle around the centres of galaxies in conformally transformed Palatini f (R) gravity.In section IV, we derive the rotational velocity of the test particles in the recently proposed model of f (R) gravity as mentioned above.In section V, we fit the predicted rotational velocity with observations as a test of the new f (R) gravity model in galactic scales.Further, in section VI, we derive the Tully-Fisher relation for the model.Finally, in section VII we conclude and discuss the results of our work.Throughout our work we use the metric signature (-, +, +, +).

II. PALATINI f (R) GRAVITY AND WEYL GEOMETRY A. Field Equations
The action that defines f (R) theories of gravity [30] has the generic form: where f (R) is a function of the Ricci scalar R = g µν R µν and κ 2 = 8πGc −4 = 1/M 2 pl .M pl is the (reduced) Planck mass ∼ 2 × 10 18 GeV.S m is the matter action that is independent of the connection but depends on the metric g µν and the matter field Φ.Here, we will apply the Palatini formalism [54,70].As mentioned earlier, unlike the metric formalism, in the Palatini approach the torsion-free connection Γ α µν and the metric g µν are treated as the dynamical variables to be independently varied.Now, varying the action (1) with respect to the metric g µν we obtain the field equations as where f R (R) is the derivative of f (R) with respect to R and the energy-momentum tensor T µν is given by Again, the variation of the action (1) with respect to the connection Γ α µν gives, It should be noted that when f (R) = R, equation (2) gives the Einstein's field equations and equation ( 4) simply becomes the definition of the Levi-Civita connection in GR.This means that in the limit f (R) = R, the Palatini approach leads to GR as expected.Trace of equation ( 2) is given as This expression shows that R = g µν R µν (Γ) can be algebraically solved in terms of the trace of energy-momentum tensor T ≡ T µ µ = g µν T µν which leads to R = R(T ) and f R = f R (T ) as functions of matter but not of the torsion-free connection.

B. Weyl transformation
Named after Hermann Weyl, the Weyl transformation [52,54] is a frame transformation from the Jordan frame to the Einstein frame.Usually this transformation is essential to get the minimally coupled scalar degree of freedom in the Einstein frame from the conventionally non-minimally coupled one in the Jordan frame in f (R) gravity [51].Moreover, in the Einstein frame the field equations in f (R) gravity can be conveniently written in the form of GR.In this transformation the Jordan frame metric g µν is related to the Einstein frame metric gµν as where Ω is the conformal factor and in the present context it is simply given as Ω = f R (R).As a result, from equation ( 4) it can be noted that the connection for the torsion-free theory is the Levi-Civita connection for the conformally related metric f R (R) g µν [54].At this point it would be appropriate to recast equations (2) for the torsion-free situation as However, in the case of torsion related environments where the connection is asymmetric in nature, the above equation will take a different form containing higher derivatives of f (R), which can be noticed in the Ref. [54].In relation to this it should be mentioned that we consider here the torsion-free Palatini approach based on the fact that our considered model is the dark energy f (R) gravity model satisfying the solar system tests and hence the model satisfies f R (R) − 1 ≪ 1 and also can be seen that f RR (R) ≈ 0 [55].If we viewed the second term on the right hand side of equations ( 7) as the term for the effective cosmological constant Λ e , i.e.
in a spacetime having mass-energy source then the equations can be rewritten as the modified Einstein's field equations in f (R) gravity in the form: where G µν is the usual Einstein tensor.
It is clear that the Palatini field equations (9) are not in the complete form of Einstein's field equations in GR.Now through the conformal transformation (6) this can be achieved by transforming equations ( 9) from the Jordan frame to the Einstein frame as where is the energy-momentum tensor and is the effective cosmological constant in the Einstein frame.Thus the field equations of f (R) gravity in Palatini formalism in the Einstein frame can be written exactly in the form of the corresponding equations in GR with the conformally transformed metric gµν .Further, it is clear from equations ( 8) and ( 11) that the effective cosmological constant in the Einstein frame is different from that in the Jordan frame depending on the conformal factor Ω = f R (R).

III. ROTATIONAL VELOCITY OF TEST PARTICLES AROUND GALAXIES
In order to proceed towards the outcome of our work, we first need to consider a test particle (say, a star) in a galaxy moving in a stable circular orbit [71].The centripetal acceleration a of this test particle is related to its orbital velocity v as Again, as the Einstein's equivalence principle can be justified for a theory of gravity that is conformally related to standard GR, hence the test particle in our case will satisfy the geodesic equation, To relate the orbital velocity of the test particle or star with its geodesic motion we need to consider that although at the center of a galaxy the velocities of stars are very high, however it turns out that these velocities are very low in comparison to the speed of light.So the condition v ≪ c is always satisfied for the motion of stars in the galactic environment.Under this condition, using the coordinates x i = (ct, r, θ, φ) with i = 0, 1, 2, 3, it follows that where x 0 = ct.Now, considering the above condition along with the weak field limit of equation ( 13) and a static spacetime (Γ 0 00 = 0), we get for the radial component, Thus from equations ( 12) and ( 15), we obtain It needs to be mentioned here that in the process of obtaining equation ( 4), the idea of the invariance under the projective transformation of the connection has been implemented [72,73].The projective transformation of the connection is defined as a class of connections related to each other as where S σ σµ ≡ A µ is a vector which generates torsion, and hence is usually set to zero from the beginning in the present study.Moreover, the different connections can define the same geodesics but are parametrized in different ways.The choice of parametrization is related to the choice of metric.In the case of Weyl transformation, the connection is a Levi-Civita connection of the conformal metric gµν .Hence, we can write the connection in equation ( 16) as [54] Γ r 00 = 1 2 grσ (g σ0,0 + g0σ,0 − g00,σ ) = − 1 2 grj ∂ j g00 .
The static spherically symmetric metric in the region exterior to the galactic baryonic mass distribution is given by the following line element: where the metric coefficients Φ and Λ are functions of the radial coordinate r only.For this metric equation ( 18) takes the form: where the prime denotes the derivative with respect r.
In terms of MTGs, the gravitational field equations can be generalized in the form [74,75]: where σ(ϕ) is a coupling factor to gravity and ϕ generally represents either curvature invariants or other fields, such as scalar fields, which adds to the dynamics of the theory.The additional tensor H µν added to the Einstein tensor G µν represents the geometrical modifications which appear as a result of MTGs.In this generalization GR can be recovered as a particular case of MTGs by considering H µν = 0 and σ(ϕ) = 1.Keeping this representation in mind, we can rewrite equation (10) as It is seen that equation ( 22) is comparable to equation ( 21) when H µν = Λe gµν and σ(ϕ) = f R (R).Moreover, in our study we consider the pressureless (p = 0) dust model of the Universe and hence we take the trace of the energy-momentum tensor T = − ρ.Thus from the field equations (22) with the metric equation ( 19) one can rewrite the connection (20) as Now it is required to have the explicit form of the metric coefficients e 2Φ and e 2Λ in terms of the galactic parameters through the f (R) gravity model.Hence from the metric equation ( 19), following the thorough calculations [75], we are able to obtain the metric coefficient e 2Λ in the familiar form as given by where the modified mass distribution is defined as Whereas the exact form of the metric coefficient e 2Φ is very complex.Hence, we employ another way of solving it as follows.
It is reasonable to assume that the metric coefficients e 2Φ and e 2Λ take the Schwarzschild form at large distances compared to the core radius of a galaxy.This form of the coefficient e 2Λ is already seen in equation (24).Hence, from the Ref. [76] we can write the metric gµν for the weak field approximation as where η µν is the metric of the Minkowski spacetime and the first order term coming from the Newtonian limit is given as Thus from equation ( 26), we can write, As we intend to study the behaviour of rotation curves of test particles in a galaxy moving in a stable circular orbit, we require to find an exact form of the connection (23) for the recently introduced model [55] of f (R) gravity.

IV. ROTATIONAL VELOCITY IN THE NEW f (R) GRAVITY MODEL
In the new f (R) gravity model [55,56], where α and β are two dimensionless constants.R c is the characteristic curvature constant with dimension similar to curvature scalar R, and one may expect R c ∽ Λ.For this model, the conformal factor is obtained as Using this conformal factor, and equations ( 24) and ( 28) we are able to derive the exact form of the connection (23) for the aforementioned modified gravity model as Finally, we arrive at the expression of rotational velocity for the new model (29) from equation ( 16) as Here we should mention that because of the complicated form of the modified mass distribution equation (25), in our work we assume a simple mass distribution equation within a galaxy in the form [44]: where M 0 (total mass of the galaxy) and r c (core radius) are the parameters to be predicted by fitting the calculated rotational velocities with observed data.R 0 is the scale length of the galaxy.And, the values of the parameter γ [41] are γ = 1 for high surface brightness (HSB) galaxies and γ = 2 for low surface brightness (LSB) and dwarf galaxies.In this section we check the validity of the new f (R) gravity model ( 29) by comparing its theoretical predictions with the observational data on the galactic rotation curves.Using equation (32), we try to obtain the flat rotation curves for a few samples of HSB, LSB and dwarf galaxies in the subsections that follow.We also investigate the rotation curve of an ultra diffuse galaxy (UDG) AGC 242019.For this analysis, we have used a well constrained set of model parameters viz.R Rc = 1.5, α = 0.005 and β = 0.044 [55].

A. Analysis of HSB galaxy sample
We fit the predictions of the new f (R) gravity model (29) as obtained from equation (32) to the data of a sample of nine HSB galaxies extracted from Ref. [77].In these fittings and subsequent ones we use the χ 2 minimization technique.Table I shows the related data for the nine HSB galaxies along with the respective best fit values of the parameters M 0 and r c and the reduced chi-squared (χ 2 red ) values.The results are depicted in Fig. 1.In Table I, the data of HI gas mass M HI (which extends well beyond the optical disk distributions) in units of 10 10 M ⊙ , the disk scale length R 0 measured in units of kpc, the adopted distance D in Mpc and the B-band luminosity in units of 10 10 L ⊙ are extracted from Ref [77].The total gas mass is given by 4  3 M HI [44,[77][78][79] and the stellar mass is obtained by subtracting the gas mass from the predicted total mass of the galaxy.We have also listed the stellar mass-to-light ratios (M/L) * .It is seen that in almost all the cases, the mass-to-light ratios inferred from the fit are reasonably close to the expected mass-to-light ratios.One should note that the expected mass-to-light ratios of galaxies lie within the range of 0.2 M ⊙ /L ⊙ < M/L < 10 M ⊙ /L ⊙ [81].Also, Fig. 1 depicts that almost all of the selected galaxies are well-fitted by the recently proposed f (R) model ( 29), except the galaxy NGC 6946 for which the χ 2 red value is very large, around 15.13.However, as one can see from Fig. 1, the shape of its rotation curve was not affected.

B. Analysis of LSB galaxy sample
It is well known that LSB galaxies are regarded as dark matter dominated galaxies and hence, they can provide a satisfactory test for a gravitational theory [77,80].To this end, we fit equation (32) to a sample of 21 LSB galaxies extracted from Ref. [77] with their maximum radial distance varying from 2.1 kpc to 18.2 kpc.Table II compiles the relevant data for these galaxies we have selected from the 111 galaxies studied in Ref. [77] along with the respective best fit values of the parameters M 0 and r c and the χ 2 red values of the fits.The results are presented in Fig. 2 and 3.It is seen from Table II that the mass-to-light ratios for few galaxy samples are much larger than the upper bound 10(M ⊙ /L ⊙ ) [81].However, as can be seen from Figs. 2 and 3, the rotation curves are well fitted with the observed data as the χ 2 red values are smaller than or equivalent to 1 for almost all of the galaxies in the sample.For the galaxies NGC 0300, NGC 3109 and UGC 6818 although the χ 2 red values are high, around 2.04, 2.78 and 1.69 respectively, these do not affect the shape of the fitted rotation curves.32) for the new f (R) gravity model (29) to the rotational velocities (in km/s) of a set of nine HSB galaxies extracted from Ref. [77] with their errors plotted as a function of radial distance (in kpc).

C. Analysis of dwarf galaxy sample
For dwarf galaxies, we have chosen a sample of nine galaxies taken from Ref. [82].Table III compiles the relevant data for the chosen galaxy sample and the results are depicted in Fig. 3.For all the galaxies the χ 2 red values are small and gives well fitted galactic rotation curves, except for two galaxies viz., UGC 7323 and UGC 7524.Furthermore, the mass-to-light ratio for most of the galaxies are within the upper bound.However, for two galaxies namely, UGC 4499 and UGC 7559, the mass-to-light ratios are notably greater than the upper bound value as expected from the population synthesis models [81].According to the population synthesis models depending on the history of the star formation and also on metallicities the blue-band mass-to-light ratio may range from a few tenths to 10.

D. Analysis of an ultra diffuse galaxy
Here, we analyze the rotation curve of an ultra diffuse galaxy (UDG), AGC 242019, which was identified by the Arecibo Legacy Fast ALFA (ALFALFA) survey of HI galaxies [83].This galaxy has an HI mass of 8.51 × 10 8 M ⊙ and a corresponding distance of 30.8 Mpc.The gas rich UDG AGC 242019 has been claimed to be like an observed LSB galaxy with a slowly rising rotation curve [66].Fig. 4 shows the fitted rotation curve with the data taken from Ref. [66].It can be seen that the curve of this galaxy is well fitted with observation.In fact, from Fig. 4 one can see that the rotation curve is slowly rising, indicating that UDG AGC 242019 can be considered as a member of a class of LSB galaxies.However, further study in this regard is still required to confirm such proclamation.From our study, the χ 2 red value is obtained as 6.84 and the fitted parameters are evaluated to be M 0 = 4.92 (10 10 M ⊙ ) and r c = 2.29 (kpc).Thus, the flat rotation curve of UDG AGC 242019 can be explained via the new model (29) of f (R) gravity.32) for the new model (29) to the rotational velocities (in km/s) of a set of 12 LSB galaxies extracted from Ref. [77] with their errors plotted as a function of radial distance (in kpc).32) for the new model (29) to the rotational velocities (in km/s) of a set of 9 LSB galaxies extracted from Ref. [77] with their errors plotted as a function of radial distance (in kpc).

VI. THE TULLY-FISHER RELATION
R. B. Tully and J. R. Fisher in 1977 published a method of determining the distances of spiral galaxies based on their empirical or observational relation, now known as the Tully-Fisher relation [67].It implies a relation between the luminosity of a galaxy and the velocity of the outermost observed point of the galaxy as where L is the observed luminosity of a galaxy in units of 10 10 L ⊙ , v out is the velocity at the outermost observed radial point of the galaxy in units of km/s, ξ is the proportionality constant and a is another constant.Taking logarithm on both sides of this equation (34), we may write  32) for the new model (29) to the rotational velocities (in km/s) of a set of nine dwarf galaxies extracted from Ref. [77] with their errors plotted as a function of radial distance (in kpc).
and the horizontal axis is the base 10 logarithm of the velocity of the outermost observed radial point in units of kms −1 .The B-band luminosity data for each galaxy sample has been extracted from Ref. [77] and are tabulated in Tables I, II and III.The solid line in the plot is the best least-square fitting of equation (35) to the observed corresponding data of galaxies as mentioned.The fitting is found to be very good with χ 2 red = 0.127 and the parameters a = 3.18 ± 0.27 and b = − 6.75 ± 0.54.Now, from the mean mass-luminosity relation M/L = constant [68] and equation (35) we can write the Tully-Fisher relation in terms of the predicted mass (in units of 10 10 M ⊙ ) and the velocity (in units of kms −1 for the new f (R) gravity model (29) as The right panel of Fig. 5 shows the best fit Tully-Fisher relation parametrized by equation (36) for the new model (29).The vertical axis denotes the base 10 logarithm of the predicted mass (in units of 10 10 M ⊙ ) and the horizontal axis denotes the base 10 logarithm of the fitted velocity (in units of kms −1 ).In this fit the parameters are found to be a = 3.22 ± 0.48 and b = − 7.62 ± 0.95 with χ 2 red = 0.397.(35) to the observed data.The right panel shows predicted Tully-Fisher relation for the masses and velocities of the outermost radial points of galaxies predicted from the new f (R) model (29).Here the solid line represents the best least-square fitting of (36) to the corresponding predicted data as mentioned.

VII. SUMMARY AND FINAL REMARKS
The behaviour of galactic rotation curves indicates the need for DM or some modifications of GR.However, till date there is no direct evidence of the existence of DM.Also, DM interacts only via gravity and hence, the question arises as to whether the effect of DM is a consequence of the modification of gravity only.Keeping this point in mind, in this work, we have employed one of the simplest MTG, the f (R) gravity to study the galactic rotation curves for some of the galaxies observed experimentally in recent times.In doing so we have mainly tried to test the viability of the recently proposed model [55,56] of f (R) gravity in the galactic scales.
In our present work, we have used the Palatini formalism along with the Weyl transformation.Thus the variation of the action for the f (R) theory first with respect to the metric and then with respect to the connection yields the field equations in the theory.We consider the Weyl transformation which is the frame transformation of the spacetime metric g µν from the Jordan frame to the Einstein frame to get the field equation in the convenient form of GR.Moreover, the connection in this torsion-free Palatini approach is the Levi-Civita connection for the conformally related metric and as such particles moving in the gravitational field follow the geodesics obtained by the connection.In the Einstein frame, we study the rotational velocity of a test particle moving in a stable circular orbit.A static spherically symmetric metric is taken into consideration and we are able to derive the co-efficients of the metric for our galactic model spacetime.
Next, we introduce the new model of f (R) gravity.As mentioned already this is a recently introduced f (R) gravity model and hence we try to test the viability of this model in the galactic scales.Making use of the metric co-efficients and the field equations we derive the analytical expression for rotational velocities of test particles in the stable circular orbits.It can be observed that the expression of the rotational velocity obtained is different from the Newtonian one.This is due to the fact that the rotational velocity obtained via our present approach has extra terms coming from the geometrical modifications in the theory.This rotational velocity expression is fitted with observed data of a few samples of galaxies.We consider nine samples of high surface brightness (HSB) galaxies, 21 samples of low surface brightness (LSB) galaxies and nine dwarf galaxies.The galactic rotation curves are well fitted with observations for the new f (R) gravity model.Although, for a few samples of galaxies the mass-to-light ratios are found to be much larger than the expected upper bound, however the χ 2 red values for most of the samples are smaller or equivalent to 1, thus indicating the well fitted rotation curves.This shows the viability of the new f (R) gravity model in the galactic scales.
In addition to these samples, we take into consideration an interesting class of galaxies i.e., the ultra diffuse galaxies (UDGs).They are fascinating objects which are either made entirely of DM or have a high content of DM.This makes them difficult to observe and analyze.We consider one such galaxy, AGC 242019, which has been claimed to be similar to low surface brightness galaxies with a slowly rising rotation curve.Our study supports this proclamation.However, further studies are necessary in this regard.
Finally, we have studied the Tully-Fisher relation in the new f (R) gravity model.We have studied the luminosity as well as total mass of galaxies as a function of velocity of the outermost observed radial point.The entire sample of galaxies have been combined for this purpose and the fits show a consistent result across the galaxies.In our work, we have made use of the SPARC catalogue (http://astroweb.cwru.edu/SPARC/)for the observational data of the sample of galaxies.
Lastly, it would be interesting to extend our work by including different MTGs with more observational data and also by comparing our results with the standard DM profiles such as Navarro-Frenk-White (NFW) [84,85] and Burkert [86] profiles.Also, in the near future, we can try to obtain much feasible mass-to-light ratios within the expected bound.

FIG. 1 :
FIG.1: Fitting of equation(32) for the new f (R) gravity model(29) to the rotational velocities (in km/s) of a set of nine HSB galaxies extracted from Ref.[77] with their errors plotted as a function of radial distance (in kpc).

FIG. 2 :
FIG.2: Fitting of equation(32) for the new model(29) to the rotational velocities (in km/s) of a set of 12 LSB galaxies extracted from Ref.[77] with their errors plotted as a function of radial distance (in kpc).

FIG. 3 :
FIG.3: Fitting of equation(32) for the new model(29) to the rotational velocities (in km/s) of a set of 9 LSB galaxies extracted from Ref.[77] with their errors plotted as a function of radial distance (in kpc).

FIG. 4 :
FIG.4: Fitting of equation(32) for the new model(29) to the rotational velocities (in km/s) of a set of nine dwarf galaxies extracted from Ref.[77] with their errors plotted as a function of radial distance (in kpc).

FIG. 5 : 62 FIG. 6 :
FIG.5:The rotational velocities of AGC 242019 fitted to equation(32) with the associated errors plotted as a function of radial distance in kpc.The data for the rotational velocities of the galaxy is taken from Ref.[66].

TABLE I :
(29)vant galaxy properties for a set of nine HSB galaxies and the best fit values for the parameters M0 and rc of these galaxies in the new f (R) gravity model(29).
V. ASTROPHYSICAL TESTS OF THE NEW f (R) GRAVITY MODEL IN GALACTIC LEVELS

TABLE II :
(29)vant galaxy properties for a set of 21 LSB galaxies and the best fit values for the parameters M0 and rc of these galaxies in the new f (R) gravity model(29).

TABLE III :
(29)vant galaxy properties for a set of nine Dwarf galaxies and the best fit values for the parameters M0 and rc of these galaxies in the new f (R) gravity model(29).