Kinematics of Radion field: A possible source of dark matter

The discrepancy between observed virial and baryonic mass in galaxy clusters have lead to the missing mass problem. To resolve this, a new, non-baryonic matter field, known as dark matter has been invoked. However, till date no possible constituents of the dark matter components are known. This has led to various models, by modifying gravity at large distances to explain the missing mass problem. The modification to gravity appears very naturally when effective field theory on a lower dimensional manifold, embedded in a higher dimensional spacetime is considered. It has been shown that in a scenario with two lower dimensional manifolds separated by a finite distance is capable to address the missing mass problem, which in turn determines the kinematics of the brane separation. Consequences for galactic rotation curves are also described.


Introduction
Recent astrophysical observations strongly suggest existence of non-baryonic dark matter at the galactic as well as extra-galactic scales (if the dark matter is baryonic in nature, the third peak in the Cosmic Microwave Background power spectrum would have been lower compared to the observed height of the spectrum [1]). These observations can be divided into two branches -(a) behavior of galactic rotation curves and (b) mass discrepancy in clusters of galaxies [2].
The first one, i.e., rotation curves of spiral galaxies, show clear evidences of problems associated with Newtonian and general relativity prescriptions [2][3][4]. In these galaxies neutral hydrogen clouds are observed much beyond the extent of luminous Baryonic matter. In Newtonian description, the equilibrium of these clouds moving in a circular orbit of radius r is obtained through equality of centrifugal and gravitational force. For cloud velocity v(r), the centrifugal force is given by v 2 /r and the gravitational force by GM (r)/r 2 , where M (r) stands for total gravitational mass within radius r. Equating these two will lead to the mass profile of the galaxy as, M (r) = rv 2 /G. This immediately posed serious problem, for at large distances from the center of the galaxy, the velocity remains nearly constant v ∼ 200 km/s, which suggests that mass inside radius r should increase monotonically with r, even though at large distance very little luminous matter can be detected [2][3][4].
The mass discrepancy of galaxy clusters also provides direct hint for existence of dark matter. The mass of galaxy clusters, which are the largest virialized structures in the universe, can be determined in two possible ways -(i) from the knowledge about motion of the member galaxies one can estimate the virial mass M V , secondly, (ii) estimating mass of individual galaxies and then summing over them in order to obtain total baryonic mass M . Almost without any exception M V turns out to be much large compared to M , typically one has M V /M ∼ 20 − 30 [2][3][4].
Recently, new methods have been developed to determine the mass of galaxy clusters, these are -(i) dynamical analysis of hot X-ray emitting gas [5] and (ii) gravitational lensing of background galaxies [6] -these methods also lead to similar results. Thus dynamical mass of galaxy clusters are always found to be in excess compared to their visible or baryonic mass. This missing mass issue can be explained through postulating that, every galaxy and galaxy cluster is embedded in a halo made up of dark matter. Thus the difference M V − M is originating from the mass of the dark matter halo, the galaxy cluster is embedded in.
The physical properties and possible candidates for dark matter can be summarized as follows: dark matter is assumed to be non-relativistic (hence cold and pressure-less), interacting only through gravity. Among many others, the most popular choice being weakly interacting massive particles. Among different models, the one with sterile neutrinos (with masses of several keV) has attracted much attention [7,8]. Despite of few successes it comes with its own limitations. In the sterile neutrino scenario the X-ray produced from their decay can enhance production of molecular hydrogen and thereby speeding up cooling of gas and early star formation [9]. Even after a decade long experimental and observational efforts no non-gravitational signature for the dark matter has ever been found. Thus a priori the possibility of breaking down of gravitational theories at galactic scale cannot be excluded [10][11][12][13][14][15][16][17][18].
A possible and viable way to modify the behaviour of gravity in our four dimensional spacetime is by introducing extra spatial dimensions. The extra dimensions were first introduced to explain the hierarchy problem (i.e., observed large difference between the weak and Planck energy scales) [19][20][21]. However the initial works did not incorporate gravity, but used large extra dimensions (and hence large volume factor) to reduce the Planck scale to TeV scale. Introduction of gravity, i.e., warped extra dimensions drastically altered the situation. In [22] it was first shown that anti-de Sitter solution in higher dimensional spacetime (henceforth referred to as bulk) leads to exponential suppression of the energy scales on the visible four dimensional embedded sub-manifold (called as brane) thereby solving the hierarchy problem. Even though this scenario of warped geometry model solves the hierarchy problem, it also introduces additional correction terms to the gravitational field equations, leading to deviations from Einstein's theory at high energy, with interesting cosmological and black hole physics applications [23][24][25][26][27][28][29][30][31][32][33]. This conclusion is not bound to Einstein's gravity alone but holds in higher curvature gravity theories 1 as well [30,31,41]. Since the gravitational field equations get modified due to intro- 1 In addition to introduction of extra dimensions we could also modify the gravity theory without invoking ghosts, which uniquely fixes the gravitational Lagrangian to be Lanczos-Lovelock Lagrangian. These Lagrangians have special thermodynamic properties and also modifies behaviour of four-dimensional gravity [34][35][36][37][38][39][40]. However in this work we shall confine ourselves exclusively within the framework of Einstein gravity and shall try to explain the missing mass problem from kinematics of the radion field. duction of extra dimensions it is legitimate to ask, whether it can solve the problem of missing mass in galaxy clusters. Several works in this direction exist and can explain the velocity profile of galaxy clusters. However they emerge through the following setup: • obtaining effective gravitational field equations on a lower dimensional hypersurface, starting from the full bulk spacetime, which involves additional contributions from the bulk Weyl tensor. The bulk Weyl tensor in spherically symmetric systems leads to a component behaving as mass and is known as "dark mass" (we should emphasize that this notion extends beyond Einstein's gravity and holds for any arbitrary dimensional reduction [30,31,41]). It has been shown in [42] that introduction of the "dark mass" term is capable to yield an effect similar to the dark matter. Some related aspects were also explored in [43][44][45][46], keeping the conclusions unchanged.
• In the second approach, the bulk spacetime is always taken to be anti-De Sitter such that bulk Weyl tensor vanishes. Unlike the previous case, which required S 1 /Z 2 orbifold symmetry, arbitrary embedding has been considered in [47] following [48]. This again introduces additional corrections to the gravitational field equations. These additional correction terms in turn lead to the observed virial mass for galaxy clusters.
However all these approaches are valid for a single brane system. In this work we generalize previous results for a two brane system. This approach not only gives a handle on the hierarchy problem at the level of Planck scale but is also capable of explaining the missing mass problem at the scale of galaxy clusters. Moreover, in this setup the additional corrections will depend on the radion field (for a comprehensive discussion see [27]), which represents the separation between the two branes. Hence in our setup the missing mass problem for galaxy clusters can also shed some light on the kinematics of the separation between the two branes. Further the same setup is also shown to explain the observed rotation curves of galaxies as well. Hence both the problems associated with dark matter, namely, the missing mass problem for galaxy clusters and the rotation curves for galaxies can be explained by the two brane system introduced in this work via the kinematics of the radion field.
The paper is organized as follows -In Section 2, after providing a brief review of the setup we have derived effective gravitational field equations on the visible brane which will involve additional correction terms originating from the radion field to modify the gravitational field equations. In Section 3 we have explored the connection between the radion field, dark matter and the mass profile of galaxy clusters using relativistic Boltzmann equations along with Section 4 describing possible applications. Then in Section 5 we have discussed the effect of our model on the rotation curve of galaxies while Section 6 deals with a few applications of our result in various contexts. Finally, we conclude with a discussion on our results.
Throughout our analysis, we have set the fundamental constant c to unity. All the Greek indices µ, ν, α, . . . run over the brane coordinates. We will also use the standard signature (− + + . . .) for the spacetime metric.

Effective Gravitational Field Equations on the Brane
The most promising candidate for getting effective gravitational field equations on the brane originates from Gauss-Codazzi equation. However these equations are valid on a lower dimensional hypersurface (i.e., on the brane) embedded in a higher dimensional bulk. Hence this works only for a single brane system. But, the brane world model, addressing hierarchy problem requires existence of two branes, where the above method is not applicable. To tackle the problem of two brane system we need to invoke the radion field (i.e., separation between two branes), which has significant role in the effective gravitational field equations. The bulk metric ansatz incorporating the above features takes the following form, respectively. Use of the expressions for (1) E µ ν and (1) K µ ν in Eq. (2) leads to the effective field equations on the visible brane (i.e., the brane on which Planck scale is exponentially suppressed) in this scenario as [49] (4) G µ ν = where The scalar field Φ(x) appearing in the above effective equation is directly connected to the radion field d 0 (x) (representing proper distance between the branes) such that ω(Φ) and Φ has the following expressions [49] We will assume d 0 (x), the brane separation to be finite and everywhere nonzero. This suggests that Φ(x) should always be greater than zero and shall never diverge. Finally we also have a differential equation satisfied by Φ from the trace of Eq. (7), which can be written as [49] where ω(Φ) has been defined in Eq. (9) and T vis and T hid stands for the trace of energy momentum tensor on the hidden and visible branes respectively. In the above expressions D µ stands for the four-dimensional covariant derivative, also D 2 Φ stands for D µ D µ Φ and (DΦ) 2 = D µ ΦD µ Φ. The above effective field equations for gravity have been obtained following [49], where no stabilization mechanism for the radion field was proposed. In this work as well we would like to emphasize that, we are working with the radion field in absence of any stabilization mechanism. However as already emphasized in [49], in order to provide a possible resolution to the gauge hierarchy problem one requires to stabilize the radion field. Even though we will not explicitly invoke a stabilization mechanism, but will outline how stabilization can be achieved and argue that it will not drastically alter the results.
In such a situation with stabilized radion field, the field Φ(x) appearing in the above equations can be thought of as fluctuations of the radion field around its stabilized value [50]. In particular stabilization of the radion field can be achieved by first introducing a bulk scalar field following [51] and then solving for it. Substitution of the solution in the action and subsequent integration over the extra spatial dimension, leads to a potential for Φ. The same will appear in the above equations through the projection of the bulk energy momentum tensor, which would involve the bulk scalar field and shall lead to an additional potential on the right hand side of the above equations, whose minima would be the stabilized value for Φ = Φ c . Choosing Φ = Φ c + Φ(x), where Φ(x) represents small fluctuations around the stabilized value, one ends up with similar equations as above with bulk terms having similar contributions as from T (vis) and T (hid) respectively. Thus the final results, to leading order, will remain unaffected by introduction of a stabilization mechanism. Even though the fact that virial mass of galaxy clusters scale with r will hold, the sub-leading correction terms in case of galactic motion will change due to presence of a stabilization mechanism due to appearance of extra bulk inherited terms in the above equations. It would be an interesting exercise to work out the above steps explicitly and obtain the relevant corrections due to stabilization mechanism, which we will pursue elsewhere.
As illustrated above for the two brane system the non-local terms get mapped to the radion field, the separation between the two branes. Hence ultimately one arrives at a system of closed field equations for a two brane system. The field equations as presented in Eq. (8) are closed since the radion field Φ satisfies its own field equation Eq. (10). Hence the problematic non-local terms in a single brane approach gets converted to the radion field in a two brane approach and makes the system of gravitational field equations at low energy closed.
We are mainly interested in spherically symmetric spacetime, in which generically the line element takes the following form: This particular form of the metric is used extensively in various physical contexts, for example in obtaining black hole solution, particle orbit, perihelion precession of planetary orbits, bending of light and in various other astrophysical phenomenon [52][53][54]. Given this metric ansatz we can compute all the derivatives of the scalar field and being a static situation, the brane separation is assumed to depend on radial coordinate only. Thus we will only have terms involving derivative with respect to r (these will be denoted by prime). First we can rewrite the scalar field equation, which will be a differential equation for Φ. We will also assume that there is no matter on the hidden brane, but only on the visible brane, which is assumed to be perfect fluid. Thus on the visible brane we have energy momentum tensor to be, T , with the trace being given by, T = −ρ + p + 2p ⊥ . From now on we will remove the label 'vis' from the energy momentum tensor, since only on the visible brane energy momentum tensor is non-zero. With these inputs and the above spherically symmetric metric ansatz we obtain the scalar field equation as, Having derived the scalar field equation, next we need to obtain the field equations for gravity with the metric ansatz given by Eq. (11). These will be differential equations for ν(r) and µ(r) respectively. We can separate out the time-time component, radial component and transverse components leading to where primes denote derivatives with respect to radial coordinate. In the above field equations along with the perfect fluid, we have contributions from the brane cosmological constant. Here we have inserted a brane energy density ρ 0 , where ρ 0 and brane cosmological constant is related via ρ 0 = Λ/8πG. Here G is the four dimensional gravitational constant. Finally, we have contribution from the radion field itself, since it appears on the right hand side of gravitational field equations. Having derived the field equations we will now proceed to determine the effect of the radion field on the kinematics of galaxy clusters and hence its implications for the missing mass problem.

Virial Theorem in Galaxy Clusters, Kinematics of The Radion Field and Dark Matter
It is well known that the galaxy clusters are the largest virialized systems in the universe [2]. We will further assume them to be isolated, spherically symmetric systems such that the spacetime metric near them can be presented by the ansatz in Eq. (11). Galaxies within the galaxy cluster are treated as identical, point particles satisfying general relativistic collision-less Boltzmann equation.
The Boltzmann equation requires setting up appropriate phase space for a multi-particle system along with the corresponding distribution function f (x, p), where x is the position of the particles in the spacetime manifold with its four-momentum p ∈ T x , where T x is the tangent space at x. Further the distribution function is assumed to be continuous, non-negative and describing a state of the system. The distribution function is defined on the phase space, yielding the number dN of the particles of the system, within a volume dV located at x and have fourmomentum p within a three surface element d − → p in momentum space. All the observables can be constructed out of various moments of the distribution function. Further details can be found in [42].
For the static and spherically symmetric line element as in Eq. (11) the distribution function can depend on the radial coordinate only and hence the relativistic Boltzmann equation reduces to the following form [42], The spherical symmetry of the problem requires the coefficient of cot θ to identically vanish. Hence the distribution function can be a function of r, u r and u 2 θ + u 2 φ only. Multiplying the above equation by mu r du, where m stands for galaxy mass and du is the velocity space element we find after integrating over the cluster [42], where R stands for the radius of the galaxy cluster. Using the distribution function, the energy momentum tensor of the matter becomes, which leads to the following expressions for energy density and pressure as, Using these expressions for energy density and pressure in the gravitational field equations presented in Eq. (13), Eq. (14) and Eq. (15) and finally adding all of them together we arrived at To obtain Eq. (20), we have used the expression for trace of the energy momentum tensor. We also remember that ρ 0 stands for vacuum energy density. At this stage it is useful to introduce certain assumptions, since actually we are interested in a post-Newtonian formulation of the effective gravitational field equations. The two assumptions are -(a) ν and λ are small so that any quadratic expressions constructed out of them can be neglected in comparison to the linear one. Secondly, (b) the velocity of the galaxies are assumed to be much smaller compared to the velocity of light, which suggests, such that all the pressure terms can be neglected in comparison to the energy density. Applying all these approximation schemes Eq. (20) can be rewritten as, We can also perform the same schemes of approximation to Eq. (17), which leads to, where K stands for the total kinetic energy of the galaxies within the galaxy cluster and has the following expression, The mass within a small volume of radial extent dr has the expression dM (r) = 4πr 2 ρdr, where in this and subsequent expressions ρ will indicate ρ(r). Thus total mass of the system can be given by integral of dM (r) over the full size of the galaxy. The main contribution comes from mass of intra-cluster gas and stars along with other particles, e.g., massive neutrinos. We can also define the gravitational potential energy Ω of the cluster as, Finally multiplying Eq. (21) by r 2 and integrating from 0 to r, we arrive at, where we have defined: This object captures all the effect of the radion field on the gravitational mass distribution of galaxy clusters and thus may be called as the "radion mass". Note that the "radion mass" defined in this work, is a completely different construct compared to the "dark mass" used in the literature. The dark mass appears from non-local effects of the bulk, specifically through the bulk Weyl tensor in the effective field equation formalism. However, in this work, we have used effective equation formalism for a two brane system as developed in [49], where the correction to gravitational field equations originate from radion dynamics. Pursuing this effective equations further, through virial theorem we have shown that the effect of radion dynamics can be summarized by introducing a radion mass as in Eq. (26). Hence conceptually and structurally the dark mass of [42] is completely different from our "radion mass". Further, the total baryonic mass of the galaxy cluster within a radius r can be obtained by integrating the energy density over the size of the galaxy cluster, which leads to, M (r) = 4π r 0 r 2 ρ(r)dr, using which we finally arrive at the following form for Eq. (25): Earlier we have defined the gravitational potential associated with M , the baryonic mass. We can define an identical object using the radion mass as well, leading to an potential term Ω Φ . Given the potentials we can introduce three radius -(a) R V , the virial radius, obtained using total baryonic potential and baryonic mass, (b) R I , the inertia radius, obtained from moment of inertia of the galaxy cluster and finally (c) R Φ , the radion radius obtained from the radion mass (for detailed expressions see Appendix A). Using these expressions and the definition for virial mass, M V = 2KR V /G, yields the following expression, For most of the clusters, the virial mass M V is three times compared to the baryonic mass M and thus for all practical purposes the first term inside the square root, which is of order unity can be neglected with respect to the other two. The second term yields the contribution from the brane cosmological constant, which is several orders of magnitude smaller compared to the observed mass and thus can also be neglected. Finally, the virial mass turns out to be, Among the various terms in the above expression, virial mass M V is determined from the study of velocity dispersion of galaxies within the cluster and is much large than the visible mass. The above expression shows that if the radion field kinematics is such that M tot is equal to M Φ . Then that in turn will lead to the correct virial mass of the galaxy clusters. The effect of radion field and hence of extra dimension can also be probed through gravitational lensing.
To see that, let us explore the differential equation for Φ, which has not yet been considered. Solving that will lead to some leading order behaviour of the radion field Φ, which in turn would affect M Φ . Thus crucial thing is whether M Φ behaves as r at large distance from the core of the cluster. In which case from the above equation, we readily observe that the galaxy virial mass would also scale as M V ∼ r explaining the issue of dark matter and galaxy rotation curve. To answer all these let us start by using the differential equation for Φ. There we will work under same approximation schemes, i.e., will be neglecting all the quadratic terms, e.g., ν Φ , Φ 2 , will set e λ ∼ 1 and shall neglect vacuum energy contribution ρ 0 to obtain (for general expression see Appendix A), Multiplying both sides by r 2 and integrating twice we obtain (noting that Φ 2 should not contribute) Here M (r) stands for the mass of the baryonic matter within radius r and we know from observations that the density of the baryonic matter falls as ρ c (r c /r) β , where 3 > β > 2 and r c stands for the core radius of the cluster. Thus it is straightforward to compute the mass profile, which goes as ∼ r 3−β , except for some constant contribution. Hence finally after integration we obtain the radion field to vary with the radial distance as r 2−β . However note that the mass of the radion field, i.e., M Φ under these approximations (matter is non-relativistic and field is weak) can be obtained as, Thus the radion mass indeed scales linearly with radial distance which would correctly reproduce the observed virial mass of the galaxy cluster. Due to the linear nature of the virial mass, the velocity profile does not die down at large r as expected. Hence the radion field kinematics can explain the kinematics of the galaxy cluster very well and thus the missing mass problem can be described without invoking any additional matter component. Before concluding the section, let us briefly mention about the connection of the above formalism with the gauge hierarchy problem. The separation between the two branes is denoted by d, which varies with the radion field Φ, logarithmically (see Eq. (9)). The radion field except for a constant contribution varies weakly with radial distance and hence leads to very small corrections to the distance d between the branes. Thus the graviton mass scale for the visible brane will be suppressed by a similar exponential factor as in the original scenario of Randall and Sundrum [22,55], leading to a possible resolution of the gauge hierarchy problem. Thus as advertised earlier, existence of an extra spatial dimension leads to a radion field, producing a possible explanation for the dark matter in galaxy clusters along with solving the gauge hierarchy problem.

Application: Cluster mass profiles
In the previous section we have discussed galaxy clusters by assuming them to be bound gravitational systems, with approximate spherical symmetry and virialized, i.e., in hydrostatic equilibrium. With these reasonable set of assumptions we have shown that, the mass of clusters receive additional contribution from the kinematics of the radion field and provides an alternative to the missing mass problem. In this section we will discuss one application of the above formalism, namely the mass profile of galaxy clusters and possible experimental consequences. We again start from collisionless Boltzmann equation in spherical symmetry and in hydrostatic equilibrium to read, Here V (r) stands for the gravitational potential of the cluster, σ r and σ θ,φ are the mass weighted velocity dispersions in the radial and tangential directions, respectively, with ρ gas being the gas density. For spherically symmetric systems, σ r = σ θ,φ and the pressure profile becomes, P (r) = σ 2 r ρ gas (r). Further if the velocity dispersion is assumed to have originated from thermal fluctuations, for a gas sphere with temperature profile T (r), the velocity dispersion becomes, σ 2 r = k B T (r)/µm p , where k B is the Boltzmann constant, µ 0.609 is the mean mass and m p being the proton mass. Thus Eq. (33) can be rewritten as, The potential can be divided into two parts, the Newtonian potential and the potential due to radion field. As multiplied by (4/3)r 2 /G, the Newtonian potential leads to the Newtonian mass M N , which includes mass of gas, galaxies and in particular of the CD galaxies. Thus finally we obtain the mass profile of a virialized galaxy cluster to be, Thus one needs two experimental input, the observed gas density profile, ρ gas and the observed temperature profile T (r). Gas density can be obtained from the characteristic properties of observed X-ray surface brightness profiles, similarly from X-ray spectral analysis one obtains the radial profile of temperature. Thus from X-ray analysis one can model the galaxy distribution and obtain the baryonic contribution to the mass of the galaxy cluster. From the difference between virial mass and the above estimate one can obtain the contribution due to radion field. At the leading order the radion mass scales linearly with the radial distance with its coefficients being O( /κ 2 ). Thus an estimate of the radion mass will lead to a possible value for /κ 2 . Assuming the bulk gravitational constant to be at the Planck scale one can possibly constrain the bulk curvature scale.

Effect on galaxy rotation curves
Having described a possible resolution of the missing mass problem in connection with galaxy clusters, let us now concentrate on the rotation curves of galaxies. To perform the same we would invoke some general Lie groups of transformation on a vacuum brane spacetime. In particular we will assume the metric to be static and spherically symmetric (i.e., expressed as in Eq. (11)), such that, £ ξ g µν = ψ(r)g µν , where the vector field ξ µ can be time dependent. These are known as conformally symmetric vacuum brane model and we consider angular velocity of a test particle in visible (i.e., negative tension) brane, which can be determined in terms of the conformal factor ψ(r). The above essentially amounts to the assumption that each brane is conformally mapped onto itself along the vector field ξ µ [56][57][58]. It turns out that both the metric and the vector field ξ µ has the following expressions upon solving the relation £ ξ g µν = ψ(r)g µν as, where k is a separation constant and B and C are integration constants. Substitution of these metric functions in the gravitational field equations presented in Eq. (8) leads to, where 'prime' denotes differentiation by radial coordinate r. Multiplying Eq. (40) by 2 and adding it to Eq. (39) one can readily equate it to Eq. (38) resulting into the following differential equation satisfied by ψ(r), The above differential equation can be readily solved, yielding r = r(ψ) [56][57][58]. However the solution depends on the mutual dependence of k on B. We will use galaxy rotation curves as the benchmark to determine the region of interest in the (k, B) plane. In connection to rotation curves, the motion of a particle on a circular orbit and its tangential velocity is of importance. For the static and spherically symmetric spacetime the tangential velocity of a particle in circular orbit corresponds to, where the last equality follows from Eq. (37). The above relation further shows the fact that rotational velocity is determined by g rr component alone. Since v tg is determined by ψ, it is possible to write all the expressions derived earlier in terms of the tangential velocity, e.g.,

From Eq. (42) it is clear that asymptotic limits exist only if
k ∈ (−2B 2 , 2B 2 ) [58]. In this case the solution to Eq. (41) corresponds to, Use of this solution leads to the following asymptotic expression for the tangential velocity, Note that for the following choices, B = 1.00000034 and k = 0.9 the limiting tangential velocity is given by v tg,∞ ∼ 216.3km/s, which is of the same order as the observed galactic rotational velocities. Thus behavior of all the metric coefficients in the solutions depend on two arbitrary constants of integration, namely, k and B. In order to obtain numerical estimate for these parameters we assume that there exist some radius r 0 beyond which baryonic matter density ρ B is negligible. Requiring exp(λ) = 1 − (2GM B /r 0 ), with M B = 4π r0 0 drr 2 ρ B , we readily obtain, Hence the ratio k 2 /B 4 can be determined observationally through the tangential velocity. It follows that around and outside r 0 the radion field will dominate and hence one can introduce a "radion mass" in an identical manner. This on use of conformal symmetry and the ratio k 2 /B 4 from the above equation immediately reads, However the tangential velocity v tg is non-relativistic, i.e., much smaller than unity (in c = 1 units) and hence the radion mass turns out have a scaling relation, The above result explicitly shows that, the "radion mass", will scale linearly with the radial distance, which stops the velocity profile from dieing down at large r. However note that, the linear behavior of radion mass is only the leading order behaviour. If we had kept higher order terms, we would have corrections over and above the linear term, leading to, where C 1 and C 2 are constants depending on κ 2 / and 1 , 2 both are strictly less than unity. Given the mass profile, the corresponding velocity profile can be obtained by dividing the mass profile by r and some suitable numerical factor. The coefficients and powers of the velocity profile (and hence the mass profile) can be determined by fitting the velocity profile with the observed one. We should emphasize that the linear term alone cannot lead to good fit, its effect is to make the velocity profile flat at large distances. Thus at smaller distances the additional correction terms in Eq. (48) are absolutely essential. Hence the effect of radion field can only be felt at large distances, preventing the velocity profile from decaying and the sub-leading factors in Eq. (48) are important for matching with experimental data. In particular, from Fig. 1 it turns out that all the four curves are consistent with the following choices of the power law behavior: 1 0.1 and 2 0.4 respectively. The coefficients C 1 and C 2 turn out to have the following numerical estimates: C 1 = −25.56 ± 4.3 and C 2 = 1.75 ± 0.08 respectively. Thus at small enough values of r the dominant contribution comes from the term C 1 r 1 , while for somewhat larger values of r, C 2 r 2 dominates. Finally at large values of r the linear term, i.e., the contribution from radion field becomes dominating, leading to flat velocity profile for the galaxies. Hence the correction terms are of quite significance in order to obtain a good fit with observational data.

Application to other scenarios
In this work we have used a two brane model with brane separation being represented by the radion field Φ. We have also assumed that our universe corresponds to the visible brane. In such a setup the effective gravitational field equations on the brane, written in a spherically symmetric context, depends on the radion field and its derivatives. Use of collisionless Boltzmann equation leads to the result: virial mass of the galaxy clusters scales linearly with radial distance. Thus without any dark matter we can reproduce the virial mass of galaxy clusters by invoking extra dimensions. However in order to become a realistic model we should apply our results to other situations and look for consistency. There are mainly three issues which we want to address -(i) advantage over other modified gravity models, (ii) reproducing the correct cosmology and (iii) connection with local gravity tests, in particular fifth force proposal. We address all these issues below.
• In present day particle physics an important and long standing problem is the gauge hierarchy problem, which originates due to the large energy separation between the weak scale and Planck scale. In our model the branes are separated by a distance d, such that the energy scale on our universe gets suppressed by M vis ∼ M Pl e −2kd , with k being related to brane tension. Thus a proper choice of k (such that kd ∼ 10 ) leads to M vis ∼ M weak and hence solves the hierarchy problem. Along with the missing mass problem, i.e., producing linear Virial mass our model has the potential of resolving the gauge hierarchy problem as well. This is a major advantage over modified gravity models, where the modifications in gravitational field equations are due to modifying the action for gravity. These models, though can explain the missing mass problem, usually does not address the gauge hierarchy problem.
• The next hurdle comes from local gravity tests. This should place some constraints on the behaviour of the radion field. The analysis using a spherically symmetric metric ansatz has been performed in [61] assuming dark matter to be a perfect fluid which is a perturbation Observed velocity in km/s

Radius in arc second
Velocity Profile for UGC 477 Figure 1: Best fit curves for four chosen low surface brightness galaxies, NGC 959, NGC 7137, UGC 11820, UGC 477, respectively [59,60]. On the vertical axis we have plotted observed velocity in km/s and the horizontal axis illustrates the radius measured in arc second. Good fit shows that the assumption of spherical symmetry is a good one, also the fact that baryonic matter plus radion field explains the galactic rotation curves fairly well. It also depicts the need for the sub-leading terms in Eq. (48), over the Schwarzschild solution. We can repeat the same analysis with our radion field mass function, which is a perturbation over the vacuum Schwarzschild solution. We then can compute the correction to the perihelion precession of mercury due to dark matter which leads to the following constraint on the bulk curvature radius [61][62][63] where ∆δφ = 0.004 ± 0.0006 arc second per century, corresponds to excess in the perihelion precession of Mercury [64]. Among others a is the semi-major axis, e stands for eccentricity, T M and T E are periods of revolution of Mercury and Earth respectively.
• Let us now briefly comment on the relation between existence of fifth force and dark matter. In all these models the generic feature corresponds to the existence of a scalar field which couples to dark matter and in turn couples weakly (or strongly) to standard model particles [65][66][67]. In our model this feature comes quite naturally, since the radion field Φ which plays the role of dark matter can also be thought of as a scalar field, couples to standard model particles through the matter energy momentum tensor with coupling parameter ∼ κ 2 / (3 + 2ω) −1 . Thus effectively we require a fifth force to accommodate modifications of gravity at small scales. There exist stringent constraints on the fifth force from various experimental and observational results (see for example [68][69][70]). We can apply these constraints on the fifth force for scalar tensor theories of gravity and that leads to the following bound on the composite object: (κ 2 /12πG )(1 + Φ) < 2.5 × 10 −5 . Hence for compatibility of the radion field presented in this work with fifth force constraints, the bulk curvature , the bulk gravitational constant κ 2 and radion field must satisfy the above mentioned inequality.
• Finally we address some cosmological implications of our work. In cosmology one averages over all the matter contributions at the scale of galaxy clusters and assumes all the matter components to be perfect fluid. The same applies to our model as well, in which the effect of radion field Φ at the galactic scale is to generate an effective dark matter density profile, with a given mass function. Since the mass function obeys the observed dark matter profile, therefore on the average in the cosmological scale it reproduces the standard dark matter content and hence the standard cosmological models.
Thus the radion field model proposed in this work not only matches with the virial mass profile of galaxy clusters but also fits well into other scenarios. The model has the advantage over other alternative gravity models since it can also address the hierarchy problem by exponentially suppressing Planck scale on our universe. Secondly, local gravity tests and fifth force phenomenology provides constraints on bulk curvature radius which is consistent with virial mass profile. Still, the results can change depending on the stabilization mechanism for radion field, which would be an interesting future avenue to explore.

Discussion
Brane world models can address some of the long standing puzzles in theoretical physics, namely -(a) the hierarchy problem and (b) the cosmological constant problem. To solve the hierarchy problem we need two branes, with warped five dimensional geometry such that energy scale on the visible brane gets suppressed exponentially leading to TeV scale physics. For the cosmological constant brane tension plays a crucial role. Two brane models naturally inherit an additional field, the separation between the branes (known as the radion field). Radion field is also very important in both macroscopic and microscopic physics, for it can have possible signatures in inflationary scenario [25][26][27], black hole physics [71,72], collider searches [73], etc. Along with the gauge hierarchy and cosmological constant problem, another very important problem in physics, is the missing mass problem. This appears since baryonic and virial mass of a galaxy cluster do not coincide. In this work using a two brane setup we have shown that, along with gauge hierarchy and cosmological constant problem, this model is also capable of addressing the missing mass problem through the kinematics of the brane separation, i.e., radion field. Due to the presence of this additional field, the gravitational field equations on the brane gets modified and yields additional correction terms on top of Einstein's field equations. By considering relativistic Boltzmann equation we have derived the virial mass of galaxy clusters, which depends on an effective additional mass constructed out of radion field. Moreover these correction terms modifies the structure of gravity and hence the motion under its influence at large distance, thereby producing a linear increase in the virial mass of the galaxy clusters. This in turn leads to the appropriate velocity law for galaxies within a galaxy cluster, solving the missing mass problem.
We start with various derivatives of the scalar field Φ in this spherically symmetric coordinate: Using which the field equation for Φ takes the following form: Let us now turn our attention to the gravitational field equations. We will start with the temporal component such that: Then the radial component: The above expression on rearrangement leads to Eq. (28). Now we need to solve for Φ, which satisfies the following differential equation, which under integration leads to, Integrating once again we will immediately obtain Eq. (31).