Gravity stabilizes itself

We show that a possible resolution to the stabilization of an extra spatial dimension (radion) can be obtained solely in the context of gravitational dynamics itself without the necessity of introducing any external stabilizing field. In this scenario the stabilized value of the radion field gets determined in terms of the parameters appearing in the higher curvature gravitational action. Furthermore, the mass of the radion field and its coupling to the standard model fields are found to be in the weak scale implying possible signatures in the TeV scale colliders. Some resulting implications are also discussed.


Introduction
Gravity has become the stumbling block in our search for a unified theory, which probably will lead to an understanding of the origin of our universe to the late time cosmic acceleration. On the other hand, even though the standard model of strong and electroweak interactions can explain a vast landscape of experimental results, it continues to have some longstanding unresolved issues, which strongly suggests to look for physics beyond the standard model. One of the major drawback of the standard model is the necessity of fine tuning, which originates from the large hierarchy between the electroweak and the Planck scale, known as the gauge hierarchy problem. It is remarkable that gravity provides a very novel solution to this fine tuning problem through the existence of extra dimensions. Such a gravity based resolution of the gauge hierarchy problem was elegantly described by Randall and Sundrum, where a single extra dimension with manifold structure S 1 /Z 2 was assumed, resulting in two branes (hypersurfaces of (3 + 1) dimensions) located at orbifold fixed points with positive and negative tensions. Subsequently, starting from the Einstein's equations in the bulk (higher dimensional spacetime) with a negative cosmological constant, they could show that physical mass of a field confined on the negative tension brane is in the weak scale, due to an exponential suppression, whose origin traces back to gravity. There have been numerous works later on to clarify some of the disadvantages of this model, as well as in generalizations to more complex settings (for a representative class of works see [1][2][3][4][5][6][7][8][9][10][11][12][13]). One of the key feature of the Randall-Sundrum model is the appearance of an additional four dimensional massless scalar field having no dynamics. This is an undesirable feature, since without a stabilization mechanism one cannot arrive at the desired exponential suppression.
Unfortunately, gravity sector alone could not cure this problem. One had to introduce an additional scalar field in the bulk, whose action when integrated over the extra spatial dimension, provided the potential necessary for stabilization. Any fluctuation about this stabilized value leads to a scalar degree of freedom, called the radion field. There have been numerous studies later on, regarding the details of the stabilization mechanism, corresponding collider signatures and of course, various generalizations, e.g., time dependent stabilization of the radion field (for a small sample of works see [14][15][16][17][18][19][20][21][22][23], and also the references therein).
Even though the above scheme of solving the gauge hierarchy problem looked promising, but there is one aspect, namely the introduction of a bulk scalar which is put in by hand. It would be really intriguing if the whole solution, i.e., the exponential warping as well as the stabilization can come from gravity alone. Further excitement will follow if the radion field so obtained has observable consequences at the collider experiments. In this work, we explore the above possibility and demonstrate explicitly that one can have the desired warping as well as can stabilize the radion field using only gravitational interactions! This is achieved by introducing higher curvature corrections to the Einstein-Hilbert action, which is expected, since the bulk spacetime is governed by Planck scale physics. We further delineate on the phenomenology of the associated radion field, whose potential is being supplied by the higher curvature corrections and demonstrate the significance for collider physics. The phenomenological study enables one to probe the properties of gravitational physics, in particular that of higher curvature gravity, using colliders, leading to new avenues of exploration.
We have organized the paper as follows: We start with a brief introduction to higher curvature gravity and the particular model we will be interested in. Proceeding further we demonstrate how one can have both the exponential warping and radion stabilization in this scenario, the main theme of this work. Finally we discuss the radion phenomenology and comment on possible collider signatures of our model before pointing out future directions of exploration.

Background: Higher curvature gravity
It is generally believed that at high energies (or, small length scales) the Einstein-Hilbert action must be supplemented with higher curvature corrections respecting the diffeomorphism invariance of the action. There are several possibilities for the same, two such candidates being f (R) theories of gravity and Lanczos-Lovelock models of gravity. The Lanczos-Lovelock models are more complicated, but is free of ghosts due to its quasi-linear structure [24][25][26][27][28][29]. On the other hand the f (R) models need special care and must satisfy few conditions in order to ensure its ghost free behavior. The success of f (R) models lie in its excellent match with observations as far as the cosmological arena is considered. Further, f (R) models with a certain constraint on its parameters can also evade the solar system tests as well [30][31][32][33][34][35][36][37][38][39][40][41][42][43][44]. In this work, it will be sufficient for our purpose to focus on the f (R) theories of gravity, satisfying a couple of constraints ensuring its ghost free behavior.
We will work with a five dimensional spacetime consisting of a single extra spacelike coordinate y. The extra dimension will assumed to be compact with a S 1 /Z 2 orbifold structure. Alike the Randall-Sundrum scenario, two branes are located at the orbifold fixed points y = 0, π respectively, with y and −y identified. The bulk gravitational action is assumed to be of the following form where, κ 5 is the five dimensional gravitational constant with mass dimension −3/2, Λ being the negative bulk cosmological constant with mass dimension 5, α and β are the higher curvature couplings having mass dimensions of −2 and −6 respectively. The structure of the above Lagrangian has been inferred from the ghost free criterion, which reads, f ′ (R) > 0, f ′ (R) < 1 and f ′′ (R) > 0 [30]. For convenience we will switch to units where Planck mass has been set to unity, using which the above conditions lead to α > 0 as well as α > |β|. Note that the model must satisfy these criterion at all curvature scales. An exact warped geometric solution to the above gravitational action has been derived recently in [45], which reads Here v is a constant of mass dimension 3/2, we have also defined b 0 = (9κ 2 5 |β|/32α 2 ) and k = −Λκ 2 5 /6 for convenience. For completeness we also present the form of the function f (y) appearing in Eq. (2), having the following structure Note that the solution derived above has the desired exponential warping, through the term kr c y in the warp factor. Such that on the visible brane at y = π, the physical mass of any field will be suppressed by exp(−kr c π), leading to weak scale behavior. This is alike the Randall-Sundrum scenario, where introducing gravity alone leads to the desired weak scale phenomenology with a choice of kr c ≃ 12. We would like to emphasize the fact that the effect of higher curvature gravity is through some combinations of α and β both, keeping only α or β is not sufficient to get the desired warping. The next hurdle is to provide a stabilization mechanism for the radion, which will drive the radion field to its value r c , compatible with the exponential suppression. This is what we will elaborate on in the next section.
3 Stabilizing the radion using higher curvature gravity In the original Randall-Sundrum scenario, one needs to introduce a bulk scalar field in order to achieve the stabilization of the radion field. The reason being, the Lagrangian for the Randall-Sundrum scenario, which is the Ricci scalar, has no potential term for the radion field. The bulk field introduces such a potential, and thereby satisfies both the suppression of the Planck scale and the stabilization of the radion field. We will show that, since our action in Eq. (1) has higher curvature terms, if one evaluates the same for the metric given in Eq. (2) and Eq. (3), it will naturally lead to a potential for the radion field and hence one can stabilize the same without ever introducing a bulk scalar field. The additional degree of freedom originating from the higher curvature terms actually plays the role of a stabilizing field. As a consequence, the stabilized value would depend upon the parameters α and β appearing as the couplings of the higher curvature terms. It will turn out that the condition kr c ≃ 12 is essentially a condition on these higher curvature couplings and the bulk cosmological constant. As we have already laid out the principles involved, we will proceed directly to the computation and will work exclusively with the higher curvature action presented in Eq. (1). Given the metric ansatz in Eq. (2) one can evaluate the Ricci scalar in a straightforward manner. Then one has to substitute the Ricci scalar in the gravitational part of the action and obtain the corresponding f (R) Lagrangian, which to leading orders in the coupling parameters reads (for details see Appendix A), Given this form of the f (R) Lagrangian, one can substitute the same in the bulk action, i.e., Eq. (1) and then integrate out the extra dimensional coordinate y over the interval [0, π], thanks to the orbifold symmetry. In this integration it will turn out that the contribution coming from the lower limit y = 0 is independent of the radion field r c and hence adds a constant contribution to the radion potential, while the contribution from y = π does have dependence on the radion field and shall serve as the radion potential. Introducing a new field φ(r c ) = Φ exp[−kr c π], where Φ = 24/kκ 2 5 , we finally obtain the potential for the radion field (or, equivalently for the new field φ) to yield, where, δ = (2b 0 /kκ 2 5 ) = (9 |β|/16kα 2 ), is a dimensionless constant. The usefulness of this quantity φ follows from the fact that one can upgrade it immediately to the status of a four dimensional scalar field with r c → r(x), spacetime dependent brane separation. Such that the vacuum expectation value of the same is given by r c and thus φ(r c ) will denote the vacuum expectation value of φ [14,15]. Note that the field φ is merely a constant depending on the radion field. If the field is being upgraded to depend on the spacetime coordinates, then the potential structure will remain identical, however it will inherit the canonical kinetic term in the action as well (see, for example [15]).
Note that the above potential is very much similar to the one obtained in [14,15] using bulk scalar field, with δ identified as (m 2 /4k 2 ), where m is the mass of the bulk scalar. Thus the higher curvature terms act as a source for the radion mass as we will explicitly illustrate later. Further in the above scheme the condition δ < 1 is identically satisfied, since the couplings to higher and higher curvature terms are more and more suppressed. Thus the scenario with higher curvature gravity leads to an identical situation as that of introducing a bulk scalar field, but follows from the gravity sector alone. The stabilized value of the radion field can be obtained by finding out the minima of the potential in Eq. (5), which is a solution of the equation ∂V /∂φ = 0, leading to the following expression for r c (see Eq. (20) in Appendix A), Thus with κ 5 v ∼ 40 and |β|/kα 2 ≃ 1/20, the logarithmic term becomes of order unity and then one readily obtains kr c ≃ 12, the value desired for exponential warping. Hence starting from a pure gravitational action, with higher curvature corrections, one can produce an exponential warping as well as a proper stabilized value for the brane separation without ever introducing any additional structure. Further the desired warping to address the hierarchy problem leads to a relation between the higher curvature couplings and the bulk cosmological constant. This completes what we set out to prove, i.e., derivation of the exponential warping leading to weak scale physics and a proper stabilization mechanism for the brane separation, both from the gravitational dynamics alone. For the sake of completeness we would like to discuss on the choice of f (R) gravity and its role in radion stabilization. Since we are working in the higher curvature regime, where bulk effects are important, it makes sense to add terms like R n to the Einstein-Hilbert action, with n positive. The first such choice corresponds to adding a R 2 term, which would lead to the original Goldberger-Wise scenario in the scalar tensor representation. The next leading order term in a bulk spacetime with negative cosmological constant, free of ghosts correspond to the one presented in Eq. (1). Interestingly for this situation one can solve for the full scalar coupled Einstein's equations as depicted in [45] and the situation will differ considerably from the Goldberger-Wise scenario. In principle one can add more higher order terms to the Lagrangian, however in those scenarios one would not be able to solve the full back-reacted problem in scalar-tensor representation. Moreover, these terms will be further suppressed and will contribute insignificant corrections over and above Eq. (6). Thus the scenario presented in this work captures all the essential features and is simple enough to be solved in an exact manner. This motivates the choice presented in Eq. (1). Given the above, it will be worthwhile to spend some time discussing the corresponding scenario in the scalar-tensor representation [45], which will bring out the difference of our approach with the existing ones. This is what we elaborate in the next section.

Stabilization in the Einstein frame: Scalar-tensor representation
It is well known that any f (R) gravity model is mathematically equivalent to a dual scalar-tensor representation [21,[36][37][38][46][47][48][49][50][51][52]. The mathematical equivalence follows from the transformation of the Jordan frame action to the Einstein frame aka conformal transformation. Surprisingly, this equivalence holds in lower dimensions as well, viz., if one starts from a higher dimensional action and projects on to a lower dimensional hypersurface the field equations are still connected by conformal transformation, provided one exercise caution about the boundary contributions. Despite the mathematical equivalence, there are situations where the two scenarios are not physically equivalent, e.g., in cosmological scenarios the f (R) frame may lead to late time acceleration, while the Einstein frame advocates late time deceleration [48,53,54]. The issue of physical inequivalence becomes important if the spacetime inherits a singularity or is undergoing a quick evolution phase. In particular, using reconstruction technique [55][56][57][58] it is possible to generate f (R) theories explaining the early inflationary phase to late time accelerating phase of the universe, all of which ultimately results into finite time future singularity. The existence of future singularity often breaks the equivalence with scalar-tensor representation [48,59]. However in absence of singularity [60] equivalence of f (R) theories with scalar tensor theories does exist. The situation discussed in this work has no such singularity in the spacetime structure, as evident from Eq. (2) and Eq. (3) respectively. Thus one may safely use the equivalence between f (R) and scalar tensor theories. Given this input, it will be worthwhile to explore the corresponding situation in this context, namely how the radion stabilization is affected as one considers the dual picture in the Einstein frame and contrast the same with the stabilization already discussed in Section 3. For the f (R) action under our consideration, the corresponding action in the Einstein frame becomes, where ψ is the scalar field in the dual picture defined as, κ 5 ψ = (2/ √ 3) ln(1 + f ′ ) and the corresponding potential becomes The expression for the potential brings out the key difference between the Goldberger-Wise stabilization mechanism and the one advocated here -the potential for the scalar field in Goldberger-Wise mechanism lacks the quartic term present in our analysis. Incidentally, the presence of this quartic term helps to solve the full back-reacted problem, while the original stabilization proposal was without incorporating the backreaction of the scalar field. Hence the stabilized value of the radion field derived above incorporates the effect of the quartic potential, as well as the back-reaction of the scalar field on the spacetime geometry and differs from the standard scenarios. The above structure of the potential also shows the reason for neglecting further higher curvature terms (e.g., R 6 ) in the action 1 . The scalar tensor representation with such higher curvature terms will involve more complicated potentials and hence cannot be solved in full generality by incorporating the back-reaction as well. Further the R 4 term in the action leads to the leading order departure from the Goldberger-Wise action, which we have studied in this paper. Additional higher curvature terms would lead to further sub-leading corrections and thus can be neglected. Given the action in the Einstein frame, one can invoke the corresponding solution (derived in Appendix A) and integrate out the extra dimensional part present in the action. This will result in a potential for the radion field r c , alike the Goldberger-Wise mechanism, whose minima would yield the following stabilized value of the radion field, where ψ 0 is the value of the dual scalar field in the y = 0 brane. Comparison with Eq. (6) reveals that the leading order behavior (i.e., the term outside logarithm) of the stabilized value of the radion field is identical in both Jordan and Einstein frame. This observation explicitly demonstrates that, in both these frames the radion is stabilized to the desired value necessary for exponential suppression of the Planck scale. Thus unlike various scenarios with either singularity or a quick evolution (where the two frames are physically inequivalent) in this particular situation the physical equivalence between the two frames is manifest.
At this stage, we may point out another alternative possibility of stabilizing the radion field, by incorporating quantum effects of the bulk scalar field at nonzero temperature [61,62]. In particular, one can think of the resulting thermal fluctuations to generate a modulus potential which may inherit a minimum, thereby stabilizing the brane separation. The brane separation necessary to solve the hierarchy problem involves considering a low temperature limit of the free energy associated with the bulk quantum field, which may have connections with the AdS/CFT correspondence [61,62]. However, in this context as well, one generally neglects the effect of back-reaction and treats the bulk quantum field to be siting on the fixed Anti-de Sitter background, unlike the scenario we have depicted. Thus after explicitly establishing the differences between our approach and the existing ones, we now concentrate on the phenomenology of the radion field, viz. mass of the radion field and its interaction with the standard model fields in the next section.

Some applications of the stabilized radion field
In this section we will briefly discuss two possible applications of the stabilized radion field in two diverse physical contexts. The first application will discuss the implications of this stabilized radion from the perspective of particle physics, while the other will address the effect of the radion field on the inflationary paradigm. In both these contexts we will demonstrate that the stabilization of the radion field (which has its origin in higher curvature gravity) plays a crucial role and leads to interesting additional structure in the corresponding situations having possible observational consequences.

Phenomenology of the stabilized radion
Given the potential V (φ), one can immediately obtain the mass of the φ excitation by computing ∂ 2 V /∂φ 2 and then expressing the same using r c , given in Eq. (6). Performing the same, one arrives at the following expression for radion mass, where, φ(r c ) can be obtained from Eq. (6). It is evident from the expression for m φ , that the there is an exponential suppression of the radion mass, which leads to a weak scale value from a Planck scale quantity. Note that this expression is very much similar to the result obtained in [14,15], with the identification of δ with m 2 /k 2 . But with one difference, which is caused by the (φ(r c )/Φ) 2δ term in Eq. (10). This results in a decrease in the radion mass as compared to the Goldberger-Wise scenario described in [15]. However at the same time the choice of κ 5 v also becomes important. To see the difference in a quantitative manner, consider the following situation: κ 5 v ∼ 40, δ ∼ 1/32, using which one obtains kr c ∼ 12, leading to, m φ ∼ 0.02k 2 e −2krcπ . While for the standard Goldberger-Wise scenario one would have obtained m 2 φ ∼ 0.05k 2 e −2krcπ . This explicitly shows that the radion mass in our approach is lighter (in this case two times) compared to the one obtained in [15], depicting quantitatively the difference between these two approaches. Further note that the results in [14,15] were derived in the context of a bulk scalar field, here we derive the same but from a purely gravitational standpoint. Similar exponential suppression will drive the masses of the low lying Kaluza-Klein excitations to TeV scale [63][64][65][66][67][68][69][70][71][72][73][74][75][76]. In the context of radion mass however there is one further suppressing factor, namely δ 2 and thus radion mass will be a bit smaller compared to the low lying Kaluza-Klein excitations of the bulk field. Since the radion mass, as presented in this work is completely of higher curvature origin, this suggests that detection of the radion field φ may be a hint not only of higher dimensions but also of higher curvature gravity.
In order to find out coupling of the radion field to standard model particles, note that in this scenario radion appears as a gravitational degree of freedom and since gravity couples to all kinds of matter, so will the radion. Since the standard model fields are confined on the brane hypersurfaces, it is clear that they will couple to the induced metric on the brane. On the brane located at y = 0, the warp factor becomes unity and the induced metric is proportional to η µν . Thus fields confined to y = 0 hypersurface does not couple to the radion field. On the other hand, on the visible brane (the hypersurface y = π) the induced metric is (φ/Φ) 2 η µν and thus radion field will couple to the standard model fields. Here, φ = φ(r c ) + δφ, where δφ is the fluctuations around the minimum value φ(r c ). In the case of a scalar field, say the standard model Higgs h(x), one can reabsorb the factors of Φ by redefining h → (Φ/φ(r c ))h(x), such that physical Higgs mass becomes m 0 e −krcπ . For a Planck scale bare mass m 0 and kr c ≃ 12 one immediately obtains the physical mass in the weak scale. Further the corresponding interaction term of a standard model field with the radion will involve where T µ µ stands for the trace of the energy momentum tensor of the standard model field and Λ φ defines the coupling of the radion to the standard model fields. (Since radion is a gravitational degree of freedom it has to couple to some combination of the matter energy momentum tensor). Note that for large Λ Φ , the coupling becomes small and as a consequence radion will couple weakly to the standard model fields.
In particular recent bounds on both radion mass and radion coupling strength shows that they are not independent, if radion mass is smaller the coupling will be larger and vice versa. To get a numerical estimate, one must provide an estimate for the value of k in Planck units. For k ∼ 0.1, the coupling satisfies the following stringent bound Λ φ > 14.3 TeV [76][77][78][79][80][81]. This leads to very weak coupling between radion and standard model fields, resulting in non-detectability of the radion field (even though the mass can be as low as 200GeV). At the same time for k ∼ 1, in Planck units, the coupling can become Λ φ ∼ 5 TeV. Even though in this case the coupling is not very weak, but the radion mass becomes higher ∼ 1TeV. Thus coupling and mass of the radion field are inversely related, making it difficult to detect in the current generation collider experiments.
For completeness, let us discuss the implications of our model in the context of electroweak precession measurements [82]. As emphasized earlier, the radion will linearly couple to standard model fields, e.g., the gauge bosons W ± and Z with a coupling Λ −1 φ , which is (TeV) −1 . These interaction terms will produce three point as well as four point interactions between radion and the gauge bosons due to gauge fixing terms as well as higher loop effects in the effective Lagrangian. Fixing the mass of the standard model Higgs boson at 125 GeV [77], affects the electroweak parameters such that they become insensitive to the radion mass, since the radion mass is further suppressed by its own vacuum expectation value [82]. Despite above, the fact that the mass of the radion field as well as its coupling with standard model fields are in the TeV scale makes the phenomenology of the radion field a nice testbed for higher dimensional as well as higher curvature physics in the next generation colliders.

Cosmology and radion stabilization
Having already discussed the imprints of the stabilization mechanism of the radion field on the phenomenological side, we will presently address the corresponding situation in a cosmological setting. In a more general context, following [65] one should have made the radion field dynamical by considering arbitrary fluctuations around the stabilized value and hence study the dynamics in an arbitrary background. However this will be a complicated exercise, due to presence of higher curvature terms in the gravitational Lagrangian. Thus we concentrate on a specific situation with the radion field depending on time alone, leading to the stabilized value in the cosmological context. Note that this situation has already been analyzed in the context of general relativity in [83] and has been elaborated further in the context of radion stabilization in [23]. Thus one may try to find out the cosmology on the brane in presence of higher curvature terms, whose detailed analysis for various cosmological epochs will be presented elsewhere. In this section we will try to answer this question in the context of inflationary paradigm alone. However for completeness, we will also provide some general computations. The metric ansatz suitable for our purpose corresponds to, where, both the branes are located at y = 0 and π respectively and they are assumed to be expanding with the scale factor a(t). Further, the radion field is assumed to be dynamical, such that, r(t) = r c + δr(t).
Here r c corresponds to the stabilized value of the radion field derived in the earlier sections and δr(t) corresponds to the fluctuations around the stabilized value. Further A(y) should take care of all the extra dimensional dependent quantities, which will behave as ky to the leading order. The Ricci scalar derived for the above ansatz, on a y = constant hypersurface becomes, Having obtained the Ricci scalar one can derive the Lagrangian with ease which is given in Eq. (1). Derivation of the Lagrangian enables one to obtain the corresponding field equations for a and r, by varying the scale factor and the radion field respectively. The gravity theory being f (R) in nature, the field equations will definitely inherit higher than second derivatives of the scale factor and the radion field, which reads, In general it looks sufficiently complicated, however in the case of an exponential expansion (with a(t) = e Ht ), i.e., for inflationary scenario the above equation simplifies a lot. In particular it is possible to approximately solve for the time dependence of the radion field explicitly, which turns out to be decreasing with time, similar with the corresponding situation in general relativity [23]. Thus as the universe expands exponentially, the radion field decreases with time, finally attaining the stabilized value as the inflation ends. Hence the scenario depicted above can also explain a dynamical procedure for stabilization of the radion field, modulo inflationary paradigm. This provides yet another application of the radion field aka higher curvature gravity in the present context.

Discussion
Gauge hierarchy problem is a very serious fine tuning problem in standard model physics. One avatar of the extra dimensional physics as depicted by Randall and Sundrum has the capability of addressing the gauge hierarchy problem by suppressing the Planck scale to weak scale gravitationally. Unfortunately, to stabilize the above scenario one needs to introduce an additional field. In this work, we have shown that the introduction of such a stabilizing field is unnecessary and one can achieve both the suppression to weak scale as well as a stabilization mechanism starting from a higher curvature gravitational action alone. In this sense gravity stabilizes itself! Starting from a higher curvature gravitational action we have derived an exact warped geometric solution with the extra spatial dimension having S 1 /Z 2 orbifold symmetry and two 3-branes located at orbifold fixed points y = 0, π respectively. We would like to emphasize that the above solution is exact with effects from the higher curvature terms duly accounted for, unlike the original Goldberger-Wise solution where back reaction of the stabilizing field was neglected. In this respect our work is more in tune with [63], where also solutions have been derived with inclusion of back-reaction as well. The warp factor again has exponential suppression and thus the Planck scale physics will be reduced to weak scale phenomenon on the visible brane located at y = π. But, surprisingly the higher curvature terms provide a potential for the radion field as well, whose minima leads to the stabilized value of the radion field gravitationally. Further, the stabilized value depends on the gravitational couplings present in the higher curvature action as well as the bulk cosmological constant, such that a particular choice of these parameters results to kr c ≃ 12, leading to the desired warping. Added to the excitement is the result that the radion mass has similar suppression leading to TeV scale physics, which is smaller compared to the low-lying Kaluza-Klein spectrum of the bulk scalar field. Being a gravitational degree of freedom, radion field couples to the standard model fields through the trace of the energy momentum tensor with a coupling having strength TeV −1 . Once again, we reiterate on the fact that the above analysis has been performed completely in a gravitational physics framework.
The above results open up a broad spectrum of further avenues to explore. From the observational point of view, the above result brings down the dynamics of the higher curvature gravity to a TeV scale phenomenon, which may become accessible in near future collider. In particular the radion mass which depends exclusively on the couplings present in the higher curvature theory may provide the first hint towards observational signatures of higher curvature gravity besides that of higher dimensions. A more careful analysis in this direction can be performed following [67][68][69][70][71], where similar analysis for an additional stabilizing scalar field has been carried out. From the theoretical hindsight, it will be worthwhile to understand the phenomenon of generating a potential for the radion field and its subsequent stabilization in the context of Lanczos-Lovelock (or, Einstein-Gauss-Bonnet for simplicity) gravity. It will also be of interest to explore the consequences of making the radion field dynamical, in which case the kinetic term will also contribute to the gravitational field equations resulting in distinctive cosmological consequences of this higher curvature brane world scenario, e.g., imprints on the Cosmic Microwave Background, inflationary scenario and so on.

Acknowledgement
Research of S.C. is supported by the SERB-NPDF grant (PDF/2016/001589) from Government of India.

A Appendix: Calculational Details
In this section, we provide all the calculational details to supplement the results presented in the main text. In order to make the calculation simple, we will introduce the following definitions: We further divide the appendix in two parts, the first done depicts the situation with f (R) gravity, while the other illustrates the dual scalar-tensor description.

A.1 Stabilization in f(R) gravity
The Ricci scalar computed for the metric presented in Eq.
(2) turns out to be where 'prime' denotes differentiation with respect to the extra spatial coordinate y. Evaluation of the Ricci scalar, given the functions A(y) and f (y) can be performed to leading orders of exp(−A 2 r c y) yielding Let us now write down the higher curvature Lagrangian of Eq. (1) using the above expression for the Ricci scalar, Thus the bulk gravitational action reads, Integrating out the extra spatial coordinate, leads to a potential term for the radion field, which contributes only at y = π and becomes, Such that the minima of the potential, leading to stabilized value of r c can be obtained by equating ∂V /∂r c = 0 resulting in The above expression can be enormously simplified by keeping terms to the leading order. In particular, one will have P = −20k 2 , Q = −(40/3)k 2 A 3 and R = 80kA 1 A 2 leading to Eq. (5), where we have assumed A 1 > k > A 2 and A 1 > 1. It turns out that in order to have the stabilized value of kr c ≃ 12, the above conditions are identically satisfied. Introduce the field φ, such that the above potential becomes, This is the expression, which along with the previous identifications has been used to arrive at Eq. (6). Given the potential V (φ), one can compute, as well as, Thus the minimum of the potential corresponds to ∂V /∂φ = 0, Thus the radion mass becomes, This helps to obtain the radion mass as in Eq. (10).

A.2 Stabilization in scalar-tensor representation
Since any f (R) theory of gravity has a dual description it would be interesting to understand the situation presented above in the dual scalar-tensor theory as well. This will provide a similar setup to the original Goldberger-Wise scheme but with back-reaction included. In this case the gravitational action is given by Eq. (7), with the following solution, ds 2 = e −2A(y) η µν dx µ dx ν + r 2 c dy 2 ; A(y) = A 0 + kr c y + κ Note that the correspondence between Jordan and Einstein frame relates the parameter A 2 appearing here with that earlier. The Ricci scalar for the above metric can be evaluated leading to, where, a 0 is a constant to be identified with 3k, while, b 0 = κ 2 5 A 2 /4. Further one can obtain, Such that the action becomes, In this case as well the minima of the potential originates from the following algebraic equation, Leading to a solution for the stabilized potential, which only differs within the logarithm and hence have identical leading order behavior. The mass as well as the coupling with standard model fields work out in an identical fashion. Thus the two frames lead to very much similar expressions for the stabilized value of the radion field as well as its mass and couplings with standard model fields. Hence at least in this context the two frames have nearly identical physical behavior.