Nested warped geometry in a non-flat braneworld scenario

We generalize nested multiply warped braneworld models by incorporating non-zero brane curvature caused by an effective cosmological constant Ω\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Omega $$\end{document} induced on the 3-branes. Starting with the doubly warped model, we first analyze the case where the maximally warped brane is identified as the visible brane. For Ω<0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Omega <0$$\end{document}, resolution of the gauge hierarchy problem imposes a small upper bound on |Ω|\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|\Omega |$$\end{document}, and can possibly lead to positivity of all the 3-brane tensions. For Ω>0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Omega >0$$\end{document}, the latter is not possible but the tuning of the cosmological constant to its tiny observed value is linked to the tuning of the extra dimensional moduli close to the inverse Planck length, justifying the original flat brane approximation. In both regimes, we study the dependence of the scale-clustering of the pair of TeV-branes on the brane cosmological constant and hence its potential role in generating a fermion mass hierarchy between these branes. Identifying the near-maximally warped brane as the visible brane instead opens up regions in the parameter space that allow positive 3-brane tensions for both anti de Sitter and de Sitter branes, subject to non-trivial constraints on the warping parameters. We conclude by generalizing the key results to arbitrary n-fold nested warped braneworlds.


Introduction
The five dimensional warped braneworld model [1] proposed by Randall and Sundrum provides, at first glance, an attractive theoretical framework within which the gauge hierarchy problem of the Standard Model can be naturally resolved without extreme fine tuning of parameters.Equipped with a modulus stabilization mechanism which is either based on BSM fields introduced in the bulk [2][3][4][5][6][7] or purely gravitational in origin [8][9][10][11], the low-energy phenomenology of the model predicts a scalar radion typically with TeV-scale mass.The lowest-lying KK-mode of the five-dimensional (5D) graviton is expected to be similarly massive, with enhanced coupling (compared to the massless zero-mode) to SM fields on the visible brane due to the same warping which generates the large mass hierarchy [12][13][14][15][16][17][18][19][20][21].The fact that these theoretical predictions fell within experimentally testable domains accessible to present generation colliders largely motivated the shared interest of both communities in the RS1 model following its conception.However, the non-detection of graviton KK-modes through a variety of channels at the LHC till date [34][35][36][37][38][39] has increasingly constrained the parameter space of the five dimensional paradigm.In particular, a small hierarchy of magnitude m H /M 5 ∼ 10 −2 between the fundamental Higgs mass (m H ) and the 5D Planck scale (M 5 ) is currently required to explain the null results.This so-called "little hierarchy" implies possible emergence of new physics nearly two orders of magnitude below M 5 .This is particularly disconcerting since it requires the orbifold radius r c to be enlarged by a factor of O( 102 ).As the exponential dependence of the warp factor on r c makes the former sensitive to minor changes in the latter, this possibility becomes severely restricted, which in turn diminishes the efficacy of the 5D model in its originally proposed form.Motivated largely by the aforementioned reasons, several higher dimensional generalizations of the original RS model have been proposed, most of which rely upon additional orbifolds of the form S 1 /Z 2 [40][41][42][43][44][45][46][47][48][49].In this regard, the possibility of a nested doubly warped flat braneworld model [50] is particularly interesting, where the singly warped five dimensional metric from [1] is further warped along the direction of a second orbifold, thereby resulting in a "brane-box" configuration with topology M(1, 3) × S 1 /Z 2 × S 1 /Z 2 .The bulk is a slice of AdS 6 and the "walls" are formed by four 4-branes, which intersect at the "vertices" of the box and give rise to four 3-branes warped to different extents, hence having distinct physical mass scales.In absence of any large hierarchy between the orbifold radii, simultaneous large warping along both directions is forbidden, which causes the four mass scales to be clustered in a near-TeV scale pair and a near-Planck scale pair.Such a setup has certain notable advantages over the 5D scenario.Firstly, the doubly warped structure of the metric renders the first KK mode of the 6D graviton heavier than that of the 5D scenario, while further suppressing its coupling to the SM fields confined to the visible brane [51].Unlike the 5D case, these features together help the doubly warped model survive current collider constraints without necessitating any little hierarchy between m H and M 6 .At the same time, significant regions of the parameter space for the doubly warped model are accessible to the LHC and its near-future upgrades and can be probed in upcoming runs, thus marking it as a model of direct experimental interest [52].Secondly, it offers a natural explanation for the mass hierarchy observed among SM fermions if the latter are described by five dimensional fields which are allowed to extend into the bulk.Following standard KK decomposition of such a fermionic wavefunction, it can be shown that boundary kinetic terms localized at the maximally and near-maximally warped 3-branes can alter the effective 4D fermion-scalar Yukawa couplings, thus leading to a splitting among the fermion masses [50,53].This splitting is quite small as the physical mass scales of this pair of 3-branes are clustered around O(TeV).Other phenomenological aspects of multiply warped flat braneworld models, e.g.moduli stabilization issues and dynamics of bulk matter and gauge fields, have been studied in [54][55][56][57][58][59], shedding light on a host of interesting properties.In a recent work [60], it has also been shown how the mass discrepancy of the W -boson with the SM prediction as observed by the CDF Collaboration [61] can be explained in a doubly warped background.
On the other hand, the negative tension of the TeV-brane has long been an unattractive feature of the 5D RS model, as such branes are known to suffer generically from instability issues in both classical Einstein gravity and beyond [62][63][64][65].Although there exist some dynamical schemes [66,67] which make the effective 4D theory on the negative tension brane consistent with general relativity in conjunction with a separate modulus stabilization mechanism (as in [2]), it would be interesting if the TeV-brane could somehow be imparted positive tension instead and the entire issue avoided.The latter is particularly desirable if one considers string motivated braneworld scenarios, where the corresponding stable D-branes turn out to have positive tensions.Interestingly, this can be achieved for a minimalistic generalization of [1] to the case of curved 3-branes as analyzed in [68].The origin of non-zero brane curvature can be traced to an effective 4D cosmological constant (Ω) induced on the 3-branes, leading to either an AdS-Schwarzschild (Ω < 0) or a dS-Schwarzschild (Ω > 0) geometry.For Ω < 0, the tension of the TeV-brane can be rendered positive.Additionally, one obtains a tiny upper bound (∼ 10 −32 ) on the admissible magnitude of |Ω|r 2 c .Besides successfully averting the issue of the negative tension visible brane, the model is thus interesting from a cosmological viewpoint as well.On the other hand, Ω > 0 permits a metastable minimum in the modulus potential whose form is determined solely by the brane vacuum energy, hence obviating the need to invoke any external bulk scalar or to modify the gravitational sector for the purpose of modulus stabilization [69].This possibility is arguably more attractive than the conventional bulk field prescription generalized to the curved scenario [70].Some other non-trivial aspects of the 5D curved braneworld, e.g.dynamics of bulk fields and radion-driven inflation, have been studied in [71][72][73].
Motivated by the features of the 5D non-flat and the 6D flat cases separately, it is natural to take the next logical step -generalizing the geometric structure of higher dimensional mul-tiply warped spacetimes (i.e.6D and beyond) by incorporating non-zero brane curvature.As we show in this work, such models have a variety of interesting properties which set them apart from both individual progenitor models.

The 6d Einstein equations
Assuming matter fields to be absent, the total gravitational action comprised of the bulk (S 6 ) and 4-brane (S 5 ) contributions is S = S 6 +S 5 , where S 6 is the Einstein-Hilbert action in presence of the bulk cosmological constant Λ 6 , and S 5 contains the brane-localized terms arising out of intrinsic vacuum energy densities of the 4-branes, i.e., the brane tensions.
where R is the six-dimensional Ricci scalar, M is the fundamental (six-dimensional) Planck mass, and g 5 and ḡ5 are the induced metrics on the corresponding 4-branes.To solve the resulting Einstein equations, we first assume a nested metric ansatz of the following form.
where R y and r z are the (constant) radii of the orbifolds.This ansatz differs from the familiar flat brane case in the consideration of an arbitrary g µν (instead of η µν ) on the 3-branes.Plugging this ansatz into the total action and extremizing it leads to the following set of Einstein equations. µν-component: yy-component: zz-component: where G µν is the four dimensional Einstein tensor, R = g µν R µν is the four dimensional Ricci scalar, and the primes and dots denote derivatives wrt y and z respectively.
In the bulk away from the turning points, the δ-functions can be safely ignored and the brane tension terms no longer contribute.Dividing both sides of Eq. (2.4) by g µν for any pair of µ and ν and rearranging terms, it becomes clear that one side of Eq. (2.4) depends on x µ alone, while the other side is a function only of the compact coordinates y and z.This allows us to separate the variables as where Ω is a global constant, which, as evident from Eq. (2.7), plays the role of an induced four-dimensional cosmological constant.This equation further implies a constant 4D curvature R = 4Ω, which can be used to eliminate R from Eq. (2.5) and Eq.(2.6) in terms of Ω. Plugging this result into Eq.(2.5) and using Eq.(2.8) to simplify, we obtain Putting this back into Eq.(2.6) and invoking Eq. (2.5) to separate the variables y and z finally gives us the following uncoupled ODEs in a(y) and b(z): where α is a dimensionless constant of separation.Solving these ODEs allows us to determine the two warp factors explicitly, followed by calculations of the brane tensions by incorporating the boundary terms in Eq. (2.5) and Eq.(2.6).But the natures of the a(y) and b(z) solutions depend crucially on the positivity (dS-Schwarzschild 3-branes)/negativity (AdS-Schwarzschild 3-branes) of Ω.In the following sections, we focus on the two possibilities separately.
3 AdS-Schwarzschild 3-branes (Ω < 0) In this regime, the general solutions of Eq. (2.10) and Eq.(2.11) are given by As in the flat brane case, we must choose an appropriate normalization for the warp factors such that their values are bounded by (0, 1].It is straightforward to see that a(0) = b(π) = 1 is the most natural choice.For b(π) = 1, one obtains ω1 = 1, and c1 becomes a superfluous constant which can be set to unity to simplify Eq. (3.2) to the following form.
Apparently, the solution of b(z) is identical to the solution in the flat brane limit (i.e., with Ω → 0), and the relation between the moduli α and β is exactly identical to the equality which relates their flat limit counterparts c and k.In hindsight, this makes sense because of the structure of the metric in Eq. (2.3).Due to nested warping, the information about the curvature of the 3-branes is encapsulated principally in a(y), and affects b(z) implicitly through the relation between α and β.Next, one can evaluate c 1 from a(0) = 1 as where we have taken only the positive discriminant solution, as in the limit ω 1 → 0 one needs to ensure c 1 → 2 in order to reduce Eq.(3.1) to the well-known flat result a(y) → e −αy .The admissible range of As the final step, the solutions a(y) and b(z) need to be plugged in Eq. (2.8) so that β can be obtained in terms of the fundamental parameters: which, like Eq. (3.3), expectedly mimicks the flat brane result and requires Λ 6 < 0 to be physically meaningful, i.e., the bulk should be AdS 6 .
We can now substitute the warp factors directly in Eq. (2.4), use Eq.(2.7) and Eq.(2.8) to eliminate the bulk part, and integrate over infinitesimal ε-intervals across each boundary point to obtain the corresponding brane tension.Owing to Z 2 symmetry of the orbifolds, one must take care to replace y with |y| in the expression of a(y) once the boundaries are taken into account.
V 3 (y) = 0 (3.8) where we have exploited the normalization a(0) = 1 (alongside the fact ω 1 < c 1 ) and the α − β relation from Eq. (3.3) to simplify.The explicitly coordinate-dependent term of each brane tension is essentially the same as that of the flat case.Only the constant coefficients are modified due to the presence of non-zero brane curvature.It is obvious that in the limit ω 1 → 0, the entire set Eq. (3.6)−Eq.(3.9) reduces to the results derived in [50] for the flat doubly warped model.
Owing to the relation between α and β in Eq. (3.3), it can be argued that α and β cannot both be large if, based on naturalness, one assumes R y ∼ r z , i.e., there is no considerable hierarchy between the two extra dimensional moduli.Instead, one must have either α > β or α << β, which respectively cause warping predominantly along either y or z direction.The maximally warped 3-brane is the (y = π, z = 0) brane, which we identify as the visible brane for now, based on the conservative assumption that no other brane has a physical mass scale lower than that of the visible brane.The total warp factor is given by a(π)b(0), which can be factorized as where x 1 = απ and x 2 = βπ, and n 1 and n 2 are positive constants quantifying the extents of warping along the y and z direction respectively.For successful resolution of the gauge hierarchy problem, one requires n 1 + n 2 16.We show that in absence of a large hierarchy between R y and r z , n 1 and n 2 cannot be of similar magnitudes.The exact solutions of Eq. (3.10) are from which one immediately obtains the constraint ω 2 1 ≤ 10 −2n1 .Now, if possible, let n 1 and n 2 both be large, e.g.n 1 ∼ n 2 ∼ 8. Then ω 2  1 is extremely small and renders c 1 ≈ 2. We can write ω 2 1 = 10 −N , with the minimum allowed value of N being N min = 2n 1 .For N → ∞, the solution of Eq. (3.11) reduces to the familiar flat result x 1 = n 1 ln10, whereas for N = N min = 2n 1 we have x 1 ≈ n 1 ln10 + ln2, i.e., the flat result plus a minor correction due to the brane curvature.However, for N >> N min , one obtains two distinct values of x 1 corresponding to the two roots of Eq. (3.11), both of which are physically valid solutions ≈ (N − n 1 )ln10 + ln4 (3.13)where x (1) 1 is practically the flat result, and x 1 is slightly larger than x 1 .Hence, in presence of brane curvature, a given degree of large warping along y direction can owe its origin to two distinct values of α.Being of similar magnitudes, these two values do not give rise to any new hierarchy that we need to worry about.On the other hand, Eq. (3.12) admits only one solution for large n 2 : For close values of n 1 and n 2 (e.g.n 1 ∼ n 2 ∼ 8), it is clear that α and β must be of similar magnitudes (e.g.α ∼ β ∼ 6), which leads back to the large R y /r z hierarchy that we are trying to avoid.The only way out is to consider dominant warping along either y or z, corresponding to either α > β or α << β as argued in the following sections.

Dominant warping along y (α > β)
It suffices to consider α ∼ 10 and β ∼ 0.1 to avoid a significantly large R y /r z ratio.These values correspond roughly to n 1 ∼ 15 and n 2 ∼ 1.The largeness of n 1 makes the two solutions of α given in Eq. (3.13) valid in this regime.As n 2 is not quite as large, the approximate solution from Eq. (3.14) does not hold as good, and one needs to consider the exact version in Eq. (3.12) instead for a precise value of β corresponding to a given n 2 .Thus, for Ω < 0 and dominant warping along the direction of the first orbifold, the gauge hierarchy problem can be resolved for two distinct values of the orbifold modulus, one of which (x In Figure 1a, using the equation a(π) = 10 −n1 from Eq. (3.10), we have plotted N as a function of x 1 for a few representative values of n 1 which result in no large radius hierarchy.The analytical form of N (x 1 ) near N min corresponds to the first solution x (1) 1 from Eq. (3.13), which is the dominant solution near the minimum.For each value of n 1 , N diverges and we obtain ω 2 1 = 0 at the limiting value of x 1 = n 1 ln10, and a global minimum exists at N min = 2n 1 .The linear relation between N min and n 1 is also evident from the set of curves shown in the plots.
Figure 1b shows the behaviour far from the minimum, i.e., for N >> N min , where the second solution x (2) 1 from Eq. (3.13) dominates and makes N (x 1 ) approximately linear.Of particular interest in this regime are the 3-brane tensions, which are given by pairwise algebraic sums of the 4-brane tensions at the corresponding boundary points.As ω 1 < c 1 , the forms of Eq. (3.6)−Eq.(3.9) imply that all the 4-brane tensions are decidedly non-negative except V 2 (z), which is difficult to tell at first glance.If V 2 (z) is positive, all the corner branes have positive tensions, which is a satisfactory feature as negative tension branes are often intrinsically unstable.The crucial sign-determining term of V 2 (z) is tanh ln 1 , and to +1 for 1 .Consequently, one ends up with two possible tensions of the y = π brane, corresponding to the two allowed values of the modulus x 1 that solve the hierarchy problem: On the other hand, the tension V 1 (z) of the y = 0 brane is independent of x 1 , and is approximately equal to V 2 (z) for small ω 1 .For 1 , all the 4-brane tensions are thus positive, with the sole exception of V 3 = 0.The tensions of the four 3-branes at the corners are therefore strictly positive.Importantly, if we identify the maximally warped (y = π, z = 0) brane as the visible brane, then its tension is simply equal to V 2 (0).This demonstrates that the second solution x

Dominant warping along z (α << β)
In the opposite regime, we have β 10, which implies an extremely tiny upper bound of say α ∼ 10 −12 so as to eliminate the unwanted R y /r z hierarchy.Thus, the contribution of α becomes negligibly small in the warp factor a(π), which approaches unity and corresponds to an infinitesimally small n 1 << 1.Hence, warping occurs predominantly along z, with the extent of shared warping between the two directions being many orders of magnitude smaller than the α > β case.
This regime is less interesting than the previous one, as it effectively reduces to a singly warped six-dimensional geometry with a nearly unwarped five-dimensional submanifold.However, it still allows the possibility that all 3-brane tensions can be positive.As before, the determining feature is the positivity of the hyperbolic term in V 2 (z), for which we require where we have used Eq.(3.3) to arrive at the second inequality.As the contribution of α is negligible, we need β 12 for an overall warping of order 10 −17 which can resolve the hierarchy problem.This leads to the following minimum value of the R y /r z ratio compatible with the condition that V 2 (z) is positive.
First, note that a small value of n 1 allows ω 1 to be correspondingly close to unity according to the constraint placed by Eq. (3.11).Conversely, for a negligibly small α, one can treat ω 2 1 as a free parameter and tune it sufficiently close to 1 to generate a very small warping n 1 along y.According to Eq. (3.18) and Eq.(3.19), such a fine tuning can render V 2 (z) positive while simultaneously ensuring no large hierarchy comes into play between R y and r z .In Figure 2, we have shown the roughly linear dependence of (R y /r z ) min on the tuning of |ω 1 | close to unity.As an example, for (R y /r z ) min to be no larger than 10 2 , this tuning has to be very delicate with (1 − |ω 1 |) 10 −23 .This reintroduces the fine tuning problem in a new guise, and renders the prospect of having a positive tension (π, 0) brane in this regime problematic, albeit theoretically possible.Note that although ω 2 1 needs to be tuned very close to 1, the smallness of α nonetheless manages to keep ΩR 2 y bounded above by a very small value.For example, α ∼ 10 −14 roughly gives |Ω|R 2 y 10 −28 , which is similar to the upper bound obtained earlier for α > β.
In the small α regime, the curved braneworld scenario exhibits another interesting property.Due to unequal warpings along y and z, a salient feature of the doubly warped model (both curved and flat) is the clustering of pairs of 3-branes around the fundamental scale and the maximally warped scale, e.g. for α << β, the (0, π) and (π, π) branes cluster around the fundamental scale M , while the (0, 0) and (π, 0) branes cluster around the maximally warped (TeV) scale (see [50]).As the effect of the induced cosmological constant shows up only through the subdominant warp factor a(y), we expect Ω to directly affect the degree of clustering between the mass scales of the (0, 0) and (π, 0) branes.The total warp factors are a(0)b(0) and a(π)b(0) respectively.Since a(0) = 1, the ratio of the maximally and near-maximally warped scales is simply equal to a(π).Owing to the smallness of α, we can expand a(π) up to first order in α to obtain its dependence on ω 2 1 .
Mass scale of maximally warped (π, 0) brane Mass scale of near-maximally warped (0, 0 For all 0 < ω 2 1 ≤ 1, this ratio is larger than (1 − απ) as obtained for the flat case with ω 1 = 0.In other words, the scales of the TeV-branes are clustered more closely when they are negatively curved than when they are flat.For α ∼ 10 −14 , the ratio virtually tends to unity for the aforementioned tuning of ω 2 1 close to 1 which is necessary for a positive tension (π, 0) brane.If the latter condition is relaxed, the ratio can be comparatively smaller than unity.For all the cases though, deviations from the flat limit are imperceptibly small from a practical point of view, as the smallness of α allows any difference to show up only at the 14th decimal place or beyond.
It should be emphasized, however, that these results depend crucially on how stringently we wish to avoid a hierarchy between R y and r z , which, in turn, dictates how small an α needs to be chosen.For example, if one ceases to bother about the former, one can choose significantly larger values of α instead, leading to larger differences from the corresponding flat results.Of course, such choices are largely not compatible with the naturalness assumption and simply bring back the hierarchy problem in a new avatar, hence are not helpful from a physical point of view.To drive home the point, we tabulate in Table 1 a few combinations of α and ω 1 which affect TeV scale clustering in a pronounced manner, alongside the orders of corresponding R y /r z ratios.α a(π) in flat limit ∼ R y /r z ω 1 a(π) for given ω 1 of the TeV scale (0, 0) and (π, 0) branes, as compared to the flat limit result.However, they are all accompanied by dangerously large radii ratios which spoil naturalness and introduce a new hierarchy in the model.For all the combinations, β 12 has been fixed beforehand to obtain the desired amount of dominant warping along z.
4 dS-Schwarzschild 3-branes (Ω > 0) In this regime, the general solutions of Eq. (2.10) and Eq.(2.11) are given by As before, the usual normalization of b(π) = 1 yields ω2 = 1 and renders c2 a superfluous constant which can be set to unity.This makes Eq. (3.3) hold good.Further, normalizing with a(0) = 1 gives the following solution of c 2 where we consider only the solution which allows c 2 → 2 so that a(y) → e −αy in the limit ω 2 → 0. Unlike the AdS case, Eq. ( 4.3) itself does not require ω 2 2 to be bounded above by unity.However, in order to ensure positivity of the warp factor a(π), the relation ω 2 /c 2 < e −απ needs to hold, furnishing an upper limit for ω 2 corresponding to any given value of α.The relation between β and r z remains the same as in Eq. (3.5).Evaluating the brane tensions proceeds in the same manner as described earlier and yields: Like in the AdS regime, one cannot have large simultaneous warping along both directions if one precludes a large R y /r z hierarchy.This can be shown by setting a(π) = 10 −n1 and b(0) = 10 −n2 as before, and solving for α and β provided n 1 ∼ n 2 and n 1 + n 2 16.The allowed regimes remain the same, i.e., either n 1 ∼ 15 corresponding to α ∼ 10 and β ∼ 0.1, or n 2 ≈ 16 with β ∼ 12 and α ∼ 10 −16 .The analysis is not repeated here as it is exactly the same as in the AdS case.Instead, we move directly on to the two possible scenarios.

Dominant warping along y (α > β)
The crucial difference with the AdS case is the fact that for α > β, the large warping n 1 ∼ 15 along y can be generated only by a unique value of α, unlike the two degenerate solutions found in Eq. (3.13).This can be seen by simply solving a(π) = 10 −n1 explicitly, which gives the solution where x 1 = απ as defined earlier.The other solution with a negative discriminant renders e −x1 < 0, hence must be discarded.
It is straightforward to see that all the brane tensions from Eq. (4.4)−Eq.(4.7) are strictly non-negative, with the sole exception of V 2 (z) whose nature depends on the sign of the constant hyperbolic coefficient as follows.
coth ln (4.9) The numerator is clearly positive, and so is the denominator by virtue of the positivity of the warp factor and Eq.(4.8).This makes the coth term positive, which renders V 2 (z) from Eq. (4.5) negative.Thus, it is impossible to have a positive tension (π, 0) brane in this case.
An interesting feature that generalizes readily from the singly warped counterpart is the variation of N (defined here via ω 2 2 = 10 −N ) with x 1 .Away from the divergence at the flat limit x 1 = n 1 ln10, the value of N decreases sharply as shown in Figure 3, leading to rapid increase of ω 2 2 .Thus, the observed tiny value of the cosmological constant requires x 1 (hence α) to be very close to the flat limit.Based on the relation between α and β from Eq. (3.3) and that between β and r z from Eq. (3.5), this naturally explains the tuning of r z close to the Planck length when the bulk cosmological constant Λ is of order −M 6 .

Dominant warping along z (α << β)
Like its AdS counterpart, this regime admits a vanishingly small α ∼ 10 −14 − 10 −12 and a unique β 12 to generate the large mass hierarchy without any unnaturally large R y /r z ratio.Consequently, the metric of the five-dimensional submanifold is left almost unwarped along y.
In accordance with Eq. (4.5), the criterion for V 2 (z) to be positive is given by coth ln Owing to the positivity of the warp factor a(π), this condition cannot be satisfied by any value of ω 2 2 or x 1 .So it is not possible to have a positive tension (π, 0) brane in this case either.Similar to the AdS case, here also one can figure out how the clustering of the O(TeV) mass scales depends on Ω.The ratio between the scales of the maximally and near-maximally warped 3-branes remains equal to a(π), which in this case yields the following result: Mass scale of maximally warped (π, 0) brane Mass scale of near-maximally warped (0, 0) brane = a(π) For ω 2 = 0 this boils down to (1 − απ) as expected.To arrive at non-trivial results, let us first recall that ω 2 in this regime can be safely larger than unity as long as it satisfies ω 2 /c 2 < e −απ as meant to ensure the positivity of a(π).This gives rise to three distinct possibilities: 1.One may choose α ∼ 10 −12 and 0 < ω 2 ≤ 1, for which the deviation of the ratio from its flat limit is too minute to be phenomenologically significant.
2. As in the Ω < 0 case, the deviation in case of the previous possibility is more pronounced for larger values of α, but they are all accompanied by extremely large R y /r z ratios which render them unappealing from a physical standpoint.
3.Even for a conservative α ∼ 10 −12 which results in no large R y /r z hierarchy, choosing ω 2 >> 1 (which still satisfies the inequality written above) can produce a significantly smaller value of a(π) compared to its flat counterpart.
The first two possibilities are shared by the AdS case as seen earlier, while the third is unique to dS 3-branes and hinges on the possibility that ω 2 can be larger than unity.In Table 2, we list a few representative combinations of α and ω 2 for both cases (2) and (3).The latter case implies 10 −6 ΩR 2 y 10 −2 for the values that have been chosen, but this magnitude can be lowered for lesser values of ω 2 (i.e. at the cost of closer clustering).Note that the second and third columns are the same as those of Table 1, as the same flat limit is recovered for both ω 1 = 0 and ω 2 = 0, and the functional dependence of the R y /r z ratio on α and β is the same in the Ω > 0 case as for Ω < 0. Furthermore, the clustering ratio a(π) is smaller than its flat limit for any non-zero ω 2 , i.e. the scales of the maximally and near-maximally warped branes are farther apart from each other when the branes are positively curved than when they are flat -which is the opposite of the corresponding AdS behaviour.
The clustering ratio a(π) can be readily related to the mass splitting ratio (r) of two fermions localized on the proximal branes, which is an interesting feature of multiply warped models.Following the method shown in [53], it is easy to see that in our case r = a(π) 3/2 .The fact that the degree of clustering, hence the proximity of fermion masses, increases with decreasing ω 2 is interesting from a cosmological standpoint.First, let us consider the case of any two SM fermions located on the proximal O(TeV) branes.Depending on the choice of SM fermions, the observed magnitude of r can range from 10 −2 (between the 172 GeV top quark and the 4.2 GeV bottom quark) to 10 −12 (between the top quark and the ∼ 1 eV neutrino mass scale).Now, we firstly need α 10 −10 to preclude any R y /r z hierarchy larger than 10 4 .Even this threshold value leads to the requirement of a large ω 2 of order at least 10 9 if we wish to generate r ∼ 0.01 (smaller values of r require even larger values of ω 2 tuned sufficiently close to its allowed maximum value).Assuming r z ∼ M −1 and hence R y ∼ 10 4 M −1 , one ends up with Ω/M 2 ∼ 10 −10 for such a high value of ω 2 (it can easily be shown that w 2 ∼ 10 n roughly corresponds to Ω/M 2 ∼ 10 2n−30 in the small α regime), which is far too large compared to its observed value of about 10 −120 .At first glance, this rules out using ω 2 as a regulator to consistently generate the presently observed mass hierarchy between any two distinct SM fermions on the clustered branes.However, in conjunction with other independent mechanisms (see for example [74][75][76][77][78][79]) which attempt to resolve the discrepancy between the theoretically predicted large value of Ω and its observed small value, the cosmological constant could still be made to function as a tuning parameter for r.However, the need to invoke such external mechanisms renders the entire setup non-minimal in nature and less appealing.
Alternatively and perhaps more interestingly, instead of readily generating such a hierarchy, the model offers room for an intriguing hypothesis.Conventional wisdom suggests that a much larger cosmological constant may have prevailed during the very early stages of the universe, with the end of inflation marking its switchover to the present scale.As we have just seen, such a large value of Ω can produce r 0.1 without difficulty.On the other hand, when the inflaton has started decaying into the SM spectrum during the process of preheating, and the cosmological constant has not yet been reduced all the way down (say Ω/M 2 ∼ 10 −12 ), the number of degrees of freedom is naturally expected to be quite high.Among the particles produced under such conditions, one may hypothesize the existence of multiple species of fermions with nearly identical physical properties but masses separated by orders of 10 −1 or 10 −2 .In our model, such a mass splitting can be interpreted as a result of the large value of Ω, providing an unified framework for both ideas.By the end of the reheating period, when Ω has reached its currently observed value, this hierarchy becomes extremely close to unity (up to O(10 −x ) for α ∼ 10 −x with x ∼ O(10)), and the different species can effectively appear as a single species to low energy observers for the rest of cosmic history.Thus, an SM fermion observed today could actually be a family of distinct but nearly-identical fermions with an extremely close-knit mass spectrum which may be probed directly only at the GUT scale or beyond, whereas in the early universe they might have had a far more pronounced splitting induced by the larger cosmological constant.In the doubly warped setting, this works for only two fermions, while higher dimensional extensions allow the inclusion of more species due to a larger number of near-maximally warped 3-branes.A full treatment of this idea requires a study of the cosmological dynamics of the curved doubly-warped braneworld model, which falls outside our current scope.So we mention this here only as an interesting possibility offered by this model, and defer such a treatment to a future work.

Prospect of near-maximally warped visible brane
So far we have identified the maximally warped (π, 0) brane as the visible brane, based on the conservative assumption that the latter should have the least possible mass scale among all the existing 3-branes.As shown in the preceding sections, this identification allows a positive tension visible brane only for Ω < 0. Relaxing this assumption, we may identify the near-α a(π) in flat limit ∼ R y /r z ω 2 a(π) for given ω 2 the TeV-branes in the α << β regime, as compared to the flat limit result.The first three sets of α, corresponding to case (2) from Section 4.2, are accompanied by dangerously large radii ratios which introduce a new hierarchy in the model.The last three sets, corresponding to case (3), are free from this issue and make use of ω 2 >> 1 allowed in the dS regime.For all the combinations, β 12 has been fixed beforehand to obtain the desired amount of dominant warping along z. maximally warped (π, π) brane as the visible brane instead, but only at the cost of allowing a brane with a slightly lower physical mass scale to exist in our system.In this section, we briefly focus on the implications of such a possibility for the visible brane tension, for each of the four cases considered so far.

Ω < 0
For α << β, the near-maximally warped brane lies at (0, 0).Its tension is equal to V 1 (0) from Eq. (3.6), which is positive by default for all values of ω 1 .Thus, its positivity does not require any fine tuning of ω 1 , unlike that of the (π, 0) brane.This aspect makes it preferable to choose the (0, 0) brane as the visible brane in this case, which of course comes at the cost of allowing a slightly lower mass-scale (π, 0) brane that can have negative tension in absence of such fine tuning.Also, to resolve the hierarchy problem at (0, 0), we must have n 2 16 (corresponding to β 12), as a(0) = 1 contributes nothing to the warping along y.The upper bound on |Ω|R 2 y from Section 3.2 holds good, and in fact, is allowed to be even smaller as ω 1 does not need to be tuned close to 1 anymore.
For α > β, the near-maximally warped brane is the (π, π) brane.Its tension is equal to V 2 (π) + V 4 (π), which is given explicitly by For the second solution x (2) 1 from Eq. (3.14), the coefficient term of sech(x 2 ) approaches +1 as shown in Eq. (3.15), which renders the tension of this brane positive.For the other solution x (1) 1 , it approaches −1 and the setup mimicks the flat limit.In that case, the tension can still be positive for any x 2 > x 2r , where x 2r is the root of the equation 4tanh(x 2r ) = 3sech(x 2r ), given by x 2r ≈ 0.694 (or equivalently β r ≈ 0.221).On the other hand, in order to steer clear of a large R y /r z ratio, x 2 cannot be arbitrarily larger than this lower limit.For example, x 2 can comfortably range from 0.69 to 2.50 (say), where the latter is an upper limit assumed to be set by some modulus stabilization mechanism.The key takeaway is the fact that a near-maximally warped visible brane can have positive tension for both the solutions of x 1 that resolve the hierarchy problem, unlike only the second solution x (2) 1 as in Section 3.1.

Ω > 0
As before, for α << β, the tension of the near-maximally warped (0, 0) brane is given by V 1 (0) which is unambiguously positive.This marks a notable departure from Section 4.2, where it was impossible to have a positive tension visible brane.
For α > β, the near-maximally warped (π, π) brane has tension Using Eq. ( 4.1) along with a(π) = 10 −n1 and exploiting the positivity of the warp factor, the condition for the positivity of the tension can be recast as This sets an upper bound to the magnitude of ω 2 .Firstly, for a real value of this upper bound, one needs x 2 0.694.This is quite expectedly identical to the value of x 2r calculated earlier, as in the Ω → 0 limit Eq. (5.1) and Eq.(5.2) become identical, and any x 2 > x 2r can satisfy Eq. ( 5.3) trivially and give a positive tension.But x 2 cannot be arbitrarily large as argued earlier, and needs to be bounded above by an O(1) upper limit as well.This, in turn, safely ensures that the scales of the maximally warped (π, 0) brane and (π, π) are sufficiently close.Secondly, to resolve the hierarchy problem on (π, π), we require n 1 = 16, as b(π) = 1 gives zero warping along z.Together, these conditions furnish a small upper bound of roughly ω 2 10 −15 .Note that this is consistent with the ω 2 /c 2 < e −απ criterion necessary for the positivity of a(π), which was observed earlier in Section 4.
Thus, in the Ω > 0 regime, the tension of the near-maximally warped brane can be positive for both α << β and α > β, with the latter case additionally entailing a small upper bound on the magnitude of the brane cosmological constant as a requisite criterion.These features differ substantially from the results obtained in Sections 4.1 and 4.2, and render the possibility of a near-maximally warped visible brane in this regime far more interesting.

Higher dimensional extensions
The analysis so far can be readily extended to spacetimes with n-fold warping over a series of nested S 1 /Z 2 orbifolds, equipped with (n + 4)-dimensional metrics of the form where R i s are the orbifold radii, y i s are the compact angular coordinates, and a i s are the corresponding warp factors.Evidently, the metric elements gAB are of the forms: where gjj are the extra dimensional components.The bulk cosmological constant is denoted by Λ n+4 .The setup consists of 2n number of p-branes where p = n + 2, with one pair of branes at every set of fixed points y i = {0, π}.The intersections of these branes produce 2 n number of 3-branes at the vertices of the resulting "brane-box" configuration.Since the effect of curvature enters principally through the first level of warping, we expect the geometric properties of these extended models to closely mimick those of the curved doubly warped scenario.We demonstrate the validity of this claim by first analyzing the triply warped case, i.e., n = 3, for which the EFEs take the following forms (where primes denote derivatives wrt respective variables): µν-component: y 1 y 1 -component: y 2 y 2 -component: y 3 y 3 -component: In the bulk, we can separate the µν and extra dimensional variables from Eq. (6.3) and write G µν + Ωg µν = 0, where the effective cosmological constant (Ω) is defined as Using R = 4Ω from the 4D equation and plugging Eq. (6.7) for Ω in Eq. (6.4), the bulk cosmological constant can be expressed as The steps to uncouple the warp factors are similar to the doubly warped case: (i) Plug Eq. (6.8) into Eq.(6.6) to eliminate Λ 7 .(ii) Compare Eq. (6.4) and Eq.(6.6) to express a 3 a 3 −a 2 3 in terms of a 1 and a 2 − a result which when inserted into the previous equation gives the first uncoupled ODE for a 1 (y 1 ).(iii) Now compare Eq. (6.5) and Eq.(6.6) to eliminate a 3 in terms of a 2 alone, and insert into the same equation to obtain the second uncoupled ODE for a 2 (y 2 ).(iv) In the eliminating relation found and used in the previous step, rearrange terms to construct the final ODE for a 3 (y 3 ).The equations thus obtained have remarkably simple forms.
; a 2 − α 2 2 a 2 = 0 (6.10) where α 1 and α 2 are two dimensionless constants of separation.With appropriate normalization, the final solutions of the warp factors are given by: where α 3 /R 3 = −Λ 7 /(60M 5 ) analogous to Eq. (3.5).Once again, it is clear that the effect of Ω appears explicitly only in a 1 (y 1 ).The equations (and solutions) of the other warp factors appear thereafter in a hierarchical manner, with nested functional relations among the warping parameters and the orbifold radii.Equipped with these solutions, the brane tensions can be found by integrating Eq. ( 6.3) across each fixed point.The non-zero tensions are: Note that only the compact coordinates of higher orbifolds contribute to the coordinate dependence of each brane tension, which is a consequence of nested warping.
In order to avoid a large hierarchy among the radii, the regimes allowed are either α 1 ∼ 10 alongside α 2 ∼ α 3 ∼ 0.1, or any one of α 2 or α 3 close to 12 with the other two being vanishingly small (∼ 10 −15 ).Regardless of the choice of regime, maximum possible warping always occurs at {y 1 = π, y 2 = 0, y 3 = 0}.Identifying it with the visible brane leads to visible brane tension equal to V 2 (π), which, owing to its form, gives nothing quantitatively different from Sections 3 and 4. The case of a non-maximally warped visible brane is arguably more interesting, for which one can proceed as in Section 5 and analyze various possible configurations.As an example, consider the case Ω < 0 and α 1 ∼ 10, for which (π, π, π) is a near-maximally warped brane with tension V 2 (π, π) + V 4 (π) + V 6 .Recall that the hierarchy problem can be resolved for two distinct values α (refer back to Eq. (3.13)), which render the tanh term in V 2 respectively −1 and +1 .The latter identically makes the tension positive for all α 2 and α 3 , whereas for the former, the positivity condition can be written explicitly as Clearly, for any given α 3 > 0, the lower bound on α 2 is smaller than the value of 0.221 obtained in Section 5.1.In fact, it can be readily checked that for α 3 0.182, any α 2 > 0 can satisfy Eq. (6.15).The triply warped model thus allows greater flexibility over the doubly warped one in rendering the tension positive.The cases of the other near-maximally warped branes can be analyzed individually in a similar fashion, thereby opening up different regions of parameter space consistent with a positive tension visible brane.In Figure 4, we show the allowed regions in the {α 2 , α 3 } parameter space leading to positive tensions of each of the four near-maximally warped branes for Ω < 0 and α 1 ∼ 10, corresponding to the first solution α The generalization to n extra dimensions for the metric in Eq. (6.1) is straightforward at this point, for which the governing set of equations turns out to be , a j − α 2 j a j = 0 (6.16) where j ranges from 1 to n − 1. Depending on the sign of Ω, the corresponding solutions of a 1 (y 1 ) are identical to the doubly (and triply) warped results.All the subsequent warp factors are independent of Ω and take the forms a j (y j ) = cosh(α j y j ) cosh(α j π) , α j−1 = R j−1 α j R j cosh(α j π) (6.18) from Eq. (3.13).The red (bounded) and blue (unbounded) regions in the first three cases correspond to negative and positive tensions respectively.Note that for the maximally warped (π, 0, 0) brane, there is no blue region for α 1 = α (1) 1 as the tension is identically negative, whereas for α 1 = α (2) 1 the reverse situation arises.
for j = 2 to n, with α n /R n = −Λ n+4 /(20nM n+2 ).The odd-labeled p-branes (with p = n + 2) located at each y j = 0 for j = 2 to n have vanishing tension V 2j−1 = 0, while the remaining branes at y 1 = {0, π} and y j = π have tensions as follows: where ω 1 , c 1 , ω 2 , and c 2 are as defined earlier.To avoid a large hierarchy among the R j s, the allowed regimes are either α 1 ∼ 12 with α j ∼ 0.1 for all j = 2 to n, or any particular α j ∼ 12 with all the others (including α 1 ) vanishingly small (∼ 10 −15 ).Maximum warping occurs uniquely at (π, 0, ..., 0), with 2 n−1 number of near-maximally warped 3-branes (with a closely spaced cluster of mass scales) arising in each of the allowed regimes.Proceeding as earlier, a case-by-case analysis of the tension of each near-maximal brane can be done for any given n.For example, in the Ω < 0 and α 1 ∼ 10 regime with α 1 = α 1 , there exists a non-trivial region in the (n − 1)dimensional {α 2 , ... , α n } parameter space for which the tension of the (π, π, ... , π) brane can be rendered positive (analogous to the top-left plot from Figure 4).Using Eq. (6.19) and Eq. ( 6.21), one can deduce that this region corresponds to the following inequality: Constraints for the other branes can be obtained similarly.Since V 2j−1 = 0 for all odd j > 1, note that the inequality for a 3-brane produced by the intersection of any y j = 0 p-brane will be independent of α j .While in principle unbounded above, such regions in the parameter space cannot be arbitrarily large if we are to avoid -as noted earlier -an unnaturally large hierarchy among the various R i s.An upper bound might be established, for instance, via some higher-dimensional analogue of the modulus stabilization mechanism explored in [69], a rigorous discussion of which lies beyond the scope of this work.This is an interesting question that we plan to address separately in near future.
As for the dependence of the clustering of scales between the maximally warped and any nearmaximally warped brane on the induced cosmological constant, the results for an arbitrary n-fold warped metric should not differ appreciably from the six dimensional case as Ω makes its appearance solely through a 1 (y 1 ) and in the same functional form.So the values of a 1 (π) from Tables 1 and 2 continue to hold, with other relevant a j (0) terms (independent of Ω) appearing multiplicatively and providing subleading corrections to the overall scale ratio.

Conclusion
In summary, we have generalized six and higher dimensional braneworld models with nested warping to include the effects of a non-zero cosmological constant (Ω) induced on the corner 3branes, which results in global brane curvature.We have first analyzed the doubly warped case in detail, for which the solutions of the warp factors from the Einstein field equations permit both negative and positive values of Ω, corresponding respectively to AdS-Schwarzschild and dS-Schwarzschild geometries for the 3-branes.An important consequence of nested warping is that the effect of Ω appears explicitly only in the warp factor associated with the first orbifold.On the other hand, simultaneous large warping along both orbifolds is forbidden if the radii R y and r z of the compact spaces are of close magnitudes, which is a natural assumption to make.This is a salient feature of nested warped braneworld models that is seen to hold for both flat and curved scenarios.Together, these two properties result in 2 × 2 = 4 distinct possibilites which we have analyzed in detail, for which we have first assumed the maximally warped (y = π, z = 0) brane to be the visible brane accommodating our four-dimensional universe.
• For Ω < 0, dominant warping along the first orbifold (y) imposes a small upper bound of order 10 −32 on the value of |Ω|R 2 y , which is interesting from a cosmological perspective.For values of the cosmological constant sufficiently smaller than this allowed maximum, the gauge hierarchy problem can be resolved for two distinct values of the warp factor α associated with this orbifold (as opposed to a unique value in the flat case).One of these solutions effectively reproduces the flat result, while the other one is distinct and non-trivial.Unlike the former, the latter can render the tension of the maximally warped 3-brane positive.On the other hand, dominant warping along the second orbifold (z) requires an extreme fine tuning of the cosmological constant if one wishes to simultaneously avoid a large R y /r z hierarchy and have a positive tension (π, 0) brane.One also observes that the physical mass scales of the maximally warped (0, 0) brane and the near-maximally warped (π, 0) brane are clustered more closely in this regime than in the case of flat branes.If one strictly chooses R y /r z 10 2 , this difference is imperceptibly small and unlikely to play any significant phenomenological role.
• For Ω > 0, positivity of the warp factor along y imposes a similar upper bound to the magnitude of Ω.Assuming dominant warping along y further renders it quite small.But unlike the AdS case, there is only one solution of α that can generate the desired large mass hierarchy on the visible brane.The magnitude of Ω rises exponentially for small deviations of the warping parameter from its flat value, thereby justifying a flat braneworld approximation from a physical perspective.Furthermore, the tension of the (π, 0) brane cannot be rendered positive for any set of values of the model parameters.Dominant warping along z cannot provide a positive tension visible brane either.As for the clustering of the (0, 0) and (π, 0) branes' scales, they are found to be comparatively farther apart than their flat limit counterparts -which is opposite of the analogous AdS behaviour.Moreover, unlike the AdS case, large allowed values of the cosmological constant (e.g. 10 −6 ΩR 2 y 10 −2 ) can significantly modify the clustering ratio without jeopardizing a conservative R y /r z ratio.This, however, is inconsistent with the much smaller value of Ω measured today, and prevents explaining the currently observed SM fermion mass hierarchy (along the lines of [50] and [53]) by using the induced cosmological constant as a regulator, unless some separate mechanism (e.g.[74][75][76][77][78][79]) to resolve the discrepancy between the theoretically predicted and observed values of Ω is invoked.Alternatively, one might hypothesize that end-of-inflation (p)reheating could have naturally led to the production of families of nearly-identical fermions each with such closely-spaced mass spectra as associated with a large cosmological constant (consistent with inflationary dynamics).Thereafter, the lowered value of Ω could have greatly reduced this mass hierarchy and rendered it extremely close to unity (up to O(10 −12 ) or even smaller for conservative R y /r z ratios), practically removing the distinction between different species to any low energy observer throughout the rest of cosmic history.This implies that any SM fermion observed today could actually be a tuple of two (for the doubly warped model) or more (for higher dimensional extensions) nearly-identical, fundamental fermions with extremely close masses that might only be distinguished at the GUT scale or higher.From a theoretical standpoint, this provides a unified framework which potentially accommodates both the natural consideration of additional degrees of freedom in the post-inflationary period (besides explaining why they are no longer distinctly detectable today) and the large cosmological constant briefly dominant during that period.In order to fully address this question however, one needs to study the cosmological aspects of this braneworld model in detail.This falls outside the scope of the present work, and we plan to address it in a future one.
A few of these features strikingly resemble those of the curved singly warped scenario analyzed in [68].As mentioned earlier, this is but a consequence of secondary warping of the alreadywarped 5D metric from [1] along a higher S 1 /Z 2 orbifold.Some further novelties concerning the visible brane tension can emerge by identifying, instead of the conventionally chosen maximally warped brane, the near-maximally warped 3-brane with our 4D universe.Recall that such a choice is possible only in nested braneworld models with at least two levels of warping.In the Ω < 0 regime, dominant warping along z gives the near-maximal (0, 0) brane a positive tension by default.Dominant warping along y, on the other hand, can potentially render the tension of the (π, π) brane positive for both solutions of α that generate the desired large mass hierarchy, as opposed to only the second solution which works for the maximally warped brane.In the opposite Ω > 0 regime, dominant warping along z similarly results in a tension which is always positive.warping along y can do the same as long as there is a small upper threshold to the value of Ω.As readily apparent, these features differ significantly from what one obtains for a maximally warped visible brane.
The curved doubly warped model can be readily generalized to seven and higher dimensions with a series of nested warpings over successive orbifolds.The equations for the warp factors, together with the relations among the warping parameters and the orbifold radii, follow a definite hierarchical pattern as seen from Eq. (6.16)−Eq.(6.18).Since the effect of the induced cosmological constant enters only through the first level of warping, the dependence of scale-clustering between the maximal and any near-maximally warped brane on Ω remains virtually identical to the 6D case.Moreover, identifying near-maximal TeV-branes (instead of the conventionally chosen maximal one) as the visible brane can potentially open up non-trivial regions in the warping parameter space that can render the corresponding brane tension positive.
Based on the results derived in this work, one can explore various phenomenological aspects of multiply warped non-flat braneworld models.The issue of moduli stabilization in curved higher dimensional scenarios is a particularly interesting avenue.Previously, the authors of the present work have demonstrated stabilization of the two moduli of the 6D flat warped braneworld model via a straightforward extension of the Goldberger-Wise mechanism [59].While a non-flat generalization of such a bulk field approach is expected to work well, it is tempting to check if modified gravity effects and/or non-zero brane curvature alone can produce a stabilizing potential without the need for any external field.This can indeed be realized in the non-flat 5D model, where the resulting radion has interesting phenomenological properties and cosmological implications [9,11,69].If such effective potentials governed solely by the gravitational sector can be shown to exist in multiply warped geometries as well, the phenomenology of the associated 4D scalar radions and their cosmological dynamics (e.g.their role in driving multi-field inflation) would be an area of considerable interest -one which we plan to explore in a future project.The inclusion of higher curvature terms in the gravitational action at sufficiently high energy scales (motivated here by the large bulk cosmological constant Λ 6 ∼ −M 6 ) is also likely to result in non-trivial modifications to the model.The effective dynamics of bulk matter and gauge fields against a curved warped brane background is another important area which needs to be focused on in greater detail.The issue of the mass discrepancy of the Standard Model W -boson [60], for example, falls within the purview of this domain.These questions lie beyond our current scope as well, and we plan to address them separately in future works.

( 1 ) 1 ) 1 10 − 30 . 2 y
is very close to the value obtained in the flat brane scenario and the other one (x (2) 1 ) slightly larger.As noted earlier, n 1 ∼ 15 puts a very small upper bound on the induced cosmological constant by ensuring ω 2 For α ∼ 10, Eq. (3.1) implies that this bound corresponds to |Ω|R 10 −28 .

( 2 ) 1 ,
which arises solely in the context of non-flat geometry, is of key importance in giving us a positive tension visible brane -a feature impossible in the flat case.

14 -Figure 2 :
Figure 2: For negative brane cosmological constant and α << β, logarithmic plot of the magnitude of the minimum admissible radius ratio (Ry/rz) min consistent with positive tension of the y = π brane, versus the required degree of fine tuning of ω 2 1 close to unity.

Table 1 :
A few representative combinations of α and ω 1 which are shown to considerably affect the clustering

Table 2 :
A few representative combinations of α and ω 2 which are shown to considerably affect the clustering of