From sine-Gordon to vacuumless systems in flat and curved spacetimes

In this work we start from the Higgs prototype model to introduce a new model, which makes a smooth transition between systems with well-located minima and systems that support no minima at all. We implement this possibility using the deformation procedure, which allows the obtaining a sine-Gordon-like model, controlled by a real parameter that gives rise to a family of models, reproducing the sine-Gordon and the so-called vacuumless models. We also study the thick brane scenarios associated with these models and investigate their stability and renormalization group flow. In particular, it is shown how gravity can change from the 5-dimensional warped geometry with a single extra dimension of infinite extent to the conventional 5-dimensional Minkowski geometry.


Introduction
Topological solutions in field theory are related to many phenomena in physics [1][2][3]. In particular, when we have models in 1+1 dimensions involving a potential written in terms of a scalar field and with a set of degenerate minima, such solutions represent transitions between consecutive minima and are called kinks. The minima define the possible vacuum states at the quantum level, and they are distributed over the many values for the field. Such distribution can appear in the most diverse ways, for instance, in the φ 4 -model one has only one topological sector, defined between two consecutive minima, while in the sine-Gordon model [4] one has an infinite copy of the same sector, always between two consecutive and well-localized minima. A system that is quite different in this perspective is the vacuumless model [5,6]. In this model we still have a topological sector, which can be interpreted as connecting two minima of the scalar potential, but now they are located at infinity. In this case, the field solution is asymptotically divergent and has infinite amplia e-mail: bazeia@fisica.ufpb.br tude, but keeps its topological character well behaved [7]. It is worth mentioning that vacuumless systems appear in a diversity of contexts in high energy physics [5][6][7][8][9][10][11].
The aim of this paper is to construct a model that, with the proper variation of a given parameter, makes the transition between systems with well-located minima and systems that support no minima at all. The model presented here has as limit cases the sine-Gordon and the vacuumless models. Using the deformation procedure developed in Ref. [12], we find a field transformation that takes us to a new system with a double sine-Gordon-like behavior, which has two infinitely degenerate sets of solutions. One of these sets transits between the sine-Gordon kink and the vacuumless solution. The second set coincides with the previous one in one of its limits and, except from a phase, in the other limit case it is destroyed, remaining only the zero-energy solutions. Thus, in one of the limits of the model the two phases coincide, generating infinite copies of the same topological sector, and in the other limit only one topological sector survives.
Models in field theory motivate generalizations of the Randall-Sundrum model [13,14] in the presence of scalar fields [15][16][17][18][19][20][21][22][23][24], which are known as thick branes. For this reason, in this work we also study the thick brane scenario generated by the scalar field model which we first introduce and study in the flat spacetime. In particular, we find a brane that transits between the sine-Gordon brane [18] and a flat 5dimensional spacetime with two zero-energy solutions for the scalar field. The idea here is similar to the case investigated in [25], in which the authors propose a braneworld scenario where the brane changes from a thick to a thin behavior. In the current study, however, we describe a mechanism in which a single parameter can be used to control the brane profile, contributing to change the 5-dimensional warped geometry into a flat geometry.
Although the braneworld model that we explore below is more involved, it also supports analytical solutions. Thus, in the 5-dimensional spacetime with a single extra dimen-sion of infinite extent, we also analyze the stability of the braneworld scenario against tensorial perturbations. In the sense of gauge/gravity duality [26][27][28][29][30], where the extra dimension can be identified with the energy scale of the holographic dual field theory, we also study its implications for the renormalization group flow (the RG flow) in the dual field theory, since stability of the gravitational sector has relevant information as regards the dual model [31].
The subject to be explored in the current work is organized as follows. In Sect. 2 we review several aspects of the firstorder formalism for a single real scalar field and its relationship with the Bogomol'nyi-Prasad-Somerfield (BPS) solutions [32,33]. Then we discuss the deformation procedure and show how it acts to generate the new model in the 2dimensional spacetime. We go on and study in Sect. 3 the new model and its topological solutions, including how the energy and energy density behave. In addition, we also show that the solutions obtained are stable. Moreover, in Sect. 4 we analyze the properties of the thick branes that can be constructed from the model. We also investigate stability of the brane against metric fluctuations and study implications of the RG flow for the dual field theory. We end the work in Sect. 5, adding some comments and conclusions.

First-order formalism
The behavior of a scalar field is usually encoded in a Lagrangian density having the general form Here φ is the scalar field and V (φ) is the potential of the model, which determines how the field behaves. In the flat spacetime with (1, 1) spacetime dimensions, the metric tensor becomes η μν = diag(1, −1), so the scalar field only depends on the two coordinates x and t, i.e., φ = φ(x, t). For simplicity, we also work with dimensionless fields and coordinates. The equation of motion for the scalar field derived from the Lagrangian (1) is given by As the Lagrangian (1) is Lorentz invariant, we can focus on static solutions, since traveling waves can be obtained from a Lorentz boost. For static configurations Eq. (2) becomes a second-order differential equation given by Another important quantity we are interested in is the energymomentum tensor In particular, its T 00 component provides the energy density of the solution we are looking for. We represent it by ρ(x), which is explicitly given by Note that, for positive-definite energy, the potential must be non-negative, i.e., V (φ) ≥ 0. A powerful tool in the treatment of these models is the use of an auxiliary function, denoted by W (φ), which is introduced as follows: where W φ = dW /dφ. In this case the expression for the energy density (5a) becomes and Eq. (3) can now be given as the first-order differential equations It implies that the quadratic term in (7) disappears and, as a consequence, the energy density is only related to the xderivative of W . Thus, the energy of the model is determined only by the asymptotic behavior of the function W in the coordinate space; that is, one can write In this case the energy (9) is called BPS energy [32,33]. Relevant phenomena occur when the system under analysis presents a set of degenerate minima. In these situations each pair of consecutive minima form distinct topological sectors that, in turn, have different solutions. These solutions are called kinks. The simplest case of models having such properties is the well-known Higgs prototype or φ 4 -model, defined by the potential which has two degenerate minima at φ = ±1 and a topological sector having a kink solution given explicitly by φ(x) = tanh(x).
An interesting way to characterize a kink is the existence of a topological current. Here, we define it by [7] j μ = μν ∂ ν W (φ(x)), (11) where μν is the antisymmetric symbol in two dimensions with 01 = 1. Associated with this current (11) we have a topological charge given by Despite the similarities in the values of the charge (12) and the BPS energy (9), they have fundamental differences. While the BPS energy is associated with a continuous symmetry and can be identified from the Noether theorem, the charge (12) is associated with the topology of each solution, and results from the transition of the topological solution in between two minima of the potential.

The deformation procedure
Once we know the behavior of a given model, with its characteristics and general behavior, it is interesting to look for new well-behaved models. In this sense, the deformation method [12] is a powerful tool to find new models in field theory. The method consists of choosing a theory with the Lagrangian L(φ, ∂ μ φ) having the form (1) and then perform a transformation of the type We then get where with the potential U (χ ) given by This is a field redefinition, but if we consider the model (15) described by the potential (16) as a new model, and call it the deformed model, in this case there is aW function such that the first-order differential equations withW χ = W φ (φ → f (χ ))/ f χ , provide solutions for the field χ(x). In particular, among the characteristics of the new model, we can highlight that the total energy of the solution is since the deformed model also has a first-order structure. However, the deformed model engenders another important property: the transformation (13) allows us to find the solution for the field χ(x) without the need to deal directly with Eq. (17). What we have is that the solution to the new field is obtained by inverting the transformation (13). Thus, one finds where φ(x) represents a solution for the previous model. In this sense, Eq. (19) represents the main link between the φmodel and the deformed χ -model. For more details, see [12].
Although we initially present the model with two parameters, we are interested in situations where only one parameter is necessary to describe the changes in the model. So we assume that θ = θ(λ). We also impose on θ(λ) the conditions θ(0) = 1 and θ(1) = 0 in order to find the vacuumless model [7] for λ = 1 and the sine-Gordon model [4] for λ = 0. As a consequence we obtain the particular cases W (χ , 1) = 4 tan −1 (e χ ) − π and W (χ , 0) = sin(2χ)/2, so the sine-Gordon and vacuumless models appear in the system as The parameter λ has an interesting behavior: when it increases from 0 to unity, the model changes from the sine-Gordon to the vacuumless model. Physically, it transforms the periodic sine-Gordon potential (25) to a non-periodic one, the hyperbolic potential (25b) which defines the vacuumless model. The distinct solutions appear in Fig. 2 Fig. 3 one illustrates how the energy density becomes more and more diffuse, as λ increases from zero to unity.

and in
Once the model has been presented, we must solve the first-order equation  Here we observe how the small kink tends to disappear in the limit λ → 1 Equation (26) has two infinite set of solutions. The first one we call large kinks and they are given by The behavior of the solution (27) is depicted in Fig. 2 for n = 0. It asymptotically approaches . We note that, as θ → 0, χ L ,λ (±∞) tends to diverge and χ L ,λ (x 0) remains well behaved. For θ → 1 we do not have pathologies in the solution. Thus, the large kink describes, as a function of the λ parameter, systems that transit between the sine-Gordon kink and the vacuumless solution. Particularly, we have The second set of solutions we find from Eq. (26) is and we call it a small kink. Its shape is depicted in Fig. 3 for n = 0, except for the phase K λ . The solution (29) asymptotically approaches The small kink has a phase K λ , which causes its topological sector to move away from the center of the potential (23) to infinity as λ → 1. With a suitable choice for θ(λ) we can destroy the topological sectors associated with small kinks when λ = 1. Thus, only the topological sector of the vacuumless solution remains at that point. If we drop the phase in the small kink (29), it is possible to note that the zero-energy solutions also remain, but at infinity. A simple choice in this direction that also obeys the conditions (θ(0), θ (1) With this choice for θ(λ) we can explicitly rewrite the solutions of large and small kinks as So now we have a double sine-Gordon-like model with two field solutions coming from the two manifest topological sectors. Both kinks retrieve the sine-Gordon solution when λ = 0, but as λ grows, such solutions have distinct properties. Large kink becomes more diffuse until reaching the solution of the vacuumless model, which has divergent amplitude. It implies that at this point all topological sectors of the model are sent to infinity, except for the sector that is at the center of the potential (23). A small kink also becomes more diffuse as λ → 1, but it is destroyed when λ = 1. As we shall see later, it happens because as the associated topological sector moves away from the center of the potential its energy approaches zero. Moreover, one can show that the mass of the meson in the minima of the potential is given by where χ min is a minimum of (23). Here, we choose to write the model in terms of that quantity, whenever possible.
We can now perform the analysis of the energy densities of the model. The energy density of the large kink (27) is given by Fig. 4 The energy density (31) for some values of λ. Here one notes how ρ L ,λ behaves as λ → 1, becoming more diffuse but keeping its localized profile and its shape is depicted in Fig. 4. It makes the transition between the curves ρ L , . We do not have information for the case λ = 1 in the general asymptotic expansion for the energy density, but a direct approach in ρ L ,1 (x) leads to the expression ρ L ,1 (x → ∞) 1/x 2 + O(x −4 ), which decays much slower than the exponential. This change in the asymptotic behavior is due to the mass scale of the quantum meson which, for λ = 1, is zero. On the other hand the behavior of (31) in the vicinity of x = 0 is given by , which shows that the central portions of the energy density increases and becomes a bit more concentrated as λ → 1.
For the small kink we have the energy density given by and its behavior is depicted in Fig. 5. Equation (32) transits between the curves ρ S,1 (x) = 0 and ρ L ,0 (x) = sech 2 (2x).
, revealing that we cannot distinguish the energy densities of the fields (27) and (29) . It shows that, despite the similar asymptotic behavior (for λ = 1), the evolution of (31) and (32) in terms of the parameter λ is very different. The height and width of the energy density (31) grow with λ, becoming more diffuse but still with a localized pro- Fig. 5 The energy density (32) for some values of λ. Here one notes how ρ S,λ behaves as λ → 1, becoming more delocalized, disappearing at λ = 1 file. Meanwhile, as λ grows, the width of the energy density (32) increases, but its central portion decreases and becomes less concentrated, which implies that as the solution evolves it becomes more diffuse, but also delocalized. As a consequence, the area under the energy density (32) become smaller until it disappears, at λ = 1.
Integrating (31) and (32) we find the energy of large and small kinks, given by and respectively. The corresponding behaviors are depicted in Fig. 6. Here one observes that the energy of the solutions (27) and (29)

Linear stability
In this section we analyze the stability of the solutions of the models presented so far. The usual procedure is to take a timedependent perturbation around the static solution written in the form χ(x, t) = χ(x) + n η n (x) cos(ω n t), for small η n (x), and then insert χ(x, t) into (2). The procedure gives Fig. 6 The energy behavior as a function of λ. Here we observe that which is a Schrödinger-like equation with a stability potential given by Inserting the large kink solution (27) which is depicted in Fig. 7. For λ = 1 we have v L ,1 (x) = 4 8x 2 − 1 / 4x 2 + 1 2 , which is a volcano potential, and for λ = 0 we have v L ,0 (x) = 4 − 8sech 2 (2x), which is a reflectionless potential. Equation (37) has a global minimum at v L ,λ (0) = (1 + λ) λ 2 + λ − 4 ; this minimum increases or decreases, depending on λ being above or below the point atλ = 1 3 √ 13 − 2 . It shows that the potential (37), as λ grows, becomes deeper and after λ =λ it comes back to its initial depth. The change in the shape of the stability potential, in this case, is due to the behavior of the meson mass in the minima of (23). As λ → 1 we have m 2 λ → 0. Note that v L ,λ (0) < 0 for all λ and translational invariance requires the existence of at least one bound state, so the transition between the reflectionless and the volcano shapes in this case describes the tran- Fig. 7 The stability potential v L ,λ (x) given by (37) for some values of λ. When λ = 0, we have a reflectionless potential and as λ increases the shape of the potential changes to become of the volcano type sition from systems with massive meson to systems having a massless meson. This is another relevant physical behavior induced by the parameter λ, which will lead to distinct possibilities when used to describe braneworld scenarios, as we discuss in Sect. 4. Now, inserting the small kink solution (29) in the general expression for the stability potential (36) which is depicted in Fig. 8. For λ = 1 we have v S,1 = 0 and for λ = 0 we have v S,0 = 4 − 8sech 2 (2x), as expected.
The translational invariance of the solutions we presented so far implies the existence of at least one bound state for each topological sector, which is given by the zero mode of Eq. (35). In the formula, if we represent the zero mode by η 0 (x), it is the derivative of the field solution, For large and small kinks the shapes are depicted in Fig. 9, where one can observe how (37) holds the nice behavior of its zero mode, and how the zero mode of (38) gets smaller until it disappears when λ = 1.

Braneworld
Models described by scalar fields have direct applications in gravitation, providing braneworld scenarios for thick branes. In this context, the scalar field acts as a source of gravity around the brane, and thus describes how gravity behaves throughout the bulk. The system we are interested consists of a 3-brane embedded in a (4+1) spacetime with an extra dimension of infinite extent. The background geometry can Fig. 9 The zero modes for the large kink (27) (right) and the small kink (29) (left), depicted for some values of λ. Both modes behave adequately for λ = 1. However, for λ = 1, we have η S,λ=1 = 0, meaning that the small topological sector disappears, while η L ,λ=1 (x ∞) 1/x 2 , which is a consequence of the presence of the massless meson be written in terms of a static warped metric given by Here a, b = 0, ..., 4, μ, ν = 0, ..., 3, ds 2 4 = η μν dx μ dx ν , and the y-coordinate describes the extra spatial dimension. The functions A(y) and e A(y) are called warp function and warp factor, respectively, and they are assumed to depend only on the extra dimension.
In this braneworld scenario, we are interested in models that can be described by the action where is the Lagrangian for the scalar field and, for simplicity, we assume 4π G 5 = 1. We also assume that φ = φ(y), i.e., the scalar field only depends of the extra dimension.
The Einstein equations that follow from the action (41) are where G ab is the Einstein tensor and T ab is the energymomentum tensor, similar to Eq. (4). The 00 and the 44 components of (42) are given by respectively, with the prime denoting derivation with respect to the coordinate y. With Eq. (43) at hand we can subtract the first equation from the second to obtain Equation (44) provides a way to rewrite the system in terms of first-order equations. For this purpose, we introduce a function W (φ(y)) in the equations through the relation As a consequence, the equation providing the solution for the scalar field is now In order to solve the equations of motion, these two first-order equations require that the potential obeys Equations (45) and (46) constitute the first-order framework and can be used to analyze possible scenarios of thick brane that can be generated by the models presented in the previous section.
The analysis of the thick branes scenario generated by the large kink (27) is similar to the case of small kinks (29), so we concentrate on the brane generated by the small kink. To simplify the investigation, let us perform a shift in the field, φ → φ − K λ . In this case, one gets Due to the shift performed over φ, the field solution of (46) becomes Its shape can be observed in Fig. 3, when we make the changes x → y and φ λ (y) = χ S,λ (y) − K λ . Note that, for λ = 1 the field approaches a constant, φ λ=1 (±∞) = ±csch −1 √ 2 . With these ingredients we analytically solve Eq. (45), together with the boundary conditions A(0) = A (0) = 0, to find where Here Li 2 (y) is the polylogarithmic function. The shape of e 2 A(y) is depicted in Fig. 10. We can observe that in the vicinity of the center the warp factor behaves like A(|y| 0) (1 − λ)y 2 + O(y 4 ), which implies that, as λ → 1, the brane becomes less and less localized, becoming effectively delocalized for λ = 1. Meanwhile, in the asymptotic regime |y| → ∞ the thick brane approaches AdS vacua with cosmological constant 5 ∼ −W (φ(±∞)) 2 , and we have A(|y| → ∞) − 1 3 E S (λ)|y|, where we used Eqs. (18) and (34). As E S (λ) decreases as λ grows (starting at λ = 0, where we have the sine-Gordon brane) the thick branes generated by Fig. 10 The warp factor that follows from (50) for some values of λ.
Here we see how it evolves from the sine-Gordon case when λ = 0 to a 5-dimensional flat spacetime when λ = 1 this model approaches AdS vacua with smaller cosmological constant as λ → 1. In particular, for λ = 1 the brane seems to fill the entire space, so we have no graviton localization anymore at this point, with all the modes dispersed through the extra dimension. The braneworld scenario is destroyed in the limit λ → 1.
From Eq. (43) we find that the potential has to obey (47) and one notes that it keeps the correct form for vacuum stability in gravitational theories [34,35]. In particular, for λ = 1 the potential U (φ) vanish, and the spacetime solution reaches a flat 5-dimensional geometry. At this point, only the dynamical term ∼ ∂ a φ∂ a φ survives in the action (41), having as solution the constants φ λ=1 (±∞). The presence of two values possible for the scalar field in the flat background is due the original Z 2 -symmetry, which act as a memory related to the kind of thick brane system from which it is derived.
The energy density of the model is given by and its shape is depicted in Fig. 11. It is well known that models derived from the first-order equations (45) and (46) have zero energy. It happens because we can rewrite the energy density as a total derivative, ρ(y) = d dy W e 2 A . As W (φ(±∞)) is finite and asymptotically the warp function falls off as e −constant 2 |y| , the integral of ρ(y) over all space must vanish. Here we see that when λ = 0, the energy density is well concentrated around the origin and, as λ increases, it becomes more and more diffuse, finally disappearing at λ = 1.
The above results describe an interesting scenario, in which the parameter λ may be used to control the physical properties of the 5-dimensional spacetime. If λ increases from 0 to unity, it may change the spacetime from a braneworld model with a single extra spatial dimension of Fig. 11 The energy density for some values of λ. Here we see how it becomes delocalized as λ → 1. For λ = 1, there is no more energy concentration around the origin and the brane disappears from the background geometry infinite extent to a 5-dimensional Minkowski spacetime with no graviton localization.

Metric fluctuations
In this section we analyze the stability of the gravitational sector. For this purpose, we perform a redefinition of variable dy 2 → e 2 A(z) dz 2 in (40), which allows us to rewrite the metric in a conformally flat scenariog ab = e 2 A(z) η ab . With a linear perturbation the metric becomes In the transverse-traceless gauge (∂ μ h μν = 0 and h μ μ = 0) the conformal Einstein tensor isḠ ab = − 1 2 ∂ c ∂ c h ab and the linearized Einstein tensor is given by In this way the μν-components of G (1) ab are and the linearized energy-momentum tensor becomes where the prime denotes the derivative with respect to the variable z. When using linearized Einstein equations, G (1) μν , we obtain the equation for h μν , which is −∂ c ∂ c h μν + 3A h μν = 0. At last, the redefinition H μν = e −i px e 3A/2 h μν allows us to rewrite the equation for h μν as Fig. 12 The stability potential (58) for some values of λ. Here we observe that the stability potential of the gravitational sector considerably decreases when λ → 1, although it always maintains its volcano shape, implying stability for gravity localization. For λ = 1, Eq. (58) vanishes Note that Eq. (57) has the form of a supersymmetric Schrödinger equation, where the stability potential is It is depicted in Fig. 12. We can observe that (58) keep its volcano shape, as usual, but as λ grows it becomes less pronounced, and we have no stability potential for λ = 1. Equation (57) has the factorized form S + S − ψ = p 2 ψ, where S ± = ±∂ z + 3A /2 . In particular, the zero-energy solution of (57) is given by ψ 0 (y) = e 3A(y)/2 , where the asymptotic behavior of A(y) ensures gravity localization around the brane for 0 ≤ λ < 1. Since the Hermitian operator S + S − is non-negative, we have no normalizable negative gravitons modes, and it ensures system stability. At λ = 1, we find a flat space, so the stability analysis as presented here fails. In this case, the stability of this space is guaranteed by the positive mass theorems on asymptotically flat spacetimes, which are valid for dimensions ≤ 7 [36,37].

RG flow
As pointed out in [30,31], by rewriting the metric (40) in the form where u = e A(y) , we are allowed to interpret the function u(y) as the renormalization scale of some quantum field theory in the sense of the gauge/gravity duality [26][27][28][29]. Note that, due to the boundary conditions on the deformation factor, u is limited to the range [0, 1] = [e A(±∞) , e A(0) ]. Thus, domain wall solutions such as those presented here naturally lead to confining regimes with a UV-cutoff located at u = 1 in the dual field theory. In this context, the solution φ(y(u)) is identified with the running coupling of the system (see [10,38,39] and the references therein), so that the β-function defined by describes the RG flow in the dual theory.
Denoting the critical points of (47) by φ λ,∞ , we can identify them with the zeros of W φ , which leads to the AdS (or flat) vacua. Expanding the β-function around its critical point we have Since in the critical points we have β(φ = φ λ,∞ ) = 0, we can find the following expression for the running coupling: where φ λ,∞ = sc −1 1 √ 1+λ , λ . Note that (62) is a monotonically increasing function of λ and we have no divergences in the running coupling since u ∈ [0, 1], but for any λ the IR regime (u → 0) at φ = φ λ,∞ is well defined in the dual field theory, even for the 5d Minkowski setup. As a consequence, none of the solutions has a conformal dual field theory.

Comments
In this work we studied a sine-Gordon-like model, which is controlled by a real parameter that continuously connects the sine-Gordon and the vacuumless models. The model appears as a deformation of the φ 4 model, and the real parameter is λ: for λ = 0 one gets the standard sine-Gordon model and for λ = 1 it reproduces the so-called vacuumless model. However, for λ in the interval (0, 1) one gets a double sine-Gordon model, which contains two distinct topological sectors, the large and the small sectors, which give rise to large and small kinks, respectively.
As shown, in the 2-dimensional spacetime, the energy of the large kink in the large sector varies from E L (λ = 0) = 1 to E L (λ = 1) = 2π , and in the case of the small sector one gets E S (λ = 0) = 1 and E S (λ = 1) = 0. We then see that the small sector, which is degenerate with the large sector at λ = 0, disappears as λ = 1, with the large sector becoming the topological sector of the vacuumless model.
In the 5-dimensional case, we considered a warped geometry with a single extra dimension of infinite extent and studied the new braneworld scenario described in the small sector. In this scenario, the brane energy density is such that the brane energy vanishes, independently of the value of λ. If we see the model with λ increasing from zero to unity, it then nicely describes a way to change a 5-dimensional warped geometry which is asymptotically AdS into a 5-dimensional Minkowski geometry. However, if λ is supposed to run in the reverse sense, decreasing from unity to zero, the model could do the reverse, changing the 5-dimensional Minkowski geometry into a braneworld scenario with a warped geometry which is asymptotically AdS 5 .
As we have shown, the model is stable under tensorial fluctuations in the metric and of current interest, so one should now investigate how it modifies Newton's law, and how fermion and gauge fields can be entrapped into the brane as λ varies in the interval [0, 1]. Another issue of current interest concerns the variation of λ: the present investigation cannot tell the value of λ, so one should search for this considering other arguments. An interesting possibility is to investigate the conformational entropy associated with the current braneworld model, and see how it behaves as λ varies in the interval [0, 1]. This has been recently investigated in other contexts in [40][41][42][43] and in the references therein, and it may provide an important guide towards the physical realization of gravity localization in the present model, since the entropy could perhaps suggest a better way λ should vary, increasing or decreasing in the interval [0, 1].